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

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.

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.

Reconstruction of ovine axonal cytoarchitecture enables more accurate models of brain biomechanics

There is an increased need and focus to understand how local brain microstructure affects the transport of drug molecules directly administered to the brain tissue, for example in convection-enhanced delivery procedures. This study reports a systematic attempt to characterize the cytoarchitecture of commissural, long association and projection fibres, namely the corpus callosum, the fornix and the corona radiata, with the specific aim to map different regions of the tissue and provide essential information for the development of accurate models of brain biomechanics. Ovine samples are imaged using scanning electron microscopy combined with focused ion beam milling to generate 3D volume reconstructions of the tissue at subcellular spatial resolution. Focus is placed on the characteristic cytological feature of the white matter: the axons and their alignment in the tissue. For each tract, a 3D reconstruction of relatively large volumes, including a significant number of axons, is performed and outer axonal ellipticity, outer axonal cross-sectional area and their relative perimeter are measured. The study of well-resolved microstructural features provides useful insight into the fibrous organization of the tissue, whose micromechanical behaviour is that of a composite material presenting elliptical tortuous tubular axonal structures embedded in the extra-cellular matrix. Drug flow can be captured through microstructurally-based models using 3D volumes, either reconstructed directly from images or generated in silico using parameters extracted from the database of images, leading to a workflow to enable physically-accurate simulations of drug delivery to the targeted tissue.

Imaging and reconstruction of sheep brain axonal cytoarchitecture provides insight for brain biomechanics models that simulate drug delivery and other biological processes governed by interstitial fluid flow and transport.

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

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

Correlation Tensor MRI deciphers underlying kurtosis sources in stroke

Noninvasively detecting and characterizing modulations in cellular scale micro-architecture remains a desideratum for contemporary neuroimaging. Diffusion MRI (dMRI) has become the mainstay methodology for probing microstructure, and, in ischemia, its contrasts have revolutionized stroke management. Diffusion kurtosis imaging (DKI) has been shown to significantly enhance the sensitivity of stroke detection compared to its diffusion tensor imaging (DTI) counterparts. However, the interpretation of DKI remains ambiguous as its contrast may arise from competing kurtosis sources related to the anisotropy of tissue components, diffusivity variance across components, and microscopic kurtosis (e.g., arising from cross-sectional variance, structural disorder, and restriction). Resolving these sources may be fundamental for developing more specific imaging techniques for stroke management, prognosis, and understanding its pathophysiology. In this study, we apply Correlation Tensor MRI (CTI) - a double diffusion encoding (DDE) methodology recently introduced for deciphering kurtosis sources based on the unique information captured in DDE's diffusion correlation tensors - to investigate the underpinnings of kurtosis measurements in acute ischemic lesions. Simulations for the different kurtosis sources revealed specific signatures for cross-sectional variance (representing neurite beading), edema, and cell swelling. Ex vivo CTI experiments at 16.4 T were then performed in an experimental photothrombotic stroke model 3 h post-stroke (N = 10), and successfully separated anisotropic, isotropic, and microscopic non-Gaussian diffusion sources in the ischemic lesions. Each of these kurtosis sources provided unique contrasts in the stroked area. Particularly, microscopic kurtosis was shown to be a primary "driver" of total kurtosis upon ischemia; its large increases, coupled with decreases in anisotropic kurtosis, are consistent with the expected elevation in cross-sectional variance, likely linked to beading effects in small objects such as neurites. In vivo experiments at 9.4 T at the same time point (3 h post ischemia, N = 5) demonstrated the stability and relevance of the findings and showed that fixation is not a dominant confounder in our findings. In future studies, the different CTI contrasts may be useful to address current limitations of stroke imaging, e.g., penumbra characterization, distinguishing lesion progression form tissue recovery, and elucidating pathophysiological correlates.

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.

Comparative analysis of signal models for microscopic fractional anisotropy estimation using q-space trajectory encoding

Microscopic diffusion anisotropy imaging using diffusion-weighted MRI and multidimensional diffusion encoding is a promising method for quantifying clinically and scientifically relevant microstructural properties of neural tissue. Several methods for estimating microscopic fractional anisotropy (µFA), a normalized measure of microscopic diffusion anisotropy, have been introduced but the differences between the methods have received little attention thus far. In this study, the accuracy and precision of µFA estimation using q-space trajectory encoding and different signal models were assessed using imaging experiments and simulations. Three healthy volunteers and a microfibre phantom were imaged with five non-zero b-values and gradient waveforms encoding linear and spherical b-tensors. Since the ground-truth µFA was unknown in the imaging experiments, Monte Carlo random walk simulations were performed using axon-mimicking fibres for which the ground truth was known. Furthermore, parameter bias due to time-dependent diffusion was quantified by repeating the simulations with tuned waveforms, which have similar power spectra, and with triple diffusion encoding, which, unlike q-space trajectory encoding, is not based on the assumption of time-independent diffusion. The truncated cumulant expansion of the powder-averaged signal, gamma-distributed diffusivities assumption, and q-space trajectory imaging, a generalization of the truncated cumulant expansion to individual signals, were used to estimate µFA. The gamma-distributed diffusivities assumption consistently resulted in greater µFA values than the second order cumulant expansion, 0.1 greater when averaged over the whole brain. In the simulations, the generalized cumulant expansion provided the most accurate estimates. Importantly, although time-dependent diffusion caused significant overestimation of µFA using all the studied methods, the simulations suggest that the resulting bias in µFA is less than 0.1 in human white matter.

Evidence for microscopic kurtosis in neural tissue revealed by correlation tensor MRI

PURPOSE

The impact of microscopic diffusional kurtosis (µK), arising from restricted diffusion and/or structural disorder, remains a controversial issue in contemporary diffusion MRI (dMRI). Recently, correlation tensor imaging (CTI) was introduced to disentangle the sources contributing to diffusional kurtosis, without relying on a-priori multi-gaussian component (MGC) or other microstructural assumptions. Here, we investigated µK in in vivo rat brains and assessed its impact on state-of-the-art methods ignoring µK.

THEORY AND METHODS

CTI harnesses double diffusion encoding (DDE) experiments, which were here improved for speed and minimal bias using four different sets of acquisition parameters. The robustness of the improved CTI protocol was assessed via simulations. In vivo CTI acquisitions were performed in healthy rat brains using a 9.4T pre-clinical scanner equipped with a cryogenic coil, and targeted the estimation of µK, anisotropic kurtosis, and isotropic kurtosis.

RESULTS

The improved CTI acquisition scheme substantially reduces scan time and importantly, also minimizes higher-order-term biases, thus enabling robust µK estimation, alongside K_aniso and K_iso metrics. Our CTI experiments revealed positive µK both in white and gray matter of the rat brain in vivo; µK is the dominant kurtosis source in healthy gray matter tissue. The non-negligible µK substantially were found to bias prior MGC analyses of K_iso and K_aniso .

CONCLUSIONS

Correlation Tensor MRI offers a more accurate and robust characterization of kurtosis sources than its predecessors. µK is non-negligible in vivo in healthy white and gray matter tissues and could be an important biomarker for future studies. Our findings thus have both theoretical and practical implications for future dMRI research.

Design and characterization of a 3D-printed axon-mimetic phantom for diffusion MRI

PURPOSE

To introduce and characterize inexpensive and easily produced 3D-printed axon-mimetic diffusion MRI phantoms in terms of pore geometry and diffusion kurtosis imaging metrics.

METHODS

Phantoms were 3D-printed with a composite printing material that, after the dissolution of the polyvinyl alcohol, exhibits microscopic fibrous pores. Confocal microscopy and synchrotron phase-contrast micro-CT imaging were performed to visualize and assess the pore sizes. Diffusion MRI scans of four identical phantoms and phantoms with varying print parameters in water were performed at 9.4 T. Diffusion kurtosis imaging was fit to both data sets and used to assess the reproducibility between phantoms and effects of print parameters on diffusion kurtosis imaging metrics. Identical scans were performed 25 and 76 days later, to test their stability.

RESULTS

Segmentation of pores in three microscopy images yielded a mean, median, and SD of equivalent pore diameters of 7.57 μm, 3.51 μm, and 12.13 μm, respectively. Phantoms had T₁ /T₂ = 2 seconds/180 ms, and those with identical parameters showed a low coefficient of variation (~10%) in mean diffusivity (1.38 × 10^-3 mm² /s) and kurtosis (0.52) metrics and radial diffusivity (1.01 × 10^-3 mm² /s) and kurtosis (1.13) metrics. Printing temperature and speed had a small effect on diffusion kurtosis imaging metrics (< 16%), whereas infill density had a larger and more variable effect (> 16%). The stability analysis showed small changes over 2.5 months (< 7%).

CONCLUSION

Three-dimension-printed axon-mimetic phantoms can mimic the fibrous structure of axon bundles on a microscopic scale, serving as complex, anisotropic diffusion MRI phantoms.

Recent Advances in Parameter Inference for Diffusion MRI Signal Models

In this paper, fundamentals and recent progress for obtaining biological features quantitatively by using diffusion MRI are reviewed. First, a brief description of diffusion MRI history, application, and development was presented. Then, well-known parametric models including diffusion tensor imaging (DTI), diffusional kurtosis imaging (DKI), and neurite orientation dispersion diffusion imaging (NODDI) are introduced with several classifications in various viewpoints with other modeling schemes. In addition, this review covers mathematical generalization and examples of methodologies for the model parameter inference from conventional fitting to recent machine learning approaches, which is called Q-space learning (QSL). Finally, future perspectives on diffusion MRI parameter inference are discussed with the aspects of imaging modeling and simulation.

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.

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.