Abstract
Background:
Lewy body disease (LBD) is a progressive neurodegenerative a-synucleinopathy, whereas isolated REM sleep behavior disorder (iRBD) is recognized as a prodromal stage of LBD. Although growing evidence implicates the gut–brain axis in neurodegeneration, the taxonomic and functional roles of the gut microbiome across the prodromal-to-symptomatic LBD continuum remain poorly defined.
Methods:
Here, we performed shotgun metagenomic sequencing on stool samples from 25 patients with LBD (10 mild cognitive impairment due to LBD [MCI-LB] and 15 dementia with Lewy bodies [DLB]), 10 individuals with iRBD, and their household matched cohabitant controls to characterize disease-associated microbial alterations while minimizing environmental confounding.
Results:
Despite no significant differences in global microbial diversity, we identified convergent shifts in microbial taxa, metabolic pathways, and gene families across disease stages. Both LBD and iRBD showed increased abundance of microbial taxa potentially associated with gut barrier disruption, as well as higher abundance of functional pathways related to lipopolysaccharide biosynthesis. LBD showed lower abundance of pathways related to complex carbohydrate fermentation, and both groups showed lower abundance of pathways associated with neurotransmitter-related metabolism. In particular, pathways and gene families associated with starch degradation were reduced in LBD, and those associated with histidine-to-glutamate/ GABA metabolism were reduced in both groups.
Discussion:
These exploratory findings represent the first high-resolution, shotgun metagenomic characterization of gut microbiome alterations across the LBD continuum, highlighting functional patterns that may serve as candidate markers of disease progression in future longitudinal and mechanistic studies.
Introduction
Lewy body disease (LBD) refers to a spectrum of clinical conditions associated with Lewy body pathology, including mild cognitive impairment due to Lewy body disease (MCI-LB) and dementia with Lewy bodies (DLB) (). It is a progressive neurodegenerative disease caused by the accumulation of abnormally phosphorylated α-synuclein (). Clinically, LBD presents with a broad range of symptoms, including progressive cognitive decline, fluctuations in attention, recurrent visual hallucinations, autonomic dysfunction, parkinsonism, and REM sleep behavior disorder (RBD) (; Walker et al., 2015; ). In current diagnostic criteria, the timing of cognitive and motor symptom onset helps distinguish LBD subtypes from other synucleinopathies, such as Parkinson’s disease (PD). Among these distinguishing clinical features, RBD is particularly prominent and occurs in the vast majority of individuals with LBD (; ; ; ; ). When RBD occurs in the absence of a diagnosed neurodegenerative disease, it is termed isolated RBD (iRBD), which is widely considered a prodromal marker of α-synucleinopathies (; ; ). Supporting this concept, one study reported that few cases clinically diagnosed with iRBD have shown evidence for LBD pathology at post-mortem examination (). Taken together, these findings highlight the close link between iRBD and LBD and support their classification along a common disease continuum ().
The gut-brain axis is a bidirectional communication system between the central nervous system and the gastrointestinal tract (; ). This system operates through multiple interconnected pathways, including immune signaling, vagal and enteric neural circuits, endocrine mechanisms, and microbiota-derived metabolites and neurotransmitters (; ; ). Growing evidence suggests that the gut-brain axis plays a critical role in neurodegenerative diseases, with the gut microbiome increasingly implicated across several conditions (; ). In PD, for example, multiple studies have consistently reported a decreased abundance of short-chain fatty acids (SCFAs) producing bacteria and an increased prevalence of pro-inflammatory microbial taxa (; ; ; ).
Despite this growing body of work in PD, relatively few studies have examined the gut microbiome in the cognitive manifestations of LBD (i.e., DLB and MCI-LB) or in its prodromal stage, iRBD. Nishiwaki et al. analyzed stool samples from 28 patients with DLB, 26 patients with iRBD, and 147 controls (). They reported decreased abundance in SCFAs producing bacterial genera, and increased levels of Ruminococcus torques and Collinsella in DLB compared with controls. More recently, Teigen et al. evaluated 27 LBD patients, 11 iRBD patients, and their cohabitant controls, and found no major differences in overall microbiome diversity, although the iRBD group showed a notably lower abundance of Bacteroides (Teigen et al., 2024).
To date, microbiome studies in LBD have relied on 16S rRNA gene amplicon sequencing, which targets a single marker gene and provides limited taxonomic resolution (primarily at the genus level) and minimal functional insight. As a result, this approach cannot fully characterize the metabolic capabilities of the microbiome or its potential contributions to disease pathology. In contrast, whole genome shotgun (WGS) metagenomic sequencing enables functional potential profiling of the microbiome by identifying microbial metabolic pathways and gene families encoded within the community (). Such functional information is particularly important for LBD, as gut microbes produce SCFAs, bile acids, neurotransmitters, and other metabolites that modulate host immune responses, intestinal barrier integrity, and neural signaling (; Strandwitz, 2018; ). Altered microbial metabolism may therefore influence processes relevant to neurodegeneration, including neuroinflammation and oxidative stress (; ). Therefore, studying these functional alterations provides an opportunity to identify mechanistic links between the gut microbiome and LBD pathogenesis.
In microbiome research, a high degree of inter-individual variation is commonly observed (The-Human-Microbiome-Project-Consortium, 2012; ), driven largely by environmental factors including diet, lifestyle, and medication use (Turnbaugh et al., 2009; Zhernakova et al., 2016; ). Importantly, individuals who share a household tend to have more similar gut microbiota than unrelated individuals (Song et al., 2013; ; Tavalire et al., 2021). For this reason, the inclusion of household controls is critical for distinguishing true disease-associated microbial signatures from differences attributable to environmental or dietary variation. Without such controls, observed microbial differences may reflect lifestyle or environmental disparities rather than disease processes. Notably, only one prior LBD study—using 16S rRNA sequencing—has incorporated household-matched controls (Teigen et al., 2024), limiting the interpretability and comparability of findings from studies that did not account for this essential source of confounding.
In this study, we performed WGS metagenomic sequencing on stool samples from patients with LBD and iRBD, together with their household cohabitants. We assessed differences in microbiome diversity, taxonomic composition, and metabolic functional potential across these groups, and examined associations between microbiome features and clinical measures of cognition and motor function. By integrating high-resolution functional profiling with a household-matched design, our study aimed to identify disease-associated taxonomic and metabolic shifts across the LBD continuum while minimizing confounding from shared environmental factors.
Materials and methods
Participant enrollment
Participant enrollment followed procedures described in our previous work (Teigen et al., 2024). This study was approved by the Mayo Clinic Institutional Review Board (IRB: 18-011778), and informed consent was obtained from all participants. All procedures involving human participants were conducted in accordance with the ethical standards of the institutional research committee and with the Declaration of Helsinki and its later amendments. Participants were recruited through the Mayo Clinic Functional Genomics of LBD program, Mayo Clinic Alzheimer’s Disease Research Center and the North American Prodromal Synucleinopathy (NAPS) cohort. The study population included 25 patients with LBD who met research criteria for MCI-LB (n = 10) or DLB (n = 15), 10 patients with polysomnography-confirmed iRBD, as well as their household-matched (cohabitant) controls. All patients underwent a comprehensive clinical interview, neurological examination, and neuropsychological testing by sleep and behavioral neurology subspecialists and neuropsychologists. Following evaluation, a consensus clinical diagnosis of LBD or iRBD was assigned based on established clinical criteria. Since consensus diagnostic criteria for iRBD are not established, iRBD in this study was confirmed by polysomnography in the setting of normal neuropsychological testing and absence of parkinsonism. Demographic and clinical information was collected for all participants if applicable, including age, sex, dementia severity (Clinical Dementia Rating [CDR] global and sum of boxes [SB]), bedside cognitive testing score (Kokmen Short Test of Mental Status [STMS], Montreal Cognitive Assessment [MoCA]), and parkinsonism severity (Movement Disorders Society-Unified Parkinson’s Disease Rating Scale Part III [MDS-UPDRS III]).
Stool sample collection, DNA extraction, and shotgun metagenomic sequencing
Fecal samples were collected using standardized stool specimen collection kits at the time of clinical assessment. Samples were promptly transferred to the University of Minnesota Genomics Center for DNA extraction and sequencing. Fecal DNA was extracted with the DNeasy 96 PowerSoil Pro QIAcube HT Kit (QIAGEN, Germantown, MD, USA) using the inhibitor removal technology (IRT) protocol according to the manufacturer’s instructions. DNA concentration and purity were measured with a NanoDrop-8000 UV-Vis spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA) and a PicoGreen assay. Sequencing libraries were prepared and run on an Illumina NovaSeq 6000 using 2 × 150 bp paired-end chemistry, targeting approximately 8 million read pairs per sample.
Quality filtration of sequenced reads
Metagenomic reads were processed using an in-house quality control pipeline. Specifically, Trimmomatic (version 0.39) () was used to remove adapter and overrepresented sequences (ILLUMINACLIP: adaptor_overrepresented.fa:2:30:10:2:keepBothReads), trim low-quality bases from the start and end of reads (LEADING:3 and TRAILING:3), and discard reads shorter than 60 bp (MINLEN:60). Host-derived reads were subsequently removed by aligning to the Human Genome Reference Database (GRCh38) using Bowtie2 (version 2.5.0) (), followed by filtering with SAMtools (version 1.18) ().
Taxonomic and functional profiling of stool metagenomes
Taxonomic profiling was performed using the MetaPhlAn4 (version 4.1.15) () phylogenetic clade identification pipeline with default parameters. Briefly, MetaPhlAn4 identifies microbial taxonomy and composition by detecting unique clade-specific marker genes. The marker database (mpa_vJun23_CHOCOPhlAnSGB_202403) is derived from a curated collection of 1.01 million prokaryotic reference and metagenome-assembled genomes, encompassing 26,970 species-level genome bins, 4,992 of which are taxonomically unidentified at the species level. Microbial taxa of viral origin and those labeled as unclassified or unknown were excluded from further analyses. Microbiome profiles were then normalized using total sum scaling (TSS) to obtain the relative abundances of microbial taxa. The relative abundance profiles of all detected microbial taxa are provided in Supplementary Table 1.
Functional profiling of microbial communities was conducted using HUMAnN3 (version 3.96) () with default parameters. Briefly, HUMAnN3 first maps quality-controlled metagenomic reads to species-specific pangenomes to quantify gene families and subsequently aligns unmapped reads to the UniRef90 database to improve detection sensitivity. Gene family abundances are then used to reconstruct microbial metabolic pathways based on the MetaCyc database. Pathways classified as unmapped (i.e., reads that did not align to any gene family) or unintegrated (i.e., gene families not assigned to any known pathway) were excluded from downstream analyses. The relative abundances of the remaining pathways were normalized using TSS via the “humann_renorm_table --units relab --special n” utility in HUMAnN3. In parallel, the relative abundances of gene families were also characterized for downstream analyses. The relative abundance profiles of all detected microbial pathways are provided in Supplementary Table 2.
Gut microbiome diversity analysis
Overall gut microbiome composition was evaluated by calculating both α-diversity (species-level Shannon index and richness) and β-diversity (Bray–Curtis distances between all sample pairs). The diversity() function from R package vegan (version 2.7.3) () was used to calculate Shannon index and species richness based on untransformed relative abundances of microbial species in each stool metagenome. To compare α-diversity between LBD or iRBD patients and their cohabitant controls, mixed-effects linear regression models were constructed in R using the lmer() function from the lmerTest package (version 3.2.1) (), with household ID included as a random effect. P-values less than 0.05 were considered statistically significant. For β-diversity analysis, principal coordinate analysis (PCoA) was performed using Bray–Curtis dissimilarities calculated on arcsine square-root-transformed relative abundances of microbial species and pathways identified by MetaPhlAn4 and HUMAnN3. This analysis was implemented using the R packages ade4 (version 1.7.24) () and vegan (version 2.7.3) (). All statistical analyses described above were conducted in R (version 4.5.3).
PERMANOVA on taxonomic composition of microbial communities
A permutational multivariate analysis of variance (PERMANOVA) was conducted on the Bray–Curtis dissimilarity matrix using the adonis2() function from the R package vegan (version 2.7.3) (). P-values for the pseudo-F statistic were calculated based on 999 permutations to assess the proportion of variance in gut microbial community composition explained by case-control status and patient characteristics such as age, sex, and BMI. For comparisons between LBD or iRBD patients and their cohabitant controls, permutations were constrained within households using the strata option to account for paired study design. The output from PERMANOVA was used to estimate the percentage of variance in microbiome composition attributable to case-control status and its statistical significance. P-values less than 0.05 were considered statistically significant. All statistical analyses described above were conducted in R (version 4.5.3).
Differential abundance and prevalence analysis
Microbial taxonomic and functional features (species, pathways, and gene families) were preprocessed using two filtering steps prior to differential analysis. First, to reduce the influence of extremely low-abundance features that may represent noise or technical artifacts, relative abundance values below a predefined threshold were set to zero based on rank-abundance plots (Supplementary Figures 1A–C). Specifically, presence/absence thresholds for each feature type were determined empirically by ranking all features by their relative abundance and identifying the elbow point of the resulting rank-abundance curve, beyond which the distribution drops sharply toward the detection limit. This approach yielded thresholds of 10–4.7,10–4.5, and 10–7 for microbial species, metabolic pathways, and gene families, respectively (Supplementary Figure 1). Second, to reduce extreme data sparsity inherent to microbiome data, microbial features were filtered using a 10% prevalence cutoff. That is, taxonomic or functional features detected in fewer than 10% of samples within each group comparison were excluded prior to downstream differential analyses. The 10% prevalence threshold was chosen based on established practice in microbiome studies (; Yang and Chen, 2022), where features present in fewer than 10% of samples are considered too sparse to yield reliable statistical estimates. Mixed-effects linear regression models were applied to arcsine square-root-transformed relative abundance data to identify differentially abundant features (species, pathways, and gene families) between LBD or iRBD patients and their cohabitant controls. The models included age and BMI as fixed effects and household ID as a random effect. Age and BMI were included as covariates given their well-documented influence on gut microbiome composition (Turnbaugh and Gordon, 2009; ). Although antibiotic use and constipation are also recognized modulators of the gut microbiome (Zhu et al., 2014; ), these variables could not be included in the models due to data availability constraints. Antibiotic use data were missing for six participants across the LBD and control groups, and constipation status was only available for LBD and iRBD participants but not for their cohabitant controls. All models were fitted in R using the lmer() function from the lmerTest package (version 3.2.1) (). To quantify the magnitude of microbial species and pathway relative abundance differences between patients and their cohabitant controls, Cohen’s d was calculated for each microbial species and pathway feature. All analyses were performed in R using the cohens_d() function from the effectsize package (version 1.0.2) (). Differentially prevalent microbial features were identified using Fisher’s exact test. All P-values reported here were unadjusted for multiple comparisons. Unadjusted P-values less than 0.05 were considered statistically significant. All statistical analyses described above were conducted in R (version 4.5.3).
Correlations between clinical characteristics and microbial features
To investigate correlation between microbial feature abundance and clinical measurements (functional impairment [CDR-SB], cognitive performance [MoCA and STMS], and motor severity [MDS-UPDRS III]), Spearman partial correlation analyses were performed separately for the LBD (n = 25) and iRBD (n = 10) groups. To reduce sparsity and improve the robustness of the analysis, untransformed relative abundances of microbial features (species and pathways) were filtered by prevalence; only features present in more than 50% of samples within each group were retained for downstream analysis.
For each retained feature, the Spearman partial correlation coefficient ρ was estimated against each clinical measurement, adjusting for age and BMI as covariates using the pcor.test() function from the ppcor package (version 1.1) () in R. Features with |ρ| > 0.4 and unadjusted P < 0.05 were considered statistically significant and retained for further analysis.
Given the small sample sizes, particularly in the iRBD group, a resampling-based validation was applied to all features meeting the above criteria to evaluate the robustness of each observed correlation. An empirical P-value was derived from a permutation null distribution constructed as follows: in each of 1,000 iterations, a random subsample of up to n−1 observations (80–96% for LBD; 80–90% for iRBD) of the available observations was drawn without replacement, the clinical measurement labels were randomly permuted within the subsample to break any true correlation, and a Spearman partial correlation adjusting for age and BMI was computed. The empirical P-value was defined as the proportion of null |ρ| values that equaled or exceeded the observed |ρ| from the full sample. All P-values reported here are unadjusted for multiple comparisons. Features with an unadjusted empirical P < 0.05 were considered to show robust correlations with the clinical characteristic of interest. All statistical analyses described above were conducted in R (version 4.5.3).
Results
Participant demographics
This study includes stool samples from 25 patients with LBD, 10 patients with iRBD, and their cohabitant controls (Figure 1A). Participant demographic and clinical characteristics are shown in Table 1. LBD patients ranged in age from 46 to 73 years (median = 68), and iRBD patients from 56 to 76 years (median = 68). Both groups were slightly older than their cohabitant controls; however, the differences were not statistically significant (P = 0.099 and 0.149, respectively; Mann–Whitney U test). Male sex accounted for 84% and 90% of LBD and iRBD patients, respectively, indicating a strong male predominance in both groups (P ≤ 0.001 for both; Fisher’s exact test). Although five iRBD patients had a global CDR score of 0.5 based on reported degree of cognitive impairment, their neuropsychometric performance was within normal limits and they were therefore classified as having iRBD rather than MCI-LB.
Figure 1
Table 1
| Patient characteristics | LBDa | iRBDb | ||||
|---|---|---|---|---|---|---|
| Patients (n = 25) | Cohabitant controls (n = 25) | P-value# | Patients (n = 10) | Cohabitant controls (n = 10) | P-value# | |
| Age | ||||||
| Median [Q1, Q3]c | 68.0 [64.0, 73.0] | 64.0 [59.0, 71.0] | 0.099 | 68.0 [63.5, 69.8] | 63.0 [57.5, 65.5] | 0.149 |
| Sex | ||||||
| Female, n (%) | 4 (16) | 23 (92) | < 0.001 | 1 (10) | 9 (90) | 0.001 |
| Male, n (%) | 21 (84) | 2 (8) | 9 (90) | 1 (10) | ||
| BMI | ||||||
| Median [Q1, Q3] | 26.4 [25.1, 31.2] | 26.6 [25.1, 29.3] | 0.288 | 28.9 [26.9, 29.4] | 27.5 [24.2, 29.5] | 0.533 |
| Not available | 2 | 10 | 0 | 3 | ||
| Global CDRd, n (%) | ||||||
| 0 | 0 (0) | N/A | 5 (50) | N/A | ||
| 0.5 | 16 (64) | 5 (50) | ||||
| 1 | 8 (32) | 0 (0) | ||||
| 2 | 0 (0) | 0 (0) | ||||
| 3 | 1 (4) | 0 (0) | ||||
| CDR-SBe, n (%) | ||||||
| 0 | 0 (0) | N/A | 5 (50) | N/A | ||
| 0.5–4.0 | 15 (60) | 5 (50) | ||||
| 4.5–9.0 | 9 (36) | 0 (0) | ||||
| 9.5–15.5 | 0 (0) | 0 (0) | ||||
| 16.0–18.0 | 1 (4) | 0 (0) | ||||
| MDS-UPDRS IIIf | ||||||
| Median [Q1, Q3] | 19.0 [1.0, 26.0] | N/A | 0.0 [0.0, 2.5] | N/A | ||
| MoCAg | ||||||
| Median [Q1, Q3] | 21.0 [16.8, 27.3] | N/A | 25.5 [25.0, 28.5] | N/A | ||
| Not available | 1 | 0 | ||||
| STMSh | ||||||
| Median [Q1, Q3] | 31.0 [27.0, 34.0] | N/A | 36.5 [36.0, 37.8] | N/A | ||
| Visual hallucinations, n (%) | 9 (36) | N/A | 0 (0) | N/A | ||
| Cognitive fluctuations, n (%) | 10 (40) | N/A | 0 (0) | N/A | ||
| Parkinsonism, n (%) | 18 (72) | N/A | 1^ (10) | N/A | ||
| REM sleep behavior disorder, n (%) | 24 (96) | N/A | 10 (100) | N/A | ||
| Antibiotic use in the last 6 months, n (%) | ||||||
| n (%) | 8 (32) | 4 (16) | 1 (4) | 2 (8) | ||
| Not available | 4 | 2 | 1 | 1 | ||
| Constipation in the last 6 months, n (%) | ||||||
| n (%) | 20 (80) | N/A | 5 (50) | N/A | ||
| Not available | 1 | 0 | ||||
Demographic and clinical characteristics of the study participants.
aLBD, Lewy body disease; biRBD, isolated REM sleep behavior disorder; cUpper and lower quartiles; dGlobal CDR, Global Clinical Dementia Rating; eCDR-SB, Clinical Dementia Rating Scale Sum of Boxes; fMDS-UPDRS III, Movement Disorder Society-Unified Parkinson’s Disease Rating Scale Part III score; gMoCA, Montreal Cognitive Assessment test score; hSTMS, Short Test of Mental Status score; #Mann-Whitney U test and Fisher’s exact test was used to test for the statistical significance of age and sex, respectively; and mixed-effects linear regression was used to test for the statistical significance of BMI accounting for household as a random effect; N/A, not applicable or not available. ^One iRBD participant had mild bradykinesia suggestive of very mild parkinsonism but without other cardinal manifestations that would meet criteria for PD.
Overall comparisons of gut microbiome communities between LBD or iRBD patients and their cohabitant controls
Species-level Shannon index and richness did not show significant differences between LBD or iRBD patients and their cohabitant controls (Figures 1B, C). This suggests that overall microbial α-diversity at the species level was similar between disease and control groups. For β-diversity, PERMANOVA analysis showed that disease diagnosis explained 2.2% and 5.3% of the total variance in gut microbial community composition in LBD and iRBD cohorts, respectively; however, the variance explained by case-control status was not significant in either cohort (P = 0.067 and 0.156; Figures 1D, E). α- and β-diversity analysis based on metabolic pathway relative abundance also showed no significant differences between disease groups and their cohabitant controls (Supplementary Figures 2A–D).
We also performed PERMANOVA analysis to evaluate other patient characteristics that might contribute to variance in gut microbial communities of LBD and iRBD cohorts. In seperate models, age, sex, and BMI explained 2.2%, 2.0%, and 2.6%, respectively, of the total variance in gut microbial communities among individuals from the LBD cohort, and 5.0%, 4.7%, and 5.3%, respectively, among those from the iRBD cohort (Table 2). However, none of these variables explained a significant amount of variance in either cohort (Table 2).
Table 2
| Cohort | Patient characteristics | Variance explained (%) | P-value# |
|---|---|---|---|
| LBD (n = 25) | Case-control status | 2.2 | 0.067 |
| Age | 2.2 | 0.076 | |
| Sex | 2.0 | 0.114 | |
| BMI | 2.6 | 0.235 | |
| iRBD (n = 10) | Case-control status | 5.3 | 0.156 |
| Age | 5.0 | 0.567 | |
| Sex | 4.7 | 0.339 | |
| BMI | 5.3 | 0.378 |
Patient characteristics contributing to the variance in gut microbial community composition.
#Permutational Multivariate Analysis of Variance (PERMANOVA) was used to test for the association between the corresponding patient characteristic and compositional variance within the gut microbiome.
The relative abundance of microbial taxonomic families in the gut microbiomes of individuals with LBD or iRBD and their cohabitant controls is shown in Figures 1F, G. Lachnospiraceae was the most abundant family in all four groups (LBD, LBD controls, iRBD, and iRBD controls) (median relative abundance = 0.304, 0.312, 0.273, 0.307, respectively). Oscillospiraceae and Bacteroidaceae were also abundant across all groups. The median relative abundance of Oscillospiraceae was 0.160, 0.168, 0.181, and 0.232 in the LBD, LBD controls, iRBD, and iRBD controls, respectively; while Bacteroidaceae had median values of 0.050, 0.087, 0.000, and 0.085 in the same groups. Some families, such as Akkermansiaceae and Bifidobacteriaceae, were present in high abundance in a few individuals but absent in others. Overall, the microbial profiles at the family level indicate clear inter-individual heterogeneity.
Differentially abundant and prevalent gut microbiome features between LBD patients and their cohabitant controls
Although no significant differences in α- and β-diversity were observed between LBD and their cohabitant controls, differential abundance and prevalence analysis of microbial taxa and metabolic pathways revealed several compositional differences between disease and cohabitant control groups. Filtering with the prevalence cut-off (described in Materials and Methods) yielded a total of 492 microbial species and 433 metabolic pathways for downstream differential analyses. Using a mixed-effects linear regression model adjusted for age and BMI, with household ID as a random effect, we identified that four species were significantly more abundant in LBD, while nine were elevated in controls (unadjusted P < 0.05; median |Cohen’s d| = 0.533; range: 0.368–0.692; Supplementary Figure 3A; Supplementary Table 3). Microbial species that were more abundant in LBD included Oscillospiraceae bacterium CLA-AA-H250, Dorea sp. AF36-15AT, GGB9522 SGB14921 (Firmicutes), and Hungatella hathewayi, whereas GGB9480 SGB14874 (Firmicutes), Roseburia hominis, Eubacterium ramulus, GGB2653 SGB3574 (Firmicutes), and Intestinimonas butyriciproducens were more abundant in controls (Figure 2A). Applying the same model to microbial metabolic pathways, we identified four pathways with significantly higher relative abundance in LBD and six pathways elevated in controls (unadjusted P < 0.05; median |Cohen’s d| = 0.451; range: 0.397–0.593; Supplementary Figure 3B, Supplementary Table 4). Pathways enriched in LBD included L-glutamine biosynthesis III, the superpathway of L-methionine biosynthesis by sulfhydrylation, ADP-L-glycero-β-D-manno-heptose biosynthesis, and heme b biosynthesis I (aerobic), whereas pathways enriched in the control group included starch degradation III, β-(1,4)-mannan degradation, L-glutamate and L-glutamine biosynthesis, L-histidine degradation III, and purine nucleobase degradation II (anaerobic) (Figure 2B). Full results for differentially abundant microbial species and pathways, including FDR-adjusted P-values, are provided in Supplementary Tables 3 and 4.
Figure 2
In addition to differentially abundant microbial features, we found that six microbial species and two metabolic pathways were significantly more prevalent in LBD, while nine microbial species were more prevalent in controls (unadjusted P < 0.05, Fisher’s exact test; Figure 2C; Supplementary Tables 5 and 6). When viewed together with the differential abundance results, several microbial species showed consistent patterns. Oscillospiraceae bacterium CLA-AA-H250 was not only more abundant in LBD but also present in a greater number of LBD individuals. Similarly, species such as Eubacterium rectale, GGB2653 SGB3574 (Firmicutes), GGB9480 SGB14874 (Firmicutes), and Longicatena caecimuris exhibited both higher relative abundance and prevalence in the control group, which indicates these taxa may serve as robust microbial signatures distinguishing LBD from controls. Full results for differentially prevalent microbial species and pathways, including FDR-adjusted P-values, are provided in Supplementary Tables 5 and 6.
Differentially abundant and prevalent gut microbiome features between iRBD patients and their cohabitant controls
After applying the prevalence cut-off (described in Materials and Methods), 589 microbial species and 439 metabolic pathways remained for further analysis. With a mixed-effects linear regression model adjusted for age, BMI, and household ID, we found 23 microbial species that were significantly more abundant in iRBD; and four species showing higher abundance in controls (unadjusted P < 0.05; median |Cohen’s d| = 0.696; range: 0.432–1.278; Supplementary Figure 3C; Supplementary Table 7). Microbial species that were more abundant in iRBD included Collinsella aerofaciens, Blautia luti, Mediterraneibacter butyricigenes, Luoshenia tenuis, and Anaerovorax odorimutans, whereas the control group showed higher abundance of Clostridiaceae bacterium, GGB3653 SGB4964 (Firmicutes), Clostridia bacterium UC5.1-1D1, and Enterocloster bolteae (Figure 3A). Metabolic pathway analysis using the same approach identified 26 pathways that had significantly higher relative abundance in iRBD, as well as 39 pathways being more abundant in controls (unadjusted P < 0.05; median |Cohen’s d| = 0.755; range: 0.008–1.984; Supplementary Figure 3D; Supplementary Table 8). The top five pathways enriched in iRBD included O-antigen building blocks biosynthesis (E. coli), Lactose and galactose degradation I, Pyrimidine deoxyribonucleotide phosphorylation, (S)-propane-1,2-diol degradation, and Bile acids epimerization, while the control group showed greater abundance of pathways such as 4-amino-2-methyl-5-diphosphomethylpyrimidine biosynthesis II, Superpathway of pyridoxal 5’-phosphate biosynthesis and salvage, Pyridoxal 5’-phosphate biosynthesis I, L-histidine degradation III, and Purine nucleobases degradation II (anaerobic) (Figure 3B). Full results for differential abundance analysis of microbial species and pathways, including FDR-adjusted P-values, are provided in Supplementary Tables 7 and 8.
Figure 3
In addition to differentially abundant microbial features, we identified three microbial species and one metabolic pathway that were significantly more prevalent in iRBD patients compared with their cohabitant controls (unadjusted P < 0.05, Fisher’s exact test; Figure 3C; Supplementary Tables 9 and 10). When integrating results from abundance and prevalence analyses, several features demonstrated consistent group-wise differences. Specifically, Luoshenia tenuis, Clostridiaceae unclassified SGB15090, and the Bile acids epimerization pathway were not only significantly more abundant in the iRBD group but also present in a greater proportion of iRBD individuals, which suggests these features as potential microbial signatures associated with iRBD. Full results for differential prevalence analysis of microbial species and pathways, including FDR-adjusted P-values, are provided in Supplementary Tables 9 and 10.
Microbial features characterizing the disease continuum of Lewy body disease
Several microbial features showed consistent directional changes in both LBD and iRBD compared with their cohabitant controls, although statistical significance was reached only in the LBD cohort (Supplementary Figure 4). For example, Oscillospiraceae bacterium CLA-AA-H250, Akkermansia muciniphila, Actinomyces oris, GGB3730 SGB5060 (Firmicutes), and GGB9627 SGB15081 (Firmicutes) were more prevalent in both disease groups, but these increases were statistically significant only in LBD. Conversely, Clostridia bacterium UC5.1-1D1, Longicatena caecimuris, and GGB9719 SGB15272 (Firmicutes) were more commonly detected in cohabitant controls, again with significance limited to LBD. Similar patterns were observed for two metabolic pathways, Superpathway of L-citrulline metabolism and L-Nδ-acetylornithine biosynthesis, which were more prevalent in both disease groups but reached statistical significance only in LBD. Together, these findings suggest that several microbial alterations may be present during the prodromal phase (iRBD) but intensify or become more consistently detectable as the disease progresses to overt LBD. This progression-like pattern supports the concept of a microbial signature that evolves along the Lewy body disease continuum.
Shared microbial features between LBD and iRBD compared with cohabitant controls
To identify gut microbial alterations common across the LBD disease spectrum, we examined features that were altered in the same direction in both LBD and iRBD relative to their cohabitant controls. At the species level, Clostridiaceae unclassified SGB15090 was significantly more prevalent in both LBD and iRBD groups (unadjusted P < 0.05; Figures 2C, 3C). At the genus level, Bacteroides was significantly less abundant in both LBD and iRBD, consistent with patterns previously reported in α-synucleinopathies (; ; Teigen et al., 2024) (unadjusted P < 0.05; Supplementary Figure 5; Supplementary Tables 11 and 12). At the pathway level, the relative abundances of L-histidine degradation III, Purine nucleobases degradation II (anaerobic), and L-glutamate and L-glutamine biosynthesis were consistently lower in both disease groups compared with controls (unadjusted P < 0.05). Together, these shared alterations suggest common taxonomic and metabolic disruptions that may be present in the prodromal stage (iRBD) through the symptomatic stage of LBD.
Decreased metabolic potential for starch degradation in LBD compared with cohabitant controls
Among pathways that were less abundant in LBD compared with controls, Starch degradation III displayed the strongest statistical significance (unadjusted P = 0.012, mixed-effects linear regression; Figure 2B). Its reduced abundance in LBD suggests a diminished microbial capacity for starch breakdown, which could contribute to lower SCFAs production.
Using HUMAnN3, we quantified species-attributed contributions to the Starch degradation III pathway in LBD patients and their cohabitant controls (Figure 4A). Specifically, this represents the proportion of total pathway abundance that can be assigned to individual microbial species based on their encoded gene families, thereby quantifying how much each species functionally contributes to the pathway. Substantial inter-individual variability was observed in the type of taxa contributing to the pathway, as well as in their relative contributions; however, Blautia obeum and Blautia wexlerae consistently accounted for the largest proportion of pathway activity across groups.
Figure 4
When comparing total species-attributed contributions within each household, the controls exhibited higher overall pathway contribution in 15 households (60%), whereas LBD patients showed higher contribution in only 10 households (40%) (Figure 4B). This household-level pattern highlights a consistent shift toward reduced starch-degrading capacity in LBD, even when controlling for shared environmental and dietary factors. At the individual species level, Eubacterium ramulus was the only species that contributed significantly more to this pathway in controls than in LBD patients (unadjusted P = 0.012; Figure 4C).
Finally, mapping gene families to specific enzymatic steps within the pathway indicated that three of the four steps had a greater number of gene families with higher abundance in controls than in LBD patients (Figure 4D). This step-level reduction, when considered alongside lower pathway abundance and diminished contributions from key carbohydrate-utilizing taxa, such as Eubacterium ramulus, points to a coordinated loss of functional capacity across multiple layers of the pathway. Together, these multi-layered analyses (pathway abundance, species-level contributions, and gene family distributions) provide convergent evidence for a reduction in microbial starch degradation potential in LBD, which may contribute to impaired fermentation of complex carbohydrates and reduced SCFAs availability.
Decreased metabolic potential for L-histidine degradation in both LBD and iRBD compared with cohabitant controls
The L-histidine degradation III pathway was one of the few metabolic pathways consistently reduced in both LBD and iRBD compared with their cohabitant controls. Given its role in producing L-glutamate, which is a microbial precursor for γ-aminobutyric Acid (GABA) synthesis, we further examined its contributing taxa and associated gene families to better characterize its functional decline. This pathway was significantly more abundant in cohabitant controls than in both LBD and iRBD groups (unadjusted P = 0.027 and unadjusted P = 5.47×10–5, respectively). Household-level stacked bar plots revealed substantial inter-individual variability in contributing species; however, Bacteroides uniformis emerged as the most consistently detected contributor across all groups (Figure 5A).
Figure 5
At the species level, several Bacteroides taxa showed significantly reduced contributions to L-histidine degradation III in disease groups. Contributions from Bacteroides eggerthii, Bacteroides uniformis, and Bacteroides dorei were significantly lower in LBD (unadjusted P < 0.05; Figures 5B–D), while Bacteroides thetaiotaomicron showed a reduction in iRBD (unadjusted P < 0.05; Figure 5E). These patterns parallel the decreased abundance of multiple Bacteroides species observed in both LBD and iRBD, suggesting that loss of Bacteroides-mediated metabolism may underlie the overall reduction in L-histidine degradation capacity.
The L-histidine degradation III pathway comprises several enzymatic steps that generate intermediates feeding into glutamate and downstream GABA production (Figure 5F). Mapping gene families to these steps demonstrated a broad reduction in functional potential across the pathway; at nearly every enzymatic step, a greater number of gene families showed higher abundance in controls than in either disease group. In iRBD, reductions were evident in earlier steps of the pathway, whereas in LBD, deficits extended further downstream, particularly within folate-dependent transformations (e.g., MTHFD1). This progressive erosion of gene family abundance from iRBD to LBD suggests a stepwise decline in microbial histidine catabolism along the disease continuum, potentially reducing microbial contributions to glutamate and GABA availability in the gut.
One of the final products of the L-histidine degradation III pathway is L-glutamate, which can enter the GABA synthesis pathway (GABA shunt) and be converted to GABA, which is one of the main inhibitory neurotransmitters of the brain. The key enzyme catalyzing this conversion is glutamate decarboxylase (EC 4.1.1.15). To investigate the microbial potential for this enzymatic step, we compared the abundance of gene families mapped to glutamate decarboxylase between LBD or iRBD and their cohabitant controls. In LBD, two differentially abundant gene families mapped to this enzyme showed higher abundance in controls, whereas in iRBD, one of three differentially abundant gene families showed higher abundance in controls (Figures 6A, B). Taken together, these findings suggest a progressive reduction in microbial glutamate decarboxylase potential from the prodromal stage (iRBD) to the symptomatic stage (LBD), which may reflect declining microbial capacity to support gut-derived GABA production along the disease continuum.
Figure 6
Correlation between clinical measures and microbial features
We calculated Spearman partial correlation coefficients (ρ) and corresponding P-values between four clinical features (CDR-SB, MoCA, STMS, and MDS-UPDRS III) and the abundances of microbial species and metabolic pathways (Supplementary Tables 13–16). Only features meeting both the primary significance criteria (unadjusted P < 0.05 and |ρ| > 0.4) and resampling-based validation (unadjusted empirical P < 0.05) are reported (Figures 7A, B; Supplementary Tables 17 and 18). Among these clinical scores, CDR-SB, MoCA, and STMS reflect dementia severity and cognitive performance: higher CDR-SB values indicate more severe dementia, whereas lower MoCA and STMS scores indicate greater cognitive impairment. MDS-UPDRS III measures motor symptom severity, with higher scores indicating worse motor function.
Figure 7
In LBD patients, several microbial metabolic pathways were positively associated with better cognitive performance and less severe dementia. The Superpathway of L-alanine biosynthesis and the Bifidobacterium shunt were positively correlated with both MoCA and STMS scores. Homolactic fermentation, Glycolysis I (from fructose 6-phosphate), and Glycolysis II (from fructose 6-phosphate) were negatively correlated with CDR-SB, and 4-amino-2-methyl-5-diphosphomethylpyrimidine biosynthesis II was positively correlated with MoCA (Figure 7A). All pathways described above also showed negative correlations with MDS-UPDRS III scores (Supplementary Table 15). At the species level, Erysipelatoclostridium ramosum was positively correlated with CDR-SB, while Evtepia gabavorous was negatively correlated with CDR-SB. Lentihominibacter faecis was positively correlated with STMS (Figure 7A). Unlike the pathway-level findings, these species associations were each specific to a single clinical measure.
In iRBD patients, several microbial features also correlated with better cognitive outcomes and less severe dementia. Clostridium leptum, Neglectibacter timonensis, and Starch degradation III were negatively correlated with CDR-SB (Figure 7B). Additionally, (S)-propane-1,2-diol degradation was positively correlated with both MoCA and STMS (Figure 7B).
Discussion
In this study, we identified alterations in gut microbiome taxonomy and functional potential across the LBD continuum, spanning the prodromal stage of iRBD through clinical manifestation of LBD. By integrating shotgun metagenomic sequencing with a household-matched study design, our analysis moves beyond community-level descriptions to reveal convergent taxonomic and metabolic potential shifts that are not apparent from global diversity metrics alone. Collectively, our findings suggest that subtle but consistent changes in microbial functions are shared across both the prodromal and symptomatic stages of the LBD continuum and may be relevant to disease progression through gut-brain axis mechanisms. These findings are presented as exploratory and hypothesis-generating, and are intended to provide specific, testable targets for future validation studies in larger cohorts.
At the level of overall community structure, neither α-diversity nor β-diversity differed significantly between patients and cohabitant controls. These findings are consistent with prior studies that utilized household-matched designs and support that broad measures of diversity may be insensitive to disease-associated changes when environmental exposure and dietary confounders are carefully controlled (Song et al., 2013; ; Teigen et al., 2024). Importantly, however, the absence of global diversity shifts does not imply functional equivalence. Rather, it highlights a potential strength of our study: the ability to detect fine-grained, functionally relevant microbial alterations that may occur within otherwise stable community structures.
To highlight the most informative microbial alterations, the discussion focuses on microbial species and pathways that showed both statistical significance (unadjusted P < 0.05) and meaningful effect size (|Cohen’s d| > 0.4). Microbial features lacking species-level characterization are not interpreted in depth. In some cases, features with more modest rankings are included because they show convergent signals across the LBD continuum or have established relevance to the disease.
One finding from our differential abundance analyses is the increased abundance of microbial taxa potentially associated with impaired gut barrier integrity. We observed higher abundance of Hungatella hathewayi (unadjusted P = 0.015, Cohen’s d = 0.477) in LBD and Collinsella aerofaciens (unadjusted P < 0.001, Cohen’s d = 1.202) in iRBD compared with their cohabitant controls. H. hathewayi has been shown to degrade glycosaminoglycans (GAGs) (), which are key structural components supporting gut epithelial barrier integrity (). Experimental studies have suggested that excessive depletion of GAG-rich layers increases tissue permeability and weakens barrier function (). In this context, increased abundance of this species may therefore be associated with altered gut integrity and low-grade inflammation. Similarly, C. aerofaciens and the Collinsella genus have previously been reported to be more abundant in RBD (), early PD (), and DLB (). Beyond neurodegenerative disease, Collinsella has also been found to be more abundant in rheumatoid arthritis, where it has been shown to stimulate IL-17A production, which is a cytokine linked to increased intestinal permeability and promotes chronic inflammation (). Therefore, our findings suggest that the increased abundance of these taxa in LBD and iRBD may be associated with gut barrier dysfunction and systemic inflammatory signaling, though this requires functional confirmation in future studies.
In addition, we observed higher abundance of microbial pathways associated with putative LPS biosynthesis in both disease groups. Specifically, the ADP-L-glycero-β-D-manno-heptose biosynthesis pathway (unadjusted P = 0.046, Cohen’s d = 0.461) was enriched in LBD, while the O-antigen building blocks biosynthesis pathway (unadjusted P = 0.002, Cohen’s d = 1.407) was elevated in iRBD. Both pathways contribute essential components to LPS assembly: ADP-L-glycero-β-D-manno-heptose serves as a precursor for the inner core (; Whitfield et al., 2020), and O-antigen building blocks form the outer polysaccharide domain of LPS (; Whitfield et al., 2020). Elevated microbial potential for LPS biosynthesis may be relevant to α-synucleinopathies through several proposed mechanisms. First, LPS has been reported to be associated with compromised gut epithelial integrity and increased intestinal permeability (; ; ), thereby potentially facilitating exposure of the enteric nervous system to microbial antigens (van Kessel and El Aidy, 2019). Moreover, in vitro studies suggest that LPS can promote α-synuclein aggregation in the gut (; Bhattacharyya and Bhunia, 2021). Overall, the increased abundance of LPS biosynthesis related pathways observed in both LBD and iRBD may be relevant to mechanisms that have been associated with heightened intestinal permeability and potential α-synuclein aggregation. One hypothesis is that this aggregation in the gut may facilitate “prion-like” propagation of α-synuclein through the enteric nervous system and ascend to the brainstem through the vagus nerve, as proposed in Braak’s staging schema (Visanji et al., 2013). However, whether this reflects active LPS production or downstream consequences for neurodegeneration remains to be directly tested.
We next examined microbial taxa and functional pathways associated with putative SCFA production, which are particularly sensitive to dietary and environmental influences. SCFAs have been reported to play important roles in maintaining gut barrier integrity (), exerting anti-inflammatory effects (), and regulating neurotransmission (; ; ). We identified a significant reduction in the starch degradation III pathway (unadjusted P = 0.012, Cohen’s d = –0.417) in LBD, accompanied by decreased abundance of multiple gene families across several enzymatic steps, suggesting reduced microbial capacity for starch degradation. The terminal products of this pathway feed into glycolysis and are metabolized to pyruvate, which can be further converted to acetate (; ). Complementing this finding, the β-(1,4)-mannan degradation pathway (unadjusted P = 0.013, Cohen’s d = –0.593) also had lower abundance in LBD group compared with their cohabitant controls. β-(1,4)-mannan is a hemicellulose that functions as a dietary fiber (). The end products of this pathway also feed into glycolysis and can be converted to acetate (). Acetate has been reported to cross the blood-brain barrier and may serve as a preferred energy source for astrocytes with proposed roles in supporting glutamate and GABA homeostasis (Waniewski and Martin, 1998; ). Within the brain, acetate is converted to acetyl-CoA, a key substrate for histone acetylation that has been linked to memory-related gene expression (; ). Thus, reduced complex carbohydrate degradation potential in LBD may be associated with reduced acetate availability with possible downstream implications for compromised metabolic and epigenetic processes along the gut–brain axis, although this remains to be directly tested. In addition, several known SCFA-producing taxa, including Roseburia hominis (unadjusted P = 0.001, Cohen’s d = –0.692), Eubacterium ramulus (unadjusted P = 0.006, Cohen’s d = –0.569), and Intestinimonas butyriciproducens (unadjusted P = 0.013, Cohen’s d = –0.533), were decreased in LBD, while Bacteroides were reduced in both LBD (unadjusted P = 0.033, Cohen’s d = –0.491) and iRBD (unadjusted P = 0.001, Cohen’s d = –1.305). The consistent loss of these beneficial commensals across disease stages suggests a shared microbial pattern that may be associated with weakened gut barrier integrity and disrupted gut-brain signaling from the prodromal through symptomatic phases.
Finally, we identified differences in microbial pathways related to putative neurotransmitter metabolism, particularly the L-histidine degradation III pathway, which was reduced in both LBD (unadjusted P = 0.027, Cohen’s d = –0.498) and iRBD (unadjusted P < 0.001, Cohen’s d = –1.921) relative to controls. Multiple gene families involved in this pathway showed lower abundance in disease groups. In this pathway, L-histidine is converted to L-glutamate, which can enter the GABA shunt pathway and be converted to GABA by glutamate decarboxylase. Among the significant gene families mapped to glutamate decarboxylase, both in LBD and one of three in iRBD showed higher abundance in controls; this suggests a possible progressive decline in microbial potential to convert L-glutamate to GABA from iRBD to LBD. Increasing evidence suggests that gut-derived GABA can influence circulating and central nervous system GABA levels, suggesting gut microbiome as a potential contributor to GABA-related neural signaling (Strandwitz, 2018; ). Consistent with this, previous studies have reported increased abundance of GABA-consuming taxa and enrichment of GABA degradation pathways in Parkinson’s Disease (Wallen et al., 2022; ). Additionally, a recent study has found that probiotic and prebiotic intervention restored brain GABA levels and reduced motor impairment and neuronal loss in animal models (Zhong et al., 2023). Furthermore, impaired GABA neurotransmission has been implicated in REM sleep without atonia, the neurophysiologic substrate of RBD (). Taken together, the reduced abundance of microbial pathways and gene families related to L-glutamate production and its conversion to GABA in both iRBD and LBD may link to alterations in microbial functional potential associated with GABA-related metabolic processes, which could have implications for disrupted excitatory-inhibitory balance along the gut-brain axis. These changes could be associated with RBD symptom development or neurodegenerative processes across the LBD spectrum, although further mechanistic studies are needed to clarify these relationships.
Several limitations should be acknowledged when interpreting the results. First, we were unable to adjust for sex in the statistical models due to strong collinearity between sex and case-control status, which is a common challenge in studies of LBD and iRBD given their marked male predominance. As a result, including sex as a covariate in mixed-effects models would not meaningfully disentangle sex effects from disease status and could obscure true disease-associated microbial signals. Second, although cohabitant controls were included to minimize environmental confounding, we did not directly quantify the degree of shared dietary intake within households. Consequently, residual dietary influences may not have been fully accounted for. Nevertheless, the use of cohabitant controls remains a major strength of this study, as it substantially reduces unmeasured environmental and lifestyle variability compared with designs relying on unrelated controls. Third, the relatively small sample size, particularly for the iRBD group and their matched controls (n = 10), limits the statistical power to detect small or moderate effects and is unlikely to be fully representative of the broader population. Future studies with larger, multi-center cohorts will be necessary to validate and extend these observations. Despite these constraints, the consistent directional patterns observed across both LBD and iRBD relative to their household-matched controls lend confidence that the identified alterations reflect biologically relevant signals rather than chance findings. Fourth, we did not apply multiple testing correction in the primary analysis, as very few results remained statistically significant after Benjamini–Hochberg false discovery rate adjustment given the large number of simultaneous comparisons across species, pathways, and gene families. Consequently, the reported results are subject to an increased risk of false-positive findings and should be interpreted as hypothesis-generating and exploratory rather than biologically conclusive (). Nevertheless, our results can provide a focused set of microbial species and pathways that may serve as candidates for validation in larger, independent cohorts (). Fifth, a key limitation of this study is the compositional nature of relative abundance data in microbiome analyses (). Total sum scaling forces all microbial features to sum to a constant value (i.e., 100%), meaning that changes in one taxon, pathway, or gene family influences the relative abundance of others. As a result, some observed differences may reflect compositional artifacts rather than true biological changes. In this study, we applied an arcsine square-root transformation to improve data distribution and stabilize variance for statistical analyses. While this transformation helps meet statistical assumptions, it does not resolve the underlying compositional issue, and our results should therefore be interpreted with appropriate caution. Finally, this study is cross-sectional, which limits the ability to infer causal relationships between gut microbial alterations and disease status. The observed associations with functions such as gut barrier integrity, neuroinflammation, and α-synuclein aggregation should therefore be interpreted with caution. Nevertheless, the consistent alterations identified across two stages of the LBD continuum provide a framework for future longitudinal and mechanistic studies. To this point, the candidate signals reported here define a clear and specific agenda for validation in larger, independent cohorts using complementary methodologies such as targeted qPCR, metatranscriptomics, and metabolomics.
Conclusions
Our findings suggest that the LBD disease continuum may be associated with changes in specific microbial taxa and metabolic pathways, rather than global disruption of the overall gut microbiome. These changes include reduced abundance of pathways related to complex carbohydrate fermentation in LBD and neurotransmitter-related metabolism in both LBD and iRBD, as well as increased abundance of taxa potentially associated with gut barrier disruption and pathways related to LPS biosynthesis in both disease groups. Overall, our findings support the possibility that gut microbial changes may accompany neurodegenerative processes along the disease continuum. By focusing on biologically interpretable taxonomic and metabolic features rather than broad diversity metrics alone, this study provides a framework for future longitudinal and mechanistic studies investigating potential links between the gut microbiome and LBD pathogenesis and progression.
Statements
Data availability statement
Sequencing data for stool metagenomes used in this study have been deposited at NCBI’s SRA data repository (PRJNA1393457) and can be downloaded without any restrictions at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1393457/. The deposited sequences include .fastq files for 70 stool metagenomes collected from 25 LBD patients, 10 iRBD patients, as well as their 35 cohabitant controls. Human reads were identified and removed prior to data upload. Code used for data analysis is available in the following GitHub repository: https://github.com/xiaowei-zhao-1111/LBD_Gut_Microbiome_2026.
Ethics statement
The studies involving humans were approved by Mayo Clinic Institutional Review Board 18-011778. 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. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
XZ: Data curation, Formal analysis, Methodology, Investigation, Visualization, Writing – original draft, Writing – review & editing. SJM: Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Investigation, Writing – review & editing. VKG: Investigation, Methodology, Writing – review & editing. KMG: Writing – review & editing. EKS: Writing – review & editing. KK: Writing – review & editing. RS: Writing – review & editing. MH: Writing – review & editing. HEV: Writing – review & editing. CS: Writing – review & editing. BFB: Funding acquisition, Project administration, Supervision, Writing – review & editing. OAR: Funding acquisition, Project administration, Supervision, Writing – review & editing. LMT: Conceptualization, Investigation, Project administration, Resources, Supervision, Writing – review & editing. JS: Conceptualization, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study was funded by NIH P30AG062677, U01NS100620, R34AG056639, U19AG071754, the University of Minnesota Data Science Initiative, the Mangurian Foundation for Lewy body disease research, The Little Family Foundation, the Turner Foundation (BFB), Ted Turner and Family and the Robert H. and Clarice Smith and Abigail van Buren Alzheimer Disease Research Program.
Acknowledgments
The authors wish to thank dedicated Mayo Clinic Alzheimer’s Disease Research Center (ADRC) staff for their valuable assistance. We are grateful to our patients and caregivers for their participation in our detailed annual assessments.
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.
The author/s ES, RS, CS declared that they were an editorial board member of Frontiers at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declared that generative AI was used in the creation of this manuscript. During the preparation of this work, the authors utilized ChatGPT (GPT5, OpenAI) during the final editing stages to improve word choice and readability. After using this tool, the authors reviewed and edited the content as needed and take full responsibility for the content of the publication.
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/frmbi.2026.1834726/full#supplementary-material
References
1
AbelesS. R.JonesM. B.Santiago-RodriguezT. M.LyM.KlitgordN.YoosephS.et al. (2016). Microbial diversity in individuals and their household contacts following typical antibiotic courses. Microbiome4, 39. doi: 10.1186/s40168-016-0187-9
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. doi: 10.1186/s13024-021-00427-6
3
BeghiniF.McIverL. J.Blanco-MíguezA.DuboisL.AsnicarF.MaharjanS.et al. (2021). Integrating taxonomic, functional, and strain-level profiling of diverse microbial communities with bioBakery 3. eLife10, e65088. doi: 10.7554/eLife.65088
4
Ben-ShacharM. S.LüdeckeD.MakowskiD. (2020). effectsize: Estimation of effect size indices and standardized parameters. J. Open Source Software5, 2815. doi: 10.21105/joss.02815
5
BertaniB.RuizN. (2018). Function and biogenesis of lipopolysaccharides. EcoSal Plus8, 10.1128. doi: 10.1128/ecosalplus.ESP-0001-2018
6
BhattacharyyaD.MohiteG. M.KrishnamoorthyJ.GayenN.MehraS.NavalkarA.et al. (2019). Lipopolysaccharide from gut microbiota modulates α-synuclein aggregation and alters its biological function. ACS Chem. Neurosci.10, 2229–2236. doi: 10.1021/acschemneuro.8b00733
7
BhattacharyyaD.MohiteG. M.KrishnamoorthyJ.GayenN.MehraS.NavalkarA.et al. (2019). Lipopolysaccharide from gut microbiota modulates α-synuclein aggregation and alters its biological function. ACS Chem. Neurosci.10, 2229–2236. doi: 10.1021/acschemneuro.8b00733
8
Blanco-MíguezA.BeghiniF.CumboF.McIverL. J.ThompsonK. N.ZolfoM.et al. (2023). Extending and improving metagenomic taxonomic profiling with uncharacterized species using MetaPhlAn 4. Nat. Biotechnol.41, 1633–1644. doi: 10.1038/s41587-023-01688-w
9
BoeveB. F.SilberM. H.FermanT. J.KokmenE.SmithG. E.IvnikR. J.et al. (1998). REM sleep behavior disorder and degenerative dementia. Neurology51, 363–370. doi: 10.1212/WNL.51.2.363
10
BoeveB. F.SilberM. H.ParisiJ. E.DicksonD. W.FermanT. J.BenarrochE. E.et al. (2003). Synucleinopathy pathology and REM sleep behavior disorder plus dementia or parkinsonism. Neurology61, 40–45. doi: 10.1212/01.WNL.0000073619.94467.B0
11
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics30, 2114–2120. doi: 10.1093/bioinformatics/btu170
12
BörnigenD.MorganX. C.FranzosaE. A.RenB.XavierR. J.GarrettW. S.et al. (2013). Functional profiling of the gut microbiome in disease-associated inflammation. Genome Med.5, 65. doi: 10.1186/gm469
13
BragaJ. D.ThongngamM.KumrungseeT. (2024). Gamma-aminobutyric acid as a potential postbiotic mediator in the gut–brain axis. NPJ Sci. Food8, 16. doi: 10.1038/s41538-024-00253-2
14
BrooksP. L.PeeverJ. H. (2011). Impaired GABA and glycine transmission triggers cardinal features of rapid eye movement sleep behavior disorder in mice. J. Neurosci.31, 7111–7121. doi: 10.1523/jneurosci.0347-11.2011
15
CaspiR.BillingtonR.FulcherC. A.KeselerI. M.KothariA.KrummenackerM.et al. (2018). The MetaCyc database of metabolic pathways and enzymes. Nucleic Acids Res.46, D633–D639. doi: 10.1093/nar/gkx935
16
ChenJ.WrightK.DavisJ. M.JeraldoP.MariettaE. V.MurrayJ.et al. (2016). An expansion of rare lineage intestinal microbes characterizes rheumatoid arthritis. Genome Med.8, 43. doi: 10.1186/s13073-016-0299-7
17
ChesselD.DufourA. B.ThioulouseJ. (2004). The ade4 package-I-One-table methods. R. News4, 5–10.
18
CryanJ. F.O'RiordanK. J.CowanC. S. M.SandhuK. V.BastiaanssenT. F. S.BoehmeM.et al. (2019). The microbiota-gut-brain axis. Physiol. Rev.99, 1877–2013. doi: 10.1152/physrev.00018.2018
19
DeleuS.MachielsK.RaesJ.VerbekeK.VermeireS. (2021). Short chain fatty acids and its producingorganisms: An overlooked therapy for IBD? eBioMedicine66, 103293. doi: 10.1016/j.ebiom.2021.103293
20
DixonP. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci.14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
21
DonaghyP. C.McKeithI. G. (2014). The clinical characteristics of dementia with Lewy bodies and a consideration of prodromal diagnosis. Alzheimer's Res. Ther.6, 46. doi: 10.1186/alzrt274
22
ErnyD.Hrabě de AngelisA. L.PrinzM. (2017). Communicating systems in the body: how microbiota and microglia cooperate. Immunology150, 7–15. doi: 10.1111/imm.12645
23
FacchinS.BertinL.BonazziE.LorenzonG.De BarbaC.BarberioB.et al. (2024). Short-chain fatty acids and human health: From metabolic pathways to current therapeutic implications. Life14, 559. doi: 10.3390/life14050559
24
FermanT. J.BoeveB. F.SmithG. E.SilberM. H.KokmenE.PetersenR. C.et al. (1999). REM sleep behavior disorder and dementia. Neurology52, 951. doi: 10.1212/WNL.52.5.951
25
FlintH. J.ScottK. P.DuncanS. H.LouisP.ForanoE. (2012). Microbial degradation of complex carbohydrates in the gut. Gut Microbes3, 289–306. doi: 10.4161/gmic.19897
26
FrancisK. L.ZhengH. B.SuskindD. L.MurphreeT. A.PhanB. A.QuahE.et al. (2024). Characterizing the human intestinal chondroitin sulfate glycosaminoglycan sulfation signature in inflammatory bowel disease. Sci. Rep.14, 11839. doi: 10.1038/s41598-024-60959-x
27
FrostG.SleethM. L.Sahuri-ArisoyluM.LizarbeB.CerdanS.BrodyL.et al. (2014). The short-chain fatty acid acetate reduces appetite via a central homeostatic mechanism. Nat. Commun.5, 3611. doi: 10.1038/ncomms4611
28
FujiwaraH.HasegawaM.DohmaeN.KawashimaA.MasliahE.GoldbergM. S.et al. (2002). α-Synuclein is phosphorylated in synucleinopathy lesions. Nat. Cell Biol.4, 160–164. doi: 10.1038/ncb748
29
GhoshT. S.DasM.JefferyI. B.O'TooleP. W. (2020). Adjusting for age improves identification of gut microbiome alterations in multiple diseases. Elife9, e50240. doi: 10.7554/elife.50240
30
GloorG. B.MacklaimJ. M.Pawlowsky-GlahnV.EgozcueJ. J. (2017). Microbiome datasets are compositional: And this is not optional. Front. Microbiol.8, 2224. doi: 10.3389/fmicb.2017.02224
31
GropperS. S.SmithJ. L.GroffJ. L. (2009). Advanced Nutrition and Human Metabolism (5th ed.). (Belmont, CA, USA: Wadsworth Cengage Learning), pp. 110.
32
GuoS.Al-SadiR.SaidH. M.MaT. Y. (2013). Lipopolysaccharide causes an increase in intestinal tight junction permeability in vitro and in vivo by inducing enterocyte membrane expression and localization of TLR-4 and CD14. Am. J. Pathol.182, 375–387. doi: 10.1016/j.ajpath.2012.10.014
33
GuoS.NighotM.Al-SadiR.AlhmoudT.NighotP.MaT. Y. (2015). Lipopolysaccharide regulation of intestinal tight junction permeability is mediated by TLR4 signal transduction pathway activation of FAK and MyD88. J. Immunol.195, 4999–5010. doi: 10.4049/jimmunol.1402598
34
HöglB.StefaniA.VidenovicA. (2018). Idiopathic REM sleep behaviour disorder and neurodegeneration — an update. Nat. Rev. Neurol.14, 40–55. doi: 10.1038/nrneurol.2017.157
35
HuM. T. (2020). REM sleep behavior disorder (RBD). Neurobiol. Dis.143, 104996. doi: 10.1016/j.nbd.2020.104996
36
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. doi: 10.1038/s41467-023-38248-4
37
HurstR. E.Van GordonS.TylerK.KroppB.TownerR.LinH.et al (2016). In the absence of overt urothelial damage, chondroitinase ABC digestion of the GAG layer increases bladder permeability in ovariectomized female rats. Am. J. Physiol. Renal Physiol.310, F1074–F1080. doi: 10.1152/ajprenal.00566.2015
38
IizumiT.BattagliaT.RuizV.PerezG. I. P. (2017). Gut microbiome and antibiotics. Arch. Med. Res.48, 727–734. doi: 10.1016/j.arcmed.2017.11.004
39
IranzoA.SantamariaJ.TolosaE. (2009). The clinical and pathophysiological relevance of REM sleep behavior disorder in neurodegenerative diseases. Sleep Med. Rev.13, 385–401. doi: 10.1016/j.smrv.2008.11.003
40
KimS. (2015). ppcor: an R package for a fast calculation to semi-partial correlation coefficients. Commun. For. Stat. Appl. Methods22, 665. doi: 10.5351/csam.2015.22.6.665
41
KneidingerB.MaroldaC.GraningerM.ZamyatinaA.McArthurF.KosmaP.et al. (2002). Biosynthesis pathway of ADP-L-glycero-beta-D-manno-heptose in Escherichia coli. J. Bacteriol.184, 363–369. doi: 10.1128/jb.184.2.363-369.2002
42
KuznetsovaA.BrockhoffP. B.ChristensenR. H. B. (2017). lmerTest package: tests in linear mixed effects models. J. Stat. Software82, 1–26. doi: 10.18637/jss.v082.i13
43
KwonD.ZhangK.PaulK. C.FolleA. D.Del RosarioI.JacobsJ. P.et al. (2024). Diet and the gut microbiome in patients with Parkinson's disease. NPJ Parkinson's Dis.10, 89. doi: 10.1038/s41531-024-00681-7
44
LangmeadB.SalzbergS. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods9, 357–359. doi: 10.1038/nmeth.1923
45
LiH.HandsakerB.WysokerA.FennellT.RuanJ.HomerN.et al. (2009). The Sequence Alignment/Map format and SAMtools. Bioinformatics25, 2078–2079. doi: 10.1093/bioinformatics/btp352
46
Lloyd-PriceJ.Abu-AliG.HuttenhowerC. (2016). The healthy human microbiome. Genome Med.8, 51. doi: 10.1186/s13073-016-0307-y
47
LohJ. S.MakW. Q.TanL. K. S.NgC. X.ChanH. H.YeowS. H.et al. (2024). Microbiota–gut–brain axis and its therapeutic applications in neurodegenerative diseases. Signal Transduction Targeted Ther.9, 37. doi: 10.1038/s41392-024-01743-1
48
LouisP.FlintH. J. (2017). Formation of propionate and butyrate by the human colonic microbiota. Environ. Microbiol.19, 29–41. doi: 10.1111/1462-2920.13589
49
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. doi: 10.3389/fnagi.2022.875261
50
MarzoukN. H.RashwanH. H.El-HadidiM.RamadanR.MysaraM. (2025). Proinflammatory and GABA eating bacteria in Parkinson's disease gut microbiome from a meta-analysis perspective. NPJ Parkinson's Dis.11, 145. doi: 10.1038/s41531-025-00950-z
51
MayàG.IranzoA.GaigC.Sánchez-ValleR.SerradellM.Molina-PorcelL.et al. (2024). Post-mortem neuropathology of idiopathic rapid eye movement sleep behaviour disorder: a case series. Lancet Neurol.23, 1238–1251. doi: 10.1016/s1474-4422(24)00402-2
52
MayerE. A. (2011). Gut feelings: the emerging biology of gut–brain communication. Nat. Rev. Neurosci.12, 453–466. doi: 10.1038/nrn3071
53
MayerE. A.TillischK.GuptaA. (2015). Gut/brain axis and the microbiota. J. Clin. Invest.125, 926–938. doi: 10.1172/jci76304
54
McKeithI. G.DicksonD. W.LoweJ.EmreM.O'BrienJ. T.FeldmanH.et al. (2005). Diagnosis and management of dementia with Lewy bodies. Neurology65, 1863–1872. doi: 10.1212/01.wnl.0000187889.17253.b1
55
McKeithI. G.BoeveB. F.DicksonD. W.HallidayG.TaylorJ. P.WeintraubD.et al. (2017). Diagnosis and management of dementia with Lewy bodies: Fourth consensus report of the DLB Consortium. Neurology89, 88–100. doi: 10.1212/wnl.0000000000004058
56
McKeithI. G.FermanT. J.ThomasA. J.BlancF.BoeveB. F.FujishiroH.et al. (2020). Research criteria for the diagnosis of prodromal dementia with Lewy bodies. Neurology94, 743–755. doi: 10.1212/wnl.0000000000009323
57
MewsP.DonahueG.DrakeA. M.LuczakV.AbelT.BergerS. L. (2017). Acetyl-CoA synthetase regulates histone acetylation and hippocampal memory. Nature546, 381–386. doi: 10.1038/nature22405
58
MinatoT.MaedaT.FujisawaY.TsujiH.NomotoK.OhnoK.et al. (2017). Progression of Parkinson's disease is associated with gut dysbiosis: Two-year follow-up study. PloS One12, e0187307. doi: 10.1371/journal.pone.0187307
59
NearingJ. T.DouglasG. M.HayesM. G.MacDonaldJ.DesaiD. K.AllwardN.et al. (2022). Microbiome differential abundance methods produce different results across 38 datasets. Nat. Commun.13, 342. doi: 10.1038/s41467-022-28034-z
60
NighotM.Al-SadiR.GuoS.RawatM.NighotP.WattersonM. D.et al. (2017). Lipopolysaccharide-Induced Increase in Intestinal Epithelial Tight Permeability Is Mediated by Toll-Like Receptor 4/Myeloid Differentiation Primary Response 88 (MyD88) Activation of Myosin Light Chain Kinase Expression. Am. J. of Pathol.187, 2698–2710. doi: 10.1016/j.ajpath.2017.08.005
61
NishiwakiH.UeyamaJ.KashiharaK.ItoM.HamaguchiT.MaedaT.et al. (2022). Gut microbiota in dementia with Lewy bodies. NPJ Parkinson's Dis.8, 169. doi: 10.1038/s41531-022-00428-2
62
OlsonE. J.BoeveB. F.SilberM. H. (2000). Rapid eye movement sleep behaviour disorder: demographic, clinical and laboratory findings in 93 cases. Brain123, 331–339. doi: 10.1093/brain/123.2.331
63
ParkerA.FonsecaS.CardingS. R. (2020). Gut microbes and metabolites as modulators of blood-brain barrier integrity and brain health. Gut Microbes11, 135–157. doi: 10.1080/19490976.2019.1638722
64
PengL.LiZ. R.GreenR. S.HolzmanI. R.LinJ. (2009). Butyrate enhances the intestinal barrier by facilitating tight junction assembly via activation of AMP-activated protein kinase in Caco-2 cell monolayers. J. Nutr.139, 1619–1625. doi: 10.3945/jn.109.104638
65
PetrovV.SaltykovaI.ZhukovaI.AlifirovaV.ZhukovaN.DorofeevaY. B.et al (2017). Analysis of gut microbiota in patients with Parkinson's disease. Bull. Exp. Biol. Med.162, 734–737. doi: 10.1007/s10517-017-3700-7
66
PostumaR. B.BergD. (2019). Prodromal parkinson's disease: the decade past, the decade to come. Movement Disord.34, 665–675. doi: 10.1002/mds.27670
67
RawatP. S.LiY.ZhangW.MengX.LiuW. (2022). Hungatella hathewayi, an efficient glycosaminoglycan-degrading Firmicutes from human gut and its chondroitin ABC exolyase with high activity and broad substrate specificity. Appl. Environ. Microbiol.88, e0154622. doi: 10.1128/aem.01546-22
68
Rios-CovianD.SalazarN.GueimondeM.de los Reyes-GavilanC. G. (2017). Shaping the metabolism of intestinal Bacteroides population through diet to improve human health. Front. Microbiol.8. doi: 10.3389/fmicb.2017.00376
69
RothmanK. J. (1990). No adjustments are needed for multiple comparisons. Epidemiology1, 43–46. doi: 10.1097/00001648-199001000-00010
70
RothschildD.WeissbrodO.BarkanE.KurilshikovA.KoremT.ZeeviD.et al. (2018). Environment dominates over host genetics in shaping human gut microbiota. Nature555, 210–215. doi: 10.1038/nature25973
71
SampsonT. R.MazmanianS. K. (2015). Control of brain development, function, and behavior by the microbiome. Cell. Host Microbe17, 565–576. doi: 10.1016/j.chom.2015.04.011
72
Savulescu-FiedlerI.BeneaS. N.CăruntuC.NancoffA. S.HomentcovschiC.BucuricaS. (2025). Rewiring the Brain Through the Gut: Insights into Microbiota-Nervous System Interactions. Curr. Issues Mol. Biol.47, 489. doi: 10.3390/cimb47070489
73
SegainJ. P.Raingeard de la BlétièreD.BourreilleA.LerayV.GervoisN.RosalesC.et al. (2000). Butyrate inhibits inflammatory responses through NFkappaB inhibition: implications for Crohn's disease. Gut47, 397–403. doi: 10.1136/gut.47.3.397
74
SharonG.SampsonT. R.GeschwindD. H.MazmanianS. K. (2016). The central nervous system and the gut microbiome. Cell.167, 915–932. doi: 10.1016/j.cell.2016.10.027
75
SolimanM. L.RosenbergerT. A. (2011). Acetate supplementation increases brain histone acetylation and inhibits histone deacetylase activity and expression. Mol. Cell. Biochem.352, 173–180. doi: 10.1007/s11010-011-0751-3
76
SongS. J.LauberC.CostelloE. K.LozuponeC. A.HumphreyG.Berg-LyonsD.et al. (2013). Cohabiting family members share microbiota with one another and with their dogs. Elife2, e00458. doi: 10.7554/eLife.00458
77
StrandwitzP. (2018). Neurotransmitter modulation by the gut microbiota. Brain Res.1693, 128–133. doi: 10.1016/j.brainres.2018.03.015
78
TavalireH. F.ChristieD. M.LeveL. D.TingN.CreskoW. A.BohannanB. J. M. (2021). Shared Environment and Genetics Shape the Gut Microbiome after Infant Adoption. mBio12, e00548-21. doi: 10.1128/mBio.00548-21
79
TeigenL. M.McCarterS. J.ZiegertZ.StaleyC.GrantK. M.GuptaV. K.et al. (2024). Taxonomic intestinal microbiota differences in Lewy body spectrum disease and cohabitant controls. Parkinsonism Related Disord.129, 107176. doi: 10.1016/j.parkreldis.2024.107176
80
The-Human-Microbiome-Project-Consortium (2012). Structure, function and diversity of the healthy human microbiome. Nature486, 207–214. doi: 10.1038/nature11234
81
TurnbaughP. J.HamadyM.YatsunenkoT.CantarelB. L.DuncanA.LeyR. E.et al. (2009). A core gut microbiome in obese and lean twins. Nature457, 480–484. doi: 10.1038/nature07540
82
TurnbaughP. J.GordonJ. I. (2009). The core gut microbiome, energy balance and obesity. J. Physiol.587, 4153–4158. doi: 10.1113/jphysiol.2009.174136
83
van KesselS. P.El AidyS. (2019). Bacterial metabolites mirror altered gut microbiota composition in patients with Parkinson's disease. J. Parkinsons Dis.9, S359–S370. doi: 10.3233/jpd-191780
84
VisanjiN. P.BrooksP. L.HazratiL. N.LangA. E. (2013). The prion hypothesis in Parkinson's disease: Braak to the future. Acta Neuropathol. Commun.1, 2. doi: 10.1186/2051-5960-1-2
85
WalkerZ.PossinK. L.BoeveB. F.AarslandD. (2015). Lewy body dementias. Lancet386, 1683–1697. doi: 10.1016/S0140-6736(15)00462-6
86
WallenZ. D.DemirkanA.TwaG.CohenG.DeanM. N.StandaertD. G.et al. (2022). Metagenomics of Parkinson's disease implicates the gut microbiome in multiple disease mechanisms. Nat. Commun.13, 6958. doi: 10.1038/s41467-022-34667-x
87
WaniewskiR. A.MartinD. L. (1998). Preferential utilization of acetate by astrocytes is attributable to transport. J. Neurosci.18, 5225–5233. doi: 10.1523/jneurosci.18-14-05225.1998
88
WhitfieldC.WilliamsD. M.KellyS. D. (2020). Lipopolysaccharide O-antigens-bacterial glycans made to measure. J. Biol. Chem.295, 10593–10609. doi: 10.1074/jbc.REV120.009402
89
YangL.ChenJ. (2022). A comprehensive evaluation of microbial differential abundance analysis methods: current status and potential solutions. Microbiome10, 130. doi: 10.1186/s40168-022-01320-0
90
ZhernakovaA.KurilshikovA.BonderM. J.TigchelaarE. F.SchirmerM.VatanenT.et al. (2016). Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science352, 565–569. doi: 10.1126/science.aad3369
91
ZhongH.-J.WangS.-Q.ZhangR.-X.ZhuangY.-P.LiL.YiS.-Z.et al. (2023). Supplementation with high-GABA-producing Lactobacillus plantarum L5 ameliorates essential tremor triggered by decreased gut bacteria-derived GABA. Trans. Neurodegeneration12, 58. doi: 10.1186/s40035-023-00391-9
92
ZhuL.LiuW.AlkhouriR.BakerR. D.BardJ. E.QuigleyE. M.et al. (2014). Structural changes in the gut microbiome of constipated patients. Physiol. Genomics46, 679–686. doi: 10.1152/physiolgenomics.00082.2014
Summary
Keywords
dementia, gut microbiome, gut-brain axis, isolated REM sleep behavior disorder (iRBD), Lewy body disease (LBD), mild cognitive impairment (MCI), shotgun metagenomic sequencing, α-synucleinopathy
Citation
Zhao X, McCarter SJ, Gupta VK, Grant KM, St. Louis EK, Kantarci K, Savica R, Hill M, Vuong HE, Staley C, Boeve BF, Ross OA, Teigen LM and Sung J (2026) Shotgun metagenomic analysis reveals taxonomic and functional alterations in the gut microbiome across prodromal and symptomatic Lewy body disease. Front. Microbiomes 5:1834726. doi: 10.3389/frmbi.2026.1834726
Received
19 March 2026
Revised
12 June 2026
Accepted
16 June 2026
Published
15 July 2026
Volume
5 - 2026
Edited by
Abhijeet Singh, Swedish University of Agricultural Sciences, Sweden
Reviewed by
Qiwen Cheng, Shandong University, China
Ovidiu Dumitru Ilie, “Ion Ionescu de la Brad” Iasi University of Life Sciences, Romania
Updates
Copyright
© 2026 Zhao, McCarter, Gupta, Grant, St. Louis, Kantarci, Savica, Hill, Vuong, Staley, Boeve, Ross, Teigen and Sung.
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: Jaeyun Sung, sung.jaeyun@mayo.edu
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.