The importance of axonal undulation in diffusion MR measurements: a Monte Carlo simulation study

From structure to mechanics: exploring the role of axons and interconnections in anisotropic behavior of brain white matter

According to various experimental studies, the role of axons in the brain's white matter (WM) is still a subject of debate: Is the role of axons in brain white matter (WM) limited to their functional significance, or do they also play a pivotal mechanical role in defining its anisotropic behavior? Micromechanics and computational models provide valuable tools for scientists to comprehend the underlying mechanisms of tissue behavior, taking into account the contribution of microstructures. In this review, we delve into the consideration of strain level, strain rates, and injury threshold to determine when WM should be regarded as anisotropic, as well as when the assumption of isotropy can be deemed acceptable. Additionally, we emphasize the potential mechanical significance of interconnections between glial cells-axons and glial cells-vessels. Moreover, we elucidate the directionality of WM stiffness under various loading conditions and define the possible roles of microstructural components in each scenario. Ultimately, this review aims to shed light on the significant mechanical contributions of axons in conjunction with glial cells, paving the way for the development of future multiscale models capable of predicting injuries and facilitating the discovery of applicable treatments.

Predicting Mesoscopic Larmor Frequency Shifts in White Matter With Diffusion MRI-A Monte Carlo Study in Axonal Phantoms

Magnetic susceptibility MRI offers potential insights into the chemical composition and microstructural organization of tissue. However, estimating magnetic susceptibility in white matter is challenging due to anisotropic subvoxel Larmor frequency shifts caused by axonal microstructure relative to the B0 field orientation. Recent biophysical models have analytically described how axonal microstructure influences the Larmor frequency shifts, relating these shifts to a mesoscopically averaged magnetic field that depends on the axons' fiber orientation distribution function (fODF), typically estimated using diffusion MRI. This study is aimed at validating the use of MRI to estimate mesoscopic magnetic fields and determining whether diffusion MRI can faithfully estimate the orientation dependence of the Larmor frequency shift in realistic axonal microstructure. To achieve this, we developed a framework for performing Monte Carlo simulations of MRI signals in mesoscopically sized white matter axon substrates segmented with electron microscopy. Our simulations demonstrated that with careful experimental design, it is feasible to estimate mesoscopic magnetic fields. Additionally, the fODF estimated by the standard model of diffusion in white matter could predict the orientation dependence of the mesoscopic Larmor frequency shift. We also found that incorporating the intra-axonal axial kurtosis into the standard model could explain a significant amount of signal variance, thereby improving the estimation of the Larmor frequency shift. This factor should not be neglected when fitting the standard model.

Bridging the 3D geometrical organisation of white matter pathways across anatomical length scales and species

We used diffusion MRI and x-ray synchrotron imaging on monkey and mice brains to examine the organisation of fibre pathways in white matter across anatomical scales. We compared the structure in the corpus callosum and crossing fibre regions and investigated the differences in cuprizone-induced demyelination in mouse brains versus healthy controls. Our findings revealed common principles of fibre organisation that apply despite the varying patterns observed across species; small axonal fasciculi and major bundles formed laminar structures with varying angles, according to the characteristics of major pathways. Fasciculi exhibited non-straight paths around obstacles like blood vessels, comparable across the samples of varying fibre complexity and demyelination. Quantifications of fibre orientation distributions were consistent across anatomical length scales and modalities, whereas tissue anisotropy had a more complex relationship, both dependent on the field-of-view. Our study emphasises the need to balance field-of-view and voxel size when characterising white matter features across length scales.

Magnetic resonance spectroscopy in cervical spondylotic myelopathy: a systematic review of metabolite changes and clinical correlations

Introduction

Cervical spondylotic myelopathy (CSM) is a common degenerative condition characterized by narrowing of the cervical spinal canal, leading to progressive spinal cord injury and functional decline. While magnetic resonance imaging (MRI) is the gold standard for diagnosing CSM, it has limitations in predicting clinical outcomes. Magnetic resonance spectroscopy (MRS) offers metabolic insights that may enhance diagnostic and prognostic capabilities in CSM.

Methods

We conducted a systematic review following the PRISMA guidelines. Comprehensive literature searches were performed in PubMed, OVID, EMBASE, Web of Science, and the Cochrane Central Register of Controlled Trials up to June 2023. Studies included human subjects with CSM, a cohort of at least 10 patients, and reported primary data on cervical spine MRS findings correlated with clinical scales such as the modified Japanese Orthopaedic Association (mJOA) scale, both pre- and post-operatively.

Results

Six prospective studies involving 123 patients (average age 45.8 to 63 years) met the inclusion criteria. Common symptoms were neck pain, radicular upper-limb pain, paresthesia, and motor impairment. MRS findings indicated that symptomatic CSM patients had reduced N-acetyl aspartate to creatine (NAA/Cr) ratios and elevated choline to creatine (Cho/Cr) and choline to NAA (Cho/NAA) ratios compared to healthy controls. Lactate peaks were detected in a significant proportion of symptomatic patients, suggesting hypoxic or inflammatory injury. Decreased NAA/Cr and increased Cho/NAA ratios correlated with lower mJOA scores, indicating more severe myelopathy. Post-operative increases in NAA/Cr ratios and decreases in Cho/NAA ratios were associated with improved mJOA scores, highlighting the prognostic value of these metabolites.

Conclusion

MRS provides valuable metabolic information correlating with clinical severity and functional outcomes in CSM. Reduced NAA/Cr and elevated Cho/Cr and Cho/NAA ratios are associated with more severe disease and may predict surgical recovery. MRS shows promise as a non-invasive tool for enhancing the diagnosis and management of CSM. Further research is needed to standardize protocols, validate findings in larger cohorts, and integrate MRS into clinical practice.

Susceptibility-induced internal gradients reveal axon morphology and cause anisotropic effects in the diffusion-weighted MRI signal

Diffusion-weighted MRI is our most promising method for estimating microscopic tissue morphology in vivo. The signal acquisition is based on scanner-generated external magnetic gradients. However, it will also be affected by susceptibility-induced internal magnetic gradients caused by interactions between the tissue and the static magnetic field of the scanner. With 3D in silico experiments, we show how internal gradients cause morphology-, compartment-, and orientation-dependence of spin-echo and pulsed-gradient spin-echo experiments in myelinated axons. These effects surpass those observed with previous 2D modelling corresponding to straight cylinders. For an ex vivo monkey brain, we observe the orientation-dependence generated only when including non-circular cross-sections in the in silico morphological configurations, and find orientation-dependent deviation of up to 17% for diffusion tensor metrics. Interestingly, we find that the orientation-dependence not only biases the signal across different brain regions, but also carries a sensitivity to the morphology of axonal cross-sections which is not attainable by the idealised theoretical diffusion-weighted MRI signal.

Comparisons of MR and EM inferred tissue microstructure properties using a human autopsy corpus callosum sample

Degeneration of white matter (WM) microstructure in the central nervous system is characteristic of many neurodegenerative conditions. Previous research indicates that axonal degeneration visible in ex vivo electron microscopy (EM) photomicrographs precede the onset of clinical symptoms. Measuring WM microstructural features, such as axon diameter and packing fraction, currently require these highly invasive methods of analysis and it is therefore of great importance to develop methods for in vivo measurements. Diffusion weighted Magnetic Resonance Imaging (MRI) is a non-invasive method which can be used in conjunction with temporal diffusion spectroscopy (TDS) and an oscillating gradient spin echo (OGSE) pulse sequence to probe micron-scale structures within neural tissue. The current experiment aims to compare axon diameter measurements, mean effective axon diameter (AxD¯), and packing fractions calculated from EM histopathological analysis and inferred values from MR images. Mathematical models of axon diameters used for analysis include the ActiveAx Frequency-Dependent Extra-Axonal Diffusion (AAD) model and the AxCaliber Frequency-Dependent Extra-Axonal Diffusion (ACD) model using ROI (Region of Interest) based analysis (RBA) and voxel-based analysis (VBA), respectively. Overall, it was observed that MRI inferred WM microstructural parameters overestimate those calculated from EM. This may be attributable to tissue shrinkage during EM dehydration, the sensitivity of MR pulse sequences to larger diameter axons, and/or inaccurate model assumptions. The results of the current study provide a means to characterize the precision and accuracy of RBA-ACD and VBA-AAD OGSE-TDS and highlight the need for further research investigating the relationship between ex vivo MRI and EM, with the goal of reaching in vivo MRI.

A comprehensive review on modeling aspects of infusion-based drug delivery in the brain

Brain disorders represent an ever-increasing health challenge worldwide. While conventional drug therapies are less effective due to the presence of the blood-brain barrier, infusion-based methods of drug delivery to the brain represent a promising option. Since these methods are mechanically controlled and involve multiple physical phases ranging from the neural and molecular scales to the brain scale, highly efficient and precise delivery procedures can significantly benefit from a comprehensive understanding of drug-brain and device-brain interactions. Behind these interactions are principles of biophysics and biomechanics that can be described and captured using mathematical models. Although biomechanics and biophysics have received considerable attention, a comprehensive mechanistic model for modeling infusion-based drug delivery in the brain has yet to be developed. Therefore, this article reviews the state-of-the-art mechanistic studies that can support the development of next-generation models for infusion-based brain drug delivery from the perspective of fluid mechanics, solid mechanics, and mathematical modeling. The supporting techniques and database are also summarized to provide further insights. Finally, the challenges are highlighted and perspectives on future research directions are provided. STATEMENT OF SIGNIFICANCE: Despite the immense potential of infusion-based drug delivery methods for bypassing the blood-brain barrier and efficiently delivering drugs to the brain, achieving optimal drug distribution remains a significant challenge. This is primarily due to our limited understanding of the complex interactions between drugs and the brain that are governed by principles of biophysics and biomechanics, and can be described using mathematical models. This article provides a comprehensive review of state-of-the-art mechanistic studies that can help to unravel the mechanism of drug transport in the brain across the scales, which underpins the development of next-generation models for infusion-based brain drug delivery. More broadly, this review will serve as a starting point for developing more effective treatments for brain diseases and mechanistic models that can be used to study other soft tissue and biomaterials.

Generalized Langevin subdiffusion in channels: The bath always wins

We consider subdiffusive motion, modeled by the generalized Langevin equation in an equilibrium setting, of tracer particles in channels of indefinite length in the x direction: the channels of varying width and the channels with sinusoidally meandering midline. The subdiffusion in the x direction is not affected by constraints put by the channel. This is especially astonishing for meandering channels whose centerline might be quite long. The same behavior is seen in a holonomic model of a bead on a sinusoidal and meandering wire, where some analytic insights are possible.

Exploring white matter dynamics and morphology through interactive numerical phantoms: the White Matter Generator

Brain white matter is a dynamic environment that continuously adapts and reorganizes in response to stimuli and pathological changes. Glial cells, especially, play a key role in tissue repair, inflammation modulation, and neural recovery. The movements of glial cells and changes in their concentrations can influence the surrounding axon morphology. We introduce the White Matter Generator (WMG) tool to enable the study of how axon morphology is influenced through such dynamical processes, and how this, in turn, influences the diffusion-weighted MRI signal. This is made possible by allowing interactive changes to the configuration of the phantom generation throughout the optimization process. The phantoms can consist of myelinated axons, unmyelinated axons, and cell clusters, separated by extra-cellular space. Due to morphological flexibility and computational advantages during the optimization, the tool uses ellipsoids as building blocks for all structures; chains of ellipsoids for axons, and individual ellipsoids for cell clusters. After optimization, the ellipsoid representation can be converted to a mesh representation which can be employed in Monte-Carlo diffusion simulations. This offers an effective method for evaluating tissue microstructure models for diffusion-weighted MRI in controlled bio-mimicking white matter environments. Hence, the WMG offers valuable insights into white matter's adaptive nature and implications for diffusion-weighted MRI microstructure models, and thereby holds the potential to advance clinical diagnosis, treatment, and rehabilitation strategies for various neurological disorders and injuries.

Fast multi-compartment Microstructure Fingerprinting in brain white matter

We proposed two deep neural network based methods to accelerate the estimation of microstructural features of crossing fascicles in the white matter. Both methods focus on the acceleration of a multi-dictionary matching problem, which is at the heart of Microstructure Fingerprinting, an extension of Magnetic Resonance Fingerprinting to diffusion MRI. The first acceleration method uses efficient sparse optimization and a dedicated feed-forward neural network to circumvent the inherent combinatorial complexity of the fingerprinting estimation. The second acceleration method relies on a feed-forward neural network that uses a spherical harmonics representation of the DW-MRI signal as input. The first method exhibits a high interpretability while the second method achieves a greater speedup factor. The accuracy of the results and the speedup factors of several orders of magnitude obtained on in vivo brain data suggest the potential of our methods for a fast quantitative estimation of microstructural features in complex white matter configurations.

Diffusion MRI marks progressive alterations in fiber integrity in the zQ175DN mouse model of Huntington's disease

Huntington's disease (HD) is a progressive neurodegenerative disease affecting motor and cognitive abilities. Multiple studies have found white matter anomalies in HD-affected humans and animal models of HD. The identification of sensitive white-matter-based biomarkers in HD animal models will be important in understanding disease mechanisms and testing the efficacy of therapeutic interventions. Here we investigated the progression of white matter deficits in the knock-in zQ175DN heterozygous (HET) mouse model of HD at 3, 6 and 11 months of age (M), reflecting different states of phenotypic progression. We compared findings from traditional diffusion tensor imaging (DTI) and advanced fixel-based analysis (FBA) diffusion metrics for their sensitivity in detecting white matter anomalies in the striatum, motor cortex, and segments of the corpus callosum. FBA metrics revealed progressive and widespread reductions of fiber cross-section and fiber density in myelinated bundles of HET mice. The corpus callosum genu was the most affected structure in HET mice at 6 and 11 M based on the DTI and FBA metrics, while the striatum showed the earliest progressive differences starting at 3 M based on the FBA metrics. Overall, FBA metrics detected earlier and more prominent alterations in myelinated fiber bundles compared to the DTI metrics. Luxol fast blue staining showed no loss in myelin density, indicating that diffusion anomalies could not be explained by myelin reduction but diffusion anomalies in HET mice were accompanied by increased levels of neurofilament light chain protein at 11 M. Altogether, our findings reveal progressive alterations in myelinated fiber bundles that can be measured using diffusion MRI, representing a candidate noninvasive imaging biomarker to study phenotype progression and the efficacy of therapeutic interventions in zQ175DN mice. Moreover, our study exposed higher sensitivity of FBA than DTI metrics, suggesting a potential benefit of adopting these advanced metrics in other contexts, including biomarker development in humans.

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 phantom study of fiber crossings at varied angles reconstructed with ODF-Fingerprinting

White matter fiber reconstructions based on seeking local maxima of Orientation Distribution Functions (ODFs) typically fail to identify fibers crossing at narrow angles below 45°. ODF-Fingerprinting (ODF-FP) replaces the ODF maxima localization mechanism with pattern matching, allowing the use of all information stored in ODFs. In this work, we study the ability of ODF-FP to reconstruct fibers crossing at varied angles spanning 10°-90° in physical diffusion phantoms composed of textile tubes with 0.8μm diameter, approaching the anatomical scale of axons. Our results show that ODF-FP is able to correctly identify 80 ± 8% of the crossing fibers regardless of the crossing angle and provide the highest average reconstruction accuracy.

Probing brain tissue microstructure with MRI: principles, challenges, and the role of multidimensional diffusion-relaxation encoding

Diffusion MRI uses the random displacement of water molecules to sensitize the signal to brain microstructure and to properties such as the density and shape of cells. Microstructure modeling techniques aim to estimate these properties from acquired data by separating the signal between virtual tissue 'compartments' such as the intra-neurite and the extra-cellular space. A key challenge is that the diffusion MRI signal is relatively featureless compared with the complexity of brain tissue. Another challenge is that the tissue microstructure is wildly different within the gray and white matter of the brain. In this review, we use results from multidimensional diffusion encoding techniques to discuss these challenges and their tentative solutions. Multidimensional encoding increases the information content of the data by varying not only the b-value and the encoding direction but also additional experimental parameters such as the shape of the b-tensor and the echo time. Three main insights have emerged from such encoding. First, multidimensional data contradict common model assumptions on diffusion and T2 relaxation, and illustrates how the use of these assumptions cause erroneous interpretations in both healthy brain and pathology. Second, many model assumptions can be dispensed with if data are acquired with multidimensional encoding. The necessary data can be easily acquired in vivo using protocols optimized to minimize Cramér-Rao lower bounds. Third, microscopic diffusion anisotropy reflects the presence of axons but not dendrites. This insight stands in contrast to current 'neurite models' of brain tissue, which assume that axons in white matter and dendrites in gray matter feature highly similar diffusion. Nevertheless, as an axon-based contrast, microscopic anisotropy can differentiate gray and white matter when myelin alterations confound conventional MRI contrasts.

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.

Monte Carlo simulation of water diffusion through cardiac tissue models

Monte Carlo diffusion simulations are commonly used to establish a reliable ground truth of tissue microstructure, including for the validation of diffusion-weighted MRI. However, selecting simulation parameters is challenging and affects validity and reproducibility. We conducted experiments to investigate critical conditions in Monte Carlo simulations, such as tissue representation complexity, simulated molecules, update duration, and compartment size. Results show significant changes in microstructure characteristics when parameters are altered, emphasizing the importance of careful control for a reliable ground truth.

CACTUS: a computational framework for generating realistic white matter microstructure substrates

Monte-Carlo diffusion simulations are a powerful tool for validating tissue microstructure models by generating synthetic diffusion-weighted magnetic resonance images (DW-MRI) in controlled environments. This is fundamental for understanding the link between micrometre-scale tissue properties and DW-MRI signals measured at the millimetre-scale, optimizing acquisition protocols to target microstructure properties of interest, and exploring the robustness and accuracy of estimation methods. However, accurate simulations require substrates that reflect the main microstructural features of the studied tissue. To address this challenge, we introduce a novel computational workflow, CACTUS (Computational Axonal Configurator for Tailored and Ultradense Substrates), for generating synthetic white matter substrates. Our approach allows constructing substrates with higher packing density than existing methods, up to 95% intra-axonal volume fraction, and larger voxel sizes of up to 500μm3 with rich fibre complexity. CACTUS generates bundles with angular dispersion, bundle crossings, and variations along the fibres of their inner and outer radii and g-ratio. We achieve this by introducing a novel global cost function and a fibre radial growth approach that allows substrates to match predefined targeted characteristics and mirror those reported in histological studies. CACTUS improves the development of complex synthetic substrates, paving the way for future applications in microstructure imaging.

Unraveling multi-fixel microstructure with tractography and angular weighting

Recent advances in MRI technology have enabled richer multi-shell sequences to be implemented in diffusion MRI, allowing the investigation of both the microscopic and macroscopic organization of the brain white matter and its complex network of neural fibers. The emergence of advanced diffusion models has enabled a more detailed analysis of brain microstructure by estimating the signal received from a voxel as the combination of responses from multiple fiber populations. However, disentangling the individual microstructural properties of different macroscopic white matter tracts where those pathways intersect remains a challenge. Several approaches have been developed to assign microstructural properties to macroscopic streamlines, but often present shortcomings. ROI-based heuristics rely on averages that are not tract-specific. Global methods solve a computationally-intensive global optimization but prevent the use of microstructural properties not included in the model and often require restrictive hypotheses. Other methods use atlases that might not be adequate in population studies where the shape of white matter tracts varies significantly between patients. We introduce UNRAVEL, a framework combining the microscopic and macroscopic scales to unravel multi-fixel microstructure by utilizing tractography. The framework includes commonly-used heuristics as well as a new algorithm, estimating the microstructure of a specific white matter tract with angular weighting. Our framework grants considerable freedom as the inputs required, a set of streamlines defining a tract and a multi-fixel diffusion model estimated in each voxel, can be defined by the user. We validate our approach on synthetic data and in vivo data, including a repeated scan of a subject and a population study of children with dyslexia. In each case, we compare the estimation of microstructural properties obtained with angular weighting to other commonly-used approaches. Our framework provides estimations of the microstructure at the streamline level, volumetric maps for visualization and mean microstructural values for the whole tract. The angular weighting algorithm shows increased accuracy, robustness to uncertainties in its inputs and maintains similar or better reproducibility compared to commonly-used analysis approaches. UNRAVEL will provide researchers with a flexible and open-source tool enabling them to study the microstructure of specific white matter pathways with their diffusion model of choice.

The biomechanical implications of neck position in cervical contusion animal models of SCI

Large animal contusion models of spinal cord injury are an essential precursor to developing and evaluating treatment options for human spinal cord injury. Reducing variability in these experiments has been a recent focus as it increases the sensitivity with which treatment effects can be detected while simultaneously decreasing the number of animals required in a study. Here, we conducted a detailed review to explore if head and neck positioning in a cervical contusion model of spinal cord injury could be a factor impacting the biomechanics of a spinal cord injury, and thus, the resulting outcomes. By reviewing existing literature, we found evidence that animal head/neck positioning affects the exposed level of the spinal cord, morphology of the spinal cord, tissue mechanics and as a result the biomechanics of a cervical spinal cord injury. We posited that neck position could be a hidden factor contributing to variability. Our results indicate that neck positioning is an important factor in studying biomechanics, and that reporting these values can improve inter-study consistency and comparability and that further work needs to be done to standardize positioning for cervical spinal cord contusion injury models.

Novel insights into axon diameter and myelin content in late childhood and adolescence

White matter microstructural development in late childhood and adolescence is driven predominantly by increasing axon density and myelin thickness. Ex vivo studies suggest that the increase in axon diameter drives developmental increases in axon density observed with pubertal onset. In this cross-sectional study, 50 typically developing participants aged 8-18 years were scanned using an ultra-strong gradient magnetic resonance imaging scanner. Microstructural properties, including apparent axon diameter $({d}_a)$, myelin content, and g-ratio, were estimated in regions of the corpus callosum. We observed age-related differences in ${d}_a$, myelin content, and g-ratio. In early puberty, males had larger ${d}_a$ in the splenium and lower myelin content in the genu and body of the corpus callosum, compared with females. Overall, this work provides novel insights into developmental, pubertal, and cognitive correlates of individual differences in apparent axon diameter and myelin content in the developing human brain.

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.

Temporal Diffusion Ratio (TDR) for imaging restricted diffusion: Optimisation and pre-clinical demonstration

Temporal Diffusion Ratio (TDR) is a recently proposed dMRI technique (Dell'Acqua et al., proc. ISMRM 2019) which provides contrast between areas with restricted diffusion and areas either without restricted diffusion or with length scales too small for characterisation. Hence, it has a potential for informing on pore sizes, in particular the presence of large axon diameters or other cellular structures. TDR employs the signal from two dMRI acquisitions obtained with the same, large, b-value but with different diffusion gradient waveforms. TDR is advantageous as it employs standard acquisition sequences, does not make any assumptions on the underlying tissue structure and does not require any model fitting, avoiding issues related to model degeneracy. This work for the first time introduces and optimises the TDR method in simulation for a range of different tissues and scanner constraints and validates it in a pre-clinical demonstration. We consider both substrates containing cylinders and spherical structures, representing cell soma in tissue. Our results show that contrasting an acquisition with short gradient duration, short diffusion time and high gradient strength with an acquisition with long gradient duration, long diffusion time and low gradient strength, maximises the TDR contrast for a wide range of pore configurations. Additionally, in the presence of Rician noise, computing TDR from a subset (50% or fewer) of the acquired diffusion gradients rather than the entire shell as proposed originally further improves the contrast. In the last part of the work the results are demonstrated experimentally on rat spinal cord. In line with simulations, the experimental data shows that optimised TDR improves the contrast compared to non-optimised TDR. Furthermore, we find a strong correlation between TDR and histology measurements of axon diameter. In conclusion, we find that TDR has great potential and is a very promising alternative (or potentially complement) to model-based approaches for informing on pore sizes and restricted diffusion in general.

Axial and radial axonal diffusivities and radii from single encoding strongly diffusion-weighted MRI

We enable the estimation of the per-axon axial diffusivity from single encoding, strongly diffusion-weighted, pulsed gradient spin echo data. Additionally, we improve the estimation of the per-axon radial diffusivity compared to estimates based on spherical averaging. The use of strong diffusion weightings in magnetic resonance imaging (MRI) allows to approximate the signal in white matter as the sum of the contributions from only axons. At the same time, spherical averaging leads to a major simplification of the modeling by removing the need to explicitly account for the unknown distribution of axonal orientations. However, the spherically averaged signal acquired at strong diffusion weightings is not sensitive to the axial diffusivity, which cannot therefore be estimated although needed for modeling axons - especially in the context of multi-compartmental modeling. We introduce a new general method for the estimation of both the axial and radial axonal diffusivities at strong diffusion weightings based on kernel zonal modeling. The method could lead to estimates that are free from partial volume bias with gray matter or other isotropic compartments. The method is tested on publicly available data from the MGH Adult Diffusion Human Connectome project. We report reference values of axonal diffusivities based on 34 subjects, and derive estimates of axonal radii from only two shells. The estimation problem is also addressed from the angle of the required data preprocessing, the presence of biases related to modeling assumptions, current limitations, and future possibilities.

Post mortem mapping of connectional anatomy for the validation of diffusion MRI

Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.

Insights into Infusion-Based Targeted Drug Delivery in the Brain: Perspectives, Challenges and Opportunities

Targeted drug delivery in the brain is instrumental in the treatment of lethal brain diseases, such as glioblastoma multiforme, the most aggressive primary central nervous system tumour in adults. Infusion-based drug delivery techniques, which directly administer to the tissue for local treatment, as in convection-enhanced delivery (CED), provide an important opportunity; however, poor understanding of the pressure-driven drug transport mechanisms in the brain has hindered its ultimate success in clinical applications. In this review, we focus on the biomechanical and biochemical aspects of infusion-based targeted drug delivery in the brain and look into the underlying molecular level mechanisms. We discuss recent advances and challenges in the complementary field of medical robotics and its use in targeted drug delivery in the brain. A critical overview of current research in these areas and their clinical implications is provided. This review delivers new ideas and perspectives for further studies of targeted drug delivery in the brain.

White matter markers and predictors for subject-specific rTMS response in major depressive disorder

Repetitive transcranial magnetic stimulation (rTMS) has established therapeutic efficacy for major depressive disorder (MDD). While translational research has focused primarily on understanding the mechanism of action of TMS on functional activation and connectivity, the effects on structural connectivity remain largely unknown especially when rTMS is applied using subject-specific brain targets. This study aims to use novel diffusion magnetic resonance imaging (dMRI) analysis to examine microstructural changes related to rTMS treatment response using a unique cohort of 21 patients with MDD treated using rTMS with subject-specific targets. White matter dMRI microstructural measures and clinical scores were captured before and after the full course of treatment. We defined disease-relevant fiber bundles connected to different subregions of the left prefrontal cortex and analyzed changes in diffusion properties as well as correlations between the changes of dMRI measures and the changes in Hamilton Depression Rating Scale (HAMD). No significant changes were observed in tracts connected to the TMS targets. rTMS significantly increased the extra-axonal free-water volume, fractional anisotropy and decreased the radial diffusivity in anterior-medial prefrontal fiber bundles but did not lead to raw changes in lateral prefrontal tracts. That said, the microstructural changes in the lateral prefrontal white matter were significantly correlated with treatment response. Moreover, pre-rTMS dMRI measures of the dorsal anterior cingulate cortex and lateral prefrontal cortex connections are correlated with changes in HAMD scores. Microstructural changes in the anterior-medial and lateral prefrontal white matter are potentially involved in treatment response to TMS, though further investigation is needed using larger datasets.

The direction-dependence of apparent water exchange rate in human white matter

Transmembrane water exchange is a potential biomarker in the diagnosis and understanding of cancers, brain disorders, and other diseases. Filter-exchange imaging (FEXI), a special case of diffusion exchange spectroscopy adapted for clinical applications, has the potential to reveal different physiological water exchange processes. However, it is still controversial whether modulating the diffusion encoding gradient direction can affect the apparent exchange rate (AXR) measurements of FEXI in white matter (WM) where water diffusion shows strong anisotropy. In this study, we explored the diffusion-encoding direction dependence of FEXI in human brain white matter by performing FEXI with 20 diffusion-encoding directions on a clinical 3T scanner in-vivo. The results show that the AXR values measured when the gradients are perpendicular to the fiber orientation (0.77 ± 0.13 s ^- 1, mean ± standard deviation of all the subjects) are significantly larger than the AXR estimates when the gradients are parallel to the fiber orientation (0.33 ± 0.14 s ^- 1, p < 0.001) in WM voxels with coherently-orientated fibers. In addition, no significant correlation is found between AXRs measured along these two directions, indicating that they are measuring different water exchange processes. What's more, only the perpendicular AXR rather than the parallel AXR shows dependence on axonal diameter, indicating that the perpendicular AXR might reflect transmembrane water exchange between intra-axonal and extra-cellular spaces. Further finite difference (FD) simulations having three water compartments (intra-axonal, intra-glial, and extra-cellular spaces) to mimic WM micro-environments also suggest that the perpendicular AXR is more sensitive to the axonal water transmembrane exchange than parallel AXR. Taken together, our results show that AXR measured along different directions could be utilized to probe different water exchange processes in WM.

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.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

SPHERIOUSLY? The challenges of estimating sphere radius non-invasively in the human brain from diffusion MRI

The Soma and Neurite Density Imaging (SANDI) three-compartment model was recently proposed to disentangle cylindrical and spherical geometries, attributed to neurite and soma compartments, respectively, in brain tissue. There are some recent advances in diffusion-weighted MRI signal encoding and analysis (including the use of multiple so-called 'b-tensor' encodings and analysing the signal in the frequency-domain) that have not yet been applied in the context of SANDI. In this work, using: (i) ultra-strong gradients; (ii) a combination of linear, planar, and spherical b-tensor encodings; and (iii) analysing the signal in the frequency domain, three main challenges to robust estimation of sphere size were identified: First, the Rician noise floor in magnitude-reconstructed data biases estimates of sphere properties in a non-uniform fashion. It may cause overestimation or underestimation of the spherical compartment size and density. This can be partly ameliorated by accounting for the noise floor in the estimation routine. Second, even when using the strongest diffusion-encoding gradient strengths available for human MRI, there is an empirical lower bound on the spherical signal fraction and radius that can be detected and estimated robustly. For the experimental setup used here, the lower bound on the sphere signal fraction was approximately 10%. We employed two different ways of establishing the lower bound for spherical radius estimates in white matter. The first, examining power-law relationships between the DW-signal and diffusion weighting in empirical data, yielded a lower bound of 7μm, while the second, pure Monte Carlo simulations, yielded a lower limit of 3μm and in this low radii domain, there is little differentiation in signal attenuation. Third, if there is sensitivity to the transverse intra-cellular diffusivity in cylindrical structures, e.g., axons and cellular projections, then trying to disentangle two diffusion-time-dependencies using one experimental parameter (i.e., change in frequency-content of the encoding waveform) makes spherical radii estimates particularly challenging. We conclude that due to the aforementioned challenges spherical radii estimates may be biased when the corresponding sphere signal fraction is low, which must be considered.

Impact of the acquisition protocol on the sensitivity to demyelination and axonal loss of clinically feasible DWI techniques: a simulation study

OBJECTIVE

To evaluate: (a) the specific effect that the demyelination and axonal loss have on the DW signal, and (b) the impact of the sequence parameters on the sensitivity to damage of two clinically feasible DWI techniques, i.e. DKI and NODDI.

METHODS

We performed a Monte Carlo simulation of water diffusion inside a novel synthetic model of white matter in the presence of axonal loss and demyelination, with three compartments with permeable boundaries between them. We compared DKI and NODDI in their ability to detect and assess the damage, using several acquisition protocols. We used the F test statistic as an index of the sensitivity for each DWI parameter to axonal loss and demyelination, respectively.

RESULTS

DKI parameters significantly changed with increasing axonal loss, but, in most cases, not with demyelination; all the NODDI parameters showed sensitivity to both the damage processes (at p < 0.01). However, the acquisition protocol strongly affected the sensitivity to damage of both the DKI and NODDI parameters and, especially for NODDI, the parameter absolute values also.

DISCUSSION

This work is expected to impact future choices for investigating white matter microstructure in focusing on specific stages of the disease, and for selecting the appropriate experimental framework to obtain optimal data quality given the purpose of the experiment.

Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study

This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI.

Diffusional Kurtosis Imaging in the Diffusion Imaging in Python Project

Diffusion-weighted magnetic resonance imaging (dMRI) measurements and models provide information about brain connectivity and are sensitive to the physical properties of tissue microstructure. Diffusional Kurtosis Imaging (DKI) quantifies the degree of non-Gaussian diffusion in biological tissue from dMRI. These estimates are of interest because they were shown to be more sensitive to microstructural alterations in health and diseases than measures based on the total anisotropy of diffusion which are highly confounded by tissue dispersion and fiber crossings. In this work, we implemented DKI in the Diffusion in Python (DIPY) project-a large collaborative open-source project which aims to provide well-tested, well-documented and comprehensive implementation of different dMRI techniques. We demonstrate the functionality of our methods in numerical simulations with known ground truth parameters and in openly available datasets. A particular strength of our DKI implementations is that it pursues several extensions of the model that connect it explicitly with microstructural models and the reconstruction of 3D white matter fiber bundles (tractography). For instance, our implementations include DKI-based microstructural models that allow the estimation of biophysical parameters, such as axonal water fraction. Moreover, we illustrate how DKI provides more general characterization of non-Gaussian diffusion compatible with complex white matter fiber architectures and gray matter, and we include a novel mean kurtosis index that is invariant to the confounding effects due to tissue dispersion. In summary, DKI in DIPY provides a well-tested, well-documented and comprehensive reference implementation for DKI. It provides a platform for wider use of DKI in research on brain disorders and in cognitive neuroscience.

Scan-rescan repeatability of axonal imaging metrics using high-gradient diffusion MRI and statistical implications for study design

Axon diameter mapping using diffusion MRI in the living human brain has attracted growing interests with the increasing availability of high gradient strength MRI systems. A systematic assessment of the consistency of axon diameter estimates within and between individuals is needed to gain a comprehensive understanding of how such methods extend to quantifying differences in axon diameter index between groups and facilitate the design of neurobiological studies using such measures. We examined the scan-rescan repeatability of axon diameter index estimation based on the spherical mean technique (SMT) approach using diffusion MRI data acquired with gradient strengths up to 300 mT/m on a 3T Connectom system in 7 healthy volunteers. We performed statistical power analyses using data acquired with the same protocol in a larger cohort consisting of 15 healthy adults to investigate the implications for study design. Results revealed a high degree of repeatability in voxel-wise restricted volume fraction estimates and tract-wise estimates of axon diameter index derived from high-gradient diffusion MRI data. On the region of interest (ROI) level, across white matter tracts in the whole brain, the Pearson’s correlation coefficient of the axon diameter index estimated between scan and rescan experiments was r = 0.72 with an absolute deviation of 0.18 μm. For an anticipated 10% effect size in studies of axon diameter index, most white matter regions required a sample size of less than 15 people to observe a measurable difference between groups using an ROI-based approach. To facilitate the use of high-gradient strength diffusion MRI data for neuroscientific studies of axonal microstructure, the comprehensive multi-gradient strength, multi-diffusion time data used in this work will be made publicly available, in support of open science and increasing the accessibility of such data to the greater scientific community.

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.

The variability of MR axon radii estimates in the human white matter

The noninvasive quantification of axonal morphology is an exciting avenue for gaining understanding of the function and structure of the central nervous system. Accurate non-invasive mapping of micron-sized axon radii using commonly applied neuroimaging techniques, that is, diffusion-weighted MRI, has been bolstered by recent hardware developments, specifically MR gradient design. Here the whole brain characterization of the effective MR axon radius is presented and the inter- and intra-scanner test-retest repeatability and reproducibility are evaluated to promote the further development of the effective MR axon radius as a neuroimaging biomarker. A coefficient-of-variability of approximately 10% in the voxelwise estimation of the effective MR radius is observed in the test-retest analysis, but it is shown that the performance can be improved fourfold using a customized along-tract analysis.

Elevated Intraindividual Variability in Executive Functions and Associations with White Matter Microstructure in Veterans with Mild Traumatic Brain Injury

OBJECTIVE

We examined whether intraindividual variability (IIV) across tests of executive functions (EF-IIV) is elevated in Veterans with a history of mild traumatic brain injury (mTBI) relative to military controls (MCs) without a history of mTBI. We also explored relationships among EF-IIV, white matter microstructure, and posttraumatic stress disorder (PTSD) symptoms.

METHOD

A total of 77 Veterans (mTBI = 43, MCs = 34) completed neuropsychological testing, diffusion tensor imaging (DTI), and PTSD symptom ratings. EF-IIV was calculated as the standard deviation across six tests of EF, along with an EF-Mean composite. DSI Studio connectometry analysis identified white matter tracts significantly associated with EF-IIV according to generalized fractional anisotropy (GFA).

RESULTS

After adjusting for EF-Mean and PTSD symptoms, the mTBI group showed significantly higher EF-IIV than MCs. Groups did not differ on EF-Mean after adjusting for PTSD symptoms. Across groups, PTSD symptoms significantly negatively correlated with EF-Mean, but not with EF-IIV. EF-IIV significantly negatively correlated with GFA in multiple white matter pathways connecting frontal and more posterior regions.

CONCLUSIONS

Veterans with mTBI demonstrated significantly greater IIV across EF tests compared to MCs, even after adjusting for mean group differences on those measures as well as PTSD severity. Findings suggest that, in contrast to analyses that explore effects of mean performance across tests, discrepancy analyses may capture unique variance in neuropsychological performance and more sensitively capture cognitive disruption in Veterans with mTBI histories. Importantly, findings show that EF-IIV is negatively associated with the microstructure of white matter pathways interconnecting cortical regions that mediate executive function and attentional processes.

Effects of Spaceflight Stressors on Brain Volume, Microstructure, and Intracranial Fluid Distribution

Astronauts are exposed to elevated CO2 levels onboard the International Space Station. Here, we investigated structural brain changes in 11 participants following 30-days of head-down tilt bed rest (HDBR) combined with 0.5% ambient CO2 (HDBR + CO2) as a spaceflight analog. We contrasted brain changes observed in the HDBR + CO2 group with those of a previous HDBR sample not exposed to elevated CO2. Both groups exhibited a global upward shift of the brain and concomitant intracranial free water (FW) redistribution. Greater gray matter changes were seen in the HDBR + CO2 group in some regions. The HDBR + CO2 group showed significantly greater FW decrements in the posterior cerebellum and the cerebrum than the HDBR group. In comparison to the HDBR group, the HDBR + CO2 group exhibited greater diffusivity increases. In half of the participants, the HDBR + CO2 intervention resulted in signs of Spaceflight Associated Neuro-ocular Syndrome (SANS), a constellation of ocular structural and functional changes seen in astronauts. We therefore conducted an exploratory comparison compared between subjects that did and did not develop SANS and found asymmetric lateral ventricle enlargement in the SANS group. These results enhance our understanding of the underlying mechanisms of spaceflight-induced brain changes, which is critical for promoting astronaut health and performance.

Beyond the diffusion standard model in fixed rat spinal cord with combined linear and planar encoding

Information about tissue on the microscopic and mesoscopic scales can be accessed by modelling diffusion MRI signals, with the aim of extracting microstructure-specific biomarkers. The standard model (SM) of diffusion, currently the most broadly adopted microstructural model, describes diffusion in white matter (WM) tissues by two Gaussian components, one of which has zero radial diffusivity, to represent diffusion in intra- and extra-axonal water, respectively. Here, we reappraise these SM assumptions by collecting comprehensive double diffusion encoded (DDE) MRI data with both linear and planar encodings, which was recently shown to substantially enhance the ability to estimate SM parameters. We find however, that the SM is unable to account for data recorded in fixed rat spinal cord at an ultrahigh field of 16.4 T, suggesting that its underlying assumptions are violated in our experimental data. We offer three model extensions to mitigate this problem: first, we generalize the SM to accommodate finite radii (axons) by releasing the constraint of zero radial diffusivity in the intra-axonal compartment. Second, we include intracompartmental kurtosis to account for non-Gaussian behaviour. Third, we introduce an additional (third) compartment. The ability of these models to account for our experimental data are compared based on parameter feasibility and Bayesian information criterion. Our analysis identifies the three-compartment description as the optimal model. The third compartment exhibits slow diffusion with a minor but non-negligible signal fraction (∼12%). We demonstrate how failure to take the presence of such a compartment into account severely misguides inferences about WM microstructure. Our findings bear significance for microstructural modelling at large and can impact the interpretation of biomarkers extracted from the standard model of diffusion.

The sensitivity of diffusion MRI to microstructural properties and experimental factors

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

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

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

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.

Axon diameter index estimation independent of fiber orientation distribution using high-gradient diffusion MRI

Axon diameter mapping using high-gradient diffusion MRI has generated great interest as a noninvasive tool for studying trends in axonal size in the human brain. One of the main barriers to mapping axon diameter across the whole brain is accounting for complex white matter fiber configurations (e.g., crossings and fanning), which are prevalent throughout the brain. Here, we present a framework for generalizing axon diameter index estimation to the whole brain independent of the underlying fiber orientation distribution using the spherical mean technique (SMT). This approach is shown to significantly benefit from the use of real-valued diffusion data with Gaussian noise, which reduces the systematic bias in the estimated parameters resulting from the elevation of the noise floor when using magnitude data with Rician noise. We demonstrate the feasibility of obtaining whole-brain orientationally invariant estimates of axon diameter index and relative volume fractions in six healthy human volunteers using real-valued diffusion data acquired on a dedicated high-gradient 3-Tesla human MRI scanner with 300 mT/m maximum gradient strength. The trends in axon diameter index are consistent with known variations in axon diameter from histology and demonstrate the potential of this generalized framework for revealing coherent patterns in axonal structure throughout the living human brain. The use of real-valued diffusion data provides a viable solution for eliminating the Rician noise floor and should be considered for all spherical mean approaches to microstructural parameter estimation.

Direct and specific assessment of axonal injury and spinal cord microenvironments using diffusion correlation imaging

We describe a practical two-dimensional (2D) diffusion MRI framework to deliver specificity and improve sensitivity to axonal injury in the spinal cord. This approach provides intravoxel distributions of correlations of water mobilities in orthogonal directions, revealing sub-voxel diffusion components. Here we use it to investigate water diffusivities along axial and radial orientations within spinal cord specimens with confirmed, tract-specific axonal injury. First, we show using transmission electron microscopy and immunohistochemistry that tract-specific axonal beading occurs following Wallerian degeneration in the cortico-spinal tract as direct sequelae to closed head injury. We demonstrate that although some voxel-averaged diffusion tensor imaging (DTI) metrics are sensitive to this axonal injury, they are non-specific, i.e., they do not reveal an underlying biophysical mechanism of injury. Then we employ 2D diffusion correlation imaging (DCI) to improve discrimination of different water microenvironments by measuring and mapping the joint water mobility distributions perpendicular and parallel to the spinal cord axis. We determine six distinct diffusion spectral components that differ according to their microscopic anisotropy and mobility. We show that at the injury site a highly anisotropic diffusion component completely disappears and instead becomes more isotropic. Based on these findings, an injury-specific MR image of the spinal cord was generated, and a radiological-pathological correlation with histological silver staining % area was performed. The resulting strong and significant correlation (r=0.70,p < 0.0001) indicates the high specificity with which DCI detects injury-induced tissue alterations. We predict that the ability to selectively image microstructural changes following axonal injury in the spinal cord can be useful in clinical and research applications by enabling specific detection and increased sensitivity to injury-induced microstructural alterations. These results also encourage us to translate DCI to higher spatial dimensions to enable assessment of traumatic axonal injury, and possibly other diseases and disorders in the brain.

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.

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

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

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.

Optimization and numerical evaluation of multi-compartment diffusion MRI using the spherical mean technique for practical multiple sclerosis imaging

BACKGROUND

The multi-compartment diffusion MRI using the spherical mean technique (SMT) has been suggested to enhance the pathological specificity to tissue injury in multiple sclerosis (MS) imaging, but its accuracy and precision have not been comprehensively evaluated.

METHODS

A Cramer-Rao Lower Bound method was used to optimize an SMT protocol for MS imaging. Finite difference computer simulations of spins in packed cylinders were then performed to evaluate the influences of five realistic pathological features in MS lesions: axon diameter, axon density, free water fraction, axonal crossing, dispersion, and undulation.

RESULTS

SMT derived metrics can be biased by some confounds of pathological variations, such as axon size and free water fraction. However, SMT in general provides valuable information to characterize pathological features in MS lesions with a clinically feasible protocol.

CONCLUSION

SMT may be used as a practical MS imaging method and should be further improved in clinical MS imaging.

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

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