On the potential for mapping apparent neural soma density via a clinically viable diffusion MRI protocol

Improving Microstructural Estimation in Time-Dependent Diffusion MRI With a Bayesian Method

BACKGROUND

Accurately fitting diffusion-time-dependent diffusion MRI (td-dMRI) models poses challenges due to complex and nonlinear formulas, signal noise, and limited clinical data acquisition.

PURPOSE

Introduce a Bayesian methodology to refine microstructural fitting within the IMPULSED (Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion) model and optimize the prior distribution within the Bayesian framework.

STUDY TYPE

Retrospective.

POPULATION

Involving 69 pediatric patients (median age 6 years, interquartile range [IQR] 3-9 years, 61% male) with 41 low-grade and 28 high-grade gliomas, of which 76.8% were identified within the brainstem or cerebellum.

FIELD STRENGTH/SEQUENCE

3 T, oscillating gradient spin-echo (OGSE) and pulsed gradient spin-echo (PGSE).

ASSESSMENT

The Bayesian method's performance in fitting cell diameter (d $$ d $$ ), intracellular volume fraction (f in $$ {f}_{in} $$ ), and extracellular diffusion coefficient (D ex $$ {D}_{ex} $$ ) was compared against the NLLS method, considering simulated and experimental data. The tumor region-of-interest (ROI) were manually delineated on the b0 images. The diagnostic performance in distinguishing high- and low-grade gliomas was assessed, and fitting accuracy was validated against H&E-stained pathology.

STATISTICAL TESTS

T-test, receiver operating curve (ROC), area under the curve (AUC) and DeLong's test were conducted. Significance considered at P < 0.05.

RESULTS

Bayesian methodology manifested increased accuracy with robust estimates in simulation (RMSE decreased by 29.6%, 40.9%, 13.6%, and STD decreased by 29.2%, 43.5%, and 24.0%, respectively ford $$ d $$ ,f in $$ {f}_{in} $$ , andD ex $$ {D}_{ex} $$ compared to NLLS), indicating fewer outliers and reduced error. Diagnostic performance for tumor grade was similar in both methods, however, Bayesian method generated smoother microstructural maps (outliers ratio decreased by 45.3% ± 19.4%) and a marginal enhancement in correlation with H&E staining result (r = 0.721 forf in $$ {f}_{in} $$ compared to r = 0.698 using NLLS, P = 0.5764).

DATA CONCLUSION

The proposed Bayesian method substantially enhances the accuracy and robustness of IMPULSED model estimation, suggesting its potential clinical utility in characterizing cellular microstructure.

EVIDENCE LEVEL

3 TECHNICAL EFFICACY: Stage 1.

Fetal brain growth and infant autistic traits

BACKGROUND

Structural differences exist in the brains of autistic individuals. To date only a few studies have explored the relationship between fetal brain growth and later infant autistic traits, and some have used fetal head circumference (HC) as a proxy for brain development. These findings have been inconsistent. Here we investigate whether fetal subregional brain measurements correlate with autistic traits in toddlers.

METHODS

A total of 219 singleton pregnancies (104 males and 115 females) were recruited at the Rosie Hospital, Cambridge, UK. 2D ultrasound was performed at 12-, 20- and between 26 and 30 weeks of pregnancy, measuring head circumference (HC), ventricular atrium (VA) and transcerebellar diameter (TCD). A total of 179 infants were followed up at 18-20 months of age and completed the quantitative checklist for autism in toddlers (Q-CHAT) to measure autistic traits.

RESULTS

Q-CHAT scores at 18-20 months of age were positively associated with TCD size at 20 weeks and with HC at 28 weeks, in univariate analyses, and in multiple regression models which controlled for sex, maternal age and birth weight.

LIMITATIONS

Due to the nature and location of the study, ascertainment bias could also have contributed to the recruitment of volunteer mothers with a higher than typical range of autistic traits and/or with a significant interest in the neurodevelopment of their children.

CONCLUSION

Prenatal brain growth is associated with toddler autistic traits and this can be ascertained via ultrasound starting at 20 weeks gestation.

Differentiation of white matter histopathology using b-tensor encoding and machine learning

Diffusion-weighted magnetic resonance imaging (DW-MRI) is a non-invasive technique that is sensitive to microstructural geometry in neural tissue and is useful for the detection of neuropathology in research and clinical settings. Tensor-valued diffusion encoding schemes (b-tensor) have been developed to enrich the microstructural data that can be obtained through DW-MRI. These advanced methods have proven to be more specific to microstructural properties than conventional DW-MRI acquisitions. Additionally, machine learning methods are particularly useful for the study of multidimensional data sets. In this work, we have tested the reach of b-tensor encoding data analyses with machine learning in different histopathological scenarios. We achieved this in three steps: 1) We induced different levels of white matter damage in rodent optic nerves. 2) We obtained ex vivo DW-MRI data with b-tensor encoding schemes and calculated quantitative metrics using Q-space trajectory imaging. 3) We used a machine learning model to identify the main contributing features and built a voxel-wise probabilistic classification map of histological damage. Our results show that this model is sensitive to characteristics of microstructural damage. In conclusion, b-tensor encoded DW-MRI data analyzed with machine learning methods, have the potential to be further developed for the detection of histopathology and neurodegeneration.

Microstructural Changes in the Corpus Callosum in Systemic Lupus Erythematous

Central nervous system (CNS) involvement in childhood-onset systemic lupus erythematosus (cSLE) occurs in more than 50% of patients. Structural magnetic resonance imaging (MRI) has identified global cerebral atrophy, as well as the involvement of the corpus callosum and hippocampus, which is associated with cognitive impairment. In this cross-sectional study we included 71 cSLE (mean age 24.7 years (SD 4.6) patients and a disease duration of 11.8 years (SD 4.8) and two control groups: (1) 49 adult-onset SLE (aSLE) patients (mean age of 33.2 (SD 3.7) with a similar disease duration and (2) 58 healthy control patients (mean age of 29.9 years (DP 4.1)) of a similar age. All of the individuals were evaluated on the day of the MRI scan (Phillips 3T scanner). We reviewed medical charts to obtain the clinical and immunological features and treatment history of the SLE patients. Segmentation of the corpus callosum was performed through an automated segmentation method. Patients with cSLE had a similar mid-sagittal area of the corpus callosum in comparison to the aSLE patients. When compared to the control groups, cSLE and aSLE had a significant reduction in the mid-sagittal area in the posterior region of the corpus callosum. We observed significantly lower FA values and significantly higher MD, RD, and AD values in the total area of the corpus callosum and in the parcels B, C, D, and E in cSLE patients when compared to the aSLE patients. Low complement, the presence of anticardiolipin antibodies, and cognitive impairment were associated with microstructural changes. In conclusion, we observed greater microstructural changes in the corpus callosum in adults with cSLE when compared to those with aSLE. Longitudinal studies are necessary to follow these changes, however they may explain the worse cognitive function and disability observed in adults with cSLE when compared to aSLE.

IMPULSED model based cytological feature estimation with U-Net: Application to human brain tumor at 3T

PURPOSE

This work introduces and validates a deep-learning-based fitting method, which can rapidly provide accurate and robust estimation of cytological features of brain tumor based on the IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) model fitting with diffusion-weighted MRI data.

METHODS

The U-Net was applied to rapidly quantify extracellular diffusion coefficient (D_ex ), cell size (d), and intracellular volume fraction (v_in ) of brain tumor. At the training stage, the image-based training data, synthesized by randomizing quantifiable microstructural parameters within specific ranges, was used to train U-Net. At the test stage, the pre-trained U-Net was applied to estimate the microstructural parameters from simulated data and the in vivo data acquired on patients at 3T. The U-Net was compared with conventional non-linear least-squares (NLLS) fitting in simulations in terms of estimation accuracy and precision.

RESULTS

Our results confirm that the proposed method yields better fidelity in simulations and is more robust to noise than the NLLS fitting. For in vivo data, the U-Net yields obvious quality improvement in parameter maps, and the estimations of all parameters are in good agreement with the NLLS fitting. Moreover, our method is several orders of magnitude faster than the NLLS fitting (from about 5 min to <1 s).

CONCLUSION

The image-based training scheme proposed herein helps to improve the quality of the estimated parameters. Our deep-learning-based fitting method can estimate the cell microstructural parameters fast and accurately.

Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange

Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging anisotropic compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along neuronal and glial processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Kärger model of two exchanging compartments. Using datasets acquired at multiple diffusion weightings (b) and diffusion times (t) in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the gray matter, compared to the other two possible contributions. Our results for the diffusion time window 20-45 ms show minimal diffusivity time-dependence and more pronounced kurtosis decay with time, which is well fit by the exchange model. Moreover, we observe lower signal for longer diffusion times at high b. In light of these observations, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for gray matter microstructure mapping. We finally highlight multi-b multi-t acquisition protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 15 - 60 ms in the rat cortex and hippocampus in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions. This suggests water exchange as an essential component for interpreting diffusion MRI measurements in gray matter.

Soma and Neurite Density MRI (SANDI) of the in-vivo mouse brain and comparison with the Allen Brain Atlas

Diffusion MRI (dMRI) provides unique insights into the neural tissue milieu by probing interactions between diffusing molecules and tissue microstructure. Most dMRI techniques focus on white matter (WM) tissues, nevertheless, interest in gray matter characterizations is growing. The Soma and Neurite Density MRI (SANDI) methodology harnesses a model incorporating water diffusion in spherical objects (assumed to be associated with cell bodies) and in impermeable "sticks" (assumed to represent neurites), which potentially enables the characterization of cellular and neurite densities. Recognising the importance of rodents in animal models of development, aging, plasticity, and disease, we here employ SANDI for in-vivo preclinical imaging and provide a first validation of the methodology by comparing SANDI metrics with cellular density reflected by the Allen mouse brain atlas. SANDI was implemented on a 9.4T scanner equipped with a cryogenic coil, and in-vivo experiments were carried out on N = 6 mice. Pixelwise, ROI-based, and atlas comparisons were performed, magnitude vs. real-valued analyses were compared, and shorter acquisitions with reduced the number of b-value shells were investigated. Our findings reveal good reproducibility of the SANDI parameters, including the sphere and stick fractions, as well as sphere size (CoV < 7%, 12% and 3%, respectively). Additionally, we find a very good rank correlation between SANDI-driven sphere fraction and Allen mouse brain atlas contrast that represents cellular density. We conclude that SANDI is a viable preclinical MRI technique that can greatly contribute to research on brain tissue microstructure.

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.

Training data distribution significantly impacts the estimation of tissue microstructure with machine learning

PURPOSE

Supervised machine learning (ML) provides a compelling alternative to traditional model fitting for parameter mapping in quantitative MRI. The aim of this work is to demonstrate and quantify the effect of different training data distributions on the accuracy and precision of parameter estimates when supervised ML is used for fitting.

METHODS

We fit a two- and three-compartment biophysical model to diffusion measurements from in-vivo human brain, as well as simulated diffusion data, using both traditional model fitting and supervised ML. For supervised ML, we train several artificial neural networks, as well as random forest regressors, on different distributions of ground truth parameters. We compare the accuracy and precision of parameter estimates obtained from the different estimation approaches using synthetic test data.

RESULTS

When the distribution of parameter combinations in the training set matches those observed in healthy human data sets, we observe high precision, but inaccurate estimates for atypical parameter combinations. In contrast, when training data is sampled uniformly from the entire plausible parameter space, estimates tend to be more accurate for atypical parameter combinations but may have lower precision for typical parameter combinations.

CONCLUSION

This work highlights that estimation of model parameters using supervised ML depends strongly on the training-set distribution. We show that high precision obtained using ML may mask strong bias, and visual assessment of the parameter maps is not sufficient for evaluating the quality of the estimates.