Simulations on the influence of myelin water in diffusion-weighted imaging

Multi-parametric MRI can detect enhanced myelination in the Gli1 -/- mouse brain

This study investigated the potential of combining multiple MR parameters to enhance the characterization of myelin in the mouse brain. We collected ex vivo multi-parametric MR data at 7 Tesla from control and Gli1 -/- mice; the latter exhibit enhanced myelination at postnatal day 10 (P10) in the corpus callosum and cortex. The MR data included relaxivity, magnetization transfer, and diffusion measurements, each targeting distinct myelin properties. This analysis was followed by and compared to myelin basic protein (MBP) staining of the same samples. Although a majority of the MR parameters included in this study showed significant differences in the corpus callosum between the control and Gli1 -/- mice, only T 2 , T 1 /T 2, and radial diffusivity (RD) demonstrated a significant correlation with MBP values. Based on data from the corpus callosum, partial least square regression suggested that combining T 2 , T 1 /T 2 , and inhomogeneous magnetization transfer ratio could explain approximately 80% of the variance in the MBP values. Myelin predictions based on these three parameters yielded stronger correlations with the MBP values in the P10 mouse brain corpus callosum than any single MR parameter. In the motor cortex, combining T 2 , T 1 /T 2, and radial kurtosis could explain over 90% of the variance in the MBP values at P10. This study demonstrates the utility of multi-parametric MRI in improving the detection of myelin changes in the mouse brain.

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.

T2 relaxation time of the normal-appearing white matter is related to the cognitive status in cerebral small vessel disease

Previous diffusion tensor imaging (DTI) studies indicate that impaired microstructural integrity of the normal-appearing white matter (NAWM) is related to cognitive impairment in cerebral small vessel disease (SVD). This study aimed to investigate whether quantitative T₂ relaxometry is a suitable imaging biomarker for the assessment of tissue changes related to cognitive abnormalities in patients with SVD. 39 patients and 18 age-matched healthy control subjects underwent 3 T magnetic resonance imaging (MRI) with T2-weighted multiple spin echo sequences for T₂ relaxometry and DTI sequences, as well as comprehensive cognitive assessment. Averaged quantitative T₂, fractional anisotropy (FA) and mean diffusivity (MD) were determined in the NAWM and related to cognitive parameters controlling for age, normalized brain volume, white matter hyperintensity volume and other conventional SVD markers. In SVD patients, quantitative T₂ values were significantly increased compared to controls (p = 0.002) and significantly negatively correlated with the global cognitive performance (r= -0.410, p = 0.014) and executive function (r= -0.399, p = 0.016). DTI parameters did not correlate with cognitive function. T₂ relaxometry of the NAWM seems to be sensitive to microstructural tissue damage associated with cognitive impairment in SVD and might be a promising imaging biomarker for evaluation of disease progression and possible effects of therapeutic interventions.

A simple estimate of axon size with diffusion MRI

Noninvasive estimation of mean axon diameter presents a new opportunity to explore white matter plasticity, development, and pathology. Several diffusion-weighted MRI (DW-MRI) methods have been proposed to measure the average axon diameter in white matter, but they typically require many diffusion encoding measurements and complicated mathematical models to fit the signal to multiple tissue compartments, including intra- and extra-axonal spaces. Here, Monte Carlo simulations uncovered a straightforward DW-MRI metric of axon diameter: the change in radial apparent diffusion coefficient estimated at different effective diffusion times, ΔD_⊥. Simulations indicated that this metric increases monotonically within a relevant range of effective mean axon diameter while being insensitive to changes in extra-axonal volume fraction, axon diameter distribution, g-ratio, and influence of myelin water. Also, a monotonic relationship was found to exist for signals coming from both intra- and extra-axonal compartments. The slope in ΔD_⊥ with effective axon diameter increased with the difference in diffusion time of both oscillating and pulsed gradient diffusion sequences.

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 .

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.

Diffusion time dependency along the human corpus callosum and exploration of age and sex differences as assessed by oscillating gradient spin-echo diffusion tensor imaging

Conventional diffusion imaging uses pulsed gradient spin echo (PGSE) waveforms with diffusion times of tens of milliseconds (ms) to infer differences of white matter microstructure. The combined use of these long diffusion times with short diffusion times (<10 ​ms) enabled by oscillating gradient spin echo (OGSE) waveforms can enable more sensitivity to changes of restrictive boundaries on the scale of white matter microstructure (e.g. membranes reflecting the axon diameters). Here, PGSE and OGSE images were acquired at 4.7 ​T from 20 healthy volunteers aged 20-73 years (10 males). Mean, radial, and axial diffusivity, as well as fractional anisotropy were calculated in the genu, body and splenium of the corpus callosum (CC). Monte Carlo simulations were also conducted to examine the relationship of intra- and extra-axonal radial diffusivity with diffusion time over a range of axon diameters and distributions. The results showed elevated diffusivities with OGSE relative to PGSE in the genu and splenium (but not the body) in both males and females, but the OGSE-PGSE difference was greater in the genu for males. Females showed positive correlations of OGSE-PGSE diffusivity difference with age across the CC, whereas there were no such age correlations in males. Simulations of radial diffusion demonstrated that for axon sizes in human brain both OGSE and PGSE diffusivities were dominated by extra-axonal water, but the OGSE-PGSE difference nonetheless increased with area-weighted outer-axon diameter. Therefore, the lack of OGSE-PGSE difference in the body is not entirely consistent with literature that suggests it is composed predominantly of axons with large diameter. The greater OGSE-PGSE difference in the genu of males could reflect larger axon diameters than females. The OGSE-PGSE difference correlation with age in females could reflect loss of smaller axons at older ages. The use of OGSE with short diffusion times to sample the microstructural scale of restriction implies regional differences of axon diameters along the corpus callosum with preliminary results suggesting a dependence on age and sex.

Monte Carlo Simulations of Water Exchange Through Myelin Wraps: Implications for Diffusion MRI

Diffusion magnetic resonance imaging (dMRI) yields parameters sensitive to brain tissue microstructure. A structurally important aspect of this microstructure is the myelin wrapping around the axons. This paper investigated the forward problem concerning whether water exchange via the spiraling structure of the myelin can meaningfully contribute to the signal in dMRI. Monte Carlo simulations were performed in a system with intra-axonal, myelin, and extra-axonal compartments. Diffusion in the myelin was simulated as a spiral wrapping the axon, with a custom number of wraps. Exchange (or intra-axonal residence) times were analyzed for various number of wraps and axon diameters. Pulsed gradient sequences were employed to simulate the dMRI signal, which was analyzed using different methods. Diffusional kurtosis imaging analysis yielded the radial diffusivity (RD) and radial kurtosis (RK), while the two-compartment Kärger model yielded estimates the intra-axonal volume fraction ( ν _ic ) and exchange time ( τ ). Results showed that τ was on the sub-second level for geometries with axon diameters below 1.0 μ m and less than eight wraps. Otherwise, exchange was negligible compared to typical experimental durations, with τ of seconds or longer. In situations where exchange influenced the signal, estimates of RK and ν _ic increased with the number of wraps, while RD decreased. τ estimates from simulated signals were in agreement with predicted ones. In conclusion, exchange through spiraling myelin permits sub-second τ for small diameters and low number of wraps. Such conditions may arise in the developing brain or in neurodegenerative disease, and thus the results could aid the interpretation of dMRI studies.

Standardization of T1w/T2w Ratio Improves Detection of Tissue Damage in Multiple Sclerosis

Normal appearing white matter (NAWM) damage develops early in multiple sclerosis (MS) and continues in the absence of new lesions. The ratio of T1w and T2w (T1w/T2w ratio), a measure of white matter integrity, has previously shown reduced intensity values in MS NAWM. We evaluate the validity of a standardized T1w/T2w ratio (sT1w/T2w ratio) in MS and whether this method is sensitive in detecting MS-related differences in NAWM. T1w and T2w scans were acquired at 3 Tesla in 47 patients with relapsing-remitting MS and 47 matched controls (HC). T1w/T2w and sT1w/T2w ratios were then calculated. We compared between-group variability between T1w/T2w and sT1w/T2w ratio in HC and MS and assessed for group differences. We also evaluated the relationship between the T1w/T2w and sT1w/T2w ratios and clinically relevant variables. Compared to the classic T1w/T2w ratio, the between-subject variability in sT1w/T2w ratio showed a significant reduction in MS patients (p < 0.001) and HC (p < 0.001). However, only sT1w/T2w ratio values were reduced in patients compared to HC (p < 0.001). The sT1w/T2w ratio intensity values were significantly influenced by age, T2 lesion volume and group status (MS vs. HC) (adjusted R² = 0.30, p < 0.001). We demonstrate the validity of the sT1w/T2w ratio in MS and that it is more sensitive to MS-related differences in NAWM compared to T1w/T2w ratio. The sT1w/T2w ratio shows promise as an easily-implemented measure of NAWM in MS using readily available scans and simple post-processing methods.

Challenges in diffusion MRI tractography - Lessons learned from international benchmark competitions

Diffusion MRI (dMRI) fiber tractography has become a pillar of the neuroimaging community due to its ability to noninvasively map the structural connectivity of the brain. Despite widespread use in clinical and research domains, these methods suffer from several potential drawbacks or limitations. Thus, validating the accuracy and reproducibility of techniques is critical for sound scientific conclusions and effective clinical outcomes. Towards this end, a number of international benchmark competitions, or "challenges", has been organized by the diffusion MRI community in order to investigate the reliability of the tractography process by providing a platform to compare algorithms and results in a fair manner, and evaluate common and emerging algorithms in an effort to advance the state of the field. In this paper, we summarize the lessons from a decade of challenges in tractography, and give perspective on the past, present, and future "challenges" that the field of diffusion tractography faces.

Towards microstructure fingerprinting: Estimation of tissue properties from a dictionary of Monte Carlo diffusion MRI simulations

Many closed-form analytical models have been proposed to relate the diffusion-weighted magnetic resonance imaging (DW-MRI) signal to microstructural features of white matter tissues. These models generally make assumptions about the tissue and the diffusion processes which often depart from the biophysical reality, limiting their reliability and interpretability in practice. Monte Carlo simulations of the random walk of water molecules are widely recognized to provide near groundtruth for DW-MRI signals. However, they have mostly been limited to the validation of simpler models rather than used for the estimation of microstructural properties. This work proposes a general framework which leverages Monte Carlo simulations for the estimation of physically interpretable microstructural parameters, both in single and in crossing fascicles of axons. Monte Carlo simulations of DW-MRI signals, or fingerprints, are pre-computed for a large collection of microstructural configurations. At every voxel, the microstructural parameters are estimated by optimizing a sparse combination of these fingerprints. Extensive synthetic experiments showed that our approach achieves accurate and robust estimates in the presence of noise and uncertainties over fixed or input parameters. In an in vivo rat model of spinal cord injury, our approach provided microstructural parameters that showed better correspondence with histology than five closed-form models of the diffusion signal: MMWMD, NODDI, DIAMOND, WMTI and MAPL. On whole-brain in vivo data from the human connectome project (HCP), our method exhibited spatial distributions of apparent axonal radius and axonal density indices in keeping with ex vivo studies. This work paves the way for microstructure fingerprinting with Monte Carlo simulations used directly at the modeling stage and not only as a validation tool.

Physical and numerical phantoms for the validation of brain microstructural MRI: A cookbook

Phantoms, both numerical (software) and physical (hardware), can serve as a gold standard for the validation of MRI methods probing the brain microstructure. This review aims to provide guidelines on how to build, implement, or choose the right phantom for a particular application, along with an overview of the current state-of-the-art of phantoms dedicated to study brain microstructure with MRI. For physical phantoms, we discuss the essential requirements and relevant characteristics of both the (NMR visible) liquid and (NMR invisible) phantom materials that induce relevant microstructural features detectable via MRI, based on diffusion, intra-voxel incoherent motion, magnetization transfer or magnetic susceptibility weighted contrast. In particular, for diffusion MRI, many useful phantoms have been proposed, ranging from simple liquids to advanced biomimetic phantoms consisting of hollow or plain microfibers and capillaries. For numerical phantoms, the focus is on Monte Carlo simulations of random walk, for which the basic principles, along with useful criteria to check and potential pitfalls are reviewed, in addition to a literature overview highlighting recent advances. While many phantoms exist already, the current review aims to stimulate further research in the field and to address remaining needs.

White matter intercompartmental water exchange rates determined from detailed modeling of the myelin sheath

PURPOSE

Magnetization exchange (ME) between hydrogen protons of water and large molecules (semisolids [SS]) in lipid bilayers is an important factor in MRI signal generation and can be exploited to study white matter pathology. Current models used to quantify ME in white matter generally consider water to reside in 1 or 2 distinct compartments, ignoring the complexities of the myelin sheath's multicompartment structure of alternating myelin SS and myelin water (MW) layers. Here, we investigated the effect of this by fitting ME data obtained from human brain at 7 T with a multilayer model of myelin.

METHODS

A multi-echo acquisition for a T₂ ^* -based separation of MW from other water signals was combined with various preparation pulses to change the (relative) state of the SS and water pools and analyzed by fitting with a multilayer exchange model.

RESULTS

The estimated lifetime within a single MW layer was 260 µs, corresponding to a lipid bilayer permeability of 6.7 µm/s. The magnetization lifetime of the aggregate of all MW was estimated at 13 ms, shorter than previously reported values in the range of 40 to 140 ms.

CONCLUSION

Contrary to expectations and previous reports, ME between protons in myelin SS and water is not limited by the myelin sheath but rather by the exchange between SS and water protons. The analysis of ME contrast should account for the relatively short MW lifetime and affects the interpretation of tissue compartmentalization from MRI contrasts such as T₁ - and diffusion-weighting.

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

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

Effect of myelin water exchange on DTI-derived parameters in diffusion MRI: Elucidation of TE dependence

PURPOSE

Water exchange exists between different neuronal compartments of brain tissue but is often ignored in most diffusion models. The goal of the current study was to demonstrate the dependence of diffusion measurements on echo time (TE) in the human brain and to investigate the underlying effects of myelin water exchange.

METHODS

Five healthy subjects were examined with single-shot pulsed-gradient spin-echo echo-planar imaging with fixed duration (δ) and separation (Δ) of diffusion gradient pulses and a set of varying TEs. The effects of water exchange and intrinsic T₂ difference in cellular environments were investigated with Monte Carlo simulations.

RESULTS

Both in vivo measurements and simulations showed that fractional anisotropy (FA) and axial diffusivity (AD) had positive correlations with TE, while radial diffusivity (RD) showed a negative correlation, which is consistent with a previous study. The simulation results further indicated the sensitivity of TE dependence to the change of g-ratio.

CONCLUSION

The exchange between myelin and intra/extra-axonal water pools often plays a non-negligible role in the observed TE dependence of diffusion parameters, which may accompany or alter the effect of intrinsic T₂ in causing such dependence. The TE dependence may potentially serve as a biomarker for demyelination processes (e.g., in multiple sclerosis and Alzheimer's disease). Magn Reson Med 79:1650-1660, 2018. © 2017 International Society for Magnetic Resonance in Medicine.