Microstructural organization of axons in the human corpus callosum quantified by diffusion-weighted magnetic resonance spectroscopy of N-acetylaspartate and post-mortem histology

Time to bury the chisel: a continuous dorsal association tract system

The arcuate fasciculus may be subdivided into a tract directly connecting frontal and temporal lobes and a pair of indirect subtracts in which the fronto-temporal connection is mediated by connections to the inferior parietal lobe. This tripartition has been advanced as an improvement over the centuries-old consensus that the lateral dorsal association fibers form a continuous system with no discernible discrete parts. Moreover, it has been used as the anatomical basis for functional hypotheses regarding linguistic abilities. Ex hypothesi, damage to the indirect subtracts leads to deficits in the repetition of multi-word sequences, whereas damage to the direct subtract leads to deficits in the immediate reproduction of single multisyllabic words. We argue that this partitioning of the dorsal association tract system enjoys no special anatomical status, and the search for the anatomical substrates of linguistic abilities should not be constrained by it. Instead, the merit of any postulated partitioning should primarily be judged on the basis of whether it enlightens or obfuscates our understanding of the behavior of patients in which individual subtracts are damaged.

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-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

A translational MRI approach to validate acute axonal damage detection as an early event in multiple sclerosis

Axonal degeneration is a central pathological feature of multiple sclerosis and is closely associated with irreversible clinical disability. Current noninvasive methods to detect axonal damage in vivo are limited in their specificity and clinical applicability, and by the lack of proper validation. We aimed to validate an MRI framework based on multicompartment modeling of the diffusion signal (AxCaliber) in rats in the presence of axonal pathology, achieved through injection of a neurotoxin damaging the neuronal terminal of axons. We then applied the same MRI protocol to map axonal integrity in the brain of multiple sclerosis relapsing-remitting patients and age-matched healthy controls. AxCaliber is sensitive to acute axonal damage in rats, as demonstrated by a significant increase in the mean axonal caliber along the targeted tract, which correlated with neurofilament staining. Electron microscopy confirmed that increased mean axonal diameter is associated with acute axonal pathology. In humans with multiple sclerosis, we uncovered a diffuse increase in mean axonal caliber in most areas of the normal-appearing white matter, preferentially affecting patients with short disease duration. Our results demonstrate that MRI-based axonal diameter mapping is a sensitive and specific imaging biomarker that links noninvasive imaging contrasts with the underlying biological substrate, uncovering generalized axonal damage in multiple sclerosis as an early event.

Diffusion-weighted SPECIAL improves the detection of J-coupled metabolites at ultrahigh magnetic field

PURPOSE

To improve the detection and subsequent estimation of the diffusion properties of strongly J-coupled metabolites in diffusion-weighted MRS (DWS).

METHODS

A new sequence for single-voxel diffusion-weighted 1 H MR spectroscopy, named DW-SPECIAL, is proposed. It combines the semi-adiabatic SPECIAL sequence with a stimulated echo diffusion block. Acquisitions with DW-SPECIAL and STE-LASER, the current gold standard for rodent DWS experiments at high fields, were performed at 14.1T on phantoms and in vivo on the rat brain. The apparent diffusion coefficient and intra-stick diffusivity (Callaghan's model, randomly-oriented sticks) were fitted and compared between the sequences for glutamate, glutamine, myo-inositol, taurine, total NAA, total Cho, total Cr, and the macromolecules.

RESULTS

The shorter TE achieved with DW-SPECIAL (18 ms against 33 ms with STE-LASER) substantially limited the metabolites' signal loss caused by J-evolution. In addition, DW-SPECIAL preserved the main advantages of STE-LASER: absence of cross-terms, diffusion time during a stimulated echo, and limited sensitivity to B1 inhomogeneities. In vivo, compared to STE-LASER, DW-SPECIAL yielded the same spectral quality and reduced the Cramer Rao Lower Bounds for J-coupled metabolites, irrespective of the b-value. DW-SPECIAL also reduced the SD of the metabolites' diffusion estimates based on individual animal fitting without loss of accuracy compared to the fit on the averaged decay.

CONCLUSION

We conclude that due to its reduced TE, DW-SPECIAL can serve as an alternative to STE-LASER when strongly J-coupled metabolites like glutamine are investigated, thereby extending the range of accessible metabolites in the context of DWS acquisitions.

Practical considerations of diffusion-weighted MRS with ultra-strong diffusion gradients

Introduction

Diffusion-weighted magnetic resonance spectroscopy (DW-MRS) offers improved cellular specificity to microstructure-compared to water-based methods alone-but spatial resolution and SNR is severely reduced and slow-diffusing metabolites necessitate higher b-values to accurately characterize their diffusion properties. Ultra-strong gradients allow access to higher b-values per-unit time, higher SNR for a given b-value, and shorter diffusion times, but introduce additional challenges such as eddy-current artefacts, gradient non-uniformity, and mechanical vibrations.

Methods

In this work, we present initial DW-MRS data acquired on a 3T Siemens Connectom scanner equipped with ultra-strong (300 mT/m) gradients. We explore the practical issues associated with this manner of acquisition, the steps that may be taken to mitigate their impact on the data, and the potential benefits of ultra-strong gradients for DW-MRS. An in-house DW-PRESS sequence and data processing pipeline were developed to mitigate the impact of these confounds. The interaction of TE, b-value, and maximum gradient amplitude was investigated using simulations and pilot data, whereby maximum gradient amplitude was restricted. Furthermore, two DW-MRS voxels in grey and white matter were acquired using ultra-strong gradients and high b-values.

Results

Simulations suggest T2-based SNR gains that are experimentally confirmed. Ultra-strong gradient acquisitions exhibit similar artefact profiles to those of lower gradient amplitude, suggesting adequate performance of artefact mitigation strategies. Gradient field non-uniformity influenced ADC estimates by up to 4% when left uncorrected. ADC and Kurtosis estimates for tNAA, tCho, and tCr align with previously published literature.

Discussion

In conclusion, we successfully implemented acquisition and data processing strategies for ultra-strong gradient DW-MRS and results indicate that confounding effects of the strong gradient system can be ameliorated, while achieving shorter diffusion times and improved metabolite SNR.

Histology-informed multiscale modeling of human brain white matter

In this study, we propose a novel micromechanical model for the brain white matter, which is described as a heterogeneous material with a complex network of axon fibers embedded in a soft ground matrix. We developed this model in the framework of RVE-based multiscale theories in combination with the finite element method and the embedded element technique for embedding the fibers. Microstructural features such as axon diameter, orientation and tortuosity are incorporated into the model through distributions derived from histological data. The constitutive law of both the fibers and the matrix is described by isotropic one-term Ogden functions. The hyperelastic response of the tissue is derived by homogenizing the microscopic stress fields with multiscale boundary conditions to ensure kinematic compatibility. The macroscale homogenized stress is employed in an inverse parameter identification procedure to determine the hyperelastic constants of axons and ground matrix, based on experiments on human corpus callosum. Our results demonstrate the fundamental effect of axon tortuosity on the mechanical behavior of the brain's white matter. By combining histological information with the multiscale theory, the proposed framework can substantially contribute to the understanding of mechanotransduction phenomena, shed light on the biomechanics of a healthy brain, and potentially provide insights into neurodegenerative processes.

Multiscale label-free imaging of myelin in human brain tissue with polarization-sensitive optical coherence tomography and birefringence microscopy

The combination of polarization-sensitive optical coherence tomography (PS-OCT) and birefringence microscopy (BRM) enables multiscale assessment of myelinated axons in postmortem brain tissue, and these tools are promising for the study of brain connectivity and organization. We demonstrate label-free imaging of myelin structure across the mesoscopic and microscopic spatial scales by performing serial-sectioning PS-OCT of a block of human brain tissue and periodically sampling thin sections for high-resolution imaging with BRM. In co-registered birefringence parameter maps, we observe good correspondence and demonstrate that BRM enables detailed validation of myelin (hence, axonal) organization, thus complementing the volumetric information content of PS-OCT.

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.

Larmor frequency shift from magnetized cylinders with arbitrary orientation distribution

The magnetic susceptibility of tissue can provide valuable information about its chemical composition and microstructural organization. However, the relation between the magnetic microstructure and the measurable Larmor frequency shift is understood only for a few idealized cases. Here we analyze the microstructure formed by magnetized, NMR-invisible infinite cylinders suspended in an NMR-reporting fluid. Through simulations, we scrutinize various geometries of mesoscopic Lorentz cavities and inclusions, and show that the cavity size should be approximately one order of magnitude larger than the width of the inclusions. We also analytically derive the Larmor frequency shift for a population of cylinders with arbitrary orientation dispersion and show that it is determined by the l = 2 Laplace expansion coefficients p 2 m of the cylinders' orientation distribution function. Our work underscores the need to account for microstructural organization when estimating magnetic tissue properties.

MP-PCA denoising for diffusion MRS data: promises and pitfalls

Diffusion-weighted (DW) magnetic resonance spectroscopy (MRS) suffers from a lower signal to noise ratio (SNR) compared to conventional MRS owing to the addition of diffusion attenuation. This technique can therefore strongly benefit from noise reduction strategies. In the present work, Marchenko-Pastur principal component analysis (MP-PCA) denoising is tested on Monte Carlo simulations and on in vivo DW-MRS data acquired at 9.4 T in rat brain and at 3 T in human brain. We provide a descriptive study of the effects observed following different MP-PCA denoising strategies (denoising the entire matrix versus using a sliding window), in terms of apparent SNR, rank selection, noise correlation within and across b-values and quantification of metabolite concentrations and fitted diffusion coefficients. MP-PCA denoising yielded an increased apparent SNR, a more accurate B0 drift correction between shots, and similar estimates of metabolite concentrations and diffusivities compared to the raw data. No spectral residuals on individual shots were observed but correlations in the noise level across shells were introduced, an effect which was mitigated using a sliding window, but which should be carefully considered.

White Matter Metabolite Relaxation and Diffusion Abnormalities in First-Episode Psychosis: A Longitudinal Study

Microstructural abnormalities in the white matter (WM) are implicated in the pathophysiology of psychosis. In vivo magnetic resonance spectroscopy (MRS) can probe the brain's intracellular microenvironment through the measurement of transverse relaxation and diffusion of neurometabolites and possibly provide cell-specific information. In our previous studies, we observed differential metabolite signal abnormalities in first episode and chronic stages of psychosis. In the present work, longitudinal data were presented for the first time on white matter cell-type specific abnormalities using a combination of diffusion tensor spectroscopy (DTS), T2 MRS, and diffusion tensor imaging (DTI) from a group of 25 first episode psychosis patients and nine matched controls scanned at baseline and one and two years of follow-up. We observed significantly reduced choline ADC in the year 1 of follow-up (0.194 µm2/ms) compared to baseline (0.229 µm2/ms), followed by a significant increase in NAA ADC in the year 2 follow-up (0.258 µm2/ms) from baseline (0.222 µm2/ms) and year 1 follow-up (0.217 µm2/ms). In contrast, NAA T2 relaxation, reflecting a related but different aspect of microenvironment from diffusion, was reduced at year 1 follow-up (257 ms) compared to baseline (278 ms). These abnormalities were observed in the absence of any abnormalities in water relaxation and diffusion at any timepoint. These findings indicate that abnormalities are seen in in glial-enriched (choline) signals in early stages of psychosis, followed by the subsequent emergence of neuronal-enriched (NAA) diffusion abnormalities, all in the absence of nonspecific water signal abnormalities.

High fidelity fiber orientation density functions from fiber ball imaging

The fiber orientation density function (fODF) in white matter is a primary physical quantity that can be estimated with diffusion MRI. It has often been employed for fiber tracking and microstructural modeling. Requirements for the construction of high fidelity fODFs, in the sense of having good angular resolution, adequate data to avoid sampling errors, and minimal noise artifacts, are described for fODFs calculated with fiber ball imaging. A criterion is formulated for the number of diffusion encoding directions needed to achieve a given angular resolution. The advantages of using large b-values (≥6000 s/mm² ) are also discussed. For the direct comparison of different fODFs, a method is developed for defining a local frame of reference tied to each voxel's individual axonal structure. The Matusita anisotropy axonal is proposed as a scalar fODF measure for quantifying angular variability. Experimental results, obtained at 3 T from human volunteers, are used as illustrations.

Does powder averaging remove dispersion bias in diffusion MRI diameter estimates within real 3D axonal architectures?

Noninvasive estimation of axon diameter with diffusion MRI holds the potential to investigate the dynamic properties of the brain network and pathology of neurodegenerative diseases. Recent studies use powder averaging to account for complex white matter architectures, but these have not been validated for real axonal geometries from regions that contain fibre crossings. Here, we present 120-304μm long segmented axons from X-ray nano-holotomography volumes of a splenium and crossing fibre region of a vervet monkey brain. We show that the axons in the complex crossing fibre region, which contains callosal, association, and corticospinal connections, are larger and exhibit a wider distribution than those of the splenium region. To accurately estimate the axon diameter in these regions, therefore, sensitivity to a wide range of diameters is required. We demonstrate how the q-value, b-value, signal-to-noise ratio and the assumed intra-axonal parallel diffusivity influence the range of measurable diameters with powder average approaches. Furthermore, we show how Gaussian distributed noise results in a wider range of measurable diameter at high b-values than Rician distributed noise, even at high signal-to-noise ratios of 100. The number of gradient directions is also shown to impose a lower bound on measurable diameter. Our results indicate that axon diameter estimation can be performed with only few b-shells, and that additional shells do not improve the accuracy of the estimate. For strong gradients available on human Connectom and preclinical scanners, Monte Carlo simulations of diffusion confirm that powder averaging techniques succeed in providing accurate estimates of axon diameter across a range of sequence parameters and diffusion times, even in complex white matter architectures. At relatively low b-values, the diameter estimate becomes sensitive to axonal microdispersion and the intra-axonal parallel diffusivity shows time dependency at both in vivo and ex vivo intrinsic diffusivities.

Cytosolic diffusivity and microscopic anisotropy of N-acetyl aspartate in human white matter with diffusion-weighted MRS at 7 T

Metabolite diffusion measurable in humans in vivo with diffusion-weighted spectroscopy (DW-MRS) provides a window into the intracellular morphology and state of specific cell types. Anisotropic diffusion in white matter is governed by the microscopic properties of the individual cell types and their structural units (axons, soma, dendrites). However, anisotropy is also markedly affected by the macroscopic orientational distribution over the imaging voxel, particularly in DW-MRS, where the dimensions of the volume of interest (VOI) are much larger than those typically used in diffusion-weighted imaging. One way to address the confound of macroscopic structural features is to average the measurements acquired with uniformly distributed gradient directions to mimic a situation where fibers present in the VOI are orientationally uniformly distributed. This situation allows the extraction of relevant microstructural features such as transverse and longitudinal diffusivities within axons and the related microscopic fractional anisotropy. We present human DW-MRS data acquired at 7 T in two different white matter regions, processed and analyzed as described above, and find that intra-axonal diffusion of the neuronal metabolite N-acetyl aspartate is in good correspondence to simple model interpretations, such as multi-Gaussian diffusion from disperse fibers where the transverse diffusivity can be neglected. We also discuss the implications of our approach for current and future applications of DW-MRS for cell-specific measurements.

Rapid acquisition diffusion MR spectroscopy of metabolites in human brain

Few studies have focused on metabolite diffusion in the human brain using ¹ H-MRS due to significant technical challenges. Moreover, such studies have required lengthy acquisition times and are therefore impractical to implement clinically. By first characterizing and then minimizing the effects of linear and oscillating eddy currents, which arise from the diffusion gradients, and by implementing phase-cycle and slice-order strategies, as well as introducing a new phase-alignment methodology, we report a method that allows data acquisition requiring 20 seconds per spectrum. This remained feasible, even for b-values >8000 s/mm² , with a rapid acquisition diffusion MRS methodology. It has allowed the nonlinear characterization of signal intensity with multiple b-values, and has improved the measurement of rotationally invariant diffusion parameters via six-direction, six b-value diffusion tensor spectroscopy (DTS) in 12 minutes at 4.7 T. The shorter DTS acquisition will enable its application to white matter regions not aligned with the gradients and permit clinical studies in a feasible time.

Compartmental diffusion and microstructural properties of human brain gray and white matter studied with double diffusion encoding magnetic resonance spectroscopy of metabolites and water

Double diffusion encoding (DDE) of the water signal offers a unique ability to separate the effect of microscopic anisotropic diffusion in structural units of tissue from the overall macroscopic orientational distribution of cells. However, the specificity in detected microscopic anisotropy is limited as the signal is averaged over different cell types and across tissue compartments. Performing side-by-side water and metabolite DDE spectroscopic (DDES) experiments provides complementary measures from which intracellular and extracellular microscopic fractional anisotropies (μFA) and diffusivities can be estimated. Metabolites are largely confined to the intracellular space and therefore provide a benchmark for intracellular μFA and diffusivities of specific cell types. By contrast, water DDES measurements allow examination of the separate contributions to water μFA and diffusivity from the intra- and extracellular spaces, by using a wide range of b values to gradually eliminate the extracellular contribution. Here, we aimed to estimate tissue and compartment specific human brain microstructure by combining water and metabolites DDES experiments. We performed our DDES measurements in two brain regions that contain widely different amounts of white matter (WM) and gray matter (GM): parietal white matter (PWM) and occipital gray matter (OGM) in a total of 20 healthy volunteers at 7 Tesla. Metabolite DDES measurements were performed at b = 7199 s/mm², while water DDES measurements were performed with a range of b values from 918 to 7199 s/mm². The experimental framework we employed here resulted in a set of insights pertaining to the morphology of the intracellular and extracellular spaces in both gray and white matter. Results of the metabolite DDES experiments in both PWM and OGM suggest a highly anisotropic intracellular space within neurons and glia, with the possible exception of gray matter glia. The water μFA obtained from the DDES results at high b values in both regions converged with that of the metabolite DDES, suggesting that the signal from the extracellular space is indeed effectively suppressed at the highest b value. The μFA measured in the OGM significantly decreased at lower b values, suggesting a considerably lower anisotropy of the extracellular space in GM compared to WM. In PWM, the water μFA remained high even at the lowest b value, indicating a high degree of organization in the interstitial space in WM. Tortuosity values in the cytoplasm for water and tNAA, obtained with correlation analysis of microscopic parallel diffusivity with respect to GM/WM tissue fraction in the volume of interest, are remarkably similar for both molecules, while exhibiting a clear difference between gray and white matter, suggesting a more crowded cytoplasm and more complex cytomorphology of neuronal cell bodies and dendrites in GM than those found in long-range axons in WM.

Probing tissue microstructure by diffusion skewness tensor imaging

Probing the cellular structure of in vivo biological tissue is a fundamental problem in biomedical imaging and medical science. This work introduces an approach for analyzing diffusion magnetic resonance imaging data acquired by the novel tensor-valued encoding technique for characterizing tissue microstructure. Our approach first uses a signal model to estimate the variance and skewness of the distribution of apparent diffusion tensors modeling the underlying tissue. Then several novel imaging indices, such as weighted microscopic anisotropy and microscopic skewness, are derived to characterize different ensembles of diffusion processes that are indistinguishable by existing techniques. The contributions of this work also include a theoretical proof that shows that, to estimate the skewness of a diffusion tensor distribution, the encoding protocol needs to include full-rank tensor diffusion encoding. This proof provides a guideline for the application of this technique. The properties of the proposed indices are illustrated using both synthetic data and in vivo data acquired from a human brain.

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.

Assessment of the structural complexity of diffusion MRI voxels using 3D electron microscopy in the rat brain

Validation and interpretation of diffusion magnetic resonance imaging (dMRI) requires detailed understanding of the actual microstructure restricting the diffusion of water molecules. In this study, we used serial block-face scanning electron microscopy (SBEM), a three-dimensional electron microscopy (3D-EM) technique, to image seven white and grey matter volumes in the rat brain. SBEM shows excellent contrast of cellular membranes, which are the major components restricting the diffusion of water in tissue. Additionally, we performed 3D structure tensor (3D-ST) analysis on the SBEM volumes and parameterised the resulting orientation distributions using Watson and angular central Gaussian (ACG) probability distributions as well as spherical harmonic (SH) decomposition. We analysed how these parameterisations described the underlying orientation distributions and compared their orientation and dispersion with corresponding parameters from two dMRI methods, neurite orientation dispersion and density imaging (NODDI) and constrained spherical deconvolution (CSD). Watson and ACG parameterisations and SH decomposition captured well the 3D-ST orientation distributions, but ACG and SH better represented the distributions due to its ability to model asymmetric dispersion. The dMRI parameters corresponded well with the 3D-ST parameters in the white matter volumes, but the correspondence was less evident in the more complex grey matter. SBEM imaging and 3D-ST analysis also revealed that the orientation distributions were often not axially symmetric, a property neatly captured by the ACG distribution. Overall, the ability of SBEM to image diffusion barriers in intricate detail, combined with 3D-ST analysis and parameterisation, provides a step forward toward interpreting and validating the dMRI signals in complex brain tissue microstructure.

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.

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.

Machine learning based white matter models with permeability: An experimental study in cuprizone treated in-vivo mouse model of axonal demyelination

The intra-axonal water exchange time (τ_i), a parameter associated with axonal permeability, could be an important biomarker for understanding and treating demyelinating pathologies such as Multiple Sclerosis. Diffusion-Weighted MRI (DW-MRI) is sensitive to changes in permeability; however, the parameter has so far remained elusive due to the lack of general biophysical models that incorporate it. Machine learning based computational models can potentially be used to estimate such parameters. Recently, for the first time, a theoretical framework using a random forest (RF) regressor suggests that this is a promising new approach for permeability estimation. In this study, we adopt such an approach and for the first time experimentally investigate it for demyelinating pathologies through direct comparison with histology. We construct a computational model using Monte Carlo simulations and an RF regressor in order to learn a mapping between features derived from DW-MRI signals and ground truth microstructure parameters. We test our model in simulations, and find strong correlations between the predicted and ground truth parameters (intra-axonal volume fraction f: R² =0.99, τ_i: R² =0.84, intrinsic diffusivity d: R² =0.99). We then apply the model in-vivo, on a controlled cuprizone (CPZ) mouse model of demyelination, comparing the results from two cohorts of mice, CPZ (N=8) and healthy age-matched wild-type (WT, N=8). We find that the RF model estimates sensible microstructure parameters for both groups, matching values found in literature. Furthermore, we perform histology for both groups using electron microscopy (EM), measuring the thickness of the myelin sheath as a surrogate for exchange time. Histology results show that our RF model estimates are very strongly correlated with the EM measurements (ρ = 0.98 for f, ρ = 0.82 for τ_i). Finally, we find a statistically significant decrease in τ_i in all three regions of the corpus callosum (splenium/genu/body) of the CPZ cohort (<τ_i>=310ms/330ms/350ms) compared to the WT group (<τ_i>=370ms/370ms/380ms). This is in line with our expectations that τ_i is lower in regions where the myelin sheath is damaged, as axonal membranes become more permeable. Overall, these results demonstrate, for the first time experimentally and in vivo, that a computational model learned from simulations can reliably estimate microstructure parameters, including the axonal permeability .

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.

The impact of realistic axonal shape on axon diameter estimation using diffusion MRI

To study axonal microstructure with diffusion MRI, axons are typically modeled as straight impermeable cylinders, whereby the transverse diffusion MRI signal can be made sensitive to the cylinder’s inner diameter. However, the shape of a real axon varies along the axon direction, which couples the longitudinal and transverse diffusion of the overall axon direction. Here we develop a theory of the intra-axonal diffusion MRI signal based on coarse-graining of the axonal shape by 3-dimensional diffusion. We demonstrate how the estimate of the inner diameter is confounded by the diameter variations (beading), and by the local variations in direction (undulations) along the axon. We analytically relate diffusion MRI metrics, such as time-dependent radial diffusivity D_⊥(t) and kurtosis K_⊥(t), to the axonal shape, and validate our theory using Monte Carlo simulations in synthetic undulating axons with randomly positioned beads, and in realistic axons reconstructed from electron microscopy images of mouse brain white matter. We show that (i) In the narrow pulse limit, the inner diameter from D_⊥(t) is overestimated by about twofold due to a combination of axon caliber variations and undulations (each contributing a comparable effect size); (ii) The narrow-pulse kurtosis K_⊥∣_t→∞ deviates from that in an ideal cylinder due to caliber variations; we also numerically calculate the fourth-order cumulant for an ideal cylinder in the wide pulse limit, which is relevant for inner diameter overestimation; (iii) In the wide pulse limit, the axon diameter overestimation is mainly due to undulations at low diffusion weightings b; and (iv) The effect of undulations can be considerably reduced by directional averaging of high-b signals, with the apparent inner diameter given by a combination of the axon caliber (dominated by the thickest axons), caliber variations, and the residual contribution of undulations.

A time-dependent diffusion MRI signature of axon caliber variations and beading

MRI provides a unique non-invasive window into the brain, yet is limited to millimeter resolution, orders of magnitude coarser than cell dimensions. Here, we show that diffusion MRI is sensitive to the micrometer-scale variations in axon caliber or pathological beading, by identifying a signature power-law diffusion time-dependence of the along-fiber diffusion coefficient. We observe this signature in human brain white matter and identify its origins by Monte Carlo simulations in realistic substrates from 3-dimensional electron microscopy of mouse corpus callosum. Simulations reveal that the time-dependence originates from axon caliber variation, rather than from mitochondria or axonal undulations. We report a decreased amplitude of time-dependence in multiple sclerosis lesions, illustrating the potential sensitivity of our method to axonal beading in a plethora of neurodegenerative disorders. This specificity to microstructure offers an exciting possibility of bridging across scales to image cellular-level pathology with a clinically feasible MRI technique.

Strong diffusion gradients allow the separation of intra- and extra-axonal gradient-echo signals in the human brain

The quantification of brain white matter properties is a key area of application of Magnetic Resonance Imaging (MRI), with much effort focused on using MR techniques to quantify tissue microstructure. While diffusion MRI probes white matter (WM) microstructure by characterising the sensitivity of Brownian motion of water molecules to anisotropic structures, susceptibility-based techniques probe the tissue microstructure by observing the effect of interaction between the tissue and the magnetic field. Here, we unify these two complementary approaches by combining ultra-strong (300 mT/m) gradients with a novel Diffusion-Filtered Asymmetric Spin Echo (D-FASE) technique. Using D-FASE we can separately assess the evolution of the intra- and extra-axonal signals under the action of susceptibility effects, revealing differences in the behaviour in different fibre tracts. We observed that the effective relaxation rate of the ASE signal in the corpus callosum decreases with increasing b-value in all subjects (from 17.1 ± 0.7 s⁻¹ at b = 0 s/mm² to 14.6 ± 0.7 s⁻¹ at b = 4800 s/mm²), while this dependence on b in the corticospinal tract is less pronounced (from 12.0± 1.1 s⁻¹ at b = 0s/mm² to 10.7 ± 0.5 s⁻¹ at b = 4800 s/mm²). Voxelwise analysis of the signal evolution with respect to b-factor and acquisition delay using a microscopic model demonstrated differences in gradient echo signal evolution between the intra- and extra-axonal pools.

Optimization of q-space sampling for mean apparent propagator MRI metrics using a genetic algorithm

Mean Apparent Propagator (MAP) MRI is a recently introduced technique to estimate the diffusion probability density function (PDF) robustly. Using the estimated PDF, MAP MRI then calculates zero-displacement and non-Gaussianity metrics, which might better characterize tissue microstructure compared to diffusion tensor imaging or diffusion kurtosis imaging. However, intensive q-space sampling required for MAP MRI limits its widespread adoption. A reduced q-space sampling scheme that maintains the accuracy of the derived metrics would make it more practical. A heuristic approach for acquiring MAP MRI with fewer q-space samples has been introduced earlier with scan duration of less than 10 minutes. However, the sampling scheme was not optimized systematically to preserve the accuracy of the model metrics. In this work, a genetic algorithm is implemented to determine optimal q-space subsampling schemes for MAP MRI that will keep total scan time under 10 min. Results show that the metrics derived from the optimized schemes more closely match those computed from the full set, especially in dense fiber tracts such as the corpus callosum.

Along-axon diameter variation and axonal orientation dispersion revealed with 3D electron microscopy: implications for quantifying brain white matter microstructure with histology and diffusion MRI

Tissue microstructure modeling of diffusion MRI signal is an active research area striving to bridge the gap between macroscopic MRI resolution and cellular-level tissue architecture. Such modeling in neuronal tissue relies on a number of assumptions about the microstructural features of axonal fiber bundles, such as the axonal shape (e.g., perfect cylinders) and the fiber orientation dispersion. However, these assumptions have not yet been validated by sufficiently high-resolution 3-dimensional histology. Here, we reconstructed sequential scanning electron microscopy images in mouse brain corpus callosum, and introduced a random-walker (RaW)-based algorithm to rapidly segment individual intra-axonal spaces and myelin sheaths of myelinated axons. Confirmed by a segmentation based on human annotations initiated with conventional machine-learning-based carving, our semi-automatic algorithm is reliable and less time-consuming. Based on the segmentation, we calculated MRI-relevant estimates of size-related parameters (inner axonal diameter, its distribution, along-axon variation, and myelin g-ratio), and orientation-related parameters (fiber orientation distribution and its rotational invariants; dispersion angle). The reported dispersion angle is consistent with previous 2-dimensional histology studies and diffusion MRI measurements, while the reported diameter exceeds those in other mouse brain studies. Furthermore, we calculated how these quantities would evolve in actual diffusion MRI experiments as a function of diffusion time, thereby providing a coarse-graining window on the microstructure, and showed that the orientation-related metrics have negligible diffusion time-dependence over clinical and pre-clinical diffusion time ranges. However, the MRI-measured inner axonal diameters, dominated by the widest cross sections, effectively decrease with diffusion time by ~ 17% due to the coarse-graining over axonal caliber variations. Furthermore, our 3d measurement showed that there is significant variation of the diameter along the axon. Hence, fiber orientation dispersion estimated from MRI should be relatively stable, while the "apparent" inner axonal diameters are sensitive to experimental settings, and cannot be modeled by perfectly cylindrical axons.

Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation

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

Glutamate diffusion in the rat brain in vivo under light and deep anesthesia conditions

PURPOSE

Glutamate (Glu) is the most abundant neurotransmitter in the human central nervous system and glutamatergic neurotransmission has been implicated in many common and severe neuropsychiatric disorders. In vivo MRS techniques have been developed to measure brain Glu concentration to investigate the pathophysiology of various brain disorders. However, it is difficult to interpret Glu signal changes because Glu plays multiple roles in the brain and is found in multiple microenvironments including cytosolic, vesicular, and extracellular.

METHODS

In vivo diffusion-weighted MRS (DW-MRS) with low to very high b-values was performed on the rat prefrontal cortex at 9.4T under both light and deep anesthetic conditions to examine Glu diffusion properties.

RESULTS

Significant alterations in Glu diffusion as well as reduced Glu concentration were observed under deep anesthesia compared with superficial anesthesia in the absence of similar changes in NAA or creatine.

CONCLUSION

The modifications in Glu diffusion under deep anesthesia might reflect changes in Glu microenvironment. The present work shows that Glu DW-MRS could be an important tool to explore Glu physiology with changing levels of neuronal activity and synaptic function.

Larmor frequency in heterogeneous media

High-precision signal phase measurements ignited an ongoing discussion of the microstructural correlates of the Larmor frequency shift of water in biological tissues. In a broader context, this is the question about the averaged precession frequency in magnetically heterogenous, in particular, porous media. In this study, the Larmor frequency shift is found analytically for water filling connected pore space between NMR-invisible magnetized inclusions with a constant magnetic susceptibility tensor. The magnetic microstructure that encompasses the inclusions' shape and their spatial arrangement is arbitrary as well as the inclusions' magnetic susceptibility tensor. The result is limited to the case of effectively fast diffusion of water molecules. In this limit, the effect of magnetic microstructure on the Larmor frequency shift is represented by only five relevant parameters in the general case and by a single parameter in the case of axially symmetric microstructure. This single parameter enters the dependence of the Larmor frequency on the sample orientation relative to the main magnetic field in the previously performed experiments. The result can help interpreting known experimental data and developing realistic models of biological tissues.

Evidence of early microstructural white matter abnormalities in multiple sclerosis from multi-shell diffusion MRI

Irreversible white matter (WM) damage, including severe demyelination and axonal loss, is a main determinant of long-term disability in multiple sclerosis (MS). Non-invasive detection of changes in microstructural WM integrity in the disease is challenging since commonly used imaging metrics lack the necessary sensitivity, especially in the early phase of the disease. This study aims at assessing microstructural WM abnormalities in early-stage MS by using ultra-high gradient strength multi-shell diffusion MRI and the restricted signal fraction (FR) from the Composite Hindered and Restricted Model of Diffusion (CHARMED), a metric sensitive to the volume fraction of axons. In 22 early MS subjects (disease duration ≤5 years) and 15 age-matched healthy controls, restricted fraction estimates were obtained through the CHARMED model along with conventional Diffusion Tensor Imaging (DTI) metrics. All imaging parameters were compared cross-sectionally between the MS subjects and controls both in WM lesions and normal-appearing white matter (NAWM). We found a significant reduction in FR focally in WM lesions and widespread in the NAWM in MS patients relative to controls (corrected p < .05). Signal fraction changes in NAWM were not driven by perilesional tissue, nor were they influenced by proximity to the ventricles, challenging the hypothesis of an outside-in pathological process driven by CSF-mediated immune cytotoxic factors. No significant differences were found in conventional DTI parameters. In a cross-validated classification task, FR showed the largest effect size and outperformed all other diffusion imaging metrics in discerning lesions from contralateral NAWM. Taken together, our data provide evidence for the presence of widespread microstructural changes in the NAWM in early MS stages that are, at least in part, unrelated to focal demyelinating lesions. Interestingly, these pathological changes were not yet detectable by conventional diffusion imaging at this early disease stage, highlighting the sensitivity and value of multi-shell diffusion imaging for better characterizing axonal microstructure in MS.

Evidence for Progressive Microstructural Damage in Early Multiple Sclerosis by Multi-Shell Diffusion Magnetic Resonance Imaging

In multiple sclerosis (MS), it would be of clinical value to be able to track the progression of axonal pathology, especially before the manifestation of clinical disability. However, non-invasive evaluation of short-term longitudinal progression of white matter integrity is challenging. This study aims at assessing longitudinal changes in the restricted (i.e. intracellular) diffusion signal fraction (FR) in early-stage MS by using ultra-high gradient strength multi-shell diffusion magnetic resonance imaging. In 11 early MS subjects (disease duration ≤5 years), FR was obtained at two timepoints (one year apart) through the Composite Hindered and Restricted Model of Diffusion, along with conventional Diffusion Tensor Imaging metrics. At follow-up, no statistically significant change was detected in clinical variables, while all imaging metrics showed statistically significant longitudinal changes (p < 0.01, corrected for multiple comparisons) in widespread regions in normal-appearing white matter (NAWM). The most extensive longitudinal changes were observed in FR, including areas known to include a large fraction of crossing fibers. Furthermore, FR was also the only metric showing significant longitudinal changes in lesions that were present at both time points (p = 0.007), with no significant differences found for conventional diffusion metrics. Finally, FR was the only diffusion metric (as compared to Diffusion Tensor Imaging) that revealed pre-lesional changes already present at baseline. Taken together, our data provide evidence for progressive microstructural damage in the NAWM of early MS cases detectable already at 1-year follow-up. Our study highlights the value of multi-shell diffusion imaging for sensitive tracking of disease evolution in MS before any clinical changes are observed. This article is part of a Special Issue entitled: SI: MRI and Neuroinflammation.

Diffusion time dependence of microstructural parameters in fixed spinal cord

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

On the scaling behavior of water diffusion in human brain white matter

Development of therapies for neurological disorders depends on our ability to non-invasively diagnose and monitor the progression of underlying pathologies at the cellular level. Physics and physiology limit the resolution of human MRI to be orders of magnitude coarser than cell dimensions. Here we identify and quantify the MRI signal coming from within micrometer-thin axons in human white matter tracts in vivo, by utilizing the sensitivity of diffusion MRI to Brownian motion of water molecules restricted by cell walls. We study a specific power-law scaling of the diffusion MRI signal with the diffusion weighting, predicted for water confined to narrow axons, and quantify axonal water fraction and orientation dispersion.

Diffusional kurtosis imaging and white matter microstructure modeling in a clinical study of major depressive disorder

Major depressive disorder (MDD) is a globally prevalent psychiatric disorder that results from disruption of multiple neural circuits involved in emotional regulation. Although previous studies using diffusion tensor imaging (DTI) found smaller values of fractional anisotropy (FA) in the white matter, predominantly in the frontal lobe, of patients with MDD, studies using diffusion kurtosis imaging (DKI) are scarce. Here, we used DKI whole-brain analysis with tract-based spatial statistics (TBSS) to investigate the brain microstructural abnormalities in MDD. Twenty-six patients with MDD and 42 age- and sex-matched control subjects were enrolled. To investigate the microstructural pathology underlying the observations in DKI, a compartment model analysis was conducted focusing on the corpus callosum. In TBSS, the patients with MDD showed significantly smaller values of FA in the genu and frontal portion of the body of the corpus callosum. The patients also had smaller values of mean kurtosis (MK) and radial kurtosis (RK), but MK and RK abnormalities were distributed more widely compared with FA, predominantly in the frontal lobe but also in the parietal, occipital, and temporal lobes. Within the callosum, the regions with smaller MK and RK were located more posteriorly than the region with smaller FA. Model analysis suggested significantly smaller values of intra-neurite signal fraction in the body of the callosum and greater fiber dispersion in the genu, which were compatible with the existing literature of white matter pathology in MDD. Our results show that DKI is capable of demonstrating microstructural alterations in the brains of patients with MDD that cannot be fully depicted by conventional DTI. Though the issues of model validation and parameter estimation still remain, it is suggested that diffusion MRI combined with a biophysical model is a promising approach for investigation of the pathophysiology of MDD.

Ultra-high-b radial diffusion-weighted imaging (UHb-rDWI) of human cervical spinal cord

BACKGROUND

Injury in the cervical spinal cord (CSC) can lead to varying degrees of neurologic deficit and persistent disability. Diffusion tensor imaging (DTI) is a promising method to evaluate white matter integrity and pathology. However, the conventional DTI results are limited with respect to the specific details of neuropathology and microstructural architecture. In this study we used ultrahigh-b radial-DWI (UHb-rDWI) with b-values ranging from 0 to ∼7500 s/mm² and calculated decay constant (D_H ) at the high b-values, which gives much deeper insight about the microscopic environment of CSC white matter.

PURPOSE

To evaluate a novel diffusion MRI, UHb-rDWI technique for imaging of the CSC.

STUDY TYPE

Longitudinal.

SUBJECTS

Four healthy controls, each scanned twice.

FIELD STRENGTH/SEQUENCE

3T/2D single shot diffusion-weighted stimulated echo planar imaging with reduced field of view.

ASSESSMENT

The signal from each pixel of b₀ (b = 0) and b-value (b ≠ 0) images were fitted to a biexponential function and normalized. The signal-b curve is obtained by dividing the latter curve by the former. D_H was obtained from the curve at b >4000 s/mm² . A Monte-Carlo Simulation (MCS) was performed to investigate how D_H changes upon the increased water-exchange at the CSC.

RESULTS

The signal-b curves plotted at multiple levels of healthy CSC are almost identical on two successive scans and show a biexponential decay behavior: fast exponential decay at lower b-values and much slower decay at UHb-values. The mean values of D_H were measured as (0.0607 ± 0.02531) ×10^-3 and (0.0357 ± 0.02072) ×10^-3 s/mm² at the lateral funiculus and posterior column, respectively. MCS of diffusion MRI shows that the D_H is elevated by increased water exchange between the intra- and extraaxonal spaces.

DATA CONCLUSION

UHb-rDWI signal-b plots of the normal CSC were highly reproducible on successive scans and their biexponential decay behavior can be used to characterize normal spinal white matter.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:204-211.

Rotationally-invariant mapping of scalar and orientational metrics of neuronal microstructure with diffusion MRI

We develop a general analytical and numerical framework for estimating intra- and extra-neurite water fractions and diffusion coefficients, as well as neurite orientational dispersion, in each imaging voxel. By employing a set of rotational invariants and their expansion in the powers of diffusion weighting, we analytically uncover the nontrivial topology of the parameter estimation landscape, showing that multiple branches of parameters describe the measurement almost equally well, with only one of them corresponding to the biophysical reality. A comprehensive acquisition shows that the branch choice varies across the brain. Our framework reveals hidden degeneracies in MRI parameter estimation for neuronal tissue, provides microstructural and orientational maps in the whole brain without constraints or priors, and connects modern biophysical modeling with clinical MRI.

Confirmation of a gyral bias in diffusion MRI fiber tractography

Diffusion MRI fiber tractography has been increasingly used to map the structural connectivity of the human brain. However, this technique is not without limitations; for example, there is a growing concern over anatomically correlated bias in tractography findings. In this study, we demonstrate that there is a bias for fiber tracking algorithms to terminate preferentially on gyral crowns, rather than the banks of sulci. We investigate this issue by comparing diffusion MRI (dMRI) tractography with equivalent measures made on myelin-stained histological sections. We begin by investigating the orientation and trajectories of axons near the white matter/gray matter boundary, and the density of axons entering the cortex at different locations along gyral blades. These results are compared with dMRI orientations and tract densities at the same locations, where we find a significant gyral bias in many gyral blades across the brain. This effect is shown for a range of tracking algorithms, both deterministic and probabilistic, and multiple diffusion models, including the diffusion tensor and a high angular resolution diffusion imaging technique. Additionally, the gyral bias occurs for a range of diffusion weightings, and even for very high-resolution datasets. The bias could significantly affect connectivity results using the current generation of tracking algorithms.

Can increased spatial resolution solve the crossing fiber problem for diffusion MRI?

It is now widely recognized that voxels with crossing fibers or complex geometrical configurations present a challenge for diffusion MRI (dMRI) reconstruction and fiber tracking, as well as microstructural modeling of brain tissues. This "crossing fiber" problem has been estimated to affect anywhere from 30% to as many as 90% of white matter voxels, and it is often assumed that increasing spatial resolution will decrease the prevalence of voxels containing multiple fiber populations. The aim of this study is to estimate the extent of the crossing fiber problem as we progressively increase the spatial resolution, with the goal of determining whether it is possible to mitigate this problem with higher resolution spatial sampling. This is accomplished using ex vivo MRI data of the macaque brain, followed by histological analysis of the same specimen to validate these measurements, as well as to extend this analysis to resolutions not yet achievable in practice with MRI. In both dMRI and histology, we find unexpected results: the prevalence of crossing fibers increases as we increase spatial resolution. The problem of crossing fibers appears to be a fundamental limitation of dMRI associated with the complexity of brain tissue, rather than a technical problem that can be overcome with advances such as higher fields and stronger gradients.

Design and validation of diffusion MRI models of white matter

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

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

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