Whole brain in-vivo g-ratio mapping using neurite orientation dispersion and density imaging (NODDI) and GRE myelin water imaging (GRE-MWI)
Woojin Jung1, Yoonho Nam2, Hui Zhang3, and Jongho Lee1

1Laboratory for Imaging Science and Technology, Department of Electrical and Computer Engineering, Seoul National University, Seoul, Korea, Republic of, 2Department of Radiology, Seoul St. Mary's Hospital, College of Medicine, The Catholic University of Korea, Seoul, Korea, Republic of, 3Department of Computer Science & Centre for Medical Image Computing, University College London, London, United Kingdom

Synopsis

A new in-vivo g-ratio mapping method that combined neurite orientation dispersion and density imaging (NODDI) and GRE myelin water imaging (GRE-MWI) is proposed. The method is substantially fast, taking 17 min for a 2 mm isotropic resolution whole brain g-ratio mapping. The resulting map reveals a reasonable range of g-ratio that has been reported in histology studies.

Introduction

The g-ratio of a nerve fiber, defined as the ratio of axonal diameter (excluding myelin thickness) to fiber diameter (including myelin thickness), is an important biophysical parameter that determines conduction velocity of the fiber [1]. Recently, an MRI method that measures an aggregated g-ratio of a voxel has been proposed [2] using quantitative magnetization transfer (qMT) imaging [3] and neurite orientation dispersion and density imaging (NODDI) [4] to estimate the parameters for the g-ratio formula: myelin volume fraction (MVF) from qMT and axon volume fraction (AVF) from NODDI. However, the qMT imaging requires a substantial scan time (over 30 min) [3] and the resulting fractional pool size (F) requires a scaling coefficient which may need complex histological investigation to translate it to MVF [2]. Recently, GRE myelin water imaging (GRE-MWI) that utilizes multi-echo GRE data to estimate a fast decaying myelin water signal has been proposed as an approach to acquire myelin water fraction (MWF) in a relatively short scan time (<10 min) [5-6]. The resulting MWF is a quantitative measure of myelin water volume and can be converted to MVF using available information (described below). Hence combining GRE-MWI and NODDI can reduce the total scan time and may improve quantification for the g-ratio imaging. In this study, we demonstrate that GRE-MWI and NODDI can deliver a high quality g-ratio map of the whole brain in a reasonable scan time.

Methods

[Data acquisition] Multi-echo GRE data for GRE-MWI and multi-shell diffusion imaging data for NODDI were acquired in four healthy subjects (IRB approved) at 3T. Scan parameters for 3D multi-echo GRE were: TR=67ms, TE=1.52:2.03:31.97ms, flip angle=30°, FOV=256×256×100mm3, resolution=2×2×2mm3, number of echoes=16, and scan time=7min. For NODDI, a three-shell diffusion imaging (b=300s/mm2 with 8 directions; b=700s/mm2 with 32 directions; b=2000s/mm2 with 64 directions; b=0s/mm2 with 13 averages) was acquired using a multi-band (band=2) SE-EPI diffusion sequence. The same resolution, slice thickness and slice number as the GRE sequence were used with TR/TE=4000/95ms, and scan time=9.75min. [Data processing] From the multi-echo GRE data, the signal intensities of myelin water and axonal/interstitial water were estimated by using the GRE-MWI algorithm described in [6]. The resulting signal fraction was translated into the volume fraction with the following considerations. Myelin sheath contains alternating layers of lipid (60% of total volume) and water (40% of total volume) [7-8] and, therefore, the myelin water signal needs to be scaled by 1/0.4 to represent the total myelin volume. The proton density of axon and interstitial space is approximately 0.85 [9] and, therefore, the axonal/interstitial water signal needs to be scaled by 0.85 to represent the axonal/interstitial water volume. As a result, the MVF can be estimated by ([myelin water signal]/0.4)/([axon and interstitial water signal]/0.85+[myelin water signal]/0.4). For AVF estimation, the fraction of free, restricted, and hindered diffusion compartments was calculated using NODDI. Among them, only the restricted diffusion compartment contributes to AVF. Since myelin signal was not detected in NODDI because of long TE of DTI, the restricted diffusion compartment $$$((1-\nu_{iso})\nu_{ic})$$$ did not contain MVF. Hence, the restricted compartment was scaled by (1-MVF) and AVF became $$$((1-MVF)(1-\nu_{iso})\nu_{ic})$$$ [10]. Finally, the aggregated g-ratio was alculated by $$$g=\sqrt{1/(1+MVF/AVF)}$$$ [2]. [Data analysis] ROIs were chosen in optic radiation (OR), corticospinal tract (CST), superior longitudinal fasciculus (SLF), and splenium (SPL) to calculate g-ratios and to test the consistency of our method across the subjects.

Results

The MVF, AVF, and g-ratio maps are shown in Figure 1. The subject averaged ROI values are listed in Table 1. A high g-ratio was observed in CST because of low MVF and high AVF whereas a low g-ratio was found in OR because of high MVF and low AVF. These results are partially supported by histological evidences that CST has large axons [11-12] whereas OR is populated with thick myelin sheath [13]. The mean g-ratios of the ROIs ranged from 0.68 to 0.86. These values are close to an optimum g-ratio of 0.77 for signal conduction [14]. Figure 2 demonstrates consistent g-ratio values across the subjects.

Discussion and conclusion

In this work we demonstrated an approach of generating a g-ratio map using GRE-MWI and NODDI. The scan time was less than 20 min for the whole brain coverage in a 2 mm isotropic voxel size. Further reduction in the DTI scan time is possible with a powerful gradient system [15]. Compared to a recent result using qMT and NODDI, which revealed a relatively uniform g-ratio map [2], our map shows larger g-ratio variations among fiber bundles. This may be in line with the histological measurements showing variations in g-ratio [14,16].

Acknowledgements

This research was supported by the Brain Research Program through the National Research Foundation of Korea(NRF) funded by the Ministry of Science, ICT & Future Planning (NRF-2015M3C7A1031969).

References

[1] RUSHTON WA. A theory of the effects of fibre size in medullated nerve. J Physiol. 1951;115(1):101-22.

[2] Stikov N, Campbell JS, Stroh T, et al. In vivo histology of the myelin g-ratio with magnetic resonance imaging. Neuroimage. 2015;118:397-405.

[3] Sled JG, Pike GB. Quantitative imaging of magnetization transfer exchange and relaxation properties in vivo using MRI. Magn Reson Med. 2001;46(5):923-31.

[4] Zhang H, Schneider T, Wheeler-Kingshott CA, et al. NODDI: Practical in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage. 2012;61(4):1000-16.

[5] Du YP1, Chu R, Hwang D, et al. Fast multislice mapping of the myelin water fraction using multicompartment analysis of T2* decay at 3T: a preliminary postmortem study. Magn Reson Med. 2007;58(5):865-70.

[6] Nam Y, Lee J, Hwang D, et al. Improved estimation of myelin water fraction using complex model fitting . Neuroimage. 2015;116:214-21.

[7] Wharton S, Bowtell R. Fiber orientation-dependent white matter contrast in gradient echo MRI. Proc Natl Acad Sci U S A. 2012;109(45):18559-64.

[8] P Morell, RH Quarles, WT Norton. Formation, structure, and biochemistry of myelin. Basic neurochemistry, 1989, Raven Press, New York.

[9] Sati P, van Gelderen P, Silva AC, et al. Micro-compartment specific T2* relaxation in the brain. Neuroimage. 2013 Aug 15;77:268-78.

[10] Campbell, J.S., Stikov, N., Dougherty, R.F., et al. Combined NODDI and qMT for full-brain g-ratio mapping with complex subvoxel microstructure. Proc. 22nd Scientific Meeting, International Society for Magnetic Resonance in Medicine, Milan, 2014; p. 393.

[11] Schramm U, Graf von Keyserlingk D. Studien u¨ berLatissimusbo¨ gen des Oberarmes. Anat. Anz. 1984. 157:97.

[12] Siawoosh Mohammadi, Daniel Carey, Fred Dick et al. Whole-brain in-vivo measurements of the axonal g-ratio in a group of 19 healthy volunteers. Proc. 23rd Power Pitch Session, International Society for Magnetic Resonance in Medicine, Toronto, 2015.

[13] Kitajima M, Korogi Y, Takahashi M, et al. MR signal intensity of the optic radiation. AJNR Am J Neuroradiol. 1996 (7):1379-83.

[14] Taylor Chomiak, Bin Hu. What is the optimal value of the g-ratio for myelinated fibers in the rat CNS? A theoretical approach. PLoS One. 2009;4(11):e7754.

[15] Sotiropoulos SN, Jbabdi S, Xu J, et al. Advances in diffusion MRI acquisition and processing in the Human Connectome Project. Neuroimage. 2013;80:125-43.

[16] Gibson EM, Purger D, Mount CW et al. Neuronal activity promotes oligodendrogenesis and adaptive myelination in the mammalian brain. Science. 2014;344(6183):1252304.

Figures

Figure 1. Axon Volume Fraction (AVF) from NODDI, Myelin Volume Fraction (MVF) from GRE-MWI and g-ratio map. Arrows in g-ratio map indicate target ROIs: OR (Slice 21), SPL (Slice 26), and CST for left arrow and SLF for right arrow (Slice 32).

Table 1. Average Myelin Volume Fraction (MVF), Axon Volume Fraction (AVF), and g-ratio in the four subjects.

Figure 2. G-ratio of each subject in the four ROIs. Each ROI shows a fairly uniform g-ratios across the subjects.



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