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 B0 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 B0 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.

Diffusion time dependence, power-law scaling, and exchange in gray matter

Characterizing neural tissue microstructure is a critical goal for future neuroimaging. Diffusion MRI (dMRI) provides contrasts that reflect diffusing spins' interactions with myriad microstructural features of biological systems. However, the specificity of dMRI remains limited due to the ambiguity of its signals vis-à-vis the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) taken spherical compartments into account (the SANDI model) in attempts to represent cell soma. The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. To validate these assumptions experimentally, observing their unique, functional properties, such as the b-1/2 power-law associated with one-dimensional diffusion, has emerged as a fruitful strategy. The absence of this signature in GM, in turn, has been explained by neurite water exchange, non-linear morphology, and/or by obscuring soma signal contributions. Here, we present diffusion simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law behavior in impermeable neurites, but also that their signal is drowned by the soma signal under typical experimental conditions. Nevertheless, by studying the GM dMRI signal's behavior as a function of diffusion weighting as well as time, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a signal behavior opposite to that which would be expected from restricted diffusion, thereby providing a functional signature that disambiguates the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field of 16.4T and observe a time dependence that is consistent with substantial exchange but also with a GM stick power-law. The first finding suggests significant water exchange between neurites and the extracellular space while the second suggests a small sub-population of impermeable neurites. To quantify these observations, we harness the Kärger exchange model and incorporate the corresponding signal time dependence in the SM and SANDI models.

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.

Test-retest reproducibility of in vivo oscillating gradient and microscopic anisotropy diffusion MRI in mice at 9.4 Tesla

BACKGROUND AND PURPOSE

Microstructure imaging with advanced diffusion MRI (dMRI) techniques have shown increased sensitivity and specificity to microstructural changes in various disease and injury models. Oscillating gradient spin echo (OGSE) dMRI, implemented by varying the oscillating gradient frequency, and microscopic anisotropy (μA) dMRI, implemented via tensor valued diffusion encoding, may provide additional insight by increasing sensitivity to smaller spatial scales and disentangling fiber orientation dispersion from true microstructural changes, respectively. The aims of this study were to characterize the test-retest reproducibility of in vivo OGSE and μA dMRI metrics in the mouse brain at 9.4 Tesla and provide estimates of required sample sizes for future investigations.

METHODS

Twelve adult C57Bl/6 mice were scanned twice (5 days apart). Each imaging session consisted of multifrequency OGSE and μA dMRI protocols. Metrics investigated included μA, linear diffusion kurtosis, isotropic diffusion kurtosis, and the diffusion dispersion rate (Λ), which explores the power-law frequency dependence of mean diffusivity. The dMRI metric maps were analyzed with mean region-of-interest (ROI) and whole brain voxel-wise analysis. Bland-Altman plots and coefficients of variation (CV) were used to assess the reproducibility of OGSE and μA metrics. Furthermore, we estimated sample sizes required to detect a variety of effect sizes.

RESULTS

Bland-Altman plots showed negligible biases between test and retest sessions. ROI-based CVs revealed high reproducibility for most metrics (CVs < 15%). Voxel-wise CV maps revealed high reproducibility for μA (CVs ~ 10%), but low reproducibility for OGSE metrics (CVs ~ 50%).

CONCLUSION

Most of the μA dMRI metrics are reproducible in both ROI-based and voxel-wise analysis, while the OGSE dMRI metrics are only reproducible in ROI-based analysis. Given feasible sample sizes (10-15), μA metrics and OGSE metrics may provide sensitivity to subtle microstructural changes (4-8%) and moderate changes (> 6%), respectively.

Double diffusion encoding and applications for biomedical imaging

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

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

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

Challenges for biophysical modeling of microstructure

The biophysical modeling efforts in diffusion MRI have grown considerably over the past 25 years. In this review, we dwell on the various challenges along the journey of bringing a biophysical model from initial design to clinical implementation, identifying both hurdles that have been already overcome and outstanding issues. First, we describe the critical initial task of selecting which features of tissue microstructure can be estimated using a model and which acquisition protocol needs to be implemented to make the estimation possible. The model performance should necessarily be tested in realistic numerical simulations and in experimental data - adapting the fitting strategy accordingly, and parameter estimates should be validated against complementary techniques, when/if available. Secondly, the model performance and validity should be explored in pathological conditions, and, if appropriate, dedicated models for pathology should be developed. We build on examples from tumors, ischemia and demyelinating diseases. We then discuss the challenges associated with clinical translation and added value. Finally, we single out four major unresolved challenges that are related to: the availability of a microstructural ground truth, the validation of model parameters which cannot be accessed with complementary techniques, the development of a generalized standard model for any brain region and pathology, and the seamless communication between different parties involved in the development and application of biophysical models of diffusion.

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.

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.

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.

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.

Microscopic anisotropy misestimation in spherical-mean single diffusion encoding MRI

Purpose

Microscopic fractional anisotropy (µFA) can disentangle microstructural information from orientation dispersion. While double diffusion encoding (DDE) MRI methods are widely used to extract accurate µFA, it has only recently been proposed that powder‐averaged single diffusion encoding (SDE) signals, when coupled with the diffusion standard model (SM) and a set of constraints, could be used for µFA estimation. This study aims to evaluate µFA as derived from the spherical mean technique (SMT) set of constraints, as well as more generally for powder‐averaged SM signals.

Methods

SDE experiments were performed at 16.4 T on an ex vivo mouse brain (Δ/δ = 12/1.5 ms). The µFA maps obtained from powder‐averaged SDE signals were then compared to maps obtained from DDE‐MRI experiments (Δ/τ/δ = 12/12/1.5 ms), which allow a model‐free estimation of µFA. Theory and simulations that consider different types of heterogeneity are presented for corroborating the experimental findings.

Results

µFA, as well as other estimates derived from powder‐averaged SDE signals produced large deviations from the ground truth in both gray and white matter. Simulations revealed that these misestimations are likely a consequence of factors not considered by the underlying microstructural models (such as intercomponent and intracompartmental kurtosis).

Conclusion

Powder‐averaged SMT and (2‐component) SM are unable to accurately report µFA and other microstructural parameters in ex vivo tissues. Improper model assumptions and constraints can significantly compromise parameter specificity. Further developments and validations are required prior to implementation of these models in clinical or preclinical research.

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.

Microstructural imaging in the spinal cord and validation strategies

In vivo histology using magnetic resonance imaging (MRI) is a newly emerging research field that aims to non-invasively characterize tissue microstructure. The implications of in vivo histology are many, from discovering novel biomarkers to studying human development, to providing tools for disease diagnosis and monitoring the effects of novel treatments on tissue. This review focuses on quantitative MRI (qMRI) techniques that are used to map spinal cord microstructure. Opening with a rationale for non-invasive imaging of the spinal cord, this article continues with a brief overview of the existing MRI techniques for axon and myelin imaging, followed by the specific challenges and potential solutions for acquiring and processing such data. The final part of this review focuses on histological validation, with suggested tissue preparation, acquisition and processing protocols for large-scale microscopy.

Distinguishing neuronal from astrocytic subcellular microstructures using in vivo Double Diffusion Encoded 1H MRS at 21.1 T

Measuring cellular microstructures non-invasively and achieving specificity towards a cell-type population within an interrogated in vivo tissue, remains an outstanding challenge in brain research. Magnetic Resonance Spectroscopy (MRS) provides an opportunity to achieve cellular specificity via the spectral resolution of metabolites such as N-Acetylaspartate (NAA) and myo-Inositol (mI), which are considered neuronal and astrocytic markers, respectively. Yet the information typically obtained with MRS describes metabolic concentrations, diffusion coefficients or relaxation rates rather than microstructures. Understanding how these metabolites are compartmentalized is a challenging but important goal, which so far has been mainly addressed using diffusion models. Here, we present direct in vivo evidence for the confinement of NAA and mI within sub-cellular components, namely, the randomly oriented process of neurons and astrocytes, respectively. Our approach applied Relaxation Enhanced MRS at ultrahigh (21.1 T) field, and used its high ¹H sensitivity to measure restricted diffusion correlations for NAA and mI using a Double Diffusion Encoding (DDE) filter. While very low macroscopic anisotropy was revealed by spatially localized Diffusion Tensor Spectroscopy, DDE displayed characteristic amplitude modulations reporting on confinements in otherwise randomly oriented anisotropic microstructures for both metabolites. This implies that for the chosen set of parameters, the DDE measurements had a biased sensitivity towards NAA and mI sited in the more confined environments of neurites and astrocytic branches, than in the cell somata. These measurements thus provide intrinsic diffusivities and compartment diameters, and revealed subcellular neuronal and astrocytic morphologies in normal in vivo rat brains. The relevance of these measurements towards human applications—which could in turn help understand CNS plasticity as well as diagnose brain diseases—is discussed.

Promise and pitfalls of g-ratio estimation with MRI

The fiber g-ratio is the ratio of the inner to the outer diameter of the myelin sheath of a myelinated axon. It has a limited dynamic range in healthy white matter, as it is optimized for speed of signal conduction, cellular energetics, and spatial constraints. In vivo imaging of the g-ratio in health and disease would greatly increase our knowledge of the nervous system and our ability to diagnose, monitor, and treat disease. MRI based g-ratio imaging was first conceived in 2011, and expanded to be feasible in full brain white matter with preliminary results in 2013. This manuscript reviews the growing g-ratio imaging literature and speculates on future applications. It details the methodology for imaging the g-ratio with MRI, and describes the known pitfalls and challenges in doing so.

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.

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.

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.

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.

Size Distribution Imaging by Non-Uniform Oscillating-Gradient Spin Echo (NOGSE) MRI

Objects making up complex porous systems in Nature usually span a range of sizes. These size distributions play fundamental roles in defining the physicochemical, biophysical and physiological properties of a wide variety of systems - ranging from advanced catalytic materials to Central Nervous System diseases. Accurate and noninvasive measurements of size distributions in opaque, three-dimensional objects, have thus remained long-standing and important challenges. Herein we describe how a recently introduced diffusion-based magnetic resonance methodology, Non-Uniform-Oscillating-Gradient-Spin-Echo (NOGSE), can determine such distributions noninvasively. The method relies on its ability to probe confining lengths with a (length)6 parametric sensitivity, in a constant-time, constant-number-of-gradients fashion; combined, these attributes provide sufficient sensitivity for characterizing the underlying distributions in μm-scaled cellular systems. Theoretical derivations and simulations are presented to verify NOGSE's ability to faithfully reconstruct size distributions through suitable modeling of their distribution parameters. Experiments in yeast cell suspensions - where the ground truth can be determined from ancillary microscopy - corroborate these trends experimentally. Finally, by appending to the NOGSE protocol an imaging acquisition, novel MRI maps of cellular size distributions were collected from a mouse brain. The ensuing micro-architectural contrasts successfully delineated distinctive hallmark anatomical sub-structures, in both white matter and gray matter tissues, in a non-invasive manner. Such findings highlight NOGSE's potential for characterizing aberrations in cellular size distributions upon disease, or during normal processes such as development.

Diffusion-MRI in neurodegenerative disorders

The ability to image the whole brain through ever more subtle and specific methods/contrasts has come to play a key role in understanding the basis of brain abnormalities in several diseases. In magnetic resonance imaging (MRI), "diffusion" (i.e. the random, thermally-induced displacements of water molecules over time) represents an extraordinarily sensitive contrast mechanism, and the exquisite structural detail it affords has proven useful in a vast number of clinical as well as research applications. Since diffusion-MRI is a truly quantitative imaging technique, the indices it provides can serve as potential imaging biomarkers which could allow early detection of pathological alterations as well as tracking and possibly predicting subtle changes in follow-up examinations and clinical trials. Accordingly, diffusion-MRI has proven useful in obtaining information to better understand the microstructural changes and neurophysiological mechanisms underlying various neurodegenerative disorders. In this review article, we summarize and explore the main applications, findings, perspectives as well as challenges and future research of diffusion-MRI in various neurodegenerative disorders including Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, Huntington's disease and degenerative ataxias.

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.

Measuring small compartment dimensions by probing diffusion dynamics via Non-uniform Oscillating-Gradient Spin-Echo (NOGSE) NMR

Noninvasive measurements of microstructure in materials, cells, and in biological tissues, constitute a unique capability of gradient-assisted NMR. Diffusion-diffraction MR approaches pioneered by Callaghan demonstrated this ability; Oscillating-Gradient Spin-Echo (OGSE) methodologies tackle the demanding gradient amplitudes required for observing diffraction patterns by utilizing constant-frequency oscillating gradient pairs that probe the diffusion spectrum, D(ω). Here we present a new class of diffusion MR experiments, termed Non-uniform Oscillating-Gradient Spin-Echo (NOGSE), which dynamically probe multiple frequencies of the diffusion spectral density at once, thus affording direct microstructural information on the compartment's dimension. The NOGSE methodology applies N constant-amplitude gradient oscillations; N-1 of these oscillations are spaced by a characteristic time x, followed by a single gradient oscillation characterized by a time y, such that the diffusion dynamics is probed while keeping (N-1)x+y≡TNOGSE constant. These constant-time, fixed-gradient-amplitude, multi-frequency attributes render NOGSE particularly useful for probing small compartment dimensions with relatively weak gradients - alleviating difficulties associated with probing D(ω) frequency-by-frequency or with varying relaxation weightings, as in other diffusion-monitoring experiments. Analytical descriptions of the NOGSE signal are given, and the sequence's ability to extract small compartment sizes with a sensitivity towards length to the sixth power, is demonstrated using a microstructural phantom. Excellent agreement between theory and experiments was evidenced even upon applying weak gradient amplitudes. An MR imaging version of NOGSE was also implemented in ex vivo pig spinal cords and mouse brains, affording maps based on compartment sizes. The effects of size distributions on NOGSE are also briefly analyzed.

Diffusion weighted MRI by spatiotemporal encoding: analytical description and in vivo validations

Diffusion-weighted (DW) MRI is a powerful modality for studying microstructure in normal and pathological tissues. The accuracy derived from DW MRI depends on the acquisition of quality images, and on a precise assessment of the b-values involved. Conventional DW MRI tends to be of limited use in regions suffering from large magnetic field or chemical shift heterogeneities, which severely distort the MR images. In this study we propose novel sequences based on SPatio-temporal ENcoding (SPEN), which overcome such shortcomings owing to SPEN's inherent robustness to offsets. SPEN, however, relies on the simultaneous application of gradients and radiofrequency-swept pulses, which may impart different diffusion weightings along the spatial axes. These will be further complicated in DW measurements by the diffusion-sensitizing gradients, and will in general lead to complex, spatially-dependent b-values. This study presents a formalism for analyzing these diffusion-weighted SPEN (dSPEN) data, which takes into account the concomitant effects of adiabatic pulses, of the imaging as well as diffusion gradients, and of the cross-terms between them. These analytical b-values derivations are subject to experimental validations in phantom systems and ex vivo spinal cords. Excellent agreement is found between the theoretical predictions and these dSPEN experiments. The ensuing methodology is then demonstrated by in vivo mapping of diffusion in human breast - organs where conventional k-space DW acquisition methods are challenged by both field and chemical shift heterogeneities. These studies demonstrate the increased robustness of dSPEN vis-à-vis comparable DW echo planar imaging, and demonstrate the value of this new methodology for medium- or high-field diffusion measurements in heterogeneous systems.

The CONNECT project: Combining macro- and micro-structure

In recent years, diffusion MRI has become an extremely important tool for studying the morphology of living brain tissue, as it provides unique insights into both its macrostructure and microstructure. Recent applications of diffusion MRI aimed to characterize the structural connectome using tractography to infer connectivity between brain regions. In parallel to the development of tractography, additional diffusion MRI based frameworks (CHARMED, AxCaliber, ActiveAx) were developed enabling the extraction of a multitude of micro-structural parameters (axon diameter distribution, mean axonal diameter and axonal density). This unique insight into both tissue microstructure and connectivity has enormous potential value in understanding the structure and organization of the brain as well as providing unique insights to abnormalities that underpin disease states. The CONNECT (Consortium Of Neuroimagers for the Non-invasive Exploration of brain Connectivity and Tracts) project aimed to combine tractography and micro-structural measures of the living human brain in order to obtain a better estimate of the connectome, while also striving to extend validation of these measurements. This paper summarizes the project and describes the perspective of using micro-structural measures to study the connectome.

In vivo detection of microscopic anisotropy using quadruple pulsed-field gradient (qPFG) diffusion MRI on a clinical scanner

We report our design and implementation of a quadruple pulsed-field gradient (qPFG) diffusion MRI pulse sequence on a whole-body clinical scanner and demonstrate its ability to non-invasively detect restriction-induced microscopic anisotropy in human brain tissue. The microstructural information measured using qPFG diffusion MRI in white matter complements that provided by diffusion tensor imaging (DTI) and exclusively characterizes diffusion of water trapped in microscopic compartments with unique measures of average cell geometry. We describe the effect of white matter fiber orientation on the expected MR signal and highlight the importance of incorporating such information in the axon diameter measurement using a suitable mathematical framework. Integration of qPFG diffusion-weighted images (DWI) with fiber orientations measured using high-resolution DTI allows the estimation of average axon diameters in the corpus callosum of healthy human volunteers. Maps of inter-hemispheric average axon diameters reveal an anterior-posterior variation in good topographical agreement with anatomical measurements reported in previous post-mortem studies. With further technical refinements and additional clinical validation, qPFG diffusion MRI could provide a quantitative whole-brain histological assessment of white and gray matter, enabling a wide range of neuroimaging applications for improved diagnosis of neurodegenerative pathologies, monitoring neurodevelopmental processes, and mapping brain connectivity.