The Effects of Bowel Preparation on Microbiota-Related Metrics Differ in Health and in Inflammatory Bowel Disease and for the Mucosal and Luminal Microbiota Compartments

Rima M Shobar, Suresh Velineni, Ali Keshavarzian, Garth Swanson, Mark T DeMeo, Joshua E Melson, John Losurdo, Philip A Engen, Yan Sun, Lars Koenig, Ece A Mutlu, Rima M Shobar, Suresh Velineni, Ali Keshavarzian, Garth Swanson, Mark T DeMeo, Joshua E Melson, John Losurdo, Philip A Engen, Yan Sun, Lars Koenig, Ece A Mutlu

Abstract

Objectives: Bowel preparations (BPs) taken before colonoscopy may introduce a confounding effect on the results of gastrointestinal microbiota studies. This study aimed to determine the effect of bowel preparation on the mucosa-associated and luminal colonic microbiota in healthy subjects (HC) and inflammatory bowel disease (IBD) patients.

Methods: Biopsy samples (n=36) and fecal samples (n=30) were collected from 10 HC and 8 IBD subjects pre- and post-BP. 16S rRNA gene was pyrosequenced using 454 Titanium protocols. We compared the differences between the pre- and post-BP samples (i.e., comparisons-across-bowel-prep); we examined the effect of BP on the expected separation of the mucosal vs. the luminal compartments (i.e., comparisons-across-compartments). Last, we compared the baseline differences between the HC vs. IBD groups (a secondary analysis), and examined whether the differences between the HC vs. IBD changed after BP.

Results: In comparisons-across-bowel-prep, the Shannon's index (SI) decreased only in the biopsy samples of IBD subjects post-BP (P=0.025) and phylogenetic diversity-whole tree (PD-WT) metric decreased in biopsy samples of HC subjects post-BP (P=0.021). In secondary comparisons, the subtle differences between the fecal samples of the HC vs. IBD groups, in terms of evenness and the SI, were not apparent post-BP. In terms of β-diversity, in comparisons-across-bowel-prep, the proportion of shared operational taxonomic units (OTUs) in pre- and post-BP samples was low (~30%) and unweighted Unifrac distances between pre- and post-BP specimens ranged from 0.52 to 0.66. HC biopsies were affected more than IBD biopsies with BP (P=0.004). In comparisons-across-compartments, the proportion of shared OTUs between biopsy and fecal samples increased and Unifrac distances decreased post-BP in IBD subjects, reducing the differences between the mucosal and luminal compartments of the gut microbiota. Interindividual differences in Unifrac distances were preserved even with BP effects, although the effects were greater on weighted Unifrac distances. Bacteroidetes and its subtypes increased post-BP in both the luminal and mucosal compartments.

Conclusions: Bowel preparations affect the composition and diversity of the fecal and luminal microbiota in the short term, introducing potential bias into experiments examining the gut microbiota. The magnitude of the effect of BP is not greater than that of interindividual variation. Both the luminal and mucosal compartments of the gut microbiota get affected, and samples from controls and IBD subjects may get affected differently. Studies of the colonic microbiota should take into account the direction and the magnitude of the change introduced by BP during the design stage of the experiments, and consider sample sizes so that potential bias is minimized.

Figures

Figure 1
Figure 1
Operational taxonomic unit (OTU)-based diversity indices of pre- and post-bowel preparation (BP) samples. (a) Richness of OTUs pre- and post-BP. Each marker represents a single sample. The lines represent the mean value for each sample group. The color red represents biopsy samples of healthy controls (HCs). The color green represents fecal samples of HCs. The color blue represents biopsy samples of inflammatory bowel disease (IBD) subjects. The color yellow represents the fecal samples of IBD. In comparisons-across-bowel-prep, a trend towards reduced richness in samples collected post-BP is seen in both HC and IBD groups, but the differences are not statistically significant (P>0.05). (b) The evenness pre- and post-BP. There was a difference in the evenness between pre-BP fecal samples of HC vs. IBD subjects (P=0.033), and this difference was not apparent after BP. (c) The Shannon's diversity index pre- and post-BP. The baseline (pre-BP) Shannon's index for the fecal samples of IBD was significantly lower compared with those for HCs (P=0.04). This difference was not apparent after BP. Although there was a trend towards reduced Shannon's indices in the post-BP samples (compared with pre-BP samples) in general, this reduction was only significant in the biopsy samples of the IBD subjects (P=0.025, paired t-test). (d) The phylogenetic diversity-whole tree (PD-WT) index pre- and post-BP. The effect of BP was more pronounced in the biopsy samples, especially in the HC group (P=0.021, related-samples sign test). The apparent lesser effect on fecal samples may be explained by the greater inconsistency in the direction of response to BP among the fecal samples.
Figure 2
Figure 2
The proportion of shared operational taxonomic units (OTUs) between samples. (a) Shared OTUs between the pre- and post-bowel prep (BP) samples in the healthy control (HC) and inflammatory bowel disease (IBD) groups. Each marker plotted represents the proportion of shared OTUs between paired samples (before and after BP) of the same subject. Red markers represent shared OTUs in HCs. Blue markers represent shared OTUs in the IBD subjects. The lines represent the mean values for the proportion of shared OTUs. The proportion of shared OTUs between pre- and post-bowel prep samples in all sample sets was low (~20–30%). (b) Shared OTUs between the mucosal and luminal compartments (i.e., biopsy and fecal samples) collected at the same visit in the HC and IBD groups. Each marker plotted represents the proportion of shared OTUs between paired samples (fecal and biopsy samples in the same visit) of the same subject. Red markers represent shared OTUs in HCs. Blue markers represent shared OTUs in the IBD subjects. There was an increase in the proportion of shared OTUs between the biopsy and fecal samples after BP (P=0.016).
Figure 3
Figure 3
Unifrac distances between paired samples. (a) The unweighted Unifrac distances between the pre- and post-bowel prep (BP) samples in the healthy control (HC) and inflammatory bowel disease (IBD) groups. Each marker plotted represents the Unifrac distance between paired pre- and post-BP samples of the same type for the same subject. Red markers represent Unifrac distances in the HCs. Blue markers represent Unifrac distances in the IBD subjects. The lines represent the mean Unifrac distances. The unweighted Unifrac distances between the pre- and post-BP biopsy samples was significantly greater in the HC group, when compared with the IBD group (P=0.04). (b) The weighted Unifrac distances between the pre- and post-BP samples in the HC and IBD groups. (c) The unweighted Unifrac distances between the mucosal and luminal compartments (i.e., biopsy and fecal samples) collected at the same visit in the HC and IBD groups. Overall, the mean-unweighted Unifrac distance between the biopsy and fecal samples decreased after bowel prep (P=0.003). This reduction in unweighted Unifrac distances between sample types because of bowel prep was significant in the IBD group (P=0.03), but not in the HC group (P=0.068). (d) The weighted Unifrac distances between the mucosal and luminal compartments (i.e., biopsy and fecal samples) collected at the same visit in the HC and IBD groups. Although there was a trend towards reduced weighted Unifrac distances between the mucosal and luminal compartments after BP, this reduction was only significant in the IBD group (P=0.036).
Figure 4
Figure 4
Unifrac distances between samples from the same subject (i.e., intra-individual differences) vs. those from other subjects (i.e., inter-individual differences). (a) The unweighted Unifrac distances between samples of the same subject (intra-individual distances) vs. the unweighted Unifrac distances between samples of different subjects (inter-individual distances) in the healthy control (HC) group (P⩽0.0001). (b) The weighted Unifrac distances between samples of the same subject (intra-individual distances) vs. the weighted Unifrac distances between samples of different subjects (inter-individual distances) in the HC group (P=0.0047). (c) The unweighted Unifrac distances between samples of the same subject (intra-individual distances) vs. the unweighted Unifrac distances between samples of different subjects (inter-individual distances) in the inflammatory bowel disease (IBD) group (P≤0.0001). (d) The weighted Unifrac distances between samples of the same subject (intra-individual distances) vs. the weighted Unifrac distances between samples of different subjects (inter-individual distances) in the IBD group (P≤0.0001). The inter-individual distances were significantly greater than the intra-individual distances for all comparisons. The line in each box represents the median and the orange marker represents the mean.
Figure 5
Figure 5
Principal coordinate analysis (PCoA) plots of samples by Unifrac distances. (a) Plots of samples collected from both groups (healthy control (HC) and inflammatory bowel disease (IBD)) colored by disease state. Each marker denotes a single sample. PCoA plots for the unweighted Unifrac distances are shown on the left, whereas the PCoA plots for the weighted Unifrac distances are shown on the right. (b) Plots of all the samples colored by sample type (biopsy vs. fecal). (c) Plots of paired pre- and post-bowel prep samples of the HC group only. (d) Plots of paired pre- and post-bowel prep samples of IBD group only. The lines in (c and d) connect paired pre- and post-bowel prep samples of the same type for the same individual subject.
Figure 6
Figure 6
Relative abundances of significant indicator phyla and genera in healthy subjects. Each marker represents the relative (%) abundance of the denoted phylum or genus in one sample. Each panel shows the relative (%) abundance in the pre- and post-bowel prep samples collected from healthy subjects. The lines represent the mean relative abundance. Black circles represent the biopsy samples and blue circles represent the fecal samples. (a) The relative abundance of the phylum Actinobacteria (an indicator of the pre-bowel prep (BP) state). (b) The relative abundance of the phylum Bacteroidetes (an indicator of the post-BP state). (c) The relative abundance of the phylum Tenericutes (an indicator of the pre-BP state). (d) The relative abundance of an unclassified genus within the Bifidobacteriales order (an indicator of the pre-BP state). (e) The relative abundance of an unclassified genus within the Bacteroidaceae family (an indicator of the post-BP state). (f) The relative abundance of an unclassified genus within the Erysipelotrichaceae family (an indicator of the pre-BP state). (g) The relative abundance of an unclassified genus within the Clostridiales order (an indicator of the pre-BP state). ISA, indicator species analysis.
Figure 7
Figure 7
Relative abundances of significant indicator phyla and genera in inflammatory bowel disease (IBD) subjects. Each marker represents the relative (%) abundance of the denoted phylum or genus in one sample. Each panel shows the relative (%) abundance in the pre- and post-bowel prep samples collected from IBD subjects. The lines represent the mean relative abundance. Black circles represent the biopsy samples and blue circles represent the fecal samples. (a) The relative abundance of the phylum Bacteroidetes (the only significant indicator phylum in IBD subjects). (b) The relative abundance of an unclassified genus within the Bacteroidaceae family (an indicator of the post-bowel prep (BP) state). (c) The relative abundance of an unclassified genus within Streptococcaceae family (an indicator of the pre-BP state). (d) The relative abundance of an unclassified genus within the Clostridiaceae family (an indicator of the pre-BP state). ISA, indicator species analysis.
Figure 8
Figure 8
Bacteroidetes/Firmicutes ratio in samples. The bar chart shows the ratio of Bacteroidetes to Firmicutes in each sample group and compares the pre-bowel prep samples to the post-bowel prep samples. The average Bacteroidetes/Firmicutes ratio increased after bowel prep in the fecal samples of healthy control (HC) subjects (P=0.023, related-samples t-test). IBD, inflammatory bowel disease.

References

    1. Sears CL. A dynamic partnership: celebrating our gut flora. Anaerobe 2005; 11: 247–251.
    1. Eckburg PB, Bik EM, Bernstein CN et al. Diversity of the human intestinal microbial flora. Science (New York, NY) 2005; 308: 1635–1638.
    1. Backhed F, Ley RE, Sonnenburg JL et al. Host–bacterial mutualism in the human intestine. Science (New York, NY) 2005; 307: 1915–1920.
    1. Rajilic-Stojanovic M, Smidt H, de Vos WM. Diversity of the human gastrointestinal tract microbiota revisited. Environ Microbiol 2007; 9: 2125–2136.
    1. Rappe MS, Giovannoni SJ. The uncultured microbial majority. Annu Rev Microbiol 2003; 57: 369–394.
    1. Lewis K, Epstein S, D'Onofrio A et al. Uncultured microorganisms as a source of secondary metabolites. J Antibiot 2010; 63: 468–476.
    1. Turnbaugh PJ, Ley RE, Hamady M et al. The human microbiome project. Nature 2007; 449: 804–810.
    1. Dethlefsen L, Relman DA. Incomplete recovery and individualized responses of the human distal gut microbiota to repeated antibiotic perturbation. Proc Natl Acad Sci USA 2011; 108 (Suppl 1): 4554–4561.
    1. Croucher LJ, Bury JP, Williams EA et al. Commonly used bowel preparations have significant and different effects upon cell proliferation in the colon: a pilot study. BMC Gastroenterol 2008; 8: 54.
    1. Bucher P, Gervaz P, Egger JF et al. Morphologic alterations associated with mechanical bowel preparation before elective colorectal surgery: a randomized trial. Dis Colon Rectum 2006; 49: 109–112.
    1. Morotomi M, Guillem JG, Pocsidio J et al. Effect of polyethylene glycol-electrolyte lavage solution on intestinal microflora. Appl Environ Microbiol 1989; 55: 1026–1028.
    1. Goossens DAM, Jonkers DMAE, Russel GMVM et al. Bowel cleansing with subsequent intake of Lactobacillus plantarum 299v does not change the composition of the faecal flora. Microb Ecol Health Dis 2006; 18: 139–146.
    1. Mai V, Greenwald B, Morris JGJr et al. Effect of bowel preparation and colonoscopy on post-procedure intestinal microbiota composition. Gut 2006; 55: 1822–1823.
    1. Jain A, Smukalla SM, Ni J et al. Variable recovery of the major phyla of the gut microbiome after colonoscopy bowel preparation. Gastroenterology 2011; 140: S-362.
    1. O'Brien CL, Allison GE, Grimpen F et al. Impact of colonoscopy bowel preparation on intestinal microbiota. PLoS One 2013; 8: e62815.
    1. Bibiloni R, Tandon P, Vargas-Voracka F et al. Differential clustering of bowel biopsy-associated bacterial profiles of specimens collected in Mexico and Canada: what do these profiles represent? J Med Microbiol 2008; 57: 111–117.
    1. Momozawa Y, Deffontaine V, Louis E et al. Characterization of bacteria in biopsies of colon and stools by high throughput sequencing of the V2 region of bacterial 16S rRNA gene in human. PLoS One 2011; 6: e16952.
    1. Aronchick CA, Lipshutz WH, Wright SH et al. Validation of an instrument to assess colon cleansing. Am J Gastroenterol 1999;94:2667.
    1. Harvey RF, Bradshaw JM. A simple index of Crohn's-disease activity. Lancet 1980; 1: 514.
    1. Powell-Tuck J, Bown RL, Lennard-Jones JE. A comparison of oral prednisolone given as single or multiple daily doses for active proctocolitis. Scand J Gastroenterol 1978; 13: 833–837.
    1. Powell-Tuck J, Day DW, Buckell NA et al. Correlations between defined sigmoidoscopic appearances and other measures of disease activity in ulcerative colitis. Dig Dis Sci 1982; 27: 533–537.
    1. Satsangi J, Silverberg MS, Vermeire S et al. The Montreal classification of inflammatory bowel disease: controversies, consensus, and implications. Gut 2006; 55: 749–753.
    1. Silverberg MS, Satsangi J, Ahmad T et al. Toward an integrated clinical, molecular and serological classification of inflammatory bowel disease: report of a Working Party of the 2005 Montreal World Congress of Gastroenterology. Canad J Gastroenterol 2005;19 (Suppl A): 5A–36A.
    1. Smith DM, Snow DE, Rees E et al. Evaluation of the bacterial diversity of pressure ulcers using bTEFAP pyrosequencing. BMC Med Genom 2010; 3: 41.
    1. Mutlu EA, Keshavarzian A, Losurdo J et al. A compositional look at the human gastrointestinal microbiome and immune activation parameters in HIV infected subjects. PLoS Pathogens 2014; 10: e1003829.
    1. Caporaso JG, Kuczynski J, Stombaugh J et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 2010; 7: 335–336.
    1. Kuczynski J, Stombaugh J, Walters WA et al. Using QIIME to analyze 16S rRNA gene sequences from microbial communities. Current Protocols in Bioinformatics 2011; 36:10.7: 10.7.1–10.7.20.
    1. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010; 26: 2460–2461.
    1. Edgar RC, Haas BJ, Clemente JC et al. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 2011; 27: 2194–2200.
    1. Caporaso JG, Bittinger K, Bushman FD et al. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 2010; 26: 266–267.
    1. DeSantis TZ, Hugenholtz P, Larsen N et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol 2006; 72: 5069–5072.
    1. Wang Q, Garrity GM, Tiedje JM et al. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol 2007; 73: 5261–5267.
    1. Kovach W. MVSP (Multivariate Statistical Package), Version 3.13m. Wales, UK: Kovach Computing Service, 2004.
    1. Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 2005; 71: 8228–8235.
    1. Lozupone C, Hamady M, Knight R. UniFrac—an online tool for comparing microbial community diversity in a phylogenetic context. BMC Bioinform 2006; 7: 371.
    1. Lozupone C, Lladser ME, Knights D et al. UniFrac: an effective distance metric for microbial community comparison. ISME J 2011; 5: 169–172.
    1. Lozupone CA, Hamady M, Kelley ST et al. Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol 2007; 73: 1576–1585.
    1. Dufrene M, Legendre P. Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol Monogr 1997; 67: 345–366.
    1. Zoetendal EG, von Wright A, Vilpponen-Salmela T et al. Mucosa-associated bacteria in the human gastrointestinal tract are uniformly distributed along the colon and differ from the community recovered from feces. Appl Environ Microbiol 2002; 68: 3401–3407.
    1. Ley RE, Bäckhed F, Turnbaugh P et al. Obesity alters gut microbial ecology. Proc Natl Acad Sci USA 2005; 102: 11070–11075.
    1. Ley RE, Turnbaugh PJ, Klein S et al. Microbial ecology: human gut microbes associated with obesity. Nature 2006; 444: 1022–1023.
    1. Larsen N, Vogensen FK, van den Berg FW et al. Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PLoS One 2010; 5: e9085.
    1. Haahtela T, Holgate S, Pawankar R et al. The biodiversity hypothesis and allergic disease: world allergy organization position statement. World Allergy Organ J 2013; 6: 3.
    1. Levy S. Reduced bacterial biodiversity is associated with increased allergy. Environ Health Perspect 2012; 120: a304.
    1. Scher JU, Ubeda C, Artacho A et al. Decreased bacterial diversity characterizes the altered gut microbiota in patients with psoriatic arthritis, resembling dysbiosis in inflammatory bowel disease. Arthritis Rheumatol (Hoboken, NJ) 2015; 67: 128–139.
    1. Ott SJ, Musfeldt M, Wenderoth DF et al. Reduction in diversity of the colonic mucosa associated bacterial microflora in patients with active inflammatory bowel disease. Gut 2004; 53: 685–693.
    1. Ott SJ, Plamondon S, Hart A et al. Dynamics of the mucosa-associated flora in ulcerative colitis patients during remission and clinical relapse. J Clin Microbiol 2008; 46: 3510–3513.
    1. Martinez C, Antolin M, Santos J et al. Unstable composition of the fecal microbiota in ulcerative colitis during clinical remission. Am J Gastroenterol 2008; 103: 643–648.
    1. Dicksved J, Halfvarson J, Rosenquist M et al. Molecular analysis of the gut microbiota of identical twins with Crohn's disease. ISME J 2008; 2: 716–727.
    1. Fujimoto T, Imaeda H, Takahashi K et al. Decreased abundance of Faecalibacterium prausnitzii in the gut microbiota of Crohn's disease. J Gastroenterol Hepatol 2013; 28: 613–619.
    1. Walker AW, Sanderson JD, Churcher C et al. High-throughput clone library analysis of the mucosa-associated microbiota reveals dysbiosis and differences between inflamed and non-inflamed regions of the intestine in inflammatory bowel disease. BMC Microbiol 2011; 11: 7.
    1. Li Q, Wang C, Tang C et al. Molecular–phylogenetic characterization of the microbiota in ulcerated and non-ulcerated regions in the patients with Crohn's disease. PLoS One 2012; 7: e34939.
    1. Chen L, Wang W, Zhou R et al. Characteristics of fecal and mucosa-associated microbiota in Chinese patients with inflammatory bowel disease. Medicine (Baltimore, MD) 2014;93:e51.
    1. Hedin C, van der Gast CJ, Rogers GB et al. Siblings of patients with Crohn's disease exhibit a biologically relevant dysbiosis in mucosal microbial metacommunities. Gut 2015. doi:10.1136/gutjnl-2014-308896.
    1. Rehman A, Rausch P, Wang J et al. Geographical patterns of the standing and active human gut microbiome in health and IBD. Gut 2016; 65: 238–248.
    1. Menees S, Higgins P, Korsnes S et al. Does colonoscopy cause increased ulcerative colitis symptoms? Inflamm Bowel Dis 2007; 13: 12–18.
    1. Harrell L, Wang Y, Antonopoulos D et al. Standard colonic lavage alters the natural state of mucosal-associated microbiota in the human colon. PLoS One 2012;7:e32545.
    1. Jalanka J, Salonen A, Salojärvi J et al. Effects of bowel cleansing on the intestinal microbiota. Gut 2014; 64: 1562–1568.

Source: PubMed

3
Prenumerera