ORIGINAL RESEARCH article

Front. Neurol., 09 July 2025

Sec. Applied Neuroimaging

Volume 16 - 2025 | https://doi.org/10.3389/fneur.2025.1605753

Investigating white matter alterations in Parkinson’s disease using multi-shell free-water DTI and NODDI: insights into neurodegeneration and levodopa effects

  • 1. Barrow Neuroimaging Innovation Center, Barrow Neurological Institute, Phoenix, AZ, United States

  • 2. School of Life Sciences, Arizona State University, Tempe, AZ, United States

  • 3. The Muhammad Ali Parkinson's Center, Barrow Neurological Institute, Phoenix, AZ, United States

Abstract

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

ControlONOFF
GroupHCPDPD
Age66.88 (7.44)66.95 (8.86)
Shapiro Wilk testW = 0.96; p = 0.6889W = 0.92; p = 0.103
t-testt = −0.028; p = 0.978
# (%F)16 (43.75%)20 (25.00%)
χ2χ2 = 38.45; p < 0.001*
MDS-UPDRS-III (total)26.55 (11.87)14.30 (6.94)
Shapiro Wilk testW = 0.95; p = 0.327W = 0.93; p = 0.142
Paired t-testt = 6.717; p < 0.001*
Hoehn and Yahr Stage1.70 (0.73)1.65 (0.59)
Shapiro Wilk testW = 0.84; p = 0.004*W = 0.74; p < 0.001*
Wilcoxon Signed-Rank TestV = 4.0; p = 0.773
Motion during dMRI acquisitions
GroupHCPDPD
ABS Motion (mm)0.62 (0.28)0.52 (0.20)0.59 (0.22)
Shapiro Wilk testW = 0.81; p = 0.004*W = 0.96; p = 0.484W = 0.95; p = 0.359
Wilcoxon Rank-Sum Test - HC vs.PDon: W = 152; p = 0.811PDoff: W = 187; p = 0.399
paired t-test (PDoff vs. PDon)t = −1.566; p = 0.134
REL Motion (mm)0.28 (0.08)0.28 (0.11)0.29 (0.14)
Shapiro Wilk testW = 0.86; p = 0.018*W = 0.92; p = 0.090W = 0.86; p = 0.007*
Wilcoxon Rank-Sum Test - HC vs.PDon: W = 170; p = 0.762PDoff: W = 158; p = 0.962
Wilcoxon Signed-Rank Test (PDoff vs. PDon)V = 97; p = 0.784
Outlier (%)0.06 (0.25)0.08 (0.27)0.06 (0.28)
Shapiro Wilk testW = 0.27; p < 0.001*W = 0.236; p < 0.001*W = 0.24; p < 0.001*
Wilcoxon Rank-Sum Test - HC vs.W = 154; p = 0.715PDon: W = 161.5; p = 0.936PDoff: W = 154; p = 0.715
Wilcoxon Signed-Rank Test (PDoff vs. PDon)V = 2.0; p = 1.000
MDS-UPDRS-III (sections)ONOFFWilcoxon Signed–Rank test
Axial0.62 (0.80)0.35 (0.61)V = 313.5; p < 0.001; pcorr < 0.001
Bradykinesia1.08 (0.94)0.61 (0.78)V = 7,921; p < 0.001; pcorr < 0.001
Rigidity0.84 (0.83)0.42 (0.64)V = 991.5; p < 0.001; pcorr < 0.001
Tremor0.53 (0.75)0.26 (0.56)V = 1,400; p < 0.001; pcorr < 0.001

Demographic, clinical, and motion characteristics of healthy controls (HC) and Parkinson’s disease (PD) patients in OFF and ON states.

Results of statistical analyses, including Shapiro–Wilk tests for normality, t-tests, and Wilcoxon rank-sum and signed-rank tests, are presented for each variable. Data are shown as mean (standard deviation, SD) or percentages, with significant findings indicated (*p < 0.05). Motor symptom evaluation of Parkinson’s disease patients in OFF and ON states based on MDS-UPDRS-III sections. Mean values (standard deviations, SD) for each motor symptom are presented alongside results of the Wilcoxon signed-rank test for paired comparisons. Significant differences between OFF and ON states are indicated by the Wilcoxon Signed-Rank Test (pcorr: p corrected for multiple comparisons by Bonferroni’s).

HC, healthy control; PD, Parkinson’s disease; F, females, χ2, chi-squared test; ABS Motion, absolute motion; REL motion, relative motion.

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.

Table 2

Disease duration# subjectsMean LEDDNumber meds (avg)
<5 yrs7188.11.7
5–10 yrs12239.92.2
>10 yrs1154.03.0

Medication specifics and disease duration of the PD cohort.

LEDD, Levodopa Equivalent Daily Dose; meds, medications.

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

Table 3

(a)
fw-FAHC<PD (OFF)HC>PD (OFF)
t-testsEffect-sizeTwo-sample t-testsEffect-size
vol (%)tvol (%)gvol (%)tvol (%)g
JHU Atlas
Anterior Thalamic Radiation L0.753.5701.331.1360.40−4.0570.41−1.236
Anterior Thalamic Radiation R2.373.4551.851.190
Cortical spinal tract L0.321.077
Cortical spinal tract R0.271.102
Cingulum Hippo L1.154.1561.541.167
Forceps Major0.321.1440.27−1.143
Inferior fronto-occipital fasc L0.241.1114.37−3.2521.32−1.183
Inferior Longitudinal fasc L2.40−3.1300.70−1.137
Inferior Longitudinal fasc R0.324.0510.571.186
Superior Longitudinal fasc L1.913.6970.971.1560.33−3.3960.38−1.149
Uncinate fasc L0.321.1230.80−3.5700.50−1.246
Sup Longitudinal fasc temporal L3.283.7421.611.1380.60−3.301
ICBM81 Atlas
Middle cerebellar peduncle0.231.133
Body of corpus callosum0.363.9590.531.177
Corticospinal tract L2.631.075
Inferior cerebellar peduncle L3.821.233
Cerebral peduncle R3.343.6806.281.180
Cerebral peduncle L4.613.7567.071.167
Anterior limb of internal capsule R9.083.5325.131.204
Anterior limb of internal capsule L2.623.4004.241.121
Posterior limb of internal capsule R2.533.1780.371.044
Retrolenticular part of internal capsule L22.56−3.3454.82−1.174
Sagittal stratum L11.61−2.869
External capsule L2.29−2.824
Fornix (cres)/Stria terminalis L0.53−2.725
fw-indexHC<PD (OFF)
Two–sample t-testsEffect-size
vol (%)tvol (%)g
JHU Atlas
Anterior Thalamic Radiation L1.162.7740.521.126
Anterior Thalamic Radiation R0.534.1020.561.309
Cingulum cingulate gyrus L1.223.3041.031.118
Cingulum cingulate gyrus R17.532.9204.651.128
Cingulum Hippo R0.283.9810.251.258
Forceps Major0.383.1230.481.093
Forceps Minor1.753.2150.681.303
Inferior fronto-occipital fasc L1.922.7440.261.069
Inferior Longitudinal fasc L2.793.1281.041.086
Inferior Longitudinal fasc R0.851.085
Superior Longitudinal fasc L6.412.7790.501.098
Superior Longitudinal fasc R0.553.4960.551.185
Uncinate fasc L1.732.776
Uncinate fasc R0.303.7530.331.100
Sup Longitudinal fasc temporal L6.302.8880.531.097
ICBM81 Atlas
Middle cerebellar peduncle0.224.0390.811.128
Pontine crossing tract (a part of MCP)2.731.071
Genu of corpus callosum0.442.252
Body of corpus callosum22.842.6514.461.143
Splenium of corpus callosum1.232.7580.531.079
Medial lemniscus L11.591.089
Cerebral peduncle L0.703.5260.441.134
Anterior limb of internal capsule L0.892.433
Posterior limb of internal capsule R1.573.7201.231.125
Posterior limb of internal capsule L1.633.6080.911.119
Anterior corona radiata R2.322.677
Anterior corona radiata L7.442.6780.851.088
Superior corona radiata R11.892.5550.641.134
Superior corona radiata L9.202.404
Posterior corona radiata R0.323.118
Posterior corona radiata L3.072.9000.941.055
Posterior thalamic radiation R0.303.180
Posterior thalamic radiation L2.562.6650.281.057
External capsule L2.952.951
Cingulum (cingulate gyrus) R14.772.8293.761.089
Cingulum (cingulate gyrus) L0.473.362
Superior longitudinal fasciculus R0.303.056
Superior longitudinal fasciculus L5.982.8400.891.129
Superior fronto-occipital fasciculus L1.582.345
(b)
fw-FAPD (OFF) < PD (ON)
t-testsEffect-size
JHU Atlasvol (%)tvol (%)g
Inferior Longitudinal fasc L0.41−3.475
fw-indexPD (OFF) < PD (ON)
t-testsEffect-size
vol (%)tvol (%)g
JHU Atlas
Forceps Major0.22−3.893
Sup Longitudinal fasc temporal L0.21−3.726
ICBM81 Atlas
Splenium of corpus callosum0.57−3.904
Superior cerebellar peduncle R9.88−0.693
Posterior limb of internal capsule R0.53−3.383
Posterior thalamic radiation R2.84−0.686
Sagittal stratum R2.69−0.675
Superior longitudinal fasciculus L0.35−3.750

(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.

t, t-value; g, Hedges’ g (effect-size); vol, volume.

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 2ac, Table 4, and Supplementary Figure 2 present the voxel-based differences between HC and PD (OFF) patients for these metrics.

Figure 2

Table 4

AtlasNODDI - FWFNODDI - ODINODDI - NDI
HC<PD (OFF)HC<PD (OFF)HC>PD (OFF)HC>PD (OFF)
t-testsEffect-sizet-testsEffect-sizet-testst-testsEffect-size
vol (%)tvol (%)gvol (%)tvol (%)gvol (%)tvol (%)tvol (%)g
JHU Atlas
Anterior Thalamic Radiation L1.422.6820.793.2510.641.1960.34−3.4090.55−3.1010.09−1.156
Anterior Thalamic Radiation R0.612.4520.773.3450.271.097
Cortical spinal tract L2.252.6460.223.412
Cortical spinal tract R0.562.4790.273.7540.651.188
Cingulum cingulate gyrus R4.312.857
Forceps Major0.60−3.202
Forceps Minor4.532.8260.451.4623.013.4301.011.2220.25−3.3680.73−3.1500.13−1.168
Inferior fronto-occipital fasc L3.643.3211.831.1880.23−3.4620.45−3.150
Inferior fronto-occipital fasc R0.363.282
Inferior Longitudinal fasc L0.372.9611.313.1710.431.1720.18−3.493
Inferior Longitudinal fasc R0.161.0940.21−4.365
Superior Longitudinal fasc L0.773.0760.161.1360.393.4420.341.1450.15−1.164
Superior Longitudinal fasc R0.682.7280.393.8870.801.175
Uncinate fasc L1.723.5721.431.1720.22−3.4430.51−3.222
Sup Longitudinal fasc temporal L0.593.3310.181.1340.523.314
Sup Longitudinal fasc temporal R0.612.9370.383.6900.921.211
ICBM81 Atlas
Middle cerebellar peduncle1.662.6850.251.2011.693.508
Pontine crossing tract (a part of MCP)2.472.683
Genu of corpus callosum1.45−3.493
Body of corpus callosum4.782.6780.681.127
Splenium of corpus callosum1.48−3.172
Corticospinal tract R6.022.455
Corticospinal tract L12.042.700
Medial lemniscus L3.433.128
Inferior cerebellar peduncle L0.723.126
Anterior limb of internal capsule L0.86−3.470
Posterior limb of internal capsule L1.12−3.346
Retrolenticular part of internal capsule L14.743.5697.051.224
Anterior corona radiata R0.792.516
Anterior corona radiata L0.54−3.479
Superior corona radiata R4.472.4960.441.166
Superior corona radiata L0.592.905
Posterior thalamic radiation R1.36−4.341
Sagittal stratum L8.562.9130.631.0392.33−3.493
External capsule L3.173.2671.751.130
Cingulum (cingulate gyrus) R1.672.623

Results of the voxel-based analysis for NODDI metrics between HC and PD (OFF).

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.

t, t-value; g, Hedges’ g (effect-size); vol, volume; FWF, Free Water Fraction; ODI, Orientation Dispersion Index; NDI, Neurite Density Index; HC, healthy control; PD, Parkinson’s disease.

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 2df, 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.

Table 5

AtlasNODDI - FWFNODDI - ODINODDI - NDI
PD (ON) > PD (OFF)PD (ON) < PD (OFF)PD (ON) > PD (OFF)PD (ON) > PD (OFF)
t-testsEffect-sizet-testsEffect-sizet-testst-testsEffect-size
vol (%)tvol (%)gvol (%)tvol (%)gvol (%)tvol (%)tvol (%)g
JHU Atlas
Anterior Thalamic Radiation R0.62−0.668
Cortical spinal tract R0.41−3.6340.60−0.661
Forceps Minor0.61−0.726
Inferior fronto-occipital fasc L0.40−3.4450.72−0.732
Inferior fronto-occipital fasc R0.56−0.7070.11−4.2201.00−3.5610.92−0.812
Inferior Longitudinal fasc L0.26−0.678
Inferior Longitudinal fasc R0.19−4.1460.44−0.6910.25−4.1692.36−3.6081.89−0.724
Superior Longitudinal fasc R0.253.7980.180.712
Uncinate fasc L0.86−3.4321.21−0.740
Uncinate fasc R0.22−0.715
Sup Longitudinal fasc temporal R0.613.8810.460.724
ICBM81 Atlas
Genu of corpus callosum0.69−0.755
Superior cerebellar peduncle R9.68−0.661
Posterior limb of internal capsule R1.76−3.6521.47−0.658
Retrolenticular part of internal capsule R0.60−0.660
Retrolenticular part of internal capsule L2.27−0.690
Anterior corona radiata R0.82−0.724
Posterior thalamic radiation R0.55−4.2521.59−0.6830.93−3.3400.52−0.814
Sagittal stratum R0.76−4.0851.48−0.7092.02−4.1988.93−3.6015.33−0.631
Sagittal stratum L0.63−0.658
External capsule R0.46−3.2730.73−0.727

Results of the voxel-based analysis for NODDI metrics between PD (OFF) and PD (ON).

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.

t, t-value; g, Hedges’ g (effect-size); vol, volume; FWF, Free Water Fraction; ODI, Orientation Dispersion Index; NDI, Neurite Density Index; HC, healthy control; PD, Parkinson’s disease.

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

Figure 4

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

AtlasMDS-UPDRS IIIHoehn and Yahr
fw-FA - PD (OFF)fw-FA - PD (ON)Δ vol% PD (ON) – PD (OFF)fw-FA - PD (OFF)fw-FA - PD (ON)Δ vol% PD (ON) – PD (OFF)
vol (%)tρvol (%)tρΔvol (%)vol (%)tρvol (%)tρΔvol (%)
JHU Atlas
Anterior Thalamic Radiation L4.58−2.472−0.4660.15−2.831−0.481−4.43
Anterior Thalamic Radiation R0.36−3.321−0.568−0.366.12−2.484−0.4710.41−3.607−0.579−5.72
Cortical spinal tract L0.47−3.777−0.495−0.470.57−2.416−0.4610.93−3.781−0.5340.36
Cortical spinal tract R0.41−3.169−0.6330.05−3.862−0.657−0.372.26−3.359−0.5082.26
Cingulum cingulate gyrus L2.90−2.884−0.566−2.9011.16−2.399−0.4590.04−3.032−0.589−11.13
Cingulum cingulate gyrus R1.55−2.909−0.356−1.550.53−1.997−0.4290.57−2.580−0.4950.05
Cingulum Hippo L8.94−3.268−0.4880.71−3.760−0.598−8.233.94−2.329−0.4844.41−2.970−0.5040.47
Cingulum Hippo R8.94−3.186−0.425−8.942.62−2.735−0.5321.83−2.812−0.566−0.79
Forceps Major7.64−3.837−0.5670.43−4.262−0.763−7.2214.67−2.556−0.49918.40−3.525−0.5433.73
Forceps Minor0.96−3.382−0.549−0.966.07−2.388−0.490−6.07
Inferior fronto-occipital fasc L0.22−3.877−0.7450.2211.17−2.346−0.4822.03−3.058−0.542−9.14
Inferior fronto-occipital fasc R0.79−3.592−0.6210.14−4.536−0.769−0.6516.76−2.322−0.4728.09−2.931−0.507−8.68
Inferior Longitudinal fasc L0.28−3.726−0.7310.288.96−2.240−0.4672.60−3.026−0.534−6.36
Inferior Longitudinal fasc R1.06−3.673−0.6140.07−4.364−0.745−0.9920.60−2.521−0.49810.08−2.931−0.528−10.52
Superior Longitudinal fasc L1.18−3.745−0.5050.06−3.725−0.708−1.125.79−2.440−0.4722.45−3.146−0.509−3.33
Superior Longitudinal fasc R0.49−3.873−0.627−0.494.21−2.637−0.5071.90−2.848−0.487−2.31
Uncinate fasc L9.60−2.505−0.507−9.60
Uncinate fasc R7.03−2.176−0.452−7.03
Sup Longitudinal fasc temporal L0.64−3.518−0.5580.05−3.891−0.723−0.594.92−2.403−0.4742.42−3.352−0.504−2.50
Sup Longitudinal fasc temporal R1.52−3.909−0.632−1.5211.48−2.690−0.5145.89−2.851−0.482−5.58
ICBM81 Atlas
Middle cerebellar peduncle0.42−3.912−0.6690.422.29−3.932−0.6312.29
Genu of corpus callosum0.32−3.420−0.554−0.325.51−2.062−0.455−5.51
Body of corpus callosum1.95−3.277−0.523−1.951.15−2.411−0.4960.24−2.517−0.488−0.91
Splenium of corpus callosum24.03−3.686−0.567−24.0334.17−2.673−0.51043.38−3.465−0.5389.22
Corticospinal tract R0.44−2.762−0.5000.44
Medial lemniscus R10.14−3.815−0.59910.14
Medial lemniscus L0.29−3.389−0.6200.29
Inferior cerebellar peduncle R4.65−3.136−0.54117.87−3.781−0.67513.22
Inferior cerebellar peduncle L6.51−3.824−0.6806.51
Cerebral peduncle R11.55−3.206−0.46011.55
Cerebral peduncle L1.80−4.059−0.584−1.8
Anterior limb of internal capsule R23.93−2.701−0.478−23.93
Anterior limb of internal capsule L10.34−2.580−0.393−10.34
Posterior limb of internal capsule R1.65−2.069−0.391−1.65
Posterior limb of internal capsule L0.24−3.461−0.443−0.240.21−2.074−0.2530.64−3.519−0.4740.43
Retrolenticular part of internal capsule R18.37−2.293−0.45219.17−2.950−0.4400.80
Retrolenticular part of internal capsule L0.85−2.884−0.5051.22−3.699−0.4760.36
Anterior corona radiata R5.37−3.368−0.547−5.3726.28−2.263−0.470−26.28
Anterior corona radiata L17.28−2.544−0.509−17.28
Superior corona radiata R4.21−3.252−0.549−4.211.67−2.446−0.487−1.67
Superior corona radiata L1.04−3.477−0.440−1.041.66−2.314−0.4660.44−3.497−0.512−1.23
Posterior corona radiata R0.72−3.168−0.643−0.721.98−2.445−0.5192.44−2.522−0.4580.46
Posterior corona radiata L1.97−3.242−0.461−1.979.88−2.347−0.4735.41−3.392−0.514−4.47
Posterior thalamic radiation R1.61−3.827−0.666−1.6132.23−2.272−0.47623.79−2.800−0.534−8.43
Posterior thalamic radiation L24.38−2.164−0.4528.12−2.923−0.528−16.26
Sagittal stratum R3.50−3.472−0.618−3.557.00−2.548−0.51229.67−2.822−0.537−27.33
Sagittal stratum L18.20−2.622−0.5316.59−3.163−0.555−11.61
External capsule R13.30−2.287−0.463−13.30
External capsule L1.24−3.580−0.465−1.242.63−2.423−0.4710.16−3.394−0.420−2.47
Cingulum (cingulate gyrus) R3.20−3.014−0.551−3.25.29−2.088−0.4695.55−3.246−0.6010.26
Cingulum (cingulate gyrus) L6.43−2.964−0.589−6.436.80−2.026−0.4170.91−2.666−0.499−5.89
Cingulum (hippocampus) R4.21−2.839−0.420−4.210.89−1.803−0.4103.56−3.142−0.6182.67
Cingulum (hippocampus) L2.60−2.841−0.372−2.60.43−1.913−0.3551.21−2.575−0.4890.78
Fornix (cres)/Stria terminalis R0.71−3.283−0.452−0.7113.70−1.996−0.403−13.70
Superior longitudinal fasciculus R1.15−3.861−0.615−1.154.25−2.798−0.5181.86−2.802−0.454−2.39
Superior longitudinal fasciculus L1.53−3.606−0.542−1.535.81−2.196−0.4492.33−3.248−0.517−3.48
Uncinate fasciculus R2.37−1.896−0.396−2.37
Tapetum R50.34−2.431−0.48753.86−3.102−0.5443.52
Tapetum L39.67−2.469−0.51649.00−3.321−0.5979.33
AtlasMDS-UPDRS IIIHoehn and Yahr
fw-index - PD (off)fw-index - PD (on)Δ vol% PD (on) - PD (off)fw-index - PD (off)fw-index - PD (on)Δ vol% PD (on) - PD (off)
vol (%)tρvol (%)tρΔvol (%)vol (%)tρvol (%)tρΔvol (%)
JHU Atlas
Anterior Thalamic Radiation L8.203.0100.516−8.2020.672.3740.46911.903.0590.539−8.77
Anterior Thalamic Radiation R11.772.6510.5012.753.6360.674−9.0218.012.2550.4575.663.1350.586−12.34
Cortical spinal tract L0.382.1900.262−0.383.412.7020.5095.723.2750.6162.32
Cortical spinal tract R1.532.3130.3720.173.4470.663−1.352.572.2520.4246.683.0590.6014.11
Cingulum cingulate gyrus L25.772.6750.462−25.7717.292.2280.4500.742.7960.343−16.55
Cingulum cingulate gyrus R7.372.5650.544−7.377.572.1390.449−7.57
Cingulum Hippo R0.000.0000.0002.501.9510.443−2.50
Forceps Major2.793.5410.582−2.7910.582.3770.4501.574.4860.650−9.01
Forceps Minor14.662.6150.5211.753.2710.602−12.9118.162.2830.4562.783.1570.557−15.38
Inferior fronto-occipital fasc L8.153.1710.5550.434.1000.671−7.7220.722.4510.4821.822.9250.413−18.89
Inferior fronto-occipital fasc R16.302.7940.5092.803.4550.668−13.5022.982.2710.4500.723.6490.608−22.27
Inferior Longitudinal fasc L2.143.3980.5750.494.0910.663−1.6510.012.2890.469−10.01
Inferior Longitudinal fasc R7.133.2210.5600.043.9630.736−7.0916.342.2750.455−16.34
Superior Longitudinal fasc L3.152.8120.4790.024.7380.642−3.136.742.2480.454−6.74
Superior Longitudinal fasc R10.802.6600.4720.183.7960.583−10.626.202.1800.4450.363.2110.378−5.84
Uncinate fasc L13.343.1590.5560.000.0000.000−13.3423.592.5210.4864.322.9300.407−19.27
Uncinate fasc R20.882.7360.5225.803.5030.666−15.0836.102.3210.4680.163.5600.547−35.94
Sup Longitudinal fasc temporal L3.862.8100.469−3.867.882.1800.445−7.88
Sup Longitudinal fasc temporal R12.042.7110.4600.033.5390.513−12.0110.112.0600.4280.283.1850.355−9.83
ICBM81 Atlas
Middle cerebellar peduncle1.532.7540.5422.292.9130.4830.76
Pontine crossing tract (a part of MCP)0.732.5200.48754.403.2690.62853.67
Genu of corpus callosum30.092.5900.5405.063.1940.569−25.0323.822.1770.45810.903.1670.576−12.91
Body of corpus callosum49.812.8490.539−49.8131.542.1660.446−31.54
Splenium of corpus callosum7.852.6330.517−7.858.042.3020.4593.864.4590.649−4.18
Corticospinal tract R8.302.7270.4794.702.8060.568−3.60
Corticospinal tract L16.862.9080.5403.363.1680.638−13.50
Medial lemniscus R67.683.1010.62267.68
Medial lemniscus L2.723.1210.58668.103.0990.63465.38
Inferior cerebellar peduncle R10.232.7440.54710.23
Inferior cerebellar peduncle L1.142.6600.5718.572.6830.5837.44
Superior cerebellar peduncle R5.444.5440.7455.4426.313.4140.59726.31
Superior cerebellar peduncle L23.292.9950.58723.29
Cerebral peduncle R0.181.7990.42714.052.8900.59213.87
Cerebral peduncle L8.082.8680.57330.423.2210.59822.34
Anterior limb of internal capsule R11.662.7120.540−11.6626.772.1890.47713.773.0520.534−13.00
Anterior limb of internal capsule L3.352.3380.300−3.3528.931.9640.42218.823.1030.438−10.11
Posterior limb of internal capsule R2.643.3970.6582.646.981.8150.40810.973.0350.5054.00
Posterior limb of internal capsule L2.451.9520.4335.602.8690.4923.14
Retrolenticular part of internal capsule R5.492.9850.448−5.4913.402.0320.3832.352.8490.589−11.05
Retrolenticular part of internal capsule L8.793.3310.575−8.7922.682.1130.4400.162.5840.586−22.52
Anterior corona radiata R36.782.5880.4738.993.4710.659−27.7940.442.3030.4430.102.6610.365−40.34
Anterior corona radiata L30.012.9990.517−30.0132.792.6980.4994.312.7750.376−28.49
Superior corona radiata R38.952.9840.450−38.9527.362.2370.423−27.36
Superior corona radiata L34.122.6360.385−34.1232.782.1520.387−32.78
Posterior corona radiata R12.202.4930.423−12.205.281.8840.357−5.28
Posterior corona radiata L3.393.2430.578−3.396.522.3740.478−6.52
Posterior thalamic radiation R4.413.1370.537−4.4125.332.1440.406−25.33
Posterior thalamic radiation L8.623.3950.5730.404.0480.706−8.2218.072.4970.472−18.07
Sagittal stratum R20.783.0970.586−20.7831.152.4100.4980.853.5540.718−30.30
Sagittal stratum L3.723.7190.684−3.7226.452.4900.518−26.45
External capsule R26.522.6830.4570.573.9000.725−25.9526.912.2310.4565.933.3550.482−20.98
External capsule L0.553.0020.551−0.5522.592.2010.4440.452.6640.357−22.14
Cingulum (cingulate gyrus) R9.742.4060.523−9.7412.342.0780.449−12.34
Cingulum (cingulate gyrus) L34.102.7810.477−34.1022.792.2410.460−22.79
Fornix (cres)/Stria terminalis R12.632.2110.4160.183.9320.637−12.46
Fornix (cres)/Stria terminalis L10.311.8970.41113.073.1580.6502.76
Superior longitudinal fasciculus R19.982.6980.4590.093.6700.521−19.8912.371.9750.422−12.37
Superior longitudinal fasciculus L6.462.7480.492−6.4612.902.3010.460−12.90
Superior fronto-occipital fasciculus R31.362.4170.479−31.3646.942.0670.457−46.94
Superior fronto-occipital fasciculus L27.222.7170.422−27.2234.522.1110.411−34.52
Uncinate fasciculus R38.163.4250.609−38.1672.632.4180.500−72.63
Uncinate fasciculus L6.121.8440.395−6.12
Tapetum R10.073.4380.628−10.0711.742.3180.465−11.74
Tapetum L12.173.8280.631−12.178.332.8700.508−8.33

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.

Results are expressed as t-values and Spearman’s rank correlation coefficient (ρ). The percentage (%) indicates the cluster volume within the respective white matter region based on the atlas. The column ‘Δ vol% PD (ON) - PD (OFF)’ represents the difference in the percentage of significant cluster volumes between the ON and OFF states.

t, t-value; g, Hedges’ g (effect-size); vol, volume; FWF, Free Water Fraction; ODI, Orientation Dispersion Index; NDI, Neurite Density Index; HC, healthy control; PD, Parkinson’s disease.

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

Figure 6

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

AtlasMDS-UPDRS IIIHoehn and Yahr
FWF - PD (OFF)FWF - PD (ON)Δ vol% PD (ON)- PD (OFF)FWF - PD (OFF)FWF - PD (ON)Δ vol% PD (ON) – PD (OFFf)
vol (%)tρvol (%)tρΔvol (%)vol (%)tρvol (%)tρΔvol (%)
JHU Atlas
Anterior Thalamic Radiation L4.283.0440.5270.673.3600.607−3.6115.362.3000.44419.882.9380.5484.51
Anterior Thalamic Radiation R7.972.6150.5114.323.8120.621−3.6514.022.2780.44612.442.8020.555−1.58
Cortical spinal tract L1.123.2700.628−1.125.052.2890.48911.332.8070.5326.29
Cortical spinal tract R3.062.3530.47613.312.8690.56310.25
Cingulum cingulate gyrus L3.642.6310.5090.054.3010.657−3.595.022.3110.451−5.02
Cingulum cingulate gyrus R1.192.6080.5840.264.5050.666−0.931.022.4130.477−1.02
Forceps Major1.174.2190.6510.044.2620.683−1.133.382.6180.4830.944.7150.622−2.44
Forceps Minor7.852.6940.5573.013.2990.575−4.839.292.2410.433−9.29
Inferior fronto-occipital fasc L5.373.0370.5351.193.6630.630−4.188.982.3990.4460.572.5180.375−8.42
Inferior fronto-occipital fasc R8.412.6780.5195.403.7650.624−3.029.842.2980.4400.032.5580.649−9.81
Inferior Longitudinal fasc L0.313.8760.6830.614.1930.6300.301.422.2970.441−1.42
Inferior Longitudinal fasc R0.304.2320.7294.9770.801−0.291.222.6580.491−1.22
Superior Longitudinal fasc L0.932.5640.5040.064.6020.582−0.871.612.5160.512−1.61
Superior Longitudinal fasc R3.882.7930.5450.293.8140.550−3.580.902.6070.507−0.90
Uncinate fasc L10.702.9410.5171.753.3390.597−8.9514.772.3660.4391.522.5410.372−13.25
Uncinate fasc R10.612.5860.5199.903.7340.602−0.7016.532.3200.451−16.53
Sup Longitudinal fasc temporal L0.952.5920.501−0.951.142.4470.504−1.14
Sup Longitudinal fasc temporal R4.612.8070.5420.143.7160.508−4.460.392.3790.458−0.39
ICBM81 Atlas
Middle cerebellar peduncle4.142.2620.4871.902.4320.449−2.24
Pontine crossing tract (a part of MCP)2.003.3240.642−2.001.132.1780.46491.002.8030.57389.87
Genu of corpus callosum14.942.6460.5675.173.0980.566−9.768.112.0890.441−8.11
Body of corpus callosum36.452.8370.5560.294.2320.627−36.157.772.3380.500−7.77
Splenium of corpus callosum1.602.3340.508−1.600.982.2090.4372.294.6920.6231.31
Fornix (column and body of fornix)19.272.0050.416−19.27
Corticospinal tract R7.862.1660.48010.722.7620.5582.86
Corticospinal tract L3.283.4300.691−3.2826.722.2270.48410.002.6080.568−16.72
Medial lemniscus R83.622.7570.57083.62
Medial lemniscus L8.872.5070.49572.102.8540.58863.23
Inferior cerebellar peduncle R0.102.1780.4899.302.6410.5279.19
Inferior cerebellar peduncle L7.232.0580.46919.832.4980.51212.60
Superior cerebellar peduncle R34.172.9710.57234.17
Superior cerebellar peduncle L27.722.8400.53227.72
Cerebral peduncle R11.062.9270.51158.912.8640.54247.85
Cerebral peduncle L7.332.4200.48356.852.9250.54349.52
Anterior limb of internal capsule R2.332.7700.557−2.3310.292.1490.4733.663.5010.570−6.63
Anterior limb of internal capsule L0.302.3450.466−0.308.081.9910.43326.812.8110.46918.72
Posterior limb of internal capsule R9.032.2240.4819.032.7910.506
Posterior limb of internal capsule L3.522.3810.49515.882.7200.43612.37
Retrolenticular part of internal capsule R4.102.5170.5664.10
Retrolenticular part of internal capsule L3.124.1080.684−3.128.552.5560.4811.582.6700.491−6.97
Anterior corona radiata R22.732.5640.47913.993.8090.601−8.7520.602.2800.410−20.60
Anterior corona radiata L18.492.9640.5150.963.2780.567−17.5319.962.5520.4420.762.5970.312−19.21
Superior corona radiata R17.952.7930.433−17.952.602.2390.4690.032.2730.605−2.57
Superior corona radiata L11.312.4300.402−11.315.692.1370.398−5.69
Posterior corona radiata R1.694.4070.6790.594.1210.530−1.100.892.6990.4740.352.7230.660−0.54
Posterior corona radiata L0.032.3770.5770.03
Posterior thalamic radiation L1.083.8940.7040.733.9990.650−0.35
Sagittal stratum L2.061.8630.363−2.06
External capsule R4.872.6940.506−4.878.632.3250.476−8.63
External capsule L0.302.4720.573−0.309.842.0300.4082.132.5170.387−7.71
Cingulum (cingulate gyrus) R0.342.2760.5640.344.4530.6591.372.6790.510−1.37
Cingulum (cingulate gyrus) L3.892.6310.547−3.893.242.5770.517−3.24
Fornix (cres)/Stria terminalis R3.202.6350.455−3.20
Fornix (cres)/Stria terminalis L12.622.1020.39821.873.0140.6349.24
Superior longitudinal fasciculus R6.332.7960.5330.363.7110.515−5.960.032.2140.501−0.03
Superior longitudinal fasciculus L1.712.4370.499−1.711.852.3270.493−1.85
Superior fronto-occipital fasciculus R8.482.3890.513−8.48
Superior fronto-occipital fasciculus L12.622.5260.473−12.621.382.0680.394−1.38
Uncinate fasciculus R17.892.6270.536−17.89
Uncinate fasciculus L1.601.8390.406−1.60
AtlasMDS-UPDRS IIIHoehn and Yahr
ODI - PD (OFF)ODI - PD (ON)Δ vol% PD (ON)- PD (OFF)ODI - PD (OFF)ODI - PD (ON)Δ vol% PD (ON) – PD (OFFf)
vol (%)tρvol (%)tρΔvol (%)vol (%)tρvol (%)tρΔvol (%)
JHU Atlas
Anterior Thalamic Radiation L0.243.1980.6390.642−0.244.172.5320.50410.862.9510.5926.69
Anterior Thalamic Radiation R0.203.5010.6870.942.8510.6290.735.012.4110.5038.122.8500.5653.11
Cortical spinal tract L0.573.0260.6310.742.9500.6580.172.112.4650.5235.783.2550.6053.67
Cortical spinal tract R0.022.9950.6381.352.7540.6181.330.962.3010.50511.292.9370.58110.33
Cingulum cingulate gyrus L0.553.1310.628−0.553.242.9720.552−3.24
Cingulum cingulate gyrus R0.703.1070.431−0.700.222.3690.489−0.22
Cingulum Hippo L4.803.3960.576−4.801.182.5100.53110.832.8870.6059.65
Cingulum Hippo R1.933.6460.648−1.931.612.6490.54925.222.8670.62323.61
Forceps Major7.483.8550.6390.344.3860.744−7.1411.392.6300.51013.643.6400.6142.25
Forceps Minor0.113.5970.712−0.115.152.4970.494−5.15
Inferior fronto-occipital fasc L0.114.5470.7270.113.342.5650.5190.422.7120.562−2.91
Inferior fronto-occipital fasc R0.433.7520.6830.374.1150.694−0.0510.192.3770.474−10.19
Inferior Longitudinal fasc L0.304.3700.6800.300.567
Inferior Longitudinal fasc R0.504.0630.6510.504.0380.67910.662.5430.4830.433.7120.647−10.23
Superior Longitudinal fasc L0.134.2760.6380.183.9370.6420.052.622.6730.534−2.62
Superior Longitudinal fasc R0.173.9570.5970.174.2550.606−0.012.882.6990.503−2.88
Uncinate fasc L7.882.5860.515−7.88
Uncinate fasc R0.523.5160.706−0.526.522.4970.519−6.52
Sup Longitudinal fasc temporal L0.024.0440.6690.114.2990.6280.100.812.6360.508−0.81
Sup Longitudinal fasc temporal R0.583.9570.5980.353.9000.573−0.237.912.6680.494−7.91
ICBM81 Atlas
Middle cerebellar peduncle0.583.3030.6132.803.1000.6342.224.272.3220.5027.513.2380.5013.24
Pontine crossing tract (a part of MCP)0.202.9400.6501.872.8770.6651.670.802.1940.48977.333.0800.62576.53
Genu of corpus callosum0.182.1910.470−0.18
Body of corpus callosum1.823.2910.558−1.820.143.0410.4880.322.7040.6250.18
Splenium of corpus callosum18.343.7970.6390.314.1210.729−18.0420.332.7270.54036.743.5110.60716.40
Corticospinal tract R10.792.7180.61210.792.062.1500.4892.722.9320.5830.66
Corticospinal tract L0.952.9660.6263.213.0110.6762.267.152.4040.5164.313.1490.582−2.85
Medial lemniscus R4.062.9740.6364.0677.973.3040.62177.97
Medial lemniscus L2.723.0520.618−2.7212.302.1670.48669.963.2730.62657.65
Inferior cerebellar peduncle R12.193.3250.58521.183.7810.6908.9910.852.2480.48031.303.4480.59020.45
Inferior cerebellar peduncle L12.913.2840.6379.403.1700.672−3.513.002.1150.45527.693.5640.66024.69
Superior cerebellar peduncle R6.352.6630.6266.3514.822.9540.51314.82
Superior cerebellar peduncle L26.612.8800.55426.61
Cerebral peduncle R42.982.7840.57242.98
Cerebral peduncle L13.483.1860.64913.48
Anterior limb of internal capsule R9.112.6330.515−9.11
Anterior limb of internal capsule L3.482.7050.527−3.48
Posterior limb of internal capsule R0.242.5590.4464.052.7980.6123.81
Posterior limb of internal capsule L0.453.4720.5980.45
Retrolenticular part of internal capsule R3.582.7170.5453.58
Retrolenticular part of internal capsule L1.663.3250.5741.66
Anterior corona radiata R1.143.4990.672−1.1425.522.4320.498−25.52
Anterior corona radiata L10.302.6610.536−10.30
Superior corona radiata R3.963.3920.651−3.96
Superior corona radiata L0.474.5030.656−0.470.592.8680.4510.373.7310.599−0.21
Posterior corona radiata R0.323.0250.603−0.320.592.5080.519−0.59
Posterior corona radiata L0.222.3020.4011.783.8340.5901.56
Posterior thalamic radiation R0.153.7840.6430.633.8720.7040.4816.312.2400.446−16.31
Sagittal stratum R0.633.8940.6690.6318.222.3160.4540.723.5140.631−17.50
Sagittal stratum L3.052.5670.5393.05
External capsule L0.144.3220.577−0.140.072.1860.4820.393.6010.5480.32
Cingulum (cingulate gyrus) R2.013.0230.583−2.012.562.3540.5122.602.6910.6080.04
Cingulum (cingulate gyrus) L2.293.2470.619−2.290.912.2650.4680.362.3420.510−0.55
Cingulum (hippocampus) R7.043.1090.596−7.040.081.9170.44333.332.8770.64033.25
Cingulum (hippocampus) L3.643.1950.612−3.642.252.5810.5347.452.8240.6035.19
Fornix (cres)/Stria terminalis L44.532.8450.59444.53
Superior longitudinal fasciculus R0.453.9700.5890.453.9340.5683.192.6880.493−3.19
Superior longitudinal fasciculus L0.094.2730.706−0.092.032.5110.494−2.03
AtlasMDS-UPDRS IIIHoehn and Yahr
NDI - PD (OFF)NDI - PD (ON)Δ vol% PD (ON)- PD (OFF)NDI - PD (OFF)NDI - PD (ON)Δ vol% PD (ON) – PD (OFFf)
vol (%)tρvol (%)tρΔvol (%)vol (%)tρvol (%)tρΔvol (%)
JHU Atlas
Anterior Thalamic Radiation L5.86−2.274−0.4446.22−2.688−0.3460.37
Anterior Thalamic Radiation R0.25−3.158−0.6420.250.62−2.042−0.4200.96−2.465−0.2660.34
Cortical spinal tract L0.39−2.776−0.5090.86−2.975−0.6520.476.53−2.219−0.4749.52−2.621−0.3782.99
Cortical spinal tract R2.81−2.965−0.6192.814.52−2.115−0.4499.76−2.549−0.3285.24
Cingulum cingulate gyrus L12.08−2.280−0.4142.41−2.366−0.388−9.67
Cingulum cingulate gyrus R0.64−4.260−0.487−0.645.17−2.308−0.4087.08−2.337−0.2321.91
Cingulum Hippo L0.41−3.675−0.476−0.418.26−3.014−0.5456.54−2.589−0.340−1.72
Cingulum Hippo R1.20−2.092−0.4446.83−2.325−0.3345.63
Forceps Major1.63−3.667−0.493−1.6313.53−2.520−0.46815.00−2.884−0.4311.48
Forceps Minor2.25−2.197−0.4474.76−2.430−0.3682.51
Inferior fronto-occipital fasc L0.46−3.340−0.606−0.4615.20−2.394−0.45820.49−2.849−0.3535.29
Inferior fronto-occipital fasc R16.69−2.377−0.46112.13−2.838−0.361−4.56
Inferior Longitudinal fasc L0.39−3.383−0.622−0.3913.79−2.364−0.46721.17−2.762−0.3567.38
Inferior Longitudinal fasc R15.74−2.212−0.44011.97−2.609−0.331−3.77
Superior Longitudinal fasc L4.46−3.057−0.458−4.4612.96−2.289−0.45618.74−2.740−0.3475.77
Superior Longitudinal fasc R9.15−2.202−0.43412.42−2.666−0.3903.28
Uncinate fasc L0.10−3.028−0.441−0.109.16−2.306−0.42911.52−2.676−0.2972.36
Uncinate fasc R0.17−2.618−0.5483.64−2.613−0.3493.47
Sup Longitudinal fasc temporal L4.15−3.048−0.442−4.1516.38−2.219−0.44224.70−2.733−0.3368.32
Sup Longitudinal fasc temporal R16.96−2.128−0.43414.26−2.578−0.397−2.70
ICBM81 Atlas
Middle cerebellar peduncle3.78−3.123−0.6153.781.43−2.631−0.546−1.43
Pontine crossing tract (a part of MCP)2.47−3.195−0.6532.470.27−2.497−0.547−0.27
Genu of corpus callosum8.12−2.207−0.45011.31−2.534−0.3823.19
Body of corpus callosum1.39−3.693−0.420−1.396.86−2.312−0.42116.43−2.766−0.3459.58
Splenium of corpus callosum5.07−3.667−0.491−5.0725.48−2.486−0.44533.42−2.808−0.3847.94
Corticospinal tract R23.42−2.948−0.61323.420.15−2.345−0.504−0.15
Corticospinal tract L0.88−2.680−0.6450.884.45−2.573−0.541−4.45
Medial lemniscus R7.68−3.252−0.6477.68
Medial lemniscus L1.43−2.702−0.6291.431.29−2.399−0.525−1.29
Inferior cerebellar peduncle R13.22−3.488−0.66613.22
Inferior cerebellar peduncle L5.17−2.937−0.6425.17
Anterior limb of internal capsule R1.69−1.916−0.3855.74−2.688−0.3334.05
Anterior limb of internal capsule L1.46−2.946−0.326−1.463.84−1.921−0.3634.08−2.379−0.1890.23
Posterior limb of internal capsule R1.23−1.786−0.3884.98−2.395−0.2513.76
Posterior limb of internal capsule L0.80−3.092−0.296−0.801.81−1.893−0.3335.41−2.375−0.1653.60
Retrolenticular part of internal capsule R34.16−2.364−0.43039.76−3.116−0.3475.61
Retrolenticular part of internal capsule L13.69−2.274−0.40722.96−2.768−0.2489.28
Anterior corona radiata R6.21−2.169−0.43910.53−2.624−0.3364.32
Anterior corona radiata L1.58−3.071−0.445−1.5822.64−2.166−0.41720.80−2.506−0.355−1.84
Superior corona radiata R3.09−3.661−0.466−3.0927.56−2.234−0.43056.44−2.894−0.33728.88
Superior corona radiata L3.77−2.994−0.401−3.7719.81−2.120−0.42835.44−2.837−0.40715.64
Posterior corona radiata R35.43−2.302−0.46833.23−2.715−0.393−2.20
Posterior corona radiata L0.24−2.904−0.592−0.2436.59−2.453−0.50636.19−2.744−0.426−0.40
Posterior thalamic radiation R50.13−2.461−0.49431.32−2.930−0.409−18.81
Posterior thalamic radiation L3.04−3.365−0.624−3.0432.88−2.678−0.53030.89−2.857−0.434−1.99
Sagittal stratum R57.27−2.184−0.43047.62−2.580−0.287−9.65
Sagittal stratum L36.89−2.401−0.46336.40−3.202−0.374−0.49
External capsule R24.42−2.041−0.40317.95−2.506−0.284−6.47
External capsule L3.99−3.033−0.363−3.9910.01−2.141−0.39519.71−2.795−0.2749.70
Cingulum (cingulate gyrus) R1.45−4.255−0.465−1.456.15−2.306−0.3918.41−2.500−0.2932.26
Cingulum (cingulate gyrus) L12.21−2.085−0.3942.11−2.365−0.314−10.11
Cingulum (hippocampus) R5.18−2.315−0.3455.18
Cingulum (hippocampus) L10.91−3.264−0.5623.90−2.210−0.302−7.01
Fornix (cres)/Stria terminalis R4.89−2.001−0.41314.50−2.502−0.3399.61
Fornix (cres)/Stria terminalis L15.91−2.316−0.44953.78−2.835−0.30137.87
Superior longitudinal fasciculus R9.19−2.096−0.40517.75−2.675−0.3978.57
Superior longitudinal fasciculus L10.20−3.089−0.445−10.2020.36−2.167−0.43436.56−2.756−0.33216.20
Superior fronto-occipital fasciculus R22.88−1.898−0.40534.91−2.443−0.19812.03
Superior fronto-occipital fasciculus L13.41−1.930−0.4311.78−2.176−0.460−11.64
Uncinate fasciculus R20.79−2.236−0.29820.79
Uncinate fasciculus L31.91−2.119−0.44128.99−2.961−0.479−2.93
Tapetum R60.40−2.804−0.55672.99−3.262−0.52012.58
Tapetum L4.00−3.247−0.660−4.0050.33−2.511−0.54057.50−3.253−0.5977.17

Voxel-based correlations between NODDI metrics and clinical scores (MDS-UPDRS III and Hoehn and Yahr) in PD groups during OFF and ON states.

Results are expressed as t-values and Spearman’s rank correlation coefficient (ρ). The percentage (%) indicates the cluster volume within the respective white matter region based on the atlas. The column ‘Δ vol% PD (on) - PD (off)’ represents the difference in the percentage of significant cluster volumes between the ON and OFF states.

t, t-value; g, Hedges’ g (effect-size); vol, volume; FWF, Free Water Fraction; ODI, Orientation Dispersion Index; NDI, Neurite Density Index; HC, healthy control; PD, Parkinson’s disease.

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

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.

Statements

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

References

  • 1.

    TysnesO-BStorsteinA. Epidemiology of Parkinsons disease. J Neural Transm. (2017) 124:9015. doi: 10.1007/s00702-017-1686-y

  • 2.

    DorseyERBloemBR. The Parkinson pandemic—a call to action. JAMA Neurol. (2018) 75:910. doi: 10.1001/jamaneurol.2017.3299

  • 3.

    KaliaLVLangAE. Parkinsons disease. Lancet. (2015) 386:896912. doi: 10.1016/S0140-6736(14)61393-3

  • 4.

    JankovicJ. Parkinsons disease: clinical features and diagnosis. J Neurol Neurosurg Psychiatry. (2008) 79:36876. doi: 10.1136/jnnp.2007.131045

  • 5.

    SimuniTChahineLMPostonKBrummMBuracchioTCampbellMet al. A biological definition of neuronal α-synuclein disease: towards an integrated staging system for research. Lancet Neurol. (2024) 23:17890. doi: 10.1016/S1474-4422(23)00405-2

  • 6.

    SchapiraAHVChaudhuriKRJennerP. Non-motor features of Parkinson disease. Nat Rev Neurosci. (2017) 18:43550. doi: 10.1038/nrn.2017.62

  • 7.

    OlanowCWStocchiF. Levodopa: a new look at an old friend. Mov Disord. (2018) 33:85966. doi: 10.1002/mds.27216

  • 8.

    KageyamaTNakamuraMMatsuoAYamasakiYTakakuraYHashidaMet al. The 4F2hc/LAT1 complex transports L-DOPA across the blood-brain barrier. Brain Res. (2000) 879:11521. doi: 10.1016/s0006-8993(00)02758-x

  • 9.

    NishijimaHTomiyamaM. 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.

    BrooksDJ. Imaging approaches to Parkinson disease. J Nucl Med. (2010) 51:596609. doi: 10.2967/jnumed.108.059998

  • 11.

    MenkeRAJbabdiSMillerKLMatthewsPMZareiM. Connectivity-based segmentation of the substantia nigra in human and its implications in Parkinsons disease. NeuroImage. (2010) 52:117580. doi: 10.1016/j.neuroimage.2010.05.086

  • 12.

    BasserPJMattielloJLeBihanD. MR diffusion tensor spectroscopy and imaging. Biophys J. (1994) 66:25967. doi: 10.1016/S0006-3495(94)80775-1

  • 13.

    BergaminoMKeelingEGMishraVRStokesAMWalshRR. 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.

    OforiEPasternakOPlanettaPJLiHBurciuRGSnyderAFet al. Longitudinal changes in free-water within the substantia nigra of Parkinsons disease. Brain. (2015) 138:232231. doi: 10.1093/brain/awv136

  • 15.

    ZhangYWuIWBuckleySCoffeyCSFosterEMendickSet al. Diffusion tensor imaging of the nigrostriatal fibers in Parkinsons disease. Mov Disord. (2015) 30:122936. doi: 10.1002/mds.26251

  • 16.

    Atkinson-ClementCPintoSEusebioACoulonO. Diffusion tensor imaging in Parkinsons disease: review and meta-analysis. NeuroImage Clin. (2017) 16:98110. doi: 10.1016/j.nicl.2017.07.011

  • 17.

    ChungJWBurciuRGOforiEShuklaPOkunMSHessCWet al. Parkinsons disease diffusion MRI is not affected by acute antiparkinsonian medication. Neuroimage Clin. (2017) 14:41721. doi: 10.1016/j.nicl.2017.02.012

  • 18.

    BergmannØHenriquesRWestinCFPasternakO. 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.

    OforiEPasternakOPlanettaPJBurciuRSnyderAFeboMet al. Increased free water in the substantia nigra of Parkinsons disease: a single-site and multi-site study. Neurobiol Aging. (2015) 36:1097104. doi: 10.1016/j.neurobiolaging.2014.10.029

  • 20.

    OforiEKrismerFBurciuRGPasternakOMcCrackenJLLewisMMet al. Free water improves detection of changes in the substantia nigra in Parkinsonism: a multisite study. Mov Disord. (2017) 32:145764. doi: 10.1002/mds.27100

  • 21.

    BurciuRGOforiEArcherDBWuSSPasternakOMcFarlandNRet al. Progression marker of Parkinsons disease: a 4-year multi-site imaging study. Brain. (2017) 140:218392. doi: 10.1093/brain/awx146

  • 22.

    ZhangHSchneiderTWheeler-KingshottCAAlexanderDC. NODDI: practical in vivo neurite orientation dispersion and density imaging of the human brain. NeuroImage. (2012) 61:100016. doi: 10.1016/j.neuroimage.2012.03.072

  • 23.

    KamagataKHatanoTAokiS. What is NODDI and what is its role in Parkinsons assessment?Expert Rev Neurother. (2016) 16:2413. doi: 10.1586/14737175.2016.1142876

  • 24.

    MitchellTArcherDBChuWTCoombesSALaiSWilkesBJet al. Neurite orientation dispersion and density imaging (NODDI). and free-water imaging in Parkinsonism. Hum Brain Mapp. (2019) 40:5094107. doi: 10.1002/hbm.24760

  • 25.

    KamagataKHatanoTOkuzumiAMotoiYAbeOShimojiKet al. Neurite orientation dispersion and density imaging in the substantia nigra in idiopathic Parkinson disease. Eur Radiol. (2016) 26:256777. doi: 10.1007/s00330-015-4066-8

  • 26.

    WeiXWangSZhangMYanYWangZWeiWet 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.

    GoetzCGTilleyBCShaftmanSRStebbinsGTFahnSMartinez-MartinPet 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:212970. doi: 10.1002/mds.22340

  • 28.

    HoehnMMYahrMD. Parkinsonism: onset, progression and mortality. Neurology. (1967) 17:42742. doi: 10.1212/wnl.17.5.427

  • 29.

    TomlinsonCLStoweRPatelSRickCGrayRClarkeCE. Systematic review of levodopa dose equivalency reporting in Parkinson's disease. Mov Disord. (2010) 25:264953. doi: 10.1002/mds.23429

  • 30.

    TournierJ-DSmithRRaffeltDTabbaraRDhollanderTPietschMet 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.

    JenkinsonMBeckmannCFBehrensTEJWoolrichMWSmithSM. FSL. NeuroImage. (2012) 62:78290. doi: 10.1016/j.neuroimage.2011.09.015

  • 32.

    VeraartJNovikovDSChristiaensDAdes-aronBSijbersJFieremansE. Denoising of diffusion MRI using random matrix theory. NeuroImage. (2016) 142:394406. doi: 10.1016/j.neuroimage.2016.08.016

  • 33.

    AnderssonJLRSotiropoulosSN. An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. NeuroImage. (2016) 125:106378. doi: 10.1016/j.neuroimage.2015.10.019

  • 34.

    BastianiMCottaarMFitzgibbonSPSuriSAlfaro-AlmagroFSotiropoulosSNet 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:80112. doi: 10.1016/j.neuroimage.2018.09.073

  • 35.

    AnderssonJLRGrahamMSZsoldosESotiropoulosSN. Incorporating outlier detection and replacement into a non-parametric framework for movement and distortion correction of diffusion MR images. NeuroImage. (2016) 141:55672. doi: 10.1016/j.neuroimage.2016.06.058

  • 36.

    DhollanderTRaffeltDConnellyA. 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.

    AvantsBBEpsteinCLGrossmanMGeeJC. Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain. Med Image Anal. (2008) 12:2641. doi: 10.1016/j.media.2007.06.004

  • 38.

    PasternakOSochenNGurYIntratorNAssafY. Free water elimination and mapping from diffusion MRI. Magn Reson Med. (2009) 62:71730. doi: 10.1002/mrm.22055

  • 39.

    HoyARKoayCGKecskemetiSRAlexanderAL. Optimization of a free water elimination two-compartment model for diffusion tensor imaging. NeuroImage. (2014) 103:32333. doi: 10.1016/j.neuroimage.2014.09.053

  • 40.

    GaryfallidisEBrettMAmirbekianBRokemAvan der WaltSDescoteauxMet al. Dipy, a library for the analysis of diffusion MRI data. Front Neuroinform. (2014) 8:8. doi: 10.3389/fninf.2014.00008

  • 41.

    DaducciACanales-RodríguezEJZhangHDyrbyTBAlexanderDCThiranJP. Accelerated Microstructure Imaging via Convex Optimization (AMICO). from diffusion MRI data. NeuroImage. (2015) 105:3244. doi: 10.1016/j.neuroimage.2014.10.026

  • 42.

    WakanaSJiangHNagae-PoetscherLMvan ZijlPCMMoriS. Fiber tract-based atlas of human white matter anatomy. Radiology. (2004) 230:7787. doi: 10.1148/radiol.2301021640

  • 43.

    HuaKZhangJWakanaSJiangHLiXReichDSet al. Tract probability maps in stereotaxic spaces: analyses of white matter anatomy and tract-specific quantification. NeuroImage. (2008) 39:33647. doi: 10.1016/j.neuroimage.2007.07.053

  • 44.

    SmithSMNicholsTE. Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. NeuroImage. (2009) 44:8398. doi: 10.1016/j.neuroimage.2008.03.061

  • 45.

    BenjaminiYHochbergY. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc Ser B (Methodological). (1995) 57:289300. doi: 10.1111/j.2517-6161.1995.tb02031.x

  • 46.

    HedgesLV. Distribution theory for Glasss estimator of effect size and related estimators. J Educ Stat. (1981) 6:10728. doi: 10.3102/10769986006002107

  • 47.

    SpearmanC. The proof and measurement of association between two things. Am J Psychol. (1904) 15:72101. doi: 10.2307/1412159

  • 48.

    Motor fluctuations and dyskinesias (diagnosis and management). Ann Indian Acad Neurol. (2011) 14:S135.

  • 49.

    KamiyaKHoriMAokiS. NODDI in clinical research. J Neurosci Methods. (2020) 346:108908. doi: 10.1016/j.jneumeth.2020.108908

  • 50.

    IngalhalikarMSmithAParkerDSatterthwaiteTDElliottMARuparelKet al. Sex differences in the structural connectome of the human brain. Proc Natl Acad Sci USA. (2014) 111:8238. doi: 10.1073/pnas.1316909110

  • 51.

    BergaminoMKeelingEGBaxterLCSiscoNJWalshRRStokesAM. Sex differences in Alzheimer's disease revealed by free-water diffusion tensor imaging and voxel-based morphometry. J Alzheimers Dis. (2022) 85:395414. doi: 10.3233/JAD-210406

  • 52.

    KoshiyamaDFukunagaMOkadaNMoritaKNemotoKUsuiKet al. White matter microstructural alterations across four major psychiatric disorders: mega-analysis study in 2937 individuals. Mol Psychiatry. (2020) 25:88395. doi: 10.1038/s41380-019-0553-7

  • 53.

    GattellaroGMinatiLGrisoliMMarianiCCarellaFOsioMet al. White matter involvement in idiopathic Parkinson disease: a diffusion tensor imaging study. Am J Neuroradiol. (2009) 30:12226. doi: 10.3174/ajnr.A1556

  • 54.

    WycocoVShroffMSudhakarSLeeW. White matter anatomy: what the radiologist needs to know. Neuroimaging Clin N Am. (2013) 23:197216. doi: 10.1016/j.nic.2012.12.002

  • 55.

    PlanettaPJOforiEPasternakOBurciuRGShuklaPDeSimoneJCet al. Free-water imaging in Parkinsons disease and atypical parkinsonism. Brain. (2016) 139:495508. doi: 10.1093/brain/awv361

  • 56.

    MoleJPSubramanianLBrachtTMorrisHMetzler-BaddeleyCLindenDEJ. Increased fractional anisotropy in the motor tracts of Parkinsons disease suggests compensatory neuroplasticity or selective neurodegeneration. Eur Radiol. (2016) 26:332735. doi: 10.1007/s00330-015-4178-1

  • 57.

    LenfeldtNLarssonANybergLBirganderRForsgrenL. Fractional anisotropy in the substantia nigra in Parkinsons disease: a complex picture. Eur J Neurol. (2015) 22:140814. doi: 10.1111/ene.12760

  • 58.

    ZhangYBurockMA. Diffusion tensor imaging in Parkinsons disease and parkinsonian syndrome: a systematic review. Front Neurol. (2020) 11:531993. doi: 10.3389/fneur.2020.531993

  • 59.

    ThomasCSadeghiNNayakATreflerASarllsJBakerCIet al. Impact of time-of-day on diffusivity measures of brain tissue derived from diffusion tensor imaging. NeuroImage. (2018) 173:2534. doi: 10.1016/j.neuroimage.2018.02.026

Summary

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

Received

03 April 2025

Accepted

10 June 2025

Published

09 July 2025

Volume

16 - 2025

Edited by

Kyle Carson Kern, University of California, Los Angeles, United States

Reviewed by

Luoyu Wang, Hangzhou First People’s Hospital, China

Jacopo Pasquini, University of Pisa, Italy

Updates

Copyright

*Correspondence: Ashley M. Stokes,

†These authors have contributed equally to this work

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics