- 1Department of General Surgery, Suzhou Ninth People’s Hospital, Suzhou, China
- 2Inner Mongolia Medical University, Hohhot, China
- 3Inner Mongolia Autonomous Region People’s Hospital, Hohhot, China
- 4Suzhou Hospital of Anhui Medical University, Anhui, China
- 5Inner Mongolia Autonomous Region Health Management Service Center, Hohhot, China
Background: Recent studies have shown that the gut microbiota is closely related to the pathogenesis of Inflammatory Bowel Disease (IBD), but the causal nature is largely unknown. The purpose of this study was to assess the causal relationship between intestinal bacteria and IBD and to identify specific pathogenic bacterial taxa via the Mendelian randomization (MR) analysis.
Materials and Methods: MR analysis was performed on genome-wide association study (GWAS) summary statistics of gut microbiota and IBD. Specifically, the TwinsUK microbiota GWAS (N = 1,126 twin pairs) was used as exposure. The UK inflammatory bowel disease (UKIBD) and the Understanding Social Program (USP) study GWAS (N = 48,328) was used as discovery outcome, and the British IBD study (N = 35,289) was used as replication outcome. SNPs associated with bacteria abundance at the suggestive significance level (α = 1.0 × 10–5) were used as instrumental variables. Bacteria were grouped into families and genera.
Results: In the discovery sample, a total of 30 features were available for analysis, including 15 families and 15 genera. Three features were nominally significant, including one family (Verrucomicrobiaceae, 2 IVs, beta = −0.04, p = 0.05) and two genera (Akkermansia, 2 IVs, beta = 0.04, p = 0.05; Dorea, 2 IVs, beta = −0.07, p = 0.04). All of them were successfully replicated in the replication sample (Verrucomicrobiaceae and Akkermansia Preplication = 0.02, Dorea Preplication = 0.01) with consistent effect direction.
Conclusion: We identified specific pathogenic bacteria features that were causally associated with the risk of IBD, thus offering new insights into the prevention and diagnosis of IBD.
Introduction
Inflammatory bowel disease (IBD) is a chronic non-specific inflammatory disease that invades colonic mucosa without gender advantage (Matsuoka et al., 2018). The peak age of IBD is between 20 and 40 years old (Loftus, 2004; Bernstein et al., 2006; Cosnes et al., 2011). The main symptoms of IBD are abdominal pain, diarrhea, mucous bloody stool, as well as extra-intestinal symptoms. IBD is mostly common in developed countries including North America, Europe, Australia, and New Zealand, with incidence rate as high as 20–100 per million people. It is estimated that as many as 1 million Americans suffer from IBD (Cohen et al., 2010; Magro et al., 2012). In recent decades, the incidence of IBD has been rising all over the world, especially in East Asian (Loftus et al., 2007; Bengtson et al., 2009; Cosnes et al., 2011; Molodecky et al., 2012).
The pathogenesis of IBD has not been fully elucidated. It has a strong genetic determinant. For instance, first-degree relatives of patients with IBD are 4 to 20 times more likely to develop IBD (Kevans et al., 2016). Recent genome-wide association studies (GWASs) have identified more than 200 responsible genomic loci associated with IBD (Turpin et al., 2018). Despite these fruitful findings, its pathogenic mechanism has not been fully understood yet. On the other hand, gut microbiota may be related to the pathogenesis of IBD (Nishida et al., 2018). Imbalance of gut microbiota coupled with impaired intestinal bacterial clearance could enhance the invasiveness of pathogens, disrupt intestinal immune response, accelerate intestinal inflammation, and eventually lead to IBD. In a recent controlled trial, patients in the fecal microbiota transplantation group showed significant clinical improvement, indicating that high-dose fecal microbiota transplantation is an effective method for the treatment of active IBD (Paramsothy et al., 2017). Another study indicates that the low abundance of Phascolarctobacterium is positively correlated with the occurrence of IBD (Bajer et al., 2017).
Although previous extensive studies have established observational associations between gut microbiota and IBD developing risk, the causal nature is largely unclear. Mendelian randomization (MR) analysis is a statistical approach that aims to infer causal relationship from observational association results (Lee and Lim, 2019). With the rapidly increasing genetic data at both microbiota and complex disease sides, MR has been widely applied in recent years. MR approach has three essential assumptions: (1) Instrumental variable (IV) is strongly associated with exposure; (2) IV is not associated with any confounders of exposure; and (3) The association of IV with outcome is only through exposure. It has been used to infer the causal relationship from gut microbiota to type 2 diabetes, neurodegenerative diseases, and bone density (Burgess et al., 2013; Bowden et al., 2015; Goodrich et al., 2016; Quigley, 2017; Verbanck et al., 2018).
In the present study, in order to explore the causal relationship from gut microbiota to IBD, and to identify specific pathogenic bacteria taxa, we conducted a two-sample MR study based on GWAS summary data. In brief, summary data from the microbiota GWAS (MGWAS) of the TwinsUK study were used as exposure, and GWAS summary statistics from two IBD GWAS were used as discovery and replication outcomes.
Materials and Methods
GWAS Summary Statistics
The MR analysis was performed on GWAS summary statistics of both microbiota and IBD. All data were retrieved from previously published studies that were released to the public.
The microbiota GWAS summary statistics from the TwinsUK study (Goodrich et al., 2016) served as exposure. In brief, The TwinsUK study collected 3,261 fecal samples from 1,126 twin pairs from the TwinsUK Registry in British. Microbiota 16S rRNA was sequenced using Illumina Miseq 2 × 250 bp sequencing platform, followed by host genome genotyping using Illumina HumanHap610 Quad Chip. For genotype imputation, the 1,000 Genomes project (Phase 3) reference panel was used. Sixty-one bacteria taxa were found to be associated with 307 host SNPs with p-values ranging from 7.33 × 10–5 to 4.94 × 10–9 (Supplementary Table 1).
The discovery outcome sample UK IBD and Understanding Social Program (UKIBD and USP) is a GWAS study based on a general prospective population cohort of European ancestry with 12,924 cases and 35,391 controls. Host genome was genotyped by the HumanCyto SNP-12 BeadChip and the Immunochip arrays, and was imputed into the UK IBD Genetics Consortium and UK10K Consortium reference panel (Burgess et al., 2013). A total of 38 genomic loci were identified at the genome-wide significance level (p < 5.0 × 10–8), increasing the number of known IBD risk sites to 200.
The replication British IBD sample was the GWAS of 16,452 IBD British cases and 18,837 controls. Participants were genotyped on the Human Core Exome v12.1, the Affymetrix 500K, or the Affymetrix 6.0 genotyping array.
Instrumental Variable Selection
The same criteria were used for IV selection in both discovery and replication samples. IVs were grouped at family or genus level. Specifically, a bacterial feature was defined as a family or genus. SNPs associated with bacterial taxa in one feature were grouped together for that feature. As a QC step, palindrome SNPs, which are defined as SNPs with ambiguous strand information (e.g., A/T or G/C polymorphisms), were removed. SNP-feature association threshold was set to be 1.0 × 10–5. In order to eliminate the effect of linkage disequilibrium (LD), SNPs within each feature were clumped with PLINK (v1.9). The LD threshold was set to be r2 < 0.1, and the clustering window was set to be 500 kb. LD was estimated on the 1,000 Genome Project sequencing data (Phase 3).
In order to minimize the effect of horizontal pleiotropy. MR-PRESSO global test and outlier test were applied (Verbanck et al., 2018). The MR-PRESSO outlier test calculates the p-value for the significance of pleiotropy for each SNP, while the MR-PRESSO Global test calculates the p-value for the overall level of pleiotropy. Each individual SNP was deleted in turn and the MR-PRESSO Outlier test was applied to the set of remaining SNPs (Verbanck et al., 2018). All significant SNPs were removed. A MR-PRESSO Global test was finally performed to monitor the overall pleiotropic effect. Non-significant SNPs were used for subsequent MR analysis.
MR Analysis
Upon the selection of qualified SNPs, MR analysis was then performed for a causal relationship from microbiota feature to IBD risk. Specifically, each microbiota feature was tested for its association. For features with multiple IVs, the inverse-variance weighted (IVW) test (Burgess et al., 2013) was applied. For features with only one IV, the Wald ratio test was applied. The results of IVW were also cross-validated by three alternative tests, including the MR-Egger regression (Bowden et al., 2015), the weighted median estimator (Bowden et al., 2016) and the MR-PRESSO (Verbanck et al., 2018).
Nominally significant results identified in the discovery sample were subjected to be replicated in the replication sample, with the same analysis procedures.
The horizontal heterogeneity effect was examined by the IVW test and the MR-Egger regression. Meanwhile, a leave-one-out sensitivity analysis was performed to monitor if significant associations were dominated by a single SNP.
All the above analyses (including sensitivity analysis and MR analysis) were implemented within the R packages TwoSampleMR1 (Hemani et al., 2018) and MRPRESSO2 (Sanna et al., 2019).
Results
The flow chart of the present study is displayed in Figure 1. In the discovery sample, there are 237 host SNPs that are associated with gut microbiota features at the significance threshold p < 1.0 × 10–5. After clumping, 168 and 80 SNPs are left for 15 families and 15 genera, respectively (Supplementary Table 2). Two families with the largest number of SNPs are Lachnospiraceae (51 SNPs) and Ruminococcaceae (51 SNPs), followed by Bacteroidaceae (36 SNPs). There are five families, Barnesiellaceae, Enterobacteriaceae, Rikenellaceae, Streptococcaceae, and Veillonellaceae, each having only one SNP. At the genus level, the genus with the largest number of SNPs is Bacteroides (36 SNPs), followed by Faecalibacterium (9 SNPs) and Coprococcus (6 SNPs). There are four genera each having only one SNP, Anaerostipes, Dorea, Streptococcus, and Veillonella. Of note, genus is a child taxon of family, therefore SNPs contained in both features may overlap. For example, the genus Faecalibacterium is a child taxon of the family Ruminococcaceae. The SNPs in them are partly identical.
Figure 1. Diagrammatic description of MR analysis in the discovery and replication. (A) The whole workflow of MR analysis. (B) The main results and the change in the number of SNPs.
For features containing multiple IVs, no outliers were detected using the MR-PRESSO outlier test and no evidence of horizontal pleiotropy (both MR-PRESSO Global test p > 0.05/15 = p > 3.3 × 10–3 and MR-Egger regression p > 0.05) was observed.
MR Analysis
In the discovery sample, after removing potentially pleiotropic SNPs, one family and two genera are significant at the nominal level (p < = 0.05): family Verrucomicrobiaceae (2 IVs, beta = 0.04, p = 0.05), genus Akkermansia (2 IVs, beta = 0.04, p = 0.05) and genus Dorea (1 IV, beta = −0.07, p = 0.04).
In sum, three features (one family ++ two genera) are causally associated with IBD in the discovery sample. These three features were replicated in the British IBD replication sample. The same IVs are available in the replication sample. Using the same IVW test, the replication p-value is significant (p = 0.02) and the effect direction is consistent for family Verrucomicrobiaceae and genus Akkermansia (Table 1). For the other genus Dorea, only one SNP rs10743315 is qualified as the IV. Using the Wald ratio test, the MR p-value is 0.01, again with consistent effect direction. Moreover, there is no evidence of heterogeneity at the three identified features in both discovery and replication sample. Detailed information of the 3 IVs is listed in Table 2.
Discussion
In this study, we used MR analysis to evaluate the causal relationship between gut microbiota and IBD. Using the summary statistics of one microbiota GWAS and 2 IBD GWASs, we identified and replicated three bacterial taxa, one family Verrucomicrobiaceae and two genera Akkermansia and Dorea, that may have causal relationship with IBD. Our study confirmed that gut microbiota can aggravate IBD, suggesting that gut microbiota plays a regulatory role in IBD.
The gut microbiota is an intricate and dynamic collective of ecological microbial communities that are colonized in the human gut, even called a “forgotten organ” (O’Hara and Shanahan, 2006; Backhed et al., 2015). Gut microbiota is not only an important part of immune and metabolic health, but also regulate central nervous system and relevant disorders, including movement disorders, neurodegenerative diseases, behavioral disorders, neuroimmune-mediated diseases, and Cerebrovascular accident (Strandwitz, 2018). More than 90% of the gut microbiota that maintain intestinal health and balance in adults consist of four phylums of Firmicutes, Bacteroides, Actinobacteria, and Proteobacteria (Matsuoka and Kanai, 2015). The large intestine comprises the densest and metabolism-active microorganism in healthy individuals, which is predominated by anaerobic microbiota, two phyla Firmicutes and Bacteroidetes, apart from Actinobacteria, Proteobacteria, and Verrucomicrobia (Eckburg et al., 2005).
The Dorea identified in this study belongs to the Lachnospiraceae family, which mainly exists in the gut microbiota of mammals and humans. One previous study has established a link between Lachnospiraceae and IBD (Lee et al., 2020). Another recent studies has also confirmed that the level of Lachnospiraceae and butyric acid gets decreased in IBD patients (Sasaki et al., 2019). The genus Akkermansia is present abundantly in the human gastrointestinal tract where it is believed to be a key symbiont member of the microbiota (Collado et al., 2007; Derrien et al., 2008; van Passel et al., 2011; Clarke et al., 2014; Guo et al., 2016). Extensive studies demonstrate that the lower level of Akkermansia is found in patients with IBD and other metabolic disorders, suggesting that Akkermansia may have potential anti-inflammatory properties (Zhang et al., 2016).
Previous studies have shown that the imbalance of gut microbiota is one of the pathogenic factors of IBD, but the specific regulatory mechanism is yet poorly understood. One possible mechanism, among others, is that the anti-inflammatory activity of IBD model is related to the regulation of inflammatory cytokines such as iNOS, MPO, IL-4, IL-10, EGF, MUC2, IL-6 and so on (Ma et al., 2018). However, this needs to be confirmed by further functional studies, which is beyond the scope of this study.
Mendelian randomization analysis is an effective method to explore causality from exposure to outcome while controlling confounding factors. The MR analysis in this study has the following advantages. First, it is a new attempt to speculate the causal relationship from gut microbiota to IBD, which provides a theoretical basis for the follow-up study of the regulation mechanism of single strain on IBD. Second, it is based on publicly available large-scale GWAS summary statistics, so it provides an effective choice for mining reliable genetic information without additional experimental cost.
Obviously, our study has certain limitations. First, due to limited sample size, the genetic loci identified in gut microbiota GWAS are still limited, which limits the statistical power of MR analysis. Second, MR analysis based on one single IV is less robust, which may bias the interpretation of our findings.
In conclusion, we evaluated the causal relationship from gut microbiota to IBD and identified specific bacterial taxa that may affect the pathogenesis of IBD by a two-sample MR analysis using publicly available GWAS summary statistics. Our results provide a basis for revealing the causal relationship from gut microbiota to IBD, and thus offer new insights into its development and treatment.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: ftp://ftp.sanger.ac.uk/pub/project/humgen/summary_statistics/human/2016-11-07/. The accession code EGAS00001000924. Key data were supplied in the Supplementary Tables 1, 2.
Ethics Statement
Ethical review and approval was not required for the study on human participants in accordance with the Local Legislation and Institutional Requirements. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
Z-JZ: conceptualization, formal analysis, data curation, and writing – original draft. H-LQ: validation and writing – original draft. NZ: formal analysis and writing – original draft. JW: resources and writing – original draft. X-YW: data curation and writing – original draft. BL: data curation, supervision, methodology, writing, and revising. RH: conceptualization, methodology, software, data curation, writing, review, editing, and supervision. All authors contributed to the article and approved the submitted version.
Funding
This study was partially supported by the funding from National Natural Science Foundation of China (81460223 to RH) and a project funded by the Natural Science Foundation of Inner Mongolia Autonomous Region (2016MS0879 to JW). BL was partially supported by grants from Jiangsu Provincial Medical Youth Talent (No. QNRC2016249), Suzhou Science and Technology Bureau (Nos. SYSD2017041 and SYS201788), Gusu Health Top-Notch Youth Talent of Suzhou Health Commission (No. GSWS2019086), and Wujiang District Health Commission (Nos. WWK201609 and WWK201806).
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.
Acknowledgments
We appreciate all the volunteers who participated in this study. We are grateful to the TwinsUK study for releasing the gut microbiota GWAS summary statistics, the UKIBD and USP study and the British IBD study for releasing the IBD GWAS summary statistics.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.631061/full#supplementary-material
Abbreviations
FDR, false discovery rate; GWAS, genome-wide association study; IBD, inflammatory bowel disease; IV, instrumental variable; IVW, inverse-variance weighted; LD, linkage disequilibrium; MGWAS, microbiome genome-wide association study; MR, Mendelian randomization; OTU, operational taxonomic unit; UKB, UK Biobank; USP, understanding social program.
Footnotes
References
Backhed, F., Roswall, J., Peng, Y., Feng, Q., Jia, H., Kovatcheva-Datchary, P., et al. (2015). Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe 17, 852. doi: 10.1016/j.chom.2015.05.012
Bajer, L., Kverka, M., Kostovcik, M., Macinga, P., Dvorak, J., Stehlikova, Z., et al. (2017). Distinct gut microbiota profiles in patients with primary sclerosing cholangitis and ulcerative colitis. World J. Gastroenterol. 23, 4548–4558. doi: 10.3748/wjg.v23.i25.4548
Bengtson, M. B., Solberg, C., Aamodt, G., Sauar, J., Jahnsen, J., Moum, B., et al. (2009). Familial aggregation in Crohn’s disease and ulcerative colitis in a Norwegian population-based cohort followed for ten years. J. Crohns Colitis 3, 92–99. doi: 10.1016/j.crohns.2008.11.002
Bernstein, C. N., Wajda, A., Svenson, L. W., MacKenzie, A., Koehoorn, M., Jackson, M., et al. (2006). The epidemiology of inflammatory bowel disease in Canada: a population-based study. Am. J. Gastroenterol. 101, 1559–1568.
Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44, 512–525. doi: 10.1093/ije/dyv080
Bowden, J., Davey Smith, G., Haycock, P. C., and Burgess, S. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40, 304–314. doi: 10.1002/gepi.21965
Burgess, S., Butterworth, A., and Thompson, S. G. (2013). Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 37, 658–665. doi: 10.1002/gepi.21758
Clarke, S. F., Murphy, E. F., O’Sullivan, O., Lucey, A. J., Humphreys, M., Hogan, A., et al. (2014). Exercise and associated dietary extremes impact on gut microbial diversity. Gut 63, 1913–1920. doi: 10.1136/gutjnl-2013-306541
Cohen, R. D., Yu, A. P., Wu, E. Q., Xie, J., Mulani, P. M., and Chao, J. (2010). Systematic review: the costs of ulcerative colitis in Western countries. Aliment. Pharmacol. Ther. 31, 693–707. doi: 10.1111/j.1365-2036.2010.04234.x
Collado, M. C., Derrien, M., Isolauri, E., de Vos, W. M., and Salminen, S. (2007). Intestinal integrity and Akkermansia muciniphila, a mucin-degrading member of the intestinal microbiota present in infants, adults, and the elderly. Appl. Environ. Microbiol. 73, 7767–7770. doi: 10.1128/aem.01477-07
Cosnes, J., Gower-Rousseau, C., Seksik, P., and Cortot, A. (2011). Epidemiology and natural history of inflammatory bowel diseases. Gastroenterology 140, 1785–1794. doi: 10.1053/j.gastro.2011.01.055
Derrien, M., Collado, M. C., Ben-Amor, K., Salminen, S., and de Vos, W. M. (2008). The Mucin degrader Akkermansia muciniphila is an abundant resident of the human intestinal tract. Appl. Environ. Microbiol. 74, 1646–1648. doi: 10.1128/aem.01226-07
Eckburg, P. B., Bik, E. M., Bernstein, C. N., Purdom, E., Dethlefsen, L., Sargent, M., et al. (2005). Diversity of the human intestinal microbial flora. Science 308, 1635–1638. doi: 10.1126/science.1110591
Goodrich, J. K., Davenport, E. R., Beaumont, M., Jackson, M. A., Knight, R., Ober, C., et al. (2016). Genetic determinants of the gut microbiome in UK twins. Cell Host Microbe 19, 731–743. doi: 10.1016/j.chom.2016.04.017
Guo, X., Zhang, J., Wu, F., Zhang, M., Yi, M., and Peng, Y. (2016). Different subtype strains of Akkermansia muciniphila abundantly colonize in southern China. J. Appl. Microbiol. 120, 452–459. doi: 10.1111/jam.13022
Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., et al. (2018). The MR-Base platform supports systematic causal inference across the human phenome. eLife 7, e34408.
Kevans, D., Silverberg, M. S., Borowski, K., Griffiths, A., Xu, W., Onay, V., et al. (2016). Genetic risk profile in healthy first-degree relatives of Crohn’s disease patients. J Crohns Colitis. 10, 209–215. doi: 10.1093/ecco-jcc/jjv197
Lee, A. A., Rao, K., Limsrivilai, J., Gillilland, M., Malamet, B., Briggs, E., et al. (2020). Temporal gut microbial changes predict recurrent clostridiodes difficile infection in patients with and without ulcerative colitis. Inflamm. Bowel Dis. 26, 1748–1758. doi: 10.1093/ibd/izz335
Lee, K., and Lim, C. Y. (2019). Mendelian randomization analysis in observational epidemiology. J. Lipid Atheroscler. 8, 67–77. doi: 10.12997/jla.2019.8.2.67
Loftus, C. G., Loftus, E. V. Jr., Harmsen, W. S., Zinsmeister, A. R., Tremaine, W. J., and Melton, L. J. III, et al. (2007). Update on the incidence and prevalence of Crohn’s disease and ulcerative colitis in Olmsted County, Minnesota, 1940-2000. Inflamm. Bowel Dis. 13, 254–261. doi: 10.1002/ibd.20029
Loftus, E. V. Jr. (2004). Clinical epidemiology of inflammatory bowel disease: incidence, prevalence, and environmental influences. Gastroenterology 126, 1504–1517. doi: 10.1053/j.gastro.2004.01.063
Ma, X., Hu, Y., Li, X., Zheng, X., Wang, Y., Zhang, J., et al. (2018). Periplaneta americana ameliorates dextran sulfate sodium-induced ulcerative colitis in rats by Keap1/Nrf-2 activation, intestinal barrier function, and gut microbiota regulation. Front. Pharmacol. 9:944. doi: 10.3389/fphar.2018.00944
Magro, F., Rodrigues, A., Vieira, A. I., Portela, F., Cremers, I., Cotter, J., et al. (2012). Review of the disease course among adult ulcerative colitis population-based longitudinal cohorts. Inflamm. Bowel Dis. 18, 573–583. doi: 10.1002/ibd.21815
Matsuoka, K., and Kanai, T. (2015). The gut microbiota and inflammatory bowel disease. Semin. Immunopathol. 37, 47–55.
Matsuoka, K., Kobayashi, T., Ueno, F., Matsui, T., Hirai, F., Inoue, N., et al. (2018). Evidence-based clinical practice guidelines for inflammatory bowel disease. J. Gastroenterol. 53, 305–353.
Molodecky, N. A., Soon, I. S., Rabi, D. M., Ghali, W. A., Ferris, M., Chernoff, G., et al. (2012). Increasing incidence and prevalence of the inflammatory bowel diseases with time, based on systematic review. Gastroenterology 142, 46–54.e42; quiz e30.
Nishida, A., Inoue, R., Inatomi, O., Bamba, S., Naito, Y., and Andoh, A. (2018). Gut microbiota in the pathogenesis of inflammatory bowel disease. Clin. J. Gastroenterol. 11, 1–10.
O’Hara, A. M., and Shanahan, F. (2006). The gut flora as a forgotten organ. EMBO Rep. 7, 688–693. doi: 10.1038/sj.embor.7400731
Paramsothy, S., Kamm, M. A., Kaakoush, N. O., Walsh, A. J., van den Bogaerde, J., Samuel, D., et al. (2017). Multidonor intensive faecal microbiota transplantation for active ulcerative colitis: a randomised placebo-controlled trial. Lancet 389, 1218–1228. doi: 10.1016/s0140-6736(17)30182-4
Quigley, E. M. M. (2017). Microbiota-brain-gut axis and neurodegenerative diseases. Curr. Neurol. Neurosci. Rep. 17, 94.
Sanna, S., van Zuydam, N. R., Mahajan, A., Kurilshikov, A., Vich Vila, A., Võsa, U., et al. (2019). Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases. Nat. Genet. 51, 600–605. doi: 10.1038/s41588-019-0350-x
Sasaki, K., Inoue, J., Sasaki, D., Hoshi, N., Shirai, T., Fukuda, I., et al. (2019). Construction of a model culture system of human colonic microbiota to detect decreased lachnospiraceae abundance and butyrogenesis in the feces of ulcerative colitis patients. Biotechnol. J. 14, e1800555.
Strandwitz, P. (2018). Neurotransmitter modulation by the gut microbiota. Brain Res. 1693, 128–133. doi: 10.1016/j.brainres.2018.03.015
Turpin, W., Goethel, A., Bedrani, L., and Croitoru Mdcm, K. (2018). Determinants of IBD heritability: genes, bugs, and more. Inflamm. Bowel Dis. 24, 1133–1148. doi: 10.1093/ibd/izy085
van Passel, M. W., Kant, R., Zoetendal, E. G., Plugge, C. M., Derrien, M., Malfatti, S. A., et al. (2011). The genome of Akkermansia muciniphila, a dedicated intestinal mucin degrader, and its use in exploring intestinal metagenomes. PLoS One 6:e16876. doi: 10.1371/journal.pone.0016876
Verbanck, M., Chen, C. Y., Neale, B., and Do, R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50, 693–698. doi: 10.1038/s41588-018-0099-7
Keywords: mendelian randomization, gut microbiota, inflammatory bowel disease, ulcerative colitis, causal relationship
Citation: Zhang Z-J, Qu H-L, Zhao N, Wang J, Wang X-Y, Hai R and Li B (2021) Assessment of Causal Direction Between Gut Microbiota and Inflammatory Bowel Disease: A Mendelian Randomization Analysis. Front. Genet. 12:631061. doi: 10.3389/fgene.2021.631061
Received: 19 November 2020; Accepted: 04 January 2021;
Published: 18 February 2021.
Edited by:
Lei Zhang, Soochow University, ChinaCopyright © 2021 Zhang, Qu, Zhao, Wang, Wang, Hai and Li. 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: Bin Li, binli10tigarres@sina.com; Rong Hai, aodaili12@hotmail.com