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

On modeling

Mapping tissue microstructure with MRI holds great promise as a noninvasive window into tissue organization at the cellular level. Having originated within the realm of diffusion NMR in the late 1970s, this field is experiencing an exponential growth in the number of publications. At the same time, model-based approaches are also increasingly incorporated into advanced MRI acquisition and reconstruction techniques. However, after about two decades of intellectual and financial investment, microstructural mapping has yet to find a single commonly accepted clinical application. Here, we suggest that slow progress in clinical translation may signify unresolved fundamental problems. We outline such problems and related practical pitfalls, as well as review strategies for developing and validating tissue microstructure models, to provoke a discussion on how to bridge the gap between our scientific aspirations and the clinical reality. We argue for recalibrating the efforts of our community toward a more systematic focus on fundamental research aimed at identifying relevant degrees of freedom affecting the measured MR signal. Such a focus is essential for realizing the truly revolutionary potential of noninvasive three-dimensional in vivo microstructural mapping.

Different patterns of cortical maturation before and after 38 weeks gestational age demonstrated by diffusion MRI in vivo

Human cortical development during the third trimester is characterised by macro- and microstructural changes which are reflected in alterations in diffusion MRI (dMRI) measures, with significant decreases in cortical mean diffusivity (MD) and fractional anisotropy (FA). This has been interpreted as reflecting increased cellular density and dendritic arborisation. However, the fall in FA stops abruptly at 38 weeks post-menstrual age (PMA), and then tends to plateau, while MD continues to fall, suggesting a more complex picture and raising the hypothesis that after this age development is dominated by continuing increase in neural and organelle density rather than alterations in the geometry of dendritic trees. To test this, we used neurite orientation dispersion and density imaging (NODDI), acquiring multi-shell, high angular resolution dMRI and measures of cortical volume and mean curvature in 99 preterm infants scanned between 25 and 47 weeks PMA. We predicted that increased neurite and organelle density would be reflected in increases in neurite density index (NDI), while a relatively unchanging geometrical structure would be associated with constant orientation dispersion index (ODI). As dendritic arborisation is likely to be one of the drivers of gyrification, we also predicted that measures of cortical volume and curvature would correlate with ODI and show slower growth after 38 weeks.

We observed a decrease of MD throughout the period, while cortical FA decreased from 25 to 38 weeks PMA and then increased. ODI increased up to 38 weeks and then plateaued, while NDI rose after 38 weeks. The evolution of ODI correlated with cortical volume and curvature. Regional analysis of cortical microstructure revealed a heterogenous pattern with increases in FA and NDI after 38 weeks confined to primary motor and sensory regions. These results support the interpretation that cortical development between 25 and 38 weeks PMA shows a predominant increase in dendritic arborisation and neurite growth, while between 38 and 47 weeks PMA it is dominated by increasing cellular and organelle density.

•

DTI and NODDI cortical measures between 25 and 47 weeks GA

•

Early cortical changes consistent with dendritic arborisation and neurite growth

•

After 38 weeks cortical changes consistent with increasing cellular density

•

Cortical curvature evolves in parallel with dendritic arborisation

Future Directions for Examination of Brain Networks in Neurodevelopmental Disorders

Neurodevelopmental disorders are associated with atypical development and maturation of brain networks. A recent focus on human connectomics research and the growing popularity of open science initiatives has created the ideal climate in which to make real progress toward understanding the neurobiology of disorders affecting youth. Here we outline future directions for neuroscience researchers examining brain networks in neurodevelopmental disorders, highlighting gaps in the current literature. We emphasize the importance of leveraging large neuroimaging and phenotypic data sets recently made available to the research community, and we suggest specific novel methodological approaches, including analysis of brain dynamics and structural connectivity, that have the potential to produce the greatest clinical insight. Transdiagnostic approaches will also become increasingly necessary as the Research Domain Criteria framework put forth by the National Institute of Mental Health permeates scientific discourse. During this exciting era of big data and increased computational sophistication of analytic tools, the possibilities for significant advancement in understanding neurodevelopmental disorders are limitless.

Delineation of early brain development from fetuses to infants with diffusion MRI and beyond

Dynamic macrostructural and microstructural changes take place from the mid-fetal stage to 2 years after birth. Delineating structural changes of the brain during early development provides new insights into the complicated processes of both typical development and the pathological mechanisms underlying various psychiatric and neurological disorders including autism, attention deficit hyperactivity disorder and schizophrenia. Decades of histological studies have identified strong spatial and functional maturation gradients in human brain gray and white matter. The recent improvements in magnetic resonance imaging (MRI) techniques, especially diffusion MRI (dMRI), relaxometry imaging, and magnetization transfer imaging (MTI) have provided unprecedented opportunities to non-invasively quantify and map the early developmental changes at whole brain and regional levels. Here, we review the recent advances in understanding early brain structural development during the second half of gestation and the first two postnatal years using modern MR techniques. Specifically, we review studies that delineate the emergence and microstructural maturation of white matter tracts, as well as dynamic mapping of inhomogeneous cortical microstructural organization unique to fetuses and infants. These imaging studies converge into maturational curves of MRI measurements that are distinctive across different white matter tracts and cortical regions. Furthermore, contemporary models offering biophysical interpretations of the dMRI-derived measurements are illustrated to infer the underlying microstructural changes. Collectively, this review summarizes findings that contribute to charting spatiotemporally heterogeneous gray and white matter structural development, offering MRI-based biomarkers of typical brain development and setting the stage for understanding aberrant brain development in neurodevelopmental disorders.

Working Memory And Brain Tissue Microstructure: White Matter Tract Integrity Based On Multi-Shell Diffusion MRI

Working memory is a complex cognitive process at the intersection of sensory processing, learning, and short-term memory and also has a general executive attention component. Impaired working memory is associated with a range of neurological and psychiatric disorders, but very little is known about how working memory relates to underlying white matter (WM) microstructure. In this study, we investigate the association between WM microstructure and performance on working memory tasks in healthy adults (right-handed, native English speakers). We combine compartment specific WM tract integrity (WMTI) metrics derived from multi-shell diffusion MRI as well as diffusion tensor/kurtosis imaging (DTI/DKI) metrics with Wechsler Adult Intelligence Scale-Fourth Edition (WAIS-IV) subtests tapping auditory working memory. WMTI is a novel tool that helps us describe the microstructural characteristics in both the intra- and extra-axonal environments of WM such as axonal water fraction (AWF), intra-axonal diffusivity, extra-axonal axial and radial diffusivities, allowing a more biophysical interpretation of WM changes. We demonstrate significant positive correlations between AWF and letter-number sequencing (LNS), suggesting that higher AWF with better performance on complex, more demanding auditory working memory tasks goes along with greater axonal volume and greater myelination in specific regions, causing efficient and faster information process.

Imaging structural and functional brain development in early childhood

In humans, the period from term birth to ∼2 years of age is characterized by rapid and dynamic brain development and plays an important role in cognitive development and risk of disorders such as autism and schizophrenia. Recent imaging studies have begun to delineate the growth trajectories of brain structure and function in the first years after birth and their relationship to cognition and risk of neuropsychiatric disorders. This Review discusses the development of grey and white matter and structural and functional networks, as well as genetic and environmental influences on early-childhood brain development. We also discuss initial evidence regarding the usefulness of early imaging biomarkers for predicting cognitive outcomes and risk of neuropsychiatric disorders.

Neurite imaging reveals microstructural variations in human cerebral cortical gray matter

We present distinct patterns of neurite distribution in the human cerebral cortex using diffusion magnetic resonance imaging (MRI). We analyzed both high-resolution structural (T1w and T2w images) and diffusion MRI data in 505 subjects from the Human Connectome Project. Neurite distributions were evaluated using the neurite orientation dispersion and density imaging (NODDI) model, optimized for gray matter, and mapped onto the cortical surface using a method weighted towards the cortical mid-thickness to reduce partial volume effects. The estimated neurite density was high in both somatosensory and motor areas, early visual and auditory areas, and middle temporal area (MT), showing a strikingly similar distribution to myelin maps estimated from the T1w/T2w ratio. The estimated neurite orientation dispersion was particularly high in early sensory areas, which are known for dense tangential fibers and are classified as granular cortex by classical anatomists. Spatial gradients of these cortical neurite properties revealed transitions that colocalize with some areal boundaries in a recent multi-modal parcellation of the human cerebral cortex, providing mutually supportive evidence. Our findings indicate that analyzing the cortical gray matter neurite morphology using diffusion MRI and NODDI provides valuable information regarding cortical microstructure that is related to but complementary to myeloarchitecture.

Multidimensional correlation of nuclear relaxation rates and diffusion tensors for model-free investigations of heterogeneous anisotropic porous materials

Despite their widespread use in non-invasive studies of porous materials, conventional MRI methods yield ambiguous results for microscopically heterogeneous materials such as brain tissue. While the forward link between microstructure and MRI observables is well understood, the inverse problem of separating the signal contributions from different microscopic pores is notoriously difficult. Here, we introduce an experimental protocol where heterogeneity is resolved by establishing 6D correlations between the individual values of isotropic diffusivity, diffusion anisotropy, orientation of the diffusion tensor, and relaxation rates of distinct populations. Such procedure renders the acquired signal highly specific to the sample's microstructure, and allows characterization of the underlying pore space without prior assumptions on the number and nature of distinct microscopic environments. The experimental feasibility of the suggested method is demonstrated on a sample designed to mimic the properties of nerve tissue. If matched to the constraints of whole body scanners, this protocol could allow for the unconstrained determination of the different types of tissue that compose the living human brain.

Post-mortem inference of the human hippocampal connectivity and microstructure using ultra-high field diffusion MRI at 11.7 T

The human hippocampus plays a key role in memory management and is one of the first structures affected by Alzheimer's disease. Ultra-high magnetic resonance imaging provides access to its inner structure in vivo. However, gradient limitations on clinical systems hinder access to its inner connectivity and microstructure. A major target of this paper is the demonstration of diffusion MRI potential, using ultra-high field (11.7 T) and strong gradients (750 mT/m), to reveal the extra- and intra-hippocampal connectivity in addition to its microstructure. To this purpose, a multiple-shell diffusion-weighted acquisition protocol was developed to reach an ultra-high spatio-angular resolution with a good signal-to-noise ratio. The MRI data set was analyzed using analytical Q-Ball Imaging, Diffusion Tensor Imaging (DTI), and Neurite Orientation Dispersion and Density Imaging models. High Angular Resolution Diffusion Imaging estimates allowed us to obtain an accurate tractography resolving more complex fiber architecture than DTI models, and subsequently provided a map of the cross-regional connectivity. The neurite density was akin to that found in the histological literature, revealing the three hippocampal layers. Moreover, a gradient of connectivity and neurite density was observed between the anterior and the posterior part of the hippocampus. These results demonstrate that ex vivo ultra-high field/ultra-high gradients diffusion-weighted MRI allows the mapping of the inner connectivity of the human hippocampus, its microstructure, and to accurately reconstruct elements of the polysynaptic intra-hippocampal pathway using fiber tractography techniques at very high spatial/angular resolutions.

Deciphering the microstructure of hippocampal subfields with in vivo DTI and NODDI: Applications to experimental multiple sclerosis

The hippocampus contains distinct populations of neurons organized into separate anatomical subfields and layers with differential vulnerability to pathological mechanisms. The ability of in vivo neuroimaging to pinpoint regional vulnerability is especially important for better understanding of hippocampal pathology at the early stage of neurodegenerative disorders and for monitoring future therapeutic strategies. This is the case for instance in multiple sclerosis whose neurodegenerative component can affect the hippocampus from the early stage. We challenged the capacity of two models, i.e. the classical diffusion tensor imaging (DTI) model and the neurite orientation dispersion and density imaging (NODDI) model, to compute quantitative diffusion MRI that could capture microstructural alterations in the individual hippocampal layers of experimental-autoimmune encephalomyelitis (EAE) mice, the animal model of multiple sclerosis. To achieve this, the hippocampal anatomy of a healthy mouse brain was first explored ex vivo with high resolution DTI and NODDI. Then, 18 EAE mice and 18 control mice were explored 20 days after immunization with in vivo diffusion MRI prior to sacrifice for the histological quantification of neurites and glial markers in each hippocampal layer. Fractional anisotropy (FA), axial diffusivity (AD), radial diffusivity (RD) and mean diffusivity (MD) maps were computed from the DTI model while the orientation dispersion index (ODI), the neurite density index (NDI) and the volume fraction of isotropic diffusivity (isoVF) maps were computed from the NODDI model. We first showed in control mice that color-coded FA and ODI maps can delineate three main hippocampal layers. The quantification of FA, AD, RD, MD, ODI, NDI and isoVF presented differences within these 3 layers, especially within the molecular layer of the dentate gyrus which displayed a specific signature based on a combination of AD (or MD), ODI and NDI. Then, the comparison between EAE and control mice showed a decrease of AD (p = 0.036) and of MD (p = 0.033) selectively within the molecular layer of EAE mice while NODDI indices did not present any difference between EAE and control mice in any layer. Histological analyses confirmed the differential vulnerability of the molecular layer of EAE mice that exhibited decreased dendritic length and decreased dendritic complexity together with activated microglia. Dendritic length and intersections within the molecular layer were independent contributors to the observed decrease of AD (R² = 0.37 and R² = 0.40, p < 0.0001) and MD (R² = 0.41 and R² = 0.42, p < 0.0001). We therefore identified that NODDI maps can help to highlight the internal microanatomy of the hippocampus but NODDI still presents limitations in grey matter as it failed to capture selective dendritic alterations occurring at early stages of a neurodegenerative disease such as multiple sclerosis, whereas DTI maps were significantly altered.

Fixel-based analysis reveals alterations is brain microstructure and macrostructure of preterm-born infants at term equivalent age

Preterm birth causes significant disruption in ongoing brain development, frequently resulting in adverse neurodevelopmental outcomes. Brain imaging using diffusion MRI may provide valuable insight into microstructural properties of the developing brain. The aim of this study was to establish whether the recently introduced fixel-based analysis method, with its associated measures of fibre density (FD), fibre bundle cross-section (FC), and fibre density and bundle cross-section (FDC), is suitable for the investigation of the preterm infant brain at term equivalent age. High-angular resolution diffusion weighted images (HARDI) of 55 preterm-born infants and 20 term-born infants, scanned around term-equivalent age, were included in this study (3 T, 64 directions, b = 2000 s/mm²). Postmenstrual age at the time of MRI, and intracranial volume (FC and FDC only), were identified as confounding variables. Gestational age at birth was correlated with all fixel measures in the splenium of the corpus callosum. Compared to term-born infants, preterm infants showed reduced FD, FC, and FDC in a number of regions, including the corpus callosum, anterior commissure, cortico-spinal tract, optic radiations, and cingulum. Preterm infants with minimal macroscopic brain abnormality showed more extensive reductions than preterm infants without any macroscopic brain abnormality; however, little differences were observed between preterm infants with no and with minimal brain abnormality. FC showed significant reductions in preterm versus term infants outside regions identified with FD and FDC, highlighting the complementary role of these measures. Fixel-based analysis identified both microstructural and macrostructural abnormalities in preterm born infants, providing a more complete picture of early brain development than previous diffusion tensor imaging (DTI) based approaches.

•

Gestational age at birth associated with measurements in corpus callosum splenium.

•

Preterms without macroscopic brain abnormality show differences to term infants.

•

No differences between preterms with minimal versus without abnormality detected.

NODDI-DTI: Estimating Neurite Orientation and Dispersion Parameters from a Diffusion Tensor in Healthy White Matter

The NODDI-DTI signal model is a modification of the NODDI signal model that formally allows interpretation of standard single-shell DTI data in terms of biophysical parameters in healthy human white matter (WM). The NODDI-DTI signal model contains no CSF compartment, restricting application to voxels without CSF partial-volume contamination. This modification allowed derivation of analytical relations between parameters representing axon density and dispersion, and DTI invariants (MD and FA) from the NODDI-DTI signal model. These relations formally allow extraction of biophysical parameters from DTI data. NODDI-DTI parameters were estimated by applying the proposed analytical relations to DTI parameters estimated from the first shell of data, and compared to parameters estimated by fitting the NODDI-DTI model to both shells of data (reference dataset) in the WM of 14 in vivo diffusion datasets recorded with two different protocols, and in simulated data. The first two datasets were also fit to the NODDI-DTI model using only the first shell (as for DTI) of data. NODDI-DTI parameters estimated from DTI, and NODDI-DTI parameters estimated by fitting the model to the first shell of data gave similar errors compared to two-shell NODDI-DTI estimates. The simulations showed the NODDI-DTI method to be more noise-robust than the two-shell fitting procedure. The NODDI-DTI method gave unphysical parameter estimates in a small percentage of voxels, reflecting voxelwise DTI estimation error or NODDI-DTI model invalidity. In the course of evaluating the NODDI-DTI model, it was found that diffusional kurtosis strongly biased DTI-based MD values, and so, making assumptions based on healthy WM, a novel heuristic correction requiring only DTI data was derived and used to mitigate this bias. Since validations were only performed on healthy WM, application to grey matter or pathological WM would require further validation. Our results demonstrate NODDI-DTI to be a promising model and technique to interpret restricted datasets acquired for DTI analysis in healthy white matter with greater biophysical specificity, though its limitations must be borne in mind.

Optimizing Filter-Probe Diffusion Weighting in the Rat Spinal Cord for Human Translation

Diffusion tensor imaging (DTI) is a promising biomarker of spinal cord injury (SCI). In the acute aftermath, DTI in SCI animal models consistently demonstrates high sensitivity and prognostic performance, yet translation of DTI to acute human SCI has been limited. In addition to technical challenges, interpretation of the resulting metrics is ambiguous, with contributions in the acute setting from both axonal injury and edema. Novel diffusion MRI acquisition strategies such as double diffusion encoding (DDE) have recently enabled detection of features not available with DTI or similar methods. In this work, we perform a systematic optimization of DDE using simulations and an in vivo rat model of SCI and subsequently implement the protocol to the healthy human spinal cord. First, two complementary DDE approaches were evaluated using an orientationally invariant or a filter-probe diffusion encoding approach. While the two methods were similar in their ability to detect acute SCI, the filter-probe DDE approach had greater predictive power for functional outcomes. Next, the filter-probe DDE was compared to an analogous single diffusion encoding (SDE) approach, with the results indicating that in the spinal cord, SDE provides similar contrast with improved signal to noise. In the SCI rat model, the filter-probe SDE scheme was coupled with a reduced field of view (rFOV) excitation, and the results demonstrate high quality maps of the spinal cord without contamination from edema and cerebrospinal fluid, thereby providing high sensitivity to injury severity. The optimized protocol was demonstrated in the healthy human spinal cord using the commercially-available diffusion MRI sequence with modifications only to the diffusion encoding directions. Maps of axial diffusivity devoid of CSF partial volume effects were obtained in a clinically feasible imaging time with a straightforward analysis and variability comparable to axial diffusivity derived from DTI. Overall, the results and optimizations describe a protocol that mitigates several difficulties with DTI of the spinal cord. Detection of acute axonal damage in the injured or diseased spinal cord will benefit the optimized filter-probe diffusion MRI protocol outlined here.

Multipulse sodium magnetic resonance imaging for multicompartment quantification: Proof-of-concept

We present a feasibility study of sodium quantification in a multicompartment model of the brain using sodium (²³Na) magnetic resonance imaging. The proposed method is based on a multipulse sequence acquisition and simulation at 7 T, which allows to differentiate the ²³Na signals emanating from three compartments in human brain in vivo: intracellular (compartment 1), extracellular (compartment 2), and cerebrospinal fluid (compartment 3). The intracellular sodium concentration C₁ and the volume fractions α₁, α₂, and α₃ of all respective three brain compartments can be estimated. Simulations of the sodium spin 3/2 dynamics during a 15-pulse sequence were used to optimize the acquisition sequence by minimizing the correlation between the signal evolutions from the three compartments. The method was first tested on a three-compartment phantom as proof-of-concept. Average values of the ²³Na quantifications in four healthy volunteer brains were α₁ = 0.54 ± 0.01, α₂ = 0.23 ± 0.01, α₃ = 1.03 ± 0.01, and C₁ = 23 ± 3 mM, which are comparable to the expected physiological values \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{\alpha }}}_{{\bf{1}}}^{{\boldsymbol{theory}}}$$\end{document}α1theory ∼ 0.6, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{\alpha }}}_{{\bf{2}}}^{{\boldsymbol{theory}}}$$\end{document}α2theory ∼ 0.2, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{\alpha }}}_{{\bf{3}}}^{{\boldsymbol{theory}}}$$\end{document}α3theory ∼ 1, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{C}}}_{{\bf{1}}}^{{\boldsymbol{theory}}}$$\end{document}C1theory ∼ 10–30 mM. The proposed method may allow a quantitative assessment of the metabolic role of sodium ions in cellular processes and their malfunctions in brain in vivo.

Design and validation of diffusion MRI models of white matter

Diffusion MRI is arguably the method of choice for characterizing white matter microstructure in vivo. Over the typical duration of diffusion encoding, the displacement of water molecules is conveniently on a length scale similar to that of the underlying cellular structures. Moreover, water molecules in white matter are largely compartmentalized which enables biologically-inspired compartmental diffusion models to characterize and quantify the true biological microstructure. A plethora of white matter models have been proposed. However, overparameterization and mathematical fitting complications encourage the introduction of simplifying assumptions that vary between different approaches. These choices impact the quantitative estimation of model parameters with potential detriments to their biological accuracy and promised specificity. First, we review biophysical white matter models in use and recapitulate their underlying assumptions and realms of applicability. Second, we present up-to-date efforts to validate parameters estimated from biophysical models. Simulations and dedicated phantoms are useful in assessing the performance of models when the ground truth is known. However, the biggest challenge remains the validation of the "biological accuracy" of estimated parameters. Complementary techniques such as microscopy of fixed tissue specimens have facilitated direct comparisons of estimates of white matter fiber orientation and densities. However, validation of compartmental diffusivities remains challenging, and complementary MRI-based techniques such as alternative diffusion encodings, compartment-specific contrast agents and metabolites have been used to validate diffusion models. Finally, white matter injury and disease pose additional challenges to modeling, which are also discussed. This review aims to provide an overview of the current state of models and their validation and to stimulate further research in the field to solve the remaining open questions and converge towards consensus.

Tissue microstructure estimation using a deep network inspired by a dictionary-based framework

Diffusion magnetic resonance imaging (dMRI) captures the anisotropic pattern of water displacement in the neuronal tissue and allows noninvasive investigation of the complex tissue microstructure. A number of biophysical models have been proposed to relate the tissue organization with the observed diffusion signals, so that the tissue microstructure can be inferred. The Neurite Orientation Dispersion and Density Imaging (NODDI) model has been a popular choice and has been widely used for many neuroscientific studies. It models the diffusion signal with three compartments that are characterized by distinct diffusion properties, and the parameters in the model describe tissue microstructure. In NODDI, these parameters are estimated in a maximum likelihood framework, where the nonlinear model fitting is computationally intensive. Therefore, efforts have been made to develop efficient and accurate algorithms for NODDI microstructure estimation, which is still an open problem. In this work, we propose a deep network based approach that performs end-to-end estimation of NODDI microstructure, which is named Microstructure Estimation using a Deep Network (MEDN). MEDN comprises two cascaded stages and is motivated by the AMICO algorithm, where the NODDI microstructure estimation is formulated in a dictionary-based framework. The first stage computes the coefficients of the dictionary. It resembles the solution to a sparse reconstruction problem, where the iterative process in conventional estimation approaches is unfolded and truncated, and the weights are learned instead of predetermined by the dictionary. In the second stage, microstructure properties are computed from the output of the first stage, which resembles the weighted sum of normalized dictionary coefficients in AMICO, and the weights are also learned. Because spatial consistency of diffusion signals can be used to reduce the effect of noise, we also propose MEDN+, which is an extended version of MEDN. MEDN+ allows incorporation of neighborhood information by inserting a stage with learned weights before the MEDN structure, where the diffusion signals in the neighborhood of a voxel are processed. The weights in MEDN or MEDN+ are jointly learned from training samples that are acquired with diffusion gradients densely sampling the q-space. We performed MEDN and MEDN+ on brain dMRI scans, where two shells each with 30 gradient directions were used, and measured their accuracy with respect to the gold standard. Results demonstrate that the proposed networks outperform the competing methods.

White matter biomarkers from fast protocols using axially symmetric diffusion kurtosis imaging

White matter tract integrity (WMTI) can characterize brain microstructure in areas with highly aligned fiber bundles. Several WMTI biomarkers have now been validated against microscopy and provided promising results in studies of brain development and aging, as well as in a number of brain disorders. Currently, WMTI is mostly used in dedicated animal studies and clinical studies of slowly progressing diseases, and has not yet emerged as a routine clinical tool. To this end, a less data intensive experimental method would be beneficial by enabling high resolution validation studies, and ease clinical applications by speeding up data acquisition compared with typical diffusion kurtosis imaging (DKI) protocols utilized as part of WMTI imaging. Here, we evaluate WMTI based on recently introduced axially symmetric DKI, which has lower data demand than conventional DKI. We compare WMTI parameters derived from conventional DKI with those calculated analytically from axially symmetric DKI. We employ numerical simulations, as well as data from fixed rat spinal cord (one sample) and in vivo human (three subjects) and rat brain (four animals). Our analysis shows that analytical WMTI based on axially symmetric DKI with sparse data sets (19 images) produces WMTI metrics that correlate strongly with estimates based on traditional DKI data sets (60 images or more). We demonstrate the preclinical potential of the proposed WMTI technique in in vivo rat brain (300 μm isotropic resolution with whole brain coverage in a 1 h acquisition). WMTI parameter estimates are subject to a duality leading to two solution branches dependent on a sign choice, which is currently debated. Results from both of these branches are presented and discussed throughout our analysis. The proposed fast WMTI approach may be useful for preclinical research and e.g. clinical evaluation of patients with traumatic white matter injuries or symptoms of neurovascular or neuroinflammatory disorders.

Diffusion MRI microstructure models with in vivo human brain Connectome data: results from a multi-group comparison

A large number of mathematical models have been proposed to describe the measured signal in diffusion-weighted (DW) magnetic resonance imaging (MRI). However, model comparison to date focuses only on specific subclasses, e.g. compartment models or signal models, and little or no information is available in the literature on how performance varies among the different types of models. To address this deficiency, we organized the 'White Matter Modeling Challenge' during the International Symposium on Biomedical Imaging (ISBI) 2015 conference. This competition aimed to compare a range of different kinds of models in their ability to explain a large range of measurable in vivo DW human brain data. Specifically, we assessed the ability of models to predict the DW signal accurately for new diffusion gradients and b values. We did not evaluate the accuracy of estimated model parameters, as a ground truth is hard to obtain. We used the Connectome scanner at the Massachusetts General Hospital, using gradient strengths of up to 300 mT/m and a broad set of diffusion times. We focused on assessing the DW signal prediction in two regions: the genu in the corpus callosum, where the fibres are relatively straight and parallel, and the fornix, where the configuration of fibres is more complex. The challenge participants had access to three-quarters of the dataset and their models were ranked on their ability to predict the remaining unseen quarter of the data. The challenge provided a unique opportunity for a quantitative comparison of diverse methods from multiple groups worldwide. The comparison of the challenge entries reveals interesting trends that could potentially influence the next generation of diffusion-based quantitative MRI techniques. The first is that signal models do not necessarily outperform tissue models; in fact, of those tested, tissue models rank highest on average. The second is that assuming a non-Gaussian (rather than purely Gaussian) noise model provides little improvement in prediction of unseen data, although it is possible that this may still have a beneficial effect on estimated parameter values. The third is that preprocessing the training data, here by omitting signal outliers, and using signal-predicting strategies, such as bootstrapping or cross-validation, could benefit the model fitting. The analysis in this study provides a benchmark for other models and the data remain available to build up a more complete comparison in the future.

Mapping White Matter Microstructure in the One Month Human Brain

White matter microstructure, essential for efficient and coordinated transmission of neural communications, undergoes pronounced development during the first years of life, while deviations to this neurodevelopmental trajectory likely result in alterations of brain connectivity relevant to behavior. Hence, systematic evaluation of white matter microstructure in the normative brain is critical for a neuroscientific approach to both typical and atypical early behavioral development. However, few studies have examined the infant brain in detail, particularly in infants under 3 months of age. Here, we utilize quantitative techniques of diffusion tensor imaging and neurite orientation dispersion and density imaging to investigate neonatal white matter microstructure in 104 infants. An optimized multiple b-value diffusion protocol was developed to allow for successful acquisition during non-sedated sleep. Associations between white matter microstructure measures and gestation corrected age, regional asymmetries, infant sex, as well as newborn growth measures were assessed. Results highlight changes of white matter microstructure during the earliest periods of development and demonstrate differential timing of developing regions and regional asymmetries. Our results contribute to a growing body of research investigating the neurobiological changes associated with neurodevelopment and suggest that characteristics of white matter microstructure are already underway in the weeks immediately following birth.

Detailing neuroanatomical development in late childhood and early adolescence using NODDI

Diffusion tensor imaging (DTI) studies have provided much evidence of white and subcortical gray matter changes during late childhood and early adolescence that suggest increasing myelination, axon density, and/or fiber coherence. Neurite orientation dispersion and density imaging (NODDI) can be used to further characterize development in white and subcortical grey matter regions in the brain by improving specificity of the MRI signal compared to conventional DTI. We used measures from NODDI and DTI to examine white and subcortical gray matter development in a group of 27 healthy participants aged 8–13 years. Neurite density index (NDI) was strongly correlated with age in nearly all regions, and was more strongly associated with age than fractional anisotropy (FA). No significant correlations were observed between orientation dispersion index (ODI) and age. This suggests that white matter and subcortical gray matter changes during late childhood and adolescence are dominated by changes in neurite density (i.e., axon density and myelination), rather than increasing coherence of axons. Within brain regions, FA was correlated with both ODI and NDI while mean diffusivity was only related to neurite density, providing further information about the structural variation across individuals. Data-driven clustering of the NODDI parameters showed that microstructural profiles varied along layers of white matter, but that that much of the white and subcortical gray matter matured in a similar manner. Clustering highlighted isolated brain regions with decreasing NDI values that were not apparent in region-of-interest analysis. Overall, these results help to more specifically understand patterns of white and gray matter development during late childhood and early adolescence.

Promise and pitfalls of g-ratio estimation with MRI

The fiber g-ratio is the ratio of the inner to the outer diameter of the myelin sheath of a myelinated axon. It has a limited dynamic range in healthy white matter, as it is optimized for speed of signal conduction, cellular energetics, and spatial constraints. In vivo imaging of the g-ratio in health and disease would greatly increase our knowledge of the nervous system and our ability to diagnose, monitor, and treat disease. MRI based g-ratio imaging was first conceived in 2011, and expanded to be feasible in full brain white matter with preliminary results in 2013. This manuscript reviews the growing g-ratio imaging literature and speculates on future applications. It details the methodology for imaging the g-ratio with MRI, and describes the known pitfalls and challenges in doing so.

Evaluating kurtosis-based diffusion MRI tissue models for white matter with fiber ball imaging

In order to quantify well-defined microstructural properties of brain tissue from diffusion MRI (dMRI) data, tissue models are typically employed that relate biological features, such as cell morphology and cell membrane permeability, to the diffusion dynamics. A variety of such models have been proposed for white matter, and their validation is a topic of active interest. In this paper, three different tissue models are tested by comparing their predictions for a specific microstructural parameter to a value measured independently with a recently proposed dMRI method known as fiber ball imaging (FBI). The three tissue models are all constructed with the diffusion and kurtosis tensors, and they are hence compatible with diffusional kurtosis imaging. Nevertheless, the models differ significantly in their details and predictions. For voxels with fractional anisotropies (FAs) exceeding 0.5, all three are reasonably consistent with FBI. However, for lower FA values, one of these, called the white matter tract integrity (WMTI) model, is found to be in much better accord with FBI than the other two, suggesting that the WMTI model has a broader range of applicability.

Robust and fast nonlinear optimization of diffusion MRI microstructure models

Advances in biophysical multi-compartment modeling for diffusion MRI (dMRI) have gained popularity because of greater specificity than DTI in relating the dMRI signal to underlying cellular microstructure. A large range of these diffusion microstructure models have been developed and each of the popular models comes with its own, often different, optimization algorithm, noise model and initialization strategy to estimate its parameter maps. Since data fit, accuracy and precision is hard to verify, this creates additional challenges to comparability and generalization of results from diffusion microstructure models. In addition, non-linear optimization is computationally expensive leading to very long run times, which can be prohibitive in large group or population studies. In this technical note we investigate the performance of several optimization algorithms and initialization strategies over a few of the most popular diffusion microstructure models, including NODDI and CHARMED. We evaluate whether a single well performing optimization approach exists that could be applied to many models and would equate both run time and fit aspects. All models, algorithms and strategies were implemented on the Graphics Processing Unit (GPU) to remove run time constraints, with which we achieve whole brain dataset fits in seconds to minutes. We then evaluated fit, accuracy, precision and run time for different models of differing complexity against three common optimization algorithms and three parameter initialization strategies. Variability of the achieved quality of fit in actual data was evaluated on ten subjects of each of two population studies with a different acquisition protocol. We find that optimization algorithms and multi-step optimization approaches have a considerable influence on performance and stability over subjects and over acquisition protocols. The gradient-free Powell conjugate-direction algorithm was found to outperform other common algorithms in terms of run time, fit, accuracy and precision. Parameter initialization approaches were found to be relevant especially for more complex models, such as those involving several fiber orientations per voxel. For these, a fitting cascade initializing or fixing parameter values in a later optimization step from simpler models in an earlier optimization step further improved run time, fit, accuracy and precision compared to a single step fit. This establishes and makes available standards by which robust fit and accuracy can be achieved in shorter run times. This is especially relevant for the use of diffusion microstructure modeling in large group or population studies and in combining microstructure parameter maps with tractography results.

•

Evaluate robustness of fit, accuracy, precision for diffusion microstructure models

•

Test three optimization algorithms and three parameter initialization strategies

•

GPU implementation removes run time constraints; whole brain fit within minutes

•

Powell conjugate-direction algorithm has superior fit, accuracy, precision

•

Initialization approaches are important for crossing fiber microstructure models

The potential of neuroimaging for identifying predictors of adolescent alcohol use initiation and misuse

BACKGROUND AND AIMS

Dysfunction in brain regions underlying impulse control, reward processing and executive function have been associated previously with adolescent alcohol misuse. However, identifying pre-existing neurobiological risk factors, as distinct from changes arising from early alcohol-use, is difficult. Here, we outline how neuroimaging data can identify the neural predictors of adolescent alcohol-use initiation and misuse by using prospective longitudinal studies to follow initially alcohol-naive individuals over time and by neuroimaging adolescents with inherited risk factors for alcohol misuse.

METHOD

A comprehensive narrative of the literature regarding neuroimaging studies published between 2010 and 2016 focusing on predictors of adolescent alcohol use initiation and misuse.

FINDINGS

Prospective, longitudinal neuroimaging studies have identified pre-existing differences between adolescents who remained alcohol-naive and those who transitioned subsequently to alcohol use. Both functional and structural grey matter differences were observed in temporal and frontal regions, including reduced brain activity in the superior frontal gyrus and temporal lobe, and thinner temporal cortices of future alcohol users. Interactions between brain function and genetic predispositions have been identified, including significant association found between the Ras protein-specific guanine nucleotide releasing factor 2 (RASGRF2) gene and reward-related striatal functioning.

CONCLUSIONS

Neuroimaging predictors of alcohol use have shown modest utility to date. Future research should use out-of-sample performance as a quantitative measure of a predictor's utility. Neuroimaging data should be combined across multiple modalities, including structural information such as volumetrics and cortical thickness, in conjunction with white-matter tractography. A number of relevant neurocognitive systems should be assayed; particularly, inhibitory control, reward processing and executive functioning. Combining a rich magnetic resonance imaging data set could permit the generation of neuroimaging risk scores, which could potentially yield targeted interventions.

Diffusion imaging of reversible and irreversible microstructural changes within the corticospinal tract in idiopathic normal pressure hydrocephalus

The symptoms of idiopathic normal pressure hydrocephalus (iNPH) can be improved by shunt surgery, but prediction of treatment outcome is not established. We investigated changes of the corticospinal tract (CST) in iNPH before and after shunt surgery by using diffusion microstructural imaging, which infers more specific tissue properties than conventional diffusion tensor imaging. Two biophysical models were used: neurite orientation dispersion and density imaging (NODDI) and white matter tract integrity (WMTI). In both methods, the orientational coherence within the CSTs was higher in patients than in controls, and some normalization occurred after the surgery in patients, indicating axon stretching and recovery. The estimated axon density was lower in patients than in controls but remained unchanged after the surgery, suggesting its potential as a marker for irreversible neuronal damage. In a Monte-Carlo simulation that represented model axons as undulating cylinders, both NODDI and WMTI separated the effects of axon density and undulation. Thus, diffusion MRI may distinguish between reversible and irreversible microstructural changes in iNPH. Our findings constitute a step towards a quantitative image biomarker that reflects pathological process and treatment outcomes of iNPH.

•

NODDI and WMTI provide markers of reversible and irreversible changes in iNPH.

•

Measures of axon orientation indicated recovery from stretching after surgery.

•

Axon density remained low after surgery, suggesting chronic neuronal damage.

•

Axon stretching in simulations differentially affected diffusion metrics.

A time-efficient acquisition protocol for multipurpose diffusion-weighted microstructural imaging at 7 Tesla

PURPOSE

Several diffusion-weighted MRI techniques have been developed and validated during the past 2 decades. While offering various neuroanatomical inferences, these techniques differ in their proposed optimal acquisition design, preventing clinicians and researchers benefiting from all potential inference methods, particularly when limited time is available. This study reports an optimal design that enables for a time-efficient diffusion-weighted MRI acquisition scheme at 7 Tesla. The primary audience of this article is the typical end user, interested in diffusion-weighted microstructural imaging at 7 Tesla.

METHODS

We tested b-values in the range of 700 to 3000 s/mm² with different number of angular diffusion-encoding samples, against a data-driven "gold standard."

RESULTS

The suggested design is a protocol with b-values of 1000 and 2500 s/mm² , with 25 and 50 samples, uniformly distributed over two shells. We also report a range of protocols in which the results of fitting microstructural models to the diffusion-weighted data had high correlation with the gold standard.

CONCLUSION

We estimated minimum acquisition requirements that enable diffusion tensor imaging, higher angular resolution diffusion-weighted imaging, neurite orientation dispersion, and density imaging and white matter tract integrity across whole brain with isotropic resolution of 1.8 mm in less than 11 min. Magn Reson Med 78:2170-2184, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Restriction Spectrum Imaging As a Potential Measure of Cortical Neurite Density in Autism

Autism postmortem studies have shown various cytoarchitectural anomalies in cortical and limbic areas including increased cell packing density, laminar disorganization, and narrowed minicolumns. However, there is little evidence on dendritic and axonal organization in ASD. Recent imaging techniques have the potential for non-invasive, in vivo studies of small-scale structure in the human brain, including gray matter. Here, Restriction Spectrum Imaging (RSI), a multi-shell diffusion-weighted imaging technique, was used to examine gray matter microstructure in 24 children with ASD (5 female) and 20 matched typically developing (TD) participants (2 female), ages 7–17 years. RSI extends the spherical deconvolution model to multiple length scales to characterize neurite density (ND) and organization. Measures were examined in 48 cortical regions of interest per hemisphere. To our knowledge, this is the first time that a multi-compartmental diffusion model has been applied to cortical gray matter in ASD. The ND measure detected robust age effects showing a significant positive relationship to age in all lobes except left temporal when groups were combined. Results were also suggestive of group differences (ASD<TD) in anterior cingulate, right superior temporal lobe and much of the parietal lobes, but these fell short of statistical significance. For MD, significant group differences (ASD>TD) in bilateral parietal regions as well as widespread age effects were detected. Our findings support the value of multi-shell diffusion imaging for assays of cortical gray matter. This approach has the potential to add to postmortem literature, examining intracortical organization, intracortical axonal content, myelination, or caliber. Robust age effects further support the validity of the ND metric for in vivo examination of gray matter microstructure in ASD and across development. While diffusion MRI does not approach the precision of histological studies, in vivo imaging measures of microstructure can complement postmortem studies, by allowing access to large sample sizes, a whole-brain field of view, longitudinal designs, and combination with behavioral and functional assays. This makes multi-shell diffusion imaging a promising technique for understanding the underlying cytoarchitecture of the disorder.

Analysis of the effects of noise, DWI sampling, and value of assumed parameters in diffusion MRI models

Purpose

This study was a systematic evaluation across different and prominent diffusion MRI models to better understand the ways in which scalar metrics are influenced by experimental factors, including experimental design (diffusion‐weighted imaging [DWI] sampling) and noise.

Methods

Four diffusion MRI models—diffusion tensor imaging (DTI), diffusion kurtosis imaging (DKI), mean apparent propagator MRI (MAP‐MRI), and neurite orientation dispersion and density imaging (NODDI)—were evaluated by comparing maps and histogram values of the scalar metrics generated using DWI datasets obtained in fixed mouse brain with different noise levels and DWI sampling complexity. Additionally, models were fit with different input parameters or constraints to examine the consequences of model fitting procedures.

Results

Experimental factors affected all models and metrics to varying degrees. Model complexity influenced sensitivity to DWI sampling and noise, especially for metrics reporting non‐Gaussian information. DKI metrics were highly susceptible to noise and experimental design. The influence of fixed parameter selection for the NODDI model was found to be considerable, as was the impact of initial tensor fitting in the MAP‐MRI model.

Conclusion

Across DTI, DKI, MAP‐MRI, and NODDI, a wide range of dependence on experimental factors was observed that elucidate principles and practical implications for advanced diffusion MRI. Magn Reson Med 78:1767–1780, 2017. © 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‐NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

MRI and M/EEG studies of the White Matter Development in Human Fetuses and Infants: Review and Opinion

Already during the last trimester of gestation, functional responses are recorded in foetuses and preterm newborns, attesting an already complex cerebral architecture. Then throughout childhood, anatomical connections are further refined but at different rates and over asynchronous periods across functional networks. Concurrently, infants gradually achieve new psychomotor and cognitive skills. Only the recent use of non-invasive techniques such as magnetic resonance imaging (MRI) and magneto- and electroencephalography (M/EEG) has opened the possibility to understand the relationships between brain maturation and skills development in vivo. In this review, we describe how these techniques have been applied to study the white matter maturation. At the structural level, the early architecture and myelination of bundles have been assessed with diffusion and relaxometry MRI, recently integrated in multi-compartment models and multi-parametric approaches. Nevertheless, technical limitations prevent us to map major developmental mechanisms such as fibers growth and pruning, and the progressive maturation at the bundle scale in case of mixing trajectories. At the functional level, M/EEG have been used to record different visual, somatosensory and auditory evoked responses. Because the conduction velocity of neural impulses increases with the myelination of connections, major changes in the components latency are observed throughout development. But so far, only a few studies have related structural and functional markers of white matter myelination. Such multi-modal approaches will be a major challenge in future research, not only to understand normal development, but also to characterize early mechanisms of pathologies and the influence of fetal and perinatal interventions on later outcome.

Microstructure Imaging of Crossing (MIX) White Matter Fibers from diffusion MRI

Diffusion MRI (dMRI) reveals microstructural features of the brain white matter by quantifying the anisotropic diffusion of water molecules within axonal bundles. Yet, identifying features such as axonal orientation dispersion, density, diameter, etc., in complex white matter fiber configurations (e.g. crossings) has proved challenging. Besides optimized data acquisition and advanced biophysical models, computational procedures to fit such models to the data are critical. However, these procedures have been largely overlooked by the dMRI microstructure community and new, more versatile, approaches are needed to solve complex biophysical model fitting problems. Existing methods are limited to models assuming single fiber orientation, relevant to limited brain areas like the corpus callosum, or multiple orientations but without the ability to extract detailed microstructural features. Here, we introduce a new and versatile optimization technique (MIX), which enables microstructure imaging of crossing white matter fibers. We provide a MATLAB implementation of MIX, and demonstrate its applicability to general microstructure models in fiber crossings using synthetic as well as ex-vivo and in-vivo brain data.

Assessing White Matter Microstructure in Brain Regions with Different Myelin Architecture Using MRI

OBJECTIVE

We investigate how known differences in myelin architecture between regions along the cortico-spinal tract and frontal white matter (WM) in 19 healthy adolescents are reflected in several quantitative MRI parameters that have been proposed to non-invasively probe WM microstructure. In a clinically feasible scan time, both conventional imaging sequences as well as microstructural MRI parameters were assessed in order to quantitatively characterise WM regions that are known to differ in the thickness of their myelin sheaths, and in the presence of crossing or parallel fibre organisation.

RESULTS

We found that diffusion imaging, MR spectroscopy (MRS), myelin water fraction (MWF), Magnetization Transfer Imaging, and Quantitative Susceptibility Mapping were myelin-sensitive in different ways, giving complementary information for characterising WM microstructure with different underlying fibre architecture. From the diffusion parameters, neurite density (NODDI) was found to be more sensitive than fractional anisotropy (FA), underlining the limitation of FA in WM crossing fibre regions. In terms of sensitivity to different myelin content, we found that MWF, the mean diffusivity and chemical-shift imaging based MRS yielded the best discrimination between areas.

CONCLUSION

Multimodal assessment of WM microstructure was possible within clinically feasible scan times using a broad combination of quantitative microstructural MRI sequences. By assessing new microstructural WM parameters we were able to provide normative data and discuss their interpretation in regions with different myelin architecture, as well as their possible application as biomarker for WM disorders.

Characteristics of Diffusional Kurtosis in Chronic Ischemia of Adult Moyamoya Disease: Comparing Diffusional Kurtosis and Diffusion Tensor Imaging

BACKGROUND AND PURPOSE

Detecting microstructural changes due to chronic ischemia potentially enables early identification of patients at risk of cognitive impairment. In this study, diffusional kurtosis imaging and diffusion tensor imaging were used to investigate whether the former provides additional information regarding microstructural changes in the gray and white matter of adult patients with Moyamoya disease.

MATERIALS AND METHODS

MR imaging (diffusional kurtosis imaging and DTI) was performed in 23 adult patients with Moyamoya disease and 23 age-matched controls. Three parameters were extracted from diffusional kurtosis imaging (mean kurtosis, axial kurtosis, and radial kurtosis), and 4, from DTI (fractional anisotropy, radial diffusivity, mean diffusivity, and axial diffusivity). Voxelwise analysis for these parameters was performed in the normal-appearing brain parenchyma. The association of these parameters with neuropsychological performance was also evaluated.

RESULTS

Voxelwise analysis revealed the greatest differences in fractional anisotropy, followed, in order, by radial diffusivity, mean diffusivity, and mean kurtosis. In patients, diffusional kurtosis imaging parameters were decreased in the dorsal deep white matter such as the corona radiata and superior longitudinal fasciculus (P < .01), including areas without DTI abnormality. Superior longitudinal fasciculus fiber-crossing areas showed weak correlations between diffusional kurtosis imaging and DTI parameters compared with tissues with a single-fiber direction (eg, the corpus callosum). Diffusional kurtosis imaging parameters were associated with general intelligence and frontal lobe performance.

CONCLUSIONS

Although DTI revealed extensive white matter changes, diffusional kurtosis imaging additionally demonstrated microstructural changes in ischemia-prone deep white matter with abundant fiber crossings. Thus, diffusional kurtosis imaging may be a useful adjunct for detecting subtle chronic ischemic injuries.

Anisotropic Sampling Shape of White Matter Microstructure Cannot Cheat Diffusional Kurtosis

Diffusion kurtosis imaging is a newly developed diffusion magnetic resonance imaging technique, which is becoming increasingly valuable in clinical practice. Although low-resolution sampling is commonly used to compensate the unsteadiness of kurtosis estimation, the influence of the sampling shape has not been investigated. In this study, by using two different acquisition protocols, isotropic and anisotropic sampling voxels were acquired and their influence on various white matter structures was observed. Fiber tracking, T-tests, and correlation analysis were used to quantify the difference between the anisotropic and isotropic sampling. A significant difference (p<0.01) was found in the fractional anisotropic level but not in kurtosis. The results presented here can provide a basis for higher resolution as well as higher quality kurtosis mapping, which may be of great significance in clinical examinations.

Age-dependent differences in brain tissue microstructure assessed with neurite orientation dispersion and density imaging

Human aging is accompanied by progressive changes in executive function and memory, but the biological mechanisms underlying these phenomena are not fully understood. Using neurite orientation dispersion and density imaging, we sought to examine the relationship between age, cellular microstructure, and neuropsychological scores in 116 late middle-aged, cognitively asymptomatic participants. Results revealed widespread increases in the volume fraction of isotropic diffusion and localized decreases in neurite density in frontal white matter regions with increasing age. In addition, several of these microstructural alterations were associated with poorer performance on tests of memory and executive function. These results suggest that neurite orientation dispersion and density imaging is capable of measuring age-related brain changes and the neural correlates of poorer performance on tests of cognitive functioning, largely in accordance with published histological findings and brain-imaging studies of people of this age range. Ultimately, this study sheds light on the processes underlying normal brain development in adulthood, knowledge that is critical for differentiating healthy aging from changes associated with dementia.

Alterations of Diffusion Kurtosis and Neurite Density Measures in Deep Grey Matter and White Matter in Parkinson's Disease

In Parkinson’s disease (PD), pathological microstructural changes occur and such changes might be detected using diffusion magnetic resonance imaging (dMRI). However, it is unclear whether dMRI improves PD diagnosis or helps differentiating between phenotypes, such as postural instability gait difficulty (PIGD) and tremor dominant (TD) PD. We included 105 patients with PD and 44 healthy controls (HC), all of whom underwent dMRI as part of the prospective Swedish BioFINDER study. Diffusion kurtosis imaging (DKI) and neurite density imaging (NDI) analyses were performed using regions of interest in the basal ganglia, the thalamus, the pons and the midbrain as well as tractography of selected white matter tracts. In the putamen, the PD group showed increased mean diffusivity (MD) (p = .003), decreased fractional anisotropy (FA) (p = .001) and decreased mean kurtosis (MK), compared to HC (p = .024). High MD and a low MK in the putamen were associated with more severe motor and cognitive symptomatology (p < .05). Also, patients with PIGD exhibited increased MD in the putamen compared to the TD patients (p = .009). In the thalamus, MD was increased (p = .001) and FA was decreased (p = .032) in PD compared to HC. Increased MD and decreased FA correlated negatively with motor speed and balance (p < .05). In the superior longitudinal fasciculus (SLF), MD (p = .019) and f_iso were increased in PD compared to HC (p = .03). These changes correlated negatively with motor speed (p < .002) and balance (p < .037). However, most of the observed changes in PD were also present in cases with either multiple system atrophy (n = 11) or progressive supranuclear palsy (n = 10). In conclusion, PD patients exhibit microstructural changes in the putamen, the thalamus, and the SLF, which are associated with worse disease severity. However, the dMRI changes are not sufficiently specific to improve the diagnostic work-up of PD. Longitudinal studies should evaluate whether dMRI measures can be used to track disease progression.

Microstructural maturation of white matter tracts in encephalopathic neonates

PURPOSE

This study aims to apply neurite orientation dispersion and density imaging (NODDI) to measure white matter microstructural features during early development.

METHODS

NODDI parameters were measured in twelve newborns and thirteen 6-month infants, all with perinatal clinical encephalopathy.

RESULTS

Between 0 and 6 months, there were significant differences in fractional anisotropy (FA) for all tracts; in neurite density for internal capsules, optic radiations, and splenium; and in orientation dispersion for anterior limb of internal capsule and optic radiations. There were no appreciable differences in NODDI parameters related to outcome.

CONCLUSION

NODDI may allow more detailed characterization of microstructural maturation than FA.

Abnormalities in Diffusional Kurtosis Metrics Related to Head Impact Exposure in a Season of High School Varsity Football

The purpose of this study was to determine whether the effects of cumulative head impacts during a season of high school football produce changes in diffusional kurtosis imaging (DKI) metrics in the absence of clinically diagnosed concussion. Subjects were recruited from a high school football team and were outfitted with the Head Impact Telemetry System (HITS) during all practices and games. Biomechanical head impact exposure metrics were calculated, including: total impacts, summed acceleration, and Risk Weighted Cumulative Exposure (RWE). Twenty-four players completed pre- and post-season magnetic resonance imaging, including DKI; players who experienced clinical concussion were excluded. Fourteen subjects completed pre- and post-season Immediate Post-Concussion Assessment and Cognitive Testing (ImPACT). DKI-derived metrics included mean kurtosis (MK), axial kurtosis (K axial), and radial kurtosis (K radial), and white matter modeling (WMM) parameters included axonal water fraction, tortuosity of the extra-axonal space, extra-axonal diffusivity (D_e axial and radial), and intra-axonal diffusivity (D_a). These metrics were used to determine the total number of abnormal voxels, defined as 2 standard deviations above or below the group mean. Linear regression analysis revealed a statistically significant relationship between RWE combined probability (RWE_CP) and MK. Secondary analysis of other DKI-derived and WMM metrics demonstrated statistically significant linear relationships with RWE_CP after covariate adjustment. These results were compared with the results of DTI-derived metrics from the same imaging sessions in this exact same cohort. Several of the DKI-derived scalars (D_a, MK, K axial, and K radial) explained more variance, compared with RWE_CP, suggesting that DKI may be more sensitive to subconcussive head impacts. No significant relationships between DKI-derived metrics and ImPACT measures were found. It is important to note that the pathological implications of these metrics are not well understood. In summary, we demonstrate a single season of high school football can produce DKI measurable changes in the absence of clinically diagnosed concussion.

Quantification of normal-appearing white matter tract integrity in multiple sclerosis: a diffusion kurtosis imaging study

Our aim was to characterize the nature and extent of pathological changes in the normal-appearing white matter (NAWM) of patients with multiple sclerosis (MS) using novel diffusion kurtosis imaging-derived white matter tract integrity (WMTI) metrics and to investigate the association between these WMTI metrics and clinical parameters. Thirty-two patients with relapsing-remitting MS and 19 age- and gender-matched healthy controls underwent MRI and neurological examination. Maps of mean diffusivity, fractional anisotropy and WMTI metrics (intra-axonal diffusivity, axonal water fraction, tortuosity and axial and radial extra-axonal diffusivity) were created. Tract-based spatial statistics analysis was performed to assess for differences in the NAWM between patients and controls. A region of interest analysis of the corpus callosum was also performed to assess for group differences and to evaluate correlations between WMTI metrics and measures of disease severity. Mean diffusivity and radial extra-axonal diffusivity were significantly increased while fractional anisotropy, axonal water fraction, intra-axonal diffusivity and tortuosity were decreased in MS patients compared with controls (p values ranging from <0.001 to <0.05). Axonal water fraction in the corpus callosum was significantly associated with the expanded disability status scale score (ρ = -0.39, p = 0.035). With the exception of the axial extra-axonal diffusivity, all metrics were correlated with the symbol digits modality test score (p values ranging from 0.001 to <0.05). WMTI metrics are thus sensitive to changes in the NAWM of MS patients and might provide a more pathologically specific, clinically meaningful and practical complement to standard diffusion tensor imaging-derived metrics.

In vivo quantification of demyelination and recovery using compartment-specific diffusion MRI metrics validated by electron microscopy

There is a need for accurate quantitative non-invasive biomarkers to monitor myelin pathology in vivo and distinguish myelin changes from other pathological features including inflammation and axonal loss. Conventional MRI metrics such as T2, magnetization transfer ratio and radial diffusivity have proven sensitivity but not specificity. In highly coherent white matter bundles, compartment-specific white matter tract integrity (WMTI) metrics can be directly derived from the diffusion and kurtosis tensors: axonal water fraction, intra-axonal diffusivity, and extra-axonal radial and axial diffusivities. We evaluate the potential of WMTI to quantify demyelination by monitoring the effects of both acute (6weeks) and chronic (12weeks) cuprizone intoxication and subsequent recovery in the mouse corpus callosum, and compare its performance with that of conventional metrics (T2, magnetization transfer, and DTI parameters). The changes observed in vivo correlated with those obtained from quantitative electron microscopy image analysis. A 6-week intoxication produced a significant decrease in axonal water fraction (p<0.001), with only mild changes in extra-axonal radial diffusivity, consistent with patchy demyelination, while a 12-week intoxication caused a more marked decrease in extra-axonal radial diffusivity (p=0.0135), consistent with more severe demyelination and clearance of the extra-axonal space. Results thus revealed increased specificity of the axonal water fraction and extra-axonal radial diffusivity parameters to different degrees and patterns of demyelination. The specificities of these parameters were corroborated by their respective correlations with microstructural features: the axonal water fraction correlated significantly with the electron microscopy derived total axonal water fraction (ρ=0.66; p=0.0014) but not with the g-ratio, while the extra-axonal radial diffusivity correlated with the g-ratio (ρ=0.48; p=0.0342) but not with the electron microscopy derived axonal water fraction. These parameters represent promising candidates as clinically feasible biomarkers of demyelination and remyelination in the white matter.

Heritability of complex white matter diffusion traits assessed in a population isolate

INTRODUCTION

Diffusion weighted imaging (DWI) methods can noninvasively ascertain cerebral microstructure by examining pattern and directions of water diffusion in the brain. We calculated heritability for DWI parameters in cerebral white (WM) and gray matter (GM) to study the genetic contribution to the diffusion signals across tissue boundaries.

METHODS

Using Old Order Amish (OOA) population isolate with large family pedigrees and high environmental homogeneity, we compared the heritability of measures derived from three representative DWI methods targeting the corpus callosum WM and cingulate gyrus GM: diffusion tensor imaging (DTI), the permeability-diffusivity (PD) model, and the neurite orientation dispersion and density imaging (NODDI) model. These successively more complex models represent the diffusion signal modeling using one, two, and three diffusion compartments, respectively.

RESULTS

We replicated the high heritability of the DTI-based fractional anisotropy (h(2)  = 0.67) and radial diffusivity (h(2)  = 0.72) in WM. High heritability in both WM and GM tissues were observed for the permeability-diffusivity index from the PD model (h(2)  = 0.64 and 0.84), and the neurite density from the NODDI model (h(2)  = 0.70 and 0.55). The orientation dispersion index from the NODDI model was only significantly heritable in GM (h(2)  = 0.68).

CONCLUSION

DWI measures from multicompartmental models were significantly heritable in WM and GM. DWI can offer valuable phenotypes for genetic research; and genes thus identified may reveal mechanisms contributing to mental and neurological disorders in which diffusion imaging anomalies are consistently found. Hum Brain Mapp 37:525-535, 2016. © 2015 Wiley Periodicals, Inc.

Bingham-NODDI: Mapping anisotropic orientation dispersion of neurites using diffusion MRI

This paper presents Bingham-NODDI, a clinically-feasible technique for estimating the anisotropic orientation dispersion of neurites. Direct quantification of neurite morphology on clinical scanners was recently realised by a diffusion MRI technique known as neurite orientation dispersion and density imaging (NODDI). However in its current form NODDI cannot estimate anisotropic orientation dispersion, which is widespread in the brain due to common fanning and bending of neurites. This work proposes Bingham-NODDI that extends the NODDI formalism to address this limitation. Bingham-NODDI characterises anisotropic orientation dispersion by utilising the Bingham distribution to model neurite orientation distribution. The new model estimates the extent of dispersion about the dominant orientation, separately along the primary and secondary dispersion orientations. These estimates are subsequently used to estimate the overall dispersion about the dominant orientation and the dispersion anisotropy. We systematically evaluate the ability of the new model to recover these key parameters of anisotropic orientation dispersion with standard NODDI protocol, both in silico and in vivo. The results demonstrate that the parameters of the proposed model can be estimated without additional acquisition requirements over the standard NODDI protocol. Thus anisotropic dispersion can be determined and has the potential to be used as a marker for normal brain development and ageing or in pathology. We additionally find that the original NODDI model is robust to the effects of anisotropic orientation dispersion, when the quantification of anisotropic dispersion is not of interest.

Including diffusion time dependence in the extra-axonal space improves in vivo estimates of axonal diameter and density in human white matter

Axonal density and diameter are two fundamental properties of brain white matter. Recently, advanced diffusion MRI techniques have made these two parameters accessible in vivo. However, the techniques available to estimate such parameters are still under development. For example, current methods to map axonal diameters capture relative trends over different structures, but consistently over-estimate absolute diameters. Axonal density estimates are more accessible experimentally, but different modeling approaches exist and the impact of the experimental parameters has not been thoroughly quantified, potentially leading to incompatibility of results obtained in different studies using different techniques. Here, we characterise the impact of diffusion time on axonal density and diameter estimates using Monte Carlo simulations and STEAM diffusion MRI at 7 T on 9 healthy volunteers. We show that axonal density and diameter estimates strongly depend on diffusion time, with diameters almost invariably overestimated and density both over and underestimated for some commonly used models. Crucially, we also demonstrate that these biases are reduced when the model accounts for diffusion time dependency in the extra-axonal space. For axonal density estimates, both upward and downward bias in different situations are removed by modeling extra-axonal time-dependence, showing increased accuracy in these estimates. For axonal diameter estimates, we report increased accuracy in ground truth simulations and axonal diameter estimates decreased away from high values given by earlier models and towards known values in the human corpus callosum when modeling extra-axonal time-dependence. Axonal diameter feasibility under both advanced and clinical settings is discussed in the light of the proposed advances.

Degeneracy in model parameter estimation for multi-compartmental diffusion in neuronal tissue

The ultimate promise of diffusion MRI (dMRI) models is specificity to neuronal microstructure, which may lead to distinct clinical biomarkers using noninvasive imaging. While multi-compartment models are a common approach to interpret water diffusion in the brain in vivo, the estimation of their parameters from the dMRI signal remains an unresolved problem. Practically, even when q space is highly oversampled, nonlinear fit outputs suffer from heavy bias and poor precision. So far, this has been alleviated by fixing some of the model parameters to a priori values, for improved precision at the expense of accuracy. Here we use a representative two-compartment model to show that fitting fails to determine the five model parameters from over 60 measurement points. For the first time, we identify the reasons for this poor performance. The first reason is the existence of two local minima in the parameter space for the objective function of the fitting procedure. These minima correspond to qualitatively different sets of parameters, yet they both lie within biophysically plausible ranges. We show that, at realistic signal-to-noise ratio values, choosing between the two minima based on the associated objective function values is essentially impossible. Second, there is an ensemble of very low objective function values around each of these minima in the form of a pipe. The existence of such a direction in parameter space, along which the objective function profile is very flat, explains the bias and large uncertainty in parameter estimation, and the spurious parameter correlations: in the presence of noise, the minimum can be randomly displaced by a very large amount along each pipe. Our results suggest that the biophysical interpretation of dMRI model parameters crucially depends on establishing which of the minima is closer to the biophysical reality and the size of the uncertainty associated with each parameter.

Whole-Brain In-vivo Measurements of the Axonal G-Ratio in a Group of 37 Healthy Volunteers

The g-ratio, quantifying the ratio between the inner and outer diameters of a fiber, is an important microstructural characteristic of fiber pathways and is functionally related to conduction velocity. We introduce a novel method for estimating the MR g-ratio non-invasively across the whole brain using high-fidelity magnetization transfer (MT) imaging and single-shell diffusion MRI. These methods enabled us to map the MR g-ratio in vivo across the brain's prominent fiber pathways in a group of 37 healthy volunteers and to estimate the inter-subject variability. Effective correction of susceptibility-related distortion artifacts was essential before combining the MT and diffusion data, in order to reduce partial volume and edge artifacts. The MR g-ratio is in good qualitative agreement with histological findings despite the different resolution and spatial coverage of MRI and histology. The MR g-ratio holds promise as an important non-invasive biomarker due to its microstructural and functional relevance in neurodegeneration.

Nonconventional MRI and microstructural cerebral changes in multiple sclerosis

MRI has become the most important paraclinical tool for diagnosing and monitoring patients with multiple sclerosis (MS). However, conventional MRI sequences are largely nonspecific in the pathology they reveal, and only provide a limited view of the complex morphological changes associated with MS. Nonconventional MRI techniques, such as magnetization transfer imaging (MTI), diffusion-weighted imaging (DWI) and susceptibility-weighted imaging (SWI) promise to complement existing techniques by revealing more-specific information on microstructural tissue changes. Past years have witnessed dramatic advances in the acquisition and analysis of such imaging data, and numerous studies have used these tools to probe tissue alterations associated with MS. Other MRI-based techniques-such as myelin-water imaging, (23)Na imaging, magnetic resonance elastography and magnetic resonance perfusion imaging-might also shed new light on disease-associated changes. This Review summarizes the rapid technical progress in the use of MRI in patients with MS, with a focus on nonconventional structural MRI. We critically discuss the present utility of nonconventional MRI in MS, and provide an outlook on future applications, including clinical practice. This information should allow appropriate selection of advanced MRI techniques, and facilitate their use in future studies of this disease.

Detection of acute nervous system injury with advanced diffusion-weighted MRI: a simulation and sensitivity analysis

Diffusion-weighted imaging (DWI) is a powerful tool to investigate the microscopic structure of the central nervous system (CNS). Diffusion tensor imaging (DTI), a common model of the DWI signal, has a demonstrated sensitivity to detect microscopic changes as a result of injury or disease. However, DTI and other similar models have inherent limitations that reduce their specificity for certain pathological features, particularly in tissues with complex fiber arrangements. Methods such as double pulsed field gradient (dPFG) and q-vector magic angle spinning (qMAS) have been proposed to specifically probe the underlying microscopic anisotropy without interference from the macroscopic tissue organization. This is particularly important for the study of acute injury, where abrupt changes in the microscopic morphology of axons and dendrites manifest as focal enlargements known as beading. The purpose of this work was to assess the relative sensitivity of DWI measures to beading in the context of macroscopic fiber organization and edema. Computational simulations of DWI experiments in normal and beaded axons demonstrated that, although DWI models can be highly specific for the simulated pathologies of beading and volume fraction changes in coherent fiber pathways, their sensitivity to a single idealized pathology is considerably reduced in crossing and dispersed fibers. However, dPFG and qMAS have a high sensitivity for beading, even in complex fiber tracts. Moreover, in tissues with coherent arrangements, such as the spinal cord or nerve fibers in which tract orientation is known a priori, a specific dPFG sequence variant decreases the effects of edema and improves specificity for beading. Collectively, the simulation results demonstrate that advanced DWI methods, particularly those which sample diffusion along multiple directions within a single acquisition, have improved sensitivity to acute axonal injury over conventional DTI metrics and hold promise for more informative clinical diagnostic use in CNS injury evaluation.

Diffusion kurtosis imaging probes cortical alterations and white matter pathology following cuprizone induced demyelination and spontaneous remyelination

Although MRI is the gold standard for the diagnosis and monitoring of multiple sclerosis (MS), current conventional MRI techniques often fail to detect cortical alterations and provide little information about gliosis, axonal damage and myelin status of lesioned areas. Diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) provide sensitive and complementary measures of the neural tissue microstructure. Additionally, specific white matter tract integrity (WMTI) metrics modelling the diffusion in white matter were recently derived. In the current study we used the well-characterized cuprizone mouse model of central nervous system demyelination to assess the temporal evolution of diffusion tensor (DT), diffusion kurtosis tensor (DK) and WMTI-derived metrics following acute inflammatory demyelination and spontaneous remyelination. While DT-derived metrics were unable to detect cuprizone induced cortical alterations, the mean kurtosis (MK) and radial kurtosis (RK) were found decreased under cuprizone administration, as compared to age-matched controls, in both the motor and somatosensory cortices. The MK remained decreased in the motor cortices at the end of the recovery period, reflecting long lasting impairment of myelination. In white matter, DT, DK and WMTI-derived metrics enabled the detection of cuprizone induced changes differentially according to the stage and the severity of the lesion. More specifically, the MK, the RK and the axonal water fraction (AWF) were the most sensitive for the detection of cuprizone induced changes in the genu of the corpus callosum, a region less affected by cuprizone administration. Additionally, microgliosis was associated with an increase of MK and RK during the acute inflammatory demyelination phase. In regions undergoing severe demyelination, namely the body and splenium of the corpus callosum, DT-derived metrics, notably the mean diffusion (MD) and radial diffusion (RD), were among the best discriminators between cuprizone and control groups, hence highlighting their ability to detect both acute and long lasting changes. Interestingly, WMTI-derived metrics showed the aptitude to distinguish between the different stages of the disease. Both the intra-axonal diffusivity (Da) and the AWF were found to be decreased in the cuprizone treated group, Da specifically decreased during the acute inflammatory demyelinating phase whereas the AWF decrease was associated to the spontaneous remyelination and the recovery period. Altogether our results demonstrate that DKI is sensitive to alterations of cortical areas and provides, along with WMTI metrics, information that is complementary to DT-derived metrics for the characterization of demyelination in both white and grey matter and subsequent inflammatory processes associated with a demyelinating event.

N-acetyl-aspartate levels correlate with intra-axonal compartment parameters from diffusion MRI

Diffusion MRI combined with biophysical modeling allows for the description of a white matter (WM) fiber bundle in terms of compartment specific white matter tract integrity (WMTI) metrics, which include intra-axonal diffusivity (Daxon), extra-axonal axial diffusivity (De||), extra-axonal radial diffusivity (De┴), axonal water fraction (AWF), and tortuosity (α) of extra-axonal space. Here we derive these parameters from diffusion kurtosis imaging to examine their relationship to concentrations of global WM N-acetyl-aspartate (NAA), creatine (Cr), choline (Cho) and myo-Inositol (mI), as measured with proton MR spectroscopy ((1)H-MRS), in a cohort of 25 patients with mild traumatic brain injury (MTBI). We found statistically significant (p<0.05) positive correlations between NAA and Daxon, AWF, α, and fractional anisotropy; negative correlations between NAA and De,┴ and the overall radial diffusivity (D┴). These correlations were supported by similar findings in regional analysis of the genu and splenium of the corpus callosum. Furthermore, a positive correlation in global WM was noted between Daxon and Cr, as well as a positive correlation between De|| and Cho, and a positive trend between De|| and mI. The specific correlations between NAA, an endogenous probe of the neuronal intracellular space, and WMTI metrics related to the intra-axonal space, combined with the specific correlations of De|| with mI and Cho, both predominantly present extra-axonally, corroborate the overarching assumption of many advanced modeling approaches that diffusion imaging can disentangle between the intra- and extra-axonal compartments in WM fiber bundles. Our findings are also generally consistent with what is known about the pathophysiology of MTBI, which appears to involve both intra-axonal injury (as reflected by a positive trend between NAA and Daxon) as well as axonal shrinkage, demyelination, degeneration, and/or loss (as reflected by correlations between NAA and De┴, AWF, and α).