Abstract
Objective:
This study aimed to develop and validate a predictive model for the decline in treatment efficacy among Parkinson’s disease (PD) patients based on clinical characteristics and biological markers, providing a basis for early risk identification and personalized therapeutic strategies.
Methods:
A retrospective study was conducted on 500 PD patients admitted to our hospital between January 2021 and December 2024. The patients were randomly divided into a training set (n = 350) and a validation set (n = 150) at a 7:3 ratio. Demographic characteristics, clinical rating scales, and biological markers were collected for all patients. In the training set, univariate analysis was performed to screen variables associated with treatment efficacy decline. After variable selection using LASSO regression, multivariate logistic regression analysis was performed to identify independent predictors. Predictive models, including random forest (RF), support vector machine (SVM), and gradient boosting, were constructed using Python 3.8.5 and the scikit-learn library. Model performance was evaluated using the area under the receiver operating characteristic curve (AUC), and the optimal model was selected based on key predictor importance.
Results:
No significant differences in baseline characteristics were observed between the training and validation sets (all p > 0.05). A multivariate logistic regression analysis identified the total MDS-UPDRS score, MDS-UPDRS II (activities of daily living), MDS-UPDRS IV (motor complications), PDQ-39 score, E. coli/Lactobacillus ratio, fecal lactoferrin, and fecal calprotectin as independent risk factors (all p < 0.05), while total fecal bacterial count was an independent protective factor (all p < 0.05). The RF model demonstrated superior predictive performance (AUC = 0.874, 95%CI: 0.831–0.917) compared to SVM (AUC = 0.806, 95%CI: 0.753–0.859) and gradient boosting (AUC = 0.842, 95%CI: 0.794–0.889).
Conclusion:
The RF model incorporating clinical and biological markers effectively predicts decline in treatment efficacy among PD patients, with fecal calprotectin, fecal lactoferrin, and the E. coli/Lactobacillus ratio serving as key predictors.
Introduction
Parkinson’s disease (PD) is a prevalent neurodegenerative disorder characterized by motor symptoms, including bradykinesia, tremor, and rigidity, as well as various non-motor symptoms that diminish quality of life and significantly impair patients’ survival (1). Current clinical management primarily relies on pharmacotherapy, including levodopa; however, the majority of patients eventually fail to maintain a stable therapeutic response and develop wearing-off phenomena, characterized by symptom fluctuations and motor complications, which pose a major challenge in disease management (2, 3). Conventional efficacy assessments, based on clinical observations and rating scales, are limited by subjectivity and delayed detection, hindering early prediction of therapeutic response trajectories.
Recent studies indicate that wearing-off in PD is closely associated with gut microbiota dysbiosis, chronic inflammatory, and the progression of non-motor symptoms (4). For instance, reduced Lactobacillus and overgrowth of Escherichia coli in feces may influence neuroinflammation and dopaminergic neuron function via the gut–brain axis, while elevated fecal inflammatory markers (e.g., lactoferrin and calprotectin) suggest intestinal barrier disruption and systemic inflammation, potentially exacerbating therapeutic decline (5, 6). However, single biomarkers exhibit limited predictive value, and constructing a robust wearing-off prediction model by integrating multidimensional indicators remains challenging.
Machine learning algorithms, with their capacity for multi-source data integration, have demonstrated significant advantages in disease prognosis and therapeutic prediction. This study aims to develop a predictive model for wearing-off in PD by combining clinical scales, quality-of-life metrics, and fecal biomarkers via machine learning, identifying key determinants to guide early intervention and therapeutic optimization.
Materials and methods
Study participants
A retrospective cohort of 500 PD patients (determined by power analysis with 80% power, α = 0.05, detecting a 15% group difference in treatment response rates) diagnosed at our neurology department between January 2021 and December 2024 was enrolled. The inclusion criteria were as follows: (1) diagnosis conforming to the Clinical Diagnostic Criteria for Parkinson’s Disease (7), (2) ≥ 6 months of levodopa-based therapy, (3) complete clinical and follow-up data, including pre- and post-treatment efficacy assessments, and (4) availability of fecal microbiota and inflammatory marker profiles. The exclusion criteria were as follows: (1) secondary parkinsonism or Parkinson-plus syndromes; (2) severe cardiac, hepatic, or renal dysfunction; (3) intestinal infections, malignancies, or other conditions affecting fecal biomarkers; and (4) lost to follow-up or incomplete data. The patients were randomly divided into a training set (n = 350) and a validation set (n = 150) at a ratio of 7:3. This study was approved by the hospital’s ethics committee. All patients signed informed consent forms.
Data collection
The following data from the electronic medical record system and the laboratory database were collected. Demographics: age and sex; Rating Scales: Movement Disorder Society-sponsored Unified Parkinson’s Disease Rating Scale (MDS-UPDRS; total score and Parts I–IV assessing non-motor symptoms, activities of daily living, motor symptoms, and motor complications, respectively), Parkinson’s Disease Fatigue Scale (PFS) (8), Non-Motor Symptoms Scale (NMSS), 39-item Parkinson’s Disease Questionnaire (PDQ-39), Wexner Constipation Scale, Geriatric Depression Scale (GDS), Parkinson’s Anxiety Scale (PAS), Lille Apathy Rating Scale (LAS), Parkinson’s Disease Sleep Scale (PDSS), Montreal Cognitive Assessment (MoCA), and Bristol Stool Form Scale.
Biological Markers: Defecation frequency (times/week), total fecal bacterial count, E. coli/Lactobacillus ratio, fecal lactoferrin, and fecal calprotectin (all fecal samples were collected within 1 week after enrollment and before standardized levodopa treatment), C-reactive protein (CRP), interleukin-6 (IL-6), plasma neurofilament light chain (NfL), serum ferritin, transferrin saturation, tumor necrosis factor-α (TNF-α), interferon-γ (IFN-γ), serum superoxide dismutase (SOD), serum dopamine, and serum glutathione (GSH).
Outcome definition
The Movement Disorder Society Consensus on Therapeutic Response in PD (2021) explicitly designated short-term (≤ 6 months) as the core window for observing therapeutic efficacy (9). The Chinese Guidelines for PD Treatment (4th edition) recommend regular follow-ups every 3–6 months to assess symptom progression, manage drug side effects, and adjust treatment plans. Thus, 6 months is a widely recognized initial screening window for therapeutic efficacy in clinical practice (10). Patients were stratified into two groups based on 6-month follow-up data:
-
Stable-response group: patients in the stable-response group met all of the following criteria at 6 months: (1) ≥ 10% reduction or ≤5% fluctuation in MDS-UPDRS Part III score versus baseline; (2) absence of motor complications (e.g., wearing-off and dyskinesias), confirmed independently by two neurologists with ≥5 years’ experience using MDS criteria; and (3) ≤ 10% increase in levodopa equivalent daily dose (LEDD).
-
Suboptimal-response group: patients in the suboptimal-response group met any of the following criteria: (1) ≥ 10% increase in MDS-UPDRS Part III score, (2) presence of motor complications (same confirmation protocol), and (3) ≥ 20% LEDD escalation to maintain symptom control.
Statistical analysis
Analyses were performed using SPSS 26.0, Python 3.8.5, and R 4.2.3. Normally distributed continuous variables were expressed as mean ± SD and compared via t-tests; non-normal data were reported as median (IQR) and analyzed using the Mann–Whitney U-tests. Categorical variables were presented as counts (%) and compared using χ2 tests. In the training set, univariate screening followed by LASSO regression and multivariate logistic regression modeling identified independent predictors, and their odds ratios (ORs) and 95% confidence intervals (CIs) were calculated. Multicollinearity was assessed via the variance inflation factor (VIF) with a threshold of VIF < 5 prior to multivariate modeling. Random forest (RF), support vector machine (SVM), and gradient boosting models were built using Python 3.8.5 (scikit-learn). The receiver operating characteristic (ROC) curve was plotted, and the area under the curve (AUC) value was calculated. A p-value of < 0.05 was considered statistically significant.
Results
Baseline characteristics of training and validation sets
A total of 500 PD patients were included. The training set (n = 350) included 198 (56.73%) stable-response patients, and the validation set (n = 150) included 84 (56%) suboptimal-response patients. There was no significant difference in the baseline characteristics between the two sets (all p > 0.05) (Table 1).
Table 1
| Variable | Training set (n = 350) | Validation set (n = 150) | t/χ2 | P |
|---|---|---|---|---|
| Age (years) | 65.23 ± 8.45 | 64.89 ± 7.98 | 0.419 | 0.675 |
| Sex (male/female) | 201/149 | 88/62 | 0.066 | 0.797 |
| Hoehn & Yahr stage (I/II/III/IV/V) | 42/135/110/55/10 | 17/58/46/23/6 | 0.088 | 0.993 |
| Medication status (On/Off) | 298/52 | 125/25 | 0.264 | 0.608 |
| LEDD (mg/day) | 523.67 ± 125.34 | 518.90 ± 119.87 | 0.395 | 0.693 |
| MDS-UPDRS (total score) | 56.32 ± 15.67 | 55.89 ± 14.98 | 0.285 | 0.776 |
| MDS-UPDRS I | 12.45 ± 4.23 | 12.12 ± 3.98 | 0.814 | 0.416 |
| MDS-UPDRS II | 18.67 ± 5.34 | 18.23 ± 4.98 | 0.861 | 0.390 |
| MDS-UPDRS III | 20.34 ± 6.78 | 19.89 ± 6.23 | 0.697 | 0.486 |
| MDS-UPDRS IV | 4.86 ± 2.12 | 4.67 ± 1.98 | 0.936 | 0.340 |
| PFS score | 15.67 ± 5.23 | 15.23 ± 4.91 | 0.878 | 0.381 |
| NMSS score | 32.45 ± 10.34 | 31.89 ± 9.87 | 0.563 | 0.574 |
| PDQ-39 score | 45.67 ± 12.31 | 44.88 ± 11.85 | 0.665 | 0.506 |
| Wexner constipation score | 8.23 ± 3.12 | 8.03 ± 3.02 | 0.663 | 0.508 |
| GDS score | 7.67 ± 2.89 | 7.45 ± 2.74 | 0.749 | 0.454 |
| PAS score | 10.34 ± 3.66 | 10.12 ± 3.45 | 0.627 | 0.521 |
| LAS score | 12.59 ± 4.18 | 12.17 ± 3.83 | 1.005 | 0.292 |
| PDSS score | 142.34 ± 25.61 | 140.77 ± 24.75 | 0.633 | 0.527 |
| MoCA score | 26.32 ± 3.25 | 26.16 ± 3.11 | 0.511 | 0.610 |
| Defecation frequency (times/week) | 4.23 ± 1.29 | 4.08 ± 1.13 | 1.235 | 0.217 |
| Total fecal bacteria (CFU/g) | 6.91 ± 1.34 | 6.78 ± 1.23 | 1.018 | 0.309 |
| E. coli/Lactobacillus ratio | 1.89 ± 0.67 | 1.82 ± 0.59 | 1.109 | 0.268 |
| CRP (mg/L) | 8.20 ± 3.17 | 8.04 ± 2.88 | 0.531 | 0.596 |
| IL-6 (pg/mL) | 15.62 ± 5.44 | 15.26 ± 5.15 | 0.689 | 0.491 |
| NfL (pg/mL) | 25.34 ± 8.67 | 24.82 ± 8.21 | 0.624 | 0.533 |
| Simplified stool form scale score | 3.67 ± 1.23 | 3.56 ± 1.12 | 0.941 | 0.347 |
| Serum ferritin (ng/mL) | 235.33 ± 68.19 | 230.75 ± 65.47 | 0.696 | 0.487 |
| Transferrin saturation (%) | 32.45 ± 8.34 | 31.90 ± 8.08 | 0.682 | 0.496 |
| TNF-α (pg/mL) | 12.38 ± 4.16 | 12.14 ± 3.93 | 0.601 | 0.548 |
| IFN-γ (pg/mL) | 9.69 ± 3.32 | 9.54 ± 3.06 | 0.434 | 0.636 |
| SOD activity (U/mL) | 85.34 ± 15.67 | 84.89 ± 15.35 | 0.296 | 0.767 |
| Fecal lactoferrin (μg/g) | 5.67 ± 2.12 | 5.55 ± 2.02 | 0.588 | 0.557 |
| Serum dopamine (ng/mL) | 125.31 ± 35.56 | 123.80 ± 34.23 | 0.440 | 0.660 |
| GSH (μmol/L) | 4.86 ± 1.34 | 4.72 ± 1.20 | 1.104 | 0.270 |
| Fecal calprotectin (μg/g) | 128.64 ± 45.43 | 125.95 ± 42.11 | 0.620 | 0.536 |
Baseline characteristics of training and validation cohorts.
Univariate analysis of factors influencing the attenuation of treatment efficacy in Parkinson’s disease
In the training set, the univariate analysis showed that statistically significant differences were observed between the stable-response group and the efficacy-attenuation group in the following indicators: MDS-UPDRS, MDS-UPDRS I, MDS-UPDRS II, MDS-UPDRS III, MDS-UPDRS IV, NMSS, PDQ-39, PDSS, MoCA, total fecal bacterial count, Escherichia coli/Lactobacillus ratio, fecal lactoferrin, and fecal calprotectin (all p < 0.05) (Table 2).
Table 2
| Variable | Stable-response group (n = 152) | Efficacy-attenuation group (n = 198) | t/χ2 | P |
|---|---|---|---|---|
| Age (years) | 65.01 ± 7.98 | 66.14 ± 8.76 | 1.243 | 0.215 |
| Sex (male/female) | 142/110 | 59/39 | 0.429 | 0.513 |
| Hoehn & Yahr stage (I/II/III/IV/V) | 17/68/40/21/4 | 25/67/70/24/6 | 5.086 | 0.166 |
| Medication status (On/Off) | 215/37 | 83/15 | 0.022 | 0.883 |
| LEDD (mg/day) | 520.33 ± 120.65 | 531.61 ± 130.29 | 0.829 | 0.408 |
| MDS-UPDRS (total score) | 53.39 ± 12.48 | 58.72 ± 13.46 | 3.789 | 0.001 |
| MDS-UPDRS I | 11.82 ± 3.89 | 12.95 ± 4.73 | 2.389 | 0.017 |
| MDS-UPDRS II | 17.85 ± 5.44 | 19.36 ± 5.81 | 2.477 | 0.014 |
| MDS-UPDRS III | 19.49 ± 5.88 | 21.03 ± 6.74 | 2.238 | 0.026 |
| MDS-UPDRS IV | 4.23 ± 2.02 | 5.38 ± 2.27 | 4.926 | 0.001 |
| PFS score | 15.24 ± 4.83 | 15.96 ± 5.32 | 1.306 | 0.193 |
| NMSS score | 30.77 ± 9.84 | 33.08 ± 11.29 | 2.005 | 0.046 |
| PDQ-39 score | 43.79 ± 11.24 | 48.30 ± 13.45 | 3.335 | 0.001 |
| Wexner constipation score | 8.05 ± 2.98 | 8.61 ± 3.41 | 1.608 | 0.109 |
| GDS score | 7.23 ± 2.79 | 7.73 ± 3.12 | 1.555 | 0.121 |
| PAS score | 9.95 ± 3.56 | 10.72 ± 4.44 | 1.749 | 0.081 |
| LAS score | 12.06 ± 4.11 | 12.68 ± 4.53 | 1.321 | 0.187 |
| PDSS score | 145.67 ± 24.34 | 138.05 ± 26.92 | 2.735 | 0.007 |
| MoCA score | 26.56 ± 3.02 | 25.89 ± 2.87 | 2.116 | 0.035 |
| Defecation frequency (times/week) | 4.34 ± 1.35 | 3.89 ± 1.22 | 1.814 | 0.071 |
| Total fecal bacteria (CFU/g) | 7.32 ± 1.58 | 6.47 ± 1.43 | 5.266 | 0.001 |
| E. coli/Lactobacillus ratio | 1.76 ± 0.52 | 2.15 ± 0.78 | 5.322 | 0.001 |
| CRP (mg/L) | 8.12 ± 2.98 | 8.56 ± 3.23 | 1.306 | 0.192 |
| IL-6 (pg/mL) | 15.09 ± 4.24 | 16.03 ± 5.30 | 1.600 | 0.111 |
| NfL (pg/mL) | 24.56 ± 7.89 | 26.11 ± 9.23 | 1.657 | 0.098 |
| Simplified stool form scale score | 3.78 ± 1.20 | 3.57 ± 1.18 | 1.638 | 0.102 |
| Serum ferritin (ng/mL) | 238.67 ± 65.34 | 229.89 ± 62.12 | 1.281 | 0.201 |
| Transferrin saturation (%) | 32.12 ± 7.98 | 31.56 ± 8.23 | 0.639 | 0.523 |
| TNF-α (pg/mL) | 11.94 ± 3.51 | 12.68 ± 4.10 | 1.770 | 0.078 |
| IFN-γ (pg/mL) | 9.63 ± 2.98 | 10.16 ± 3.45 | 1.510 | 0.134 |
| SOD activity (U/mL) | 85.56 ± 14.98 | 83.34 ± 15.67 | 1.339 | 0.182 |
| Fecal lactoferrin (μg/g) | 5.12 ± 1.97 | 6.56 ± 2.33 | 6.122 | 0.001 |
| Serum dopamine (ng/mL) | 126.35 ± 32.45 | 124.78 ± 28.67 | 0.437 | 0.663 |
| GSH (μmol/L) | 4.98 ± 1.23 | 4.76 ± 1.12 | 1.745 | 0.082 |
| Fecal calprotectin (μg/g) | 112.79 ± 35.88 | 146.78 ± 52.34 | 6.863 | 0.001 |
Univariate analysis of factors influencing the attenuation of therapeutic efficacy in Parkinson’s disease.
Multivariate logistic regression analysis of factors influencing the decline in treatment efficacy for Parkinson’s disease
The post-treatment efficacy status of patients was set as the dependent variable (1 = suboptimal-response group, 0 = stable-response group). Variables demonstrating statistical significance in the univariate analysis were incorporated into LASSO regression for variable selection (Supplementary Table 1). The screening criterion lambda.1se was applied to identify variables, as illustrated in Supplementary Figures 1, 2.
Variables with appropriate predictive capacity were subsequently included in the multivariate logistic regression analysis. The results revealed that MDS-UPDRS, MDS-UPDRS II, MDS-UPDRS IV, PDQ-39, Escherichia coli/Lactobacillus ratio, fecal lactoferrin, and fecal calprotectin were independent risk factors for the decline in treatment efficacy in Parkinson’s disease patients (all p < 0.05). Conversely, total fecal bacterial count was identified as an independent protective factor against treatment efficacy decline (all p < 0.05) (Figure 1).
Figure 1

Results of the multivariate logistic regression analysis (X1: MDS-UPDRS, X3L: MDS-UPDRS II, X5: MDS-UPDRS IV, X7: PDQ-39, X10: total fecal bacterial count X11: Escherichia coli/Lactobacillus ratio, X12: fecal lactoferrin, X13: fecal calprotectin).
Prediction performance of machine learning models in training and validation sets
The random forest model, support vector machine model, and gradient boosting model were applied to both the training and validation sets. In the training set, the AUC values were 0.874 (95%CI: 0.831–0.917) for random forest, 0.806 (95%CI: 0.753–0.859) for support vector machine, and 0.842 (95%CI: 0.794–0.889) for gradient boosting. In the validation set, the corresponding AUC values were 0.832 (95%CI: 0.774–0.890), 0.768 (95%CI: 0.703–0.833), and 0.801 (95%CI: 0.738–0.864). The random forest model achieved the highest AUC values in both sets and was selected as the optimal model (Figure 2).
Figure 2

Area under the receiver operating characteristic curve of machine learning models (A: ROC curve for training set and B: ROC curve for validation set).
Construction of the predictive model for diminished therapeutic efficacy in Parkinson’s disease
The importance scores of independent influencing factors for diminished therapeutic efficacy in Parkinson’s disease were calculated using the random forest model. The descending order of feature importance was as follows: fecal calprotectin, fecal lactoferrin, Escherichia coli/Lactobacillus ratio, PDQ-39, total fecal bacterial count, MDS-UPDRS, MDS-UPDRS IV, and MDS-UPDRS II (Figures 3, 4).
Figure 3

Trend of mean out-of-bag error rate with an increasing number of decision trees.
Figure 4

Variable importance ranking in the random forest model (note: X1: MDS-UPDRS, X3: MDS-UPDRS II, X5: MDS-UPDRS IV, X7: PDQ-39, X10: total fecal bacterial count, X11: Escherichia coli/Lactobacillus ratio, X12: fecal lactoferrin, and X13: fecal calprotectin).
Discussion
The individual variability in therapeutic efficacy and the risk of efficacy decline during PD treatment remain core challenges in clinical management. Identifying the risk of efficacy decline early through objective indicators and optimizing treatment strategies are crucial for improving patients’ quality of life. This study integrated gut microbiota features, clinical rating scales, and inflammatory markers from 500 PD patients. Variables were screened using LASSO regression, and the multivariate logistic regression analysis identified the Escherichia coli/Lactobacillus ratio, fecal lactoferrin, and fecal calprotectin as independent risk factors for efficacy decline, while total bacterial count was an independent protective factor. A random forest model demonstrated the highest predictive performance (AUC = 0.874). Fecal calprotectin, lactoferrin, and the E. coli/Lactobacillus ratio ranked highest in feature importance, suggesting that gut microbiota characteristics and related inflammatory markers may serve as core biomarkers for predicting PD treatment efficacy. These findings provide a novel theoretical basis for precise clinical evaluation and personalized intervention strategies.
In this study, gut microbiota-related indicators played a pivotal role in predicting efficacy decline, which can be mechanistically explained by the multi-pathway regulatory network of the “gut-brain axis.” This involves the synergistic effects of gut microbial dysbiosis, intestinal barrier damage, and systemic inflammatory responses. An elevated E. coli/Lactobacillus ratio was a strong risk factor for efficacy decline (OR = 3.139), which was directly linked to disrupted gut–brain axis homeostasis due to microbial imbalance. Under normal physiological conditions, Lactobacillus, as a dominant species, maintains intestinal barrier integrity by secreting short-chain fatty acids (e.g., butyrate), suppresses pro-inflammatory cytokine release, and promotes the synthesis of neurotransmitter precursors such as γ-aminobutyric acid (GABA), thereby indirectly modulating central nervous system function (11, 12). In contrast, E. coli overgrowth may exacerbate therapeutic decline through multiple pathways, including the production of lipopolysaccharides (LPS), which activate intestinal mucosal immune cells (e.g., macrophages and dendritic cells) and trigger local inflammation (13). Additionally, E. coli disrupts tight junction proteins (e.g., occludin and claudin), thereby increasing intestinal permeability and allowing LPS and pro-inflammatory factors to enter systemic circulation, inducing low-grade systemic inflammation (14). Furthermore, E. coli may affect the substantia nigra via the vagus nerve or bloodstream, exacerbating oxidative stress damage in dopaminergic neurons and reducing drug sensitivity (15). The positive correlation between this ratio and MDS-UPDRS IV (motor complications) scores further supports the notion that microbial dysbiosis may accelerate efficacy decline by promoting motor complications.
Total bacterial count emerged as an independent protective factor (OR = 0.650), with higher levels indicating sufficient microbial abundance, reflecting the positive role of gut microbiota diversity in maintaining treatment response. A diverse gut microbiota is fundamental to metabolic homeostasis, with sufficient microbial abundance conferring multiple benefits (16). It enhances the synthesis of neuroactive compounds such as short-chain fatty acids and tryptophan metabolites, strengthening the suppression of central neuroinflammation (17). Additionally, it improves colonization resistance against pathogens, reducing the risk of intestinal infections and chronic inflammation (18). Moreover, sufficient microbial abundance helps stabilize the intestinal metabolism of drugs such as levodopa, regulating bioavailability and ensuring consistent therapeutic effects (19). Previous studies have confirmed that PD patients commonly exhibit reduced gut microbiota diversity, which correlates with disease severity. Our findings further suggest a positive association between microbial abundance and drug efficacy (20).
Fecal lactoferrin (a neutrophil activation marker) and calprotectin (a gut-specific inflammatory indicator) serve as direct evidence of intestinal barrier damage and immune activation. Both were independent risk factors for efficacy decline and ranked highest in feature importance in the random forest model, indicating that local intestinal inflammation may serve as a core predictive signal. The underlying mechanisms include the following: (1) gut inflammation spreading to the central nervous system via immune pathways of the gut–brain axis, activating microglia and promoting pro-inflammatory cytokine (e.g., IL-1β and TNF-α) release in the substantia nigra, thereby accelerating dopaminergic neuron degeneration (21); (2) intestinal motility dysfunction under inflammatory conditions, which affects drug absorption stability and leads to fluctuating plasma drug concentrations, thereby reducing therapeutic efficacy (22); and (3) chronic gut inflammation exacerbating neuroendocrine stress responses via the hypothalamic–pituitary–adrenal axis, thereby worsening motor fluctuations and drug tolerance (23). Notably, the correlation between these markers and MDS-UPDRS II (activities of daily living) scores suggests that gut inflammation may indirectly reduce treatment compliance by impairing daily functioning, creating a vicious cycle.
In this study, clinical indicators such as MDS-UPDRS total score, MDS-UPDRS II/IV, and PDQ-39 scores were incorporated into the optimal model alongside gut microbiota features, reflecting the multidimensional relationship between disease severity, functional status, and gut microbial ecology. For instance, MDS-UPDRS IV scores directly reflect motor complications (e.g., wearing-off and dyskinesia), which are not only dose-dependent but may also be accelerated by neuroinflammation triggered by microbial dysbiosis. The PDQ-39, a quality-of-life metric, positively correlated with fecal calprotectin levels, suggesting that gut inflammation may indirectly reduce subjective treatment satisfaction by exacerbating non-motor symptoms (e.g., pain and anxiety) (24, 25). The synergistic effects of these clinical and microbial indicators confirm the holistic regulatory role of the gut–brain axis in PD treatment response and highlight the limitations of relying solely on traditional scales, necessitating the integration of biomarkers for improved predictive accuracy.
The random forest model achieved a significantly higher AUC (0.874) than support vector machines and gradient boosting models, owing to its superior ability to integrate multidimensional data. It effectively captures non-linear relationships (e.g., interactions between the E. coli/Lactobacillus ratio and MDS-UPDRS III scores), reduces prediction variability caused by single-marker fluctuations through bootstrap sampling and decision-tree ensembles (particularly beneficial for fecal markers susceptible to detection variability), and provides intuitive feature importance rankings, highlighting the central role of gut inflammatory markers (calprotectin and lactoferrin) in guiding clinical monitoring. The model’s translational value lies in its potential to identify high-risk patients early via non-invasive fecal testing and simplified scale assessments, supporting personalized adjustments in drug dosing, probiotic supplementation (e.g., Lactobacillus), or anti-inflammatory interventions.
However, there are some limitations to this study. First, this single-center retrospective study may have selection bias due to the inherent limitations of retrospective data collection, which may affect the external validity of the findings. Second, the lack of external validation with multi-center datasets restricts the generalizability of the established random forest model to broader populations. Third, all fecal biomarkers were collected at a single time point (pre-treatment), which may not capture the dynamic changes of gut microbiota during PD progression; future studies should consider serial sampling to improve predictive accuracy. The 6-month follow-up period is relatively short considering the long-term progressive nature of PD, limiting our ability to assess the dynamic associations between gut microbiota and treatment efficacy over 3–5 years. Fourth, the use of probiotics or macrobiotics by patients was not systematically recorded in this study, which may act as confounding variables affecting gut microbiota composition and treatment response. Furthermore, gut microbiota analysis lacked metagenomic sequencing, potentially omitting key microbial features (e.g., species-specific metabolic pathways). The study did not explore gut microbiota modulation (e.g., probiotics) as an intervention, relying on indirect mechanistic evidence. Future multicenter prospective studies integrating metagenomics and metabolomics are needed to elucidate molecular pathways linking gut microbiota to drug efficacy, and external validation of the established random forest model will be conducted using multi-center datasets to enhance its clinical applicability. Subsequently, randomized controlled trials validating targeted interventions (e.g., probiotic supplementation) are warranted.
Conclusion
Gut microbiota features are key predictors of PD treatment efficacy decline. The random forest model incorporating these markers demonstrated excellent predictive performance, with fecal calprotectin and lactoferrin serving as clinically accessible biomarkers. This study provides novel insights into optimizing PD treatment strategies from a gut–brain axis perspective, offering potential for improving therapeutic responses through targeted microbial modulation.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by Yantai Yantaishan Hospital Medical Ethics Committee (no. 2020079). 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
CL: Conceptualization, Investigation, Data curation, Methodology, Resources, Visualization, Writing – original draft, Validation, Formal analysis. TY: Investigation, Project administration, Resources, Supervision, Writing – review & editing, Methodology, Conceptualization.
Funding
The author(s) declared that financial support was not received for this work and/or its publication.
Conflict of interest
The author(s) declared that this work 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 author(s) declared that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.1686882/full#supplementary-material
References
1.
Abbasinezhad-Moud F Shirazinia M Mirzabeyki R Einafshar E Alavi M Shaban S et al . Neuroprotective effects of thymol-loaded selenium nanoparticles against 6-OHDA-induced apoptosis and oxidative stress in an in vitro Parkinson's disease model. Curr Pharm Des. (2025) 32. doi: 10.2174/0113816128380006250630103711,
2.
Alaeddin S Chatterjee A Roberts TL Steiner-Lim GZ Jensen SO Gyengesi E et al . Exploring the effects of faecal microbiota transplantation on cognitive function: a review of clinical trials. Brain Behav Immun. (2025) 48:101049. doi: 10.1016/j.bbih.2025.101049,
3.
Basaia S Sarasso E Sciancalepore F Balestrino R Musicco S Pisano S et al . Multi-Center 3D CNN for Parkinson's disease diagnosis and prognosis using clinical and T1-weighted MRI data. Neuroimage Clin. (2025) 48:103859. doi: 10.1016/j.nicl.2025.103859,
4.
Carrossa G Misenti V Faggin S Giron MC Antonini A . The small intestinal microbiota and the gut-brain Axis in Parkinson's disease: a narrative review. Biomedicine. (2025) 13:13. doi: 10.3390/biomedicines13071769,
5.
Clarence T Masse N Porras C Nm P Berretta S Haroutunian V et al . A transcriptomic roadmap of Parkinson's disease progression at single cell resolution. medRxiv. (2025). doi: 10.1101/2025.07.30.25332436
6.
De Sciscio M Bryant RV Haylock-Jacobs S Day AS Pitchers W Iansek R et al . Faecal microbiota transplant in Parkinson's disease: pilot study to establish safety & tolerability. NPJ Parkinsons Dis. (2025) 11:203. doi: 10.1038/s41531-025-01061-5,
7.
Du Y Wang L Cui Y Xu X Zhang M Li Y et al . Effect of probiotics supplementation on REM sleep behavior disorder and motor symptoms in Parkinson's disease: a pilot study. CNS Neurosci Ther. (2025) 31:e70541. doi: 10.1111/cns.70541,
8.
Ducrocq H Puisieux S Hopes L Frismand S Colnat-Coulbois S Kyheng M et al . Non-motor subtypes in candidates for subthalamic deep brain stimulation for Parkinson's disease. Parkinsonism Relat Disord. (2025) 138:107945. doi: 10.1016/j.parkreldis.2025.107945,
9.
Francisco M Grau R . Biofilm proficient Bacillus subtilis prevents neurodegeneration in Caenorhabditis elegans Parkinson's disease models via PMK-1/p38 MAPK and SKN-1/Nrf2 signaling. Sci Rep. (2025) 15:9864. doi: 10.1038/s41598-025-93737-4,
10.
Gao J Yao L Yin L Wang Q . The impact of touching-communication-relief therapy on Parkinson's disease patients: a spiritual care perspective study. BMC Neurol. (2025) 25:306. doi: 10.1186/s12883-025-04300-7,
11.
Hein ZM Arbain MFF Kumar S Mehat MZ Hamid HA Che Ramli MD et al . Intermittent fasting as a neuroprotective strategy: gut-brain Axis modulation and metabolic reprogramming in neurodegenerative disorders. Nutrients. (2025) 17:2266. doi: 10.3390/nu17142266,
12.
Huang YB Tian LL Zhu ZW Zhou K-G Lai X Peng Y-Z et al . Apigenin enhances Nrf2-induced chaperone-mediated autophagy and mitigates α-synuclein pathology: implications for Parkinson's disease therapy. Phytomedicine. (2025) 141:156652. doi: 10.1016/j.phymed.2025.156652,
13.
Li J Zhao G Yu C Qu Y Shen X Zhao Y et al . Ginsenoside re exerts neuroprotective in MPTP mice: potential links to gut microbiota and serum metabolism. Neuropharmacology. (2025):110596. doi: 10.1016/j.neuropharm.2025.110596
14.
Li X Zhou J Peng W Zhao R Sun Q Liu Z et al . Mapping the global burden of early-onset Parkinson's disease: socioeconomic and regional inequalities from the global burden of disease study 2021. Front Public Health. (2025) 13:1618533. doi: 10.3389/fpubh.2025.1618533,
15.
Lin SH Lin RJ Chu CL Chen Y-L Fu S-C . Associations between gut microbiota composition and impulse control disorders in Parkinson's disease. Int J Mol Sci. (2025) 26:6146. doi: 10.3390/ijms26136146,
16.
Ling PZ Wong KH Ho YS Cheng WY Chang RCC . The role of gut-brain Axis in modulating the impact of sterile inflammation on neuroimmune responses in neurodegenerative diseases—Alzheimer's disease and Parkinson's disease. Neuroimmunomodulation. (2025):1–22. doi: 10.1159/000547746,
17.
Nabil Y Helal MM Qutob IA Dawoud AIA Allam S Haddad R et al . Efficacy and safety of fecal microbiota transplantation in the management of parkinson's disease: a systematic review. BMC Neurol. (2025) 25:291. doi: 10.1186/s12883-025-04105-8,
18.
Ohira J Sawamura M Kawano K Sato R Taguchi T Ishimoto T et al . Impact of α-synuclein fibril structure on seeding activity in experimental models of Parkinson's disease. NPJ Parkinsons Dis. (2025) 11:224. doi: 10.1038/s41531-025-01080-2,
19.
Rajkovaca Latic I Popovic Z Mijatovic K Sahinovic I Pekic V Vucic D et al . Association of intestinal inflammation and permeability markers with clinical manifestations of Parkinson's disease. Parkinsonism Relat Disord. (2024) 123:106948. doi: 10.1016/j.parkreldis.2024.106948,
20.
Rashidi M Arima S Stetco AC Coppola C Musarò D Greco M et al . Prediction of Parkinson disease using long-term, short-term acoustic features based on machine learning. Brain Sci. (2025) 15:739. doi: 10.3390/brainsci15070739,
21.
Sopel A Szczuciński W Gosiewski T Salamon D . The role of fecal microbiota transplantation in selected neurodegenerative diseases and neurodevelopmental disorders. Prz Gastroenterol. (2025) 20:127–41. doi: 10.5114/pg.2024.146118,
22.
Stinson SE Shadrin AA Rahman Z Rødevand L Broce IJ Persson K et al . Distinct metabolic signatures of Alzheimer's and Parkinson's disease revealed through genetic overlap with metabolic markers. medRxiv. (2025). doi: 10.1101/2025.07.31.25332114,
23.
Ugwu OP Okon MB Alum EU Ugwu CN Anyanwu EG Mariam B et al . Unveiling the therapeutic potential of the gut microbiota-brain axis: novel insights and clinical applications in neurological disorders. Medicine (Baltimore). (2025) 104:e43542. doi: 10.1097/MD.0000000000043542,
24.
Villette R Novikova PV Laczny CC Mollenhauer B May P Wilmes P . Human gut microbiome gene co-expression network reveals a loss in taxonomic and functional diversity in Parkinson's disease. NPJ Biofilms Microbiomes. (2025) 11:142. doi: 10.1038/s41522-025-00780-0,
25.
Wang J Hou J Du M Hou X Wang Y Cheng H . The mechanism of alpha-lipoic acid regulating dopaminergic neuronal damage in Parkinson's disease through the Nrf2/HMOX1 pathway. Neurosci Lett. (2025) 865:138338. doi: 10.1016/j.neulet.2025.138338,
Summary
Keywords
gut microbiota, Parkinson’s disease, fecal calprotectin, fecal lactoferrin, E. coli/Lactobacillus ratio
Citation
Liu C and Yu T (2025) A gut microbiota-based predictive model for the treatment efficacy of Parkinson’s disease. Front. Neurol. 16:1686882. doi: 10.3389/fneur.2025.1686882
Received
16 August 2025
Revised
18 November 2025
Accepted
02 December 2025
Published
18 December 2025
Volume
16 - 2025
Edited by
Ying Li, Jishou University, China
Reviewed by
Ruggero Bacchin, Azienda Provinciale per i Servizi Sanitari (APSS), Italy
Kanathip Singsai, Khon Kaen University, Thailand
Updates
Copyright
© 2025 Liu and Yu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tianxia Yu, YTX7894512313465@163.com
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.