Altered gut microbial networks and metabolic pathways in multiple system atrophy: a comparative 16S rRNA study

Abstract

Introduction:

The alterations in the gut microbial network in multiple system atrophy (MSA) remain poorly understood. This study aimed to identify key gut microbial interaction networks in MSA through comprehensive multimodal analyses.

Methods:

Demographic information and frozen fecal specimens were collected from 119 participants [MSA, n = 26; Parkinson’s disease (PD), n = 66; healthy control (HC), n = 27]. Raw amplicons of the bacterial 16S rRNA V3–V4 gene region were processed using two methods: DADA2-denoising and clustering into operational taxonomic units. We conducted univariate and multivariable analyses to assess the differential abundance of bacterial genera and predicted metabolic pathways using four statistical methods: ANCOM, ANCOM-BC, ALDEx2, and MaAsLin 2. Interbacterial interactions were assessed using four correlation and two network analyses.

Results:

We consistently observed lower levels of Fusicatenibacter in MSA patients and lower levels of Butyricicoccus in PD patients compared with HCs (q < 0.05), both before and after adjusting for comorbidities, diet, and constipation status. The random forest classifiers effectively differentiated between MSA and PD, achieving high AUCs (0.75–0.78) in 5-fold cross-validation. A significant positive interbacterial interaction between Ruminococcus gnavus group and Erysipelatoclostridium was uniquely observed in MSA patients. Additionally, we identified an increase in the ARGORNPROST-PWY pathway (L-arginine degradation, q = 0.003) and a decrease in the PWY-6478 pathway (GDP-D-glycero-α-D-manno-heptose biosynthesis, q = 0.015) in MSA patients compared with HCs.

Conclusion:

Future studies are warranted to determine whether fecal microbiome-derived signatures can serve as reliable biomarkers for MSA.

1 Introduction

Multiple system atrophy (MSA) is a rare adult-onset neurodegenerative disorder characterized by glial cytoplasmic inclusions made of misfolded α-synuclein (Schrag et al., 1999; Spillantini et al., 1998). Although both MSA and Parkinson’s disease (PD) are classified as α-synucleinopathies, MSA is characterized by levodopa-unresponsive parkinsonism, autonomic dysfunction, and cerebellar ataxia (Krismer and Wenning, 2017). In contrast to PD, MSA progresses more rapidly, leading to severe disability within approximately 5 years and mortality within a decade of symptom onset (Watanabe et al., 2002; Wenning et al., 2013). While current epidemiological studies have failed to establish consistent associations between external environmental risk factors and MSA, internal environmental factors, particularly the gut microbiome, may contribute to the pathogenesis of MSA (Vanacore, 2005).

The gut microbiome plays a significant role in human health and disease through the gut-brain axis, including stimulation of immune responses and production of short-chain fatty acids (SCFAs) (Wang and Kasper, 2014). Numerous studies, including cross-sectional, longitudinal, and meta-analyses, have identified consistent alterations in the gut microbiome associated with PD, such as the enrichment of Bifidobacterium and the depletion of SCFA-producing bacteria like Lachnospiraceae and Faecalibacterium (Romano et al., 2021). Similar to PD, gut microbial alterations may contribute to the pathogenesis of MSA, as patients with MSA exhibit a proinflammatory colonic microbiome, impaired intestinal barrier integrity, and elevated levels of endotoxin-related intestinal inflammatory biomarkers (Engen et al., 2017). Additionally, both MSA and PD exhibit α-synuclein over-deposition (Brück et al., 2016). Dysbiosis of the gut microbiome, together with gut inflammation and hyperpermeability, can lead to α-synuclein seeding and propagation from the enteric nervous system to the central nervous system (CNS) via the vagus nerve (Klann et al., 2021). It is plausible that the gut microbiome contributes to MSA pathogenesis through the gut-brain axis.

In contrast to the growing body of evidence supporting the gut microbiome’s role in the pathogenesis and progression of PD (Menozzi et al., 2025; Tan et al., 2022), investigations in MSA have proven to be challenging because of its rare incidence and prevalence. Only few studies with limited sample sizes have examined the gut microbial composition in MSA. Furthermore, recent advancements in computational tools and bioinformatic methodologies, coupled with the lack of standardized protocols, can lead to inconsistent results in microbiome research (Bharti and Grimm, 2021). This study aimed to identify distinctive gut microbial signatures in MSA through comprehensive multimodal analyses. We employed two amplicon-preprocessing approaches and four statistical methods to enhance the robustness and reproducibility of our findings. Two statistical models were utilized to account for potential confounders, including comorbidities, diet, and constipation status. Through the identification of MSA-specific microbial fingerprints, interaction networks, and relevant metabolic pathways, we sought to advance the capabilities of differential diagnosis and enhance our understanding of the pathogenesis of MSA.

2 Materials and methods

2.1 Study subjects and collection of fecal specimens

This hospital-based case-control study was carried out at the National Taiwan University Cancer Center (NTUCC), Taipei, Taiwan, from August 2020 to December 2022. Patients meeting the criteria for clinically probable MSA according to the 2022 Movement Disorder Society criteria, as well as patients diagnosed with PD based on the United Kingdom Parkinson’s Disease Society Brain Bank clinical diagnostic criteria, were enrolled by two specialists in movement disorders (S.Y. Cheng, M.C. Kuo) (Hughes et al., 1992; Wenning et al., 2022). Each patient underwent a minimum follow-up period of 1 year. Healthy controls (HCs) were recruited from healthy spouses or caregivers of patients, or asymptomatic volunteers at the hospital. Individuals under 40 years old, with concurrent malignancies, cerebrovascular diseases, or known neurodegenerative diseases, were excluded from the study. All subjects completed a face-to-face interview with a structured, semi-quantitative questionnaire adapted from the previous study (Kuo et al., 2022). The Unified Parkinson’s Disease Rating Scale Part III (UPDRS III) score was assessed in PD patients. (Martínez-Martín et al., 2015).

Fecal specimens were collected with informed consent and with the approval of the National Taiwan University Hospital (NTUH) Research Ethics Committee (REC) (No. 202006060RINB). The specimens were stored in a sterile tube at −80°C until processing. All procedures involving human participants in this study were conducted in accordance with the ethical standards of the NTUH REC and the Declaration of Helsinki. The study was approved by the NTUCC REC.

2.2 DNA extraction, PCR amplification, and sequencing

Fecal DNA was extracted from 200 mg frozen samples using the QIAamp Fast DNA Stool Mini Kit (Qiagen, Hilden, Germany). The DNA samples were then stored at −20°C. The 341f forward primer and 805r reverse primer were designed to include sequences complementary to the upstream and downstream regions of the 16S rRNA V3–V4 gene segment, along with the Illumina overhang adapter sequences. The dual-index barcodes and the Illumina sequencing adapters were attached to the PCR products using the Nextera XT Index Kit (Illumina, CA, United States). The libraries were cleaned, normalized, pooled, and sequenced using the MiSeq System (Illumina, CA, United States) with V3 reagent for paired-end sequencing (2 × 300 bps).

2.3 Preprocessing of amplicons and taxonomic classification

The analyses of microbial composition were primarily conducted utilizing QIIME 2 2022.11 unless otherwise specified (Bolyen et al., 2019). Two distinct approaches were employed to preprocess the amplicons of 16S rRNA V3-V4 gene segment: denoising using DADA2 and de novo clustering into operational taxonomic units (OTUs) at a 99% identity threshold (Callahan et al., 2016; Rognes et al., 2016). The reference taxonomy was SILVA release SSU Ref NR 138.1¹ (Quast et al., 2013). Because SILVA does not curate sequences to the species level, species-level information was excluded, and genus-level assignments may mask biologically divergent species. We assembled ecologically informed prior class weights and built a bespoke taxonomy classifier (Bokulich et al., 2018; Kaehler et al., 2019).

2.4 Differential abundance analyses

Multiple algorithms designed for compositional data were used to conduct differential abundance analyses across the study groups and between subjects with and without constipation. These algorithms included: (1) ANalysis of COmposition of Microbiomes (ANCOM) (Mandal et al., 2015), (2) ANalysis of COmposition of Microbiomes with Bias Correction (ANCOM-BC) (Lin and Peddada, 2020), (3) ANOVA-Like Differential Expression 2 (ALDEx2) (Fernandes et al., 2013), and (4) Microbiome Multivariable Associations with Linear Models 2 (MaAsLin 2) (Mallick et al., 2021).

ANCOM-BC and MaAsLin 2 support multivariable differential abundance analyses. We estimated the differential abundances of genera using two statistical models. Model 1 incorporated constipation status, since it is associated with gut microbial dysbiosis and can potentially lead to spurious associations in microbiome studies of human diseases, particularly in Parkinsonism (Aho et al., 2019; Vujkovic-Cvijin et al., 2020). Model 2 incorporated potential confounders that differed significantly across the study groups, including age, sex, hypertension, constipation status, medications for constipation, and usage of probiotics.

2.5 Correlation between MSA-associated and PD-associated genera

The ranked relative differentials estimated by Songbird were visualized using Qurro (Fedarko et al., 2020; Morton et al., 2019). Log ratios were computed between the top and bottom 25% of differentially abundant genera in comparisons of MSA vs. HC and PD vs. HC. The Student’s t-test was used to determine if there was a significant difference between the disease groups and HC. To investigate whether the genera that were more abundant in MSA were also more abundant in PD, and vice versa, log ratios of a genus’s mean relative abundance in MSA vs. HC were plotted against PD vs. HC.

2.6 Using random forest classifiers to predict disease and constipation status

To assess the potential of microbial composition as a predictor of disease status, random forest classifiers were trained. Each analysis focused on either MSA or PD, along with HC, utilizing nested 5-fold cross-validation. Each genus was assigned an importance score by the scikit-learn learning estimator (Pedregosa et al., 2011). The area under the curve (AUC) and confusion matrix were used to evaluate the predictive performance of the classifiers.

2.7 Microbial correlations and networks

Four statistical methods were used to calculate the microbial correlation matrices within each study group: (1) the Spearman’s correlation on centered log-ratio (CLR) transformed data, (2) SparCC, (Friedman and Alm, 2012) and Sparse Estimation of Correlations among Microbiomes (SECOM) using (3) linear correlation and (4) distance correlation (Lin et al., 2022). To ensure sufficient representation for correlation analyses while limiting the number of multiple comparisons, genera with a prevalence ≥ 50% and an average abundance ≥ 20 within each study group were included in the correlation analyses.

The inference of microbial networks was conducted using Sparse Cooccurrence Network Investigation for Compositional data (SCNIC) and SParse InversE Covariance Estimation for Ecological Association Inference (SPIEC-EASI) (Kurtz et al., 2015; Shaffer et al., 2023). Genera with an average abundance < 2 in each study group were excluded.

2.8 Differential pathway abundance analyses

Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2) was used to predict the functional abundances of Enzyme Commission (EC) numbers and to regroup EC numbers into MetaCyc pathways (Caspi et al., 2016; Douglas et al., 2020). The outputs from PICRUSt2 were visualized using STAMP (Parks et al., 2014). Univariate and multivariable differential abundance analyses of MetaCyc pathways were conducted using four statistical methods (ANCOM, ANCOM-BC, ALDEx2, and MaAsLin 2).

2.9 Statistical analyses

The categorical variables are presented as count (percentage), while the numeric variables are presented as mean ± standard deviation (SD). The Fisher’s exact test was used to compare categorical variables across the study groups. Comparisons of age across the study groups were performed using the ANOVA with the Tukey’s post-hoc test. The Student’s t-test was used to compare age of onset and disease duration between MSA and PD. P-values < 0.05 were considered statistical significance.

The Supplementary material provided details on the microbial composition analyses. The features were collapsed at the genus level before subsequent analyses. To control false discovery rates, q values derived from the Benjamini-Hochberg procedure were used in differential abundance analyses (ANCOM, ALDEx2, and MaAsLin 2), except for ANCOM-BC, which adopted the Holm-Bonferroni method as its default approach. Q values < 0.05 were considered statistical significance. The analyses were carried out using R packages² and SAS version 9.4 (SAS Institute Inc., Cary, NC, United States) (R Core Team, 2023).

3 Results

3.1 Clinicodemographic characteristics

A total of 119 participants (MSA, n = 26; PD, n = 66; HC, n = 27) were enrolled. After quality control, 24 MSA, 52 PD, and 27 HC samples remained after DADA2-denoising, while 24 MSA, 50 PD, and 27 HC samples remained after OTU-clustering. Table 1 presents the clinicodemographic characteristics of the study subjects with qualified fecal specimens using both approaches. The PD patients were older than the MSA patients and HCs (p < 0.01). In addition, a greater proportion of MSA and PD patients reported experiencing constipation, taking medications for constipation, and using probiotics compared with HCs (p < 0.01). The PD patients displayed mild motor symptoms according to the UPDRS III score (15.2 ± 12.0). Among PD patients, those with constipation had higher UPDRS III scores (p = 0.011, Supplementary Table 1).

3.2 Alpha diversities and beta diversities

There were no significant differences in the observed number of features, the Shannon’s diversity index, or the Faith’s phylogenetic diversity across the study groups using either DADA2-denoising or OTU-clustering. The structures of the microbial composition are depicted in Supplementary Figure 1. A notable distinction in beta diversities, as indicated by unweighted and weighted UniFrac distances, was observed between MSA and HC in PERmutational Multivariate ANalysis of Variance (PERMANOVA) (q < 0.05, Supplementary Table 2; Anderson, 2001). On the other hand, not all beta diversity measurements between PD and HC, or between MSA and PD, were significantly different.

3.3 Univariate differential abundance analyses

As shown in Table 2, Fusicatenibacter consistently exhibited a significantly lower abundance in MSA than in HC across all four algorithms using both DADA2-denoising and OTU-clustering. Furthermore, Butyricicoccus was significantly less abundant in MSA than in HC when ANCOM-BC was used. On the other hand, a significant decrease in Butyricicoccus abundance was observed in PD compared with HC across all algorithms, except ALDEx2. When comparing MSA and PD, only ANCOM-BC identified Limosilactobacillus as exhibiting a significantly greater abundance in PD using DADA2-denoising (q = 0.016).

3.4 Multivariable differential abundance analyses

The differential abundances of genera were estimated after adjusting for potential confounders in two models. Model 1 included constipation status, while Model 2 included age, sex, hypertension, constipation status, medications for constipation, and usage of probiotics. Fusicatenibacter and Butyricicoccus remained significantly less abundant in MSA and PD, respectively, than in HC (p < 0.01, Supplementary Table 3). The effect size exhibited consistency with those observed in the univariate analyses.

In multivariable ANCOM-BC, Phascolarctobacterium exhibited a lower abundance in MSA than in HC in Model 1 using both DADA2-denoising and OTU-clustering (q < 0.05, Supplementary Table 4). In Model 2, both multivariable ANCOM-BC and multivariable MaAsLin 2 identified Butyricicoccus as being less abundant and Intestinibacter as being more abundant in MSA compared with HC (q < 0.05, Supplementary Tables 4, 5).

3.5 MSA-associated genera and PD-associated genera were correlated

The ranked relative differentials are illustrated in Figure 1A. Compared with those in HC, genera exhibiting greater abundances in MSA tended to be more abundant in PD, and vice versa. At the same time, genera with a lower abundance in MSA generally displayed a reduced abundance in PD, and vice versa. Figure 1B displays the log ratios between the top and bottom 25% of differentially abundant genera between MSA and HC, as well as between PD and HC, across the three study groups. The genera that differentiated MSA from HC were capable of differentiating PD from HC, and vice versa (p < 0.005). Both DADA2-denoising and OTU-clustering yielded similar results (Supplementary Figure 2). Figure 1C displays the log ratios of the mean relative abundances of genera (MSA vs. HC against PD vs. HC). The Pearson correlation coefficients of the genera identified through DADA2-denoising and OTU-clustering were both 0.66 (p < 0.001).

FIGURE 1

3.6 Using random forest classifiers to predict disease status and identify genera with high importance scores

The AUCs of MSA-HC classifier and PD-HC classifier were 0.74 and 0.78 using DADA2 denoising (Supplementary Figure 3). When using OTU-clustering, the AUCs of MSA-HC classifier and PD-HC classifier were 0.67 and 0.78. The random forest classifiers built for differentiating between MSA and PD had AUCs of 0.78 and 0.75 using DADA2-denoising and OTU-clustering, respectively. Supplementary Table 6 presents the confusion matrices, which illustrate the sensitivity and specificity of the predictive performance.

Table 3 lists the genera with importance scores ≥ 0.05 in the random forest classifiers. Fusicatenibacter had the highest importance score in MSA-HC classifiers, while Butyricicoccus had the highest importance score in PD-HC classifiers. Additionally, using both DADA2-denoising and OTU-clustering, Faecalibacterium and Blautia were found to have importance scores ≥ 0.05 in MSA-HC classifiers and PD-HC classifiers, respectively. In MSA-PD classifiers, the identified genus was Collinsella.

3.7 Microbial correlations

Table 4 presents the significant interbacterial correlations identified through the correlation analyses. Only correlations discovered using both DADA2-denoising and OTU-clustering are shown. No interbacterial correlation was found when employing the Spearman’s correlation on CLR-transformed data or SparCC after controlling for false discovery rates. On the other hand, SECOM using linear correlation and distance correlation revealed three significant correlations: (1) Ruminococcus torques group * Erysipelotrichaceae UCG-003 (in MSA), (2) Eggerthella * Anaerostipes (in HC), and (3) Butyricicoccus * Anaerostipes (in HC). The microbial correlation matrices of MSA using DADA2-denoising are displayed in Figure 2. The Spearman’s correlation and SparCC exhibited similar patterns of interbacterial correlations, while SECOM produced sparse correlation matrices. SECOM revealed an additional interbacterial correlation using distance correlation compared with linear correlation in MSA. Supplementary Figures 4-1 to 4-5 present the remaining correlation matrices.

FIGURE 2

3.8 Microbial networks

Significant interbacterial interactions were identified in the network analyses using both DADA2-denoising and OTU-clustering. In MSA, the positive interbacterial interaction of Ruminococcus gnavus group * Erysipelatoclostridium was consistently observed using both SCNIC and SPIEC-EASI (Figure 3 and Supplementary Table 7). This positive interaction was found to be unique to MSA, as it was not detected in the SCNIC or SPEC-EASI networks of PD and HC using either DADA2-denoising and OTU-clustering (Supplementary Figures 5–1, 5–2).

FIGURE 3

3.9 Differential pathway abundance analyses

Univariate differential abundance analyses of MetaCyc pathways were conducted using four statistical methods (ANCOM, ANCOM-BC, ALDEx2, and MaAsLin 2). The ARGORNPROST-PWY pathway, a metabolic pathway leading to L-arginine degradation (Caspi et al., 2017), was significantly more abundant in MSA than in HC using DADA2-denoising across all four algorithms (Figures 4A, B and Supplementary Tables 8, 9). A similar trend was observed using OTU-clustering, although the results did not reach statistical significance. On the other hand, the PWY-6478 pathway, a metabolic pathway responsible for GDP-D-glycero-α-D-manno-heptose biosynthesis (Caspi et al., 2017), was significantly less abundant in MSA than in HC in certain algorithms (Figures 4C, D and Supplementary Tables 8, 9). When comparing MSA with PD, no pathway consistently exhibited a significantly differential abundance.

FIGURE 4

The differential abundances of MetaCyc pathways were further assessed after adjusting for potential confounders in two models. A significantly increased abundance of the ARGORNPROST-PWY pathway and a decreased abundance of the PWY-6478 pathway in MSA compared with HC were still identified using multivariable ANCOM-BC and multivariable MaAsLin 2, respectively (q < 0.05, Supplementary Tables 10, 11).

3.10 Constipation status was not a confounder

Because constipation may serve as a primary confounder in this study, we examined the associations between constipation status and microbial signatures. There were no significant differences in the alpha diversities between subjects with and without constipation. PERMANOVA revealed no significant differences in the beta diversities associated with constipation status (Supplementary Figures 6–1, 6–2). Furthermore, no genera exhibited a significantly differential abundance between subjects with and without constipation using either DADA2-denoising or OTU-clustering. The sensitivity analyses, including only HCs for comparison, yielded similar negative results.

Distinct genera associated with constipation status were identified in different study groups using random forest classifiers (Supplementary Table 12). The random forest classifiers for predicting constipation status in HC had AUCs of 0.77 and 0.78 using DADA2-denoising and OTU-clustering, respectively. However, when applied to MSA or PD, these classifiers exhibited limited ability to differentiate constipation status. The AUCs ranged from 0.27 to 0.57, suggesting that constipation status was associated with a distinct gut microbial shift in the disease groups compared with HC (Supplementary Table 13). In summary, the changes in gut microbiome in MSA and PD could not be attributed to constipation status.

4 Discussion

Previous studies have failed to consistently identify the gut microbial signatures associated with MSA and did not take potential confounders into consideration. In MSA, only higher levels of Akkermansia and lower levels of Blautia and Faecalibacterium were observed in more than one study (Barichella et al., 2019; Engen et al., 2017; Wan et al., 2019). Besides, only one of them directly compared the gut microbial composition and connections between MSA and PD using a simple co-abundance network analysis (Barichella et al., 2019). In this study, comprehensive multimodal analyses were conducted to reveal the altered gut microbial compositions in MSA compared with PD and HC.

To minimize the false-positive findings caused by discrepancies between the statistical methods, univariate multimodal analyses were conducted. Consistent reductions in the abundances of Fusicatenibacter and Butyricicoccus were observed in MSA and PD, respectively, in comparison with HC. Similar to other SCFA-producing bacteria, reduced levels of Fusicatenibacter and Butyricicoccus have been associated with PD (Romano et al., 2021). Studies have also shown a positive correlation between the abundance of Butyricicoccus and the levels of SCFAs in PD (Aho et al., 2021). Taken together, the mechanism underlying the relationship between Fusicatenibacter and MSA likely involves the anti-inflammatory effects of SCFAs, which is partially consistent with PD.

Only Limosilactobacillus exhibited a notable difference in abundance between MSA and PD in the univariate analyses. Limosilactobacillus reuteri is a well-studied probiotic bacterium known for its antimicrobial and immunomodulatory effects (Abuqwider et al., 2022). It can stabilize the blood-brain barrier (BBB) dysfunction associated with maternal immune activation in rodents and strengthening the gut epithelial barrier in adults (Kaur et al., 2023; Lu et al., 2023). Furthermore, a randomized controlled trial demonstrated the efficacy of a probiotic mixture containing Limosilactobacillus reuteri in alleviating constipation in PD (Tan et al., 2021). The lower levels of Limosilactobacillus in MSA compared with PD may accelerate the immune response from the gut-blood barrier to the BBB.

In multivariable analyses, Phascolarctobacterium exhibited lower abundances, while Intestinibacter exhibited higher abundances in MSA than in HC. Phascolarctobacterium is capable of producing propionate by fermenting succinate (Watanabe et al., 2012). Given the diverse functions of succinate, including serving as a metabolite for cross-feeding and regulating intestinal homeostasis (Fernández-Veledo and Vendrell, 2019), the associations between Phascolarctobacterium and neurodegenerative diseases have shown inconsistent results. On the other hand, Intestinibacter was found to be more abundant in children with neurodevelopmental disorders (Bojović et al., 2020). However, the relationship between Intestinibacter and PD has not been consistent across studies.

Although α-synuclein deposits predominate in neurons in PD, in MSA they mainly accumulate in oligodendrocytes (Spillantini et al., 1998). Distinct microbial communities in MSA versus PD may generate different signaling mediators, such as bacterial amyloidogenic proteins or metabolites like SCFAs, contributing to cell type-specific α-synuclein aggregation (Schmitt et al., 2023). SCFAs produced by the gut microbiome can modulate neuronal function and microglial maturation (Silva et al., 2020). Of note, among plasma SCFAs derived from the gut microbiome, acetic and propionic acids were decreased in MSA but increased in PD (Chen et al., 2022; He et al., 2021; Yang et al., 2022). Consequently, MSA and PD may be characterized by distinct patterns of gut microbial dysbiosis that lead to divergent plasma SCFA levels, ultimately affecting different neuronal populations and driving disease-specific α-synuclein pathology.

We further conducted an ensemble learning method using random forest classifiers to leverage the comprehensive gut microbial composition to differentiate between MSA, PD, and HC. The AUCs of the random forest classifiers ranged from 0.75 to 0.78, suggesting the potential utility of the gut microbiome in differentiating MSA from PD. Collinsella, which had importance scores > 0.05 in MSA-PD classifiers, is capable of producing SCFAs (Kageyama and Benno, 2015). Collinsella exhibits a greater abundance in numerous neurodegenerative diseases, including Alzheimer’s disease (AD) (Cammann et al., 2023), dementia with Lewy bodies (DLB) (Nishiwaki et al., 2022), and early PD (Huang et al., 2023). Notably, Collinsella was also found to increase with the stage and duration of PD (Zhang et al., 2020). Although the random forest classifiers based on the gut microbiome could differentiate MSA from PD, it was difficult to trace exactly how a prediction was made, limiting their direct clinical application. Furthermore, external validation in independent cohorts is needed to confirm common gut microbial alterations across populations.

It is debatable whether a single bacterial genus can significantly impact the host’s health status. Therefore, we utilized correlation analyses and network analyses to investigate the underlying structure of gut microbial compositions across the study groups. In MSA, a positive correlation between Ruminococcus torques group and Erysipelotrichaceae UCG-003 was identified in SECOM using both linear and distance correlation. Additionally, a positive interaction between Ruminococcus gnavus group and Erysipelatoclostridium was observed in SCNIC and SPIEC-EASI networks of MSA. Ruminococcus torques and Ruminococcus gnavus are mucosa-associated bacteria that can produce SCFAs by degrading mucin (Crost et al., 2023). These bacteria may contribute to the pathogenesis of inflammatory bowel disease by either providing substrates for non-mucolytic bacteria or producing proinflammatory polysaccharides (Wiredu Ocansey et al., 2023). Furthermore, an increased abundance of Ruminococcus torques has been observed in DLB and multiple sclerosis (Nishiwaki et al., 2022; Thirion et al., 2023). On the other hand, Erysipelotrichaceae UCG-003 and Erysipelatoclostridium have been shown to exhibit increased abundances in PD (Lubomski et al., 2022; Rosario et al., 2021). Erysipelatoclostridium was also positively related to the severity of PD (Papić et al., 2022). It is plausible that an interwoven metabolic interaction between these coupled bacteria is responsible for the observed positive correlations in MSA.

Notably, functional inference using PICRUSt2 predicted a significant increase in the abundance of the ARGORNPROST-PWY pathway and a significant decrease in the abundance of the PWY-6478 pathway in MSA compared with HC. The ARGORNPROST-PWY pathway leads to the degradation of L-arginine, the primary substrate for nitric oxide (NO) synthesis (Caspi et al., 2017). NO is an important regulator of the CNS, promoting optimal cerebral blood flow, regulating synaptic plasticity, and modulating neurosecretion (Virarkar et al., 2013). A previous study reported decreased nitrate levels, the degradation product of NO, in the cerebrospinal fluid of MSA patients, indicating reduced CNS production of NO (Kuiper et al., 1994). L-arginine itself can act as an immunosuppressive agent. In a rat model of cerebral ischemia-reperfusion injury, L-arginine inhibits the microglial inflammatory response, thereby exerting a neuroprotective effect (Chen et al., 2020). In PD and AD, agmatine, a metabolite of L-arginine, can reduce oxidative stress and decreases neuronal apoptosis by inhibiting excitatory amino acid-induced neurotoxicity (Xu et al., 2018). L-arginine also inhibits the aggregation of amyloidogenic β-sheet structures by interacting with the hydrophobic regions of proteins (Das et al., 2007). Taken together, the accelerated degradation of L-arginine may contribute to the susceptibility to α-synuclein aggregation in MSA.

On the other hand. the PWY-6478 pathway is responsible for the GDP-D-glycero-α-D-manno-heptose biosynthesis (Caspi et al., 2017). Heptoses are common components of proinflammatory lipopolysaccharides (LPSs) present on the surface of gram-negative bacterial cells (Holst et al., 1991). Gut microbial dysbiosis can lead to excessive leakage of LPS from the gut lining. LPSs primarily bind to toll-like receptor 4 (TLR4) on immune cells, causing excessive production of proinflammatory cytokines. Chronic neuroinflammation and neuronal cell death elicited via the gut-brain axis are key components of neurodegenerative diseases (Kalyan et al., 2022). MSA patients have been shown to exhibit increased TLR4 expression in their colonic mucosa (Engen et al., 2017). Despite TLR4’s role in the secretion of inflammatory mediators, it may also promote intestinal homeostasis and is essential for protection against epithelial injury and bacterial invasion (Craig et al., 2023). In our study, the decreased abundance of the PWY-6478 pathway in MSA may be associated with downregulation of intestinal TLR4 and a reduced reparative response to intestinal injury. More studies are needed to elucidate the roles of LPSs and TLR4 in the gut environment of MSA. Despite a moderate correlation between the genera associated with MSA and PD, investigations into the underlying metabolic mechanisms deserve more attention.

Our study has several limitations. First, the sample size was relatively small due to the rarity of MSA, which limited the statistical power to detect additional genera and metabolic pathways associated with MSA and PD. Furthermore, the limited number of subjects precluded the use of statistical matching techniques; therefore, we conducted multivariable analyses to control for potential confounders. Second, the associations between diseases and changes in the gut microbiome were established through the study’s cross-sectional design. Further longitudinal studies are needed to determine whether the identified microbial shifts are upstream, downstream, or incidental to MSA pathogenesis. Although the comprehensive multimodal analyses lend our study relatively robust internal validity, open datasets for external validation are needed to confirm these gut microbiome-derived signatures. We carefully controlled for potential confounders in the statistical analyses. Nonetheless, information on constipation status and daily intake of vegetables and yogurt was self-reported and assessed using a subjective true/false question, which may have introduced bias. Finally, while statistical tests did not show a significant influence of constipation on the observed microbial differences in MSA and PD, constipation may still exert residual confounding effects, as certain dysbiosis could be exclusively related to constipation in MSA patients.

In summary, our study provides the first evidence of pivotal alterations in gut microbial networks and metabolic pathways specific to MSA. Through comprehensive multimodal analyses, consistent findings were identified, particularly the decreased abundance of Fusicatenibacter in MSA. Although the genera associated with MSA and PD showed a moderate correlation, our correlation and network analyses revealed unique interbacterial interactions specific to MSA, including a novel positive correlation between Ruminococcus gnavus group and Erysipelatoclostridium. Further studies with larger cohorts are needed to establish the clinical utility of fecal microbiome-derived signatures that incorporate interbacterial interactions in MSA. Furthermore, we identified significantly altered metabolic pathways in MSA, such as an increase in the ARGORNPROST-PWY pathway and a decrease in the PWY-6478 pathway. Experimental validation, such as shotgun metagenomic sequencing, should be performed to validate these functional predictions. This provides new directions for future studies on microbial metabolic interactions. Our findings may contribute to the differential diagnosis of MSA and enhance our understanding of its underlying pathogenesis.

References

• 1

  AbuqwiderJ.AltamimiM.MaurielloG. (2022). Limosilactobacillus reuteri in health and disease.Microorganisms10:522. 10.3390/microorganisms10030522

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AhoV. T. E.HouserM. C.PereiraP. A. B.ChangJ.RudiK.PaulinL.et al (2021). Relationships of gut microbiota, short-chain fatty acids, inflammation, and the gut barrier in Parkinson’s disease.Mol. Neurodegener.16:6. 10.1186/s13024-021-00427-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  AhoV. T. E.PereiraP. A. B.VoutilainenS.PaulinL.PekkonenE.AuvinenP.et al (2019). Gut microbiota in Parkinson’s disease: Temporal stability and relations to disease progression.EBioMedicine44691–707. 10.1016/j.ebiom.2019.05.064

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  AndersonM. J. (2001). A new method for non-parametric multivariate analysis of variance.Aus. Ecol.2632–46. 10.1111/j.1442-9993.2001.01070.pp.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BarichellaM.SevergniniM.CiliaR.CassaniE.BolliriC.CaronniS.et al (2019). Unraveling gut microbiota in Parkinson’s disease and atypical parkinsonism.Mov. Disord.34396–405. 10.1002/mds.27581

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BhartiR.GrimmD. G. (2021). Current challenges and best-practice protocols for microbiome analysis.Brief. Bioinform.22178–193. 10.1093/bib/bbz155

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BojovićK.IgnjatovićÐI.Soković BajićS.Vojnović MilutinovićD.TomićM.GolićN.et al (2020). Gut microbiota dysbiosis associated with altered production of short chain fatty acids in children with neurodevelopmental disorders.Front. Cell Infect. Microbiol.10:223. 10.3389/fcimb.2020.00223

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  BokulichN. A.KaehlerB. D.RideoutJ. R.DillonM.BolyenE.KnightR.et al (2018). Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin.Microbiome6:90. 10.1186/s40168-018-0470-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  BolyenE.RideoutJ. R.DillonM. R.BokulichN. A.AbnetC. C.Al-GhalithG. A.et al (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2.Nat. Biotechnol.37852–857. 10.1038/s41587-019-0209-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  BrückD.WenningG. K.StefanovaN.FellnerL. (2016). Glia and alpha-synuclein in neurodegeneration: A complex interaction.Neurobiol. Dis.85262–274. 10.1016/j.nbd.2015.03.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  CallahanB. J.McMurdieP. J.RosenM. J.HanA. W.JohnsonA. J.HolmesS. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data.Nat. Methods13581–583. 10.1038/nmeth.3869

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  CammannD.LuY.CummingsM. J.ZhangM. L.CueJ. M.DoJ.et al (2023). Genetic correlations between Alzheimer’s disease and gut microbiome genera.Sci. Rep.13:5258. 10.1038/s41598-023-31730-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  CaspiR.BillingtonR.FerrerL.FoersterH.FulcherC. A.KeselerI. M.et al (2016). The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases.Nucleic Acids Res.44D471–D480. 10.1093/nar/gkv1164

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  CaspiR.BillingtonR.FulcherC. A.KeselerI. M.KothariA.KrummenackerM.et al (2017). The MetaCyc database of metabolic pathways and enzymes.Nucleic Acids Res.46D633–D639. 10.1093/nar/gkx935

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  ChenS. F.PanM. X.TangJ. C.ChengJ.ZhaoD.ZhangY.et al (2020). Arginine is neuroprotective through suppressing HIF-1α/LDHA-mediated inflammatory response after cerebral ischemia/reperfusion injury.Mol. Brain13:63. 10.1186/s13041-020-00601-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  ChenS. J.ChenC. C.LiaoH. Y.LinY. T.WuY. W.LiouJ. M.et al (2022). Association of fecal and plasma levels of short-chain fatty acids with gut microbiota and clinical severity in patients with Parkinson disease.Neurology98e848–e858. 10.1212/WNL.0000000000013225

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  CraigC. F.FinkelsteinD. I.McQuadeR. M.DiwakarlaS. (2023). Understanding the potential causes of gastrointestinal dysfunctions in multiple system atrophy.Neurobiol. Dis.187:106296. 10.1016/j.nbd.2023.106296

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  CrostE. H.ColettoE.BellA.JugeN. (2023). Ruminococcus gnavus: Friend or foe for human health.FEMS Microbiol. Rev.47:fuad014. 10.1093/femsre/fuad014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  DasU.HariprasadG.EthayathullaA. S.ManralP.DasT. K.PashaS.et al (2007). Inhibition of protein aggregation: Supramolecular assemblies of arginine hold the key.PLoS One2:e1176. 10.1371/journal.pone.0001176

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  DouglasG. M.MaffeiV. J.ZaneveldJ. R.YurgelS. N.BrownJ. R.TaylorC. M.et al (2020). PICRUSt2 for prediction of metagenome functions.Nat. Biotechnol.38685–688. 10.1038/s41587-020-0548-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  EngenP. A.DodiyaH. B.NaqibA.ForsythC. B.GreenS. J.VoigtR. M.et al (2017). The potential role of gut-derived inflammation in multiple system atrophy.J. Parkinsons Dis.7331–346. 10.3233/JPD-160991

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  FedarkoM. W.MartinoC.MortonJ. T.GonzálezA.RahmanG.MarotzC. A.et al (2020). Visualizing ’omic feature rankings and log-ratios using Qurro.NAR Genom. Bioinform.2:lqaa023. 10.1093/nargab/lqaa023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  FernandesA. D.MacklaimJ. M.LinnT. G.ReidG.GloorG. B. (2013). ANOVA-like differential expression (ALDEx) analysis for mixed population RNA-Seq.PLoS One8:e67019. 10.1371/journal.pone.0067019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  Fernández-VeledoS.VendrellJ. (2019). Gut microbiota-derived succinate: Friend or foe in human metabolic diseases?Rev. Endocrine Metab. Disord.20439–447. 10.1007/s11154-019-09513-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  FriedmanJ.AlmE. J. (2012). Inferring correlation networks from genomic survey data.PLoS Comput. Biol.8:e1002687. 10.1371/journal.pcbi.1002687

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  HeX.QianY.XuS.ZhangY.MoC.GuoW.et al (2021). Plasma short-chain fatty acids differences in multiple system atrophy from Parkinson’s disease.J. Parkinsons Dis.111167–1176. 10.3233/JPD-212604

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  HolstO.ZähringerU.BradeH.ZamojskiA. (1991). Structural analysis of the heptose/hexose region of the lipopolysaccharide from Escherichia coli K-12 strain W3100.Carbohydr. Res.215323–335. 10.1016/0008-6215(91)84031-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  HuangB.ChauS. W. H.LiuY.ChanJ. W. Y.WangJ.MaS. L.et al (2023). Gut microbiome dysbiosis across early Parkinson’s disease, REM sleep behavior disorder and their first-degree relatives.Nat. Commun.14:2501. 10.1038/s41467-023-38248-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  HughesA. J.DanielS. E.KilfordL.LeesA. J. (1992). Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: A clinico-pathological study of 100 cases.J. Neurol. Neurosurg. Psychiatry55181–184. 10.1136/jnnp.55.3.181

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  KaehlerB. D.BokulichN. A.McDonaldD.KnightR.CaporasoJ. G.HuttleyG. A. (2019). Species abundance information improves sequence taxonomy classification accuracy.Nat. Commun.10:4643. 10.1038/s41467-019-12669-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  KageyamaA.BennoY. (2015). Collinsella. Bergey’s manual of systematics of archaea and bacteria.Hoboken, NJ: Wiley.

  • Google Scholar
  • View reference in article

• 32

  KalyanM.TousifA. H.SonaliS.VichitraC.SunandaT.PraveenrajS. S.et al (2022). Role of endogenous lipopolysaccharides in neurological disorders.Cells11:4038. 10.3390/cells11244038

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  KaurS.SharmaP.MayerM. J.NeuertS.NarbadA.KaurS. (2023). Beneficial effects of GABA-producing potential probiotic Limosilactobacillus fermentum L18 of human origin on intestinal permeability and human gut microbiota.Microb. Cell Fact.22:256. 10.1186/s12934-023-02264-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  KlannE. M.DissanayakeU.GurralaA.FarrerM.ShuklaA. W.Ramirez-ZamoraA.et al (2021). The gut-brain axis and its relation to Parkinson’s disease: A review.Front. Aging Neurosci.13:782082. 10.3389/fnagi.2021.782082

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  KrismerF.WenningG. K. (2017). Multiple system atrophy: Insights into a rare and debilitating movement disorder.Nat. Rev. Neurol.13232–243. 10.1038/nrneurol.2017.26

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  KuiperM. A.VisserJ. J.BergmansP. L.ScheltensP.WoltersE. C. (1994). Decreased cerebrospinal fluid nitrate levels in Parkinson’s disease, Alzheimer’s disease and multiple system atrophy patients.J. Neurol. Sci.12146–49. 10.1016/0022-510x(94)90155-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  KuoM. C.LuY. C.TaiC. H.SoongB. W.HuF. C.ChenM. L.et al (2022). COQ2 and SNCA polymorphisms interact with environmental factors to modulate the risk of multiple system atrophy and subtype disposition.Eur. J. Neurol.292956–2966. 10.1111/ene.15475

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  KurtzZ. D.MüllerC. L.MiraldiE. R.LittmanD. R.BlaserM. J.BonneauR. A. (2015). Sparse and compositionally robust inference of microbial ecological networks.PLoS Comput. Biol.11:e1004226. 10.1371/journal.pcbi.1004226

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  LinH.PeddadaS. D. (2020). Analysis of compositions of microbiomes with bias correction.Nat. Commun.11:3514. 10.1038/s41467-020-17041-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  LinH.EggesbøM.PeddadaS. D. (2022). Linear and nonlinear correlation estimators unveil undescribed taxa interactions in microbiome data.Nat. Commun.13:4946. 10.1038/s41467-022-32243-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  LuJ.FanX.LuL.YuY.MarkiewiczE.LittleJ. C.et al (2023). Limosilactobacillus reuteri normalizes blood-brain barrier dysfunction and neurodevelopment deficits associated with prenatal exposure to lipopolysaccharide.Gut Microbes15:2178800. 10.1080/19490976.2023.2178800

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  LubomskiM.XuX.HolmesA. J.MullerS.YangJ. Y. H.DavisR. L.et al (2022). The gut microbiome in Parkinson’s disease: A longitudinal study of the impacts on disease progression and the use of device-assisted therapies.Front. Aging Neurosci.14:875261. 10.3389/fnagi.2022.875261

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  MallickH.RahnavardA.McIverL. J.MaS.ZhangY.NguyenL. H.et al (2021). Multivariable association discovery in population-scale meta-omics studies.PLoS Comput. Biol.17:e1009442. 10.1371/journal.pcbi.1009442

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  MandalS.Van TreurenW.WhiteR. A.EggesbøM.KnightR.PeddadaS. D. (2015). Analysis of composition of microbiomes: A novel method for studying microbial composition.Microb. Ecol. Health Dis.26:27663. 10.3402/mehd.v26.27663

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  Martínez-MartínP.Rodríguez-BlázquezC.Mario Alvarez, ArakakiT.ArilloV. C.ChanáP.et al (2015). Parkinson’s disease severity levels and MDS-Unified Parkinson’s disease rating scale.Parkinson. Relat. Disord.2150–54. 10.1016/j.parkreldis.2014.10.026

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  MenozziE.SchapiraA. H. V.BorghammerP. (2025). The gut-brain axis in Parkinson disease: Emerging concepts and therapeutic implications.Mov. Disord. Clin. Pract.10.1002/mdc3.70029[Epub ahead of print].

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  MortonJ. T.MarotzC.WashburneA.SilvermanJ.ZaramelaL. S.EdlundA.et al (2019). Establishing microbial composition measurement standards with reference frames.Nat. Commun.10:2719. 10.1038/s41467-019-10656-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  NishiwakiH.UeyamaJ.KashiharaK.ItoM.HamaguchiT.MaedaT.et al (2022). Gut microbiota in dementia with Lewy bodies.NPJ Parkinsons. Dis.8:169. 10.1038/s41531-022-00428-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  PapićE.RačkiV.HeroM.TomićZ.Starčević-ČižmarevićN.KovandaA.et al (2022). The effects of microbiota abundance on symptom severity in Parkinson’s disease: A systematic review.Front. Aging Neurosci.14:1020172. 10.3389/fnagi.2022.1020172

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  ParksD. H.TysonG. W.HugenholtzP.BeikoR. G. (2014). STAMP: Statistical analysis of taxonomic and functional profiles.Bioinformatics303123–3124. 10.1093/bioinformatics/btu494

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  PedregosaF.VaroquauxG.GramfortA.MichelV.ThirionB.GriselO.et al (2011). Scikit-learn: Machine learning in Python.J. Mach. Learn. Res.122825–2830.

  • Google Scholar
  • View reference in article

• 52

  QuastC.PruesseE.YilmazP.GerkenJ.SchweerT.YarzaP.et al (2013). The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools.Nucleic Acids Res.41D590–D596. 10.1093/nar/gks1219

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  R Core Team (2023). R: A language and environment for statistical computing. R foundation for statistical computing.Vienna: R Core Team.

  • Google Scholar

• 54

  RognesT.FlouriT.NicholsB.QuinceC.MahéF. (2016). VSEARCH: A versatile open source tool for metagenomics.PeerJ4:e2584. 10.7717/peerj.2584

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  RomanoS.SavvaG. M.BedarfJ. R.CharlesI. G.HildebrandF.NarbadA. (2021). Meta-analysis of the Parkinson’s disease gut microbiome suggests alterations linked to intestinal inflammation.NPJ Parkinsons Dis.7:27. 10.1038/s41531-021-00156-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  RosarioD.BidkhoriG.LeeS.BedarfJ.HildebrandF.Le ChatelierE.et al (2021). Systematic analysis of gut microbiome reveals the role of bacterial folate and homocysteine metabolism in Parkinson’s disease.Cell Rep.34:108807. 10.1016/j.celrep.2021.108807

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  SchmittV.MasanetzR. K.WeidenfellerM.EbbinghausL. S.SüßP.RosshartS. P.et al (2023). Gut-to-brain spreading of pathology in synucleinopathies: A focus on molecular signalling mediators.Behav. Brain Res.452:114574. 10.1016/j.bbr.2023.114574

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  SchragA.Ben-ShlomoY.QuinnN. P. (1999). Prevalence of progressive supranuclear palsy and multiple system atrophy: A cross-sectional study.Lancet3541771–1775. 10.1016/s0140-6736(99)04137-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  ShafferM.ThurimellaK.SterrettJ. D.LozuponeC. A. (2023). SCNIC: Sparse correlation network investigation for compositional data.Mol. Ecol. Resour.23312–325. 10.1111/1755-0998.13704

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  SilvaY. P.BernardiA.FrozzaR. L. (2020). The role of short-chain fatty acids from gut microbiota in gut-brain communication.Front. Endocrinol.11:25. 10.3389/fendo.2020.00025

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  SpillantiniM. G.CrowtherR. A.JakesR.CairnsN. J.LantosP. L.GoedertM. (1998). Filamentous alpha-synuclein inclusions link multiple system atrophy with Parkinson’s disease and dementia with Lewy bodies.Neurosci. Lett.251205–208. 10.1016/s0304-3940(98)00504-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  TanA. H.LimS. Y.LangA. E. (2022). The microbiome-gut-brain axis in Parkinson disease - from basic research to the clinic.Nat. Rev. Neurol.18476–495. 10.1038/s41582-022-00681-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  TanA. H.LimS. Y.ChongK. K.A ManapM. A. A.HorJ. W.LimJ. L.et al (2021). Probiotics for constipation in parkinson disease: A randomized placebo-controlled study.Neurology96e772–e782. 10.1212/WNL.0000000000010998

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  ThirionF.SellebjergF.FanY.LyuL.HansenT. H.PonsN.et al (2023). The gut microbiota in multiple sclerosis varies with disease activity.Genome Med.15:1. 10.1186/s13073-022-01148-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  VanacoreN. (2005). Epidemiological evidence on multiple system atrophy.J. Neural Transm.1121605–1612. 10.1007/s00702-005-0380-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  VirarkarM.AlappatL.BradfordP. G.AwadA. B. (2013). L-arginine and nitric oxide in CNS function and neurodegenerative diseases.Crit. Rev. Food Sci. Nutr.531157–1167. 10.1080/10408398.2011.573885

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67

  Vujkovic-CvijinI.SklarJ.JiangL.NatarajanL.KnightR.BelkaidY. (2020). Host variables confound gut microbiota studies of human disease.Nature587448–454. 10.1038/s41586-020-2881-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  WanL.ZhouX.WangC.ChenZ.PengH.HouX.et al (2019). Alterations of the gut microbiota in multiple system atrophy patients.Front. Neurosci.13:1102. 10.3389/fnins.2019.01102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  WangY.KasperL. H. (2014). The role of microbiome in central nervous system disorders.Brain Behav. Immun.381–12. 10.1016/j.bbi.2013.12.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  WatanabeH.SaitoY.TeraoS.AndoT.KachiT.MukaiE.et al (2002). Progression and prognosis in multiple system atrophy: An analysis of 230 Japanese patients.Brain125(Pt 5), 1070–1083. 10.1093/brain/awf117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  WatanabeY.NagaiF.MorotomiM. (2012). Characterization of Phascolarctobacterium succinatutens sp. nov., an asaccharolytic, succinate-utilizing bacterium isolated from human feces.Appl. Environ. Microbiol.78511–518. 10.1128/AEM.06035-11

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  WenningG. K.GeserF.KrismerF.SeppiK.DuerrS.BoeschS.et al (2013). The natural history of multiple system atrophy: A prospective European cohort study.Lancet Neurol.12264–274. 10.1016/S1474-4422(12)70327-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  WenningG. K.StankovicI.VignatelliL.FanciulliA.Calandra-BuonauraG.SeppiK.et al (2022). The movement disorder society criteria for the diagnosis of multiple system atrophy.Mov. Disord.371131–1148. 10.1002/mds.29005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  Wiredu OcanseyD. K.HangS.YuanX.QianH.ZhouM.Valerie OlovoC.et al (2023). The diagnostic and prognostic potential of gut bacteria in inflammatory bowel disease.Gut Microbes15:2176118. 10.1080/19490976.2023.2176118

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  XuW.GaoL.LiT.ShaoA.ZhangJ. (2018). Neuroprotective role of agmatine in neurological diseases.Curr. Neuropharmacol.161296–1305. 10.2174/1570159X15666170808120633

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  YangX.AiP.HeX.MoC.ZhangY.XuS.et al (2022). Parkinson’s disease is associated with impaired gut-blood barrier for short-chain fatty acids.Mov. Disord.371634–1643. 10.1002/mds.29063

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77

  ZhangF.YueL.FangX.WangG.LiC.SunX.et al (2020). Altered gut microbiota in Parkinson’s disease patients/healthy spouses and its association with clinical features.Parkinson. Relat. Disord.8184–88. 10.1016/j.parkreldis.2020.10.034

  • Pubmed Abstract
  • CrossRef
  • Google Scholar