Validation and noise robustness assessment of microscopic anisotropy estimation with clinically feasible double diffusion encoding MRI

Comparison of two types of microscopic diffusion anisotropy in mouse brain

Two distinct types of microscopic diffusion anisotropy (MA) are compared in brain for both normal control and transgenic (3xTg-AD) mice, which develop Alzheimer's disease pathology. The first type of MA is the commonly used microscopic fractional anisotropy (μFA), and the second is a new MA measure referred to as μFA'. These two MA parameters have different symmetry properties that are central to their physical interpretations. Specifically, μFA is invariant with respect to local rotations of compartmental diffusion tensors while μFA' is invariant with respect to global diffusion tensor deformations. A key distinction between μFA and μFA' is that μFA is affected by the same type of orientationally coherent diffusion anisotropy as the conventional fractional anisotropy (FA) while μFA' is not. Furthermore, μFA can be viewed as having independent contributions from FA and μFA', as is quantified by an equation relating all three anisotropies. The normal control and transgenic mice are studied at ages ranging from 2 to 15 months, with double diffusion encoding MRI being used to estimate μFA and μFA'. μFA and μFA' are nearly identical in low FA brain regions, but they show notable differences when FA is large. In particular, μFA and FA are found to be strongly correlated in the fimbria, but μFA' and FA are not. In addition, both μFA and μFA' are seen to increase with age in the corpus callosum and external capsule, and modest differences between normal control and transgenic mice are observed for μFA and μFA' in the corpus callosum and for μFA in the fimbria. The triad of FA, μFA, and μFA' is proposed as a useful combination of parameters for assessing diffusion anisotropy in brain.

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.

Rapid microscopic fractional anisotropy imaging via an optimized linear regression formulation

Water diffusion anisotropy in the human brain is affected by disease, trauma, and development. Microscopic fractional anisotropy (μFA) is a diffusion MRI (dMRI) metric that can quantify water diffusion anisotropy independent of neuron fiber orientation dispersion. However, there are several different techniques to estimate μFA and few have demonstrated full brain imaging capabilities within clinically viable scan times and resolutions. Here, we present an optimized spherical tensor encoding (STE) technique to acquire μFA directly from the 2nd order cumulant expansion of the powder averaged dMRI signal obtained from direct linear regression (i.e. diffusion kurtosis) which requires fewer powder-averaged signals than other STE fitting techniques and can be rapidly computed. We found that the optimal dMRI parameters for white matter μFA imaging were a maximum b-value of 2000 s/mm² and a ratio of STE to LTE tensor encoded acquisitions of 1.7 for our system specifications. We then compared two implementations of the direct regression approach to the well-established gamma model in 4 healthy volunteers on a 3 Tesla system. One implementation used mean diffusivity (D) obtained from a 2nd order fit of the cumulant expansion, while the other used a linear estimation of D from the low b-values. Both implementations of the direct regression approach showed strong linear correlations with the gamma model (ρ = 0.97 and ρ = 0.90) but mean biases of -0.11 and - 0.02 relative to the gamma model were also observed, respectively. All three μFA measurements showed good test-retest reliability (ρ ≥ 0.79 and bias = 0). To demonstrate the potential scan time advantage of the direct approach, 2 mm isotropic resolution μFA was demonstrated over a 10 cm slab using a subsampled data set with fewer powder-averaged signals that would correspond to a 3.3-min scan. Accordingly, our results introduce an optimization procedure that has enabled nearly full brain μFA in only several minutes.

In vivo tensor-valued diffusion MRI of focal demyelination in white and deep grey matter of rodents

•

We performed in-vivo tensor-valued diffusion MRI in demyelinating rodents.

•

Lysolecithin was injected in white and deep grey matter to cause focal demyelination.

•

Focal demyelination reduced microscopic fractional anisotropy (µFA).

•

Isotropic kurtosis may be particularly sensitive to deep grey matter lesions.

Background

Multiple sclerosis (MS) is a chronic inflammatory demyelinating disease leading to damage of white matter (WM) and grey matter (GM). Magnetic resonance imaging (MRI) is the modality of choice to assess brain damage in MS, but there is an unmet need in MRI for achieving higher sensitivity and specificity to MS-related microstructural alterations in WM and GM.

Objective

To explore whether tensor-valued diffusion MRI (dMRI) can yield sensitive microstructural read-outs for focal demyelination in cerebral WM and deep GM (DGM).

Methods

Eight rats underwent L-α-Lysophosphatidylcholine (LPC) injections in the WM and striatum to introduce focal demyelination. Multimodal MRI was performed at 7 Tesla after 7 days. Tensor-valued dMRI was complemented by diffusion tensor imaging, quantitative MRI and proton magnetic resonance spectroscopy (MRS).

Results

Quantitative MRI and MRS confirmed that LPC injections caused inflammatory demyelinating lesions in WM and DGM. Tensor-valued dMRI illustrated a significant decline of microscopic fractional anisotropy (µFA) in both LPC-treated WM and DGM (P < 0.005) along with a marked increase of isotropic kurtosis (MK_I) in DGM (P < 0.0001).

Conclusion

Tensor-valued dMRI bears considerable potential for microstructural imaging in MS, suggesting a regional µFA decrease may be a sensitive indicator of MS lesions, while a regional MK_I increase may be particularly sensitive in detecting DGM lesions of MS.

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

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

Double diffusion encoding and applications for biomedical imaging

Diffusion Magnetic Resonance Imaging (dMRI) is one of the most important contemporary non-invasive modalities for probing tissue structure at the microscopic scale. The majority of dMRI techniques employ standard single diffusion encoding (SDE) measurements, covering different sequence parameter ranges depending on the complexity of the method. Although many signal representations and biophysical models have been proposed for SDE data, they are intrinsically limited by a lack of specificity. Advanced dMRI methods have been proposed to provide additional microstructural information beyond what can be inferred from SDE. These enhanced contrasts can play important roles in characterizing biological tissues, for instance upon diseases (e.g. neurodegenerative, cancer, stroke), aging, learning, and development. In this review we focus on double diffusion encoding (DDE), which stands out among other advanced acquisitions for its versatility, ability to probe more specific diffusion correlations, and feasibility for preclinical and clinical applications. Various DDE methodologies have been employed to probe compartment sizes (Section 3), decouple the effects of microscopic diffusion anisotropy from orientation dispersion (Section 4), probe displacement correlations, study exchange, or suppress fast diffusing compartments (Section 6). DDE measurements can also be used to improve the robustness of biophysical models (Section 5) and study intra-cellular diffusion via magnetic resonance spectroscopy of metabolites (Section 7). This review discusses all these topics as well as important practical aspects related to the implementation and contrast in preclinical and clinical settings (Section 9) and aims to provide the readers a guide for deciding on the right DDE acquisition for their specific application.

Brain White-Matter Degeneration Due to Aging and Parkinson Disease as Revealed by Double Diffusion Encoding

Microstructure imaging by means of multidimensional diffusion encoding is increasingly applied in clinical research, with expectations that it yields a parameter that better correlates with clinical disability than current methods based on single diffusion encoding. Under the assumption that diffusion within a voxel can be well described by a collection of diffusion tensors, several parameters of this diffusion tensor distribution can be derived, including mean size, variance of sizes, orientational dispersion, and microscopic anisotropy. The information provided by multidimensional diffusion encoding also enables us to decompose the sources of the conventional fractional anisotropy and mean kurtosis. In this study, we explored the utility of the diffusion tensor distribution approach for characterizing white-matter degeneration in aging and in Parkinson disease by using double diffusion encoding. Data from 23 healthy older subjects and 27 patients with Parkinson disease were analyzed. Advanced age was associated with greater mean size and size variances, as well as smaller microscopic anisotropy. By analyzing the parameters underlying diffusion kurtosis, we found that the reductions of kurtosis in aging and Parkinson disease reported in the literature are likely driven by the reduction in microscopic anisotropy. Furthermore, microscopic anisotropy correlated with the severity of motor impairment in the patients with Parkinson disease. The present results support the use of multidimensional diffusion encoding in clinical studies and are encouraging for its future clinical implementation.

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.