- 1Department of Cardiology, Second Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, China
- 2Department of Cardiology, The First People's Hospital of Jiashan, Jiaxing, China
- 3Cardiovascular Key Laboratory of Zhejiang Province, Hangzhou, China
Background: Observational studies have shown gut microbiomes were associated with cardiovascular diseases (CVDs), but their roles remain controversial, and these associations have not yet been established causally.
Methods: Two-sample Mendelian randomization (MR) was used to investigate whether gut microbiome had a causal effect on the risk of CVDs. To obtain comprehensive results, we performed two sets of MR analyses, one with single nucleotide polymorphisms (SNPs) that smaller than the genome-wide statistical significance threshold (5 × 10−8) as instrumental variables, and the other with SNPs that lower than the locus-wide significance level (1 × 10−5). Summary-level statistics for CVDs, including coronary artery disease (CAD), myocardial infarction, heart failure, atrial fibrillation, stroke and its subtypes were collected. The ME estimation was performed using the inverse-variance weighted and Wald ratio methods. Sensitivity analysis was performed using the weighted median, MR-Egger, leave-one-out analysis, MR pleiotropy residual sum and outlier and MR Steiger.
Results: Based on the locus-wide significance level, genetically predicted genus Oxalobacter was positively associated with the risk of CAD (odds ratio (OR) = 1.06, 95% confidence interval (CI), 1.03 – 1.10, P = 1.67 × 10−4), family Clostridiaceae_1 was negatively correlated with stroke risk (OR = 0.83,95% CI, 0.75–0.93, P = 7.76 × 10−4) and ischemic stroke risk (OR = 0.823,95% CI, 0.74–0.92, P = 4.15 × 10−4). There was no causal relationship between other genetically predicted gut microbiome components and CVDs risk. Based on the genome-wide statistical significance threshold, the results showed that the gut microbiome had no causal relationship with CVDs risk.
Conclusion: Our findings reveal that there are beneficial or adverse causal effects of gut microbiome components on CVDs risk and provide novel insights into strategies for the prevention and management of CVDs through the gut microbiome.
Key messages
Several studies have shown significant alternations in the structure and composition of gut microbiome in cardiovascular diseases (CVDs) patients. However, the causality between the gut microbiome and CVDs remains unclear. In this 2-sample Mendelian randomization analysis, we found that genetically predicted genus Oxalobacter was positively associated with coronary artery disease risk and family Clostridiaceae_1 was associated with decrease risk of stroke and ischemic stroke.
Introduction
Cardiovascular diseases (CVDs) remain the leading cause of mortality and morbidity in the world, though the treatment of CVDs has advanced (1). CVDs cause 17.7 million deaths (31% of global deaths) every year, a figure equivalent to a third of all deaths in US and a quarter of all deaths in Europe (2). The causes of CVDs are not fully understood through substantial progress in prevention and control (3). Therefore, it remains crucial to identify protective or causative factors for CVDs.
The human gut microbiome is a complex ecosystem that provides essential functions to its host. Recently, several studies have shown significant alternations in the structure and composition of gut microbiome in cardiovascular diseases (CVDs) patients. Emoto et al. (4) reported the changes in gut microbiome composition in patients with coronary artery disease (CAD), such as the reduced abundance of the Bacteroidetes and the grown abundance of the Lactobacillales. Zuo et al. (5) found that Ruminococcus, Streptococcus and Enterococcus were overgrown in patients with atrial fibrillation (AF), and Faecalibacterium, Alistipes, Oscillibacter, and Bilophila were reduced. Research has provided further evidence of a link between gut microbiome and CVDs susceptibility through direct gut microbial transplantation. Gregory et al. (6) shown that atherosclerosis susceptibility can be transmitted via transplantation of gut microbiota. However, the causal association between the gut microbiome and CVDs remains unclear, as many other factors such as age, gender and ethnicity can influence not only gut microbiome but also CVDs development, which complicating this matter. What's more, evidence from traditional epidemiological studies fails to address the confusion caused by various biases and reverse causality and is limited by small sample sizes.
In this context, Mendelian randomization (MR) provides a way to explore causality between exposures and outcomes without any potentially detrimental intervention (7). In this study, we performed a 2-sample MR study to elucidate the potential impact of genetically predicted gut microbiome on 9 CVDs: CAD, myocardial infarction (MI), AF, heart failure (HF), and stroke and its subtypes. We also performed multivariable Mendelian randomization (MVMR) to assess the potential mediating effects of blood pressure on the identified causal associations, as numbers of studies have demonstrated that blood pressure is a risk factor for CVDs (8, 9).
Materials and methods
Study design
We conducted a 2-sample MR study using data obtained from the publicly available GWAS catalog to investigate the causality between gut microbiome and CVDs (http://www.ebi.ac.uk/gwas). Ethical approval and consent to participate were given in the original publications. Figure 1 shown an overview of the study design.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971376/fcvm-09-971376-HTML-r1/image_m/fcvm-09-971376-g001.jpg)
Figure 1. Study flow diagram. Dashed lines indicate potential pleiotropic or direct causal effects between variables that may violate MR assumptions. IV, instrumental variable; CAD, coronary artery disease; MI, myocardial infarction; HF, heart failure; AF, atrial fibrillation; IS, ischemic stroke; CES, cardioembolic stroke; LAS, large-artery atherosclerotic stroke; SVS, small-vessel stroke; IVW, multiplicative random effects inverse-variance weighted; MR, Mendelian randomization; WME, weighted-median estimator; MR-PRESSO, MR pleiotropy residual sum and outlier.
Selection of genetic instrumental variables
Single nucleotide polymorphisms (SNPs) associated with the composition of human gut microbiome were selected as instrumental variables (IVs), which from a large-scale, multiethnic GWAS study involving 18,473 individuals from various countries with 122,110 loci of variation (10). To obtain more comprehensive results, this study collected two groups of SNPs (11). One group was lower than the genome-wide statistical significance threshold (5 × 10−8) and the other group was smaller than the locus-wide significance level (1 × 10−5) (11). All SNPs were required to independently (linkage disequilibrium [LD] (9), r2 ≤ 0.01) predicted human microbiome composition. Only the one with the lowest P-value was selected if there were SNPs with highly linkage disequilibrium. To prevent potential pleiotropy, we further searched these SNPs using PhenoScanner V2 (http://www.phenoscanner.medschl.cam.ac.uk/) to assess whether the IVs were potentially related to confounders or risk factors for CVDs (Supplementary Table 12) (12, 13). The IVs were excluded from the analysis once they were potentially related to confounders or risk factors for CVDs, such body mass index, past tobacco smoking, low density lipoprotein or other factors that have been reported (14–16).
Outcome data sources
Summary-level data for CAD were extracted from a large-scale meta-analysis GWAS including 122,733 cases and 424,528 controls from the CARDIoGRAMplusC4D consortium and UK Biobank (17). Summary statistics for MI were also came from a large-scale GWAS meta-analysis of 48 studies from CARDIoGRAMplusC4D consortium with 43,676 cases and 128,199 controls (18). Genetic associations with AF were derived from a large-scale GWAS that comprised 55,114 cases and 482,295 controls of European ancestry (19). Summary data for HF obtained from a GWAS meta-analysis of 26 studies from HERMES consortium, including 47,309 cases and 930,014 controls of European ancestry (20). Aggregated data for stroke were came from a large-scale meta-analysis GWAS conducted by MEGASTROKE consortium, which comprised 40,585 cases and 406,111 controls of European ancestry (21). Among these cases, 34,217 patients were ischemic stroke, which was further divided into three subtypes, including 7,193 cardioembolic stroke cases, 4,373 large-artery atherosclerotic stroke cases and 5,386 small-vessel stroke cases. There was no overlap between the exposures and outcomes GWASs population.
Instrument strength
The variance (R2) in the MR studies stands for the proportion of the variability of the exposure explained by each genetic instrument (22). Based on previously study, the R2 for the gut microbiomes was calculated as the following formula: R2 = 2 × EAF × (1−EAF) × beta2/[2 × EAF × (1−EAF) × beta2+2 × EAF × (1−EAF) × N × se2], where EAF means effect allele frequency, beta and se means the estimated effect and its standard error of SNP on certain gut microbiome, and N means the sample size (22). Furthermore, we used the following formula to calculate the F-statistics to evaluate the weak instrument bias: F = R2 × (N−2)/(1−R2), where N refers to the sample size (23).
Statistical analysis
In this 2-sample MR, we harmonized the effect of gut microbiomes and CVD datasets, which comprised comprehensive information on SNPs, especially effect allele, standard error, beta-coefficient, P-value and sample size. When a specific gut microbiome-associated SNPs were missing from the outcome datasets, proxy SNPs (r2 > 0.8) were applied. SNPs without suitable proxies were excluded from the analyses. The multiplicative random effects inverse variance-weighted (IVW) method was applied for the primary MR analysis, which meta-analyzed the SNP-specific Wald estimates with the assumption of balanced pleiotropy (24). The Wald ratio method was performed when the MR estimate contained only one single SNP.
Sensitivity analysis
To examine the existence of horizontal pleiotropy that violated the main MR assumptions, this study performed several statistical tests. Cochran Q statistic was calculated to quantify the heterogeneity in effect sizes produced from the selected genetic IVs. An MR pleiotropy residual sum and outlier (MR-PRESSO) analysis was also be applied to check and adjust for horizontal pleiotropy by removing outliers (25).In addition, an MR-Egger regression and weighted-median estimator (WME) were performed for sensitivity analyses. We estimated the deviation of the MR-Egger intercept to detect the horizontal pleiotropy, and the difference from 0 indicating potential bias in the MR estimates (26). The WME method was supplemented to generate robust and consistent estimates of the effect, even though up to 50% of the weight came from invalid IVs (27). We also applied a leave-one-out analysis to detect for any pleiotropy affected by a single SNP. Besides, the MR Steiger test was performed to evaluate the potential effect of reverse causality of CVDs on gut microbiome (28).
To consider multiple-testing correction, the significance threshold for various taxa levels was set as P = 0.05/n, n means the taxa size. For example, the significance threshold for locus-wide significance level group was as following: class P = 3.13 × 10−3= 0.05/16, family P = 1.43 × 10−3= 0.05/35, genus p = 3.82 × 10−4 = 0.05/131, order P = 2.50 × 10−3= 0.05/20, and phylum P = 5.56 × 10−3= 0.05/9. In MVMR, a cutoff of P = 0.05 was given. All these analyses were implemented using the “TwoSampleMR” package in R Version 4.1.2.
Results
Genetic instruments for gut microbiomes
After removing SNPs that had LD effects and possibly related to confounders or risk factors for CVD, a total of 2,789 (P < 1 × 10−5) and 17 (P < 5 × 10−8) SNPs were selected as IVs (Supplementary Table 1). These SNPs were categorized according to five biological categories, including phylum, class, order, family and genus. For instance, a total of 2,672 SNPs is associated with CAD in locus wide significance level which then categorized into 9 phyla (115SNPs), 16 classes (216 SNPs), 20 orders (265 SNPs), 35 families (470 SNPs), 131 genera (1606 SNPs), and A total of 14 SNPs is associated with CAD in genome-wide statistical significance threshold. The main information of SNPs was collected systematically for further analysis, including effect allele, other allele, beta-coefficient, standard error, P-value and EAF.
Locus-wide significance level
Causal effects of gut microbiomes on CVDs
We found evidence that genus Oxalobacter was positively related to CAD risk (odds ratio (OR) = 1.06, 95% confidence interval (CI), 1.03 – 1.10, P = 1.67 × 10−4) but were not associated with the other 8 CVDs. We observed that family Clostridiaceae_1 was associated with deceased risk of stroke (OR = 0.83, 95% CI, 0.75–0.93, P = 7.76 × 10−4) and ischemic stroke (OR = 0.82, 95% CI, 0.74–0.92, P = 4.15 × 10−4), whereas the associations with the other CVDs were not significant (Table 1). In addition, no causal relationship was found between other genetically determined gut microbiome components and CVDs risk.
We then performed MVMR to adjust for CVD related traits to reassess the causality observed in our primary analysis and to explore the potential mediators (Supplementary Tables 10, 11). The causal association between genus Oxalobacter and CAD attenuated to null after adjusting for systolic blood pressure (SBP, OR = 1.01, 95% CI, 0.97–1.05, P = 0.617), diastolic blood pressure (DBP, OR = 1.00,95% CI, 0.96–1.04, P = 0.922), or both (OR = 1.04,95% CI, 0.99–1.09, P = 0.119). The causal relationship between family Clostridiaceae_1 and stroke also weakened to null after adjusting for SBP (OR = 1.03, 95% CI, 0.97–1.09, P = 0.319), DBP (OR = 1.02, 95% CI, 0.96–1.08, P = 0.480), but causality remained when adjusted for both SBP and DBP (OR = 1.07, 95% CI, 1.00–1.13, P = 0.045). The causality between family Clostridiaceae_1 and ischemic stroke also weakened to null after adjusting for SBP (OR = 1.01, 95% CI, 0.95–1.08, P =0.670), DBP (OR = 1.02, 95% CI, 0.95–1.09, P = 0.525), or both (OR = 1.05, 95% CI, 0.98–1.12, P = 0.205). These results suggest that gut microbiome might affect CVDs risk through blood pressure.
Sensitivity analyses
In the Cochran's Q statistic, heterogeneity was detected in several diseases. After using the random effects model to estimate IVW, these results did not change significantly (Supplementary Tables 3, 6). The MR-Egger regression results showed that there was no horizontal pleiotropy between the genus Oxalobacter and CAD (P = 0.896), family Clostridiaceae_1 and stroke (P = 0.145), family Clostridiaceae_1 and ischemic stroke (P = 0.851). The MR-PRESSO analysis showed that there are no outliers in the analysis of genus Oxalobacter, family Clostridiaceae_1. The results obtained by WME method were consistent with those achieved with the IVW method (Table 1). There was no significant change in the risk estimations for genetically predicted in leave-one-out analysis, proving that the causal relationship was not driven by specific SNPs (Supplementary Figures 1–3). There were no weak instrumental variables bias as the F statistics of the SNPs were all >10 (Table 2; Supplementary Table 1). The MR Steiger test indicated that there was no reverse causality (Supplementary Table 5).
Genome-wide statistical significance threshold
When the MR analysis on gut microbiome was performed as an entire, the results of IVW, MR Egger, WME and weighted mode showed that gut microbiome was not related with any CVD risk (Supplementary Table 7). The Cochrane Q statistics results showed no significant heterogeneity except for HF (P = 0.049) (Supplementary Table 7). In addition, the F statistics of all SNPs were >10. After using the random effects model to estimate IVW, the result for HF did not change significantly (Supplementary Table 7). There was no evidence of horizontal pleiotropy between IVs and CVDs as shown by MR-Egger regression analysis (Supplementary Table 7). The results of gut microbiome classification also indicated that no causal relationship between gut microbiome and CVDs (Supplementary Table 8). Heterogeneity and horizontal pleiotropy were failed to be examined due to the limited number of included SNPs.
Discussion
In this study, we used MR to investigate the potential causal relationship between gut microbiome and 9 CVDs. Our findings indicate that genetically predicted level of genus Oxalobacter might be the risk factor for CAD, and family Clostridiaceae_1 was related to a reduced risk of stroke and ischemic stroke, probably acting through blood pressure. We failed to find evidence to support a causal association between other gut microbiome and CVDs in both locus-wide significance level and genome-wide statistical significance threshold.
A compelling finding of this study is that genus Oxalobacter might increase the risk of CAD. The type species of genus Oxalobacter is Oxalobacter formigenes, which has been widely studies in nephrolithiasis (29). Oxalobacter was thought to prevent calcium nephrolithiasis through two different mechanisms: degradation of oxalate in the gut lumen with reduction of mucosal absorption and promotion of endogenous oxalate secretion by the gut mucosa (30). However, studies on the role of Oxalobacter in CAD are limited. Emoto et al. (4). analyzed the gut microbiota composition in CAD patients for the first time but did not find a link between the genus Oxalobacter and CAD. Recently, Zheng and her colleagues (31) found that 28 genera, including Oxalobacter, were significantly increased in CAD patients, which was consistent with our findings. It is a pity that Zheng et al. did not conduct further research and discussion on this finding. Thus, it is necessary to further study the possible role of Oxalobacter in CAD development.
Although many studies have examined changes in Clostridiaceae in stroke patients, the relationship between them remains unclear. Previous studies have found an overabundance of Clostridial species in post-stroke mice, which as a part of stroke-induced shift in microbiological composition (32). On the contrast, Lee et al. (33) found that Clostridiaceae were enrich in microbiota of young stroke mice compared to aged and associated with improved outcomes. In human, Xia et al. (34) found that stroke patients had decreased abundance of Clostridiaceae compared to health control, which was accord with our findings. They also found that decreased abundance of Clostridiaceae was strongly associated with more severe brain injury and a greater likelihood of unfavorable outcomes. It should be noted that the results of this study suggested that Clostridiaceae was related to both stroke and ischemic stroke, but ischemic stroke account for 85% of all stroke cases. This suggests the causality of Clostridiaceae and stroke may be due to its causal relationship with ischemic stroke to some extent.
The mechanisms of the neuroprotective effects of Clostridiaceae on stroke are not fully understood. One possibility is the protective effect of Clostridiaceae-driven short-chain fatty acids (SCFAs) on ischemic stroke. Studies have demonstrated that higher risk of stroke was correlated with lessened levels of butyrate-producing bacteria in the Gut microbiome (35). Sun et al. (36) observed that Clostridium butyricum, the type species of Clostridiaceae, has a neuroprotective effect against cerebral ischemia/reperfusion injury mice, and this neuroprotective effect may be involved to its ability to reverse the decrease of butyric acid content in the brain. Furthermore, researchers have shown that poor stroke outcome in older mice can be reversed by poststroke bacteriotherapy supplemented with 4 SCFA-producers, including Clostridium symbiosum (33).
In the initial studies, trimethylamine N-oxide (TMAO), a metabolite formed after diet, has attracted extensive attention as a potential causal link between gut microbiome and CVDs. Wang et al. (37) first identified TMAO as a gut microbiota-derived factor that has been shown to predict risk for CVDs. Studies have observed that high TMAO levels were positively associated with the risk of major cardiovascular events (38).A meta-analysis of date from 19 studies shown that participants with higher levels of TMAO had a 62% increased risk for the development of cardiovascular events (39). However, researchers failed to demonstrate a significant association of genetically predicted higher levels of TMAO and its predecessor with cardiometabolic disease through MR (40).
Our study has several advantages. The main advantage is the MR design, which enables us to estimate the causality between gut microbiome and 9 CVDs without disturbance from residual confounding or reverse causal relationship. We strictly screen the related SNPs using Plink clumping and Phenoscanner before the MR analysis. As far as we know, this is the first MR analysis on this topic.
There are also several limitations in this study. First, the available data we used were not individual-level statistics, so it might lead to inevitable biases in our results. Second, due to the lack of demographic data in the original research, we were unable to conduct the subgroup analysis, such as gender-specific causal association between gut microbiome and the risk of CVDs. Third, gut microbiome gene regulation can be greatly influenced by epigenetic phenomena (e.g., methylation) and developmental compensation mechanism (41), which may also influence the association between gut microbiome and CVDs, but such effects cannot be assessed because these are inherent defects of MR.
In conclusion, our MR study supports that there are beneficial or adverse causal effects of gut microbiome components on CVDs risk. We find suggestive evidence that genus Oxalobacter are causally association with higher risk of CAD and family Clostridiaceae_1 are causally related to lower risk of stroke and ischemic stroke. Our findings provide novel insights into strategies for the prevention and management of CVDs through the gut microbiome.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary materials, further inquiries can be directed to the corresponding author.
Ethics statement
Ethical approval was not provided for this study on human participants because ethical approval and consent to participate were carried out in the original publications. Written informed consent was not provided because ethical approval and consent to participate were carried out in the original publications.
Author contributions
YZ, DC, JL, and JJ contributed to data analysis. YZ, XZ, DC, QG, and JF structured the manuscript giving contributed to table, figures, and text editing. JJ and YZ revisited the article implementing the final manuscript form. All authors contributed to the article and approved the submitted version.
Funding
JJ was supported by grant from the National Natural Science Foundation of China (No. 82170332), Primary Research and Development Plan of Zhejiang Province (No. 2020C03016), and Chinese Cardiovascular Association-V.G fund (No. 2017-CCA-VG-006). XZ was supported by grant from the National Natural Science Foundation for young scientists of China (No. 82100346).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2022.971376/full#supplementary-material
References
1. Diseases GBD, and Injuries C, Global Global burden of 369 diseases and injuries in 204 countries and territories 1990-2019: 1990-2019: a systematic analysis for the global burden of disease study. (2019). Lancet. (2020) 396:1204–22. doi: 10.1016/S0140-6736(20)30925-9
2. Benjamin EJ, Virani SS, Callaway CW, Chamberlain AM, Chang AR, Cheng S, et al. Heart disease and stroke statistics-2018 update: a report from the American heart association. Circulation. (2018) 137:e67–492. doi: 10.1161/CIR.0000000000000573
3. Ezzati M, Obermeyer Z, Tzoulaki I, Mayosi BM, Elliott P, Leon DA, Contributions Contributions of risk factors and medical care to cardiovascular mortality trends. Nat Rev Cardiol. (2015) 12:508–30. doi: 10.1038/nrcardio.2015.82
4. Emoto T, Yamashita T, Sasaki N, Hirota Y, Hayashi T, So A, et al. Analysis of gut microbiota in coronary artery disease patients: a possible link between gut microbiota and coronary artery disease. J Atheroscler Thromb. (2016) 23:908–21. doi: 10.5551/jat.32672
5. Zuo K, Li J, Li K, Hu C, Gao Y, Chen M, et al. Disordered gut microbiota and alterations in metabolic patterns are associated with atrial fibrillation. Gigascience. (2019) 8:58. doi: 10.1093/gigascience/giz058
6. Gregory JC, Buffa JA, Org E, Wang Z, Levison BS, Zhu W, et al. Transmission of atherosclerosis susceptibility with gut microbial transplantation. J Biol Chem. (2015) 290:5647–60. doi: 10.1074/jbc.M114.618249
7. Lawlor DA, Harbord RM, Sterne JAC, Timpson N, Smith GD. Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology. Stat Med. (2008) 27:1133–63. doi: 10.1002/sim.3034
8. Fuchs FD, Whelton PK. High blood pressure and cardiovascular disease. Hypertension. (2020) 75:285–92. doi: 10.1161/HYPERTENSIONAHA.119.14240
9. Wan EYF, Fung WT, Schooling CM, Au Yeung SL, Kwok MK, Yu EYT, et al. Blood pressure and risk of cardiovascular disease in UK biobank: a mendelian randomization study. Hypertension. (2021) 77:367–75. doi: 10.1161/HYPERTENSIONAHA.120.16138
10. Kurilshikov A, Medina-Gomez C, Bacigalupe R, Radjabzadeh D, Wang J, Demirkan A, et al. Large-scale association analyses identify host factors influencing human gut microbiome composition. Nat Genet. (2021) 53:156–65. doi: 10.1038/s41588-020-00763-1
11. Sanna S, van Zuydam NR, Mahajan A, Kurilshikov A, Vich Vila A, Vosa U, et al., Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases. Nat Genet. (2019) 51:600–5. doi: 10.1038/s41588-019-0350-x
12. Staley JR, Blackshaw J, Kamat MA, Ellis S, Surendran P, Sun BB, et al. PhenoScanner: a database of human genotype-phenotype associations. Bioinformatics. (2016) 32:3207–9. doi: 10.1093/bioinformatics/btw373
13. Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, et al. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. (2019) 35:4851–53. doi: 10.1093/bioinformatics/btz469
14. Kim MS, Kim WJ, Khera AV, Kim JY, Yon DK, Lee SW, et al. Association between adiposity and cardiovascular outcomes: an umbrella review and meta-analysis of observational and Mendelian randomization studies. Eur Heart J. (2021) 42:3388–403. doi: 10.1093/eurheartj/ehab454
15. Rosoff DB, Davey Smith G, Mehta N, Clarke TK, and Lohoff FW, Evaluating the relationship between alcohol consumption, tobacco use, and cardiovascular disease: a multivariable Mendelian randomization study. PLoS Med. (2020) 17:e1003410. doi: 10.1371/journal.pmed.1003410
16. Burgess S, Ference BA, Staley JR, Freitag DF, Mason AM, Nielsen SF, et al. Association of LPA Variants with risk of coronary disease and the implications for lipoprotein(a)-lowering therapies: a mendelian randomization analysis. JAMA Cardiol. (2018) 3:619–27. doi: 10.1001/jamacardio.2018.1470
17. van der Harst P, Verweij N. Identification of 64 novel genetic loci provides an expanded view on the genetic architecture of coronary artery disease. Circ Res. (2018) 122:433–43. doi: 10.1161/CIRCRESAHA.117.312086
18. Nikpay M, Goel A, Won HH, Hall LM, Willenborg C, Kanoni S, et al. A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet. (2015) 47:1121–30. doi: 10.1038/ng.3396
19. Roselli C, Chaffin MD, Weng LC, Aeschbacher S, Ahlberg G, Albert CM, et al. Multi-ethnic genome-wide association study for atrial fibrillation. Nat Genet. (2018) 50:1225–33. doi: 10.1038/s41588-018-0133-9
20. Shah S, Henry A, Roselli C, Lin H, Sveinbjornsson G, Fatemifar G, et al. Genome-wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure. Nat Commun. (2020) 11:163. doi: 10.1038/s41467-019-13690-5
21. Malik R, Chauhan G, Traylor M, Sargurupremraj M, Okada Y, Mishra A, et al. Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. Nat Gen. (2018) 50:524–37. doi: 10.1038/s41588-018-0058-3
22. Shim H, Chasman DI, Smith JD, Mora S, Ridker PM, Nickerson DA, et al. A multivariate genome-wide association analysis of 10 LDL subfractions, and their response to statin treatment, in 1,868 Caucasians. PLoS ONE. (2015) 10:e0120758. doi: 10.1371/journal.pone.0120758
23. Burgess S, Thompson SG, Collaboration CCG, Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. (2011) 40:755–64. doi: 10.1093/ije/dyr036
24. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J, A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. (2017) 36:1783–802. doi: 10.1002/sim.7221
25. Verbanck M, Chen CY, Neale B, and Do R, Detection Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
26. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. (2015) 44:512–25. doi: 10.1093/ije/dyv080
27. Bowden J, Davey Smith G, Haycock PC, and Burgess S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965
28. Hemani G, Tilling K, and Davey Smith G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. (2017) 13:e1007081. doi: 10.1371/journal.pgen.1007081
29. Dawson KA, Allison MJ, Hartman PA. Isolation and some characteristics of anaerobic oxalate-degrading bacteria from the Rumen. Appl Environ Microb. (1980) 40:833–9. doi: 10.1128/aem.40.4.833-839.1980
30. Arvans D, Jung YC, Antonopoulos D, Koval J, Granja I, Bashir M, et al., Oxalobacter formigenes-derived bioactive factors stimulate oxalate transport by intestinal epithelial cells. J Am Soc Nephrol. (2017) 28:876–87. doi: 10.1681/ASN.2016020132
31. Zheng YY, Wu TT, Liu ZQ, Li A, Guo QQ, Ma YY, et al. Gut microbiome-based diagnostic model to predict coronary artery disease. J Agr Food Chem. (2020) 68:3548–57. doi: 10.1021/acs.jafc.0c00225
32. Stanley D, Moore RJ, Wong CHY, An insight into intestinal mucosal microbiota disruption after stroke. Sci Rep-Uk. (2018) 8:8904. doi: 10.1038/s41598-017-18904-8
33. Lee J, d'Aigle J, Atadja L, Quaicoe V, Honarpisheh P, Ganesh BP, et al. Gut microbiota-derived short-chain fatty acids promote poststroke recovery in aged mice. Circ Res. (2020) 127:453–65. doi: 10.1161/CIRCRESAHA.119.316448
34. Xia GH, You C, Gao XX, Zeng XL, Zhu JJ, Xu KY, et al. Stroke dysbiosis index. (SDI) in gut microbiome are associated with brain injury and prognosis of stroke. Front Neurol. (2019) 10:397. doi: 10.3389/fneur.2019.00397
35. Zeng X, Gao X, Peng Y, Wu Q, Zhu J, Tan C, et al. Higher risk of stroke is correlated with increased opportunistic pathogen load and reduced levels of butyrate-producing bacteria in the gut. Front Cell Infect Microbiol. (2019) 9:4. doi: 10.3389/fcimb.2019.00004
36. Sun J, Ling Z, Wang F, Chen W, Li H, Jin J, et al. Clostridium butyricum pretreatment attenuates cerebral ischemia/reperfusion injury in mice via anti-oxidation and anti-apoptosis. Neurosci Lett. (2016) 613:30–5. doi: 10.1016/j.neulet.2015.12.047
37. Wang Z, Klipfell E, Bennett BJ, Koeth R, Levison BS, Dugar B, et al. Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature. (2011) 472:57–63. doi: 10.1038/nature09922
38. Koeth RA, Wang Z, Levison BS, Buffa JA, Org E, Sheehy BT, et al. Intestinal microbiota metabolism of L-carnitine, a nutrient in red meat, promotes atherosclerosis. Nat Med. (2013) 19:576–85. doi: 10.1038/nm.3145
39. Heianza Y, Ma W, Manson JE, Rexrode KM, 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:4947. doi: 10.1161/JAHA.116.004947
40. Jia J, Dou P, Gao M, Kong X, Li C, Liu Z, et al. Assessment of causal direction between gut microbiota-dependent metabolites and cardiometabolic health: a bidirectional mendelian randomization analysis. Diabetes. (2019) 68:1747–55. doi: 10.2337/db19-0153
Keywords: Mendelian randomization, gut microbiome, cardiovascular disease, causality, Oxalobacter, Clostridiaceae
Citation: Zhang Y, Zhang X, Chen D, Lu J, Gong Q, Fang J and Jiang J (2022) Causal associations between gut microbiome and cardiovascular disease: A Mendelian randomization study. Front. Cardiovasc. Med. 9:971376. doi: 10.3389/fcvm.2022.971376
Received: 30 June 2022; Accepted: 15 August 2022;
Published: 30 August 2022.
Edited by:
Daniele Pastori, Sapienza University of Rome, ItalyReviewed by:
Hanbin Cui, Ningbo First Hospital, ChinaKazuyuki Kasahara, Nanyang Technological University, Singapore
Copyright © 2022 Zhang, Zhang, Chen, Lu, Gong, Fang and Jiang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jun Jiang, jiang-jun@zju.edu.cn
†ORCID: Jun Jiang orcid.org/0000-0001-6926-9516