Numerical simulations of short-mixing-time double-wave-vector diffusion-weighting experiments with multiple concatenations on whole-body MR systems

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.

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.

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.

NMR-based diffusion pore imaging

Nuclear magnetic resonance (NMR) diffusion experiments offer a unique opportunity to study boundaries restricting the diffusion process. In a recent Letter [Phys. Rev. Lett. 107, 048102 (2011)], we introduced the idea and concept that such diffusion experiments can be interpreted as NMR imaging experiments. Consequently, images of closed pores, in which the spins diffuse, can be acquired. In the work presented here, an in-depth description of the diffusion pore imaging technique is provided. Image artifacts due to gradient profiles of finite duration, field inhomogeneities, and surface relaxation are considered. Gradients of finite duration lead to image blurring and edge enhancement artifacts. Field inhomogeneities have benign effects on diffusion pore images, and surface relaxation can lead to a shrinkage and shift of the pore image. The relation between boundary structure and the imaginary part of the diffusion weighted signal is analyzed, and it is shown that information on pore coherence can be obtained without the need to measure the phase of the diffusion weighted signal. Moreover, it is shown that quite arbitrary gradient profiles can be used for diffusion pore imaging. The matrices required for numerical calculations are stated and provided as supplemental material.

Towards compartment size estimation in vivo based on double wave vector diffusion weighting

It has been shown that double wave vector diffusion weighting, which employs two gradient pulse pairs of independent directions, can provide information about tissue structure that is not easily available otherwise, such as cell size or shape in a tissue sample. One approach to measure cell size is based on the signal difference between parallel and antiparallel gradient orientations at small mixing times between the two diffusion weightings. A major difficulty for in vivo application is the small size of the signal difference if clinical MR systems with limited gradient hardware are employed. In this study, the method is applied to human brain tissue in vivo, using whole-body gradients. Data are reported for the corticospinal tracts. The characteristics of the observed signal difference between parallel and antiparallel gradient orientations are consistent with both analytical and numerical predictions. As an estimate of pore size, the resulting mean squared radius of gyration of the pores amounts to approximately 4 µm(2) . An analysis that accounts for the finite values of gradient pulse duration and diffusion time yields a volume contribution-weighted mean pore diameter of 13 μm if a cylindrical pore shape is assumed. The results demonstrate that the technique can be applied in vivo.

Overcoming apparent susceptibility-induced anisotropy (aSIA) by bipolar double-pulsed-field-gradient NMR

Double-Pulsed-Field-Gradient (d-PFG) MR is emerging as a powerful new means for obtaining unique microstructural information in opaque porous systems that cannot be obtained by conventional single-PFG (s-PFG) methods. The angular d-PFG MR methodology is particularly important since it can utilize the effects of microscopic anisotropy (μA) and compartment shape anisotropy (csA) in the E(ψ) profile at the different t(m) regimes to provide detailed information on compartment size and eccentricity. An underlying assumption is that the PFGs that are imparted to weigh diffusion are the only gradients present; however, in realistic systems and especially where there are randomly oriented anisotropic pores, susceptibility effects may induce strong internal gradients. In this study, the effects of such internal gradients on E(ψ) plots obtained from angular d-PFG MR and on microstructural information that can be obtained from s-PFG and d-PFG MR were investigated. First, it was found that internal gradients induce a bias in the s-PFG MR results, thus creating an anisotropy that is not related to microstructure, termed apparent-Susceptibility-Induced-Anisotropy (aSIA). We then show that aSIA effects are also manifest in different ways in the angular d-PFG MR experiment in controlled phantoms and in realistic systems such as quartz sand, emulsions, and biological systems. The effects of aSIA in some cases completely masked the effects of μA and csA; however, we subsequently show that by introducing bipolar gradients to the d-PFG MR (bp-d-PFG), the effects of aSIA can be largely suppressed, restoring the E(ψ) plots that are expected from the theory along with the microstructural information that it conveys. We conclude that when specimens are characterized by strong internal gradients, the novel information on μA and csA that is manifest in the E(ψ) plots can indeed be inferred when bp-d-PFG MR is used, i.e. when bipolar gradients are applied.