The Larmor frequency shift of a white matter magnetic microstructure model with multiple sources.

Magnetic susceptibility imaging may provide valuable information about chemical composition and microstructural organization of tissue. However, its estimation from the MRI signal phase is particularly difficult as it is sensitive to magnetic tissue properties ranging from the molecular to the macroscopic scale. The MRI Larmor frequency shift measured in white matter (WM) tissue depends on the myelinated axons and other magnetizable sources such as iron-filled ferritin. We have previously derived the Larmor frequency shift arising from a dense medium of cylinders with scalar susceptibility and arbitrary orientation dispersion. Here, we extend our model to include microscopic WM susceptibility anisotropy as well as spherical inclusions with scalar susceptibility to represent subcellular structures, biologically stored iron, and so forth. We validate our analytical results with computer simulations and investigate the feasibility of estimating susceptibility using simple iterative linear least squares without regularization or preconditioning. This is done in a digital brain phantom synthesized from diffusion MRI measurements of an ex vivo mouse brain at ultra-high field.

Rapid microstructural plasticity in the cortical semantic network following a short language learning session.

Despite the clear importance of language in our life, our vital ability to quickly and effectively learn new words and meanings is neurobiologically poorly understood. Conventional knowledge maintains that language learning-especially in adulthood-is slow and laborious. Furthermore, its structural basis remains unclear. Even though behavioural manifestations of learning are evident near instantly, previous neuroimaging work across a range of semantic categories has largely studied neural changes associated with months or years of practice. Here, we address rapid neuroanatomical plasticity accompanying new lexicon acquisition, specifically focussing on the learning of action-related language, which has been linked to the brain's motor systems. Our results show that it is possible to measure and to externally modulate (using transcranial magnetic stimulation (TMS) of motor cortex) cortical microanatomic reorganisation after mere minutes of new word learning. Learning-induced microstructural changes, as measured by diffusion kurtosis imaging (DKI) and machine learning-based analysis, were evident in prefrontal, temporal, and parietal neocortical sites, likely reflecting integrative lexico-semantic processing and formation of new memory circuits immediately during the learning tasks. These results suggest a structural basis for the rapid neocortical word encoding mechanism and reveal the causally interactive relationship of modal and associative brain regions in supporting learning and word acquisition.

MP-PCA denoising of fMRI time-series data can lead to artificial activation "spreading".

MP-PCA denoising has become the method of choice for denoising MRI data since it provides an objective threshold to separate the signal components from unwanted thermal noise components. In rodents, thermal noise in the coils is an important source of noise that can reduce the accuracy of activation mapping in fMRI. Further confounding this problem, vendor data often contains zero-filling and other post-processing steps that may violate MP-PCA assumptions. Here, we develop an approach to denoise vendor data and assess activation "spreading" caused by MP-PCA denoising in rodent task-based fMRI data. Data was obtained from N = 3 mice using conventional multislice and ultrafast fMRI acquisitions (1 s and 50 ms temporal resolution, respectively), using a visual stimulation paradigm. MP-PCA denoising produced SNR gains of 64% and 39%, and Fourier Spectral Amplitude (FSA) increases in BOLD maps of 9% and 7% for multislice and ultrafast data, respectively, when using a small [2 2] denoising window. Larger windows provided higher SNR and FSA gains with increased spatial extent of activation that may or may not represent real activation. Simulations showed that MP-PCA denoising can incur activation "spreading" with increased false positive rate and smoother functional maps due to local "bleeding" of principal components, and that the optimal denoising window for improved specificity of functional mapping, based on Dice score calculations, depends on the data's tSNR and functional CNR. This "spreading" effect applies also to another recently proposed low-rank denoising method (NORDIC), although to a lesser degree. Our results bode well for enhancing spatial and/or temporal resolution in future fMRI work, while taking into account the sensitivity/specificity trade-offs of low-rank denoising methods.

Tensor denoising of multidimensional MRI data.

PURPOSE: To develop a denoising strategy leveraging redundancy in high-dimensional data. THEORY AND METHODS: The SNR fundamentally limits the information accessible by MRI. This limitation has been addressed by a host of denoising techniques, recently including the so-called MPPCA: principal component analysis of the signal followed by automated rank estimation, exploiting the Marchenko-Pastur distribution of noise singular values. Operating on matrices comprised of data patches, this popular approach objectively identifies noise components and, ideally, allows noise to be removed without introducing artifacts such as image blurring, or nonlocal averaging. The MPPCA rank estimation, however, relies on a large number of noise singular values relative to the number of signal components to avoid such ill effects. This condition is unlikely to be met when data patches and therefore matrices are small, for example due to spatially varying noise. Here, we introduce tensor MPPCA (tMPPCA) for the purpose of denoising multidimensional data, such as from multicontrast acquisitions. Rather than combining dimensions in matrices, tMPPCA uses each dimension of the multidimensional data's inherent tensor-structure to better characterize noise, and to recursively estimate signal components. RESULTS: Relative to matrix-based MPPCA, tMPPCA requires no additional assumptions, and comparing the two in a numerical phantom and a multi-TE diffusion MRI data set, tMPPCA dramatically improves denoising performance. This is particularly true for small data patches, suggesting that tMPPCA can be especially beneficial in such cases. CONCLUSIONS: The MPPCA denoising technique can be extended to high-dimensional data with improved performance for smaller patch sizes.

To mask or not to mask? Investigating the impact of accounting for spatial frequency distributions and susceptibility sources on QSM quality.

PURPOSE: Estimating magnetic susceptibility using MRI depends on inverting a forward relationship between the susceptibility and measured Larmor frequency. However, an often-overlooked constraint in susceptibility fitting is that the Larmor frequency is only measured inside the sample, and after successful background field removal, susceptibility sources should only reside inside the same sample. Here, we test the impact of accounting for these constraints in susceptibility fitting. THEORY AND METHODS: Two different digital brain phantoms with scalar susceptibility were examined. We used the MEDI phantom, a simple phantom with no background fields, to examine the effect of the imposed constraints for various levels of SNR. Next, we considered the QSM reconstruction challenge 2.0 phantom with and without background fields. We estimated the parameter accuracy of openly-available QSM algorithms by comparing fitting results to the ground truth. Next, we implemented the mentioned constraints and compared to the standard approach. RESULTS: Including the spatial distribution of frequencies and susceptibility sources decreased the RMS-error compared to standard QSM on both brain phantoms when background fields were absent. When background field removal was unsuccessful, as is presumably the case in most in vivo conditions, it is better to allow sources outside the brain. CONCLUSION: Informing QSM algorithms about the location of susceptibility sources and where Larmor frequency was measured improves susceptibility fitting for realistic SNR levels and efficient background field removal. However, the latter remains the bottleneck of the algorithm. Allowing for external sources regularizes unsuccessful background field removal and is currently the best strategy in vivo.

A user independent denoising method for x-nuclei MRI and MRS.

PURPOSE: X-nuclei (also called non-proton MRI) MRI and spectroscopy are limited by the intrinsic low SNR as compared to conventional proton imaging. Clinical translation of x-nuclei examination warrants the need of a robust and versatile tool improving image quality for diagnostic use. In this work, we compare a novel denoising method with fewer inputs to the current state-of-the-art denoising method. METHODS: Denoising approaches were compared on human acquisitions of sodium (23 Na) brain, deuterium (2 H) brain, carbon (13 C) heart and brain, and simulated dynamic hyperpolarized 13 C brain scans, with and without additional noise. The current state-of-the-art denoising method Global-local higher order singular value decomposition (GL-HOSVD) was compared to the few-input method tensor Marchenko-Pastur principal component analysis (tMPPCA). Noise-removal was quantified by residual distributions, and statistical analyses evaluated the differences in mean-square-error and Bland-Altman analysis to quantify agreement between original and denoised results of noise-added data. RESULTS: GL-HOSVD and tMPPCA showed similar performance for the variety of x-nuclei data analyzed in this work, with tMPPCA removing ˜5% more noise on average over GL-HOSVD. The mean ratio between noise-added and denoising reproducibility coefficients of the Bland-Altman analysis when compared to the original are also similar for the two methods with 3.09 ± 1.03 and 2.83 ± 0.79 for GL-HOSVD and tMPPCA, respectively. CONCLUSION: The strength of tMPPCA lies in the few-input approach, which generalizes well to different data sources. This makes the use of tMPPCA denoising a robust and versatile tool in x-nuclei imaging improvements and the preferred denoising method.

Larmor frequency shift from magnetized cylinders with arbitrary orientation distribution.

The magnetic susceptibility of tissue can provide valuable information about its chemical composition and microstructural organization. However, the relation between the magnetic microstructure and the measurable Larmor frequency shift is understood only for a few idealized cases. Here we analyze the microstructure formed by magnetized, NMR-invisible infinite cylinders suspended in an NMR-reporting fluid. Through simulations, we scrutinize various geometries of mesoscopic Lorentz cavities and inclusions, and show that the cavity size should be approximately one order of magnitude larger than the width of the inclusions. We also analytically derive the Larmor frequency shift for a population of cylinders with arbitrary orientation dispersion and show that it is determined by the l = 2 Laplace expansion coefficients p 2 m of the cylinders' orientation distribution function. Our work underscores the need to account for microstructural organization when estimating magnetic tissue properties.

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.

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.

Beyond the diffusion standard model in fixed rat spinal cord with combined linear and planar encoding.

Information about tissue on the microscopic and mesoscopic scales can be accessed by modelling diffusion MRI signals, with the aim of extracting microstructure-specific biomarkers. The standard model (SM) of diffusion, currently the most broadly adopted microstructural model, describes diffusion in white matter (WM) tissues by two Gaussian components, one of which has zero radial diffusivity, to represent diffusion in intra- and extra-axonal water, respectively. Here, we reappraise these SM assumptions by collecting comprehensive double diffusion encoded (DDE) MRI data with both linear and planar encodings, which was recently shown to substantially enhance the ability to estimate SM parameters. We find however, that the SM is unable to account for data recorded in fixed rat spinal cord at an ultrahigh field of 16.4 T, suggesting that its underlying assumptions are violated in our experimental data. We offer three model extensions to mitigate this problem: first, we generalize the SM to accommodate finite radii (axons) by releasing the constraint of zero radial diffusivity in the intra-axonal compartment. Second, we include intracompartmental kurtosis to account for non-Gaussian behaviour. Third, we introduce an additional (third) compartment. The ability of these models to account for our experimental data are compared based on parameter feasibility and Bayesian information criterion. Our analysis identifies the three-compartment description as the optimal model. The third compartment exhibits slow diffusion with a minor but non-negligible signal fraction (∼12%). We demonstrate how failure to take the presence of such a compartment into account severely misguides inferences about WM microstructure. Our findings bear significance for microstructural modelling at large and can impact the interpretation of biomarkers extracted from the standard model of diffusion.

Microstructural changes in the brain after long-term post-concussion symptoms: A randomized trial.

A recent randomized controlled trial in young patients with long-term post-concussion symptoms showed that a novel behavioral intervention "Get going After concussIoN" is superior to enhanced usual care in terms of symptom reduction. It is unknown whether these interventional effects are associated with microstructural brain changes. The aim of this study was to examine whether diffusion-weighted MRI indices, which are sensitive to the interactions between cellular structures and water molecules' Brownian motion, respond differently to the interventions of the above-mentioned trial and whether such differences correlate with the improvement of post-concussion symptoms. Twenty-three patients from the intervention group (mean age 22.8, 18 females) and 19 patients from the control group (enhanced usual care) (mean age 23.9, 14 females) were enrolled. The primary outcome measure was the mean kurtosis tensor, which is sensitive to the microscopic complexity of brain tissue. The mean kurtosis tensor was significantly increased in the intervention group (p = 0.003) in the corpus callosum but not in the thalamus (p = 0.78) and the hippocampus (p = 0.34). An increase in mean kurtosis tensor in the corpus callosum tended to be associated with a reduction in symptoms, but this association did not reach significance (p = 0.059). Changes in diffusion tensor imaging metrics did not differ between intervention groups and were not associated with symptoms. The current study found different diffusion-weighted MRI responses from the microscopic cellular structures of the corpus callosum between patients receiving a novel behavioral intervention and patients receiving enhanced usual care. Correlations with improvement of post-concussion symptoms were not evident.

Evidence for microscopic kurtosis in neural tissue revealed by correlation tensor MRI.

PURPOSE: The impact of microscopic diffusional kurtosis (µK), arising from restricted diffusion and/or structural disorder, remains a controversial issue in contemporary diffusion MRI (dMRI). Recently, correlation tensor imaging (CTI) was introduced to disentangle the sources contributing to diffusional kurtosis, without relying on a-priori multi-gaussian component (MGC) or other microstructural assumptions. Here, we investigated µK in in vivo rat brains and assessed its impact on state-of-the-art methods ignoring µK. THEORY AND METHODS: CTI harnesses double diffusion encoding (DDE) experiments, which were here improved for speed and minimal bias using four different sets of acquisition parameters. The robustness of the improved CTI protocol was assessed via simulations. In vivo CTI acquisitions were performed in healthy rat brains using a 9.4T pre-clinical scanner equipped with a cryogenic coil, and targeted the estimation of µK, anisotropic kurtosis, and isotropic kurtosis. RESULTS: The improved CTI acquisition scheme substantially reduces scan time and importantly, also minimizes higher-order-term biases, thus enabling robust µK estimation, alongside Kaniso and Kiso metrics. Our CTI experiments revealed positive µK both in white and gray matter of the rat brain in vivo; µK is the dominant kurtosis source in healthy gray matter tissue. The non-negligible µK substantially were found to bias prior MGC analyses of Kiso and Kaniso . CONCLUSIONS: Correlation Tensor MRI offers a more accurate and robust characterization of kurtosis sources than its predecessors. µK is non-negligible in vivo in healthy white and gray matter tissues and could be an important biomarker for future studies. Our findings thus have both theoretical and practical implications for future dMRI research.

Toward more robust and reproducible diffusion kurtosis imaging.

PURPOSE: The general utility of diffusion kurtosis imaging (DKI) is challenged by its poor robustness to imaging artifacts and thermal noise that often lead to implausible kurtosis values. THEORY AND METHODS: A robust scalar kurtosis index can be estimated from powder-averaged diffusion-weighted data. We introduce a novel DKI estimator that uses this scalar kurtosis index as a proxy for the mean kurtosis to regularize the fit. RESULTS: The regularized DKI estimator improves the robustness and reproducibility of the kurtosis metrics and results in parameter maps with enhanced quality and contrast. CONCLUSION: Our novel DKI estimator promotes the wider use of DKI in clinical research and potentially diagnostics by improving the reproducibility and precision of DKI fitting and, as such, enabling enhanced visual, quantitative, and statistical analyses of DKI parameters.

Correlation tensor magnetic resonance imaging.

Diffusional Kurtosis Magnetic Resonance Imaging (DKI) quantifies the extent of non-Gaussian water diffusion, which has been shown to be a sensitive biomarker for microstructure in health and disease. However, DKI is not specific to any microstructural property per se since kurtosis may emerge from several different sources. Q-space trajectory encoding schemes have been proposed for decoupling kurtosis arising from the variance of mean diffusivities (isotropic kurtosis) from kurtosis driven by microscopic anisotropy (anisotropic kurtosis). Still, these methods assume that the system is comprised of multiple Gaussian diffusion components with vanishing intra-compartmental kurtosis (associated with restricted diffusion). Here, we develop a more general framework for resolving the underlying kurtosis sources without relying on the multiple Gaussian diffusion approximation. We introduce Correlation Tensor MRI (CTI) - an approach harnessing the versatility of double diffusion encoding (DDE) and its sensitivity to displacement correlation tensors capable of explicitly decoupling isotropic and anisotropic kurtosis components from intra-compartmental kurtosis effects arising from restricted (and time-dependent) diffusion. Additionally, we show that, by subtracting these isotropic and anisotropic kurtosis components from the total diffusional kurtosis, 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 experiments in ex vivo mouse brains at ultrahigh field of 16.4 T, 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 synthetic microenvironments show, however, that our current CTI protocol for estimating intra-compartmental kurtosis is limited by higher order terms that were not taken into account in this study. CTI measurements were then extended to in vivo settings and used to map heathy rat brains at 9.4 T. These in vivo CTI results were found to be consistent with our ex vivo findings. Although future studies are still required to assess and mitigate the higher order effects on the intra-compartmental kurtosis index, our results show that CTI's more general estimates of anisotropic and isotropic kurtosis contributions are already ripe for future in vivo studies, which can have significant impact our understanding of the mechanisms underlying diffusion metrics extracted in health and disease.

The impact of realistic axonal shape on axon diameter estimation using diffusion MRI.

To study axonal microstructure with diffusion MRI, axons are typically modeled as straight impermeable cylinders, whereby the transverse diffusion MRI signal can be made sensitive to the cylinder's inner diameter. However, the shape of a real axon varies along the axon direction, which couples the longitudinal and transverse diffusion of the overall axon direction. Here we develop a theory of the intra-axonal diffusion MRI signal based on coarse-graining of the axonal shape by 3-dimensional diffusion. We demonstrate how the estimate of the inner diameter is confounded by the diameter variations (beading), and by the local variations in direction (undulations) along the axon. We analytically relate diffusion MRI metrics, such as time-dependent radial diffusivity D⊥(t)and kurtosis K⊥(t),to the axonal shape, and validate our theory using Monte Carlo simulations in synthetic undulating axons with randomly positioned beads, and in realistic axons reconstructed from electron microscopy images of mouse brain white matter. We show that (i) In the narrow pulse limit, the inner diameter from D⊥(t)is overestimated by about twofold due to a combination of axon caliber variations and undulations (each contributing a comparable effect size); (ii) The narrow-pulse kurtosis K⊥|t→∞deviates from that in an ideal cylinder due to caliber variations; we also numerically calculate the fourth-order cumulant for an ideal cylinder in the wide pulse limit, which is relevant for inner diameter overestimation; (iii) In the wide pulse limit, the axon diameter overestimation is mainly due to undulations at low diffusion weightings b; and (iv) The effect of undulations can be considerably reduced by directional averaging of high-b signals, with the apparent inner diameter given by a combination of the axon caliber (dominated by the thickest axons), caliber variations, and the residual contribution of undulations.

Diffusion Kurtosis Imaging maps neural damage in the EAE model of multiple sclerosis.

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 associated 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 with varied grades of 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 reveal a strong relationship between disability and measured microstructural parameters in normal appearing white matter and gray matter. Relationships between disability and mean of the kurtosis tensor, radial kurtosis, radial diffusivity were similar to what has been found in other hypomyelinating MS models, and in patients. However, the changes 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.

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.

Resolving degeneracy in diffusion MRI biophysical model parameter estimation using double diffusion encoding.

PURPOSE: Biophysical tissue models are increasingly used in the interpretation of diffusion MRI (dMRI) data, with the potential to provide specific biomarkers of brain microstructural changes. However, it has been shown recently that, in the general Standard Model, parameter estimation from dMRI data is ill-conditioned even when very high b-values are applied. We analyze this issue for the Neurite Orientation Dispersion and Density Imaging with Diffusivity Assessment (NODDIDA) model and demonstrate that its extension from single diffusion encoding (SDE) to double diffusion encoding (DDE) resolves the ill-posedness for intermediate diffusion weightings, producing an increase in accuracy and precision of the parameter estimation. METHODS: We analyze theoretically the cumulant expansion up to fourth order in b of SDE and DDE signals. Additionally, we perform in silico experiments to compare SDE and DDE capabilities under similar noise conditions. RESULTS: We prove analytically that DDE provides invariant information non-accessible from SDE, which makes the NODDIDA parameter estimation injective. The in silico experiments show that DDE reduces the bias and mean square error of the estimation along the whole feasible region of 5D model parameter space. CONCLUSIONS: DDE adds additional information for estimating the model parameters, unexplored by SDE. We show, as an example, that this is sufficient to solve the previously reported degeneracies in the NODDIDA model parameter estimation.

Evaluation of principal component analysis image denoising on multi-exponential MRI relaxometry.

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 ≈3×, 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.

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.