Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation

Bias-reduced neural networks for parameter estimation in quantitative MRI

PURPOSE

To develop neural network (NN)-based quantitative MRI parameter estimators with minimal bias and a variance close to the Cramér-Rao bound.

THEORY AND METHODS

We generalize the mean squared error loss to control the bias and variance of the NN's estimates, which involves averaging over multiple noise realizations of the same measurements during training. Bias and variance properties of the resulting NNs are studied for two neuroimaging applications.

RESULTS

In simulations, the proposed strategy reduces the estimates' bias throughout parameter space and achieves a variance close to the Cramér-Rao bound. In vivo, we observe good concordance between parameter maps estimated with the proposed NNs and traditional estimators, such as nonlinear least-squares fitting, while state-of-the-art NNs show larger deviations.

CONCLUSION

The proposed NNs have greatly reduced bias compared to those trained using the mean squared error and offer significantly improved computational efficiency over traditional estimators with comparable or better accuracy.

The Diffusion Exchange Ratio (DEXR): A minimal sampling of diffusion exchange spectroscopy to probe exchange, restriction, and time-dependence

Water exchange is increasingly recognized as an important biological process that can affect the study of biological tissue using diffusion MR. Methods to measure exchange, however, remain immature as opposed to those used to characterize restriction, with no consensus on the optimal pulse sequence (s) or signal model (s). In general, the trend has been towards data-intensive fitting of highly parameterized models. We take the opposite approach and show that a judicious sub-sample of diffusion exchange spectroscopy (DEXSY) data can be used to robustly quantify exchange, as well as restriction, in a data-efficient manner. This sampling produces a ratio of two points per mixing time: (i) one point with equal diffusion weighting in both encoding periods, which gives maximal exchange contrast, and (ii) one point with the same total diffusion weighting in just the first encoding period, for normalization. We call this quotient the Diffusion EXchange Ratio (DEXR). Furthermore, we show that it can be used to probe time-dependent diffusion by estimating the velocity autocorrelation function (VACF) over intermediate to long times (∼2-500ms). We provide a comprehensive theoretical framework for the design of DEXR experiments in the case of static or constant gradients. Data from Monte Carlo simulations and experiments acquired in fixed and viable ex vivo neonatal mouse spinal cord using a permanent magnet system are presented to test and validate this approach. In viable spinal cord, we report the following apparent parameters from just 6 data points: τk=17±4ms, fNG=0.72±0.01, Reff=1.05±0.01μm, and κeff=0.19±0.04μm/ms, which correspond to the exchange time, restricted or non-Gaussian signal fraction, an effective spherical radius, and permeability, respectively. For the VACF, we report a long-time, power-law scaling with ≈t-2.4, which is approximately consistent with disordered domains in 3-D. Overall, the DEXR method is shown to be highly efficient, capable of providing valuable quantitative diffusion metrics using minimal MR data.

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.

Exploring white matter dynamics and morphology through interactive numerical phantoms: the White Matter Generator

Brain white matter is a dynamic environment that continuously adapts and reorganizes in response to stimuli and pathological changes. Glial cells, especially, play a key role in tissue repair, inflammation modulation, and neural recovery. The movements of glial cells and changes in their concentrations can influence the surrounding axon morphology. We introduce the White Matter Generator (WMG) tool to enable the study of how axon morphology is influenced through such dynamical processes, and how this, in turn, influences the diffusion-weighted MRI signal. This is made possible by allowing interactive changes to the configuration of the phantom generation throughout the optimization process. The phantoms can consist of myelinated axons, unmyelinated axons, and cell clusters, separated by extra-cellular space. Due to morphological flexibility and computational advantages during the optimization, the tool uses ellipsoids as building blocks for all structures; chains of ellipsoids for axons, and individual ellipsoids for cell clusters. After optimization, the ellipsoid representation can be converted to a mesh representation which can be employed in Monte-Carlo diffusion simulations. This offers an effective method for evaluating tissue microstructure models for diffusion-weighted MRI in controlled bio-mimicking white matter environments. Hence, the WMG offers valuable insights into white matter's adaptive nature and implications for diffusion-weighted MRI microstructure models, and thereby holds the potential to advance clinical diagnosis, treatment, and rehabilitation strategies for various neurological disorders and injuries.

In vivo mapping of cellular resolution neuropathology in brain ischemia with diffusion MRI

Noninvasive mapping of cellular pathology can provide critical diagnostic and prognostic information. Recent advances in diffusion magnetic resonance imaging enabled in vivo examination of tissue microstructures well beyond the imaging resolution. Here, we proposed to use diffusion time-dependent diffusion kurtosis imaging (tDKI) to simultaneously assess cellular morphology and transmembrane permeability in hypoxic-ischemic (HI) brain injury. Through numerical simulations and organoid imaging, we demonstrated the feasibility of capturing effective size and permeability changes using tDKI. In vivo MRI of HI-injured mouse brains detected a shift of the tDKI peak to longer diffusion times, suggesting swelling of the cellular processes. Furthermore, we observed a faster decrease of the tDKI tail, reflecting increased transmembrane permeability associated with up-regulated water exchange or necrosis. Such information, unavailable from a single diffusion time, can predict salvageable tissues. Preliminary applications of tDKI in patients with ischemic stroke suggested increased transmembrane permeability in stroke regions, illustrating tDKI's potential for detecting pathological changes in the clinics.

Tuned exchange imaging: Can the filter exchange imaging pulse sequence be adapted for applications with thin slices and restricted diffusion?

Filter exchange imaging (FEXI) is a double diffusion-encoding (DDE) sequence that is specifically sensitive to exchange between sites with different apparent diffusivities. FEXI uses a diffusion-encoding filtering block followed by a detection block at varying mixing times to map the exchange rate. Long mixing times enhance the sensitivity to exchange, but they pose challenges for imaging applications that require a stimulated echo sequence with crusher gradients. Thin imaging slices require strong crushers, which can introduce significant diffusion weighting and bias exchange rate estimates. Here, we treat the crushers as an additional encoding block and consider FEXI as a triple diffusion-encoding sequence. This allows the bias to be corrected in the case of multi-Gaussian diffusion, but not easily in the presence of restricted diffusion. Our approach addresses challenges in the presence of restricted diffusion and relies on the ability to independently gauge sensitivities to exchange and restricted diffusion for arbitrary gradient waveforms. It follows two principles: (i) the effects of crushers are included in the forward model using signal cumulant expansion; and (ii) timing parameters of diffusion gradients in filter and detection blocks are adjusted to maintain the same level of restriction encoding regardless of the mixing time. This results in the tuned exchange imaging (TEXI) protocol. The accuracy of exchange mapping with TEXI was assessed through Monte Carlo simulations in spheres of identical sizes and gamma-distributed sizes, and in parallel hexagonally packed cylinders. The simulations demonstrate that TEXI provides consistent exchange rates regardless of slice thickness and restriction size, even with strong crushers. However, the accuracy depends on b-values, mixing times, and restriction geometry. The constraints and limitations of TEXI are discussed, including suggestions for protocol adaptations. Further studies are needed to optimize the precision of TEXI and assess the approach experimentally in realistic, heterogeneous substrates.

Investigating apparent differences between standard DKI and axisymmetric DKI and its consequences for biophysical parameter estimates

PURPOSE

The purpose of the study is to identify differences between axisymmetric diffusion kurtosis imaging (DKI) and standard DKI, their consequences for biophysical parameter estimates, and the protocol choice influence on parameter estimation.

METHODS

Noise-free and noisy, synthetic diffusion MRI human brain data is simulated using standard DKI for a standard and the fast "199" acquisition protocol. First the noise-free "baseline" difference between both DKI models is estimated and the influence of fiber complexity is investigated. Noisy data is used to establish the signal-to-noise ratio at which the baseline difference exceeds noise variability. The influence of protocol choices and denoising is investigated. The five axisymmetric DKI tensor metrics (AxTM), the parallel and perpendicular diffusivity and kurtosis and mean of the kurtosis tensor are used to compare both DKI models. Additionally, the baseline difference is also estimated for the five parameters of the WMTI-Watson model.

RESULTS

The parallel and perpendicular kurtosis and all of the WMTI-Watson parameters had large baseline differences. Using a Westin or FA mask reduced the number of voxels with large baseline difference, that is, by selecting voxels with less complex fibers. For the noisy data, precision was worsened by the fast "199" protocol but adaptive denoising can help counteract these effects.

CONCLUSION

For the diffusivities and mean of the kurtosis tensor, axisymmetric DKI with a standard protocol delivers similar results as standard DKI. Fiber complexity is one main driver of the baseline differences. Using the "199" protocol worsens precision in noisy data but adaptive denoising mitigates these effects.

Signatures of microstructure in gradient-echo and spin-echo signals

PURPOSE

To determine whether the spatial scale and magnetic susceptibility of microstructure can be evaluated robustly from the decay of gradient-echo and spin-echo signals.

THEORY AND METHODS

Gradient-echo and spin-echo images were acquired from suspensions of spherical polystyrene microbeads of 10, 20, and 40 μm nominal diameter. The sizes of the beads and their magnetic susceptibility relative to the medium were estimated from the signal decay curves, using a lookup table generated from Monte Carlo simulations and an analytic model based on the Gaussian phase approximation.

RESULTS

Fitting Monte Carlo predictions to spin-echo data yielded acceptable estimates of microstructural parameters for the 20 and 40 μm microbeads. Using gradient-echo data, the Monte Carlo lookup table provided satisfactory parameter estimates for the 20 μm beads but unstable results for the diameter of the largest beads. Neither spin-echo nor gradient-echo data allowed accurate parameter estimation for the smallest beads. The analytic model performed poorly over all bead sizes.

CONCLUSIONS

Microstructural sources of magnetic susceptibility produce distinctive non-exponential signatures in the decay of gradient-echo and spin-echo signals. However, inverting the problem to extract microstructural parameters from the signals is nontrivial and, in certain regimes, ill-conditioned. For microstructure with small characteristic length scales, parameter estimation is hampered by the difficulty of acquiring accurate data at very short echo times. For microstructure with large characteristic lengths, the gradient-echo signal approaches the static-dephasing regime, where it becomes insensitive to size. Applicability of the analytic model was further limited by failure of the Gaussian phase approximation for all but the smallest beads.

Inflammation-Triggered Enlargement of Choroid Plexus in Subacute COVID-19 Patients with Neurological Symptoms

OBJECTIVE

To investigate whether choroid plexus volumes in subacute coronavirus disease 2019 (COVID-19) patients with neurological symptoms could indicate inflammatory activation or barrier dysfunction and assess their association with clinical data.

METHODS

Choroid plexus volumes were measured in 28 subacute COVID-19 patients via cerebral magnetic resonance imaging (MRI), compared with those in infection-triggered non-COVID-19 encephalopathy patients (n = 25), asymptomatic individuals after COVID-19 (n = 21), and healthy controls (n = 21). Associations with inflammatory serum markers (peak counts of leukocytes, C-reactive protein [CRP], interleukin 6), an MRI-based marker of barrier dysfunction (CSF volume fraction [V-CSF]), and clinical parameters like olfactory performance and cognitive scores (Montreal Cognitive Assessment) were investigated.

RESULTS

COVID-19 patients showed significantly larger choroid plexus volumes than control groups (p < 0.001, η2 = 0.172). These volumes correlated significantly with peak leukocyte levels (p = 0.001, Pearson's r = 0.621) and V-CSF (p = 0.009, Spearman's rho = 0.534), but neither with CRP nor interleukin 6. No significant correlations were found with clinical parameters.

INTERPRETATION

In patients with subacute COVID-19, choroid plexus volume is a marker of central nervous system inflammation and barrier dysfunction in the presence of neurologic symptoms. The absence of plexus enlargement in infection-triggered non-COVID-19 encephalopathy suggests a specific severe acute respiratory syndrome coronavirus 2 effect. This study also documents an increase in choroid plexus volume for the first time as a parainfectious event. ANN NEUROL 2024.

Callosal Interhemispheric Communication in Mild Traumatic Brain Injury: A Mediation Analysis on WM Microstructure Effects

BACKGROUND AND PURPOSE

Because the corpus callosum connects the left and right hemispheres and a variety of WM bundles across the brain in complex ways, damage to the neighboring WM microstructure may specifically disrupt interhemispheric communication through the corpus callosum following mild traumatic brain injury. Here we use a mediation framework to investigate how callosal interhemispheric communication is affected by WM microstructure in mild traumatic brain injury.

MATERIALS AND METHODS

Multishell diffusion MR imaging was performed on 23 patients with mild traumatic brain injury within 1 month of injury and 17 healthy controls, deriving 11 diffusion metrics, including DTI, diffusional kurtosis imaging, and compartment-specific standard model parameters. Interhemispheric processing speed was assessed using the interhemispheric speed of processing task (IHSPT) by measuring the latency between word presentation to the 2 hemivisual fields and oral word articulation. Mediation analysis was performed to assess the indirect effect of neighboring WM microstructures on the relationship between the corpus callosum and IHSPT performance. In addition, we conducted a univariate correlation analysis to investigate the direct association between callosal microstructures and IHSPT performance as well as a multivariate regression analysis to jointly evaluate both callosal and neighboring WM microstructures in association with IHSPT scores for each group.

RESULTS

Several significant mediators in the relationships between callosal microstructure and IHSPT performance were found in healthy controls. However, patients with mild traumatic brain injury appeared to lose such normal associations when microstructural changes occurred compared with healthy controls.

CONCLUSIONS

This study investigates the effects of neighboring WM microstructure on callosal interhemispheric communication in healthy controls and patients with mild traumatic brain injury, highlighting that neighboring noncallosal WM microstructures are involved in callosal interhemispheric communication and information transfer. Further longitudinal studies may provide insight into the temporal dynamics of interhemispheric recovery following mild traumatic brain injury.

Diffusion MRI of the Hippocampus

The hippocampus is a brain structure that plays key roles in a variety of cognitive processes. Critically, a wide range of neurological disorders are associated with degeneration of the hippocampal microstructure, defined as neurons, dendrites, glial cells, and more. Thus, the hippocampus is a key target for methods that are sensitive to these microscale properties. Diffusion MRI is one such method, which can noninvasively probe neural architecture. Here we review the extensive use of diffusion MRI to capture hippocampal microstructure in both health and disease. The results of these studies indicate that (1) diffusion tensor imaging is sensitive but not specific to the hippocampal microstructure; (2) biophysical modeling of diffusion MRI signals is a promising avenue to capture more specific aspects of the hippocampal microstructure; (3) use of ultra-short diffusion times have shown unique laminar-specific microstructure and response to hippocampal injury; (4) dispersion of microstructure is likely abundant in the hippocampus; and (5) the angular richness of the diffusion MRI signal can be leveraged to improve delineation of the internal hippocampal circuitry. Overall, extant findings suggest that diffusion MRI offers a promising avenue for characterizing hippocampal microstructure.

Cerebral microstructural alterations in Post-COVID-condition are related to cognitive impairment, olfactory dysfunction and fatigue

After contracting COVID-19, a substantial number of individuals develop a Post-COVID-Condition, marked by neurologic symptoms such as cognitive deficits, olfactory dysfunction, and fatigue. Despite this, biomarkers and pathophysiological understandings of this condition remain limited. Employing magnetic resonance imaging, we conduct a comparative analysis of cerebral microstructure among patients with Post-COVID-Condition, healthy controls, and individuals that contracted COVID-19 without long-term symptoms. We reveal widespread alterations in cerebral microstructure, attributed to a shift in volume from neuronal compartments to free fluid, associated with the severity of the initial infection. Correlating these alterations with cognition, olfaction, and fatigue unveils distinct affected networks, which are in close anatomical-functional relationship with the respective symptoms.

Deep learning with diffusion MRI as in vivo microscope reveals sex-related differences in human white matter microstructure

Biological sex is a crucial variable in neuroscience studies where sex differences have been documented across cognitive functions and neuropsychiatric disorders. While gross statistical differences have been previously documented in macroscopic brain structure such as cortical thickness or region size, less is understood about sex-related cellular-level microstructural differences which could provide insight into brain health and disease. Studying these microstructural differences between men and women paves the way for understanding brain disorders and diseases that manifest differently in different sexes. Diffusion MRI is an important in vivo, non-invasive methodology that provides a window into brain tissue microstructure. Our study develops multiple end-to-end classification models that accurately estimates the sex of a subject using volumetric diffusion MRI data and uses these models to identify white matter regions that differ the most between men and women. 471 male and 560 female healthy subjects (age range, 22-37 years) from the Human Connectome Project are included. Fractional anisotropy, mean diffusivity and mean kurtosis are used to capture brain tissue microstructure characteristics. Diffusion parametric maps are registered to a standard template to reduce bias that can arise from macroscopic anatomical differences like brain size and contour. This study employ three major model architectures: 2D convolutional neural networks, 3D convolutional neural networks and Vision Transformer (with self-supervised pretraining). Our results show that all 3 models achieve high sex classification performance (test AUC 0.92-0.98) across all diffusion metrics indicating definitive differences in white matter tissue microstructure between males and females. We further use complementary model architectures to inform about the pattern of detected microstructural differences and the influence of short-range versus long-range interactions. Occlusion analysis together with Wilcoxon signed-rank test is used to determine which white matter regions contribute most to sex classification. The results indicate that sex-related differences manifest in both local features as well as global features / longer-distance interactions of tissue microstructure. Our highly consistent findings across models provides new insight supporting differences between male and female brain cellular-level tissue organization particularly in the central white matter.

Resiliency to Alzheimer's disease neuropathology can be distinguished from dementia using cortical astrogliosis imaging

Despite the presence of significant Alzheimer's disease (AD) pathology, characterized by amyloid β (Aβ) plaques and phosphorylated tau (pTau) tangles, some cognitively normal elderly individuals do not inevitably develop dementia. These findings give rise to the notion of cognitive 'resilience', suggesting maintained cognitive function despite the presence of AD neuropathology, highlighting the influence of factors beyond classical pathology. Cortical astroglial inflammation, a ubiquitous feature of symptomatic AD, shows a strong correlation with cognitive impairment severity, potentially contributing to the diversity of clinical presentations. However, noninvasively imaging neuroinflammation, particularly astrogliosis, using MRI remains a significant challenge. Here we sought to address this challenge and to leverage multidimensional (MD) MRI, a powerful approach that combines relaxation with diffusion MR contrasts, to map cortical astrogliosis in the human brain by accessing sub-voxel information. Our goal was to test whether MD-MRI can map astroglial pathology in the cerebral cortex, and if so, whether it can distinguish cognitive resiliency from dementia in the presence of hallmark AD neuropathological changes. We adopted a multimodal approach by integrating histological and MRI analyses using human postmortem brain samples. Ex vivo cerebral cortical tissue specimens derived from three groups comprised of non-demented individuals with significant AD pathology postmortem, individuals with both AD pathology and dementia, and non-demented individuals with minimal AD pathology postmortem as controls, underwent MRI at 7 T. We acquired and processed MD-MRI, diffusion tensor, and quantitative T 1 and T 2 MRI data, followed by histopathological processing on slices from the same tissue. By carefully co-registering MRI and microscopy data, we performed quantitative multimodal analyses, leveraging targeted immunostaining to assess MD-MRI sensitivity and specificity towards Aβ, pTau, and glial fibrillary acidic protein (GFAP), a marker for astrogliosis. Our findings reveal a distinct MD-MRI signature of cortical astrogliosis, enabling the creation of predictive maps for cognitive resilience amid AD neuropathological changes. Multiple linear regression linked histological values to MRI changes, revealing that the MD-MRI cortical astrogliosis biomarker was significantly associated with GFAP burden (standardized β=0.658, pFDR<0.0001), but not with Aβ (standardized β=0.009, p FDR =0.913) or pTau (standardized β=-0.196, p FDR =0.051). Conversely, none of the conventional MRI parameters showed significant associations with GFAP burden in the cortex. While the extent to which pathological glial activation contributes to neuronal damage and cognitive impairment in AD is uncertain, developing a noninvasive imaging method to see its affects holds promise from a mechanistic perspective and as a potential predictor of cognitive outcomes.

Optimisation of quantitative brain diffusion-relaxation MRI acquisition protocols with physics-informed machine learning

Diffusion-relaxation MRI aims to extract quantitative measures that characterise microstructural tissue properties such as orientation, size, and shape, but long acquisition times are typically required. This work proposes a physics-informed learning framework to extract an optimal subset of diffusion-relaxation MRI measurements for enabling shorter acquisition times, predict non-measured signals, and estimate quantitative parameters. In vivo and synthetic brain 5D-Diffusion-T1-T2∗-weighted MRI data obtained from five healthy subjects were used for training and validation, and from a sixth participant for testing. One fully data-driven and two physics-informed machine learning methods were implemented and compared to two manual selection procedures and Cramér-Rao lower bound optimisation. The physics-informed approaches could identify measurement-subsets that yielded more consistently accurate parameter estimates in simulations than other approaches, with similar signal prediction error. Five-fold shorter protocols yielded error distributions of estimated quantitative parameters with very small effect sizes compared to estimates from the full protocol. Selected subsets commonly included a denser sampling of the shortest and longest inversion time, lowest echo time, and high b-value. The proposed framework combining machine learning and MRI physics offers a promising approach to develop shorter imaging protocols without compromising the quality of parameter estimates and signal predictions.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency,ω $$ \omega $$ , in addition to the diffusion tensor,D $$ \mathbf{D} $$ , and relaxation,R 1 $$ {R}_1 $$ ,R 2 $$ {R}_2 $$ , correlations. AD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on theirD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

A detailed spatio-temporal atlas of the white matter tracts for the fetal brain

This study presents the construction of a comprehensive spatiotemporal atlas detailing the development of white matter tracts in the fetal brain using diffusion magnetic resonance imaging (dMRI). Our research leverages data collected from fetal MRI scans conducted between 22 and 37 weeks of gestation, capturing the dynamic changes in the brain's microstructure during this critical period. The atlas includes 60 distinct white matter tracts, including commissural, projection, and association fibers. We employed advanced fetal dMRI processing techniques and tractography to map and characterize the developmental trajectories of these tracts. Our findings reveal that the development of these tracts is characterized by complex patterns of fractional anisotropy (FA) and mean diffusivity (MD), reflecting key neurodevelopmental processes such as axonal growth, involution of the radial-glial scaffolding, and synaptic pruning. This atlas can serve as a useful resource for neuroscience research and clinical practice, improving our understanding of the fetal brain and potentially aiding in the early diagnosis of neurodevelopmental disorders. By detailing the normal progression of white matter tract development, the atlas can be used as a benchmark for identifying deviations that may indicate neurological anomalies or predispositions to disorders.

The connection of motor improvement after deep brain stimulation in Parkinson's disease and microstructural integrity of the substantia nigra and subthalamic nucleus

BACKGROUND

Nigrostriatal microstructural integrity has been suggested as a biomarker for levodopa response in Parkinson's disease (PD), which is a strong predictor for motor response to deep brain stimulation (DBS) of the subthalamic nucleus (STN). This study aimed to explore the impact of microstructural integrity of the substantia nigra (SN), STN, and putamen on motor response to STN-DBS using diffusion microstructure imaging.

METHODS

Data was collected from 23 PD patients (mean age 63 ± 7, 6 females) who underwent STN-DBS, had preoperative 3 T diffusion magnetic resonance imaging including multishell diffusion-weighted MRI with b-values of 1000 and 2000 s/mm2 and records of motor improvement available.

RESULTS

The association between a poorer DBS-response and increased free interstitial fluid showed notable effect sizes (rho > |0.4|) in SN and STN, but not in putamen. However, this did not reach significance after Bonferroni correction and controlling for sex and age.

CONCLUSION

Microstructural integrity of SN and STN are potential biomarkers for the prediction of therapy efficacy following STN-DBS, but further studies are required to confirm these associations.

Developing an AI-empowered head-only ultra-high-performance gradient MRI system for high spatiotemporal neuroimaging

Recent advances in neuroscience requires high-resolution MRI to decipher the structural and functional details of the brain. Developing a high-performance gradient system is an ongoing effort in the field to facilitate high spatial and temporal encoding. Here, we proposed a head-only gradient system NeuroFrontier, dedicated for neuroimaging with an ultra-high gradient strength of 650 mT/m and 600 T/m/s. The proposed system features in 1) ultra-high power of 7MW achieved by running two gradient power amplifiers using a novel paralleling method; 2) a force/torque balanced gradient coil design with a two-step mechanical structure that allows high-efficiency and flexible optimization of the peripheral nerve stimulation; 3) a high-density integrated RF system that is miniaturized and customized for the head-only system; 4) an AI-empowered compressed sensing technique that enables ultra-fast acquisition of high-resolution images and AI-based acceleration in q-t space for diffusion MRI (dMRI); and 5) a prospective head motion correction technique that effectively corrects motion artifacts in real-time with 3D optical tracking. We demonstrated the potential advantages of the proposed system in imaging resolution, speed, and signal-to-noise ratio for 3D structural MRI (sMRI), functional MRI (fMRI) and dMRI in neuroscience applications of submillimeter layer-specific fMRI and dMRI. We also illustrated the unique strength of this system for dMRI-based microstructural mapping, e.g., enhanced lesion contrast at short diffusion-times or high b-values, and improved estimation accuracy for cellular microstructures using diffusion-time-dependent dMRI or for neurite microstructures using q-space approaches.

Bias-Reduced Neural Networks for Parameter Estimation in Quantitative MRI

Purpose

To develop neural network (NN)-based quantitative MRI parameter estimators with minimal bias and a variance close to the Cramér-Rao bound.

Theory and Methods

We generalize the mean squared error loss to control the bias and variance of the NN's estimates, which involves averaging over multiple noise realizations of the same measurements during training. Bias and variance properties of the resulting NNs are studied for two neuroimaging applications.

Results

In simulations, the proposed strategy reduces the estimates' bias throughout parameter space and achieves a variance close to the Cramér-Rao bound. In vivo, we observe good concordance between parameter maps estimated with the proposed NNs and traditional estimators, such as non-linear least-squares fitting, while state-of-the-art NNs show larger deviations.

Conclusion

The proposed NNs have greatly reduced bias compared to those trained using the mean squared error and offer significantly improved computational efficiency over traditional estimators with comparable or better accuracy.

Microstructural white matter alterations associated with social anxiety disorders: A systematic review

BACKGROUND

Social anxiety disorder (SAD) is a psychiatric condition characterized by impaired social functioning that negatively impacts individuals' quality of life. Previous neuroimaging studies have revealed morphological and functional changes in various brain regions associated with SAD. Recent advances in diffusion tensor imaging (DTI) and diffusion-weighted imaging (DWI) have enabled the investigation of microstructural white matter (WM) alterations in SAD. This study aims to provide an overview of DTI/DWI studies exploring WM microstructure changes in SAD.

METHODS

A systematic search on PubMed, Scopus, Web of Science, and PsycINFO was conducted for relevant records on July 8, 2023. An exploratory meta-analysis was also performed.

RESULTS

Eight studies were reviewed. Consistent findings indicated reduced fractional anisotropy and increased diffusivity measures in different WM tracts in SAD patients compared to healthy controls. These alterations were mostly observed within regions of the fronto-limbic network, like uncinate fasciculus (UF) and superior and inferior longitudinal fasciculi (SLF and ILF). Finally, our exploratory meta-analysis on four studies utilizing a voxel-wise analytic approach yielded no significant differences between SAD patients and controls.

LIMITATIONS

Limited number of studies, small sample sizes, and heterogeneity in analysis methods.

CONCLUSIONS

Patients with SAD exhibited altered WM integrity, particularly in the UF, SLF, and ILF, compared to healthy controls. However, due to the limited number of included studies, our meta-analysis yielded no significant differences between SAD patients and controls. Therefore, future research is crucial to unravel the link between altered WM structure and the involvement of other limbic and cortical structures in SAD pathogenesis.

Blurred streamlines: A novel representation to reduce redundancy in tractography

Tractography is a powerful tool to study brain connectivity in vivo, but it is well known to suffer from an intrinsic trade-off between sensitivity and specificity. A critical - but usually underrated - parameter to choose that can heavily impact the quality of the estimates is the number of streamlines to be reconstructed for a given data set. In fact, sensitivity can be improved by generating more and more streamlines, as all real anatomical connections are likely reconstructed, but lots of false positives are inevitably introduced, too. Consequently, so-called tractography filtering techniques have become increasingly popular to get rid of these false positives and improve specificity. However, increasing number of streamlines introduces redundancy in tractography reconstructions, which may negatively impact the performance of filtering algorithms, especially those based on linear formulations. To address this problem, we introduce a novel streamlines representation, called "blurred streamlines", which drastically reduces the redundancy among streamlines by (i) clustering similar trajectories and (ii) spatially blurring the corresponding signal contributions. We tested the effectiveness of the blurred streamlines both on synthetic and in vivo data. Our results clearly show that this new representation is as accurate as state-of-the-art methods despite using only 5% of the input streamlines, thus significantly decreasing the computational complexity of filtering algorithms as well as storage requirements of the resulting reconstructions.

The effects of axonal beading and undulation on axonal diameter estimation from diffusion MRI: Insights from simulations in human axons segmented from three-dimensional electron microscopy

The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings b , where the deviation from the expected 1 / b scaling in white matter yields a finite transverse diffusivity, which is then translated into an axon diameter estimate. While axons are usually modeled as perfectly straight, impermeable cylinders, local variations in diameter (caliber variation or beading) and direction (undulation) are known to influence axonal diameter estimates and have been observed in microscopy data of human axons. In this study, we performed Monte Carlo simulations of diffusion in axons reconstructed from three-dimensional electron microscopy of a human temporal lobe specimen using simulated sequence parameters matched to the maximal gradient strength of the next-generation Connectome 2.0 human MRI scanner ( ≲ 500 mT/m). We show that axon diameter estimation is accurate for nonbeaded, nonundulating fibers; however, in fibers with caliber variations and undulations, the axon diameter is heavily underestimated due to caliber variations, and this effect overshadows the known overestimation of the axon diameter due to undulations. This unexpected underestimation may originate from variations in the coarse-grained axial diffusivity due to caliber variations. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Facilitating diffusion tensor imaging of the brain during continuous gross head motion with first and second order motion compensating diffusion gradients

PURPOSE

To demonstrate the feasibility of motion compensating diffusion gradient schemes in the acquisition of quality diffusion tensor images (DTI) of the brain during continuous gross head motion.

METHODS

Five healthy subjects were scanned using a clinical 3 T MRI with and without continuous head motion. For one volunteer, DTI data was acquired using standard (M0) diffusion-weighted (DW) gradients, and first (M1) and second (M2) order gradient schemes that were previously developed for use in cardiac DTI. In four additional volunteers, DTI data was acquired with M0 and M2 gradients. DTI parameters were calculated and compared with established retrospective motion corrections.

RESULTS

In the absence of motion, DTI parameters calculated from M0, M1, and M2 data were consistent. In the presence of motion, up to 44% of DW images acquired with M0 gradients were corrupted by signal dropout, compared to 0% of the M2 images. In voxelwise comparisons, DTI parameters calculated using motion-M0 data were elevated compared to reference data. Retrospective corrections for extreme motion applied to motion-M0 data did not improve consistency with reference data in cases where motion corrupted >15% of DW images. In contrast, DTI parameters calculated with motion-M2 data were consistent with reference data.

CONCLUSION

This proof-of-principle study demonstrates that motion compensating diffusion gradients can mitigate artifacts because of continuous motion in DTI of the brain and offers promise for improved DTI accessibility. Further study will be necessary to determine the robustness of the approach in patient populations with high susceptibility to head motion.

Neural orientation distribution fields for estimation and uncertainty quantification in diffusion MRI

Inferring brain connectivity and structure in-vivo requires accurate estimation of the orientation distribution function (ODF), which encodes key local tissue properties. However, estimating the ODF from diffusion MRI (dMRI) signals is a challenging inverse problem due to obstacles such as significant noise, high-dimensional parameter spaces, and sparse angular measurements. In this paper, we address these challenges by proposing a novel deep-learning based methodology for continuous estimation and uncertainty quantification of the spatially varying ODF field. We use a neural field (NF) to parameterize a random series representation of the latent ODFs, implicitly modeling the often ignored but valuable spatial correlation structures in the data, and thereby improving efficiency in sparse and noisy regimes. An analytic approximation to the posterior predictive distribution is derived which can be used to quantify the uncertainty in the ODF estimate at any spatial location, avoiding the need for expensive resampling-based approaches that are typically employed for this purpose. We present empirical evaluations on both synthetic and real in-vivo diffusion data, demonstrating the advantages of our method over existing approaches.

Superficial white matter microstructural imaging method based on time-space fractional-order diffusion

Objective. Microstructure imaging based on diffusion magnetic resonance signal is an advanced imaging technique that enablesin vivomapping of the brain's microstructure. Superficial white matter (SWM) plays an important role in brain development, maturation, and aging, while fewer microstructure imaging methods address the SWM due to its complexity. Therefore, this study aims to develop a diffusion propagation model to investigate the microstructural characteristics of the SWM region.Approach. In this paper, we hypothesize that the effect of cell membrane permeability and the water exchange between soma and dendrites cannot be neglected for typical clinical diffusion times (20 ms Collapse

Key Words

• anomalous diffusion
• diffusion magnetic resonance imaging
• fractional diffusion equation
• microstructure imaging
• superficial white matter

Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Jianglin He School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, People's Republic of China
Yuanjun Wang School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, People's Republic of China

Collapse

Insights and improvements in correspondence between axonal volume fraction measured with diffusion-weighted MRI and electron microscopy

Biophysical diffusion-weighted imaging (DWI) models are increasingly used in neuroscience to estimate the axonal water fraction (f AW ), which in turn is key for noninvasive estimation of the axonal volume fraction (f A ). These models require thorough validation by comparison with a reference method, for example, electron microscopy (EM). While EM studies often neglect the unmyelinated axons and solely report the fraction of myelinated axons, in DWI both myelinated and unmyelinated axons contribute to the DWI signal. However, DWI models often include simplifications, for example, the neglect of differences in the compartmental relaxation times or fixed diffusivities, which in turn might affect the estimation off AW . We investigate whether linear calibration parameters (scaling and offset) can improve the comparability between EM- and DWI-based metrics off A . To this end, we (a) used six DWI models based on the so-called standard model of white matter (WM), including two models with fixed compartmental diffusivities (e.g., neurite orientation dispersion and density imaging, NODDI) and four models that fitted the compartmental diffusivities (e.g., white matter tract integrity, WMTI), and (b) used a multimodal data set including ex vivo diffusion DWI and EM data in mice with a broad dynamic range of fibre volume metrics. We demonstrated that the offset is associated with the volume fraction of unmyelinated axons and the scaling factor is associated with different compartmentalT 2 and can substantially enhance the comparability between EM- and DWI-based metrics off A . We found that DWI models that fitted compartmental diffusivities provided the most accurate estimates of the EM-basedf A . Finally, we introduced a more efficient hybrid calibration approach, where only the offset is estimated but the scaling is fixed to a theoretically predicted value. Using this approach, a similar one-to-one correspondence to EM was achieved for WMTI. The method presented can pave the way for use of validated DWI-based models in clinical research and neuroscience.

Restricted diffusion characteristics in oscillating gradient spin echo with mesoscopic phantom

In diffusion magnetic resonance imaging, oscillating gradient spin echo (OGSE) has an extremely short diffusion time if motion probing gradient (MPG) is applied to the waveform. Further, it can detect microstructural specificity. OGSE changes sensitivity to spin displacement velocity based on the MPG phase. The current study aimed to investigate the restricted diffusion characteristics of each OGSE waveform using the capillary phantom with various b-values, frequencies, and MPG phases. We performed OGSE (b-value = 300, 500, 800, 1200, 1600, and 2000 s/mm2) for the sine and cosine waveforms using the capillary phantom (6, 12, 25, 50, and 100 μm and free water) with a 9.4-T experimental magnetic resonance imaging system and a solenoid coil. We evaluated the axial and radial diffusivity (AD, RD) of each structure size. The output current of the MPG was assessed with an oscilloscope and analyzed with the gradient modulation power spectra by fast Fourier transform. In sine, the sidelobe spectrum was enhanced with increasing frequency, and the central spectrum slightly increased. The difference in RD was detected at 6 and 12 μm; however, it did not depend on the structure scale at 50 or 100 μm and free water. In cosine, the diffusion spectrum was enhanced, whereas the central spectrum decreased with increasing frequency. Both AD and RD in cosine had a frequency dependence, and AD and RD increased with a higher frequency regardless of structure size. AD and RD in either sine or cosine had no evident b-value dependence. We evaluated the OGSE-restricted diffusion characteristics. The measurements obtained diffusion information similar to the pulsed gradient spin echo. Hence, the cosine measurements indicated that a higher frequency could capture faster diffusion within the diffusion phenomena.

Assessment of Precision and Accuracy of Brain White Matter Microstructure using Combined Diffusion MRI and Relaxometry

Joint modeling of diffusion and relaxation has seen growing interest due to its potential to provide complementary information about tissue microstructure. For brain white matter, we designed an optimal diffusion-relaxometry MRI protocol that samples multiple b-values, B-tensor shapes, and echo times (TE). This variable-TE protocol (27 min) has as subsets a fixed-TE protocol (15 min) and a 2-shell dMRI protocol (7 min), both characterizing diffusion only. We assessed the sensitivity, specificity and reproducibility of these protocols with synthetic experiments and in six healthy volunteers. Compared with the fixed-TE protocol, the variable-TE protocol enables estimation of free water fractions while also capturing compartmental T 2 relaxation times. Jointly measuring diffusion and relaxation offers increased sensitivity and specificity to microstructure parameters in brain white matter with voxelwise coefficients of variation below 10%.

Incorporating the effect of white matter microstructure in the estimation of magnetic susceptibility in ex vivo mouse brain

PURPOSE

To extend quantitative susceptibility mapping to account for microstructure of white matter (WM) and demonstrate its effect on ex vivo mouse brain at 16.4T.

THEORY AND METHODS

Previous studies have shown that the MRI measured Larmor frequency also depends on local magnetic microstructure at the mesoscopic scale. Here, we include effects from WM microstructure using our previous results for the mesoscopic Larmor frequencyΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ of cylinders with arbitrary orientations. We scrutinize the validity of our model and QSM in a digital brain phantom includingΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ from a WM susceptibility tensor and biologically stored iron with scalar susceptibility. We also apply susceptibility tensor imaging to the phantom and investigate how the fitted tensors are biased fromΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ . Last, we demonstrate how to combine multi-gradient echo and diffusion MRI images of ex vivo mouse brains acquired at 16.4T to estimate an apparent scalar susceptibility without sample rotations.

RESULTS

Our new model improves susceptibility estimation compared to QSM for the brain phantom. Applying susceptibility tensor imaging to the phantom withΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ from WM axons with scalar susceptibility produces a highly anisotropic susceptibility tensor that mimics results from previous susceptibility tensor imaging studies. For the ex vivo mouse brain we find theΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ due to WM microstructure to be substantial, changing susceptibility in WM up to 25% root-mean-squared-difference.

CONCLUSION

Ω ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ impacts susceptibility estimates and biases susceptibility tensor imaging fitting substantially. Hence, it should not be neglected when imaging structurally anisotropic tissue such as brain WM.

Advanced diffusion imaging reveals microstructural characteristics of primary CNS lymphoma, allowing differentiation from glioblastoma

Background

Primary CNS lymphoma (PCNSL) and glioblastoma (GBM) both represent frequent intracranial malignancies with differing clinical management. However, distinguishing PCNSL from GBM with conventional MRI can be challenging when atypical imaging features are present. We employed advanced dMRI for noninvasive characterization of the microstructure of PCNSL and differentiation from GBM as the most frequent primary brain malignancy.

Methods

Multiple dMRI metrics including Diffusion Tensor Imaging, Neurite Orientation Dispersion and Density Imaging, and Diffusion Microstructure Imaging were extracted from the contrast-enhancing tumor component in 10 PCNSL and 10 age-matched GBM on 3T MRI. Imaging findings were correlated with cell density and axonal markers obtained from histopathology.

Results

We found significantly increased intra-axonal volume fractions (V-intra and intracellular volume fraction) and microFA in PCNSL compared to GBM (all P < .001). In contrast, mean diffusivity (MD), axial diffusivity (aD), and microADC (all P < .001), and also free water fractions (V-CSF and V-ISO) were significantly lower in PCNSL (all P < .01). Receiver-operating characteristic analysis revealed high predictive values regarding the presence of a PCNSL for MD, aD, microADC, V-intra, ICVF, microFA, V-CSF, and V-ISO (area under the curve [AUC] in all >0.840, highest for MD and ICVF with an AUC of 0.960). Comparative histopathology between PCNSL and GBM revealed a significantly increased cell density in PCNSL and the presence of axonal remnants in a higher proportion of samples.

Conclusions

Advanced diffusion imaging enables the characterization of the microstructure of PCNSL and reliably distinguishes PCNSL from GBM. Both imaging and histopathology revealed a relatively increased cell density and a preserved axonal microstructure in PCNSL.

Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain

The dependence of the diffusion MRI signal on the diffusion time carries signatures of restricted diffusion and exchange. Here we seek to highlight these signatures in the human brain by performing experiments using free gradient waveforms designed to be selectively sensitive to the two effects. We examine six healthy volunteers using both strong and ultra-strong gradients (80, 200 and 300 mT/m). In an experiment featuring a large set of 150 gradient waveforms with different sensitivities to restricted diffusion and exchange, our results reveal unique and different time-dependence signatures in grey and white matter. Grey matter was characterised by both restricted diffusion and exchange and white matter predominantly by restricted diffusion. Exchange in grey matter was at least twice as fast as in white matter, across all subjects and all gradient strengths. The cerebellar cortex featured relatively short exchange times (115 ms). Furthermore, we show that gradient waveforms with tailored designs can be used to map exchange in the human brain. We also assessed the feasibility of clinical applications of the method used in this work and found that the exchange-related contrast obtained with a 25-minute protocol at 300 mT/m was preserved in a 4-minute protocol at 300 mT/m and a 10-minute protocol at 80 mT/m. Our work underlines the utility of free waveforms for detecting time dependence signatures due to restricted diffusion and exchange in vivo, which may potentially serve as a tool for studying diseased tissue.

A simulation-driven supervised learning framework to estimate brain microstructure using diffusion MRI

We propose a framework to train supervised learning models on synthetic data to estimate brain microstructure parameters using diffusion magnetic resonance imaging (dMRI). Although further validation is necessary, the proposed framework aims to seamlessly incorporate realistic simulations into dMRI microstructure estimation. Synthetic data were generated from over 1,000 neuron meshes converted from digital neuronal reconstructions and linked to their neuroanatomical parameters (such as soma volume and neurite length) using an optimized diffusion MRI simulator that produces intracellular dMRI signals from the solution of the Bloch-Torrey partial differential equation. By combining random subsets of simulated neuron signals with a free diffusion compartment signal, we constructed a synthetic dataset containing dMRI signals and 40 tissue microstructure parameters of 1.45 million artificial brain voxels. To implement supervised learning models we chose multilayer perceptrons (MLPs) and trained them on a subset of the synthetic dataset to estimate some microstructure parameters, namely, the volume fractions of soma, neurites, and the free diffusion compartment, as well as the area fractions of soma and neurites. The trained MLPs perform satisfactorily on the synthetic test sets and give promising in-vivo parameter maps on the MGH Connectome Diffusion Microstructure Dataset (CDMD). Most importantly, the estimated volume fractions showed low dependence on the diffusion time, the diffusion time independence of the estimated parameters being a desired property of quantitative microstructure imaging. The synthetic dataset we generated will be valuable for the validation of models that map between the dMRI signals and microstructure parameters. The surface meshes and microstructures parameters of the aforementioned neurons have been made publicly available.

Reconstructing the somatotopic organization of the corticospinal tract remains a challenge for modern tractography methods

The corticospinal tract (CST) is a critically important white matter fiber tract in the human brain that enables control of voluntary movements of the body. The CST exhibits a somatotopic organization, which means that the motor neurons that control specific body parts are arranged in order within the CST. Diffusion magnetic resonance imaging (MRI) tractography is increasingly used to study the anatomy of the CST. However, despite many advances in tractography algorithms over the past decade, modern, state-of-the-art methods still face challenges. In this study, we compare the performance of six widely used tractography methods for reconstructing the CST and its somatotopic organization. These methods include constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD-Stream) methods, unscented Kalman filter (UKF) tractography methods including multi-fiber (UKF2T) and single-fiber (UKF1T) models, the generalized q-sampling imaging (GQI) based deterministic tractography method, and the TractSeg method. We investigate CST somatotopy by dividing the CST into four subdivisions per hemisphere that originate in the leg, trunk, hand, and face areas of the primary motor cortex. A quantitative and visual comparison is performed using diffusion MRI data (N = 100 subjects) from the Human Connectome Project. Quantitative evaluations include the reconstruction rate of the eight anatomical subdivisions, the percentage of streamlines in each subdivision, and the coverage of the white matter-gray matter (WM-GM) interface. CST somatotopy is further evaluated by comparing the percentage of streamlines in each subdivision to the cortical volumes for the leg, trunk, hand, and face areas. Overall, UKF2T has the highest reconstruction rate and cortical coverage. It is the only method with a significant positive correlation between the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex. However, our experimental results show that all compared tractography methods are biased toward generating many trunk streamlines (ranging from 35.10% to 71.66% of total streamlines across methods). Furthermore, the coverage of the WM-GM interface in the largest motor area (face) is generally low (under 40%) for all compared tractography methods. Different tractography methods give conflicting results regarding the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex, indicating that there is generally no clear relationship, and that reconstruction of CST somatotopy is still a large challenge. Overall, we conclude that while current tractography methods have made progress toward the well-known challenge of improving the reconstruction of the lateral projections of the CST, the overall problem of performing a comprehensive CST reconstruction, including clinically important projections in the lateral (hand and face areas) and medial portions (leg area), remains an important challenge for diffusion MRI tractography.

Outliers in diffusion-weighted MRI: Exploring detection models and mitigation strategies

Diffusion-weighted MRI (dMRI) is a medical imaging method that can be used to investigate the brain microstructure and structural connections between different brain regions. The method, however, requires relatively complex data processing frameworks and analysis pipelines. Many of these approaches are vulnerable to signal dropout artefacts that can originate from subjects moving their head during the scan. To combat these artefacts and eliminate such outliers, researchers have proposed two approaches: to replace outliers or to downweight outliers during modelling and analysis. With the rising interest in dMRI for clinical research, these types of corrections are increasingly important. Therefore, we set out to investigate the differences between outlier replacement and weighting approaches to help the dMRI community to select the best tool for their data processing pipelines. We evaluated dMRI motion correction registration and single tensor model fit pipelines using Gaussian Process and Spherical Harmonic based replacement approaches and outlier downweighting using highly realistic whole-brain simulations. As a proof of concept, we applied these approaches to dMRI infant data sets that contained varying numbers of dropout artefacts. Based on our results, we concluded that the Gaussian Process based outlier replacement provided similar tensor fit results to Gaussian Process based outlier detection and downweighting. Therefore, if only the least-squares estimate of the single tensor model is of interest, our recommendation is to use outlier replacement. However, outlier downweighting can potentially provide a more accurate estimate of the model precision which could be relevant for applications such as probabilistic tractoraphy.

Advanced Diffusion-Weighted MRI for Cancer Microstructure Assessment in Body Imaging, and Its Relationship With Histology

Diffusion-weighted magnetic resonance imaging (DW-MRI) aims to disentangle multiple biological signal sources in each imaging voxel, enabling the computation of innovative maps of tissue microstructure. DW-MRI model development has been dominated by brain applications. More recently, advanced methods with high fidelity to histology are gaining momentum in other contexts, for example, in oncological applications of body imaging, where new biomarkers are urgently needed. The objective of this article is to review the state-of-the-art of DW-MRI in body imaging (ie, not including the nervous system) in oncology, and to analyze its value as compared to reference colocalized histology measurements, given that demonstrating the histological validity of any new DW-MRI method is essential. In this article, we review the current landscape of DW-MRI techniques that extend standard apparent diffusion coefficient (ADC), describing their acquisition protocols, signal models, fitting settings, microstructural parameters, and relationship with histology. Preclinical, clinical, and in/ex vivo studies were included. The most used techniques were intravoxel incoherent motion (IVIM; 36.3% of used techniques), diffusion kurtosis imaging (DKI; 16.7%), vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT; 13.3%), and imaging microstructural parameters using limited spectrally edited diffusion (IMPULSED; 11.7%). Another notable category of techniques relates to innovative b-tensor diffusion encoding or joint diffusion-relaxometry. The reviewed approaches provide histologically meaningful indices of cancer microstructure (eg, vascularization/cellularity) which, while not necessarily accurate numerically, may still provide useful sensitivity to microscopic pathological processes. Future work of the community should focus on improving the inter-/intra-scanner robustness, and on assessing histological validity in broader contexts. LEVEL OF EVIDENCE: NA TECHNICAL EFFICACY: Stage 2.

Probing brain tissue microstructure with MRI: principles, challenges, and the role of multidimensional diffusion-relaxation encoding

Diffusion MRI uses the random displacement of water molecules to sensitize the signal to brain microstructure and to properties such as the density and shape of cells. Microstructure modeling techniques aim to estimate these properties from acquired data by separating the signal between virtual tissue 'compartments' such as the intra-neurite and the extra-cellular space. A key challenge is that the diffusion MRI signal is relatively featureless compared with the complexity of brain tissue. Another challenge is that the tissue microstructure is wildly different within the gray and white matter of the brain. In this review, we use results from multidimensional diffusion encoding techniques to discuss these challenges and their tentative solutions. Multidimensional encoding increases the information content of the data by varying not only the b-value and the encoding direction but also additional experimental parameters such as the shape of the b-tensor and the echo time. Three main insights have emerged from such encoding. First, multidimensional data contradict common model assumptions on diffusion and T2 relaxation, and illustrates how the use of these assumptions cause erroneous interpretations in both healthy brain and pathology. Second, many model assumptions can be dispensed with if data are acquired with multidimensional encoding. The necessary data can be easily acquired in vivo using protocols optimized to minimize Cramér-Rao lower bounds. Third, microscopic diffusion anisotropy reflects the presence of axons but not dendrites. This insight stands in contrast to current 'neurite models' of brain tissue, which assume that axons in white matter and dendrites in gray matter feature highly similar diffusion. Nevertheless, as an axon-based contrast, microscopic anisotropy can differentiate gray and white matter when myelin alterations confound conventional MRI contrasts.

Subject classification and cross-time prediction based on functional connectivity and white matter microstructure features in a rat model of Alzheimer's using machine learning

BACKGROUND

The pathological process of Alzheimer's disease (AD) typically takes decades from onset to clinical symptoms. Early brain changes in AD include MRI-measurable features such as altered functional connectivity (FC) and white matter degeneration. The ability of these features to discriminate between subjects without a diagnosis, or their prognostic value, is however not established.

METHODS

The main trigger mechanism of AD is still debated, although impaired brain glucose metabolism is taking an increasingly central role. Here, we used a rat model of sporadic AD, based on impaired brain glucose metabolism induced by an intracerebroventricular injection of streptozotocin (STZ). We characterized alterations in FC and white matter microstructure longitudinally using functional and diffusion MRI. Those MRI-derived measures were used to classify STZ from control rats using machine learning, and the importance of each individual measure was quantified using explainable artificial intelligence methods.

RESULTS

Overall, combining all the FC and white matter metrics in an ensemble way was the best strategy to discriminate STZ rats, with a consistent accuracy over 0.85. However, the best accuracy early on was achieved using white matter microstructure features, and later on using FC. This suggests that consistent damage in white matter in the STZ group might precede FC. For cross-timepoint prediction, microstructure features also had the highest performance while, in contrast, that of FC was reduced by its dynamic pattern which shifted from early hyperconnectivity to late hypoconnectivity.

CONCLUSIONS

Our study highlights the MRI-derived measures that best discriminate STZ vs control rats early in the course of the disease, with potential translation to humans.

Implications of fitting a two-compartment model in single-shell diffusion MRI

It is becoming increasingly common for studies to fit single-shell diffusion MRI data to a two-compartment model, which comprises a hindered cellular compartment and a freely diffusing isotropic compartment. These studies consistently find that the fraction of the isotropic compartment (f) is sensitive to white matter (WM) conditions and pathologies, although the actual biological source of changes infhas not been validated. In this work we put aside the biological interpretation offand study the sensitivity implications of fitting single-shell data to a two-compartment model. We identify a nonlinear transformation between the one-compartment model (diffusion tensor imaging, DTI) and a two-compartment model in which the mean diffusivities of both compartments are effectively fixed. While the analytic relationship implies that fitting this two-compartment model does not offer any more information than DTI, it explains why metrics derived from a two-compartment model can exhibit enhanced sensitivity over DTI to certain types of WM processes, such as age-related WM differences. The sensitivity enhancement should not be viewed as a substitute for acquiring multi-shell data. Rather, the results of this study provide insight into the consequences of choosing a two-compartment model when only single-shell data is available.

Multishell diffusion MRI reveals whole-brain white matter changes in HIV

Diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) have been previously used to explore white matter related to human immunodeficiency virus (HIV) infection. While DTI and DKI suffer from low specificity, the Combined Hindered and Restricted Model of Diffusion (CHARMED) provides additional microstructural specificity. We used these three models to evaluate microstructural differences between 35 HIV-positive patients without neurological impairment and 20 healthy controls who underwent diffusion-weighted imaging using three b-values. While significant group effects were found in all diffusion metrics, CHARMED and DKI analyses uncovered wider involvement (80% vs. 20%) of all white matter tracts in HIV infection compared with DTI. In restricted fraction (FR) analysis, we found significant differences in the left corticospinal tract, middle cerebellar peduncle, right inferior cerebellar peduncle, right corticospinal tract, splenium of the corpus callosum, left superior cerebellar peduncle, left superior cerebellar peduncle, pontine crossing tract, left posterior limb of the internal capsule, and left/right medial lemniscus. These are involved in language, motor, equilibrium, behavior, and proprioception, supporting the functional integration that is frequently impaired in HIV-positivity. Additionally, we employed a machine learning algorithm (XGBoost) to discriminate HIV-positive patients from healthy controls using DTI and CHARMED metrics on an ROIwise basis, and unique contributions to this discrimination were examined using Shapley Explanation values. The CHARMED and DKI estimates produced the best performance. Our results suggest that biophysical multishell imaging, combining additional sensitivity and built-in specificity, provides further information about the brain microstructural changes in multimodal areas involved in attentive, emotional and memory networks often impaired in HIV patients.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω , in addition to the diffusion tensor, D , and relaxation, R 1 , R 2 , correlations. A D ( ω ) - R 1 - R 2 clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D ( ω ) - R 1 - R 2 distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Mapping the individual human cortex using multidimensional MRI and unsupervised learning

Human evolution has seen the development of higher-order cognitive and social capabilities in conjunction with the unique laminar cytoarchitecture of the human cortex. Moreover, early-life cortical maldevelopment has been associated with various neurodevelopmental diseases. Despite these connections, there is currently no noninvasive technique available for imaging the detailed cortical laminar structure. This study aims to address this scientific and clinical gap by introducing an approach for imaging human cortical lamina. This method combines diffusion-relaxation multidimensional MRI with a tailored unsupervised machine learning approach that introduces enhanced microstructural sensitivity. This new imaging method simultaneously encodes the microstructure, the local chemical composition and importantly their correlation within complex and heterogenous tissue. To validate our approach, we compared the intra-cortical layers obtained using our ex vivo MRI-based method with those derived from Nissl staining of postmortem human brain specimens. The integration of unsupervised learning with diffusion-relaxation correlation MRI generated maps that demonstrate sensitivity to areal differences in cytoarchitectonic features observed in histology. Significantly, our observations revealed layer-specific diffusion-relaxation signatures, showing reductions in both relaxation times and diffusivities at the deeper cortical levels. These findings suggest a radial decrease in myelin content and changes in cell size and anisotropy, reflecting variations in both cytoarchitecture and myeloarchitecture. Additionally, we demonstrated that 1D relaxation and high-order diffusion MRI scalar indices, even when aggregated and used jointly in a multimodal fashion, cannot disentangle the cortical layers. Looking ahead, our technique holds the potential to open new avenues of research in human neurodevelopment and the vast array of disorders caused by disruptions in neurodevelopment.

Computational mesoscale framework for biological clustering and fractal aggregation

Hierarchical clustering due to diffusion and reaction is a widespread occurrence in natural phenomena, displaying fractal behavior with non-integer size scaling. The study of this phenomenon has garnered interest in both biological systems such as morphogenesis and blood clotting, and synthetic systems such as colloids and polymers. The modeling of biological clustering can be difficult, as it can occur on a variety of scales and involve multiple mechanisms, necessitating the use of various methods to capture its behavior. Here, we propose a novel framework, the generalized-mesoscale-clustering (GMC), for the study of complex hierarchical clustering phenomena in biological systems. The GMC framework incorporates the effects of hydrodynamic interactions, bonding, and surface tension, and allows for the analysis of both static and dynamic states of cluster development. The framework is applied to a range of biological clustering mechanisms, with a focus on blood-related clustering from fibrin network formation to platelet aggregation. Our study highlights the importance of a comprehensive characterization of the structural properties of the cluster, including fractal dimension, pore-scale diffusion, initiation time, and consolidation time, in fully understanding the behavior of biological clustering systems. The GMC framework also provides the potential to investigate the temporal evolution and mechanical properties of the clusters by tracking bond density and including hydrodynamic interactions.

Fractal relaxation model with a nonlinear diffusion coefficient for fitting anomalous diffusion data in magnetic resonance imaging

In this paper a relaxation model of signal attenuation in diffusion-weighted magnetic resonance imaging (dMRI) with varying diffusion coefficient in terms of fractal derivative is proposed, in which the diffusion coefficient is a power law of the effective diffusion time. The relaxation model provides measures of diffusion constant, fractal dimension of diffusive trajectory of water molecule and the time power-law behavior of the diffusion coefficient. The proposed model was used to describe the magnetic resonance attenuation signal of the bullfrog sciatic nerve, and the corresponding spectral entropy was calculated to detect the environmental complexity in bullfrog sciatic nerve for water molecular diffusion. The results showed that the fractal derivative relaxation model (the VDC model) can accurately depict the diffusion pattern of water molecules in complex heterogeneous biological media at large b values. The VDC model provides an alternative theoretical reference for biological tissue detection based on time-dependent diffusion of water molecules.

Anomalous diffusion, nonergodicity, non-Gaussianity, and aging of fractional Brownian motion with nonlinear clocks

How do nonlinear clocks in time and/or space affect the fundamental properties of a stochastic process? Specifically, how precisely may ergodic processes such as fractional Brownian motion (FBM) acquire predictable nonergodic and aging features being subjected to such conditions? We address these questions in the current study. To describe different types of non-Brownian motion of particles-including power-law anomalous, ultraslow or logarithmic, as well as superfast or exponential diffusion-we here develop and analyze a generalized stochastic process of scaled-fractional Brownian motion (SFBM). The time- and space-SFBM processes are, respectively, constructed based on FBM running with nonlinear time and space clocks. The fundamental statistical characteristics such as non-Gaussianity of particle displacements, nonergodicity, as well as aging are quantified for time- and space-SFBM by selecting different clocks. The latter parametrize power-law anomalous, ultraslow, and superfast diffusion. The results of our computer simulations are fully consistent with the analytical predictions for several functional forms of clocks. We thoroughly examine the behaviors of the probability-density function, the mean-squared displacement, the time-averaged mean-squared displacement, as well as the aging factor. Our results are applicable for rationalizing the impact of nonlinear time and space properties superimposed onto the FBM-type dynamics. SFBM offers a general framework for a universal and more precise model-based description of anomalous, nonergodic, non-Gaussian, and aging diffusion in single-molecule-tracking observations.

The impact of head orientation with respect to B0 on diffusion tensor MRI measures

Diffusion tensor MRI (DT-MRI) remains the most commonly used approach to characterise white matter (WM) anisotropy. However, DT estimates may be affected by tissue orientation w.r.t. B → 0 due to local gradients and intrinsic T 2 orientation dependence induced by the microstructure. This work aimed to investigate whether and how diffusion tensor MRI-derived measures depend on the orientation of the head with respect to the static magnetic field, B → 0 . By simulating WM as two compartments, we demonstrated that compartmental T 2 anisotropy can induce the dependence of diffusion tensor measures on the angle between WM fibres and the magnetic field. In in vivo experiments, reduced radial diffusivity and increased axial diffusivity were observed in white matter fibres perpendicular to B → 0 compared to those parallel to B → 0 . Fractional anisotropy varied by up to 20 % as a function of the angle between WM fibres and the orientation of the main magnetic field. To conclude, fibre orientation w.r.t. B → 0 is responsible for up to 7 % variance in diffusion tensor measures across the whole brain white matter from all subjects and head orientations. Fibre orientation w.r.t. B → 0 may introduce additional variance in clinical research studies using diffusion tensor imaging, particularly when it is difficult to control for (e.g., fetal or neonatal imaging, or when the trajectories of fibres change due to, e.g., space occupying lesions).

Microstructural Alterations in Tract Development in College Football and Volleyball Players: A Longitudinal Diffusion MRI Study

BACKGROUND AND OBJECTIVES

Repeated impacts in high-contact sports such as American football can affect the brain's microstructure, which can be studied using diffusion MRI. Most imaging studies are cross-sectional, do not include low-contact players as controls, or lack advanced tract-specific microstructural metrics. We aimed to investigate longitudinal changes in high-contact collegiate athletes compared with low-contact controls using advanced diffusion MRI and automated fiber quantification.

METHODS

We examined brain microstructure in high-contact (football) and low-contact (volleyball) collegiate athletes with up to 4 years of follow-up. Inclusion criteria included university and team enrollment. Exclusion criteria included history of neurosurgery, severe brain injury, and major neurologic or substance abuse disorder. We investigated diffusion metrics along the length of tracts using nested linear mixed-effects models to ascertain the acute and chronic effects of subconcussive and concussive impacts, and associations between diffusion changes with clinical, behavioral, and sports-related measures.

RESULTS

Forty-nine football and 24 volleyball players (271 total scans) were included. Football players had significantly divergent trajectories in multiple microstructural metrics and tracts. Longitudinal increases in fractional anisotropy and axonal water fraction, and decreases in radial/mean diffusivity and orientation dispersion index, were present in volleyball but absent in football players (all findings |T-statistic|> 3.5, p value <0.0001). This pattern was present in the callosum forceps minor, superior longitudinal fasciculus, thalamic radiation, and cingulum hippocampus. Longitudinal differences were more prominent and observed in more tracts in concussed football players (n = 24, |T|> 3.6, p < 0.0001). An analysis of immediate postconcussion scans (n = 12) demonstrated a transient localized increase in axial diffusivity and mean/radial kurtosis in the uncinate and cingulum hippocampus (|T| > 3.7, p < 0.0001). Finally, within football players, those with high position-based impact risk demonstrated increased intracellular volume fraction longitudinally (T = 3.6, p < 0.0001).

DISCUSSION

The observed longitudinal changes seen in football, and especially concussed athletes, could reveal diminished myelination, altered axonal calibers, or depressed pruning processes leading to a static, nondecreasing axonal dispersion. This prospective longitudinal study demonstrates divergent tract-specific trajectories of brain microstructure, possibly reflecting a concussive and repeated subconcussive impact-related alteration of white matter development in football athletes.

Imaging local diffusion in microstructures using NV-based pulsed field gradient NMR

Understanding diffusion in microstructures plays a crucial role in many scientific fields, including neuroscience, medicine, or energy research. While magnetic resonance (MR) methods are the gold standard for diffusion measurements, spatial encoding in MR imaging has limitations. Here, we introduce nitrogen-vacancy (NV) center-based nuclear MR (NMR) spectroscopy as a powerful tool to probe diffusion within microscopic sample volumes. We have developed an experimental scheme that combines pulsed gradient spin echo (PGSE) with optically detected NV-NMR spectroscopy, allowing local quantification of molecular diffusion and flow. We demonstrate correlated optical imaging with spatially resolved PGSE NV-NMR experiments probing anisotropic water diffusion within an individual model microstructure. Our optically detected PGSE NV-NMR technique opens up prospects for extending the current capabilities of investigating diffusion processes with the future potential of probing single cells, tissue microstructures, or ion mobility in thin film materials for battery applications.

In vivo Mapping of Cellular Resolution Neuropathology in Brain Ischemia by Diffusion MRI

Non-invasive mapping of cellular pathology can provide critical diagnostic and prognostic information. Recent developments in diffusion MRI have produced new tools for examining tissue microstructure at a level well below the imaging resolution. Here, we report the use of diffusion time ( t )-dependent diffusion kurtosis imaging ( t DKI) to simultaneously assess the morphology and transmembrane permeability of cells and their processes in the context of pathological changes in hypoxic-ischemic brain (HI) injury. Through Monte Carlo simulations and cell culture organoid imaging, we demonstrate feasibility in measuring effective size and permeability changes based on the peak and tail of t DKI curves. In a mouse model of HI, in vivo imaging at 11.7T detects a marked shift of the t DKI peak to longer t in brain edema, suggesting swelling and beading associated with the astrocytic processes and neuronal neurites. Furthermore, we observed a faster decrease of the t DKI tail in injured brain regions, reflecting increased membrane permeability that was associated with upregulated water exchange upon astrocyte activation at acute stage as well as necrosis with disrupted membrane integrity at subacute stage. Such information, unavailable with conventional diffusion MRI at a single t, can predict salvageable tissues. For a proof-of-concept, t DKI at 3T on an ischemic stroke patient suggested increased membrane permeability in the stroke region. This work therefore demonstrates the potential of t DKI for in vivo detection of the pathological changes in microstructural morphology and transmembrane permeability after ischemic injury using a clinically translatable protocol.

Comparison of different neurite density metrics with brain asymmetry evaluation

The standard diffusion MRI model with intra- and extra-axonal water pools offers a set of microstructural parameters describing brain white matter architecture. However, non-linearities in the standard model and diffusion data contamination by noise and imaging artefacts make estimation of diffusion metrics challenging. In order to develop reliable diffusion approaches and to avoid computational model degeneracy, additional theoretical assumptions allowing stable numerical implementations are required. Advanced diffusion approaches allow for estimation of intra-axonal water fraction (AWF), describing a key structural characteristic of brain tissue. AWF can be interpreted as an indirect measure or proxy of neurite density and has a potential as useful clinical biomarker. Established diffusion approaches such as white matter tract integrity, neurite orientation dispersion and density imaging (NODDI), and spherical mean technique provide estimates of AWF within their respective theoretical frameworks. In the present study, we estimated AWF metrics using different diffusion approaches and compared measures of brain asymmetry between the different metrics in a sub-sample of 182 subjects from the UK Biobank. Multivariate decomposition by mean of linked independent component analysis revealed that the various AWF proxies derived from the different diffusion approaches reflect partly non-overlapping variance of independent components, with distinct anatomical distributions and sensitivity to age. Further, voxel-wise analysis revealed age-related differences in AWF-based brain asymmetry, indicating less apparent left-right hemisphere difference with higher age. Finally, we demonstrated that NODDI metrics suffer from a quite strong dependence on used numerical algorithms and post-processing pipeline. The analysis based on AWF metrics strongly depends on the used diffusion approach and leads to poorly reproducible results.