Skip to main content

Impact of the gut microbiota and associated metabolites on cardiometabolic traits, chronic diseases and human longevity: a Mendelian randomization study

Abstract

Features of the gut microbiota have been associated with several chronic diseases and longevity in preclinical models as well as in observational studies. Whether these relations underlie causal effects in humans remains to be established. We aimed to determine whether the gut microbiota influences cardiometabolic traits as well as the risk of chronic diseases and human longevity using a comprehensive 2-Sample Mendelian randomization approach. We included as exposures 10 gut-associated metabolites and pathways and 57 microbial taxa abundance. We included as outcomes nine cardiometabolic traits (fasting glucose, fasting insulin, systolic blood pressure, diastolic blood pressure, HDL cholesterol, LDL cholesterol, triglycerides, estimated glomerular filtration rate, body mass index [BMI]), eight chronic diseases previously linked with the gut microbiota in observational studies (Alzheimer’s disease, depression, type 2 diabetes, non-alcoholic fatty liver disease, coronary artery disease (CAD), stroke, osteoporosis and chronic kidney disease), as well as parental lifespan and longevity. We found 7 associations with evidence of causality before and after sensitivity analyses, but not after multiple testing correction (1198 tests). Most effect sizes (4/7) were small. The two largest exposure-outcome effects were markedly attenuated towards the null upon inclusion of BMI or alcohol intake frequency in multivariable MR analyses. While finding robust genetic instruments for microbiota features is challenging hence potentially inflating type 2 errors, these results do not support a large causal impact of human gut microbita features on cardiometabolic traits, chronic diseases or longevity. These results also suggest that the previously documented associations between gut microbiota and human health outcomes may not always underly causal relations.

Introduction

The human gut microbiota is the microbial symbiotic organ residing in the gut. It is involved in key metabolic and immunological processes including host immunity, food digestion, intestinal endocrine function and intestinal permeability [47]. Several observational studies revealed that the gut microbiota is associated with a wide range of cardiometabolic risk factors and human diseases [32]. The systemic effects of the gut microbiota is partly mediated through its by-products. These microbial metabolites can reach the peripheral circulation via the portal vein [6], or diffuse readily and be taken up by the gut mucosa [24], where they can reach organs and act as substrates or signalling molecules. This bidirectional crosstalk between the gut microbiota and different organs occurs via the gut-liver axis [2], the gut-brain axis [15], the gut-bone axis [66], the gut-kidney axis [22], the gut-lung axis [48] and the gut-heart axis [4]. Specific classes of microbiota-derived metabolites, notably short-chain fatty acids [19], branched-chain amino acids [3], trimethylamine N-oxide [58], and derivatives of tryptophan [5] have been implicated in the pathogenesis of metabolic disorders [1], lifespan [72], neurological and cardiovascular diseases [49].

Over the past few years, the gut microbiota emerged as a therapeutic target of great interest to prevent and/or treat chronic diseases and improve human lifespan and healthspan. An overwhelming amount of supportive evidence from preclinical models contributed to the widely accepted view that a large number of diseases and pathological processes could be influenced by the microbiome, from early metabolic perturbations to full-blown diseases and premature mortality. Fecal transplantation studies in rodents have provided promising results for the treatment of obesity [51], type 2 diabetes (T2D) [71], depression and chronic stress [42], liver injury [44], myocarditis [33] and aging [16]. Human microbiota-associated (HMA) studies, consisting of the transplantation of feces from human patients into germ-free mice while control mice receive feces from healthy humans, further supported these associations. A systematic review conducted in 2019 on the HMA method to study the impact of the microbiota on chronic diseases reported that 95% of such studies (36/38) concluded that fecal transplantation from a sick human donor resulted in at least one worsened symptom compared to healthy controls [69].

This finding was deemed “implausible” by the authors of this systematic review [69]. According to Walter et al., in the vast majority of cases, these studies lacked adequate replication and they had statistical and methodological flaws that artificially inflated the odds of obtaining positive findings [69]. Together with the “file-drawer effect” (whereby positive studies are more likely to be published compared to negative studies), these caveats may distort the odds of translating findings from preclinical models into microbiota-targeting therapies to prevent or treat human diseases. Observational studies in humans with various diseases have identified relevant differences in intestinal microbiota composition [59]. However, they are subject to biases such as reverse causality and confounding (through unmeasured confounders) and cannot, by design, assess causality. Obesity, pharmacotherapy, diet, alcohol intake and many other factors appear to be important confounders in the microbiota-health relationships [67]. Given these limitations, Walter et al. suggested that novel and innovative methods such as Mendelian randomization (MR) should be used to investigate the causal role of the gut microbiota in human disease etiology.

MR is an epidemiological approach that mitigates many of the biases of observational studies such as reverse causality or confounding. Under specific assumptions, it has the potential to evaluate potential causal effects between multiple exposures (gut microbiota features) and outcomes (cardiometabolic traits, chronic diseases and human longevity). Briefly, MR uses genetic variants strongly associated with an exposure (gut microbiota features) to infer causality with an outcome (cardiometabolic traits, chronic diseases or human longevity). Twin studies have shown that heritability of the abundance of different bacterial taxa is on average 20%, although some variation exists between taxa [26, 27]. This is consistent with the view that genes play a non-negligible role in determining gut microbiota composition, making MR a valuable tool to assess the potential causal role of the gut microbiota in human diseases.

Here, we used a 2-sample MR study design to investigate the potential causal links between gut microbiota features and nine cardiometabolic traits (fasting glucose, fasting insulin, diastolic blood, systolic blood pressure, HDL cholesterol, LDL cholesterol, triglycerides, estimated glomerular filtration rate, body mass index [BMI]) eight chronic disease outcomes encompassing different body systems (coronary artery disease [CAD], T2D, ischemic stroke [IS], nonalcoholic fatty liver disease [NAFLD], chronic kidney disease [CKD], osteoporosis, Alzheimer disease [AD] and depression), and human lifespan (as defined by parental lifespan and living beyond the 90th percentile). In this large-scale MR study, we first investigated the potential causal effect of micriobiota associated metabolites on diseases, metabolic risk factors and lifespan. Second, we leveraged summary statistics from two large genome-wide association studies (GWAS) of gut microbiota abundance to investigate the causal effect of genetically predicted taxa abundance on chronic diseases and human longevity.

Results

The conceptual framework of this MR analysis as well as the datasets used to derive the study exposures and outcomes are presented in Fig. 1 and Additional file 1: Table S1. Briefly, the objective of this MR analysis was to test the hypothesis that the gut microbiome causally impacts chronic diseases and longevity and to provide estimates for each exposure/outcome association. We performed two sample MR on microbiota features as exposures and relevant outcomes of cardiometabolic risk factors, chronic diseases and human longevity. We used publicly available genome-wide association study (GWAS) summary statistics to extract 67 traits related to the microbiome including 10 fecal and plasma metabolites associated with the gut microbiota, microbial abundance of 57 taxa partly under genetic control, nine cardiometabolic risk factors and 10 disease-related outcomes and human longevity (see Methods). The ten metabolites were selected based on the existence of taxa that metabolize them and based on their previous association with chronic diseases, as described in Additional file 1: Table S2. Analyses were restricted to participants from European ancestry except for a fraction of participants included in the study of microbial abundance (MiBioGen) consortium and CAD (CARDIoGRAMplusC4D). Samples from exposures and outcomes overlapped to a minor extent and at different degrees depending of data sources. We selected only exposures that had at least three independent (r2 < 0.01) genetic instruments at minimum p-value < 1e-5 (the threshold differed between data source depending on the availability of genetic instruments [Additional file 1: Table S3]) with mean F statistics  > 10, resulting in 67 microbiota-related exposures available for MR. These criteria were chosen to minimize weak instrument bias and allow the use of sensitivity analysis to assess the validity of the MR assumptions. The harmonized dataset of the associations between genetic variant and exposure, and between genetic variant and outcome is presented in Additional file 1: Table S4.

Fig. 1
figure 1

Overview of the Mendelian randomization framework used to investigate the causal effect of gut microbiome features (blood and gut-derived metabolites and microbial taxa abundance) on cardiometabolic traits, chronic diseases and human longevity

Effect of gut microbiome-related metabolites on chronic diseases and longevity

We first sought to determine whether 10 blood metabolites associated with the gut microbiota and fecal microbial metabolites or their functional pathways could influence cardiometabolic risk factors, chronic diseases and longevity. We selected genetic instruments for short-chain fatty acids such as fecal propionate, a gut metabolite linked with T2D in a recent MR study, from a GWAS of 898 participants [56] and acetate [36]. We included the microbial pathway involved in 4-aminobutanoate (GABA) degradation (PWY-5022 pathway) acting as a proxy for butyrate production by the gut [56]. We also included trimethylamine N-oxide [54], branched chain amino acids (leucine, isoleucine and valine) [46] and derivatives of tryptophan produced by the gut microbiota, namely indole 3-propionate [54], serotonin and kynurenin [46]. We performed IVW-MR for each of the health outcomes under study (Fig. 2). A total of 190 exposure-outcome associations were tested. The mean absolute effect size was 0.04, meaning that a 1-SD deviation increase in a microbiota associated metabolites would only increase by 0.04 SD a risk factor or by ~ 4% the risk of a disease. Seven associations passed a nominal p-value significance threshold of 0.05 (including the propionate-T2D association) but none of the gut microbiota metabolites were associated with chronic diseases and longevity after multiple testing correction. Figure 2 also reports the association of LDL cholesterol, used as positive control, with the outcomes of interest. As expected, each SD increase in LDL cholesterol was positively associated with cardiovascular diseases (OR = 1.5 95% CI 1.4–1.6, P = 1.4e-31) and negatively associated with human longevity (OR = 0.64 95% CI 0.56–0.74, P = 8.4e-11). Altogether, this analysis identified some blood and gut-derived metabolites that may be associated with cardiometabolic traits, chronic diseases and human longevity but their effect sizes were small. Spurious associations also cannot entirely be ruled out since none of the metabolites were associated with outcomes of interest after correction for multiple testing.

Fig. 2
figure 2

Balloon plot of the association of microbial fecal metabolites, microbial pathway and plasma metabolites with all 19 health outcomes. LDL cholesterol is included as positive control. Non-available (NA) associations stem from a lack of overlapping SNPs or proxies between exposure and outcome data resulting in fewer than three genetic instruments in the harmonized data set. Associations at P-value  > 0.05 are depicted with crosses. For readability, the effect of LDL on LDL was forced to be non-available

Effect of gut microbial abundance on chronic diseases and longevity

We explored the impact of different taxa abundance on health-related outcomes. We included available genetic information on 57 microbial taxa abundance from the two recent GWAS studies of Kurilshikov et al., and Rühlemann et al. [40, 55]. We then performed IVW-MR on each of the 19 health outcomes (Fig. 3). Similar to microbiota associated metabolites, the mean absolute effect size was 0.04. Out of 1008 exposure-outcome tests, 69 passed a significance threshold of 0.05, but no association remained after Benjamini-Hochberg correction for multiple testing. Altogether, this analysis identified some microbes that may be associated with cardiometabolic traits, chronic diseases and human longevity but, as observed with blood and gut-derived metabolites, reported effect sizes were small.

Fig. 3
figure 3

Balloon plot of the association of gut microbial taxa abundance with all 19 health outcomes. LDL cholesterol is included as positive control. Non-available (NA) associations stem from a lack of overlapping SNPs or proxies between exposure and outcome data resulting in fewer than three genetic instruments in the harmonized data set. Associations at P-value  > 0.05 are depicted with crosses. For readability, the effect of LDL on LDL was forced to be non-available

Exploration of promising findings and tests for pleiotropy

One of the key assumptions underlying MR is that genetic instruments affecting gut microbiota features do not affect diseases or longevity by other mechanisms than the one associated with gut microbiota features (Davies, Holmes, and Davey Smith 2018). This phenomenon is known as horizontal pleiotropy [43]. We performed robust MR analysis on all 76 potentially causal relationships (primary analysis P-value  < 0.05) to estimate the robustness of our primary causal estimates to pleiotropy. We used four different methods that make different assumptions about the nature of the underlying pleiotropy: the weighted median, the MR-PRESSO or when it could not be performed the MR-Radial and the contamination mixture approaches. Consistency across the estimates of the methods provides support to causality. The MR-Egger intercept was also used to assess robustness to horizontal pleiotropy. We removed the MR results that were not supported by robust MR analyses (weighted median method P > 0.05, MR-Egger intercept P < 0.05, MR PRESSO outliers-adjusted test or MR-Radial test (P > 0.05), contamination mixture (P > 0.05). Of the 76 associations tested, 7 remained as they were consistent with a true causal effect unlikely to be confounded by pleiotropy (Fig. 4 and Additional file 1: Table S6). Genetic instruments of these associations were located outside the ABO, HLA and APOE gene window (± 1 Mb), which may represent genetic regions potentially causing pleiotropy and hence were less likely to be pleiotropic. Genetic instruments did not display evidence of reverse causality as indicated by Steiger filtering tests. Most (4/7) causal effects were low (abs(b)  < 0.1). Notably, the causal effect of circulating serotonin levels on CAD and the causal effect of the order lactobacillales on T2D were strong (OR = 1.23 95% CI 1.07–1.42, P = 4.0e-03) and (1.26 95% CI 1.11–1.44, P = 3.6e-04) respectively. This analysis identified some gut microbiota associated features that may be worth exploring in further studies.

Fig. 4
figure 4

Forest plot of the associations that were consistent across robust MR analyses. Dichotomous traits are reported on a log(OR) scale. Continuous are reported on 1-SD scale. Dots depicts the point estimate. Horizontal bars depicts 95% confidence interval (CI)

Exploration of BMI and alcohol intake as potential confounding factors

Obesity and alcohol intake frequency were recently identified as major confounding factors in microbiome-disease associations [67]. We performed multivariable MR on the 7 promising associations to determine if the causal effects were robust to the inclusion of obesity and alcohol intake frequency. This analysis provided evidence that the two aforementioned largest reported associations might, to a certain extent, be confounded by BMI or alcohol consumption (Fig. 5 and Additional file 1: Table S7). For example, the causal effect of serotonin plasma level on CAD attenuated towards the null to (1.13 95% CI 0.94–1.36, P = 2.0e-01) upon inclusion of alcohol intake frequency as covariate and the causal effect of the order Lactobacillales on type 2 diabetes risk attenuated towards the null to (OR = 1.14 95% CI 0.95–1.37, P = 1.6e-01) upon inclusion as BMI as covariate.

Fig. 5
figure 5

IVW-MR results before and after correcting for BMI and alcohol intake frequency using multivariable MR framework

Discussion

In order to determine whether previously reported studies using preclinical models or observational study designs in humans were consistent with causal effects, we assessed the roles of a wide range of microbial factors and nine cardiometabolic risk factors, eight chronic diseases as well as human longevity using MR. We found 7 associations with evidence of causality before and after sensitivity analyses, but not after multiple testing correction.Most effect sizes were small. The causal effect of serotonin levels on CAD and of the Lactobacillales order abundance on T2D were strong, but their effect substantially decreased upon inclusion of BMI and alcohol intake frequency in multivariable MR analyses. Altogether, results of this study suggest that previously reported associations between the human gut microbiome and human disease might have been due to biases such as reverse causality or confounding and that the impact of the gut microbiota on cardiometabolic traits, chronic diseases and human longevity may not be as prominent as previously suggested.

Comparisons with other studies

Our results generally contrast with those from previous observational studies. Microbial metabolites have been associated with health and disease such as neurological disease, NAFLD, cardiovascular disease, survival and T2D (Agus, Clément, and Sokol 2020; [49, 72]. In a prospective study of apparently healthy participants with eight-year follow-up, elevated plasma TMAO (4th quartile vs. 1st quartile) predicted CAD even after adjustments for traditional risk factors (adjusted odds ratio (OR) 1.58 [95% confidence interval (CI) 1.21–2.06], P < 0.001) [62]. This result was similar in another study cohort [29]. In mice, increasing through dietary supplementation the level of TMAO accelerated atherosclerosis [39]. By contrast, our MR analysis does not support causality of TMAO on CAD (0.99 95% CI 0.98–1, P = 9.9e-02), despite precisely estimating the effect size, as previously suggested in another MR investigation on CAD [34]. Dietary factors could arguably act as confounding factors, since meat intake increases TMAO levels [39].

Several differences in the microbial composition of diseased and healthy individuals have been identified, but causality remains to be elucidated. RCT of fecal microbiome transfer (FMT) in humans are currently employed to establish causality between microbiome and health, but few have been attempted, and even fewer have been conclusive [20, 50]. To date, most successful randomized control study of FMT on humans has been applied to the treatment of recurrent or refractory Clostridioides difficile infections [73] and some to ulcerative colitis [17]. Mice FMTs are a valuable exploratory tool, but inference to human subjects is hazardous. Particularly, a substantial proportion of species in the human gut are not present in mice [37]. For example, several FMT in mice from lean to obese mice resulted in improved cardiometabolic profile [41, 76], but these findings failed to replicate in humans. A systematic review of all three randomized placebo-controlled studies to treat obesity published to date found no impact of FMT on obesity, fasting plasma glucose, hepatic insulin sensitivity, or cholesterol markers across all included studies [75]. For human observational studies, multiple confounding factors could create spurious correlation between microbiome and chronic diseases, including antibiotic use, age, sex, diet, geography, BMI and alcohol intake [38]. Moreover, alteration of the gut microbiota could potentially be a consequence of disease states rather than a causal factor [14].

Altogether, the large proportion of null findings (i.e., human gut-related traits may not cause chronic disease) is in line with recent literature showing an overwhelming positive publication bias in the microbiome literature [69]. Publication bias can occur if studies are not adequately corrected for multiple testing and can be identified with attempt and failure to replicate (i-e winners curse bias) [57]. For example, it was originally published that individuals with obesity were more likely to have lower bacterial diversity and relative abundances of the phylum Bacteroidetes [63], but this result failed to replicate in 9 independent cohorts [23, 61, 70]. This highlight the importance of triangulating with different methods such as MR to address a causal research question.

Strengths and limitations

An important strength of this study is the use of the MR design with the largest publicly available GWAS datasets. Because alleles are randomly assigned and fixed at conception, biases due to confounding and reverse causality are mitigated in an MR analysis [18]. A further strength is that sample was mostly restricted to individuals of European ancestry to reduce bias due to population stratification. However, it also restricts the generalizability of the results to this ethnic group.

Our study, however, has limitations. Robust genetic instruments for microbial species are challenging to find. First, microbiome heterogeneity and interindividual variability are high, substantially reducing the statistical power of microbiome GWAS analyses. Second, the phenotype is distal from individual genes making it a complex polymorphic trait with many variants of small effect size which could be prone to pleiotropy. Lastly, the twin heritability for gut microbiota taxa abundance is only on average 20% [26, 27]. This heritability estimate represents the upper ceiling that variance explained by genetic instruments can attain, reducing power. These factors all contributed to the fact that the number of mbQTLs identified to date is rather modest. For these reasons, we used a less stringent p-value cut-off to include a greater number of genetic variants to allow the use of sensitivity analysis and increase power. However, a less stringent p-value cut-off has the trade-off to potentially increase the chance of including false-positive effect variants which induce biases. The most important bias it introduces, the winner’s curse bias, refers to the fact that the genetic hits in discovery samples are more likely to be false positive, adding noise to the analysis which will typically bias MR results towards the null. Second in importance, the weak instrument bias occurs when the variance explained by the instrument and the sample size are low [13]. In the two sample MR setting, it will bias towards the null [11]. Third in importance, invalid instruments are pleiotropic variants that affect the outcome via another pathway than the one going through the exposure. Since genetic variants are derived from the host, instruments for the microbiome as exposure have a high chance of being pleiotropic. In general, pleiotropy is more likely to bias estimates away from the null [60]. Therefore, our null results are generally more robust to horizontal pleiotropy. We minimized weak instrument biases by including only exposures with genetic instruments with an average F-statistic above 10. We minimized the propensity of the results to be biased by pleiotropic variants by systematically including sensitivity analyses.

The Steiger test we performed can be biased in the presence of a difference in measurement error between the exposure and the outcome [30]. Typically non-differential measurement error will decrease variance explained hence will bias the Steiger test towards lower measurement error to higher measurement error directionality [30]. Since microbiota features (our exposures) are arguably measured less precisely than the health phenotype (our outcomes), one can assume that Steiger filtering will be biased towards the reverse direction. This bias should have little impact on our results since no genetic instruments were removed by Steiger filtering.A second limitation is that the microbiome GWAS included in the current analysis did not target the entire 16S RNA gene, which greatly diminished their ability to achieve a sufficient taxonomic resolution to identify potential therapeutic targets. The meta-analysis by Ruhlemann targeted the V1-V2 subregion while the meta-analysis by Kurilshikov et al. included mostly the V4 subregion and to a lesser extent the V1-V3 subregion. Targeting only 16S subregions such as V4 leads to lower taxonomic resolution achieved compared to sequencing the full V1-V9 16S RNA gene [35]. Indeed, using a variable region as a surrogate for the entire 16S RNA gene only allows for the identification of taxa at the genus level or above [35]. Being confident at the genus level provides little to inform disease treatment. Indeed, within a high-level taxon such as a phylum, some species may have a positive correlation with a disease, but some neutral or negative. For example, in our study, the phylum Actinobacteria was potentially protective for CAD, while the subsequent level, the class Actinobacteria was a potential risk factor for the same disease. Employing third-generation technologies has the potential to allow the sequencing of the full 16S RNA gene in a high throughput manner, and improve taxonomic discrimination.

A third limitation is the difference in ancestry from data in the MiBioGen and CARDIoGRAMplusC4D consortium that could potentially violate the independence assumption through population stratification. Population stratification occurs in the advent that population subgroups have different disease rates or different distributions of continuous traits and have different frequencies of alleles. This phenomenon could violate the independence MR assumption. However, this bias is likely to be minimal as principal components were included and participants were mostly from European descent.

Conclusions

Using MR, an approach less subject to reverse causality and confounding factors in comparison to traditional methods, we showed that several features of human gut microbioata including plasma metabolites and microbial taxa abundance had no evidence of causal effect on nine cardiometabolic traits, eight chronic diseases and human longevity. While finding robust genetic instruments for microbiota features is challenging potentially inflating type 2 errors, these results do not support large causal effects of the human gut microbiota and microbial metabolites on human chronic diseases and longevity. As the microbiome field matures, the use of larger microbiome GWAS study taking advantage of discriminatory potential of the full 16S RNA gene is warranted to fully elucidate the association of the human gut microbiota in the etiology of chronic diseases.

Methods

Study exposures (gut microbiota-derived metabolites)

We derived our ten gut microbiota-derived exposures of interest from five different publicly available data sources (Additional file 1: Table S1). We selected independent (r2 ≤ 0.01) SNPs (for all studies: P-value  < 1e-5; except Lotta et al.: P-value  < 1e-6) as genetic instruments. Genetic instruments for fecal propionate and PWY-5022 were obtained from a GWAS on 952 normoglycemic participants of the LifeLines-DEEP cohort, a population-based cohort from northern Netherlands (age ranges 18–84 years) [56]. Fecal propionate levels were measured by gas chromatography-mass spectrometry (GCMS). The functional pathway PWY-5022 was obtained with HUMAnN2 (v 0.4.0) [25] and MetaCyc metabolic-pathway database [64]. Genotyping was carried out with two Illumina arrays, HumanCytoSNP-12 BeadChip and ImmunoChip.

Genetic instruments for plasma TMAO and indole-3-propionate were extracted from a GWAS conducted in 2076 participants from European ancestry from the Framingham Heart Study (FHS) offspring cohort [54]. The FHS offspring cohort is a prospective community-based cohort from Framingham, Massachusetts, USA. Children of the spouse of the FHS study were recruited in 1971. Metabolites profiling was performed by liquid chromatography-mass spectrometry 63. Genotyping was conducted using the Affymetrix 500 K mapping array and the Affymetrix 50 K gene-focused MIP array. The participants all provided their informed consent and the study was approved by the Boston University Medical Center.

Genetic instruments for acetate were extracted from a meta-analysis of GWAS conducted on 10 European cohorts totalizing 24,925 individuals [36]. Human blood metabolites were quantified with quantitative high-throughput NMR metabolomics platform. Genetic instruments or plasma branched-chain amino acids (leucine, isoleucine and valine) and tryptophan derived metabolites kynurenine and serotonin were extracted from a meta-analysis of seven cohorts on up to 86,401 participants [46]. Human blood metabolites were quantified with quantitative high-throughput NMR metabolomics platform.

Study exposures (gut-microbiata abundance)

We first identified 4 recent GWAS on gut microbe abundance with available summary statistics [40, 45, 53, 55]. We filtered all microbiota quantitative trait loci (mbQTLs) with P-value  < 1.0e-6 for all taxa abundance present in their analysis and kept only exposures with at least three shared mbQTLs with mean F statistics  > 10. The study of Lopera-Maya et al., and Qin et al., were removed as none of the exposures satisfied our criteria (≥ 3 independent mbQTLs at a P-value  < 1e-6, with mean F-statistics  > 10). In total, we derived our microbiota taxa abundance exposures of interest from two different publicly available data sources (Additional file 1: Table S1). Genetic instruments for bacterial taxon were extracted from a GWAS of bacterial taxon abundances of 8,956 German individuals from the PopGen (population-based cohort), the FoCus (population registry based), the KORA FF4 (population-based adult cohort initiated in 1984) and the SHIP cohort (longitudinal population-based cohort)67. Human Genotyping and fecal microbial 16S rRNA gene surveys were performed using multiple arrays. Additionally, other genetic instruments for bacterial taxon were extracted from a meta-analysis conducted by the MiBioGen consortium on 16S fecal microbiome data from 18,340 individuals (24 cohorts) [40]. All cohorts implemented the standardized 16S processing pipeline that uses SILVA as a reference database, with truncation of the taxonomic resolution of the database to genus level. Cohorts were Middle Eastern, East Asian, American Hispanic/Latin, African American and admixed, although the majority of the sample (more than 72%) came from European descent. Notably, The two studies shared the PopGen (n = 721), the FOCUS (n = 960), and the SHIP (n = 1901) cohorts for a total of 3582 individuals or approximately half the samplesize of the study by ruhlemann et al. For this reason, we only sourced genetic instruments in the study by Ruhlemann et al., for exposures that were absent from Kurilshikov.

Study outcomes

We used publicly available GWAS summary statistics of the largest GWAS of nine cardiometabolic traits, eight chronic diseases, parental lifespan and longevity. Relevant information on the GWAS summary statistics are presented in Additional file 1: Table S1.

Selection of genetic variants and variants harmonization

We first identified all SNPs associated with exposures. Summary parameters for genetic instrument selection can be found at Additional file 1: Table S7. These SNPs were then clumped using the 1000Genomes Project Phase 3 European LD reference panel to make sure instrumental variables were independent with a 10 Mb window and pairwise linkage disequilibrium (LD) r2 < 0.01. This step was implemented with the gwwasvcf package in R [21]. Instrument strength was quantified using the F-statistic [13], and the variance explained was quantified using the r2 [52]. Variant harmonization was performed by aligning the betas of different studies on the same effect allele with the TwoSampleMR package [31]. When a particular exposure SNP was not present in the outcome dataset, we used proxy SNPs instead (r2 > 0.8). We used the LD matrix of the 1000 Genomes Project-European sample of the Utah residents from North and West Europe. We kept only the results based on at least three independent shared SNPs with mean F statistics  > 10 to reduce weak instrument bias and allow for robust MR analyses.

Primary Mendelian randomization analyses

We conducted primary MR analysis on each outcome and exposure association. As primary method for causal inference, we performed the IVW method with multiplicative random effects with a standard error correction for under dispersion as recommended by recent MR guidelines [10]. The IVW-MR combines the ratio estimates from each genetic instrument in a meta-analysis model by giving more weight to the ratio estimates with lower variance [12]. A total of 1198 primary analyses were performed: 67 exposures * 19 outcomes—75 exposures/outcomes with fewer than three overlapping genetic instruments or proxies. We applied a Benjamini Hochberg correction for multiple testing using a false discovery rate of 5% to reduce the propensity of false positive finding. For dichotomous traits, we transformed ORs and CIs to effect sizes and standard error when it was not already done. Lifespan, blood pressure, fasting glucose, fasting insulin, and glomerular filtration rate were reported in years, mmHG, log(pmol/L), log(pmol/L) and log(eGFR) respectively. For better interpretability and comparability, these summary statistics were transformed to a one standard deviation scale using the sdY.est function in the coloc package [68]. The other continuous variables were already inverse-rank normal transformed in the GWAS.

Sensitivity analyses

For associations with IVW-MR P-value  < 0.05, we performed sensitivity analyses to estimate the robustness of our primary causal estimate. We used 5 different robust methods that make different assumptions about the nature of the underlying pleiotropy. As a general test of the presence of pleiotropy, we used the intercept term from MR-Egger regression), MR-Egger is similar to the IVW method except the regression model estimates an intercept [7]. An intercept significantly different from zero gives indication of pleiotropy. The contamination mixture provides consistent estimate under the plurality valid assumption [60]59. It assumes that not all genetic variants are valid IVs and runs a likelihood function to categorize genetic instruments as valid or invalid. As a general test of robustness of the IVW-MR estimates, we used the weighted median. The weighted median estimates an unbiased causal effect if the “majority valid” assumption is upheld, that is if up to 50% of the weights comes from variants that are valid IVs [8]. Finally, as a general test to the presence of outliers, we used the outlier-robust method MR-PRESSO, which is a simulation approach where genetic variants are removed based on their contributions to heterogeneity [65]. This method provides consistent estimates under the same assumptions as the IVW-MR method for the set of genetic variants that are not identified as outliers [60]. When the MR-PRESSO could not be performed because there was less than four instruments, we performed MR-Radial as an outlier robust test instead. The MR-Radial uses a simulation-based approach to detect and remove outlier variants to re-estimate the exposure-outcome relation [9].

As additional sensitivity analysis, we excluded variants from known pleiotropic gene regions and performed Steiger filtering. Because of their known association with pleiotropic pathways, we excluded from the analysis all SNPs of the HLA, ABO and APOE genetic regions. We also performed Steiger filtering to remove variants with evidence of a stronger association with the outcome than its association with the exposure. The Steiger test provides a p-value under the null hypothesis that the difference in variance explained is null [30].

Obesity and alcohol intake frequency were recently identified as major confounding factors in the gut-disease associations [67] because they are both to some extent associated with the health outcome under study while potentially simultaneously influencing microbiome composition. Although other confounding factors may exist, adding BMI and alcohol intake frequency are the most important predictors of microbiota composition and health and adding them as covariates in linear mixed-effect models reduced the numbers of spurious microbiome health associations [67]. To account for this, we performed multivariable MR as a sensitivity analysis to correct for measured confounders [28]. BMI and alcohol intake frequency GWAS from the UK biobank were obtained from publicly accessible source (Additional file 1: Table S1). Multivariable MR IVW estimates were computed using the MendelianRandomization package [74].

Availability of data and materials

All data used in this study are in the public domain. Additional file 1: Table S1 describes the data used and relevant information to retrieve the summary statistics. Code was performed in the R V.4.0.0 computing environment using publicly accessible functions from the TwoSampleMR V.0.5.5 https://github.com/MRCIEU/TwoSampleMR, the MendelianRandomization V.0.5.1 https://cran.rproject.org/web/packages/MendelianRandomization/index.html and the data.table V.1.14.0 https://github.com/Rdatatable/data.table packages. The tidyverse V.1.3.1 collection of R packages was also used. Code to reproduce the analysis can be found on https://github.com/gagelo01/Dysbiose_project.

References

  1. Agus A, Clément K, Sokol H. Gut microbiota-derived metabolites as central regulators in metabolic disorders. Gut December Gutjnl. 2020. https://doi.org/10.1136/gutjnl-2020-323071.

    Article  Google Scholar 

  2. Albillos A, de Gottardi A, Rescigno M. The gut-liver axis in liver disease: pathophysiological basis for therapy. J Hepatol. 2020;72(3):558–77. https://doi.org/10.1016/j.jhep.2019.10.003.

    Article  CAS  Google Scholar 

  3. Arany Z, Neinast M. Branched Chain amino acids in metabolic disease. Curr DiabRep. 2018;18(10):76. https://doi.org/10.1007/s11892-018-1048-7.

    Article  CAS  Google Scholar 

  4. Bartolomaeus H, McParland V, Wilck N. Gut-heart axis : how gut bacteria influence cardiovascular diseases. Herz. 2020;45(2):134–41. https://doi.org/10.1007/s00059-020-04897-0.

    Article  Google Scholar 

  5. Beaumont M, Neyrinck AM, Olivares M, Rodriguez J, de Rocca A, Serra MR, Bindels LB, et al. The gut microbiota metabolite indole alleviates liver inflammation in mice. FASEB J. 2018;32(12):6681–93. https://doi.org/10.1096/fj.201800544.

    Article  CAS  Google Scholar 

  6. den Besten G, van Eunen K, Groen AK, Venema K, Reijngoud D-J, Bakker BM. The role of short-Chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. J Lipid Res. 2013;54(9):2325–40. https://doi.org/10.1194/jlr.R036012.

    Article  CAS  Google Scholar 

  7. Bowden J, Smith GD, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int J Epidemiol. 2015;44(2):512–25. https://doi.org/10.1093/ije/dyv080.

    Article  Google Scholar 

  8. Bowden J, Smith GD, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14. https://doi.org/10.1002/gepi.21965.

    Article  Google Scholar 

  9. Bowden J, Spiller W, Fabiola Del Greco M, Sheehan N, Thompson J, Minelli C, Smith GD. Improving the visualization, interpretation and analysis of two-sample summary data mendelian randomization via the radial plot and radial regression. Int J Epidemiol. 2018;47(4):1264–78. https://doi.org/10.1093/ije/dyy101.

    Article  Google Scholar 

  10. Stephen B, Smith GD, Davies NM, Dudbridge F, Dipender G, Maria Glymour M, Hartwig FP, et al. Guidelines for performing mendelian randomization investigations. Wellcome Open Res. 2020. https://doi.org/10.12688/wellcomeopenres.15555.2.

    Article  Google Scholar 

  11. Burgess S, Davies NM, Thompson SG. Bias due to participant overlap in two-sample mendelian randomization. Genet Epidemiol. 2016;40(7):597–608. https://doi.org/10.1002/gepi.21998.

    Article  Google Scholar 

  12. Burgess S, Small DS, Thompson SG. A review of instrumental variable estimators for mendelian randomization. Stat Methods Med Res. 2017;26(5):2333–55. https://doi.org/10.1177/0962280215597579.

    Article  Google Scholar 

  13. Burgess S, Thompson SG, Genetics Collaboration CRPCHD. Avoiding bias from weak instruments in mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64. https://doi.org/10.1093/ije/dyr036.

    Article  Google Scholar 

  14. Cani PD. Human gut microbiome: hopes, threats and promises. Gut. 2018;67(9):1716–25. https://doi.org/10.1136/gutjnl-2018-316723.

    Article  CAS  Google Scholar 

  15. Carabotti M, Scirocco A, Maselli MA, Severi C. The gut-brain axis: interactions between enteric microbiota, central and enteric nervous systems. Ann Gastroenterol Quart Publ Hellenic Soc Gastroenterol. 2015;28(2):203–9.

    Google Scholar 

  16. Yinfeng C, Zhang S, Zeng B, Zhao J, Yang M, Zhang M, Li Y, Ni Q, De W, Li Y. Transplant of microbiota from long-living people to mice reduces aging-related indices and transfers beneficial bacteria. Aging (Albany NY). 2020;12(6):4778–93. https://doi.org/10.18632/aging.102872.

    Article  Google Scholar 

  17. Costello SP, Waters O, Bryant RV, Katsikeros R, Makanyanga J, Schoeman M, Mountifield RE, et al. Short duration, low intensity, pooled fecal microbiota transplantation induces remission in patients with mild-moderately active ulcerative colitis: a randomised controlled trial. Gastroenterology. 2017;152(5):S198-9. https://doi.org/10.1016/S0016-5085(17)30969-1.

    Article  Google Scholar 

  18. Davies NM, Holmes MV, Smith GD. Reading mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018. https://doi.org/10.1136/bmj.k601.

    Article  Google Scholar 

  19. Vadder De, Filipe P-D, Goncalves D, Vinera J, Zitoun C, Duchampt A, Bäckhed F, Mithieux G. Microbiota-generated metabolites promote metabolic benefits via gut-brain neural circuits. Cell. 2014;156(1–2):84–96. https://doi.org/10.1016/j.cell.2013.12.016.

    Article  CAS  Google Scholar 

  20. Depommier C, Everard A, Druart C, Plovier H, Van Hul M, Vieira-Silva S, Falony G, et al. Supplementation with akkermansia muciniphila in overweight and obese human volunteers: a proof-of-concept exploratory study. Nat Med. 2019;25(7):1096–103. https://doi.org/10.1038/s41591-019-0495-2.

    Article  CAS  Google Scholar 

  21. Elsworth B, Lyon M, Alexander T, Liu Y, Matthews P, Hallett J, Bates P, et al. The MRC IEU opengwas data infrastructure. BioRxiv. 2020. https://doi.org/10.1101/2020.08.10.244293.

    Article  Google Scholar 

  22. Evenepoel P, Poesen R, Meijers B. The gut-kidney axis. Pediatr Nephrol. 2017;32(11):2005–14. https://doi.org/10.1007/s00467-016-3527-x.

    Article  Google Scholar 

  23. Finucane MM, Sharpton TJ, Laurent TJ, Pollard KS. A taxonomic signature of obesity in the microbiome? getting to the guts of the matter. PLoS ONE. 2014;9(1):e84689. https://doi.org/10.1371/journal.pone.0084689.

    Article  CAS  Google Scholar 

  24. Flint HJ. Gut microbial metabolites in health and disease. Gut Microbes. 2016;7(3):187–8. https://doi.org/10.1080/19490976.2016.1182295.

    Article  Google Scholar 

  25. Franzosa EA, McIver LJ, Rahnavard G, Thompson LR, Schirmer M, Weingart G, Lipson KS, et al. Species-level functional profiling of metagenomes and metatranscriptomes. Nat Method. 2018;15(11):962–8. https://doi.org/10.1038/s41592-018-0176-y.

    Article  CAS  Google Scholar 

  26. Goodrich JK, Davenport ER, Beaumont M, Jackson MA, Knight R, Ober C, Spector TD, Bell JT, Clark AG, Ley RE. Genetic determinants of the gut microbiome in UK Twins. Cell Host Microbe. 2016;19(5):731–43. https://doi.org/10.1016/j.chom.2016.04.017.

    Article  CAS  Google Scholar 

  27. Goodrich JK, Waters JL, Poole AC, Sutter JL, Koren O, Blekhman R, Beaumont M, et al. Human genetics shape the gut microbiome. Cell. 2014;159(4):789–99. https://doi.org/10.1016/j.cell.2014.09.053.

    Article  CAS  Google Scholar 

  28. Gormley M, Dudding T, Sanderson E, Martin RM, Thomas S, Tyrrell J, Ness AR, et al. A multivariable mendelian randomization analysis investigating smoking and alcohol consumption in oral and oropharyngeal cancer. Nat Commun. 2020;11(1):6071. https://doi.org/10.1038/s41467-020-19822-6.

    Article  CAS  Google Scholar 

  29. Guo F, Zhou J, Li Z, Zaixin Yu, Ouyang D. The association between trimethylamine N-oxide and its predecessors choline, L-carnitine, and betaine with coronary artery disease and artery stenosis. Cardiol Res Pract. 2020;2020:5854919. https://doi.org/10.1155/2020/5854919.

    Article  Google Scholar 

  30. Hemani G, Tilling K, Smith GD. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13(11):e1007081. https://doi.org/10.1371/journal.pgen.1007081.

    Article  CAS  Google Scholar 

  31. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, Laurin C, et al. The MR-base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408. https://doi.org/10.7554/eLife.34408.

    Article  Google Scholar 

  32. Hooks KB, O’Malley MA. Dysbiosis and its discontents. MBio. 2017. https://doi.org/10.1128/mBio.01492-17.

    Article  Google Scholar 

  33. Hu X-F, Zhang W-Y, Wen Q, Chen W-J, Wang Z-M, Chen J, Zhu F, et al. Fecal microbiota transplantation alleviates myocardial damage in myocarditis by restoring the microbiota composition. Pharmacol Res. 2019;139:412–21. https://doi.org/10.1016/j.phrs.2018.11.042.

    Article  Google Scholar 

  34. Jia J, Dou P, Gao M, Kong X, Li C, Liu Z, Huang T. Assessment of causal direction between gut microbiota-dependent metabolites and cardiometabolic health: a bidirectional mendelian randomization analysis. Diabetes. 2019;68(9):1747–55. https://doi.org/10.2337/db19-0153.

    Article  CAS  Google Scholar 

  35. Johnson JS, Spakowicz DJ, Hong B-Y, Petersen LM, Demkowicz P, Chen L, Leopold SR, et al. Evaluation of 16S RRNA gene sequencing for species and strain-level microbiome analysis. Nat Commun. 2019;10(1):5029. https://doi.org/10.1038/s41467-019-13036-1.

    Article  CAS  Google Scholar 

  36. Kettunen J, Demirkan A, Würtz P, Draisma HHM, Haller T, Rawal R, Vaarhorst A, et al. Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat Commun. 2016;7:11122. https://doi.org/10.1038/ncomms11122.

    Article  CAS  Google Scholar 

  37. Kieser S, Zdobnov EM, Trajkovski M. Comprehensive mouse microbiota genome catalog reveals major difference to its human counterpart. PLoS Comput Biol. 2022;18(3):e1009947. https://doi.org/10.1371/journal.pcbi.1009947.

    Article  CAS  Google Scholar 

  38. Kim D, Hofstaedter CE, Zhao C, Mattei L, Tanes C, Clarke E, Lauder A, et al. Optimizing methods and dodging pitfalls in microbiome research. Microbiome. 2017. https://doi.org/10.1186/s40168-017-0267-5.

    Article  Google Scholar 

  39. Koeth RA, Wang Z, Levison BS, Buffa JA, Org E, Sheehy BT, Britt EB, et al. Intestinal microbiota metabolism of L-carnitine, a nutrient in red meat, promotes atherosclerosis. Nat Med. 2013;19(5):576–85. https://doi.org/10.1038/nm.3145.

    Article  CAS  Google Scholar 

  40. Kurilshikov A, Medina-Gomez C, Bacigalupe R, Radjabzadeh D, Wang J, Demirkan A, Le Roy CI, et al. Large-scale association analyses identify host factors influencing human gut microbiome composition. Nature Genet. 2021. https://doi.org/10.1038/s41588-020-00763-1.

    Article  Google Scholar 

  41. Lai Z-L, Tseng C-H, Ho HJ, Cheung CKY, Lin J-Y, Chen Y-J, Cheng F-C, et al. Fecal microbiota transplantation confers beneficial metabolic effects of diet and exercise on diet-induced obese mice. Sci Rep. 2018;8(1):1–11. https://doi.org/10.1038/s41598-018-33893-y.

    Article  CAS  Google Scholar 

  42. Langgartner D, Vaihinger CA, Haffner-Luntzer M, Kunze JF, Weiss A-L, Foertsch S, Bergdolt S, Ignatius A, Reber SO. The role of the intestinal microbiome in chronic psychosocial stress-induced pathologies in male mice. Front Behav Neurosci. 2018. https://doi.org/10.3389/fnbeh.2018.00252.

    Article  Google Scholar 

  43. DA Lawlor, K Wade, MC Borges, T Palmer, FP Hartwig, G Hemani, J Bowden. 2019. ‘A Mendelian Randomization Dictionary: Useful Definitions and Descriptions for Undertaking, Understanding and Interpreting Mendelian Randomization Studies’. Preprint. Open Science Framework. https://doi.org/10.31219/osf.io/6yzs7.

  44. Liu Y, Fan L, Cheng Z, Lei Yu, Cong S, Yaxin Hu, Zhu L, et al. Fecal transplantation alleviates acute liver injury in mice through regulating Treg/Th17 cytokines balance. Sci Rep. 2021;11(1):1611. https://doi.org/10.1038/s41598-021-81263-y.

    Article  CAS  Google Scholar 

  45. Lopera-Maya EA, Kurilshikov A, van der Graaf A, Hu S, Andreu-Sánchez S, Chen L, Vich Vila A, et al. Effect of host genetics on the gut microbiome in 7,738 participants of the dutch microbiome project. BioRxiv. 2020. https://doi.org/10.1101/2020.12.09.417642.

    Article  Google Scholar 

  46. Lotta LA, Pietzner M, Stewart ID, Wittemans LBL, Li C, Bonelli R, Raffler J, et al. A cross-platform approach identifies genetic regulators of human metabolism and health. Nat Genet. 2021;53(1):54–64. https://doi.org/10.1038/s41588-020-00751-5.

    Article  CAS  Google Scholar 

  47. Lynch SV, Pedersen O. The human intestinal microbiome in health and disease. N Engl J Med. 2016;375(24):2369–79. https://doi.org/10.1056/NEJMra1600266.

    Article  CAS  Google Scholar 

  48. Marsland BJ, Trompette A, Gollwitzer ES. The gut-lung axis in respiratory disease. Ann Am Thorac Soc. 2015;12(Suppl 2):S150-156. https://doi.org/10.1513/AnnalsATS.201503-133AW.

    Article  Google Scholar 

  49. Martinez KB, Leone V, Chang EB. Microbial metabolites in health and disease: navigating the unknown in search of function. J Biol Chem. 2017;292(21):8553–9. https://doi.org/10.1074/jbc.R116.752899.

    Article  CAS  Google Scholar 

  50. Mocanu V, Zhang Z, Deehan EC, Kao DH, Hotte N, Karmali S, Birch DW, Samarasinghe KK, Walter J, Madsen KL. Fecal microbial transplantation and fiber supplementation in patients with severe obesity and metabolic syndrome: a randomized double-blind, placebo-controlled phase 2 trial. Nat Med. 2021;27(7):1272–9. https://doi.org/10.1038/s41591-021-01399-2.

    Article  CAS  Google Scholar 

  51. Pérez-Matute P, Íñiguez M, de Toro M, Recio-Fernández E, Oteo JA. Autologous fecal transplantation from a lean state potentiates caloric restriction effects on body weight and adiposity in obese mice. Sci Rep. 2020;10(1):9388. https://doi.org/10.1038/s41598-020-64961-x.

    Article  CAS  Google Scholar 

  52. Pierce BL, Ahsan H, VanderWeele TJ. Power and instrument strength requirements for mendelian randomization studies using multiple genetic variants. Int J Epidemiol. 2011;40(3):740–52. https://doi.org/10.1093/ije/dyq151.

    Article  Google Scholar 

  53. Qin Y, Havulinna AS, Liu Y, Jousilahti P, Ritchie SC, Tokolyi A, Sanders JG, et al. Combined effects of host genetics and diet on human gut microbiota and incident disease in a single population cohort. MedRxiv. 2020. https://doi.org/10.1101/2020.09.12.20193045.

    Article  Google Scholar 

  54. Rhee EP, Ho JE, Chen M-H, Shen D, Cheng S, Larson MG, Ghorbani A, et al. A genome-wide association study of the human metabolome in a community-based cohort. Cell Metab. 2013;18(1):130–43. https://doi.org/10.1016/j.cmet.2013.06.013.

    Article  CAS  Google Scholar 

  55. Rühlemann MC, Hermes BM, Bang C, Doms S, Moitinho-Silva L, Thingholm LB, Frost F, et al. Genome-wide association study in 8,956 german individuals identifies influence of ABO histo-blood groups on gut microbiome. Nature Genet. 2021. https://doi.org/10.1038/s41588-020-00747-1.

    Article  Google Scholar 

  56. Sanna S, van Zuydam NR, Mahajan A, Kurilshikov A, Vila AV, Võsa U, Mujagic Z, et al. Causal relationships between gut microbiome, short-chain fatty acids and metabolic diseases. Nat Genet. 2019;51(4):600–5. https://doi.org/10.1038/s41588-019-0350-x.

    Article  CAS  Google Scholar 

  57. Schloss PD. Identifying and overcoming threats to reproducibility, replicability, robustness, and generalizability in microbiome research. MBio. 2018. https://doi.org/10.1128/mBio.00525-18.

    Article  Google Scholar 

  58. Shan Z, Sun T, Huang H, Chen S, Chen L, Luo C, Yang W, et al. Association between microbiota-dependent metabolite trimethylamine-N-oxide and type 2 diabetes. Am J Clin Nutr. 2017;106(3):888–94. https://doi.org/10.3945/ajcn.117.157107.

    Article  CAS  Google Scholar 

  59. Shreiner AB, Kao JY, Young VB. The gut microbiome in health and in disease. Curr Opin Gastroenterol. 2015;31(1):69–75. https://doi.org/10.1097/MOG.0000000000000139.

    Article  CAS  Google Scholar 

  60. Slob EAW, Burgess S. A comparison of robust mendelian randomization methods using summary data. Genet Epidemiol. 2020;44(4):313–29. https://doi.org/10.1002/gepi.22295.

    Article  Google Scholar 

  61. Sze MA, Schloss PD. Looking for a signal in the noise: revisiting obesity and the microbiome. MBio. 2016. https://doi.org/10.1128/mBio.01018-16.

    Article  Google Scholar 

  62. Tang WHW, Li XS, Yuping Wu, Wang Z, Khaw K-T, Wareham NJ, Max Nieuwdorp S, Boekholdt M, Hazen SL. Plasma trimethylamine N-oxide (TMAO) levels predict future risk of coronary artery disease in apparently healthy individuals in the EPIC-norfolk prospective population study. Am Heart J. 2021;236:80–6. https://doi.org/10.1016/j.ahj.2021.01.020.

    Article  CAS  Google Scholar 

  63. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457(7228):480–4. https://doi.org/10.1038/nature07540.

    Article  CAS  Google Scholar 

  64. Vatanen T, Kostic AD, d’Hennezel E, Siljander H, Franzosa EA, Yassour M, Kolde R, et al. Variation in microbiome LPS immunogenicity contributes to autoimmunity in humans. Cell. 2016;165(4):842–53. https://doi.org/10.1016/j.cell.2016.04.007.

    Article  CAS  Google Scholar 

  65. Verbanck M, Chen C-Y, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8. https://doi.org/10.1038/s41588-018-0099-7.

    Article  CAS  Google Scholar 

  66. Villa CR, Ward WE, Comelli EM. Gut microbiota-bone axis. Crit Rev Food Sci Nutr. 2017;57(8):1664–72. https://doi.org/10.1080/10408398.2015.1010034.

    Article  Google Scholar 

  67. Vujkovic-Cvijin I, Sklar J, Jiang L, Natarajan L, Knight R, Belkaid Y. Host variables confound gut microbiota studies of human disease. Nature. 2020;587(7834):448–54. https://doi.org/10.1038/s41586-020-2881-9.

    Article  CAS  Google Scholar 

  68. Wallace C. Eliciting priors and relaxing the single causal variant assumption in colocalisation analyses. PLoS Genet. 2020;16(4):e1008720. https://doi.org/10.1371/journal.pgen.1008720.

    Article  CAS  Google Scholar 

  69. Walter J, Armet AM, Brett Finlay B, Shanahan F. Establishing or exaggerating causality for the gut microbiome: lessons from human microbiota-associated rodents. Cell. 2020;180(2):221–32. https://doi.org/10.1016/j.cell.2019.12.025.

    Article  CAS  Google Scholar 

  70. Walters WA, Zech Xu, Knight R. Meta-analyses of human gut microbes associated with obesity and IBD. FEBS Lett. 2014;588(22):4223–33. https://doi.org/10.1016/j.febslet.2014.09.039.

    Article  CAS  Google Scholar 

  71. Wang H, Yuan Lu, Yan Y, Tian S, Zheng D, Leng D, Wang C, Jiao J, Wang Z, Bai Y. Promising treatment for type 2 diabetes: fecal microbiota transplantation reverses insulin resistance and impaired islets. Front Cell Infect Microbiol. 2020. https://doi.org/10.3389/fcimb.2019.00455.

    Article  Google Scholar 

  72. Wilmanski T, Diener C, Rappaport N, Patwardhan S, Wiedrick J, Lapidus J, Earls JC, et al. Gut microbiome pattern reflects healthy ageing and predicts survival in humans. Nat Metab. 2021;3(2):274–86. https://doi.org/10.1038/s42255-021-00348-0.

    Article  CAS  Google Scholar 

  73. Wortelboer K, Nieuwdorp M, Herrema H. Fecal microbiota transplantation beyond clostridioides difficile infections. EBioMedicine. 2019;44:716–29. https://doi.org/10.1016/j.ebiom.2019.05.066.

    Article  Google Scholar 

  74. Yavorska OO, Burgess S. MendelianRandomization: An R package for performing mendelian randomization analyses using summarized data. Int J Epidemiol. 2017;46(6):1734–9. https://doi.org/10.1093/ije/dyx034.

    Article  Google Scholar 

  75. Zhang Z, Mocanu V, Cai C, Dang J, Slater L, Deehan EC, Walter J, Madsen KL. Impact of Fecal Microbiota Transplantation on Obesity and Metabolic Syndrome—A Systematic Review. Nutrients. 2019;11(10):2291. https://doi.org/10.3390/nu11102291.

    Article  CAS  Google Scholar 

  76. Zoll J, Read MN, Heywood SE, Estevez E, Marshall JPS, Kammoun HL, Allen TL, Holmes AJ, Febbraio MA, Henstridge DC. Fecal Microbiota Transplantation from High Caloric-Fed Donors Alters Glucose Metabolism in Recipient Mice, Independently of Adiposity or Exercise Status. Am J Physiol Endocrinol Metab. 2020;319(1):E203-16. https://doi.org/10.1152/ajpendo.00037.2020.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank all study participants as well as all investigators of the studies that were used throughout the course of this investigation.

Funding

EG holds a masters research award from the Fonds de recherche du Québec: Santé. (FRQS). BJA and ST hold junior scholar awards from the FRQS. Part of this study was supported by the European Union through the European Regional Development fund. The work of Estonian Genome Center, Univ. of Tartu has been supported by the European Regional Development Fund and grants No. GENTRANSMED (2014-2020.4.01.15-0012), MOBERA5 (Norface Network project no 462.16.107) and 2014-2020.4.01.16-0125. This study was also funded by the European Union through Horizon 2020 research and innovation programme under grant no 810645 and through the European Regional Development Fund project no. MOBEC008 and Estonian Research Council Grant PUT1660.

Author information

Authors and Affiliations

Authors

Contributions

Conception and design EG, PMit, BJA. Data acquisition and analysis ÉG, HM, EA, NT, TE, NG, ST, PMat. Drafting of the work EG, PMit, BJA. Revision for critical important content ÉG, HM, EA, NT, TE, NG, ST, PMat. All authors approved the final version of the manuscript.

Corresponding author

Correspondence to Benoit J. Arsenault.

Ethics declarations

Ethics approval and consent to participate

All data used in this study are in the public domain. All participants provided informed consent and study protocols were approved by their respective local ethical committees 29–32,64,67,86–94,94–97. This project was approved by the Institutional Review Board of the Quebec Heart and Lung Institute.

Competing interests

BJA is a consultant for Novartis and Silence Therapeutics and has received research contracts from Pfizer, Ionis Pharmaceuticals and Silence Therapeutics.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Description of the datasets used. Table S2. Description of the metabolites selection rational. Table S3. Summary of the instrument selection criteria for each data source. LD R2, pvalue treshold and Fstatistics. Table S4. Harmonised dataset for each exposure outcome. Table S5. Primary MR results and associated statistics. Table S6. Robust MR results and other sensitivity analyses results. Table S7. Multivariable MR results and associated statistics.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gagnon, E., Mitchell, P.L., Manikpurage, H.D. et al. Impact of the gut microbiota and associated metabolites on cardiometabolic traits, chronic diseases and human longevity: a Mendelian randomization study. J Transl Med 21, 60 (2023). https://doi.org/10.1186/s12967-022-03799-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-022-03799-5