Abstract
Background:
Many observational studies have investigated the link between the gut microbiota and Alzheimer’s disease (AD), but the causality remains uncertain.
Objective:
This study aimed to evaluate the causal impact of gut microbiota on AD.
Methods:
A two-sample Mendelian randomization (MR) study was conducted employing summary data. Summary statistics for AD were from the latest genome-wide association study (cases and proxy cases: 85,934; controls: 401,577). Summary data for gut microbiota were acquired from MiBioGen consortium. Causal effect estimations primarily relied on the inverse variance weighting method along with the sensitivity analyses for testing for pleiotropy and heterogeneity. Additionally, reverse MR analyses were performed to examine potential reverse causality.
Results:
Seven gut microbiota were identified as associated with AD risk. Order
Conclusions:
This study unveiled a causal link between certain gut microbiota and AD, offering new insights for advancing clinical treatments.
INTRODUCTION
Alzheimer’s disease (AD) is a neurodegenerative disorder that accounts for 60–80% of cases of cognitive decline [1]. Regrettably, the etiology of AD is poorly understood [2]. Because the AD biomarkers testing is radioactive and invasive, the diagnosis of AD mainly relies on the assessment of clinical manifestations. Regarding treatment, despite substantial investments in the field, there remains a deficiency in effective medications capable of altering the pathological progression of AD [3]. Hence, further elucidating the pathological mechanisms of AD, exploring novel non-invasive early biomarkers, and developing effective treatment strategies are imperative.
Gut microbiota are the microorganisms in the digestive tracts of animals that play a crucial role in metabolic processes, immune modulation, and intestinal epithelial upkeep [4]. Lately, considerable attention has been directed towards the influence of gut microbiota on central nervous system function, commonly known as the gut-brain axis [5]. Growing evidence suggests that changes in the gut microbiome are linked to the initiation and progression of AD [6–9]. In an observational study, the gut microbiome composition in individuals with AD was found to be significantly different from that of asymptomatic controls. By characterizing the bacterial taxonomic makeup of fecal samples from subjects with and without diagnosed AD, Vogt et al. identified differences in bacterial abundance from phylum to genus in the microbiomes of AD participants, including Firmicutes Phylum decreased, Bacteroidetes increased, and Bifidobacteria decreased [10]. It was also discovered that some gut microbiota in individuals with AD differed at the taxonomic level from controls, such as Bacteroidetes, Actinobacteria, Ruminococcus, Lachnospiraceae, and Lunamonales, indicating potential involvement of gut microbiota in AD pathogenesis [11]. Furthermore, neuroinflammation in AD patients may be related to amyloidosis and peripheral inflammation caused by altered gut microbiota. Specifically, AD patients exhibited elevated levels of species of endobacteria with pro-inflammatory properties such as Escherichia/Shigella, alongside reduced levels of the anti-inflammatory microbial group E. rectale. [12]. In addition, gut dysbiosis can contribute to the occurrence and onset of AD through oxidative stress and insulin resistance [13]. However, the debate persists over whether variations in gut microbiota composition are a causative factor or a consequence of AD, with meta-analyses of observational studies yielding inconclusive findings [11]. The results in these studies may be affected by unmeasured confounding factors (e.g., health conditions), potentially leading to inaccuracies in reflecting the true causal effects.
Mendelian randomization (MR) is a methodological study design, comparable to a randomized controlled trial, that estimates causal effects by exploiting the random allocation of genetic variants [14]. By leveraging strongly associated single nucleotide polymorphisms (SNPs) in the role of instrumental variables (IVs) for inference, MR produces unbiased estimates even when there is unmeasured confounding factors. While the MR approach has been applied across various studies to investigate potential causal links between the gut microbiota and various diseases [15–18], evidence linking the gut microbiota to a causal role in AD is limited.
In this case, we conducted a two-sample bidirectional MR investigation employing summary level statistics obtained from the most extensive and current genome-wide association studies (GWASs) of common gut microbiota species and AD. By using multiple MR analysis methods and performing comprehensive sensitivity analyses, we found potential causal effects of the presence of 7 gut microbiota species on AD proxied by 1,947 SNPs. Specifically, Order
MATERIALS AND METHODS
MR and assumptions
This study utilized two-sample MR with genetic instruments to explore potential causal links between gut microbiota and AD, leveraging data from the latest GWAS. SNPs associated with the composition of microorganisms inhabiting the gastrointestinal tract of humans were chosen as IVs. We employed a bidirectional design to test the association of gut microbiota with AD and to examine whether the gut microbiota composition difference among individuals is a downstream effect of AD.
MR relies on three assumptions: 1) Relevance: the IVs exhibit strong association with the exposure under scrutiny; 2) Effective random assignment: the genetic variants utilized as IVs are free from associations with the unmeasured potential confounders; 3) Exclusion restriction: the IVs solely influence the outcome by means of the exposure under investigation [19].
Genetic associations with gut microbiota
The genetic data regarding gut microbiota were sourced via the recent GWAS conducted by the international consortium MiBioGen [20, 21]. This comprehensive study incorporated sequencing data of 16 S ribosomal RNA genes and genotyping information from a total of 18,340 individuals across 24 cohorts from diverse nations, including Finland, Germany, the UK, etc. It aimed to investigate the connection involving autosomal human genetic variations and the gut microbiome, encompassing a comprehensive set of 211 taxa. Of these, 15 microbial taxa without particular species designation were excluded, leaving a total of 196 taxa for analysis.
The data source for AD
The summary statistics pertaining to AD were acquired via a recent GWAS carried out by the European Alzheimer and Dementia Biobank (EADB) consortium [22]. It encompasses a total of 39,106 individuals clinically diagnosed with AD, 46,828 individuals classified as proxy-ADD cases (individuals whose parent received a diagnosis of AD or dementia [23]) and 401,577 control individuals of European descent.
Genetic instruments selection
This study utilized SNPs that strongly and independently predicted exposures as genetic instruments. To broaden the scope of correlation findings, we chose to employ genetic instruments with a threshold of
Instrumental strength for the association between SNPs and exposure was assessed by F-statistic. F value is calculated by
Statistical analysis
Mendelian randomization analysis
In the main analysis, we evaluated the causal relationships between gut microbiota and AD by using the inverse variance weighting (IVW) estimation [29]. Besides, to enhance the robustness of the research, we utilized four alternative approaches, including MR-Egger regression, weighted median (WM), maximum likelihood estimation (MLE), as well as weighted mode (WMODE) method.
Pleiotropy denotes the occurrence where genetic variants impact multiple phenotypes, potentially influencing result reliability. IVW estimation involves meta-analyzing variant-specific Wald ratios for each variant, assuming balanced pleiotropy [28]. The four alternative approaches mentioned above can address horizontal pleiotropy—an infringement of the exclusion restriction assumption—under varying assumptions about the causal structure of the pleiotropy. The WM estimator requires that more than 50% SNPs are unaffected by pleiotropy [30]. Besides, in MR-Egger analysis, it is assumed that the magnitude of the pleiotropic effect is independent of the strength of the association between the genetic variant and the exposure across all instruments [28]. The maximum likelihood estimation stands as a classical tool for point estimation with lower standard error than the IVW method [31]. Finally, the weighted mode method allows a relaxation of the instrumental variable assumptions even if the majority of instruments are invalid [32]. As IVW estimation is reputed to be more robust when multiple IVs are used [30], this study primarily relies on the IVW method, complemented by the other four approaches as supplementary analyses.
Reverse Mendelian randomization
Bidirectional MR involves conducting two analyses on the same pair of phenotypes, with the exposure and the outcome variables swapped. In the reverse MR context, AD was regarded as the exposure and SNPs that reach significance
Sensitivity analyses
Given the challenge of verifying the exclusion restriction assumption, extensive sensitivity analyses were carried out to ensure the consistency and reliability of the inference. We employed the MR-Egger intercept test [31] and the MR-pleiotropy residual sum and outlier (MR-PRESSO) global test [33] to detect data points that deviate significantly from the rest of the dataset and assess horizontal pleiotropy within the IVs. In both tests, a
Results were considered significant based on the following criteria [34, 35]: 1) the IVW method demonstrated a statistically significant outcome (
We depicted the detailed schematic overview of this MR study in Fig. 1. All statistical procedures were executed utilizing R software version 4.3.2 and the R packages in the data availability statement section (‘TwosampleMR’[33], ‘Mendelian randomization’ [36], and ‘MR-PRESSO’ [37]).

Workflow for performing the two-sample bidirectional Mendelian randomization study. GWAS, genome-wide association study; MR, Mendelian randomization.
RESULTS
Collection of instrumental variable
After conducting quality control and removing exposures with fewer than 3 associated SNPs, we collected 1,947 SNPs jointly associated with AD and 183 taxa of gut microbiota, from fine to rough, comprising 111 genera, 29 families, 19 orders, 15 classes, and 9 phyla. Details of the SNPs and gut microbiota included are in Supplementary Table 1. The F-statistics of the IVs are all larger than 16.97, indicating a low likelihood of weak IVs bias.
Causal relationships between genetically proxied gut microbiota and AD
Figure 2 visually illustrates the preliminary findings regarding the relationships between the genetically proxied gut microbiota and the genetically proxied AD. Seven taxa of gut microbiota were identified as having a causal association with AD by applying the IVW method at the significance level 0.05 (Fig. 3). And the other four methods showed consistent directions with the IVW method.

The causal effects of gut microbiota on Alzheimer’s disease were detected by five MR analysis methods. The numbers in the outermost circle are the IDs of each gut microbiota, and their correspondence with the names of the microbiota taxa is shown in Supplementary Table 2. Each red dot represents the estimated causal impact of the gut microbiota on AD using the IVW method, that is, the odds ratio (OR). The dashed line indicates OR = 1. OR < 1 indicates that the microbiota may have a protective effect on AD, otherwise they may be harmful. The circular heat map shows the

Forest plot of MR estimate for the causal effect of seven gut microbiota on the risk of Alzheimer’s disease. OR, odds ratio; 95% CI, 95% confidence interval; IVW, inverse variance weighting estimation; WM, weighted median; MLE, maximum likelihood estimation; WMODE, weighted mode.
On the one hand, the IVW method showed three gut microbiota taxa had a positive causal effect on AD risk: Order
Sensitivity analyses of pleiotropy and heterogeneity
The comprehensive outcomes of the sensitivity assessments are presented in Supplementary Table 2. According to the MR-Egger regression, there was no indication that the IVs associated with the seven significant gut microbiota taxa had any pleiotropic effect on AD, which was coincided with the analysis of MR-PRESSO. Moreover, Cochran’s Q tests also showed no evidence for any heterogeneity in the analysis (Table 1). These results further showed the robustness of the significant causal effects of the seven specific gut microbiota species on AD.
Sensitivity analyses of MR results between gut microbiota and Alzheimer’s disease risk
IVW, inverse variance weighting estimation; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier.
Reverse Mendelian randomization
In order to rule out the possibility of reverse causal effects and further validate the previous results, we performed a reverse MR analysis employing the IVW method. After removing the SNPs in linkage disequilibrium, we identified 141 SNPs possessing F-statistics exceeding 10, as showed in Supplementary Table 3.
No statistically significant reverse causality was found between AD and each of these seven specific gut microbiota species, including Class
DISCUSSION
Our bidirectional MR study revealed potential causal effects of the presence of seven gut microbiota species on AD. Specifically, Order
Weighted median estimation is more robust than the IVW estimation if there exist pleiotropic IVs in MR studies. Based on this fact, we only focus on
To our knowledge, this study represents the initial exploration into the causality of gut microbiota and AD risk through MR approach. We uncovered a causal association between some common gut microbiota species and AD, with certain gut microbiota actively contributing to disease progression. A notable strength of our study lies in the rigorous implementation of MR methods, adeptly inferring the direction of causation, thereby yielding more precise conclusions about causal inferences.
While the MR design is generally less prone to confounding factors compared to other observational studies, there are limitations. First, relying on summary level statistics in GWAS databases instead of raw clinical data limits the capacity for further analyses on different forms of the disease. Second, because of the limited gut microbiota data sources, the same dataset is used for both selecting SNPs and estimating the impact size of the SNPs, resulting in potential selection bias. Third, the inclusion of primarily individuals with European descent in GWAS studies may lead to bias, as changes in the composition of the dominant intestinal microbiota, which are affected by different exposures such as dietary habits, may not be fully represented. Fourth, considering that AD typically manifests in older age, participants would need to have survived the exposure in order to be recruited, which may introduce bias. Finally, because the pathobiology of AD is complex, stringent multiple testing corrections might inadvertently overlook gut microbiota causally linked to AD. To address this concern, we opted not to apply multiple corrections.
In summary, this study unveiled a causal relationship between particular gut microbiota and the risk of AD. These findings have clinical implications, because they may provide novel avenues for advancing therapeutic approaches in clinical settings for AD and strategies to prevent AD. Whether AD has a causal impact on the composition of the gut microbiota needs more study.
AUTHOR CONTRIBUTIONS
Anqi Chen (Conceptualization; Formal analysis; Methodology; Writing – original draft); Yuquan Wang (Data curation; Formal analysis; Writing – original draft); Yue-Qing Hu (Conceptualization; Project administration; Supervision; Writing – review & editing).
Footnotes
ACKNOWLEDGMENTS
FUNDING
This research was supported partially by the National Key R&D Program of China [2023YFF1205101].
CONFLICT OF INTEREST
The authors have no conflict of interest to report.
DATA AVAILABILITY
The R package TwosampleMR is available at https://github.com/MRCIEU/TwoSampleMR. The R package MR-PRESSO is available at
. The source codes to reproduce the analyses can be accessed in the Supplementary Material.
