Reduced Field-of-View Diffusion-Weighted Imaging of the Lumbosacral Enlargement: A Pilot In Vivo study of the healthy spinal cord using a clinical 3T MR system
Marios C Yiannakas1, Polymnia Louka1, Francesco Grussu1, Ferran Prados1,2, Rebecca S Samson1, Marco Battiston1, Sebastien Ourselin2, David H Miller1,3, and Claudia Angela Michela Gandini Wheeler-Kingshott1,4

1NMR Research Unit, Queen Square MS Centre, Department of Neuroinflammation, UCL Institute of Neurology, University College London, London, United Kingdom, 2Translational Imaging Group, Medical Physics and Biomedical Engineering, University College London, London, United Kingdom, 3NIHR Biomedical Research Centre, UCL-UCLH, London, United Kingdom, 4Brain Connectivity Center, C. Mondino National Neurological Institute, Pavia, Italy

Synopsis

The use of imaging methods to study the lower spinal cord has been hindered by a number of technical challenges; hence the relative contribution of pathology affecting the lower spinal cord to the observed clinical symptoms remains largely unexplored. In this pilot study, we investigate the feasibility of obtaining tissue-specific (grey matter and white matter) diffusion tensor imaging metrics within the lumbosacral enlargement in vivo in healthy volunteers using reduced field-of-view echo planar imaging on a clinical 3T MRI system. Preliminary results show that such measures may be obtained reliably and within clinically acceptable scan times.

INTRODUCTION:

Imaging methods that can be used to study the lower spinal cord have potential to provide new insights in to the pathophysiology of symptoms such as bladder and sexual dysfunction, which are often associated with neurological conditions including multiple sclerosis1, spinal cord injury2 and multiple system atrophy3. Using magnetic resonance imaging (MRI), it has recently been possible to depict grey matter (GM) and white matter (WM) in the lumbosacral enlargement (LSE) using clinical systems4. Quantitative MRI methods such as diffusion tensor imaging (DTI) could provide measures that reflect tissue microstructure at that level, hitherto unreported. In this pilot study, we investigated the feasibility of obtaining tissue-specific (i.e. GM and WM) DTI metrics within the LSE in vivo in healthy volunteers using reduced field-of-view echo planar imaging (ZOOM-EPI5,6) on a clinical 3T MRI system.

METHOD:

A) Participants: Six healthy subjects were recruited (mean age 27 years, 5 female). The local institutional review board approved the study and informed consent was obtained from all participants. B) MR Imaging: A 3T Philips Achieva system with RF dual-transmit technology (Philips Healthcare, Best, Netherlands) and the product 15-channel SENSE spine coil were used. A sagittal T2-weighted image of the lumbar spine was first obtained and used to facilitate prescription of subsequent scans perpendicular to the cord at the T11 - L1 level, to ensure coverage of the LSE4. For identifying the LSE, a 3D fast field-echo (3D-FFE) sequence was used with fat suppression: TR = 22 ms, TE = 4.4 ms, flip angle a =10°, FOV = 180 x 180 mm2, voxel size = 0.5 x 0.5 x 5 mm3, NEX = 8, slices = 10, acquisition time ~ 9, 5 mins. For calculating the DTI metrics, a cardiac-gated ZOOM-EPI was used with identical slice geometry as the 3D-FFE, 60 diffusion directions at b = 1000 s/mm2 interleaved with 7 b = 0 measurements: TR = 6000 ms, TE = 40 ms, flip angle a = 90°, number of slices = 10; FOV = 64 x 48 mm2, voxel size = 1 x 1 x 5 mm3, acquisition time ~ 15 mins (depending on heart rate). C) Image analysis: Using the 3D-FFE, three slices (i.e. 15 mm) through the widest section of the lumbar cord (i.e. the LSE) were identified, using the active surface model (ASM) segmentation method in JIM 6.0 (http://www.xinapse.com)4. The diffusion-weighted (DW) images were firstly corrected for motion using slice-wise linear registration implemented in FSL (http://www.fmrib.ox.ac.uk/fsl/), with registration transformations estimated among non-DW images7. DTI fitting was performed using CAMINO (http://cmic.cs.ucl.ac.uk/camino/) and maps of axial, radial and mean diffusivity (AD/RD/MD) and fractional anisotropy (FA) were obtained. Using the mean b = 0 volume, the three slices corresponding to the LSE, previously identified from the 3D-FFE, were segmented using ASM to obtain the whole cord outline. GM was manually outlined on an image obtained by averaging DW images7,8. Binary masks were subsequently created and eroded prior to their application to the DTI maps. To examine the impact of the number of diffusion directions on the quality of DTI indices, the DTI model was retrospectively fitted to gradually reduced sets of measurements optimally spread over the sphere (60 to 10, with steps of 10), extracted with CAMINO. D) Reproducibility assessment: All volunteers had a repeated scan on a different occasion, and the reproducibility of all DTI metrics obtained with all diffusion protocols within GM, WM and the whole cord was assessed by calculating the intra-class correlation coefficient (ICC) and the coefficient of variation (%COV).

RESULTS:

High-resolution images through the LSE and corresponding DTI maps are shown in Figure 1. Tissue-specific values (mean ± SD) of AD, RD, FA and MD in 6 healthy subjects are shown in Table 1. ICC and %COV results from the reproducibility assessment are shown in Table 2 and Table 3, respectively. The effect of using a different number of diffusion directions on the %COV of the DTI metrics within each tissue-type is demonstrated in Figure 2.

DISCUSSION AND CONCLUSION:

This study has shown that tissue-specific DTI metrics within the LSE can be obtained reliably using a commercially available 3T MR system. This allows characterisation of microstructural features at lumbar level, which show some differences compared to other levels, such as lower FA in WM than in the upper cervical spinal cord. Future studies will be focused on developing the most time-efficient acquisition protocol (i.e. minimum number of diffusion directions), refining the image segmentation method further and assessing the value of the final protocol in the investigation of neurological disease.

Acknowledgements

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

References

1) McCombe PA, Gordon TP, Jackson MW. Bladder dysfunction in multiple sclerosis. Expert Rev Neurother. 2009; 9: 331-340 2) Anderson KD, Borisoff JF, Johnson RD et al. The impact of spinal cord injury on sexual function: concerns of the general population. Spinal Cord 2007; 45: 328-337 3) Sakakibara R, Hattori T, Tojo M et al. Micturitional disturbance in multiple system atrophy. Jpn J Psychiatry Neurol. 1993; 47: 591-598. 4) Yiannakas MC, Kakar P, Hoy LR et al. The Use of the Lumbosacral Enlargement as an Intrinsic Imaging Biomarker: Feasibility of Grey Matter and White Matter Cross-Sectional Area Measurements Using MRI at 3T. PLoS One 2014; 9: e105544, 5) Wilm BJ, Gamper U, Henning A et al. Diffusion-weighted imaging of the entire spinal cord. NMR Biomed. 2009; 22:174-81, 6) Wheeler-Kingshott CA, Hickman SJ, Parker GJ et al. Investigating cervical spinal cord structure using axial diffusion tensor imaging. NeuroImage 2002; 16:93–102, 7) Grussu F, Schneider T, Zhang H et al. Neurite orientation dispersion and density imaging of the healthy cervical cord in vivo. Neuroimage 2015; 111: 590-601 8) Kearney H, Schneider T, Yiannakas MC et al. Spinal cord grey matter abnormalities are associated with secondary progression and physical disability in multiple sclerosis. J Neurol Neurosurg Psychiatry. 2015 Jun; 86: 608-14.

Figures

Table 1. DTI metrics obtained within the lumbosacral enlargement (LSE); mean (SD) values in grey matter (GM), white matter (WM) and the whole cord across six healthy volunteers.

Table 2. The intra-class correlation coefficient (ICC) of the DTI metrics is reported separately for grey matter (GM), white matter (WM) and the whole cord. Here the 1-ICC is reported instead, which provides an estimate of the fraction of variability due to measurement error (within-subject) over the total variation.

Table 3. The percentage coefficient of variation (%COV) of the DTI metrics is reported separately for grey matter (GM), white matter (WM) and whole cord.

Figure 1. High resolution image of the lumbosacral enlargement acquired with the fast field echo (FFE) sequence and corresponding maps of axial diffusivity (AD), radial diffusivity (RD), fractional anisotropy (FA) and mean diffusivity (MD).

Figure 2. The effect of acquiring data using a different number of diffusion directions on the percentage coefficient of variation (%COV) of the DTI metrics, here plotted separately for grey matter (GM), white matter (WM) and the whole cord.



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