Rapid microscopic fractional anisotropy imaging via an optimized linear regression formulation

Investigating microstructural changes between in vivo and perfused ex vivo marmoset brains using oscillating gradient and b-tensor encoded diffusion MRI at 9.4 T

PURPOSE

To investigate microstructural alterations induced by perfusion fixation in brain tissues using advanced diffusion MRI techniques and estimate their potential impact on the application of ex vivo models to in vivo microstructure.

METHODS

We used oscillating gradient spin echo (OGSE) and b-tensor encoding diffusion MRI to examine in vivo and ex vivo microstructural differences in the marmoset brain. OGSE was used to shorten effective diffusion times, whereas b-tensor encoding allowed for the differentiation of isotropic and anisotropic kurtosis. Additionally, we performed Monte Carlo simulations to estimate the potential microstructural changes in the tissues.

RESULTS

We report large changes (˜50%-60%) in kurtosis frequency dispersion (OGSE) and in both anisotropic and isotropic kurtosis (b-tensor encoding) after perfusion fixation. Structural MRI showed an average volume reduction of about 10%. Monte Carlo simulations indicated that these alterations could likely be attributed to extracellular fluid loss possibly combined with axon beading and increased dot compartment signal fraction. Little evidence was observed for reductions in axonal caliber.

CONCLUSION

Our findings shed light on advanced MRI parameter changes that are induced by perfusion fixation and potential microstructural sources for these changes. This work also suggests that caution should be exercised when applying ex vivo models to infer in vivo tissue microstructure, as significant differences may arise.

Developmental Mismatch Across Brain Modalities in Young Children

Background: Brain development during the preschool period is complex and extensive and underlies ongoing behavioral and cognitive maturation. Increasing understanding of typical brain maturation during this time is critical to early identification of atypical development and could inform treatments and interventions. Previous studies have suggested mismatches between brain structural and functional development in later childhood and adolescence. The current study aimed to delineate the developmental matches and mismatches between brain measures from multiple magnetic resonance imaging modalities in young children. Methods: Brain volume, cortical thickness, fractional anisotropy, cerebral blood flow (CBF), amplitude of low-frequency fluctuations (ALFF), regional homogeneity (ReHo), and eigenvector centrality mapping (ECM) were included. Multi-modal neuroimages for 159 datasets from 67 typically developing preschoolers (2.0-7.6 years old) were collected and analyzed. Results: Functional measures (CBF, ECM, ReHo, ALFF) had similar developmental trajectories across regions, whereas development trajectories for brain volumes and cortical thickness were more heterogeneous. Furthermore, within individuals, brain volumes and cortical thickness were very good at predicting individual scans from prior longitudinal scans. Conclusions: These findings provide a more detailed characterization of the complex interplay of different types of brain development in the early years, laying the foundation for future studies on the impact of environmental factors and neurodevelopmental disorders on the development matches/mismatches patterns between brain areas and modalities.

Nonparametric distributions of tensor-valued Lorentzian diffusion spectra for model-free data inversion in multidimensional diffusion MRI

Magnetic resonance imaging (MRI) is the method of choice for noninvasive studies of micrometer-scale structures in biological tissues via their effects on the time- and frequency-dependent (restricted) and anisotropic self-diffusion of water. While new designs of time-dependent magnetic field gradient waveforms have enabled disambiguation between different aspects of translational motion that are convolved in traditional MRI methods relying on single pairs of field gradient pulses, data analysis for complex heterogeneous materials remains a challenge. Here, we propose and demonstrate nonparametric distributions of tensor-valued Lorentzian diffusion spectra, or "D(ω) distributions," as a general representation with sufficient flexibility to describe the MRI signal response from a wide range of model systems and biological tissues investigated with modulated gradient waveforms separating and correlating the effects of restricted and anisotropic diffusion.

Optimal experimental design and estimation for q-space trajectory imaging

Tensor-valued diffusion encoding facilitates data analysis by q-space trajectory imaging. By modeling the diffusion signal of heterogeneous tissues with a diffusion tensor distribution (DTD) and modulating the encoding tensor shape, this novel approach allows disentangling variations in diffusivity from microscopic anisotropy, orientation dispersion, and mixtures of multiple isotropic diffusivities. To facilitate the estimation of the DTD parameters, a parsimonious acquisition scheme coupled with an accurate and precise estimation of the DTD is needed. In this work, we create two precision-optimized acquisition schemes: one that maximizes the precision of the raw DTD parameters, and another that maximizes the precision of the scalar measures derived from the DTD. The improved precision of these schemes compared to a naïve sampling scheme is demonstrated in both simulations and real data. Furthermore, we show that the weighted linear least squares (WLLS) estimator that uses the squared reciprocal of the noisy signal as weights can be biased, whereas the iteratively WLLS estimator with the squared reciprocal of the predicted signal as weights outperforms the conventional unweighted linear LS and nonlinear LS estimators in terms of accuracy and precision. Finally, we show that the use of appropriate constraints can considerably increase the precision of the estimator with only a limited decrease in accuracy.

A longitudinal microstructural MRI dataset in healthy C57Bl/6 mice at 9.4 Tesla

Multimodal microstructural MRI has shown increased sensitivity and specificity to changes in various brain disease and injury models in the preclinical setting. Here, we present an in vivo longitudinal dataset, including a subset of ex vivo data, acquired as control data and to investigate microstructural changes in the healthy mouse brain. The dataset consists of structural T2-weighted imaging, magnetization transfer ratio and saturation imaging, and advanced quantitative diffusion MRI (dMRI) methods. The dMRI methods include oscillating gradient spin echo (OGSE) dMRI and microscopic anisotropy (μA) dMRI, which provide additional insight by increasing sensitivity to smaller spatial scales and disentangling fiber orientation dispersion from true microstructural changes, respectively. The technical skills required to analyze microstructural MRI data are complex and include MRI sequence development, acquisition, and computational neuroimaging expertise. Here, we share unprocessed and preprocessed data, and scalar maps of quantitative MRI metrics. We envision utility of this dataset in the microstructural MRI field to develop and test biophysical models, methods that model temporal brain dynamics, and registration and preprocessing pipelines.

Estimation of free water-corrected microscopic fractional anisotropy

Water diffusion anisotropy MRI is sensitive to microstructural changes in the brain that are hallmarks of various neurological conditions. However, conventional metrics like fractional anisotropy are confounded by neuron fiber orientation dispersion, and the relatively low resolution of diffusion-weighted MRI gives rise to significant free water partial volume effects in many brain regions that are adjacent to cerebrospinal fluid. Microscopic fractional anisotropy is a recent metric that can report water diffusion anisotropy independent of neuron fiber orientation dispersion but is still susceptible to free water contamination. In this paper, we present a free water elimination (FWE) technique to estimate microscopic fractional anisotropy and other related diffusion indices by implementing a signal representation in which the MRI signal within a voxel is assumed to come from two distinct sources: a tissue compartment and a free water compartment. A two-part algorithm is proposed to rapidly fit a set of diffusion-weighted MRI volumes containing both linear- and spherical-tensor encoding acquisitions to the representation. Simulations and in vivo acquisitions with four healthy volunteers indicated that the FWE method may be a feasible technique for measuring microscopic fractional anisotropy and other indices with greater specificity to neural tissue characteristics than conventional methods.

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.