Association of gut microbiota with post-operative clinical course in Crohn’s disease
© Dey et al.; licensee BioMed Central Ltd. 2013
Received: 8 January 2013
Accepted: 6 August 2013
Published: 22 August 2013
The gut microbiome is altered in Crohn’s disease. Although individual taxa have been correlated with post-operative clinical course, global trends in microbial diversity have not been described in this context.
We collected mucosal biopsies from the terminal ileum and ascending colon during surgery and post-operative colonoscopy in 6 Crohn’s patients undergoing ileocolic resection (and 40 additional Crohn’s and healthy control patients undergoing either surgery or colonoscopy). Using next-generation sequencing technology, we profiled the gut microbiota in order to identify changes associated with remission or recurrence of inflammation.
We performed 16S ribosomal profiling using 101 base-pair single-end sequencing on the Illumina GAIIx platform with deep coverage, at an average depth of 1.3 million high quality reads per sample. At the time of surgery, Crohn’s patients who would remain in remission were more similar to controls and more species-rich than Crohn’s patients with subsequent recurrence. Patients remaining in remission also exhibited greater stability of the microbiota through time.
These observations permitted an association of gut microbial profiles with probability of recurrence in this limited single-center study. These results suggest that profiling the gut microbiota may be useful in guiding treatment of Crohn’s patients undergoing surgery.
KeywordsCrohn's disease Gut microbiome Next-generation sequencing Microbial profiling 16S rRNA gene
Crohn’s disease is an inflammatory bowel disease (IBD) with courses of recurrent inflammation, the etiology of which is likely rooted in the interplay of gut microbial imbalances (dysbiosis) and host immunity . Culture-independent studies of the gut microbiota in Crohn’s disease have documented that at least a subset of patients harbor microbial communities distinct from those in individuals without IBD [2–8]. The unpredictability of Crohn’s disease renders sampling of patients prior to the onset of inflammation inherently challenging. Here we addressed this by studying patients undergoing surgery and following them through the initial post-operative colonoscopies. Patients were thereby effectively “synchronized” in the history of their disease by starting at a common event. This design also permitted us to study the onset of inflammation after a state of known health (i.e., following the resection of all diseased bowel). Individual taxa, such as Faecalibacterium prausnitzii, have been associated with post-operative clinical course . This finding raises the question that we examined here: can profiling the gut microbiota at large be a useful clinical tool?
We characterized the gut microbiota in 6 patients with Crohn’s disease undergoing surgical resection of the terminal ileum to remove all evident inflammation, a procedure known as ileocolic resection. We collected samples during surgery and the first post-operative colonoscopy. We also characterized 5 Crohn’s patients undergoing surgery alone, 9 patients without IBD undergoing ileocolic resection (due to a colonic polyp or cancer), and 26 patients receiving only colonoscopy (15 with Crohn’s for the purpose of assessing disease activity as per the discretion of the treating gastroenterologist, 11 without IBD receiving colorectal cancer screening) (Additional file 1: Table S1). All patients were recruited via continuous enrollment from a single medical center in the San Francisco Bay Area.
Crohn’s patients undergoing ileocolic resection are at risk for post-operative recurrence of inflammation in the newly terminal ileum: 65% have visible inflammation at colonoscopy after 1 year . Powerful anti-inflammatory treatments can prevent or reduce this but are typically reserved for patients at high risk for recurrence [10, 11]. Gastroenterologists use clinical history and colonoscopy (initially performed approximately six months after surgery) to inform treatment decisions. At colonoscopy, mucosal health is graded with the Rutgeerts score (0–1 indicates remission; 2–4, recurrence) , which informs prognosis. The aim of the present study is to determine whether profiling the gut microbiota can further inform prognosis in this setting. We asked whether post-operative recurrence could be correlated with microbiota of the resected surgical tissue, either alone or in conjunction with specimens obtained subsequently. Mucosal biopsies were taken from the terminal ileum and ascending colon, distinguishing between inflamed and healthy-appearing mucosa. All biopsies from an individual patient at a given time point were obtained with the same surgical or endoscopic instruments, consistent with previously published protocols [3, 5, 13–15] but leading to the possibility of cross-contamination. Fecal samples were obtained one day prior to procedures when feasible. We performed microbial profiling of samples based on deep sequencing of a 16S ribosomal gene region. Based on observed associations between species abundance and post-operative recurrence, we propose that microbial profiling of patients with Crohn’s disease undergoing surgery may serve as a prognostic indicator for subsequent clinical course.
Patients scheduled for surgery and/or colonoscopy were recruited to the study from the University of California, San Francisco (UCSF), between November 2008 and August 2010. No patients underwent a procedure specifically for this study. Informed consent was obtained of each patient by one of the authors (ND). The use of human subjects was approved by the Committee on Human Research (CHR) at UCSF (study 10–01519). All samples analyzed were acquired specifically for this study; none were taken from a tissue bank.
Stool samples (250 mg) were homogenized in PowerBead tubes (MoBio, Carlsbad, CA) using a Mini-BeadBeater-8 (BioSpec) at the highest setting for 1 minute. DNA was extracted using the MoBio Power Soil DNA Kit (Carlsbad, CA) according to the manufacturer’s instructions. Frozen biopsies (30 mg) were suspended in 600 μl of Buffer RLT (Qiagen, Valencia, CA) that contained β-Mercaptoethanol (1 μl per 100 μl of buffer), and transferred to silica-bead-containing tubes (Lysing Matrix B, MP Biomedicals, Solon, OH). Tubes were placed in a Mini-BeadBeater-8 (BioSpec) and homogenized at the highest setting for 3 minutes. Samples were further thawed in a 37°C water bath for 5 minutes. DNA was subsequently extracted using the column-based Qiagen AllPrep kit (Valencia, CA) as per manufacturer instructions.
16S rRNA amplification
A fragment of the 16S ribosomal genes (containing the V1-V3 hypervariable regions) was amplified by PCR with 9F-containing (5′-CAAGCAGAAGACGGCATACGAGATNNNNNNGTGACTGGAGCGAGTTTGATC [AC]TGGCTCAG-3′) and 529R-containing (5′-AATGATCGGCGACCACCGAGATCTACACTCTTTCACCGCGGC[GT]GCTGGC-3′) primers at a final concentration of 400 nM using Phusion high-fidelity PCR Master Mix (Finnzymes). (Adapter sequences are underlined. Sample-specific 6-mer barcodes are represented by N’s. 9F and 529R sequences are italicized.) 10 barcodes were used, in varying combinations: AAGCTA, ACATCG, ATTGGC, CACTGT, CGTGAT, CTGATC, GATCTG, GCCTAA, TCAAGT, and TGGTCA. PCR conditions were as follows: 98°C initial denaturation for 30 seconds; 10 cycles of 98°C for 10 seconds, 53°C for 30 seconds, and 72°C for 30 seconds; a final extension at 72°C for 5 minutes. Following PCR, the sample was run on a 1.5% agarose gel (Invitrogen UltraPure agarose, Carlsbad, CA) at 80 Volts for 15 minutes. The amplified fragment was then excised and purified using GeneCatcher tips (the Gel Company, San Francisco, CA). Technical triplicates were performed for each sample.
Equimolar amounts of amplicons from 6–7 samples were pooled for multiplexing on each flow cell lane. The V3 hypervariable region was sequenced on an Illumina Genome Analyzer IIx using standard reagents. A custom sequencing primer (5′-GATCTACACTCTTTCACCGCGGC[GT]GCTGGC-3′), which contained 529R (italicized), was used to initiate sequencing beyond the 529R primer region, effectively increasing sequencing read length of the region of interest. Despite this strategy, conservation in 16S sequence beyond the 529R primer region required low cluster density on flow cells for successful sequencing. Matrix and phasing values generated from a library control lane were applied to the lanes containing the amplified 16S fragments. We generated 23.3 gigabases of 101 base pair single-end reads of partial 16S ribosomal sequence data.
16S rRNA data processing and analysis
Ribosomal sequence data was processed using tools from the QIIME pipeline v1.2.1 . Sequences were de-replicated at 100% sequence identity. No chimeric sequences were seen using ChimeraSlayer . Reads were retained if they were present in at least two samples and represented in at least ten times in the combined dataset. No read trimming or filtering was performed based on Q scores. Reads of presumed human DNA were identified and removed if they both bore 100% sequence identity to human DNA (reference human genome downloaded from http://www.ncbi.nlm.nih.gov/genome/ on May 19, 2011) and were not seen in the Greengenes ribosomal database (downloaded from http://greengenes.lbl.gov/Download/Sequence_Data/Greengenes_format/ on March 17, 2011). Contaminant human DNA constituted 15% of sequences. No contaminant fungal DNA was seen (18S fungal genes downloaded from NCBI on May 19, 2011). 18.4 gigabases of ribosomal sequence data were considered to be high quality and included in our analysis. This total represented a mean coverage of 1.3 million reads per sample ± 0.86 million [s.d.], with a range of 6,400-5,300,000 reads per sample. 14 of the total 139 samples had fewer than 250,000 reads (9 Crohn’s, 5 non-IBD). The mean number of reads per sample was not significantly different in Crohn’s compared to non-IBD controls. The range for samples from the six highlighted Crohn’s patients was 22,000-4,100,000 reads per sample. Reads were aligned using PyNAST . Phylogenetic trees were constructed using FastTree . The unique sequences were assigned to known taxa, to the extent possible, using the RDP classifier . The relative proportions of sequences represented by each taxa were averaged by patient per time point prior to group comparisons, which were made with the Student’s t-test (Welch two sample t-test) using values from a single time point per patient as indicated in the text.
Alpha diversity was measured using the Chao1 and Faith phylogenetic diversity (PD) metrics [21, 22]. Beta diversity was measured using the weighted UniFrac metric . UniFrac distance matrices generated via QIIME were used to compute average distances within and between various groups of samples (e.g. Crohn’s versus control biopsies). The arithmetic mean of coordinates in each dimension of the ordination space (i.e., the PCoA space, a multi-dimensional Euclidean space) was taken to determine the centroid. Autocorrelation, or pseudo-replication, was accounted for by taking the centroid in ordination space of all biopsies collected from a given individual at a given time. For example, in comparing Crohn’s patients to control patients, the centroids of biopsies collected at patients’ initial procedures within this study were used to determine average UniFrac distances. For comparisons that involved multiple pairwise UniFrac distances for an individual, such as inflamed versus healthy-appearing biopsies within a Crohn’s patient, the average was used for the overall calculation. Principal coordinates analysis (PCoA) plots were generated in R  using the ggplot2 package . The Student’s t-test with 1,000 Monte Carlo simulations was used to identify statistically significant differences between groups [26, 27].
Heat map generation showing likelihood of recurrence or remission
The two-sided Student’s t-test was used to compare sets of weighted UniFrac distances between surgical biopsies from Crohn’s patients and the following: surgical control biopsies, biopsies from the same patient at the time of the first post-operative colonoscopy, or all control biopsies (i.e., surgery and colonoscopy) (Additional file 2: Figure S4A). In addition, for each Crohn’s patient undergoing surgery, intra-individual comparisons of surgical biopsies were made (Additional file 2: Figure S4B). Likelihood heat maps were generated by overlaying distributions of weighted UniFrac distances in remission and recurrence, then calculating the proportion of UniFrac distances (in bins of width 0.1) seen in remission or recurrence. Leave-one-out cross-validation analyses were performed by generating a predictive model using the aggregate of UniFrac distances from all but one patient, and then “predicting” the clinical outcome of the excluded patient. We generated these predictions by computing the average probability of recurrence arising from each UniFrac distance. Statistical calculations were performed in R . Heat maps were generated using the “heatmap.2” function in the gplots package .
Deep sequencing of a 16S ribosomal RNA gene region
We amplified 16S ribosomal genes with 9F and 529R primers. We used a 529R-containing custom sequencing primer to initiate sequencing beyond the conserved primer sequence region and generate 101 base-pair single-end reads of the V3 hypervariable region using the Illumina GAIIx platform. A total of 139 samples were sequenced from 46 subjects. After filtering, our dataset was comprised of 18.4 gigabases of high-quality reads (mean = 1.3 million reads per sample ± 0.86 million [s.d.]; range 6,400-5,300,000 reads) (see Methods for additional details). The combined dataset contained 300,389 unique sequences, or 57,603 operational taxonomic units (OTUs) when clustering by 97% nucleotide sequence identity. This is a higher number of OTUs than previously reported to be associated with human gut [26, 29–32], but we suspected this to be a consequence of sequencing error combined with the deep level of coverage . The region of 16S selected evolves at a rate similar to the 16S gene overall , so this is not likely to be a major source of difference. We confirmed this by simulating deep sequencing of the gut microbiota. Using previously published datasets of 9,920  and 7,208  near-full-length 16S sequences, we generated large datasets in silico of 200 million reads (comparable to our combined dataset of 180 million reads), adding mismatches at rates of 0.01%, 0.1%, or 0.5% per base. Where the original samples contained 451 and 204 OTUs, respectively, the simulated deeply sequenced samples with a 0.1% error rate contained 32,868 and 26,622 OTUs, confirming the prior observation that sequencing noise on the Illumina platform can inflate OTU estimates . With an error rate of 0.5% per base, we observed 554,073 and 431,544 OTUs, respectively, suggesting that – even if the 57,603 OTUs we observed do not represent novel taxa revealed by deep sequencing but instead can largely be accounted for by sequencing error – the true sequencing error rate in our dataset was closer to 0.1% per base. We also performed the inverse experiment, subsampling sets of 10,000 sequences from each of the gut samples we collected, in order to estimate OTU counts that might have been seen with shallower sequencing. These subsamples contained 57–245 (95% confidence interval) OTUs, similar to OTU counts at this sequencing depth seen in prior studies [26, 29].
Inter-individual heterogeneity in Crohn’s disease
Still, several significant differences existed. The two most common phyla in both Crohn’s and control samples were Bacteroidetes and Firmicutes (Figures 1C-D). All biopsies considered, Crohn’s patients had greater relative abundance of Proteobacteria (P<10-4) and Fusobacteria (P<0.05), and lower relative abundance of Bacteroidetes (P<10-4), compared to non-IBD patients (Additional file 1: Table S3). Similar to several previous reports [9, 36, 37], our data showed greater relative abundance of Faecalibacterium in healthy patients compared to Crohn’s patients, and Crohn’s remission patients compared to Crohn’s recurrence patients; however, these differences were non-significant. The role of F. prausnitzii in the pathogenesis of Crohn’s disease is yet unclear . Consistent with previous observations, the family Enterobacteriaceae was more abundant in Crohn’s patients in our study (P<10-3) [4, 7, 15]. The genus Ruminococcus was less abundant in surgical biopsies from Crohn’s patients relative to non-IBD surgical controls (P<0.05), in line with a study reporting lower levels of Ruminococcus gnavus in Crohn’s disease . In the context of mixed reports [9, 36, 40–42], our data showed lower abundance of Bacteroides in Crohn’s disease patients (P<10-2). The presence of discrepancies between studies may represent species- or strain-level differences, or be a consequence of the inter-individual heterogeneity observed here and elsewhere [2, 15, 43].
Dysbiosis was neither localized nor specific to inflamed mucosa. The average UniFrac distance between inflamed biopsies and adjacent healthy-appearing biopsies from the same individual taken at the same time (“inflamed vs healthy-appearing biopsies” in Additional file 2: Figure S1) was significantly less than the distance between inflamed biopsies from different individuals (“inflamed vs inflamed [between persons]”) (P<0.005), indicating that measured microbial community composition is fairly consistent around the ileocolic junction within a patient, regardless of local inflammation. As noted above, this may be due in part to cross-contamination. This held true for both comparisons within the same region (ileum or colon) or between regions (ileum versus colon) (Additional file 2: Figure S1). Prior studies have reported mixed results to this end [5, 15, 44, 45]. These data also suggest that inflamed biopsies in different patients vary in microbial composition, as supported by the comparison of the average distance between inflamed biopsies from different individuals (“inflamed vs inflamed [between persons]” in Figure 1B) with the average distance between inflamed biopsies from an individual taken at a given time point (“inflamed vs inflamed [within individual]”) (P<0.05). Indeed, inflamed biopsies from different patients did not cluster together, but rather varied across PCoA space (Figure 1A), suggesting multiple varieties of dysbiosis. Finally, stool and biopsies acquired concomitantly from a single individual were no more similar to each other than biopsies (or stool) from different individuals (Additional file 2: Figures S1, S2).
Surgical biopsies from Crohn’s patients who remain in remission are more similar to controls
Given these patterns, we asked whether a single surgical biopsy from a Crohn’s patient could be used to estimate probability of recurrence or remission. We took two approaches to studying this. First, we randomly selected a single surgical biopsy from each of the six Crohn’s patients and determined the UniFrac distance to a randomly selected control biopsy. We then compared the distributions of three distances from remission patients versus three distances from recurrence patients. A different control biopsy was used for each comparison, in order to test whether any observed differences were robust to variation among controls. Repeating this exercise 10,000 times, we found that these distributions were significantly different in 34% of cases, as defined by a P value of less than 0.05 (Figure 2D). Second, we determined UniFrac distances from each Crohn’s surgical biopsy to the centroid in ordination space of all control biopsies. We compared distributions in remission versus recurrence, examining all 486 combinations of pairwise comparisons between each surgical sample from the six Crohn’s patients and the centroid of control biopsies (Figure 2E). The distributions of distances in remission and recurrence were significantly different in 56% of comparisons. Both of these calculations suggest that use of a single biopsy specimen for prognostic purposes may not be reliably useful.
We therefore examined the utility of combining information from multiple biopsies. For each of the six Crohn’s patients, we estimated outcome based on likelihood models constructed from the other five patients (Additional file 2: Figures S5, S6). In this leave-one-patient-out cross-validation analysis, pairwise UniFrac distances between surgical biopsies of the five-patient “training set” and all controls were used to generate heat maps reflecting probabilities of recurrence given a single UniFrac distance. We then computed the probability of recurrence by averaging the probabilities generated by each surgical biopsy from the “test patient.” This scheme correctly predicted recurrence or remission in all cases. In two cases, the correct prediction was made with low confidence, based on a probability of 53-55%; in the remaining cases, the signal was substantially stronger, such that the mean probability associated with all predictions in this leave-one-patient-out analysis was 79% ± 9%. Multiple biopsies appear to improve signal despite the general trend that biopsies from a patient are most similar to each other. This may relate to the inconsistency seen between biopsies from surgery in patients with recurrence, discussed below.
Crohn’s patients with post-operative recurrence have less consistency between surgical biopsies and less stability through time
Reduced species richness in Crohn’s pronounced in recurrence
Consistent with prior reports [2, 4, 15, 46], we observed reduced alpha diversity (i.e., species richness) in the gut microbiota of Crohn’s patients. In light of variations in sequencing depth, we subsampled sets of 10,000 sequences from each sample, excluding the single sample with fewer total reads (a colonic biopsy from a patient without IBD). We observed 143 +/− 5 (s.e.m.) OTUs in patients with Crohn’s disease, as compared to 176 +/− 6 (s.e.m.) OTUs in non-IBD controls. The Chao1 nonparametric estimator and the Faith phylogenetic diversity (PD) metric also indicated reduced diversity in Crohn's patients compared to controls (Additional file 2: Figures S9, S10) [21, 22]. The reduction in alpha diversity in Crohn’s was significant in patients with recurrence with both Chao1 and PD metrics relative to non-IBD controls (P<0.05). These data are consistent with the hypothesis that a rich gut microbiota may protect against disease.
We identified a limited number of significant taxonomic differences in samples from Crohn’s patients with recurrence relative to those who stayed in remission. At surgery, patients who would go on to experience recurrence had a significantly lower relative abundance of the families Lachnospiraceae (P<10-2) and Erysipelotrichaceae (P<0.05), as well as an unidentified genus within Clostridia (P<10-2), all members of the phylum Firmicutes (Additional file 1: Table S3). Recurrence was associated with a greater abundance of Rhodobacteraceae (of class Alphaproteobacteria; P<0.05) and an unknown Proteobacteria (P<10-2). At post-operative colonoscopy, the genus Rhizobium (also Alphaproteobacteria) was the lone significant difference between Crohn’s patients experiencing recurrence and those remaining in remission, more abundant in patients with recurrence (P<0.05).
Understanding the gut microbiome in IBD is of interest to both clinicians and patients, who are eager to predict clinical course of the disease and identify appropriate treatments to improve health outcomes, including probiotics and diet . Studying the etiologic role of the gut microbiome in Crohn’s disease is challenging, due to multiple issues including unpredictability of disease course and patient selection. Many studies have characterized patients at a single time point and attributed to Crohn’s disease the observed microbial differences from controls [2, 4, 7, 15, 36, 37, 42–46]. Global trends in diversity have been underexplored. Here, we followed carefully phenotyped patients over time to correlate microbial profiles with clinical course, synchronizing the time course with ileocolic resection, and controlling for biologic medication usage and geographic impact. Consistent with prior reports suggesting translational utility of gut microbial profiling in Crohn’s disease [9, 48], we show here that patients experiencing recurrence of inflammation after surgical resection of all diseased bowel have distinct diversity profiles relative to patients who remain in remission. A clinical test characterizing multiple biopsies from a surgically resected specimen may ultimately be useful for informing prognosis.
While this study addresses several challenges of study design pertaining to variability and unpredictability in disease course, it is also subject to certain limitations. The findings presented here represent a small number of patients on biologic therapies at a single medical center. Conceivably, these findings may therefore only apply to a limited subset of Crohn’s patients, whether it be patients on related medications or those in the study’s catchment population. Given differences in study design, patient characteristics, sequencing strategies, and other experimental variables, it is difficult to directly evaluate patients described in previous studies in the context of our findings. The development of a practical test based on these findings will require validation on a larger and more diverse patient population involving multiple medical centers.
There were several questions we were unable to address. We wished to collect fecal and biopsy samples at all time points in order to examine the relationships of both to outcomes. Unfortunately, we were only able collect fecal samples from already-enrolled Crohn’s patients concurrent with their post-operative colonoscopy. Subsequent investigations examining the predictive capacity of stool may build upon current knowledge of fecal biomarkers , perhaps by focusing on patients with inflammation of the distal gut or in patients who are not otherwise scheduled for an invasive evaluation (e.g. colonoscopy).
We also could not interpret the observed post-operative “migration” of samples on the PCoA plot towards the healthy controls in cases of both recurrence and remission. It is unclear whether this constituted some degree of recovery from dysbiosis. Pre-operative data on the Crohn’s patients or longitudinal data on the non-IBD surgical controls may have enabled a better understanding of this migration. Future studies should evaluate pharmacologic or dietary interventions in the context of microbial shifts through time. Recent reviews of probiotics and prebiotics in Crohn’s disease have underscored a lack of knowledge in this area [50–53].
In summary, we demonstrated in this limited cohort at a single medical center that analyses of surgical biopsies using the UniFrac metric are associated with subsequent clinical course. Whether the severity of dysbiosis or the instability of the gut microbiota (or both) is the critical determinant of recurrence in Crohn’s patients undergoing surgery will need to be further explored. Whether these influence inflammation or are the result thereof remains unclear. Dysbiosis in Crohn’s has many varieties and does not appear to be specific to inflamed mucosa. With apologies to Tolstoy, “unhappy guts” in Crohn’s appear to be “unhappy” in different ways . Taken together with the broad spectrum of genetic mutations associated with Crohn’s disease [55, 56], it is not surprising that a specific cause has not been identified, and that IBD patients require individualized care. Further studies are warranted to validate these conclusions in a larger set of patients and other clinical contexts.
Inflammatory bowel disease
Operational taxonomic unit
Principal coordinates analysis.
The authors wish to thank all the patients who generously agreed to participate in this study. We are grateful to Edward Kim, Averil Ma, Uma Mahadevan, Jonathan Terdiman, Madhulika Varma, and Fernando Velayos for facilitating patient recruitment and biopsy acquisition.
We have obtained a letter of certification to submit the combined 16S rRNA sequence dataset to the NCBI database of Genotypes and Phenotypes (dbGaP). Data are being submitted at the time of publication.
National Academies Keck Futures Initiative (Geno11), National Institutes of Health (T32 Institutional Training grant DK00700733; HG4872), and Marie Curie International Outgoing Fellowship (PIOF-GA-2009-237751).
- Abraham C, Cho JH: Inflammatory bowel disease. N Engl J Med. 2009, 361: 2066-2078. 10.1056/NEJMra0804647.View ArticlePubMedPubMed CentralGoogle Scholar
- 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-13785. 10.1073/pnas.0706625104.View ArticlePubMedPubMed CentralGoogle Scholar
- Neut C, Bulois P, Desreumaux P, Membré J-M, Lederman E, Gambiez L, Cortot A, Quandalle P, van Kruiningen H, Colombel J-F: Changes in the bacterial flora of the neoterminal ileum after ileocolonic resection for Crohn’s disease. Am J Gastroenterol. 2002, 97: 939-946. 10.1111/j.1572-0241.2002.05613.x.View ArticlePubMedGoogle Scholar
- Dicksved J, Halfvarson J, Rosenquist M, Järnerot G, Tysk C, Apajalahti J, Engstrand L, Jansson JK: Molecular analysis of the gut microbiota of identical twins with Crohn’s disease. ISME J. 2008, 2: 716-727. 10.1038/ismej.2008.37.View ArticlePubMedGoogle Scholar
- Bibiloni R, Mangold M, Madsen KL, Fedorak RN, Tannock GW: The bacteriology of biopsies differs between newly diagnosed, untreated, Crohn’s disease and ulcerative colitis patients. J Med Microbiol. 2006, 55 (Pt 8): 1141-1149.View ArticlePubMedGoogle Scholar
- Swidsinski A, Weber J, Loening-Baucke V, Hale LP, Lochs H: Spatial organization and composition of the mucosal flora in patients with inflammatory bowel disease. J Clin Microbiol. 2005, 43: 3380-3389. 10.1128/JCM.43.7.3380-3389.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Willing BP, Dicksved J, Halfvarson J, Andersson AF, Lucio M, Zheng Z, Järnerot G, Tysk C, Jansson JK, Engstrand L: A pyrosequencing study in twins shows that gastrointestinal microbial profiles vary with inflammatory bowel disease phenotypes. Gastroenterology. 2010, 139: 1844-1854. 10.1053/j.gastro.2010.08.049.View ArticlePubMedGoogle Scholar
- Frank DN, Robertson CE, Hamm CM, Kpadeh Z, Zhang T, Chen H, Zhu W, Sartor RB, Boedeker EC, Harpaz N, Pace NR, Li E: Disease phenotype and genotype are associated with shifts in intestinal‒associated microbiota in inflammatory bowel diseases. Inflamm Bowel Dis. 2011, 17: 179-184. 10.1002/ibd.21339.View ArticlePubMedGoogle Scholar
- Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermúdez-Humarán LG, Gratadoux J-J, Blugeon S, Bridonneau C, Furet J-P, Corthier G, Grangette C, Vasquez N, Pochart P, Trugnan G, Thomas G, Blottière HM, Doré J, Marteau P, Seksik P, Langella P: Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci USA. 2008, 105: 16731-16736. 10.1073/pnas.0804812105.View ArticlePubMedPubMed CentralGoogle Scholar
- Terdiman JP: Prevention of postoperative recurrence in Crohn’s disease. Clin Gastroenterol Hepatol. 2008, 6: 616-620. 10.1016/j.cgh.2007.08.020.View ArticlePubMedGoogle Scholar
- Regueiro M, Schraut W, Baidoo L, Kip KE, Sepulveda AR, Pesci M, Harrison J, Plevy SE: Infliximab prevents Crohn’s disease recurrence after ileal resection. Gastroenterology. 2009, 136: 441-450. 10.1053/j.gastro.2008.10.051.View ArticlePubMedGoogle Scholar
- Rutgeerts P, Geboes K, Vantrappen G, Kerremans R, Coenegrachts JL, Coremans G: Natural history of recurrent Crohn’s disease at the ileocolonic anastomosis after curative surgery. Gut. 1984, 25: 665-672. 10.1136/gut.25.6.665.View ArticlePubMedPubMed CentralGoogle Scholar
- Dave M, Johnson LA, Walk ST, Young VB, Stidham RW, Chaudhary MN, Funnell J, Higgins PDR: A randomised trial of sheathed versus standard forceps for obtaining uncontaminated biopsy specimens of microbiota from the terminal ileum. Gut. 2011, 60: 1043-1049. 10.1136/gut.2010.224337.View ArticlePubMedPubMed CentralGoogle Scholar
- Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, Gill SR, Nelson KE, Relman DA: Diversity of the human intestinal microbial flora. Science. 2005, 308: 1635-1638. 10.1126/science.1110591.View ArticlePubMedPubMed CentralGoogle Scholar
- Walker AW, Sanderson JD, Churcher C, Parkes GC, Hudspith BN, Rayment N, Brostoff J, Parkhill J, Dougan G, Petrovska L: 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-10.1186/1471-2180-11-7.View ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7: 335-336. 10.1038/nmeth.f.303.View ArticlePubMedPubMed CentralGoogle Scholar
- Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E, Methé B, DeSantis TZ, Petrosino JF, Knight R, Birren BW: Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011, 21: 494-504. 10.1101/gr.112730.110.View ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R: PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010, 26: 266-267. 10.1093/bioinformatics/btp636.View ArticlePubMedGoogle Scholar
- Price MN, Dehal PS, Arkin AP: FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009, 26: 1641-1650. 10.1093/molbev/msp077.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Q, Garrity GM, Tiedje JM, Cole JR: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007, 73: 5261-5267. 10.1128/AEM.00062-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Chao A: Non-parametric estimation of the number of classes in a population. Scand J Stat. 1987, 11: 265-270.Google Scholar
- Faith DP: Phylogenetic pattern and the quantification of organismal biodiversity. Philos Trans R Soc Lond B Biol Sci. 1994, 345: 45-58. 10.1098/rstb.1994.0085.View ArticlePubMedGoogle Scholar
- Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005, 71: 8228-8235. 10.1128/AEM.71.12.8228-8235.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2011, Vienna: R Foundation for Statistical ComputingGoogle Scholar
- Wickham H: Ggplot2: Elegant Graphics for Data Analysis. 2009, Springer New York, http://cran.r-project.org/web/packages/ggplot2/,View ArticleGoogle Scholar
- Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, Gordon JI: A core gut microbiome in obese and lean twins. Nature. 2009, 457: 480-484. 10.1038/nature07540.View ArticlePubMedGoogle Scholar
- Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP, Heath AC, Warner B, Reeder J, Kuczynski J, Caporaso JG, Lozupone CA, Lauber C, Clemente JC, Knights D, Knight R, Gordon JI: Human gut microbiome viewed across age and geography. Nature. 2012, 486: 222-227.PubMedPubMed CentralGoogle Scholar
- Warnes GR, Bolker B, Bonebakker L, Gentleman R, Liaw WHA, Lumley T, Maechler M, Magnusson A, Moeller S, Schwartz M, Venables B: gplots: Various R programming tools for plotting data. 2010Google Scholar
- Dethlefsen L, Huse S, Sogin ML, Relman DA: The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing. PLoS Biol. 2008, 6: e280-10.1371/journal.pbio.0060280.View ArticlePubMedPubMed CentralGoogle Scholar
- Ley RE, Turnbaugh PJ, Klein S, Gordon JI: Microbial ecology: human gut microbes associated with obesity. Nature. 2006, 444: 1022-1023. 10.1038/4441022a.View ArticlePubMedGoogle Scholar
- Bäckhed F, Ley RE, Sonnenburg JL, Peterson DA, Gordon JI: Host-bacterial mutualism in the human intestine. Science. 2005, 307: 1915-1920. 10.1126/science.1104816.View ArticlePubMedGoogle Scholar
- Muegge BD, Kuczynski J, Knights D, Clemente JC, González A, Fontana L, Henrissat B, Knight R, Gordon JI: Diet Drives Convergence in Gut Microbiome Functions Across Mammalian Phylogeny and Within Humans. Science. 2011, 332: 970-974. 10.1126/science.1198719.View ArticlePubMedPubMed CentralGoogle Scholar
- Huse SM, Welch DM, Morrison HG, Sogin ML: Ironing out the wrinkles in the rare biosphere through improved OTU clustering. Environ Microbiol. 2010, 12: 1889-1898. 10.1111/j.1462-2920.2010.02193.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Soergel DAW: Computational Methods for Evaluating Microbial Diversity (Ph.D. dissertation, University of California, Berkeley). 2010, Berkeley: Ph.D. dissertation. University of CaliforniaGoogle Scholar
- Degnan PH, Ochman H: Illumina-based analysis of microbial community diversity. ISME J. 2011, 6: 183-194.View ArticlePubMedPubMed CentralGoogle Scholar
- Kang S, Denman SE, Morrison M, Yu Z, Dore J, Leclerc M, McSweeney CS: Dysbiosis of fecal microbiota in Crohn’s disease patients as revealed by a custom phylogenetic microarray. Inflamm Bowel Dis. 2010, 16: 2034-2042. 10.1002/ibd.21319.View ArticlePubMedGoogle Scholar
- Mondot S, Kang S, Furet JP, Aguirre de Carcer D, McSweeney C, Morrison M, Marteau P, Doré J, Leclerc M: Highlighting new phylogenetic specificities of Crohn’s disease microbiota. Inflammatory Bowel Diseases. 2011, 17: 185-192. 10.1002/ibd.21436.View ArticlePubMedGoogle Scholar
- Jia W, Whitehead RN, Griffiths L, Dawson C, Waring RH, Ramsden DB, Hunter JO, Cole JA: Is the abundance of Faecalibacterium prausnitzii relevant to Crohn’s disease?. FEMS Microbiol Lett. 2010, 310: 138-144. 10.1111/j.1574-6968.2010.02057.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Joossens M, Huys G, Cnockaert M, De Preter V, Verbeke K, Rutgeerts P, Vandamme P, Vermeire S: Dysbiosis of the faecal microbiota in patients with Crohn’s disease and their unaffected relatives. Gut. 2011, 60: 631-637. 10.1136/gut.2010.223263.View ArticlePubMedGoogle Scholar
- Linskens RK, Huijsdens XW, Savelkoul PH, Vandenbroucke-Grauls CM, Meuwissen SG: The bacterial flora in inflammatory bowel disease: current insights in pathogenesis and the influence of antibiotics and probiotics. Scand J Gastroenterol Suppl. 2001, 234: 29-40.View ArticlePubMedGoogle Scholar
- Andoh A, Tsujikawa T, Sasaki M, Mitsuyama K, Suzuki Y, Matsui T, Matsumoto T, Benno Y, Fujiyama Y: Faecal microbiota profile of Crohn’s disease determined by terminal restriction fragment length polymorphism analysis. Aliment Pharmacol Ther. 2009, 29: 75-82. 10.1111/j.1365-2036.2008.03860.x.View ArticlePubMedGoogle Scholar
- Ott SJ, Musfeldt M, Wenderoth DF, Hampe J, Brant O, Fölsch UR, Timmis KN, Schreiber S: Reduction in diversity of the colonic mucosa associated bacterial microflora in patients with active inflammatory bowel disease. Gut. 2004, 53: 685-693. 10.1136/gut.2003.025403.View ArticlePubMedPubMed CentralGoogle Scholar
- Martinez-Medina M, Aldeguer X, Gonzalez-Huix F, Acero D, Garcia-Gil LJ: Abnormal microbiota composition in the ileocolonic mucosa of Crohn’s disease patients as revealed by polymerase chain reaction-denaturing gradient gel electrophoresis. Inflamm Bowel Dis. 2006, 12: 1136-1145. 10.1097/01.mib.0000235828.09305.0c.View ArticlePubMedGoogle Scholar
- Gophna U, Sommerfeld K, Gophna S, Doolittle WF, Veldhuyzen van Zanten SJO: Differences between tissue-associated intestinal microfloras of patients with Crohn’s disease and ulcerative colitis. J Clin Microbiol. 2006, 44: 4136-4141. 10.1128/JCM.01004-06.View ArticlePubMedPubMed CentralGoogle Scholar
- Keighley MR, Arabi Y, Dimock F, Burdon DW, Allan RN, Alexander-Williams J: Influence of inflammatory bowel disease on intestinal microflora. Gut. 1978, 19: 1099-1104. 10.1136/gut.19.12.1099.View ArticlePubMedPubMed CentralGoogle Scholar
- Manichanh C, Rigottier-Gois L, Bonnaud E, Gloux K, Pelletier E, Frangeul L, Nalin R, Jarrin C, Chardon P, Marteau P, Roca J, Dore J: Reduced diversity of faecal microbiota in Crohn’s disease revealed by a metagenomic approach. Gut. 2006, 55: 205-211. 10.1136/gut.2005.073817.View ArticlePubMedPubMed CentralGoogle Scholar
- Hedin CRH, Mullard M, Sharratt E, Jansen C, Sanderson JD, Shirlaw P, Howe LC, Djemal S, Stagg AJ, Lindsay JO, Whelan K: Probiotic and prebiotic use in patients with inflammatory bowel disease: a case–control study. Inflamm Bowel Dis. 2010, 16: 2099-2108. 10.1002/ibd.21286.View ArticlePubMedGoogle Scholar
- Scanlan PD, Shanahan F, O’Mahony C, Marchesi JR: Culture-independent analyses of temporal variation of the dominant fecal microbiota and targeted bacterial subgroups in Crohn’s disease. J Clin Microbiol. 2006, 44: 3980-3988. 10.1128/JCM.00312-06.View ArticlePubMedPubMed CentralGoogle Scholar
- Judd TA, Day AS, Lemberg DA, Turner D, Leach ST: Update of fecal markers of inflammation in inflammatory bowel disease. J Gastroenterol Hepatol. 2011, 26: 1493-1499. 10.1111/j.1440-1746.2011.06846.x.View ArticlePubMedGoogle Scholar
- Butterworth AD, Thomas AG, Akobeng AK: Probiotics for induction of remission in Crohn’s disease. Cochrane Database Syst Rev. 2008, 3: CD006634-PubMedGoogle Scholar
- Doherty G, Bennett G, Patil S, Cheifetz A, Moss AC: Interventions for prevention of post-operative recurrence of Crohn’s disease. Cochrane Database Syst Rev. 2009, 4: CD006873-PubMedGoogle Scholar
- Rolfe VE, Fortun PJ, Hawkey CJ, Bath-Hextall F: Probiotics for maintenance of remission in Crohn’s disease. Cochrane Database Syst Rev. 2006, 4: CD004826-PubMedGoogle Scholar
- Benjamin JL, Hedin CRH, Koutsoumpas A, Ng SC, McCarthy NE, Hart AL, Kamm MA, Sanderson JD, Knight SC, Forbes A, Stagg AJ, Whelan K, Lindsay JO: Randomised, double-blind, placebo-controlled trial of fructo-oligosaccharides in active Crohn’s disease. Gut. 2011, 60: 923-929. 10.1136/gut.2010.232025.View ArticlePubMedGoogle Scholar
- Tolstoy L: Anna Karenina. 2004, Penguin ClassicsGoogle Scholar
- Barrett JC, Hansoul S, Nicolae DL, Cho JH, Duerr RH, Rioux JD, Brant SR, Silverberg MS, Taylor KD, Barmada MM, Bitton A, Dassopoulos T, Datta LW, Green T, Griffiths AM, Kistner EO, Murtha MT, Regueiro MD, Rotter JI, Schumm LP, Steinhart AH, Targan SR, Xavier RJ, Libioulle C, Sandor C, Lathrop M, Belaiche J, Dewit O, Gut I, Heath S, et al: Genome-wide association defines more than 30 distinct susceptibility loci for Crohn’s disease. Nat Genet. 2008, 40: 955-962. 10.1038/ng.175.View ArticlePubMedPubMed CentralGoogle Scholar
- Franke A, McGovern DPB, Barrett JC, Wang K, Radford-Smith GL, Ahmad T, Lees CW, Balschun T, Lee J, Roberts R, Anderson CA, Bis JC, Bumpstead S, Ellinghaus D, Festen EM, Georges M, Green T, Haritunians T, Jostins L, Latiano A, Mathew CG, Montgomery GW, Prescott NJ, Raychaudhuri S, Rotter JI, Schumm P, Sharma Y, Simms LA, Taylor KD, Whiteman D, et al: Genome-wide meta-analysis increases to 71 the number of confirmed Crohn’s disease susceptibility loci. Nat Genet. 2010, 42: 1118-1125. 10.1038/ng.717.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-230X/13/131/prepub