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×100mm
3,
resolution=2×2×2mm
3, number
of echoes=16, and scan time=7min.
For NODDI, a three-shell diffusion
imaging (b=300s/mm
2 with 8 directions; b=700s/mm
2 with 32
directions; b=2000s/mm2 with 64 directions; b=0s/mm
2 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.