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

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

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.

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/mm2, while water DDES measurements were performed with a range of b values from 918 to 7199 s/mm2. 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.

Q-space trajectory imaging with positivity constraints (QTI+).

Q-space trajectory imaging (QTI) enables the estimation of useful scalar measures indicative of the local tissue structure. This is accomplished by employing generalized gradient waveforms for diffusion sensitization alongside a diffusion tensor distribution (DTD) model. The first two moments of the underlying DTD are made available by acquisitions at low diffusion sensitivity (b-values). Here, we show that three independent conditions have to be fulfilled by the mean and covariance tensors associated with distributions of symmetric positive semidefinite tensors. We introduce an estimation framework utilizing semi-definite programming (SDP) to guarantee that these conditions are met. Applying the framework on simulated signal profiles for diffusion tensors distributed according to non-central Wishart distributions demonstrates the improved noise resilience of QTI+ over the commonly employed estimation methods. Our findings on a human brain data set also reveal pronounced improvements, especially so for acquisition protocols featuring few number of volumes. Our method's robustness to noise is expected to not only improve the accuracy of the estimates, but also enable a meaningful interpretation of contrast in the derived scalar maps. The technique's performance on shorter acquisitions could make it feasible in routine clinical practice.

Mapping prostatic microscopic anisotropy using linear and spherical b-tensor encoding: A preliminary study.

PURPOSE: Tensor-valued diffusion encoding provides more specific information than conventional diffusion-weighted imaging (DWI), but has mainly been applied in neuroimaging studies. This study aimed to assess its potential for the imaging of prostate cancer (PCa). METHODS: Seventeen patients with histologically proven PCa were enrolled. DWI of the prostate was performed with linear and spherical tensor encoding using a maximal b-value of 1.5 ms/µm2 and a voxel size of 3 × 3 × 4 mm3 . The gamma-distribution model was used to estimate the mean diffusivity (MD), the isotropic kurtosis (MKI ), and the anisotropic kurtosis (MKA ). Regions of interest were placed in MR-defined cancerous tissues, as well as in apparently healthy tissues in the peripheral and transitional zones (PZs and TZs). RESULTS: DWI with linear and spherical encoding yielded different image contrasts at high b-values, which enabled the estimation of MKA and MKI . Compared with healthy tissue (PZs and TZs combined) the cancers displayed a significantly lower MD (P < .05), higher MKI (P < 10-5 ), and lower MKA (P < .05). Compared with the TZ, tissue in the PZ showed lower MD (P < 10-3 ) and higher MKA (P < 10-3 ). No significant differences were found between cancers of different Gleason scores, possibly because of the limited sample size. CONCLUSION: Tensor-valued diffusion encoding enabled mapping of MKA and MKI in the prostate. The elevated MKI in PCa compared with normal tissues suggests an elevated heterogeneity in the cancers. Increased in-plane resolution could improve tumor delineation in future studies.

Cytosolic diffusivity and microscopic anisotropy of N-acetyl aspartate in human white matter with diffusion-weighted MRS at 7 T.

Metabolite diffusion measurable in humans in vivo with diffusion-weighted spectroscopy (DW-MRS) provides a window into the intracellular morphology and state of specific cell types. Anisotropic diffusion in white matter is governed by the microscopic properties of the individual cell types and their structural units (axons, soma, dendrites). However, anisotropy is also markedly affected by the macroscopic orientational distribution over the imaging voxel, particularly in DW-MRS, where the dimensions of the volume of interest (VOI) are much larger than those typically used in diffusion-weighted imaging. One way to address the confound of macroscopic structural features is to average the measurements acquired with uniformly distributed gradient directions to mimic a situation where fibers present in the VOI are orientationally uniformly distributed. This situation allows the extraction of relevant microstructural features such as transverse and longitudinal diffusivities within axons and the related microscopic fractional anisotropy. We present human DW-MRS data acquired at 7 T in two different white matter regions, processed and analyzed as described above, and find that intra-axonal diffusion of the neuronal metabolite N-acetyl aspartate is in good correspondence to simple model interpretations, such as multi-Gaussian diffusion from disperse fibers where the transverse diffusivity can be neglected. We also discuss the implications of our approach for current and future applications of DW-MRS for cell-specific measurements.

Tensor-valued diffusion MRI in under 3 minutes: an initial survey of microscopic anisotropy and tissue heterogeneity in intracranial tumors.

PURPOSE: To evaluate the feasibility of a 3-minutes protocol for assessment of the microscopic anisotropy and tissue heterogeneity based on tensor-valued diffusion MRI in a wide range of intracranial tumors. METHODS: B-tensor encoding was performed in 42 patients with intracranial tumors (gliomas, meningiomas, adenomas, and metastases). Microscopic anisotropy and tissue heterogeneity were evaluated by estimating the anisotropic kurtosis (MKA ) and isotropic kurtosis (MKI ), respectively. An extensive imaging protocol was compared with a 3-minutes protocol. RESULTS: The fast imaging protocol yielded parameters with characteristics in terms of bias and precision similar to the full protocol. Glioblastomas had lower microscopic anisotropy than meningiomas (MKA = 0.29 ± 0.06 vs. 0.45 ± 0.08, P = 0.003). Metastases had higher tissue heterogeneity (MKI = 0.57 ± 0.07) than both the glioblastomas (0.44 ± 0.06, P < 0.001) and meningiomas (0.46 ± 0.06, P = 0.03). CONCLUSION: Evaluation of the microscopic anisotropy and tissue heterogeneity in intracranial tumor patients is feasible in clinically relevant times frames.

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

PURPOSE: Double diffusion encoding (DDE) MRI enables the estimation of microscopic diffusion anisotropy, yielding valuable information on tissue microstructure. A recent study proposed that the acquisition of rotationally invariant DDE metrics, typically obtained using a spherical "5-design," could be greatly simplified by assuming Gaussian diffusion, facilitating reduced acquisition times that are more compatible with clinical settings. Here, we aim to validate the new minimal acquisition scheme against the standard DDE 5-design, and to quantify the proposed method's noise robustness to facilitate future clinical use. THEORY AND METHODS: DDE MRI experiments were performed on both ex vivo and in vivo rat brains at 9.4 T using the 5-design and the proposed minimal design and taking into account the difference in the number of acquisitions. The ensuing microscopic fractional anisotropy (µFA) maps were compared over a range of b-values up to 5000 s/mm2 . Noise robustness was studied using analytical calculations and numerical simulations. RESULTS: The minimal protocol quantified µFA at an accuracy comparable to the estimates obtained by means of the more theoretically robust DDE 5-design. µFA's sensitivity to noise was found to strongly depend on compartment anisotropy and tensor magnitude in a nonlinear manner. When µFA < 0.75 or when mean diffusivity is particularly low, very high signal-to-noise ratio is required for precise quantification of µFA. CONCLUSION: Our work supports using DDE for quantifying microscopic diffusion anisotropy in clinical settings but raises hitherto overlooked precision issues when measuring µFA with DDE and typical clinical signal-to-noise ratio.

In vivo demonstration of microscopic anisotropy in the human kidney using multidimensional diffusion MRI.

PURPOSE: To demonstrate the feasibility of multidimensional diffusion MRI to probe and quantify microscopic fractional anisotropy (µFA) in human kidneys in vivo. METHODS: Linear tensor encoded (LTE) and spherical tensor encoded (STE) renal diffusion MRI scans were performed in 10 healthy volunteers. Respiratory triggering and image registration were used to minimize motion artefacts during the acquisition. Kidney cortex-medulla were semi-automatically segmented based on fractional anisotropy (FA) values. A model-free analysis of LTE and STE signal dependence on b-value in the renal cortex and medulla was performed. Subsequently, µFA was estimated using a single-shell approach. Finally, a comparison of conventional FA and µFA is shown. RESULTS: The hallmark effect of µFA (divergence of LTE and STE signal with increasing b-value) was observed in all subjects. A statistically significant difference between LTE and STE signal was found in the cortex and medulla, starting from b = 750 s/mm2 and b = 500 s/mm2 , respectively. This difference was maximal at the highest b-value sampled (b = 1000 s/mm2 ) which suggests that relatively high b-values are required for µFA mapping in the kidney compared to conventional FA. Cortical and medullary µFA were, respectively, 0.53 ± 0.09 and 0.65 ± 0.05, both respectively higher than conventional FA (0.19 ± 0.02 and 0.40 ± 0.02). CONCLUSION: The feasibility of combining LTE and STE diffusion MRI to probe and quantify µFA in human kidneys is demonstrated for the first time. By doing so, we show that novel microstructure information-not accessible by conventional diffusion encoding-can be probed by multidimensional diffusion MRI. We also identify relevant technical limitations that warrant further development of the technique for body MRI.

Accurate estimation of microscopic diffusion anisotropy and its time dependence in the mouse brain.

Microscopic diffusion anisotropy (µA) has been recently gaining increasing attention for its ability to decouple the average compartment anisotropy from orientation dispersion. Advanced diffusion MRI sequences, such as double diffusion encoding (DDE) and double oscillating diffusion encoding (DODE) have been used for mapping µA, usually using measurements from a single b shell. However, the accuracy of µA estimation vis-à-vis different b-values was not assessed. Moreover, the time-dependence of this metric, which could offer additional insights into tissue microstructure, has not been studied so far. Here, we investigate both these concepts using theory, simulation, and experiments performed at 16.4T in the mouse brain, ex-vivo. In the first part, simulations and experimental results show that the conventional estimation of microscopic anisotropy from the difference of D(O)DE sequences with parallel and orthogonal gradient directions yields values that highly depend on the choice of b-value. To mitigate this undesirable bias, we propose a multi-shell approach that harnesses a polynomial fit of the signal difference up to third order terms in b-value. In simulations, this approach yields more accurate µA metrics, which are similar to the ground-truth values. The second part of this work uses the proposed multi-shell method to estimate the time/frequency dependence of µA. The data shows either an increase or no change in µA with frequency depending on the region of interest, both in white and gray matter. When comparing the experimental results with simulations, it emerges that simple geometric models such as infinite cylinders with either negligible or finite radii cannot replicate the measured trend, and more complex models, which, for example, incorporate structure along the fibre direction are required. Thus, measuring the time dependence of microscopic anisotropy can provide valuable information for characterizing tissue microstructure.

Neurite density imaging versus imaging of microscopic anisotropy in diffusion MRI: A model comparison using spherical tensor encoding.

In diffusion MRI (dMRI), microscopic diffusion anisotropy can be obscured by orientation dispersion. Separation of these properties is of high importance, since it could allow dMRI to non-invasively probe elongated structures such as neurites (axons and dendrites). However, conventional dMRI, based on single diffusion encoding (SDE), entangles microscopic anisotropy and orientation dispersion with intra-voxel variance in isotropic diffusivity. SDE-based methods for estimating microscopic anisotropy, such as the neurite orientation dispersion and density imaging (NODDI) method, must thus rely on model assumptions to disentangle these features. An alternative approach is to directly quantify microscopic anisotropy by the use of variable shape of the b-tensor. Along those lines, we here present the 'constrained diffusional variance decomposition' (CODIVIDE) method, which jointly analyzes data acquired with diffusion encoding applied in a single direction at a time (linear tensor encoding, LTE) and in all directions (spherical tensor encoding, STE). We then contrast the two approaches by comparing neurite density estimated using NODDI with microscopic anisotropy estimated using CODIVIDE. Data were acquired in healthy volunteers and in glioma patients. NODDI and CODIVIDE differed the most in gray matter and in gliomas, where NODDI detected a neurite fraction higher than expected from the level of microscopic diffusion anisotropy found with CODIVIDE. The discrepancies could be explained by the NODDI tortuosity assumption, which enforces a connection between the neurite density and the mean diffusivity of tissue. Our results suggest that this assumption is invalid, which leads to a NODDI neurite density that is inconsistent between LTE and STE data. Using simulations, we demonstrate that the NODDI assumptions result in parameter bias that precludes the use of NODDI to map neurite density. With CODIVIDE, we found high levels of microscopic anisotropy in white matter, intermediate levels in structures such as the thalamus and the putamen, and low levels in the cortex and in gliomas. We conclude that accurate mapping of microscopic anisotropy requires data acquired with variable shape of the b-tensor.

Double oscillating diffusion encoding and sensitivity to microscopic anisotropy.

PURPOSE: To introduce a novel diffusion pulse sequence, namely double oscillating diffusion encoding (DODE), and to investigate whether it adds sensitivity to microscopic diffusion anisotropy (µA) compared to the well-established double diffusion encoding (DDE) methodology. METHODS: We simulate measurements from DODE and DDE sequences for different types of microstructures exhibiting restricted diffusion. First, we compare the effect of varying pulse sequence parameters on the DODE and DDE signal. Then, we analyse the sensitivity of the two sequences to the microstructural parameters (pore diameter and length) which determine µA. Finally, we investigate specificity of measurements to particular substrate configurations. RESULTS: Simulations show that DODE sequences exhibit similar signal dependence on the relative angle between the two gradients as DDE sequences, however, the effect of varying the mixing time is less pronounced. The sensitivity analysis shows that in substrates with elongated pores and various orientations, DODE sequences increase the sensitivity to pore diameter, while DDE sequences are more sensitive to pore length. Moreover, DDE and DODE sequence parameters can be tailored to enhance/suppress the signal from a particular range of substrates. CONCLUSIONS: A combination of DODE and DDE sequences maximize sensitivity to µA, compared to using just the DDE method. Magn Reson Med 78:550-564, 2017. © 2016 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 License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

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

The structural heterogeneity of tumor tissue can be probed by diffusion MRI (dMRI) in terms of the variance of apparent diffusivities within a voxel. However, the link between the diffusional variance and the tissue heterogeneity is not well-established. To investigate this link we test the hypothesis that diffusional variance, caused by microscopic anisotropy and isotropic heterogeneity, is associated with variable cell eccentricity and cell density in brain tumors. We performed dMRI using a novel encoding scheme for diffusional variance decomposition (DIVIDE) in 7 meningiomas and 8 gliomas prior to surgery. The diffusional variance was quantified from dMRI in terms of the total mean kurtosis (MKT), and DIVIDE was used to decompose MKT into components caused by microscopic anisotropy (MKA) and isotropic heterogeneity (MKI). Diffusion anisotropy was evaluated in terms of the fractional anisotropy (FA) and microscopic fractional anisotropy (µFA). Quantitative microscopy was performed on the excised tumor tissue, where structural anisotropy and cell density were quantified by structure tensor analysis and cell nuclei segmentation, respectively. In order to validate the DIVIDE parameters they were correlated to the corresponding parameters derived from microscopy. We found an excellent agreement between the DIVIDE parameters and corresponding microscopy parameters; MKA correlated with cell eccentricity (r=0.95, p<10-7) and MKI with the cell density variance (r=0.83, p<10-3). The diffusion anisotropy correlated with structure tensor anisotropy on the voxel-scale (FA, r=0.80, p<10-3) and microscopic scale (µFA, r=0.93, p<10-6). A multiple regression analysis showed that the conventional MKT parameter reflects both variable cell eccentricity and cell density, and therefore lacks specificity in terms of microstructure characteristics. However, specificity was obtained by decomposing the two contributions; MKA was associated only to cell eccentricity, and MKI only to cell density variance. The variance in meningiomas was caused primarily by microscopic anisotropy (mean±s.d.) MKA=1.11±0.33 vs MKI=0.44±0.20 (p<10-3), whereas in the gliomas, it was mostly caused by isotropic heterogeneity MKI=0.57±0.30 vs MKA=0.26±0.11 (p<0.05). In conclusion, DIVIDE allows non-invasive mapping of parameters that reflect variable cell eccentricity and density. These results constitute convincing evidence that a link exists between specific aspects of tissue heterogeneity and parameters from dMRI. Decomposing effects of microscopic anisotropy and isotropic heterogeneity facilitates an improved interpretation of tumor heterogeneity as well as diffusion anisotropy on both the microscopic and macroscopic scale.

Q-space trajectory imaging for multidimensional diffusion MRI of the human brain.

This work describes a new diffusion MR framework for imaging and modeling of microstructure that we call q-space trajectory imaging (QTI). The QTI framework consists of two parts: encoding and modeling. First we propose q-space trajectory encoding, which uses time-varying gradients to probe a trajectory in q-space, in contrast to traditional pulsed field gradient sequences that attempt to probe a point in q-space. Then we propose a microstructure model, the diffusion tensor distribution (DTD) model, which takes advantage of additional information provided by QTI to estimate a distributional model over diffusion tensors. We show that the QTI framework enables microstructure modeling that is not possible with the traditional pulsed gradient encoding as introduced by Stejskal and Tanner. In our analysis of QTI, we find that the well-known scalar b-value naturally extends to a tensor-valued entity, i.e., a diffusion measurement tensor, which we call the b-tensor. We show that b-tensors of rank 2 or 3 enable estimation of the mean and covariance of the DTD model in terms of a second order tensor (the diffusion tensor) and a fourth order tensor. The QTI framework has been designed to improve discrimination of the sizes, shapes, and orientations of diffusion microenvironments within tissue. We derive rotationally invariant scalar quantities describing intuitive microstructural features including size, shape, and orientation coherence measures. To demonstrate the feasibility of QTI on a clinical scanner, we performed a small pilot study comparing a group of five healthy controls with five patients with schizophrenia. The parameter maps derived from QTI were compared between the groups, and 9 out of the 14 parameters investigated showed differences between groups. The ability to measure and model the distribution of diffusion tensors, rather than a quantity that has already been averaged within a voxel, has the potential to provide a powerful paradigm for the study of complex tissue architecture.

Model-based estimation of microscopic anisotropy using diffusion MRI: a simulation study.

Non-invasive estimation of cell size and shape is a key challenge in diffusion MRI. This article presents a model-based approach that provides independent estimates of pore size and eccentricity from diffusion MRI data. The technique uses a geometric model of finite cylinders with gamma-distributed radii to represent pores of various sizes and elongations. We consider both macroscopically isotropic substrates and substrates of semi-coherently oriented anisotropic pores and we use Monte Carlo simulations to generate synthetic data. We compare the sensitivity of single and double diffusion encoding (SDE and DDE) sequences to the size distribution and eccentricity, and further analyse different protocols of DDE sequences with parallel and/or perpendicular pairs of gradients. We show that explicitly accounting for size distribution is necessary for accurate microstructural parameter estimates, and a model that assumes a single size yields biased eccentricity values. We also find that SDE sequences support estimates, although DDE sequences with mixed parallel and perpendicular gradients enhance accuracy. In the case of macroscopically anisotropic substrates, this model-based approach can be extended to a rotationally invariant framework to provide features of pore shape (specifically eccentricity) in the presence of size distribution and orientation dispersion.

Quantification of microscopic diffusion anisotropy disentangles effects of orientation dispersion from microstructure: applications in healthy volunteers and in brain tumors.

The anisotropy of water diffusion in brain tissue is affected by both disease and development. This change can be detected using diffusion MRI and is often quantified by the fractional anisotropy (FA) derived from diffusion tensor imaging (DTI). Although FA is sensitive to anisotropic cell structures, such as axons, it is also sensitive to their orientation dispersion. This is a major limitation to the use of FA as a biomarker for "tissue integrity", especially in regions of complex microarchitecture. In this work, we seek to circumvent this limitation by disentangling the effects of microscopic diffusion anisotropy from the orientation dispersion. The microscopic fractional anisotropy (µFA) and the order parameter (OP) were calculated from the contrast between signal prepared with directional and isotropic diffusion encoding, where the latter was achieved by magic angle spinning of the q-vector (qMAS). These parameters were quantified in healthy volunteers and in two patients; one patient with meningioma and one with glioblastoma. Finally, we used simulations to elucidate the relation between FA and µFA in various micro-architectures. Generally, µFA was high in the white matter and low in the gray matter. In the white matter, the largest differences between µFA and FA were found in crossing white matter and in interfaces between large white matter tracts, where µFA was high while FA was low. Both tumor types exhibited a low FA, in contrast to the µFA which was high in the meningioma and low in the glioblastoma, indicating that the meningioma contained disordered anisotropic structures, while the glioblastoma did not. This interpretation was confirmed by histological examination. We conclude that FA from DTI reflects both the amount of diffusion anisotropy and orientation dispersion. We suggest that the µFA and OP may complement FA by independently quantifying the microscopic anisotropy and the level of orientation coherence.

Double-pulsed diffusional kurtosis imaging for the in vivo assessment of human brain microstructure.

We have recently extended conventional single-pulsed-field-gradient (s-PFG) diffusional kurtosis imaging (DKI) to double-pulsed-field-gradient (d-PFG) diffusion MRI sequences, with a method known as double-pulsed DKI (DP-DKI). By virtue of a six-dimensional (6D) formulation for q-space, many of the results and insights of s-PFG DKI are generalized to those of DP-DKI. Owing to the fact that DP-DKI isolates the second order contributions to the d-PFG signal (i.e. second order in b-value), the 6D diffusional kurtosis encodes information beyond what is available from s-PFG sequences. Previously, we have demonstrated DP-DKI for in vivo mouse brain at 7 T, and it is the objective of this study to demonstrate the feasibility of DP-DKI at 3 T for the in vivo assessment of human brain microstructure. In addition, an example is given of how to utilize the additional information obtained from DP-DKI for the purpose of biophysical modeling. The relationship between a specific microscopic anisotropy metric estimated from DP-DKI and other recently proposed measures is also discussed.

Detecting compartmental non-Gaussian diffusion with symmetrized double-PFG MRI.

Diffusion in tissue and porous media is known to be non-Gaussian and has been used for clinical indications of stroke and other tissue pathologies. However, when conventional NMR techniques are applied to biological tissues and other heterogeneous materials, the presence of multiple compartments (pores) with different Gaussian diffusivities will also contribute to the measurement of non-Gaussian behavior. Here we present symmetrized double PFG (sd-PFG), which can separate these two contributions to non-Gaussian signal decay as having distinct angular modulation frequencies. In contrast to prior angular d-PFG methods, sd-PFG can unambiguously extract kurtosis as an oscillation from samples with isotropic or uniformly oriented anisotropic pores, and can generally extract a combination of compartmental anisotropy and kurtosis. The method further fixes its sensitivity with respect to the time dependence of the apparent diffusion coefficient. We experimentally demonstrate the measurement of the fourth cumulant (kurtosis) of diffusion and find it consistent with theoretical predictions. By enabling the unambiguous identification of contributions of compartmental kurtosis to the signal, sd-PFG has the potential to help identify the underlying micro-structural changes corresponding to current kurtosis based diagnostics, and act as a novel source of contrast to better resolve tissue micro-structure.

Double-pulsed diffusional kurtosis imaging.

Diffusional kurtosis imaging (DKI) is extended to double-pulsed-field-gradient (d-PFG) diffusion MRI sequences. This gives a practical approach for acquiring and analyzing d-PFG data. In particular, the leading d-PFG effects, beyond what conventional single-pulsed field gradient (s-PFG) provides, are interpreted in terms of the kurtosis for a diffusion displacement probability density function (dPDF) in a six-dimensional (6D) space. The 6D diffusional kurtosis encodes the unique information provided by d-PFG sequences up to second order in the b-value. This observation leads to a compact expression for the signal magnitude, and it suggests novel data acquisition and analysis methods. Double-pulsed DKI (DP-DKI) is demonstrated for in vivo mouse brain with d-PFG data obtained at 7 T. Copyright © 2014 John Wiley & Sons, Ltd.

Orientationally invariant metrics of apparent compartment eccentricity from double pulsed field gradient diffusion experiments.

Pulsed field gradient diffusion sequences (PFG) with multiple diffusion encoding blocks have been indicated to offer new microstructural tissue information, such as the ability to detect nonspherical compartment shapes in macroscopically isotropic samples, i.e. samples with negligible directional signal dependence on diffusion gradients in standard diffusion experiments. However, current acquisition schemes are not rotationally invariant in the sense that the derived metrics depend on the orientation of the sample, and are affected by the interplay of sampling directions and compartment orientation dispersion when applied to macroscopically anisotropic systems. Here we propose a new framework, the d-PFG 5-design, to enable rotationally invariant estimation of double wave vector diffusion metrics (d-PFG). The method is based on the idea that an appropriate orientational average of the signal emulates the signal from a powder preparation of the same sample, where macroscopic anisotropy is absent by construction. Our approach exploits the theory of exact numerical integration (quadrature) of polynomials on the rotation group, and we exemplify the general procedure with a set consisting of 60 pairs of diffusion wave vectors (the d-PFG 5-design) facilitating a theoretically exact determination of the fourth order Taylor or cumulant expansion of the orientationally averaged signal. The d-PFG 5-design is evaluated with numerical simulations and ex vivo high field diffusion MRI experiments in a nonhuman primate brain. Specifically, we demonstrate rotational invariance when estimating compartment eccentricity, which we show offers new microstructural information, complementary to that of fractional anisotropy (FA) from diffusion tensor imaging (DTI). The imaging observations are supported by a new theoretical result, directly relating compartment eccentricity to FA of individual pores.