One diffusion acquisition and different white matter models: how does microstructure change in human early development based on WMTI and NODDI?

The impact of edema and fiber crossing on diffusion MRI metrics assessed in an ex vivo nerve phantom: Multi-tensor model vs. diffusion orientation distribution function

Diffusion tensor imaging (DTI) has been employed for over 2 decades to noninvasively quantify central nervous system diseases/injuries. However, DTI is an inadequate simplification of diffusion modeling in the presence of coexisting inflammation, edema and crossing nerve fibers. We employed a tissue phantom using fixed mouse trigeminal nerves coated with various amounts of agarose gel to mimic crossing fibers in the presence of vasogenic edema. Diffusivity measures derived by DTI and diffusion basis spectrum imaging (DBSI) were compared at increasing levels of simulated edema and degrees of fiber crossing. Furthermore, we assessed the ability of DBSI, diffusion kurtosis imaging (DKI), generalized q-sampling imaging (GQI), q-ball imaging (QBI) and neurite orientation dispersion and density imaging to resolve fiber crossing, in reference to the gold standard angles measured from structural images. DTI-computed diffusivities and fractional anisotropy were significantly confounded by gel-mimicked edema and crossing fibers. Conversely, DBSI calculated accurate diffusivities of individual fibers regardless of the extent of simulated edema and degrees of fiber crossing angles. Additionally, DBSI accurately and consistently estimated crossing angles in various conditions of gel-mimicked edema when compared with the gold standard (r2 = 0.92, P = 1.9 × 10-9 , bias = 3.9°). Small crossing angles and edema significantly impact the diffusion orientation distribution function, making DKI, GQI and QBI less accurate in detecting and estimating fiber crossing angles. Lastly, we used diffusion tensor ellipsoids to demonstrate that DBSI resolves the confounds of edema and crossing fibers in the peritumoral edema region from a patient with lung cancer metastasis, while DTI failed. In summary, DBSI is able to separate two crossing fibers and accurately recover their diffusivities in a complex environment characterized by increasing crossing angles and amounts of gel-mimicked edema. DBSI also indicated better angular resolution compared with DKI, QBI and GQI.

Nonlinear associations of neurite density and myelin content with age revealed using multicomponent diffusion and relaxometry magnetic resonance imaging

Most magnetic resonance imaging (MRI) studies investigating the relationship between regional brain myelination or axonal density and aging have relied upon nonspecific methods to probe myelin and axonal content, including diffusion tensor imaging and relaxation time mapping. While these studies have provided pivotal insights into changes in cerebral architecture with aging and pathology, details of the underlying microstructural alterations have not been fully elucidated. In the current study, we used the BMC-mcDESPOT analysis, a direct and specific multicomponent relaxometry method for imaging of myelin water fraction (MWF), a marker of myelin content, and NODDI, an emerging multicomponent diffusion technique, for neurite density index (NDI) imaging, a proxy of axonal density. We investigated age-related differences in MWF and NDI in several white matter brain regions in a cohort of cognitively unimpaired participants over a wide age range. Our results indicate a quadratic, inverted U-shape, relationship between MWF and age in all brain regions investigated, suggesting that myelination continues until middle age followed by a decrease at older ages, in agreement with previous work. We found a similarly complex regional association between NDI and age, with several cerebral structures also exhibiting a quadratic, inverted U-shape, relationship. This novel observation suggests an increase in axonal density until the fourth decade of age followed by a rapid loss at older ages. We also observed that these age-related differences in MWF and NDI vary across different brain regions, as expected. Finally, our study indicates no significant association between MWF and NDI in most cerebral structures investigated, although this association approached significance in a limited number of brain regions, indicating the complementary nature of their information and encouraging further investigation. Overall, we find evidence of nonlinear associations between age and myelin or axonal density in a sample of well-characterized adults, using direct myelin and axonal content imaging methods.

Probing Tissue Microarchitecture of the Baby Brain via Spherical Mean Spectrum Imaging

During the first years of life, the human brain undergoes dynamic spatially-heterogeneous changes, invo- lving differentiation of neuronal types, dendritic arbori- zation, axonal ingrowth, outgrowth and retraction, synaptogenesis, and myelination. To better quantify these changes, this article presents a method for probing tissue microarchitecture by characterizing water diffusion in a spectrum of length scales, factoring out the effects of intra-voxel orientation heterogeneity. Our method is based on the spherical means of the diffusion signal, computed over gradient directions for a set of diffusion weightings (i.e., b -values). We decompose the spherical mean profile at each voxel into a spherical mean spectrum (SMS), which essentially encodes the fractions of spin packets undergoing fine- to coarse-scale diffusion proce- sses, characterizing restricted and hindered diffusion stemming respectively from intra- and extra-cellular water compartments. From the SMS, multiple orientation distribution invariant indices can be computed, allowing for example the quantification of neurite density, microscopic fractional anisotropy ( μ FA), per-axon axial/radial diffusivity, and free/restricted isotropic diffusivity. We show that these indices can be computed for the developing brain for greater sensitivity and specificity to development related changes in tissue microstructure. Also, we demonstrate that our method, called spherical mean spectrum imaging (SMSI), is fast, accurate, and can overcome the biases associated with other state-of-the-art microstructure models.

White matter microstructure across the adult lifespan: A mixed longitudinal and cross-sectional study using advanced diffusion models and brain-age prediction

The macro- and microstructural architecture of human brain white matter undergoes substantial alterations throughout development and ageing. Most of our understanding of the spatial and temporal characteristics of these lifespan adaptations come from magnetic resonance imaging (MRI), including diffusion MRI (dMRI), which enables visualisation and quantification of brain white matter with unprecedented sensitivity and detail. However, with some notable exceptions, previous studies have relied on cross-sectional designs, limited age ranges, and diffusion tensor imaging (DTI) based on conventional single-shell dMRI. In this mixed cross-sectional and longitudinal study (mean interval: 15.2 months) including 702 multi-shell dMRI datasets, we combined complementary dMRI models to investigate age trajectories in healthy individuals aged 18 to 94 years (57.12% women). Using linear mixed effect models and machine learning based brain age prediction, we assessed the age-dependence of diffusion metrics, and compared the age prediction accuracy of six different diffusion models, including diffusion tensor (DTI) and kurtosis imaging (DKI), neurite orientation dispersion and density imaging (NODDI), restriction spectrum imaging (RSI), spherical mean technique multi-compartment (SMT-mc), and white matter tract integrity (WMTI). The results showed that the age slopes for conventional DTI metrics (fractional anisotropy [FA], mean diffusivity [MD], axial diffusivity [AD], radial diffusivity [RD]) were largely consistent with previous research, and that the highest performing advanced dMRI models showed comparable age prediction accuracy to conventional DTI. Linear mixed effects models and Wilk's theorem analysis showed that the 'FA fine' metric of the RSI model and 'orientation dispersion' (OD) metric of the NODDI model showed the highest sensitivity to age. The results indicate that advanced diffusion models (DKI, NODDI, RSI, SMT mc, WMTI) provide sensitive measures of age-related microstructural changes of white matter in the brain that complement and extend the contribution of conventional DTI.

Challenges for biophysical modeling of microstructure

The biophysical modeling efforts in diffusion MRI have grown considerably over the past 25 years. In this review, we dwell on the various challenges along the journey of bringing a biophysical model from initial design to clinical implementation, identifying both hurdles that have been already overcome and outstanding issues. First, we describe the critical initial task of selecting which features of tissue microstructure can be estimated using a model and which acquisition protocol needs to be implemented to make the estimation possible. The model performance should necessarily be tested in realistic numerical simulations and in experimental data - adapting the fitting strategy accordingly, and parameter estimates should be validated against complementary techniques, when/if available. Secondly, the model performance and validity should be explored in pathological conditions, and, if appropriate, dedicated models for pathology should be developed. We build on examples from tumors, ischemia and demyelinating diseases. We then discuss the challenges associated with clinical translation and added value. Finally, we single out four major unresolved challenges that are related to: the availability of a microstructural ground truth, the validation of model parameters which cannot be accessed with complementary techniques, the development of a generalized standard model for any brain region and pathology, and the seamless communication between different parties involved in the development and application of biophysical models of diffusion.

Evaluation of the diffusion MRI white matter tract integrity model using myelin histology and Monte-Carlo simulations

Quantitative evaluation of brain myelination has drawn considerable attention. Conventional diffusion-based magnetic resonance imaging models, including diffusion tensor imaging and diffusion kurtosis imaging (DKI),¹ have been used to infer the microstructure and its changes in neurological diseases. White matter tract integrity (WMTI) was proposed as a biophysical model to relate the DKI-derived metrics to the underlying microstructure. Although the model has been validated on ex vivo animal brains, it was not well evaluated with ex vivo human brains. In this study, histological samples (namely corpus callosum) from postmortem human brains have been investigated based on WMTI analyses on a clinical 3T scanner and comparisons with gold standard myelin staining in proteolipid protein and Luxol fast blue. In addition, Monte Carlo simulations were conducted to link changes from ex vivo to in vivo conditions based on the microscale parameters of water diffusivity and permeability. The results show that WMTI metrics, including axonal water fraction AWF, radial extra-axonal diffusivity D_e^⊥, and intra-axonal diffusivity D_awere needed to characterize myelin content alterations. Thus, WMTI model metrics are shown to be promising candidates as sensitive biomarkers of demyelination.

Abnormal neurite orientation dispersion and density imaging of white matter in children with primary nocturnal enuresis

•

NODDI probes white matter (WM) microstructure in PNE children.

•

Different patterns of WM microstructural abnormalities were found in PNE children.

•

NDI of anterior thalamic radiation was correlated with abnormal arousal in PNE children.

Several lines of evidence indicate that multiple abnormalities of gray matter are related to the pathogenesis of primary nocturnal enuresis (PNE); however, few studies have been conducted with respect to abnormalities in white matter (WM) of children with PNE. The present work investigated the microstructure of WM in children with PNE using a neurite orientation dispersion and density imaging (NODDI) method. NODDI data were obtained from 29 children with PNE (age = 9.8 ± 1.2 years, 59% males) and 34 healthy controls (age = 10.3 ± 1.6 years, 56% males) in this study. Multi-b-value diffusion-weighted imaging data were acquired with a 3 T MR system, and the orientation dispersion index (ODI) and neurite density index (NDI) maps were calculated. Tract-Based Spatial Statistics analyses of WM tracts were performed with ODI and NDI maps in children with PNE and controls. Children with PNE had lower ODIs in WM fiber tracts of the bilateral superior longitudinal fasciculus (SLF) and higher ODIs in the bilateral internal capsule (IC) and right anterior thalamic radiation (ATR) than controls. PNE children also had lower NDIs in the bilateral IC and the cingulum and higher NDIs in the bilateral SLF. These changes in NODDI indices, which indicated abnormal neural maturation of the WM microstructures, may be related to abnormal sleep and enuresis in children with PNE.

NODDI in clinical research

Diffusion MRI (dMRI) has proven to be a useful imaging approach for both clinical diagnosis and research investigating the microstructures of nervous tissues, and it has helped us to better understand the neurophysiological mechanisms of many diseases. Though diffusion tensor imaging (DTI) has long been the default tool to analyze dMRI data in clinical research, acquisition with stronger diffusion weightings beyond the DTI regimen is now possible with modern clinical scanners, potentially enabling even more detailed characterization of tissue microstructures. To take advantage of such data, neurite orientation dispersion and density imaging (NODDI) has been proposed as a way to relate the dMRI signal to tissue features via biophysically inspired modeling. The number of reports demonstrating the potential clinical utility of NODDI is rapidly increasing. At the same time, the pitfalls and limitations of NODDI, and general challenges in microstructure modeling, are becoming increasingly recognized by clinicians. dMRI microstructure modeling is a rapidly evolving field with great promise, where people from different scientific backgrounds, such as physics, medicine, biology, neuroscience, and statistics, are collaborating to build novel tools that contribute to improving human healthcare. Here, we review the applications of NODDI in clinical research and discuss future perspectives for investigations toward the implementation of dMRI microstructure imaging in clinical practice.

Abnormal Microstructural Development of the Cerebral Cortex in Neonates With Congenital Heart Disease Is Associated With Impaired Cerebral Oxygen Delivery

Background Abnormal macrostructural development of the cerebral cortex has been associated with hypoxia in infants with congenital heart disease ( CHD ). Animal studies have suggested that hypoxia results in cortical dysmaturation at the cellular level. New magnetic resonance imaging techniques offer the potential to investigate the relationship between cerebral oxygen delivery and cortical microstructural development in newborn infants with CHD . Methods and Results We measured cortical macrostructural and microstructural properties in 48 newborn infants with serious or critical CHD and 48 age-matched healthy controls. Cortical volume and gyrification index were calculated from high-resolution structural magnetic resonance imaging. Neurite density and orientation dispersion indices were modeled using high-angular-resolution diffusion magnetic resonance imaging. Cerebral oxygen delivery was estimated in infants with CHD using phase contrast magnetic resonance imaging and preductal pulse oximetry. We used gray matter-based spatial statistics to examine voxel-wise group differences in cortical microstructure. Microstructural development of the cortex was abnormal in 48 infants with CHD , with regions of increased fractional anisotropy and reduced orientation dispersion index compared with 48 healthy controls, correcting for gestational age at birth and scan (family-wise error corrected for multiple comparisons at P<0.05). Regions of reduced cortical orientation dispersion index in infants with CHD were related to impaired cerebral oxygen delivery ( R²=0.637; n=39). Cortical orientation dispersion index was associated with the gyrification index ( R²=0.589; P<0.001; n=48). Conclusions This study suggests that the primary component of cerebral cortex dysmaturation in CHD is impaired dendritic arborization, which may underlie abnormal macrostructural findings reported in this population, and that the degree of impairment is related to reduced cerebral oxygen delivery.

Diffusion MRI biomarkers of white matter microstructure vary nonmonotonically with increasing cerebral amyloid deposition

Beta amyloid (Aβ) accumulation is the earliest pathological marker of Alzheimer's disease (AD), but early AD pathology also affects white matter (WM) integrity. We performed a cross-sectional study including 44 subjects (23 healthy controls and 21 mild cognitive impairment or early AD patients) who underwent simultaneous PET-MR using 18F-Florbetapir, and were categorized into 3 groups based on Aβ burden: Aβ- [mean mSUVr ≤1.00], Aβi [1.00 < mSUVr <1.17], Aβ+ [mSUVr ≥1.17]. Intergroup comparisons of diffusion MRI metrics revealed significant differences across multiple WM tracts. Aβi group displayed more restricted diffusion (higher fractional anisotropy, radial kurtosis, axonal water fraction, and lower radial diffusivity) than both Aβ- and Aβ+ groups. This nonmonotonic trend was confirmed by significant continuous correlations between mSUVr and diffusion metrics going in opposite direction for 2 cohorts: pooled Aβ-/Aβi and pooled Aβi/Aβ+. The transient period of increased diffusion restriction may be due to inflammation that accompanies rising Aβ burden. In the later stages of Aβ accumulation, neurodegeneration is the predominant factor affecting diffusion.

Microstructural Alterations in Hippocampal Subfields Mediate Age-Related Memory Decline in Humans

Aging, even in the absence of clear pathology of dementia, is associated with cognitive decline. Neuroimaging, especially diffusion-weighted imaging, has been highly valuable in understanding some of these changes in live humans, non-invasively. Traditional tensor techniques have revealed that the integrity of the fornix and other white matter tracts significantly deteriorates with age, and that this deterioration is highly correlated with worsening cognitive performance. However, traditional tensor techniques are still not specific enough to indict explicit microstructural features that may be responsible for age-related cognitive decline and cannot be used to effectively study gray matter properties. Here, we sought to determine whether recent advances in diffusion-weighted imaging, including Neurite Orientation Dispersion and Density Imaging (NODDI) and Constrained Spherical Deconvolution, would provide more sensitive measures of age-related changes in the microstructure of the medial temporal lobe. We evaluated these measures in a group of young (ages 20-38 years old) and older (ages 59-84 years old) adults and assessed their relationships with performance on tests of cognition. We found that the fiber density (FD) of the fornix and the neurite density index (NDI) of the fornix, hippocampal subfields (DG/CA3, CA1, and subiculum), and parahippocampal cortex, varied as a function of age in a cross-sectional cohort. Moreover, in the fornix, DG/CA3, and CA1, these changes correlated with memory performance on the Rey Auditory Verbal Learning Test (RAVLT), even after regressing out the effect of age, suggesting that they were capturing neurobiological properties directly related to performance in this task. These measures provide more details regarding age-related neurobiological properties. For example, a change in fiber density could mean a reduction in axonal packing density or myelination, and the increase in NDI observed might be explained by changes in dendritic complexity or even sprouting. These results provide a far more comprehensive view than previously determined on the possible system-wide processes that may be occurring because of healthy aging and demonstrate that advanced diffusion-weighted imaging is evolving into a powerful tool to study more than just white matter properties.

Sex differences in the effects of gonadal hormones on white matter microstructure development in adolescence

Adolescence is characterized by rapid brain development in white matter (WM) that is attributed in part to surges in gonadal hormones. To date, however, there have been few longitudinal investigations relating changes in gonadal hormones and WM development in adolescents. We acquired diffusion-weighted MRI to estimate mean fractional anisotropy (FA) from 10 WM tracts and salivary testosterone from 51 females and 29 males (ages 9-14 years) who were matched on pubertal stage and followed, on average, for 2 years. We tested whether interactions between sex and changes in testosterone levels significantly explained changes in FA. We found positive associations between changes in testosterone and changes in FA within the corpus callosum, cingulum cingulate, and corticospinal tract in females (all ps<0.05, corrected) and non-significant associations in males. We also collected salivary estradiol from females and found that increases in estradiol were associated with increases in FA in the left uncinate fasciculus (p = 0.04, uncorrected); however, this effect was no longer significant after accounting for changes in testosterone. Our findings indicate there are sex differences in how changes in testosterone relate to changes in WM microstructure of tracts that support impulse control and emotion regulation across the pubertal transition.

Magnitude and timing of major white matter tract maturation from infancy through adolescence with NODDI

White matter maturation is a nonlinear and heterogeneous phenomenon characterized by axonal packing, increased axon caliber, and a prolonged period of myelination. While current in vivo diffusion MRI (dMRI) methods, like diffusion tensor imaging (DTI), have successfully characterized the gross structure of major white matter tracts, these measures lack the specificity required to unravel the distinct processes that contribute to microstructural development. Neurite orientation dispersion and density imaging (NODDI) is a dMRI approach that probes tissue compartments and provides biologically meaningful measures that quantify neurite density index (NDI) and orientation dispersion index (ODI). The purpose of this study was to characterize the magnitude and timing of major white matter tract maturation with NODDI from infancy through adolescence in a cross-sectional cohort of 104 subjects (0.6–18.8 years). To probe the regional nature of white matter development, we use an along-tract approach that partitions tracts to enable more fine-grained analysis. Major white matter tracts showed exponential age-related changes in NDI with distinct maturational patterns. Overall, analyses revealed callosal fibers developed before association fibers. Our along-tract analyses elucidate spatially varying patterns of maturation with NDI that are distinct from those obtained with DTI. ODI was not significantly associated with age in the majority of tracts. Our results support the conclusion that white matter tract maturation is heterochronous process and, furthermore, we demonstrate regional variability in the developmental timing within major white matter tracts. Together, these results help to disentangle the distinct processes that contribute to and more specifically define the time course of white matter maturation.

Neonatal morphometric similarity mapping for predicting brain age and characterizing neuroanatomic variation associated with preterm birth

Multi-contrast MRI captures information about brain macro- and micro-structure which can be combined in an integrated model to obtain a detailed "fingerprint" of the anatomical properties of an individual's brain. Inter-regional similarities between features derived from structural and diffusion MRI, including regional volumes, diffusion tensor metrics, neurite orientation dispersion and density imaging measures, can be modelled as morphometric similarity networks (MSNs). Here, individual MSNs were derived from 105 neonates (59 preterm and 46 term) who were scanned between 38 and 45 weeks postmenstrual age (PMA). Inter-regional similarities were used as predictors in a regression model of age at the time of scanning and in a classification model to discriminate between preterm and term infant brains. When tested on unseen data, the regression model predicted PMA at scan with a mean absolute error of 0.70  ±  0.56 weeks, and the classification model achieved 92% accuracy. We conclude that MSNs predict chronological brain age accurately; and they provide a data-driven approach to identify networks that characterise typical maturation and those that contribute most to neuroanatomic variation associated with preterm birth.

Diffusion Kurtosis Imaging maps neural damage in the EAE model of multiple sclerosis

Diffusion kurtosis imaging (DKI) is an imaging modality that yields novel disease biomarkers and in combination with nervous tissue modeling, provides access to microstructural parameters. Recently, DKI and subsequent estimation of microstructural model parameters has been used for assessment of tissue changes in neurodegenerative diseases and associated animal models. In this study, mouse spinal cords from the experimental autoimmune encephalomyelitis (EAE) model of multiple sclerosis (MS) were investigated for the first time using DKI in combination with biophysical modeling to study the relationship between microstructural metrics and degree of animal dysfunction. Thirteen spinal cords were extracted from animals with varied grades of disability and scanned in a high-field MRI scanner along with five control specimen. Diffusion weighted data were acquired together with high resolution T₂* images. Diffusion data were fit to estimate diffusion and kurtosis tensors and white matter modeling parameters, which were all used for subsequent statistical analysis using a linear mixed effects model. T₂* images were used to delineate focal demyelination/inflammation. Our results reveal a strong relationship between disability and measured microstructural parameters in normal appearing white matter and gray matter. Relationships between disability and mean of the kurtosis tensor, radial kurtosis, radial diffusivity were similar to what has been found in other hypomyelinating MS models, and in patients. However, the changes in biophysical modeling parameters and in particular in extra-axonal axial diffusivity were clearly different from previous studies employing other animal models of MS. In conclusion, our data suggest that DKI and microstructural modeling can provide a unique contrast capable of detecting EAE-specific changes correlating with clinical disability.

Imaging G-Ratio in Multiple Sclerosis Using High-Gradient Diffusion MRI and Macromolecular Tissue Volume

BACKGROUND AND PURPOSE

Remyelination represents an area of great therapeutic interest in multiple sclerosis but currently lacks a robust imaging marker. The purpose of this study was to use high-gradient diffusion MRI and macromolecular tissue volume imaging to obtain estimates of axonal volume fraction, myelin volume fraction, and the imaging g-ratio in patients with MS and healthy controls and to explore their relationship to neurologic disability in MS.

MATERIALS AND METHODS

Thirty individuals with MS (23 relapsing-remitting MS, 7 progressive MS) and 19 age-matched healthy controls were scanned on a 3T MRI scanner equipped with 300 mT/m maximum gradient strength using a comprehensive multishell diffusion MRI protocol. Macromolecular tissue volume imaging was performed to quantify the myelin volume fraction. Diffusion data were fitted to a 3-compartment model of white matter using a spheric mean approach to yield estimates of axonal volume fraction. The imaging g-ratio was calculated from the ratio of myelin volume fraction and axonal volume fraction. Imaging metrics were compared between groups using 2-sided t tests with a Bonferroni correction.

RESULTS

The mean g-ratio was significantly elevated in lesions compared with normal-appearing WM (0.74 vs 0.67, P < .001). Axonal volume fraction (0.17 vs 0.23, P < .001) and myelin volume fraction (0.17 vs 0.25, P < .001) were significantly lower in lesions than normal-appearing WM. Myelin volume fraction was lower in normal-appearing WM compared with that in healthy controls (0.25 vs 0.27, P = .009). Disability, as measured by the Expanded Disability Status Scale, was significantly associated with myelin volume fraction (β = -40.5, P = .001) and axonal volume fraction (β = -41.0, P = .016) in normal-appearing WM.

CONCLUSIONS

The imaging g-ratio may serve as a biomarker for the relative degree of axonal and myelin loss in MS.

Population-based Bayesian regularization for microstructural diffusion MRI with NODDIDA

PURPOSE

Information on the brain microstructure can be probed by Diffusion Magnetic Resonance Imaging (dMRI). Neurite Orientation Dispersion and Density Imaging with Diffusivities Assessment (NODDIDA) is one of the simplest microstructural model proposed. However, the estimation of the NODDIDA parameters from clinically plausible dMRI acquisition is ill-posed, and different parameter sets can describe the same measurements equally well. A few approaches to resolve this problem focused on developing better optimization strategies for this non-convex optimization. However, this fundamentally does not resolve ill-posedness. This article introduces a Bayesian estimation framework, which is regularized through knowledge from an extensive dMRI measurement set on a population of healthy adults (henceforth population-based prior).

METHODS

We reformulate the problem as a Bayesian maximum a posteriori estimation, which includes as a special case previous approach using non-informative uniform priors. A population-based prior is estimated from 35 subjects of the MGH Adult Diffusion data (Human Connectome Project), acquired with an extensive acquisition protocol including high b-values. The accuracy and robustness of different approaches with and without the population-based prior is tested on subsets of the MGH dataset, and an independent dataset from a clinically comparable scanner, with only clinically plausible dMRI measurements.

RESULTS

The population-based prior produced substantially more accurate and robust parameter estimates, compared to the conventional uniform priors, for clinically feasible protocols, without introducing any evident bias.

CONCLUSIONS

The use of the proposed Bayesian population-based prior can lead to clinically feasible and robust estimation of NODDIDA parameters without changing the acquisition protocol.

Optimizing the intrinsic parallel diffusivity in NODDI: An extensive empirical evaluation

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 μm²⋅ms⁻¹, though the effects of this assumption have not been extensively explored. This work investigates the optimality of d_∥ = 1.7 μm²⋅ms⁻¹ 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 μm²⋅ms⁻¹. 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 μm²⋅ms⁻¹ is appropriate for adult brain white matter but it is suboptimal for gray matter with optimal values being significantly lower. d_∥ = 1.7 μm²⋅ms⁻¹ 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 μm²⋅ms⁻¹ is suboptimal in gray matter and infant brains.

White Matter Tract Integrity: An Indicator of Axonal Pathology after Mild Traumatic Brain Injury

We seek to elucidate the underlying pathophysiology of injury sustained after mild traumatic brain injury (mTBI) using multi-shell diffusion magnetic resonance imaging, deriving compartment-specific white matter tract integrity (WMTI) metrics. WMTI allows a more biophysical interpretation of white matter (WM) changes by describing microstructural characteristics in both intra- and extra-axonal environments. Thirty-two patients with mTBI within 30 days of injury and 21 age- and sex-matched controls were imaged on a 3 Tesla magnetic resonance scanner. Multi-shell diffusion acquisition was performed with five b-values (250-2500 sec/mm²) along 6-60 diffusion encoding directions. Tract-based spatial statistics (TBSS) was used with family-wise error (FWE) correction for multiple comparisons. TBSS results demonstrated focally lower intra-axonal diffusivity (D_axon) in mTBI patients in the splenium of the corpus callosum (sCC; p < 0.05, FWE-corrected). The area under the curve value for D_axon was 0.76 with a low sensitivity of 46.9% but 100% specificity. These results indicate that D_axon may be a useful imaging biomarker highly specific for mTBI-related WM injury. The observed decrease in D_axon suggests restriction of the diffusion along the axons occurring shortly after injury.

Diffusion MRI microstructural models in the cervical spinal cord - Application, normative values, and correlations with histological analysis

Multi-compartment tissue modeling using diffusion magnetic resonance imaging has proven valuable in the brain, offering novel indices sensitive to the tissue microstructural environment in vivo on clinical MRI scanners. However, application, characterization, and validation of these models in the spinal cord remain relatively under-studied. In this study, we apply a diffusion "signal" model (diffusion tensor imaging, DTI) and two commonly implemented "microstructural" models (neurite orientation dispersion and density imaging, NODDI; spherical mean technique, SMT) in the human cervical spinal cord of twenty-one healthy controls. We first provide normative values of DTI, SMT, and NODDI indices in a number of white matter ascending and descending pathways, as well as various gray matter regions. We then aim to validate the sensitivity and specificity of these diffusion-derived contrasts by relating these measures to indices of the tissue microenvironment provided by a histological template. We find that DTI indices are sensitive to a number of microstructural features, but lack specificity. The microstructural models also show sensitivity to a number of microstructure features; however, they do not capture the specific microstructural features explicitly modelled. Although often regarded as a simple extension of the brain in the central nervous system, it may be necessary to re-envision, or specifically adapt, diffusion microstructural models for application to the human spinal cord with clinically feasible acquisitions - specifically, adjusting, adapting, and re-validating the modeling as it relates to both theory (i.e. relevant biology, assumptions, and signal regimes) and parameter estimation (for example challenges of acquisition, artifacts, and processing).

Resolving degeneracy in diffusion MRI biophysical model parameter estimation using double diffusion encoding

PURPOSE

Biophysical tissue models are increasingly used in the interpretation of diffusion MRI (dMRI) data, with the potential to provide specific biomarkers of brain microstructural changes. However, it has been shown recently that, in the general Standard Model, parameter estimation from dMRI data is ill-conditioned even when very high b-values are applied. We analyze this issue for the Neurite Orientation Dispersion and Density Imaging with Diffusivity Assessment (NODDIDA) model and demonstrate that its extension from single diffusion encoding (SDE) to double diffusion encoding (DDE) resolves the ill-posedness for intermediate diffusion weightings, producing an increase in accuracy and precision of the parameter estimation.

METHODS

We analyze theoretically the cumulant expansion up to fourth order in b of SDE and DDE signals. Additionally, we perform in silico experiments to compare SDE and DDE capabilities under similar noise conditions.

RESULTS

We prove analytically that DDE provides invariant information non-accessible from SDE, which makes the NODDIDA parameter estimation injective. The in silico experiments show that DDE reduces the bias and mean square error of the estimation along the whole feasible region of 5D model parameter space.

CONCLUSIONS

DDE adds additional information for estimating the model parameters, unexplored by SDE. We show, as an example, that this is sufficient to solve the previously reported degeneracies in the NODDIDA model parameter estimation.

White matter plasticity and maturation in human cognition

White matter plasticity likely plays a critical role in supporting cognitive development. However, few studies have used the imaging methods specific to white matter tissue structure or experimental designs sensitive to change in white matter necessary to elucidate these relations. Here we briefly review novel imaging approaches that provide more specific information regarding white matter microstructure. Furthermore, we highlight recent studies that provide greater clarity regarding the relations between changes in white matter and cognition maturation in both healthy children and adolescents and those with white matter insult. Finally, we examine the hypothesis that white matter is linked to cognitive function via its impact on neural synchronization. We test this hypothesis in a population of children and adolescents with recurrent demyelinating syndromes. Specifically, we evaluate group differences in white matter microstructure within the optic radiation; and neural phase synchrony in visual cortex during a visual task between 25 patients and 28 typically developing age-matched controls. Children and adolescents with demyelinating syndromes show evidence of myelin and axonal compromise and this compromise predicts reduced phase synchrony during a visual task compared to typically developing controls. We investigate one plausible mechanism at play in this relationship using a computational model of gamma generation in early visual cortical areas. Overall, our findings show a fundamental connection between white matter microstructure and neural synchronization that may be critical for cognitive processing. In the future, longitudinal or interventional studies can build upon our knowledge of these exciting relations between white matter, neural communication, and cognition.

Neurite orientation dispersion and density imaging of mouse brain microstructure

Advanced biophysical models like neurite orientation dispersion and density imaging (NODDI) have been developed to estimate the microstructural complexity of voxels enriched in dendrites and axons for both in vivo and ex vivo studies. NODDI metrics derived from high spatial and angular resolution diffusion MRI using the fixed mouse brain as a reference template have not yet been reported due in part to the extremely long scan time required. In this study, we modified the three-dimensional diffusion-weighted spin-echo pulse sequence for multi-shell and undersampling acquisition to reduce the scan time. This allowed us to acquire several exhaustive datasets that would otherwise not be attainable. NODDI metrics were derived from a complex 8-shell diffusion (1000-8000 s/mm²) dataset with 384 diffusion gradient-encoding directions at 50 µm isotropic resolution. These provided a foundation for exploration of tradeoffs among acquisition parameters. A three-shell acquisition strategy covering low, medium, and high b values with at least angular resolution of 64 is essential for ex vivo NODDI experiments. The good agreement between neurite density index (NDI) and the orientation dispersion index (ODI) with the subsequent histochemical analysis of myelin and neuronal density highlights that NODDI could provide new insight into the microstructure of the brain. Furthermore, we found that NDI is sensitive to microstructural variations in the corpus callosum using a well-established demyelination cuprizone model. The study lays the ground work for developing protocols for routine use of high-resolution NODDI method in characterizing brain microstructure in mouse models.

A deep network for tissue microstructure estimation using modified LSTM units

Diffusion magnetic resonance imaging (dMRI) offers a unique tool for noninvasively assessing tissue microstructure. However, accurate estimation of tissue microstructure described by complicated signal models can be challenging when a reduced number of diffusion gradients are used. Deep learning based microstructure estimation has recently been developed and achieved promising results. In particular, optimization-based learning, where deep network structures are constructed by unfolding the iterative processes performed for solving optimization problems, has demonstrated great potential in accurate microstructure estimation with a reduced number of diffusion gradients. In this work, using the optimization-based learning strategy, we propose a deep network structure that is motivated by the use of historical information in iterative optimization for tissue microstructure estimation, and such incorporation of historical information has not been previously explored in the design of deep networks for microstructure estimation. We assume that (1) diffusion signals can be sparsely represented by a dictionary and its coefficients jointly in the spatial and angular domain, and (2) tissue microstructure can be computed from the sparse representation. Following these assumptions, our network comprises two cascaded stages. The first stage takes image patches as input and computes the spatial-angular sparse representation of the input with learned weights. Specifically, the network structure in the first stage is constructed by unfolding an iterative process for solving sparse reconstruction problems, where historical information is incorporated. The components in this network can be shown to correspond to modified long short-term memory (LSTM) units. In the second stage, fully connected layers are added to compute the mapping from the sparse representation to tissue microstructure. The weights in the two stages are learned jointly by minimizing the mean squared error of microstructure estimation. Experiments were performed on dMRI scans with a reduced number of diffusion gradients. For demonstration, we evaluated the estimation of tissue microstructure described by three signal models: the neurite orientation dispersion and density imaging (NODDI) model, the spherical mean technique (SMT) model, and the ensemble average propagator (EAP) model. The results indicate that the proposed approach outperforms competing methods.

Imaging brain microstructure with diffusion MRI: practicality and applications

This article gives an overview of microstructure imaging of the brain with diffusion MRI and reviews the state of the art. The microstructure-imaging paradigm aims to estimate and map microscopic properties of tissue using a model that links these properties to the voxel scale MR signal. Imaging techniques of this type are just starting to make the transition from the technical research domain to wide application in biomedical studies. We focus here on the practicalities of both implementing such techniques and using them in applications. Specifically, the article summarizes the relevant aspects of brain microanatomy and the range of diffusion-weighted MR measurements that provide sensitivity to them. It then reviews the evolution of mathematical and computational models that relate the diffusion MR signal to brain tissue microstructure, as well as the expanding areas of application. Next we focus on practicalities of designing a working microstructure imaging technique: model selection, experiment design, parameter estimation, validation, and the pipeline of development of this class of technique. The article concludes with some future perspectives on opportunities in this topic and expectations on how the field will evolve in the short-to-medium term.

Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation

We review, systematize and discuss models of diffusion in neuronal tissue, by putting them into an overarching physical context of coarse-graining over an increasing diffusion length scale. From this perspective, we view research on quantifying brain microstructure as occurring along three major avenues. The first avenue focusses on transient, or time-dependent, effects in diffusion. These effects signify the gradual coarse-graining of tissue structure, which occurs qualitatively differently in different brain tissue compartments. We show that transient effects contain information about the relevant length scales for neuronal tissue, such as the packing correlation length for neuronal fibers, as well as the degree of structural disorder along the neurites. The second avenue corresponds to the long-time limit, when the observed signal can be approximated as a sum of multiple nonexchanging anisotropic Gaussian components. Here, the challenge lies in parameter estimation and in resolving its hidden degeneracies. The third avenue employs multiple diffusion encoding techniques, able to access information not contained in the conventional diffusion propagator. We conclude with our outlook on future directions that could open exciting possibilities for designing quantitative markers of tissue physiology and pathology, based on methods of studying mesoscopic transport in disordered systems.

A review of diffusion MRI of typical white matter development from early childhood to young adulthood

Understanding typical, healthy brain development provides a baseline from which to detect and characterize brain anomalies associated with various neurological or psychiatric disorders and diseases. Diffusion MRI is well suited to study white matter development, as it can virtually extract individual tracts and yield parameters that may reflect alterations in the underlying neural micro-structure (e.g. myelination, axon density, fiber coherence), though it is limited by its lack of specificity and other methodological concerns. This review summarizes the last decade of diffusion imaging studies of healthy white matter development spanning childhood to early adulthood (4-35 years). Conclusions about anatomical location, rates, and timing of white matter development with age are discussed, as well as the influence of image acquisition, analysis, age range/sample size, and statistical model. Despite methodological variability between studies, some consistent findings have emerged from the literature. Specifically, diffusion studies of neurodevelopment overwhelmingly demonstrate regionally varying increases of fractional anisotropy and decreases of mean diffusivity during childhood and adolescence, some of which continue into adulthood. While most studies use linear fits to model age-related changes, studies with sufficient sample sizes and age range provide clear evidence that white matter development (as indicated by diffusion) is non-linear. Several studies further suggest that maturation in association tracts with frontal-temporal connections continues later than commissural and projection tracts. The emerging contributions of more advanced diffusion methods are also discussed, as they may reveal new aspects of white matter development. Although non-specific, diffusion changes may reflect increases of myelination, axonal packing, and/or coherence with age that may be associated with changes in cognition.

Comparison of NODDI and spherical mean signal for measuring intra-neurite volume fraction

PURPOSE

Neurite orientation dispersion and density imaging (NODDI) is a clinically feasible approach to measure intra-neurite volume fraction (fin). However, the sophisticated fitting procedure takes several hours. And the NODDI model relied on several questionable assumptions. Recent analytical work demonstrated that fin could be simply calculated from the spherical mean signal (MEANS) averaged over all gradient directions with a more solid theoretical foundation. The current study aims to compare NODDI and MEANS for measuring fin in human brain and investigate the potential of MEANS as a fast approach in clinics.

METHODS

NODDI fin and MEANS fin were measured and compared on the same dataset. NODDI fin was obtained using the NODDI MATLAB Toolbox. MEANS fin is the product of the spherical mean signal and 2bD/π, where D is the intra-neurite intrinsic diffusivity.

RESULTS

NODDI fin and MEANS fin maps are similar. The voxel-by-voxel correlation suggests that NODDI fin and MEANS fin are approximately equivalent to each other.

CONCLUSION

MEANS may have potential to serve a fast and simple approach to estimate fin in clinics.

Application of neurite orientation dispersion and density imaging to characterize brain microstructural abnormalities in type-2 diabetics with mild cognitive impairment

BACKGROUND

Diffusion-tensor-imaging (DTI) is sensitive in detecting white matter changes in type-2 diabetes mellitus (T2DM). However, DTI indices can be affected by either neurite density or spatial variation. A novel diffusion MRI technique, termed neurite orientation dispersion and density imaging (NODDI), can provide distinct indices of fiber density and dispersion.

PURPOSE

To characterize brain microstructural alterations in T2DM patients with mild cognitive impairment (MCI) using the NODDI model.

STUDY TYPE

Cross-sectional.

SUBJECTS

Twenty T2DM patients with (DM-MCI group), 18 age- and gender-matched T2DM patients with normal cognition (DM-NC group), and 28 euglycemic healthy controls (HC).

FIELD STRENGTH/SEQUENCE

3T/NODDI.

ASSESSMENT

Diffusion data were analyzed using tract-based-spatial-statistics (TBSS) analysis in white matter and voxel-based analysis in both white and gray matter. NODDI indices, including intracellular volume fraction (Vic) and orientation dispersion index (ODI), were estimated from multiple regions and compared among these groups.

STATISTICAL TESTS

Differences between groups were compared by Student's t-test, Pearson chi-square test, or analysis of variance when appropriate. Correlation analyses were performed to investigate the relationship between NODDI variables and clinical measurements.

RESULTS

Whole-brain TBSS revealed that 2.29% and 2.02% of the white matter regions exhibited decreased fractional anisotropy and Vic, respectively, between the DM-NC and HC, while considerably larger white matter areas showed decreased fractional anisotropy (38.38%) and Vic (34.64%) between the DM-MCI and HC (Student's t-test, P < 0.05). However, the angular variation of neurites, characterized by ODI, exhibited very little (0.1%, P < 0.05) or no difference (P > 0.05) between either the DM-MCI or DM-NC groups and HC. Decreased Vic values in the genu of the corpus callosum (R = 0.580, 0.551 and 0.586, P < 0.01) and thalamus (R = 0.570, 0.616, and 0.595, P < 0.05) correlated with glycosylated hemoglobin A1c level, disease duration, and neuropsychological scores, respectively.

DATA CONCLUSION

T2DM patients with cognitive decline had reduced Vic, which indicated decreased density of axons and dendrites. NODDI might be able to help probe microstructural changes in white and gray matter and provide information on diabetic encephalopathy, including those with cognitive impairment.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:889-898.

IMAge/enGINE: a freely available software for rapid computation of high-dimensional quantification

BACKGROUND

High-dimensional image data including diffusion weighted imaging, diffusion tensor imaging and dynamic imaging are important in exploring the connectivity, cellularity, pharmacokinetic and blood supply. IMAge/enGINE is software especially designed for high-dimensional medical image computing.

METHODS

IMAge/enGINE is implemented based on open-source and cross-platform tools such as Qt, ITK and VTK. It processes the high-dimensional image data in a slice-by-slice computation mechanism. For computational efficiency, C++ is used for implementing IMAge/enGINE and multi-thread computing is handled in the scale of voxels. The architecture of IMAge/enGINE is modularized for easier extension.

RESULTS

IMAge/enGINE has following features: (I) IMAge/enGINE is free for research use; (II) it has an easy-to-use graphic user interface designed for clinical users without programming or engineering background; (III) its frame work is open-source and extensible. Developers can implement algorithms as modules and integrate them into IMAge/enGINE or generate their own application.

CONCLUSIONS

The source of IMAge/enGINE is hosted at https://github.com/VusionMed/IMAge-enGINE. Multiple diffusion and perfusion models are implemented and integrated into IMAge/enGINE and its binaries can be downloaded freely at http://www.vusion.com.cn/?page_id=14971.

Applying microstructural models to understand the role of white matter in cognitive development

Diffusion MRI (dMRI) holds great promise for illuminating the biological changes that underpin cognitive development. The diffusion of water molecules probes the cellular structure of brain tissue, and biophysical modeling of the diffusion signal can be used to make inferences about specific tissue properties that vary over development or predict cognitive performance. However, applying these models to study development requires that the parameters can be reliably estimated given the constraints of data collection with children. Here we collect repeated scans using a typical multi-shell diffusion MRI protocol in a group of children (ages 7-12) and use two popular modeling techniques to examine individual differences in white matter structure. We first assess scan-rescan reliability of model parameters and show that axon water faction can be reliably estimated from a relatively fast acquisition, without applying spatial smoothing or de-noising. We then investigate developmental changes in the white matter, and individual differences that correlate with reading skill. Specifically, we test the hypothesis that previously reported correlations between reading skill and diffusion anisotropy in the corpus callosum reflect increased axon water fraction in poor readers. Both models support this interpretation, highlighting the utility of these approaches for testing specific hypotheses about cognitive development.

Microscopic anisotropy misestimation in spherical-mean single diffusion encoding MRI

Purpose

Microscopic fractional anisotropy (µFA) can disentangle microstructural information from orientation dispersion. While double diffusion encoding (DDE) MRI methods are widely used to extract accurate µFA, it has only recently been proposed that powder‐averaged single diffusion encoding (SDE) signals, when coupled with the diffusion standard model (SM) and a set of constraints, could be used for µFA estimation. This study aims to evaluate µFA as derived from the spherical mean technique (SMT) set of constraints, as well as more generally for powder‐averaged SM signals.

Methods

SDE experiments were performed at 16.4 T on an ex vivo mouse brain (Δ/δ = 12/1.5 ms). The µFA maps obtained from powder‐averaged SDE signals were then compared to maps obtained from DDE‐MRI experiments (Δ/τ/δ = 12/12/1.5 ms), which allow a model‐free estimation of µFA. Theory and simulations that consider different types of heterogeneity are presented for corroborating the experimental findings.

Results

µFA, as well as other estimates derived from powder‐averaged SDE signals produced large deviations from the ground truth in both gray and white matter. Simulations revealed that these misestimations are likely a consequence of factors not considered by the underlying microstructural models (such as intercomponent and intracompartmental kurtosis).

Conclusion

Powder‐averaged SMT and (2‐component) SM are unable to accurately report µFA and other microstructural parameters in ex vivo tissues. Improper model assumptions and constraints can significantly compromise parameter specificity. Further developments and validations are required prior to implementation of these models in clinical or preclinical research.

Early assessment of lateralization and sex influences on the microstructure of the white matter corticospinal tract in healthy term neonates

We assessed the sex and the lateralization differences in the corticospinal tract (CST) during the early postnatal period. Twenty-five healthy term neonates (13 girls, aged 39.2 ± 1.2 weeks, and 12 boys aged 38.6 ± 3.0 weeks) underwent Diffusion Tensor Imaging (DTI). Fiber tracking was performed to extract bilaterally the CST pathways and to quantify the parallel (E₁ ) and perpendicular (E₂₃ ) diffusions, the apparent diffusion coefficient (ADC), and fractional anisotropy (FA). The measurements were performed on the entire CST fibers and on four segments: base of the pons (CST-Po), cerebral peduncles (CST-CP), posterior limb of the internal capsule (CST-PLIC), and corona-radiata (CST-CR). Significantly higher E₁ , lower E_23, and higher FA in the right compared to the left were noted in the CST-PLIC of the girls. Significantly lower E₂₃ and lower ADC with higher FA in the right compared to left were observed in the CST-CP of the boys. Moreover, the CST-PLIC of the boys had significantly higher E₁ in the right compared to the left. There was a significant increase in left CST E₁ of boys when compared with girls. Girls had a significantly lower E₁ , lower E₂₃ and, lower ADC in the left CST-CP compared with boys. In addition, girls had a significantly lower E₂₃ and higher FA in the right CST-PLIC compared with boys. Sex differences and lateralization in structure-based segments of the CST were found in healthy term infants during early postnatal period. These findings are vital to understanding motor development of healthy term born neonates to better interpret newborn infants with abnormal neurodevelopment.

Diffusional kurtosis imaging of the corpus callosum in autism

Background

The corpus callosum is implicated in the pathophysiology of autism spectrum disorder (ASD). However, specific structural deficits and underlying mechanisms are yet to be well defined.

Methods

We employed diffusional kurtosis imaging (DKI) metrics to characterize white matter properties within five discrete segments of the corpus callosum in 17 typically developing (TD) adults and 16 age-matched participants with ASD without co-occurring intellectual disability (ID). The DKI metrics included axonal water fraction (f_axon) and intra-axonal diffusivity (D_axon), which reflect axonal density and caliber, and extra-axonal radial (RD_extra) and axial (AD_extra) diffusivities, which reflect myelination and microstructural organization of the extracellular space. The relationships between DKI metrics and processing speed, a cognitive feature known to be impaired in ASD, were also examined.

Results

ASD group had significantly decreased callosal f_axon and D_axon (p = .01 and p = .045), particularly in the midbody, isthmus, and splenium. Regression analysis showed that variation in DKI metrics, primarily in the mid and posterior callosal regions explained up to 70.7% of the variance in processing speed scores for TD (p = .001) but not for ASD (p > .05).

Conclusion

Decreased DKI metrics suggested that ASD may be associated with axonal deficits such as reduced axonal caliber and density in the corpus callosum, especially in the mid and posterior callosal areas. These data suggest that impaired interhemispheric connectivity may contribute to decreased processing speed in ASD participants.

Electronic supplementary material

The online version of this article (10.1186/s13229-018-0245-1) contains supplementary material, which is available to authorized users.

Genomic kinship construction to enhance genetic analyses in the human connectome project data

Imaging genetic analyses quantify genetic control over quantitative measurements of brain structure and function using coefficients of relationship (CR) that code the degree of shared genetics between subjects. CR can be inferred through self-reported relatedness or calculated empirically using genome-wide SNP scans. We hypothesized that empirical CR provides a more accurate assessment of shared genetics than self-reported relatedness. We tested this in 1,046 participants of the Human Connectome Project (HCP) (480 M/566 F) recruited from the Missouri twin registry. We calculated the heritability for 17 quantitative traits drawn from four categories (brain diffusion and structure, cognition, and body physiology) documented by the HCP. We compared the heritability and genetic correlation estimates calculated using self-reported and empirical CR methods Kinship-based INference for GWAS (KING) and weighted allelic correlation (WAC). The polygenetic nature of traits was assessed by calculating the empirical CR from chromosomal SNP sets. The heritability estimates based on whole-genome empirical CR were higher but remained significantly correlated (r ∼0.9) with those obtained using self-reported values. Population stratification in the HCP sample has likely influenced the empirical CR calculations and biased heritability estimates. Heritability values calculated using empirical CR for chromosomal SNP sets were significantly correlated with the chromosomal length (r 0.7) suggesting a polygenic nature for these traits. The chromosomal heritability patterns were correlated among traits from the same knowledge domains; among traits with significant genetic correlations; and among traits sharing biological processes, without being genetically related. The pedigree structures generated in our analyses are available online as a web-based calculator (www.solar-eclipse-genetics.org/HCP).

Diffusion time dependence of microstructural parameters in fixed spinal cord

Biophysical modelling of diffusion MRI is necessary to provide specific microstructural tissue properties. However, estimating model parameters from data with limited diffusion gradient strength, such as clinical scanners, has proven unreliable due to a shallow optimization landscape. On the other hand, estimation of diffusion kurtosis (DKI) parameters is more robust, and its parameters may be connected to microstructural parameters, given an appropriate biophysical model. However, it was previously shown that this procedure still does not provide sufficient information to uniquely determine all model parameters. In particular, a parameter degeneracy related to the relative magnitude of intra-axonal and extra-axonal diffusivities remains. Here we develop a model of diffusion in white matter including axonal dispersion and demonstrate stable estimation of all model parameters from DKI in fixed pig spinal cord. By employing the recently developed fast axisymmetric DKI, we use stimulated echo acquisition mode to collect data over a two orders of magnitude diffusion time range with very narrow diffusion gradient pulses, enabling finely resolved measurements of diffusion time dependence of both net diffusion and kurtosis metrics, as well as model intra- and extra-axonal diffusivities, and axonal dispersion. Our results demonstrate substantial time dependence of all parameters except volume fractions, and the additional time dimension provides support for intra-axonal diffusivity to be larger than extra-axonal diffusivity in spinal cord white matter, although not unambiguously. We compare our findings for the time-dependent compartmental diffusivities to predictions from effective medium theory with reasonable agreement.

The development of brain white matter microstructure

Throughout infancy, childhood, and adolescence, our brains undergo remarkable changes. Processes including myelination and synaptogenesis occur rapidly across the first 2-3 years of life, and ongoing brain remodeling continues into young adulthood. Studies have sought to characterize the patterns of structural brain development, and early studies predominately relied upon gross anatomical measures of brain structure, morphology, and organization. MRI offers the ability to characterize and quantify a range of microstructural aspects of brain tissue that may be more closely related to fundamental neurodevelopmental processes. Techniques such as diffusion, magnetization transfer, relaxometry, and myelin water imaging provide insight into changing cyto- and myeloarchitecture, neuronal density, and structural connectivity. In this review, we focus on the growing body of literature exploiting these MRI techniques to better understand the microstructural changes that occur in brain white matter during maturation. Our review focuses on studies of normative brain development from birth to early adulthood (∼25 years), and places particular emphasis on longitudinal studies and newer techniques that are being used to study microstructural white matter development. All imaging methods demonstrate consistent, rapid microstructural white matter development over the first 3 years of life, suggesting increased myelination and axonal packing. Diffusion studies clearly demonstrate continued white matter maturation during later childhood and adolescence, though the lack of consistent findings in other modalities suggests changes may be mainly due to axonal packing. An emerging literature details differential microstructural development in boys and girls, and connects developmental trajectories to cognitive abilities, behaviour, and/or environmental factors, though the nature of these relationships remains unclear. Future research will need to focus on newer imaging techniques and longitudinal studies to provide more detailed information about microstructural white matter development, particularly in the childhood years.

Prenatal stress and the developing brain: Risks for neurodevelopmental disorders

The prenatal period is increasingly considered as a crucial target for the primary prevention of neurodevelopmental and psychiatric disorders. Understanding their pathophysiological mechanisms remains a great challenge. Our review reveals new insights from prenatal brain development research, involving (epi)genetic research, neuroscience, recent imaging techniques, physical modeling, and computational simulation studies. Studies examining the effect of prenatal exposure to maternal distress on offspring brain development, using brain imaging techniques, reveal effects at birth and up into adulthood. Structural and functional changes are observed in several brain regions including the prefrontal, parietal, and temporal lobes, as well as the cerebellum, hippocampus, and amygdala. Furthermore, alterations are seen in functional connectivity of amygdalar-thalamus networks and in intrinsic brain networks, including default mode and attentional networks. The observed changes underlie offspring behavioral, cognitive, emotional development, and susceptibility to neurodevelopmental and psychiatric disorders. It is concluded that used brain measures have not yet been validated with regard to sensitivity, specificity, accuracy, or robustness in predicting neurodevelopmental and psychiatric disorders. Therefore, more prospective long-term longitudinal follow-up studies starting early in pregnancy should be carried out, in order to examine brain developmental measures as mediators in mediating the link between prenatal stress and offspring behavioral, cognitive, and emotional problems and susceptibility for disorders.

Modeling white matter tract integrity in aging with diffusional kurtosis imaging

Myelin breakdown and neural fiber loss occur in aging. This study used white matter tract integrity metrics derived from biophysical modeling using Diffusional Kurtosis Imaging to assess loss of myelin (i.e., extraaxonal diffusivity, radial direction, De,⊥) and axonal density (i.e., axonal water fraction) in cognitively unimpaired older adults. Tract-based spatial statistics and region of interest analyses sought to identify ontogenic differences and age-related changes in white matter tracts using cross-sectional and longitudinal data analyzed with general linear and mixed-effects models. In addition to pure diffusion parameters (i.e., fractional anisotropy, mean diffusivity, mean kurtosis), we found that white matter tract integrity metrics significantly differentiated early- from late-myelinating tracts, correlated with age in spatially distinct regions, and identified primarily extraaxonal changes over time. Percent metric changes were |0.3-0.9|% and |0.0-1.9|% per year using cross-sectional data and longitudinal data, respectively. There was accelerated decline in some late- versus early-myelinating tracts in older age. These results demonstrate that these metrics may inform further study of the transition from age-related changes to neurodegenerative decline.

Relevance of time-dependence for clinically viable diffusion imaging of the spinal cord

Purpose

Time‐dependence is a key feature of the diffusion‐weighted (DW) signal, knowledge of which informs biophysical modelling. Here, we study time‐dependence in the human spinal cord, as its axonal structure is specific and different from the brain.

Methods

We run Monte Carlo simulations using a synthetic model of spinal cord white matter (WM) (large axons), and of brain WM (smaller axons). Furthermore, we study clinically feasible multi‐shell DW scans of the cervical spinal cord (b = 0; b = 711 s mm⁻²; b = 2855 s mm⁻²), obtained using three diffusion times (Δ of 29, 52 and 76 ms) from three volunteers.

Results

Both intra‐/extra‐axonal perpendicular diffusivities and kurtosis excess show time‐dependence in our synthetic spinal cord model. This time‐dependence is reflected mostly in the intra‐axonal perpendicular DW signal, which also exhibits strong decay, unlike our brain model. Time‐dependence of the total DW signal appears detectable in the presence of noise in our synthetic spinal cord model, but not in the brain. In WM in vivo, we observe time‐dependent macroscopic and microscopic diffusivities and diffusion kurtosis, NODDI and two‐compartment SMT metrics. Accounting for large axon calibers improves fitting of multi‐compartment models to a minor extent.

Conclusions

Time‐dependence of clinically viable DW MRI metrics can be detected in vivo in spinal cord WM, thus providing new opportunities for the non‐invasive estimation of microstructural properties. The time‐dependence of the perpendicular DW signal may feature strong intra‐axonal contributions due to large spinal axon caliber. Hence, a popular model known as “stick” (zero‐radius cylinder) may be sub‐optimal to describe signals from the largest spinal axons.

Filter-probe diffusion imaging improves spinal cord injury outcome prediction

OBJECTIVE

Diffusion-weighted imaging (DWI) is a powerful tool for investigating spinal cord injury (SCI), but has limited specificity for axonal damage, which is the most predictive feature of long-term functional outcome. In this study, a technique designed to detect acute axonal injury, filter-probe double diffusion encoding (FP-DDE), is compared with standard DWI for predicting long-term functional and cellular outcomes.

METHODS

This study extends FP-DDE to predict long-term functional and histological outcomes in a rat SCI model of varying severities (n = 58). Using a 9.4T magnetic resonance imaging (MRI) system, a whole-cord FP-DDE spectroscopic voxel was acquired in 3 minutes at the lesion site and compared to DWI at 48 hours postinjury. Relationships with chronic (30-day) locomotor and histological outcomes were evaluated with linear regression.

RESULTS

The FP-DDE measure of parallel diffusivity (ADC|| ) was significantly related to chronic hind limb locomotor functional outcome (R2  = 0.63, p < 0.0001), and combining this measurement with acute functional scores demonstrated prognostic benefit versus functional testing alone (p = 0.0007). Acute ADC|| measurements were also more closely related to the number of injured axons measured 30 days after the injury than standard DWI. Furthermore, acute FP-DDE images showed a clear and easily interpretable pattern of injury that closely corresponded with chronic MRI and histology observations.

INTERPRETATION

Collectively, these results demonstrate FP-DDE benefits from greater specificity for acute axonal damage in predicting functional and histological outcomes with rapid acquisition and fully automated analysis, improving over standard DWI. FP-DDE is a promising technique compatible with clinical settings, with potential research and clinical applications for evaluation of spinal cord pathology. Ann Neurol 2018;83:37-50.

Characterizing intra-axonal water diffusion with direction-averaged triple diffusion encoding MRI

For large diffusion weightings, the direction-averaged diffusion MRI (dMRI) signal from white matter is typically dominated by the contribution of water confined to axons. This fact can be exploited to characterize intra-axonal diffusion properties, which may be valuable for interpreting the biophysical meaning of diffusion changes associated with pathology. However, using just the classic Stejskal-Tanner pulse sequence, it has proven challenging to obtain reliable estimates for both the intrinsic intra-axonal diffusivity and the intra-axonal water fraction. Here we propose to apply a modification of the Stejskal-Tanner sequence designed for achieving such estimates. The key feature of the sequence is the addition of a set of extra diffusion encoding gradients that are orthogonal to the direction of the primary gradients, which corresponds to a specific type of triple diffusion encoding (TDE) MRI sequence. Given direction-averaged dMRI data for this TDE sequence, it is shown how the intra-axonal diffusivity and the intra-axonal water fraction can be determined by applying simple, analytic formulae. The method is illustrated with numerical simulations, which suggest that it should be accurate for b-values of about 4000 s/mm² or higher.

A biomimetic tumor tissue phantom for validating diffusion-weighted MRI measurements

PURPOSE

To develop a biomimetic tumor tissue phantom which more closely reflects water diffusion in biological tissue than previously used phantoms, and to evaluate the stability of the phantom and its potential as a tool for validating diffusion-weighted (DW) MRI measurements.

METHODS

Coaxial-electrospraying was used to generate micron-sized hollow polymer spheres, which mimic cells. The bulk structure was immersed in water, providing a DW-MRI phantom whose apparent diffusion coefficient (ADC) and microstructural properties were evaluated over a period of 10 months. Independent characterization of the phantom's microstructure was performed using scanning electron microscopy (SEM). The repeatability of the construction process was investigated by generating a second phantom, which underwent high resolution synchrotron-CT as well as SEM and MR scans.

RESULTS

ADC values were stable (coefficients of variation (CoVs) < 5%), and varied with diffusion time, with average values of 1.44 ± 0.03 µm2 /ms (Δ = 12 ms) and 1.20 ± 0.05 µm2 /ms (Δ = 45 ms). Microstructural parameters showed greater variability (CoVs up to 13%), with evidence of bias in sphere size estimates. Similar trends were observed in the second phantom.

CONCLUSION

A novel biomimetic phantom has been developed and shown to be stable over 10 months. It is envisaged that such phantoms will be used for further investigation of microstructural models relevant to characterizing tumor tissue, and may also find application in evaluating acquisition protocols and comparing DW-MRI-derived biomarkers obtained from different scanners at different sites. Magn Reson Med 80:147-158, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Empirical reproducibility, sensitivity, and optimization of acquisition protocol, for Neurite Orientation Dispersion and Density Imaging using AMICO

Neurite Orientation Dispersion and Density Imaging (NODDI) has been gaining prominence for estimating multiple diffusion compartments from MRI data acquired in a clinically feasible time. To establish a pathway for adoption of NODDI in clinical studies, it is important to understand the sensitivity and reproducibility of NODDI metrics on empirical data in the context of acquisition protocol and brain anatomy. Previous studies addressed reproducibility across the 3 T scanners and within session and between subject reproducibility at 1.5 T and 3 T. However, empirical reproducibility on the performance of NODDI metrics based on b-value and diffusion-sensitized directions has not yet been addressed. In this study, we investigate a high angular resolution dataset with 11 repeats of a study with five b-values shells (1000, 1500, 2000, 2500 and 3000 s/mm2) and 96 directions per shell on a single subject. We validated the findings with a dataset from second subject with 10 repeats and 3 b-value shells (1000, 2000, 3000 s/mm2). The NODDI model was estimated using Accelerated Microstructure Imaging via Convex Optimization (AMICO) for different b-values and gradient directions on two-shell High Angular Resolution Density Imaging (HARDI) data fixing the lower shell at b = 1000 s/mm2. NODDI model applied to all acquired imaging data was used as a baseline gold standard for comparison. Additionally, we characterize orientation dispersion index (ODI) reproducibility using single-shell data. The experimental findings confirmed the sensitivity of intracellular volume fraction (Vic) with the choice of outer shell b-value more than with the choice of gradient directions. On the other hand, ODI is more sensitive to the number of gradient directions compared to b-value selection. Single-shell results for ODI are more comparable to 2-shell data at lower b-values than higher b-values. Recommended settings by region of interest and acquisition time are reported for the researchers considering using NODDI in human studies and/or comparing results across acquisition protocols.

White matter biomarkers from diffusion MRI

As part of an issue celebrating 2 decades of Joseph Ackerman editing the Journal of Magnetic Resonance, this paper reviews recent progress in one of the many areas in which Ackerman and his lab has made significant contributions: NMR measurement of diffusion in biological media, specifically in brain tissue. NMR diffusion signals display exquisite sensitivity to tissue microstructure, and have the potential to offer quantitative and specific information on the cellular scale orders of magnitude below nominal image resolution when combined with biophysical modeling. Here, I offer a personal perspective on some recent advances in diffusion imaging, from diffusion kurtosis imaging to microstructural modeling, and the connection between the two. A new result on the estimation accuracy of axial and radial kurtosis with axially symmetric DKI is presented. I moreover touch upon recently suggested generalized diffusion sequences, promising to offer independent microstructural information. We discuss the need and some methods for validation, and end with an outlook on some promising future directions.