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 T
1 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).