Optimizing the intrinsic parallel diffusivity in NODDI: An extensive empirical evaluation
Authors:
Jose M. Guerrero aff001; Nagesh Adluru aff002; Barbara B. Bendlin aff003; H. Hill Goldsmith aff002; Stacey M. Schaefer aff005; Richard J. Davidson aff005; Steven R. Kecskemeti aff002; Hui Zhang aff006; Andrew L. Alexander aff001
Authors place of work:
Department of Medical Physics, University of Wisconsin - Madison, Madison, WI, United States of America
aff001; Waisman Center, University of Wisconsin - Madison, Madison, WI, United States of America
aff002; Department of Medicine, University of Wisconsin - Madison, Madison, WI, United States of America
aff003; Department of Psychology, University of Wisconsin - Madison, Madison, WI, United States of America
aff004; Center for Healthy Minds, University of Wisconsin - Madison, Madison, WI, United States of America
aff005; Department of Computer Science, University College London, London, United Kingdom
aff006
Published in the journal:
PLoS ONE 14(9)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0217118
Summary
Purpose
NODDI is widely used in parameterizing microstructural brain properties. The model includes three signal compartments: intracellular, extracellular, and free water. The neurite compartment intrinsic parallel diffusivity (d∥) is set to 1.7 μm2⋅ms−1, though the effects of this assumption have not been extensively explored. This work investigates the optimality of d∥ = 1.7 μm2⋅ms−1 under varying imaging protocol, age groups, sex, and tissue type in comparison to other biologically plausible values of d∥.
Methods
Model residuals were used as the optimality criterion. The model residuals were evaluated in function of d∥ over the range from 0.5 to 3.0 μm2⋅ms−1. This was done with respect to tissue type (i.e., white matter versus gray matter), sex, age (infancy to late adulthood), and diffusion-weighting protocol (maximum b-value). Variation in the estimated parameters with respect to d∥ was also explored.
Results
Results show d∥ = 1.7 μm2⋅ms−1 is appropriate for adult brain white matter but it is suboptimal for gray matter with optimal values being significantly lower. d∥ = 1.7 μm2⋅ms−1 was also suboptimal in the infant brain for both white and gray matter with optimal values being significantly lower. Minor optimum d∥ differences were observed versus diffusion protocol. No significant sex effects were observed. Additionally, changes in d∥ resulted in significant changes to the estimated NODDI parameters.
Conclusion
The default (d∥) of 1.7 μm2⋅ms−1 is suboptimal in gray matter and infant brains.
Keywords:
Neonates – Age groups – Diffusion tensor imaging – Neuroimaging – Central nervous system – Data acquisition – Neurites – Diffusion magnetic resonance imaging
Introduction
In diffusion weighted magnetic resonance imaging (dMRI), biophysical models are used for relating the dMRI signal to microstructural properties in white and gray matter [1–7]. Neurite orientation dispersion and density imaging (NODDI) [7], separates the brain tissue microstructure landscape into three compartments: intracellular space or neurites (axons, dendrites), extracellular tissue matrix, and a free water compartment. In spite of its shortcomings, much like the case of other techniques such as diffusion tensor imaging (DTI), NODDI offers useful information and has been widely used in the investigation of brain tissue microstructure as a function of early development, cognitive function and aging as well as a number of neurological conditions [8–13].
Biophysical modeling relies on simplifying assumptions about the tissue properties. Besides the separation of tissue into three compartments, the NODDI model is characterized by the following features or assumptions. Each compartment is represented by its own normalized signal and volume fraction. Water exchange between compartments is assumed negligible. Neurites are modeled as sticks (cylinders of zero radius) for capturing highly anisotropic architecture of neuronal tissue. Diffusion inside the neurites is described by a diffusivity parallel to the sticks, which is referred to as the intrinsic diffusivity, d∥, and zero diffusivity perpendicular to them. The orientation distribution function (ODF) of the sticks at each voxel is modeled by an axially symmetric Watson distribution, W [14], which itself is characterized by a concentration parameter κ and mean orientation μ. Highly aligned sticks like those seen in white matter bundles are reflected by high κ values, while highly dispersed sticks like those seen in gray matter fibers are reflected by low κ. The extra-neurite compartment is directionally correlated with the neurite ODF, and modeled as a Gaussian anisotropic compartment.
The local parallel diffusivity of the extracellular space is set equal to the intra-neurite intrinsic diffusivity, d∥, whereas the perpendicular diffusivity d⊥ is related to the neurite water fraction, fic, and d∥ by the mean-field tortuosity model [15] as d⊥ = (1 − fic)d∥. The free-water compartment is modeled as having isotropic diffusion with free diffusivity diso = 3 μm2 ⋅ ms−1 and volume fraction fiso. The intrinsic diffusivity d∥ for NODDI is assumed to be 1.7 μm2⋅ms−1. This is selected to be a biologically reasonable value, which approximates the mean parallel diffusivity from DTI in a healthy coherent white matter region [1]. The parameters that are estimated from acquired data using non-linear gradient descent and heuristic initializations are the water fraction of the neurite compartment fic, the concentration (κ) and mean orientation (μ) of the Watson distribution. The signal S(b, g) from the unit diffusion gradient direction g for sticks oriented along unit vector n and b-matrix (bggt) is given by
Aic, Aec, and Aiso, are the intra-cellular, extra-cellular, and free-water isotropic compartments signal contributions respectively. W(n, μ; κ) is the Watson distribution with κ concentration and oriented along μ. S0 is the un-attenuated signal i.e. S(0, 1), and De(n) = ficd∥nnt + (1 − fic) d∥I3 is the axially symmetric extra-cellular apparent diffusion tensor.
Recently, the model assumptions have become a topic of discussion in the field. The more relevant discussions have focused around the fixed parallel intrinsic diffusivity and equality between parallel intrinsic diffusivity of the extra- and intra-cellular compartment [16]. Of the two, the equality assumption is the more difficult to assess, but has been explored in several reports. While no consensus has been reached, most reports suggest that the intra-cellular parallel intrinsic diffusivity is larger than the extra-cellular one [17–20]. Yet, this may depend on tissue type [21] and most studies have focused on white matter. Also, some sustain that the differences may not be substantial and independent validation experiments are needed [16].
With respect to the fixed diffusivity assumption, [22] proposed a framework for relaxing the fixed constraints. The study reported that microscopic parallel diffusivities varied across the brain, and that white matter values where considerably larger than that assumed by NODDI. It is important to note, however, that the ability to “estimate intrinsic diffusivity” in [22] comes at a cost, which is the reduction to two-compartment model. In this sense, then, the model in [22] is not fully comparable to the model in NODDI, since the former gives up on estimating the CSF volume fraction. Others [23, 24] have also relaxed the fixed diffusivity constraint and made it a free parameter. However, this resulted in unwanted effects on the other parameters in the form of unstable and degenerate estimates. Originally, it was considered unlikely that variation in d∥ across regions and subjects was significant enough to remove trends in the estimated parameters [1]. Additionally, the fixing of d∥ is necessary for stability in the parameter estimates and for speeding up convergence of the fitting procedure. Plus, the value that was chosen was the value that minimized the fitting errors for voxels in the midsagittal plane of the corpus callosum [1].
Taking into consideration the non-consensus on the equality assumption and the still widespread use of the technique, here we choose to build on earlier work [25] which investigated the assumption of fixed diffusivity. This consisted of looking at optimal values of the parallel intrinsic diffusivity according to the model residuals. Results suggested that the default value was reasonable in white matter, but it was sub-optimal in gray matter. While recent publications have found our method useful [26, 27], this earlier work only considered a single axial slice from three age matched participants and dMRI data acquired with the same imaging protocol. For this reason, we propose a more extensive investigation that considers a diverse array of data in terms of age populations, imaging protocols, and is conducted across the full brain.
We choose to do this only for the case of the original NODDI technique and not for its variants [28, 29], as the vast majority of applications have implemented the original version.
Materials and methods
Data
Datasets acquired with multiple b-value sequences (suitable for implementing the NODDI technique [7]) were readily available for use in this work from a number of existing neuroimaging studies. These include imaging data from individuals with a broad range of ages and acquired with imaging protocols that vary in regards to number and magnitude of b values as well as number of diffusion encoding directions. dMRI sets include infants, adolescents, young adults, adults, and aging adults. All dMRI sets were collected on a 3T MR750 Discovery scanner (General Electric, Waukesha, WI). A brief description of each study is provided below and details are summarized in Table 1. All procedures for the included studies were approved by the University of Wisconsin—Madison Institutional Review Board.
Neonates study (Neonates)
Participants are from a study of neonatal white matter microstructure. Diffusion scans contain 6 non-diffusion weighted volumes and diffusion encoded along 63 directions. Other imaging parameters include: TR/TE = 8400/94ms, 2mm isotropic resolution.
Teen study (Teen-I)
Participants in this cohort were drawn from a study of emotion in adolescents. Diffusion scans contain 6 non-diffusion weighted volumes and diffusion encoded along 62 non-collinear directions. Other imaging parameters include: TR/TE = 8400/94 ms and 2 mm isotropic resolution.
Twin teen study (Teen-II)
Participants are from a cohort of 130 adolescent twins. Diffusion scans contain 6 non-diffusion weighted volumes and diffusion encoded along 57 directions. Other parameters include 2.0 mm isotropic resolution and TR/TE = 8000/66.2 ms.
Midlife meditation study (Midlife-I)
Participants in this cohort were drawn from a study of emotion regulation, asthma, and sleep part of the National Center for Complementary and Alternative Medicine (NCCAM). Diffusion scans contain 6 non-diffusion weighted volumes and diffusion encoded along 57 directions. Other parameters include 2.0 mm isotropic resolution and TR/TE = 8000/66.2 ms.
Preclinical Alzheimer’s disease risk study (AD-Risk)
Participants were cognitively unimpaired individuals with and without increased risk for Alzheimer’s disease recruited from the Wisconsin Registry for Alzheimer’s Prevention and Wisconsin Alzheimer’s Disease Research Center. Diffusion scans contain 7 non-diffusion weighted volumes and diffusion encoded along 105 non-collinear directions. Other imaging parameters include: TR/TE = 6500/102 ms, sagittal slices 3mm thick, and in-plane resolution of 2.5 mm × 2.5 mm.
Intrinsic diffusivity optimization
Optimality of d∥ = 1.7 μm2⋅ms−1 was considered by minimizing the model residuals as in [1]. Other biologically plausible values were considered for comparison in the interval [0.5, 3.0] μm2⋅ms−1 in increments of 0.1 μm2⋅ms−1. For each of the 26 values, the model was fitted to the measured dMRI signal voxel by voxel using the Matlab (The MathWorks, Inc., Natick, MA) NODDI toolbox (http://nitrc.org/projects/noddi_toolbox). Predictions of the signal were then calculated at each voxel from the estimated parameters. With the measured and predicted signals for each d∥ setting, the root mean squared (RMS) residual was computed at each voxel. A linear search across the 26 different points was then performed for locating the value of d∥ corresponding to the lowest RMS residual value per voxel. This was done in order to generate a brain map of optimum d∥ and for looking at the optimality of d∥ = 1.7 μm2⋅ms−1 across brain regions.
Tissue type segmentations
White matter (WM) and gray matter (GM) masks were obtained for each individual in order to probe the influence of tissue type on the fitting residuals. This was conducted by running FSL’s [30] FAST tool [31] with meand diffusivity (MD) and fractional anisotropy (FA) maps as input channels. FA and MD maps were obtained from tensor fits using a weighted least squares method. For the AD-risk study, the shells with b values of 4.8 and 7.5 ms⋅μm−2 were excluded in the tensor fitting.
Influences of age, sex, and protocol
The availability of data from the various studies allowed for selection of several subgroups that were organized according to age, sex, and protocol. With the data sets organized this way, the residual analysis was performed for the following three cases:
Groups for age analysis
Subgroups of 16 participants (roughly half male and half female) were selected from three studies as follows: One group of 16 subjects age approximately one month from the Neonates study. One group of 16 subjects ages between 10 and 19 from the Teen-II study. Six groups, 16 subjects each, extracted from the Midlife-I study, for the six age categories of: 20-29, 30-39, 40-49, 50-59, and 60-65 years. Note that, except for the neonates, these data sets have matching protocols so that the main difference per category was age. In order to help disambiguate protocol from age influences, two additional scans were obtained for one adult: one with the infant protocol and one with the adult protocol.
Groups for sex analysis
From the Teen-I study, two subgroups one of 30 females and one of 30 males were selected. The two groups were matched by age (13 years old), so that the main difference between the groups was sex.
Groups for protocol analysis
Two groups of 16 subjects (roughly half females and half males) with ages ranging from 50-59 years were selected, one from the Midlife-I study and one from the AD-risk study. In this case, the assumed main difference between the groups was the acquisition protocol.
Results
The results are organized as follows. (1) We first show how variation in d∥ translates to variability in the estimated parameters. (2) Then, the model RMS residuals, with respect to d∥ are shown to differ between tissue types. (3) This is followed by the presentation of voxel-wise optimized d∥ maps and the ways in which the optimality of d∥ = 1.7 μm2⋅ms−1 is influenced by age, sex, protocol and tissue type.
Estimated model parameters and d∥
Upon completion of the various model fits, the dependence of the estimated model parameters to variations in d∥ was explored. For all model parameter maps, mean values were calculated over WM and GM regions. Fig 1 shows these values plotted with respect to d∥. This analysis reveals a dependence on d∥ for all three parameters irrespective of the study as well as variation in the comparison of parameters among the studies. For example, for gray matter values of d∥ that are lower than the assumed value would weaken variation of the neurite density across the teen and adult subjects. On the other hand, lower values of d∥ in gray matter would enhance differences in the ODF concentration parameter across all studies.
Model residuals with respect to d∥
The values of d∥ that result in the closest agreement between the measured and predicted signals as dictated by the RMS residuals were explored next. For each of the resulting 26 RMS residual maps, mean values across WM and GM were calculated. These are plotted with respect to d∥ in Fig 2. These plots reveal that d∥ values in GM that achieve minimum RMS residuals deviate from the default setting (1.7 μm2⋅ ms−1) for all studies. For WM, the lowest values in the RMS residual curves occur in the neighborhood of the default setting. Notably, most WM curves, with the exception of the Neonate study, exhibit broad ranges of lowest values as compared to the majority of GM curves. The better defined minima in WM for the infants could be related to a maximum b value that better matches the characteristics of the young brain tissue (i.e. longer T2, higher water content) such that diffusion weighting in the signal is more adequate. This is in line with the AD-risk study, which used a max b-value of 7.5 ms⋅μm−2 and the WM RMS residual curves are noticeably more convex. The remaining studies have maximum b values that are likely on the low end of the optimal range for capturing effects of more restrictive intra-neurite environment, which could help explain the shallower curves in WM.
Optimized d∥ maps
Optimum intrinsic diffusivity whole brain maps were created by selecting at each voxel the value that corresponded to the smallest RMS residual. Resulting optimal d∥ maps were median filtered using a box kernel (size 3x3x3 in voxels). The filtering helps to enhance the underlying structure in the distribution of values between white and gray matter. The pattern is spatially consistent before filtering, but it is more difficult to appreciate due to the shallowness of the residual curves for white matter for some of the studies (i.e. Teen-I, Teen-II, Midlife-I). Fig 3 shows optimum d∥ maps for one subject selected randomly from each of the six studies. These maps reveal moderate to substantial contrast between WM and GM regions. The non-uniformly distributed d∥ in these maps suggests that d∥ = 1.7 μm2⋅ms−1 may not be appropriate for all brain regions and all populations.
Optimized d∥ and age
Optimal d∥ maps were computed for the cohort organized by age group. These maps were further masked into WM and GM regions and average optimal d∥ values were obtained for each region. Fig 4A shows the distributions of average optimal d∥ values according to age group. These plots show distinct distributions between WM and GM average optimal d∥ for all age groups greater than 10 years. The majority of WM optimal d∥ values are distributed around the default operating point (1.7 μm2⋅ms−1), while all GM optimal d∥ values are reduced by at least 0.4 μm2⋅ms−1. These trends are fairly consistent for all distributions corresponding to ages 10 years and above. For the group of less than 1 year (i.e. infants) there is a greater degree of closeness between the WM and GM distributions of average optimal d∥ in comparison to the rest of the age groups. In this case, optimal d∥ values fall approximately between 1.4 and 1.5 μm2⋅ms−1 for WM and 1.2 and 1.3 μm2⋅ms−1 for GM. For each age group, a pairwise t-test was conducted in order to assess statistical significance of the tissue-wise difference in average optimal d∥. The testing showed that for all groups the optimum d∥ for GM and WM were significantly different (p < 0.01). A multiple group test revealed that average optimal d∥ is significantly different between the infant and the rest of the older age groups in both WM and GM, while no significant differences were found between any of the other groups. The mean optimum d∥ values for the two additional scans on one adult, Fig 4B, are in agreement with those values from same age group for both the infant and adult protocols, pointing to the fact that the observed trends are more a result of differences in age rather than in protocol.
Optimized d∥ and sex
Optimal d∥ maps were also computed for the cohort organized according to sex. Average optimal d∥ values were obtained across WM and GM regions. The distributions of average optimal d∥ values according to sex category revealed significantly different values between WM and GM with ranges that are consistent with the same age group (10-19 years) from the age-dependence analysis. Yet, no significant effects of sex were observed, a result that is compatible with the age-dependent analysis, which also showed no obvious split in optimal d∥ between the male and female participants.
Optimized d∥ and acquisition protocol
Finally, optimal d∥ maps were also computed for the cohort of subjects with data acquired under differing imaging protocols. Based on the observation that the age dependence analysis revealed no obvious age effects for ages 10 and above, data from the Teen-I study was also included in this cohort despite the unmatched age. This resulted in 3 protocol categories. Fig 5 shows the distribution of WM and GM average optimal d∥ values according to imaging protocol.
The data sets from the groups with the highest b value protocol show optimal d∥ values that are lower than d∥ = 1.7 μm2⋅ms−1. In GM, this analysis reveals a seemingly decaying trend in optimal d∥ distributions with respect to maximum b value. Pair-wise t-tests revealed all distributions in GM are significantly shifted down compared to WM distributions, consistent with the observed trend in the previous age and sex comparisons.
Discussion
In this work we studied the implications of diverse multi-shell dMRI data on the optimality of the NODDI parallel intrinsic diffusivity d∥ = 1.7 μm2⋅ms−1. The results suggests model assumptions for d∥ may be suboptimal for specific ages (i.e., infants) and also in gray matter. Although not examined, the optimality of d∥ = 1.7 μm2⋅ms−1 may also vary with pathology. We also observed that suboptimal d∥ leads to biases in the estimated NODDI parameters. Of particular interest would be a drop of neurite density in gray matter, a result that is consistent with findings in a recent study [32].
For gray matter, the optimal d∥ is significantly lower than 1.7 μm2⋅ms−1. In white matter of the adult brain, values of the optimal d∥ hover around the default d∥ = 1.7 μm2⋅ms−1 and just below the range [1.9, 2.2] μm2⋅ms−1 of intra-axonal diffusivities in white matter reported elsewhere [33], though, further analysis (see below) suggested high FA regions in the adult brain contained average optimal d∥ that falls in this range. It is important to note, however, that the ranges of residual minima in white matter were broad and shallow.
Further, a finer grain analysis indicates that protocol and age also have an impact on the optimality of d∥ = 1.7 μm2⋅ms−1, both in white and gray matter. The age-dependence analysis revealed that the newborn brain optimum d∥ in white and gray matter are closer in value compared to those in the adult brains. Both WM and GM values of optimum d∥ are different, however, from that used in recent studies [24, 34] that have implemented NODDI in the infant brain. The value in these studies was set to 2.0 μm2⋅ms−1, likely because average DTI axial diffusivity in high FA regions (see below) of newborns is close to this number. Interestingly, at this setting, and using the 1.7 μm2⋅ms−1 for the adult brain, nearly any difference between the infants ODF concentration parameter and that of the older age brains would be removed in gray matter. Using the optimal setting for d∥, would result in appreciable differences in ODF concentration parameter between the adults and the infants. On the other hand, using the optimal settings for d∥, would weaken the differences in intra-cellular volume fraction between the infant and the older subjects.
This analysis also showed that in the adult brain optimum intrinsic diffusivity values do not vary appreciably with age. However, optimum d∥ values in GM are much lower than those in WM and different from the default d∥ = 1.7 μm2⋅ms−1. With regards to imaging protocol, high b value and more diffusion weighted volumes appeared to yield less noisy and more stable optimal intrinsic diffusivity and NODDI parameter estimates.
In hindsight, the sub-optimality of the assumed d∥ = 1.7 μm2⋅ms−1 in gray matter is not surprising since this value was originally estimated in the adult corpus callosum [1]. Also, sub-optimality of the current state of the model in gray matter might be related to the idea that the impermeable ‘stick’ representation of neurites is only adequate for myelinated axons but not for dendrites or non-myelinated axons, as others have suggested [35]. In general, however, the variation of optimal intrinsic diffusivity across tissue types is in agreement with findings of axial diffusivity variation across the brain reported in [32].
Studies have reported decreasing DTI axial diffusivity with age [36–38]. Thus, the trend of increasing optimum d∥ with age in WM seen in Fig 4A prompted further investigation. For comparison, averages of DTI axial diffusivity over WM and GM were computed for all subjects in all age groups, Fig 6A and 6B. The resulting axial diffusivity age trajectories are in agreement with previous studies [36–38]. However, while these numbers pertain to the whole of white matter, regional differences in developmental trajectories of DTI quantities in the neonate brain have been observed [39]. In the infants, a further look into high FA (>0.5) regions, which reduce to portions of the corpus callosum and the internal capsule, revealed that average optimal d∥ in these regions is comparable to that seen in the adult global WM. These regions in the infant are thought to be myelinated by one month after birth and to have higher fiber coherence than other white matter areas [39]. The lower FA regions (not shown) in the infant brain, which presumably reflect less or not-yet myelinated axons and or lower fiber coherence, exhibit values of average optimal d∥ that are similar to those of whole WM. For the older age groups, the axial diffusivity distributions in gray matter mimic those of the optimal d∥. For the infants, this is true for both the WM and GM distributions. Also, the optimal d∥ distribution separation between WM and GM is less for the infants than for the rest of the older age groups. Based on all this, it could be speculated that the neonatal gray matter neurites and white matter neurites are more similar than they are in the adults. Therefore, the model fit for less coherent, non-myelinated fibers in neonatal white matter would be more similar to the fit in the neonatal gray matter than to the fit in the adult whole WM, as it is illustrated in Fig 4A.
Limitations
Assumed equal intra- and extra-cellular d∥
As mentioned in the introduction, another important assumption of the model is that of equal d∥ in the intra- and extra-cellular compartments. Thus, one of the limitations of this work is that it was carried out while maintaining this and other assumptions of the model.
In order to glimpse at the appropriateness of this assumption as it pertains to this work, a similar model residual optimization was done for the case where the extra- to intra-cellular parallel diffusivity ratio took on values different than 1. In this case, the model was adjusted so that the extra-cellular diffusivity was expressed as a fraction of the intra-cellular diffusivity value. The ratios ranged from 0.1 to 1.3 in 0.1 increments. In this case the number of fits increases dramatically for each subject (26x13 = 338), as do the memory and time requirements. Therefore, the analysis was restricted to two subjects, one infant and one adult, and for a single axial slice. Additionally, in order to circumvent the long fitting times using the Matlab tool box, for this part of the analysis the AMICO NODDI toolbox [40] was used instead.
Model RMS residuals were calculated for each of the 26 intra-cellular d∥ values in [0.5 μm2⋅ms−1, 3.0 μm2⋅ms−1] and each of the 13 extra- to intra-cellular d∥ ratio values in [0.1,1.3]. Average RMS residuals over WM and GM were plotted with respect to both, the intra-cellular d∥ and the ratio of extra- to intra-cellular d∥. These results are shown by the contour plots in Fig 7. Both in white and in gray matter, the regions of minimum residual values extend over several values in the two dimensions of the graphs. These poorly defined minima point to a multiplicity of solutions when constraints on the model diffusivities are not imposed. Similar results have been presented by other reports [23, 32], which show that unconstrained multi-compartment biophysical models lead to issues in parameter estimation. Particularly, the shape of the lowest residual regions in these contour plots is evocative the pipe-like structures for the fitting cost function landscapes of non-constrained multi-compartment models reported in [32] and [23].
Generalizability
Great effort was made in order to make this as an exhaustive analysis as possible in terms of the diversity of the data that was used. Yet, we acknowledge it is not fully generalizable to the wider scope of neuroimaging biophysical modeling diffusion research, for which it should consider, among others, conditions of pathology and ex-vivo experiments. Nonetheless, we believe that these results are highly informative considering the broad range of ages and imaging protocols investigated. Finally, this analysis was performed for Watson-NODDI only, not for other flavors of the technique which include Bingham-NODDI [28] or NODDIx [29].
Conclusion
In this work, dependence of the estimated NODDI parameters on the parallel intrinsic diffusivity d∥ was observed. Optimum d∥ in white matter of the adult brain is similar to the currently used value d∥ = 1.7 μm2⋅ms−1 but significantly lower in gray matter. Optimal d∥ is also lower than the default value for the newborn brain in white and gray matter. Effects of imaging protocol on the optimum d∥ were also observed. Finally, it is important to consider that, despite its limitations, recent analysis suggests that NODDI metrics provide information that is congruent with histologically equivalent metrics [41].
Zdroje
1. Alexander DC, Hubbard PL, Hall MG, Moore EA, Ptito M, Parker GJ, et al. Orientationally invariant indices of axon diameter and density from diffusion MRI. Neuroimage. 2010;52(4):1374–1389. doi: 10.1016/j.neuroimage.2010.05.043 20580932
2. Assaf Y, Blumenfeld-Katzir T, Yovel Y, Basser PJ. AxCaliber: a method for measuring axon diameter distribution from diffusion MRI. Magnetic resonance in medicine. 2008;59(6):1347–1354. doi: 10.1002/mrm.21577 18506799
3. Assaf Y, Basser PJ. Composite hindered and restricted model of diffusion (CHARMED) MR imaging of the human brain. Neuroimage. 2005;27(1):48–58. doi: 10.1016/j.neuroimage.2005.03.042 15979342
4. Fieremans E, Jensen JH, Helpern JA. White matter characterization with diffusional kurtosis imaging. Neuroimage. 2011;58(1):177–188. doi: 10.1016/j.neuroimage.2011.06.006 21699989
5. Jespersen SN, Kroenke CD, Østergaard L, Ackerman JJ, Yablonskiy DA. Modeling dendrite density from magnetic resonance diffusion measurements. Neuroimage. 2007;34(4):1473–1486. doi: 10.1016/j.neuroimage.2006.10.037 17188901
6. Stanisz GJ, Wright GA, Henkelman RM, Szafer A. An analytical model of restricted diffusion in bovine optic nerve. Magnetic Resonance in Medicine. 1997;37(1):103–111. doi: 10.1002/mrm.1910370115 8978638
7. Zhang H, Schneider T, Wheeler-Kingshott CA, Alexander DC. NODDI: practical in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage. 2012;61(4):1000–1016. doi: 10.1016/j.neuroimage.2012.03.072 22484410
8. Billiet T, Vandenbulcke M, Mädler B, Peeters R, Dhollander T, Zhang H, et al. Age-related microstructural differences quantified using myelin water imaging and advanced diffusion MRI. Neurobiology of aging. 2015;36(6):2107–2121. doi: 10.1016/j.neurobiolaging.2015.02.029 25840837
9. Billiet T, Mädler B, D’Arco F, Peeters R, Deprez S, Plasschaert E, et al. Characterizing the microstructural basis of “unidentified bright objects” in neurofibromatosis type 1: A combined in vivo multicomponent T2 relaxation and multi-shell diffusion MRI analysis. NeuroImage: Clinical. 2014;4:649–658. doi: 10.1016/j.nicl.2014.04.005
10. Caverzasi E, Papinutto N, Castellano A, Zhu AH, Scifo P, Riva M, et al. Neurite orientation dispersion and density imaging color maps to characterize brain diffusion in neurologic disorders. Journal of Neuroimaging. 2016;26(5):494–498. doi: 10.1111/jon.12359 27214558
11. Colgan N, Siow B, O’Callaghan JM, Harrison IF, Wells JA, Holmes HE, et al. Application of neurite orientation dispersion and density imaging (NODDI) to a tau pathology model of Alzheimer’s disease. NeuroImage. 2016;125:739–744. doi: 10.1016/j.neuroimage.2015.10.043 26505297
12. Kamagata K, Hatano T, Okuzumi A, Motoi Y, Abe O, Shimoji K, et al. Neurite orientation dispersion and density imaging in the substantia nigra in idiopathic Parkinson disease. European radiology. 2016;26(8):2567–2577. doi: 10.1007/s00330-015-4066-8 26515546
13. Owen JP, Chang YS, Pojman NJ, Bukshpun P, Wakahiro ML, Marco EJ, et al. Aberrant white matter microstructure in children with 16p11. 2 deletions. Journal of Neuroscience. 2014;34(18):6214–6223. doi: 10.1523/JNEUROSCI.4495-13.2014 24790192
14. Mardia KV, Jupp PE. Directional statistics. vol. 494. John Wiley & Sons; 2009.
15. Szafer A, Zhong J, Gore JC. Theoretical model for water diffusion in tissues. Magnetic resonance in medicine. 1995;33(5):697–712. doi: 10.1002/mrm.1910330516 7596275
16. Jelescu IO, Budde MD. Design and validation of diffusion MRI models of white matter. Frontiers in physics. 2017;5:61. doi: 10.3389/fphy.2017.00061
17. Skinner NP, Kurpad SN, Schmit BD, Tugan Muftuler L, Budde MD. Rapid in vivo detection of rat spinal cord injury with double-diffusion-encoded magnetic resonance spectroscopy. Magnetic resonance in medicine. 2017;77(4):1639–1649. doi: 10.1002/mrm.26243 27080726
18. Szczepankiewicz F, Lasič S, van Westen D, Sundgren PC, Englund E, Westin CF, et al. Quantification of microscopic diffusion anisotropy disentangles effects of orientation dispersion from microstructure: applications in healthy volunteers and in brain tumors. NeuroImage. 2015;104:241–252. doi: 10.1016/j.neuroimage.2014.09.057 25284306
19. Jespersen SN, Olesen JL, Hansen B, Shemesh N. Diffusion time dependence of microstructural parameters in fixed spinal cord. Neuroimage. 2018;182:329–342. doi: 10.1016/j.neuroimage.2017.08.039 28818694
20. Veraart J, Novikov DS, Fieremans E. TE dependent Diffusion Imaging (TEdDI) distinguishes between compartmental T2 relaxation times. Neuroimage. 2018;182:360–369. doi: 10.1016/j.neuroimage.2017.09.030 28935239
21. Novikov DS, Veraart J, Jelescu IO, Fieremans E. Mapping orientational and microstructural metrics of neuronal integrity with in vivo diffusion MRI. arXiv preprint arXiv:160909144. 2016;10(11).
22. Kaden E, Kelm ND, Carson RP, Does MD, Alexander DC. Multi-compartment microscopic diffusion imaging. NeuroImage. 2016;139:346–359. doi: 10.1016/j.neuroimage.2016.06.002 27282476
23. Jelescu IO, Veraart J, Fieremans E, Novikov DS. Degeneracy in model parameter estimation for multi-compartmental diffusion in neuronal tissue. NMR in Biomedicine. 2016;29(1):33–47. doi: 10.1002/nbm.3450 26615981
24. Jelescu IO, Veraart J, Adisetiyo V, Milla SS, Novikov DS, Fieremans E. One diffusion acquisition and different white matter models: how does microstructure change in human early development based on WMTI and NODDI? Neuroimage. 2015;107:242–256. doi: 10.1016/j.neuroimage.2014.12.009 25498427
25. Guerrero J, Adluru N, Kecskemeti S, Davidson R, Alexander A. Investigating the effects of intrinsic diffusivity on neurite orientation dispersion and density imaging (NODDI). International Society for Magnetic Resonance in Medicine (ISMRM) Singapore. 2016;.
26. Genç E, Fraenz C, Schlüter C, Friedrich P, Hossiep R, Voelkle MC, et al. Diffusion markers of dendritic density and arborization in gray matter predict differences in intelligence. Nature communications. 2018;9(1):1905. doi: 10.1038/s41467-018-04268-8 29765024
27. Fukutomi H, Glasser MF, Zhang H, Autio JA, Coalson TS, Okada T, et al. Neurite imaging reveals microstructural variations in human cerebral cortical gray matter. Neuroimage. 2018;182:488–499. doi: 10.1016/j.neuroimage.2018.02.017 29448073
28. Tariq M, Schneider T, Alexander DC, Wheeler-Kingshott CAG, Zhang H. Bingham–NODDI: mapping anisotropic orientation dispersion of neurites using diffusion MRI. NeuroImage. 2016;133:207–223. doi: 10.1016/j.neuroimage.2016.01.046 26826512
29. Farooq H, Xu J, Nam JW, Keefe DF, Yacoub E, Georgiou T, et al. Microstructure imaging of crossing (MIX) white matter fibers from diffusion MRI. Scientific reports. 2016;6:38927. doi: 10.1038/srep38927 27982056
30. Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM. Fsl. Neuroimage. 2012;62(2):782–790. doi: 10.1016/j.neuroimage.2011.09.015 21979382
31. Zhang Y, Brady M, Smith S. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE transactions on medical imaging. 2001;20(1):45–57. doi: 10.1109/42.906424 11293691
32. Novikov DS, Veraart J, Jelescu IO, Fieremans E. Rotationally-invariant mapping of scalar and orientational metrics of neuronal microstructure with diffusion MRI. NeuroImage. 2018;174:518–538. doi: 10.1016/j.neuroimage.2018.03.006 29544816
33. Veraart J, Fieremans E, Novikov DS. Universal power-law scaling of water diffusion in human brain defines what we see with MRI. arXiv preprint arXiv:160909145. 2016;.
34. Kunz N, Zhang H, Vasung L, O’brien KR, Assaf Y, Lazeyras F, et al. Assessing white matter microstructure of the newborn with multi-shell diffusion MRI and biophysical compartment models. Neuroimage. 2014;96:288–299. doi: 10.1016/j.neuroimage.2014.03.057 24680870
35. Veraart J, Fieremans E, Rudrapatna U, Jones D, Novikov DS. Biophysical modeling of the gray matter: does the “stick” model hold? Proceedings of the 27th Annual Meeting of ISMRM, Paris, France 2018;.
36. Mukherjee P, Miller JH, Shimony JS, Philip JV, Nehra D, Snyder AZ, et al. Diffusion-tensor MR imaging of gray and white matter development during normal human brain maturation. American Journal of Neuroradiology. 2002;23(9):1445–1456. 12372731
37. Faria AV, Zhang J, Oishi K, Li X, Jiang H, Akhter K, et al. Atlas-based analysis of neurodevelopment from infancy to adulthood using diffusion tensor imaging and applications for automated abnormality detection. Neuroimage. 2010;52(2):415–428. doi: 10.1016/j.neuroimage.2010.04.238 20420929
38. Yoshida S, Oishi K, Faria AV, Mori S. Diffusion tensor imaging of normal brain development. Pediatric radiology. 2013;43(1):15–27. doi: 10.1007/s00247-012-2496-x 23288475
39. Dean D, Planalp E, Wooten W, Adluru N, Kecskemeti S, Frye C, et al. Mapping white matter microstructure in the one month human brain. Scientific reports. 2017;7(1):9759. doi: 10.1038/s41598-017-09915-6 28852074
40. Daducci A, Canales-Rodríguez EJ, Zhang H, Dyrby TB, Alexander DC, Thiran JP. Accelerated microstructure imaging via convex optimization (AMICO) from diffusion MRI data. NeuroImage. 2015;105:32–44. doi: 10.1016/j.neuroimage.2014.10.026 25462697
41. Grussu F, Schneider T, Tur C, Yates RL, Tachrount M, Ianuş A, et al. Neurite dispersion: a new marker of multiple sclerosis spinal cord pathology? Annals of clinical and translational neurology. 2017;4(9):663–679. doi: 10.1002/acn3.445 28904988
Článek vyšel v časopise
PLOS One
2019 Číslo 9
- Pomohou nám s detekcí karcinomu plic i včely?
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Umělá inteligence v mozkovém implantátu pomáhá muži po mrtvici s dvojjazyčnou konverzací
Nejčtenější v tomto čísle
- Graviola (Annona muricata) attenuates behavioural alterations and testicular oxidative stress induced by streptozotocin in diabetic rats
- CH(II), a cerebroprotein hydrolysate, exhibits potential neuro-protective effect on Alzheimer’s disease
- Comparison between Aptima Assays (Hologic) and the Allplex STI Essential Assay (Seegene) for the diagnosis of Sexually transmitted infections
- Assessment of glucose-6-phosphate dehydrogenase activity using CareStart G6PD rapid diagnostic test and associated genetic variants in Plasmodium vivax malaria endemic setting in Mauritania