Structure of the gut microbiome following colonization with human feces determines colonic tumor burden
© Baxter et al.; licensee BioMed Central Ltd. 2014
Received: 5 March 2014
Accepted: 4 June 2014
Published: 17 June 2014
A growing body of evidence indicates that the gut microbiome plays a role in the development of colorectal cancer (CRC). Patients with CRC harbor gut microbiomes that are structurally distinct from those of healthy individuals; however, without the ability to track individuals during disease progression, it has not been possible to observe changes in the microbiome over the course of tumorigenesis. Mouse models have demonstrated that these changes can further promote colonic tumorigenesis. However, these models have relied upon mouse-adapted bacterial populations and so it remains unclear which human-adapted bacterial populations are responsible for modulating tumorigenesis.
We transplanted fecal microbiota from three CRC patients and three healthy individuals into germ-free mice, resulting in six structurally distinct microbial communities. Subjecting these mice to a chemically induced model of CRC resulted in different levels of tumorigenesis between mice. Differences in the number of tumors were strongly associated with the baseline microbiome structure in mice, but not with the cancer status of the human donors. Partitioning of baseline communities into enterotypes by Dirichlet multinomial mixture modeling resulted in three enterotypes that corresponded with tumor burden. The taxa most strongly positively correlated with increased tumor burden were members of the Bacteroides, Parabacteroides, Alistipes, and Akkermansia, all of which are Gram-negative. Members of the Gram-positive Clostridiales, including multiple members of Clostridium Group XIVa, were strongly negatively correlated with tumors. Analysis of the inferred metagenome of each community revealed a negative correlation between tumor count and the potential for butyrate production, and a positive correlation between tumor count and the capacity for host glycan degradation. Despite harboring distinct gut communities, all mice underwent conserved structural changes over the course of the model. The extent of these changes was also correlated with tumor incidence.
Our results suggest that the initial structure of the microbiome determines susceptibility to colonic tumorigenesis. There appear to be opposing roles for certain Gram-negative (Bacteroidales and Verrucomicrobia) and Gram-positive (Clostridiales) bacteria in tumor susceptibility. Thus, the impact of community structure is potentially mediated by the balance between protective, butyrate-producing populations and inflammatory, mucin-degrading populations.
KeywordsColorectal cancer Community structure Germ-free Gut microbiome Humanized mice Microbiota
Colorectal cancer (CRC) is the second leading cause of cancer-related death in the United States each year . Recent evidence suggests that the community of microbes inhabiting the gastrointestinal tract plays an important role in the development and progression of CRC [2–4]. This community, termed the gut microbiome, is known to influence cancer-related functions, including cell proliferation, angiogenesis, and apoptosis, and it is strongly linked to diet, obesity, and inflammation, which are known risk factors of CRC [5–9]. Using a mouse model of CRC, we have shown that structural changes to the microbiome occur during tumorigenesis and result in a gut microbiome with an increased tumorigenic capacity . These findings demonstrate that the gut microbiome has a causal role in the development and progression of CRC.
Several survey-based studies have shown that CRC patients harbor microbial communities that are structurally distinct from those of healthy individuals [11–15]. However, there has been no consensus among these studies as to which bacterial populations are important. In mouse models, several gut commensals have been shown to promote tumorigenesis in the colon. Both enterotoxigenic Bacteroides fragilis (ETBF) and strains of Escherichia coli that carry the pks pathogenicity island can promote tumorigenesis by the production of toxins [3, 4]. Fusobacterium nucleatum has also been shown to potentiate tumorigenesis in mouse models and cell culture experiments by stimulating inflammation via myeloid cell recruitment or activation of β-catenin signaling [2, 16]. Fusobacterium was also found to be enriched in a subset of human colon adenomas . Although there is increasing evidence that Fusobacterium is involved in CRC cases, it was detected in less than half of adenomas, which suggests that other bacterial populations are capable of potentiating tumorigenesis . In fact, it may be that CRC is a polymicrobial disease requiring combinations of these or other populations to influence tumorigenesis.
While individual bacterial species have been associated with some human CRC cases, in other cases the capacity of the microbiome to modulate tumorigenesis could be determined by the structure of the community as a whole rather than the presence or absence of individual populations [4, 17]. The potentially polymicrobial influence of the gut microbiome on this disease necessitates the disentangling of the complex interactions between bacterial populations in the gut. Understanding these interactions requires investigation of the relationship between the microbiome and tumorigenesis under a diverse set of community structures. Unfortunately, mechanistic studies typically rely on experiments with conventionally reared inbred mice living in homogenous, controlled environments, leading to relatively little variation in microbiome structure between individual animals. Although, experiments in conventional mice are useful for understanding the mechanisms by which the microbiome modulates tumorigenesis, they are limited by investigating only those bacterial strains found in laboratory mice, many of which are absent in human beings. It is reasonable to expect that incorporating human-associated microbial populations into these experiments would increase the ability to translate results to human beings.
To investigate the role of microbiome structure in tumorigenesis, we combined the advantages of the high interpersonal variation among human beings and the convenience of a mouse model. We inoculated germ-free mice with microbiota from human subjects harboring distinct microbiomes. This technique enabled us to test the effect of different baseline microbiome communities with variation beyond what is seen in conventionally reared mice. The transfer of human microbiota to germ-free mice, sometimes referred to as ‘humanization’, has been employed to study the microbiome in the context of several other diseases. In studies of diabetes, obesity, and malnutrition, colonization with human feces has been reported to recapitulate the phenotype of the human donors in the recipient mice [18–21]. Thus, in addition to searching for tumor-modulating community structures, we sought to determine whether this strategy could be used to recapitulate the tumor-promoting capacity of CRC patients’ microbiota in mice.
DNA extraction and 16S rRNA gene sequencing
Mouse fecal samples were collected throughout the experiment and frozen at -20°C. Genomic DNA from samples collected on days 0 and 73 and the human inocula were isolated using the PowerSoil-htp 96 Well Soil DNA isolation kit (MO BIO, Carlsbad, CA, USA) using an epMotion 5075 automated pipetting system. The V4 region of the 16S rRNA gene was amplified using custom barcoded primers and sequenced as described previously using an Illumina MiSeq sequencer . All fastq files and the MIMARKS spreadsheet are available online .
Sequence curation and analysis
The 16S rRNA gene sequences were curated using the mothur software package, as described previously [25, 27]. Briefly, paired end reads were assembled into contigs and aligned to the SILVA 16S rRNA sequence database . Sequences that failed to align or were flagged as possible chimeras by UCHIME were removed . Each sequence was classified using a naive Bayesian classifier trained against a 16S rRNA gene training set provided by the Ribosomal Database Project [30, 31]. Finally, sequences were grouped based on their taxonomic classification or clustered into operational taxonomic units (OTUs) based on a 97% similarity cutoff. The number of sequences in each sample was rarefied to 3,306 sequences per sample, to minimize the effects of uneven sampling. Parallel sequencing and processing of a mock community indicated that the error rate of the curated sequences was 0.085%.
The dissimilarity in community structure between samples was calculated using the ΘYC metric . The ΘYC distances between samples were used for ordination analysis by nonmetric dimensional scaling (NMDS) in two dimensions. Ten iterations were performed and the resulting ordination that had the lowest stress was used for data visualization. Dirichlet multinomial mixture models were generated to group samples into enterotypes based on the abundance of bacterial genera in each sample . To identify conserved changes that occurred over the course of the AOM/DSS model, the samples from each mouse on day 0, and the samples collected at the end of the model were grouped into ‘baseline’ and ‘endpoint’ categories, respectively. The R randomForest package was used to identify the OTUs that best distinguished between the two categories based on their importance for the classification model [34, 35].
The Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) software package was used to infer the metagenomic content of each sample, based on the taxonomy and abundance of each OTU . Although this method is limited by the number of available genomes, it has been shown to replicate metagenomes to a high degree of accuracy, especially for human-adapted bacterial communities. The weighted nearest sequenced taxon index (NSTI) for our samples was 0.056 ± 0.01. In general, NSTI values below 0.06 suggest that closely related reference genomes were available for the dataset . From the inferred metagenomes, we identified KEGG orthologs that could be used as markers for butyrate production or host glycan degradation. Because either butyrate kinase or butyryl-CoA:acetate CoA-transferase is required for butyrate production in the gut, the KEGG orthologs chosen as markers for butyrate production were K00929 (butyrate kinase (EC:188.8.131.52)), K01034 (acetate CoA-transferase α subunit (EC:184.108.40.206)), K01035 (acetate CoA-transferase β subunit (EC:220.127.116.11)) . To choose markers for glycan degradation, we found all of the KEGG orthologs annotated as sialidases, fucosidases, sulfatases, or members of the glycoside hydrolase family 18, as these classes of enzymes are necessary, and moderately specific for host glycan degradation [39, 40]. Ten such KEGG orthologs were found in the metagenomes and used as markers: K01138 (uncharacterized sulfatase (EC:3.1.6.-)), K01130 (arylsufatase (EC:18.104.22.168)), K01135 (arylsufatase B (EC:22.214.171.124)), K01137 (N-acetylglucosamine-6-sulfatase (EC:126.96.36.199)), K01134 (arylsufatase A (EC:188.8.131.52)), K01186 (sialidase-1 (EC:184.108.40.206)), K01206 (α-L-fucosidase (EC:3.2.151)), K01183 (1,4-β-poly-N-acetylglucosaminidase (EC:220.127.116.11)), K01205 (α-N-acetylglucosaminidase (EC:18.104.22.168)), and K05970 (sialate O-acetylesterase (EC:22.214.171.124)). Finally, we calculated the Spearman correlation coefficients between tumor counts and these KEGG orthologs.
Differences in tumor counts between Dirichlet multinomial mixture (DMM) partitions were examined using a Wilcoxon rank-sum test. To test whether there was a significant difference in tumor counts between groups that received healthy or cancer-associated inocula, we rank transformed the tumor counts to correct for heteroscedasticity and performed a nested analysis of variance (ANOVA). Differences in community structure were examined using analysis of molecular variance (AMOVA) in mothur .
Colonization of germ-free mice with human microbiota
Tumor incidence is linked to initial community structure
To determine which OTUs were driving this trend, we generated a biplot using the NMDS axes generated from the ΘYC distances between samples collected at the time of AOM injection (day 0; Figure 3B). Among the OTUs most strongly correlated with high tumor counts were two OTUs from the genus Bacteroides (OTUs 1 and 4). More detailed characterization of these OTUs indicated that OTU 1 was closely affiliated with B. uniformis and OTU 4 was affiliated with a mixture of Bacteroides species, including B. fragilis, B. ovatus, B. xylanisolvens, and B. thetaiotaomicron. Both of these OTUs were found in all six cohorts of mice and their initial abundances were positively correlated with tumor counts (ρ = 0.47 and 0.49, respectively; both P < 0.005; Spearman correlation). Interestingly, all samples were PCR-negative for the ETBF toxin gene, suggesting that OTU 4 was not ETBF. Other OTUs associated with high tumor counts were affiliated with the genera Parabacteroides (OTU 18) and Alistipes (OTU 19), as well as an OTU affiliated with the species Akkermansia muciniphila (OTU 11). In addition, several OTUs associated with Clostridium Group XIVa (OTUs 7, 9, 15, and 17), Clostridium Group IV (OTU 49), and unclassified members of the Lachnospiraceae (OTUs 67 and 13) were correlated with lower tumor counts. These results indicate that the relative abundance of specific OTUs in the starting community could be associated with tumor counts.
Changes in the microbiome during the AOM/DSS model
To determine the extent to which the microbiomes of each group changed over the course of the AOM/DSS model, we calculated the ΘYC distances between the communities in mice at the time of AOM injection and at the end of the experiment. Interestingly, the two groups with the highest tumor counts (H2, C3) changed very little over time (ΘYC = 0.12 and 0.14), while the microbiomes of the three groups with the lowest tumor counts (C2, H3, C1) changed substantially (ΘYC = 0.73, 0.76, 0.83) (Figure 2A). Thus, the closer the initial community of each group was to the tumor-associated endpoint community, the more tumors those mice developed.
Tumor incidence is linked to butyrate production and host glycan degradation
Our experiments suggested that Clostridiales, Bacteroidales, and Akkermansia played a role in modulating tumorigenesis. Members of the Clostridiales, especially Clostridium Group XIVa, are the predominate producers of intestinal butyrate, an important anti-inflammatory and anti-tumorigenic metabolite in the gut [38, 43, 44]. Bacteroides and Akkermansia, on the other hand, are known to break down host-derived glycans, especially mucin. Mucin degradation has been linked to intestinal inflammation and can facilitate colonization of intestinal pathogens [45–47]. To test whether the genomic potential for these metabolic activities is linked to tumor incidence, we used the PICRUSt software package to predict the metagenomic content for each sample at the time of AOM injection. Butyrate production in the gut requires either butyryl-CoA:acetate CoA-transferase or butyrate kinase . KEGG orthologs (KOs) of the α and β subunits of butyryl-CoA:acetate CoA-transferase were negatively correlated with tumor incidence (ρ < -0.35, P < 0.05). Butyrate kinase had the same trend, but the correlation was not statistically significant (ρ = -0.30, P = 0.08). Next, we identified KOs for sialidases, fucosidases, sulfatases, and N-acetylglucosaminidases, which are indicative of host glycan degradation [39, 40]. Of the ten such KOs found in our metagenomes, seven (two arylsulfatases, an uncharacterized sulfatase, α-L-fucosidase, sialate O-acetylesterase, α-N-acetylglucosaminidase, 1,4-β-poly-N-acetylglucosaminidase) were positively correlated with tumor count (ρ > 0.47, P <0.01). None of the three remaining KOs correlated with tumors. Together, these data suggest that the correlation between tumor incidence and the microbiome might be dependent on metabolic activity rather than bacterial phylogeny.
The results of this study demonstrate that the structure of the gut microbiome is important for determining susceptibility to inflammation-associated tumorigenesis. We observed strong correlations between the initial community structure of the gut microbiome and tumor multiplicity. This relationship is driven primarily by two distinct groups of bacteria. In general, we found that members of the Bacteroidales (Bacteroides, Parabacteroides, Alistipes, and Porphyromonodaceae) were associated with a higher rate of tumorigenesis, while members of the Clostridiales, especially Clostridium Group XIVa, were associated with a decreased rate of tumorigenesis. There were exceptions to this pattern, however, as a few OTUs associated with Clostridiales (OTUs associated with Roseburia, Blautia, and Subdoligranulum) were enriched in the groups with higher tumor counts (Figure 3B). However, these OTUs were less abundant (<0.7% mean abundance) than those Clostridiales that were negatively correlated with tumors (≈2% mean abundance). Therefore, the data generally support a model in which susceptibility to colonic tumorigenesis is determined by the balance between the abundance of members of Bacteroidales and Clostridiales. One limitation of this study is that we only assayed the fecal communities. While this was necessary for correlating baseline community structure with the numbers of tumors that developed, characterization of the mucosal microbiota could potentially yield additional associations with tumor burden. It is also important to note that, although we observed variation in the number of tumors within inoculum groups, we were unable to correlate these differences with any differences in their microbiomes.
Based on our predicted metagenomic analysis, the roles of Clostridiales and Bacteroidales could be dependent on specific metabolic activities. Members of Clostridium Group XIVa are the predominant producers of butyrate in the gut . Given the anti-inflammatory and anti-tumorigenic properties of intestinal butyrate, its production by members of Clostridium Group XIVa could explain the association with lower susceptibility to colon tumorigenesis [43, 44]. This hypothesis is supported by our predicted metagenomic data, which correlated the increased potential for butyrate production with decreased tumorigenesis. Bacteroides and Akkermansia were the two genera most strongly correlated with higher rates of tumorigenesis. Both are known mucin degraders, and several genes linked to mucin degradation were positively correlated with tumor incidence. Additionally, previous studies have linked mucin degradation by Bacteroides and Akkermansia with intestinal inflammation [45–47]. It is possible that an overabundance of these or other mucin degraders could undermine the integrity of the mucosal barrier, leading to increased inflammation. Such a mechanism could be an alternative to the ETBF-based model of tumorigenesis, as we were unable to detect the gene for the ETBF toxin in any of our samples. While we cannot exclude the possibility of a novel toxin in the Bacteroides populations in our experiment, the additional correlation with Akkermansia muciniphila supports a model in which inflammation is induced by mucin degradation. If further experiments confirm this model, blocking mucin degradation could be used as a therapeutic for preventing or slowing the progression of tumorigenesis.
In this study, we observed a relationship between tumor multiplicity and the extent to which the microbiome shifted over the course of the model. The gut community of mice with high tumor counts changed very little over the course of the model, while the microbiome of groups with low tumor counts changed drastically. Thus, the more similar the baseline community was to the endpoint community, the more tumors the host developed. We hypothesize that the microbiome of these mice was not significantly altered by the AOM/DSS model since it was already in a state of dysbiosis. Therefore, there was a greater exposure to a tumorigenic microbiome. Similarly, in a previous study, we colonized germ-free mice with the feces of conventional mice that had already gone through the model . These mice developed more tumors than germ-free mice colonized with feces from normal mice. Thus, in addition to needing a dysbiotic community to exacerbate tumorigenesis, the length of exposure to that community is important to tumor formation.
In contrast with earlier studies where human feces were used to colonize germ-free mice, we were unable to recapitulate the structures of the human microbiota donors, as numerous members of the donor community failed to colonize the recipients and others colonized in different abundances. For example, one of the donor communities (C1) was dominated by Fusobacterium species (58% relative abundance). Another inoculum (C3), contained F. nucleatum at 2% relative abundance . However, we did not recover any sequences from the phylum Fusobacteria in the recipient mice. We were also unable to culture it from the original human stool sample, suggesting that it might not have survived the freezing and thawing of the sample or was never alive in the stool. Colonizing germ-free mice with human feces and recovering a similar microbiome is probably an unreasonable expectation. The mouse host certainly selects for specific populations of bacteria based on its immune system and metabolic profile [48, 49]. In addition, some bacteria will only colonize after other bacteria have suitably prepared the environment. In the oral cavity, Fusobacterium only colonizes after streptococcal populations have first attached to the tooth surface . Although we did not fully recapitulate the community structure or phenotype of the human donors, colonizing mice with human fecal communities did serve as a useful tool for generating novel community structures to test the influence of specific bacterial populations on tumorigenesis. This strategy also allowed us to investigate the role of human microbiota, which should be more clinically relevant, while maintaining the tractability of a mouse model.
In this study, we found that the process of colonizing germ-free mice with human fecal communities did not recapitulate the phenotype of the human donors in this particular mouse model of CRC. Nonetheless, our findings demonstrate the importance of the initial microbiome structure in determining the rate of tumorigenesis. Furthermore, we identified several bacterial populations correlated with tumor incidence in the context of six distinct gut communities. Multiple OTUs associated with the order Bacteroidales and the species Akkermansia muciniphila were correlated with exacerbated tumorigenesis, while several OTUs associated with Clostridium Group XIVa and other Clostridiales were correlated with protection. Based on inferred metagenomes of the baseline communities, we provided evidence that the positive correlations between Akkermansia and Bacteroidales and tumor incidence could be a result of their ability to degrade mucin, and the negative correlation between the Clostridiales and tumor incidence could be due to the production of butyrate. The results are consistent with a model in which susceptibility is determined by the balance between mucin degradation and short chain fatty acid production. More studies are needed to confirm these results and to test the mechanisms by which these or other bacterial populations influence colon tumorigenesis. A better understanding of microbiome structures with a propensity to promote or inhibit tumorigenesis could lead to the development of prebiotic or probiotic therapies to prevent or slow the development and progression of CRC.
analysis of molecular variance
analysis of variance
Dirichlet multinomial mixture
dextran sulfate sodium
enterotoxigenic Bacteroides fragilis
nonmetric dimensional scaling
nearest sequenced taxon index
operational taxonomic unit
polymerase chain reaction
Phylogenetic Investigation of Communities by Reconstruction of Unobserved States.
This work was supported by grants from the National Institutes for Health to PDS (R01GM099514, R01HG005975, P30DK034933, and University of Michigan GI SPORE) and GYC (University of Michigan GI SPORE and ARRA Supplement P30CA4659-22S3). The funding agencies had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
- American Cancer Society: Colorectal Cancer Facts & Figures 2011-2013. 2013, AtlantaGoogle Scholar
- Kostic AD, Chun E, Robertson L, Glickman JN, Gallini CA, Michaud M, Clancy TE, Chung DC, Lochhead P, Hold GL, El-Omar EM, Brenner D, Fuchs CS, Meyerson M, Garrett WS: Fusobacterium nucleatum potentiates intestinal tumorigenesis and modulates the tumor-immune microenvironment. Cell Host Microbe. 2013, 14 (2): 207-215.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu S, Rhee KJ, Albesiano E, Rabizadeh S, Wu X, Yen HR, Huso DL, Brancati FL, Wick E, McAllister F, Housseau F, Pardoll DM, Sears CL: A human colonic commensal promotes colon tumorigenesis via activation of T helper type 17 T cell responses. Nat Med. 2009, 15 (9): 1016-1022.View ArticlePubMedPubMed CentralGoogle Scholar
- Arthur JC, Perez-Chanona E, Mühlbauer M, Tomkovich S, Uronis JM, Fan TJ, Campbell BJ, Abujamel T, Dogan B, Rogers AB, Rhodes JM, Stintzi A, Simpson KW, Hansen JJ, Keku TO, Fodor AA, Jobin C: Intestinal inflammation targets cancer-inducing activity of the microbiota. Science. 2012, 338 (6103): 120-123.View ArticlePubMedPubMed CentralGoogle Scholar
- Cheesman SE, Neal JT, Mittge E, Seredick BM, Guillemin K: Epithelial cell proliferation in the developing zebrafish intestine is regulated by the Wnt pathway and microbial signaling via Myd88. Proc Natl Acad Sci USA. 2011, 108 (Suppl 1): 4570-4577.View ArticlePubMedPubMed CentralGoogle Scholar
- Dolara P, Caderni G, Salvadori M, Morozzi G, Fabiani R, Cresci A, Orpianesi C, Trallori G, Russo A, Palli D: Fecal levels of short-chain fatty acids and bile acids as determinants of colonic mucosal cell proliferation in humans. Nutr Cancer. 2002, 42 (2): 186-190.View ArticlePubMedGoogle Scholar
- Stappenbeck TS, Hooper LV, Gordon JI: Developmental regulation of intestinal angiogenesis by indigenous microbes via Paneth cells. Proc Natl Acad Sci USA. 2002, 99 (24): 15451-15455.View ArticlePubMedPubMed CentralGoogle Scholar
- Huxley RR, Ansary-Moghaddam A, Clifton P, Czernichow S, Parr CL, Woodward M: The impact of dietary and lifestyle risk factors on risk of colorectal cancer: a quantitative overview of the epidemiological evidence. Int J Cancer. 2009, 125 (1): 171-180.View ArticlePubMedGoogle Scholar
- Chambers WM, Warren BF, Jewell DP, Mortensen NJM: Cancer surveillance in ulcerative colitis. Br J Surg. 2005, 92 (8): 928-936.View ArticlePubMedGoogle Scholar
- Zackular JP, Baxter NT, Iverson KD, Sadler WD, Petrosino JF, Chen GY, Schloss PD: The gut microbiome modulates colon tumorigenesis. mBio. 2013, 4 (6): e00692-e00613.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang T, Cai G, Qiu Y, Fei N, Zhang M, Pang X, Jia W, Cai S, Zhao L: Structural segregation of gut microbiota between colorectal cancer patients and healthy volunteers. ISME J. 2012, 6 (2): 320-329.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen HM, Yu YN, Wang JL, Lin YW, Kong X, Yang CQ, Yang L, Liu ZJ, Yuan YZ, Liu F, Wu JX, Zhong L, Fang DC, Zou W, Fang JY: Decreased dietary fiber intake and structural alteration of gut microbiota in patients with advanced colorectal adenoma. Am J Clin Nutr. 2013, 97 (5): 1044-1052.View ArticlePubMedGoogle Scholar
- Chen W, Liu F, Ling Z, Tong X, Xiang C: Human intestinal lumen and mucosa-associated microbiota in patients with colorectal cancer. PLoS One. 2012, 7 (6): e39743-View ArticlePubMedPubMed CentralGoogle Scholar
- Shen XJ, Rawls JF, Randall T, Burcal L, Mpande CN, Jenkins N, Jovov B, Abdo Z, Sandler RS, Keku TO: Molecular characterization of mucosal adherent bacteria and associations with colorectal adenomas. Gut Microbes. 2010, 1 (3): 138-147.View ArticlePubMedPubMed CentralGoogle Scholar
- Kostic AD, Gevers D, Pedamallu CS, Michaud M, Duke F, Earl AM, Ojesina AI, Jung J, Bass AJ, Tabernero J, Baselga J, Liu C, Shivdasani RA, Ogino S, Birren BW, Huttenhower C, Garrett WS, Meyerson M: Genomic analysis identifies association of Fusobacterium with colorectal carcinoma. Genome Res. 2012, 22 (2): 292-298.View ArticlePubMedPubMed CentralGoogle Scholar
- Rubinstein MR, Wang X, Liu W, Hao Y, Cai G, Han YW: Fusobacterium nucleatum promotes colorectal carcinogenesis by modulating E-cadherin/β-catenin signaling via its FadA adhesin. Cell Host Microbe. 2013, 14 (2): 195-206.View ArticlePubMedPubMed CentralGoogle Scholar
- Toprak NU, Yagci A, Gulluoglu BM, Akin ML, Demirkalem P, Celenk T, Soyletir G: A possible role of Bacteroides fragilis enterotoxin in the aetiology of colorectal cancer. Clin Microbiol Infec. 2006, 12 (8): 782-786.View ArticleGoogle Scholar
- Alpert C, Sczesny S, Gruhl B, Blaut M: Long-term stability of the human gut microbiota in two different rat strains. Curr Issues Mol Biol. 2008, 10 (1–2): 17-24.PubMedGoogle 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 (7122): 1027-1031.View ArticlePubMedGoogle Scholar
- Koren O, Goodrich JK, Cullender TC, Spor A, Laitinen K, Bäckhed HK, Gonzalez A, Werner JJ, Angenent LT, Knight R, Bäckhed F, Isolauri E, Salminen S, Ley RE: Host remodeling of the gut microbiome and metabolic changes during pregnancy. Cell. 2012, 150 (3): 470-480.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith MI, Yatsunenko T, Manary MJ, Trehan I, Mkakosya R, Cheng J, Kau AL, Rich SS, Concannon P, Mychaleckyj JC, Liu J, Houpt E, Li JV, Holmes E, Nicholson J, Knights D, Ursell LK, Knight R, Gordon JI: Gut microbiomes of Malawian twin pairs discordant for kwashiorkor. Science. 2013, 339 (6119): 548-554.View ArticlePubMedPubMed CentralGoogle Scholar
- Franco AA, Cheng RK, Goodman A, Sears CL: Modulation of bft expression by the Bacteroides fragilis pathogenicity island and its flanking region. Mol Microbiol. 2002, 45 (4): 1067-1077.View ArticlePubMedGoogle Scholar
- Chen GY, Shaw MH, Redondo G, Nunez G: The innate immune receptor Nod1 protects the intestine from inflammation-induced tumorigenesis. Cancer Res. 2008, 68 (24): 10060-10067.View ArticlePubMedPubMed CentralGoogle Scholar
- Tanaka T, Kohno H, Suzuki R, Yamada Y, Sugie S, Mori H: A novel inflammation-related mouse colon carcinogenesis model induced by azoxymethane and dextran sodium sulfate. Cancer Sci. 2003, 94 (11): 965-973.View ArticlePubMedGoogle Scholar
- Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD: Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013, 79 (17): 5112-5120.View ArticlePubMedPubMed CentralGoogle Scholar
- Schloss P: Data File Repository for the Manuscript, ‘Structure of the Gut Microbiome Following Colonization with Human Feces Determines Colonic Tumor Burden’ by Niel Baxter, Joseph Zackular and Colleagues.http://www.mothur.org/human_mouse_aomdss,
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009, 75 (23): 7537-7541.View ArticlePubMedPubMed CentralGoogle Scholar
- Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glockner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007, 35 (21): 7188-7196.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011, 27 (16): 2194-2200.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 (16): 5261-5267.View ArticlePubMedPubMed CentralGoogle Scholar
- RDP Classifier.http://sourceforge.net/projects/rdp-classifier/,
- Yue JC, Clayton MK: A similarity measure based on species proportions. Commun Stat Theory Methods. 2005, 34 (11): 2123-2131.View ArticleGoogle Scholar
- Holmes I, Harris K, Quince C: Dirichlet multinomial mixtures: generative models for microbial metagenomics. Plos One. 2012, 7 (2): e30126-View ArticlePubMedPubMed CentralGoogle Scholar
- Breiman L: Random forests. Mach Learn. 2001, 45 (1): 5-32.View ArticleGoogle Scholar
- Liaw A, Wiener M: Classification and regression by randomForest. R News. 2002, 2: 18-22.Google Scholar
- Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, Beiko RG, Huttenhower C: Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013, 31 (9): 814-821.View ArticlePubMedPubMed CentralGoogle Scholar
- Quality Control of PICRUST Predictions.http://picrust.github.io/picrust/tutorials/quality_control.html,
- Pryde SE, Duncan SH, Hold GL, Stewart CS, Flint HJ: The microbiology of butyrate formation in the human colon. FEMS Microbiol Lett. 2002, 217 (2): 133-139.View ArticlePubMedGoogle Scholar
- Corfield AP, Wagner SA, Clamp JR, Kriaris MS, Hoskins LC: Mucin degradation in the human colon: production of sialidase, sialate O-acetylesterase, N-acetylneuraminate lyase, arylesterase, and glycosulfatase activities by strains of fecal bacteria. Infect Immun. 1992, 60 (10): 3971-3978.PubMedPubMed CentralGoogle Scholar
- Martens EC, Chiang HC, Gordon JI: Mucosal glycan foraging enhances fitness and transmission of a saccharolytic human gut bacterial symbiont. Cell Host Microbe. 2008, 4 (5): 447-457.View ArticlePubMedPubMed CentralGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 131 (2): 479-491.PubMedPubMed CentralGoogle Scholar
- De Robertis M, Massi E, Poeta ML, Carotti S, Morini S, Cecchetelli L, Signori E, Fazio VM: The AOM/DSS murine model for the study of colon carcinogenesis: from pathways to diagnosis and therapy studies. J Carcinog. 2011, 10: 9-View ArticlePubMedGoogle 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 NFκB inhibition: implications for Crohn’s disease. Gut. 2000, 47 (3): 397-403.View ArticlePubMedPubMed CentralGoogle Scholar
- Hague A, Elder DJ, Hicks DJ, Paraskeva C: Apoptosis in colorectal tumour cells: induction by the short chain fatty acids butyrate, propionate and acetate and by the bile salt deoxycholate. Int J Cancer. 1995, 60 (3): 400-406.View ArticlePubMedGoogle Scholar
- Bloom SM, Bijanki VN, Nava GM, Sun L, Malvin NP, Donermeyer DL, Dunne WM, Allen PM, Stappenbeck TS: Commensal Bacteroides species induce colitis in host-genotype-specific fashion in a mouse model of inflammatory bowel disease. Cell Host Microbe. 2011, 9 (5): 390-403.View ArticlePubMedPubMed CentralGoogle Scholar
- Ng KM, Ferreyra JA, Higginbottom SK, Lynch JB, Kashyap PC, Gopinath S, Naidu N, Choudhury B, Weimer BC, Monack DM, Sonnenburg JL: Microbiota-liberated host sugars facilitate post-antibiotic expansion of enteric pathogens. Nature. 2013, 502 (7469): 96-99.View ArticlePubMedPubMed CentralGoogle Scholar
- Ganesh BP, Klopfleisch R, Loh G, Blaut M: Commensal Akkermansia muciniphila exacerbates gut inflammation in Salmonella Typhimurium-infected gnotobiotic mice. Plos One. 2013, 8 (9): e74963-View ArticlePubMedPubMed CentralGoogle Scholar
- Rawls JF, Mahowald MA, Ley RE, Gordon JI: Reciprocal gut microbiota transplants from zebrafish and mice to germ-free recipients reveal host habitat selection. Cell. 2006, 127 (2): 423-433.View ArticlePubMedPubMed CentralGoogle Scholar
- Chung H, Pamp SJ, Hill JA, Surana NK, Edelman SM, Troy EB, Reading NC, Villablanca EJ, Wang S, Mora JR, Umesaki Y, Mathis D, Benoist C, Relman DA, Kasper DL: Gut immune maturation depends on colonization with a host-specific microbiota. Cell. 2012, 149 (7): 1578-1593.View ArticlePubMedPubMed CentralGoogle Scholar
- Kolenbrander PE, Palmer RJ, Periasamy S, Jakubovics NS: Oral multispecies biofilm development and the key role of cell-cell distance. Nat Rev Microbiol. 2010, 8 (7): 471-480.View ArticlePubMedGoogle 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.