Imaging brain microstructure with diffusion MRI: practicality and applications

Secondary Degeneration of White Matter After Focal Sensorimotor Cortical Ischemic Stroke in Rats

Ischemic lesions could lead to secondary degeneration in remote regions of the brain. However, the spatial distribution of secondary degeneration along with its role in functional deficits is not well understood. In this study, we explored the spatial and connectivity properties of white matter (WM) secondary degeneration in a focal unilateral sensorimotor cortical ischemia rat model, using advanced microstructure imaging on a 14 T MRI system. Significant axonal degeneration was observed in the ipsilateral external capsule and even remote regions including the contralesional external capsule and corpus callosum. Further fiber tractography analysis revealed that only fibers having direct axonal connections with the primary lesion exhibited a significant degeneration. These results suggest that focal ischemic lesions may induce remote WM degeneration, but limited to fibers tied to the primary lesion. These “direct” fibers mainly represent perilesional, interhemispheric, and subcortical axonal connections. At last, we found that primary lesion volume might be the determining factor of motor function deficits.

Longitudinal Reproducibility of Neurite Orientation Dispersion and Density Imaging (NODDI) Derived Metrics in the White Matter

Diffusion-weighted magnetic resonance imaging (DWI) is undergoing constant evolution with the ambitious goal of developing in-vivo histology of the brain. A recent methodological advancement is Neurite Orientation Dispersion and Density Imaging (NODDI), a histologically validated multi-compartment model to yield microstructural features of brain tissue such as geometric complexity and neurite packing density, which are especially useful in imaging the white matter. Since NODDI is increasingly popular in clinical research and fields such as developmental neuroscience and neuroplasticity, it is of vast importance to characterize its reproducibility (or reliability). We acquired multi-shell DWI data in 29 healthy young subjects twice over a rescan interval of 4 weeks to assess the within-subject coefficient of variation (CV_WS), between-subject coefficient of variation (CV_BS) and the intraclass correlation coefficient (ICC), respectively. Using these metrics, we compared regional and voxel-by-voxel reproducibility of the most common image analysis approaches (tract-based spatial statistics [TBSS], voxel-based analysis with different extents of smoothing ["VBM-style"], ROI-based analysis). We observed high test-retest reproducibility for the orientation dispersion index (ODI) and slightly worse results for the neurite density index (NDI). Our findings also suggest that the choice of analysis approach might have significant consequences for the results of a study. Collectively, the voxel-based approach with Gaussian smoothing kernels of ≥4 mm FWHM and ROI-averaging yielded the highest reproducibility across NDI and ODI maps (CV_WS mostly ≤3%, ICC mostly ≥0.8), respectively, whilst smaller kernels and TBSS performed consistently worse. Furthermore, we demonstrate that image quality (signal-to-noise ratio [SNR]) is an important determinant of NODDI metric reproducibility. We discuss the implications of these results for longitudinal and cross-sectional research designs commonly employed in the neuroimaging field.

Towards in vivo g-ratio mapping using MRI: Unifying myelin and diffusion imaging

BACKGROUND

The g-ratio, quantifying the comparative thickness of the myelin sheath encasing an axon, is a geometrical invariant that has high functional relevance because of its importance in determining neuronal conduction velocity. Advances in MRI data acquisition and signal modelling have put in vivo mapping of the g-ratio, across the entire white matter, within our reach. This capacity would greatly increase our knowledge of the nervous system: how it functions, and how it is impacted by disease.

NEW METHOD

This is the second review on the topic of g-ratio mapping using MRI.

RESULTS

This review summarizes the most recent developments in the field, while also providing methodological background pertinent to aggregate g-ratio weighted mapping, and discussing pitfalls associated with these approaches.

COMPARISON WITH EXISTING METHODS

Using simulations based on recently published data, this review reveals caveats to the state-of-the-art calibration methods that have been used for in vivo g-ratio mapping. It highlights the need to estimate both the slope and offset of the relationship between these MRI-based markers and the true myelin volume fraction if we are really to achieve the goal of precise, high sensitivity g-ratio mapping in vivo. Other challenges discussed in this review further evidence the need for gold standard measurements of human brain tissue from ex vivo histology.

CONCLUSIONS

We conclude that the quest to find the most appropriate MRI biomarkers to enable in vivo g-ratio mapping is ongoing, with the full potential of many novel techniques yet to be investigated.

The sensitivity of diffusion MRI to microstructural properties and experimental factors

Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.

Axon morphology is modulated by the local environment and impacts the noninvasive investigation of its structure-function relationship

Axonal conduction velocity, which ensures efficient function of the brain network, is related to axon diameter. Noninvasive, in vivo axon diameter estimates can be made with diffusion magnetic resonance imaging, but the technique requires three-dimensional (3D) validation. Here, high-resolution, 3D synchrotron X-ray nano-holotomography images of white matter samples from the corpus callosum of a monkey brain reveal that blood vessels, cells, and vacuoles affect axonal diameter and trajectory. Within single axons, we find that the variation in diameter and conduction velocity correlates with the mean diameter, contesting the value of precise diameter determination in larger axons. These complex 3D axon morphologies drive previously reported 2D trends in axon diameter and g-ratio. Furthermore, we find that these morphologies bias the estimates of axon diameter with diffusion magnetic resonance imaging and, ultimately, impact the investigation and formulation of the axon structure-function relationship.

Along-tract statistics of neurite orientation dispersion and density imaging diffusion metrics to enhance MR tractography quantitative analysis in healthy controls and in patients with brain tumors

Along‐tract statistics analysis enables the extraction of quantitative diffusion metrics along specific white matter fiber tracts. Besides quantitative metrics derived from classical diffusion tensor imaging (DTI), such as fractional anisotropy and diffusivities, new parameters reflecting the relative contribution of different diffusion compartments in the tissue can be estimated through advanced diffusion MRI methods as neurite orientation dispersion and density imaging (NODDI), leading to a more specific microstructural characterization. In this study, we extracted both DTI‐ and NODDI‐derived quantitative microstructural diffusion metrics along the most eloquent fiber tracts in 15 healthy subjects and in 22 patients with brain tumors. We obtained a robust intraprotocol reference database of normative along‐tract microstructural metrics, and their corresponding plots, from healthy fiber tracts. Each diffusion metric of individual patient's fiber tract was then plotted and statistically compared to the normative profile of the corresponding metric from the healthy fiber tracts. NODDI‐derived metrics appeared to account for the pathological microstructural changes of the peritumoral tissue more accurately than DTI‐derived ones. This approach may be useful for future studies that may compare healthy subjects to patients diagnosed with other pathological conditions.

Realistic Microstructure Simulator (RMS): Monte Carlo simulations of diffusion in three-dimensional cell segmentations of microscopy images

BACKGROUND

Monte Carlo simulations of diffusion are commonly used as a model validation tool as they are especially suitable for generating the diffusion MRI signal in complicated tissue microgeometries.

NEW METHOD

Here we describe the details of implementing Monte Carlo simulations in three-dimensional (3d) voxelized segmentations of cells in microscopy images. Using the concept of the corner reflector, we largely reduce the computational load of simulating diffusion within and exchange between multiple cells. Precision is further achieved by GPU-based parallel computations.

RESULTS

Our simulation of diffusion in white matter axons segmented from a mouse brain demonstrates its value in validating biophysical models. Furthermore, we provide the theoretical background for implementing a discretized diffusion process, and consider the finite-step effects of the particle-membrane reflection and permeation events, needed for efficient simulation of interactions with irregular boundaries, spatially variable diffusion coefficient, and exchange.

COMPARISON WITH EXISTING METHODS

To our knowledge, this is the first Monte Carlo pipeline for MR signal simulations in a substrate composed of numerous realistic cells, accounting for their permeable and irregularly-shaped membranes.

CONCLUSIONS

The proposed RMS pipeline makes it possible to achieve fast and accurate simulations of diffusion in realistic tissue microgeometry, as well as the interplay with other MR contrasts. Presently, RMS focuses on simulations of diffusion, exchange, and T₁ and T₂ NMR relaxation in static tissues, with a possibility to straightforwardly account for susceptibility-induced T₂^* effects and flow.

CRAFT (Cerclage after full dilatation caesarean section): protocol of a mixed methods study investigating the role of previous in-labour caesarean section in preterm birth risk

BACKGROUND

Full dilatation caesarean sections are associated with recurrent early spontaneous preterm birth and late miscarriage. The risk following first stage caesarean sections, are less well defined, but appears to be increased in late-first stage of labour. The mechanism for this increased risk of late miscarriage and early spontaneous preterm birth in these women is unknown and there are uncertainties with regards to clinical management. Current predictive models of preterm birth (based on transvaginal ultrasound and quantitative fetal fibronectin) have not been validated in these women and it is unknown whether the threshold to define a short cervix (≤25 mm) is reliable in predicting the risk of preterm birth. In addition the efficacy of standard treatments or whether benefit may be derived from prophylactic interventions such as a cervical cerclage is unknown.

METHODS

There are three distinct components to the CRAFT project (CRAFT-OBS, CRAFT-RCT and CRAFT-IMG).

CRAFT-OBS

Observational Study; To evaluate subsequent pregnancy risk of preterm birth in women with a prior caesarean section in established labour. This prospective study of cervical length and quantitative fetal fibronectin data will establish a predictive model of preterm birth.

CRAFT-RCT

Randomised controlled trial arm; To assess treatment for short cervix in women at high risk of preterm birth following a fully dilated caesarean section.

CRAFT-IMG

Imaging sub-study; To evaluate the use of MRI and transvaginal ultrasound imaging of micro and macrostructural cervical features which may predispose to preterm birth in women with a previous fully dilated caesarean section, such as scar position and niche.

DISCUSSION

The CRAFT project will quantify the risk of preterm birth or late miscarriage in women with previous in-labour caesarean section, define the best management and shed light on pathological mechanisms so as to improve the care we offer to women and their babies.

TRIAL REGISTRATION

CRAFT was prospectively registered on 25th November 2019 with the ISRCTN registry ( https://doi.org/10.1186/ISRCTN15068651 ).

Diffusion MRI simulation of realistic neurons with SpinDoctor and the Neuron Module

The diffusion MRI signal arising from neurons can be numerically simulated by solving the Bloch-Torrey partial differential equation. In this paper we present the Neuron Module that we implemented within the Matlab-based diffusion MRI simulation toolbox SpinDoctor. SpinDoctor uses finite element discretization and adaptive time integration to solve the Bloch-Torrey partial differential equation for general diffusion-encoding sequences, at multiple b-values and in multiple diffusion directions. In order to facilitate the diffusion MRI simulation of realistic neurons by the research community, we constructed finite element meshes for a group of 36 pyramidal neurons and a group of 29 spindle neurons whose morphological descriptions were found in the publicly available neuron repository NeuroMorpho.Org. These finite elements meshes range from having 15,163 nodes to 622,553 nodes. We also broke the neurons into the soma and dendrite branches and created finite elements meshes for these cell components. Through the Neuron Module, these neuron and cell components finite element meshes can be seamlessly coupled with the functionalities of SpinDoctor to provide the diffusion MRI signal attributable to spins inside neurons. We make these meshes and the source code of the Neuron Module available to the public as an open-source package. To illustrate some potential uses of the Neuron Module, we show numerical examples of the simulated diffusion MRI signals in multiple diffusion directions from whole neurons as well as from the soma and dendrite branches, and include a comparison of the high b-value behavior between dendrite branches and whole neurons. In addition, we demonstrate that the neuron meshes can be used to perform Monte-Carlo diffusion MRI simulations as well. We show that at equivalent accuracy, if only one gradient direction needs to be simulated, SpinDoctor is faster than a GPU implementation of Monte-Carlo, but if many gradient directions need to be simulated, there is a break-even point when the GPU implementation of Monte-Carlo becomes faster than SpinDoctor. Furthermore, we numerically compute the eigenfunctions and the eigenvalues of the Bloch-Torrey and the Laplace operators on the neuron geometries using a finite elements discretization, in order to give guidance in the choice of the space and time discretization parameters for both finite elements and Monte-Carlo approaches. Finally, we perform a statistical study on the set of 65 neurons to test some candidate biomakers that can potentially indicate the soma size. This preliminary study exemplifies the possible research that can be conducted using the Neuron Module.

NOise reduction with DIstribution Corrected (NORDIC) PCA in dMRI with complex-valued parameter-free locally low-rank processing

Diffusion-weighted magnetic resonance imaging (dMRI) has found great utility for a wide range of neuroscientific and clinical applications. However, high-resolution dMRI, which is required for improved delineation of fine brain structures and connectomics, is hampered by its low signal-to-noise ratio (SNR). Since dMRI relies on the acquisition of multiple different diffusion weighted images of the same anatomy, it is well-suited for denoising methods that utilize correlations across the image series to improve the apparent SNR and the subsequent data analysis. In this work, we introduce and quantitatively evaluate a comprehensive framework, NOise Reduction with DIstribution Corrected (NORDIC) PCA method for processing dMRI. NORDIC uses low-rank modeling of g-factor-corrected complex dMRI reconstruction and non-asymptotic random matrix distributions to remove signal components which cannot be distinguished from thermal noise. The utility of the proposed framework for denoising dMRI is demonstrated on both simulations and experimental data obtained at 3 Tesla with different resolutions using human connectome project style acquisitions. The proposed framework leads to substantially enhanced quantitative performance for estimating diffusion tractography related measures and for resolving crossing fibers as compared to a conventional/state-of-the-art dMRI denoising method.

The connections of the insular VEN area in great apes: A histologically-guided ex vivo diffusion tractography study

We mapped the connections of the insular von Economo neuron (VEN) area in ex vivo brains of a bonobo, an orangutan and two gorillas with high angular resolution diffusion MRI imaging acquired in 36 h imaging sessions for each brain. The apes died of natural causes without neurological disorders. The localization of the insular VEN area was based on cresyl violet-stained histological sections from each brain that were coregistered with structural and diffusion images from the same individuals. Diffusion MRI tractography showed that the insular VEN area is connected with olfactory, gustatory, visual and other sensory systems, as well as systems for the mediation of appetite, reward, aversion and motivation. The insular VEN area in apes is most strongly connected with frontopolar cortex, which could support their capacity to choose voluntarily among alternative courses of action particularly in exploring for food resources. The frontopolar cortex may also support their capacity to take note of potential resources for harvesting in the future (prospective memory). All of these faculties may support insight and volitional choice when contemplating courses of action as opposed to rule-based decision-making.

Microglial Density Alters Measures of Axonal Integrity and Structural Connectivity

Diffusion tensor imaging (DTI) has fundamentally transformed how we interrogate diseases and disorders of the brain in neuropsychiatric illness. DTI and recently developed multicompartment diffusion-weighted imaging (MC-DWI) techniques, such as NODDI (neurite orientation dispersion and density imaging), measure diffusion anisotropy presuming a static neuroglial environment; however, microglial morphology and density are highly dynamic in psychiatric illness, and how alterations in microglial density might influence intracellular measures of diffusion anisotropy in DTI and MC-DWI brain microstructure is unknown. To address this question, DTI and MC-DWI studies of murine brains depleted of microglia were performed, revealing significant alterations in axonal integrity and fiber tractography in DTI and in commonly used MC-DWI models. With accumulating evidence of the role of microglia in neuropsychiatric illness, our findings uncover the unexpected contribution of microglia to measures of axonal integrity and structural connectivity and provide unanticipated insights into the potential influence of microglia in diffusion imaging studies of neuropsychiatric disease.

Double diffusion encoding and applications for biomedical imaging

Diffusion Magnetic Resonance Imaging (dMRI) is one of the most important contemporary non-invasive modalities for probing tissue structure at the microscopic scale. The majority of dMRI techniques employ standard single diffusion encoding (SDE) measurements, covering different sequence parameter ranges depending on the complexity of the method. Although many signal representations and biophysical models have been proposed for SDE data, they are intrinsically limited by a lack of specificity. Advanced dMRI methods have been proposed to provide additional microstructural information beyond what can be inferred from SDE. These enhanced contrasts can play important roles in characterizing biological tissues, for instance upon diseases (e.g. neurodegenerative, cancer, stroke), aging, learning, and development. In this review we focus on double diffusion encoding (DDE), which stands out among other advanced acquisitions for its versatility, ability to probe more specific diffusion correlations, and feasibility for preclinical and clinical applications. Various DDE methodologies have been employed to probe compartment sizes (Section 3), decouple the effects of microscopic diffusion anisotropy from orientation dispersion (Section 4), probe displacement correlations, study exchange, or suppress fast diffusing compartments (Section 6). DDE measurements can also be used to improve the robustness of biophysical models (Section 5) and study intra-cellular diffusion via magnetic resonance spectroscopy of metabolites (Section 7). This review discusses all these topics as well as important practical aspects related to the implementation and contrast in preclinical and clinical settings (Section 9) and aims to provide the readers a guide for deciding on the right DDE acquisition for their specific application.

Magnetic resonance diffusion tensor tractography of a midbrain auditory circuit in Alligator

Knowledge of brain circuitry is critical for understanding the organization, function, and evolution of central nervous systems. Most commonly, brain connections have been elucidated using histological and experimental methods that require animal sacrifice. On the other hand, magnetic resonance diffusion tensor imaging and associated tractography have emerged as a preferred method to noninvasively visualize brain white matter tracts. However, existing studies have primarily examined large, heavily myelinated fiber tracts. Whether tractography can visualize fiber bundles that contain thin and poorly myelinated axons is uncertain. To address this question, the midbrain auditory pathway to the thalamus was investigated in Alligator. This species was chosen because of its evolutionary importance as it is the reptilian group most closely related to birds and because its brain contains many thin and poorly myelinated tracts. Furthermore, this auditory pathway is well documented in other reptiles, including a related crocodilian. Histological observations and experimental determination of anterograde connections confirmed this path in Alligator. Tractography identified these tracts in Alligator and provided a 3-dimensional picture that accurately identified the neural elements of this circuit. In addition, tractography identified one possible unrecognized pathway. These results demonstrate that tractography can visualize circuits containing thin, poorly myelinated fibers. These findings open the door for future studies to examine these types of pathways in other vertebrates.

In vivo microstructural heterogeneity of white matter lesions in healthy elderly and Alzheimer's disease participants using tissue compositional analysis of diffusion MRI data

White matter hyperintensities (WMH) are regions of high signal intensity typically identified on fluid attenuated inversion recovery (FLAIR). Although commonly observed in elderly individuals, they are more prevalent in Alzheimer's disease (AD) patients. Given that WMH appear relatively homogeneous on FLAIR, they are commonly partitioned into location- or distance-based classes when investigating their relevance to disease. Since pathology indicates that such lesions are often heterogeneous, probing their microstructure in vivo may provide greater insight than relying on such arbitrary classification schemes. In this study, we investigated WMH in vivo using an advanced diffusion MRI method known as single-shell 3-tissue constrained spherical deconvolution (SS3T-CSD), which models white matter microstructure while accounting for grey matter and CSF compartments. Diffusion MRI data and FLAIR images were obtained from AD (n = 48) and healthy elderly control (n = 94) subjects. WMH were automatically segmented, and classified: (1) as either periventricular or deep; or (2) into three distance-based contours from the ventricles. The 3-tissue profile of WMH enabled their characterisation in terms of white matter-, grey matter-, and fluid-like characteristics of the diffusion signal. Our SS3T-CSD findings revealed substantial heterogeneity in the 3-tissue profile of WMH, both within lesions and across the various classes. Moreover, this heterogeneity information indicated that the use of different commonly used WMH classification schemes can result in different disease-based conclusions. We conclude that future studies of WMH in AD would benefit from inclusion of microstructural information when characterising lesions, which we demonstrate can be performed in vivo using SS3T-CSD.

Super-Resolved q-Space deep learning with uncertainty quantification

Diffusion magnetic resonance imaging (dMRI) provides a noninvasive method for measuring brain tissue microstructure. q-Space deep learning(q-DL) methods have been developed to accurately estimate tissue microstructure from dMRI scans acquired with a reduced number of diffusion gradients. In these methods, deep networks are trained to learn the mapping directly from diffusion signals to tissue microstructure. However, the quality of tissue microstructure estimation can be limited not only by the reduced number of diffusion gradients but also by the low spatial resolution of typical dMRI acquisitions. Therefore, in this work we extend q-DL to super-resolved tissue microstructure estimation and propose super-resolvedq-DL (SR-q-DL), where deep networks are designed to map low-resolution diffusion signals undersampled in the q-space to high-resolution tissue microstructure. Specifically, we use a patch-based strategy, where a deep network takes low-resolution patches of diffusion signals as input and outputs high-resolution tissue microstructure patches. The high-resolution patches are then combined to obtain the final high-resolution tissue microstructure map. Motivated by existing q-DL methods, we integrate the sparsity of diffusion signals in the network design, which comprises two functional components. The first component computes sparse representation of diffusion signals for the low-resolution input patch, and the second component maps the low-resolution sparse representation to high-resolution tissue microstructure. The weights in the two components are learned jointly and the trained network performs end-to-end tissue microstructure estimation. In addition to SR-q-DL, we further propose probabilistic SR-q-DL, which can quantify the uncertainty of the network output as well as achieve improved estimation accuracy. In probabilistic SR-q-DL, a deep ensemble strategy is used. Specifically, the deep network for SR-q-DL is revised to produce not only tissue microstructure estimates but also the uncertainty of the estimates. Then, multiple deep networks are trained and their results are fused for the final prediction of high-resolution tissue microstructure and uncertainty quantification. The proposed method was evaluated on two independent datasets of brain dMRI scans. Results indicate that our approach outperforms competing methods in terms of estimation accuracy. In addition, uncertainty measures provided by our method correlate with estimation errors, which indicates potential application of the proposed uncertainty quantification method in brain studies.

The present and the future of microstructure MRI: From a paradigm shift to normal science

The aspiration of imaging tissue microstructure with MRI is to uncover micrometer-scale tissue features within millimeter-scale imaging voxels, in vivo. This kind of super-resolution has fueled a paradigm shift within the biomedical imaging community. However, what feels like an ongoing revolution in MRI, has been conceptually experienced in physics decades ago; from this point of view, our current developments can be seen as Thomas Kuhn's "normal science" stage of progress. While the concept of model-based quantification below the nominal imaging resolution is not new, its possibilities in neuroscience and neuroradiology are only beginning to be widely appreciated. This disconnect calls for communicating the progress of tissue microstructure MR imaging to its potential users. Here, a number of recent research developments are outlined in terms of the overarching concept of coarse-graining the tissue structure over an increasing diffusion length. A variety of diffusion models and phenomena are summarized on the phase diagram of diffusion MRI, with the unresolved problems and future directions corresponding to its unexplored domains.

Low-frequency dominant electrical conductivity imaging of in vivo human brain using high-frequency conductivity at Larmor-frequency and spherical mean diffusivity without external injection current

Diffusion weighted imaging based on random Brownian motion of water molecules within a voxel provides information on the micro-structure of biological tissues through water molecule diffusivity. As the electrical conductivity is primarily determined by the concentration and mobility of ionic charge carriers, the macroscopic electrical conductivity of biological tissues is also related to the diffusion of electrical ions. This paper aims to investigate the low-frequency electrical conductivity by relying on a pre-defined biological model that separates the brain into the intracellular (restricted) and extracellular (hindered) compartments. The proposed method uses B1 mapping technique, which provides a high-frequency conductivity distribution at Larmor frequency, and the spherical mean technique, which directly estimates the microscopic tissue structure based on the water molecule diffusivity and neurite orientation distribution. The total extracellular ion concentration, which is separated from the high-frequency conductivity, is recovered using the estimated diffusivity parameters and volume fraction in each compartment. We propose a method to reconstruct the low-frequency dominant conductivity tensor by taking into consideration the extracted extracellular diffusion tensor map and the reconstructed electrical parameters. To demonstrate the reliability of the proposed method, we conducted two phantom experiments. The first one used a cylindrical acrylic cage filled with an agar in the background region and four anomalies for the effect of ion concentration on the electrical conductivity. The other experiment, in which the effect of ion mobility on the conductivity was verified, used cell-like materials with thin insulating membranes suspended in an electrolyte. Animal and human brain experiments were conducted to visualize the low-frequency dominant conductivity tensor images. The proposed method using a conventional MRI scanner can predict the internal current density map in the brain without directly injected external currents.

Neuroanatomical underpinning of diffusion kurtosis measurements in the cerebral cortex of healthy macaque brains

PURPOSE

To investigate the neuroanatomical underpinning of healthy macaque brain cortical microstructure measured by diffusion kurtosis imaging (DKI), which characterizes non-Gaussian water diffusion.

METHODS

High-resolution DKI was acquired from 6 postmortem macaque brains. Neurofilament density (ND) was quantified based on structure tensor from neurofilament histological images of a different macaque brain sample. After alignment of DKI-derived mean kurtosis (MK) maps to the histological images, MK and histology-based ND were measured at corresponding regions of interests characterized by distinguished cortical MK values in the prefrontal/precentral-postcentral and temporal cortices. Pearson correlation was performed to test significant correlation between these cortical MK and ND measurements.

RESULTS

Heterogeneity of cortical MK across different cortical regions was revealed, with significantly and consistently higher MK measurements in the prefrontal/precentral-postcentral cortex compared to those in the temporal cortex across all six scanned macaque brains. Corresponding higher ND measurements in the prefrontal/precentral-postcentral cortex than in the temporal cortex were also found. The heterogeneity of cortical MK is associated with heterogeneity of histology-based ND measurements, with significant correlation between cortical MK and corresponding ND measurements (P < .005).

CONCLUSION

These findings suggested that DKI-derived MK can potentially be an effective noninvasive biomarker quantifying underlying neuroanatomical complexity inside the cerebral cortical mantle for clinical and neuroscientific research.

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

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

Lattice Boltzmann method for simulation of diffusion magnetic resonance imaging physics in multiphase tissue models

We report an implementation of the lattice Boltzmann method (LBM) to integrate the Bloch-Torrey equation, which describes the evolution of the transverse magnetization vector and the fate of the signal of diffusion magnetic resonance imaging (dMRI). Motivated by the need to interpret dMRI experiments in biological tissues, and to offset the small time-step limitation of classical LBM, a hybrid LBM scheme is introduced and implemented to solve the Bloch-Torrey equation. A membrane boundary condition is presented which is able to accurately represent the effects of thin curvilinear membranes typically found in biological tissues. As implemented, the hybrid LBM scheme accommodates piece-wise uniform transport, dMRI parameters, periodic and mirroring outer boundary conditions, and finite membrane permeabilities on non-boundary-conforming inner boundaries. By comparing with analytical solutions of limiting cases, we demonstrate that the hybrid LBM scheme is more accurate than the classical LBM scheme. The proposed explicit LBM scheme maintains second-order spatial accuracy, stability, and first-order temporal accuracy for a wide range of parameters. The parallel implementation of the hybrid LBM code in a multi-CPU computer system, as well as on GPUs, is straightforward and efficient. Along with offering certain advantages over finite element or Monte Carlo schemes, the proposed hybrid LBM constitutes a flexible scheme that can by easily adapted to model more complex interfacial conditions and physics in heterogeneous multiphase tissue models and to accommodate sophisticated dMRI sequences.

Characterizing Intra-soma Diffusion with Spherical Mean Spectrum Imaging

Most brain microstructure models are dedicated to the quantification of white matter microstructure, using for example sticks, cylinders, and zeppelins to model intra- and extra-axonal environments. Gray matter presents unique micro-architecture with cell bodies (somas) exhibiting diffusion characteristics that differ from axons in white matter. In this paper, we introduce a method to quantify soma microstructure, giving measures such as volume fraction, diffusivity, and kurtosis. Our method captures a spectrum of diffusion patterns and scales and does not rely on restrictive model assumptions. We show that our method yields unique and meaningful contrasts that are in agreement with histological data. We demonstrate its application in the mapping of the distinct spatial patterns of soma density in the cortex.

Challenges for biophysical modeling of microstructure

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

Probing myelin content of the human brain with MRI: A review

Rapid and efficient transmission of electric signals among neurons of vertebrates is ensured by myelin-insulating sheaths surrounding axons. Human cognition, sensation, and motor functions rely on the integrity of these layers, and demyelinating diseases often entail serious cognitive and physical impairments. Magnetic resonance imaging radically transformed the way these disorders are monitored, offering an irreplaceable tool to noninvasively examine the brain structure. Several advanced techniques based on MRI have been developed to provide myelin-specific contrasts and a quantitative estimation of myelin density in vivo. Here, the vast offer of acquisition strategies developed to date for this task is reviewed. Advantages and pitfalls of the different approaches are compared and discussed.

Diffusion-weighted MRI in neurodegenerative and psychiatric animal models: Experimental strategies and main outcomes

Preclinical MRI approaches constitute a key tool to study a wide variety of neurological and psychiatric illnesses, allowing a more direct investigation of the disorder substrate and, at the same time, the possibility of back-translating such findings to human subjects. However, the lack of consensus on the optimal experimental scheme used to acquire the data has led to relatively high heterogeneity in the choice of protocols, which can potentially impact the comparison between results obtained by different groups, even using the same animal model. This is especially true for diffusion-weighted MRI data, where certain experimental choices can impact not only on the accuracy and precision of the extracted biomarkers, but also on their biological meaning. With this in mind, we extensively examined preclinical imaging studies that used diffusion-weighted MRI to investigate neurodegenerative, neurodevelopmental and psychiatric disorders in rodent models. In this review, we discuss the main findings for each preclinical model, with a special focus on the analysis and comparison of the different acquisition strategies used across studies and their impact on the heterogeneity of the findings.

NODDI in clinical research

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

ConFiG: Contextual Fibre Growth to generate realistic axonal packing for diffusion MRI simulation

This paper presents Contextual Fibre Growth (ConFiG), an approach to generate white matter numerical phantoms by mimicking natural fibre genesis. ConFiG grows fibres one-by-one, following simple rules motivated by real axonal guidance mechanisms. These simple rules enable ConFiG to generate phantoms with tuneable microstructural features by growing fibres while attempting to meet morphological targets such as user-specified density and orientation distribution. We compare ConFiG to the state-of-the-art approach based on packing fibres together by generating phantoms in a range of fibre configurations including crossing fibre bundles and orientation dispersion. Results demonstrate that ConFiG produces phantoms with up to 20% higher densities than the state-of-the-art, particularly in complex configurations with crossing fibres. We additionally show that the microstructural morphology of ConFiG phantoms is comparable to real tissue, producing diameter and orientation distributions close to electron microscopy estimates from real tissue as well as capturing complex fibre cross sections. Signals simulated from ConFiG phantoms match real diffusion MRI data well, showing that ConFiG phantoms can be used to generate realistic diffusion MRI data. This demonstrates the feasibility of ConFiG to generate realistic synthetic diffusion MRI data for developing and validating microstructure modelling approaches.

•

We present ConFiG, a biologically motivated numerical phantom generator for white matter.

•

ConFiG produces phantoms with state-of-the-art density and realistic microstructure.

•

Diffusion MRI simulations in ConFiG phantoms are comparable to real dMRI signals.

A new method for accurate in vivo mapping of human brain connections using microstructural and anatomical information

Diffusion magnetic resonance imaging is a noninvasive imaging modality that has been extensively used in the literature to study the neuronal architecture of the brain in a wide range of neurological conditions using tractography. However, recent studies highlighted that the anatomical accuracy of the reconstructions is inherently limited and challenged its appropriateness. Several solutions have been proposed to tackle this issue, but none of them proved effective to overcome this fundamental limitation. In this work, we present a novel processing framework to inject into the reconstruction problem basic prior knowledge about brain anatomy and its organization and evaluate its effectiveness using both simulated and real human brain data. Our results indicate that our proposed method dramatically increases the accuracy of the estimated brain networks and, thus, represents a major step forward for the study of connectivity.

Diffusion Imaging in the Post HCP Era

Diffusion imaging is a critical component in the pursuit of developing a better understanding of the human brain. Recent technical advances promise enabling the advancement in the quality of data that can be obtained. In this review the context for different approaches relative to the Human Connectome Project are compared. Significant new gains are anticipated from the use of high-performance head gradients. These gains can be particularly large when the high-performance gradients are employed together with ultrahigh magnetic fields. Transmit array designs are critical in realizing high accelerations in diffusion-weighted (d)MRI acquisitions, while maintaining large field of view (FOV) coverage, and several techniques for optimal signal-encoding are now available. Reconstruction and processing pipelines that precisely disentangle the acquired neuroanatomical information are established and provide the foundation for the application of deep learning in the advancement of dMRI for complex tissues. Level of Evidence: 3 Technical Efficacy Stage: Stage 3.

Disentangling white-matter damage from physiological fibre orientation dispersion in multiple sclerosis

Multiple sclerosis leads to diffuse damage of the central nervous system, affecting also the normal-appearing white matter. Demyelination and axonal degeneration reduce regional fractional anisotropy in normal-appearing white matter, which can be routinely mapped with diffusion tensor imaging. However, the standard fractional anisotropy metric is also sensitive to physiological variations in orientation dispersion of white matter fibres. This complicates the detection of disease-related damage in large parts of cerebral white matter where microstructure physiologically displays a high degree of fibre dispersion. To resolve this ambiguity, we employed a novel tensor-valued encoding method for diffusion MRI, which yields a microscopic fractional anisotropy metric that is unaffected by regional variations in orientation dispersion. In 26 patients with relapsing-remitting multiple sclerosis, 14 patients with primary-progressive multiple sclerosis and 27 age-matched healthy controls, we compared standard fractional anisotropy mapping with the novel microscopic fractional anisotropy mapping method, focusing on normal-appearing white matter. Mean microscopic fractional anisotropy and standard fractional anisotropy of normal-appearing white matter were significantly reduced in both patient groups relative to healthy controls, but microscopic fractional anisotropy yielded a better reflection of disease-related white-matter alterations. The reduction in mean microscopic fractional anisotropy showed a significant positive linear relationship with physical disability, as reflected by the expanded disability status scale. Mean reduction of microscopic fractional anisotropy in normal-appearing white matter also scaled positively with individual cognitive dysfunction, as measured with the symbol digit modality test. Mean microscopic fractional anisotropy reduction in normal-appearing white matter also showed a positive relationship with total white-matter lesion load as well as lesion load in specific tract systems. None of these relationships between normal-appearing white-matter microstructure and clinical, cognitive or structural measures emerged when using mean fractional anisotropy. Together, the results provide converging evidence that microscopic fractional anisotropy mapping substantially advances the assessment of cerebral white matter in multiple sclerosis by disentangling microstructure damage from variations in physiological fibre orientation dispersion at the stage of data acquisition. Since tensor-valued encoding can be implemented in routine diffusion MRI, microscopic fractional anisotropy mapping bears considerable potential for the future assessment of disease progression in normal-appearing white matter in both relapsing-remitting and progressive forms of multiple sclerosis as well as other white-matter-related brain diseases.

Leveraging multi-shell diffusion for studies of brain development in youth and young adulthood

•

Multi-shell imaging sequences may improve sensitivity to developmental effects.

•

Models that leverage multi-shell information are often less sensitive to the confounding effects of motion.

•

Multi-shell sequences and models that leverage this data may be of particular utility for studying the developing brain.

Diffusion weighted imaging (DWI) has advanced our understanding of brain microstructure evolution over development. Recently, the use of multi-shell diffusion imaging sequences has coincided with advances in modeling the diffusion signal, such as Neurite Orientation Dispersion and Density Imaging (NODDI) and Laplacian-regularized Mean Apparent Propagator MRI (MAPL). However, the relative utility of recently-developed diffusion models for understanding brain maturation remains sparsely investigated. Additionally, despite evidence that motion artifact is a major confound for studies of development, the vulnerability of metrics derived from contemporary models to in-scanner motion has not been described. Accordingly, in a sample of 120 youth and young adults (ages 12–30) we evaluated metrics derived from diffusion tensor imaging (DTI), NODDI, and MAPL for associations with age and in-scanner head motion at multiple scales. Specifically, we examined mean white matter values, white matter tracts, white matter voxels, and connections in structural brain networks. Our results revealed that multi-shell diffusion imaging data can be leveraged to robustly characterize neurodevelopment, and demonstrate stronger age effects than equivalent single-shell data. Additionally, MAPL-derived metrics were less sensitive to the confounding effects of head motion. Our findings suggest that multi-shell imaging data and contemporary modeling techniques confer important advantages for studies of neurodevelopment.

SANDI: A compartment-based model for non-invasive apparent soma and neurite imaging by diffusion MRI

This work introduces a compartment-based model for apparent cell body (namely soma) and neurite density imaging (SANDI) using non-invasive diffusion-weighted MRI (DW-MRI). The existing conjecture in brain microstructure imaging through DW-MRI presents water diffusion in white (WM) and gray (GM) matter as restricted diffusion in neurites, modelled by infinite cylinders of null radius embedded in the hindered extra-neurite water. The extra-neurite pool in WM corresponds to water in the extra-axonal space, but in GM it combines water in the extra-cellular space with water in soma. While several studies showed that this microstructure model successfully describe DW-MRI data in WM and GM at b ​≤ ​3,000 ​s/mm² (or 3 ​ms/μm²), it has been also shown to fail in GM at high b values (b≫3,000 ​s/mm² or 3 ​ms/μm²). Here we hypothesise that the unmodelled soma compartment (i.e. cell body of any brain cell type: from neuroglia to neurons) may be responsible for this failure and propose SANDI as a new model of brain microstructure where soma of any brain cell type is explicitly included. We assess the effects of size and density of soma on the direction-averaged DW-MRI signal at high b values and the regime of validity of the model using numerical simulations and comparison with experimental data from mouse (b_max ​= ​40,000 ​s/mm², or 40 ​ms/μm²) and human (b_max ​= ​10,000 ​s/mm², or 10 ​ms/μm²) brain. We show that SANDI defines new contrasts representing complementary information on the brain cyto- and myelo-architecture. Indeed, we show maps from 25 healthy human subjects of MR soma and neurite signal fractions, that remarkably mirror contrasts of histological images of brain cyto- and myelo-architecture. Although still under validation, SANDI might provide new insight into tissue architecture by introducing a new set of biomarkers of potential great value for biomedical applications and pure neuroscience.

The role of diffusion tractography in refining glial tumor resection

Primary brain tumors are notoriously hard to resect surgically. Due to their infiltrative nature, finding the optimal resection boundary without damaging healthy tissue can be challenging. One potential tool to help make this decision is diffusion-weighted magnetic resonance imaging (dMRI) tractography. dMRI exploits the diffusion of water molecule along axons to generate a 3D modelization of the white matter bundles in the brain. This feature is particularly useful to visualize how a tumor affects its surrounding white matter and plan a surgical path. This paper reviews the different ways in which dMRI can be used to improve brain tumor resection, its benefits and also its limitations. We expose surgical tools that can be paired with dMRI to improve its impact on surgical outcome, such as loading the 3D tractography in the neuronavigation system and direct electrical stimulation to validate the position of the white matter bundles of interest. We also review articles validating dMRI findings using other anatomical investigation techniques, such as postmortem dissections, manganese-enhanced MRI, electrophysiological stimulations, and phantom studies with known ground truth. We will be discussing the areas of the brain where dMRI performs well and where the future challenges are. We will conclude this review with suggestions and take home messages for neurosurgeons, tractographers, and vendors for advancing the field and on how to benefit from tractography's use in clinical practice.

Double Diffusion Encoding for Probing Radiation-Induced Microstructural Changes in a Tumor Model: A Proof-of-Concept Study With Comparison to the Apparent Diffusion Coefficient and Histology

BACKGROUND

Microstructure analyses are gaining interest in cancer MRI as an alternative to the conventional apparent diffusion coefficient (ADC), of which the determinants remain unclear.

PURPOSE

To assess the sensitivity of parameters calculated from a double diffusion encoding (DDE) sequence to changes in a tumor's microstructure early after radiotherapy and to compare them with ADC and histology.

STUDY TYPE

Cohort study on experimental tumors.

ANIMAL MODEL

Sixteen WAG/Rij rats grafted with one rhabdomyosarcoma fragment in each thigh. Thirty-one were imaged at days 1 and 4, of which 17 tumors received a 20 Gy radiation dose after the first imagery.

FIELD STRENGTH/SEQUENCE

3T. Diffusion-weighted imaging, DDE with flow compensated, and noncompensated measurements.

ASSESSMENTS

1) To compare, after irradiation, DDE-derived parameters (intracellular fraction, cell size, and cell density) to their histological counterparts (fraction of stained area, minimal Feret diameter, and nuclei count, respectively). 2) To compare percentage changes in DDE-derived parameters and ADC. 3) To evaluate the evolution of DDE-derived parameters describing perfusion.

STATISTICAL TESTS

Wilcoxon rank sum test.

RESULTS

1) Intracellular fraction, cell size, and cell density were respectively lower (-24%, P < 0.001), higher (+7.5%, P < 0.001) and lower (-38%, P < 0.001) in treated tumors as compared to controls. Fraction of stained area, minimal Feret diameter, and nuclei count were respectively lower (-20%, P < 0.001), higher (+28%, P < 0.001), and lower (-34%, P < 0.001) in treated tumors. 2) The magnitude of ADC's percentage change due to irradiation (16.4%) was superior to the one of cell size (8.4%, P < 0.01) but inferior to those of intracellular fraction (35.5%, P < 0.001) and cell density (42%, P < 0.001). 3) After treatment, the magnitude of the vascular fraction's decrease was higher than the increase of flow velocity (33.3%, vs. 13.3%, P < 0.001).

DATA CONCLUSION

The DDE sequence allows quantitatively monitoring the effects of radiotherapy on a tumor's microstructure, whereas ADC only reveals global changes.

EVIDENCE LEVEL

2.

TECHNICAL EFFICACY

Stage 4. J. Magn. Reson. Imaging 2020;52:941-951.

Towards unconstrained compartment modeling in white matter using diffusion-relaxation MRI with tensor-valued diffusion encoding

PURPOSE

To optimize diffusion-relaxation MRI with tensor-valued diffusion encoding for precise estimation of compartment-specific fractions, diffusivities, and T₂ values within a two-compartment model of white matter, and to explore the approach in vivo.

METHODS

Sampling protocols featuring different b-values (b), b-tensor shapes (b_Δ ), and echo times (TE) were optimized using Cramér-Rao lower bounds (CRLB). Whole-brain data were acquired in children, adults, and elderly with white matter lesions. Compartment fractions, diffusivities, and T₂ values were estimated in a model featuring two microstructural compartments represented by a "stick" and a "zeppelin."

RESULTS

Precise parameter estimates were enabled by sampling protocols featuring seven or more "shells" with unique b/b_Δ /TE-combinations. Acquisition times were approximately 15 minutes. In white matter of adults, the "stick" compartment had a fraction of approximately 0.5 and, compared with the "zeppelin" compartment, featured lower isotropic diffusivities (0.6 vs. 1.3 μm² /ms) but higher T₂ values (85 vs. 65 ms). Children featured lower "stick" fractions (0.4). White matter lesions exhibited high "zeppelin" isotropic diffusivities (1.7 μm² /ms) and T₂ values (150 ms).

CONCLUSIONS

Diffusion-relaxation MRI with tensor-valued diffusion encoding expands the set of microstructure parameters that can be precisely estimated and therefore increases their specificity to biological quantities.

Self-navigation for 3D multishot EPI with data-reference

PURPOSE

In this study, we sought to develop a self-navigation strategy for improving the reconstruction of diffusion weighted 3D multishot echo planar imaging (EPI). We propose a method for extracting the phase correction information from the acquisition itself, eliminating the need for a 2D navigator, further accelerating the acquisition.

METHODS

In-vivo acquisitions at 3T with 0.9 mm and 1.5 mm isotropic resolutions were used to evaluate the performance of the self-navigation strategy. Sensitivity to motion was tested using a large difference in pitch position of the head. Using a multishell diffusion weighted acquisition, tractography results were obtained at (0.9 mm)³ to validate the quality with conventional acquisition.

RESULTS

The use of 3D multislab EPI with self-navigation enables 3D diffusion-weighted spin echo EPI acquisitions that have the same efficiency as 2D single-shot acquisition. For matched acquisition time the image signal-to-noise ratio (SNR) between 3D and 2D acquisition is shown to be comparable for whole-brain coverage with (1.5 mm)³ resolution and for (0.9 mm)³ resolution the 3D acquisition has higher SNR than what can be obtained with 2D acquisitions using current state-of-art multiband techniques. The self-navigation technique was shown to be stable under inter-volume motion. In tractography analysis, the higher resolution afforded by our technique enabled clear delineation of the tapetum and posterior corona radiata.

CONCLUSION

The proposed self-navigation approach utilized a self-consistent phase in 3D diffusion weighted acquisitions. Its efficiency and stability were demonstrated for a plurality of common acquisitions. The proposed self-navigation approach allows for faster acquisition of 3D multishot EPI desirable for large field of view and/or higher resolution.

Time-dependent diffusion in undulating thin fibers: Impact on axon diameter estimation

Diffusion MRI may enable non-invasive mapping of axonal microstructure. Most approaches infer axon diameters from effects of time-dependent diffusion on the diffusion-weighted MR signal by modeling axons as straight cylinders. Axons do not, however, propagate in straight trajectories, and so far the impact of the axonal trajectory on diameter estimation has been insufficiently investigated. Here, we employ a toy model of axons, which we refer to as the undulating thin fiber model, to analyze the impact of undulating trajectories on the time dependence of diffusion. We study time-dependent diffusion in the frequency domain and characterize the diffusion spectrum by its height, width, and low-frequency behavior (power law exponent). Results show that microscopic orientation dispersion of the thin fibers is the main parameter that determines the characteristics of the diffusion spectra. At lower frequencies (longer diffusion times), straight cylinders and undulating thin fibers can have virtually identical spectra. If the straight-cylinder assumption is used to interpret data from undulating thin axons, the diameter is overestimated by an amount proportional to the undulation amplitude and microscopic orientation dispersion of the fibers. At higher frequencies (shorter diffusion times), spectra from cylinders and undulating thin fibers differ. The low-frequency behavior of the spectra from the undulating thin fibers may also differ from that of cylinders, because the power law exponent of undulating fibers can reach values below 2 for experimentally relevant frequency ranges. In conclusion, we argue that the non-straight nature of axonal trajectories should not be overlooked when analyzing and interpreting diffusion MRI data.

Revisiting double diffusion encoding MRS in the mouse brain at 11.7T: Which microstructural features are we sensitive to?

Brain metabolites, such as N-acetylaspartate or myo-inositol, are constantly probing their local cellular environment under the effect of diffusion. Diffusion-weighted NMR spectroscopy therefore presents unparalleled potential to yield cell-type specific microstructural information. Double diffusion encoding (DDE) consists in applying two diffusion blocks, where gradient's direction in the second block is varied during the course of the experiment. Unlike single diffusion encoding, DDE measurements at long mixing time display some angular modulation of the signal amplitude which reflects microscopic anisotropy (μA), while requiring relatively low gradient strength. This angular dependence has been formerly used to quantify cell fiber diameter using a model of isotropically oriented infinite cylinders. However, how additional features of the cell microstructure (such as cell body diameter, fiber length and branching) may also influence the DDE signal has been little explored. Here, we used a cryoprobe as well as state-of-the-art post-processing to perform DDE acquisitions with high accuracy and precision in the mouse brain at 11.7 ​T. We then compared our results to simulated DDE datasets obtained in various 3D cell models in order to pinpoint which features of cell morphology may influence the most the angular dependence of the DDE signal. While the infinite cylinder model poorly fits our experimental data, we show that incorporating branched fiber structure in our model allows more realistic interpretation of the DDE signal. Lastly, data acquired in the short mixing time regime suggest that some sensitivity to cell body diameter might be retrieved, although additional experiments would be required to further support this statement.

Constant gradient FEXSY: A time-efficient method for measuring exchange

Filter-Exchange NMR Spectroscopy (FEXSY) is a method for measurement of apparent transmembranal water exchange rates. The experiment is comprised of two co-linear sequential pulsed-field gradient (PFG) blocks, separated by a mixing period in which exchange takes place. The first block remains constant and serves as a diffusion-based filter that removes signal coming from fast-diffusing water. The mixing time and the gradient area (q-value) of the second block are varied on repeated iterations to produce a 2D data set that is analyzed using a bi-compartmental model which assumes that intra- and extra-cellular water are slow and fast diffusing, respectively. Here we suggest a variant of the FEXSY method in which measurements for different mixing times are taken at a constant gradient. This Constant Gradient FEXSY (CG-FEXSY) allows for the determination of the exchange rate by using a smaller 1D data set, so that the same information can be gathered during a considerably shorter scan time. Furthermore, in the limit of high diffusion weighting, such that the extra-cellular water signal is removed while the intra-cellular signal is retained, CG-FEXSY also allows for determination of the intra-cellular mean residence time (MRT). The theoretical results are validated on a living yeast cells sample and on a fixed porcine optic nerve, where the values obtained from the two methods are shown to be in agreement.

An improved deep network for tissue microstructure estimation with uncertainty quantification

Deep learning based methods have improved the estimation of tissue microstructure from diffusion magnetic resonance imaging (dMRI) scans acquired with a reduced number of diffusion gradients. These methods learn the mapping from diffusion signals in a voxel or patch to tissue microstructure measures. In particular, it is beneficial to exploit the sparsity of diffusion signals jointly in the spatial and angular domains, and the deep network can be designed by unfolding iterative processes that adaptively incorporate historical information for sparse reconstruction. However, the number of network parameters is huge in such a network design, which could increase the difficulty of network training and limit the estimation performance. In addition, existing deep learning based approaches to tissue microstructure estimation do not provide the important information about the uncertainty of estimates. In this work, we continue the exploration of tissue microstructure estimation using a deep network and seek to address these limitations. First, we explore the sparse representation of diffusion signals with a separable spatial-angular dictionary and design an improved deep network for tissue microstructure estimation. The procedure for updating the sparse code associated with the separable dictionary is derived and unfolded to construct the deep network. Second, with the formulation of sparse representation of diffusion signals, we propose to quantify the uncertainty of network outputs with a residual bootstrap strategy. Specifically, because of the sparsity constraint in the signal representation, we perform a Lasso bootstrap strategy for uncertainty quantification. Experiments were performed on brain dMRI scans with a reduced number of diffusion gradients, where the proposed method was applied to two representative biophysical models for describing tissue microstructure and compared with state-of-the-art methods of tissue microstructure estimation. The results show that our approach compares favorably with the competing methods in terms of estimation accuracy. In addition, the uncertainty measures provided by our method correlate with estimation errors and produce reasonable confidence intervals; these results suggest potential application of the proposed uncertainty quantification method in brain studies.

Mild traumatic brain injury impacts associations between limbic system microstructure and post-traumatic stress disorder symptomatology

BACKGROUND

Post-traumatic stress disorder (PTSD) is a psychiatric disorder that afflicts many individuals, yet the neuropathological mechanisms that contribute to this disorder remain to be fully determined. Moreover, it is unclear how exposure to mild traumatic brain injury (mTBI), a condition that is often comorbid with PTSD, particularly among military personnel, affects the clinical and neurological presentation of PTSD. To address these issues, the present study explores relationships between PTSD symptom severity and the microstructure of limbic and paralimbic gray matter brain regions, as well as the impact of mTBI comorbidity on these relationships.

METHODS

Structural and diffusion MRI data were acquired from 102 male veterans who were diagnosed with current PTSD. Diffusion data were analyzed with free-water imaging to quantify average CSF-corrected fractional anisotropy (FA) and mean diffusivity (MD) in 18 limbic and paralimbic gray matter regions. Associations between PTSD symptom severity and regional average dMRI measures were examined with repeated measures linear mixed models. Associations were studied separately in veterans with PTSD only, and in veterans with PTSD and a history of military mTBI.

RESULTS

Analyses revealed that in the PTSD only cohort, more severe symptoms were associated with higher FA in the right amygdala-hippocampus complex, lower FA in the right cingulate cortex, and lower MD in the left medial orbitofrontal cortex. In the PTSD and mTBI cohort, more severe PTSD symptoms were associated with higher FA bilaterally in the amygdala-hippocampus complex, with higher FA bilaterally in the nucleus accumbens, with lower FA bilaterally in the cingulate cortex, and with higher MD in the right amygdala-hippocampus complex.

CONCLUSIONS

These findings suggest that the microstructure of limbic and paralimbic brain regions may influence PTSD symptomatology. Further, given the additional associations observed between microstructure and symptom severity in veterans with head trauma, we speculate that mTBI may exacerbate the impact of brain microstructure on PTSD symptoms, especially within regions of the brain known to be vulnerable to chronic stress. A heightened sensitivity to the microstructural environment of the brain could partially explain why individuals with PTSD and mTBI comorbidity experience more severe symptoms and poorer illness prognoses than those without a history of brain injury. The relevance of these microstructural findings to the conceptualization of PTSD as being a disorder of stress-induced neuronal connectivity loss is discussed.

Self-Immolative Activation of β-Galactosidase-Responsive Probes for In Vivo MR Imaging in Mouse Models

Our lab has developed a new series of self-immolative MR agents for the rapid detection of enzyme activity in mouse models expressing β-galactosidase (β-gal). We investigated two molecular architectures to create agents that detect β-gal activity by modulating the coordination of water to GdIII . The first is an intermolecular approach, wherein we designed several structural isomers to maximize coordination of endogenous carbonate ions. The second involves an intramolecular mechanism for q modulation. We incorporated a pendant coordinating carboxylate ligand with a 2, 4, 6, or 8 carbon linker to saturate ligand coordination to the GdIII ion. This renders the agent ineffective. We show that one agent in particular (6-C pendant carboxylate) is an extremely effective MR reporter for the detection of enzyme activity in a mouse model expressing β-gal.

Age-Related Changes of Peak Width Skeletonized Mean Diffusivity (PSMD) Across the Adult Lifespan: A Multi-Cohort Study

Parameters of water diffusion in white matter derived from diffusion-weighted imaging (DWI), such as fractional anisotropy (FA), mean, axial, and radial diffusivity (MD, AD, and RD), and more recently, peak width of skeletonized mean diffusivity (PSMD), have been proposed as potential markers of normal and pathological brain ageing. However, their relative evolution over the entire adult lifespan in healthy individuals remains partly unknown during early and late adulthood, and particularly for the PSMD index. Here, we gathered and analyzed cross-sectional diffusion tensor imaging (DTI) data from 10 population-based cohort studies in order to establish the time course of white matter water diffusion phenotypes from post-adolescence to late adulthood. DTI data were obtained from a total of 20,005 individuals aged 18.1 to 92.6 years and analyzed with the same pipeline for computing skeletonized DTI metrics from DTI maps. For each individual, MD, AD, RD, and FA mean values were computed over their FA volume skeleton, PSMD being calculated as the 90% peak width of the MD values distribution across the FA skeleton. Mean values of each DTI metric were found to strongly vary across cohorts, most likely due to major differences in DWI acquisition protocols as well as pre-processing and DTI model fitting. However, age effects on each DTI metric were found to be highly consistent across cohorts. RD, MD, and AD variations with age exhibited the same U-shape pattern, first slowly decreasing during post-adolescence until the age of 30, 40, and 50 years, respectively, then progressively increasing until late life. FA showed a reverse profile, initially increasing then continuously decreasing, slowly until the 70s, then sharply declining thereafter. By contrast, PSMD constantly increased, first slowly until the 60s, then more sharply. These results demonstrate that, in the general population, age affects PSMD in a manner different from that of other DTI metrics. The constant increase in PSMD throughout the entire adult life, including during post-adolescence, indicates that PSMD could be an early marker of the ageing process.

Modified connectivity of vulnerable brain nodes in multiple sclerosis, their impact on cognition and their discriminative value

Brain structural network modifications in multiple sclerosis (MS) seem to be clinically relevant. The discriminative ability of those changes to identify MS patients or their cognitive status remains unknown. Therefore, this study aimed to investigate connectivity changes in MS patients related to their cognitive status, and to define an automatic classification method to classify subjects as patients and healthy volunteers (HV) or as cognitively preserved (CP) and impaired (CI) patients. We analysed structural brain connectivity in 45 HV and 188 MS patients (104 CP and 84 CI). A support vector machine with k-fold cross-validation was built using the graph metrics features that best differentiate the groups (p < 0.05). Local efficiency (LE) and node strength (NS) network properties showed the largest differences: 100% and 69.7% of nodes had reduced LE and NS in CP patients compared to HV. Moreover, 55.3% and 57.9% of nodes had decreased LE and NS in CI compared to CP patients, in associative multimodal areas. The classification method achieved an accuracy of 74.8–77.2% to differentiate patients from HV, and 59.9–60.8% to discriminate CI from CP patients. Structural network integrity is widely reduced and worsens as cognitive function declines. Central network properties of vulnerable nodes can be useful to classify MS patients.

A comprehensive approach for correcting voxel-wise b-value errors in diffusion MRI

PURPOSE

In diffusion MRI, the actual b-value played out on the scanner may deviate from the nominal value due to magnetic field imperfections. A simple image-based correction method for this problem is presented.

METHODS

The apparent diffusion constant (ADC) of a water phantom was measured voxel-wise along 64 diffusion directions at b = 1000 s/mm² . The true diffusion constant of water was estimated, considering the phantom temperature. A voxel-wise correction factor, providing an effective b-value including any magnetic field deviations, was determined for each diffusion direction by relating the measured ADC to the true diffusion constant. To test the method, the measured b-value map was used to calculate the corrected voxel-wise ADC for additionally acquired diffusion data sets on the same water phantom and data sets acquired on a small water phantom at three different positions. Diffusion tensor was estimated by applying the measured b-value map to phantom and in vivo data sets.

RESULTS

The b-value-corrected ADC maps of the phantom showed the expected spatial uniformity as well as a marked improvement in consistency across diffusion directions. The b-value correction for the brain data resulted in a 5.8% and 5.5% decrease in mean diffusivity and angular differences of the primary diffusion direction of 2.71° and 0.73° inside gray and white matter, respectively.

CONCLUSION

The actual b-value deviates significantly from its nominal setting, leading to a spatially variable error in the common diffusion outcome measures. The suggested method measures and corrects these artifacts.

On the need for bundle-specific microstructure kernels in diffusion MRI

Probing microstructure with diffusion magnetic resonance imaging (dMRI) on a scale orders of magnitude below the imaging resolution relies on biophysical modelling of the signal response in the tissue. The vast majority of these biophysical models of diffusion in white matter assume that the measured dMRI signal is the sum of the signals emanating from each of the constituent compartments, each of which exhibits a distinct behaviour in the b-value and/or orientation domain. Many of these models further assume that the dMRI behaviour of the oriented compartments (e.g. the intra-axonal space) is identical between distinct fibre populations, at least at the level of a single voxel. This implicitly assumes that any potential biological differences between fibre populations are negligible, at least as far as is measurable using dMRI. Here, we validate this assumption by means of a voxel-wise, model-free signal decomposition that, under the assumption above and in the absence of noise, is shown to be rank-1. We evaluate the effect size of signal components beyond this rank-1 representation and use permutation testing to assess their significance. We conclude that in the healthy adult brain, the dMRI signal is adequately represented by a rank-1 model, implying that biologically more realistic, but mathematically more complex fascicle-specific microstructure models do not capture statistically significant or anatomically meaningful structure, even in extended high-b diffusion MRI scans.

Single- and double-Diffusion encoding MRI for studying ex vivo apparent axon diameter distribution in spinal cord white matter

Mapping average axon diameter (AAD) and axon diameter distribution (ADD) in neuronal tissues non-invasively is a challenging task that may have a tremendous effect on our understanding of the normal and diseased central nervous system (CNS). Water diffusion is used to probe microstructure in neuronal tissues, however, the different water populations and barriers that are present in these tissues turn this into a complex task. Therefore, it is not surprising that recently we have witnessed a burst in the development of new approaches and models that attempt to obtain, non-invasively, detailed microstructural information in the CNS. In this work, we aim at challenging and comparing the microstructural information obtained from single diffusion encoding (SDE) with double diffusion encoding (DDE) MRI. We first applied SDE and DDE MR spectroscopy (MRS) on microcapillary phantoms and then applied SDE and DDE MRI on an ex vivo porcine spinal cord (SC), using similar experimental conditions. The obtained diffusion MRI data were fitted by the same theoretical model, assuming that the signal in every voxel can be approximated as the superposition of a Gaussian-diffusing component and a series of restricted components having infinite cylindrical geometries. The diffusion MRI results were then compared with histological findings. We found a good agreement between the fittings and the experimental data in white matter (WM) voxels of the SC in both diffusion MRI methods. The microstructural information and apparent AADs extracted from SDE MRI were found to be similar or somewhat larger than those extracted from DDE MRI especially when the diffusion time was set to 40 ms. The apparent ADDs extracted from SDE and DDE MRI show reasonable agreement but somewhat weaker correspondence was observed between the diffusion MRI results and histology. The apparent subtle differences between the microstructural information obtained from SDE and DDE MRI are briefly discussed.