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

Modelling white matter microstructure using diffusion OGSE MRI: Model and analysis choices

Alterations in white matter (WM) microstructure of the central nervous system have been shown to be pathophysiological presentations of various neurodegenerative disorders. Current methods for measuring such WM features require ex vivo tissue samples analyzed using electron microscopy. Magnetic Resonance Imaging (MRI) diffusion-weighted pulse sequences provide a non-invasive tool for estimating such microstructural features in vivo. The current project investigated the use of two methods of analysis, including the ROI-based (Region of Interest, RBA) and voxel-based analysis (VBA), as well as four mathematical models of WM microstructure, including the ActiveAx Frequency-Independent Extra-Axonal Diffusion (AAI), ActiveAx Frequency-Dependent Extra-Axonal Diffusion (AAD), AxCaliber Frequency-Independent Extra-Axonal Diffusion (ACI), and AxCaliber Frequency-Dependent Extra-Axonal Diffusion (ACD) models. Two mice samples imaged at 7 T and 15.2 T were analyzed. Both the AAI and AAD models provide a single value for each of the fit parameters, including mean effective axon diameter AxD¯, packing fraction fin, intra-cellular and Din and extra-cellular Dex diffusion coefficients, as well as the frequency dependence of Dex, βex for the AAD model. The ACI and ACD models provide this, in addition to a distribution of axon diameters for a chosen ROI. VBA extends this, providing a parameter value for each voxel within the selected ROI, at the cost of increased computational load and analysis time. Overall, RBA-ACD and VBA-AAD were found to be optimal for parameter fitting to physically relevant values in a reasonable time frame. A full comparison of each combination of RBA and VBA with AAI, AAD, ACI, and ACD is provided to give the reader sufficient information to make an informed decision of which model is best for their own experiments.

Tracing volitional recovery in post-stroke akinetic mutism using longitudinal microstructure imaging: Insights from a single case study

Lesions in the frontal-subcortical circuitry can lead to akinetic mutism (AM) characterized by diminished volition. However, the microstructural changes in the damaged network underlying its recovery remain unknown. Clinical examination and neuropsychological assessment were performed on a patient with post-stroke AM. Multimodal MRI scans were performed at baseline and follow-ups. We used diffusion MRI and biophysical models, specifically utilizing neurite orientation dispersion and density imaging for assessing gray matter microstructure, and fixel-based analysis for the evaluation of white matter. Longitudinal comparisons were performed between the patient and healthy controls. Pronounced recovery of volition was observed after dopamine agonist therapy combined with physical therapy. In addition to infarcts in the bilateral medial cortex, microstructure imaging detected reduced neurite density in extensive areas, specifically in temporal areas and subcortical nuclei, and decreased fiber density of white matter tracts (TFCE-corrected p < .05). Microstructural degeneration in the anterior cingulate cortex and cingulum was relatively persistent (Bonferroni-corrected p < .05). However, most tracts within the frontal-subcortical circuitry showed increased fiber density during the recovery stage. Microstructure of an extensive network may contribute to the disruption and recovery of volition. Fiber density within the frontal-subcortical circuitry could be a promising biomarker indicating volitional recovery.

Advanced diffusion MRI provides evidence for altered axonal microstructure and gradual peritumoral infiltration in GBM in comparison to brain metastases

PURPOSE

In contrast to peritumoral edema in metastases, GBM is histopathologically characterized by infiltrating tumor cells within the T2 signal alterations. We hypothesized that depending on the distance from the outline of the contrast-enhancing tumor we might reveal imaging evidence of gradual peritumoral infiltration in GBM and predominantly vasogenic edema around metastases. We thus investigated the gradual change of advanced diffusion metrics with the peritumoral zone in metastases and GBM.

METHODS

In 30 patients with GBM and 28 with brain metastases, peritumoral T2 hyperintensity was segmented in 33% partitions based on the total volume beginning at the enhancing tumor margin and divided into inner, middle and outer zones. Diffusion Tensor Imaging (DTI)-derived fractional anisotropy and mean diffusivity as well as Diffusion Microstructure Imaging (DMI)-based parameters Dax-intra, Dax-extra, V‑CSF and V-intra were employed to assess group-wise differences between inner and outer zones as well as within-group gradients between the inner and outer zones.

RESULTS

In metastases, fractional anisotropy and Dax-extra were significantly reduced in the inner zone compared to the outer zone (FA p = 0.01; Dax-extra p = 0.03). In GBM, we noted a reduced Dax-extra and significantly lower intraaxonal volume fraction (Dax-extra p = 0.008, V‑intra p = 0.006) accompanied by elevated axial intraaxonal diffusivity in the inner zone (p = 0.035). Between-group comparison of the outer to the inner zones revealed significantly higher gradients in metastases over GBM for FA (p = 0.04) as well as the axial diffusivity in the intra- (p = 0.02) and extraaxonal compartment (p < 0.001).

CONCLUSION

Our findings provide evidence of gradual alterations within the peritumoral zone of brain tumors. These are compatible with predominant (vasogenic) edema formation in metastases, whereas our findings in GBM are in line with an axonal destructive component in the immediate peritumoral area and evidence of tumor cell infiltration with accentuation in the tumor's vicinity.

Predictors of cognition after glioma surgery: connectotomy, structure-function phenotype, plasticity

Determining preoperatively the maximal extent of resection that would preserve cognitive functions is the core challenge of brain tumour surgery. Over the past decade, the methodological framework to achieve this goal has been thoroughly renewed: the population-level topographically-focused voxel-based lesion-symptom mapping has been progressively overshadowed by machine learning (ML) algorithmics, in which the problem is framed as predicting cognitive outcomes in a patient-specific manner from a typically large set of variables. However, the choice of these predictors is of utmost importance, as they should be both informative and parsimonious. In this perspective, we first introduce the concept of connectotomy: instead of parameterizing resection topography through the status (intact/resected) of a huge number of voxels (or parcels) paving the whole brain in the Cartesian 3D-space, the connectotomy models the resection in the connectivity space, by computing a handful number of networks disconnection indices, measuring how the structural connectivity sustaining each network of interest was hit by the resection. This connectivity-informed reduction of dimensionality is a necessary step for efficiently implementing ML tools, given the relatively small number of patient-examples in available training datasets. We further argue that two other major sources of interindividual variability must be considered to improve the accuracy with which outcomes are predicted: the underlying structure-function phenotype and neuroplasticity, for which we provide an in-depth review and propose new ways of determining relevant predictors. We finally discuss the benefits of our approach for precision surgery of glioma.

Updates in diagnostic tools for diagnosing nerve injury and compressions

Predicting prognosis after nerve injury and compression can be challenging, even for the experienced clinician. Although thorough clinical assessment can aid diagnosis, we cannot always be precise about long-term functional recovery of either motor or sensory nerves. To evaluate the severity of nerve injury, surgical exploration remains the gold standard, particularly after iatrogenic injury and major nerve injury from trauma, such as brachial plexus injury. Recently, advances in imaging techniques (ultrasound, magnetic resonance imaging [MRI] and MR neurography) along with multimodality assessment, including electrodiagnostic testing, have allowed us to have a better preoperative understanding of nerve continuity and prediction of nerve health and possible recovery. This article outlines the current and potential roles for clinical assessment, exploratory surgery, electrodiagnostic testing ultrasound and MRI in entrapment neuropathies, inflammatory neuritis and trauma. Emphasis is placed on those modalities that are improving in diagnostic accuracy of nerve assessment before any surgical intervention.

The effects of axonal beading and undulation on axonal diameter estimation from diffusion MRI: Insights from simulations in human axons segmented from three-dimensional electron microscopy

The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings b , where the deviation from the expected 1 / b scaling in white matter yields a finite transverse diffusivity, which is then translated into an axon diameter estimate. While axons are usually modeled as perfectly straight, impermeable cylinders, local variations in diameter (caliber variation or beading) and direction (undulation) are known to influence axonal diameter estimates and have been observed in microscopy data of human axons. In this study, we performed Monte Carlo simulations of diffusion in axons reconstructed from three-dimensional electron microscopy of a human temporal lobe specimen using simulated sequence parameters matched to the maximal gradient strength of the next-generation Connectome 2.0 human MRI scanner ( ≲ 500 mT/m). We show that axon diameter estimation is accurate for nonbeaded, nonundulating fibers; however, in fibers with caliber variations and undulations, the axon diameter is heavily underestimated due to caliber variations, and this effect overshadows the known overestimation of the axon diameter due to undulations. This unexpected underestimation may originate from variations in the coarse-grained axial diffusivity due to caliber variations. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Diffusion kurtosis MRI tracks gray matter myelin content in the primate cerebral cortex

Diffusion magnetic resonance imaging (dMRI) has been widely employed to model the trajectory of myelinated fiber bundles in white matter. Increasingly, dMRI is also used to assess local tissue properties throughout the brain. In the cerebral cortex, myelin content is a critical indicator of the maturation, regional variation, and disease related degeneration of gray matter tissue. Gray matter myelination can be measured and mapped using several non-diffusion MRI strategies; however, first order diffusion statistics such as fractional anisotropy (FA) show only weak spatial correlation with cortical myelin content. Here we show that a simple higher order diffusion parameter, the mean diffusion kurtosis (MK), is strongly correlated with the laminar and regional variation of myelin in the primate cerebral cortex. We carried out ultra-high resolution, multi-shelled dMRI in ex vivo marmoset monkey brains and compared dMRI parameters from a number of higher order models (diffusion kurtosis, NODDI and MAP MRI) to the distribution of myelin obtained using histological staining, and via Magnetization Transfer Ratio MRI (MTR), a non-diffusion MRI method. In contrast to FA, MK closely matched the myelin content assessed by histology and by MTR in the same sample. The parameter maps from MAP-MRI and NODDI also showed good correspondence with cortical myelin content. The results demonstrate that dMRI can be used to assess the variation of local myelin content in the primate cortical cortex, which may be of great value for assessing tissue integrity and tracking disease in living human patients.

Advanced diffusion imaging reveals microstructural characteristics of primary CNS lymphoma, allowing differentiation from glioblastoma

Background

Primary CNS lymphoma (PCNSL) and glioblastoma (GBM) both represent frequent intracranial malignancies with differing clinical management. However, distinguishing PCNSL from GBM with conventional MRI can be challenging when atypical imaging features are present. We employed advanced dMRI for noninvasive characterization of the microstructure of PCNSL and differentiation from GBM as the most frequent primary brain malignancy.

Methods

Multiple dMRI metrics including Diffusion Tensor Imaging, Neurite Orientation Dispersion and Density Imaging, and Diffusion Microstructure Imaging were extracted from the contrast-enhancing tumor component in 10 PCNSL and 10 age-matched GBM on 3T MRI. Imaging findings were correlated with cell density and axonal markers obtained from histopathology.

Results

We found significantly increased intra-axonal volume fractions (V-intra and intracellular volume fraction) and microFA in PCNSL compared to GBM (all P < .001). In contrast, mean diffusivity (MD), axial diffusivity (aD), and microADC (all P < .001), and also free water fractions (V-CSF and V-ISO) were significantly lower in PCNSL (all P < .01). Receiver-operating characteristic analysis revealed high predictive values regarding the presence of a PCNSL for MD, aD, microADC, V-intra, ICVF, microFA, V-CSF, and V-ISO (area under the curve [AUC] in all >0.840, highest for MD and ICVF with an AUC of 0.960). Comparative histopathology between PCNSL and GBM revealed a significantly increased cell density in PCNSL and the presence of axonal remnants in a higher proportion of samples.

Conclusions

Advanced diffusion imaging enables the characterization of the microstructure of PCNSL and reliably distinguishes PCNSL from GBM. Both imaging and histopathology revealed a relatively increased cell density and a preserved axonal microstructure in PCNSL.

Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain

The dependence of the diffusion MRI signal on the diffusion time carries signatures of restricted diffusion and exchange. Here we seek to highlight these signatures in the human brain by performing experiments using free gradient waveforms designed to be selectively sensitive to the two effects. We examine six healthy volunteers using both strong and ultra-strong gradients (80, 200 and 300 mT/m). In an experiment featuring a large set of 150 gradient waveforms with different sensitivities to restricted diffusion and exchange, our results reveal unique and different time-dependence signatures in grey and white matter. Grey matter was characterised by both restricted diffusion and exchange and white matter predominantly by restricted diffusion. Exchange in grey matter was at least twice as fast as in white matter, across all subjects and all gradient strengths. The cerebellar cortex featured relatively short exchange times (115 ms). Furthermore, we show that gradient waveforms with tailored designs can be used to map exchange in the human brain. We also assessed the feasibility of clinical applications of the method used in this work and found that the exchange-related contrast obtained with a 25-minute protocol at 300 mT/m was preserved in a 4-minute protocol at 300 mT/m and a 10-minute protocol at 80 mT/m. Our work underlines the utility of free waveforms for detecting time dependence signatures due to restricted diffusion and exchange in vivo, which may potentially serve as a tool for studying diseased tissue.

Frequency-dependent diffusion kurtosis imaging in the human brain using an oscillating gradient spin echo sequence and a high-performance head-only gradient

Measuring the time/frequency dependence of diffusion MRI is a promising approach to distinguish between the effects of different tissue microenvironments, such as membrane restriction, tissue heterogeneity, and compartmental water exchange. In this study, we measure the frequency dependence of diffusivity (D) and kurtosis (K) with oscillating gradient diffusion encoding waveforms and a diffusion kurtosis imaging (DKI) model in human brains using a high-performance, head-only MAGNUS gradient system, with a combination of b-values, oscillating frequencies (f), and echo time that has not been achieved in human studies before. Frequency dependence of diffusivity and kurtosis are observed in both global and local white matter (WM) and gray matter (GM) regions and characterized with a power-law model ∼Λ*fθ. The frequency dependences of diffusivity and kurtosis (including changes between fmin and fmax, Λ, and θ) vary over different WM and GM regions, indicating potential microstructural differences between regions. A trend of decreasing kurtosis over frequency in the short-time limit is successfully captured for in vivo human brains. The effects of gradient nonlinearity (GNL) on frequency-dependent diffusivity and kurtosis measurements are investigated and corrected. Our results show that the GNL has prominent scaling effects on the measured diffusivity values (3.5∼5.5% difference in the global WM and 6∼8% difference in the global cortex) and subsequently affects the corresponding power-law parameters (Λ, θ) while having a marginal influence on the measured kurtosis values (<0.05% difference) and power-law parameters (Λ, θ). This study expands previous OGSE studies and further demonstrates the translatability of frequency-dependent diffusivity and kurtosis measurements to human brains, which may provide new opportunities to probe human brain microstructure in health and disease.

The influence of axonal beading and undulation on axonal diameter mapping

We consider the effect of non-cylindrical axonal shape on axonal diameter mapping with diffusion MRI. Practical sensitivity to axon diameter is attained at strong diffusion weightings b , where the deviation from the 1 / b scaling yields the finite transverse diffusivity, which is then translated into axon diameter. While axons are usually modeled as perfectly straight, impermeable cylinders, the local variations in diameter (caliber variation or beading) and direction (undulation) have been observed in microscopy data of human axons. Here we quantify the influence of cellular-level features such as caliber variation and undulation on axon diameter estimation. For that, we simulate the diffusion MRI signal in realistic axons segmented from 3-dimensional electron microscopy of a human brain sample. We then create artificial fibers with the same features and tune the amplitude of their caliber variations and undulations. Numerical simulations of diffusion in fibers with such tunable features show that caliber variations and undulations result in under- and over-estimation of axon diameters, correspondingly; this bias can be as large as 100%. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain

The dependence of the diffusion MRI signal on the diffusion time carries signatures of restricted diffusion and exchange. Here we seek to highlight these signatures in the human brain by performing experiments using free gradient waveforms that are selectively sensitive to the two effects. We examine six healthy volunteers using both strong and ultra-strong gradients (80, 200 and 300 mT/m). In an experiment featuring a large set of gradient waveforms with different sensitivities to restricted diffusion and exchange (150 samples), our results reveal unique time-dependence signatures in grey and white matter, where the former is characterised by both restricted diffusion and exchange and the latter predominantly exhibits restricted diffusion. Furthermore, we show that gradient waveforms with independently varying sensitivities to restricted diffusion and exchange can be used to map exchange in the human brain. We consistently find that exchange in grey matter is at least twice as fast as in white matter, across all subjects and all gradient strengths. The shortest exchange times observed in this study were in the cerebellar cortex (115 ms). We also assess the feasibility of future clinical applications of the method used in this work, where we find that the grey-white matter exchange contrast obtained with a 25-minute 300 mT/m protocol is preserved by a 4-minute 300 mT/m and a 10-minute 80 mT/m protocol. Our work underlines the utility of free waveforms for detecting time-dependence signatures due to restricted diffusion and exchange in vivo, which may potentially serve as a tool for studying diseased tissue.

Assessment of complementary white matter microstructural changes and grey matter atrophy in a preclinical model of Alzheimer's disease

Alzheimer's disease (AD) has been associated with amyloid and tau pathology, as well as neurodegeneration. Beyond these hallmark features, white matter microstructural abnormalities have been observed using MRI. The objective of this study was to assess grey matter atrophy and white matter microstructural changes in a preclinical mouse model of AD (3xTg-AD) using voxel-based morphometry (VBM) and free-water (FW) diffusion tensor imaging (FW-DTI). Compared to controls, lower grey matter density was observed in the 3xTg-AD model, corresponding to the small clusters in the caudate-putamen, hypothalamus, and cortex. DTI-based fractional anisotropy (FA) was decreased in the 3xTg model, while the FW index was increased. Notably, the largest clusters for both FW-FA and FW index were in the fimbria, with other regions including the anterior commissure, corpus callosum, forebrain septum, and internal capsule. Additionally, the presence of amyloid and tau in the 3xTg model was confirmed with histopathology, with significantly higher levels observed across many regions of the brain. Taken together, these results are consistent with subtle neurodegenerative and white matter microstructural changes in the 3xTg-AD model that manifest as increased FW, decreased FW-FA, and decreased grey matter density.

Linear multi-scale modeling of diffusion MRI data: A framework for characterization of oriented structures across length scales

Diffusion-weighted magnetic resonance imaging (DW-MRI) has evolved to provide increasingly sophisticated investigations of the human brain's structural connectome in vivo. Restriction spectrum imaging (RSI) is a method that reconstructs the orientation distribution of diffusion within tissues over a range of length scales. In its original formulation, RSI represented the signal as consisting of a spectrum of Gaussian diffusion response functions. Recent technological advances have enabled the use of ultra-high b-values on human MRI scanners, providing higher sensitivity to intracellular water diffusion in the living human brain. To capture the complex diffusion time dependence of the signal within restricted water compartments, we expand upon the RSI approach to represent restricted water compartments with non-Gaussian response functions, in an extended analysis framework called linear multi-scale modeling (LMM). The LMM approach is designed to resolve length scale and orientation-specific information with greater specificity to tissue microstructure in the restricted and hindered compartments, while retaining the advantages of the RSI approach in its implementation as a linear inverse problem. Using multi-shell, multi-diffusion time DW-MRI data acquired with a state-of-the-art 3 T MRI scanner equipped with 300 mT/m gradients, we demonstrate the ability of the LMM approach to distinguish different anatomical structures in the human brain and the potential to advance mapping of the human connectome through joint estimation of the fiber orientation distributions and compartment size characteristics.

Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange

Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging anisotropic compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along neuronal and glial processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Kärger model of two exchanging compartments. Using datasets acquired at multiple diffusion weightings (b) and diffusion times (t) in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the gray matter, compared to the other two possible contributions. Our results for the diffusion time window 20-45 ms show minimal diffusivity time-dependence and more pronounced kurtosis decay with time, which is well fit by the exchange model. Moreover, we observe lower signal for longer diffusion times at high b. In light of these observations, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for gray matter microstructure mapping. We finally highlight multi-b multi-t acquisition protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 15 - 60 ms in the rat cortex and hippocampus in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions. This suggests water exchange as an essential component for interpreting diffusion MRI measurements in gray matter.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

Diffusion time dependence, power-law scaling, and exchange in gray matter

Characterizing neural tissue microstructure is a critical goal for future neuroimaging. Diffusion MRI (dMRI) provides contrasts that reflect diffusing spins' interactions with myriad microstructural features of biological systems. However, the specificity of dMRI remains limited due to the ambiguity of its signals vis-à-vis the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) taken spherical compartments into account (the SANDI model) in attempts to represent cell soma. The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. To validate these assumptions experimentally, observing their unique, functional properties, such as the b-1/2 power-law associated with one-dimensional diffusion, has emerged as a fruitful strategy. The absence of this signature in GM, in turn, has been explained by neurite water exchange, non-linear morphology, and/or by obscuring soma signal contributions. Here, we present diffusion simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law behavior in impermeable neurites, but also that their signal is drowned by the soma signal under typical experimental conditions. Nevertheless, by studying the GM dMRI signal's behavior as a function of diffusion weighting as well as time, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a signal behavior opposite to that which would be expected from restricted diffusion, thereby providing a functional signature that disambiguates the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field of 16.4T and observe a time dependence that is consistent with substantial exchange but also with a GM stick power-law. The first finding suggests significant water exchange between neurites and the extracellular space while the second suggests a small sub-population of impermeable neurites. To quantify these observations, we harness the Kärger exchange model and incorporate the corresponding signal time dependence in the SM and SANDI models.

In vivo Correlation Tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner

Diffusion MRI (dMRI) has become one of the most important imaging modalities for noninvasively probing tissue microstructure. Diffusional Kurtosis MRI (DKI) quantifies the degree of non-gaussian diffusion, which in turn has been shown to increase sensitivity towards, e.g., disease and orientation mapping in neural tissue. However, the specificity of DKI is limited as different sources can contribute to the total intravoxel diffusional kurtosis, including: variance in diffusion tensor magnitudes (K_iso), variance due to diffusion anisotropy (K_aniso), and microscopic kurtosis (μK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, μK is typically ignored in diffusion MRI signal modeling as it is assumed to be negligible in neural tissues. However, recently, Correlation Tensor MRI (CTI) based on Double-Diffusion-Encoding (DDE) was introduced for kurtosis source separation, revealing non negligible μK in preclinical imaging. Here, we implemented CTI for the first time on a clinical 3T scanner and investigated the sources of total kurtosis in healthy subjects. A robust framework for kurtosis source separation in humans is introduced, followed by estimation of μK (and the other kurtosis sources) in the healthy brain. Using this clinical CTI approach, we find that μK significantly contributes to total diffusional kurtosis both in gray and white matter tissue but, as expected, not in the ventricles. The first μK maps of the human brain are presented, revealing that the spatial distribution of μK provides a unique source of contrast, appearing different from isotropic and anisotropic kurtosis counterparts. Moreover, group average templates of these kurtosis sources have been generated for the first time, which corroborated our findings at the underlying individual-level maps. We further show that the common practice of ignoring μK and assuming the multiple gaussian component approximation for kurtosis source estimation introduces significant bias in the estimation of other kurtosis sources and, perhaps even worse, compromises their interpretation. Finally, a twofold acceleration of CTI is discussed in the context of potential future clinical applications. We conclude that CTI has much potential for future in vivo microstructural characterizations in healthy and pathological tissue.

Anomalous diffusion, aging, and nonergodicity of scaled Brownian motion with fractional Gaussian noise: overview of related experimental observations and models

How does a systematic time-dependence of the diffusion coefficient $D (t)$ affect the ergodic and statistical characteristics of fractional Brownian motion (FBM)? Here, we examine how the behavior of the...

Macroscopic Structural and Connectome Mapping of the Mouse Brain Using Diffusion Magnetic Resonance Imaging

Translational work in rodents elucidates basic mechanisms that drive complex behaviors relevant to psychiatric and neurological conditions. Nonetheless, numerous promising studies in rodents later fail in clinical trials, highlighting the need for improving the translational utility of preclinical studies in rodents. Imaging of small rodents provides an important strategy to address this challenge, as it enables a whole-brain unbiased search for structural and dynamic changes that can be directly compared to human imaging. The functional significance of structural changes identified using imaging can then be further investigated using molecular and genetic tools available for the mouse. Here, we describe a pipeline for unbiased search and characterization of structural changes and network properties, based on diffusion MRI data covering the entire mouse brain at an isotropic resolution of 100 µm. We first used unbiased whole-brain voxel-based analyses to identify volumetric and microstructural alterations in the brain of adult mice exposed to unpredictable postnatal stress (UPS), which is a mouse model of complex early life stress (ELS). Brain regions showing structural abnormalities were used as nodes to generate a grid for assessing structural connectivity and network properties based on graph theory. The technique described here can be broadly applied to understand brain connectivity in other mouse models of human disorders, as well as in genetically modified mouse strains. Graphic abstract: Pipeline for characterizing structural connectome in the mouse brain using diffusion magnetic resonance imaging. Scale bar = 1 mm.

Test-retest reproducibility of in vivo oscillating gradient and microscopic anisotropy diffusion MRI in mice at 9.4 Tesla

BACKGROUND AND PURPOSE

Microstructure imaging with advanced diffusion MRI (dMRI) techniques have shown increased sensitivity and specificity to microstructural changes in various disease and injury models. Oscillating gradient spin echo (OGSE) dMRI, implemented by varying the oscillating gradient frequency, and microscopic anisotropy (μA) dMRI, implemented via tensor valued diffusion encoding, may provide additional insight by increasing sensitivity to smaller spatial scales and disentangling fiber orientation dispersion from true microstructural changes, respectively. The aims of this study were to characterize the test-retest reproducibility of in vivo OGSE and μA dMRI metrics in the mouse brain at 9.4 Tesla and provide estimates of required sample sizes for future investigations.

METHODS

Twelve adult C57Bl/6 mice were scanned twice (5 days apart). Each imaging session consisted of multifrequency OGSE and μA dMRI protocols. Metrics investigated included μA, linear diffusion kurtosis, isotropic diffusion kurtosis, and the diffusion dispersion rate (Λ), which explores the power-law frequency dependence of mean diffusivity. The dMRI metric maps were analyzed with mean region-of-interest (ROI) and whole brain voxel-wise analysis. Bland-Altman plots and coefficients of variation (CV) were used to assess the reproducibility of OGSE and μA metrics. Furthermore, we estimated sample sizes required to detect a variety of effect sizes.

RESULTS

Bland-Altman plots showed negligible biases between test and retest sessions. ROI-based CVs revealed high reproducibility for most metrics (CVs < 15%). Voxel-wise CV maps revealed high reproducibility for μA (CVs ~ 10%), but low reproducibility for OGSE metrics (CVs ~ 50%).

CONCLUSION

Most of the μA dMRI metrics are reproducible in both ROI-based and voxel-wise analysis, while the OGSE dMRI metrics are only reproducible in ROI-based analysis. Given feasible sample sizes (10-15), μA metrics and OGSE metrics may provide sensitivity to subtle microstructural changes (4-8%) and moderate changes (> 6%), respectively.

Diffusion MRI tractography for neurosurgery: the basics, current state, technical reliability and challenges

Diffusion magnetic resonance imaging (dMRI) tractography is currently the only imaging technique that allows for non-invasive delineation and visualisation of white matter (WM) tractsin vivo,prompting rapid advances in related fields of brain MRI research in recent years. One of its major clinical applications is for pre-surgical planning and intraoperative image guidance in neurosurgery, where knowledge about the location of WM tracts nearby the surgical target can be helpful to guide surgical resection and optimise post-surgical outcomes. Surgical injuries to these WM tracts can lead to permanent neurological and functional deficits, making the accuracy of tractography reconstructions paramount. The quality of dMRI tractography is influenced by many modifiable factors, ranging from MRI data acquisition through to the post-processing of tractography output, with the potential of error propagation based on decisions made at each and subsequent processing steps. Research over the last 25 years has significantly improved the anatomical accuracy of tractography. An updated review about tractography methodology in the context of neurosurgery is now timely given the thriving research activities in dMRI, to ensure more appropriate applications in the clinical neurosurgical realm. This article aims to review the dMRI physics, and tractography methodologies, highlighting recent advances to provide the key concepts of tractography-informed neurosurgery, with a focus on the general considerations, the current state of practice, technical challenges, potential advances, and future demands to this field.

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.