Search

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

Diffusional Kurtosis Imaging (DKI) is a sensitive biomarker for microstructure in health and disease. However, DKI is not specific to any microstructural property since it may emerge from several different sources. Q-space trajectory encoding has been proposed for decoupling isotropic from anisotropic kurtosis. Still, this method assume that the system is comprised of multiple Gaussian diffusion components. Here, we develop a more general framework for resolving the underlying kurtosis sources. We introduce Correlation Tensor MRI (CTI) - an approach harnessing the versatility of double diffusion encoding (DDE) and capable of explicitly decoupling isotropic and anisotropic kurtosis components from intra-compartmental kurtosis effects arising from restricted diffusion. Additionally, CTI provides an index that is potentially sensitive to intra-compartmental kurtosis. The theoretical foundations of CTI, as well as the first proof-of-concept CTI ex vivo experiments in mouse brains at a field of 16.4T, are presented. We find that anisotropic and isotropic kurtosis can decouple microscopic anisotropy from substantial partial volume effects between tissue and free water. Our intra-compartmental kurtosis index exhibited positive values in both white and grey matter tissues. Simulations in different microenvironments show, however, that our current CTI protocol for estimating intra-compartmental kurtosis is limited by higher-order terms that were not considered in this study. CTI measurements were then extended to in vivo settings and used to map healthy rat brains at 9.4T. These in vivo CTI results were consistent with our ex vivo findings. Although future studies are still required to mitigate the higher-order effects on the intra-compartmental kurtosis index, results show that CTI's more general estimates of anisotropic and isotropic kurtosis contributions are already ripe for future studies., Comment: 37 pages (main text), 8 main figures + 4 supplementary figures

This work introduces a compartment-based model for apparent soma and neurite density imaging (SANDI) using non-invasive diffusion-weighted MRI (DW-MRI). The existing conjecture in brain microstructure imaging trough DW-MRI presents water diffusion in white (WM) and grey (GM) matter as restricted diffusion in neurites, modelled by infinite cylinders of null radius embedded in the hindered extra-neurite water. The extra-neurite pool in WM corresponds to water in the extra-axonal space, but in GM it combines water in the extra-cellular space with water in soma. While several studies showed that this microstructure model successfully describe DW-MRI data in WM and GM at b<3 ms/{\mum^2}, it has been also shown to fail in GM at high b values (b>>3 ms/{\mum^2}). Here we hypothesize that the unmodelled soma compartment may be responsible for this failure and propose SANDI as a new model of brain microstructure where soma (i.e. cell body of any brain cell type: from neuroglia to neurons) is explicitly included. We assess the effects of size and density of soma on the direction-averaged DW-MRI signal at high b values and the regime of validity of the model using numerical simulations and comparison with experimental data from mouse (bmax = 40 ms/{/mum^2}) and human (bmax = 10 ms/{\mum^2}) brain. We show that SANDI defines new contrasts representing new complementary information on the brain cyto- and myelo-architecture. Indeed, we show for the first-time maps from 25 healthy human subjects of MR soma and neurite signal fractions, that remarkably mirror contrasts of histological images of brain cyto- and myelo-architecture. Although still under validation, SANDI might provide new insight into tissue architecture by introducing a new set of biomarkers of potential great value for biomedical applications and pure neuroscience.

Purpose: Most MR spectroscopy (MRS) pulse sequences rely on broadband excitation with water saturation and typically focus on upfield signals. By contrast, the downfield spectrum, which contains many potentially useful resonances, is typically not targeted because conventional water-suppressed techniques indirectly saturate the labile protons through exchange. Relaxation-enhanced MRS (RE-MRS) uses frequency-selective excitation while actively avoiding bulk water perturbation, thereby enabling high-quality downfield spectroscopy. However, RE-MRS typically requires very long (typically >40 ms) echo times (TEs) due to its localization module, which inevitably decreases sensitivity and filters shorter T2 components. Here, we overcome this limitation by combining RE-MRS and image selected in vivo spectroscopy (ISIS) localization, abbreviated iRE-MRS, which in turn allows very short TEs (5 ms using our hardware). Methods: Experiments were performed in vitro for validation as well as and in in vivo rat brains at 9.4T. Results: The new iRE-MRS methodology was validated in phantoms where good performance was noted. When the downfield spectrum was investigated at short TEs in in vivo rat brains, iRE-MRS provided very high sensitivity; the ensuing downfield spectra encompassed numerous broad peaks, as well as a broad baseline. All downfield spectral peaks were highly attenuated by increasing TEs as well as by applying water saturation, although to different extent. The signal ratios also varied between TEs, suggesting that exchange rates are different among the downfield signals. Conclusions: Short-TE iRE 1H downfield MRS opens new directions in the investigation of in vivo downfield metabolites and their role on healthy and disease processes.

Purpose: Microscopic fractional anisotropy ({\mu}FA) can disentangle microstructural information from orientation dispersion. While double diffusion encoding (DDE) MRI methods are widely used to extract accurate {\mu}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 {\mu}FA estimation. This study aims to evaluate {\mu}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 ({\Delta}/{\delta} = 12/1.5 ms). The {\mu}FA maps obtained from powder-averaged SDE signals were then compared to maps obtained from DDE-MRI experiments ({\Delta}/{\tau}/{\delta} = 12/12/1.5 ms), which allow a model-free estimation of {\mu}FA. Theory and simulations that consider different types of heterogeneity are presented for corroborating the experimental findings. Results: {\mu}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 {\mu}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.

Investigating neural activity from a global brain perspective in-vivo has been in the domain of functional Magnetic Resonance Imaging (fMRI) over the past few decades. The intricate neurovascular couplings that govern fMRI's blood-oxygenation-level-dependent (BOLD) functional contrast are invaluable in mapping active brain regions, but they also entail significant limitations, such as non-specificity of the signal to active foci. Diffusion-weighted functional MRI (dfMRI) with relatively high diffusion-weighting strives to ameliorate this shortcoming as it offers functional contrasts more intimately linked with the underlying activity. Insofar, apart from somewhat smaller activation foci, dfMRI's contrasts have not been convincingly shown to offer significant advantages over BOLD, and its contrasts relied on significant modelling. Here, we study whether dfMRI could offer a better representation of neural activity in the thalamocortical pathway compared to its (spin-echo (SE)) BOLD counterpart. Using high-end forepaw stimulation experiments in the rat at 9.4 T, and with significant sensitivity enhancements due to the use of cryocoils, we show for the first time that dfMRI signals exhibit layer specificity, and, additionally, display signals in areas devoid of SE-BOLD responses. We find that dfMRI signals in the thalamocortical pathway cohere with each other, namely, dfMRI signals in the ventral posterolateral (VPL) thalamic nucleus cohere specifically with layers IV and V in the somatosensory cortex. These activity patterns are much better correlated (compared with SE-BOLD signals) with literature-based electrophysiological recordings in the cortex as well as thalamus. All these findings suggest that dfMRI signals better represent the underlying neural activity in the pathway. In turn, this may entail significant implications towards a much more specific and accurate (...)

PURPOSE: Multi-exponential relaxometry is a powerful tool for characterizing tissue, but generally requires high image signal-to-noise ratio (SNR). This work evaluates the use of principal-component-analysis (PCA) denoising to mitigate these SNR demands and improve the precision of relaxometry measures. METHODS: PCA denoising was evaluated using both simulated and experimental MRI data. Bi-exponential transverse relaxation signals were simulated for a wide range of acquisition and sample parameters, and experimental data were acquired from three excised and fixed mouse brains. In both cases, standard relaxometry analysis was performed on both original and denoised image data, and resulting estimated signal parameters were compared. RESULTS: Denoising reduced the root-mean-square-error of parameters estimated from multi-exponential relaxometry by factors of 2 to 4x, for typical acquisition and sample parameters. Denoised images and subsequent parameter maps showed little or no signs of spatial artifact or loss of resolution. CONCLUSION: Experimental studies and simulations demonstrate that PCA denoising of MRI relaxometry data is an effective method of improving parameter precision without sacrificing image resolution. This simple yet important processing step thus paves the way for broader applicability of multi-exponential MRI relaxometry.

Purpose: Diffusion Magnetic Resonance Imaging (dMRI) is confounded by its long acquisition duration, thereby thwarting the detection of rapid microstructural changes, especially when diffusivity variations are accompanied by rapid changes in T2. The purpose of the present study is to accelerate dMRI to a single scan acquisition, and to enable a more accurate estimation of diffusivity as function of time. Methods: A general methodology termed Incomplete Initial Nutation Diffusion Imaging (INDI) capturing two diffusion contrasts in a single shot, is presented. INDI creates a longitudinal magnetization reservoir that facilitates the successive acquisition of two images separated by only a few milliseconds. INDI's theory is presented, followed by proof-of-concept ex- and in-vivo experiments at 16.4 T and 9.4 T. Results: Mean diffusivities (MDs) extracted from INDI were comparable with Diffusion Tensor Imaging (DTI) and the two-shot IDE in the water phantom. As expected in the brain tissues, DTI provided lower MD than UF-IDE and IDE, but IDE and UF-IDE were in excellent agreement. Simulations are presented for identifying the regimes where INDI is most beneficial. Conclusions: INDI accelerates dMRI acquisition to single-shot mode, which can be of great importance for mapping dynamic microstructural properties in-vivo without T2 bias.

The auditory pathway is widely distributed throughout the brain, and is perhaps one of the most interesting networks in the context of neuroplasticity. Accurate mapping of neural activity in the entire pathway, preferably noninvasively, and with high resolution, could be instrumental for understanding such longitudinal processes. Functional magnetic resonance imaging (fMRI) has clear advantages for such characterizations, as it is noninvasive, provides relatively high spatial resolution and lends itself for repetitive studies, albeit relying on an indirect neurovascular coupling to deliver its information. Indeed, fMRI has been previously used to characterize the auditory pathway in humans and in rats. In the mouse, however, the auditory pathway has insofar only been mapped using manganese-enhanced MRI. Here, we describe a novel setup specifically designed for high-resolution mapping of the mouse auditory pathway using high-field fMRI. Robust and consistent Blood-Oxygenation-Level-Dependent (BOLD) responses were documented along nearly the entire auditory pathway, from the cochlear nucleus (CN), through the superior olivary complex (SOC), nuclei of the lateral lemniscus (LL), inferior colliculus (IC) and the medial geniculate body (MGB). By contrast, clear BOLD responses were not observed in auditory cortex (AC) in this study. Diverse BOLD latencies were mapped ROI- and pixel-wise using coherence analysis, evidencing different averaged BOLD time courses in different auditory centers. Some degree of tonotopy was identified in the IC, SOC, and MGB in the pooled dataset though it could not be assessed per subject due to a lack of statistical power. Given the importance of the mouse model in plasticity studies, animal models, and optogenetics, and fMRI's potential to map dynamic responses to specific cues, this first fMRI study of the mouse auditory pathway paves the way for future (...)

Many developmental processes, such as plasticity and aging, or pathological processes such as neurological diseases are characterized by modulations of specific cellular types and their microstructures. Diffusion-weighted Magnetic Resonance Imaging (DW-MRI) is a powerful technique for probing microstructure, yet its information arises from the ubiquitous, non-specific water signal. By contrast, diffusion-weighted Magnetic Resonance Spectroscopy (DW-MRS) allows specific characterizations of tissues such as brain and muscle in vivo by quantifying the diffusion properties of MR-observable metabolites. Many brain metabolites are predominantly intracellular, and some of them are preferentially localized in specific brain cell populations, e.g., neurons and glia. Given the microstructural sensitivity of diffusion-encoding filters, investigation of metabolite diffusion properties using DW-MRS can thus provide exclusive cell and compartment-specific information. Furthermore, since many models and assumptions are used for quantification of water diffusion, metabolite diffusion may serve to generate a-priori information for model selection in DW-MRI. However, DW-MRS measurements are extremely challenging, from the acquisition to the accurate and correct analysis and quantification stages. In this review, we survey the state-of-the-art methods that have been developed for the robust acquisition, quantification and analysis of DW-MRS data and discuss the potential relevance of DW-MRS for elucidating brain microstructure in vivo. The review highlights that when accurate data on the diffusion of multiple metabolites is combined with accurate computational and geometrical modelling, DW-MRS can provide unique cell-specific information on the intracellular structure of brain tissue, in health and disease, which could serve as incentives for further application in vivo in human (...)

Mapping tissue microstructure accurately and noninvasively is one of the frontiers of biomedical imaging. Diffusion Magnetic Resonance Imaging (MRI) is at the forefront of such efforts, as it is capable of reporting on microscopic structures orders of magnitude smaller than the voxel size by probing restricted diffusion. Double Diffusion Encoding (DDE) and Double Oscillating Diffusion Encoding (DODE) in particular, are highly promising for their ability to report on microscopic fractional anisotropy ({\mu}FA), a measure of the pore anisotropy in its own eigenframe, irrespective of orientation distribution. However, the underlying correlates of {\mu}FA have insofar not been studied. Here, we extract {\mu}FA from DDE and DODE measurements at ultrahigh magnetic field of 16.4T in the aim to probe fixed rat spinal cord microstructure. We further endeavor to correlate {\mu}FA with Myelin Water Fraction (MWF) derived from multiexponential T2 relaxometry, as well as with literature-based spatially varying axonal diameters. In addition, a simple new method is presented for extracting unbiased {\mu}FA from three measurements at different b-values. Our findings reveal strong anticorrelations between {\mu}FA (derived from DODE) and axon diameter in the distinct spinal cord tracts; a moderate correlation was also observed between {\mu}FA derived from DODE and MWF. These findings suggest that axonal membranes strongly modulate {\mu}FA, which - owing to its robustness towards orientation dispersion effects - reflects axon diameter much better than its typical FA counterpart. The {\mu}FA exhibited modulations when measured via oscillating or blocked gradients, suggesting selective probing of different parallel path lengths and providing insight into how those modulate {\mu}FA metrics. Our findings thus shed light into the underlying microstructural correlates of {\mu}FA and are (...)

Multiexponential T2 (MET2) Relaxometry and Magnetization Transfer (MT) are among the most promising MRI-derived techniques for white matter (WM) characterization. Both techniques are shown to have histologically correlated sensitivity to myelin, but these correlations are not fully understood. Furthermore, MET2 and MT reports on different features of WM, thus being specific to different (patho)physiological states. Two-dimensional studies, such as those commonly used in NMR, have been rarely performed in this context. Here, off-resonance irradiation effects on MET2 components were evaluated in fixed rat spinal cord white matter at 16.4T. These 2D MT-MET2 experiments reveal that MT affects both short and long T2 in a tract-specific fashion. The spatially distinct modulations enhanced contrast between microstructurally-distinct spinal cord tracts. Two hypotheses to explain these findings were proposed: either selective elimination of a short T2 component through pre-saturation is combined with intercompartmental water exchange effects that occur on the timescale of irradiation; or other macromolecular species that exist within the tissue - other than myelin - such as neurofilaments, may be involved in the apparent microstructural segregation of the spinal cord (SC) from MET2. Though further investigation is required to elucidate the underlying mechanism, this phenomenon adds a new dimension for WM characterization.

Diffusion kurtosis imaging (DKI), is an imaging modality that yields novel disease biomarkers and in combination with nervous tissue modeling, provides access to microstructural parameters. Recently, DKI and subsequent estimation of microstructural model parameters has been used for assessment of tissue changes in neurodegenerative diseases and their animal models. In this study, mouse spinal cords from the experimental autoimmune encephalomyelitis (EAE) model of multiple sclerosis (MS) were investigated for the first time using DKI in combination with biophysical modeling to study the relationship between microstructural metrics and degree of animal dysfunction. Thirteen spinal cords were extracted from animals of variable disability and scanned in a high-field MRI scanner along with five control specimen. Diffusion weighted data were acquired together with high resolution T2* images. Diffusion data were fit to estimate diffusion and kurtosis tensors and white matter modeling parameters, which were all used for subsequent statistical analysis using a linear mixed effects model. T2* images were used to delineate focal demyelination/inflammation. Our results unveil a strong relationship between disability and measured microstructural parameters in normal appearing white matter and gray matter. The changes we found in biophysical modeling parameters and in particular in extra-axonal axial diffusivity were clearly different from previous studies employing other animal models of MS. In conclusion, our data suggest that DKI and microstructural modeling can provide a unique contrast capable of detecting EAE-specific changes correlating with clinical disability. These findings could close the gap between MRI findings and clinical presentation in patients and deepen our understanding of EAE and the MS mechanisms., Comment: 40 pages, 6 figures, 4 tables

Designing novel diffusion-weighted pulse sequences to probe tissue microstructure beyond the conventional Stejskal-Tanner family is currently of broad interest. One such technique, multidimensional diffusion MRI, has been recently proposed to afford model-free decomposition of diffusion signal kurtosis into terms originating from either ensemble variance of isotropic diffusivity or microscopic diffusion anisotropy. This ability rests on the assumption that diffusion can be described as a sum of multiple Gaussian compartments, but this is often not strictly fulfilled. The effects of nongaussian diffusion on single shot isotropic diffusion sequences were first considered in detail by de Swiet and Mitra in 1996. They showed theoretically that anisotropic compartments lead to anisotropic time dependence of the diffusion tensors, which causes the measured isotropic diffusivity to depend on gradient frame orientation. Here we show how such deviations from the multiple Gaussian compartments assumption conflates orientation dispersion with ensemble variance in isotropic diffusivity. Second, we consider additional contributions to the apparent variance in isotropic diffusivity arising due to intracompartmental kurtosis. These will likewise depend on gradient frame orientation. We illustrate the potential importance of these confounds with analytical expressions, numerical simulations in simple model geometries, and microimaging experiments in fixed spinal cord using isotropic diffusion encoding waveforms with 7.5 ms duration and 3000 mT/m maximum amplitude., Comment: 26 pages, 9 figures. Appearing in J. Magn. Reson

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