Meta-analysis dataset composition and quality control
To explore the role of the gut microbiome in the MGBA across multiple neurological disorders, we conducted a comprehensive meta-analysis of publicly available 16S rRNA gene amplicon sequencing datasets. We searched relevant studies in public data repositories, such as the NCBI Sequence Read Archive (SRA), European Nucleotide Archive (ENA), and MG-RAST. Our search criteria focused on case-control studies that compared the gut microbiome of patients with AD, PD, or ASD to healthy controls. After conducting a thorough screening and quality control process, we included 31 datasets in our meta-analysis: 6 AD studies, 11 PD studies, and 14 ASD studies (Details see Methods Study Selection). These studies encompass diverse populations from various regions, including Europe, North America, and Asia. The sample sizes ranged from 20 to 507 individuals, with a total of 4,049 samples across all datasets. Detailed information for each study—such as sample size, country of origin, 16S rRNA gene region targeted (e.g., V3-V4, V4), specific disease (AD, PD, or ASD), and project accession numbers are provided in Supplementary Table 1.
Figure 1 illustrates the number of samples included in the meta-analysis for each neurological disorder. Among the three disorders, ASD studies contributed the largest number of samples, with 14 datasets. PD studies followed closely, with 11 datasets, while AD studies comprised 6 datasets. It is important to note that the number of datasets included for each disorder may not fully reflect the extent of research in the field, as our stringent inclusion criteria may have excluded some relevant studies. Nevertheless, the selected studies represent a comprehensive collection of high-quality, case-control investigations into the gut microbiome alterations associated with these neurological disorders, providing a solid foundation for a robust meta-analysis.
To ensure consistency and comparability across datasets, we implemented a standardized bioinformatic pipeline to process the raw 16S rRNA gene sequencing data. First, we downloaded the raw fastq files from the respective repositories and performed quality control using the FastQC tool. We then employed the DADA2 pipeline[12] for denoising, filtering, and merging the paired-end reads, as well as for removing chimeric sequences. DADA2 pipeline was selected for its high accuracy in resolving amplicon sequence variants (ASVs) and its ability to generate reproducible results across diverse studies.
After processing the raw data through the DADA2 pipeline, we identified 310,758 unique ASV across all datasets. The number of ASVs per study ranged from 500 to 70,465, with a median of 5,233 ASVs. Then we aggregated the ASVs at the different taxonomic levels using the Silva v138 reference database[13]. These ASVs represented a wide range of bacterial taxa, spanning 58 phyla, 148 classes, 345 orders, 533 families, and 1,603 genera. The most prevalent phyla across all datasets were Firmicutes, Bacteroidetes, Actinobacteria, and Proteobacteria, which is consistent with previous studies characterizing the human gut microbiome[14].
Alpha and beta diversity profiles of the gut microbiome in the MGBA context
We first evaluated the differences in microbial profiles between three neurological disorders and healthy controls using alpha diversity metrics. Several commonly used indices were calculated, including observed species and Chao1 for richness[15][16], Shannon and Simpson indices for diversity[17][18], and Pielou for evenness[19]. Additionally, we examined dominance patterns and assessed the abundance of rare taxa using measures such as core abundance and rare abundance indices[20][21].
The alpha diversity results, shown in Figure 2, were analyzed through a random-effects meta-analysis pooling estimates across studies. No statistically significant differences were found in alpha diversity measures between patients and controls for AD (p = 0.47) and ASD (p = 0.59), suggesting that the overall richness and evenness of the gut microbiome may not be consistently altered in these disorders [22]. In contrast, PD studies revealed significant differences in alpha diversity compared to controls. Specifically, PD patients showed significantly higher richness, as reflected by a greater number of observed species (mean difference = 0.12, 95% CI: 0.02~0.22, p = 0.02), and higher Chao1 index (mean difference = 0.11, 95% CI: 0.01~0.21, p = 0.02). This indicates that the gut microbiome of PD patients harbors more taxa compared to healthy individuals. Moreover, the rare abundance index was significantly elevated in PD (mean difference = 0.21, 95% CI: 0.07~0.35, p = 0.003), while the core abundance index was lower (mean difference = -0.21, 95% CI: -0.39~-0.03, p = 0.02), indicating an increased prevalence of low-abundance taxa in the PD gut microbiome. These findings suggest that the gut microbiome in PD is characterized by a more diverse and even community, with a higher representation of rare species[23]. Previous studies have also reported this increased alpha diversity in PD, which may be linked to the altered gut environment and impaired gut motility associated with the disease[24].
Next, we investigated beta diversity to assess differences in microbial community composition between samples. Bray-Curtis dissimilarities were calculated and visualized using principal coordinates analysis (PCoA). Across the three disorders, both disease status and study contributed significantly to the variance in beta diversity (PERMANOVA, p < 0.001). However, the “study” factor (e.g., data collected from different labs) had a more substantial impact on the overall gut microbiome composition than disease status. This may be due to study-specific variations in the dataset, such as geographical location, sequencing platform, or DNA extraction method. The effect sizes for “study” (AD: R2 = 0.17; PD: R2 = 0.28; ASD: R2 = 0.28) were consistently larger than those for disease status (AD: R2 = 0.011; PD: R2 = 0.008; ASD: R2 = 0.003) (Figure 3). This finding highlights the importance of controlling for technical and geographical variation in microbiome studies and is consistent with previous reports on the strong influence of study-specific factors on microbiome composition[25][26].
To further explore the impact of disease status on beta diversity while controlling for study-specific effects, we performed a constrained analysis of principal coordinates (CAP) using “study” as a conditioning variable to mitigate variation introduced by different research centers. As is shown in Figure 4, we observed significant differences in microbial community composition between patients and healthy controls across all three neurological disorders (permutation test, p < 0.01). Additionally, a clear separation between patient and control groups along the first principal coordinate axis, which explained 16.56% of the variation for AD, 17.79% for PD, and 23.58% for ASD. These results suggest that neurological disorders are associated with specific alterations in gut microbiome composition. For example, in Parkinson's disease, Scheperjans et al. reported a significant reduction in the relative abundance of Prevotellaceae in PD patients compared to healthy controls, which was correlated with the severity of motor symptoms[27]. However, the modest proportions of variance explained by disease status indicate that additional factors: such as diet, medication, or host genetics--may also play a significant role in shaping the gut microbiome in these disorders. This multifactorial influence on gut microbiome composition in neurological disorders has been emphasized in recent reviews[28][29].
Consensus and specific microbial signatures for neurological diseases
We next investigated the microbial signatures associated with each neurological disorder, aiming to identify both consensus and disease-specific microbial signatures. We leveraged a two-step approach combining different statistical methods. First, we analyzed each of the 31 datasets independently using the ALDEx2 method [30], which has been recommended by previous benchmarking study[31] for its conservative and reliable performance in microbiome differential abundance testing[31]. Then, to identify consensus and disease-specific microbial signatures, we focused on taxa that were consistently differentially abundant across multiple independent studies of the same disorder. This lead to identify robust microbial alterations that were replicable across different cohorts, reducing the impact of study-specific confounding factors.
Figure 5 revealed both shared and disorder-specific microbial signatures across PD, AD, and ASD. The genera Blautia and Bacteroides exhibited differential abundance in all three disorders, suggesting potential common microbial signatures relevant to multiple neurological conditions. However, it is crucial to note that the direction of change varied across disorders. For instance, Blautia was depleted in PD but enriched in AD and ASD, reflecting the complex and disease-specific nature of gut microbiome alterations in neurological disorders[32][33].For instance, Blautia, known for producing short-chain fatty acids (SCFAs), particularly butyrate[34], was depleted in PD but enriched in AD and ASD. Butyrate's neuroprotective and anti-inflammatory properties[35] suggest that Blautia's varying impact on the microbiota-gut-brain axis may be context-dependent and disease-specific. Similarly, Bacteroides, crucial for carbohydrate metabolism and gut homeostasis[36], showed alterations linked to changes in gut permeability and systemic inflammation[37], known to affect the gut-brain axis. These findings align with previous studies, which reported alterations in Bacteroides abundance in AD and PD, respectively[32][38], underscoring the importance of these microbial changes in neurological disorders.
Several genera showed significant differences in two of the three disorders, suggesting potential shared pathways in neurological conditions. UBA1819, Lachnoclostridium, and Bifidobacterium showed significant changes in both PD and AD, potentially indicating common mechanisms in neurodegenerative processes. Bifidobacterium, extensively studied for its probiotic properties and potential neuroprotective effects[39], may play a role in these processes through the production of neuroactive compounds and immune system modulation[40]. Roseburia, Lachnospiraceae UCG-004, Lachnospiraceae ND3007 group, Lachnospira, and Agathobacter were altered in both PD and ASD, hinting at common mechanisms despite the different nature of these disorders. Notably, Roseburia, a butyrate-producing genus crucial for gut barrier integrity and with anti-inflammatory properties[41], suggests a potential shared pathway involving gut barrier function and neuroinflammation. Additionally, Clostridium sensu stricto 1 and [Ruminococcus] torques group were differentially abundant in both AD and ASD, indicating possible shared aspects in cognitive and neurodevelopmental disorders.
Among the disorders studies, PD had the most distinct microbial profile, with 19 genera showing PD-specific alterations. Notably, Akkermansia was significantly enriched in PD patients, a finding consistent with previous studies[38]. This enrichment has been associated with increased gut permeability and inflammation, which are known features of PD pathophysiology[42]. Conversely, Faecalibacterium was depleted in PD patients, potentially contributing to the pro-inflammatory status observed in PD due to its role in producing the anti-inflammatory short-chain fatty acid butyrate[43]. Additional PD-specific alterations were seen in genera such as Fusicatenibacter, Ligilactobacillus, and Porphyromonas, indicating widespread reshaping of the gut microbial community in PD.
In contrast, AD showed relatively few unique microbial alterations, with only three genera specifically altered in AD patients. Among these, Eggerthella was enriched in AD patients, contributing to the growing evidence of specific microbial alterations linked to AD pathogenesis. Eggerthella has been associated with the production of toxic metabolites and low-grade inflammation[44], which could potentially contribute to AD pathogenesis through the gut-brain axis. This finding is consistent with a study by Vogt et al., which also reported an increase in Eggerthella in AD patients[24]. Eggerthella is known for its ability to produce amyloids, which may contribute to the amyloid accumulation characteristic of AD[45]. The other two genera, Alistipes and [Ruminococcus] gnavus group, also exhibited AD-specific alterations. Alistipes has been linked to gut inflammation and mental health disorders in previous studies[46], while the [Ruminococcus] gnavus group has been associated with inflammatory conditions and alterations in gut permeability[47]. These taxa were also identified in a study by Zhuang et al., which found them to be differentially abundant in AD patients compared to healthy controls[48].
ASD exhibited an intermediate level of microbial specificity, with 10 genera uniquely altered. Notably, Escherichia-Shigella and Collinsella were found to be enriched in ASD patients, findings that are consistent with previous studies[49]. The enrichment of these bacteria, along with alterations in other genera like CAG-352, supports the growing evidence for gut microbiome involvement in ASD and suggests that microbial alterations may contribute to the heterogeneous nature of the disorder.
Consensus and specific microbial network structure for neurological diseases
To explore the complex microbial interactions and co-occurrence patterns in the context of neurological disorders, we performed a network analysis using the SpiecEasi (Sparse Inverse Covariance Estimation for Ecological Association and Statistical Inference) method[50]. We first fit networks for both the disease and control groups within each study (as illustrated in Figure 6) and then compared the consensus and unique edges for each disease to identify disease-specific alterations of microbial association.
It was shown that two edges were consistently altered across all three neurological disorders (PD, AD, and ASD): the Coprococcus -- Lachnospiraceae FCS020 group and Gemella – Streptococcus connections. The presence of these shared alterations suggests potential common mechanisms within the gut microbiome that may play a role in various neurological conditions. The first connection involves known butyrate producers, which are crucial for gut health and neuroprotection [51], while the second connection involves genera that, when altered, have been associated with inflammatory conditions [52]. Disruptions in these microbial interactions could impair butyrate production and contribute to systemic inflammation, both of which are implicated in the pathogenesis of multiple neurological disorders [53].
We summarized these edge alternations in Figure 7. PD exhibited the most extensive changes, with 199 disease-specific edges, indicating widespread disruption of microbial interactions in the PD gut microbiome. AD demonstrated 63 disease-specific edges, while ASD had 53 unique edges. These findings highlight that all three disorders involve distinct microbial network alterations compared to the controls, with PD showing the most substantial changes, followed by AD and ASD.
Interestingly, PD and AD shared the greatest number of common alterations, with 12 shared edges. This overlap suggests that there may be common microbial mechanisms or pathways between these two neurodegenerative disorders. This could reflect certain shared pathophysiological features or common gut-brain axis that play a role in both PD and AD [54].
The identification of these disease-specific and shared microbial network alterations provides insights into the intricate interplay between gut microbes in neurological disorders. The larger number of PD-specific alterations, in particular, highlights the profound impact of this disease on gut microbial interactions. Moreover, the shared microbial network changes across all three disorders, as well as the significant overlap between PD and AD, suggest the existence of common gut-brain axis pathways that may contribute to multiple neurological conditions.