- 1Barrow Neuroimaging Innovation Center, Barrow Neurological Institute, Phoenix, AZ, United States
- 2School of Life Sciences, Arizona State University, Tempe, AZ, United States
- 3The Muhammad Ali Parkinson's Center, Barrow Neurological Institute, Phoenix, AZ, United States
Introduction: Parkinson’s disease (PD) is a progressive neurodegenerative disorder characterized by motor and non-motor symptoms. Levodopa remains the primary treatment, temporarily restoring dopamine levels and improving motor symptoms. Advanced diffusion MRI techniques, such as free-water corrected diffusion tensor imaging (fw-DTI) and neurite orientation dispersion and density imaging (NODDI), provide insights into PD-related microstructural changes beyond conventional DTI.
Methods: This study investigates white matter alterations in PD using multi-shell fw-DTI and NODDI to compare voxel-wise differences between PD patients both OFF and ON levodopa, with comparison to healthy controls (HC). Effect sizes and receiver operating characteristic (ROC) analyses assessed the discriminative power of imaging metrics.
Results: PD (OFF) exhibited increased free-water, reduced neurite density (NDI), and altered orientation dispersion (ODI) in key motor pathways in comparison to HC, while fw-FA offered robust group discrimination (AUC=0.956). Levodopa (ON state) increased NDI and NODDI-FWF, suggesting acute microstructural plasticity, though this finding contrasted with minimal fw-DTI FW changes. Additionally, voxel-based correlation analyses linked free-water and neurite integrity metrics with disease severity.
Discussion: Our findings suggest that fw-DTI and NODDI provide complementary information on PD-related neurodegeneration and the transient effects of levodopa. These results underscore the potential of advanced diffusion MRI techniques as biomarkers for tracking PD progression and treatment response.
1 Introduction
Parkinson’s disease (PD) is the second most common neurodegenerative disorder after Alzheimer’s disease, affecting an estimated 6 million individuals worldwide in 2015, with projections suggesting this number may exceed 12 million by 2040 (1, 2). This progressive disease is primarily characterized by alpha-synuclein deposition in central and peripheral nervous system and degeneration of dopaminergic neurons in the substantia nigra pars compacta, leading to the hallmark motor manifestations of bradykinesia, resting tremor, rigidity, and postural instability (3–5). While the causes of PD remain unclear, the degeneration of dopaminergic neurons in the substantia nigra (SN) leads to decreased dopaminergic innervation in the striatum, which is a key factor in the pathophysiology of PD motor symptoms (4, 6).
Levodopa (L-DOPA) remains the most effective pharmacological treatment for PD, often providing marked improvement in motor function when patients are in the so-called therapeutically ON state (7). L-DOPA is able to cross the blood–brain barrier, where it is enzymatically converted into dopamine, temporarily restoring dopaminergic signaling, and improving motor symptoms (8). This mechanism directly targets the core pathophysiology of PD by replenishing dopamine levels in the striatum, a brain region essential for regulating motor function (9). The OFF and ON states experienced by patients reflect dynamic changes in dopaminergic signaling and can significantly impact a patient’s quality of life. Therefore, investigating the OFF and ON states is crucial for understanding disease progression, optimizing treatment regimens, and developing more targeted interventions.
In recent years, magnetic resonance imaging (MRI) has played an important role in studying the structural and functional changes associated with PD. Although conventional MRI sequences often fail to detect PD-specific alterations, especially in the early stages, advanced MRI techniques have shown promise in identifying microstructural and functional changes that correlate with disease severity and progression (10, 11). One advanced MRI method is diffusion MRI (dMRI), which is particularly sensitive to disruptions in the microstructural integrity of white matter by quantifying the diffusion of water molecules within brain tissue (12, 13).
The most commonly used dMRI technique is Diffusion Tensor Imaging (DTI), yielding metrics such as fractional anisotropy (FA) and mean diffusivity (MD). These biomarkers have been investigated in various regions affected by PD, including the substantia nigra, basal ganglia, and cortical areas (14, 15). Changes in these diffusion parameters have been connected to neurodegeneration and disease severity (16). While these biomarkers have proven effective in identifying chronic neurodegenerative changes, previous work has suggested a lack of dynamic differences between the OFF and ON states following L-DOPA administration (17).
Free-water DTI (fw-DTI) extends the standard DTI model to distinguish fw in extracellular spaces and water restricted within neural tissues. This differentiation is crucial for accurately assessing brain microstructure because the fw content in the brain can confound traditional DTI metrics (18). In the context of PD, fw-DTI has been used to identify microstructural white matter changes, with studies reporting elevated free-water levels in the posterior SN of PD patients that increase over time (14, 19). Furthermore, fw-DTI metrics have been shown to correlate with motor scores and cognition, as well as disease progression (20, 21).
Another advanced dMRI technique is Neurite Orientation Dispersion and Density Imaging (NODDI), which is designed to overcome some of the limitations of traditional DTI using multi-shell dMRI acquisitions to provide more specific metrics of microstructural complexity in the brain (22). NODDI involves a multi-compartment model that distinguishes intra-neurite, extra-neurite, and cerebrospinal fluid (CSF) contributions to the overall diffusion signal. Therefore, NODDI can estimate indices such as the Neurite Density Index (NDI) and Orientation Dispersion Index (ODI), which reflect the density and angular variability of neurites (axons and dendrites) within a brain voxel, respectively. Similar to FW-DTI, NODDI provides estimates of the free water fraction (FWF) inside the brain (22).
Previous studies have applied NODDI in PD to identify microstructural alterations in the substantia nigra and other nigrostriatal pathways—regions critically involved in PD pathology (23). For instance, reductions in the NDI may reflect decreased neurite density that is consistent with dopaminergic neuronal loss, while altered ODI might signify changes in neurite organization or orientation as the disease progresses (24, 25). Moreover, Wei et al. showed that NODDI metrics were significantly correlated with objective measures of gait in PD patients, particularly in the ON state (26).
In this study, multi-shell fw-DTI and NODDI were employed to investigate voxel-based microstructural differences in a cohort of 20 PD patients between their OFF and ON (levodopa) states. Additionally, voxel-wise differences were assessed between PD patients in the OFF state and a cohort of healthy controls (HC). For the PD group, voxel-based correlations were performed between clinical measures [MDS-Unified Parkinson’s Disease Rating Scale (27) and Hoehn and Yahr (H&Y) (28)] and fw-DTI/NODDI-derived metrics. This multimodal imaging approach may provide insights into how neurodegenerative processes and levodopa-induced changes influence brain microstructure, thereby enhancing the characterization of Parkinson’s disease pathology.
2 Methods
2.1 Participants
This study was performed in accordance with the local Institutional Review Board. Informed consent was obtained from all participants in this study, which was conducted in compliance with the Health Insurance Portability and Accountability Act (HIPAA). A total of 20 patients with non-demented, treated PD (25% female; mean age ± standard deviation (SD): 66.95 ± 8.86 years) and 16 healthy controls (HC) (43.75% female; mean age ± SD: 66.88 ± 7.44 years) were included in this study. All PD subjects were recruited by a movement disorders specialist, while HC were recruited through word-of-mouth referrals or our local volunteer database between September 2019 and November 2024.
Inclusion criteria required participants to be 40–90 years of age capable of providing informed consent, with PD patients maintaining stable antiparkinsonian medication and HCs being age-matched without neurological/psychiatric conditions. Exclusion criteria for both groups included MRI contraindications (bioimplants, aneurysm clips), severe claustrophobia/vertigo, marked tremor, or use of anticholinergic/dopamine-blocking drugs.
PD participants underwent motor assessment using the MDS-UPDRS-III, which evaluates motor symptoms such as rigidity, bradykinesia, tremor, gait, and postural stability through clinician-rated tasks. Each item is scored on a scale from 0 (normal) to 4 (severe impairment), with higher total scores indicating greater motor dysfunction (27). The H&Y scale was used to classify disease severity based on motor symptoms and functional impairment, ranging from mild unilateral involvement to severe disability requiring full-time assistance (28).
All study participants underwent MRI scanning, with HC subjects scanned once and PD subjects scanned twice, first in the OFF state (defined as >12 h since their last dose of levodopa) and subsequently in the ON state (practically defined as >1 h following levodopa administration). Following the OFF scans, PD patient received their typical morning dose of medication. Levodopa equivalent daily dose (LEDD) was calculated for all subjects1 (29).
2.2 MRI acquisition
Diffusion MRI data were acquired using a 3.0 T Philips Ingenia scanner. Imaging was performed using 2D Echo-Planar Imaging (EPI) with the following acquisition parameters: field of view (FOV) = 256 × 256 mm2, matrix size = 128 × 128, number of slices = 72, voxel size = 2.0 × 2.0 × 2.0 mm3, repetition time (TR) = 7,000 ms, echo time (TE) = 121 ms, and flip angle = 90°.
The multi-shell dMRI acquisition included one non-diffusion-weighted (b = 0 s/mm2, b0) image and four diffusion-weighted shells with b-values of 500, 1,000, 2,000, and 3,000 s/mm2, each acquired with 20 diffusion directions using phase-encoding along the antero-posterior (A/P) direction. Additionally, a b = 0 s/mm2 image with reversed phase-encoding polarity (posterior–anterior, P/A) was acquired to enable EPI distortion correction.
2.3 dMRI pre-processing
Diffusion MRI pre-processing was performed using MRtrix3 (version 3.0.4-145) (30), FSL (version 6.0.7.16) (31), and ANTs (version v2.4.4).2 The preprocessing pipeline included denoising using dwidenoise (MRtrix3) (32), EPI distortion correction, eddy current correction, and motion correction using TOPUP and eddy (FSL) (33). The eddy quality control tools were used to assess the quality of the DTI dataset (34). Slices with signal loss due to subject movement coinciding with diffusion encoding were identified and replaced using Gaussian process predictions (35). Subjects with more than 5% total outliers and/or an average absolute volume-to-volume head motion exceeding 3 mm were excluded.
Prior to MNI registration, the preprocessed dMRI native images in 2 mm were upsampled to 1.25 mm isotropic resolution using mrgrid (MRtrix3) (32) to reduce resolution mismatch and improve the accuracy of subsequent non-linear registration to the 1 mm MNI template using ANTs SyN. All brain extraction was performed on the upsampled b0 images using dwi2mask (MRtrix3) (36).
Subsequently, brain-extracted b0 images were non-linearly coregistered to the MNI standard space (1 × 1 × 1 mm) using the ANTs SyN registration algorithm (37). The resulting deformation fields were used to transform diffusion-derived maps into standard space for group-level analyses. Therefore, for each participant, fw-DTI and NODDI maps (see Free-Water Algorithm and NODDI sub-sections) were coregistered to MNI space by applying both the registration matrix and the warp field obtained from the b0 coregistration using ANTs. For the final voxel-based analysis, spatial smoothing with a 4 mm FWHM Gaussian kernel was applied to enhance signal-to-noise ratio and improve sensitivity to group differences, while maintaining a balance with spatial specificity relative to the native resolution.
2.4 Free water algorithm
The fw elimination model aims to mitigate the negative impact of the CSF partial volume effects on diffusion measurements (38). By distinguishing between freely moving water molecules and those that are hindered or restricted, this model provides a more accurate characterization of brain tissue microstructure. Multi-shell dMRI data were used for fw-DTI to overcome the known limitations of single-shell dMRI (39).
In this study, fw-corrected diffusion metrics, including fractional anisotropy (fw-FA) and the fw-index, were computed using the DIPY (version 1.10.0) (40) algorithm through the fwdti.FreeWaterTensorModel class and object and fwdtimodel.fit function. The multi-shell dMRI data were processed via an in-house Python (version 3.11.5) script.
2.5 Neurite orientation dispersion and density imaging (NODDI)
NODDI is an advanced dMRI model designed to provide a more detailed characterization of brain microstructure by distinguishing between different tissue compartments (22). Unlike conventional DTI, which assumes Gaussian diffusion, NODDI accounts for the complexity of neurite structures by modeling the diffusion signal as a combination of three compartments: the intra-neurite compartment, representing diffusion within axons and dendrites and characterized by restricted diffusion; the extra-neurite compartment, reflecting diffusion in the extracellular space where water diffusion is hindered by cellular structures; and the isotropic (fw) compartment, corresponding to unrestricted diffusion primarily associated with cerebrospinal fluid (CSF) and extracellular free water.
Using a multi-shell dMRI acquisition, NODDI estimates several key microstructural metrics including the neurite density index (NDI), which represents the volume fraction of the intra-neurite space and is an indicator of axonal and dendritic density; the orientation dispersion index (ODI), which quantifies the variability in neurite orientation and reflects the complexity of neural organization; and the isotropic volume fraction (FWF), which measures the fraction of freely diffusing water and is useful for detecting CSF contamination and neuroinflammation.
NODDI maps were generated using an in-house Python script implementing the Accelerated Microstructure Imaging via Convex Optimization (AMICO, version 2.0.3) framework (41), which enables efficient and robust estimation of NODDI metrics.3
2.6 Statistical analyses
The demographic and clinical characteristics, including age, sex, MDS-UPDRS-III, H&Y, motion, and outliers, are presented as means and standard deviations for each group.
Differences in age between groups were examined using a Student’s t-test, confirmed by the Shapiro–Wilk (SW) test for normality, while differences in sex were evaluated using the χ2 test.
Motor symptoms and disease severity for the PD group between OFF and ON conditions were analyzed using the paired t-test for the MDS-UPDRS-III test and the Wilcoxon Signed-Rank Test for the H&Y test (Table 1). Differences in tremor, bradykinesia, rigidity, and axial subscores between ON and OFF states were evaluated by the Wilcoxon Signed-Rank Test, with p-values reported as both uncorrected and corrected for multiple comparisons by Bonferroni’s method. Group differences in motion and outliers across groups or ‘OFF/ON’ status were analyzed using the Wilcoxon Rank-Sum Test, Wilcoxon Signed-Rank Test, and paired t-test. Complete information can be found in Table 1.

Table 1. Demographic, clinical, and motion characteristics of healthy controls (HC) and Parkinson’s disease (PD) patients in OFF and ON states.
Voxel-based analyses, including between-group comparisons [HC vs. PD (OFF)], ‘OFF/ON’ status assessments, and correlations with cognitive test scores, were restricted to voxels within the combined masks derived from the ICBM-DTI-81 white matter labels atlas and the JHU white matter tractography atlas (42, 43). Statistically significant clusters were identified and annotated based on these atlases.
To evaluate voxel-based differences between HC and PD (OFF) across all diffusion-related metrics, a Student’s t-test was performed using a linear model. Age, sex, and average absolute motion were included as covariates to control for potential confounding factors. For within-subject comparisons between PD (OFF) and PD (ON), a paired Student’s t-test was applied, with average absolute motion as covariate. Additionally, for the PD group (both OFF and ON conditions), voxel-based correlations between diffusion metrics and cognitive scores were computed using a linear model, adjusting for average absolute motion. All analyses were conducted using an in-house Python script.
For voxel-based statistical analyses, the Threshold-Free Cluster Enhancement (TFCE) method was employed to enhance cluster detection while avoiding arbitrary thresholding and addressing multiple comparisons (44). Additionally, a Family-Wise Error (FWE) correction at the 0.05 level was applied using the Benjamini-Hochberg procedure (45) (FDR < 0.05). Effect sizes were computed for all analyses using Hedges’ g (46), where |g| ≥ 1.03 and |g| ≥ 0.61 indicated a large effect size for differences between HC and PD (OFF) and differences between PD (OFF) and PD (ON), respectively (both at α = 0.05, power = 0.85). For correlations with motor symptoms and disease severity, Spearman’s correlation coefficient (ρ) (47) was used, along with the corresponding t value.
To assess the discriminative power of NODDI and fw-DTI metrics, we computed the area under the receiver operating characteristic curve (AUC) for each metric. The AUC distributions were estimated using a nonparametric bootstrap resampling approach with 10,000 iterations. This method provides robust estimates of AUC variability by resampling the data with replacement, generating empirical distributions that account for the uncertainty in AUC estimation.
For each metric, we computed the mean AUC along with 95% confidence intervals (CIs) derived from the bootstrap distribution. The confidence intervals were estimated using the bias-corrected and accelerated (BCa) bootstrap method to account for potential skewness in the AUC distribution.
Additionally, to statistically compare the performance of different metrics, pairwise differences in AUC values were computed using the bootstrap method. The differences in AUCs were calculated for all metric pairs, and corresponding 95% confidence intervals were derived. A statistically significant difference was inferred when the confidence interval did not include zero, indicating a meaningful distinction in classification performance between the metrics.
All AUC related analyses were performed using R [version 4.4.2 (2024-10-31)] and RStudio (version 2024.09.0) by using the pROC and boot packages.
3 Results
For the PD cohort, 19 participants received levodopa, with 7 of them also taking amantadine, 3 taking rasagiline, and 2 taking ropinirole. One participant was just on a combination of amantadine and rasagiline. The mean time between the last medication dose and the OFF scan was 12.8 ± 0.1 h, while the mean time between drug administration and the ON scan was 68.9 ± 6.0 min. The mean disease duration was 5.5 ± 2.9 years. Demographic and clinical characteristics of the recruited individuals are summarized in Table 1, while disease duration and medications are summarized for PD participants in Table 2.
No significant difference in age was found between HC and PD (t = −0.028, p = 0.978). However, a significant difference was observed between the two groups for sex with males more common in the PD group (χ2 = 38.45, p < 0.001).
For within-subject comparisons between PD (OFF) and PD (ON), a significant difference was found for MDS-UPDRS-III (t = 6.717, p < 0.001), while no significant difference was observed for the H&Y score (V = 4.0, p = 0.773).
No significant differences were detected in the dMRI-derived motion parameters, including absolute motion, relative motion, or total outliers, across groups or between the OFF and ON stages. In all final analyses, no subjects were excluded due to motion or outlier effects. All statistical results for cognitive scores, motion, and outliers are summarized in Table 1.
Table 1 shows tremor, bradykinesia, rigidity, and axial subscores for PD patients in ON and OFF states based on MDS-UPDRS-III sub-sections. Significant statistical differences between the OFF and ON stages were found for all subscores (p < 0.001).
3.1 fw-DTI results: HC vs. PD (OFF)
The voxel-wise analysis identified statistical differences between HC and PD (OFF) in fw-index and fw-FA metrics (Figures 1a,b, Table 3(a), and Supplementary Figure 1).

Figure 1. (a,b) Voxel-wise comparison of free-water index (fw-index) and free-water–corrected fractional anisotropy (fw-FA) between healthy controls (HC) and Parkinson’s Disease (PD) patients in the OFF medication state. (c,d) Comparison of fw-index and fw-FA between PD patients in the OFF and ON states. Statistical maps are overlaid on MNI-space anatomical images. Significant group differences (FWE-corrected p < 0.05) are displayed using red-to-yellow color scales for increased values in PD (t > 0) and blue-to-light blue for increased values in HC or PD ON (t < 0). Clusters with large effect sizes are outlined in green (g positive) and purple (g negative).

Table 3. (a) Results of the voxel-based analysis for fw-DTI between HC and PD (OFF). The analysis includes fw-FA and the fw-index. Results are reported as statistical values (t-values) and effect sizes. The percentage (%) represents the volume of the cluster within the corresponding white matter location based on the atlas. (b) Results of the voxel-based analysis for fw-DTI between PD (OFF) and PD (ON). The analysis includes fw-FA and the fw-index. Results are reported as statistical values (t-values) and effect sizes. The percentage (%) represents the volume of the cluster within the corresponding white matter location based on the atlas.
The fw-index was significantly higher in PD (OFF) than in HC, with statistical differences across large portions of the white matter. Additionally, fw-FA was also elevated in PD (OFF) in several white matter regions, including the right and left anterior thalamic radiation (t = 3.57, g = 1.14; t = 3.46, g = 1.19, respectively), left superior longitudinal fasciculus (t = 3.74, g = 1.14), right and left cerebral peduncle (t = 3.68, g = 1.18; t = 3.76, g = 1.17, respectively), and the right and left anterior limb of the internal capsule (t = 3.53, g = 1.20; t = 3.40, g = 1.12, respectively).
Conversely, fw-FA was significantly lower in PD (OFF) than HC in a large cluster within the left retrolenticular part of the internal capsule (t = −3.35, g = 1.17). The differences between the two groups were also supported by effect size (with large effect-size values) for both metrics.
3.2 fw-DTI results: PD (OFF) vs. PD (ON)
Figures 1c,d, Table 3 (b), and Supplementary Figure 1 present the voxel-based statistical differences between PD (OFF) and PD (ON). We observed higher values for both fw-related metrics in the ON state compared to the OFF state. However, these differences were limited to small clusters (cluster volume less than 1% of the corresponding white matter area). Additionally, these findings were not supported by the effect size analysis.
3.3 NODDI results: HC vs. PD (OFF)
This study analyzed three NODDI-derived metrics: FWF, ODI, and NDI. Figures 2a–c, Table 4, and Supplementary Figure 2 present the voxel-based differences between HC and PD (OFF) patients for these metrics.

Figure 2. (a–c) Voxel-wise comparison of NODDI-derived metrics between healthy controls (HC) and Parkinson’s Disease (PD) patients in the OFF medication state, including free-water fraction (FWF), orientation dispersion index (ODI), and neurite density index (NDI), respectively. (d–f) Comparison of the same NODDI metrics between PD patients in the OFF and ON states. Statistical group differences are overlaid on MNI-space anatomical images. Significant regions (p < 0.05, FWE-corrected) are shown using red-to-yellow color scales for positive t-values and blue-to-light blue for negative t-values. Green and purple contours indicate clusters with large effect sizes.
Comparisons between the two groups revealed significant differences across all metrics. FWF values were higher in PD (OFF) compared to HC in multiple white matter regions, with large effect sizes observed mainly in the forceps minor (t = 2.83, g = 1.46), left superior longitudinal fasciculus (t = 3.08, g = 1.14), middle cerebellar peduncle (t = 2.69, g = 1.20), the body of the corpus callosum (t = 2.68, g = 1.13), and right superior corona radiata (t = 2.50, g = 1.17).
For ODI, significant increases in PD (OFF) were detected in multiple white matter regions, with the largest significant clusters located in the left retrolenticular part of the internal capsule (t = 3.57, g = 1.22) and the left sagittal stratum (t = 2.91, g = 1.04). Additionally, several smaller clusters showed reduced ODI values in PD (OFF) compared to HC in different white matter locations. However, they were not confirmed by the effect size at a large effect level.
Lastly, NDI was significantly lower in PD (OFF) compared to HC, with large effect sizes confirming this difference mainly inside the left anterior thalamic radiation (t = −3.10, g = −1.16) and forceps minor (t = −3.15, g = −1.17).
3.4 NODDI: PD (OFF) vs. PD (ON)
Figures 2d–f, Table 5, and Supplementary Figure 2 show the voxel-based statistical differences in NODDI-related metrics between PD patients in the OFF and ON medication states.
FWF values were higher in the ON state compared to the OFF state across several white matter regions, with large effect sizes observed mainly in the right corticospinal tract (t = −3.63, g = −0.66), left inferior fronto-occipital fasciculus (t = −3.45, g = −0.73), left uncinate fasciculus (t = −3.43, g = −0.74), right posterior limb of the internal capsule (t = −3.65, g = −0.66), right posterior thalamic radiation (t = −4.25, g = −0.68), and right sagittal stratum (t = −4.09, g = −0.71).
Only minimal differences were detected between the OFF and ON states in both statistical directions for ODI. However, these findings were not supported by large effect sizes.
Finally, compared to the OFF state, higher NDI values were observed in the ON state within the right inferior fronto-occipital fasciculus (t = −3.56, g = −0.81), right inferior longitudinal fasciculus (t = −3.61, g = −0.72), and right sagittal stratum (t = −3.60, g = −0.63).
3.5 fw-DTI metrics and correlations with MDS-UPDRS III/H&Y scores
Figures 3, 4 show voxel-based correlations between fw-DTI metrics and MDS-UPDRS-III/H&Y in the PD group OFF and ON states.

Figure 3. Voxel-based correlations between fw-DTI metrics and MDS-UPDRS-III scores in the Parkinson’s Disease (PD) group. (a) Positive correlations with the fw-index are shown warm colors (red-yellow). (b) Negative correlations with the fw-FA index are shown cool colors (blue-light blue). The color bars represent t-values and corresponding correlation coefficients (ρ). Differences between OFF and ON states (green for OFF and red for ON) indicate a reduction in correlation strength during the ON state.

Figure 4. Voxel-based correlations between fw-DTI metrics and H&Y scores in the Parkinson’s Disease (PD) group. (a) Positive correlations with the fw-index are shown warm colors (red-yellow). (b) Negative correlations with the fw-FA index are shown cool colors (blue-light blue). The color bars represent t-values and corresponding correlation coefficients (ρ). Differences between OFF and ON states (green for OFF and red for ON) indicate the variation in correlation strength during the ON state.
Significant positive correlations were observed between the fw-index and MDS-UPDRS-III in both OFF and ON states, though these correlations were reduced during the ON state. Conversely, significant negative correlations were found between the fw-FA index and MDS-UPDRS-III, with a similar reduction in correlation strength observed in the ON state compared to OFF. For both metrics, a prominent statistical cluster was identified within the splenium of the corpus callosum.
Similarly, significant positive correlations were detected between the fw-index and H&Y scores in both OFF and ON states. In contrast, significant negative correlations were found between fw-FA and H&Y across both conditions. Differences in correlation strength were also observed between the OFF and ON states.
Table 6 provides a summary of these findings, including cluster volumes, t-values, ρ-values, and the percentage change in volume between the OFF and ON states.

Table 6. Voxel-based correlations between fw-DTI metrics and clinical scores (MDS-UPDRS-III and Hoehn and Yahr) in PD groups during the OFF and ON stages.
3.6 NODDI metrics and correlations with MDS-UPDRS III/H&Y scores
Figures 5, 6 show voxel-based correlations between NODDI metrics and MDS-UPDRS III/H&Y in the PD group across OFF and ON states.

Figure 5. Voxel-based correlations between NODDI metrics and MDS-UPDRS-III scores in the Parkinson’s Disease (PD) group. (a,b) Positive correlations with the FWF and ODI are shown warm colors (red-yellow). (c) Negative correlations with the NDI are shown cool colors (blue-light blue). The color bars represent t-values and corresponding correlation coefficients (ρ). Differences between OFF and ON states (green for OFF and red for ON) indicate a reduction in correlation strength during the ON state.

Figure 6. Voxel-based correlations between NODDI metrics and H&Y scores in the Parkinson’s Disease (PD) group. (a,b) Positive correlations with the FWF and ODI are shown warm colors (red-yellow). (c) Negative correlations with the NDI are shown cool colors (blue-light blue). The color bars represent t-values and corresponding correlation coefficients (ρ). Differences between OFF and ON states (green for OFF and red for ON) indicate a reduction in correlation strength during the ON state.
Significant positive correlations were observed between MDS-UPDRS-III and both FWF and ODI metrics, whereas significant negative correlations were found for the NDI. Additionally, differences in correlation patterns were detected between the OFF and ON states.
Similar trends were observed for the H&Y scores, with significant correlations in the same direction of MDS-UPDRS-III for all NODDI metrics.
Table 7 provides a summary of these findings, including cluster volumes, t-values, ρ-values, and the percentage change in volume between the OFF and ON states.

Table 7. Voxel-based correlations between NODDI metrics and clinical scores (MDS-UPDRS III and Hoehn and Yahr) in PD groups during OFF and ON states.
3.7 Bootstrap distribution of AUCs and differences in AUC
Figure 7 illustrates the bootstrap distributions of AUC values for fw-DTI and NODDI metrics, along with the pairwise comparisons of AUC differences. The classification algorithm was used to differentiate between HC and PD in OFF conditions. The AUC distributions were calculated within the significant clusters identified after TFCE and FDR correction (p < 0.05).

Figure 7. Bootstrap distribution of area under curve (AUC) values (using 10,000 iterations) for fw-DTI and NODDI metrics. The density plot illustrates the distribution of AUC values obtained through bootstrapping for each metric. The accompanying tables summarize the mean AUC values with 95% confidence intervals (CI) and pairwise differences in AUC between metrics.
Across all diffusion-related metrics, the highest classification performance was observed for fw-FA, with a mean AUC of 0.956 (95% CI: 0.883–1.000), followed by ODI (AUC = 0.935, 95% CI: 0.847–0.993) and FWF (AUC = 0.926, 95% CI: 0.883–0.992). The fw-index also showed strong discriminative power with an AUC of 0.907 (95% CI: 0.793–0.984), whereas NDI had the lowest AUC at 0.872 (95% CI: 0.742–0.968).
Pairwise comparisons of AUC values revealed differences across the diffusion metrics. The largest differences were found between fw-FA and NDI (<ΔAUC> = 0.085, 95% CI: −0.032 to 0.225).
4 Discussion
Today, levodopa remains the most effective treatment for Parkinson’s disease (PD), alleviating many of the motor symptoms. Prolonged use of levodopa leads to motor fluctuations between OFF (withdrawn) and ON (post-levodopa) states, which can significantly impact daily functioning and quality of life (7, 48). Advanced dMRI techniques, such as multi-shell fw-DTI (fw-DTI) (39) and NODDI (22), offer deeper insights into PD-related microstructural alterations beyond conventional DTI. In this study, we applied multi-shell fw-DTI and NODDI to investigate voxel-wise white matter differences in 20 PD patients across OFF and ON states and compared PD (OFF) with a cohort of HC participants. By correlating imaging findings with clinical scores, we aimed to enhance the understanding of PD-related neurodegeneration and levodopa-induced neural changes.
To ensure precise characterization of microstructural alterations, we implemented a robust dMRI processing pipeline and advanced statistical analyses. Traditional single-shell dMRI, which uses only one b-value, is constrained by a simplified diffusion model, making it challenging to separate free-water contamination from tissue-specific diffusion properties. In contrast, multi-shell dMRI samples diffusion at multiple b-values, enabling more precise separation of the free-water and tissue compartments (39). This improved modeling reduces biases inherent in single-shell methods, leading to more accurate estimates of white matter microstructure and better sensitivity to pathological changes in neurodegenerative diseases like PD (18). Therefore, multi-shell fw-DTI can offer improved sensitivity and specificity in detecting microstructural alterations in neurodegenerative diseases. Additionally, NODDI provides complementary information on neurite density and orientation dispersion, further refining the assessment of PD-related white matter integrity (22).
The use of multi-shell fw-DTI and NODDI was guided by pathophysiological features specific to PD. In PD, free-water accumulation in the posterior substantia nigra—measured via fw-DTI—has been strongly associated with dopaminergic neuron loss (19) and correlates with disease progression (14). Unlike conventional DTI metrics (e.g., FA, MD), which conflate extracellular fluid changes with tissue damage, fw-DTI uses compartmental modeling to isolate free-water effects, enabling earlier detection of neurodegeneration (38). For example, Ofori et al. showed that increases in nigral free-water precede detectable FA changes in early PD, underscoring fw-DTI’s sensitivity to early pathological changes (19).
NODDI, in contrast, is a complementary biophysical model that addresses limitations of traditional DTI and DKI—such as the assumption of Gaussian diffusion and nonspecific metric changes. Several studies support NODDI’s utility for PD diagnosis and tracking disease progression (23, 25, 49). Together, fw-DTI and NODDI overcome the constraints of single-shell DTI, offering enhanced specificity in distinguishing overlapping pathological processes—a key step toward developing stage-specific biomarkers for PD.
In this study, both age and sex were included as a covariate in the statistical analysis to account for potential confounding effects. Our results showed no significant age differences between HC and PD groups, though a significant sex difference was observed. Sex differences have been reported in white matter microstructure, including variations in fractional anisotropy and mean diffusivity (50, 51). By controlling for sex, we ensure that observed group differences in DTI metrics are not driven by underlying sex-related variations in brain structure. Motion parameters were comparable across groups and conditions, ensuring no motion-related bias in imaging analyses.
Significant differences in fw-DTI metrics between HC and PD (OFF) groups indicate widespread microstructural alterations in PD-related white matter pathways. The elevated fw-index in PD (OFF) suggests increased extracellular free water, potentially reflecting neuroinflammatory or neurodegenerative processes. The regions with increased fw-index values included the right cingulum (cingulate gyrus), body of the corpus callosum, and superior corona radiata—regions integral to motor control and cognitive function (52–54). Our results align with previous fw-DTI studies on PD-related white matter alterations. Ofori et al. reported increased free-water content in the posterior substantia nigra, suggesting neurodegenerative changes in this region (14), while Planetta et al. found elevated free-water levels in motor-related pathways correlated with disease severity (55). These findings suggest that fw-DTI could provide a non-invasive biomarker of neurodegenerative changes in PD, especially in white matter regions implicated in motor functions. Clinically, this could support earlier diagnosis or stratification of patients based on disease burden.
Similarly, fw-FA was significantly higher in PD (OFF) across multiple white matter tracts, including the anterior thalamic radiation (bilaterally), left superior longitudinal fasciculus, and anterior limb of the internal capsule. These findings, supported by large effect sizes, may reflect underlying pathological changes or adaptive responses in PD. In contrast, a significant decrease in fw-FA in the left retrolenticular part of the internal capsule highlighted regional variability in PD-associated microstructural changes. Increased FA values in PD compared to HC have also been reported. Mole et al. found significantly higher FA in motor pathways, including the corticospinal and thalamocortical tracts, possibly indicating compensatory neuroplasticity or selective neurodegeneration (56). Additionally, Lenfeldt et al. observed increased FA in the substantia nigra, potentially reflecting gliosis, inflammation, or surrounding fiber intrusion into the shrinking structure (57).
Voxel-based comparisons between PD (OFF) and PD (ON) states revealed higher fw-related metric values in the ON state. However, these differences were confined to small clusters (<1% of the corresponding white matter region) and were not supported by effect size analysis, suggesting limited clinical relevance. The absence of widespread free-water changes following levodopa administration may reflect the short timeframe between OFF and ON states, which is insufficient for large-scale neurochemical or microstructural alterations to manifest. These findings appear consistent with a study by Chung et al., which explored the use of single-shell fw-DTI in PD patients. The study found no significant differences in fw-FA or fw-index between the OFF and ON states following L-DOPA administration. This suggests that fw-DTI may be better suited for evaluating chronic structural changes in PD rather than acute, medication-induced fluctuations (17).
NODDI-derived metrics revealed significant microstructural differences between HC and PD (OFF). Increased FWF in PD (OFF) across multiple white matter regions—particularly the forceps minor, superior longitudinal fasciculus, middle cerebellar peduncle, body of the corpus callosum, and right superior corona radiata—suggested elevated extracellular water, indicative of neuroinflammatory or neurodegenerative processes. Large effect sizes confirmed the robustness of these findings. Mitchell et al. previously compared NODDI and fw-DTI biomarkers in PD and atypical Parkinsonism. Both techniques detected microstructural changes, with fw imaging showing higher effect sizes for extracellular changes (e.g., elevated free water in the substantia nigra). NODDI provided additional insights into intracellular volume and orientation dispersion. Both methods were able to discriminate PD from atypical Parkinsonism, with fw-DTI being more sensitive to extracellular neurodegeneration (24).
Interestingly, the FWF differences observed in the PD OFF vs. ON condition were predominantly right-lateralized. While our cohort was mostly right-handed (16 right-handed, 2 left-handed, 2 ambidextrous), information regarding symptom lateralization was not collected for this study. Given that PD often presents with asymmetric motor symptoms and lateralized neurodegeneration, this pattern may reflect underlying hemispheric differences in disease burden. Handedness and symptom asymmetry are potential modulators of such imaging findings, and future studies with detailed clinical lateralization data are warranted to further explore these associations.
The findings from both fw-DTI and NODDI highlight the utility of free-water imaging in PD. Both models identified increased extracellular water in PD (OFF), with FWF (from NODDI) and fw-index (from fw-DTI) showing consistent results across overlapping white matter regions, such as the superior longitudinal fasciculus. This agreement suggests that fw metrics are robust biomarkers for neuroinflammatory or neurodegenerative processes in PD. Additionally, FWF may offer advantages over the fw-index (from multi-shell fw-DTI) due to its ability to disentangle intracellular and extracellular contributions more precisely, providing a more nuanced understanding of tissue microstructure.
In addition to FWF changes, NODDI-based ODI was significantly increased in PD (OFF) in several white matter regions, including the left retrolenticular internal capsule and left sagittal stratum, suggesting greater neurite orientation complexity, possibly due to compensatory mechanisms or altered connectivity. However, regions with lower ODI values in PD (OFF) than HC were also identified, though without strong effect size support, indicating regional variability in disease-related changes.
NDI was significantly lower in PD (OFF), with the most significant reductions in the left anterior thalamic radiation and forceps minor, reflecting decreased neurite density consistent with known PD-related neurodegeneration. Large effect sizes supported the relevance of these findings. Similar to these findings, Kamagata et al. previously reported reduced intracellular volume fraction (similar to NDI) in the substantia nigra pars compacta (SNpc) and putamen, correlating with disease severity (25). Additionally, Wei et al. investigated white matter microstructural alterations in PD patients compared to HCs using NODDI, identifying significant changes in NDI, ODI, and ISO in tracts such as the anterior thalamic radiation and corticospinal tract. These changes correlated with gait impairments, highlighting axonal degeneration’s role in motor deficits (26).
For NODDI metrics, white matter changes were also observed between pre- and post-levodopa. Higher FWF in the ON state, particularly in the right corticospinal tract, left inferior fronto-occipital fasciculus, and right posterior thalamic radiation, suggests transient neurochemical changes or medication-induced fluid regulation shifts. Significant NDI increases in associative and visual processing pathways (e.g., right inferior fronto-occipital and inferior longitudinal fasciculus) could suggest transient neurite density enhancement, possibly due to improved neurotransmission or compensatory mechanisms. However, ODI differences between OFF and ON states were minimal and lacked strong effect sizes, suggesting a limited impact of levodopa on neurite orientation.
Within-subject comparisons revealed a significant improvement in motor function after levodopa administration (via the MDS-UPDRS-III). At the same time, the H&Y score remained unchanged, likely due to its categorical nature and lower sensitivity to short-term motor fluctuations. The voxel-based correlation analysis between fw-DTI/NODDI metrics and MDS-UPDRS-III/H&Y revealed significant associations between white matter microstructure and disease severity in both OFF and ON states. A positive correlation was observed between the fw-index and MDS-UPDRS-III scores in both conditions, suggesting that greater free-water content is associated with more severe motor impairment. Additionally, the volumes of the correlation clusters were reduced in the ON state. Similarly, negative correlations between fw-FA and MDS-UPDRS-III indicate that reduced white matter integrity is linked to worse motor function. These correlations were strongest in the splenium of the corpus callosum, a region implicated in interhemispheric communication and motor coordination. Correlations with the H&Y scores followed a similar pattern, with the fw-index showing positive correlations and fw-FA showing negative correlations in both OFF and ON states, with a reduction in correlation volumes in the ON state.
FWF showed significant positive correlations with MDS-UPDRS-III and H&Y scores in both states, indicating that increased extracellular water content is associated with greater motor impairment. This aligns with the hypothesis that neurodegenerative processes, including neuroinflammation and axonal loss, contribute to worsening PD symptoms. However, these correlations were primarily observed in the OFF state, with a dramatic decrease in significant correlation volume in the ON state. This could indicate that these biomarkers primarily reflect the more chronic state, independent of the temporary improvements induced by levodopa. Orientation dispersion index (ODI) also exhibited positive correlations with both clinical scores, further supporting the idea that increased neurite complexity is associated with disease severity. This could reflect compensatory structural reorganization or maladaptive changes in response to neurodegeneration. In contrast, neurite density index (NDI) was negatively correlated with MDS-UPDRS III and H&Y scores, with lower NDI values linked to more severe motor impairment. This suggests that reduced neurite density, indicative of axonal degeneration, is a key factor in PD progression.
These findings highlight the potential of fw-DTI and NODDI metrics as biomarkers for tracking disease severity and treatment response in PD. The distinct correlation patterns of FWF, ODI, and NDI with clinical scores suggest that multiple microstructural processes contribute to PD pathology, emphasizing the need for multimodal imaging approaches to fully understand the disease’s impact on white matter integrity.
Several studies have examined the relationship between diffusion imaging metrics—DTI, fw-DTI, and NODDI—and motor symptom severity in PD as assessed by the MDS-UPDRS-III. DTI-derived measures, such as FA and mean diffusivity (MD), have been linked to MDS-UPDRS-III scores, suggesting that white matter microstructural alterations are associated with disease severity (58). Fw-DTI studies have shown that increased free-water content in the posterior substantia nigra correlates with higher MDS-UPDRS-III scores, indicating that extracellular water accumulation may contribute to motor impairment (14, 55). NODDI-based analyses have reported significant correlations between orientation dispersion index (ODI) and neurite density index (NDI) with MDS-UPDRS-III, reflecting the role of neurite complexity and axonal degeneration in PD-related motor deficits (26).
The bootstrap analysis of AUC distributions provides important information about the discriminative power of fw-DTI and NODDI metrics in distinguishing between groups. As illustrated in Figure 7, fw-FA exhibited the highest classification performance with an AUC of 0.956, indicating its strong sensitivity in detecting microstructural differences. This was followed by ODI (AUC = 0.935), FWF (AUC = 0.926), and fw-index (AUC = 0.907), all demonstrating robust predictive capability. In contrast, NDI showed the lowest AUC (0.872), suggesting relatively weaker discriminative power in this context. These results indicate that fw-FA and related diffusion metrics have strong potential for clinical use in distinguishing PD patients from healthy individuals. Although currently limited to research settings, these advanced imaging metrics may eventually aid in diagnostic support, particularly if future studies confirm their utility in very early-stage patients and in cases where the diagnosis is uncertain.
Pairwise comparisons of AUC values further highlight differences in classification performance among the metrics. The largest observed difference was between fw-FA and NDI (<ΔAUC> = 0.085), suggesting that free-water-corrected fractional anisotropy provides superior group differentiation compared to neurite density index. However, the CI for all pairwise comparisons include zero, indicating that the observed AUC differences are not statistically significant.
Although the current study did not directly examine gray matter structures such as the SN, previous diffusion MRI studies have reported consistent microstructural alterations in the SN of individuals with PD, including increased FW content and reduced FA (14, 55). These changes are thought to reflect neurodegeneration and gliosis in this key dopaminergic region. While our analyses were restricted to white matter, the observed alterations in diffusion metrics, particularly in projection pathways and frontal tracts, may reflect downstream effects of SN pathology, consistent with the known progression of PD-related neurodegeneration along cortico-subcortical circuits.
The primary limitation of this study is the relatively small sample size, which may reduce statistical power and limit the generalizability of the results. Additionally, the OFF scans were consistently acquired before the ON scans, we cannot fully rule out the potential influence of time-of-day effects on diffusion metrics; however, existing literature suggests such circadian variation may primarily affect CSF or gray matter rather than white matter, and reported findings remain mixed regarding the directionality of these changes (59). There are also well-known limitations of voxel-based analysis that may be overcome with advanced methods such as tract-based spatial statistics (TBSS); however, TBSS may limit sensitivity to white matter changes in the peripheral and boundary white matter regions where fw-DTI and NODDI metrics often show pathological changes. While the high AUC values (see CI of fw-FA) suggest strong classification performance, the modest cohort size raises concerns about potential overfitting, where the model may capture noise rather than biologically meaningful patterns. Replication in larger, independent cohorts is essential to validate these findings, improve the reliability of voxel-based comparisons, and enable robust subgroup analyses. Additionally, while our pipeline involved an initial upsampling of the original 2 mm DTI data to 1.25 mm using mrgrid prior to registration with the 1 mm MNI template, a step taken to optimize subsequent non-linear registration accuracy with ANTs SyN, we acknowledge the theoretical limitation of introducing potential interpolation artifacts during the upsampling process. However, we employed a high-quality interpolation method within mrgrid and visually inspected registration results to mitigate this risk. Future work could explore the impact of direct registration from the native 2 mm resolution or alternative multi-resolution registration strategies. Another limitation is the application of spatial smoothing in our data. While implemented to enhance signal detection and account for inter-subject variability, this smoothing reduced the effective spatial resolution compared to the native 2 mm in-plane acquisition. Consequently, fine anatomical details may have been blurred, and highly focal microstructural differences smaller than the smoothing kernel could have been averaged out.
5 Conclusion
This study provides new insights into PD-related white matter alterations using advanced diffusion MRI techniques, including multi-shell fw-DTI and NODDI. Our findings highlight significant microstructural differences between PD patients and HCs, particularly in free-water content, neurite density, and orientation dispersion, suggesting widespread neurodegenerative changes beyond the basal ganglia. Additionally, intra-individual comparisons between OFF and ON medication states revealed regionally constrained effects of levodopa on white matter microstructure, emphasizing the complexity of dopaminergic treatment responses. Correlations between diffusion metrics and MDS-UPDRS-III/H&Y further reinforce the role of white matter integrity in motor impairment. This study underscores the potential of fw-DTI and NODDI as sensitive biomarkers for tracking disease progression and treatment effects in PD.
Data availability statement
The datasets and analysis tools generated and used in this study are not publicly available but can be made available from the corresponding author upon reasonable request.
Ethics statement
The studies involving humans were approved by Dignity Health Institutional Review Board. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
MB: Methodology, Writing – original draft, Conceptualization, Formal analysis, Software, Validation, Data curation, Supervision, Writing – review & editing. SZ: Conceptualization, Formal analysis, Writing – review & editing, Software, Data curation, Writing – original draft. HS: Visualization, Conceptualization, Writing – review & editing. AS: Validation, Funding acquisition, Supervision, Project administration, Writing – review & editing, Visualization, Conceptualization, Methodology.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Barrow Neurological Foundation, Sam & Peggy Grossman Family Foundation, Samuel P. Mandell Foundation, and NIH (NINDS 1R01NS124575).
Acknowledgments
We would like to express our deepest gratitude to Dr. Ryan Walsh, MD, PhD, FAAN, whose contributions to this work were invaluable. Although he is no longer with us, his expertise, insight, and dedication played a significant role in shaping this study. We honor his memory and are grateful for his impact on this research and the scientific community.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2025.1605753/full#supplementary-material
Footnotes
1. ^https://www.parkinsonsmeasurement.org/toolBox/levodopaEquivalentDose.htm
References
1. Tysnes, O-B, and Storstein, A. Epidemiology of Parkinsons disease. J Neural Transm. (2017) 124:901–5. doi: 10.1007/s00702-017-1686-y
2. Dorsey, ER, and Bloem, BR. The Parkinson pandemic—a call to action. JAMA Neurol. (2018) 75:9–10. doi: 10.1001/jamaneurol.2017.3299
3. Kalia, LV, and Lang, AE. Parkinsons disease. Lancet. (2015) 386:896–912. doi: 10.1016/S0140-6736(14)61393-3
4. Jankovic, J. Parkinsons disease: clinical features and diagnosis. J Neurol Neurosurg Psychiatry. (2008) 79:368–76. doi: 10.1136/jnnp.2007.131045
5. Simuni, T, Chahine, LM, Poston, K, Brumm, M, Buracchio, T, Campbell, M, et al. A biological definition of neuronal α-synuclein disease: towards an integrated staging system for research. Lancet Neurol. (2024) 23:178–90. doi: 10.1016/S1474-4422(23)00405-2
6. Schapira, AHV, Chaudhuri, KR, and Jenner, P. Non-motor features of Parkinson disease. Nat Rev Neurosci. (2017) 18:435–50. doi: 10.1038/nrn.2017.62
7. Olanow, CW, and Stocchi, F. Levodopa: a new look at an old friend. Mov Disord. (2018) 33:859–66. doi: 10.1002/mds.27216
8. Kageyama, T, Nakamura, M, Matsuo, A, Yamasaki, Y, Takakura, Y, Hashida, M, et al. The 4F2hc/LAT1 complex transports L-DOPA across the blood-brain barrier. Brain Res. (2000) 879:115–21. doi: 10.1016/s0006-8993(00)02758-x
9. Nishijima, H, and Tomiyama, M. What mechanisms are responsible for the reuptake of levodopa-derived dopamine in parkinsonian striatum? Front Neurosci. (2016) 10:575. doi: 10.3389/fnins.2016.00575
10. Brooks, DJ. Imaging approaches to Parkinson disease. J Nucl Med. (2010) 51:596–609. doi: 10.2967/jnumed.108.059998
11. Menke, RA, Jbabdi, S, Miller, KL, Matthews, PM, and Zarei, M. Connectivity-based segmentation of the substantia nigra in human and its implications in Parkinsons disease. NeuroImage. (2010) 52:1175–80. doi: 10.1016/j.neuroimage.2010.05.086
12. Basser, PJ, Mattiello, J, and LeBihan, D. MR diffusion tensor spectroscopy and imaging. Biophys J. (1994) 66:259–67. doi: 10.1016/S0006-3495(94)80775-1
13. Bergamino, M, Keeling, EG, Mishra, VR, Stokes, AM, and Walsh, RR. Assessing white matter pathology in early-stage Parkinson disease using diffusion MRI: a systematic review. Front Neurol. (2020) 11:314. doi: 10.3389/fneur.2020.00314
14. Ofori, E, Pasternak, O, Planetta, PJ, Li, H, Burciu, RG, Snyder, AF, et al. Longitudinal changes in free-water within the substantia nigra of Parkinsons disease. Brain. (2015) 138:2322–31. doi: 10.1093/brain/awv136
15. Zhang, Y, Wu, IW, Buckley, S, Coffey, CS, Foster, E, Mendick, S, et al. Diffusion tensor imaging of the nigrostriatal fibers in Parkinsons disease. Mov Disord. (2015) 30:1229–36. doi: 10.1002/mds.26251
16. Atkinson-Clement, C, Pinto, S, Eusebio, A, and Coulon, O. Diffusion tensor imaging in Parkinsons disease: review and meta-analysis. NeuroImage Clin. (2017) 16:98–110. doi: 10.1016/j.nicl.2017.07.011
17. Chung, JW, Burciu, RG, Ofori, E, Shukla, P, Okun, MS, Hess, CW, et al. Parkinsons disease diffusion MRI is not affected by acute antiparkinsonian medication. Neuroimage Clin. (2017) 14:417–21. doi: 10.1016/j.nicl.2017.02.012
18. Bergmann, Ø, Henriques, R, Westin, CF, and Pasternak, O. Fast and accurate initialization of the free-water imaging model parameters from multi-shell diffusion MRI. NMR Biomed. (2020) 33:e4219. doi: 10.1002/nbm.4219
19. Ofori, E, Pasternak, O, Planetta, PJ, Burciu, R, Snyder, A, Febo, M, et al. Increased free water in the substantia nigra of Parkinsons disease: a single-site and multi-site study. Neurobiol Aging. (2015) 36:1097–104. doi: 10.1016/j.neurobiolaging.2014.10.029
20. Ofori, E, Krismer, F, Burciu, RG, Pasternak, O, McCracken, JL, Lewis, MM, et al. Free water improves detection of changes in the substantia nigra in Parkinsonism: a multisite study. Mov Disord. (2017) 32:1457–64. doi: 10.1002/mds.27100
21. Burciu, RG, Ofori, E, Archer, DB, Wu, SS, Pasternak, O, McFarland, NR, et al. Progression marker of Parkinsons disease: a 4-year multi-site imaging study. Brain. (2017) 140:2183–92. doi: 10.1093/brain/awx146
22. Zhang, H, Schneider, T, Wheeler-Kingshott, CA, and Alexander, DC. NODDI: practical in vivo neurite orientation dispersion and density imaging of the human brain. NeuroImage. (2012) 61:1000–16. doi: 10.1016/j.neuroimage.2012.03.072
23. Kamagata, K, Hatano, T, and Aoki, S. What is NODDI and what is its role in Parkinsons assessment? Expert Rev Neurother. (2016) 16:241–3. doi: 10.1586/14737175.2016.1142876
24. Mitchell, T, Archer, DB, Chu, WT, Coombes, SA, Lai, S, Wilkes, BJ, et al. Neurite orientation dispersion and density imaging (NODDI). and free-water imaging in Parkinsonism. Hum Brain Mapp. (2019) 40:5094–107. doi: 10.1002/hbm.24760
25. Kamagata, K, Hatano, T, Okuzumi, A, Motoi, Y, Abe, O, Shimoji, K, et al. Neurite orientation dispersion and density imaging in the substantia nigra in idiopathic Parkinson disease. Eur Radiol. (2016) 26:2567–77. doi: 10.1007/s00330-015-4066-8
26. Wei, X, Wang, S, Zhang, M, Yan, Y, Wang, Z, Wei, W, et al. Gait impairment-related axonal degeneration in Parkinsons disease by neurite orientation dispersion and density imaging. npj Parkinsons Dis. (2024) 10:45. doi: 10.1038/s41531-024-00654-w
27. Goetz, CG, Tilley, BC, Shaftman, SR, Stebbins, GT, Fahn, S, Martinez-Martin, P, et al. Movement Disorder Society-sponsored revision of the Unified Parkinsons Disease Rating Scale (MDS-UPDRS): scale presentation and clinimetric testing results. Mov Disord. (2008) 23:2129–70. doi: 10.1002/mds.22340
28. Hoehn, MM, and Yahr, MD. Parkinsonism: onset, progression and mortality. Neurology. (1967) 17:427–42. doi: 10.1212/wnl.17.5.427
29. Tomlinson, CL, Stowe, R, Patel, S, Rick, C, Gray, R, and Clarke, CE. Systematic review of levodopa dose equivalency reporting in Parkinson's disease. Mov Disord. (2010) 25:2649–53. doi: 10.1002/mds.23429
30. Tournier, J-D, Smith, R, Raffelt, D, Tabbara, R, Dhollander, T, Pietsch, M, et al. MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. NeuroImage. (2019) 202:116137. doi: 10.1016/j.neuroimage.2019.116137
31. Jenkinson, M, Beckmann, CF, Behrens, TEJ, Woolrich, MW, and Smith, SM. FSL. NeuroImage. (2012) 62:782–90. doi: 10.1016/j.neuroimage.2011.09.015
32. Veraart, J, Novikov, DS, Christiaens, D, Ades-aron, B, Sijbers, J, and Fieremans, E. Denoising of diffusion MRI using random matrix theory. NeuroImage. (2016) 142:394–406. doi: 10.1016/j.neuroimage.2016.08.016
33. Andersson, JLR, and Sotiropoulos, SN. An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. NeuroImage. (2016) 125:1063–78. doi: 10.1016/j.neuroimage.2015.10.019
34. Bastiani, M, Cottaar, M, Fitzgibbon, SP, Suri, S, Alfaro-Almagro, F, Sotiropoulos, SN, et al. Automated quality control for within and between studies diffusion MRI data using a non-parametric framework for movement and distortion correction. NeuroImage. (2019) 184:801–12. doi: 10.1016/j.neuroimage.2018.09.073
35. Andersson, JLR, Graham, MS, Zsoldos, E, and Sotiropoulos, SN. Incorporating outlier detection and replacement into a non-parametric framework for movement and distortion correction of diffusion MR images. NeuroImage. (2016) 141:556–72. doi: 10.1016/j.neuroimage.2016.06.058
36. Dhollander, T, Raffelt, D, and Connelly, A. Unsupervised 3-tissue response function estimation from single-shell or multi-shell diffusion MR data without a co-registered T1 image, ISMRM workshop on breaking the barriers of diffusion MRI, ISMRM (International Society for Magnetic Resonance in Medicine). vol. 5 (2016).
37. Avants, BB, Epstein, CL, Grossman, M, and Gee, JC. Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain. Med Image Anal. (2008) 12:26–41. doi: 10.1016/j.media.2007.06.004
38. Pasternak, O, Sochen, N, Gur, Y, Intrator, N, and Assaf, Y. Free water elimination and mapping from diffusion MRI. Magn Reson Med. (2009) 62:717–30. doi: 10.1002/mrm.22055
39. Hoy, AR, Koay, CG, Kecskemeti, SR, and Alexander, AL. Optimization of a free water elimination two-compartment model for diffusion tensor imaging. NeuroImage. (2014) 103:323–33. doi: 10.1016/j.neuroimage.2014.09.053
40. Garyfallidis, E, Brett, M, Amirbekian, B, Rokem, A, van der Walt, S, Descoteaux, M, et al. Dipy, a library for the analysis of diffusion MRI data. Front Neuroinform. (2014) 8:8. doi: 10.3389/fninf.2014.00008
41. Daducci, A, Canales-Rodríguez, EJ, Zhang, H, Dyrby, TB, Alexander, DC, and Thiran, JP. Accelerated Microstructure Imaging via Convex Optimization (AMICO). from diffusion MRI data. NeuroImage. (2015) 105:32–44. doi: 10.1016/j.neuroimage.2014.10.026
42. Wakana, S, Jiang, H, Nagae-Poetscher, LM, van Zijl, PCM, and Mori, S. Fiber tract-based atlas of human white matter anatomy. Radiology. (2004) 230:77–87. doi: 10.1148/radiol.2301021640
43. Hua, K, Zhang, J, Wakana, S, Jiang, H, Li, X, Reich, DS, et al. Tract probability maps in stereotaxic spaces: analyses of white matter anatomy and tract-specific quantification. NeuroImage. (2008) 39:336–47. doi: 10.1016/j.neuroimage.2007.07.053
44. Smith, SM, and Nichols, TE. Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. NeuroImage. (2009) 44:83–98. doi: 10.1016/j.neuroimage.2008.03.061
45. Benjamini, Y, and Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc Ser B (Methodological). (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
46. Hedges, LV. Distribution theory for Glasss estimator of effect size and related estimators. J Educ Stat. (1981) 6:107–28. doi: 10.3102/10769986006002107
47. Spearman, C. The proof and measurement of association between two things. Am J Psychol. (1904) 15:72–101. doi: 10.2307/1412159
48. Motor fluctuations and dyskinesias (diagnosis and management). Ann Indian Acad Neurol. (2011) 14:S13–5.
49. Kamiya, K, Hori, M, and Aoki, S. NODDI in clinical research. J Neurosci Methods. (2020) 346:108908. doi: 10.1016/j.jneumeth.2020.108908
50. Ingalhalikar, M, Smith, A, Parker, D, Satterthwaite, TD, Elliott, MA, Ruparel, K, et al. Sex differences in the structural connectome of the human brain. Proc Natl Acad Sci USA. (2014) 111:823–8. doi: 10.1073/pnas.1316909110
51. Bergamino, M, Keeling, EG, Baxter, LC, Sisco, NJ, Walsh, RR, and Stokes, AM. Sex differences in Alzheimer's disease revealed by free-water diffusion tensor imaging and voxel-based morphometry. J Alzheimers Dis. (2022) 85:395–414. doi: 10.3233/JAD-210406
52. Koshiyama, D, Fukunaga, M, Okada, N, Morita, K, Nemoto, K, Usui, K, et al. White matter microstructural alterations across four major psychiatric disorders: mega-analysis study in 2937 individuals. Mol Psychiatry. (2020) 25:883–95. doi: 10.1038/s41380-019-0553-7
53. Gattellaro, G, Minati, L, Grisoli, M, Mariani, C, Carella, F, Osio, M, et al. White matter involvement in idiopathic Parkinson disease: a diffusion tensor imaging study. Am J Neuroradiol. (2009) 30:1222–6. doi: 10.3174/ajnr.A1556
54. Wycoco, V, Shroff, M, Sudhakar, S, and Lee, W. White matter anatomy: what the radiologist needs to know. Neuroimaging Clin N Am. (2013) 23:197–216. doi: 10.1016/j.nic.2012.12.002
55. Planetta, PJ, Ofori, E, Pasternak, O, Burciu, RG, Shukla, P, DeSimone, JC, et al. Free-water imaging in Parkinsons disease and atypical parkinsonism. Brain. (2016) 139:495–508. doi: 10.1093/brain/awv361
56. Mole, JP, Subramanian, L, Bracht, T, Morris, H, Metzler-Baddeley, C, and Linden, DEJ. Increased fractional anisotropy in the motor tracts of Parkinsons disease suggests compensatory neuroplasticity or selective neurodegeneration. Eur Radiol. (2016) 26:3327–35. doi: 10.1007/s00330-015-4178-1
57. Lenfeldt, N, Larsson, A, Nyberg, L, Birgander, R, and Forsgren, L. Fractional anisotropy in the substantia nigra in Parkinsons disease: a complex picture. Eur J Neurol. (2015) 22:1408–14. doi: 10.1111/ene.12760
58. Zhang, Y, and Burock, MA. Diffusion tensor imaging in Parkinsons disease and parkinsonian syndrome: a systematic review. Front Neurol. (2020) 11:531993. doi: 10.3389/fneur.2020.531993
Keywords: Parkinson’s disease, levodopa, diffusion MRI, white matter alterations, NODDI, free-water DTI
Citation: Bergamino M, Zhu S, Shill HA and Stokes AM (2025) Investigating white matter alterations in Parkinson’s disease using multi-shell free-water DTI and NODDI: insights into neurodegeneration and levodopa effects. Front. Neurol. 16:1605753. doi: 10.3389/fneur.2025.1605753
Edited by:
Kyle Carson Kern, University of California, Los Angeles, United StatesReviewed by:
Luoyu Wang, Hangzhou First People’s Hospital, ChinaJacopo Pasquini, University of Pisa, Italy
Copyright © 2025 Bergamino, Zhu, Shill and Stokes. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ashley M. Stokes, YXNobGV5LnN0b2tlc0BiYXJyb3duZXVyby5vcmc=
†These authors have contributed equally to this work