Optimized multi-echo gradient-echo magnetic resonance imaging for gray and white matter segmentation in the lumbosacral cord at 3 T

Silvan Büeler, Marios C Yiannakas, Zdravko Damjanovski, Patrick Freund, Martina D Liechti, Gergely David, Silvan Büeler, Marios C Yiannakas, Zdravko Damjanovski, Patrick Freund, Martina D Liechti, Gergely David

Abstract

Atrophy in the spinal cord (SC), gray (GM) and white matter (WM) is typically measured in-vivo by image segmentation on multi-echo gradient-echo magnetic resonance images. The aim of this study was to establish an acquisition and analysis protocol for optimal SC and GM segmentation in the lumbosacral cord at 3 T. Ten healthy volunteers underwent imaging of the lumbosacral cord using a 3D spoiled multi-echo gradient-echo sequence (Siemens FLASH, with 5 echoes and 8 repetitions) on a Siemens Prisma 3 T scanner. Optimal numbers of successive echoes and signal averages were investigated comparing signal-to-noise (SNR) and contrast-to-noise ratio (CNR) values as well as qualitative ratings for segmentability by experts. The combination of 5 successive echoes yielded the highest CNR between WM and cerebrospinal fluid and the highest rating for SC segmentability. The combination of 3 and 4 successive echoes yielded the highest CNR between GM and WM and the highest rating for GM segmentability in the lumbosacral enlargement and conus medullaris, respectively. For segmenting the SC and GM in the same image, we suggest combining 3 successive echoes. For SC or GM segmentation only, we recommend combining 5 or 3 successive echoes, respectively. Six signal averages yielded good contrast for reliable SC and GM segmentation in all subjects. Clinical applications could benefit from these recommendations as they allow for accurate SC and GM segmentation in the lumbosacral cord.

Conflict of interest statement

The authors declare no competing interests.

© 2022. The Author(s).

Figures

Figure 1
Figure 1
(A) Sagittal T2-weighted turbo spin echo acquisition in the lower spine used for subsequent prescription of the high-resolution axial acquisition. (B) Corresponding axial slices acquired with the 3D multi-echo gradient-echo sequence ( Siemens FLASH) in the caudal-rostral direction (slices 1–20). Highlighted are the slice in the lumbosacral enlargement (LSE) with the largest cord cross-sectional spinal cord area (defined as the "LSE slice" and shown in light blue in A and B; here: slice 15), and the most caudal slice in the conus medullaris (CM) where the gray matter still has the characteristic butterfly shape (defined as the "CM slice" and shown in red in A and B; here: slice 9). A saturation band, displayed as yellow shaded area in A, was placed anterior to the spine to suppress signal and possible artifacts arising from abdominal peristalsis.
Figure 2
Figure 2
Visual representation of echoes, echo combinations, signal averages, and image segmentations. The 3D spoiled multi-echo gradient echo sequence (Siemens FLASH) consisted of 5 echoes and was acquired with 8 individual repetitions. For each subject, a series of images was created by successively combining echoes (echo 1, 1–2, 1–3, 1–4, 1–5) and averaging across repetitions (number of signal averages (NSA): 1, 2, 3, …, 8), resulting in a total of 72 images. (A) Image series of individual echoes (NSA = 8) for a representative slice in the lumbosacral enlargement (LSE). (B) Image series of increasing number of combined echoes in the same slice as in (A) (NSA = 8). (C) Image series with increasing NSA (3 combined echoes). (D) Spinal cord (SC) and gray matter (GM) were segmented manually in each slice (here a representative LSE and conus medullaris slice are shown). A mask of cerebrospinal fluid (CSF) was drawn anterior to the SC. White matter (WM) mask was obtained by subtracting GM from the SC mask.
Figure 3
Figure 3
Quantitative comparison of the individual and combined echoes (8 signal averages) in terms of contrast-to-noise ratio (CNR) and contrast between white matter (WM) and cerebrospinal fluid (CSF) and betwen gray matter (GM) and WM, and signal-to-noise ratio (SNR) of WM and GM. Values are displayed for individual echoes (1, 2, 3, 4, and 5) and combined echoes (1–2, 1–3, 1–4 and 1–5). All measures are displayed separately for the lumbosacral enlargement (LSE) in blue and conus medullaris (CM) in red. In all subplots, error bars represent standard deviation across participants (n = 10).
Figure 4
Figure 4
Quantitative comparison of the number of signal averages (NSA) (3 combined echoes) in terms of contrast-to-noise ratio (CNR) and contrast between white matter (WM) and cerebrospinal fluid (CSF) and between gray matter (GM) and WM, and signal-to-noise ratio (SNR) of WM and GM. All measures are displayed separately for the lumbosacral enlargement (LSE) in blue and conus medullaris (CM) in red. In all subplots, error bars represent standard deviation across participants (n = 10).
Figure 5
Figure 5
Qualitative comparison of the individual and combined echoes (8 signal averages). Shown are ranking scores assigned to individual echoes (1, 2, 3, 4, 5) and combined echoes (1–2, 1–3, 1–4, 1–5) and averaged across all participants (n = 10) and raters (n = 5). Echoes/echo combinations were ranked separately for spinal cord (SC) only, gray matter (GM) only and joint SC and GM segmentability. The analysis was performed in a representative slice at the lumbosacral enlargement (LSE) in blue and conus medullaris (CM) in red. Ranking score ranges from 1 (lowest score) to 9 (highest score). Error bars represent standard deviation across participants.

References

    1. Anderson KD, Borisoff JF, Johnson RD, Stiens SA, Elliott SL. The impact of spinal cord injury on sexual function: Concerns of the general population. Spinal Cord. 2007;45:328–337. doi: 10.1038/sj.sc.3101977.
    1. Fernández O. Mechanisms and current treatments of urogenital dysfunction in multiple sclerosis. J. Neurol. 2002;249:1–8. doi: 10.1007/PL00007835.
    1. McCombe PA, Gordon TP, Jackson MW. Bladder dysfunction in multiple sclerosis. Expert Rev. Neurother. 2009;9:331–340. doi: 10.1586/14737175.9.3.331.
    1. Panicker JN, et al. Early presentation of urinary retention in multiple system atrophy: can the disease begin in the sacral spinal cord? J. Neurol. 2020;267:659–664. doi: 10.1007/s00415-019-09597-2.
    1. Watanabe T, et al. High incidence of occult neurogenic bladder dysfunction in neurologically intact patients with thoracolumbar spinal injuries. J. Urol. 1998;159:965–968. doi: 10.1016/S0022-5347(01)63786-8.
    1. Yokota K, et al. Pathological changes of distal motor neurons after complete spinal cord injury. Mol. Brain. 2019;12:1–15. doi: 10.1186/s13041-018-0422-3.
    1. Petrova N, Carassiti D, Altmann DR, Baker D, Schmierer K. Axonal loss in the multiple sclerosis spinal cord revisited. Brain Pathol. 2018;28:334–348. doi: 10.1111/bpa.12516.
    1. Konno H, Yamamoto T, Iwasaki Y, Iizuka H. Shy-Drager syndrome and amyotrophic lateral sclerosis: Cytoarchitectonic and morphometric studies of sacral autonomic neurons. J. Neurol. Sci. 1986;73:193–204. doi: 10.1016/0022-510X(86)90130-9.
    1. Cohen-Adad J, et al. Involvement of spinal sensory pathway in ALS and specificity of cord atrophy to lower motor neuron degeneration. Amyotroph. Lateral Scler. Front. Degener. 2013;14:30–38. doi: 10.3109/17482968.2012.701308.
    1. Cohen-Adad J, et al. Demyelination and degeneration in the injured human spinal cord detected with diffusion and magnetization transfer MRI. Neuroimage. 2011;55:1024–1033. doi: 10.1016/j.neuroimage.2010.11.089.
    1. Freund P, et al. Disability, atrophy and cortical reorganization following spinal cord injury. Brain. 2011;134:1610–1622. doi: 10.1093/brain/awr093.
    1. Freund P, et al. MRI investigation of the sensorimotor cortex and the corticospinal tract after acute spinal cord injury: A prospective longitudinal study. Lancet. Neurol. 2013;12:873–881. doi: 10.1016/S1474-4422(13)70146-7.
    1. Lundell H, et al. Independent spinal cord atrophy measures correlate to motor and sensory deficits in individuals with spinal cord injury. Spinal Cord. 2011;49:70–75. doi: 10.1038/sc.2010.87.
    1. Mann RS, Constantinescu CS, Tench CR. Upper cervical spinal cord cross-sectional area in relapsing remitting multiple sclerosis: Application of a new technique for measuring cross-sectional area on magnetic resonance images. J. Magn. Reson. Imaging. 2007;26:61–65. doi: 10.1002/jmri.20959.
    1. Rocca MA, et al. A multicenter assessment of cervical cord atrophy among MS clinical phenotypes. Neurology. 2011;76:2096–2102. doi: 10.1212/WNL.0b013e31821f46b8.
    1. Horsfield MA, et al. Rapid semi-automatic segmentation of the spinal cord from magnetic resonance images: Application in multiple sclerosis. Neuroimage. 2010;50:446–455. doi: 10.1016/j.neuroimage.2009.12.121.
    1. Barry RL, Smith SA. Measurement of T2* in the human spinal cord at 3 T. Magn. Reson. Med. 2019;82:743–748. doi: 10.1002/mrm.27755.
    1. Martin N, et al. Comparison of MERGE and axial T2-weighted fast spin-echo sequences for detection of multiple sclerosis lesions in the cervical spinal cord. Am. J. Roentgenol. 2012;199:157–162. doi: 10.2214/AJR.11.7039.
    1. Yiannakas MC, et al. Feasibility of grey matter and white matter segmentation of the upper cervical cord in vivo: A pilot study with application to magnetisation transfer measurements. Neuroimage. 2012;63:1054–1059. doi: 10.1016/j.neuroimage.2012.07.048.
    1. Asiri A, et al. Comparison between 2D and 3D MEDIC for human cervical spinal cord MRI at 3 T. J. Med. Radiat. Sci. 2021;68:4–12. doi: 10.1002/jmrs.433.
    1. White ML, Zhang Y, Healey K. Cervical spinal cord multiple sclerosis: Evaluation with 2D multi-echo recombined gradient echo MR imaging. J. Spinal Cord Med. 2011;34:93–98. doi: 10.1179/107902610X12911165975025.
    1. Huber E, et al. Dorsal and ventral horn atrophy is associated with clinical outcome after spinal cord injury. Neurology. 2018;90:E1510–E1522. doi: 10.1212/WNL.0000000000005361.
    1. Grabher P, et al. Voxel-based analysis of grey and white matter degeneration in cervical spondylotic myelopathy. Sci. Rep. 2016;6:1–10. doi: 10.1038/srep24636.
    1. Paquin M-EE, et al. Spinal cord gray matter atrophy in amyotrophic lateral sclerosis. Am. J. Neuroradiol. 2018;39:184–192. doi: 10.3174/ajnr.A5427.
    1. Yiannakas MC, Kakar P, Hoy LR, Miller DH, Wheeler-Kingshott CAM. 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 3 T. PLoS ONE. 2014;9:105544. doi: 10.1371/journal.pone.0105544.
    1. Yiannakas MC, et al. Gray vs. white matter segmentation of the conus medullaris: Reliability and variability in healthy volunteers. J. Neuroimaging. 2019;29:410–417. doi: 10.1111/jon.12591.
    1. David G, et al. In vivo evidence of remote neural degeneration in the lumbar enlargement after cervical injury. Neurology. 2019;92:e1367–e1377. doi: 10.1212/WNL.0000000000007137.
    1. David G, et al. Extent of cord pathology in the lumbosacral enlargement in non-traumatic versus traumatic spinal cord injury. J. Neurotrauma. 2022 doi: 10.1089/neu.2021.0389.
    1. David G, et al. Longitudinal changes of spinal cord grey and white matter following spinal cord injury. J. Neurol. Neurosurg. Psychiatry. 2021;1:1–9.
    1. Vallotton K, et al. Tracking white and gray matter degeneration along the spinal cord axis in degenerative cervical myelopathy. J. Neurotrauma. 2021 doi: 10.1089/neu.2021.0148.
    1. Sijbers J, Scheunders P, Bonnet N, Van Dyck D, Raman E. Quantification and improvement of the signal-to-noise ratio in a magnetic resonance image acquisition procedure. Magn. Reson. Imaging. 1996;14:1157–1163. doi: 10.1016/S0730-725X(96)00219-6.
    1. Cohen-Adad J, et al. Generic acquisition protocol and open-access data for quantitative MRI of the spinal cord. Nat. Protoc. 2021 doi: 10.1038/S41596-021-00588-0.
    1. Mohammadi S, Möller HE, Kugel H, Müller DK, Deppe M. Correcting eddy current and motion effects by affine whole-brain registrations: Evaluation of three-dimensional distortions and comparison with slicewise correction. Magn. Reson. Med. 2010;64:1047–1056. doi: 10.1002/mrm.22501.
    1. Smith SA, Edden RAE, Farrell JAD, Barker PB, Van Zijl PCM. Measurement of T1 and T2 in the cervical spinal cord at 3 Tesla. Magn. Reson. Med. 2008;60:213–219. doi: 10.1002/mrm.21596.
    1. Brown RW, Cheng Y-CN, Haacke EM, Thompson MR, Venkatesan R. Magnetic Resonance Imaging: Physical Principles and Sequence Design. 2. Wiley; 2014. pp. 447–510.
    1. Cohen-Adad J, et al. Comparison of multicenter MRI protocols for visualizing the spinal cord gray matter. Magn. Reson. Med. 2022;88:849–859. doi: 10.1002/mrm.29249.
    1. Ashburner J, Ridgway G. Symmetric diffeomorphic modeling of longitudinal structural MRI. Front. Neurosci. 2013;6:197. doi: 10.3389/fnins.2012.00197.
    1. Drake-Pérez M, et al. Normal values of magnetic relaxation parameters of spine components with the synthetic MRI sequence. Am. J. Neuroradiol. 2018;39:788–795. doi: 10.3174/ajnr.A5566.
    1. Islam H, et al. Dynamic per slice shimming for simultaneous brain and spinal cord fMRI. Magn. Reson. Med. 2019;81:825–838. doi: 10.1002/mrm.27388.
    1. Finsterbusch J, Eippert F, Büchel C. Single, slice-specific z-shim gradient pulses improve T2*-weighted imaging of the spinal cord. Neuroimage. 2012;59:2307–2315. doi: 10.1016/j.neuroimage.2011.09.038.
    1. Topfer R, Foias A, Stikov N, Cohen-Adad J. Real-time correction of respiration-induced distortions in the human spinal cord using a 24-channel shim array. Magn. Reson. Med. 2018;80:935–946. doi: 10.1002/mrm.27089.
    1. Verma T, Cohen-Adad J. Effect of respiration on the B0 field in the human spinal cord at 3 T. Magn. Reson. Med. 2014;72(6):1629–1636. doi: 10.1002/mrm.25075.
    1. Vannesjo SJ, Miller KL, Clare S, Tracey I. Spatiotemporal characterization of breathing-induced B0 field fluctuations in the cervical spinal cord at 7T. Neuroimage. 2018;167:191–202. doi: 10.1016/j.neuroimage.2017.11.031.
    1. Vannesjo SJ, Clare S, Kasper L, Tracey I, Miller KL. A method for correcting breathing-induced field fluctuations in T2*-weighted spinal cord imaging using a respiratory trace. Magn. Reson. Med. 2019;81:3745–3753. doi: 10.1002/mrm.27664.
    1. Magnotta, V. A., Friedman, L. & BIRN, F. Measurement of Signal-to-Noise and Contrast-to-Noise in the fBIRN Multicenter Imaging Study. J. Digit. Imaging19, 140–147 (2006).
    1. Papinutto N, Henry RG. Evaluation of intra- and interscanner reliability of MRI protocols for spinal cord gray matter and total cross-sectional area measurements. J. Magn. Reson. Imaging. 2019;49:1078–1090. doi: 10.1002/jmri.26269.
    1. Prados F, et al. Spinal cord grey matter segmentation challenge. Neuroimage. 2017;152:312–329. doi: 10.1016/j.neuroimage.2017.03.010.
    1. Yiannakas MC, et al. Reduced field-of-view diffusion-weighted imaging of the lumbosacral enlargement: A pilot in vivo study of the healthy spinal cord at 3t. PLoS ONE. 2016;11:e0164890. doi: 10.1371/journal.pone.0164890.
    1. Grussu F, Schneider T, Zhang H, Alexander DC, Wheeler-Kingshott CAM. Neurite orientation dispersion and density imaging of the healthy cervical spinal cord in vivo. Neuroimage. 2015;111:590–601. doi: 10.1016/j.neuroimage.2015.01.045.
    1. Gringel T, et al. Optimized high-resolution mapping of magnetization transfer (MT) at 3 Tesla for direct visualization of substructures of the human thalamus in clinically feasible measurement time. J. Magn. Reson. Imaging. 2009;29:1285–1292. doi: 10.1002/jmri.21756.

Source: PubMed

3
Tilaa