- Open Access
Microbial imbalance in inflammatory bowel disease patients at different taxonomic levels
Gut Pathogens volume 12, Article number: 1 (2020)
Inflammatory bowel disease (IBD), is a debilitating group of chronic diseases including Crohn’s Disease (CD) and ulcerative colitis (UC), which causes inflammation of the gut and affects millions of people worldwide. At different taxonomic levels, the structure of the gut microbiota is significantly altered in IBD patients compared to that of healthy individuals. However, it is unclear how these IBD-affected bacterial groups are related to other common bacteria in the gut, and how they are connected across different disease conditions at the global scale.
In this study, using faecal samples from patients with IBD, we show through diversity analysis of the microbial community structure based on the 16S rRNA gene that the gut microbiome of IBD patients is less diverse compared to healthy individuals. Furthermore, we have identified which bacterial groups change in abundance in both CD and UC compared to healthy controls. A substantial imbalance was observed across four major bacterial phyla including Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria, which together constitute > 98% of the gut microbiota. Next, we reconstructed a bacterial family co-abundance network based on the correlation of abundance profiles obtained from the public gut microbiome data of > 22,000 samples of faecal and gut biopsies taken from both diseased and healthy individuals. The data was compiled using the EBI metagenomics database (Mitchell et al. in Nucleic Acids Res 46:D726–D735, 2018). By mapping IBD-altered bacterial families to the network, we show that the bacterial families which exhibit an increased abundance in IBD conditions are not well connected to other groups, implying that these families generally do not coexist together with common gut organisms. Whereas, the bacterial families whose abundance is reduced or did not change in IBD conditions compared to healthy conditions are very well connected to other bacterial groups, suggesting they are highly important groups of bacteria in the gut that can coexist with other bacteria across a range of conditions.
IBD patients exhibited a less diverse gut microbiome compared to healthy individuals. Bacterial groups which changed in IBD patients were found to be groups which do not co-exist well with common commensal gut bacteria, whereas bacterial groups which did not change in patients with IBD were found to commonly co-exist with commensal gut microbiota. This gives a potential insight into the dynamics of the gut microbiota in patients with IBD.
Inflammatory bowel disease (IBD), a group of chronic intestinal disorders including Crohn’s disease (CD) and ulcerative colitis (UC), causes inflammation of the gut and affects millions of people worldwide [1,2,3,4]. Both CD and UC diseases are differentiated by their location and levels of inflammation in the gastrointestinal (GI) tract. UC mostly involves inflammation to the rectum and colon, whereas CD most often affects the terminal ileum and colon though in some cases it can affect any part of the GI tract [2, 5]. Currently, there is no full cure for IBD, but different treatments such as taking steroids, immunosuppressants, liquid diet or surgery can help in reducing the symptoms . To date, the exact cause of IBD is not understood, however, a combination of genetic variants, environmental factors, deregulated host immune system, and gut microbiota dysbiosis is associated with IBD [6,7,8,9,10,11,12,13,14,15,16].
More than 215 IBD-associated loci have been identified so far from various genome-wide association studies (GWAS) . It has been reported that a large percentage (~ 30%) of these loci are common between CD and UC, showing involvement of common biological processes in both conditions [14, 17]. Moreover, these IBD-associated loci are mostly involved with immune system deregulation, a process which the gut microbiome has also been implicated in . The gut microbiota, which has a complex community of a hundred trillion bacterial and archaeal cells comprising more than a thousand species, provides benefits to the host such as short-chain fatty acids (SCFA) and amino acids, metabolism of undigested carbohydrate, and stimulation of the immune system [19, 20]. In patients with IBD, the structure and the composition of the gut microbiota is severely altered compared to that of a healthy condition [10, 13, 18, 21]. Previous work has reported imbalances in IBD patients for the Firmicutes and Bacteroidetes at the phylum level, and Ruminococcaceae, Veillonellaceae, Christensenellaceae, Bacteroidaceae, and Rikenellaceae at the family level. However, there is a large degree of variability across studies, with many reporting contradictory findings. In particular, it is unclear what the relationship is between microbial groups when there is inflammation of the gut epithelium during various diseased states. There is furthermore, a general knowledge gap in establishing the relationships between microbial groups across different disease conditions.
The aim of the current work was to investigate the relationships between changing microbial groups in IBD. In particular, we wanted to understand which microbial groups differ during IBD, and how these groups differ in co-abundance patterns across a variety of diseases at the global scale. To do this, we initially investigated the gut microbial imbalance, at different taxonomic levels for healthy volunteers and CD and UC patients. We next reconstructed a network of the co-abundance patterns of different bacterial groups using publicly available data from a variety of studies at a global scale. Our results indicate that the bacterial groups which increase in abundance during IBD are specific to both CD and UC conditions. In comparison, bacterial groups which did not change in abundance during different disease states are well connected in global networks, giving us a better understanding of the dynamics of the microbiome in both health and disease.
Results and discussion
Gut microbial richness in IBD patients
We collected faecal samples from 30 individuals (9 CD patients, 11 UC patients and 10 healthy volunteers) and performed 16S rRNA taxonomic profiling to understand changes in community structure during disease with resulting data analysed using the DADA2 pipeline. Amplicon sequence variants (ASVs) were used as a measure of diversity. As previously reported, species richness in the gut of IBD patients (both CD and UC) was lower than that of healthy volunteers [10, 13]. Moreover, within the IBD patients groups, the gut of CD patients exhibit substantially lower species richness than that of UC patients (Fig. 1a). As a measure of diversity, we identified a total of 2261 ASVs, of which 81% belong to the Firmicutes, 9.25% Bacteroidetes, 5.13% Proteobacteria, and 3.14% ASVs are from Actinobacteria. Combined, these four phyla constitute more than 98% of the total identified ASVs [22, 23]. To understand how this compared to other studies, we investigated the global microbial species-richness in the human gut across > 22,000 samples from 113 different studies from a variety of conditions (EBI metagenomics accession numbers  in Additional file 1: Table S1). We have considered only the known species in each study and made a unique list of gut bacterial species across studies. Similar to our experimental observations, the four phyla including Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria comprised > 94% of the total diversity in the gut, with the Firmicutes being the most species-rich phylum across conditions (Fig. 1b). However, the number of species identified as belonging to the Firmicute Phyla in our dataset (81% of all identified species from 20 IBD and 10 healthy condition samples) is substantially higher than what is usually reported at a global level (38.18%). As the observation comes from 2/3 of IBD patient and 1/3 of healthy control samples, this increased level of Firmicutes species-richness was attributed to the high number of IBD samples. Analysis of global studies for species-richness of Proteobacteria demonstrated this phyla usually accounts for 31.18% of all species, however this was substantially lower in our dataset (5.13%). Finally, the richness of Bacteroidetes was also reduced in our dataset compared to the global datasets (9.24% in our study compared to 14.35% globally). We find a similar observation when considering all OTUs from individual EBI gut microbiome studies (Additional file 2: Figure S1).
After identifying the differences at different taxonomic levels for each disease condition, for future work it is crucial to understand the reasons for such dysbiosis and whether they are causative or consequential of disease. Studies suggest that metabolic dependency  and nutritional preferences  between microorganisms are driving forces in microbial community formation. For example, metabolic cooperation between bacteria is crucial to microbial assemblages and changes to this could cause shifts across the whole community. For future work, it would be interesting to investigate the microbial metabolic interactions during disease and how this compares to a healthy gut.
The gut microbial abundance at different taxonomic levels in IBD patients
Compared to the healthy controls, both IBD patient groups (CD and UC patients) demonstrated strong microbial imbalance at different taxonomic levels (Fig. 2). At the phylum level, both IBD conditions exhibit an increased abundance of Firmicutes and Actinobacteria, relative to the controls. In particular, the abundance of this phyla during UC, was far greater than CD or healthy controls. For the two other dominant bacterial phyla, the Proteobacteria and Bacteroidetes, the abundance profiles varied across disease conditions (Fig. 2a). In CD, the abundance of Bacteroidetes, which is often associated with a healthy gut, was deceased 2.4-fold, whereas the abundance of Proteobacteria, a phyla associated including wide variety of pathogens, was increased 3.8 fold. Interestingly, for patients with UC the abundance of Proteobacteria was decreased (3.4 fold) relative to controls and there was no significant differences in levels Bacteroidetes . Several studies have reported the gut microbial imbalances for IBD, however, the imbalance at the level of different phylum is variable across studies [10, 13, 18, 21, 27, 28]. This could likely be a reflection on the lack of standardisation across microbiome techniques, or perhaps due to a heterogeneity in the microbiome associated with disease.
We further investigated how different taxonomic levels belonging to each of the main four phyla, Firmicutes, Proteobacteria, Bacteroidetes, and Actinobacteria, were changed during IBD. For Firmicutes, the most abundant phylum in the gut in all conditions, we observed four classes and nine different families which changed in abundance relative to healthy controls. For CD patients, the abundance of two classes including Clostridia and Erysipelotrichia, was reduced, and three families including Ruminococcaceae, Christensenellaceae and Erysipelotrichaceae were reduced relative to healthy controls. The level of two other two classes such as Negativicutes and Bacilli (obligately aerobic) and five families including Veillonellaceae, Lactobacillaceae, Acidaminococcaceae, Streptococcaceae and Peptostreptococcaceae was increased, similar to the imbalance in their parent phylum Firmicutes. Interestingly, Lachnospiraceae, the most abundant Firmicutes family, was at a similar level to the control. For UC patients, the abundance of the Firmicute classes Clostridia, Negativicutes and Bacilli, and Firmicute families Ruminococcaceae, Lachnospiraceae, Veillonellaceae, Streptococcaceae and Peptostreptococcaceae were increased. The Erysipelotrichia class and Erysipelotrichaceae family were the same as the controls, whereas, three families including Acidaminococcaceae, Christensenellaceae and Lactobacillaceae, were reduced in abundance. For Bacteroidetes, which is the only reduced phylum in CD, we observed a reduced abundance in the Bactersoidia class and Bacteroidia families Bacteroidaceae and Rikenellaceae. The Prevotellaceae family, in particular, was increased in CD patients. For UC patients, the only families to change of the Bacteroidetes was the Rikenellaceae and Tannerellaceae families which were decreased in abundance, and the Prevotellaceae which increased in abundance as with CD. Finally, for the Proteobacteria phylum, we observed an imbalance in Enterobacteriaceae and Burkholderiaceae families, with the abundance level of Enterobacteriaceae increased in CD patients and decreased in UC patients compared to the controls. Burkholderiaceae abundance was increased for both CD and UC patients. Finally, for the Actinobacteria phylum, the abundance level of both the class Coriobacteriia and family Coriobacteriaceae was increased in both IBD conditions relative to the controls, whereas the class Actinobacteria and family Bifidobacteriaceae was reduced [13, 29]. In summary, we demonstrate that multiple families of a class, classes of a phyla differ between both the IBD conditions, and between each IBD condition and healthy control. This suggests that changes in one bacterial family has consequences for others. To investigate this further, we used co-occurrence network analysis to establish patterns of how bacterial groups increase and decrease across global studies.
Global co-abundance in the gut of different bacterial families
The gut microbiota abundance profiles from > 22,000 samples across a variety of conditions from 113 different studies was analysed to explore how different bacterial groups change across global studies. Using the Pearson’s correlation test, we built a network of significantly co-abundant (Pearson’s correlation coefficient > 0.3 and P value < 1e−10) bacterial families across a range of conditions collected from global studies (Fig. 3ai). We observed that the majority of bacterial families in the network belonged to the phyla Proteobacteria, Actinobacteria, Firmicutes and Bacteroidetes (Fig. 3aii). In the global gut bacterial family–family co-abundance network (Fig. 3ai), family nodes of 45% of the total connections are from the same phylum, compared to a random network of the same size where percentage connections were much lower (~ 30%) (Fig. 3aiii). This suggests that groups from the same phylum, which likely have similar metabolic requirements, are likely to change abundance as a collective. To understand the family level microbial imbalance during IBD, we further considered a subnetwork of the global family co-abundance network, where at least one family node was from the most abundant bacterial families in either CD, UC or healthy subjects. In this subnetwork, we highlighted bacterial families which were changed during CD or UC, compared to the healthy (Fig. 3b, c).
In CD, seven bacterial families including Coriobacteriaceae, Prevotellaceae, Burkholderiaceae, Veillonellaceae, Streptococcaceae, Pseudomonadaceae and Acidaminococcaceae have increased abundance compared to the healthy controls (Fig. 2c), however, only two families including Prevotellaceae and Veillonellaceae are connected in the global network (Fig. 3b). In contrast, the level of five other families including Erysipelotrichaceae, Christensenellaceae, Ruminococcaceae, Bacteroidaceae and Rikenellaceae, were well connected in the global network (Fig. 3b) and had reduced abundance in CD (Fig. 2c). Similarly, for UC, families with an increased abundance in UC were less well connected on the global scale (Fig. 3c). This suggests that bacterial groups which increase in abundance during IBD, are not typically associated with the healthy gut microbiome, nor do they commonly co-exist with commensals observed in the healthy gut. Furthermore, families which had increased abundance levels in healthy conditions compared to CD are very well connected, suggesting microbes in the gut of healthy individuals exist as a co-operative microbial assemblage. In particular, the connection between families such as Bacteroidaceae and Ruminococcaceae in the co-abundance network indicates that they may coexist together in across conditions, potentially due to similarities in physiology or the presence of metabolically cooperating species. For future work, it would be highly interesting to examine species of these families and investigate the relationships between these organisms.
In summary, our analysis demonstrates that IBD patients (both CD and UC) and healthy volunteers have reduced species richness, and imbalances in families, classes, and phyla, relative to healthy volunteers. Four bacterial phyla including Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria comprised > 98% of the species in this study. To understand how bacteria assemblages depend on co-operation, we reconstructed a large co-abundance network based on the public gut microbiome data of > 22,000 samples. From this we demonstrated that the bacterial families which have an increased abundance level in IBD conditions are not well connected to other bacterial groups in the global family co-abundance network. This suggests that these bacteria do not co-exist with healthy gut microbial commensals and supports the concept that healthy assemblages are dependent on metabolic co-operation, due to the high connectivity of bacterial groups found in healthy conditions across > 22,000 samples.
Sample collection and DNA extraction
Samples were collected from patients in standard 300 ml sterilin tubes and frozen immediately in −80 °C. Patients were asked to produce the first-morning sample for consistency and to avoid alcohol the previous 24 h. Samples were thawed and DNA was extracted using FastDNA Spin Kit for Soil (MPBiomedicals)  as per the manufacturer’s instructions.
16S rRNA sequencing
454 pyrosequencing using 16S universal eubacterial primers 27F and 534R  was performed by Molecular Research (MRDNA), Shallowater, Texas, using an adapted protocol developed in . Number of reads per sample ranged from 6936 to 100,972, with an average of 38,931 reads per sample.
Bioinformatic analysis of 16S rRNA sequencing data
16S rRNA high-throughput sequencing data was analysed by following the workflow from Callahan et al. . Quality checking, filtering and trimming of fastq files were performed by functions from the dada2 package in R . After filtering the reads, high-resolution Amplicon Sequence Variants (ASVs) were inferred using dada function . ASVs are a higher-resolution analogue of the traditional OTUs. Chimeric sequences were removed and taxonomy assigned to ASVs based on the naive Bayesian classifier method with silva_nr_v132_train_set.fa as the training set . Species-richness and alpha diversity (Shannon) were analysed by plot_richness function from the phyloseq package in R . To make the rarefaction species richness curve, ‘rarecurve’ function from the vegan package  in R was used.
Statistical analysis: bacterial family co-abundance network based on microbiome data
Taxonomic assignments, containing a detailed taxonomy and abundance data of OTUs or ASVs in samples, of 113 gut microbiome studies, covering more than 22,000 samples, were downloaded from the EBI metagenomics database . The data were then parsed and tables containing bacterial abundance from different phyla, classes and families were generated. Abundance at the level of phylum was then used for the global gut microbial abundance. The family-level bacterial abundance was used to construct the bacterial family–family coexistence network. For each pair of bacterial families, Pearson’s correlation test was performed. Family nodes were connected when P-value < 1e−10 and Pearson’s correlation coefficient > 0.3.
Availability of data and materials
Sequence data of the study has been submitted to European Nucleotide Archive (ENA) (https://www.ebi.ac.uk/ena/submit/sra/#studies) under the Accession Number PRJEB33711
Mitchell AL, Scheremetjew M, Denise H, Potter S, Tarkowska A, Qureshi M, et al. EBI Metagenomics in 2017: enriching the analysis of microbial communities, from sequence reads to assemblies. Nucleic Acids Res. 2018;46:D726–35.
Abraham C, Cho JH. Inflammatory bowel disease. N Engl J Med. 2009;361:2066–78.
Kaplan GG. The global burden of IBD: from 2015 to 2025. Nat Rev Gastroenterol Hepatol. 2015;12:720–7.
Mozdiak E, O’Malley J, Arasaradnam R. Inflammatory bowel disease. BMJ. 2015. h4416.
Baumgart DC, Sandborn WJ. Crohn’s disease. Lancet. 2012;380:1590–605.
McGovern DPB, Kugathasan S, Cho JH. Genetics of inflammatory bowel diseases. Gastroenterology. 2015;149(1163–76):e2.
Luo Y, de Lange KM, Jostins L, Moutsianas L, Randall J, Kennedy NA, et al. Exploring the genetic architecture of inflammatory bowel disease by whole-genome sequencing identifies association at ADCY7. Nat Genet. 2017;49:186–92.
Molodecky NA, Kaplan GG. Environmental risk factors for inflammatory bowel disease. Gastroenterol Hepatol. 2010;6:339–46.
Ananthakrishnan AN. Environmental risk factors for inflammatory bowel disease. Gastroenterol Hepatol. 2013;9:367–74.
Sheehan D, Moran C, Shanahan F. The microbiota in inflammatory bowel disease. J Gastroenterol. 2015;50:495–507.
Goldsmith JR, Sartor RB. The role of diet on intestinal microbiota metabolism: downstream impacts on host immune function and health, and therapeutic implications. J Gastroenterol. 2014;49:785–98.
Liu JZ, van Sommeren S, Huang H, Ng SC, Alberts R, Takahashi A, et al. Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations. Nat Genet. 2015;47:979–86.
Nishida A, Inoue R, Inatomi O, Bamba S, Naito Y, Andoh A. Gut microbiota in the pathogenesis of inflammatory bowel disease. Clin J Gastroenterol. 2018;11:1–10.
Khor B, Gardet A, Xavier RJ. Genetics and pathogenesis of inflammatory bowel disease. Nature. 2011;474:307–17.
Lee SH, Kwon JE, Cho M-L. Immunological pathogenesis of inflammatory bowel disease. Intest Res. 2018;16:26–42.
Shaw KA, Cutler DJ, Okou D, Dodd A, Aronow BJ, Haberman Y, et al. Genetic variants and pathways implicated in a pediatric inflammatory bowel disease cohort. Genes Immun. 2019;20:131–42.
Ek WE, D’Amato M, Halfvarson J. The history of genetics in inflammatory bowel disease. Ann Gastroenterol Hepatol. 2014;27:294–303.
Belkaid Y, Hand TW. Role of the Microbiota in Immunity and Inflammation. Cell. 2014;157:121–41.
Nishida AH, Ochman H. A great-ape view of the gut microbiome. Nat Rev Genet. 2019;20:195–206.
Cani PD. Human gut microbiome: hopes, threats and promises. Gut. 2018;67:1716–25.
Franzosa EA, Sirota-Madi A, Avila-Pacheco J, Fornelos N, Haiser HJ, Reinker S, et al. Gut microbiome structure and metabolic activity in inflammatory bowel disease. Nat Microbiol. 2019;4:293–305.
Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, et al. Diversity of the human intestinal microbial flora. Science. 2005;308:1635–8.
Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science. 2008;320:1647–51.
Zelezniak A, Andrejev S, Ponomarova O, Mende DR, Bork P, Patil KR. Metabolic dependencies drive species co-occurrence in diverse microbial communities. Proc Natl Acad Sci USA. 2015;112:6449–54.
Tramontano M, Andrejev S, Pruteanu M, Klünemann M, Kuhn M, Galardini M, et al. Nutritional preferences of human gut bacteria reveal their metabolic idiosyncrasies. Nat Microbiol. 2018;3:514–22.
Gophna U, Sommerfeld K, Gophna S, Doolittle WF, van Zanten SJO. Differences between tissue-associated intestinal microfloras of patients with Crohn’s disease and ulcerative colitis. J Clin Microbiol. 2006;44:4136–41.
Vester-Andersen MK, Mirsepasi-Lauridsen HC, Prosberg MV, Mortensen CO, Träger C, Skovsen K, et al. Increased abundance of proteobacteria in aggressive Crohn’s disease seven years after diagnosis. Sci Rep. 2019;9:13473.
Frank DN, St Amand AL, Feldman RA, Boedeker EC, Harpaz N, Pace NR. Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases. Proc Natl Acad Sci USA. 2007;104:13780–5.
LeBlanc JG, Laiño JE, del Valle MJ, Vannini V, van Sinderen D, Taranto MP, et al. B-group vitamin production by lactic acid bacteria—current knowledge and potential applications. J Appl Microbiol. 2011;111:1297–309.
King HC, Murphy A, James P, Travis E, Porter D, Sawyer J, et al. Performance of a Noninvasive test for detecting Mycobacterium bovis shedding in European Badger (Meles meles) populations. J Clin Microbiol. 2015;53:2316–23.
Tamaki H, Wright CL, Li X, Lin Q, Hwang C, Wang S, et al. Analysis of 16S rRNA amplicon sequencing options on the Roche/454 next-generation titanium sequencing platform. PLoS ONE. 2011;6:e25263.
Dowd SE, Sun Y, Wolcott RD, Domingo A, Carroll JA. Bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP) for microbiome studies: bacterial diversity in the ileum of newly weaned Salmonella-infected pigs. Foodborne Pathog Dis. 2008;5:459–72.
Callahan BJ, Sankaran K, Fukuyama JA, McMurdie PJ, Holmes SP. Bioconductor workflow for Microbiome Data analysis: from raw reads to community analyses. F1000 Res. 2016;5:1492.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217.
Oksanen J, Guillaume Blanchet F, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community Ecology Package. 2018. https://CRAN.R-project.org/package=vegan.
We thank BBSRC (BBSRC BB/L027801/1) and FAPIC (EU H2020 Proposal Number: SEP-210176939 2015-2020) for funding the project. MTA was Funded by the Warwick medical school.
EMHW and GCAA received financial support from BBSRC (BBSRC BB/L027801/1), Warwick Medical School and AM from FAPIC—Fast Assay for Pathogen Identification and Characterisation EU H2020 Proposal Number: SEP-210176939 2015-2020.
Ethics approval and consent to participate
Patients were recruited as part of the Famished Multicentre Study (Warwickshire local research Ethics No: 09/H1211/38).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of gut microbiome studies from EBI metagenomics database.
Phylum level richness in the gut microbiota from our samples (shown in green) compared to the gut microbial species richness obtained across studies, each with more than 100 OTUs (total 81 study), from a variety of conditions (shown in black).
About this article
Cite this article
Alam, M.T., Amos, G.C.A., Murphy, A.R.J. et al. Microbial imbalance in inflammatory bowel disease patients at different taxonomic levels. Gut Pathog 12, 1 (2020). https://doi.org/10.1186/s13099-019-0341-6
- Inflammatory bowel disease
- Crohn’s disease
- Ulcerative colitis
- Gut microbiota
- Microbial imbalance