1454

Making myelin water imaging “normal”: Robust acquisition for cervical cord with low peripheral nerve stimulation
Sharada Balaji1, Adam Dvorak1, Neale Wiley1, Laura Barlow2, Irene M. Vavasour3,4, Alex MacKay1,3, and Shannon H. Kolind1,3,4,5
1Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada, 2MRI Research Centre, University of British Columbia, Vancouver, BC, Canada, 3Radiology, University of British Columbia, Vancouver, BC, Canada, 4International Collaboration on Repair Discoveries, Vancouver, BC, Canada, 5Medicine, University of British Columbia, Vancouver, BC, Canada

Synopsis

The standard myelin water imaging (MWI) sequence for cervical spinal cord, 3D gradient and spin echo (GRASE), was modified to allow scanning in “normal mode” (limiting imaging parameters that may cause physiologic stress) for patients who may not tolerate peripheral nerve stimulation or tissue heating due to MR-conditional implants or medical issues. Traditional 32-echo GRASE was replaced with 48-echo GRASE in normal mode. Myelin water fraction maps from the new sequence had better repeatability than the standard sequence. 48-echo GRASE is the recommended sequence for MWI in spinal cord, particularly for subjects who cannot be scanned outside of normal mode.

Introduction

Myelin water imaging (MWI) is a T2 relaxation-based quantitative biomarker for myelin content in the central nervous system1. MWI decomposes the T2 decay into components corresponding to myelin-associated water (<40ms), intra/extra-cellular water (40-200ms) and cerebrospinal fluid (~2s). Myelin water fraction (MWF) is calculated as the ratio of the signal from myelin-associated water to the total signal.

A commonly used technique for MWI of cervical spinal cord is an 8.5 minute 32-echo 3D Gradient and Spin Echo (3D GRASE) sequence2 operating in first-level controlled mode, which allows for faster gradient switching. This mode is not suitable for scanning subjects with MR-conditional implants or susceptible physiological conditions (e.g. poor cardiac output, opioid drug use, pregnancy) as it may cause overheating of subjects or implants due to RF pulses.

Normal mode3,4 limits the RF power deposition and gradient switching rate to reduce the chance of overheating. The peripheral nerve stimulation (PNS) threshold is limited to 80%3,4 by restricting the time-varying magnetic field gradients (dB/dt). The specific absorption rate (SAR), a measure of RF power deposition per unit mass, and B1+ rms, the time average of the applied RF field, are also limited. To prevent the core body temperature from rising by more than 0.5C, normal mode requires that whole body SAR<2W/kg3–5. Here, a modified 3D GRASE sequence was implemented and evaluated in normal mode.

Methods

Acquisition: A healthy volunteer underwent scanning on a Philips Elition 3T scanner using a 16-channel head and neck coil, with scans centered at the C2/C3 level of the spinal cord. The scanning protocol consisted of:

1. Multi-echo gradient echo (GRE-ME): resolution 0.5x0.5x5mm3, TR/TE1/ΔTE=600/7.1/6.5ms, Compressed Sense (CS) factor=2, 2 signal averages, scan duration=4.5min, average of 3 acquired echoes used for anatomical segmentation and registration
2. Axial 32-echo 3D GRASE in first-level controlled mode: 32 echoes, TE=10ms, TR=1500ms, refocusing flip angle=180°, acquired resolution=0.75x0.75x5mm3, FOV=180x150x40mm3 (AP, RL, FH), reconstructed resolution=0.6x0.6x2.5mm3, SENSE factor 2 (RL), scan duration=8.3min
3. Axial 48-echo 3D GRASE in normal mode: 48 echoes, TE=8ms, TR=2000ms, refocusing flip angle=180°, acquired resolution=1.4x0.9x5mm3, FOV=170x170x40mm3 (AP, RL, FH), reconstructed resolution=0.5x0.5x2.5mm3, SENSE factor 2 (RL), scan duration=10.8min
To assess repeatability, both GRASE sequences were collected twice during the scan session without repeating any scan twice in a row. The volunteer stayed in the scanner for the full duration of the exam, without re-positioning.

Analysis: MWF maps were generated using the Julia-based DECAES software6, including corrections for stimulated echo artifacts and regularized non-negative least squares (NNLS) fitting7. Images were segmented and GRE-ME was registered with the PAM50 template using Spinal Cord Toolbox8. Regions of interest (ROIs) were obtained by registering the GRASE images with the PAM50 template aligned in GRE-ME space, then concatenating the two warp fields to transform data from GRASE space to PAM50 template space using a single interpolation step. Regionwise MWFs were extracted for six ROIs (whole spinal cord, white matter, grey matter, dorsal columns, lateral funiculi, ventral funiculi). MWFs from each scan were compared to assess repeatability.

Results

Figure 1 shows the generated MWF maps for the same slice from each scan.

Table 1 gives a summary of MWFs separated by region. The average scan-rescan coefficient of variation9 for the 32-echo sequence was 14%, and for the 48-echo sequence was 11%. The average MWF WM/GM contrast for the 32-echo sequence was 1.3 and for the 48-echo sequence was 1.5.

Table 2 compares the whole body SAR, B1+ rms, PNS threshold as a percentage of the maximum, and dB/dt (time varying gradients) from each sequence.

The fit-to-noise ratio (FNR) for the new 48-echo sequence was on average 60% higher than for the 32-echo sequence.

Discussion

The 48-echo sequence has a slightly longer acquisition duration due to the requirement of lower power deposition as well as more RF pulses, which resulted in a longer TR.

Scan parameters relating to PNS were reduced in normal mode as expected. Although SAR was well within normal mode operating limits in both sequences, a slight increase for the 48-echo sequence was noted, likely due to more RF pulses. PNS is the limiting factor for scanning subjects who would otherwise be unable to be scanned using GRASE MWI, because the chain of gradient-echo/spin-echo/gradient-echo (for an echo-planar imaging factor of three) must be acquired in a brief timespan to provide the short echo spacing required for MWI, which requires high dB/dt rates.

The 48-echo sequence compared favourably with the traditional 32-echo sequence in scan-rescan repeatability. Higher fit-to-noise ratios were seen with the 48-echo sequence, indicating better fits for the T2 decay curves. This may be due to the ability of the 48-echo sequence to better fit slower decaying signal because of the longer echo train. The shorter echo spacing on the new sequence should be able to better quantify the fast-relaxing myelin water compartment.

Conclusion

Modifications were made to the current standard for cervical spinal cord MWI, the 32-echo 3D GRASE sequence, to operate in normal mode and be compatible for participants who would otherwise be unable to be scanned safely. The new 48-echo sequence features shorter echo spacing, better fits for the T2 decay curves, similar MWF values, and better repeatability than the standard sequence.

Acknowledgements

We thank the researchers, volunteers and the MR technical staff at the UBC MRI Research Centre for their support. This research was funded by Michael Smith Health Care BC (SK), NSERC (SB, AD, SK Grant RGPIN-2018-03904).

References

1. MacKay A, Whittall K, Adler J, Li D, Paty D, Graeb D. In vivo visualization of myelin water in brain by magnetic resonance. Magnetic resonance in medicine. 1994;31(6):673-677. http://doi.wiley.com/10.1002/mrm.1910310614

2. Ljungberg E, Vavasour I, Tam R, et al. Rapid myelin water imaging in human cervical spinal cord. Magnetic Resonance in Medicine. 2017;78(4). doi:10.1002/mrm.26551

3. Loren A. Zaremba. FDA Guidelines for Magnetic Resonance Equipment Safety. https://www.aapm.org/meetings/02AM/pdf/8356-48054.pdf

4. Medicines and Healthcare Products Regulatory Agency. Safety Guidelines for Magnetic Resonance Imaging Equipment in Clinical Use.; 2021.

5. Allison J, Yanasak N. What MRI Sequences Produce the Highest Specific Absorption Rate (SAR), and Is There Something We Should Be Doing to Reduce the SAR During Standard Examinations? American Journal of Roentgenology. 2015;205(2). doi:10.2214/AJR.14.14173

6. Doucette J, Kames C, Rauscher A. DECAES – DEcomposition and Component Analysis of Exponential Signals. Zeitschrift für Medizinische Physik. 2020;30(4). doi:10.1016/j.zemedi.2020.04.001

7. Prasloski T, Mädler B, Xiang QS, MacKay A, Jones C. Applications of stimulated echo correction to multicomponent T2 analysis. Magnetic resonance in medicine. 2011;67(6):1803-1814. http://doi.wiley.com/10.1002/mrm.23157

8. de Leener B, Lévy S, Dupont SM, et al. SCT: Spinal Cord Toolbox, an open-source software for processing spinal cord MRI data. NeuroImage. 2017;145. doi:10.1016/j.neuroimage.2016.10.009

9. Meyers SM, Vavasour IM, Mädler B, et al. Multicenter measurements of myelin water fraction and geometric mean T 2 : Intra- and intersite reproducibility. Journal of Magnetic Resonance Imaging. 2013;38(6). doi:10.1002/jmri.24106

Figures

Figure 1. MWF maps from the (a) two 32-echo scans and (b) two 48-echo scans at the same slice location.

Table 1. Regionwise mean MWF values (standard deviation) from each scan-rescan test, and the scan-rescan coefficient of variation (CoV) accounting for only 2 samples

Table 2. Scan parameters from both sequences

Proc. Intl. Soc. Mag. Reson. Med. 30 (2022)
1454
DOI: https://doi.org/10.58530/2022/1454