The link between diffusion MRI and tumor heterogeneity: Mapping cell eccentricity and density by diffusional variance decomposition (DIVIDE)

An evidence-based approach to assess the accuracy of diffusion kurtosis imaging in characterization of gliomas

OBJECTIVE

Accurate and noninvasive pathologic grading of glioma patients before surgery was crucial to guiding clinicians to select appropriate treatment and improve patient prognosis. This study was performed to investigate the potential diagnostic value of diffusion kurtosis imaging (DKI) to distinguish high-grade gliomas (HGGs) from low-grade gliomas (LGGs) based on an evidence-based approach.

METHODS

Relevant articles that used DKI to distinguish HGG from LGG in Embase, PubMed, China Knowledge Resource Integrated database (CNKI), Web of Knowledge, and Cochrane Libraries databases were electronically searched to April 31, 2018 by 2 reviewers. All analysis was performed by using Meta-disc1.4 and Stata. Influence factors on the diagnostic accuracy were evaluated using meta-regression analysis.

RESULTS

Five eligible studies were included in this meta-analysis. The pooled sensitivity (SEN) and specificity (SPE) was 91% (confidence interval [CI]: 0.78-0.96; P = .02) and 91% (CI: 0.80-0.97; P = .01). The pooled data showed that diagnostic odds ratio (DOR) of DKI was 79.75 (CI: 31.57-201.45). The area under the curve (AUC) of summary receiver operating characteristic curve was 0.96. There is no evidence that our research has a threshold effect (Spearman correlation coefficient: 0.300, P = .624) and publication bias. Meta regression analysis identified that country, language, field strength, and parameter of magnetic resonance imaging had no significant effect on diagnostic performance.

CONCLUSION

The present meta-analysis shows that the mean kurtosis values derived from DKI may be useful in characterization of gliomas with high sensitivity and specificity. Taken into consideration the small sample of this study, we need to be cautious when interpreting the results of this study.

Measuring non-parametric distributions of intravoxel mean diffusivities using a clinical MRI scanner

We measure spectra of water mobilities (i.e., mean diffusivities) from intravoxel pools in brain tissues of healthy subjects with a non-parametric approach. Using a single-shot isotropic diffusion encoding (IDE) preparation, we eliminate signal confounds caused by anisotropic diffusion, including microscopic anisotropy, and acquire in vivo diffusion-weighted images (DWIs) over a wide range of diffusion sensitizations. We analyze the measured IDE signal decays using a regularized inverse laplace transform (ILT) to derive a probability distribution of mean diffusivities of tissue water in each voxel. Based on numerical simulations we assess the sensitivity and accuracy of our ILT analysis and optimize an experimental protocol for use with clinical MRI scanners. In vivo spectra of intravoxel mean diffusivities measured in healthy subjects generally show single-peak distributions throughout the brain parenchyma, with small differences in peak location and shape among white matter, cortical and subcortical gray matter, and cerebrospinal fluid. Mean diffusivity distributions (MDDs) with multiple peaks are observed primarily in voxels at tissue interfaces and are likely due to partial volume contributions. To quantify tissue-specific MDDs with improved statistical power, we average voxel-wise normalized MDDs in corresponding regions-of-interest (ROIs). This non-parametric, rotation-invariant assessment of isotropic diffusivities of tissue water may reflect important microstructural information, such as cell packing and cell size, and active physiological processes, such as water transport and exchange, which may enhance biological specificity in the clinical diagnosis and characterization of ischemic stroke, cancer, neuroinflammation, and neurodegenerative disorders and diseases.

Diffusion MRI of cancer: From low to high b-values

Following its success in early detection of cerebral ischemia, diffusion-weighted imaging (DWI) has been increasingly used in cancer diagnosis and treatment evaluation. These applications are propelled by the rapid development of novel diffusion models to extract biologically valuable information from diffusion-weighted MR signals, and significant advances in MR hardware that has enabled image acquisition with high b-values. This article reviews recent technical developments and clinical applications in cancer imaging using DWI, with a special emphasis on high b-value diffusion models. The article is organized in four sections. First, we provide an overview of diffusion models that are relevant to cancer imaging. The model parameters are discussed in relation to three tissue properties-cellularity, vascularity, and microstructures. An emphasis is placed on characterization of microstructural heterogeneity, given its novelty and close relevance to cancer. Second, we illustrate diffusion MR clinical applications in each of the following three categories: 1) cancer detection and diagnosis; 2) cancer grading, staging, and classification; and 3) cancer treatment response prediction and evaluation. Third, we discuss several practical issues, including selection of image acquisition parameters, reproducibility and reliability, motion management, image distortion, etc., that are commonly encountered when applying DWI to cancer in clinical settings. Lastly, we highlight a few ongoing challenges and provide some possible future directions, particularly in the area of establishing standards via well-organized multicenter clinical trials to accelerate clinical translation of advanced DWI techniques to improving cancer care on a large scale. Level of Evidence: 5 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:23-40.

Comparative Analysis of Diffusional Kurtosis Imaging, Diffusion Tensor Imaging, and Diffusion-Weighted Imaging in Grading and Assessing Cellular Proliferation of Meningiomas

BACKGROUND AND PURPOSE

An accurate evaluation of the World Health Organization grade and cellular proliferation is particularly important in meningiomas. Our aim was to prospectively evaluate and compare diffusional kurtosis imaging, DTI, and DWI metrics in determining the grade and cellular proliferation of meningiomas.

MATERIALS AND METHODS

Ninety-six consecutive patients with histopathologically confirmed meningiomas were included in this study. Mean kurtosis, radial kurtosis, axial kurtosis, fractional anisotropy, mean diffusivity, and ADC were semiautomatically obtained in the solid components of tumors. Each normalized diffusion value was compared between high-grade meningiomas and low-grade meningiomas using the Mann-Whitney U test. Receiver operating characteristic, multiple logistic regression, and Pearson correlation analysis were used for statistical evaluations.

RESULTS

Diffusional kurtosis imaging metrics (mean kurtosis, radial kurtosis, and axial kurtosis) were significantly higher in high-grade meningiomas than in low-grade meningiomas (P ≤ .001). Mean diffusivity and ADC were significantly lower in high-grade meningiomas than in low-grade meningiomas (P = .003 and .002). Mean kurtosis had significantly greater area the under curve values than mean diffusivity and fractional anisotropy in differentiating high-grade meningiomas from low-grade meningiomas (P = .038 and .002). Mean kurtosis was the only variable that could be used to independently differentiate high-grade meningiomas and low-grade meningiomas (P < .001). Significant correlations were found between the Ki-67 labeling index and kurtosis metrics (P < .001), as well as for mean diffusivity and ADC (P = .004, and .007).

CONCLUSIONS

Compared with other diffusion metrics, mean kurtosis may serve as an optimal parameter for evaluating and predicting the meningioma grade. Moreover, diffusion metrics may potentially reflect cellular proliferation.

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

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

Liquid crystal phantom for validation of microscopic diffusion anisotropy measurements on clinical MRI systems

PURPOSE

To develop a phantom for validating MRI pulse sequences and data processing methods to quantify microscopic diffusion anisotropy in the human brain.

METHODS

Using a liquid crystal consisting of water, detergent, and hydrocarbon, we designed a 0.5-L spherical phantom showing the theoretically highest possible degree of microscopic anisotropy. Data were acquired on the Connectome scanner using echo-planar imaging signal readout and diffusion encoding with axisymmetric b-tensors of varying magnitude, anisotropy, and orientation. The mean diffusivity, fractional anisotropy (FA), and microscopic FA (µFA) parameters were estimated.

RESULTS

The phantom was observed to have values of mean diffusivity similar to brain tissue, and relaxation times compatible with echo-planar imaging echo times on the order of 100 ms. The estimated values of µFA were at the theoretical maximum of 1.0, whereas the values of FA spanned the interval from 0.0 to 0.8 as a result of varying orientational order of the anisotropic domains within each voxel.

CONCLUSIONS

The proposed phantom can be manufactured by mixing three widely available chemicals in volumes comparable to a human head. The acquired data are in excellent agreement with theoretical predictions, showing that the phantom is ideal for validating methods for measuring microscopic diffusion anisotropy on clinical MRI systems. Magn Reson Med 79:1817-1828, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution-NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

Magnetic resonance temporal diffusion tensor spectroscopy of disordered anisotropic tissue

Molecular diffusion measured with diffusion weighted MRI (DWI) offers a probe for tissue microstructure. However, inferring microstructural properties from conventional DWI data is a complex inverse problem and has to account for heterogeneity in sizes, shapes and orientations of the tissue compartments contained within an imaging voxel. Alternative experimental means for disentangling the signal signatures of such features could provide a stronger link between the data and its interpretation. Double diffusion encoding (DDE) offers the possibility to factor out variation in compartment shapes from orientational dispersion of anisotropic domains by measuring the correlation between diffusivity in multiple directions. Time dependence of the diffusion is another effect reflecting the dimensions and distributions of barriers. In this paper we extend on DDE with a modified version of the oscillating gradient spin echo (OGSE) experiment, giving a basic contrast mechanism closely linked to both the temporal diffusion spectrum and the compartment anisotropy. We demonstrate our new method on post mortem brain tissue and show that we retrieve the correct temporal diffusion tensor spectrum in synthetic data from Monte Carlo simulations of random walks in a range of disordered geometries of different sizes and shapes.

Multidimensional correlation of nuclear relaxation rates and diffusion tensors for model-free investigations of heterogeneous anisotropic porous materials

Despite their widespread use in non-invasive studies of porous materials, conventional MRI methods yield ambiguous results for microscopically heterogeneous materials such as brain tissue. While the forward link between microstructure and MRI observables is well understood, the inverse problem of separating the signal contributions from different microscopic pores is notoriously difficult. Here, we introduce an experimental protocol where heterogeneity is resolved by establishing 6D correlations between the individual values of isotropic diffusivity, diffusion anisotropy, orientation of the diffusion tensor, and relaxation rates of distinct populations. Such procedure renders the acquired signal highly specific to the sample's microstructure, and allows characterization of the underlying pore space without prior assumptions on the number and nature of distinct microscopic environments. The experimental feasibility of the suggested method is demonstrated on a sample designed to mimic the properties of nerve tissue. If matched to the constraints of whole body scanners, this protocol could allow for the unconstrained determination of the different types of tissue that compose the living human brain.

Efficient experimental designs for isotropic generalized diffusion tensor MRI (IGDTI)

PURPOSE

We propose a new generalized diffusion tensor imaging (GDTI) experimental design and analysis framework for efficiently measuring orientationally averaged diffusion-weighted images (DWIs), which remove bulk signal modulations attributed to diffusion anisotropy and quantify isotropic higher-order diffusion tensors (HOT). We illustrate how this framework accelerates the clinical measurement of rotation-invariant tissue microstructural parameters derived from HOT, such as the HOT-Trace and the mean t-kurtosis.

THEORY AND METHODS

For a large range of b-values, we compare orientationally averaged DWIs measured with high angular resolution diffusion imaging to those obtained with the proposed isotropic GDTI (IGDTI) experimental design. We compare rotation-invariant microstructural parameters measured with IGDTI to those derived from HOTs measured explicitly with GDTI.

RESULTS

In both fixed-brain microimaging and in vivo clinical experiments, IGDTI accurately quantifies mean apparent diffusion coefficient (mADC)-weighted DWIs over a wide range of b-values and allows efficient computation of HOT-derived scalar tissue parameters from a small number of DWIs.

CONCLUSIONS

IGDTI provides direct and accurate estimates of orientationally averaged tissue water mobilities over a wide range of b-values. This efficient method may enable new, sensitive, and quantitative assessments for clinical applications in which changes in mADC can be observe,d such as detecting and characterizing stroke, cancers, and neurodegenerative diseases. Magn Reson Med 79:180-194, 2018. Published 2017. This article is a U.S. Government work and is in the public domain in the USA.

Double diffusion encoding MRI for the clinic

PURPOSE

The purpose of this study is to develop double diffusion encoding (DDE) MRI methods for clinical use. Microscopic diffusion anisotropy measurements from DDE promise greater specificity to changes in tissue microstructure compared with conventional diffusion tensor imaging, but implementation of DDE sequences on whole-body MRI scanners is challenging because of the limited gradient strengths and lengthy acquisition times.

METHODS

A custom single-refocused DDE sequence was implemented on a 3T whole-body scanner. The DDE gradient orientation scheme and sequence parameters were optimized based on a Gaussian diffusion assumption. Using an optimized 5-min DDE acquisition, microscopic fractional anisotropy (μFA) maps were acquired for the first time in multiple sclerosis patients.

RESULTS

Based on simulations and in vivo human measurements, six parallel and six orthogonal diffusion gradient pairs were found to be the minimum number of diffusion gradient pairs necessary to produce a rotationally invariant measurement of μFA. Simulations showed that optimal precision and accuracy of μFA measurements were obtained using b-values between 1500 and 3000 s/mm² . The μFA maps showed improved delineation of multiple sclerosis lesions compared with conventional fractional anisotropy and distinct contrast from T₂ -weighted fluid attenuated inversion recovery and T₁ -weighted imaging.

CONCLUSION

The μFA maps can be measured using DDE in a clinical setting and may provide new opportunities for characterizing multiple sclerosis lesions and other types of tissue degeneration. Magn Reson Med 80:507-520, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Precision and accuracy of diffusion kurtosis estimation and the influence of b-value selection

Diffusion kurtosis imaging (DKI) is an extension of diffusion tensor imaging that accounts for leading non-Gaussian diffusion effects. In DKI studies, a wide range of different gradient strengths (b-values) is used, which is known to affect the estimated diffusivity and kurtosis parameters. Hence there is a need to assess the accuracy and precision of the estimated parameters as a function of b-value. This work examines the error in the estimation of mean of the kurtosis tensor (MKT) with respect to the ground truth, using simulations based on a biophysical model for both gray (GM) and white (WM) matter. Model parameters are derived from densely sampled experimental data acquired in ex vivo rat brain and in vivo human brain. Additionally, the variability of MKT is studied using the experimental data. Prevalent fitting protocols are implemented and investigated. The results show strong dependence on the maximum b-value of both net relative error and standard deviation of error for all of the employed fitting protocols. The choice of b-values with minimum MKT estimation error and standard deviation of error was found to depend on the protocol type and the tissue. Protocols that utilize two terms of the cumulant expansion (DKI) were found to achieve minimum error in GM at b-values less than 1 ms/μm² , whereas maximal b-values of about 2.5 ms/μm² were found to be optimal in WM. Protocols including additional higher order terms of the cumulant expansion were found to provide higher accuracy for the more commonly used b-value regime in GM, but were associated with higher error in WM. Averaged over multiple voxels, a net average error of around 15% for both WM and GM was observed for the optimal b-value choice. These results suggest caution when using DKI generated metrics for microstructural modeling and when comparing results obtained using different fitting techniques and b-values.

Cancer Metabolism and Tumor Heterogeneity: Imaging Perspectives Using MR Imaging and Spectroscopy

Cancer cells reprogram their metabolism to maintain viability via genetic mutations and epigenetic alterations, expressing overall dynamic heterogeneity. The complex relaxation mechanisms of nuclear spins provide unique and convertible tissue contrasts, making magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) pertinent imaging tools in both clinics and research. In this review, we summarized MR methods that visualize tumor characteristics and its metabolic phenotypes on an anatomical, microvascular, microstructural, microenvironmental, and metabolomics scale. The review will progress from the utilities of basic spin-relaxation contrasts in cancer imaging to more advanced imaging methods that measure tumor-distinctive parameters such as perfusion, water diffusion, magnetic susceptibility, oxygenation, acidosis, redox state, and cell death. Analytical methods to assess tumor heterogeneity are also reviewed in brief. Although the clinical utility of tumor heterogeneity from imaging is debatable, the quantification of tumor heterogeneity using functional and metabolic MR images with development of robust analytical methods and improved MR methods may offer more critical roles of tumor heterogeneity data in clinics. MRI/MRS can also provide insightful information on pharmacometabolomics, biomarker discovery, disease diagnosis and prognosis, and treatment response. With these future directions in mind, we anticipate the widespread utilization of these MR-based techniques in studying in vivo cancer biology to better address significant clinical needs.

Separating blood and water: Perfusion and free water elimination from diffusion MRI in the human brain

The assessment of the free water fraction in the brain provides important information about extracellular processes such as atrophy and neuroinflammation in various clinical conditions as well as in normal development and aging. Free water estimates from diffusion MRI are assumed to account for freely diffusing water molecules in the extracellular space, but may be biased by other pools of molecules in rapid random motion, such as the intravoxel incoherent motion (IVIM) of blood, where water molecules perfuse in the randomly oriented capillary network. The goal of this work was to separate the signal contribution of the perfusing blood from that of free-water and of other brain diffusivities. The influence of the vascular compartment on the estimation of the free water fraction and other diffusivities was investigated by simulating perfusion in diffusion MRI data. The perfusion effect in the simulations was significant, especially for the estimation of the free water fraction, and was maintained as long as low b-value data were included in the analysis. Two approaches to reduce the perfusion effect were explored in this study: (i) increasing the minimal b-value used in the fitting, and (ii) using a three-compartment model that explicitly accounts for water molecules in the capillary blood. Estimation of the model parameters while excluding low b-values reduced the perfusion effect but was highly sensitive to noise. The three-compartment model fit was more stable and additionally, provided an estimation of the volume fraction of the capillary blood compartment. The three-compartment model thus disentangles the effects of free water diffusion and perfusion, which is of major clinical importance since changes in these components in the brain may indicate different pathologies, i.e., those originating from the extracellular space, such as neuroinflammation and atrophy, and those related to the vascular space, such as vasodilation, vasoconstriction and capillary density. Diffusion MRI data acquired from a healthy volunteer, using multiple b-shells, demonstrated an expected non-zero contribution from the blood fraction, and indicated that not accounting for the perfusion effect may explain the overestimation of the free water fraction evinced in previous studies. Finally, the applicability of the method was demonstrated with a dataset acquired using a clinically feasible protocol with shorter acquisition time and fewer b-shells.

A test-retest study on Parkinson's PPMI dataset yields statistically significant white matter fascicles

In this work, we propose a diffusion MRI protocol for mining Parkinson's disease diffusion MRI datasets and recover robust disease-specific biomarkers. Using advanced high angular resolution diffusion imaging (HARDI) crossing fiber modeling and tractography robust to partial volume effects, we automatically dissected 50 white matter (WM) fascicles. These fascicles connect deep nuclei (thalamus, putamen, pallidum) to different cortical functional areas (associative, motor, sensorimotor, limbic), basal forebrain and substantia nigra. Then, among these 50 candidate WM fascicles, only the ones that passed a test-retest reproducibility procedure qualified for further tractometry analysis. Leveraging the unique 2-timepoints test-retest Parkinson's Progression Markers Initiative (PPMI) dataset of over 600 subjects, we found statistically significant differences in tract profiles along the subcortico-cortical pathways between Parkinson's disease patients and healthy controls. In particular, significant increases in FA, apparent fiber density, tract-density and generalized FA were detected in some locations of the nigro-subthalamo-putaminal-thalamo-cortical pathway. This connection is one of the major motor circuits balancing the coordination of motor output. Detailed and quantifiable knowledge on WM fascicles in these areas is thus essential to improve the quality and outcome of Deep Brain Stimulation, and to target new WM locations for investigation.

Tensor estimation for double-pulsed diffusional kurtosis imaging

Double-pulsed diffusional kurtosis imaging (DP-DKI) represents the double diffusion encoding (DDE) MRI signal in terms of six-dimensional (6D) diffusion and kurtosis tensors. Here a method for estimating these tensors from experimental data is described. A standard numerical algorithm for tensor estimation from conventional (i.e. single diffusion encoding) diffusional kurtosis imaging (DKI) data is generalized to DP-DKI. This algorithm is based on a weighted least squares (WLS) fit of the signal model to the data combined with constraints designed to minimize unphysical parameter estimates. The numerical algorithm then takes the form of a quadratic programming problem. The principal change required to adapt the conventional DKI fitting algorithm to DP-DKI is replacing the three-dimensional diffusion and kurtosis tensors with the 6D tensors needed for DP-DKI. In this way, the 6D diffusion and kurtosis tensors for DP-DKI can be conveniently estimated from DDE data by using constrained WLS, providing a practical means for condensing DDE measurements into well-defined mathematical constructs that may be useful for interpreting and applying DDE MRI. Data from healthy volunteers for brain are used to demonstrate the DP-DKI tensor estimation algorithm. In particular, representative parametric maps of selected tensor-derived rotational invariants are presented.

Resolution limit of cylinder diameter estimation by diffusion MRI: The impact of gradient waveform and orientation dispersion

Diffusion MRI has been proposed as a non-invasive technique for axonal diameter mapping. However, accurate estimation of small diameters requires strong gradients, which is a challenge for the transition of the technique from preclinical to clinical MRI scanners, since these have weaker gradients. In this work, we develop a framework to estimate the lower bound for accurate diameter estimation, which we refer to as the resolution limit. We analyse only the contribution from the intra-axonal space and assume that axons can be represented by impermeable cylinders. To address the growing interest in using techniques for diffusion encoding that go beyond the conventional single diffusion encoding (SDE) sequence, we present a generalised analysis capable of predicting the resolution limit regardless of the gradient waveform. Using this framework, waveforms were optimised to minimise the resolution limit. The results show that, for parallel cylinders, the SDE experiment is optimal in terms of yielding the lowest possible resolution limit. In the presence of orientation dispersion, diffusion encoding sequences with square-wave oscillating gradients were optimal. The resolution limit for standard clinical MRI scanners (maximum gradient strength 60-80 mT/m) was found to be between 4 and 8 μm, depending on the noise levels and the level of orientation dispersion. For scanners with a maximum gradient strength of 300 mT/m, the limit was reduced to between 2 and 5 μm.

Diffusion Kurtosis Imaging of Gliomas Grades II and III - A Study of Perilesional Tumor Infiltration, Tumor Grades and Subtypes at Clinical Presentation

BACKGROUND

Diffusion kurtosis imaging (DKI) allows for assessment of diffusion influenced by microcellular structures. We analyzed DKI in suspected low-grade gliomas prior to histopathological diagnosis. The aim was to investigate if diffusion parameters in the perilesional normal-appearing white matter (NAWM) differed from contralesional white matter, and to investigate differences between glioma malignancy grades II and III and glioma subtypes (astrocytomas and oligodendrogliomas).

PATIENTS AND METHODS

Forty-eight patients with suspected low-grade glioma were prospectively recruited to this institutional review board-approved study and investigated with preoperative DKI at 3T after written informed consent. Patients with histologically proven glioma grades II or III were further analyzed (n=35). Regions of interest (ROIs) were delineated on T2FLAIR images and co-registered to diffusion MRI parameter maps. Mean DKI data were compared between perilesional and contralesional NAWM (student's t-test for dependent samples, Wilcoxon matched pairs test). Histogram DKI data were compared between glioma types and glioma grades (multiple comparisons of mean ranks for all groups). The discriminating potential for DKI in assessing glioma type and grade was assessed with receiver operating characteristics (ROC) curves.

RESULTS

There were significant differences in all mean DKI variables between perilesional and contralesional NAWM (p=<0.000), except for axial kurtosis (p=0.099). Forty-four histogram variables differed significantly between glioma grades II (n=23) and III (n=12) (p=0.003-0.048) and 10 variables differed significantly between ACs (n=18) and ODs (n=17) (p=0.011-0.050). ROC curves of the best discriminating variables had an area under the curve (AUC) of 0.657-0.815.

CONCLUSIONS

Mean DKI variables in perilesional NAWM differ significantly from contralesional NAWM, suggesting altered microstructure by tumor infiltration not depicted on morphological MRI. Histogram analysis of DKI data identifies differences between glioma grades and subtypes.

Multidimensional diffusion MRI

Principles from multidimensional NMR spectroscopy, and in particular solid-state NMR, have recently been transferred to the field of diffusion MRI, offering non-invasive characterization of heterogeneous anisotropic materials, such as the human brain, at an unprecedented level of detail. Here we revisit the basic physics of solid-state NMR and diffusion MRI to pinpoint the origin of the somewhat unexpected analogy between the two fields, and provide an overview of current diffusion MRI acquisition protocols and data analysis methods to quantify the composition of heterogeneous materials in terms of diffusion tensor distributions with size, shape, and orientation dimensions. While the most advanced methods allow estimation of the complete multidimensional distributions, simpler methods focus on various projections onto lower-dimensional spaces as well as determination of means and variances rather than actual distributions. Even the less advanced methods provide simple and intuitive scalar parameters that are directly related to microstructural features that can be observed in optical microscopy images, e.g. average cell eccentricity, variance of cell density, and orientational order - properties that are inextricably entangled in conventional diffusion MRI. Key to disentangling all these microstructural features is MRI signal acquisition combining isotropic and directional dimensions, just as in the field of multidimensional solid-state NMR from which most of the ideas for the new methods are derived.