Neuroprosthetic baroreflex controls haemodynamics after spinal cord injury.

Spinal cord injury (SCI) induces haemodynamic instability that threatens survival1-3, impairs neurological recovery4,5, increases the risk of cardiovascular disease6,7, and reduces quality of life8,9. Haemodynamic instability in this context is due to the interruption of supraspinal efferent commands to sympathetic circuits located in the spinal cord10, which prevents the natural baroreflex from controlling these circuits to adjust peripheral vascular resistance. Epidural electrical stimulation (EES) of the spinal cord has been shown to compensate for interrupted supraspinal commands to motor circuits below the injury11, and restored walking after paralysis12. Here, we leveraged these concepts to develop EES protocols that restored haemodynamic stability after SCI. We established a preclinical model that enabled us to dissect the topology and dynamics of the sympathetic circuits, and to understand how EES can engage these circuits. We incorporated these spatial and temporal features into stimulation protocols to conceive a clinical-grade biomimetic haemodynamic regulator that operates in a closed loop. This 'neuroprosthetic baroreflex' controlled haemodynamics for extended periods of time in rodents, non-human primates and humans, after both acute and chronic SCI. We will now conduct clinical trials to turn the neuroprosthetic baroreflex into a commonly available therapy for people with SCI.

Diffusion-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling.

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

Parameter estimation for WMTI-Watson model of white matter using encoder-decoder recurrent neural network.

PURPOSE: Biophysical modeling of the diffusion MRI (dMRI) signal provides estimates of specific microstructural tissue properties. Although non-linear least squares (NLLS) is the most widespread fitting method, it suffers from local minima and high computational cost. Deep learning approaches are steadily replacing NLLS, but come with the limitation that the model needs to be retrained for each acquisition protocol and noise level. In this study, a novel fitting approach was proposed based on the encoder-decoder recurrent neural network (RNN) to accelerate model estimation with good generalization to various datasets. METHODS: The white matter tract integrity (WMTI)-Watson model as an implementation of the Standard Model of diffusion in white matter derives its parameters indirectly from the diffusion and kurtosis tensors (DKI). The RNN-based solver, which estimates the WMTI-Watson model from DKI, is therefore more readily translatable to various data, irrespective of acquisition protocols as long as the DKI was pre-computed from the signal. An embedding approach was also used to render the model insensitive to potential differences in distributions between training data and experimental data. The analytical solution, NLLS, RNN-, and a multilayer perceptron (MLP)-based methods were evaluated on synthetic and in vivo datasets of rat and human brain. RESULTS: The proposed RNN solver showed highly reduced computation time over the analytical solution and NLLS, with similar accuracy but improved robustness, and superior generalizability over MLP. CONCLUSION: The RNN estimator can be easily applied to various datasets without retraining, which shows great potential for a widespread use.

Cellular Exchange Imaging (CEXI): Evaluation of a diffusion model including water exchange in cells using numerical phantoms of permeable spheres.

PURPOSE: Biophysical models of diffusion MRI have been developed to characterize microstructure in various tissues, but existing models are not suitable for tissue composed of permeable spherical cells. In this study we introduce Cellular Exchange Imaging (CEXI), a model tailored for permeable spherical cells, and compares its performance to a related Ball & Sphere (BS) model that neglects permeability. METHODS: We generated DW-MRI signals using Monte-Carlo simulations with a PGSE sequence in numerical substrates made of spherical cells and their extracellular space for a range of membrane permeability. From these signals, the properties of the substrates were inferred using both BS and CEXI models. RESULTS: CEXI outperformed the impermeable model by providing more stable estimates cell size and intracellular volume fraction that were diffusion time-independent. Notably, CEXI accurately estimated the exchange time for low to moderate permeability levels previously reported in other studies ( κ < 25 µ m / s $$ \kappa <25\kern0.3em \mu \mathrm{m}/\mathrm{s} $$ ). However, in highly permeable substrates ( κ = 50 µ m / s $$ \kappa =50\kern0.3em \mu \mathrm{m}/\mathrm{s} $$ ), the estimated parameters were less stable, particularly the diffusion coefficients. CONCLUSION: This study highlights the importance of modeling the exchange time to accurately quantify microstructure properties in permeable cellular substrates. Future studies should evaluate CEXI in clinical applications such as lymph nodes, investigate exchange time as a potential biomarker of tumor severity, and develop more appropriate tissue models that account for anisotropic diffusion and highly permeable membranes.

MP-PCA denoising for diffusion MRS data: promises and pitfalls.

Diffusion-weighted (DW) magnetic resonance spectroscopy (MRS) suffers from a lower signal to noise ratio (SNR) compared to conventional MRS owing to the addition of diffusion attenuation. This technique can therefore strongly benefit from noise reduction strategies. In the present work, Marchenko-Pastur principal component analysis (MP-PCA) denoising is tested on Monte Carlo simulations and on in vivo DW-MRS data acquired at 9.4 T in rat brain and at 3 T in human brain. We provide a descriptive study of the effects observed following different MP-PCA denoising strategies (denoising the entire matrix versus using a sliding window), in terms of apparent SNR, rank selection, noise correlation within and across b-values and quantification of metabolite concentrations and fitted diffusion coefficients. MP-PCA denoising yielded an increased apparent SNR, a more accurate B0 drift correction between shots, and similar estimates of metabolite concentrations and diffusivities compared to the raw data. No spectral residuals on individual shots were observed but correlations in the noise level across shells were introduced, an effect which was mitigated using a sliding window, but which should be carefully considered.

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.

Synchronous nonmonotonic changes in functional connectivity and white matter integrity in a rat model of sporadic Alzheimer's disease.

Brain glucose hypometabolism has been singled out as an important contributor and possibly main trigger to Alzheimer's disease (AD). Intracerebroventricular injections of streptozotocin (icv-STZ) cause brain glucose hypometabolism without systemic diabetes. Here, a first-time longitudinal study of brain glucose metabolism, functional connectivity and white matter microstructure was performed in icv-STZ rats using PET and MRI. Histological markers of pathology were tested at an advanced stage of disease. STZ rats exhibited altered functional connectivity and intra-axonal damage and demyelination in brain regions typical of AD, in a temporal pattern of acute injury, transient recovery/compensation and chronic degeneration. In the context of sustained glucose hypometabolism, these nonmonotonic trends - also reported in behavioral studies of this animal model as well as in human AD - suggest a compensatory mechanism, possibly recruiting ketone bodies, that allows a partial and temporary repair of brain structure and function. The early acute phase could thus become a valuable therapeutic window to strengthen the recovery phase and prevent or delay chronic degeneration, to be considered both in preclinical and clinical studies of AD. In conclusion, this work reveals the consequences of brain insulin resistance on structure and function, highlights signature nonmonotonic trajectories in their evolution and proposes potent MRI-derived biomarkers translatable to human AD and diabetic populations.

Intra- and extra-axonal axial diffusivities in the white matter: Which one is faster?

A two-compartment model of diffusion in white matter, which accounts for intra- and extra-axonal spaces, is associated with two plausible mathematical scenarios: either the intra-axonal axial diffusivity Da,‖ is higher than the extra-axonal De,‖ (Branch 1), or the opposite, i.e. Da,‖ < De,‖ (Branch 2). This duality calls for an independent validation of compartment axial diffusivities, to determine which of the two cases holds. The aim of the present study was to use an intracerebroventricular injection of a gadolinium-based contrast agent to selectively reduce the extracellular water signal in the rat brain, and compare diffusion metrics in the genu of the corpus callosum before and after gadolinium infusion. The diffusion metrics considered were diffusion and kurtosis tensor metrics, as well as compartment-specific estimates of the WMTI-Watson two-compartment model. A strong decrease in genu T1 and T2 relaxation times post-Gd was observed (p < 0.001), as well as an increase of 48% in radial kurtosis (p < 0.05), which implies that the relative fraction of extracellular water signal was selectively decreased. This was further supported by a significant increase in intra-axonal water fraction as estimated from the two-compartment model, for both branches (p < 0.01 for Branch 1, p < 0.05 for Branch 2). However, pre-Gd estimates of axon dispersion in Branch 1 agreed better with literature than those of Branch 2. Furthermore, comparison of post-Gd changes in diffusivity and dispersion between data and simulations further supported Branch 1 as the biologically plausible solution, i.e. Da,‖ > De,‖. This result is fully consistent with other recent measurements of compartment axial diffusivities that used entirely different approaches, such as diffusion tensor encoding.

Rotationally-invariant mapping of scalar and orientational metrics of neuronal microstructure with diffusion MRI.

We develop a general analytical and numerical framework for estimating intra- and extra-neurite water fractions and diffusion coefficients, as well as neurite orientational dispersion, in each imaging voxel. By employing a set of rotational invariants and their expansion in the powers of diffusion weighting, we analytically uncover the nontrivial topology of the parameter estimation landscape, showing that multiple branches of parameters describe the measurement almost equally well, with only one of them corresponding to the biophysical reality. A comprehensive acquisition shows that the branch choice varies across the brain. Our framework reveals hidden degeneracies in MRI parameter estimation for neuronal tissue, provides microstructural and orientational maps in the whole brain without constraints or priors, and connects modern biophysical modeling with clinical MRI.

Functional magnetic resonance microscopy at single-cell resolution in Aplysia californica.

In this work, we show the feasibility of performing functional MRI studies with single-cell resolution. At ultrahigh magnetic field, manganese-enhanced magnetic resonance microscopy allows the identification of most motor neurons in the buccal network of Aplysia at low, nontoxic Mn(2+) concentrations. We establish that Mn(2+) accumulates intracellularly on injection into the living Aplysia and that its concentration increases when the animals are presented with a sensory stimulus. We also show that we can distinguish between neuronal activities elicited by different types of stimuli. This method opens up a new avenue into probing the functional organization and plasticity of neuronal networks involved in goal-directed behaviors with single-cell resolution.

In vivo quantification of demyelination and recovery using compartment-specific diffusion MRI metrics validated by electron microscopy.

There is a need for accurate quantitative non-invasive biomarkers to monitor myelin pathology in vivo and distinguish myelin changes from other pathological features including inflammation and axonal loss. Conventional MRI metrics such as T2, magnetization transfer ratio and radial diffusivity have proven sensitivity but not specificity. In highly coherent white matter bundles, compartment-specific white matter tract integrity (WMTI) metrics can be directly derived from the diffusion and kurtosis tensors: axonal water fraction, intra-axonal diffusivity, and extra-axonal radial and axial diffusivities. We evaluate the potential of WMTI to quantify demyelination by monitoring the effects of both acute (6weeks) and chronic (12weeks) cuprizone intoxication and subsequent recovery in the mouse corpus callosum, and compare its performance with that of conventional metrics (T2, magnetization transfer, and DTI parameters). The changes observed in vivo correlated with those obtained from quantitative electron microscopy image analysis. A 6-week intoxication produced a significant decrease in axonal water fraction (p<0.001), with only mild changes in extra-axonal radial diffusivity, consistent with patchy demyelination, while a 12-week intoxication caused a more marked decrease in extra-axonal radial diffusivity (p=0.0135), consistent with more severe demyelination and clearance of the extra-axonal space. Results thus revealed increased specificity of the axonal water fraction and extra-axonal radial diffusivity parameters to different degrees and patterns of demyelination. The specificities of these parameters were corroborated by their respective correlations with microstructural features: the axonal water fraction correlated significantly with the electron microscopy derived total axonal water fraction (ρ=0.66; p=0.0014) but not with the g-ratio, while the extra-axonal radial diffusivity correlated with the g-ratio (ρ=0.48; p=0.0342) but not with the electron microscopy derived axonal water fraction. These parameters represent promising candidates as clinically feasible biomarkers of demyelination and remyelination in the white matter.

Gibbs ringing in diffusion MRI.

PURPOSE: To study and reduce the effect of Gibbs ringing artifact on computed diffusion parameters. METHODS: We reduce the ringing by extrapolating the k-space of each diffusion weighted image beyond the measured part by selecting an adequate regularization term. We evaluate several regularization terms and tune the regularization parameter to find the best compromise between anatomical accuracy of the reconstructed image and suppression of the Gibbs artifact. RESULTS: We demonstrate empirically and analytically that the Gibbs artifact, which is typically observed near sharp edges in magnetic resonance images, has a significant impact on the quantification of diffusion model parameters, even for infinitesimal diffusion weighting. We find the second order total generalized variation to be a good choice for the penalty term to regularize the extrapolation of the k-space, as it provides a parsimonious representation of images, a practically full suppression of Gibbs ringing, and the absence of staircasing artifacts typical for total variation methods. CONCLUSIONS: Regularized extrapolation of the k-space data significantly reduces truncation artifacts without compromising spatial resolution in comparison to the default option of window filtering. In particular, accuracy of estimating diffusion tensor imaging and diffusion kurtosis imaging parameters improves so much that unconstrained fits become possible. Magn Reson Med 76:301-314, 2016. © 2015 Wiley Periodicals, Inc.

Degeneracy in model parameter estimation for multi-compartmental diffusion in neuronal tissue.

The ultimate promise of diffusion MRI (dMRI) models is specificity to neuronal microstructure, which may lead to distinct clinical biomarkers using noninvasive imaging. While multi-compartment models are a common approach to interpret water diffusion in the brain in vivo, the estimation of their parameters from the dMRI signal remains an unresolved problem. Practically, even when q space is highly oversampled, nonlinear fit outputs suffer from heavy bias and poor precision. So far, this has been alleviated by fixing some of the model parameters to a priori values, for improved precision at the expense of accuracy. Here we use a representative two-compartment model to show that fitting fails to determine the five model parameters from over 60 measurement points. For the first time, we identify the reasons for this poor performance. The first reason is the existence of two local minima in the parameter space for the objective function of the fitting procedure. These minima correspond to qualitatively different sets of parameters, yet they both lie within biophysically plausible ranges. We show that, at realistic signal-to-noise ratio values, choosing between the two minima based on the associated objective function values is essentially impossible. Second, there is an ensemble of very low objective function values around each of these minima in the form of a pipe. The existence of such a direction in parameter space, along which the objective function profile is very flat, explains the bias and large uncertainty in parameter estimation, and the spurious parameter correlations: in the presence of noise, the minimum can be randomly displaced by a very large amount along each pipe. Our results suggest that the biophysical interpretation of dMRI model parameters crucially depends on establishing which of the minima is closer to the biophysical reality and the size of the uncertainty associated with each parameter.

One diffusion acquisition and different white matter models: how does microstructure change in human early development based on WMTI and NODDI?

White matter microstructural changes during the first three years of healthy brain development are characterized using two different models developed for limited clinical diffusion data: White Matter Tract Integrity (WMTI) metrics from Diffusional Kurtosis Imaging (DKI) and Neurite Orientation Dispersion and Density Imaging (NODDI). Both models reveal a non-linear increase in intra-axonal water fraction and in tortuosity of the extra-axonal space as a function of age, in the genu and splenium of the corpus callosum and the posterior limb of the internal capsule. The changes are consistent with expected behavior related to myelination and asynchrony of fiber development. The intra- and extracellular axial diffusivities as estimated with WMTI do not change appreciably in normal brain development. The quantitative differences in parameter estimates between models are examined and explained in the light of each model's assumptions and consequent biases, as highlighted in simulations. Finally, we discuss the feasibility of a model with fewer assumptions.

Effects of hypotonic stress and ouabain on the apparent diffusion coefficient of water at cellular and tissue levels in Aplysia.

There is evidence that physiological or pathological cell swelling is associated with a decrease of the apparent diffusion coefficient (ADC) of water in tissues, as measured with MRI. However the mechanism remains unclear. Magnetic resonance microscopy, performed on small tissue samples, has the potential to distinguish effects occurring at cellular and tissue levels. A three-dimensional diffusion prepared fast imaging with steady-state free precession sequence for MR microscopy was implemented on a 17.2 T imaging system and used to investigate the effect of two biological challenges known to cause cell swelling, exposure to a hypotonic solution or to ouabain, on Aplysia nervous tissue. The ADC was measured inside isolated neuronal soma and in the region of cell bodies of the buccal ganglia. Both challenges resulted in an ADC increase inside isolated neuronal soma (+31 ± 24% and +30 ± 11%, respectively) and an ADC decrease at tissue level in the buccal ganglia (-12 ± 5% and -18 ± 8%, respectively). A scenario involving a layer of water molecules bound to the inflating cell membrane surface is proposed to reconcile this apparent discrepancy.

Diffusion of brain metabolites highlights altered brain microstructure in type C hepatic encephalopathy: a 9.4 T preliminary study.

Introduction: Type C hepatic encephalopathy (HE) is a decompensating event of chronic liver disease leading to severe motor and cognitive impairment. The progression of type C HE is associated with changes in brain metabolite concentrations measured by 1H magnetic resonance spectroscopy (MRS), most noticeably a strong increase in glutamine to detoxify brain ammonia. In addition, alterations of brain cellular architecture have been measured ex vivo by histology in a rat model of type C HE. The aim of this study was to assess the potential of diffusion-weighted MRS (dMRS) for probing these cellular shape alterations in vivo by monitoring the diffusion properties of the major brain metabolites. Methods: The bile duct-ligated (BDL) rat model of type C HE was used. Five animals were scanned before surgery and 6- to 7-week post-BDL surgery, with each animal being used as its own control. 1H-MRS was performed in the hippocampus (SPECIAL, TE = 2.8 ms) and dMRS in a voxel encompassing the entire brain (DW-STEAM, TE = 15 ms, diffusion time = 120 ms, maximum b-value = 25 ms/µm2) on a 9.4 T scanner. The in vivo MRS acquisitions were further validated with histological measures (immunohistochemistry, Golgi-Cox, electron microscopy). Results: The characteristic 1H-MRS pattern of type C HE, i.e., a gradual increase of brain glutamine and a decrease of the main organic osmolytes, was observed in the hippocampus of BDL rats. Overall increased metabolite diffusivities (apparent diffusion coefficient and intra-stick diffusivity-Callaghan's model, significant for glutamine, myo-inositol, and taurine) and decreased kurtosis coefficients were observed in BDL rats compared to control, highlighting the presence of osmotic stress and possibly of astrocytic and neuronal alterations. These results were consistent with the microstructure depicted by histology and represented by a decline in dendritic spines density in neurons, a shortening and decreased number of astrocytic processes, and extracellular edema. Discussion: dMRS enables non-invasive and longitudinal monitoring of the diffusion behavior of brain metabolites, reflecting in the present study the globally altered brain microstructure in BDL rats, as confirmed ex vivo by histology. These findings give new insights into metabolic and microstructural abnormalities associated with high brain glutamine and its consequences in type C HE.

Microstructural alterations in the locus coeruleus-entorhinal cortex pathway in Alzheimer's disease and frontotemporal dementia.

INTRODUCTION: We investigated in vivo the microstructural integrity of the pathway connecting the locus coeruleus to the transentorhinal cortex (LC-TEC) in patients with Alzheimer's disease (AD) and frontotemporal dementia (FTD). METHODS: Diffusion-weighted MRI scans were collected for 21 AD, 20 behavioral variants of FTD (bvFTD), and 20 controls. Fractional anisotropy (FA), mean, axial, and radial diffusivities (MD, AxD, RD) were computed in the LC-TEC pathway using a normative atlas. Atrophy was assessed using cortical thickness and correlated with microstructural measures. RESULTS: We found (i) higher RD in AD than controls; (ii) higher MD, RD, and AxD, and lower FA in bvFTD than controls and AD; and (iii) a negative association between LC-TEC MD, RD, and AxD, and entorhinal cortex (EC) thickness in bvFTD (all p < 0.050). DISCUSSION: LC-TEC microstructural alterations are more pronounced in bvFTD than AD, possibly reflecting neurodegeneration secondary to EC atrophy. Highlights: Microstructural integrity of LC-TEC pathway is understudied in AD and bvFTD.LC-TEC microstructural alterations are present in both AD and bvFTD.Greater LC-TEC microstructural alterations in bvFTD than AD.LC-TEC microstructural alterations in bvFTD are associated to EC neurodegeneration.

Highlighting manganese dynamics in the nervous system of Aplysia californica using MEMRI at ultra-high field.

Exploring the pathways of manganese (Mn(2+)) transport in the nervous system becomes of interest as many recent studies use Mn(2+) as a neural tract tracer in mammals. In this study, we performed manganese enhanced MRI (MEMRI) at 17.2 T on the buccal ganglia of Aplysia californica. The main advantage of this model over mammalian systems is that it contains networks of large identified neurons. Using Mn(2+) retrograde transport along selected nerves, we first validated the mapping of motor neurons' axonal projections into peripheral nerves, previously obtained from optical imaging (Morton et al., 1991). This protocol was found not to alter the functional properties of the neuronal network. Second, we compared the Mn(2+) dynamics inside the ganglia in the presence or absence of chemical stimulation. We found that 2h of stimulation with the modulatory transmitter dopamine increased the extent of areas of intermediate signal enhancement caused by manganese accumulation. In the absence of dopamine, an overall decrease of the enhanced areas in favor of non-enhanced areas was found, as a result of natural Mn(2+) washout. This supports the hypothesis that, upon activation, Mn(2+) is released from labeled neurons and captured by other, initially unlabeled, neurons. However, the latter could not be clearly identified due to lack of sensitivity and multiplicity of possible pathways starting from labeled cells. Nonetheless, the Aplysia buccal ganglia remain a well-suited model for attempting to visualize Mn(2+) transport from neuron to neuron upon activation, as well as for studying dopaminergic modulation in a motor network.

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

BACKGROUND: The pathological process of Alzheimer's disease (AD) typically takes decades from onset to clinical symptoms. Early brain changes in AD include MRI-measurable features such as altered functional connectivity (FC) and white matter degeneration. The ability of these features to discriminate between subjects without a diagnosis, or their prognostic value, is however not established. METHODS: The main trigger mechanism of AD is still debated, although impaired brain glucose metabolism is taking an increasingly central role. Here, we used a rat model of sporadic AD, based on impaired brain glucose metabolism induced by an intracerebroventricular injection of streptozotocin (STZ). We characterized alterations in FC and white matter microstructure longitudinally using functional and diffusion MRI. Those MRI-derived measures were used to classify STZ from control rats using machine learning, and the importance of each individual measure was quantified using explainable artificial intelligence methods. RESULTS: Overall, combining all the FC and white matter metrics in an ensemble way was the best strategy to discriminate STZ rats, with a consistent accuracy over 0.85. However, the best accuracy early on was achieved using white matter microstructure features, and later on using FC. This suggests that consistent damage in white matter in the STZ group might precede FC. For cross-timepoint prediction, microstructure features also had the highest performance while, in contrast, that of FC was reduced by its dynamic pattern which shifted from early hyperconnectivity to late hypoconnectivity. CONCLUSIONS: Our study highlights the MRI-derived measures that best discriminate STZ vs control rats early in the course of the disease, with potential translation to humans.