Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment
- Xochitl C Morgan†1,
- Timothy L Tickle†1, 2,
- Harry Sokol†3, 4, 5,
- Dirk Gevers2,
- Kathryn L Devaney3, 4,
- Doyle V Ward2,
- Joshua A Reyes1,
- Samir A Shah6,
- Neal LeLeiko6,
- Scott B Snapper7,
- Athos Bousvaros7,
- Joshua Korzenik3, 7,
- Bruce E Sands8,
- Ramnik J Xavier2, 3, 4 and
- Curtis Huttenhower1, 2Email author
© Morgan et al.; licensee BioMed Central Ltd. 2012
Received: 16 April 2012
Accepted: 26 September 2012
Published: 26 September 2012
The inflammatory bowel diseases (IBD) Crohn's disease and ulcerative colitis result from alterations in intestinal microbes and the immune system. However, the precise dysfunctions of microbial metabolism in the gastrointestinal microbiome during IBD remain unclear. We analyzed the microbiota of intestinal biopsies and stool samples from 231 IBD and healthy subjects by 16S gene pyrosequencing and followed up a subset using shotgun metagenomics. Gene and pathway composition were assessed, based on 16S data from phylogenetically-related reference genomes, and associated using sparse multivariate linear modeling with medications, environmental factors, and IBD status.
Firmicutes and Enterobacteriaceae abundances were associated with disease status as expected, but also with treatment and subject characteristics. Microbial function, though, was more consistently perturbed than composition, with 12% of analyzed pathways changed compared with 2% of genera. We identified major shifts in oxidative stress pathways, as well as decreased carbohydrate metabolism and amino acid biosynthesis in favor of nutrient transport and uptake. The microbiome of ileal Crohn's disease was notable for increases in virulence and secretion pathways.
This inferred functional metagenomic information provides the first insights into community-wide microbial processes and pathways that underpin IBD pathogenesis.
Inflammatory bowel disease (IBD), a chronic and relapsing inflammatory condition of the gastrointestinal (GI) tract, is intimately linked to the microbial communities of the human gut. Although it is now widely accepted that IBD results from altered interactions between gut microbes and the intestinal immune system, the precise nature of the intestinal microbiota dysfunction in IBD remains to be elucidated . IBD further includes two main subtypes, ulcerative colitis (UC) and Crohn's disease (CD), which each include distinct microbial perturbations and tissue localizations. The former is confined to the colon, while the latter may affect any part of the digestive tract, with unclear implications for microbial involvement or causality . In particular, the microbial mechanisms and metabolism underlying the role of the GI microbiome in IBD onset and its alteration in the course of active treatment and recovery are still unknown.
In the last decade, advances in DNA sequencing have allowed exploration of the 40% of the gut microbiome that is still uncultured , setting the stage for investigation of the IBD microbiome. The GI microbiome of healthy humans is dominated by four major bacterial phyla: Firmicutes, Bacteroidetes, and to a lesser degree Proteobacteria and Actinobacteria [4, 5]. Many studies have observed imbalances or dysbioses in the GI microbiomes of IBD patients [6–13]; in both CD and UC patients, there is decreased biodiversity, a lower proportion of Firmicutes, and an increase in Gammaproteobacteria . In CD, proportions of the Clostridia are altered: the Roseburia and Faecalibacterium genera of the Lachnospiracae and Ruminococcaceae families are decreased, whereas Ruminococcus gnavus increases [15–17]. Specific features of UC-associated dysbiosis are less described, although increased sulfate-reducing Deltaproteobacteria have been reported [18, 19]. These studies have described typical changes in composition of the IBD gut community, but the functional roles of these organisms - or of the entirety of a dysbiotic community - remain less clear.
The normal gut microbiome exhibits tremendous functional diversity encoded by a collection of bacterial genes numbering more than 100 times the human gene set [4, 20]. Thus, the genomic potential of the human microbiome is far greater than that of its host, and treatments, diets, or medications that affect the host will also likely affect the microbiome. A primary example of the importance of the microbiome to host health is in the digestion of dietary fiber, which is used by the microbiota of the lower GI tract as their main source of energy . Fibrolytic bacteria degrade polysaccharides into smaller carbohydrates, which are then fermented into short-chain fatty acids (SCFAs) such as acetate, propionate, and butyrate. Butyrate in particular is a major source of energy for colonocytes, but all three of these have demonstrated immunomodulatory properties [22–27]. In addition to these metabolic functions, many genetic studies in IBD have highlighted the central role of host-microbe interactions in IBD pathogenesis [1, 28–30]. Specific host pathways linked to microbial response in IBD include T-cell activation, the IL-23/T helper 17 pathway, autophagy , and Paneth cell function . Together, these results support the centrality of host-microbiota crosstalk for gut homeostasis and in turn the role of dysfunctional crosstalk between the host and GI microbiome in IBD.
Little work has yet bridged the gap between IBD pathogenesis in a human host, individual microbes, and alterations in metabolism of the GI microbial community in IBD. Few studies of the IBD gut microbiome have investigated microbiome function , and these have not systematically accounted for the influences of treatments and environmental factors. We have thus analyzed the GI microbiomes of 121 CD patients, 75 UC patients, and 27 healthy controls using a novel multivariate metagenomic analysis pipeline specifically accounting for environmental factors (including treatment, age, and tobacco use). In addition to assessing microbiome composition, we have analyzed the inferred metagenome as determined from phylogenetically-associated reference genomes, including metabolic modules and pathways also associated with disease status and with environmental factors such as medications and smoking. Not only were these GI microbiomes characterized by shifts in bacterial populations during disease as previously described, but these dysbioses were highly functionally coordinated. Cross-species enrichments included mucin metabolism and redox tolerance by means of glutathione transport, cysteine biosynthesis, and riboflavin metabolism. Conversely, processes linked broadly to clades IV and XIVa Clostridia were depleted, particularly short chain fatty acid production. Dysbioses in IBD are correspondingly not simply structural changes in the gut microbiota, but are instead associated with major impairments in many fundamental microbial metabolic functions with potential impact on the host.
Characteristics of patients in this study
Female gender (n)
Age (lower 95%-upper 95%)
Active disease (n)a
Assessing significant covariation of microbiome structure with host IBD status, treatment, and environment
An important consideration that informed the remainder of our analysis, and which is often overlooked in studies of the microbiome, was the consistent covariation among disease status, aspects of subject environment, and microbiome structure. For example, the factor most associated with changes in microbiome composition was not disease but whether the sample origin was stool or biopsy. Biopsy location induced minor changes in microbiome composition (Additional files 2, 3 and 4) relative to the extreme differences between stool and biopsy communities, in agreement with previous studies [36, 37]. In this cohort, iCD was always represented by biopsy, whereas 18.4% of non-iCD and 36% of UC samples were biopsies. iCD was also associated with greater likelihood of immunosuppressant treatment: iCD, non-iCD, and UC patients were treated by immunosuppressants in 74.4%, 19.2% and 16% of samples, respectively. In contrast, non-iCD and UC cases were more likely to be treated with mesalamine or antibiotics: mesalamine was used for 30.2% of iCD samples, 69.2% of non-iCD samples, and 77.3% of UC samples, while antibiotics were used in 2.3% of iCD, 17.9% of non-iCD, and 13.3% of UC samples. These associations lead to a range of non-independent covariates. Although disease activity may influence microbiome composition, after adjusting for the other factors, it was not independently associated with a specific shift in the microbiome composition in our analysis, and there were no significant (P < 0.01) associations between microbiome composition and gender (Additional file 5).
The second largely independent factor influencing microbiome composition was age, itself negatively associated with smoking (Figure 1; Additional file 1). Twenty-four (10.4%) of the available subjects were less than 18 years of age and 26 were 60 years or older. Aging is associated with continual changes in the microbiome, primarily a gradual decrease in Bifidobacterium as observed here (Additional file 6) and by others [38, 39]. After observing these overall patterns of covariation among disease, treatment, environment, and gut microbiome composition, we continued our analysis only after assessing the significance of microbiome-disease associations in a multivariate manner to account for host environment and treatment.
Microbial clades differentially abundant specifically in IBD include Roseburia, the Ruminococcaceae, and the Enterobacteriaceae
The Ruminococcaceae, which are acetate producers , were decreased in CD, while the Leuconostocaceae, which produce acetate and lactate , were decreased in UC. The only major clade with a significant increase in abundance specific to CD was the Enterobacteriaceae, specifically Escherichia/Shigella. This family has been previously implicated in intestinal inflammation [6, 45–47].
Crohn's disease with ileal involvement presents a distinct microbiome phenotype including reduced Faecalibacterium, and Odoribacter is reduced both in iCD and in pancolonic UC
In CD patients with ileal involvement, sequences of the Ruminococcaceae family and of Faecalibacterium in particular were dramatically reduced compared to other subjects (Figure 2a), confirming previous studies [48, 49]. Faecalibacterium prausnitzii, the only cultured representative of Faecalibacterium, is able to metabolize both diet-derived polysaccharides and host-derived substrates such as N-acetyl glucosamine from intestinal mucus . It is also a major butyrate producer and exhibits anti-inflammatory effects in a colitis setting . The Ruminococcaceae represent the first step of microbiome-linked carbohydrate metabolism, as they degrade several types of polysaccharides in the lower GI tract, including starch, cellulose, and xylan . The Roseburia genus, which is significantly reduced in all IBD patients (including iCD), and the Ruminococcaceae are further functionally connected in that the latter consume hydrogen and produce acetate that can be utilized by Roseburia to produce butyrate [41, 43]. Consistent reductions in all of these clades may thus have functional consequences on the ability of the host to repair the epithelium and to regulate inflammation.
The genera Escherichia/Shigella (indistinguishable as a 16S-based phylotype) were particularly highly enriched in iCD (q < 0.2; Additional file 1) above their general overabundance in CD patients. Lipopolysaccharide produced by Gram-negative bacteria such as Escherichia coli is a canonical microbe-associated molecular pattern, known to activate toll-like receptor 4 (TLR4) signaling  and thus trigger inflammatory cascades. TLR4 expression is highly up-regulated in the intestinal epithelium of IBD patients , and mutations in TLR4 are associated with both CD and UC . Previous culture-based studies have found that E. coli, specifically E. coli exhibiting pathogen-like behaviors such as adhesion and invasiveness , are more frequently cultured from iCD biopsies, and culture-independent studies have found an enrichment in E. coli that contain virulence-associated genes in iCD . This suggests that CD-involved ileum is a favorable milieu for establishment of E. coli with pathobiont features, which may have implications for IBD exacerbations and its chronicity. An inflamed ileum may furnish a specialized niche permissive for microbes with enhanced fitness in inflamed conditions.
The most severe form of UC is pancolitis, in which UC affects the entire colon; this condition is associated with greatly increased risk of colon cancer . Patients with pancolitis did not harbor a clear specificity in their dysbiosis. However, both these patients and iCD patients had a reduced abundance of the Odoribacter genus, which belongs to the Porphyromonadaceae family and to the Bacteroidetes phylum. As Odoribacter splanchnus is a known producer of acetate, propionate, and butyrate , decreased Odoribacter may affect host inflammation via reduced SCFA availability.
Microbiome composition is also strongly associated with subject age, treatment, smoking, and sample biogeography
Critical to determining causality in links between IBD and the gut microbiome, IBD treatments were also associated with alterations in microbiome composition. Mesalamine (5-aminosalicylic acid) is a bowel-specific aminosalicylate drug. Although its exact mode of action is unknown, it is thought to act as an antioxidant and to decrease intestinal inflammation, in part by peroxisome proliferator-activated receptor-γ (PPARγ) activation and inhibition of NFκB and pro-inflammatory eicosanoid production. Here, its use was linked to strong reductions in Escherichia/Shigella (> 100% of average abundance, q < 0.04; Additional file 7), in agreement with a recent study . Both 5-aminosalicylic acid and immunosuppressant treatment were associated with modest increases in Enterococcus, the only genus perturbed in immunosuppressant-treated patients with low false discovery rate (also > 100% of average abundance, q < 0.09).
Antibiotics were among the strongest factors associated with a reduction in ecological diversity (Figure 2b). Many individual clades were greatly reduced or nearly absent after administration of antibiotics, including the Collinsella, Dorea, Butyricicoccus, Subdoligranulum, and Acetivibrio (all q < 0.2; Additional file 1). These genera are predominantly from the Clostridiales order, Gram-positive and anaerobic bacteria that are targeted by the antibiotics commonly used in IBD, such as ciprofloxacin and metronidazole.
Smoking is likely the best-known environmental factor that impacts IBD . It is associated with increased risk of CD and is conversely protective towards developing UC . The only common organism to which tobacco usage was linked in these individuals was Anaerostipes (Firmicutes phylum), which decreased (> 60% average abundance, q < 0.15; Figure 3) in current or former tobacco users, beyond any change due solely to smokers' higher average age. The Anaerostipes genera can utilize lactate to produce butyrate , which is beneficial to colonic health.
Finally, as previously mentioned, samples of the stool as opposed to mucosal biopsies differed strongly in microbiome composition (Additional file 2). More than 70 clades were significantly over- or under-enriched in stool samples relative to biopsies at q < 0.2. This effect extended to entire phyla, as the Firmicutes were approximately twofold more abundant in stool (Additional file 1). Microbial habitat dictates the composition of microbial communities ; in the GI tract, this has been suggested to occur on biogeographical scales of intestinal regions [37, 63] or even millimeters apart [64, 65], and luminal/mucosal differences may be further perturbed by bowel preparation prior to colonoscopy . The data did not suggest that the luminal and mucosal communities were independent; rather, all 14 clades significantly associated with IBD retained the same trend when stratified by sample origin (Additional file 8). The fecal microbiome appeared to convey a consistent but numerically transformed function of mucosal communities, both of which shifted in composition in association with host environment, treatment, and disease.
In a closer analysis of intestinal biogeography as reflected by biopsies drawn from distinct regions, differences in most clades were modest and correlated largely with previously described changes in pH (Additional files 2, 3 and 9) . The clades with the largest regional changes included the Roseburia and Ruminococcaceae, with lower abundance in the low-pH terminal ileum, transverse, and right colon; Alistipes, following a similar pattern; and the Fusobacteria and Enterobacteriaceae, with an opposite pattern of somewhat increased abundance in the ileum and right colon. Particularly as the former have also been associated with the colorectal cancer microenvironment in previous work [68, 69], it is of note that these variations in the microbiota with respect to biogeography and pH are similar to those we observed with respect to IBD and potentially redox status as detailed below.
The metagenomic abundances of microbial metabolic pathways are more consistently perturbed in IBD than are organismal abundances
Considering only the contrast between IBD (CD or UC) and healthy subjects, 24 of 200 (12%) total metabolic modules were differentially abundant at q < 0.2. This is in stark contrast to the microbial shifts discussed above, in which only 6 of 263 (2%) genus-level clades reached this significance threshold. Even in the absence of metagenomic or metatranscriptomic data and only leveraging the genes and pathways in reference genomes associated with these communities, changes in microbial function were more consistent than changes in community structure. This has been noted in environmental communities  and suggested with respect to obesity and other biometrics [73, 74], but to date it has not been reported for disease-linked dysbioses or IBD.
We validated these functional shifts by shotgun metagenomic sequencing of the small subset of available samples with appropriate stool DNA, seven healthy controls and four CD patients (Additional file 14). These were sequenced to a shallow depth averaging 119 meganucleotides per sample of 150-nucleotide paired-end Illumina MiSeq reads, reducing our effective limit of detection but otherwise providing close agreement with inferred metabolic shifts in the IBD metagenome. Of the modules highlighted below and in Figure 5, one (cobalamin biosynthesis) fell below the limit of detection, and the remaining six retained the expected trend of over- or under-enrichment in Crohn's disease, as did additional processes detailed below, including glycolysis and bacterial secretion.
Amino acid biosynthesis and carbohydrate metabolism are reduced in the IBD microbiome in favor of nutrient uptake
We observed that even basic GI microbiome metabolism was altered in both UC and CD. Amino acid metabolism showed major perturbation: genes for the metabolism and biosynthesis of nearly all amino acids (particularly histidine and lysine) decreased in abundance (Figure 4), while arginine, histidine, and lysine transport (Figure 5) gene abundance increased. In iCD we also observed a decrease in glutamine-related functional modules, which would lead to a lower amount of glutamate required for gamma-aminobutyric acid, ornithine, and arginine biosynthesis; abundance of all three of these modules also decreased. In marked contrast to the other amino acids, genes for metabolism of the sulfur-containing amino acid cysteine significantly increased in abundance, with even greater increase in iCD. This corresponded with an overrepresentation of genes related to sulfate transport in UC and CD (Figure 5), and in increase in sulfur and nitrogen metabolism in CD (Figure 4).
CD was associated with increased abundance of many genes related to carbohydrate transport (Figure 5). There were large increases in pentose phosphate pathway and fructose/mannose metabolism gene abundance in iCD (Figure 4), which were accompanied by increase in carbohydrate metabolism, but they were not significant in UC and CD. In addition, iCD showed increased abundance of transporter genes for glucose, hexoses, maltose, and mono-, di-, and oligosaccharides (Figure 5). We observed a decrease in both butanoate and propanoate metabolism in iCD (Figure 4), suggesting a potential decrease in SCFA production by the microbiome, possibly due to the observed decrease in Roseburia and Faecalibacterium.
We saw an increase in glutathione transport gene abundance in UC and CD (Figure 5) and an increase in glutathione metabolism gene abundance in UC. Glutathione is a tripeptide of cysteine and glutamate, synthesized by Proteobacteria and a few streptococci and enterococci , which allows bacteria to maintain homeostasis during oxidative or acid stress. Inflammatory cascades include production of highly reactive oxygen and nitrogen metabolites, which are greatly increased in active IBD . Lamina propia monocytes also release homocysteine during inflammation, which further contributes to oxidative stress; IBD is associated with higher levels of both mucosal and serum homocysteine . Thus, the increases in sulfate transport, cysteine metabolism, and glutathione metabolism may reflect a mechanism by which the gut microbiome addresses the oxidative stress caused by inflammation.
Extreme functional shifts in iCD include changes in redox metabolism, enrichment of signaling/secretion, and suggest a 'pathobiont-like' invasive metagenome
CD with ileal involvement exhibited specific dysfunction at the module level. It was associated with an increase in several modules involved in glycolysis and carbohydrate transport and metabolism (Figure 5). Conversely, iCD exhibited lower abundance of genes involved in lipid metabolism and catabolism, confirming a major imbalance in energy metabolism. We observed a global decrease in nicotinamide, purine, and pyrimidine nucleotide biosynthesis modules in iCD, CD, and UC (Figure 5).
There was a decrease in vitamin biosynthesis associated with iCD, but increases in thiamine and particularly riboflavin metabolism modules (Figure 4). Interestingly, this pathway is fed by the pentose phosphate pathway, which was also overrepresented in iCD. Riboflavin is necessary for regenerating oxidized glutathione back to its reduced form, and is thus essential for pH and oxidative stress homeostasis, as is NADPH, a product of the pentose phosphate pathway. Metabolism of the sulfur-containing amino acids cysteine and methionine was increased in iCD, in marked contrast to the IBD-associated decreases in the non-sulfur-containing amino acids such as lysine and glutamine. As homocysteine is easily convertible to methionine, this may indicate a further mechanism of maintaining redox homeostasis. Alternatively, this may be connected to the iCD-specific increase in carbohydrate metabolism, as cysteine may be metabolized to pyruvate.
Finally, genes involved in pathogenesis processes, such as secretion systems and adherence/invasion, were overrepresented in iCD (Figure 4). For example, genes involved in the shigellosis pathway were more abundant in CD, and type II secretion genes were more abundant in iCD. Type II secretion is involved in the secretion of cell wall-degrading enzymes  and the secretion of toxins such as heat-labile enterotoxin, similar to cholera toxin . These functions are typical of pathobiont adherent-invasive E. coli, which have been observed to increase in iCD in our own study and others [6, 55]. This may be associated with tissue damage, either primarily as a result of toxin secretion, or secondarily as a result of stimulated cytokine production. This tissue destruction is a likely source of metabolites for microbial overgrowth, selecting for auxotrophic specialists able to thrive in this environment and resulting in the microbiome-wide loss of basic biosynthetic processes (Figures 4 and 5). This would in turn lead to further tissue breakdown, bacterial overgrowth, and community structural and functional dysbiosis.
The GI microbiome influences dietary energy extraction, immune system development, vitamin production, and drug metabolism, yet most molecular and metabolic functions of the bacteria of the GI microbiome are uncharacterized . To gain insight into the functional consequences of IBD-associated dysbiosis, we used a novel approach pairing microbial community 16S gene sequence profiles with information from the closest available whole-genome sequences. This defined an inferred metagenome and thus complement of metabolic functional modules for each microbiome in this study. This allowed us to identify unique functional perturbations in the microbiomes of IBD patients. Interestingly, although we identified only nine changes in bacterial clades that associated with UC (of 350 total, 2.6%), we identified 21 statistically significant differences in functional pathways and metabolic modules (of 295, 7.1%); this pattern held for CD and iCD function as well. This underscores the fact that phylogenetically diverse changes in the composition of the GI microbiome can be functionally coordinated and lead to major modifications in the metabolic potential of the microbiota.
The microbial metabolic information available in this study represents only one step in the functional investigation of the IBD microbiota, as it is an accurate but approximate inference using prior knowledge of microbial genomes. The metagenomes inferred from our 16S data were supported by shotgun sequencing of a subset of samples, providing one confirmation that they were representative of community functional capability. As sequencing costs continue to fall, rich metagenomic data for dozens or hundreds of samples will further improve our ability to resolve species-level gene function in communities. Of course, a community expresses only a variable subset of its functional capability at any given time, in response to environmental stimuli. Thus, metatranscriptomic, proteomic, and metabolomic data will continue to add to our understanding of which of a community's potential functions are most strongly affecting the host during inflammatory disease.
Alternatively, the increased biosynthesis of cysteine (a precursor of glutathione) and of glutathione transport modules may speak to the microbiome's response to the oxidative stress (high levels of reactive oxygen and nitrogen species) of the inflamed IBD gut . In support of this concept, we found that riboflavin metabolism, which is required to convert glutathione between its oxidized and reduced forms, is increased in iCD. Furthermore, the pentose phosphate pathway, which produces the NADPH also required for glutathione reduction, is increased as well. Recent studies have shown that redox stress allows Salmonella to use ethanolamine as a carbon source  and allows enterohemorrhagic E. coli to use it as a nitrogen source , thus conferring a competitive advantage to these microbes. This raises the interesting possibility that E. coli or related species in IBD may be highly represented because they gain a competitive advantage from oxidative stress and are better able to compensate for it with glutathione production.
In both UC and CD, there were decreases in the biosynthesis of lysine, arginine, and histidine in favor of transport in both UC and CD; a further decrease in tryptophan metabolism was associated with iCD. The data showed additional broad decreases in many essential processes, such as cobalamin synthesis, purine and pyrimidine biosythesis, lipid catabolism, and phospholipid metabolism, as well as marked increases in transport. This overall decrease in abundance of genes for amino acid and nucleotide biosynthesis bears striking resemblance to the lifestyle of highly symbiotic bacteria that are intrinsically auxotropic and also of some pathobionts (Figure 6). One such example are segmented filamentous bacteria (SFB), a symbiont that belongs to the Candidadatus Arthromitis, a sub-group of clade I (sensu stricto) Clostridia. A recently sequenced SFB genome lacked genes for nucleotide biosynthesis as well as nearly all vitamins and amino acids [82, 83]. SFB are often abundant in the rodent terminal ileum and are responsible for the maturation of Th17 cells , which play an important role in CD-associated inflammation . To date, neither SFB nor phylogenetically related sequences have been observed in humans [82, 86]; this was also true in our data (zero 16S sequences with > 90% identity to X77814 SFB). However, a functional trend similar to SFB was observed in these IBD community metagenomes, as biosynthetic mechanisms throughout central carbon metabolism, amino acid biosynthesis, and nucleotide maintenance were all reduced (Figures 4 and 5), hinting that humans may host functional equivalents of SFB-like pathobionts that increase in IBD but are not phylogenetically close to Candidatus Arthromitis. Host tissue destruction, either inflammation-mediated or bacterially mediated, would provide a ready nutrient source (Figure 6).
The data presented here show that IBD and iCD in particular are associated with a dysbiosis characterized by changes in Firmicutes and Proteobacteria phyla. Environmental factors and, notably, treatments were also associated with independent changes in the GI microbiome; these must be taken into account during future studies of the microbiota in IBD. These perturbations in bacterial composition, although modest, were associated with major perturbations of GI microbiome function, which revolved around metabolism in the presence of oxidative stress and perturbed nutrient availability during tissue damage. Further studies, particularly including transcriptomic, proteomic, or metabolomic characterization, longitudinal data, and dietary metadata, will be needed to additionally define the consequences of the IBD-associated microbiome dysfunction on the host and the specific mechanisms by which they are carried out or regulated by the microbiota.
Materials and methods
The OSCCAR and PRISM cohorts
The Ocean State Crohn's and Colitis Area Registry (OSCCAR) is a state-based, prospective inception cohort of IBD patients that was designed to study the epidemiology of IBD, to determine the incidence of IBD in Rhode Island, and to extrapolate these rates to the general population of the United States. The diverse population of over 1 million, limited geographic range, and well-circumscribed gastroenterology community of Rhode Island were ideal circumstances for establishing a prospective inception cohort of IBD patients. All but one of the 98 gastroenterologists/colorectal surgeons in Rhode Island agreed to refer patients to OSCCAR, and 11 gastroenterologists practicing in Massachusetts just over the Rhode Island border also agreed to refer their newly diagnosed IBD patients who resided in Rhode Island. Enrollment began 1 January 2008. All Rhode Island residents with a newly confirmed diagnosis of CD, UC, or indeterminate colitis were eligible for inclusion (within 12 months from diagnosis). Ethnic background of the subjects was not available for consideration in the analysis, and indeterminate colitis patients were analyzed only for other metadata and not for IBD diagnosis. Diagnosis of CD, UC, or indeterminate colitis was made by endoscopic, pathologic, or radiographic findings according to the criteria of the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) IBD Genetics Consortium. OSCCAR research protocols were reviewed and approved by three institutional review boards (Lifespan (#0214-07), the Partners Human Research Committee (#2007-P-001705), and the Program for the Protection of Human Subjects/Mount Sinai School of Medicine (#11-01479)), and all experiments adhered to the regulations of these review boards. Informed consent and HIPAA (Health Insurance Portability and Accountability Act) authorization were obtained from each subject prior to study participation. Individuals diagnosed with IBD prior to the study start date, pregnant women, those unwilling to provide informed consent for study participation, and those who were prisoners at the time of diagnosis were not permitted to enroll.
The Prospective Registry in IBD Study at MGH (PRISM) is a referral center-based, prospective cohort of IBD patients. Enrollment began 1 January 2005. Patients aged 18 years and older with a diagnosis of CD or UC based upon standard endoscopic, radiographic, and histologic criteria were eligible to participate. Controls consisted of healthy patients aged 18 years and older, from whom biopsies were obtained during colonoscopies performed for screening purposes.
Patients were excluded from the healthy volunteer group for current acute illness, if awaiting transplant, or if chronically ill (for example, renal failure, diabetes, congestive heart failure). During routine colonoscopies, subjects were offered the opportunity to donate biopsy samples. After sampling, intestinal biopsies were stored in 5% glycerol at -80°C until DNA extraction. Stool samples were kept at 4°C for less than 24 h before storage at -80°C until DNA extraction. PRISM research protocols were reviewed and approved by the Partners Human Research Committee (#2004-P-001067), and all experiments adhered to the regulations of this review board.
DNA from stool and biopsy samples was extracted using the QIAamp DNA Stool Mini Kit (Qiagen, Inc., Valencia, CA, USA) according to manufacturer's instructions and as described previously . The manufacturer's protocol was altered to accommodate larger stool volumes and to improve homogenization using bead-beating at several steps: a) a minimum of 2 ml of Buffer ASL and 300 mg of stool was used in the protocol; b) a ratio of 700 μl of Buffer ASL per 100 mg of stool weight was used for larger volumes using no more than 1,500 mg of stool and 10.5 ml of Buffer ASL; c) following the addition of Buffer ASL to each sample (step 2), 0.70 mm Garnet Beads (MO BIO Laboratories, Inc., Carlsbad, CA, USA) were added to the suspension and vortexed for 10 seconds; d) a second bead-beating was performed following the heating of the suspension (step 3) in 0.1 mm Glass Bead Tubes (MO BIO Laboratories, Inc.), and vortexed for 10 minutes.
Amplification and 454 sequencing of the 16S gene
The 16S gene dataset consists of 454 FLX Titanium sequences spanning the V3 to V5 variable regions. Detailed protocols used for 16S amplification and sequencing are available on the Human Microbiome Project Data Analysis and Coordination Center website . In brief, genomic DNA was subjected to 16S amplifications using primers designed incorporating the FLX Titanium adapters and a sample barcode sequence, allowing directional sequencing covering variable regions V5 to partial V3 (primers: 357F 5'-CCTACGGGAGGCAGCAG-3' and 926R 5' CCGTCAATTCMTTTRAGT-3'). PCR mixtures (25 μl) contained 10 ng of template, 1× Easy A reaction buffer (Stratagene, La Jolla, CA, USA), 200 mM of each dNTP (Stratagene), 200 nM of each primer, and 1.25 U AccuPrime hifi cloning enzyme (Invitrogen, Carlsbad, CA, USA). The cycling conditions for the V3-V5 consisted of an initial denaturation of 95°C for 2 minutes, followed by 25 cycles of denaturation at 95°C for 40 s, annealing at 50°C for 30 s, extension at 72°C for 5 minutes and a final extension at 72°C for 7 minutes. Amplicons were confirmed on 1.2% Flash Gels (Lonza, Rockland, ME, USA), purified with AMPure XP DNA purification beads (Beckman Coulter, Danvers, MA, USA) according to the manufacturer, and eluted in 25 μl of 1× low TE buffer (pH 8.0). Amplicons were quantified on Agilent Bioanalyzer 2100 DNA 1000 chips (Agilent Technologies, Santa Clara, CA, USA) and pooled in equimolar concentration. Emulsion PCR and sequencing were performed according to the manufacturer's specifications.
Processing sequencing samples
Sequences were processed in a data curation pipeline implemented in MOTHUR , which removed sequences from the analysis if they were less than 200 nucleotides or greater than 600 nucleotides, had a low read quality score (< 25), contained ambiguous characters, had a non-exact barcode match, or had more than 4 mismatches to the reverse primer sequences (926R). Remaining sequences were assigned to samples based on barcode matches, and barcode and primer sequences were then trimmed. Chimeric sequences were identified using the ChimeraSlayer  algorithm, and reads were classified with the MSU RDP classifier v2.2  using the taxonomy maintained at the Ribosomal Database Project (RDP 10 database, version 6). Sequencing depth after processing averaged 2,860 (standard deviation 1,730) reads per sample.
Metagenome inference from microbiome composition
for GG tree nodes g and h separated by phylogenetic branch length dist(g, h) and annotated with KO genome vectors and .
Using this vector representation of genomes, the abundance of an individual gene family (KO) i in a community due to the presence of a specific phylotype g is the product of the corresponding gene count and the measured abundance of phylotype g. Therefore, the total relative abundance of each KO was estimated for each sample by adding the individual contributions of all phylotypes present in the sample. Using this method, we inferred the functional composition for each sampled community. The inference process's accuracy was validated by comparing inferred KO abundances in 16S datasets from the Human Microbiome Project with their metagenomically sequenced counterparts (Additional file 15).
Metabolic pathway reconstruction
Inferred per-community gene (KO) abundances were subsequently reconstructed into microbial pathway relative abundances using HUMAnN, the Human Microbiome Project metabolic reconstruction pipeline . KOs were grouped into pathways represented as gene sets using HUMAnN, which chooses pathways by maximum parsimony using MinPath  and computes each pathway's relative abundance as a smoothed average over all genes within it, taking into account outliers and gap filling. We ran HUMAnN three times to reconstruct three complementary types of pathways from these genes: small metabolic modules (using KEGG's conjunctive normal form logic), large metabolic pathways, and Gene Ontology terms (using annotation-to-KO mappings from nine well-characterized KEGG microbes: ban, cje, cpe, eco, nse, pae, sce, son, and vch). For each of these three types of pathway, HUMAnN input the inferred relative abundances of all genes in each sample, and output the relative abundances of pathways within the sample. Subsequent analysis handled these sample-by-pathway relative abundances in the same manner as sample-by-clade microbial abundances.
Significant associations of microbial clades and pathways with sample metadata
Inverse Simpson diversity, Chao1 richness (using the R fossil package), and Pielou evenness were calculated for clade abundance, KEGG pathway and module abundance, and Gene Ontology term abundance [94–97]. Next, data were pre-processed for quality control before modeling. Clinical metadata were removed when more than 10% of data were missing, or when they did not vary in value over the available samples. Clades, pathways, and features of very low abundance (< 0.001 in ≥ 90% of samples) and feature outliers outside of the lower or upper outer fence (3× interquartile range) were removed. Missing data were imputed for significance testing with the mean abundance of the sample; missing factor metadata were imputed with a 'NA' factor level using the na.gam.replace function from the R package . Unless stated otherwise, all subsequent analyses and calculations were performed using these processed data. After processing, 228 and 231 samples passed quality control for clade abundance and functional abundance analyses, respectively.
where p are the clinical metadata selected from boosting.
Within each metadatum/clade association independently, multiple comparisons over factor levels were adjusted using a Bonferonni correction; multiple hypothesis tests over all clades and metadata were adjusted to produce a final Benjamini and Hochberg false discovery rate . Unless otherwise indicated, significant association was considered below a q-value threshold of 0.25; the KEGG pathway sulfur metabolism (ko00920) had an average q-value of 0.26 for association with Crohn's disease. Multiple factor analysis was performed to visualize the relationships within heterogeneous factor data as well as with a select group of taxa found to be significantly associated with metadata (using the FactoMineR R package ). Total abundances and significant associations between metadata, taxa, and functions are listed in Additional files 1 and 11.
Sequence alignment for segmented filamentous bacteria
To determine whether SFB were present in samples, three sequences of SFB (X80834, X87244, and X77814) from three species (chicken, rat, and mouse) were aligned by blastn, using both a 20 and 15 seed word. No sequences were found with ≥ 95% identity over an alignment length of at least 100 nucleotides. The average sequence length from the study was 435 nucleotides.
Shotgun metagenomic sequencing
To provide internal validation of inferred microbial community gene and pathway compositions, stool DNA from seven healthy controls and four CD patients was subjected to metagenomic shotgun sequencing. Libraries were constructed with the Illumina Nextera XT kit and sequenced on an Illumina MiSeq using 2 × 150 bp paired-end sequencing according to the manufacturer's instructions. This resulted in sequencing depths ranging from 3.9 to 270 meganucleotides, average 119 meganucleotides, from which microbial community function was determined with HUMAnN  as described above.
Sequence accession numbers and availability
Sequences generated in this study are publicly available (NCBI BioProject ID numbers 82111 and 175224).
inflammatory bowel disease
ileal Crohn's disease
Kyoto Encyclopedia of Genes and Genomes
Ocean State Crohn's and Colitis Area Registry
Prospective Registry in IBD Study at MGH
Ribosomal Database Project
short-chain fatty acid
segmented filamentous bacteria
We thank Wendy Garrett, Patrick Veiga, Ken Croitoru, and Mark Silverberg for their valuable comments on this manuscript and study, as well as Bruce Birren and the Broad Genome Sequencing and Analysis Program. This work was funded in part by grants from the French Society of Gastroenterology, the Bettencourt Schueller Foundation, the Philippe Foundation, and the Arthur Sachs Scholarship (HS); NIH (NIHU54HG004969) and the Crohn's and Colitis Foundation of America (DG and DVW); the philanthropic support of the MacInnes and Rasmussen families; NIH 1R21HD058828-01A1 to NL; and NIH 1R01HG005969, NSF DBI-1053486, and ARO W911NF-11-1-0473 and W911NF-11-1-0429 to CH. Dr. Sands acknowledges NIH (7 R21 DK078555), Centers for Disease Control and Prevention (5 U01 DP002676, 5 U01 DP000340), and the Crohn's and Colitis Foundation of America
- Khor B, Gardet A, Xavier RJ: Genetics and pathogenesis of inflammatory bowel disease. Nature. 2011, 474: 307-317. 10.1038/nature10209.PubMedPubMed CentralView ArticleGoogle Scholar
- Cosnes J, Gower-Rousseau C, Seksik P, Cortot A: Epidemiology and natural history of inflammatory bowel diseases. Gastroenterology. 2011, 140: 1785-1794. 10.1053/j.gastro.2011.01.055.PubMedView ArticleGoogle Scholar
- Goodman AL, Kallstrom G, Faith JJ, Reyes A, Moore A, Dantas G, Gordon JI: Extensive personal human gut microbiota culture collections characterized and manipulated in gnotobiotic mice. Proc Natl Acad Sci USA. 2011, 108: 6252-6257. 10.1073/pnas.1102938108.PubMedPubMed CentralView ArticleGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Structure, function and diversity of the healthy human microbiome. Nature. 2012, 486: 207-214. 10.1038/nature11234.
- Baumgart M, Dogan B, Rishniw M, Weitzman G, Bosworth B, Yantiss R, Orsi RH, Wiedmann M, McDonough P, Kim SG, et al: Culture independent analysis of ileal mucosa reveals a selective increase in invasive Escherichia coli of novel phylogeny relative to depletion of Clostridiales in Crohn's disease involving the ileum. Isme J. 2007, 1: 403-418. 10.1038/ismej.2007.52.PubMedView ArticleGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Manichanh C, Rigottier-Gois L, Bonnaud E, Gloux K, Pelletier E, Frangeul L, Nalin R, Jarrin C, Chardon P, Marteau P, et al: Reduced diversity of faecal microbiota in Crohn's disease revealed by a metagenomic approach. Gut. 2006, 55: 205-211. 10.1136/gut.2005.073817.PubMedPubMed CentralView ArticleGoogle Scholar
- Ott SJ, Musfeldt M, Wenderoth DF, Hampe J, Brant O, Folsch 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Packey CD, Sartor RB: Commensal bacteria, traditional and opportunistic pathogens, dysbiosis and bacterial killing in inflammatory bowel diseases. Curr Opin Infect Dis. 2009, 22: 292-301. 10.1097/QCO.0b013e32832a8a5d.PubMedPubMed CentralView ArticleGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Sokol H, Lay C, Seksik P, Tannock GW: Analysis of bacterial bowel communities of IBD patients: what has it revealed?. Inflamm Bowel Dis. 2008, 14: 858-867. 10.1002/ibd.20392.PubMedView ArticleGoogle Scholar
- Li E, Hamm CM, Gulati AS, Sartor RB, Chen H, Wu X, Zhang T, Rohlf FJ, Zhu W, Gu C, et al: Inflammatory Bowel Diseases Phenotype, C. difficile and NOD2 Genotype Are Associated with Shifts in Human Ileum Associated Microbial Composition. PLoS ONE. 2012, 7: e26284-10.1371/journal.pone.0026284.PubMedPubMed CentralView ArticleGoogle Scholar
- Sokol H, Seksik P: The intestinal microbiota in inflammatory bowel diseases: time to connect with the host. Curr Opin Gastroenterol. 2010, 26: 327-331. 10.1097/MOG.0b013e328339536b.PubMedView ArticleGoogle 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.PubMedView ArticleGoogle Scholar
- Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermudez-Humaran LG, Gratadoux JJ, Blugeon S, Bridonneau C, Furet JP, Corthier G, et al: 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Sokol H, Seksik P, Furet JP, Firmesse O, Nion-Larmurier I, Beaugerie L, Cosnes J, Corthier G, Marteau P, Dore J: Low counts of Faecalibacterium prausnitzii in colitis microbiota. Inflamm Bowel Dis. 2009, 15: 1183-1189. 10.1002/ibd.20903.PubMedView ArticleGoogle Scholar
- Gibson GR, Cummings JH, Macfarlane GT: Growth and activities of sulphate-reducing bacteria in gut contents of healthy subjects and patients with ulcerative colitis. FEMS Microbiol Ecol. 1991, 9: 103-111. 10.1111/j.1574-6941.1991.tb01742.x.View ArticleGoogle Scholar
- Roediger WE, Moore J, Babidge W: Colonic sulfide in pathogenesis and treatment of ulcerative colitis. Dig Dis Sci. 1997, 42: 1571-1579. 10.1023/A:1018851723920.PubMedView ArticleGoogle Scholar
- Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, et al: A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010, 464: 59-65. 10.1038/nature08821.PubMedPubMed CentralView ArticleGoogle Scholar
- Flint HJ, Bayer EA, Rincon MT, Lamed R, White BA: Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis. Nat Rev Microbiol. 2008, 6: 121-131. 10.1038/nrmicro1817.PubMedView ArticleGoogle Scholar
- Tedelind S, Westberg F, Kjerrulf M, Vidal A: Anti-inflammatory properties of the short-chain fatty acids acetate and propionate: a study with relevance to inflammatory bowel disease. World J Gastroenterol. 2007, 13: 2826-2832.PubMedPubMed CentralGoogle Scholar
- Clausen MR, Mortensen PB: Kinetic studies on colonocyte metabolism of short chain fatty acids and glucose in ulcerative colitis. Gut. 1995, 37: 684-689. 10.1136/gut.37.5.684.PubMedPubMed CentralView ArticleGoogle Scholar
- Inan MS, Rasoulpour RJ, Yin L, Hubbard AK, Rosenberg DW, Giardina C: The luminal short-chain fatty acid butyrate modulates NF-kappaB activity in a human colonic epithelial cell line. Gastroenterology. 2000, 118: 724-734. 10.1016/S0016-5085(00)70142-9.PubMedView ArticleGoogle Scholar
- Saemann MD, Bohmig GA, Osterreicher CH, Burtscher H, Parolini O, Diakos C, Stockl J, Horl WH, Zlabinger GJ: Anti-inflammatory effects of sodium butyrate on human monocytes: potent inhibition of IL-12 and up-regulation of IL-10 production. Faseb J. 2000, 14: 2380-2382.PubMedGoogle Scholar
- Segain JP, Raingeard de la Bletiere D, Bourreille A, Leray V, Gervois N, Rosales C, Ferrier L, Bonnet C, Blottiere HM, Galmiche JP: Butyrate inhibits inflammatory responses through NFkappaB inhibition: implications for Crohn's disease. Gut. 2000, 47: 397-403. 10.1136/gut.47.3.397.PubMedPubMed CentralView ArticleGoogle Scholar
- Velazquez OC, Lederer HM, Rombeau JL: Butyrate and the colonocyte. Production, absorption, metabolism, and therapeutic implications. Adv Exp Med Biol. 1997, 427: 123-134. 10.1007/978-1-4615-5967-2_14.PubMedView ArticleGoogle Scholar
- Barrett JC, Hansoul S, Nicolae DL, Cho JH, Duerr RH, Rioux JD, Brant SR, Silverberg MS, Taylor KD, Barmada MM, 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Franke A, McGovern DP, Barrett JC, Wang K, Radford-Smith GL, Ahmad T, Lees CW, Balschun T, Lee J, Roberts R, 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.PubMedPubMed CentralView ArticleGoogle Scholar
- McGovern DP, Gardet A, Torkvist L, Goyette P, Essers J, Taylor KD, Neale BM, Ong RT, Lagace C, Li C, et al: Genome-wide association identifies multiple ulcerative colitis susceptibility loci. Nat Genet. 2010, 42: 332-337. 10.1038/ng.549.PubMedPubMed CentralView ArticleGoogle Scholar
- Hampe J, Franke A, Rosenstiel P, Till A, Teuber M, Huse K, Albrecht M, Mayr G, De La Vega FM, Briggs J, et al: A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1. Nat Genet. 2007, 39: 207-211. 10.1038/ng1954.PubMedView ArticleGoogle Scholar
- Vaishnava S, Behrendt CL, Ismail AS, Eckmann L, Hooper LV: Paneth cells directly sense gut commensals and maintain homeostasis at the intestinal host-microbial interface. Proc Natl Acad Sci USA. 2008, 105: 20858-20863. 10.1073/pnas.0808723105.PubMedPubMed CentralView ArticleGoogle Scholar
- Greenblum S, Turnbaugh PJ, Borenstein E: Metagenomic systems biology of the human gut microbiome reveals topological shifts associated with obesity and inflammatory bowel disease. Proceedings of the National Academy of Sciences of the United States of America. 2012, 109: 594-599. 10.1073/pnas.1116053109.PubMedPubMed CentralView ArticleGoogle Scholar
- Frank DN, Robertson CE, Hamm CM, Kpadeh Z, Zhang T, Chen H, Zhu W, Sartor RB, Boedeker EC, Harpaz N, et al: Disease phenotype and genotype are associated with shifts in intestinal-associated microbiota in inflammatory bowel diseases. Inflamm Bowel Dis. 2010Google Scholar
- Lepage P, Hasler R, Spehlmann ME, Rehman A, Zvirbliene A, Begun A, Ott S, Kupcinskas L, Dore J, Raedler A, Schreiber S: Twin study indicates loss of interaction between microbiota and mucosa of patients with ulcerative colitis. Gastroenterology. 2011, 141: 227-236. 10.1053/j.gastro.2011.04.011.PubMedView ArticleGoogle Scholar
- Costello EK, Lauber CL, Hamady M, Fierer N, Gordon JI, Knight R: Bacterial community variation in human body habitats across space and time. Science. 2009, 326: 1694-1697. 10.1126/science.1177486.PubMedPubMed CentralView ArticleGoogle Scholar
- Stearns JC, Lynch MDJ, Senadheera DB, Tenenbaum HC, Goldberg MB, Cvitkovitch DG, Croitoru K, Moreno-Hagelsieb G, Neufeld JD: Bacterial biogeography of the human digestive tract. Sci Rep. 2011, 1:Google Scholar
- Agans R, Rigsbee L, Kenche H, Michail S, Khamis HJ, Paliy O: Distal gut microbiota of adolescent children is different from that of adults. FEMS Microbiol Ecol. 2011, 77: 404-412. 10.1111/j.1574-6941.2011.01120.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Mariat D, Firmesse O, Levenez F, Guimaraes V, Sokol H, Dore J, Corthier G, Furet JP: The Firmicutes/Bacteroidetes ratio of the human microbiota changes with age. BMC Microbiol. 2009, 9: 123-10.1186/1471-2180-9-123.PubMedPubMed CentralView ArticleGoogle Scholar
- Atarashi K, Tanoue T, Shima T, Imaoka A, Kuwahara T, Momose Y, Cheng G, Yamasaki S, Saito T, Ohba Y, et al: Induction of colonic regulatory T cells by indigenous Clostridium species. Science. 2011, 331: 337-341. 10.1126/science.1198469.PubMedPubMed CentralView ArticleGoogle Scholar
- Duncan SH, Hold GL, Barcenilla A, Stewart CS, Flint HJ: Roseburia intestinalis sp. nov., a novel saccharolytic, butyrate-producing bacterium from human faeces. Int J Syst Evol Microbiol. 2002, 52: 1615-1620. 10.1099/ijs.0.02143-0.PubMedGoogle Scholar
- Watanabe Y, Nagai F, Morotomi M: Characterization of Phascolarctobacterium succinatutens sp. nov., an asaccharolytic, succinate-utilizing bacterium isolated from human feces. Appl Environ Microbiol. 2012, 78: 511-518. 10.1128/AEM.06035-11.PubMedPubMed CentralView ArticleGoogle Scholar
- Chassard C, Bernalier-Donadille A: H2 and acetate transfers during xylan fermentation between a butyrate-producing xylanolytic species and hydrogenotrophic microorganisms from the human gut. FEMS Microbiol Lett. 2006, 254: 116-122. 10.1111/j.1574-6968.2005.00016.x.PubMedView ArticleGoogle Scholar
- Cogan JG, Jordan KN: Metabolism of Leuconostoc bacteria. J Dairy Sci. 1994, 77: 2704-2717. 10.3168/jds.S0022-0302(94)77213-1.View ArticleGoogle Scholar
- Garrett WS, Gallini CA, Yatsunenko T, Michaud M, DuBois A, Delaney ML, Punit S, Karlsson M, Bry L, Glickman JN, et al: Enterobacteriaceae act in concert with the gut microbiota to induce spontaneous and maternally transmitted colitis. Cell Host Microbe. 2010, 8: 292-300. 10.1016/j.chom.2010.08.004.PubMedPubMed CentralView ArticleGoogle Scholar
- Kleessen B, Kroesen AJ, Buhr HJ, Blaut M: Mucosal and invading bacteria in patients with inflammatory bowel disease compared with controls. Scand J Gastroenterol. 2002, 37: 1034-1041. 10.1080/003655202320378220.PubMedView ArticleGoogle Scholar
- Mylonaki M, Rayment NB, Rampton DS, Hudspith BN, Brostoff J: Molecular characterization of rectal mucosa-associated bacterial flora in inflammatory bowel disease. Inflamm Bowel Dis. 2005, 11: 481-487. 10.1097/01.MIB.0000159663.62651.4f.PubMedView ArticleGoogle 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.PubMedView ArticleGoogle Scholar
- Mondot S, Barreau F, Al Nabhani Z, Dussaillant M, Le Roux K, Dore J, Leclerc M, Hugot JP, Lepage P: Altered gut microbiota composition in immune-impaired Nod2-/- mice. Gut. 2011Google Scholar
- Lopez-Siles M, Khan TM, Duncan SH, Harmsen HJ, Garcia-Gil LJ, Flint HJ: Cultured Representatives of Two Major Phylogroups of Human Colonic Faecalibacterium prausnitzii Can Utilize Pectin, Uronic Acids, and Host-Derived Substrates for Growth. Appl Environ Microbiol. 2012, 78: 420-428. 10.1128/AEM.06858-11.PubMedPubMed CentralView ArticleGoogle Scholar
- Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermudez-Humaran LG, Gratadoux JJ, Blugeon S, Bridonneau C, Furet JP, Corthier G, et al: 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Poltorak A, He X, Smirnova I, Liu MY, Van Huffel C, Du X, Birdwell D, Alejos E, Silva M, Galanos C, et al: Defective LPS signaling in C3H/HeJ and C57BL/10ScCr mice: mutations in Tlr4 gene. Science. 1998, 282: 2085-2088.PubMedView ArticleGoogle Scholar
- Cario E, Podolsky DK: Differential alteration in intestinal epithelial cell expression of toll-like receptor 3 (TLR3) and TLR4 in inflammatory bowel disease. Infect Immun. 2000, 68: 7010-7017. 10.1128/IAI.68.12.7010-7017.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Franchimont D, Vermeire S, El Housni H, Pierik M, Van Steen K, Gustot T, Quertinmont E, Abramowicz M, Van Gossum A, Deviere J, Rutgeerts P: Deficient host-bacteria interactions in inflammatory bowel disease? The toll-like receptor (TLR)-4 Asp299gly polymorphism is associated with Crohn's disease and ulcerative colitis. Gut. 2004, 53: 987-992. 10.1136/gut.2003.030205.PubMedPubMed CentralView ArticleGoogle Scholar
- Darfeuille-Michaud A, Boudeau J, Bulois P, Neut C, Glasser AL, Barnich N, Bringer MA, Swidsinski A, Beaugerie L, Colombel JF: High prevalence of adherent-invasive Escherichia coli associated with ileal mucosa in Crohn's disease. Gastroenterology. 2004, 127: 412-421. 10.1053/j.gastro.2004.04.061.PubMedView ArticleGoogle Scholar
- Ekbom A, Helmick C, Zack M, Adami HO: Ulcerative colitis and colorectal cancer. A population-based study. N Engl J Med. 1990, 323: 1228-1233. 10.1056/NEJM199011013231802.PubMedView ArticleGoogle Scholar
- Goker M, Gronow S, Zeytun A, Nolan M, Lucas S, Lapidus A, Hammon N, Deshpande S, Cheng JF, Pitluck S, et al: Complete genome sequence of Odoribacter splanchnicus type strain (1651/6). Stand Genomic Sci. 2011, 4: 200-209. 10.4056/sigs.1714269.PubMedPubMed CentralView ArticleGoogle Scholar
- Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI: An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006, 444: 1027-1031. 10.1038/nature05414.PubMedView ArticleGoogle Scholar
- Benjamin JL, Hedin CR, Koutsoumpas A, Ng SC, McCarthy NE, Prescott NJ, Pessoa-Lopes P, Mathew CG, Sanderson J, Hart AL, et al: Smokers with active Crohn's disease have a clinically relevant dysbiosis of the gastrointestinal microbiota. Inflamm Bowel Dis. 2011Google Scholar
- Mahid SS, Minor KS, Soto RE, Hornung CA, Galandiuk S: Smoking and inflammatory bowel disease: a meta-analysis. Mayo Clin Proc. 2006, 81: 1462-1471. 10.4065/81.11.1462.PubMedView ArticleGoogle Scholar
- Cosnes J: Tobacco and IBD: relevance in the understanding of disease mechanisms and clinical practice. Best Pract Res Clin Gastroenterol. 2004, 18: 481-496. 10.1016/j.bpg.2003.12.003.PubMedView ArticleGoogle Scholar
- Sato T, Matsumoto K, Okumura T, Yokoi W, Naito E, Yoshida Y, Nomoto K, Ito M, Sawada H: Isolation of lactate-utilizing butyrate-producing bacteria from human feces and in vivo administration of Anaerostipes caccae strain L2 and galacto-oligosaccharides in a rat model. FEMS Microbiol Ecol. 2008, 66: 528-536. 10.1111/j.1574-6941.2008.00528.x.PubMedView ArticleGoogle Scholar
- Turnbaugh PJ, Ridaura VK, Faith JJ, Rey FE, Knight R, Gordon JI: The effect of diet on the human gut microbiome: a metagenomic analysis in humanized gnotobiotic mice. Sci Transl Med. 2009, 1: 6ra14-10.1126/scitranslmed.3000322.PubMedPubMed CentralGoogle Scholar
- Hong PY, Croix JA, Greenberg E, Gaskins HR, Mackie RI: Pyrosequencing-based analysis of the mucosal microbiota in healthy individuals reveals ubiquitous bacterial groups and micro-heterogeneity. PLoS One. 2011, 6: e25042-10.1371/journal.pone.0025042.PubMedPubMed CentralView ArticleGoogle Scholar
- Li T, Mazeas L, Sghir A, Leblon G, Bouchez T: Insights into networks of functional microbes catalysing methanization of cellulose under mesophilic conditions. Environ Microbiol. 2009, 11: 889-904. 10.1111/j.1462-2920.2008.01810.x.PubMedView ArticleGoogle Scholar
- Harrell L, Wang Y, Antonopoulos D, Young V, Lichtenstein L, Huang Y, Hanauer S, Chang E: Standard colonic lavage alters the natural state of mucosal-associated microbiota in the human colon. PLoS One. 2012, 7: e32545-10.1371/journal.pone.0032545.PubMedPubMed CentralView ArticleGoogle Scholar
- Hooper LV, Littman DR, Macpherson AJ: Interactions between the microbiota and the immune system. Science. 2012, 336: 1268-1273. 10.1126/science.1223490.PubMedPubMed CentralView ArticleGoogle Scholar
- Kostic AD, Gevers D, Pedamallu CS, Michaud M, Duke F, Earl AM, Ojesina AI, Jung J, Bass AJ, Tabernero J, et al: Genomic analysis identifies association of Fusobacterium with colorectal carcinoma. Genome research. 2012, 22: 292-298. 10.1101/gr.126573.111.PubMedPubMed CentralView ArticleGoogle Scholar
- Castellarin M, Warren RL, Freeman JD, Dreolini L, Krzywinski M, Strauss J, Barnes R, Watson P, Allen-Vercoe E, Moore RA, Holt RA: Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Genome research. 2012, 22: 299-306. 10.1101/gr.126516.111.PubMedPubMed CentralView ArticleGoogle Scholar
- Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38: D355-360. 10.1093/nar/gkp896.PubMedPubMed CentralView ArticleGoogle Scholar
- Abubucker S, Segata N, Goll J, Schubert AM, Izard J, Cantarel BL, Rodriguez-Mueller B, Zucker J, Thiagarajan M, Henrissat B, et al: Metabolic reconstruction for metagenomic data and its application to the human microbiome. PLoS Comput Biol. 2012, 8: e1002358-10.1371/journal.pcbi.1002358.PubMedPubMed CentralView ArticleGoogle Scholar
- Burke C, Steinberg P, Rusch D, Kjelleberg S, Thomas T: Bacterial community assembly based on functional genes rather than species. Proc Natl Acad Sci USA. 2011, 108: 14288-14293. 10.1073/pnas.1101591108.PubMedPubMed CentralView ArticleGoogle Scholar
- Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto JM, et al: Enterotypes of the human gut microbiome. Nature. 2011, 473: 174-180. 10.1038/nature09944.PubMedPubMed CentralView ArticleGoogle Scholar
- Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, et al: A core gut microbiome in obese and lean twins. Nature. 2009, 457: 480-484. 10.1038/nature07540.PubMedPubMed CentralView ArticleGoogle Scholar
- Sherrill C, Fahey RC: Import and metabolism of glutathione by Streptococcus mutans. J Bacteriol. 1998, 180: 1454-1459.PubMedPubMed CentralGoogle Scholar
- Keshavarzian A, Banan A, Farhadi A, Komanduri S, Mutlu E, Zhang Y, Fields JZ: Increases in free radicals and cytoskeletal protein oxidation and nitration in the colon of patients with inflammatory bowel disease. Gut. 2003, 52: 720-728. 10.1136/gut.52.5.720.PubMedPubMed CentralView ArticleGoogle Scholar
- Danese S, Sgambato A, Papa A, Scaldaferri F, Pola R, Sans M, Lovecchio M, Gasbarrini G, Cittadini A, Gasbarrini A: Homocysteine triggers mucosal microvascular activation in inflammatory bowel disease. Am J Gastroenterol. 2005, 100: 886-895. 10.1111/j.1572-0241.2005.41469.x.PubMedView ArticleGoogle Scholar
- Jha G, Rajeshwari R, Sonti RV: Bacterial type two secretion system secreted proteins: double-edged swords for plant pathogens. Mol Plant Microbe Interact. 2005, 18: 891-898. 10.1094/MPMI-18-0891.PubMedView ArticleGoogle Scholar
- Sandkvist M: Type II secretion and pathogenesis. Infect Immun. 2001, 69: 3523-3535. 10.1128/IAI.69.6.3523-3535.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Thiennimitr P, Winter SE, Winter MG, Xavier MN, Tolstikov V, Huseby DL, Sterzenbach T, Tsolis RM, Roth JR, Baumler AJ: Intestinal inflammation allows Salmonella to use ethanolamine to compete with the microbiota. Proc Natl Acad Sci USA. 2011, 108: 17480-17485. 10.1073/pnas.1107857108.PubMedPubMed CentralView ArticleGoogle Scholar
- Bertin Y, Girardeau JP, Chaucheyras-Durand F, Lyan B, Pujos-Guillot E, Harel J, Martin C: Enterohaemorrhagic Escherichia coli gains a competitive advantage by using ethanolamine as a nitrogen source in the bovine intestinal content. Environ Microbiol. 2011, 13: 365-377. 10.1111/j.1462-2920.2010.02334.x.PubMedView ArticleGoogle Scholar
- Sczesnak A, Segata N, Qin X, Gevers D, Petrosino JF, Huttenhower C, Littman DR, Ivanov II: The genome of th17 cell-inducing segmented filamentous bacteria reveals extensive auxotrophy and adaptations to the intestinal environment. Cell Host Microbe. 2011, 10: 260-272. 10.1016/j.chom.2011.08.005.PubMedPubMed CentralView ArticleGoogle Scholar
- Kuwahara T, Ogura Y, Oshima K, Kurokawa K, Ooka T, Hirakawa H, Itoh T, Nakayama-Imaohji H, Ichimura M, Itoh K, et al: The lifestyle of the segmented filamentous bacterium: a non-culturable gut-associated immunostimulating microbe inferred by whole-genome sequencing. DNA Res. 2011, 18: 291-303. 10.1093/dnares/dsr022.PubMedPubMed CentralView ArticleGoogle Scholar
- Gaboriau-Routhiau V, Rakotobe S, Lecuyer E, Mulder I, Lan A, Bridonneau C, Rochet V, Pisi A, De Paepe M, Brandi G, et al: The key role of segmented filamentous bacteria in the coordinated maturation of gut helper T cell responses. Immunity. 2009, 31: 677-689. 10.1016/j.immuni.2009.08.020.PubMedView ArticleGoogle Scholar
- Brand S: Crohn's disease: Th1, Th17 or both? The change of a paradigm: new immunological and genetic insights implicate Th17 cells in the pathogenesis of Crohn's disease. Gut. 2009, 58: 1152-1167. 10.1136/gut.2008.163667.PubMedView ArticleGoogle Scholar
- Ivanov II, Atarashi K, Manel N, Brodie EL, Shima T, Karaoz U, Wei D, Goldfarb KC, Santee CA, Lynch SV, et al: Induction of intestinal Th17 cells by segmented filamentous bacteria. Cell. 2009, 139: 485-498. 10.1016/j.cell.2009.09.033.PubMedPubMed CentralView ArticleGoogle Scholar
- Papa E, Docktor M, Smillie C, Weber S, Preheim SP, Gevers D, Giannoukos G, Ciulla D, Tabbaa D, Ingram J, et al: Non-invasive mapping of the gastrointestinal microbiota identifies children with inflammatory bowel disease. PLoS ONE. 2012, 7: e39242-10.1371/journal.pone.0039242.PubMedPubMed CentralView ArticleGoogle Scholar
- The Human Microbiome Project Data Analysis and Coordination Center: Tools and Technology. [http://www.hmpdacc.org/tools_protocols/tools_protocols.php]
- Schloss PD, Handelsman J: A statistical toolbox for metagenomics: assessing functional diversity in microbial communities. BMC Bioinformatics. 2008, 9: 34-10.1186/1471-2105-9-34.PubMedPubMed CentralView ArticleGoogle Scholar
- Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E, et al: 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-Syed-Mohideen AS, McGarrell DM, Marsh T, Garrity GM, Tiedje JM: The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009, 37: D141-145. 10.1093/nar/gkn879.PubMedPubMed CentralView ArticleGoogle Scholar
- Price MN, Dehal PS, Arkin AP: FastTree 2--approximately maximum-likelihood trees for large alignments. PLoS One. 2010, 5: e9490-10.1371/journal.pone.0009490.PubMedPubMed CentralView ArticleGoogle Scholar
- Ye Y, Doak TG: A parsimony approach to biological pathway reconstruction/inference for genomes and metagenomes. PLoS Comput Biol. 2009, 5: e1000465-10.1371/journal.pcbi.1000465.PubMedPubMed CentralView ArticleGoogle Scholar
- Chao A: Nonparametric estimation of the number of classes in a population. Scandinavian Journal of Statistics. 1984, 11: 265-270.Google Scholar
- Pielou EC: Indices of diversity and evenness. Ecological Diversity. Edited by: Pielou EC. 1975, John Wiley and SonsGoogle Scholar
- Simpson EH: Measurement of diversity. Nature. 1949, 163: 688-10.1038/163688a0.View ArticleGoogle Scholar
- Vavrek MJ: fossil: Palaeoecological and palaeogeographical analysis tools. Palaeontologia Electronica. 2011, 14:Google Scholar
- Hastie T, Tibshirani R, Friedman J: The Elements of Statistical Learning. 2009, New York, NY: SpringerView ArticleGoogle Scholar
- Friedman JH: Stochastic gradient boosting. Computational Statistics and Data Analysis. 2002, 38: 367-378. 10.1016/S0167-9473(01)00065-2.View ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Statistical Society B. 1995, 57: 289-300.Google Scholar
- Lê S, Husson F: FactorMineR: An R Package for Multivariate Analysis. Journal of Statistical Software. 2008, 25: 1-18.View ArticleGoogle Scholar
- Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C: Metagenomic biomarker discovery and explanation. Genome biology. 2011, 12: R60-10.1186/gb-2011-12-6-r60.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.