In vivo T1 mapping of the spinal cord using a reduced Field-of-View Inversion Recovery sequence (IR-ZOOM-EPI)
Marco Battiston1, Torben Schneider2, Claudia Angela Michela Gandini Wheeler-Kingshott1,3, and Rebecca S Samson1

1NMR Research Unit, Queen Square MS Centre, Department of Neuroinflammation, UCL Institute of Neurology, University College London, London, United Kingdom, 2Philips Healthcare, Guildford, United Kingdom, 3Brain Connectivity Center, C. Mondino National Neurological Institute, Pavia, Italy

Synopsis

The T1 relaxation time is a fundamental quantitative Magnetic Resonance parameter widely used to characterize healthy and pathological tissue. However, investigation of quantitative T1 in the human spinal cord has been limited to date. In this work, we propose a scan time efficient protocol in the spinal cord for Inversion Recovery T1 mapping, which is considered the “gold-standard” method. The mean (± standard deviation) T1 for white matter and grey matter in the cervical spinal cord were found to be respectively 1096 (±26) ms and 1153 (±24) ms.

Purpose

To develop a reliable method for measuring T1 in the cervical spinal cord in vivo at 3T.

Introduction

The longitudinal relaxation time (T1) is related to macromolecular concentration, water binding and water content [1], and is therefore important for tissue characterisation and assessment of pathology. Furthermore, the accurate knowledge of T1 serves as the basis for several other quantitative MR methods, including in vivo spectroscopy, perfusion imaging and quantitative magnetization transfer imaging.

In the spinal cord, conventional protocols are hampered by low resolution, limited coverage and long scan times, and therefore limited in clinical practise. Currently, the only a study reporting T1 values in vivo at 3T is limited to a single slice at the C2/C3 cord level [2] and requires ~20 minutes of scan time.

We propose a T1 mapping protocol for the spinal cord that addresses spatial coverage and scan time limitations by combining reduced field-of-view (FOV) acquisition [3] with a slice-shuffling Inversion Recovery (IR) approach [4]. We show that our protocol achieves robust T1 mapping of the whole cervical spinal cord in under 9 minutes.

Materials and Methods

Sequence description:

Zonally Magnified Oblique Multi-slice Echo Planar Imaging (ZOOM-EPI) enables the acquisition of a small FOV without aliasing artefacts of surrounding tissue. However, due to the oblique excitation, a TR>>T1 is required between contiguous slice excitations for the longitudinal magnetisation to recover. Multi-slice acquisition therefore requires adjacent slices to be split into multiple packages with a long TR between each of the packages (figure 1a). The intrinsic constraint TR>>T1 can be exploited to perform an IR experiment for T1 measurement within a clinically feasible scan time, i.e. the base sequence is prepended by an inversion pulse followed by different Inversion Times (TI).

We use a non-selective inversion pulse (i.e. the same pulse is experienced by all slices within a package) to acquire data at multiple TI values in a time-efficient way by shuffling the slice acquisition order within the package and varying initial delay following the inversion pulse over different sequence repetitions (figure 1b).

The interplay between the slice shuffling mechanism and short ZOOM-EPI readout times (~40ms) improves scan time efficiency immensely, in particular for large contiguous slice coverage.

Data Acquisition:

To image 15 cm of the human cervical spinal cord, we split a stack of 30 5mm-thick slices into 5 packages, with 12 different TIs. Sequence details: FOV 64x48mm2, 1x1mm2 in-plane resolution, reconstructed to 0.5x0.5mm2, recovery time Trec=6s, TIs=50, 250, 450, 650, 850, 1050, 1250, 1450, 1650, 1850, 2050, 2250ms. Scan time was 8min16sec.

4 healthy subjects (3M, 25-28 years) were scanned using the IR-ZOOM-EPI sequence and 8 repetitions of an identical spin echo ZOOM-EPI to estimate the noise standard deviation (SD).

Data Analysis:

All data were registered slice-wise to the first TI volume with a 3 degrees-of-freedom model using FLIRT. A mono-exponential model was fitted to magnitude data using maximum likelihood estimation under the assumption of Rician noise. Mean T1 values for Grey (GM) and White Matter (WM) were extracted from regions of interest (ROIs) manually drawn on the first TI volume and applied to voxel-wise parametric maps.

Results

Figure 2a shows the 12 images acquired at different TIs in a representative slice. In figure 2b, mean WM and GM ROI signal is plotted against predicted signal when voxel-wise parameter estimates are averaged within ROIs. Figure 3 shows T1 maps at different levels of the cervical spinal cord for a single subject.

Mean T1 (±SD) for both WM and GM in all subjects are reported in figure 4. The mean (± SD) T1 for WM and GM was found to be respectively 1096 (±26) ms and 1153 (±24) ms, giving coefficients-of-variation (COVs) of 2.36% and 2.09%.

Discussion and Conclusions

We have proposed a new scan time efficient T1 mapping protocol for large coverage in the spinal cord. We demonstrate feasibility of T1 relaxation mapping over the whole cervical cord in 9 minutes using the standard IR approach with excellent inter-subject reproducibility.

Compared to previous findings in the spinal cord at 3T [2], we found higher values for both WM and GM, closer to the ranges usually reported for the brain [5], and reduced differentiation between the tissue types. However, differences in the acquisition protocols used are known to have an impact on T1 estimation [6].

Due to the short acquisition time it could easily be added to a routine protocol for spinal cord T1 quantification. In future work, we will investigate whether the acquisition time can be further decreased, e.g. by reducing the number of TIs for T1 estimation.

Acknowledgements

The UK MS Society and the UCL-UCLH Biomedical Research Centre for ongoing support.

References

[1] Gowland P. T1: the Longitudinal Relaxation Time. Quantitative MRI of the brain: Measuring changes caused by disease (2003): 111.

[2] Smith S A. Measurement of T1 and T2 in the cervical spinal cord at 3 Tesla. MRM 60:213-210(2008).

[3] Wheeler-Kingshott C AM. Investigating cervical spinal cord structure using axial diffusion tensor imaging. Neuroimage 16.1(2002):93-102.

[4] Clare S. Rapid T1 mapping using multislice echo planar imaging. MRM 45:630-634(2001).

[5] Bottomley P A. A review of normal tissue hydrogen NMR relaxation times and relaxation mechanisms from 1-100MHz: dependence on tissue type, NMR frequency, temperature, species, excision and age. Medical Physics 11.4(1984):425-448.

[6] Stikov N. On the accuracy of T1 mapping: searching for common ground. MRM 73:514-522(2015).

Figures

(a) Interleaved ZOOM-EPI acquisition. Slices belonging to the same package (identified by colours) are acquired sequentially. (b) Slice order is shuffled within packages to avoid TR lengthening in a multi-TI experiment. To ensure full recovery of magnetisation between inversion pulses of different packages, a delay Trec (>5s) is added.

(a) Images of the 12 sampled TIs during magnetisation recovery following inversion for a representative slice. (b) Left: WM (red) and GM (blue) ROIs; Right: corresponding mean signal plotted against model prediction given by mean parameters within the same ROIs.

Average images of the repeated spin echo ZOOM-EPI acquisition and corresponding T1 maps for a representative subject at different level of the spinal cord.

Mean and SD of T1 in white matter and grey matter for all subjects, with relative COVs.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
4390