Orientationally invariant indices of axon diameter and density from diffusion MRI

Microstructural Damage and Repair in the Spinal Cord of Patients With Early Multiple Sclerosis and Association With Disability at 5 Years

BACKGROUND AND OBJECTIVES

The dynamics of microstructural spinal cord (SC) damage and repair in people with multiple sclerosis (pwMS) and their clinical relevance have yet to be explored. We set out to describe patient-specific profiles of microstructural SC damage and change during the first year after MS diagnosis and to investigate their associations with disability and SC atrophy at 5 years.

METHODS

We performed a longitudinal monocentric cohort study among patients with relapsing-remitting MS: first relapse <1 year, no relapse <1 month, and high initial severity on MRI (>9 T2 lesions on brain MRI and/or initial myelitis). pwMS and age-matched healthy controls (HCs) underwent cervical SC magnetization transfer (MT) imaging at baseline and at 1 year for pwMS. Based on HC data, SC MT ratio z-score maps were computed for each person with MS. An index of microstructural damage was calculated as the proportion of voxels classified as normal at baseline and identified as damaged after 1 year. Similarly, an index of repair was also calculated (voxels classified as damaged at baseline and as normal after 1 year). Linear models including these indices and disability or SC cross-sectional area (CSA) change between baseline and 5 years were implemented.

RESULTS

Thirty-seven patients and 19 HCs were included. We observed considerable variability in the extent of microstructural SC damage at baseline (0%-58% of SC voxels). We also observed considerable variability in damage and repair indices over 1 year (0%-31% and 0%-20%), with 18 patients showing predominance of damage and 18 predominance of repair. The index of microstructural damage was associated positively with the Expanded Disability Status Scale score (r = 0.504, p = 0.002) and negatively with CSA change (r = -0.416, p = 0.02) at 5 years, independent of baseline SC lesion volume.

DISCUSSION

People with early relapsing-remitting MS exhibited heterogeneous profiles of microstructural SC damage and repair. Progression of microstructural damage was associated with disability progression and SC atrophy 5 years later. These results indicate a potential for microstructural repair in the SC to prevent disability progression in pwMS.

Early-life gut inflammation drives sex-dependent shifts in the microbiome-endocrine-brain axis

Despite recent advances in understanding the connection between the gut microbiota and the adult brain, significant knowledge gaps remain regarding how gut inflammation affects brain development. We hypothesized that gut inflammation during early life would negatively affect neurodevelopment by disrupting microbiota communication to the brain. We therefore developed a novel pediatric chemical model of inflammatory bowel disease (IBD), an incurable condition affecting millions of people worldwide. IBD is characterized by chronic intestinal inflammation, and is associated with comorbid symptoms such as anxiety, depression and cognitive impairment. Notably, 25% of patients with IBD are diagnosed during childhood, and the effects of chronic inflammation during this critical developmental period remain poorly understood. This study investigated the effects of early-life gut inflammation induced by DSS (dextran sulfate sodium) on a range of microbiota, endocrine, and behavioral outcomes, focusing on sex-specific impacts. DSS-treated mice exhibited increased intestinal inflammation and altered microbiota membership, which correlated with changes in microbiota-derived circulating metabolites. The majority of behavioral measures were unaffected, with the exception of impaired mate-seeking behaviors in DSS-treated males. DSS-treated males also showed significantly smaller seminal vesicles, lower circulating androgens, and decreased intestinal hormone-activating enzyme activity compared to vehicle controls. In the brain, DSS treatment led to chronic, sex-specific alterations in microglial morphology. These results suggest that early-life gut inflammation causes changes in gut microbiota composition, affecting short-chain fatty acid (SCFA) producers and glucuronidase (GUS) activity, correlating with altered SCFA and androgen levels. The findings highlight the developmental sensitivity to inflammation-induced changes in endocrine signalling and emphasize the long-lasting physiological and microbiome changes associated with juvenile IBD.

A simulated annealing algorithm for randomizing weighted networks

Scientific discovery in connectomics relies on network null models. The prominence of network features is conventionally evaluated against null distributions estimated using randomized networks. Modern imaging technologies provide an increasingly rich array of biologically meaningful edge weights. Despite the prevalence of weighted graph analysis in connectomics, randomization models that only preserve binary node degree remain most widely used. Here we propose a simulated annealing procedure for generating randomized networks that preserve weighted degree (strength) sequences. We show that the procedure outperforms other rewiring algorithms and generalizes to multiple network formats, including directed and signed networks, as well as diverse real-world networks. Throughout, we use morphospace representation to assess the sampling behavior of the algorithm and the variability of the resulting ensemble. Finally, we show that accurate strength preservation yields different inferences about brain network organization. Collectively, this work provides a simple but powerful method to analyze richly detailed next-generation connectomics datasets.

Susceptibility-induced internal gradients reveal axon morphology and cause anisotropic effects in the diffusion-weighted MRI signal

Diffusion-weighted MRI is our most promising method for estimating microscopic tissue morphology in vivo. The signal acquisition is based on scanner-generated external magnetic gradients. However, it will also be affected by susceptibility-induced internal magnetic gradients caused by interactions between the tissue and the static magnetic field of the scanner. With 3D in silico experiments, we show how internal gradients cause morphology-, compartment-, and orientation-dependence of spin-echo and pulsed-gradient spin-echo experiments in myelinated axons. These effects surpass those observed with previous 2D modelling corresponding to straight cylinders. For an ex vivo monkey brain, we observe the orientation-dependence generated only when including non-circular cross-sections in the in silico morphological configurations, and find orientation-dependent deviation of up to 17% for diffusion tensor metrics. Interestingly, we find that the orientation-dependence not only biases the signal across different brain regions, but also carries a sensitivity to the morphology of axonal cross-sections which is not attainable by the idealised theoretical diffusion-weighted MRI signal.

Romer-EPTI: Rotating-view motion-robust super-resolution EPTI for SNR-efficient distortion-free in-vivo mesoscale diffusion MRI and microstructure imaging

PURPOSE

To overcome the major challenges in diffusion MRI (dMRI) acquisition, including limited SNR, distortion/blurring, and susceptibility to motion artifacts.

THEORY AND METHODS

A novel Romer-EPTI technique is developed to achieve SNR-efficient acquisition while providing distortion-free imaging, minimal spatial blurring, high motion robustness, and simultaneous multi-TE imaging. It introduces a ROtating-view Motion-robust supEr-Resolution technique (Romer) combined with a distortion/blurring-free Echo Planar Time-resolved Imaging (EPTI) readout. Romer enhances SNR through simultaneous multi-thick-slice acquisition with rotating-view encoding, while providing high motion-robustness via a high-fidelity, motion-aware super-resolution reconstruction. Instead of EPI, the in-plane encoding is performed using EPTI readout to prevent geometric distortion, T2/T2*-blurring, and importantly, dynamic distortions that could introduce additional blurring/artifacts after super-resolution reconstruction due to combining volumes with inconsistent geometries. This further improves effective spatial resolution and motion robustness. Additional developments include strategies to address slab-boundary artifacts, achieve minimized TE and optimized readout for additional SNR gain, and increase robustness to strong phase variations at high b-values.

RESULTS

Using Romer-EPTI, we demonstrated distortion-free whole-brain mesoscale in-vivo dMRI at both 3T (500-μm isotropic [iso] resolution) and 7T (485-μm iso resolution) for the first time. Motion experiments demonstrated the technique's motion robustness and its ability to obtain high-resolution diffusion images in the presence of subject motion. Romer-EPTI also demonstrated high SNR gain and robustness in high b-value (b = 5000 s/mm2) and time-dependent dMRI.

CONCLUSION

The high SNR efficiency, improved image quality, and motion robustness of Romer-EPTI make it a highly efficient acquisition for high-resolution dMRI and microstructure imaging.

Synaptic Density in Early Stages of Psychosis and Clinical High Risk

Importance

Synaptic dysfunction is involved in schizophrenia pathophysiology. However, whether in vivo synaptic density is reduced in early stages of psychosis, including its high-risk states, remains unclear.

Objective

To investigate whether synaptic density (synaptic vesicle glycoprotein 2A [SV2A] binding potential) is reduced in first-episode psychosis (FEP) and in clinical high risk (CHR) and investigate the effect of cannabis use on synaptic density and examine its relationship with psychotic symptoms and gray matter microstructure across groups.

Design, Setting, and Participants

This cross-sectional study was performed in a tertiary care psychiatric hospital from July 2021 to October 2023. Participants were patients with antipsychotic-free or minimally exposed FEP or CHR and healthy controls with a clean urine drug screen (except cannabis).

Main Outcomes and Measures

Synaptic density was quantified with dynamic 90-minute [18F]SynVesT-1 positron emission tomography (PET) scans across prioritized brain regions of interest (ROIs) delineated in individual magnetic resonance images (MRIs). Cannabis use was confirmed with urine drug screens. Gray matter microstructure was assessed using diffusion-weighted MRI to estimate neurite density.

Results

A total of 49 participants were included, including 16 patients with FEP (mean [SD] age, 26.1 [4.6] years; 9 males and 7 females), 17 patients at CHR (mean [SD] age, 21.2 [3.5] years; 8 males and 9 females), and 16 healthy controls (mean [SD] age, 23.4 [3.6] years; 7 males and 9 females). Synaptic density was significantly different between groups (F2,273 = 4.02, P = .02, Cohen F = 0.17; ROI: F5,273 = 360.18, P < .01, Cohen F = 2.55) with a group × ROI interaction (F10,273 = 2.67, P < .01, Cohen F = 0.32). Synaptic density was lower in cannabis users (F1,272 = 5.31, P = .02, Cohen F = 0.14). Lower synaptic density across groups was associated with more negative symptoms (Positive and Negative Syndrome Scale negative scores: F1,81 = 4.31, P = .04, Cohen F = 0.23; Scale of Psychosis-Risk Symptoms negative scores: F1,90 = 4.12, P = .04, Cohen F = 0.21). SV2A binding potential was significantly associated with neurite density index (F1,138 = 6.76, P = .01, Cohen F = 0.22).

Conclusions and Relevance

This study found that synaptic density reductions were present during the early stages of psychosis and its risk states and associated with negative symptoms. The implications of SV2A for negative symptoms in psychosis and CHR warrant further investigation. Future studies should investigate the impact of cannabis use on synaptic density in CHR longitudinally.

Improving MR axon radius estimation in human white matter using spiral acquisition and field monitoring

PURPOSE

To compare MR axon radius estimation in human white matter using a multiband spiral sequence combined with field monitoring to the current state-of-the-art echo-planar imaging (EPI)-based approach.

METHODS

A custom multiband spiral sequence was used for diffusion-weighted imaging at ultra-highb $$ b $$ -values. Field monitoring and higher order image reconstruction were employed to greatly reduce artifacts in spiral images. Diffusion weighting parameters were chosen to match a state-of-the art EPI-based axon radius mapping protocol. The spiral approach was compared to the EPI approach by comparing the image signal-to-noise ratio (SNR) and performing a test-retest study to assess the respective variability and repeatability of axon radius mapping. Effective axon radius estimates were compared over white matter voxels and along the left corticospinal tract.

RESULTS

Increased SNR and reduced artifacts in spiral images led to reduced variability in resulting axon radius maps, especially in low-SNR regions. Test-retest variability was reduced by a factor of approximately 1.5 using the spiral approach. Reduced repeatability due to significant bias was found for some subjects in both spiral and EPI approaches, and attributed to scanner instability, pointing to a previously unknown limitation of the state-of-the-art approach.

CONCLUSION

Combining spiral readouts with field monitoring improved mapping of the effective axon radius compared to the conventional EPI approach.

Age-related alterations in human cortical microstructure across the lifespan: Insights from high-gradient diffusion MRI

The human brain undergoes age-related microstructural alterations across the lifespan. Soma and Neurite Density Imaging (SANDI), a novel biophysical model of diffusion MRI, provides estimates of cell body (soma) radius and density, and neurite density in gray matter. The goal of this cross-sectional study was to assess the sensitivity of high-gradient diffusion MRI toward age-related alterations in cortical microstructure across the adult lifespan using SANDI. Seventy-two cognitively unimpaired healthy subjects (ages 19-85 years; 40 females) were scanned on the 3T Connectome MRI scanner with a maximum gradient strength of 300mT/m using a multi-shell diffusion MRI protocol incorporating 8 b-values and diffusion time of 19 ms. Intra-soma signal fraction obtained from SANDI model-fitting to the data was strongly correlated with age in all major cortical lobes (r = -0.69 to -0.60, FDR-p < 0.001). Intra-soma signal fraction (r = 0.48-0.63, FDR-p < 0.001) and soma radius (r = 0.28-0.40, FDR-p < 0.04) were significantly correlated with cortical volume in the prefrontal cortex, frontal, parietal, and temporal lobes. The strength of the relationship between SANDI metrics and age was greater than or comparable to the relationship between cortical volume and age across the cortical regions, particularly in the occipital lobe and anterior cingulate gyrus. In contrast to the SANDI metrics, all associations between diffusion tensor imaging (DTI) and diffusion kurtosis imaging metrics and age were low to moderate. These results suggest that high-gradient diffusion MRI may be more sensitive to underlying substrates of neurodegeneration in the aging brain than DTI and traditional macroscopic measures of neurodegeneration such as cortical volume and thickness.

Comparisons of MR and EM inferred tissue microstructure properties using a human autopsy corpus callosum sample

Degeneration of white matter (WM) microstructure in the central nervous system is characteristic of many neurodegenerative conditions. Previous research indicates that axonal degeneration visible in ex vivo electron microscopy (EM) photomicrographs precede the onset of clinical symptoms. Measuring WM microstructural features, such as axon diameter and packing fraction, currently require these highly invasive methods of analysis and it is therefore of great importance to develop methods for in vivo measurements. Diffusion weighted Magnetic Resonance Imaging (MRI) is a non-invasive method which can be used in conjunction with temporal diffusion spectroscopy (TDS) and an oscillating gradient spin echo (OGSE) pulse sequence to probe micron-scale structures within neural tissue. The current experiment aims to compare axon diameter measurements, mean effective axon diameter (AxD¯), and packing fractions calculated from EM histopathological analysis and inferred values from MR images. Mathematical models of axon diameters used for analysis include the ActiveAx Frequency-Dependent Extra-Axonal Diffusion (AAD) model and the AxCaliber Frequency-Dependent Extra-Axonal Diffusion (ACD) model using ROI (Region of Interest) based analysis (RBA) and voxel-based analysis (VBA), respectively. Overall, it was observed that MRI inferred WM microstructural parameters overestimate those calculated from EM. This may be attributable to tissue shrinkage during EM dehydration, the sensitivity of MR pulse sequences to larger diameter axons, and/or inaccurate model assumptions. The results of the current study provide a means to characterize the precision and accuracy of RBA-ACD and VBA-AAD OGSE-TDS and highlight the need for further research investigating the relationship between ex vivo MRI and EM, with the goal of reaching in vivo MRI.

Restriction-induced time-dependent transcytolemmal water exchange: Revisiting the Kӓrger exchange model

The Kӓrger model and its derivatives have been widely used to incorporate transcytolemmal water exchange rate, an essential characteristic of living cells, into analyses of diffusion MRI (dMRI) signals from tissues. The Kӓrger model consists of two homogeneous exchanging components coupled by an exchange rate constant and assumes measurements are made with sufficiently long diffusion time and slow water exchange. Despite successful applications, it remains unclear whether these assumptions are generally valid for practical dMRI sequences and biological tissues. In particular, barrier-induced restrictions to diffusion produce inhomogeneous magnetization distributions in relatively large-sized compartments such as cancer cells, violating the above assumptions. The effects of this inhomogeneity are usually overlooked. We performed computer simulations to quantify how restriction effects, which in images produce edge enhancements at compartment boundaries, influence different variants of the Kӓrger-model. The results show that the edge enhancement effect will produce larger, time-dependent estimates of exchange rates in e.g., tumors with relatively large cell sizes (>10 μm), resulting in overestimations of water exchange as previously reported. Moreover, stronger diffusion gradients, longer diffusion gradient durations, and larger cell sizes, all cause more pronounced edge enhancement effects. This helps us to better understand the feasibility of the Kärger model in estimating water exchange in different tissue types and provides useful guidance on signal acquisition methods that may mitigate the edge enhancement effect. This work also indicates the need to correct the overestimated transcytolemmal water exchange rates obtained assuming the Kärger-model.

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.

DIMOND: DIffusion Model OptimizatioN with Deep Learning

Diffusion magnetic resonance imaging is an important tool for mapping tissue microstructure and structural connectivity non-invasively in the in vivo human brain. Numerous diffusion signal models are proposed to quantify microstructural properties. Nonetheless, accurate estimation of model parameters is computationally expensive and impeded by image noise. Supervised deep learning-based estimation approaches exhibit efficiency and superior performance but require additional training data and may be not generalizable. A new DIffusion Model OptimizatioN framework using physics-informed and self-supervised Deep learning entitled "DIMOND" is proposed to address this problem. DIMOND employs a neural network to map input image data to model parameters and optimizes the network by minimizing the difference between the input acquired data and synthetic data generated via the diffusion model parametrized by network outputs. DIMOND produces accurate diffusion tensor imaging results and is generalizable across subjects and datasets. Moreover, DIMOND outperforms conventional methods for fitting sophisticated microstructural models including the kurtosis and NODDI model. Importantly, DIMOND reduces NODDI model fitting time from hours to minutes, or seconds by leveraging transfer learning. In summary, the self-supervised manner, high efficacy, and efficiency of DIMOND increase the practical feasibility and adoption of microstructure and connectivity mapping in clinical and neuroscientific applications.

Romer-EPTI: rotating-view motion-robust super-resolution EPTI for SNR-efficient distortion-free in-vivo mesoscale dMRI and microstructure imaging

Purpose

To overcome the major challenges in dMRI acquisition, including low SNR, distortion/blurring, and motion vulnerability.

Methods

A novel Romer-EPTI technique is developed to provide distortion-free dMRI with significant SNR gain, high motion-robustness, sharp spatial resolution, and simultaneous multi-TE imaging. It introduces a ROtating-view Motion-robust supEr-Resolution technique (Romer) combined with a distortion/blurring-free EPTI encoding. Romer enhances SNR by a simultaneous multi-thick-slice acquisition with rotating-view encoding, while providing high motion-robustness through a motion-aware super-resolution reconstruction, which also incorporates slice-profile and real-value diffusion, to resolve high-isotropic-resolution volumes. The in-plane encoding is performed using distortion/blurring-free EPTI, which further improves effective spatial resolution and motion robustness by preventing not only T2/T2*-blurring but also additional blurring resulting from combining encoded volumes with inconsistent geometries caused by dynamic distortions. Self-navigation was incorporated to enable efficient phase correction. Additional developments include strategies to address slab-boundary artifacts, achieve minimal TE for SNR gain at 7T, and achieve high robustness to strong phase variations at high b-values.

Results

Using Romer-EPTI, we demonstrate distortion-free whole-brain mesoscale in-vivo dMRI at both 3T (500-μm-iso) and 7T (485-μm-iso) for the first time, with high SNR efficiency (e.g., 25 × ), and high image quality free from distortion and slab-boundary artifacts with minimal blurring. Motion experiments demonstrate Romer-EPTI's high motion-robustness and ability to recover sharp images in the presence of motion. Romer-EPTI also demonstrates significant SNR gain and robustness in high b-value (b=5000s/mm2) and time-dependent dMRI.

Conclusion

Romer-EPTI significantly improves SNR, motion-robustness, and image quality, providing a highly efficient acquisition for high-resolution dMRI and microstructure imaging.

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.

Reduced White Matter Fiber Density in Patients with Multiple Sclerosis

Introduction: Improved understanding of multiple sclerosis (MS) symptomatology, disease mechanisms, and clinical effectiveness can be achieved by investigating microstructural damage. The aim was to gain deeper insights into changes in white matter (WM) tracts in MS patients. Methods: Diffusion magnetic resonance imaging-based tractography was utilized to segment WM tracts into regions of interest for further quantitative analysis. However, tractography is susceptible to false-positive findings, reducing its specificity and clinical feasibility. To address these limitations, the Convex Optimization Modeling for Microstructure Informed Tractography (COMMIT) technique was used. COMMIT was used to derive measures of intracellular compartment (IC) and isotropic compartments from multishell diffusion data of 40 healthy controls (HCs) and 40 MS patients. Results: The analysis revealed a widespread pattern of significantly decreased IC values in MS patients compared with HCs across 61,581 voxels (pFWE < 0.05, threshold-free cluster enhancement [TFCE] corrected). Similar WM structures studied using the fractional anisotropy (FA) value also showed a reduction in FA among MS patients compared with HCs across 57,304 voxels (pFWE < 0.05, TFCE corrected). Out of the 61,581 voxels exhibiting lower IC, a substantial overlap of 47,251 voxels (76.72%) also demonstrated lower FA in MS patients compared with HCs. Discussion: The data suggested that lower IC values contributed to the explanation of FA reductions. In addition, IC showed promising potential for evaluating microstructural abnormalities in WM in MS, potentially being more sensitive than the frequently used FA value.

Quantitative diffusion imaging and genotype-by-sex interactions in a rat model of Alexander disease

PURPOSE

The clinical diagnosis and classification of Alexander disease (AxD) relies in part on qualitative neuroimaging biomarkers; however, these biomarkers fail to distinguish and discriminate different subtypes of AxD, especially in the presence of overlap in clinical symptoms. To address this gap in knowledge, we applied neurite orientation dispersion and density imaging (NODDI) to an innovative CRISPR-Cas9 rat genetic model of AxD to gain quantitative insights into the neural substrates and brain microstructural changes seen in AxD and to potentially identify novel quantitative NODDI biomarkers of AxD.

METHODS

Multi-shell DWI of age- and sex-matched AxD and wild-type Sprague Dawley rats (n = 6 per sex per genotype) was performed and DTI and NODDI measures calculated. A 3 × 2 × 2 analysis of variance model was used to determine the effect of genotype, biological sex, and laterality on quantitative measures of DTI and NODDI across regions of interest implicated in AxD.

RESULTS

There is a significant effect of genotype in the amygdala, hippocampus, neocortex, and thalamus in measures of both DTI and NODDI brain microstructure. A genotype by biological sex interaction was identified in DTI and NODDI measures in the corpus callosum, hippocampus, and neocortex.

CONCLUSION

We present the first application of NODDI to the study of AxD using a rat genetic model of AxD. Our analysis identifies alterations in NODDI and DTI measures to large white matter tracts and subcortical gray nuclei. We further identified genotype by sex interactions, suggesting a possible role for biological sex in the neuropathogenesis of AxD.

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.

A translational MRI approach to validate acute axonal damage detection as an early event in multiple sclerosis

Axonal degeneration is a central pathological feature of multiple sclerosis and is closely associated with irreversible clinical disability. Current noninvasive methods to detect axonal damage in vivo are limited in their specificity and clinical applicability, and by the lack of proper validation. We aimed to validate an MRI framework based on multicompartment modeling of the diffusion signal (AxCaliber) in rats in the presence of axonal pathology, achieved through injection of a neurotoxin damaging the neuronal terminal of axons. We then applied the same MRI protocol to map axonal integrity in the brain of multiple sclerosis relapsing-remitting patients and age-matched healthy controls. AxCaliber is sensitive to acute axonal damage in rats, as demonstrated by a significant increase in the mean axonal caliber along the targeted tract, which correlated with neurofilament staining. Electron microscopy confirmed that increased mean axonal diameter is associated with acute axonal pathology. In humans with multiple sclerosis, we uncovered a diffuse increase in mean axonal caliber in most areas of the normal-appearing white matter, preferentially affecting patients with short disease duration. Our results demonstrate that MRI-based axonal diameter mapping is a sensitive and specific imaging biomarker that links noninvasive imaging contrasts with the underlying biological substrate, uncovering generalized axonal damage in multiple sclerosis as an early event.

Fast and robust quantification of uncertainty in non-linear diffusion MRI models

Diffusion MRI (dMRI) allows for non-invasive investigation of brain tissue microstructure. By fitting a model to the dMRI signal, various quantitative measures can be derived from the data, such as fractional anisotropy, neurite density and axonal radii maps. We investigate the Fisher Information Matrix (FIM) and uncertainty propagation as a generally applicable method for quantifying the parameter uncertainties in linear and non-linear diffusion MRI models. In direct comparison with Markov Chain Monte Carlo (MCMC) sampling, the FIM produces similar uncertainty estimates at much lower computational cost. Using acquired and simulated data, we then list several characteristics that influence the parameter variances, including data complexity and signal-to-noise ratio. For practical purposes we investigate a possible use of uncertainty estimates in decreasing intra-group variance in group statistics by uncertainty-weighted group estimates. This has potential use cases for detection and suppression of imaging artifacts.

The human brain connectome weighted by the myelin content and total intra-axonal cross-sectional area of white matter tracts

A central goal in neuroscience is the development of a comprehensive mapping between structural and functional brain features, which facilitates mechanistic interpretation of brain function. However, the interpretability of structure-function brain models remains limited by a lack of biological detail. Here, we characterize human structural brain networks weighted by multiple white matter microstructural features including total intra-axonal cross-sectional area and myelin content. We report edge-weight-dependent spatial distributions, variance, small-worldness, rich club, hubs, as well as relationships with function, edge length, and myelin. Contrasting networks weighted by the total intra-axonal cross-sectional area and myelin content of white matter tracts, we find opposite relationships with functional connectivity, an edge-length-independent inverse relationship with each other, and the lack of a canonical rich club in myelin-weighted networks. When controlling for edge length, networks weighted by either fractional anisotropy, radial diffusivity, or neurite density show no relationship with whole-brain functional connectivity. We conclude that the co-utilization of structural networks weighted by total intra-axonal cross-sectional area and myelin content could improve our understanding of the mechanisms mediating the structure-function brain relationship.

A cross-comparative analysis of in vivo versus ex vivo MRI indices in a mouse model of concussion

BACKGROUND

We present a cross-sectional, case-matched, and pair-wise comparison of structural magnetic resonance imaging (MRI), diffusion tensor imaging (DTI), and neurite orientation dispersion and density imaging (NODDI) measures in vivo and ex vivo in a mouse model of concussion, thus aiming to establish the concordance of structural and diffusion imaging findings in living brain and after fixation.

METHODS

We allocated 28 male mice aged 3-4 months to sham injury and concussion (CON) groups. CON mice had received a single concussive impact on day 0 and underwent MRI at day 2 (n = 9) or 7 (n = 10) post-impact, and sham control mice likewise underwent imaging at day 2 (n = 5) or 7 (n = 4). Immediately after the final scanning, we collected the perfusion-fixed brains, which were stored for imaging ex vivo 6-12 months later. We then compared the structural imaging, DTI, and NODDI results between different methods.

RESULTS

In vivo to ex vivo structural and DTI/NODDI findings were in notably poor agreement regarding the effects of concussion on structural integrity of the brain.

COMPARISON WITH EXISTING METHODS

ex vivo imaging was frequently done to study the effects of diseases and treatments, but our results showed that ex vivo and in vivo imaging can detect completely opposite and contradictory results. This is also the first study that compares in vivo and ex vivo NODDI.

CONCLUSION

Our findings call for caution in extrapolating translational capabilities obtained ex vivo to physiological measurements in vivo. The divergent findings may reflect fixation artefacts and the contribution of the glymphatic system changes.

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

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.

Study of structural network connectivity using DTI tractography in insomnia disorder

Most of tractography studies on insomnia disorder (ID) have reported decreased structural connectivity between cortical and subcortical structures. Tractography based on standard diffusion tensor imaging (DTI) can generate high number of false-positive streamlines connections between gray matter regions. In the present study, we employed the convex optimization modeling for microstructure informed tractography-2 (COMMIT2) to improve the accuracy of the reconstructed whole-brain connectome and filter implausible brain connections in 28 patients with ID and compared with 27 healthy controls. Then, we used NBS-predict (a prediction-based extension to the network-based statistic method) in the COMMIT2-weighted connectome. Our results revealed decreased structural connectivity between subregions of the left somatomotor, ventral attention, frontoparietal, dorsal attention and default mode networks in the insomnia group. Moreover, there is a negative correlation between sleep efficiency and structural connectivity within the left frontoparietal, visual, default mode network, limbic, dorsal attention, right dorsal attention as well as right default mode networks. By comparing with standard connectivity analysis, we showed that by removing of false-positive streamlines connections after COMMIT2 filtering, abnormal structural connectivity was reduced in patients with ID compared to controls. Our results demonstrate the importance of improving the accuracy of tractography for understanding structural connectivity networks in ID.

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.

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.

Microstructure Fingerprinting for Heterogeneously Oriented Tissue Microenvironments

Most diffusion biophysical models capture basic properties of tissue microstructure, such as diffusivity and anisotropy. More realistic models that relate the diffusion-weighted signal to cell size and membrane permeability often require simplifying assumptions such as short gradient pulse and Gaussian phase distribution, leading to tissue features that are not necessarily quantitative. Here, we propose a method to quantify tissue microstructure without jeopardizing accuracy owing to unrealistic assumptions. Our method utilizes realistic signals simulated from the geometries of cellular microenvironments as fingerprints, which are then employed in a spherical mean estimation framework to disentangle the effects of orientation dispersion from microscopic tissue properties. We demonstrate the efficacy of microstructure fingerprinting in estimating intra-cellular, extra-cellular, and intra-soma volume fractions as well as axon radius, soma radius, and membrane permeability.

Single-subject electroencephalography measurement of interhemispheric transfer time for the in-vivo estimation of axonal morphology

Assessing axonal morphology in vivo opens new avenues for the combined study of brain structure and function. A novel approach has recently been introduced to estimate the morphology of axonal fibers from the combination of magnetic resonance imaging (MRI) data and electroencephalography (EEG) measures of the interhemispheric transfer time (IHTT). In the original study, the IHTT measures were computed from EEG data averaged across a group, leading to bias of the axonal morphology estimates. Here, we seek to estimate axonal morphology from individual measures of IHTT, obtained from EEG data acquired in a visual evoked potential experiment. Subject-specific IHTTs are computed in a data-driven framework with minimal a priori constraints, based on the maximal peak of neural responses to visual stimuli within periods of statistically significant evoked activity in the inverse solution space. The subject-specific IHTT estimates ranged from 8 to 29 ms except for one participant and the between-session variability was comparable to between-subject variability. The mean radius of the axonal radius distribution, computed from the IHTT estimates and the MRI data, ranged from 0 to 1.09 μm across subjects. The change in axonal g-ratio with axonal radius ranged from 0.62 to 0.81 μm-α . The single-subject measurement of the IHTT yields estimates of axonal morphology that are consistent with histological values. However, improvement of the repeatability of the IHTT estimates is required to improve the specificity of the single-subject axonal morphology estimates.

Tractography passes the test: Results from the diffusion-simulated connectivity (disco) challenge

Estimating structural connectivity from diffusion-weighted magnetic resonance imaging is a challenging task, partly due to the presence of false-positive connections and the misestimation of connection weights. Building on previous efforts, the MICCAI-CDMRI Diffusion-Simulated Connectivity (DiSCo) challenge was carried out to evaluate state-of-the-art connectivity methods using novel large-scale numerical phantoms. The diffusion signal for the phantoms was obtained from Monte Carlo simulations. The results of the challenge suggest that methods selected by the 14 teams participating in the challenge can provide high correlations between estimated and ground-truth connectivity weights, in complex numerical environments. Additionally, the methods used by the participating teams were able to accurately identify the binary connectivity of the numerical dataset. However, specific false positive and false negative connections were consistently estimated across all methods. Although the challenge dataset doesn't capture the complexity of a real brain, it provided unique data with known macrostructure and microstructure ground-truth properties to facilitate the development of connectivity estimation methods.

Reduced neurite density index in the prefrontal cortex of adults with autism assessed using neurite orientation dispersion and density imaging

Background

Core symptoms of autism-spectrum disorder (ASD) have been associated with prefrontal cortex abnormalities. However, the mechanisms behind the observation remain incomplete, partially due to the challenges of modeling complex gray matter (GM) structures. This study aimed to identify GM microstructural alterations in adults with ASD using neurite orientation dispersion and density imaging (NODDI) and voxel-wise GM-based spatial statistics (GBSS) to reduce the partial volume effects from the white matter and cerebrospinal fluid.

Materials and methods

A total of 48 right-handed participants were included, of which 22 had ASD (17 men; mean age, 34.42 ± 8.27 years) and 26 were typically developing (TD) individuals (14 men; mean age, 32.57 ± 9.62 years). The metrics of NODDI (neurite density index [NDI], orientation dispersion index [ODI], and isotropic volume fraction [ISOVF]) were compared between groups using GBSS. Diffusion tensor imaging (DTI) metrics and surface-based cortical thickness were also compared. The associations between magnetic resonance imaging-based measures and ASD-related scores, including ASD-spectrum quotient, empathizing quotient, and systemizing quotient were also assessed in the region of interest (ROI) analysis.

Results

After controlling for age, sex, and intracranial volume, GBSS demonstrated significantly lower NDI in the ASD group than in the TD group in the left prefrontal cortex (caudal middle frontal, lateral orbitofrontal, pars orbitalis, pars triangularis, rostral middle frontal, and superior frontal region). In the ROI analysis of individuals with ASD, a significantly positive correlation was observed between the NDI in the left rostral middle frontal, superior frontal, and left frontal pole and empathizing quotient score. No significant between-group differences were observed in all DTI metrics, other NODDI (i.e., ODI and ISOVF) metrics, and cortical thickness.

Conclusion

GBSS analysis was used to demonstrate the ability of NODDI metrics to detect GM microstructural alterations in adults with ASD, while no changes were detected using DTI and cortical thickness evaluation. Specifically, we observed a reduced neurite density index in the left prefrontal cortices associated with reduced empathic abilities.

Estimating axon radius using diffusion-relaxation MRI: calibrating a surface-based relaxation model with histology

Axon radius is a potential biomarker for brain diseases and a crucial tissue microstructure parameter that determines the speed of action potentials. Diffusion MRI (dMRI) allows non-invasive estimation of axon radius, but accurately estimating the radius of axons in the human brain is challenging. Most axons in the brain have a radius below one micrometer, which falls below the sensitivity limit of dMRI signals even when using the most advanced human MRI scanners. Therefore, new MRI methods that are sensitive to small axon radii are needed. In this proof-of-concept investigation, we examine whether a surface-based axonal relaxation process could mediate a relationship between intra-axonal T2 and T1 times and inner axon radius, as measured using postmortem histology. A unique in vivo human diffusion-T1-T2 relaxation dataset was acquired on a 3T MRI scanner with ultra-strong diffusion gradients, using a strong diffusion-weighting (i.e., b = 6,000 s/mm2) and multiple inversion and echo times. A second reduced diffusion-T2 dataset was collected at various echo times to evaluate the model further. The intra-axonal relaxation times were estimated by fitting a diffusion-relaxation model to the orientation-averaged spherical mean signals. Our analysis revealed that the proposed surface-based relaxation model effectively explains the relationship between the estimated relaxation times and the histological axon radius measured in various corpus callosum regions. Using these histological values, we developed a novel calibration approach to predict axon radius in other areas of the corpus callosum. Notably, the predicted radii and those determined from histological measurements were in close agreement.

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.

Axon fiber orientation as the source of T1 relaxation anisotropy in white matter: A study on corpus callosum in vivo and ex vivo

PURPOSE

Recent studies indicate that T1 in white matter (WM) is influenced by fiber orientation in B0 . The purpose of the study was to investigate the interrelationships between axon fiber orientation in corpus callosum (CC) and T1 relaxation time in humans in vivo as well as in rat brain ex vivo.

METHODS

Volunteers were scanned for relaxometric and diffusion MRI at 3 T and 7 T. Angular T1 plots from WM were computed using fractional anisotropy and fiber-to-field-angle maps. T1 and fiber-to-field angle were measured in five sections of CC to estimate the effects of inherently varying fiber orientations on T1 within the same tracts in vivo. Ex vivo rat-brain preparation encompassing posterior CC was rotated in B0 and T1 , and diffusion MRI images acquired at 9.4 T. T1 angular plots were determined at several rotation angles in B0 .

RESULTS

Angular T1 plots from global WM provided reference for estimated fiber orientation-linked T1 changes within CC. In anterior midbody of CC in vivo, where small axons are dominantly present, a shift in axon orientation is accompanied by a change in T1 , matching that estimated from WM T1 data. In CC, where large and giant axons are numerous, the measured T1 change is about 2-fold greater than the estimated one. Ex vivo rotation of the same midsagittal CC region of interest produced angular T1 plots at 9.4 T, matching those observed at 7 T in vivo.

CONCLUSION

These data causally link axon fiber orientation in B0 to the T1 relaxation anisotropy in WM.

CACTUS: a computational framework for generating realistic white matter microstructure substrates

Monte-Carlo diffusion simulations are a powerful tool for validating tissue microstructure models by generating synthetic diffusion-weighted magnetic resonance images (DW-MRI) in controlled environments. This is fundamental for understanding the link between micrometre-scale tissue properties and DW-MRI signals measured at the millimetre-scale, optimizing acquisition protocols to target microstructure properties of interest, and exploring the robustness and accuracy of estimation methods. However, accurate simulations require substrates that reflect the main microstructural features of the studied tissue. To address this challenge, we introduce a novel computational workflow, CACTUS (Computational Axonal Configurator for Tailored and Ultradense Substrates), for generating synthetic white matter substrates. Our approach allows constructing substrates with higher packing density than existing methods, up to 95% intra-axonal volume fraction, and larger voxel sizes of up to 500μm3 with rich fibre complexity. CACTUS generates bundles with angular dispersion, bundle crossings, and variations along the fibres of their inner and outer radii and g-ratio. We achieve this by introducing a novel global cost function and a fibre radial growth approach that allows substrates to match predefined targeted characteristics and mirror those reported in histological studies. CACTUS improves the development of complex synthetic substrates, paving the way for future applications in microstructure imaging.

Novel insights into axon diameter and myelin content in late childhood and adolescence

White matter microstructural development in late childhood and adolescence is driven predominantly by increasing axon density and myelin thickness. Ex vivo studies suggest that the increase in axon diameter drives developmental increases in axon density observed with pubertal onset. In this cross-sectional study, 50 typically developing participants aged 8-18 years were scanned using an ultra-strong gradient magnetic resonance imaging scanner. Microstructural properties, including apparent axon diameter $({d}_a)$, myelin content, and g-ratio, were estimated in regions of the corpus callosum. We observed age-related differences in ${d}_a$, myelin content, and g-ratio. In early puberty, males had larger ${d}_a$ in the splenium and lower myelin content in the genu and body of the corpus callosum, compared with females. Overall, this work provides novel insights into developmental, pubertal, and cognitive correlates of individual differences in apparent axon diameter and myelin content in the developing human brain.

Integrating multimodality magnetic resonance imaging to the Allen Mouse Brain Common Coordinate Framework

High-resolution magnetic resonance imaging (MRI) affords unique image contrasts to nondestructively probe the tissue microstructure; validation of MRI findings with conventional histology is essential to better understand the MRI contrasts. However, the dramatic difference in the spatial resolution and image contrast of these two techniques impedes accurate comparison between MRI metrics and traditional histology. To better validate various MRI metrics, we acquired whole mouse brain multigradient recalled-echo and multishell diffusion MRI datasets at 25-μm isotropic resolution. The recently developed Allen Mouse Brain Common Coordinate Framework (CCFv3) provides opportunities to integrate multimodal and multiscale datasets of the whole mouse brain in a common three-dimensional (3D) space. The T2*, quantitative susceptibility mapping, diffusion tensor imaging, and neurite orientation dispersion and density imaging parameters were compared with both serial two-photon tomography images and 3D Nissl staining images in the CCFv3 at the same spatial resolution. The correlation between MRI and Nissl staining strongly depends on different metrics and different regions of the brain. Integrating different imaging modalities to the same space may substantially improve our understanding of the complexity of the brain at different scales.

A new open, high-resolution, multishell, diffusion-weighted imaging dataset of the living squirrel monkey

Although very well adapted to brain study, Magnetic Resonance Imaging (MRI) remains limited by the facilities and capabilities required to acquire data, especially for non-human primates. Addressing the data gaps resulting from these limitations requires making data more accessible and open. In contempt of the regular use of Saimiri sciureus in neuroscience research, in vivo diffusion has yet to be openly available for this species. Here we built and made openly available a unique new resource consisting of a high-resolution, multishell diffusion-weighted dataset in the anesthetized Saimiri sciureus. The data were acquired on 11 individuals with an 11.7 T MRI scanner (isotropic resolution of 400 µm³). This paper presents an overview of our dataset and illustrates some of its possible use through example analyses. To assess the quality of our data, we analyzed long-range connections (whole-brain tractography), microstructure (Neurite Orientation Dispersion and Density Imaging), and axon diameter in the corpus callosum (ActiveAx). Constituting an essential new resource for primate evolution studies, all data are openly available.

The influence of axonal beading and undulation on axonal diameter mapping

We consider the effect of non-cylindrical axonal shape on axonal diameter mapping with diffusion MRI. Practical sensitivity to axon diameter is attained at strong diffusion weightings b , where the deviation from the 1 / b scaling yields the finite transverse diffusivity, which is then translated into axon diameter. While axons are usually modeled as perfectly straight, impermeable cylinders, the local variations in diameter (caliber variation or beading) and direction (undulation) have been observed in microscopy data of human axons. Here we quantify the influence of cellular-level features such as caliber variation and undulation on axon diameter estimation. For that, we simulate the diffusion MRI signal in realistic axons segmented from 3-dimensional electron microscopy of a human brain sample. We then create artificial fibers with the same features and tune the amplitude of their caliber variations and undulations. Numerical simulations of diffusion in fibers with such tunable features show that caliber variations and undulations result in under- and over-estimation of axon diameters, correspondingly; this bias can be as large as 100%. 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.

Temporal Diffusion Ratio (TDR) for imaging restricted diffusion: Optimisation and pre-clinical demonstration

Temporal Diffusion Ratio (TDR) is a recently proposed dMRI technique (Dell'Acqua et al., proc. ISMRM 2019) which provides contrast between areas with restricted diffusion and areas either without restricted diffusion or with length scales too small for characterisation. Hence, it has a potential for informing on pore sizes, in particular the presence of large axon diameters or other cellular structures. TDR employs the signal from two dMRI acquisitions obtained with the same, large, b-value but with different diffusion gradient waveforms. TDR is advantageous as it employs standard acquisition sequences, does not make any assumptions on the underlying tissue structure and does not require any model fitting, avoiding issues related to model degeneracy. This work for the first time introduces and optimises the TDR method in simulation for a range of different tissues and scanner constraints and validates it in a pre-clinical demonstration. We consider both substrates containing cylinders and spherical structures, representing cell soma in tissue. Our results show that contrasting an acquisition with short gradient duration, short diffusion time and high gradient strength with an acquisition with long gradient duration, long diffusion time and low gradient strength, maximises the TDR contrast for a wide range of pore configurations. Additionally, in the presence of Rician noise, computing TDR from a subset (50% or fewer) of the acquired diffusion gradients rather than the entire shell as proposed originally further improves the contrast. In the last part of the work the results are demonstrated experimentally on rat spinal cord. In line with simulations, the experimental data shows that optimised TDR improves the contrast compared to non-optimised TDR. Furthermore, we find a strong correlation between TDR and histology measurements of axon diameter. In conclusion, we find that TDR has great potential and is a very promising alternative (or potentially complement) to model-based approaches for informing on pore sizes and restricted diffusion in general.

Mapping Human Brain Pathways: Challenges and Opportunities in the Integration of Scales

The human brain is composed of a complex web of pathways. Diffusion magnetic resonance (MR) tractography is a neuroimaging technique that relies on the principle of diffusion to reconstruct brain pathways. Its tractography is broadly applicable to a range of problems as it is amenable for study in individuals of any age and from any species. However, it is well known that this technique can generate biologically implausible pathways, especially in regions of the brain where multiple fibers cross. This review highlights potential misconnections in two cortico-cortical association pathways with a focus on the aslant tract and inferior frontal occipital fasciculus. The lack of alternative methods to validate observations from diffusion MR tractography means there is a need to develop new integrative approaches to trace human brain pathways. This review discusses integrative approaches in neuroimaging, anatomical, and transcriptional variation as having much potential to trace the evolution of human brain pathways.

Linear multi-scale modeling of diffusion MRI data: A framework for characterization of oriented structures across length scales

Diffusion-weighted magnetic resonance imaging (DW-MRI) has evolved to provide increasingly sophisticated investigations of the human brain's structural connectome in vivo. Restriction spectrum imaging (RSI) is a method that reconstructs the orientation distribution of diffusion within tissues over a range of length scales. In its original formulation, RSI represented the signal as consisting of a spectrum of Gaussian diffusion response functions. Recent technological advances have enabled the use of ultra-high b-values on human MRI scanners, providing higher sensitivity to intracellular water diffusion in the living human brain. To capture the complex diffusion time dependence of the signal within restricted water compartments, we expand upon the RSI approach to represent restricted water compartments with non-Gaussian response functions, in an extended analysis framework called linear multi-scale modeling (LMM). The LMM approach is designed to resolve length scale and orientation-specific information with greater specificity to tissue microstructure in the restricted and hindered compartments, while retaining the advantages of the RSI approach in its implementation as a linear inverse problem. Using multi-shell, multi-diffusion time DW-MRI data acquired with a state-of-the-art 3 T MRI scanner equipped with 300 mT/m gradients, we demonstrate the ability of the LMM approach to distinguish different anatomical structures in the human brain and the potential to advance mapping of the human connectome through joint estimation of the fiber orientation distributions and compartment size characteristics.

Both noise-floor and tissue compartment difference in diffusivity contribute to FA dependence on b-value in diffusion MRI

Noninvasive diffusion magnetic resonance imaging (dMRI) has been widely employed in both clinical and research settings to investigate brain tissue microstructure. Despite the evidence that dMRI-derived fractional anisotropy (FA) correlates with white matter properties, the metric is not specific. Recent studies have reported that FA is dependent on the b-value, and its origin has primarily been attributed to either the influence of microstructure or the noise-floor effect. A systematic investigation into the inter-relationship of these two effects is however still lacking. This study aims to quantify contributions of the reported differences in intra- and extra-neurite diffusivity to the observed changes in FA, in addition to the noise in measurements. We used in-vivo and post-mortem human brain imaging, as well as numerical simulations and histological validation, for this purpose. Our investigations reveal that the percentage difference of FA between b-values (pdFA) has significant positive associations with neurite density index (NDI), which is derived from in-vivo neurite orientation dispersion and density imaging (NODDI), or Bielschowsky's silver impregnation (BIEL) staining sections of fixed post-mortem human brain samples. Furthermore, such an association is found to be varied with Signal-to-Noise Ratio (SNR) level, indicating a nonlinear interaction effect between tissue microstructure and noise. Finally, a multicompartment model simulation revealed that these findings can be driven by differing diffusivities of intra- and extra-neurite compartments in tissue, with the noise-floor further amplifying the effect. In conclusion, both the differences in intra- and extra-neurite diffusivity and noise-floor effects significantly contribute to the FA difference associated with the b-value.

Axial and radial axonal diffusivities and radii from single encoding strongly diffusion-weighted MRI

We enable the estimation of the per-axon axial diffusivity from single encoding, strongly diffusion-weighted, pulsed gradient spin echo data. Additionally, we improve the estimation of the per-axon radial diffusivity compared to estimates based on spherical averaging. The use of strong diffusion weightings in magnetic resonance imaging (MRI) allows to approximate the signal in white matter as the sum of the contributions from only axons. At the same time, spherical averaging leads to a major simplification of the modeling by removing the need to explicitly account for the unknown distribution of axonal orientations. However, the spherically averaged signal acquired at strong diffusion weightings is not sensitive to the axial diffusivity, which cannot therefore be estimated although needed for modeling axons - especially in the context of multi-compartmental modeling. We introduce a new general method for the estimation of both the axial and radial axonal diffusivities at strong diffusion weightings based on kernel zonal modeling. The method could lead to estimates that are free from partial volume bias with gray matter or other isotropic compartments. The method is tested on publicly available data from the MGH Adult Diffusion Human Connectome project. We report reference values of axonal diffusivities based on 34 subjects, and derive estimates of axonal radii from only two shells. The estimation problem is also addressed from the angle of the required data preprocessing, the presence of biases related to modeling assumptions, current limitations, and future possibilities.

Diffusion radiomics for subtyping and clustering in autism spectrum disorder: A preclinical study

Autism spectrum disorder (ASD) is a highly prevalent, heterogenous neurodevelopmental disorder. Neuroimaging methods such as functional, structural, and diffusion MRI have been used to identify candidate imaging biomarkers for ASD, but current findings remain non-specific and likely arise from the heterogeneity present in ASD. To account for this, efforts to subtype ASD have emerged as a potential strategy for both the study of ASD and advancement of tailored behavioral therapies and therapeutics. Towards these ends, to improve upon current neuroimaging methods, we propose combining biologically sensitive neurite orientation dispersion and density index (NODDI) diffusion MR imaging with radiomics image processing to create a new methodological approach that, we hypothesize, can sensitively and specifically capture neurobiology. We demonstrate this method can sensitively distinguish differences between four genetically distinct rat models of ASD (Fmr1, Pten, Nrxn1, Disc1). Further, we demonstrate diffusion radiomic analyses hold promise for subtyping in ASD as we show unsupervised clustering of NODDI radiomic data generates clusters specific to the underlying genetic differences between the animal models. Taken together, our findings suggest the unique application of radiomic analysis on NODDI diffusion MRI may have the capacity to sensitively and specifically disambiguate the neurobiological heterogeneity present in the ASD population.

Neurite Orientation Dispersion and Density Imaging in Psychiatric Disorders: A Systematic Literature Review and a Technical Note

While major psychiatric disorders lack signature diagnostic neuropathologies akin to dementias, classic postmortem studies have established microstructural involvement, i.e., cellular changes in neurons and glia, as a key pathophysiological finding. Advanced magnetic resonance imaging techniques allow mapping of cellular tissue architecture and microstructural abnormalities in vivo, which holds promise for advancing our understanding of the pathophysiology underlying psychiatric disorders. Here, we performed a systematic review of case-control studies using neurite orientation dispersion and density imaging (NODDI) to assess brain microstructure in psychiatric disorders and a selective review of technical considerations in NODDI. Of the 584 potentially relevant articles, 18 studies met the criteria to be included in this systematic review. We found a general theme of abnormal gray and white matter microstructure across the diagnostic spectrum. We also noted significant variability in patterns of neurite density and fiber orientation within and across diagnostic groups, as well as associations between brain microstructure and phenotypical variables. NODDI has been successfully used to detect subtle microstructure abnormalities in patients with psychiatric disorders. Given that NODDI indices may provide a more direct link to pathophysiological processes, this method may not only contribute to advancing our mechanistic understanding of disease processes, it may also be well positioned for next-generation biomarker development studies.

Improving MR cell size imaging by inclusion of transcytolemmal water exchange

The goal of the current study is to include transcytolemmal water exchange in MR cell size imaging using the IMPULSED model for more accurate characterization of tissue cellular properties (e.g., apparent volume fraction of intracellular space v in ) and quantification of indicators of transcytolemmal water exchange. We propose a heuristic model that incorporates transcytolemmal water exchange into a multicompartment diffusion-based method (IMPULSED) that was developed previously to extract microstructural parameters (e.g., mean cell size d and apparent volume fraction of intracellular space v in ) assuming no water exchange. For t diff ≤ 5 ms, the water exchange can be ignored, and the signal model is the same as the IMPULSED model. For t diff ≥ 30 ms, we incorporated the modified Kärger model that includes both restricted diffusion and exchange between compartments. Using simulations and previously published in vitro cell data, we evaluated the accuracy and precision of model-derived parameters and determined how they are dependent on SNR and imaging parameters. The joint model provides more accurate d values for cell sizes ranging from 10 to 12 microns when water exchange is fast (e.g., intracellular water pre-exchange lifetime τ in ≤ 100 ms) than IMPULSED, and reduces the bias of IMPULSED-derived estimates of v in , especially when water exchange is relatively slow (e.g., τ in &gt; 200 ms). Indicators of transcytolemmal water exchange derived from the proposed joint model are linearly correlated with ground truth τ in values and can detect changes in cell membrane permeability induced by saponin treatment in murine erythroleukemia cancer cells. Our results suggest this joint model not only improves the accuracy of IMPULSED-derived microstructural parameters, but also provides indicators of water exchange that are usually ignored in diffusion models of tissues.

Advanced spinal cord MRI in multiple sclerosis: Current techniques and future directions

Spinal cord magnetic resonance imaging (MRI) has a central role in multiple sclerosis (MS) clinical practice for diagnosis and disease monitoring. Advanced MRI sequences capable of visualizing and quantifying tissue macro- and microstructure and reflecting different pathological disease processes have been used in MS research; however, the spinal cord remains under-explored, partly due to technical obstacles inherent to imaging this structure. We propose that the study of the spinal cord merits equal ambition in overcoming technical challenges, and that there is much information to be exploited to make valuable contributions to our understanding of MS. We present a narrative review on the latest progress in advanced spinal cord MRI in MS, covering in the first part structural, functional, metabolic and vascular imaging methods. We focus on recent studies of MS and those making significant technical steps, noting the challenges that remain to be addressed and what stands to be gained from such advances. Throughout we also refer to other works that presend more in-depth review on specific themes. In the second part, we present several topics that, in our view, hold particular potential. The need for better imaging of gray matter is discussed. We stress the importance of developing imaging beyond the cervical spinal cord, and explore the use of ultra-high field MRI. Finally, some recommendations are given for future research, from study design to newer developments in analysis, and the need for harmonization of sequences and methods within the field. This review is aimed at researchers and clinicians with an interest in gaining an overview of the current state of advanced MRI research in this field and what is primed to be the future of spinal cord imaging in MS research.

Advanced diffusion MRI and image texture analysis detect widespread brain structural differences between relapsing-remitting and secondary progressive multiple sclerosis

Introduction

Disease development in multiple sclerosis (MS) causes dramatic structural changes, but the exact changing patterns are unclear. Our objective is to investigate the differences in brain structure locally and spatially between relapsing-remitting MS (RRMS) and its advanced form, secondary progressive MS (SPMS), through advanced analysis of diffusion magnetic resonance imaging (MRI) and image texture.

Methods

A total of 20 patients with RRMS and nine patients with SPMS from two datasets underwent 3T anatomical and diffusion tensor imaging (DTI). The DTI was harmonized, augmented, and then modeled, which generated six voxel- and sub-voxel-scale measures. Texture analysis focused on T2 and FLAIR MRI, which produced two phase-based measures, namely, phase congruency and weighted mean phase. Data analysis was 3-fold, i.e., histogram analysis of whole-brain normal appearing white matter (NAWM); region of interest (ROI) analysis of NAWM and lesions within three critical white matter tracts, namely, corpus callosum, corticospinal tract, and optic radiation; and along-tract statistics. Furthermore, by calculating the z-score of core-rim pathology within lesions based on diffusion measures, we developed a novel method to define chronic active lesions and compared them between cohorts.

Results

Histogram features from diffusion and all but one texture measure differentiated between RRMS and SPMS. Within-tract ROI analysis detected cohort differences in both NAWM and lesions of the corpus callosum body in three measures of neurite orientation and anisotropy. Along-tract statistics detected cohort differences from multiple measures, particularly lesion extent, which increased significantly in SPMS in posterior corpus callosum and optic radiations. The number of chronic active lesions were also significantly higher (by 5-20% over z-scores 0.5 and 1.0) in SPMS than RRMS based on diffusion anisotropy, neurite content, and diameter.

Conclusion

Advanced diffusion MRI and texture analysis may be promising approaches for thorough understanding of brain structural changes from RRMS to SPMS, thereby providing new insight into disease development mechanisms in MS.

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.

Post mortem mapping of connectional anatomy for the validation of diffusion MRI

Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.