Ensemble average propagator-based detection of microstructural alterations after stroke

Feature attention graph neural network for estimating brain age and identifying important neural connections in mouse models of genetic risk for Alzheimer's disease

Alzheimer's disease (AD) remains one of the most extensively researched neurodegenerative disorders due to its widespread prevalence and complex risk factors. Age is a crucial risk factor for AD, which can be estimated by the disparity between physiological age and estimated brain age. To model AD risk more effectively, integrating biological, genetic, and cognitive markers is essential. Here, we utilized mouse models expressing the major APOE human alleles and human nitric oxide synthase 2 to replicate genetic risk for AD and a humanized innate immune response. We estimated brain age employing a multivariate dataset that includes brain connectomes, APOE genotype, subject traits such as age and sex, and behavioral data. Our methodology used Feature Attention Graph Neural Networks (FAGNN) for integrating different data types. Behavioral data were processed with a 2D Convolutional Neural Network (CNN), subject traits with a 1D CNN, brain connectomes through a Graph Neural Network using quadrant attention module. The model yielded a mean absolute error for age prediction of 31.85 days, with a root mean squared error of 41.84 days, outperforming other, reduced models. In addition, FAGNN identified key brain connections involved in the aging process. The highest weights were assigned to the connections between cingulum and corpus callosum, striatum, hippocampus, thalamus, hypothalamus, cerebellum, and piriform cortex. Our study demonstrates the feasibility of predicting brain age in models of aging and genetic risk for AD. To verify the validity of our findings, we compared Fractional Anisotropy (FA) along the tracts of regions with the highest connectivity, the Return-to-Origin Probability (RTOP), Return-to-Plane Probability (RTPP), and Return-to-Axis Probability (RTAP), which showed significant differences between young, middle-aged, and old age groups. Younger mice exhibited higher FA, RTOP, RTAP, and RTPP compared to older groups in the selected connections, suggesting that degradation of white matter tracts plays a critical role in aging and for FAGNN's selections. Our analysis suggests a potential neuroprotective role of APOE2, relative to APOE3 and APOE4, where APOE2 appears to mitigate age-related changes. Our findings highlighted a complex interplay of genetics and brain aging in the context of AD risk modeling.

Correspondence of mean apparent propagator MRI metrics with phosphorylated tau and astrogliosis in chronic traumatic encephalopathy

Chronic traumatic encephalopathy is a neurodegenerative disease that is diagnosed and staged based on the localization and extent of phosphorylated tau pathology. Although its identification remains the primary diagnostic criteria to distinguish chronic traumatic encephalopathy from other tauopathies, the hyperphosphorylated tau that accumulates in neurofibrillary tangles in cortical grey matter and perivascular regions is often accompanied by concomitant pathology such as astrogliosis. Mean apparent propagator MRI is a clinically feasible diffusion MRI method that is suitable to characterize microstructure of complex biological media efficiently and comprehensively. We performed quantitative correlations between propagator metrics and underlying phosphorylated tau and astroglial pathology in a cross-sectional study of 10 ex vivo human tissue specimens with 'high chronic traumatic encephalopathy' at 0.25 mm isotropic voxels. Linear mixed effects analysis of regions of interest showed significant relationships of phosphorylated tau with propagator-estimated non-Gaussianity in cortical grey matter (P = 0.002) and of astrogliosis with propagator anisotropy in superficial cortical white matter (P = 0.0009). The positive correlation between phosphorylated tau and non-Gaussianity was found to be modest but significant (R2 = 0.44, P = 6.0 × 10-5) using linear regression. We developed an unsupervised clustering algorithm with non-Gaussianity and propagator anisotropy as inputs, which was able to identify voxels in superficial cortical white matter that corresponded to astrocytes that were accumulated at the grey-white matter interface. Our results suggest that mean apparent propagator MRI at high spatial resolution provides a means to not only identify phosphorylated tau pathology but also detect regions with astrocytic pathology and may therefore prove diagnostically valuable in the evaluation of concomitant pathology in cortical tissue with complex microstructure.

Viability of AMURA biomarkers from single-shell diffusion MRI in clinical studies

Diffusion Tensor Imaging (DTI) is the most employed method to assess white matter properties using quantitative parameters derived from diffusion MRI, but it presents known limitations that restrict the evaluation of complex structures. The objective of this study was to validate the reliability and robustness of complementary diffusion measures extracted with a novel approach, Apparent Measures Using Reduced Acquisitions (AMURA), with a typical diffusion MRI acquisition from a clinical context in comparison with DTI with application to clinical studies. Fifty healthy controls, 51 episodic migraine and 56 chronic migraine patients underwent single-shell diffusion MRI. Four DTI-based and eight AMURA-based parameters were compared between groups with tract-based spatial statistics to establish reference results. On the other hand, following a region-based analysis, the measures were assessed for multiple subsamples with diverse reduced sample sizes and their stability was evaluated with the coefficient of quartile variation. To assess the discrimination power of the diffusion measures, we repeated the statistical comparisons with a region-based analysis employing reduced sample sizes with diverse subsets, decreasing 10 subjects per group for consecutive reductions, and using 5,001 different random subsamples. For each sample size, the stability of the diffusion descriptors was evaluated with the coefficient of quartile variation. AMURA measures showed a greater number of statistically significant differences in the reference comparisons between episodic migraine patients and controls compared to DTI. In contrast, a higher number of differences was found with DTI parameters compared to AMURA in the comparisons between both migraine groups. Regarding the assessments reducing the sample size, the AMURA parameters showed a more stable behavior than DTI, showing a lower decrease for each reduced sample size or a higher number of regions with significant differences. However, most AMURA parameters showed lower stability in relation to higher coefficient of quartile variation values than the DTI descriptors, although two AMURA measures showed similar values to DTI. For the synthetic signals, there were AMURA measures with similar quantification to DTI, while other showed similar behavior. These findings suggest that AMURA presents favorable characteristics to identify differences of specific microstructural properties between clinical groups in regions with complex fiber architecture and lower dependency on the sample size or assessing technique than DTI.

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.

HYDI-DSI revisited: Constrained non-parametric EAP imaging without q-space re-gridding

Hybrid Diffusion Imaging (HYDI) was one of the first attempts to use multi-shell samplings of the q-space to infer diffusion properties beyond Diffusion Tensor Imaging (DTI) or High Angular Resolution Diffusion Imaging (HARDI). HYDI was intended as a flexible protocol embedding both DTI (for lower b-values) and HARDI (for higher b-values) processing, as well as Diffusion Spectrum Imaging (DSI) when the entire data set was exploited. In the latter case, the spherical sampling of the q-space is re-gridded by interpolation to a Cartesian lattice whose extent covers the range of acquired b-values, hence being acquisition-dependent. The Discrete Fourier Transform (DFT) is afterwards used to compute the corresponding Cartesian sampling of the Ensemble Average Propagator (EAP) in an entirely non-parametric way. From this lattice, diffusion markers such as the Return To Origin Probability (RTOP) or the Mean Squared Displacement (MSD) can be numerically estimated. We aim at re-formulating this scheme by means of a Fourier Transform encoding matrix that eliminates the need for q-space re-gridding at the same time it preserves the non-parametric nature of HYDI-DSI. The encoding matrix is adaptively designed at each voxel according to the underlying DTI approximation, so that an optimal sampling of the EAP can be pursued without being conditioned by the particular acquisition protocol. The estimation of the EAP is afterwards carried out as a regularized Quadratic Programming (QP) problem, which allows to impose positivity constraints that cannot be trivially embedded within the conventional HYDI-DSI. We demonstrate that the definition of the encoding matrix in the adaptive space allows to analytically (as opposed to numerically) compute several popular descriptors of diffusion with the unique source of error being the cropping of high frequency harmonics in the Fourier analysis of the attenuation signal. They include not only RTOP and MSD, but also Return to Axis/Plane Probabilities (RTAP/RTPP), which are defined in terms of specific spatial directions and are not available with the former HYDI-DSI. We report extensive experiments that suggest the benefits of our proposal in terms of accuracy, robustness and computational efficiency, especially when only standard, non-dedicated q-space samplings are available.

Diffusion-relaxation scattered MR signal representation in a multi-parametric sequence

This work focuses on obtaining a more general diffusion magnetic resonance imaging (MRI) signal representation that accounts for a longitudinal T₁ and transverse T₂^⋆ relaxations while at the same time integrating directional diffusion in the context of scattered multi-parametric acquisitions, where only a few diffusion gradient directions and b-values are available for each pair of echo and inversion times. The method is based on the three-dimensional simple harmonic oscillator-based reconstruction and estimation (SHORE) representation of the diffusion signal, which enables the estimation of the orientation distribution function and the retrieval of various quantitative indices such as the generalized fractional anisotropy or the return-to-the-origin probability while simultaneously resolving for T₁ and T₂^⋆ relaxation times. Our technique, the Relax-SHORE, has been tested on both in silico and in vivo diffusion-relaxation scattered MR data. The results show that Relax-SHORE is accurate in the context of scattered acquisitions while guaranteeing flexibility in the diffusion signal representation from multi-parametric sequences.

Beyond Diffusion Tensor MRI Methods for Improved Characterization of the Brain after Ischemic Stroke: A Review

Ischemic stroke is a worldwide problem, with 15 million people experiencing a stroke annually. MR imaging is a valuable tool for understanding and assessing brain changes after stroke and predicting recovery. Of particular interest is the use of diffusion MR imaging in the nonacute stage 1-30 days poststroke. Thousands of articles have been published on the use of diffusion MR imaging in stroke, including several recent articles reviewing the use of DTI for stroke. The goal of this work was to survey and put into context the recent use of diffusion MR imaging methods beyond DTI, including diffusional kurtosis, generalized fractional anisotropy, spherical harmonics methods, and neurite orientation and dispersion models, in patients poststroke. Early studies report that these types of beyond-DTI methods outperform DTI metrics either in being more sensitive to poststroke changes or by better predicting outcome motor scores. More and larger studies are needed to confirm the improved prediction of stroke recovery with the beyond-DTI methods.

White matter microstructural impairments in amyotrophic lateral sclerosis: A mean apparent propagator MRI study

Background

White matter (WM) impairment is a hallmark of amyotrophic lateral sclerosis (ALS). This study evaluated the capacity of mean apparent propagator magnetic resonance imaging (MAP-MRI) for detecting ALS-related WM alterations.

Methods

Diffusion images were obtained from 52 ALS patients and 51 controls. MAP-derived indices [return-to-origin/-axis/-plane probability (RTOP/RTAP/RTPP) and non-Gaussianity (NG)/perpendicular/parallel NG (NG_⊥/NG_||)] were computed. Measures from diffusion tensor/kurtosis imaging (DTI/DKI) and neurite orientation dispersion and density imaging (NODDI) were also obtained. Voxel-wise analysis (VBA) was performed to determine differences in these parameters. Relationship between MAP parameters and disease severity (assessed by the revised ALS Functional Rating Scale (ALSFRS-R)) was evaluated by Pearson’s correlation analysis in a voxel-wise way. ALS patients were further divided into two subgroups: 29 with limb-only involvement and 23 with both bulbar and limb involvement. Subgroup analysis was then conducted to investigate diffusion parameter differences related to bulbar impairment.

Results

The VBA (with threshold of P < 0.05 after family-wise error correction (FWE)) showed that ALS patients had significantly decreased RTOP/RTAP/RTPP and NG/ NG_⊥/NG_|| in a set of WM areas, including the bilateral precentral gyrus, corona radiata, posterior limb of internal capsule, midbrain, middle corpus callosum, anterior corpus callosum, parahippocampal gyrus, and medulla. MAP-MRI had the capacity to capture WM damage in ALS, which was higher than DTI and similar to DKI/NODDI. RTOP/RTAP/NG/NG_⊥/NG_|| parameters, especially in the bilateral posterior limb of internal capsule and middle corpus callosum, were significantly correlated with ALSFRS-R (with threshold of FWE-corrected P < 0.05). The VBA (with FWE-corrected P < 0.05) revealed the significant RTAP reduction in subgroup with both bulbar and limb involvement, compared with those with limb-only involvement.

Conclusions

Microstructural impairments in corticospinal tract and corpus callosum represent the consistent characteristic of ALS. MAP-MRI could provide alternative measures depicting ALS-related WM alterations, complementary to the common diffusion imaging methods.

Interpretable deep learning as a means for decrypting disease signature in multiple sclerosis

Objective.The mechanisms driving multiple sclerosis (MS) are still largely unknown, calling for new methods allowing to detect and characterize tissue degeneration since the early stages of the disease. Our aim is to decrypt the microstructural signatures of the Primary Progressive versus the Relapsing-Remitting state of disease based on diffusion and structural magnetic resonance imaging data.Approach.A selection of microstructural descriptors, based on the 3D-Simple Harmonics Oscillator Based Reconstruction and Estimation and the set of new algebraically independent Rotation Invariant spherical harmonics Features, was considered and used to feed convolutional neural networks (CNNs) models. Classical measures derived from diffusion tensor imaging, that are fractional anisotropy and mean diffusivity, were used as benchmark for diffusion MRI (dMRI). Finally, T1-weighted images were also considered for the sake of comparison with the state-of-the-art. A CNN model was fit to each feature map and layerwise relevance propagation (LRP) heatmaps were generated for each model, target class and subject in the test set. Average heatmaps were calculated across correctly classified patients and size-corrected metrics were derived on a set of regions of interest to assess the LRP contrast between the two classes.Main results.Our results demonstrated that dMRI features extracted in grey matter tissues can help in disambiguating primary progressive multiple sclerosis from relapsing-remitting multiple sclerosis patients and, moreover, that LRP heatmaps highlight areas of high relevance which relate well with what is known from literature for MS disease.Significance.Within a patient stratification task, LRP allows detecting the input voxels that mostly contribute to the classification of the patients in either of the two classes for each feature, potentially bringing to light hidden data properties which might reveal peculiar disease-state factors.

Unraveling the MRI-Based Microstructural Signatures Behind Primary Progressive and Relapsing-Remitting Multiple Sclerosis Phenotypes

Background

The mechanisms driving primary progressive and relapsing–remitting multiple sclerosis (PPMS/RRMS) phenotypes are unknown. Magnetic resonance imaging (MRI) studies support the involvement of gray matter (GM) in the degeneration, highlighting its damage as an early feature of both phenotypes. However, the role of GM microstructure is unclear, calling for new methods for its decryption.

Purpose

To investigate the morphometric and microstructural GM differences between PPMS and RRMS to characterize GM tissue degeneration using MRI.

Study Type

Prospective cross‐sectional study.

Subjects

Forty‐five PPMS (26 females) and 45 RRMS (32 females) patients.

Field Strength/Sequence

3T scanner. Three‐dimensional (3D) fast field echo T1‐weighted (T1‐w), 3D turbo spin echo (TSE) T2‐w, 3D TSE fluid‐attenuated inversion recovery, and spin echo‐echo planar imaging diffusion MRI (dMRI).

Assessment

T1‐w and dMRI data were employed for providing information about morphometric and microstructural features, respectively. For dMRI, both diffusion tensor imaging and 3D simple harmonics oscillator based reconstruction and estimation models were used for feature extraction from a predefined set of regions. A support vector machine (SVM) was used to perform patients' classification relying on all these measures.

Statistical Tests

Differences between MS phenotypes were investigated using the analysis of covariance and statistical tests (P < 0.05 was considered statistically significant).

Results

All the dMRI indices showed significant microstructural alterations between the considered MS phenotypes, for example, the mode and the median of the return to the plane probability in the hippocampus. Conversely, thalamic volume was the only morphometric feature significantly different between the two MS groups. Ten of the 12 features retained by the selection process as discriminative across the two MS groups regarded the hippocampus. The SVM classifier using these selected features reached an accuracy of 70% and a precision of 69%.

Data Conclusion

We provided evidence in support of the ability of dMRI to discriminate between PPMS and RRMS, as well as highlight the central role of the hippocampus.

Level of Evidence

2

Technical Efficacy Stage

3

Diffusion weighted imaging in acute ischemic stroke: A review of its interpretation pitfalls and advanced diffusion imaging application

Diffusion weighted imaging (DWI) is a widely used imaging technique to evaluate patients with stroke. It can detect brain ischemia within minutes of stroke onset. However, DWI has few potential pitfalls that should be recognized during interpretation. DWI lesion could be reversible in the early hours of stroke and the entire lesion may not represent ischemic core. False negative DWI could lead to diagnosis of DWI negative stroke or to a missed stroke diagnosis. Ischemic stroke mimics can occur on DWI with non-cerebrovascular neurological conditions. In this article, the history of DWI, its clinical applications, and potential pitfalls for use in acute ischemic stroke are reviewed. Advanced diffusion imaging techniques with reference to Diffusion Kurtosis Imaging and Diffusion Tensor Imaging that has been studied to evaluate ischemic core are discussed.

Efficient and accurate EAP imaging from multi-shell dMRI with micro-structure adaptive convolution kernels and dual Fourier Integral Transforms (MiSFIT)

A number of computational techniques have been lately devised to image the Ensemble Average Propagator (EAP) within the white matter of the brain, propelled by the deployment of multi-shell acquisition protocols and databases: approaches like Mean Apparent Propagator Imaging (MAP-MRI) and its Laplacian-regularized version (MAPL) aim at describing the low frequency spectrum of the EAP (limited by the maximum b-value acquired) and afterwards computing scalar indices that embed useful descriptions of the white matter, e. g. the Return-to-Origin, Plane, or Axis Probabilities (RTOP, RTPP, RTAP). These methods resort to a non-parametric, bandwidth limited representation of the EAP that implies fitting a set of 3-D basis functions in a large-scale optimization problem. We propose a semi-parametric approach inspired by signal theory: the EAP is approximated as the spherical convolution of a Micro-Structure adaptive Gaussian kernel with a non-parametric orientation histogram, which aims at representing the low-frequency response of an ensemble of coherent sets of fiber bundles at the white matter. This way, the optimization involves just the 2 to 3 parameters that describe the kernel, making our approach far more efficient than the related state of the art. We devise dual Fourier domains Integral Transforms to analytically compute RTxP-like scalar indices as moments of arbitrary orders over either the whole 3-D space, particular directions, or particular planes. The so-called MiSFIT is both time efficient (a typical multi-shell data set can be processed in roughly one minute) and accurate: it provides estimates of widely validated indices like RTOP, RTPP, and RTAP comparable to MAPL for a wide variety of white matter configurations.

Development of Brain Structural Networks Over Age 8: A Preliminary Study Based on Diffusion Weighted Imaging

Brain structural network changes provide key information about the aging process of the brain. Unfortunately, there has yet to be a detailed characterization of these structural networks across different age groups. Efforts to classify these networks have also been hampered by their reliance on technically limited traditional methods, which are unable to track multiple fiber orientations within a voxel and consequently are prone to false detection and artifacts. In this study, a newly developed Ensemble Average Propagator (EAP) based probabilistic tractography method was applied to construct a structural network, with the strength of the link between any two brain functional regions estimated according to the alignment of the EAP along connecting pathways. Age-related changes in the topological organization of human brain structural networks were thereby characterized across a broad age range (ages 8-75 years). The data from 48 healthy participants were divided into four age groups (Group 1 aged 8-15 years; Group 2 aged 25-35 years; Group 3 aged 45-55 years; and, Group 4 aged 65-75 years; N = 12 per group). We found that the brain structural network continues to strengthen during later adolescence and adulthood, through the first 20-30 years of life. Older adults, aged 65-75, had a significantly less optimized topological organization in their structural network, with decreased global efficiency and increased path lengths versus subjects in other groups. This study suggests that probabilistic tractography based on EAP provides a reliable method to construct macroscale structural connectivity networks to capture the age-associated changes of brain structures.

Micro-structure diffusion scalar measures from reduced MRI acquisitions

In diffusion MRI, the Ensemble Average diffusion Propagator (EAP) provides relevant micro-structural information and meaningful descriptive maps of the white matter previously obscured by traditional techniques like Diffusion Tensor Imaging (DTI). The direct estimation of the EAP, however, requires a dense sampling of the Cartesian q-space involving a huge amount of samples (diffusion gradients) for proper reconstruction. A collection of more efficient techniques have been proposed in the last decade based on parametric representations of the EAP, but they still imply acquiring a large number of diffusion gradients with different b-values (shells). Paradoxically, this has come together with an effort to find scalar measures gathering all the q-space micro-structural information probed in one single index or set of indices. Among them, the return-to-origin (RTOP), return-to-plane (RTPP), and return-to-axis (RTAP) probabilities have rapidly gained popularity. In this work, we propose the so-called "Apparent Measures Using Reduced Acquisitions" (AMURA) aimed at computing scalar indices that can mimic the sensitivity of state of the art EAP-based measures to micro-structural changes. AMURA drastically reduces both the number of samples needed and the computational complexity of the estimation of diffusion properties by assuming the diffusion anisotropy is roughly independent from the radial direction. This simplification allows us to compute closed-form expressions from single-shell information, so that AMURA remains compatible with standard acquisition protocols commonly used even in clinical practice. Additionally, the analytical form of AMURA-based measures, as opposed to the iterative, non-linear reconstruction ubiquitous to full EAP techniques, turns the newly introduced apparent RTOP, RTPP, and RTAP both robust and efficient to compute.

Using the Wild Bootstrap to Quantify Uncertainty in Mean Apparent Propagator MRI

Purpose: Estimation of uncertainty of MAP-MRI metrics is an important topic, for several reasons. Bootstrap derived uncertainty, such as the standard deviation, provides valuable information, and can be incorporated in MAP-MRI studies to provide more extensive insight. Methods: In this paper, the uncertainty of different MAP-MRI metrics was quantified by estimating the empirical distributions using the wild bootstrap. We applied the wild bootstrap to both phantom data and human brain data, and obtain empirical distributions for the MAP-MRI metrics return-to-origin probability (RTOP), non-Gaussianity (NG), and propagator anisotropy (PA). Results: We demonstrated the impact of diffusion acquisition scheme (number of shells and number of measurements per shell) on the uncertainty of MAP-MRI metrics. We demonstrated how the uncertainty of these metrics can be used to improve group analyses, and to compare different preprocessing pipelines. We demonstrated that with uncertainty considered, the results for a group analysis can be different. Conclusion: Bootstrap derived uncertain measures provide additional information to the MAP-MRI derived metrics, and should be incorporated in ongoing and future MAP-MRI studies to provide more extensive insight.

On the Viability of Diffusion MRI-Based Microstructural Biomarkers in Ischemic Stroke

Recent tract-based analyses provided evidence for the exploitability of 3D-SHORE microstructural descriptors derived from diffusion MRI (dMRI) in revealing white matter (WM) plasticity. In this work, we focused on the main open issues left: (1) the comparative analysis with respect to classical tensor-derived indices, i.e., Fractional Anisotropy (FA) and Mean Diffusivity (MD); and (2) the ability to detect plasticity processes in gray matter (GM). Although signal modeling in GM is still largely unexplored, we investigated their sensibility to stroke-induced microstructural modifications occurring in the contralateral hemisphere. A more complete picture could provide hints for investigating the interplay of GM and WM modulations. Ten stroke patients and ten age/gender-matched healthy controls were enrolled in the study and underwent diffusion spectrum imaging (DSI). Acquisitions at three and two time points (tp) were performed on patients and controls, respectively. For all subjects and acquisitions, FA and MD were computed along with 3D-SHORE-based indices [Generalized Fractional Anisotropy (GFA), Propagator Anisotropy (PA), Return To the Axis Probability (RTAP), Return To the Plane Probability (RTPP), and Mean Square Displacement (MSD)]. Tract-based analysis involving the cortical, subcortical and transcallosal motor networks and region-based analysis in GM were successively performed, focusing on the contralateral hemisphere to the stroke. Reproducibility of all the indices on both WM and GM was quantitatively proved on controls. For tract-based, longitudinal group analyses revealed the highest significant differences across the subcortical and transcallosal networks for all the indices. The optimal regression model for predicting the clinical motor outcome at tp3 included GFA, PA, RTPP, and MSD in the subcortical network in combination with the main clinical information at baseline. Region-based analysis in the contralateral GM highlighted the ability of anisotropy indices in discriminating between groups mainly at tp1, while diffusivity indices appeared to be altered at tp2. 3D-SHORE indices proved to be suitable in probing plasticity in both WM and GM, further confirming their viability as a novel family of biomarkers in ischemic stroke in WM and revealing their potential exploitability in GM. Their combination with tensor-derived indices can provide more detailed insights of the different tissue modulations related to stroke pathology.