Skip to main content
  • More from ADA
    • Diabetes Care
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care in Diabetes
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care
  • Subscribe
  • Log in
  • My Cart
  • Follow ada on Twitter
  • RSS
  • Visit ada on Facebook
Diabetes

Advanced Search

Main menu

  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • ADA Scientific Sessions Abstracts
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • ADA Scientific Sessions Abstracts
    • Diabetes COVID-19 Article Collection
    • Diabetes Symposium 2020
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Submit Cover Art
    • ADA Journal Policies
    • Instructions for Authors
    • ADA Peer Review
  • More from ADA
    • Diabetes Care
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care in Diabetes
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care

User menu

  • Subscribe
  • Log in
  • My Cart

Search

  • Advanced search
Diabetes
  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • ADA Scientific Sessions Abstracts
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • ADA Scientific Sessions Abstracts
    • Diabetes COVID-19 Article Collection
    • Diabetes Symposium 2020
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Submit Cover Art
    • ADA Journal Policies
    • Instructions for Authors
    • ADA Peer Review
Metabolism

Assessment of Causal Direction Between Gut Microbiota–Dependent Metabolites and Cardiometabolic Health: A Bidirectional Mendelian Randomization Analysis

  1. Jinzhu Jia1,2,
  2. Pan Dou3,
  3. Meng Gao1,
  4. Xuejun Kong4,
  5. Changwei Li5,
  6. Zhonghua Liu6 and
  7. Tao Huang1,7,8⇑
  1. 1Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China
  2. 2Department of Biostatistics, School of Public Health, Peking University, Beijing, China
  3. 3Department of Clinical Nutrition, Peking University First Hospital, Beijing, China
  4. 4Martinos Center, Massachusetts General Hospital, Charleston, MA
  5. 5Department of Epidemiology and Biostatistics, University of Georgia, Athens, GA
  6. 6Department of Statistics and Actuarial Science, The University of Hong Kong, Hong Kong, China
  7. 7Department of Global Health, School of Public Health, Peking University, Beijing, China
  8. 8Key Laboratory of Molecular Cardiovascular Sciences (Peking University), Ministry of Education, Beijing, China
  1. Corresponding author: Tao Huang, huangtaotao{at}pku.edu.cn, or Zhonghua Liu, zhhliu{at}hku.hk
Diabetes 2019 Sep; 68(9): 1747-1755. https://doi.org/10.2337/db19-0153
PreviousNext
  • Article
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF
Loading

Abstract

We examined the causal direction between gut microbiota–dependent metabolite trimethylamine N-oxide (TMAO) or its predecessors and cardiometabolic diseases, such as risk of type 2 diabetes mellitus (T2DM), coronary artery disease (CAD), myocardial infarction (MI), stroke, atrial fibrillation (AF), and chronic kidney disease (CKD). We used genetic variants as instruments to test the causal associations. Genetically predicted higher TMAO and carnitine were not associated with higher odds of T2DM, AF, CAD, MI, stroke, and CKD after Bonferroni correction (P ≤ 0.0005). However, we observed that genetically increased choline showed a suggestive association with higher risk of T2DM (odds ratio 1.84 [95% CI 1.00–3.42] per 10 units, P = 0.05). In contrast, genetically predicted higher betaine (0.68 [0.48–0.95] per 10 units, P = 0.023) was suggestively associated with a lower risk of T2DM. We observed a suggestive association of genetically increased choline with a lower level of body fat percentage (β ± SE −0.28 ± 0.11, P = 0.013) but a higher estimated glomerular filtration rate (0.10 ± 0.05, P = 0.034). We further found that T2DM (0.130 ± 0.036, P < 0.0001) and CKD (0.483 ± 0.168, P = 0.004) were causally associated with higher TMAO levels. Our Mendelian randomization findings support that T2DM and kidney disease increase TMAO levels and that observational evidence for cardiovascular diseases may be due to confounding or reverse causality.

Introduction

Gut microbiota has been recently implicated as a novel endocrine organ that plays an important role in the development of cardiovascular disease (CVD) and type 2 diabetes mellitus (T2DM) through modulating blood levels of bioactive metabolites (1,2). There is a growing appreciation that a gut microbiota–dependent metabolite, trimethylamine N-oxide (TMAO), formed from ingested choline, betaine, and carnitine in humans (3,4), is predictive of cardiovascular events (1,5). Furthermore, a greater increase in TMAO, choline, and L-carnitine has been implicated in lesser improvements in insulin sensitivity (6), adiposity, and energy metabolism (7). Animal studies have shown that the TMAO meta-organismal pathway enhances the accumulation of cholesterol in macrophages and platelet hyperreactivity (1) and exacerbates impaired glucose tolerance by blocking the hepatic insulin signaling pathway (8), which promotes CVD and diabetes in animal models.

However, the results from human studies are not always consistent (9,10). Although a recent meta-analysis of 19 observational cohorts found that elevated levels of TMAO and its precursors are associated with an increased risk of major adverse cardiovascular events (11), available evidence from observational studies mainly relies on self-reported information and is susceptible to confounding or reverse causation bias. Therefore, the causality of these observations remains unclear. Documentation of such causality can inform potential lifestyle or drug treatment targets for prevention of cardiometabolic disease.

The Mendelian randomization (MR) approach, which uses genetic variants as an instrumental variable in epidemiological study, has been widely accepted to explore the potential causal effect of exposure on disease (12–16). Such an MR method is analogous to a randomized controlled trial, where genetic alleles are randomly assorted during conception, and is less likely to be affected by confounding or reverse causation (13,17). Therefore, we performed a bidirectional MR analysis to examine the causal direction between the gut microbiota–dependent metabolite TMAO or its predecessors (e.g., choline, betaine, and carnitine) and cardiometabolic diseases (e.g., T2DM, coronary artery disease [CAD], myocardial infarction [MI], stroke, atrial fibrillation [AF], and chronic kidney disease [CKD]) and related traits using summary data from genome-wide association studies (GWAS).

Research Design and Methods

Study Design

The MR approach must satisfy three assumptions (18) (Fig. 1): 1) the genetic variant selected as the instrumental variable is associated with the gut microbiota–dependent metabolite TMAO and its predecessors; 2) the genetic variant is not associated with any unmeasured confounders of the gut microbiota–dependent metabolites and the cardiometabolic relationship; and 3) the genetic variant is associated with cardiometabolic events only through gut microbiota–dependent metabolites, not through other pathways.

Figure 1
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1

Schematic representation of an MR analysis. MR can be used to test the hypothesis that exposure (gut microbiota–related metabolites) causes T2DM, CAD, MI, AF, stroke, and CKD. The three assumptions of MR are as follows: 1) Genetic variants must be associated with gut microbiota–related metabolites; 2) genetic variants must not be associated with confounders; and 3) genetic variants must influence disease outcomes only through gut microbiota–related metabolites, not through other pathways.

Gut Microbiota–Dependent Metabolites and Instrumental Variable Selection

The gut microbiota–dependent metabolites include TMAO, choline, betaine, and carnitine. We searched PubMed for GWAS of the gut microbiota–dependent metabolites and cardiometabolic diseases and identified genetic variants for each metabolite (19) and disease (20–24). For each metabolite, we selected genetic variants (single nucleotide polymorphisms [SNPs]) at thresholds for suggestive genome-wide significance (P < 5 × 10−5) from the published GWAS (19) (Supplementary Table 1). Details on GWAS where we extracted summary-level data are presented in Table 1 and Supplementary Table 1. Contributing studies received ethical approval from their respective institutional review boards. Informed consent was obtained from all participants of contributing studies.

View this table:
  • View inline
  • View popup
Table 1

Description of modifiable factors or cardiometabolic disease

Cardiometabolic Diseases and Data Sources

For disease outcomes, summary-level data were extracted from the Diabetes Genetics Replication and Meta-analysis (DIAGRAM) Consortium for T2DM (n = 149,821) (20); the Atrial Fibrillation Consortium (AFGen) for AF (21); the Coronary Artery Disease Genomewide Replication and Meta-analysis (CARDIoGRAM) plus the Coronary Artery Disease (C4D) Genetics (CARDIoGRAMplusC4D) Consortium for CAD (60,801 patients and 123,504 control subjects) and MI (43,676 patients and 128,197 control subjects), respectively (22); the National Institute of Neurological Disorders and Stroke Genetics Network (SiGN) and International Stroke Genetics Consortium (ISGC) for stroke (37,792 patients and 397,209 control subjects) (23); and the Chronic Kidney Disease Consortium (CKDGen) (n = 133,814) for CKD (24).

Patients with T2DM were diagnosed according to 1999 World Health Organization criteria of fasting plasma glucose concentration ≥7.0 mmol/L or 2-h plasma glucose concentration ≥11.1 mmol/L, by report of diabetes medication use, or on the basis of medical record review (25). Ascertainment of AF in UK Biobank included samples with one or more of the following codes: 1) noncancer illness code, self-reported (1471, 1483); 2) operation code (1524); 3) diagnoses main/secondary (ICD-10 codes I48, I48.0–4, I48.9); 4) underlying (primary/secondary) cause of death (ICD-10 codes I48, I48.0–4, I48.9); 5) diagnoses main/secondary (ICD-9 code 4273); and 6) operative procedures main/secondary (OPCS codes K57.1, K62.1–4) (21). CKD was classified in patients with an estimated glomerular filtration rate based on creatinine (eGFRcreat) <60 mL/min/1.73 m2 (24). Cases of CAD were determined with a broad definition, including MI, acute coronary syndrome, chronic stable angina, and coronary artery stenosis >50% (22).

Cardiometabolic Traits and Data Sources

We searched PubMed for GWAS of cardiometabolic traits. Summary-level data were extracted from the Genetic Investigation of Anthropometric Traits (GIANT) Consortium for adiposity indicators such as BMI (26), body fat (27), and waist-to-hip ratio–adjusted BMI (28); the Meta-Analyses of Glucose and Insulin-Related Traits Consortium (MAGIC) (25,29) for glycemic traits such as fasting glucose, log-transformed fasting insulin, 2-h glucose, log-transformed HOMA of β-cell function (HOMA-B), log-transformed HOMA of insulin resistance (HOMA-IR), log-transformed proinsulin, and hemoglobin A1c (HbA1c); the Global Lipids Genetics Consortium (GLGC) (n = 188,577) for lipids such as HDL (30), LDL (30), total cholesterol (30), and triglycerides (30); the CKDGen (n = 133,814) for kidney function indicators such as eGFRcreat, which was estimated using the four-variable MDRD equation, and eGFR based on cystatin C (eGFRcys), which was estimated as 76.7 × (serum cystatin C) − 1.19 (24); and GWAS summary data for leptin (31) and heart rate (32).

Statistical Analysis

Linkage Disequilibrium Assessment and Pleiotropy Assessment

To verify that the SNPs selected in this study met assumptions 1 and 2, we examined the genetic association with each metabolite, measured linkage disequilibrium (LD) among all the SNPs for the same metabolite (33), and selected independent genetic variants for each metabolite (13). We chose the variant with the lowest P value for association with each metabolite if genetic variants were in LD; therefore, the SNPs selected did not violate assumption 2. All SNPs for the same metabolite showed strong association (F statistic >21.10, the strength of the instrument), thus meeting assumption 1. We used MR-Egger regression to assess the presence of pleiotropic effects on cardiometabolic outcomes (18). Using the MR-Egger method, the SNP’s effect upon each metabolite is plotted against its effect upon outcomes, and an intercept distinct from the origin provides evidence for pleiotropic effects.

MR Analysis

The estimates of the causal effect of each metabolite on outcomes were analyzed using inverse-variance weighting (IVW), which provides a combined estimate of the causal estimate from each SNP. IVW is equivalent to a two-stage least squares or allele score analysis using individual-level data and is hence considered here as conventional MR (34). Results are presented as odds ratios (ORs) (95% CIs) per 10 units of genetically predicted increase in each metabolite. In addition, results for causal effects of diseases on metabolites are presented as β ± SE per each unit higher in log odds of disease. Furthermore, complementary approaches such as simple median method, weighted median method, the mode-based estimate (MBE), and MR-Egger for multiple genetic variants were used to examine causal effect (34–37). Detailed information on these MR methods has been described previously (34–37). Finally, the MR-Egger regression test was used to evaluate the pleiotropic effects (18). Using the MR-Egger regression method, the effect of an instrumental variable on the exposure is plotted against its effect on the outcome, and an intercept distinct from the origin provides evidence for pleiotropic effects. The slope of the MR-Egger regression can provide pleiotropy-corrected causal estimates (18). In the present MR analyses, median weighted, MBE, and MR-Egger methods were considered as sensitivity analyses for MR investigations with multiple genetic variants (34–36). Power calculations for MR were conducted on the basis of mRnd (http://cnsgenomics.com/shiny/mRnd).

Analyses were performed using R 3.2.3 software (R Project for Statistical Computing). The threshold of statistical significance was P ≤ 0.0005 (0.05/100) after Bonferroni correction. P ≤ 0.05 but above the Bonferroni-corrected significance threshold was considered as suggestive of evidence for a potential association.

Results

Characteristics for Selected SNPs

The characteristics of the selected SNPs for each gut microbiota–dependent metabolite are presented in Table 1 and Supplementary Table 1. To examine assumptions 2 and 3, we tested whether any of the selected SNPs were influenced by LD. We chose the variant with the lowest P value for association with each metabolite if genetic variants were in LD. Our genetic analysis showed that 15.4% of betaine, 17.1% of carnitine, 8.0% of choline, and 9.6% of TMAO were explained by their SNPs (Supplementary Table 2). MR power calculation showed that we had 87%, 84%, 78%, and 81% power to test significant (P < 0.05) causal effects (OR 1.3) of betaine, carnitine, choline, and TMAO, respectively, on cardiometabolic events.

Gut Microbiota–Dependent Metabolites and Cardiometabolic Disease

Genetically predicted higher TMAO was not associated with a higher odds of T2DM (OR 0.96 [95% CI 0.59–1.57], P = 0.863), AF (0.99 [0.79–1.26], P = 0.961), CAD (1.00 [0.70–1.43], P = 0.986), MI (1.08 [0.73–1.60], P = 0.708), stroke (0.94 [0.53–1.65], P = 0.83), and CKD (1.08 [0.59–1.99], P = 0.794) per 10 units after Bonferroni correction (Fig. 2). Likewise, genetically predicted higher carnitine was also not associated with each cardiometabolic disease. However, we observed that genetically increased choline was suggestively associated with higher risk of T2DM (1.84 [1.00–3.42] per 10 units, P = 0.05) but not with cardiovascular events and CKD. In contrast, genetically predicted higher betaine (0.68 [0.48–0.95] per 10 units, P = 0.023) was suggestively associated with a lower risk of T2DM (Fig. 2).

Figure 2
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2

MR of gut microbiota–related metabolites and risk of T2DM, CAD, MI, AF, stroke, and CKD.

In addition, the associations were consistent in sensitivity analyses that used the simple median and weighted median methods (Supplementary Table 3). The MR-Egger regression test was used to examine the presence of a pleiotropic effect. The intercepts (SE) from MR-Egger regression testing of each outcome were centered at the origin with a CI including the null, suggesting that the results were not influenced by pleiotropy (Supplementary Table 3).

Gut Microbiota–Dependent Metabolites and Cardiometabolic Traits

Genetically predicted higher TMAO and betaine were not associated with any cardiometabolic factors after Bonferroni correction (Table 2). However, we observed a suggestive association (P < 0.05) of genetically increased choline with a lower level of body fat percentage (β ± SE −0.28 ± 0.11, P = 0.013) but a higher eGFRcys (0.10 ± 0.05, P = 0.034). Additionally, genetically predicted higher carnitine was suggestively associated with a lower level of triglycerides (−0.17 ± 0.07, P = 0.012) (Table 2). Results were similar in sensitivity analyses that used the simple median and weighted median methods (Supplementary Table 4). The results were not influenced by directional pleiotropy (P > 0.05). We further examined the genetic association of TMAO SNPs with potential confounders, such as alcohol intake and dietary factors (Supplementary Table 5). We did not find that TMAO SNPs were associated with any potential confounders, thus validating assumption 2.

View this table:
  • View inline
  • View popup
Table 2

MR analysis of gut microbiota–dependent metabolite and cardiometabolic traits

Causal Effects of Cardiometabolic Diseases on Gut-Dependent Metabolites

We further examined the causal effects of cardiometabolic diseases on gut-dependent metabolites. We found that T2DM was causally associated with lower betaine (β ± SE −0.111 ± 0.035, P = 0.002) and higher TMAO levels (0.13 ± 0.036, P < 0.0001) per 1-unit higher log odds. CKD was also causally associated with higher TMAO levels (0.483 ± 0.168, P = 0.004) per 1-unit higher log odds (Table 3 and Supplementary Table 6).

View this table:
  • View inline
  • View popup
Table 3

MR analysis of cardiometabolic diseases on gut microbiota–dependent metabolites

Discussion

In the present MR study, we did not observe a significant association of genetically predicted higher gut microbiota–dependent metabolite TMAO and its predecessor with cardiometabolic disease and related traits. However, we observed a suggestive association of genetically increased choline with a higher risk of T2DM. Furthermore, genetically increased choline showed a suggestive association with a lower body fat percentage but a higher level of eGFR, whereas genetically predicted higher carnitine was suggestively associated with a lower level of triglycerides. We further found that T2DM and CKD were causally associated with higher TMAO levels. Our findings support that T2DM and CKD increase TMAO levels and that the previously observed association of TMAO and CVD may be due to confounding or reverse causality.

Gut microbiota alteration has been implicated in the development of diabetes (38). Previous observational findings suggested that higher plasma TMAO is associated with an increased odds of newly diagnosed T2DM (39). Animal studies showed that TMAO exacerbates impaired glucose tolerance and hyperglycemia through causing inflammation or blocking the hepatic insulin signaling pathway in adipose tissue (8). Likewise, a weight loss diet intervention study demonstrated that a greater increase in TMAO is related to lesser improvements in the glycemic traits (6). However, previous studies have yielded inconsistent associations between plasma TMAO and T2DM (39–41). Our MR analysis also did not confirm the causal effect of genetically increased TMAO on risk of T2DM. Interestingly, we found that T2DM increases TMAO levels. In addition, a recent genome-wide scan suggested that human gut microbiota is influenced by a host gene that has been associated with diabetes (42), suggesting that a TMAO-T2DM association may share a genetic basis rather than a causal relationship. We cannot exclude the possibility that the suggestive genome-wide SNP (P < 5.0 × 10−5) used in the MR analysis might introduce a weak instrument bias (19), thus diluting the association. Future investigation for identifying genome-wide significant SNPs for TMAO and its causal role in the development of T2DM is required.

Furthermore, a recent prospective study indicated that a higher intake of phosphatidylcholine, the precursor to the generation of TMAO, is independently associated with an increased risk of T2DM in three U.S. populations (43). Interestingly, our MR results showed a suggestive association of genetically increased choline with a higher risk of T2DM. Our findings were supported by a previous weight loss diet intervention study that demonstrated that greater decreases in choline and L-carnitine are significantly associated with greater improvements in fasting insulin concentrations and insulin resistance (6), underscoring the importance of changes in choline and L-carnitine in improving insulin sensitivity (6).

Our findings did not corroborate the results from a previous meta-analysis of data from 19 prospective studies, which provided quantitative pooled estimates of the associations of circulating TMAO level with the incidence of cardiovascular events (11). It was estimated that participants with high TMAO levels had a 62% increased risk for the development of cardiovascular events compared with those with low TMAO levels (11). However, our MR results showed that genetically increased TMAO is not associated with a higher risk of cardiovascular events, such as CAD, MI, stroke, and AF, suggesting that observational evidence might result from confounding and reverse causality. Although we measured LD among all selected SNPs using northern Europeans from Utah samples from the 1000 Genomes Project (33), we cannot exclude the possibility that our results might be affected by unmeasured confounders. Furthermore, the MR-Egger regression showed that our results are not influenced by pleiotropy (18). Nevertheless, it is possible that the results reflect a shared genetic basis between gut microbiota–related metabolites and CVD rather than a causal relationship.

Despite strong observational evidence for an association of predecessors of TMAO, such as carnitine, choline, or betaine, with risk of CAD, MI, stroke, and AF, our results disagree with an earlier meta-analysis of observational studies showing a positive association of L‐carnitine, choline, or betaine with risk of cardiovascular events (11). These null findings suggest that our findings do not support a causal role of predecessors of TMAO in the development of CVD, but we have to acknowledge that more variants, especially suggestive genome-wide SNPs, may lead to greater pleiotropy and dilute the association in our analysis. However, we found that genetically predicted higher carnitine was suggestively associated with a lower level of triglycerides. The potential mechanism may be that dietary carnitine is metabolized in the liver by intestinal bacteria to produce TMAO by the enzyme flavin monooxygenase 3 (3,4), which is reported to be a key integrator of hepatic cholesterol and lipid metabolism (1,5).

Interestingly, an observational study has reported that TMAO, cleared by the kidney, is elevated in patients with impaired renal function (44) and portends poorer long-term survival (45), suggesting that circulating TMAO might play a role in kidney dysfunction (9,44–46). In the present MR study, our findings did not support a causal role of TMAO in kidney dysfunction. Surprisingly, we found that genetically increased choline showed a suggestive association with a higher level of eGFR. Importantly, we demonstrate that kidney disease increases TMAO levels. Given that the gut microbiota–mediated metabolism of choline and carnitine each ultimately produces TMAO (3,4), we speculate that there is a threshold in the observed association between choline and kidney function and that the optimal level of choline may be beneficial for improving kidney function. However, chronic high dietary exposure to choline or carnitine that increases TMAO appears to be toxic and directly contributes to progressive renal fibrosis and dysfunction (45); further investigation on the mechanism is thus required.

This study has several strengths. First, we have systematically assessed for the first time in our knowledge the causal role of gut microbiota–dependent metabolites in the development of T2DM, CAD, MI, stroke, AF, and CKD. Second, the present MR analysis used summary-level data from large GWAS, which might avoid bias from reverse causation and reduce confounding. Furthermore, the large sample size from GWAS summary statistics has sufficient power for reliable and lifelong causal estimation. Third, the consistent causal estimation across five methods, such as the conventional IVW, simple median, weighted median, MBE, and MR-Egger methods, suggests robustness of our findings.

However, several limitations merit consideration. First, we assumed that the associations between continuous metabolites and outcomes are linear, but our findings might suggest that there is a threshold in observed association. Therefore, further investigation is warranted on the causal role of gut microbiota–dependent metabolites in the development of cardiometabolic disease. Second, we used the suggestive genome-wide SNPs for the metabolites; thus, we cannot exclude the possibility that our findings might have been affected by weak instrument bias, which depends on the strength of the genetic instrument through the F statistic. Third, we do not have diet or lifestyle information; therefore, we cannot test whether the instrumental variables are associated with these factors. However, we used diet/lifestyle information from UK Biobank and demonstrated that the genetic variants are not associated with confounders. Fourth, we cannot exclude the possible diet-gene or gene-environment interactions on outcomes, which may influence the observed results. Finally, completely ruling out an alternative direct causal pathway is a challenge for all MR analyses, particularly for metabolites determined by both gut microbiota and multiple genetic variants.

In summary, our findings from an MR approach support that TMAO and kidney disease increase TMAO levels and that observational evidence for gut microbiota–dependent metabolites and CVDs may be due to confounding or reverse causality. The human genome influences both gut microbiota and chronic disease; therefore, we cannot exclude the possibility that such an association may have a genetic basis rather than a causal relationship. A large-scale genome-wide scan for genetic variants of gut microbiota–related metabolites and further investigation in understanding the potential role of gut microbiota–dependent metabolites in the development of cardiometabolic disease are required.

Article Information

Acknowledgments. The authors thank Dr. Robert Gerszten (Division of Cardiovascular Medicine, Beth Israel Deaconess Medical Center, Boston, MA), Dr. Ming-Huei Chen (Population Sciences Branch, National Heart, Lung, and Blood Institute’s Framingham Heart Study, Framingham, MA), and Dr. Qiong Yang (Department of Biostatistics, Boston University School of Public Health, Boston, MA) for sharing the GWAS summary statistics of gut-dependent metabolites.

Funding. The study was supported by Peking University Start-up grant BMU2018YJ002, High-performance Computing Platform of Peking University, and Beijing Technology and Business University grant 88442Y0033.

The funding organizations had no role in the preparation of the manuscript.

Duality of Interest. No potential conflicts of interest relevant to this article were reported.

Author Contributions. J.J., P.D., M.G., X.K., C.L., Z.L., and T.H. contributed to the statistical analysis, performed the data analysis, wrote the manuscript, critically reviewed the manuscript during the writing process, and approved the final version to be published. J.J., P.D., and T.H. designed the research. J.J., Z.L., and T.H. are the guarantors of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

Data and Resource Availability. All data used in the current study were obtained from GWAS summary statistics, which were publicly released by genetic consortia.

Footnotes

  • This article contains Supplementary Data online at http://diabetes.diabetesjournals.org/lookup/suppl/doi:10.2337/db19-0153/-/DC1.

  • Received February 12, 2019.
  • Accepted May 29, 2019.
  • © 2019 by the American Diabetes Association.
http://www.diabetesjournals.org/content/license

Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. More information is available at http://www.diabetesjournals.org/content/license.

References

  1. ↵
    1. Wang Z,
    2. Klipfell E,
    3. Bennett BJ, et al
    . Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature 2011;472:57–63pmid:21475195
    OpenUrlCrossRefPubMedWeb of Science
  2. ↵
    1. Forslund K,
    2. Hildebrand F,
    3. Nielsen T, et al.; MetaHIT Consortium
    . Disentangling type 2 diabetes and metformin treatment signatures in the human gut microbiota. Nature 2015;528:262–266pmid:26633628
    OpenUrlCrossRefPubMed
  3. ↵
    1. Tang WH,
    2. Wang Z,
    3. Levison BS, et al
    . Intestinal microbial metabolism of phosphatidylcholine and cardiovascular risk. N Engl J Med 2013;368:1575–1584pmid:23614584
    OpenUrlCrossRefPubMedWeb of Science
  4. ↵
    1. Bennett BJ,
    2. de Aguiar Vallim TQ,
    3. Wang Z, et al
    . Trimethylamine-N-oxide, a metabolite associated with atherosclerosis, exhibits complex genetic and dietary regulation. Cell Metab 2013;17:49–60pmid:23312283
    OpenUrlCrossRefPubMedWeb of Science
  5. ↵
    1. LeLeiko RM,
    2. Vaccari CS,
    3. Sola S, et al
    . Usefulness of elevations in serum choline and free F2)-isoprostane to predict 30-day cardiovascular outcomes in patients with acute coronary syndrome. Am J Cardiol 2009;104:638–643pmid:19699337
    OpenUrlCrossRefPubMed
  6. ↵
    1. Heianza Y,
    2. Sun D,
    3. Li X, et al
    . Gut microbiota metabolites, amino acid metabolites and improvements in insulin sensitivity and glucose metabolism: the POUNDS Lost trial. Gut 2019;68:263–270pmid:29860242
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Heianza Y,
    2. Sun D,
    3. Smith SR,
    4. Bray GA,
    5. Sacks FM,
    6. Qi L
    . Changes in gut microbiota-related metabolites and long-term successful weight loss in response to weight-loss diets: the POUNDS Lost trial. Diabetes Care 2018;41:413–419pmid:29305401
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Gao X,
    2. Liu X,
    3. Xu J,
    4. Xue C,
    5. Xue Y,
    6. Wang Y
    . Dietary trimethylamine N-oxide exacerbates impaired glucose tolerance in mice fed a high fat diet. J Biosci Bioeng 2014;118:476–481pmid:24721123
    OpenUrlCrossRefPubMed
  9. ↵
    1. Kaysen GA,
    2. Johansen KL,
    3. Chertow GM, et al
    . Associations of trimethylamine N-oxide with nutritional and inflammatory biomarkers and cardiovascular outcomes in patients new to dialysis. J Ren Nutr 2015;25:351–356
    OpenUrlCrossRefPubMed
  10. ↵
    1. Shafi T,
    2. Powe NR,
    3. Meyer TW, et al
    . Trimethylamine N-oxide and cardiovascular events in hemodialysis patients. J Am Soc Nephrol 2017;28:321–331pmid:27436853
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Heianza Y,
    2. Ma W,
    3. Manson JE,
    4. Rexrode KM,
    5. Qi L
    . Gut microbiota metabolites and risk of major adverse cardiovascular disease events and death: a systematic review and meta-analysis of prospective studies. J Am Heart Assoc 2017;6:e004947pmid:28663251
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Geng T,
    2. Smith CE,
    3. Li C,
    4. Huang T
    . Childhood BMI and adult type 2 diabetes, coronary artery diseases, chronic kidney disease, and cardiometabolic traits: a Mendelian randomization analysis. Diabetes Care 2018;41:1089–1096pmid:29483184
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Lawlor DA,
    2. Harbord RM,
    3. Sterne JAC,
    4. Timpson N,
    5. Davey Smith G
    . Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med 2008;27:1133–1163pmid:17886233
    OpenUrlCrossRefPubMed
    1. Ding M,
    2. Huang T,
    3. Bergholdt HK,
    4. Nordestgaard BG,
    5. Ellervik C,
    6. Qi L; CHARGE Consortium
    . Dairy consumption, systolic blood pressure, and risk of hypertension: Mendelian randomization study. BMJ 2017;356:j1000pmid:28302601
    OpenUrlAbstract/FREE Full Text
    1. Huang T,
    2. Ren J,
    3. Huang J,
    4. Li D
    . Association of homocysteine with type 2 diabetes: a meta-analysis implementing Mendelian randomization approach. BMC Genomics 2013;14:867pmid:24320691
    OpenUrlCrossRefPubMed
  14. ↵
    1. Benn M,
    2. Nordestgaard BG
    . From genome-wide association studies to Mendelian randomization: novel opportunities for understanding cardiovascular disease causality, pathogenesis, prevention, and treatment. Cardiovasc Res 2018;114:1192–1208pmid:29471399
    OpenUrlCrossRefPubMed
  15. ↵
    1. Didelez V,
    2. Sheehan N
    . Mendelian randomization as an instrumental variable approach to causal inference. Stat Methods Med Res 2007;16:309–330pmid:17715159
    OpenUrlCrossRefPubMedWeb of Science
  16. ↵
    1. Bowden J,
    2. Davey Smith G,
    3. Burgess S
    . Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol 2015;44:512–525pmid:26050253
    OpenUrlCrossRefPubMed
  17. ↵
    1. Rhee EP,
    2. Ho JE,
    3. Chen MH, et al
    . A genome-wide association study of the human metabolome in a community-based cohort. Cell Metab 2013;18:130–143pmid:23823483
    OpenUrlCrossRefPubMedWeb of Science
  18. ↵
    1. Mahajan A,
    2. Go MJ,
    3. Zhang W, et al.; DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium; Asian Genetic Epidemiology Network Type 2 Diabetes (AGEN-T2D) Consortium; South Asian Type 2 Diabetes (SAT2D) Consortium; Mexican American Type 2 Diabetes (MAT2D) Consortium; Type 2 Diabetes Genetic Exploration by Next-generation sequencing in multi-Ethnic Samples (T2D-GENES) Consortium
    . Genome-wide trans-ancestry meta-analysis provides insight into the genetic architecture of type 2 diabetes susceptibility. Nat Genet 2014;46:234–244pmid:24509480
    OpenUrlCrossRefPubMed
  19. ↵
    1. Roselli C,
    2. Chaffin MD,
    3. Weng LC, et al
    . Multi-ethnic genome-wide association study for atrial fibrillation. Nat Genet 2018;50:1225–1233pmid:29892015
    OpenUrlPubMed
  20. ↵
    1. Nikpay M,
    2. Goel A,
    3. Won HH, et al
    . A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet 2015;47:1121–1130pmid:26343387
    OpenUrlCrossRefPubMed
  21. ↵
    1. NINDS Stroke Genetics Network (SiGN)
    2. International Stroke Genetics Consortium (ISGC)
    . Loci associated with ischaemic stroke and its subtypes (SiGN): a genome-wide association study [published correction appears in Lancet Neurol 2016;15:241]. Lancet Neurol 2016;15:174–184pmid:26708676
    OpenUrlCrossRefPubMed
  22. ↵
    1. Pattaro C,
    2. Teumer A,
    3. Gorski M, et al.; ICBP Consortium; AGEN Consortium; CARDIOGRAM; CHARGe-Heart Failure Group; ECHOGen Consortium
    . Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function. Nat Commun 2016;7:10023pmid:26831199
    OpenUrlCrossRefPubMed
  23. ↵
    1. Scott RA,
    2. Lagou V,
    3. Welch RP, et al.; DIAbetes Genetics Replication and Meta-analysis (DIAGRAM) Consortium
    . Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways. Nat Genet 2012;44:991–1005pmid:22885924
    OpenUrlCrossRefPubMed
  24. ↵
    1. Locke AE,
    2. Kahali B,
    3. Berndt SI, et al.; LifeLines Cohort Study; ADIPOGen Consortium; AGEN-BMI Working Group; CARDIOGRAMplusC4D Consortium; CKDGen Consortium; GLGC; ICBP; MAGIC Investigators; MuTHER Consortium; MIGen Consortium; PAGE Consortium; ReproGen Consortium; GENIE Consortium; International Endogene Consortium
    . Genetic studies of body mass index yield new insights for obesity biology. Nature 2015;518:197–206pmid:25673413
    OpenUrlCrossRefPubMed
  25. ↵
    1. Lu Y,
    2. Day FR,
    3. Gustafsson S, et al
    . New loci for body fat percentage reveal link between adiposity and cardiometabolic disease risk. Nat Commun 2016;7:10495pmid:26833246
    OpenUrlCrossRefPubMed
  26. ↵
    1. Shungin D,
    2. Winkler TW,
    3. Croteau-Chonka DC, et al.; ADIPOGen Consortium; CARDIOGRAMplusC4D Consortium; CKDGen Consortium; GEFOS Consortium; GENIE Consortium; GLGC; ICBP; International Endogene Consortium; LifeLines Cohort Study; MAGIC Investigators; MuTHER Consortium; PAGE Consortium; ReproGen Consortium
    . New genetic loci link adipose and insulin biology to body fat distribution. Nature 2015;518:187–196pmid:25673412
    OpenUrlCrossRefPubMedWeb of Science
  27. ↵
    1. Dupuis J,
    2. Langenberg C,
    3. Prokopenko I, et al.; DIAGRAM Consortium; GIANT Consortium; Global BPgen Consortium; Anders Hamsten on behalf of Procardis Consortium; MAGIC Investigators
    . New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet 2010;42:105–116pmid:20081858
    OpenUrlCrossRefPubMedWeb of Science
  28. ↵
    1. Willer CJ,
    2. Schmidt EM,
    3. Sengupta S, et al.; Global Lipids Genetics Consortium
    . Discovery and refinement of loci associated with lipid levels. Nat Genet 2013;45:1274–1283pmid:24097068
    OpenUrlCrossRefPubMed
  29. ↵
    1. Kilpeläinen TO,
    2. Carli JFM,
    3. Skowronski AA, et al
    . Genome-wide meta-analysis uncovers novel loci influencing circulating leptin levels. Nat Commun 2016;7:10494pmid:26833098
    OpenUrlCrossRefPubMed
  30. ↵
    1. Eppinga RN,
    2. Hagemeijer Y,
    3. Burgess S, et al
    . Identification of genomic loci associated with resting heart rate and shared genetic predictors with all-cause mortality. Nat Genet 2016;48:1557–1563pmid:27798624
    OpenUrlCrossRefPubMed
  31. ↵
    1. Abecasis GR,
    2. Auton A,
    3. Brooks LD, et al.; 1000 Genomes Project Consortium
    . An integrated map of genetic variation from 1,092 human genomes. Nature 2012;491:56–65pmid:23128226
    OpenUrlCrossRefPubMedWeb of Science
  32. ↵
    1. Burgess S,
    2. Butterworth A,
    3. Thompson SG
    . Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol 2013;37:658–665pmid:24114802
    OpenUrlCrossRefPubMed
    1. Yavorska OO,
    2. Burgess S
    . MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data. Int J Epidemiol 2017;46:1734–1739pmid:28398548
    OpenUrlCrossRefPubMed
  33. ↵
    1. Bowden J,
    2. Davey Smith G,
    3. Haycock PC,
    4. Burgess S
    . Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol 2016;40:304–314pmid:27061298
    OpenUrlCrossRefPubMed
  34. ↵
    1. Hartwig FP,
    2. Davey Smith G,
    3. Bowden J
    . Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol 2017;46:1985–1998pmid:29040600
    OpenUrlCrossRefPubMed
  35. ↵
    1. Moon JY,
    2. Zolnik CP,
    3. Wang Z, et al
    . Gut microbiota and plasma metabolites associated with diabetes in women with, or at high risk for, HIV infection. EBioMedicine 2018;37:392–400pmid:30366816
    OpenUrlPubMed
  36. ↵
    1. Shan Z,
    2. Sun T,
    3. Huang H, et al
    . Association between microbiota-dependent metabolite trimethylamine-N-oxide and type 2 diabetes. Am J Clin Nutr 2017;106:888–894pmid:28724646
    OpenUrlAbstract/FREE Full Text
    1. Qiu G,
    2. Zheng Y,
    3. Wang H, et al
    . Plasma metabolomics identified novel metabolites associated with risk of type 2 diabetes in two prospective cohorts of Chinese adults. Int J Epidemiol 2016;45:1507–1516pmid:27694567
    OpenUrlCrossRefPubMed
  37. ↵
    1. Svingen GF,
    2. Schartum-Hansen H,
    3. Pedersen ER, et al
    . Prospective associations of systemic and urinary choline metabolites with incident type 2 diabetes. Clin Chem 2016;62:755–765pmid:26980210
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Wang J,
    2. Thingholm LB,
    3. Skiecevičienė J, et al
    . Genome-wide association analysis identifies variation in vitamin D receptor and other host factors influencing the gut microbiota. Nat Genet 2016;48:1396–1406pmid:27723756
    OpenUrlCrossRefPubMed
  39. ↵
    1. Li Y,
    2. Wang DD,
    3. Chiuve SE, et al
    . Dietary phosphatidylcholine intake and type 2 diabetes in men and women. Diabetes Care 2015;38:e13–e14pmid:25614692
    OpenUrlFREE Full Text
  40. ↵
    1. Bain MA,
    2. Faull R,
    3. Fornasini G,
    4. Milne RW,
    5. Evans AM
    . Accumulation of trimethylamine and trimethylamine-N-oxide in end-stage renal disease patients undergoing haemodialysis. Nephrol Dial Transplant 2006;21:1300–1304pmid:16401621
    OpenUrlCrossRefPubMedWeb of Science
  41. ↵
    1. Tang WHW,
    2. Wang Z,
    3. Kennedy DJ, et al
    . Gut microbiota-dependent trimethylamine N-oxide (TMAO) pathway contributes to both development of renal insufficiency and mortality risk in chronic kidney disease. Circ Res 2015;116:448–455pmid:25599331
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Stubbs JR,
    2. House JA,
    3. Ocque AJ, et al
    . Serum trimethylamine-N-oxide is elevated in CKD and correlates with coronary atherosclerosis burden. J Am Soc Nephrol 2016;27:305–313pmid:26229137
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top
Diabetes: 68 (9)

In this Issue

September 2019, 68(9)
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by Author
  • Masthead (PDF)
Sign up to receive current issue alerts
View Selected Citations (0)
Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word about Diabetes.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Assessment of Causal Direction Between Gut Microbiota–Dependent Metabolites and Cardiometabolic Health: A Bidirectional Mendelian Randomization Analysis
(Your Name) has forwarded a page to you from Diabetes
(Your Name) thought you would like to see this page from the Diabetes web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Assessment of Causal Direction Between Gut Microbiota–Dependent Metabolites and Cardiometabolic Health: A Bidirectional Mendelian Randomization Analysis
Jinzhu Jia, Pan Dou, Meng Gao, Xuejun Kong, Changwei Li, Zhonghua Liu, Tao Huang
Diabetes Sep 2019, 68 (9) 1747-1755; DOI: 10.2337/db19-0153

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Add to Selected Citations
Share

Assessment of Causal Direction Between Gut Microbiota–Dependent Metabolites and Cardiometabolic Health: A Bidirectional Mendelian Randomization Analysis
Jinzhu Jia, Pan Dou, Meng Gao, Xuejun Kong, Changwei Li, Zhonghua Liu, Tao Huang
Diabetes Sep 2019, 68 (9) 1747-1755; DOI: 10.2337/db19-0153
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Research Design and Methods
    • Results
    • Discussion
    • Article Information
    • Footnotes
    • References
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Ncor2/PPARα-Dependent Upregulation of MCUb in the Type 2 Diabetic Heart Impacts Cardiac Metabolic Flexibility and Function
  • Role of the Neutral Amino Acid Transporter SLC7A10 in Adipocyte Lipid Storage, Obesity, and Insulin Resistance
  • Heme Oxygenase-1 Regulates Ferrous Iron and Foxo1 in Control of Hepatic Gluconeogenesis
Show more Metabolism

Similar Articles

Navigate

  • Current Issue
  • Online Ahead of Print
  • Scientific Sessions Abstracts
  • Collections
  • Archives
  • Submit
  • Subscribe
  • Email Alerts
  • RSS Feeds

More Information

  • About the Journal
  • Instructions for Authors
  • Journal Policies
  • Reprints and Permissions
  • Advertising
  • Privacy Policy: ADA Journals
  • Copyright Notice/Public Access Policy
  • Contact Us

Other ADA Resources

  • Diabetes Care
  • Clinical Diabetes
  • Diabetes Spectrum
  • Scientific Sessions Abstracts
  • Standards of Medical Care in Diabetes
  • BMJ Open - Diabetes Research & Care
  • Professional Books
  • Diabetes Forecast

 

  • DiabetesJournals.org
  • Diabetes Core Update
  • ADA's DiabetesPro
  • ADA Member Directory
  • Diabetes.org

© 2021 by the American Diabetes Association. Diabetes Print ISSN: 0012-1797, Online ISSN: 1939-327X.