Compensation of concomitant field effects in double diffusion encoding by means of added oscillating gradients

Maxwell or concomitant fields imprint additional phases on the transverse magnetization. This concomitant phase may cause severe image artifacts like signal voids or distort the quantitative parameters due to the induced intravoxel dephasing. In particular, double diffusion encoding (DDE) schemes with two pairs of bipolar diffusion-weighting gradients separated by a refocusing radiofrequency (RF) pulse are prone to concomitant field-induced artifacts. In this work, a method for reducing concomitant field effects in these DDE sequences based on additional oscillating gradients is presented. These oscillating gradient pulses obtained by constrained optimization were added to the original gradient waveforms. The modified sequences reduced the accumulated concomitant phase without significant changes in the original sequence characteristics. The proposed method was applied to a DDE acquisition scheme consisting of 60 pairs of diffusion wave vectors. For phantom as well as for in vivo experiments, a considerable increase in the signal-to-noise ratio (SNR) was obtained. For phantom measurements with a diffusion weighting of b = 2000 s/mm2 for each of the gradient pairs, an SNR increase of up to 40% was observed for a transversal slice that had a distance of 5 cm from the isocenter. For equivalent slice parameters, in vivo measurements in the brain of a healthy volunteer exhibited an increase in SNR of up to 35% for b = 750 s/mm2 for each weighting. These findings are supported by corresponding simulations, which also predict a positive effect on the SNR. In summary, the presented method leads to an SNR gain without additional RF refocusing pulses.

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.

Double pulsed field gradient diffusion MRI to assess skeletal muscle microstructure

PURPOSE

Preliminary study to determine whether double pulsed field gradient (PFG) diffusion MRI is sensitive to key features of muscle microstructure related to function.

METHODS

The restricted diffusion profile of molecules in models of muscle microstructure derived from histology were systematically simulated using a numerical simulation approach. Diffusion tensor subspace imaging analysis of the diffusion signal was performed, and spherical anisotropy (SA) was calculated for each model. Linear regression was used to determine the predictive capacity of SA on the fiber area, fiber diameter, and surface area to volume ratio of the models. Additionally, a rat model of muscle hypertrophy was scanned using a single PFG and a double PFG pulse sequence, and the restricted diffusion measurements were compared with histological measurements of microstructure.

RESULTS

Excellent agreement between SA and muscle fiber area (r2  = 0.71; p < 0.0001), fiber diameter (r2  = 0.83; p < 0.0001), and surface area to volume ratio (r2  = 0.97; p < 0.0001) in simulated models was found. In a scanned rat leg, the distribution of these microstructural features measured from histology was broad and demonstrated that there is a wide variance in the microstructural features observed, similar to the SA distributions. However, the distribution of fractional anisotropy measurements in the same tissue was narrow.

CONCLUSIONS

This study demonstrates that SA-a scalar value from diffusion tensor subspace imaging analysis-is highly sensitive to muscle microstructural features predictive of function. Furthermore, these techniques and analysis tools can be translated to real experiments in skeletal muscle. The increased dynamic range of SA compared with fractional anisotropy in the same tissue suggests increased sensitivity to detecting changes in tissue microstructure.

Extra-axonal contribution to double diffusion encoding-based pore size estimates in the corticospinal tract

OBJECTIVE

To study the origin of compartment size overestimation in double diffusion encoding MRI (DDE) in vivo experiments in the human corticospinal tract. Here, the extracellular space is hypothesized to be the origin of the DDE signal. By exploiting the DDE sensitivity to pore shape, it could be possible to identify the origin of the measured signal. The signal difference between parallel and perpendicular diffusion gradient orientation can indicate if a compartment is regular or eccentric in shape. As extracellular space can be considered an eccentric compartment, a positive difference would mean a high contribution to the compartment size estimates.

MATERIALS AND METHODS

Computer simulations using MISST and in vivo experiments in eight healthy volunteers were performed. DDE experiments using a double spin-echo preparation with eight perpendicular directions were measured in vivo. The difference between parallel and perpendicular gradient orientations was analyzed using a Wilcoxon signed-rank test and a Mann-Whitney U test.

RESULTS

Simulations and MR experiments showed a statistically significant difference between parallel and perpendicular diffusion gradient orientation signals ([Formula: see text]).

CONCLUSION

The results suggest that the DDE-based size estimate may be considerably influenced by the extra-axonal compartment. However, the experimental results are also consistent with purely intra-axonal contributions in combination with a large fiber orientation dispersion.

MR Fingerprinting with b-Tensor Encoding for Simultaneous Quantification of Relaxation and Diffusion in a Single Scan

PURPOSE

Although both relaxation and diffusion imaging are sensitive to tissue microstructure, studies have reported limited sensitivity and robustness of using relaxation or conventional diffusion alone to characterize tissue microstructure. Recently, it has been shown that tensor-valued diffusion encoding and joint relaxation-diffusion quantification enable more reliable quantification of compartment-specific microstructural properties. However, scan times to acquire such data can be prohibitive. Here, we aim to simultaneously quantify relaxation and diffusion using MR fingerprinting (MRF) and b-tensor encoding in a clinically feasible time.

METHODS

We developed multidimensional MRF scans (mdMRF) with linear and spherical b-tensor encoding (LTE and STE) to simultaneously quantify T1, T2, and ADC maps from a single scan. The image quality, accuracy, and scan efficiency were compared between the mdMRF using LTE and STE. Moreover, we investigated the robustness of different sequence designs to signal errors and their impact on the maps.

RESULTS

T1 and T2 maps derived from the mdMRF scans have consistently high image quality, while ADC maps are sensitive to different sequence designs. Notably, the fast imaging steady state precession (FISP)-based mdMRF scan with peripheral pulse gating provides the best ADC maps that are free of image distortion and shading artifacts.

CONCLUSION

We demonstrated the feasibility of quantifying T1, T2, and ADC maps simultaneously from a single mdMRF scan in around 24 s/slice. The map quality and quantitative values are consistent with the reference scans.

Nonparametric 5D D-R2 distribution imaging with single-shot EPI at 21.1 T: Initial results for in vivo rat brain

In vivo human diffusion MRI is by default performed using single-shot EPI with greater than 50-ms echo times and associated signal loss from transverse relaxation. The individual benefits of the current trends of increasing B₀ to boost SNR and employing more advanced signal preparation schemes to improve the specificity for selected microstructural properties eventually may be cancelled by increased relaxation rates at high B₀ and echo times with advanced encoding. Here, initial attempts to translate state-of-the-art diffusion-relaxation correlation methods from 3 T to 21.1 T are made to identify hurdles that need to be overcome to fulfill the promises of both high SNR and readily interpretable microstructural information.

A review on investigation of the basic contrast mechanism underlying multidimensional diffusion MRI in assessment of neurological disorders

INTRODUCTION

Multidimensional diffusion MRI (MDD MRI) is a novel diffusion technique that uses advanced gradient waveforms for microstructural tissue characterization to provide information about average rate, anisotropy and orientation of the diffusion and to disentangle the signal fraction from specific cell types i.e., elongated cells, isotropic cells and free water.

AIM

To review the diagnostic potential of MDD MRI in the clinical setting for microstructural tissue characterization in patients with neurological disorders to aid in patient care and treatment.

METHOD

A scoping review on the clinical applications of MDD MRI was conducted from original articles published in PubMed and Scopus from 2015 to 2021 using the keywords "Multidimensional diffusion MRI" OR "diffusion tensor distribution" OR "Tensor-Valued Diffusion" OR "b-tensor encoding" OR "microscopic diffusion anisotropy" OR "microscopic anisotropy" OR "microscopic fractional anisotropy" OR "double diffusion encoding" OR "triple diffusion encoding" OR "double pulsed field gradients" OR "double wave vector" OR "correlation tensor imaging" AND "brain" OR "axons".

RESULTS

Initially 145 articles were screened and after applying inclusion and exclusion criteria, nine articles were included in the final analysis. In most of these studies, microscopic diffusion anisotropy within the lesion showed deviation from the normal-appearing tissue.

CONCLUSION

Multidimensional diffusion MRI can provide better quantification and visualization of tissue microstructure than conventional diffusion MRI and can be used in the clinical setting for diagnosis of neurological disorders.

Separating Glioma Hyperintensities From White Matter by Diffusion-Weighted Imaging With Spherical Tensor Encoding

Background

Tumor-related hyperintensities in high b-value diffusion-weighted imaging (DWI) are radiologically important in the workup of gliomas. However, the white matter may also appear as hyperintense, which may conflate interpretation.

Purpose

To investigate whether DWI with spherical b-tensor encoding (STE) can be used to suppress white matter and enhance the conspicuity of glioma hyperintensities unrelated to white matter.

Materials and Methods

Twenty-five patients with a glioma tumor and at least one pathology-related hyperintensity on DWI underwent conventional MRI at 3 T. The DWI was performed both with linear and spherical tensor encoding (LTE-DWI and STE-DWI). The LTE-DWI here refers to the DWI obtained with conventional diffusion encoding and averaged across diffusion-encoding directions. Retrospectively, the differences in contrast between LTE-DWI and STE-DWI, obtained at a b-value of 2,000 s/mm², were evaluated by comparing hyperintensities and contralateral normal-appearing white matter (NAWM) both visually and quantitatively in terms of the signal intensity ratio (SIR) and contrast-to-noise ratio efficiency (CNR_eff).

Results

The spherical tensor encoding DWI was more effective than LTE-DWI at suppressing signals from white matter and improved conspicuity of pathology-related hyperintensities. The median SIR improved in all cases and on average by 28%. The median (interquartile range) SIR was 1.9 (1.6 – 2.1) for STE and 1.4 (1.3 – 1.7) for LTE, with a significant difference of 0.4 (0.3 –0.5) (p < 10^–4, paired U-test). In 40% of the patients, the SIR was above 2 for STE-DWI, but with LTE-DWI, the SIR was below 2 for all patients. The CNR_eff of STE-DWI was significantly higher than of LTE-DWI: 2.5 (2 – 3.5) vs. 2.3 (1.7 – 3.1), with a significant difference of 0.4 (−0.1 –0.6) (p < 10^–3, paired U-test). The STE improved CNR_eff in 70% of the cases. We illustrate the benefits of STE-DWI in three patients, where STE-DWI may facilitate an improved radiological description of tumor-related hyperintensity, including one case that could have been missed out if only LTE-DWI was inspected.

Conclusion

The contrast mechanism of high b-value STE-DWI results in a stronger suppression of white matter than conventional LTE-DWI, and may, therefore, be more sensitive and specific for assessment of glioma tumors and DWI-hyperintensities.

In vivo Correlation Tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner

Diffusion MRI (dMRI) has become one of the most important imaging modalities for noninvasively probing tissue microstructure. Diffusional Kurtosis MRI (DKI) quantifies the degree of non-gaussian diffusion, which in turn has been shown to increase sensitivity towards, e.g., disease and orientation mapping in neural tissue. However, the specificity of DKI is limited as different sources can contribute to the total intravoxel diffusional kurtosis, including: variance in diffusion tensor magnitudes (K_iso), variance due to diffusion anisotropy (K_aniso), and microscopic kurtosis (μK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, μK is typically ignored in diffusion MRI signal modeling as it is assumed to be negligible in neural tissues. However, recently, Correlation Tensor MRI (CTI) based on Double-Diffusion-Encoding (DDE) was introduced for kurtosis source separation, revealing non negligible μK in preclinical imaging. Here, we implemented CTI for the first time on a clinical 3T scanner and investigated the sources of total kurtosis in healthy subjects. A robust framework for kurtosis source separation in humans is introduced, followed by estimation of μK (and the other kurtosis sources) in the healthy brain. Using this clinical CTI approach, we find that μK significantly contributes to total diffusional kurtosis both in gray and white matter tissue but, as expected, not in the ventricles. The first μK maps of the human brain are presented, revealing that the spatial distribution of μK provides a unique source of contrast, appearing different from isotropic and anisotropic kurtosis counterparts. Moreover, group average templates of these kurtosis sources have been generated for the first time, which corroborated our findings at the underlying individual-level maps. We further show that the common practice of ignoring μK and assuming the multiple gaussian component approximation for kurtosis source estimation introduces significant bias in the estimation of other kurtosis sources and, perhaps even worse, compromises their interpretation. Finally, a twofold acceleration of CTI is discussed in the context of potential future clinical applications. We conclude that CTI has much potential for future in vivo microstructural characterizations in healthy and pathological tissue.

Correlation Tensor MRI deciphers underlying kurtosis sources in stroke

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

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.

Motion-compensated gradient waveforms for tensor-valued diffusion encoding by constrained numerical optimization

PURPOSE

Diffusion-weighted MRI is sensitive to incoherent tissue motion, which may confound the measured signal and subsequent analysis. We propose a "motion-compensated" gradient waveform design for tensor-valued diffusion encoding that negates the effects bulk motion and incoherent motion in the ballistic regime.

METHODS

Motion compensation was achieved by constraining the magnitude of gradient waveform moment vectors. The constraint was incorporated into a numerical optimization framework, along with existing constraints that account for b-tensor shape, hardware restrictions, and concomitant field gradients. We evaluated the efficacy of encoding and motion compensation in simulations, and we demonstrated the approach by linear and planar b-tensor encoding in a healthy heart in vivo.

RESULTS

The optimization framework produced asymmetric motion-compensated waveforms that yielded b-tensors of arbitrary shape with improved efficiency compared with previous designs for tensor-valued encoding, and equivalent efficiency to previous designs for linear (conventional) encoding. Technical feasibility was demonstrated in the heart in vivo, showing vastly improved data quality when using motion compensation. The optimization framework is available online in open source.

CONCLUSION

Our gradient waveform design is both more flexible and efficient than previous methods, facilitating tensor-valued diffusion encoding in tissues in which motion would otherwise confound the signal. The proposed design exploits asymmetric encoding times, a single refocusing pulse or multiple refocusing pulses, and integrates compensation for concomitant gradient effects throughout the imaging volume.

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.

Microstructure with diffusion MRI: what scale we are sensitive to?

Diffusion-weighted MRI is the forerunner of the rapidly developed microstructure MRI (μMRI) aimed at in vivo evaluation of the cellular tissue architecture. This brief review focuses on the spatiotemporal scales of the microstructure that are accessible using different diffusion MRI techniques and the need to weight the measurability against the interpretability of results. Diffusion phenomena and models are first classified in two-dimensional space (the q-t-plane) of the measurement with narrow gradient pulses. Three-dimensional parameter space of the Stejskal-Tanner diffusion weighting adds more phenomena to this collection. Modern measurement techniques with larger number of parameters are briefly discussed under the overarching idea of diffusion weighting matching the geometry of the targeted cell species.

Tensor-valued diffusion MRI differentiates cortex and white matter in malformations of cortical development associated with epilepsy

Objective:

Delineation of malformations of cortical development (MCD) is central in presurgical evaluation of drug-resistant epilepsy. Delineation using magnetic resonance imaging (MRI) can be ambiguous, however, because the conventional T₁- and T₂-weighted contrasts depend strongly on myelin for differentiation of cortical tissue and white matter. Variations in myelin content within both cortex and white matter may cause MCD findings on MRI to change size, become undetectable, or disagree with histopathology. The novel tensor-valued diffusion MRI (dMRI) technique maps microscopic diffusion anisotropy, which is sensitive to axons rather than myelin. This work investigated whether tensor-valued dMRI may improve differentiation of cortex and white matter in the delineation of MCD.

Methods:

Tensor-valued dMRI was performed on a 7 T MRI scanner in 13 MCD patients (age = 32 ± 13 years) featuring periventricular heterotopia, subcortical heterotopia, focal cortical dysplasia, and polymicrogyria. Data analysis yielded maps of microscopic anisotropy that were compared with T₁-weighted and T₂-fluid-attenuated inversion recovery images and with the fractional anisotropy from diffusion tensor imaging.

Results:

Maps of microscopic anisotropy revealed large white matter-like regions within MCD that were uniformly cortex-like in the conventional MRI contrasts. These regions were seen particularly in the deep white matter parts of subcortical heterotopias and near the gray-white boundaries of focal cortical dysplasias and polymicrogyrias.

Significance:

By being sensitive to axons rather than myelin, mapping of microscopic anisotropy may yield a more robust differentiation of cortex and white matter and improve MCD delineation in presurgical evaluation of epilepsy.

JEDI: Joint Estimation Diffusion Imaging of macroscopic and microscopic tissue properties

PURPOSE

A new method for enhancing the sensitivity of diffusion MRI (dMRI) by combining the data from single (sPFG) and double (dPFG) pulsed field gradient experiments is presented.

METHODS

This method uses our JESTER framework to combine microscopic anisotropy information from dFPG experiments using a new method called diffusion tensor subspace imaging (DiTSI) to augment the macroscopic anisotropy information from sPFG data analyzed using our guided by entropy spectrum pathways method. This new method, called joint estimation diffusion imaging (JEDI), combines the sensitivity to macroscopic diffusion anisotropy of sPFG with the sensitivity to microscopic diffusion anisotropy of dPFG methods.

RESULTS

Its ability to produce significantly more detailed anisotropy maps and more complete fiber tracts than existing methods within both brain white matter (WM) and gray matter (GM) is demonstrated on normal human subjects on data collected using a novel fast, robust, and clinically feasible sPFG/dPFG acquisition.

CONCLUSIONS

The potential utility of this method is suggested by an initial demonstration of its ability to mitigate the problem of gyral bias. The capability of more completely characterizing the tissue structure and connectivity throughout the entire brain has broad implications for the utility and scope of dMRI in a wide range of research and clinical applications.

Single- and double-Diffusion encoding MRI for studying ex vivo apparent axon diameter distribution in spinal cord white matter

Mapping average axon diameter (AAD) and axon diameter distribution (ADD) in neuronal tissues non-invasively is a challenging task that may have a tremendous effect on our understanding of the normal and diseased central nervous system (CNS). Water diffusion is used to probe microstructure in neuronal tissues, however, the different water populations and barriers that are present in these tissues turn this into a complex task. Therefore, it is not surprising that recently we have witnessed a burst in the development of new approaches and models that attempt to obtain, non-invasively, detailed microstructural information in the CNS. In this work, we aim at challenging and comparing the microstructural information obtained from single diffusion encoding (SDE) with double diffusion encoding (DDE) MRI. We first applied SDE and DDE MR spectroscopy (MRS) on microcapillary phantoms and then applied SDE and DDE MRI on an ex vivo porcine spinal cord (SC), using similar experimental conditions. The obtained diffusion MRI data were fitted by the same theoretical model, assuming that the signal in every voxel can be approximated as the superposition of a Gaussian-diffusing component and a series of restricted components having infinite cylindrical geometries. The diffusion MRI results were then compared with histological findings. We found a good agreement between the fittings and the experimental data in white matter (WM) voxels of the SC in both diffusion MRI methods. The microstructural information and apparent AADs extracted from SDE MRI were found to be similar or somewhat larger than those extracted from DDE MRI especially when the diffusion time was set to 40 ms. The apparent ADDs extracted from SDE and DDE MRI show reasonable agreement but somewhat weaker correspondence was observed between the diffusion MRI results and histology. The apparent subtle differences between the microstructural information obtained from SDE and DDE MRI are briefly discussed.

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.

Searching for the neurite density with diffusion MRI: Challenges for biophysical modeling

In vivo mapping of the neurite density with diffusion MRI (dMRI) is a high but challenging aim. First, it is unknown whether all neurites exhibit completely anisotropic ("stick-like") diffusion. Second, the "density" of tissue components may be confounded by non-diffusion properties such as T2 relaxation. Third, the domain of validity for the estimated parameters to serve as indices of neurite density is incompletely explored. We investigated these challenges by acquiring data with "b-tensor encoding" and multiple echo times in brain regions with low orientation coherence and in white matter lesions. Results showed that microscopic anisotropy from b-tensor data is associated with myelinated axons but not with dendrites. Furthermore, b-tensor data together with data acquired for multiple echo times showed that unbiased density estimates in white matter lesions require data-driven estimates of compartment-specific T2 values. Finally, the "stick" fractions of different biophysical models could generally not serve as neurite density indices across the healthy brain and white matter lesions, where outcomes of comparisons depended on the choice of constraints. In particular, constraining compartment-specific T2 values was ambiguous in the healthy brain and had a large impact on estimated values. In summary, estimating neurite density generally requires accounting for different diffusion and/or T2 properties between axons and dendrites. Constrained "index" parameters could be valid within limited domains that should be delineated by future studies.

Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation

We review, systematize and discuss models of diffusion in neuronal tissue, by putting them into an overarching physical context of coarse-graining over an increasing diffusion length scale. From this perspective, we view research on quantifying brain microstructure as occurring along three major avenues. The first avenue focusses on transient, or time-dependent, effects in diffusion. These effects signify the gradual coarse-graining of tissue structure, which occurs qualitatively differently in different brain tissue compartments. We show that transient effects contain information about the relevant length scales for neuronal tissue, such as the packing correlation length for neuronal fibers, as well as the degree of structural disorder along the neurites. The second avenue corresponds to the long-time limit, when the observed signal can be approximated as a sum of multiple nonexchanging anisotropic Gaussian components. Here, the challenge lies in parameter estimation and in resolving its hidden degeneracies. The third avenue employs multiple diffusion encoding techniques, able to access information not contained in the conventional diffusion propagator. We conclude with our outlook on future directions that could open exciting possibilities for designing quantitative markers of tissue physiology and pathology, based on methods of studying mesoscopic transport in disordered systems.

Imaging brain microstructure with diffusion MRI: practicality and applications

This article gives an overview of microstructure imaging of the brain with diffusion MRI and reviews the state of the art. The microstructure-imaging paradigm aims to estimate and map microscopic properties of tissue using a model that links these properties to the voxel scale MR signal. Imaging techniques of this type are just starting to make the transition from the technical research domain to wide application in biomedical studies. We focus here on the practicalities of both implementing such techniques and using them in applications. Specifically, the article summarizes the relevant aspects of brain microanatomy and the range of diffusion-weighted MR measurements that provide sensitivity to them. It then reviews the evolution of mathematical and computational models that relate the diffusion MR signal to brain tissue microstructure, as well as the expanding areas of application. Next we focus on practicalities of designing a working microstructure imaging technique: model selection, experiment design, parameter estimation, validation, and the pipeline of development of this class of technique. The article concludes with some future perspectives on opportunities in this topic and expectations on how the field will evolve in the short-to-medium term.

Tensor-valued diffusion encoding for diffusional variance decomposition (DIVIDE): Technical feasibility in clinical MRI systems

Microstructure imaging techniques based on tensor-valued diffusion encoding have gained popularity within the MRI research community. Unlike conventional diffusion encoding-applied along a single direction in each shot-tensor-valued encoding employs diffusion encoding along multiple directions within a single preparation of the signal. The benefit is that such encoding may probe tissue features that are not accessible by conventional encoding. For example, diffusional variance decomposition (DIVIDE) takes advantage of tensor-valued encoding to probe microscopic diffusion anisotropy independent of orientation coherence. The drawback is that tensor-valued encoding generally requires gradient waveforms that are more demanding on hardware; it has therefore been used primarily in MRI systems with relatively high performance. The purpose of this work was to explore tensor-valued diffusion encoding on clinical MRI systems with varying performance to test its technical feasibility within the context of DIVIDE. We performed whole-brain imaging with linear and spherical b-tensor encoding at field strengths between 1.5 and 7 T, and at maximal gradient amplitudes between 45 and 80 mT/m. Asymmetric gradient waveforms were optimized numerically to yield b-values up to 2 ms/μm2. Technical feasibility was assessed in terms of the repeatability, SNR, and quality of DIVIDE parameter maps. Variable system performance resulted in echo times between 83 to 115 ms and total acquisition times of 6 to 9 minutes when using 80 signal samples and resolution 2×2×4 mm3. As expected, the repeatability, signal-to-noise ratio and parameter map quality depended on hardware performance. We conclude that tensor-valued encoding is feasible for a wide range of MRI systems-even at 1.5 T with maximal gradient waveform amplitudes of 33 mT/m-and baseline experimental design and quality parameters for all included configurations. This demonstrates that tissue features, beyond those accessible by conventional diffusion encoding, can be explored on a wide range of MRI systems.

Microstructural imaging of the human brain with a 'super-scanner': 10 key advantages of ultra-strong gradients for diffusion MRI

The key component of a microstructural diffusion MRI 'super-scanner' is a dedicated high-strength gradient system that enables stronger diffusion weightings per unit time compared to conventional gradient designs. This can, in turn, drastically shorten the time needed for diffusion encoding, increase the signal-to-noise ratio, and facilitate measurements at shorter diffusion times. This review, written from the perspective of the UK National Facility for In Vivo MR Imaging of Human Tissue Microstructure, an initiative to establish a shared 300 mT/m-gradient facility amongst the microstructural imaging community, describes ten advantages of ultra-strong gradients for microstructural imaging. Specifically, we will discuss how the increase of the accessible measurement space compared to a lower-gradient systems (in terms of Δ, b-value, and TE) can accelerate developments in the areas of 1) axon diameter distribution mapping; 2) microstructural parameter estimation; 3) mapping micro-vs macroscopic anisotropy features with gradient waveforms beyond a single pair of pulsed-gradients; 4) multi-contrast experiments, e.g. diffusion-relaxometry; 5) tractography and high-resolution imaging in vivo and 6) post mortem; 7) diffusion-weighted spectroscopy of metabolites other than water; 8) tumour characterisation; 9) functional diffusion MRI; and 10) quality enhancement of images acquired on lower-gradient systems. We finally discuss practical barriers in the use of ultra-strong gradients, and provide an outlook on the next generation of 'super-scanners'.

Preliminary evaluation of accelerated microscopic diffusional kurtosis imaging (μDKI) in a rodent model of epilepsy

PURPOSE

Our study aimed to develop accelerated microscopic diffusional kurtosis imaging (μDKI) and preliminarily evaluated it in a rodent model of chronic epilepsy.

METHODS

We investigated two μDKI acceleration schemes of reduced sampling density and angular range in a phantom and wild-type rats, and further tested μDKI method in pilocarpine-induced epilepsy rats using a 4.7 Tesla MRI. Single slice average μD_app and μK_app maps were derived, and Nissl staining was obtained.

RESULTS

The kurtosis maps from two accelerated μDKI sampling schemes (sampling density and range) are very similar to that using fully sampled data (SSIM > 0.95). For the epileptic models, μDKI showed noticeably different contrast from those obtained with conventional DKI. Specifically, the average μK_app was significantly less than that of the average of K_app (0.15 ± 0.01 vs. 0.47 ± 0.02) in the ventricle.

CONCLUSIONS

Our study demonstrated the feasibility of accelerated in vivo μDKI. Our work revealed that μDKI provides complementary information to conventional DKI method, suggesting that advanced DKI sequences are promising to elucidate tissue microstructure in neurological diseases.

Influence of the size and curvedness of neural projections on the orientationally averaged diffusion MR signal

Neuronal and glial projections can be envisioned to be tubes of infinitesimal diameter as far as diffusion magnetic resonance (MR) measurements via clinical scanners are concerned. Recent experimental studies indicate that the decay of the orientationally-averaged signal in white-matter may be characterized by the power-law, Ē(q) ∝ q^-1, where q is the wavenumber determined by the parameters of the pulsed field gradient measurements. One particular study by McKinnon et al. [1] reports a distinctively faster decay in gray-matter. Here, we assess the role of the size and curvature of the neurites and glial arborizations in these experimental findings. To this end, we studied the signal decay for diffusion along general curves at all three temporal regimes of the traditional pulsed field gradient measurements. We show that for curvy projections, employment of longer pulse durations leads to a disappearance of the q^-1 decay, while such decay is robust when narrow gradient pulses are used. Thus, in clinical acquisitions, the lack of such a decay for a fibrous specimen can be seen as indicative of fibers that are curved. We note that the above discussion is valid for an intermediate range of q-values as the true asymptotic behavior of the signal decay is Ē(q) ∝ q^-4 for narrow pulses (through Debye-Porod law) or steeper for longer pulses. This study is expected to provide insights for interpreting the diffusion-weighted images of the central nervous system and aid in the design of acquisition strategies.

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.

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.

Internal gradient distributions: A susceptibility-derived tensor delivering morphologies by magnetic resonance

Nuclear magnetic resonance is a powerful tool for probing the structures of chemical and biological systems. Combined with field gradients it leads to NMR imaging (MRI), a widespread tool in non-invasive examinations. Sensitivity usually limits MRI’s spatial resolution to tens of micrometers, but other sources of information like those delivered by constrained diffusion processes, enable one extract morphological information down to micron and sub-micron scales. We report here on a new method that also exploits diffusion – isotropic or anisotropic– to sense morphological parameters in the nm-mm range, based on distributions of susceptibility-induced magnetic field gradients. A theoretical framework is developed to define this source of information, leading to the proposition of internal gradient-distribution tensors. Gradient-based spin-echo sequences are designed to measure these new observables. These methods can be used to map orientations even when dealing with unconstrained diffusion, as is here demonstrated with studies of structured systems, including tissues.

Studying microstructure and microstructural changes in plant tissues by advanced diffusion magnetic resonance imaging techniques

As sessile organisms, plants must respond to the environment by adjusting their growth and development. Most of the plant body is formed post-embryonically by continuous activity of apical and lateral meristems. The development of lateral adventitious roots is a complex process, and therefore the development of methods that can visualize, non-invasively, the plant microstructure and organ initiation that occur during growth and development is of paramount importance. In this study, relaxation-based and advanced diffusion magnetic resonance imaging (MRI) methods including diffusion tensor (DTI), q-space diffusion imaging (QSI), and double-pulsed-field-gradient (d-PFG) MRI, at 14.1 T, were used to characterize the hypocotyl microstructure and the microstructural changes that occurred during the development of lateral adventitious roots in tomato. Better contrast was observed in relaxation-based MRI using higher in-plane resolution but this also resulted in a significant reduction in the signal-to-noise ratio of the T2-weighted MR images. Diffusion MRI revealed that water diffusion is highly anisotropic in the vascular cylinder. QSI and d-PGSE MRI showed that in the vascular cylinder some of the cells have sizes in the range of 6-10 μm. The MR images captured cell reorganization during adventitious root formation in the periphery of the primary vascular bundles, adjacent to the xylem pole that broke through the cortex and epidermis layers. This study demonstrates that MRI and diffusion MRI methods allow the non-invasive study of microstructural features of plants, and enable microstructural changes associated with adventitious root formation to be followed.

Anisotropic phantom to calibrate high-q diffusion MRI methods

A silicon oil-filled glass capillary array is proposed as an anisotropic diffusion MRI phantom. Together with a computational/theoretical pipeline these provide a gold standard for calibrating and validating high-q diffusion MRI experiments. The phantom was used to test high angular resolution diffusion imaging (HARDI) and double pulsed-field gradient (d-PFG) MRI acquisition schemes. MRI-based predictions of microcapillary diameter using both acquisition schemes were compared with results from optical microscopy. This phantom design can be used for quality control and quality assurance purposes and for testing and validating proposed microstructure imaging experiments and the processing pipelines used to analyze them.

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.

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.

Two-Dimensional Correlation of Isotropic and Directional Diffusion Using NMR

Diffusion nuclear magnetic resonance (NMR) is a powerful technique for studying porous media, but yields ambiguous results when the sample comprises multiple regions with different pore sizes, shapes, and orientations. Inspired by solid-state NMR techniques for correlating isotropic and anisotropic chemical shifts, we propose a diffusion NMR method to resolve said ambiguity. Numerical data inversion relies on sparse representation of the data in a basis of radial and axial diffusivities. Experiments are performed on a composite sample with a cell suspension and a liquid crystal.

Eddy current compensated double diffusion encoded (DDE) MRI

PURPOSE

Eddy currents might lead to image distortions in diffusion-weighted echo planar imaging. A method is proposed to reduce their effects on double diffusion encoding (DDE) MRI experiments and the thereby derived microscopic fractional anisotropy (μFA).

METHODS

The twice-refocused spin echo scheme was adapted for DDE measurements. To assess the effect of individual diffusion encodings on the image distortions, measurements of a grid of plastic rods in water were performed. The effect of eddy current compensation on μFA measurements was evaluated in the brains of six healthy volunteers.

RESULTS

The use of an eddy current compensation reduced the signal variation. As expected, the distortions caused by the second encoding were larger than those of the first encoding, entailing a stronger need to compensate for them. For an optimal result, however, both encodings had to be compensated. The artifact reduction strongly improved the measurement of the μFA in ventricles and gray matter by reducing the overestimation. An effect of the compensation on absolute μFA values in white matter was not observed.

CONCLUSION

It is advisable to compensate both encodings in DDE measurements for eddy currents. Magn Reson Med 77:328-335, 2017. © 2015 Wiley Periodicals, Inc.

Detection of acute nervous system injury with advanced diffusion-weighted MRI: a simulation and sensitivity analysis

Diffusion-weighted imaging (DWI) is a powerful tool to investigate the microscopic structure of the central nervous system (CNS). Diffusion tensor imaging (DTI), a common model of the DWI signal, has a demonstrated sensitivity to detect microscopic changes as a result of injury or disease. However, DTI and other similar models have inherent limitations that reduce their specificity for certain pathological features, particularly in tissues with complex fiber arrangements. Methods such as double pulsed field gradient (dPFG) and q-vector magic angle spinning (qMAS) have been proposed to specifically probe the underlying microscopic anisotropy without interference from the macroscopic tissue organization. This is particularly important for the study of acute injury, where abrupt changes in the microscopic morphology of axons and dendrites manifest as focal enlargements known as beading. The purpose of this work was to assess the relative sensitivity of DWI measures to beading in the context of macroscopic fiber organization and edema. Computational simulations of DWI experiments in normal and beaded axons demonstrated that, although DWI models can be highly specific for the simulated pathologies of beading and volume fraction changes in coherent fiber pathways, their sensitivity to a single idealized pathology is considerably reduced in crossing and dispersed fibers. However, dPFG and qMAS have a high sensitivity for beading, even in complex fiber tracts. Moreover, in tissues with coherent arrangements, such as the spinal cord or nerve fibers in which tract orientation is known a priori, a specific dPFG sequence variant decreases the effects of edema and improves specificity for beading. Collectively, the simulation results demonstrate that advanced DWI methods, particularly those which sample diffusion along multiple directions within a single acquisition, have improved sensitivity to acute axonal injury over conventional DTI metrics and hold promise for more informative clinical diagnostic use in CNS injury evaluation.

Microscopic diffusion anisotropy in the human brain: reproducibility, normal values, and comparison with the fractional anisotropy

Human neuroimaging of tissue microstructure, such as axonal density and integrity, is key in clinical and neuroscience research. Most studies rely on diffusion tensor imaging (DTI) and the measures derived from it, most prominently fractional anisotropy (FA). However, FA also depends on fiber orientation distribution, a more macroscopic tissue property. Recently introduced measures of so-called microscopic diffusion anisotropy, diffusion anisotropy on a cellular or microscopic level, overcome this limitation because they are independent of the orientation distributions of axons and fibers. In this study, we evaluate the feasibility of two measures of microscopic diffusion anisotropy I(MA) and MA indices, for human neuroscience and clinical research. Both indices reflect the eccentricity of the cells but while I(MA) also depends on the cell size, MA is independent of the cell size and, like FA, scaled between 0 and 1. In whole-brain measurements of a group of 19 healthy volunteers, we measured average values and variability, evaluated their reproducibility, both within and between sessions, and compared MA to FA values in selected regions-of-interest (ROIs). The within- and between-session comparison did not show substantial differences but the reproducibility was much better for the MA than I(MA) (coefficient of variation between sessions 10.5% vs. 28.9%). The reproducibility was less for MA than FA overall, but comparable in the defined ROIs and the average group sizes required for between-group comparisons was similar (about 60 participants for a relative difference of 5%). Group-averaged values of MA index were generally larger and showed less variation across white-matter brain ROIs than FA (mean ± standard deviation of seven ROIs 0.83 ± 0.10 vs. 0.58 ± 0.13). Even in some gray-matter ROIs, MA values comparable to those of white matter ROIs were observed. Furthermore, the within-group variation of the values in white matter ROIs was lower for the MA compared to the FA (mean standard deviation over volunteers 0.038 vs. 0.049) which could be due to significant variability in the distribution of fiber orientation contributing to FA. These results indicate that MA (i) should be preferred to I(MA), (ii) has a reproducibility and group-size requirements comparable to those of FA; (iii) is less sensitive to the fiber orientation distribution than FA; and (iv) could be more sensitive to differences or changes of the tissue microstructure than FA. R1.1.

High-resolution diffusion kurtosis imaging at 3T enabled by advanced post-processing

Diffusion Kurtosis Imaging (DKI) is more sensitive to microstructural differences and can be related to more specific micro-scale metrics (e.g., intra-axonal volume fraction) than diffusion tensor imaging (DTI), offering exceptional potential for clinical diagnosis and research into the white and gray matter. Currently DKI is acquired only at low spatial resolution (2–3 mm isotropic), because of the lower signal-to-noise ratio (SNR) and higher artifact level associated with the technically more demanding DKI. Higher spatial resolution of about 1 mm is required for the characterization of fine white matter pathways or cortical microstructure. We used restricted-field-of-view (rFoV) imaging in combination with advanced post-processing methods to enable unprecedented high-quality, high-resolution DKI (1.2 mm isotropic) on a clinical 3T scanner. Post-processing was advanced by developing a novel method for Retrospective Eddy current and Motion ArtifacT Correction in High-resolution, multi-shell diffusion data (REMATCH). Furthermore, we applied a powerful edge preserving denoising method, denoted as multi-shell orientation-position-adaptive smoothing (msPOAS). We demonstrated the feasibility of high-quality, high-resolution DKI and its potential for delineating highly myelinated fiber pathways in the motor cortex. REMATCH performs robustly even at the low SNR level of high-resolution DKI, where standard EC and motion correction failed (i.e., produced incorrectly aligned images) and thus biased the diffusion model fit. We showed that the combination of REMATCH and msPOAS increased the contrast between gray and white matter in mean kurtosis (MK) maps by about 35% and at the same time preserves the original distribution of MK values, whereas standard Gaussian smoothing strongly biases the distribution.