An analytical model of restricted diffusion in bovine optic nerve

Considerations and recommendations from the ISMRM Diffusion Study Group for preclinical diffusion MRI: Part 3-Ex vivo imaging: Data processing, comparisons with microscopy, and tractography

Preclinical diffusion MRI (dMRI) has proven value in methods development and validation, characterizing the biological basis of diffusion phenomena, and comparative anatomy. While dMRI enables in vivo non-invasive characterization of tissue, ex vivo dMRI is increasingly being used to probe tissue microstructure and brain connectivity. Ex vivo dMRI has several experimental advantages that facilitate high spatial resolution and high SNR images, cutting-edge diffusion contrasts, and direct comparison with histological data as a methodological validation. However, there are a number of considerations that must be made when performing ex vivo experiments. The steps from tissue preparation, image acquisition and processing, and interpretation of results are complex, with many decisions that not only differ dramatically from in vivo imaging of small animals, but ultimately affect what questions can be answered using the data. This work concludes a three-part series of recommendations and considerations for preclinical dMRI. Herein, we describe best practices for dMRI of ex vivo tissue, with a focus on image pre-processing, data processing, and comparisons with microscopy. In each section, we attempt to provide guidelines and recommendations but also highlight areas for which no guidelines exist (and why), and where future work should lie. We end by providing guidelines on code sharing and data sharing and point toward open-source software and databases specific to small animal and ex vivo imaging.

LLPS REDIFINE allows the biophysical characterization of multicomponent condensates without tags or labels

Liquid-liquid phase separation (LLPS) phenomenon plays a vital role in multiple cell biology processes, providing a mechanism to concentrate biomolecules and promote cellular reactions locally. Despite its significance in biology, there is a lack of conventional techniques suitable for studying biphasic samples in their biologically relevant form. Here, we present a label-free and non-invasive approach to characterize biomolecular condensates termed LLPS REstricted DIFusion of INvisible speciEs (REDIFINE). Relying on diffusion NMR measurements, REDIFINE exploits the exchange dynamics between molecules in the condensed and dispersed phases to determine not only diffusion constants and the fractions in both phases but also the average radius of the condensed droplets and the exchange rate between the phases. Observing proteins, RNAs, water, as well as small molecules, and even assessing the concentrations of biomolecules in both phases, REDIFINE analysis allows a rapid biophysical characterization of multicomponent condensates which is important to understand their functional roles. In comparing multiple systems, REDIFINE reveals that folded RNA-binding proteins form smaller and more dynamic droplets compared to the disordered ones.

Enhanced structural brain connectivity analyses using high diffusion-weighting strengths

Tractography is a unique modality for the in vivo measurement of structural connectivity, crucial for understanding brain networks and neurological conditions. With increasing b-value, the diffusion-weighting signal becomes primarily sensitive to the intra-axonal signal. However, it remains unclear how tractography is affected by this observation. Here, using open-source datasets, we showed that at high b-values, DWI reduces the uncertainty in estimating fiber orientations. Specifically, we found the ratio of biologically-meaningful longer-range connections increases, accompanied with downstream impact of redistribution of connectome and network metrics. However, when going beyond b = 6000 s/mm2, the loss of SNR imposed a penalty. Lastly, we showed that the data reaches satisfactory reproducibility with b-values above 1200 s/mm2. Overall, the results suggest that using b-values above 2500 s/mm2 is essential for more accurate connectome reconstruction by reducing uncertainty in fiber orientation estimation, supporting the use of higher b-value protocols in standard diffusion MRI scans and pipelines.

Imaging the structural connectome with hybrid MRI-microscopy tractography

Mapping how neurons are structurally wired into whole-brain networks can be challenging, particularly in larger brains where 3D microscopy is not available. Multi-modal datasets combining MRI and microscopy provide a solution, where high resolution but 2D microscopy can be complemented by whole-brain but lowresolution MRI. However, there lacks unified approaches to integrate and jointly analyse these multi-modal data in an insightful way. To address this gap, we introduce a data-fusion method for hybrid MRI-microscopy fibre orientation and connectome reconstruction. Specifically, we complement precise "in-plane" orientations from microscopy with "through-plane" information from MRI to construct 3D hybrid fibre orientations at resolutions far exceeding that of MRI whilst preserving microscopy's myelin specificity, resulting in superior fibre tracking. Our method is openly available, can be deployed on standard 2D microscopy, including different microscopy contrasts, and is species agnostic, facilitating neuroanatomical investigation in both animal models and human brains.

When CEST meets diffusion: Multi-echo diffusion-encoded CEST (dCEST) MRI to measure intracellular and extracellular CEST signal distributions

PURPOSE

To develop a multi-echo, diffusion-encoded chemical exchange saturation transfer (dCEST) imaging technique for estimating the intracellular and extracellular/intravascular contributions to the conventional CEST signal.

METHODS

A dCEST pulse sequence was developed to quantify the signal fractions, transverse relaxation times (T2), and apparent diffusion coefficient (ADC) of the intracellular and extracellular/intravascular water compartments. dCEST images were acquired across a wide range of TE, b-values, RF saturation strengths, and frequency offsets. The data were analyzed using a two-compartment model with distinct diffusivities and T2 values. Intracellular and extracellular fractions of conventional water-saturation spectra (Z-spectra) and corresponding amide proton transfer (APT) signals were estimated from human brain scans of healthy volunteers at 3 T.

RESULTS

The multi-echo diffusion results showed that the intracellular water fractions were significantly higher than the extracellular water fractions, whereas the intracellular T2 values were shorter than those of the extracellular/intravascular compartments. The ADC for the intracellular compartment was significantly lower than that of the extracellular compartment. The dCEST analysis showed that the average intracellular and extracellular fractions of the Z-spectra were 85 ± 7% and 15 ± 4%, respectively. The overall intracellular APT-weighted values were higher than the total (i.e., intracellular + extracellular) APT-weighted values.

CONCLUSIONS

The dCEST imaging technique provides valuable insight into the source of signals in conventional CEST MRI, offering potential utility for clinical applications.

Ex vivo massively multidimensional diffusion-relaxation correlation MRI: scan-rescan reproducibility and caveats

Massively multidimensional diffusion-relaxation correlation MRI (MMD-MRI) provides information beyond the traditional voxel-averaged metric that may better characterize the microstructural characteristics of biological tissues. MMD-MRI reproducibility has been established in clinical settings, but has yet to be thoroughly evaluated under preclinical conditions, where superior hardware and modulated gradient waveforms enhance its performance. In this study, we investigate the reproducibility of MMD-MRI on a micro-imaging system using ex vivo mouse brains. Notably, the estimated signal fractions of intra-voxel spectral components in the MD-MRI distribution, corresponding to white and gray matter, along with the frequency-dependent parameters, demonstrated high reproducibility. We identified bias between scan and rescan in some of the metrics, which we attribute to the time gap between repeated scans pointing to a long-time progressive fixation effect. We compare our results with in vivo results from clinical scanners and show the reproducibility of diffusion frequency-dependent metrics to benefit from the improved gradient hardware on our preclinical setup. Our results inform future micro-imaging ex vivo MMD-MRI studies of the reproducibility of MMD-MRI metrics and their dependence on fixation time.

Revealing membrane integrity and cell size from diffusion kurtosis time dependence

PURPOSE

The nonmonotonic dependence of diffusion kurtosis on diffusion time has been observed in biological tissues, yet its relation to membrane integrity and cellular geometry remains to be clarified. Here we establish and explain the characteristic asymmetric shape of the kurtosis peak. We also derive the relation between the peak timet peak $$ {t}_{\mathrm{peak}} $$ , when kurtosis reaches its maximum, and tissue parameters.

METHODS

The peak shape and its positiont peak $$ {t}_{\mathrm{peak}} $$ qualitatively follow from the adiabatic extension of the Kärger model onto the case of intra-cellular diffusivity time-dependence. This intuition is corroborated by the effective medium theory-based calculation, as well as by Monte Carlo simulations of diffusion and exchange in randomly and densely packed spheres for various values of permeability, cell fractions and sizes, and intrinsic diffusivity.

RESULTS

We establish thatt peak $$ {t}_{\mathrm{peak}} $$ is proportional to the geometric mean of two characteristic time scales: extra-cellular correlation time (determined by cell size) and intra-cellular residence time (determined by membrane permeability). When exchange is barrier-limited, the peak shape approaches a universal scaling form determined by the ratiot / t peak $$ t/{t}_{\mathrm{peak}} $$ .

CONCLUSION

Numerical simulations and theory provide an interpretation of a specific feature of kurtosis time-dependence, offering a potential biomarker for in vivo evaluation of pathology by disentangling the functional (permeability) and structural (cell size) integrity in tissues. This is relevant as the time-dependent diffusion cumulants are sensitive to pathological changes in membrane integrity and cellular structure in diseases, such as ischemic stroke, tumors, and Alzheimer's disease.

Investigating microstructural changes between in vivo and perfused ex vivo marmoset brains using oscillating gradient and b-tensor encoded diffusion MRI at 9.4 T

PURPOSE

To investigate microstructural alterations induced by perfusion fixation in brain tissues using advanced diffusion MRI techniques and estimate their potential impact on the application of ex vivo models to in vivo microstructure.

METHODS

We used oscillating gradient spin echo (OGSE) and b-tensor encoding diffusion MRI to examine in vivo and ex vivo microstructural differences in the marmoset brain. OGSE was used to shorten effective diffusion times, whereas b-tensor encoding allowed for the differentiation of isotropic and anisotropic kurtosis. Additionally, we performed Monte Carlo simulations to estimate the potential microstructural changes in the tissues.

RESULTS

We report large changes (˜50%-60%) in kurtosis frequency dispersion (OGSE) and in both anisotropic and isotropic kurtosis (b-tensor encoding) after perfusion fixation. Structural MRI showed an average volume reduction of about 10%. Monte Carlo simulations indicated that these alterations could likely be attributed to extracellular fluid loss possibly combined with axon beading and increased dot compartment signal fraction. Little evidence was observed for reductions in axonal caliber.

CONCLUSION

Our findings shed light on advanced MRI parameter changes that are induced by perfusion fixation and potential microstructural sources for these changes. This work also suggests that caution should be exercised when applying ex vivo models to infer in vivo tissue microstructure, as significant differences may arise.

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.

Accelerated, Physics-Inspired Inference of Skeletal Muscle Microstructure From Diffusion-Weighted MRI

Muscle health is a critical component of overall health and quality of life. However, current measures of skeletal muscle health take limited account of microstructural variations within muscle, which play a crucial role in mediating muscle function. To address this, we present a physics-inspired, machine learning-based framework for the non-invasive estimation of microstructural organization in skeletal muscle from diffusion-weighted MRI (dMRI) in an uncertainty-aware manner. To reduce the computational expense associated with direct numerical simulations of dMRI physics, a polynomial meta-model is developed that accurately represents the input/output relationships of a high-fidelity numerical model. This meta-model is used to develop a Gaussian process (GP) model that provides voxel-wise estimates and confidence intervals of microstructure organization in skeletal muscle. Given noise-free data, the GP model accurately estimates microstructural parameters. In the presence of noise, the diameter, intracellular diffusion coefficient, and membrane permeability are accurately estimated with narrow confidence intervals, while volume fraction and extracellular diffusion coefficient are poorly estimated and exhibit wide confidence intervals. A reduced-acquisition GP model, consisting of one-third the diffusion-encoding measurements, is shown to predict parameters with similar accuracy to the original model. The fiber diameter and volume fraction estimated by the reduced GP model is validated via histology, with both parameters accurately estimated, demonstrating the capability of the proposed framework as a promising non-invasive tool for assessing skeletal muscle health and function.

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 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-500 ms). 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 ± 4 m s ,f N G = 0.71 ± 0.01 ,R e f f = 1.10 ± 0.01 μ m , andκ eff = 0.21 ± 0.06 μ m / m s , 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.

Validation of conductivity tensor imaging against diffusion tensor magnetic resonance electrical impedance tomography

Diffusion tensor magnetic resonance electrical impedance tomography (DT-MREIT) and electrodeless conductivity tensor imaging (CTI) are two emerging modalities that can quantify low-frequency tissue anisotropic conductivity properties by assuming similar properties underlie ionic mobility and water diffusion. While both methods have potential applications to estimating neuro-modulation fields or formulating forward models used for electrical source imaging, a direct comparison of the two modalities has not yet been performed in-vitro or in-vivo. Therefore, the aim of this study was to test the equivalence of these two modalities. We scanned a tissue phantom and the head of human subject using DT-MREIT and CTI protocols and reconstructed conductivity tensor and effective low frequency conductivities. We found both gray and white matter conductivities recovered by each technique were equivalent within 0.05 S/m. Both DT-MREIT and CTI require multiple processing steps, and we further assess the effects of each factor on reconstructions and evaluate the extent to which different measurement mechanisms potentially cause discrepancies between the two methods. Finally, we discuss the implications for spectral models of measuring conductivity using these techniques. The study further establishes the credibility of CTI as an electrodeless non-invasive method of measuring low frequency conductivity properties.

Utilizing multicompartmental restriction spectrum magnetic resonance imaging for liver fibrosis characterization in a mouse model

BACKGROUND

Currently, an advanced imaging method may be necessary for magnetic resonance imaging (MRI) to diagnosis and quantify liver fibrosis (LF).

PURPOSE

To evaluate the feasibility of the multicompartmental restriction spectrum imaging (RSI) model to characterize LF in a mouse model.

METHODS

Thirty mice with carbon tetrachloride (CCl4)-induced LF and eight control mice were investigated using multi-b-value (ranging from 0 to 2000 s/mm2) diffusion-weighted imaging (DWI) on a 3T scanner. DWI data were processed using RSI model (2-5 compartments) with the Bayesian Information Criterion (BIC) determining the optimal model. Conventional ADC value and signal fraction of each compartment in the optimal RSI model were compared across groups. Receiver operating characteristics (ROC) curve analysis was performed to determine the diagnosis performances of different parameters, while Spearman correlation analysis was employed to investigate the correlation between different tissue compartments and the stage of LF.

RESULTS

According to BIC results, a 4-compartment RSI model (RSI4) with optimal ADCs of 0.471 × 10-3, 1.653 × 10-3, 9.487 × 10-3, and > 30 × 10-3, was the optimal model to characterize LF. Significant differences in signal contribution fraction of the C1 and C3 compartments were observed between LF and control groups (P = 0.018 and 0.003, respectively). ROC analysis showed that RSI4-C3 was the most effective single diffusion parameter for characterizing LF (AUC = 0.876, P = 0.003). Furthermore, the combination of ADC values and RSI4-C3 value increased the diagnosis performance significantly (AUC = 0.894, P = 0.002).

CONCLUSION

The 4-compartment RSI model has the potential to distinguish LF from the control group based on diffusion parameters. RSI4-C3 showed the highest diagnostic performance among all the parameters. The combination of ADC and RSI4-C3 values further improved the discrimination performance.

Diffusion tensor imaging metrics as natural markers of multiple sclerosis-induced brain disorders with a low Expanded Disability Status Scale score

Non-invasive and effective differentiation along with determining the degree of deviations compared to the healthy cohort is important in the case of various brain disorders, including multiple sclerosis (MS). Evaluation of the effectiveness of diffusion tensor metrics (DTM) in 3T DTI for recording MS-related deviations was performed using a time-acceptable MRI protocol with unique comprehensive detection of systematic errors related to spatial heterogeneity of magnetic field gradients. In a clinical study, DTMs were acquired in segmented regions of interest (ROIs) for 50 randomly selected healthy controls (HC) and 50 multiple sclerosis patients. Identical phantom imaging was performed for each clinical measurement to estimate and remove the influence of systematic errors using the b-matrix spatial distribution in the DTI (BSD-DTI) technique. In the absence of statistically significant differences due to age in healthy volunteers and patients with multiple sclerosis, the existence of significant differences between groups was proven using DTM. Moreover, a statistically significant impact of spatial systematic errors occurs for all ROIs and DTMs in the phantom and for approximately 90 % in the HC and MS groups. In the case of a single patient measurement, this appears for all the examined ROIs and DTMs. The obtained DTMs effectively discriminate healthy volunteers from multiple sclerosis patients with a low mean score on the Expanded Disability Status Scale. The magnitude of the group differences is typically significant, with an effect size of approximately 0.5, and similar in both the standard approach and after elimination of systematic errors. Differences were also observed between metrics obtained using these two approaches. Despite a small alterations in mean DTMs values for groups and ROIs (1-3 %), these differences were characterized by a huge effect (effect size ∼0.8 or more). These findings indicate the importance of determining the spatial distribution of systematic errors specific to each MR scanner and DTI acquisition protocol in order to assess their impact on DTM in the ROIs examined. This is crucial to establish accurate DTM values for both individual patients and mean values for a healthy population as a reference. This approach allows for an initial reliable diagnosis based on DTI metrics.

Spatiotemporal patterns of cortical microstructural maturation in children and adolescents with diffusion MRI

Neocortical maturation is a dynamic process that proceeds in a hierarchical manner; however, the spatiotemporal organization of cortical microstructure with diffusion MRI has yet to be fully defined. This study characterized cortical microstructural maturation using diffusion MRI (fwe-diffusion tensor imaging [DTI] and neurite orientation dispersion and density imaging [NODDI] multicompartment modeling) in a cohort of 637 children and adolescents between 8 and 21 years of age. We found spatially heterogeneous developmental patterns broadly demarcated into functional domains where NODDI metrics increased, and fwe-DTI metrics decreased with age. By applying nonlinear growth models in a vertex-wise analysis, we observed a general posterior-to-anterior pattern of maturation, where the fwe-DTI measures mean diffusivity and radial diffusivity reached peak maturation earlier than the NODDI metrics neurite density index. Using non-negative matrix factorization, we found occipito-parietal cortical regions that correspond to lower order sensory domains mature earlier than fronto-temporal higher order association domains. Our findings corroborate previous histological and neuroimaging studies that show spatially varying patterns of cortical maturation that may reflect unique developmental processes of cytoarchitectonically determined regional patterns of change.

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.

MRI with generalized diffusion encoding reveals damaged white matter in patients previously hospitalized for COVID-19 and with persisting symptoms at follow-up

There is mounting evidence of the long-term effects of COVID-19 on the central nervous system, with patients experiencing diverse symptoms, often suggesting brain involvement. Conventional brain MRI of these patients shows unspecific patterns, with no clear connection of the symptomatology to brain tissue abnormalities, whereas diffusion tensor studies and volumetric analyses detect measurable changes in the brain after COVID-19. Diffusion MRI exploits the random motion of water molecules to achieve unique sensitivity to structures at the microscopic level, and new sequences employing generalized diffusion encoding provide structural information which are sensitive to intravoxel features. In this observational study, a total of 32 persons were investigated: 16 patients previously hospitalized for COVID-19 with persisting symptoms of post-COVID condition (mean age 60 years: range 41-79, all male) at 7-month follow-up and 16 matched controls, not previously hospitalized for COVID-19, with no post-COVID symptoms (mean age 58 years, range 46-69, 11 males). Standard MRI and generalized diffusion encoding MRI were employed to examine the brain white matter of the subjects. To detect possible group differences, several tissue microstructure descriptors obtainable with the employed diffusion sequence, the fractional anisotropy, mean diffusivity, axial diffusivity, radial diffusivity, microscopic anisotropy, orientational coherence (Cc) and variance in compartment's size (CMD) were analysed using the tract-based spatial statistics framework. The tract-based spatial statistics analysis showed widespread statistically significant differences (P < 0.05, corrected for multiple comparisons using the familywise error rate) in all the considered metrics in the white matter of the patients compared to the controls. Fractional anisotropy, microscopic anisotropy and Cc were lower in the patient group, while axial diffusivity, radial diffusivity, mean diffusivity and CMD were higher. Significant changes in fractional anisotropy, microscopic anisotropy and CMD affected approximately half of the analysed white matter voxels located across all brain lobes, while changes in Cc were mainly found in the occipital parts of the brain. Given the predominant alteration in microscopic anisotropy compared to Cc, the observed changes in diffusion anisotropy are mostly due to loss of local anisotropy, possibly connected to axonal damage, rather than white matter fibre coherence disruption. The increase in radial diffusivity is indicative of demyelination, while the changes in mean diffusivity and CMD are compatible with vasogenic oedema. In summary, these widespread alterations of white matter microstructure are indicative of vasogenic oedema, demyelination and axonal damage. These changes might be a contributing factor to the diversity of central nervous system symptoms that many patients experience after COVID-19.

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.

A Method for In-Vivo Mapping of Axonal Diameter Distributions in the Human Brain Using Diffusion-Based Axonal Spectrum Imaging (AxSI)

In this paper we demonstrate a generalized and simplified pipeline called axonal spectrum imaging (AxSI) for in-vivo estimation of axonal characteristics in the human brain. Whole-brain estimation of the axon diameter, in-vivo and non-invasively, across all fiber systems will allow exploring uncharted aspects of brain structure and function relations with emphasis on connectivity and connectome analysis. While axon diameter mapping is important in and of itself, its correlation with conduction velocity will allow, for the first time, the explorations of information transfer mechanisms within the brain. We demonstrate various well-known aspects of axonal morphometry (e.g., the corpus callosum axon diameter variation) as well as other aspects that are less explored (e.g., axon diameter-based separation of the superior longitudinal fasciculus into segments). Moreover, we have created an MNI based mean axon diameter map over the entire brain for a large cohort of subjects providing the reference basis for future studies exploring relation between axon properties, its connectome representation, and other functional and behavioral aspects of the brain.

Diffusion time dependency of extracellular diffusion

PURPOSE

To quantify the variations of the power-law dependences on diffusion time t or gradient frequencyf $$ f $$ of extracellular water diffusion measured by diffusion MRI (dMRI).

METHODS

Model cellular systems containing only extracellular water were used to investigate thet / f $$ t/f $$ dependence ofD ex $$ {D}_{ex} $$ , the extracellular diffusion coefficient. Computer simulations used a randomly packed tissue model with realistic intracellular volume fractions and cell sizes. DMRI measurements were performed on samples consisting of liposomes containing heavy water(D2 O, deuterium oxide) dispersed in regular water (H2 O).D ex $$ {D}_{ex} $$ was obtained over a broadt $$ t $$ range (∼1-1000 ms) and then fit power-law equationsD ex ( t ) = D const + const · t - ϑ t $$ {D}_{ex}(t)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {t}^{-{\vartheta}_t} $$ andD ex ( f ) = D const + const · f ϑ f $$ {D}_{ex}(f)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {f}^{\vartheta_f} $$ .

RESULTS

Both simulated and experimental results suggest that no single power-law adequately describes the behavior ofD ex $$ {D}_{ex} $$ over the range of diffusion times of most interest in practical dMRI. Previous theoretical predictions are accurate over only limitedt $$ t $$ ranges; for example,θ t = θ f = - 1 2 $$ {\theta}_t={\theta}_f=-\frac{1}{2} $$ is valid only for short times, whereasθ t = 1 $$ {\theta}_t=1 $$ orθ f = 3 2 $$ {\theta}_f=\frac{3}{2} $$ is valid only for long times but cannot describe other ranges simultaneously. For the specifict $$ t $$ range of 5-70 ms used in typical human dMRI measurements,θ t = θ f = 1 $$ {\theta}_t={\theta}_f=1 $$ matches the data well empirically.

CONCLUSION

The optimal power-law fit of extracellular diffusion varies with diffusion time. The dependency obtained at short or longt $$ t $$ limits cannot be applied to typical dMRI measurements in human cancer or liver. It is essential to determine the appropriate diffusion time range when modeling extracellular diffusion in dMRI-based quantitative microstructural imaging.

Measuring water exchange on a preclinical MRI system using filter exchange and diffusion time dependent kurtosis imaging

PURPOSE

Filter exchange imaging (FEXI) and diffusion time (t)-dependent diffusion kurtosis imaging (DKI(t)) are both sensitive to water exchange between tissue compartments. The restrictive effects of tissue microstructure, however, introduce bias to the exchange rate obtained by these two methods, as their interpretation conventionally rely on the Kärger model of barrier limited exchange between Gaussian compartments. Here, we investigated whether FEXI and DKI(t) can provide comparable exchange rates in ex vivo mouse brains.

THEORY AND METHODS

FEXI and DKI(t) data were acquired from ex vivo mouse brains on a preclinical MRI system. Phase cycling and negative slice prewinder gradients were used to minimize the interferences from imaging gradients.

RESULTS

In the corpus callosum, apparent exchange rate (AXR) from FEXI correlated with the exchange rate (the inverse of exchange time, 1/τex ) from DKI(t) along the radial direction. In comparison, discrepancies between FEXI and DKI(t) were found in the cortex due to low filter efficiency and confounding effects from tissue microstructure.

CONCLUSION

The results suggest that FEXI and DKI(t) are sensitive to the same exchange processes in white matter when separated from restrictive effects of microstructure. The complex microstructure in gray matter, with potential exchange among multiple compartments and confounding effects of microstructure, still pose a challenge for FEXI and DKI(t).

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.

Adult lifespan maturation and degeneration patterns in gray and white matter: A mean apparent propagator (MAP) MRI study

The relationship between brain microstructure and aging has been the subject of intense study, with diffusion MRI perhaps the most effective modality for elucidating these associations. Here, we used the mean apparent propagator (MAP)-MRI framework, which is suitable to characterize complex microstructure, to investigate age-related cerebral differences in a cohort of cognitively unimpaired participants and compared the results to those derived using diffusion tensor imaging. We studied MAP-MRI metrics, among them the non-Gaussianity (NG) and propagator anisotropy (PA), and established an opposing pattern in white matter of higher NG alongside lower PA among older adults, likely indicative of axonal degradation. In gray matter, however, these two indices were consistent with one another, and exhibited regional pattern heterogeneity compared to other microstructural parameters, which could indicate fewer neuronal projections across cortical layers along with an increased glial concentration. In addition, we report regional variations in the magnitude of age-related microstructural differences consistent with the posterior-anterior shift in aging paradigm. These results encourage further investigations in cognitive impairments and neurodegeneration.

Spatiotemporal patterns of cortical microstructural maturation in children and adolescents with diffusion MRI

Neocortical maturation is a dynamic process that proceeds in a hierarchical manner; however, the spatiotemporal organization of cortical microstructure with diffusion MRI has yet to be fully defined. This study characterized cortical microstructural maturation using diffusion MRI (fwe-DTI and NODDI multi-compartment modeling) in a cohort of 637 children and adolescents between 8 and 21 years of age. We found spatially heterogeneous developmental patterns broadly demarcated into functional domains where NODDI metrics increased and fwe-DTI metrics decreased with age. Using non-negative matrix factorization, we found cortical regions that correspond to lower-order sensory regions mature earlier than higher-order association regions. Our findings corroborate previous histological and neuroimaging studies that show spatially-varying patterns of cortical maturation that may reflect unique developmental processes of cytoarchitectonically-determined regional patterns of change.

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.

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 MRI with pulsed and free gradient waveforms: Effects of restricted diffusion and exchange

Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm and exchange rates in the simulated range of 0 to 20 s - 1 , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.

COnstrained Reference frame diffusion TEnsor Correlation Spectroscopic (CORTECS) MRI: A practical framework for high-resolution diffusion tensor distribution imaging

High-resolution imaging studies have consistently shown that in cortical tissue water diffuses preferentially along radial and tangential orientations with respect to the cortical surface, in agreement with histology. These dominant orientations do not change significantly even if the relative contributions from microscopic water pools to the net voxel signal vary across experiments that use different diffusion times, b-values, TEs, and TRs. With this in mind, we propose a practical new framework for imaging non-parametric diffusion tensor distributions (DTDs) by constraining the microscopic diffusion tensors of the DTD to be diagonalized using the same orthonormal reference frame of the mesoscopic voxel. In each voxel, the constrained DTD (cDTD) is completely determined by the correlation spectrum of the microscopic principal diffusivities associated with the axes of the voxel reference frame. Consequently, all cDTDs are inherently limited to the domain of positive definite tensors and can be reconstructed efficiently using Inverse Laplace Transform methods. Moreover, the cDTD reconstruction can be performed using only data acquired efficiently with single diffusion encoding, although it also supports datasets with multiple diffusion encoding. In tissues with a well-defined architecture, such as the cortex, we can further constrain the cDTD to contain only cylindrically symmetric diffusion tensors and measure the 2D correlation spectra of principal diffusivities along the radial and tangential orientation with respect to the cortical surface. To demonstrate this framework, we perform numerical simulations and analyze high-resolution dMRI data from a fixed macaque monkey brain. We estimate 2D cDTDs in the cortex and derive, in each voxel, the marginal distributions of the microscopic principal diffusivities, the corresponding distributions of the microscopic fractional anisotropies and mean diffusivities along with their 2D correlation spectra to quantify the cDTD shape-size characteristics. Signal components corresponding to specific bands in these cDTD-derived spectra show high specificity to cortical laminar structures observed with histology. Our framework drastically simplifies the measurement of non-parametric DTDs in high-resolution datasets with mesoscopic voxel sizes much smaller than the radius of curvature of the underlying anatomy, e.g., cortical surface, and can be applied retrospectively to analyze existing diffusion MRI data from fixed cortical tissues.

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.

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.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

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

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

Diffusion-based microstructure models in brain tumours: Fitting in presence of a model-microstructure mismatch

•

We analyzed the performance of NODDI, SMT and DKI inside the tumoral lesion.

•

Goodness of fit were comparable to normal tissue, for DKI and NODDI.

•

Parameter precision was similar to normal tissues for all quantified metrics.

•

Parameters should not be given their healthy physiological meaning in the tumour.

•

The three models are usable as signal representations of the tumoral tissue.

Diffusion-based biophysical models have been used in several recent works to study the microenvironment of brain tumours. While the pathophysiological interpretation of the parameters of these models remains unclear, their use as signal representations may yield useful biomarkers for monitoring the treatment and the progression of this complex and heterogeneous disease. Up to now, however, no study was devoted to assessing the mathematical stability of these approaches in cancerous brain regions. To this end, we analyzed in 11 brain tumour patients the fitting results of two microstructure models (Neurite Orientation Dispersion and Density Imaging and the Spherical Mean Technique) and of a signal representation (Diffusion Kurtosis Imaging) to compare the reliability of their parameter estimates in the healthy brain and in the tumoral lesion. The framework of our between-tissue analysis included the computation of 1) the residual sum of squares as a goodness-of-fit measure 2) the standard deviation of the models’ derived metrics and 3) models’ sensitivity functions to analyze the suitability of the employed protocol for parameter estimation in the different microenvironments. Our results revealed no issues concerning the fitting of the models in the tumoral lesion, with similar goodness of fit and parameter precisions occurring in normal appearing and pathological tissues. Lastly, with the aim of highlight possible biomarkers, in our analysis we briefly discuss the correlation between the metrics of the three techniques, identifying groups of indices which are significantly collinear in all tissues and thus provide no additional information when jointly used in data-driven analyses.

Fractional Diffusion with Geometric Constraints: Application to Signal Decay in Magnetic Resonance Imaging (MRI)

We investigate diffusion in three dimensions on a comb-like structure in which the particles move freely in a plane, but, out of this plane, are constrained to move only in the perpendicular direction. This model is an extension of the two-dimensional version of the comb model, which allows diffusion along the backbone when the particles are not in the branches. We also consider memory effects, which may be handled with different fractional derivative operators involving singular and non-singular kernels. We find exact solutions for the particle distributions in this model that display normal and anomalous diffusion regimes when the mean-squared displacement is determined. As an application, we use this model to fit the anisotropic diffusion of water along and across the axons in the optic nerve using magnetic resonance imaging. The results for the observed diffusion times (8 to 30 milliseconds) show an anomalous diffusion both along and across the fibers.

Disentangling the effects of restriction and exchange with diffusion exchange spectroscopy

Diffusion exchange spectroscopy (DEXSY) is a multidimensional NMR technique that can reveal how water molecules exchange between compartments within heterogeneous media, such as biological tissue. Data from DEXSY experiments is typically processed using numerical inverse Laplace transforms (ILTs) to produce a diffusion-diffusion spectrum. A tacit assumption of this ILT approach is that the signal behavior is Gaussian - i.e., the spin echo intensity decays exponentially with the degree of diffusion weighting. The assumptions that underlie Gaussian signal behavior may be violated, however, depending on the gradient strength applied and the sample under study. We argue that non-Gaussian signal behavior due to restrictions is to be expected in the study of biological tissue using diffusion NMR. Further, we argue that this signal behavior can produce confounding features in the diffusion-diffusion spectra obtained from numerical ILTs of DEXSY data - entangling the effects of restriction and exchange. Specifically, restricted signal behavior can result in broadening of peaks and in the appearance of illusory exchanging compartments with distributed diffusivities, which pearl into multiple peaks if not highly regularized. We demonstrate these effects on simulated data. That said, we suggest the use of features in the signal acquisition domain that can be used to rapidly probe exchange without employing an ILT. We also propose a means to characterize the non-Gaussian signal behavior due to restrictions within a sample using DEXSY measurements with a near zero mixing time or storage interval. We propose a combined acquisition scheme to independently characterize restriction and exchange with various DEXSY measurements, which we term Restriction and Exchange from Equally-weighted Double and Single Diffusion Encodings (REEDS-DE). We test this method on ex vivo neonatal mouse spinal cord - a sample consisting primarily of gray matter - using a low-field, static gradient NMR system. In sum, we highlight critical shortcomings of prevailing DEXSY analysis methods that conflate the effects of restriction and exchange, and suggest a viable experimental approach to disentangle them.

The direction-dependence of apparent water exchange rate in human white matter

Transmembrane water exchange is a potential biomarker in the diagnosis and understanding of cancers, brain disorders, and other diseases. Filter-exchange imaging (FEXI), a special case of diffusion exchange spectroscopy adapted for clinical applications, has the potential to reveal different physiological water exchange processes. However, it is still controversial whether modulating the diffusion encoding gradient direction can affect the apparent exchange rate (AXR) measurements of FEXI in white matter (WM) where water diffusion shows strong anisotropy. In this study, we explored the diffusion-encoding direction dependence of FEXI in human brain white matter by performing FEXI with 20 diffusion-encoding directions on a clinical 3T scanner in-vivo. The results show that the AXR values measured when the gradients are perpendicular to the fiber orientation (0.77 ± 0.13 s ^- 1, mean ± standard deviation of all the subjects) are significantly larger than the AXR estimates when the gradients are parallel to the fiber orientation (0.33 ± 0.14 s ^- 1, p < 0.001) in WM voxels with coherently-orientated fibers. In addition, no significant correlation is found between AXRs measured along these two directions, indicating that they are measuring different water exchange processes. What's more, only the perpendicular AXR rather than the parallel AXR shows dependence on axonal diameter, indicating that the perpendicular AXR might reflect transmembrane water exchange between intra-axonal and extra-cellular spaces. Further finite difference (FD) simulations having three water compartments (intra-axonal, intra-glial, and extra-cellular spaces) to mimic WM micro-environments also suggest that the perpendicular AXR is more sensitive to the axonal water transmembrane exchange than parallel AXR. Taken together, our results show that AXR measured along different directions could be utilized to probe different water exchange processes in WM.

Axonal T2 estimation using the spherical variance of the strongly diffusion-weighted MRI signal

In magnetic resonance imaging, the application of a strong diffusion weighting suppresses the signal contributions from the less diffusion-restricted constituents of the brain's white matter, thus enabling the estimation of the transverse relaxation time T₂ that arises from the more diffusion-restricted constituents such as the axons. However, the presence of cell nuclei and vacuoles can confound the estimation of the axonal T₂, as diffusion within those structures is also restricted, causing the corresponding signal to survive the strong diffusion weighting. We devise an estimator of the axonal T₂ based on the directional spherical variance of the strongly diffusion-weighted signal. The spherical variance T₂ estimates are insensitive to the presence of isotropic contributions to the signal like those provided by cell nuclei and vacuoles. We show that with a strong diffusion weighting these estimates differ from those obtained using the directional spherical mean of the signal which contains both axonal and isotropically-restricted contributions. Our findings hint at the presence of an MRI-visible isotropically-restricted contribution to the signal in the white matter ex vivo fixed tissue (monkey) at 7T, and do not allow us to discard such a possibility also for in vivo human data collected with a clinical 3T system.

MASiVar: Multisite, multiscanner, and multisubject acquisitions for studying variability in diffusion weighted MRI

PURPOSE

Diffusion-weighted imaging allows investigators to identify structural, microstructural, and connectivity-based differences between subjects, but variability due to session and scanner biases is a challenge.

METHODS

To investigate DWI variability, we present MASiVar, a multisite data set consisting of 319 diffusion scans acquired at 3 T from b = 1000 to 3000 s/mm2 across 14 healthy adults, 83 healthy children (5 to 8 years), three sites, and four scanners as a publicly available, preprocessed, and de-identified data set. With the adult data, we demonstrate the capacity of MASiVar to simultaneously quantify the intrasession, intersession, interscanner, and intersubject variability of four common DWI processing approaches: (1) a tensor signal representation, (2) a multi-compartment neurite orientation dispersion and density model, (3) white-matter bundle segmentation, and (4) structural connectomics. Respectively, we evaluate region-wise fractional anisotropy, mean diffusivity, and principal eigenvector; region-wise CSF volume fraction, intracellular volume fraction, and orientation dispersion index; bundle-wise shape, volume, fractional anisotropy, and length; and whole connectome correlation and maximized modularity, global efficiency, and characteristic path length.

RESULTS

We plot the variability in these measures at each level and find that it consistently increases with intrasession to intersession to interscanner to intersubject effects across all processing approaches and that sometimes interscanner variability can approach intersubject variability.

CONCLUSIONS

This study demonstrates the potential of MASiVar to more globally investigate DWI variability across multiple levels and processing approaches simultaneously and suggests harmonization between scanners for multisite analyses should be considered before inference of group differences on subjects.

Microstructural Tissue Changes in a Rat Model of Mild Traumatic Brain Injury

Our study investigates the potential of diffusion MRI (dMRI), including diffusion tensor imaging (DTI), fixel-based analysis (FBA) and neurite orientation dispersion and density imaging (NODDI), to detect microstructural tissue abnormalities in rats after mild traumatic brain injury (mTBI). The brains of sham-operated and mTBI rats 35 days after lateral fluid percussion injury were imaged ex vivo in a 11.7-T scanner. Voxel-based analyses of DTI-, fixel- and NODDI-based metrics detected extensive tissue changes in directly affected brain areas close to the primary injury, and more importantly, also in distal areas connected to primary injury and indirectly affected by the secondary injury mechanisms. Histology revealed ongoing axonal abnormalities and inflammation, 35 days after the injury, in the brain areas highlighted in the group analyses. Fractional anisotropy (FA), fiber density (FD) and fiber density and fiber bundle cross-section (FDC) showed similar pattern of significant areas throughout the brain; however, FA showed more significant voxels in gray matter areas, while FD and FDC in white matter areas, and orientation dispersion index (ODI) in areas most damage based on histology. Region-of-interest (ROI)-based analyses on dMRI maps and histology in selected brain regions revealed that the changes in MRI parameters could be attributed to both alterations in myelinated fiber bundles and increased cellularity. This study demonstrates that the combination of dMRI methods can provide a more complete insight into the microstructural alterations in white and gray matter after mTBI, which may aid diagnosis and prognosis following a mild brain injury.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Axon diameter inferences in the human corpus callosum using oscillating gradient spin echo sequences

Previous methods used to infer axon diameter distributions using magnetic resonance imaging (MRI) primarily use single diffusion encoding sequences such as pulsed gradient spin echo (PGSE) and are thus sensitive to axons of diameters >5 μm. We applied oscillating gradient spin echo (OGSE) sequences to study human axons in the 1-2 μm range in the corpus callosum, which include the majority of axons constituting cortical connections. The ActiveAx model was applied to calculate the fitted mean effective diameter for axons (AxD) and was compared with values found using histology. Axon diameters from histological data were calculated using three different datasets; true diameters (minimum diameter), a combination of minimum and maximum diameters, and diameters measured across a consistent diffusion direction. The AxD estimates from MRI were 1.8 ± 0.1 μm to 2.34 ± 0.04 μm with an average of 2.0 ± 0.2 μm for the ActiveAx model. The histology AxD values were 1.43 ± 0.02 μm when using the true minimum axon diameters, 5.52 ± 0.02 μm when using the combination of minimum and maximum axon diameters, and 2.20 ± 0.02 μm when collecting measurements across a consistent diffusion direction. This experiment demonstrates the first known usage of OGSE to calculate axon diameters in the human corpus callosum on a 1-2 μm scale. The importance for the model to account for axonal orientation dispersion is indicated by histological results which more closely match the MRI model results depending on the direction of axon diameter measurements. These initial steps using this non-invasive imaging method can be applied to future methodology to develop in vivo axon diameter measurements in human brain tissue.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

SPHERIOUSLY? The challenges of estimating sphere radius non-invasively in the human brain from diffusion MRI

The Soma and Neurite Density Imaging (SANDI) three-compartment model was recently proposed to disentangle cylindrical and spherical geometries, attributed to neurite and soma compartments, respectively, in brain tissue. There are some recent advances in diffusion-weighted MRI signal encoding and analysis (including the use of multiple so-called 'b-tensor' encodings and analysing the signal in the frequency-domain) that have not yet been applied in the context of SANDI. In this work, using: (i) ultra-strong gradients; (ii) a combination of linear, planar, and spherical b-tensor encodings; and (iii) analysing the signal in the frequency domain, three main challenges to robust estimation of sphere size were identified: First, the Rician noise floor in magnitude-reconstructed data biases estimates of sphere properties in a non-uniform fashion. It may cause overestimation or underestimation of the spherical compartment size and density. This can be partly ameliorated by accounting for the noise floor in the estimation routine. Second, even when using the strongest diffusion-encoding gradient strengths available for human MRI, there is an empirical lower bound on the spherical signal fraction and radius that can be detected and estimated robustly. For the experimental setup used here, the lower bound on the sphere signal fraction was approximately 10%. We employed two different ways of establishing the lower bound for spherical radius estimates in white matter. The first, examining power-law relationships between the DW-signal and diffusion weighting in empirical data, yielded a lower bound of 7μm, while the second, pure Monte Carlo simulations, yielded a lower limit of 3μm and in this low radii domain, there is little differentiation in signal attenuation. Third, if there is sensitivity to the transverse intra-cellular diffusivity in cylindrical structures, e.g., axons and cellular projections, then trying to disentangle two diffusion-time-dependencies using one experimental parameter (i.e., change in frequency-content of the encoding waveform) makes spherical radii estimates particularly challenging. We conclude that due to the aforementioned challenges spherical radii estimates may be biased when the corresponding sphere signal fraction is low, which must be considered.

Impact of the acquisition protocol on the sensitivity to demyelination and axonal loss of clinically feasible DWI techniques: a simulation study

OBJECTIVE

To evaluate: (a) the specific effect that the demyelination and axonal loss have on the DW signal, and (b) the impact of the sequence parameters on the sensitivity to damage of two clinically feasible DWI techniques, i.e. DKI and NODDI.

METHODS

We performed a Monte Carlo simulation of water diffusion inside a novel synthetic model of white matter in the presence of axonal loss and demyelination, with three compartments with permeable boundaries between them. We compared DKI and NODDI in their ability to detect and assess the damage, using several acquisition protocols. We used the F test statistic as an index of the sensitivity for each DWI parameter to axonal loss and demyelination, respectively.

RESULTS

DKI parameters significantly changed with increasing axonal loss, but, in most cases, not with demyelination; all the NODDI parameters showed sensitivity to both the damage processes (at p < 0.01). However, the acquisition protocol strongly affected the sensitivity to damage of both the DKI and NODDI parameters and, especially for NODDI, the parameter absolute values also.

DISCUSSION

This work is expected to impact future choices for investigating white matter microstructure in focusing on specific stages of the disease, and for selecting the appropriate experimental framework to obtain optimal data quality given the purpose of the experiment.

On the use of multicompartment models of diffusion and relaxation for placental imaging

Multi-compartment models of diffusion and relaxation are ubiquitous in magnetic resonance research especially applied to neuroimaging applications. These models are increasingly making their way into the world of placental imaging. This review provides a framework for their motivation and implementation and describes some of the outstanding questions that need to be answered before they can be routinely adopted.

Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study

This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI.

Diffusion in Sephadex Gel Structures: Time Dependency Revealed by Multi-Sequence Acquisition over a Broad Diffusion Time Range

It has been increasingly reported that in biological tissues diffusion-weighted MRI signal attenuation deviates from mono-exponential decay, especially at high b-values. A number of diffusion models have been proposed to characterize this non-Gaussian diffusion behavior. One of these models is the continuous-time random-walk (CTRW) model, which introduces two new parameters: a fractional order time derivative α and a fractional order spatial derivative β. These new parameters have been linked to intravoxel diffusion heterogeneities in time and space, respectively, and are believed to depend on diffusion times. Studies on this time dependency are limited, largely because the diffusion time cannot vary over a board range in a conventional spin-echo echo-planar imaging sequence due to the accompanying T2 decays. In this study, we investigated the time-dependency of the CTRW model in Sephadex gel phantoms across a broad diffusion time range by employing oscillating-gradient spin-echo, pulsed-gradient spin-echo, and pulsed-gradient stimulated echo sequences. We also performed Monte Carlo simulations to help understand our experimental results. It was observed that the diffusion process fell into the Gaussian regime at extremely short diffusion times whereas it exhibited a strong time dependency in the CTRW parameters at longer diffusion times.

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

•

B-tensor encoding enables estimation of spherical cellular structures in the brain.

•

Spherical compartments may provide markers for apparent neural soma density.

•

Model parameters can be estimated in a fast and robust way using deep learning.

•

Practical acquisition times are achievable on widely available clinical scanners.

Diffusion MRI is a valuable tool for probing tissue microstructure in the brain noninvasively. Today, model-based techniques are widely available and used for white matter characterisation where their development is relatively mature. Conversely, tissue modelling in grey matter is more challenging, and no generally accepted models exist. With advances in measurement technology and modelling efforts, a clinically viable technique that reveals salient features of grey matter microstructure, such as the density of quasi-spherical cell bodies and quasi-cylindrical cell projections, is an exciting prospect. As a step towards capturing the microscopic architecture of grey matter in clinically feasible settings, this work uses a biophysical model that is designed to disentangle the diffusion signatures of spherical and cylindrical structures in the presence of orientation heterogeneity, and takes advantage of B-tensor encoding measurements, which provide additional sensitivity compared to standard single diffusion encoding sequences. For the fast and robust estimation of microstructural parameters, we leverage recent advances in machine learning and replace conventional fitting techniques with an artificial neural network that fits complex biophysical models within seconds. Our results demonstrate apparent markers of spherical and cylindrical geometries in healthy human subjects, and in particular an increased volume fraction of spherical compartments in grey matter compared to white matter. We evaluate the extent to which spherical and cylindrical geometries may be interpreted as correlates of neural soma and neural projections, respectively, and quantify parameter estimation errors in the presence of various departures from the modelling assumptions. While further work is necessary to translate the ideas presented in this work to the clinic, we suggest that biomarkers focussing on quasi-spherical cellular geometries may be valuable for the enhanced assessment of neurodevelopmental disorders and neurodegenerative diseases.