Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 18 February 2021
Sec. Statistical Genetics and Methodology
This article is part of the Research Topic Mendelian Randomization: Approach and Applications View all 15 articles

Assessment of Causal Direction Between Gut Microbiota and Inflammatory Bowel Disease: A Mendelian Randomization Analysis

\r\nZi-Jia Zhang,,Zi-Jia Zhang1,2,3Hong-Lei QuHong-Lei Qu4Na ZhaoNa Zhao3Jing WangJing Wang3Xiu-Yan WangXiu-Yan Wang3Rong Hai,*Rong Hai3,5*Bin Li*Bin Li1*
  • 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
www.frontiersin.org

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.

TABLE 1
www.frontiersin.org

Table 1. MR analysis of gut microbiota on IBD in both discovery and replication samples.

TABLE 2
www.frontiersin.org

Table 2. Instrumental variables used in both discovery and replication studies.

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

  1. ^ https://github.com/MRCIEU/TwoSampleMR
  2. ^ https://github.com/rondolab/MR-PRESSO

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsuoka, K., and Kanai, T. (2015). The gut microbiota and inflammatory bowel disease. Semin. Immunopathol. 37, 47–55.

Google Scholar

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.

Google Scholar

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.

Google Scholar

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.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Quigley, E. M. M. (2017). Microbiota-brain-gut axis and neurodegenerative diseases. Curr. Neurol. Neurosci. Rep. 17, 94.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

Strandwitz, P. (2018). Neurotransmitter modulation by the gut microbiota. Brain Res. 1693, 128–133. doi: 10.1016/j.brainres.2018.03.015

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Wu, X., Cao, S., Wang, L., Wang, D., Yang, H., et al. (2016). Caffeic acid ameliorates colitis in association with increased Akkermansia population in the gut microbiota of mice. Oncotarget 7, 31790–31799. doi: 10.18632/oncotarget.9306

PubMed Abstract | CrossRef Full Text | Google Scholar

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, China

Reviewed by:

Feng Zhang, Xi’an Jiaotong University, China
Sumei Xiao, Sun Yat-sen University, China

Copyright © 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

Disclaimer: 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.