Skip to main content

Composition of the adult digestive tract bacterial microbiome based on seven mouth surfaces, tonsils, throat and stool samples

Abstract

Background

To understand the relationship between our bacterial microbiome and health, it is essential to define the microbiome in the absence of disease. The digestive tract includes diverse habitats and hosts the human body's greatest bacterial density. We describe the bacterial community composition of ten digestive tract sites from more than 200 normal adults enrolled in the Human Microbiome Project, and metagenomically determined metabolic potentials of four representative sites.

Results

The microbiota of these diverse habitats formed four groups based on similar community compositions: buccal mucosa, keratinized gingiva, hard palate; saliva, tongue, tonsils, throat; sub- and supra-gingival plaques; and stool. Phyla initially identified from environmental samples were detected throughout this population, primarily TM7, SR1, and Synergistetes. Genera with pathogenic members were well-represented among this disease-free cohort. Tooth-associated communities were distinct, but not entirely dissimilar, from other oral surfaces. The Porphyromonadaceae, Veillonellaceae and Lachnospiraceae families were common to all sites, but the distributions of their genera varied significantly. Most metabolic processes were distributed widely throughout the digestive tract microbiota, with variations in metagenomic abundance between body habitats. These included shifts in sugar transporter types between the supragingival plaque, other oral surfaces, and stool; hydrogen and hydrogen sulfide production were also differentially distributed.

Conclusions

The microbiomes of ten digestive tract sites separated into four types based on composition. A core set of metabolic pathways was present across these diverse digestive tract habitats. These data provide a critical baseline for future studies investigating local and systemic diseases affecting human health.

Background

The bacterial microbiome of the human digestive tract contributes to both health and disease. In health, bacteria are key components in the development of mucosal barrier function and in innate and adaptive immune responses, and they also work to suppress establishment of pathogens [1]. In disease, with breach of the mucosal barrier, commensal bacteria can become a chronic inflammatory stimulus to adjacent tissues [2, 3] as well as a source of immune perturbation in conditions such as atherosclerosis, type 2 diabetes, non-alcoholic fatty liver disease, obesity and inflammatory bowel disease [48]. It is therefore critically important to define the microbiome of healthy persons in order to detect significant variations both in disease states and in pre-clinical conditions to understand disease onset and progression.

The Human Microbiome Project (HMP) established by the National Institutes of Health aims to characterize the microbiome of a large cohort of normal adult subjects [9], providing an unprecedented survey of the microbiome. The HMP includes over 200 subjects and has collected microbiome samples from 15 to 18 body habitats per person [10]. This unique dataset permits novel studies of the human digestive tract within a large number of subjects, allows for comparisons of microbial communities between habitats, and enables the definition of distinct metabolic niches within and among individuals. Previous studies of the healthy adult digestive tract microbiota have typically included less than 20 individuals [1121] and the studies with over 100 individuals have most often focused on a single body site [2226]. The increased throughput, the improved sensitivity of assays and the improvements in next generation sequencing technologies have enabled cataloging of microbial community membership and structure [12, 19, 27] as well as the metagenomic gene pool present in each community in large numbers of samples from large numbers of subjects. The HMP in particular includes, for each sample, both 16S rRNA gene surveys and shotgun metagenomic sequences, from a subset of the subjects recruited at two geographically distant locations in the United States. The recruitment criteria included a set of objective, composite measurements performed by healthcare professionals [10], defining this reference population and enabling this investigation to focus on defining the integrated oral, oropharyngeal, and gut microbiomes in the absence of host disease.

The focus of this study, complementary to other activities in the HMP consortium, was to measure and compare the composition, relative abundance, phylogenetic and metabolic potential of the bacterial populations inhabiting multiple sites along the digestive tract in the defined adult reference HMP subject population. The digestive tract was represented by ten microbiome samples from distinct body habitats: seven samples were from the mouth (buccal mucosa, keratinized attached gingiva, hard palate, saliva, tongue and two surfaces along the tooth); two oropharyngeal sites (back wall of the oropharynx (refered to here as throat) and the palatine tonsils); and the colon (stool). In addition to their distinct anatomic locations, these sites were chosen because sampling minimally disturbed the existing micobiota and involved minimal risk to participants. Although existing data indicate that mucosa-associated communities below the pharynx may have distinct microbiomes, these sites were not included, as sampling would have required invasive procedures [16, 17, 28].

The results show that the ten body habitats examined here formed four categories of microbial community types. These four community types included taxa typically classified as 'environmental' phyla. Genera characterized by pathogenic species and thus associated with disease were also widely distributed among the population. Most striking, each body site (within as well as between the four groups) possessed a highly distinctive community structure with moderate variability across the population, and with distinct abundances of some microbial metabolic processes within each community. The combination of high-throughput sequencing technologies and a large, well-characterized study population has thus provided quantitative and qualitative outputs that allow a comprehensive definition of the normal adult digestive tract microbiome.

Results

Microbial community structure indicates four distinct community types within the ten digestive tract sites

At all phylogenetic levels, from phylum to genus, we identified four groups of body habitats that maintain a distinct pattern of numerically dominant bacterial taxa as profiled using the 16S rRNA gene (Figure 1a), as classified by the Ribosome Database Project (RDP) [29]. While only two phyla, the Firmicutes and Bacteroidetes, dominated the communities of all ten sites, their proportions and that of nearly all taxa in the sampled body habitats form groups as follows: Group 1, buccal mucosa, keratinized gingiva, and hard palate; Group 2, saliva, tongue, tonsils, and throat (back wall of oropharynx); Group 3, sub- and supra-gingival plaque; and Group 4, stool. The microbiota of Group 1 consisted mostly of Firmicutes followed in decreasing order of relative abundance by Proteobacteria, Bacteroidetes and either Actinobacteria or Fusobacteria (Figure 1a; Additional file 1). In comparison, Group 2 had a decreased relative abundance of Firmicutes and increased levels of four phyla: Bacteroidetes, Fusobacteria, Actinobacteria and TM7. Group 3, which consisted of both tooth plaque sites, had a further decrease in Firmicutes compared to Groups 1 and 2, with a marked increase in the relative abundance of Actinobacteria. Finally, the stool microbiota (Group 4) consisted of mostly Bacteroidetes (over 60%) followed by Firmicutes, with very low relative abundances of Proteobacteria and Actinobacteria, and less than 0.01% of Fusobacteria (Additional file 1).

Figure 1
figure 1

Groups detected in the sampled digestive tract microbiome sites based on similarities in microbial composition. (a) Taxonomic composition of the microbiota in the ten digestive tract body habitats investigated based on average relative abundance of 16S rRNA pyrosequencing reads assigned to phylum (upper chart) and genus (lower chart). Microbiota from the ten habitats are grouped based on the ratio of Firmicutes to Bacteroidetes as follows: Group 1 (G1), buccal mucosa (BM), keratinized gingiva (KG) and hard palate (HP); Group 2 (G2), throat (Th), palatine tonsils (PT), tongue dorsum (TD) and saliva (Sal); Group 3 (G3), supraginval (SupP) and subgingival plaques (SubP); and Group 4 (G4), stool (Stool). Labels indicate genera at average relative abundance ≥2% in at least one body site. The remaining genera were binned together in each phylum as 'other' along with the fraction of reads that could not be assigned at the genus level as 'unclassified' (uncl.). See Additional file 1 for detailed values. (b) Circular cladogram reporting taxa consistently differential among the body habitats in at least one group detected using LEfSe. Colors indicate the group in which each differential clade was most abundant. See Additional file 3 for the detailed list of taxa whose representation was statistically different among the groups. The representation is based on RDP phylogenetic hierarchy.

These dramatic differences occurred consistently throughout the cohort, with closely juxtaposed body sites (for example, tongue dorsum (Group 2) and hard palate (Group 1)) possessing strikingly different microbial community structure even when considering the phylum level alone and independently of the structure of the tissue (Additional file 2). This supports strong local selective pressure on community membership even in the absence of disease, and these differences reach to the genus level (Figure 1a). In terms of genera, Group 1 was characterized by a very high relative abundance of Streptococcus, while Group 4 contained predominantly Bacteroides. In contrast, Groups 2 and 3, rather than having a single genus present at such high relative abundance, were characterized by a more even distribution of the most abundant genera. Streptococcus, Veillonella, Prevotella, Neisseria, Fusobacterium, Actinomyces and Leptotrichia were each present over 2% on average in Group 2. These seven genera plus Corynebacterium, Capnocytophaga, Rothia and Porphyromonas accounted for genera present at more than 2% in Group 3 (Figure 1a; Additional file 1).

Examining clade abundances at all taxonomic levels, we used the LEfSe (LDA Effect Size) system for biomarker discovery [30] to determine statistically significant biomarkers among these four groups within the digestive tract. These included both high and low abundance clades that significantly and consistently varied in abundance among and within body habitats, for example, in the three oral groups (Figure 1b; Additional file 3). For example, both the phylum Actinobacteria and individual taxa within the Actinomycetales were consistently more abundant on the tooth surfaces in Group 3 (Figure 1b; Additional file 3). When comparing Group 1 against the other three groups (a slightly more stringent setting than comparing all groups against each other as in Figure 1b and Additional file 1) two genera from the Firmicutes were identified as genus-level biomarkers: Streptococcus, from the Streptococcaceae (mean 47 ± 18% abundance in Group 1), and Gemella, from the Staphylococcaceae (mean 5.2 ± 5.1% abundance in Group 1) (Additional file 1). Although the Firmicutes phylum as a whole was most differentially abundant in Group 1, more specific taxa within the Firmicutes were detected as biomarkers for Groups 2 and 4 (Figure 1b; Additional file 3). For example, in Group 2, biomarkers, when compared to the other three groups, included Oribacterium and Catonella, members of the Lachnospiraceae, and Veillonella, a member of the Veillonellaceae (all Clostridia). The abundances of Veillonella and Prevotella overall were comparable in Group 2 (10.2 ± 5.4% versus 11.6 ± 7.3%, respectively), and both were identified as differentially abundant in this group. The other genus-level biomarkers for Group 2 detected at >2% were Porphyromonas (3.8 ± 4.2%) and Neisseria (6.6 ± 7.6%) (Additional file 1). Several genus-level biomarkers for Group 4 (stool) were also Firmicutes, mostly from the families Lachnospiraceae and Ruminococcaceae (Figure 1b; Additional file 3). These results support the overall consistency of the different microbial populations characterizing each of the four groups, and they also emphasize the need to take multiple levels of phylogenetic specificity into account when performing any analysis of the microbiome. Phylum relative abundances differentiated very distinct body habitats. As additionaly discussed below, these differences were reflected at the genus level within each body site in the healthy adult human.

The four observed groups differed significantly not only based on their specific microbial compositions, but also by several ecological summary statistics. Most strikingly, after comparing every pair of samples using the Bray-Curtis measure of beta diversity [31], within-group distance was very significantly lower (greater similarity) than between-group distance (lower similarity) for all four groups (Additional file 4; Table 1; all P < 10-20). The coarse level of species richness measurement offered by phylotype data did not distinguish strongly among any body habitats, but evenness and the resulting within-community alpha diversity ranged widely among groups as measured by the inverse Simpson index [32] (Additional file 5). For example, the Group 1 body sites together averaged below a relative diversity of 5.3, Group 2 ranged from 7.3 ± 3.0 (tonsils) to 10.6 ± 3.1 (saliva), the plaques in Group 3 had average diversities of 9.6 ± 3.1 and 9.8 ± 3.0, and Group 4 (stool) declined to a mean of 4.6 ± 2.9. The lower diversities in Group 1 are largely an effect of Streptococcus abundance, and likewise the gut microbiota's diversity is lowered by the prevalence of the Bacteroides in these data (both detailed above and below). These differences are highly statistically significant (for example, Group 1 versus 2 P < 1e-50 by t-test) and provide evidence in support of the four-group distinction at the levels of both individual bacterial clade and overall ecological structure.

Table 1 Community structure similarity is higher for samples in the same digestive tract group than for samples in different groups or outside the digestive tract

Phyla typically identified with environmental communities are part of the natural microbiota of healthy humans

Bacterial phyla originally thought to be exclusively environmental have recently been observed to possess human host-associated membership [3336]. This phenomenon was widely observed within this normal population. The phylum TM7 was highly prevalent, detected in at least one sampling site of the upper digestive tract of 85% of subjects and in the stool of 13.6% of the subjects (Additional file 6). The phyla SR1 and Synergistetes were present in at least one upper digestive tract site of 65.3% and 58.5% of the subjects and in the stool of 1.4% and 8.8% of the subjects, respectively. The phylum Verrucomicrobia, represented mainly by the genus Akkermansia [35], and the phylum Lentisphaerae, represented by the genus Victivallis [34], were present in the lower digestive tract of 41.5% and 15.0% of the subjects and in the upper digestive tract of 13.6% and 1.4% of the subjects. TM7 bacteria accounted for a mean of 3.1 ± 5.7% of the saliva population and 0.6 ± 1.2% of the bacteria found in plaque communities (Figures 1a and 2; Additional file 1). The SR1 phylum was also most abundant in saliva (mean 0.4 ± 1.2%), and both TM7 and SR1 phyla were found at trace amounts in stool. While these phyla were varyingly prevalent (Figure 2), they occurred near-uniformly at low but significantly non-zero abundances, which highlights their lack of detection in smaller studies without deep high-throughput sequencing.

Figure 2
figure 2

Noticeable relative abundance and variability of TM7, Synergistetes, and SR1 per body habitat. Representation of the relative abundances of the phyla TM7, Synergistetes (Synerg.), and SR1 among the subject population, expressed as percentage on a log scale (left). The high relative abundances of members of these phyla among the subjects, in particular for TM7, indicate a potential role in eubiosis. The body habitats and groups are labeled as in Figure 1.

Genera characterized by pathogenic members and thus associated with disease were prevalent at low abundance in the normal human microbiota

Clades populated with known bacterial oral pathogens were well represented in this reference adult cohort, typically with moderate to high population penetrance but low relative abundance in each individual. Among the Spirochetes, Treponema species are associated with periodontal and endodontic diseases [37, 38] and were present in at least one of the upper digestive tract sites of 96% of this disease-free population (and in all the nine oral sites of 6.8%). Treponema had a variable relative abundance among the oral body habitats, with highest representation in the subgingival biofilm (mean 2.2 ± 4.1%) and non-zero abundances in several other sites, for example, palatine tonsils (Figure 3; Additional file 1). In contrast, a minority of stool samples (3.4%) contained trace levels of Treponema. The previously published rarity and specificity of Brachyspira to the gut was confirmed by its detectable presence in only one stool sample (226 stool samples in total; Additional file 7) and absence from all the upper digestive tract sites (1,879 samples; Additional file 7). Other periodontal pathogens were lower in abundance. Aggregatibacter were found mostly along the tooth surfaces (Group 3; mean 0.4 ± 0.7/0.8% from supra- and sub-gingival biofilms), and Megasphaera were found mostly in Group 2 (from mean 0.4 ± 0.6% in the tonsils and tongue dorsum to 0.8 ± 0.9% in saliva). Bifidobacteriaceae, implicated in the formation of caries [39, 40], were very rare at all oral sites (means <0.03%), but possessed high prevalence (40.8%). In the stool, the genus Bifidobacterium was most represented with a low mean relative abundance of 0.08 ± 0.3%. The low abundance of Bifidobacteriaceae in the oral cavity may be a reflection of the lack of carious lesions in this healthy subject population. Porphyromonas, which includes Porphyromonas gingivalis (one of the most studied oral pathogens) and non-pathogeneic strains, was present in the upper digestive tract of all the subjects (mean 3.0 ± 3.8%, 3.8 ± 4.2%, and 3.0 ± 3.5% in the three oral groups, respectively) and in 25% of the lower digestive tract samples, though in very low abundance in the stool (Additional file 1). Tannerella, thought to incur similar host phenotypes, was present in the upper digestive tract of 97.3% and in the stool of 3.4% of the subjects. Both genera, Porphyromonas and Tannerella, were almost uniquely distributed in average abundance among individual body sites within the oral cavity, whereas the other relevant genera in the family Porphyromonadaceae (Parabacteroides, Barnesiella, Odoribacter, and Butyricimonas) predominantly colonize the stool (Figure 4).

Figure 3
figure 3

Most microbes in the digestive tract communities vary widely in relative abundance among body habitats and individuals. Genera with the lowest (top) to highest (bottom) variability among samples spanning all ten body sites, with coefficients of variation reported numerically (right column) and relative abundance colored on a log scale. The scale bar shows the color-coding of the average relative abundance expressed as percentage, from low (black) to high (red). All genera present >0.001% in at least half of the samples are reported. Prevotella, Veillonella, and Streptococcus are least variable across both body sites and individuals.

Figure 4
figure 4

Genera within the Porphyromonadaceae, Veillonellaceae and Lachnospiraceae families are differentially abundant across microbial communities between the upper and lower digestive tract. These three families were detected among all ten digestive body habitats, but genera within them showed varying patterns of niche specialization to sites along the digestive tract. All genera with at least 0.001% abundance in at least one body site are reported here. Clades showing a statistically significant difference (by LEfSe) specifically between oral and stool samples are indiocated with asterisks. Abundances are reported on a log scale as averages. The scale bar shows the color-coding of the average relative abundance expressed as percentage, from low (black) to high (red). The Porphyromonadaceae family is interesting in that its average abundances are higher in the gut than in the oral body habitats, but specific genera within the family diverge: Tannerella and Porphyromonas are predominantly present in the oral cavity, whereas Parabacteroides, Barnesiella, Odoribacter and Butyricimonas show higher relative abundances in the gut. BM, buccal mucosa; KG, keratinized gingiva; HP, hard palate; Th, throat; PT, palatine tonsils; TD, tongue dorsum; Sal, saliva; SupP, supraginval; SubP, subgingival plaques.

Genera that include important human pathogens colonizing the throat/tonsils - Streptococcus pneumoniae, Streptococcus pyogenes, Neisseria meningitidis, and Haemophilus influenzae - were all well represented in the microbiota of the upper digestive tract sites (Figure 1a; Additional file 6). The known difficulty of performing species-level identification from 16S rRNA pyrosequencing experiments [41] precluded the determination of prevalence for these specific pathogens in this cohort. The genus Moraxella, which includes the common sinus pathogen Moraxella catarrhalis, was detected in the upper digestive tract microbiota at low relative abundance, reaching a mean >0.5% only in the throat (Additional file 1). Interestingly, the high standard deviation (4.7%) of the relative abundances of Moraxella in the throat suggested variable colonization within this population.

In the lower intestinal tract, genera containing known pathogens were typically low in both prevalence and relative abundance. Helicobacter, implicated in gastrointestinal diseases, appeared in only 1.4% of stool samples in trace amounts while studies of Helicobacter pylori stool antigen prevalence in healthy European adults ranged up to 33% [42]. Enterobacteriaceae abundances were uneven among individuals in the gut and within each individual among body sites, with the most abundant genus being the Escherichia/Shigella complex (mean 0.1 ± 0.67%), which was detected in 33% of stool samples. Finally, Faecalibacterium, a genus of considerable interest due to its observed decrease in abundance in active Crohn's disease [4346], was highly represented in the stool (98% of subjects and mean 4.6 ± 5.2%) but present only at trace levels in the oral cavity (always below 0.05%), suggesting that it may be adapted to a very specific niche within the gut.

Comparison of microbial communities from the two tooth surface-associated sites

Within the oral cavity, the Group 3 sub- and supra-gingival plaque bacterial communities were most distinct and differed strongly from the other body sites, but further differences characterized each of these two sites individually. The tooth surface adjacent to the soft gingival tissues specifically comprises two distinct ecological niches, supragingival, and subgingival (Additional file 8). The supragingival region sits above the gingival margin, exposed to the oral cavity, bathed in saliva and exposed to ingested substances; the subgingival region is bathed in a serum transudate that flows from the base of the crevice outward to the oral cavity. A key known physiological difference between these two regions is the lower redox potential found subgingivally [47]. Correspondingly, we observed differences in the non-diseased plaque biofilm communities from these two regions distinguished by proportional shifts consistent with these physiological distinctions (Figure 5a; Additional file 1). Shifts at the phylum level were driven by subgingival increases in the obligate anaerobic genera Fusobacterium, Prevotella, and Treponema, and by lesser relative abundances of Dialister, Eubacterium, Selenomonas, and Parvimonas. In contrast, groups significantly increased in the supragingival plaque consisted predominantly of facultative anaerobic genera, including Streptococcus, Capnocytophaga, Neisseria, Haemophilus, Leptotrichia, Actinomyces, Rothia, Corynebacterium, and Kingella (Figure 1a; Additional file 1). This suggests that along these tooth surfaces, where direct bacterial interaction with host cells is diminished, oxygen availability - an environmental factor - may be a major driver of community composition.

Figure 5
figure 5

Niche specialization is widespread throughout the digestive tract even among adjacent body habitats. (a) Circular cladogram based on the RDP Taxonomy [29] reporting taxa significantly more abundant in supragingival (red) and subgingival plaque (green) and demonstrating the extensive specialization even at these highly related sites. At the class level, Actinobacteria, Bacilli, Gamma-proteobacteria, Beta-proteobacteria, and Flavobacteria are characteristic of the supragingival plaque, whereas Fusobacteria, Clostridia, Epsilon-proteobacteria, Spirochaetes, Bacteroidia, and unclassified Bacteroidetes are biomarkers for the subgingival plaque. (b) Circular cladogram comparing the digestive tract (red, GI) with non-mucosal body habitats (green, NON-GI: comprising samples from the anterior nares, and from the bilateral skin sites, antecubital fossae, and retroauricular creases). Only a few clades are detected as differentially present and abundant throughout the entire digestive tract, as the high degree of specialization and community variability at each body site prevents any individual community member from being representative of all ten body habitats. BM, buccal mucosa; TD, tongue dorsum; SupP, supraginval.

The oropharyngeal microbiota lacked abundant site-specific bacteria across all samples when compared to the mouth

The pharynx is the site of carriage of a number of important bacterial pathogens that impact both healthy and immunocompromised individuals. LEfSe analysis of all samples did not identify any genus-level organisms characteristic of the microbiome of the throat and/or tonsils consistently present above our limit of detection. For example, when throat and tonsil samples were compared to the mouth sites, the genera Butyrivibrio and Mogibacterium (both from the phylum Firmicutes) were identified as differentially abundant, but both were present at only low levels (mean 0.057 ± 0.09%, and 0.188 ± 0.316%, respectively, corresponding to only approximately 1 to 5 sequences/sample; Additional file 1). The palatine tonsils, located in the oropharynx, are unique among the sites sampled in this study as the only lymphoid tissue. However, the genus-level tonsil-specific biomarker when compared to the mouth, Peptococcus, was again present at very low relative abundance (mean 0.049 ± 0.079%; Additional file 1). This lack of throat- or tonsil-specific biomarkers among bacterial taxa with a relative abundance >1% likely reflects the similarity of the microbiome of these two oropharyngeal sites with those of the tongue dorsum and saliva (Group 2 in Figure 1) despite their differences in tissue type (Additional file 2). This observation is supported by the comparison of the complete Group 2 with all other groups, which revealed distinct and abundant biomarkers as discussed above (Figure 1b; Additional file 3). Microbiota composition and the path of swallowed saliva suggest a potential role of saliva as one of the host factors influencing microbiota of Group 2.

No genus-level microbial biomarkers characterize the entire digestive tract microbiota as contrasted with non-mucosal body habitats

After analyzing the microbiota of body habitats within the digestive tract, we next asked if there were bacteria whose differential abundance characterized the digestive tract as a whole. The non-mucosal sites sampled in the HMP included anterior nares, both post-auricular creases (crease behind the ear), and both antecubital fossae (inner elbow crease). Propionibacterium, Staphylococcus, and Pseudomonaswere identified as biomarkers for the non-mucosal sites, based on a LEfSe analysis of all ten digestive tract sites versus the non-mucosal sites (Figure 5b). However, no genus-level biomarkers were identified as uniquely present throughout the digestive tract microbiota. The unclassified Veillonellaceae and Porphyromonadaceae (Figure 5b) are unlikely to be true biomarkers due to their low representation. Further analysis was impaired by the lack of reference sequences for them within RDP. Members of Veillonellaceae and Porphyromonadaceae families were much less abundant at non-mucosal sites, and were essentially absent from the HMP vaginal samples, suggesting that their adaptation is to the digestive tract mucosa rather than mucosal surfaces in general .

Bacterial families common throughout the digestive tract possess variable distributions of genera distinct to upper and lower sites

Bacterial genera membership overlap in the same subject between oral and stool samples was limited when considered if present in at least 45% of the subjects. It included Bacteroides, Faecalibacterium, Parabacteroides, Eubacterium, Alistipes, Dialister, Streptococcus, Prevotella, Roseburia, Coprococcus, Veillonella, Oscilibacter, and yet-to-be-classified genera from a subset of families (Additional file 6). Interestingly, the presence of genera in a large portion of subjects was not related to a stable relative abundance in the microbial communities, as Bacteroides and Dialaster were among the four most variable genera among subjects. In contrast, Prevotella, Veillonella and Streptococcus were the genera with the most consistent presence in the subject population (Figure 3). The importance of Lachnospiraceae, Veillonellaceae, and Porphyromonadaceae families in the healthy digestive tract microbiome was indicated by their relative abundance among all body habitats and among subjects (Figures 1 and 4; Additional files 1 and 6). Bacteria of the Lachnospiraceae and Veillonellaceae families specifically were present in all subjects' oral cavities and stools (Additional file 6). Porphyromonadaceae were present in the oral cavity of all subjects and the stool of 95.9% of subjects (Additional file 6), although their relative abundance of member genera varied by body habitat (Figure 4). Porphyromonas was present primarily in the oral sites, while Parabacteroides, Barnesiella, Odoribacter and Butyricinomonas were predominant in the stool (Figure 4). The significance of this variation in genus distribution was confirmed by LEfSe conmparisons of the upper (oral) and lower (stool) digestive tract sites (Additional file 3). In contrast, Tannerella (Figure 4) was present in most oral sites, but due to a lower relative abundance specifically in the keratinized gingiva, it was not found to be statistically signicant between the oral and gut sites. The pattern of variable genus distribution between the upper and lower parts of the digestive tract holds for the Lachnospiraceae and Veillonellaceae as well (Figure 4), again suggesting a pattern of niche specialization among human body habitats extending from the bacterial family level down to specific genera.

Differential representation of microbial metabolic function among body sites using reconstruction from whole shotgun sequencing

In addition to relative abundances of bacterial organisms based on 16S rRNA genes, we examined the abundances of microbial metabolic pathways as profiled from metagenomic shotgun sequencing of a subset of the available body habitats [48]. These data from the HMP included one body site within each of the four digestive tract groups: the buccal mucosa (Group 1), the tongue dorsum (Group 2), the supragingival plaque (Group 3) and the stool (Group 4). The data analyzed below include the relative abundances of individual enzyme families (Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthologous groups (KOs) [49]) and of complete metabolic modules (KMods) (Figure 6; Additional file 9).

Figure 6
figure 6

Functional characterization of the digestive microbiota based on metabolic pathway abundances in the buccal mucosa, supragingival plaque, tongue dorsum, and stool from metagenomic shotgun sequencing. Cladogram represents the KEGG BRITE functional hierarchy, with the outermost circles representing individual metabolic modules and the innermost very broad functional categories. Pathways coloration denotes modules showing significant differential abundances in at least one of the four body habitats. Metabolic profiling was performed with HUMAnN [48], revealing a much lower degree of variability among individuals and significant specifity of many pathways' relative abundance to individual body habitats. In particular, sugar transport and metabolism varies at each of the four habitats with metagenomic data, as does iron uptake and utilization.

Bacterial cells use a wide variety of aerobic or anaerobic degradation pathways as energy sources, and this was most evident in the differences in relative abundance of specific sugar transporters when comparing the oral sites to the gut. PTS transporters for small sugars were most abundant in the oral cavity and were represented for monosaccharides by mannose (M00276) and fructose (M00304) transporters, as well as the transporter of galactosamine (M00287), derived from the breakdown of sugar-decorated glycoproteins. The supragingival plaque microbiome was enriched for threhalose (M00270, M00204), alpha-glucosides (M00201, M00200), and cellobiose (M00206) transport; in contrast, the stool microbiome was enriched for the transport of lactose/arabinose (M00199) and oligogalacturonide (M00202), and for the degradation of the larger dermatan (M00076), chondroitin (M00077) and heparin (M00078) sulfate polysaccharides. Surprisingly, while anaerobiosis-related pathways were expected throughout the digestive tract, putrescine transporters in particular were most represented in the oral cavity (M00193, M00300). This is of potential interest as concurrent production and import of putrescine is a delicate balance, and excess putrescine is one source of halitosis [50].

Consistent with what is known about the function of the colonic gut bacteria, we observed several prominent enzymes and metabolic pathways most abundantly in the stool metagenome. For instance, β-glucosidase (K05349) was specifically abundant in the gut microbiota and not at oral sites; this enzyme is critical in the pathway of cellulose breakdown to β-D-glucose. Concomitantly, given that the Embden-Meyerhoff pathway is also known to be the major route for glucose metabolism to pyruvate in the colon, the highly associated glycolysis pathway module (M00001) was also significantly enriched in stool [51]. This finding is further in agreement with the 16S rRNA gene sequencing data, which included prevalent Ruminococcus in stool that are important colonizers of plant-derived material in the gut and possess cellulolytic activity [52]. The stool bacteria were also uniquely associated with pathways for ammonia (M00028, urea cycle, and M00029, ornithine biosynthesis) and methane (M00356 and M00357, both methanogenesis) production; the prominence of these enzymes is consonant with the colonic microbiome as a significant source of ammonia production. In fact, targeting the colonic microbiome with antibiotics has been shown to be a successful therapy in acquired diseases of hyperammonemia such as encephalopathy complicating hepatic insufficiency [53]. Relatedly, compared to upper digestive tract sites, there was very high abundance of a specific multiple antibiotic resistance protein (K05595) and association with the pyruvate:ferredoxin oxidoreductase pathway, which, due to its role in converson of metronidazole to its active nitroso form, can also determine sensitivity to this antibiotic. These potential pathogenically linked behaviors are of course in addition to the expected colonic bacterial activities detected for producing energy from undigested cellulose, nitrogen-containing compounds, and vitamins and cofactors important in support of basic metabolic pathways.

Although HMP protocols were optimized for bacterial sequences, shotgun sequencing also provides an initial means of assessing the community structure of non-16S assayable microbes. As reported in Additional file 10, the fractions of Archaea (0.04% in the stool; below the detection threshold in the oral cavity) and small eukaryotes (0.34% in the buccal mucosal; <0.1% in the other body sites) detected here proved to be very small. Although this may be due in part to the HMP's specific DNA handling protocols [10], this suggests that 16S rRNA gene-based community surveys provide an accurate overview of these digestive-tract associated microbial communities. Likewise, ribosomal and shotgun sequencing in the HMP cohort have been compared elsewhere and provide consistent quantitative estimates of genus-level abundances [54] without systematic phylum-level biases.

Integration of gene/pathway abundances from metagenomic data and bacterial clades based on 16S rRNA gene data

A subset of the HMP's microbiome samples was assayed with both shotgun metagenomic and 16S rRNA gene sequencing. This allowed us to assess co-variation of microbial abundances with inferred metabolic pathways. Strong correlations between the abundances of bacterial clades (from 16S rRNA data) and gene or pathway abundances (from metagenomic data) in some cases clearly highlighted genes carried by these organisms, and in others denoted less clear pangenomic elements or metabolic dependencies. An example of the former was the arabinofuranosyltransferase genes aftA and aftB (K13686 and K13687). These genes were only present in the tooth surface habitats and are known to be encoded by Corynebacterium, a biomarker of the plaques as discussed above. This was confirmed by the genes' strong association with Corynebacterium clades in these data (Spearman correlation 0.76, P-value <1e-15; Additional file 11). Archaea were not included in our analysis, as these were not detectable by 16S rRNA gene sequencing (due to lack of the conserved sequence in the primers used) and were poorly represented in the shotgun sequencing data.

The acquisition and export of metals for bacterial homeostasis and for pathogenicity is ubiquitous throughout the human microbiota, with iron being most generally necessary. Iron transporters were widely distributed among the microbiota of all four body sites, but again, the specific mechanisms of iron uptake and sequestration differed as needed for niche specialization. One use of the iron is its incorporation in porphyrin, and there was a wide distribution of cytochrome c heme lyase (K01764), which appeared to be ubiquitous and was not strongly associated with individual organisms. Conversely, uroporphyrinogen synthase (K01719) occurred at higher relative abundance in stool, inversely associated with members of the Clostridiales (Spearman correlation -0.79, P-value <1e-15; Additional file 11). This can be contrasted to protoporphyrinogen oxidase (K00231) in the oral cavity, which is potentially linked to the Prevotella enrichment (Spearman correlation 0.71, P-value <1e-15). Within the oral cavity specifically, coproporphyrinogen oxidase (K00228) and protoporphyrinogen oxidase (K00231) were both more abundant on the tongue and in supragingival plaque than on the buccal mucosa, expected to be linked to the increased relative abundance of Porphyromonas and Prevotella on those surfaces [55] (Figure 1; Additional file 11).

Metal export and utilization were likewise ubiquitous throughout the microbiota, but differed in the prevalence of specific mechanisms. Most genes encoding exporters needed for heme tolerance [56], such as MtrCDE (K00579, K00580) and HrtAB (K09814, K09813), were present at low levels throughout the digestive tract, although MtrCDE was somewhat enriched in the more anaerobic habitats, stool and plaques. None were significantly associated with specific organisms in these data. The gene encoding hemerithryn (K07216) was detected at multiple body sites but was highly enriched in stool. This enzyme for iron utilization is most often found in members of the Methylococcaceae family [57], but these were again not detectable in this study due to their absence from the RDP 16S rRNA database (see Materials and methods). Intriguingly, the hemerithryn (K07216) gene consistently associated with members of the Clostridiales when present in the gut (Spearman correlation 0.72, P-value <1e-15; Additional file 11). Finally, other metals, including copper and zinc, are also both necessary co-factors and potential toxins, and remediation pathways and transporters for both were observed consistently (copper resistance K07245; copper homeostasis K06201, K06079; zinc resistance K07803; and also many other metal transporters).

Although recent work has provided extensive insights into the mechanisms of bacterial interaction with the host immune system in the gut, much less is known about the relationship of the microbiota with host immunity for other body habitats and cell types. Two pathways observed in both the upper and lower digestive tracts and known to be involved in immunomodulation were hydrogen (H2) and hydrogen sulfide (H2S) production. Hydrogen production has been shown to be an important byproduct of acetogenic bacteria and also has an anti-inflammatory activity [58]. Enzymes both for utilization (for example, CoM methyltransferase, K14082) and for production (for example, hydrogenase-4, K12136) of hydrogen were identified specifically in the oral cavity (nearly completely absent from the gut), with potential bacterial contributors including Veillonella and Selenomonas species genomically [59] and, in one of the strongest links between genes and organisms in these data, an unclassified Pasteurellaceae clade in the oral cavity (Spearman correlation for K12136 >0.78, P-value <1e-15 in supragingival placque and tongue dorsum).

Hydrogen sulfide gas is involved in regulation of the host response at low concentrations and in host-cell toxicity and inhibition of short chain fatty acid production, specifically in the colon, at high concentrations [6065]. H2S may thus serve different purposes among the distinct bacterial communities of the digestive tract. The potential for its production was particularly enriched in stool (for example, by cystathione-beta-lyase, K14155), and somewhat enriched in the more anaerobic habitats, stool and plaques (for example, by methionine-gamma-lyase, K01761). A possible role in host-cell toxicity was strongly suggested by K01761's close association with the Treponema and Fusobacterium genera in plaque (Spearman correlation 0.74 and 0.82, respectively, P-values <1e-15), both of which include members specifically associated with periodontal disease (Additional file 11). These genes were again, however, present at low levels among all body sites analyzed here, consistent with a low-level immunomodulatory role for H2S throughout the digestive tract.

Discussion

The large reference population of the HMP has provided, to our knowledge, the first opportunity for a comprehensive description of the human gastrointestinal microbiota, focused here on the bacterial composition and function of ten independently sampled body habitats throughout the digestive tract. Using taxonomically binned 16S rRNA gene sequences, we identified the representation and relative abundance of organisms in 2,105 samples. We used the LEfSe system for metagenomic biomarker discovery to identify clades at all taxonomic levels whose distribution varied among four classes of body habitats, and which included rare clades not expected as commensals in the human microbiome. We also observed prevalent but low abundance of genera characterized by common pathogenic species, even in this asymptomatic reference population. Finally, we performed a complementary analysis of the metabolic modules and enzymes detected in a subset of these body sites, revealing strong variation in sugar and metal utilization among the digestive tract communities.

Four distinct groups were delineated among the microbial communities from the digestive tract sites. The groups were rooted in the ratio of the relative abundances of the two major phyla, Firmicutes and Bacteroidetes (Figure 1a), and the differences extended to the genus level. In the absence of disease, these groupings suggest that it might be possible to sample one representative site from each group in future studies as a strategy to decrease sequencing costs. For example, the buccal mucosa (Group 1), tongue dorsum (Group 2), supragingival plaque (Group 3) and stool (Group 4) could be used to represent all ten sites examined here. Samples from the suggested body habitats can be obtained with minimal discomfort and risk to participants, and are likely to provide the biomass needed to yield sufficient DNA for community whole genome shotgun analysis. Since the current study includes only healthy subjects, however, additional validation would be required to investigate pre-disease and disease states at targeted sites for both local and systemic diseases.

The oral microbiome as revealed in this investigation was generally consistent with earlier studies [11, 13, 14, 22, 66, 67]. Firmicutes largely dominated the microbial communities on oral tissue surfaces and in saliva. Dental plaque taxa were more evenly distributed, dominated by Firmicutes, Bacteriodetes, Actinobacteria, Proteobacteria and Fusobacteria. The differences in the plaque communities relative to oral tissue sites are likely driven by the ability of the microbial community to accumulate on the non-shedding tooth surface and the physiological status relative to oxygen distribution in the resulting biofilm. Porphyromonas, Tannerella and Treponema, genera consisting of recognized pathogens in periodontal diseases, were highly prevalent. The presence of these genera in greater than 95% of individuals in this non-diseased population provides strong evidence that they are part of the commensal oral microbiome. These data suggest, rather than a complete absence of pathogenic organisms from the normal microbiota, the possibility of low-level carriage of potential pathogens [6870].

The stool microbiota was distinguished from the microbiota of the upper digestive tract sites (Figure 1a), as expected, and set apart by a high abundance of Bacteroidetes. A notable difference in the composition of the stool microbiome of the HMP dataset compared to existing 16S rRNA gene profiles is the increased ratio of Bacteroidetes (>60% of the sequences) to Firmicutes (≤30% of the sequences). Many previous studies of adult American populations have observed the reverse, a preponderance of Firmicutes [15, 7173], and similar observations have been reported in geographically diverse populations [74, 75] and in infant gut microbiome colonization investigations [76]. It should be noted that all HMP gut communities were assayed from stool samples, which may differ extensively from colonic biopsies. For example, using endoscopic biopsies from just two subjects, Wang et al. [77] reported 49% of 16S rRNA gene clones were from the Firmicutes and 27.7% were from Bacteroidetes. However, even this distinction is unclear, as a study of 16S rRNA sequences from regional gut biopsies and spontaneously passed stool involving three subjects similarly showed the majority of phylotypes belonged to Firmicutes (76%) compared to 16% for Bacteroidetes [15]. In a study of stool from 154 adult women (twins and their mothers), Firmicutes had a mean relative abundance of >60% using several different methods to assess the 16S rRNA gene content of stool [24]. Finally, a recently published study of fecal microbiota in 161 older subjects (≥65 years) corroborate our findings, namely a Bacteroidetes-dominant distribution (57%) compared to Firmicutes (40%) [26]. The difference in the Firmicutes:Bacteroides ratio in stool samples analyzed by 16S rRNA composition was confirmed by whole genome shotgun data from the same samples in the HMP dataset [54]. While it is possible that these differences are linked to any of geographic location, host genetics, or differences in technical procedures, further study will be critical in explaining these apparently dramatic variations in gut microbiota composition in adults.

An estimated 1011 bacterial cells per day flow from the mouth to the stomach [78, 79]. Both cultivation and molecular techniques demonstrate an overlap in the oral, pharyngeal, esophageal and intestinal microbiomes [12, 27, 28, 75, 8085]. It has thus been hypothesized that the oral microbiota might significantly contribute to distal digestive tract populations. Among HMP subjects, the genera Bacteroides, Faecalibacterium, Parabacteroides, Eubacterium, Alistipes, Dialister, Streptococcus, Prevotella, Roseburia, Coprococcus, Veillonella, and Oscilibacter were detected in both the oral cavity and stool in more than 45% of subjects. However, the short sequence reads did not permit species-level identification, leaving open both the possibility that there are distinct distributions of species of these common genera along the digestive tract, and the question of whether oral microbes seed distal sites below the stomach.

Based on the commonality of genera detected in the upper digestive tract, we postulate that saliva, via its impact on pH (as a buffer) and nutrient availability (high mucin content) [86], is a key driver of microbial composition in the habitats above the stomach. The epithelium is likely another key driver as most of the upper gastrointestinal mucosal surfaces share a common epithelial lining (nonkeratinized, stratified, squamous epithelium), with the exception of the keratinized gingiva, hard palate and parts of the tongue dorsum, which instead share a keratinized, stratified, squamous epithelium (Additional file 2). The upper digestive tract sites are also constantly exposed to both inhaled and ingested microbes. A substantial portion of the variability observed in the upper digestive tract tract microbiota might then be explained by interactions between the saliva, host cell type, and exogenous factors such as oxygen availability and oral intake.

In contrast to these potentially homogenizing effects, the throat, among the nine upper digestive tract sites sampled, is uniquely the recipient of small particles, including microbes, that are trapped in mucus and propelled by respiratory cilia up from the trachea and down from the nasal cavity en route to the stomach. This might impose an additional selective pressure on pharyngeal microbiota. However, no such effect was evident in the oropharynx, which segregated nicely into Group 2 with sites not exposed to the constant flow of respiratory tract mucus. Group 2, with the tongue, tonsils, throat and saliva, is a reminder of the important overlap between the upper segments of the digestive and respiratory tracts: the aerodigestive tract, which consists of the 'lips, mouth, tongue, nose, throat, vocal cords, and part of the esophagus and windpipe' [87]. Evidence suggests that the pool of microbes from Group 2, and other oral sites, contribute to colonization of the airways in disease. A few examples of this from the polymicrobial airway infections of cystic fibrosis follow: one of the earliest cystic fibrosis pulmonary pathogens is Haemophilus influenzae, a common colonizer of the upper aerodigestive tract [88]; members of the Streptococcus milleri group were recently implicated as cystic fibrosis pathogens [89], and are known colonizers of the oral cavity; and lastly, members of the oropharyngeal microbiome might modulate the virulence of the key cystic fibrosis pathogen Pseudomonas [90]. To explain microbial community structure throughout the aerodigestive tract and airways, one might speculatively extend the basic argument above, noting that the counterpart of saliva is mucus in regions not bathed by its flow, including sites sampled by the HMP but not investigated here (for example, the anterior nares) and habitats that require more invasive methods for sampling (for example, nasal cavity, nasopharynx, esophagus and airways).

Several 'environmental' phyla observed in human microbiota [33, 91] appear to be strongly host-associated in this study. The Synergistetes phylum, for example, has only recently been described in detailed association with the human oral cavity [36, 92], and is still considered potentially environmental due to its common occurrence in, for example, bioreactors [93, 94]. Although completely absent from all ten sites in many individuals, it conversely comprised up to 10% of the community in some samples, and tended to recur at multiple body habitats within the same individual. This property - a dichotomy of apparent niches that includes specific and potentially stable occupation of human microbiome sites - can now be extended to TM7 and SR1 based on the HMP oral cavity data. As sequencing costs drop, deeper shotgun sequencing will provide access to such organisms with higher confidence, as most of those organisms are only known through their phylogenetically conserved genes.

Conclusions

Analysis of the HMP dataset described here has provided a comprehensive characterization of the disease-free digestive tract microbiome, and will further serve as a foundation for the study of comparable disease-associated microbial communities. By surveying the HMP population, these results can be further integrated into other currently ongoing studies of the cohort's core microbiome [9] or enterotype structure [25], if any. The personalized nature of the digestive tract microbiota revealed here speaks to its potential as a therapeutic target or point of intervention in genomic medicine, particularly as future studies are able to additionally account for host genetic polymorphism. Few examples yet exist where the overall composition, relative abundances, or microbial proportions of a microbiome are conclusively causal in human disease. However, it is clear that disease states are often associated with a disruption of the microbial community, frequently resulting in one or a few pathogenic organisms emerging [95, 96]. A classic example of this is the frequent ingestion of fermentable sugars that leads to increases in the mutans streptococci, etiological agents of dental caries [97]. Similarly, in the periodontal subgingival habitat, ecological shifts in redox potential facilitate the emergence of anaerobic pathogenic microbes such as the porphyromonads, which are prevalent but in low abundance in the non-diseased state [97, 98]. It is likely that microbial biomarkers at one or more body habitats will eventually be found to be prognostic indicators of future disease status, and even this reference population could contain as-yet-undetected pre-disease states. We thus hope that this profile of the human microbiota will provide a reference for subsequent investigations of its role in the onset and alleviation of diseases along the human digestive tract.

Materials and methods

Population recruitment, sample collection, and DNA purification

Healthy adults 18 to 40 years old were recruited at two academic centers [10]. Fifteen and 18 body habitats were collected from enrolled males and females, respectively. The sites sampled included anterior nares, oropharynx (two specimens), oral cavity (seven specimens), skin (four specimens), stool, and vagina (three specimens per female) [10]. The Manual of Procedures and the Core Microbiome Sampling Protocol are available at the Data Analysis and Coordination Center for the HMP [99], as well as dbGaP [100]. Genomic DNA was isolated from the collected samples using the MO Bio PowerSoil DNA Isolation Kit (MO BIO laboratories, Inc., Carlsbad, California, USA) [10].

Sequencing and binning of 16S rRNA genes and read processing

Detailed protocols used for 16S rRNA bacterial gene amplification and sequencing, using the 454 FLX Titanium platform and kits (Roche Diagnostic, Corp., Indianapolis, Indiana, USA), are available on the HMP Data Analysis and Coordination Center website [99], and are also described elsewhere [10]. In brief, sequences were processed using a data curation pipeline implemented in mothur [10, 101] starting with quality trimmed for homopolymer runs and a minimum 50 bp window average of 35. Any sequences not aligning against the appropriate subset of the SILVA database [102] were removed, as were chimeric sequences. Resulting sequences were processed using a data curation pipeline implemented in mothur [10, 101]. Remaining sequences were classified with the MSU RDP classifier v2.2 [29] using the taxonomy maintained at the RDP (RDP 10 database, version 6). Definition of a sequence's taxonomy was determined using a pseudobootstrap threshold of 80% [10].

16S rRNA gene dataset post-processing and quality control

A table of read counts from the 16S rRNA bacterial gene pipeline was created by summing clade counts from the three regions and was further processed for removing low-coverage samples. Firstly, those taxa not supported in the whole dataset by at least two sequences in at least two samples were removed. Then, the quality control procedure compared, for each sample, the read count of the most abundant taxon t and the highest abundance value that the same taxon t achieved in the entire dataset. If the former term of the comparison is <1% of the latter, the sample was discarded. Second, third, and fourth time-point samples from the same subjects were discarded. The resulting dataset of read counts containing 2,105 samples is reported in Additional file 12, which represents 210 ± 7 samples per body site. Further analysis of the dataset was performed using the per sample normalization to relative abundances. In the text, mean values are presented with standard deviation. The number of subjects with samples in the digestive tract retained for the 16S rRNA-based analysis was 209 post-quality control, from which 147 had sample data for all 10 body sites post-quality control. Unless otherwise noted, only first visit samples were used in all analyses.

Biomarker discovery and visualization

The characterization of functional and organismal features differentiating the microbial communities specific to different body sites in the digestive tract was performed using our method for biomarker discovery and explanation called LEfSe [30]. LEfSe, publicly available [103], couples a standard test for statistical significance with a quantitative test for biological consistency, finally ranking the results by effect size. Briefly, it first uses the non-parametric factorial Kruskal-Wallis test to detect features (taxonomic clades or metabolic pathways) with abundances that differ below a significance threshold among groups of samples. Biological consistency is subsequently tested using the unpaired Wilcoxon rank-sum test among all pairs of sample groups; in our case this occurred between single body habitats. Finally, linear discriminant analysis (LDA) with bootstrapping is then used to rank differentially abundant features based on their effect sizes. A significance alpha of 0.05 and an effect size threshold of 2 were used for all biomarkers discussed in this study. Organismal and functional biomarkers are graphically represented here on hierarchical trees reflecting the RDP taxonomy [29] for 16S rRNA gene data and the KEGG BRITE hierarchy [49] on KEGG modules for metagenomic functional data.

Clustering and statistical significance of four groups of body site habitats

For assessing bacterial community structure similarities between different samples and body sites, we compared the relative abundances of every pair of samples in our dataset using the Bray-Curtis measure of beta diversity [31]. The comparisons have been summarized in terms of within- and between-group averages as reported in Table 1; moreover, statistical significance has been tested for within versus between group distances, providing strong support (all P-values <10-20) for the clustering of all four groups in distinct community structures. A multidimensional scaling analysis was then performed on the Bray-Curtis diversity matrix and the four groups were denoted with different colors for highlighting the clustering structure (Additional file 4).

Whole genome shotgun sequencing, read processing, and community metabolic profiling

Whole genome shotgun sequencing employed the Illumina GAIIx platform (Illumina, Inc.) as previously described [10]. The number of samples and nucleotide content from 98 subjects is summarized in Additional file 12. The abundances and presence (or absence) of pathways in these metagenomic data were inferred using the HUMAnN pipeline (HMP Unified Metabolic Analysis Network) [48]. Briefly, the metabolic and biomolecular potential of each sample was profiled starting from the 100 bp Illumina sequences after quality and length filtering. Reads were mapped to KEGG v54 orthologous gene families (KEGG KOs [49]) using MBLASTX (MulticoreWare, St. Louis, MO, USA), an accelerated translated BLAST implementation, using default parameters and a maximum E-value of 1. Hits were mapped to abundances of each KO using up to the 20 most significant hits, weighted by the quality of each hit (inverse blastx P-value) and normalized by the length of the hit gene. Pathway information was then recovered by assigning KO gene families to KEGG modules (representing small pathways of approximately 5 to 20 genes) using a combination of MinPath [104], filtering of pathways not consistent with the BLAST-derived taxonomic composition of the community, and gap filling of likely missing enzymes. The resulting KO and KEGG module relative abundances were used in the presented analysis. Further details of the HUMAnN methodology, its software implementation, and an extensive validation of each computational step appear in [48].

Data accessibility

The datasets used in these analyses were deposited by the NIH Common Fund Human Microbiome Consortium at the Data Analysis and Coordination Center (DACC) for the Human Microbiome Project. Specifically, the downloadable packaged datasets are the 16S rRNA gene dataset [105], phylotype-classification of the 16S rRNA gene dataset [106, 107], Human Microbiome Illumina whole genome shotgun reads [108], and the metabolic reconstruction tables [109]. The phylotype classification processed for normalization and quality control is available in Additional file 7.

Abbreviations

HMP:

Human Microbiome Project

HUMAnN:

HMP Unified Metabolic Analysis Network

KEGG:

Kyoto Encyclopedia of Genes and Genomes

KO:

KEGG Orthology

LEfSe:

LDA Effect Size

RDP:

Ribosome Database Project.

References

  1. Tlaskalová-Hogenová H, Stěpánková R, Kozáková H, Hudcovic T, Vannucci L, Tučková L, Rossmann P, Hrnčíř T, Kverka M, Zákostelská Z, Klimešová K, Přibylová J, Bártová J, Sanchez D, Fundová P, Borovská D, Srůtková D, Zídek Z, Schwarzer M, Drastich P, Funda DP: The role of gut microbiota (commensal bacteria) and the mucosal barrier in the pathogenesis of inflammatory and autoimmune diseases and cancer: contribution of germ-free and gnotobiotic animal models of human diseases. Cell Mol Immunol. 2011, 8: 110-120. 10.1038/cmi.2010.67.

    PubMed  PubMed Central  Google Scholar 

  2. Tanner A, Maiden MF, Macuch PJ, Murray LL, Kent RL: Microbiota of health, gingivitis, and initial periodontitis. J Clin Periodontol. 1998, 25: 85-98. 10.1111/j.1600-051X.1998.tb02414.x.

    PubMed  CAS  Google Scholar 

  3. Kumar PS, Leys EJ, Bryk JM, Martinez FJ, Moeschberger ML, Griffen AL: Changes in periodontal health status are associated with bacterial community shifts as assessed by quantitative 16S cloning and sequencing. J Clin Microbiol. 2006, 44: 3665-3673. 10.1128/JCM.00317-06.

    PubMed  CAS  PubMed Central  Google Scholar 

  4. Al-Attas OS, Al-Daghri NM, Al-Rubeaan K, da Silva NF, Sabico SL, Kumar S, McTernan PG, Harte AL: Changes in endotoxin levels in T2DM subjects on anti-diabetic therapies. Cardiovasc Diabetol. 2009, 8: 20-10.1186/1475-2840-8-20.

    PubMed  PubMed Central  Google Scholar 

  5. Creely SJ, McTernan PG, Kusminski CM, Fisher M, Da Silva NF, Khanolkar M, Evans M, Harte AL, Kumar S: Lipopolysaccharide activates an innate immune system response in human adipose tissue in obesity and type 2 diabetes. Am J Physiol Endocrinol Metab. 2007, 292: E740-747.

    PubMed  CAS  Google Scholar 

  6. Ott SJ, El Mokhtari NE, Musfeldt M, Hellmig S, Freitag S, Rehman A, Kühbacher T, Nikolaus S, Namsolleck P, Blaut M, Hampe J, Sahly H, Reinecke A, Haake N, Günther R, Krüger D, Lins M, Herrmann G, Fölsch UR, Simon R, Schreiber S: Detection of diverse bacterial signatures in atherosclerotic lesions of patients with coronary heart disease. Circulation. 2006, 113: 929-937. 10.1161/CIRCULATIONAHA.105.579979.

    PubMed  Google Scholar 

  7. Ley RE: Obesity and the human microbiome. Curr Opin Gastroenterol. 2010, 26: 5-11. 10.1097/MOG.0b013e328333d751.

    PubMed  Google Scholar 

  8. Pussinen PJ, Tuomisto K, Jousilahti P, Havulinna AS, Sundvall J, Salomaa V: Endotoxemia, immune response to periodontal pathogens, and systemic inflammation associate with incident cardiovascular disease events. Arterioscler Thromb Vasc Biol. 2007, 27: 1433-1439. 10.1161/ATVBAHA.106.138743.

    PubMed  CAS  Google Scholar 

  9. NIH HMP Working Group, Peterson J, Garges S, Giovanni M, McInnes P, Wang L, Schloss JA, Bonazzi V, McEwen JE, Wetterstrand KA, Deal C, Baker CC, Di Francesco V, Howcroft TK, Karp RW, Lunsford RD, Wellington CR, Belachew T, Wright M, Giblin C, David H, Mills M, Salomon R, Mullins C, Akolkar B, Begg L, Davis C, Grandison L, Humble M, Khalsa J, et al: The NIH Human Microbiome Project. Genome Res. 2009, 19: 2317-2323.

    Google Scholar 

  10. The Human Microbiome Project Consortium: A framework for human microbiome research. Nature. 2012, doi:10.1038/nature11209

    Google Scholar 

  11. Xie G, Chain PS, Lo CC, Liu KL, Gans J, Merritt J, Qi F: Community and gene composition of a human dental plaque microbiota obtained by metagenomic sequencing. Mol Oral Microbiol. 2010, 25: 391-405. 10.1111/j.2041-1014.2010.00587.x.

    PubMed  CAS  PubMed Central  Google Scholar 

  12. Andersson AF, Lindberg M, Jakobsson H, Backhed F, Nyren P, Engstrand L: Comparative analysis of human gut microbiota by barcoded pyrosequencing. PLoS One. 2008, 3: e2836-10.1371/journal.pone.0002836.

    PubMed  PubMed Central  Google Scholar 

  13. Aas JA, Paster BJ, Stokes LN, Olsen I, Dewhirst FE: Defining the normal bacterial flora of the oral cavity. J Clin Microbiol. 2005, 43: 5721-5732. 10.1128/JCM.43.11.5721-5732.2005.

    PubMed  PubMed Central  Google Scholar 

  14. Bik EM, Long CD, Armitage GC, Loomer P, Emerson J, Mongodin EF, Nelson KE, Gill SR, Fraser-Liggett CM, Relman DA: Bacterial diversity in the oral cavity of 10 healthy individuals. ISME J. 2010, 4: 962-974. 10.1038/ismej.2010.30.

    PubMed  PubMed Central  Google Scholar 

  15. 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.

    PubMed  PubMed Central  Google Scholar 

  16. Bik EM, Eckburg PB, Gill SR, Nelson KE, Purdom EA, Francois F, Perez-Perez G, Blaser MJ, Relman DA: Molecular analysis of the bacterial microbiota in the human stomach. Proc Natl Acad Sci USA. 2006, 103: 732-737. 10.1073/pnas.0506655103.

    PubMed  CAS  PubMed Central  Google Scholar 

  17. Pei Z, Bini EJ, Yang L, Zhou M, Francois F, Blaser MJ: Bacterial biota in the human distal esophagus. Proc Natl Acad Sci USA. 2004, 101: 4250-4255. 10.1073/pnas.0306398101.

    PubMed  CAS  PubMed Central  Google Scholar 

  18. Lemon KP, Klepac-Ceraj V, Schiffer HK, Brodie EL, Lynch SV, Kolter R: Comparative analyses of the bacterial microbiota of the human nostril and oropharynx. MBio. 2010, 1: e00129-00110.

    PubMed  PubMed Central  Google Scholar 

  19. 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.

    PubMed  CAS  PubMed Central  Google Scholar 

  20. Tap J, Mondot S, Levenez F, Pelletier E, Caron C, Furet JP, Ugarte E, Muñoz-Tamayo R, Paslier DL, Nalin R, Dore J, Leclerc M: Towards the human intestinal microbiota phylogenetic core. Environ Microbiol. 2009, 11: 2574-2584. 10.1111/j.1462-2920.2009.01982.x.

    PubMed  Google Scholar 

  21. Charlson ES, Chen J, Custers-Allen R, Bittinger K, Li H, Sinha R, Hwang J, Bushman FD, Collman RG: Disordered microbial communities in the upper respiratory tract of cigarette smokers. PLoS One. 2010, 5: e15216-10.1371/journal.pone.0015216.

    PubMed  CAS  PubMed Central  Google Scholar 

  22. Nasidze I, Li J, Quinque D, Tang K, Stoneking M: Global diversity in the human salivary microbiome. Genome Res. 2009, 19: 636-643. 10.1101/gr.084616.108.

    PubMed  CAS  PubMed Central  Google Scholar 

  23. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, Mende DR, Li J, Xu J, Li S, Li D, Cao J, Wang B, Liang H, Zheng H, Xie Y, Tap J, Lepage P, Bertalan M, Batto JM, Hansen T, Le Paslier D, Linneberg A, Nielsen HB, Pelletier E, Renault P, et al: A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010, 464: 59-65. 10.1038/nature08821.

    PubMed  CAS  PubMed Central  Google Scholar 

  24. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, Gordon JI: A core gut microbiome in obese and lean twins. Nature. 2009, 457: 480-484. 10.1038/nature07540.

    PubMed  CAS  PubMed Central  Google Scholar 

  25. Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto JM, Bertalan M, Borruel N, Casellas F, Fernandez L, Gautier L, Hansen T, Hattori M, Hayashi T, Kleerebezem M, Kurokawa K, Leclerc M, Levenez F, Manichanh C, Nielsen HB, Nielsen T, Pons N, Poulain J, Qin J, Sicheritz-Ponten T, Tims S, et al: Enterotypes of the human gut microbiome. Nature. 2011, 473: 174-180. 10.1038/nature09944.

    PubMed  CAS  PubMed Central  Google Scholar 

  26. Claesson MJ, Cusack S, O'Sullivan O, Greene-Diniz R, de Weerd H, Flannery E, Marchesi JR, Falush D, Dinan T, Fitzgerald G, Stanton C, van Sinderen D, O'Connor M, Harnedy N, O'Connor K, Henry C, O'Mahony D, Fitzgerald AP, Shanahan F, Twomey C, Hill C, Ross RP, O'Toole PW: Composition, variability, and temporal stability of the intestinal microbiota of the elderly. Proc Natl Acad Sci USA. 2011, 108 (Suppl 1): 4586-4591.

    PubMed  CAS  PubMed Central  Google Scholar 

  27. 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: 1-10.

    Google Scholar 

  28. Savage DC: Microbial ecology of the gastrointestinal tract. Annu Rev Microbiol. 1977, 31: 107-133. 10.1146/annurev.mi.31.100177.000543.

    PubMed  CAS  Google Scholar 

  29. 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-D145. 10.1093/nar/gkn879.

    PubMed  CAS  PubMed Central  Google Scholar 

  30. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C: Metagenomic Biomarker Discovery and Explanation. Genome Biol. 2011, 12: R60-10.1186/gb-2011-12-6-r60.

    PubMed  PubMed Central  Google Scholar 

  31. Bray JR, Curtis JT: An ordination of the upland forest communities of southern Wisconsin. Ecol Monographs. 1957, 27: 325-349. 10.2307/1942268.

    Google Scholar 

  32. Simpson EH: Measurement of diversity. Nature. 1949, 163: 1-

    Google Scholar 

  33. Brinig MM, Lepp PW, Ouverney CC, Armitage GC, Relman DA: Prevalence of bacteria of division TM7 in human subgingival plaque and their association with disease. Appl Environ Microbiol. 2003, 69: 1687-1694. 10.1128/AEM.69.3.1687-1694.2003.

    PubMed  CAS  PubMed Central  Google Scholar 

  34. Zoetendal EG, Plugge CM, Akkermans AD, de Vos WM: Victivallis vadensis gen. nov, sp nov, a sugar-fermenting anaerobe from human faeces Int J Syst Evol Microbiol. 2003, 53: 211-215.

    PubMed  Google Scholar 

  35. Derrien M, Vaughan EE, Plugge CM, de Vos WM: Akkermansia muciniphila gen. nov, sp nov, a human intestinal mucin-degrading bacterium Int J Syst Evol Microbiol. 2004, 54: 1469-1476.

    PubMed  CAS  Google Scholar 

  36. Downes J, Vartoukian SR, Dewhirst FE, Izard J, Chen T, Yu W, Sutliffe IC, Wade WG: Pyramidobacter pisciolens gen. nov, sp nov, a member of the phylum 'Synergistetes' isolated from the human oral cavity. Int J Syst Evol Microbiol. 2009, 59: 972-980.

    CAS  Google Scholar 

  37. Armitage GC, Dickinson WR, Jenderseck RS, Levine SM, Chambers DW: Relationship between the percentage of subgingival spirochetes and the severity of periodontal disease. J Periodontol. 1982, 53: 550-556. 10.1902/jop.1982.53.9.550.

    PubMed  CAS  Google Scholar 

  38. Cavrini F, Pirani C, Foschi F, Montebugnoli L, Sambri V, Prati C: Detection of Treponema denticola in root canal systems in primary and secondary endodontic infections. A correlation with clinical symptoms. New Microbiol. 2008, 31: 67-73.

    PubMed  Google Scholar 

  39. Kanasi E, Dewhirst FE, Chalmers NI, Kent R, Moore A, Hughes CV, Pradhan N, Loo CY, Tanner AC: Clonal analysis of the microbiota of severe early childhood caries. Caries Res. 2010, 44: 485-497. 10.1159/000320158.

    PubMed  CAS  PubMed Central  Google Scholar 

  40. Mantzourani M, Gilbert SC, Sulong HN, Sheehy EC, Tank S, Fenlon M, Beighton D: The isolation of bifidobacteria from occlusal carious lesions in children and adults. Caries Res. 2009, 43: 308-313. 10.1159/000222659.

    PubMed  CAS  Google Scholar 

  41. Schloss PD: The effects of alignment quality, distance calculation method, sequence filtering, and region on the analysis of 16S rRNA gene-based studies. PLoS Comput Biol. 2010, 6: e1000844-10.1371/journal.pcbi.1000844.

    PubMed  PubMed Central  Google Scholar 

  42. Ford AC, Axon ATR: Epidemiology of Helicobacter pylori infection and public health implications. Helicobacter. 2010, 15: 1-6.

    PubMed  Google Scholar 

  43. Hold GL, Schwiertz A, Aminov RI, Blaut M, Flint HJ: Oligonucleotide probes that detect quantitatively significant groups of butyrate-producing bacteria in human feces. Appl Environ Microbiol. 2003, 69: 4320-4324. 10.1128/AEM.69.7.4320-4324.2003.

    PubMed  CAS  PubMed Central  Google Scholar 

  44. 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.

    PubMed  Google Scholar 

  45. Martinez-Medina M, Aldeguer X, Gonzalez-Huix F, Acero D, Garcia-Gil LJ: Abnormal microbiota composition in the ileocolonic mucosa of Crohn's disease patients as revealed by polymerase chain reaction-denaturing gradient gel electrophoresis. Inflamm Bowel Dis. 2006, 12: 1136-1145. 10.1097/01.mib.0000235828.09305.0c.

    PubMed  Google Scholar 

  46. Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermúdez-Humarán LG, Gratadoux JJ, Blugeon S, Bridonneau C, Furet JP, Corthier G, Grangette C, Vasquez N, Pochart P, Trugnan G, Thomas G, Blottière HM, Doré J, Marteau P, Seksik P, Langella P: Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci USA. 2008, 105: 16731-16736. 10.1073/pnas.0804812105.

    PubMed  CAS  PubMed Central  Google Scholar 

  47. Kenney EB, Ash MM: Oxidation reduction potential of developing plaque, periodontal pockets and gingival sulci. J Periodontol. 1969, 40: 630-633.

    PubMed  CAS  Google Scholar 

  48. Abubucker S, Segata N, Goll J, Schubert AM, Izard J, Cantarel BL, Rodriguez-Mueller B, Zucker J, Henrissat B, White O, Kelley ST, Methé B, Schloss PD, Gevers D, Mitreva M, Huttenhower C: Metabolic reconstruction for metagenomic data and its application to the human microbiome. PLoS Comput Biol. 2012, doi 10.1371/journal.pcbi.1002358

    Google Scholar 

  49. 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-D360. 10.1093/nar/gkp896.

    PubMed  CAS  PubMed Central  Google Scholar 

  50. Barnes VM, Teles R, Trivedi HM, Devizio W, Xu T, Mitchell MW, Milburn MV, Guo L: Acceleration of purine degradation by periodontal diseases. J Dent Res. 2009, 88: 851-855. 10.1177/0022034509341967.

    PubMed  CAS  Google Scholar 

  51. Miller TL, Wolin MJ: Pathways of acetate, propionate, and butyrate formation by the human fecal microbial flora. Appl Environ Microbiol. 1996, 62: 1589-1592.

    PubMed  CAS  PubMed Central  Google Scholar 

  52. Flint HJ, Bayer EA: Plant cell wall breakdown by anaerobic microorganisms from the Mammalian digestive tract. Ann N Y Acad Sci. 2008, 1125: 280-288. 10.1196/annals.1419.022.

    PubMed  CAS  Google Scholar 

  53. Prakash R, Mullen KD: Mechanisms, diagnosis and management of hepatic encephalopathy. Nat Rev Gastroenterol Hepatol. 2010, 7: 515-525. 10.1038/nrgastro.2010.116.

    PubMed  Google Scholar 

  54. The Human Microbiome Project Consortium: Structure, function and diversity of the healthy human microbiome. Nature. 2012, doi:10.1038/nature11234

    Google Scholar 

  55. Soukos NS, Som S, Abernethy AD, Ruggiero K, Dunham J, Lee C, Doukas AG, Goodson JM: Phototargeting oral black-pigmented bacteria. Antimicrob Agents Chemother. 2005, 49: 1391-1396. 10.1128/AAC.49.4.1391-1396.2005.

    PubMed  CAS  PubMed Central  Google Scholar 

  56. Anzaldi LL, Skaar EP: Overcoming the heme paradox: heme toxicity and tolerance in bacterial pathogens. Infect Immun. 2010, 78: 4977-4989. 10.1128/IAI.00613-10.

    PubMed  CAS  PubMed Central  Google Scholar 

  57. Karlsen OA, Ramsevik L, Bruseth LJ, Larsen O, Brenner A, Berven FS, Jensen HB, Lillehaug JR: Characterization of a prokaryotic haemerythrin from the methanotrophic bacterium Methylococcus capsulatus (Bath). FEBS J. 2005, 272: 2428-2440. 10.1111/j.1742-4658.2005.04663.x.

    PubMed  CAS  Google Scholar 

  58. Kajiya M, Silva MJ, Sato K, Ouhara K, Kawai T: Hydrogen mediates suppression of colon inflammation induced by dextran sodium sulfate. Biochem Biophys Res Commun. 2009, 386: 11-15. 10.1016/j.bbrc.2009.05.117.

    PubMed  CAS  Google Scholar 

  59. Van Palenstein Helderman WH, Rosman I: Hydrogen-dependent organisms from the human gingival crevice resembling Vibrio succinogenes. Antonie Van Leeuwenhoek. 1976, 42: 107-118. 10.1007/BF00399454.

    PubMed  CAS  Google Scholar 

  60. Roediger WE, Duncan A, Kapaniris O, Millard S: Reducing sulfur compounds of the colon impair colonocyte nutrition: implications for ulcerative colitis. Gastroenterology. 1993, 104: 802-809.

    PubMed  CAS  Google Scholar 

  61. Roediger WE: The colonic epithelium in ulcerative colitis: an energy-deficiency disease?. Lancet. 1980, 2: 712-715.

    PubMed  CAS  Google Scholar 

  62. Wallace JL, Dicay M, McKnight W, Martin GR: Hydrogen sulfide enhances ulcer healing in rats. FASEB J. 2007, 21: 4070-4076. 10.1096/fj.07-8669com.

    PubMed  CAS  Google Scholar 

  63. Attene-Ramos MS, Wagner ED, Plewa MJ, Gaskins HR: Evidence that hydrogen sulfide is a genotoxic agent. Mol Cancer Res. 2006, 4: 9-14. 10.1158/1541-7786.MCR-05-0126.

    PubMed  CAS  Google Scholar 

  64. Nicholls P, Kim JK: Sulphide as an inhibitor and electron donor for the cytochrome c oxidase system. Can J Biochem. 1982, 60: 613-623. 10.1139/o82-076.

    PubMed  CAS  Google Scholar 

  65. Berglin EH, Carlsson J: Potentiation by sulfide of hydrogen peroxide-induced killing of Escherichia coli. Infect Immun. 1985, 49: 538-543.

    PubMed  CAS  PubMed Central  Google Scholar 

  66. Keijser BJ, Zaura E, Huse SM, van der Vossen JM, Schuren FH, Montijn RC, ten Cate JM, Crielaard W: Pyrosequencing analysis of the oral microflora of healthy adults. J Dent Res. 2008, 87: 1016-1020. 10.1177/154405910808701104.

    PubMed  CAS  Google Scholar 

  67. Lazarevic V, Whiteson K, Huse S, Hernandez D, Farinelli L, Osteras M, Schrenzel J, Francois P: Metagenomic study of the oral microbiota by Illumina high-throughput sequencing. J Microbiol Methods. 2009, 79: 266-271. 10.1016/j.mimet.2009.09.012.

    PubMed  CAS  PubMed Central  Google Scholar 

  68. Gandon S, Mackinnon MJ, Nee S, Read AF: Imperfect vaccines and the evolution of pathogen virulence. Nature. 2001, 414: 751-756. 10.1038/414751a.

    PubMed  CAS  Google Scholar 

  69. Lenski RE, May RM: The evolution of virulence in parasites and pathogens: reconciliation between two competing hypotheses. J Theoret Biol. 1994, 169: 253-265. 10.1006/jtbi.1994.1146.

    CAS  Google Scholar 

  70. Little TJ, Shuker DM, Colegrave N, Day T, Graham AL: The coevolution of virulence: tolerance in perspective. PLoS pathogens. 2010, 6: e1001006-10.1371/journal.ppat.1001006.

    PubMed  PubMed Central  Google Scholar 

  71. Ley RE, Turnbaugh PJ, Klein S, Gordon JI: Microbial ecology: human gut microbes associated with obesity. Nature. 2006, 444: 1022-1023. 10.1038/4441022a.

    PubMed  CAS  Google Scholar 

  72. Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, Schlegel ML, Tucker TA, Schrenzel MD, Knight R, Gordon JI: Evolution of mammals and their gut microbes. Science. 2008, 320: 1647-1651. 10.1126/science.1155725.

    PubMed  CAS  PubMed Central  Google Scholar 

  73. 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.

    PubMed  CAS  PubMed Central  Google Scholar 

  74. De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, Collini S, Pieraccini G, Lionetti P: Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci USA. 2010, 107: 14691-14696. 10.1073/pnas.1005963107.

    PubMed  PubMed Central  Google Scholar 

  75. Kurokawa K, Itoh T, Kuwahara T, Oshima K, Toh H, Toyoda A, Takami H, Morita H, Sharma VK, Srivastava TP, Taylor TD, Noguchi H, Mori H, Ogura Y, Ehrlich DS, Itoh K, Takagi T, Sakaki Y, Hayashi T, Hattori M: Comparative metagenomics revealed commonly enriched gene sets in human gut microbiomes. DNA Res. 2007, 14: 169-181. 10.1093/dnares/dsm018.

    PubMed  CAS  PubMed Central  Google Scholar 

  76. Koenig JE, Spor A, Scalfone N, Fricker AD, Stombaugh J, Knight R, Angenent LT, Ley RE: Succession of microbial consortia in the developing infant gut microbiome. Proc Natl Acad Sci USA. 2011, 108 (Suppl 1): 4578-4585.

    PubMed  CAS  PubMed Central  Google Scholar 

  77. Wang X, Heazlewood SP, Krause DO, Florin TH: Molecular characterization of the microbial species that colonize human ileal and colonic mucosa by using 16S rDNA sequence analysis. J Appl Microbiol. 2003, 95: 508-520. 10.1046/j.1365-2672.2003.02005.x.

    PubMed  CAS  Google Scholar 

  78. Socransky SS, Haffajee AD: Dental biofilms: difficult therapeutic targets. Periodontology 2000. 2002, 28: 12-55. 10.1034/j.1600-0757.2002.280102.x.

    PubMed  Google Scholar 

  79. Richardson RL, Jones M: A bacteriologic census of human saliva. J Dent Res. 1958, 37: 697-709. 10.1177/00220345580370041701.

    PubMed  CAS  Google Scholar 

  80. Seville LA, Patterson AJ, Scott KP, Mullany P, Quail MA, Parkhill J, Ready D, Wilson M, Spratt D, Roberts AP: Distribution of tetracycline and erythromycin resistance genes among human oral and fecal metagenomic DNA. Microb Drug Resist. 2009, 15: 159-166. 10.1089/mdr.2009.0916.

    PubMed  CAS  Google Scholar 

  81. Dowd SE, Callaway TR, Wolcott RD, Sun Y, McKeehan T, Hagevoort RG, Edrington TS: Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP). BMC Microbiol. 2008, 8: 125-10.1186/1471-2180-8-125.

    PubMed  PubMed Central  Google Scholar 

  82. Zhang H, DiBaise JK, Zuccolo A, Kudrna D, Braidotti M, Yu Y, Parameswaran P, Crowell MD, Wing R, Rittmann BE, Krajmalnik-Brown R: Human gut microbiota in obesity and after gastric bypass. Proc Natl Acad Sci USA. 2009, 106: 2365-2370. 10.1073/pnas.0812600106.

    PubMed  CAS  PubMed Central  Google Scholar 

  83. Dal Bello F, Hertel C: Oral cavity as natural reservoir for intestinal lactobacilli. Syst Appl Microbiol. 2006, 29: 69-76. 10.1016/j.syapm.2005.07.002.

    PubMed  CAS  Google Scholar 

  84. Maukonen J, Matto J, Suihko ML, Saarela M: Intra-individual diversity and similarity of salivary and faecal microbiota. J Med Microbiol. 2008, 57: 1560-1568. 10.1099/jmm.0.47352-0.

    PubMed  CAS  Google Scholar 

  85. Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner AC, Yu WH, Lakshmanan A, Wade WG: The human oral microbiome. J Bacteriol. 2010, 192: 5002-5017. 10.1128/JB.00542-10.

    PubMed  CAS  PubMed Central  Google Scholar 

  86. Wilson M: Bacteriology of Humans an Ecological Perspective. 2008, Blackwell Publishing Ltd

    Google Scholar 

  87. Dictionary of Cancer Terms. [http://www.cancer.gov/dictionary]

  88. Goss CH, Burns JL: Exacerbations in cystic fibrosis. 1: Epidemiology and pathogenesis. Thorax. 2007, 62: 360-367. 10.1136/thx.2006.060889.

    PubMed  PubMed Central  Google Scholar 

  89. Sibley CD, Parkins MD, Rabin HR, Duan K, Norgaard JC, Surette MG: A polymicrobial perspective of pulmonary infections exposes an enigmatic pathogen in cystic fibrosis patients. Proc Natl Acad Sci USA. 2008, 105: 15070-15075. 10.1073/pnas.0804326105.

    PubMed  CAS  PubMed Central  Google Scholar 

  90. Duan K, Dammel C, Stein J, Rabin H, Surette MG: Modulation of Pseudomonas aeruginosa gene expression by host microflora through interspecies communication. Mol Microbiol. 2003, 50: 1477-1491. 10.1046/j.1365-2958.2003.03803.x.

    PubMed  CAS  Google Scholar 

  91. Kumar PS, Griffen AL, Barton JA, Paster BJ, Moeschberger ML, Leys EJ: New bacterial species associated with chronic periodontitis. J Dent Res. 2003, 82: 338-344. 10.1177/154405910308200503.

    PubMed  CAS  Google Scholar 

  92. Zijnge V, van Leeuwen MB, Degener JE, Abbas F, Thurnheer T, Gmur R, Harmsen HJ: Oral biofilm architecture on natural teeth. PLoS One. 2010, 5: e9321-10.1371/journal.pone.0009321.

    PubMed  PubMed Central  Google Scholar 

  93. Kumar AG, Nagesh N, Prabhakar TG, Sekaran G: Purification of extracellular acid protease and analysis of fermentation metabolites by Synergistes sp. utilizing proteinaceous solid waste from tanneries. Bioresour Technol. 2008, 99: 2364-2372. 10.1016/j.biortech.2007.05.001.

    PubMed  CAS  Google Scholar 

  94. Godon JJ, Moriniere J, Moletta M, Gaillac M, Bru V, Delgenes JP: Rarity associated with specific ecological niches in the bacterial world: the 'Synergistes' example. Environ Microbiol. 2005, 7: 213-224. 10.1111/j.1462-2920.2004.00693.x.

    PubMed  CAS  Google Scholar 

  95. Nataro JP, Kaper JB: Diarrheagenic Escherichia coli. Clin Microbiol Rev. 1998, 11: 142-201.

    PubMed  CAS  PubMed Central  Google Scholar 

  96. Sakamoto M, Huang Y, Ohnishi M, Umeda M, Ishikawa I, Benno Y: Changes in oral microbial profiles after periodontal treatment as determined by molecular analysis of 16S rRNA genes. J Med Microbiol. 2004, 53: 563-571. 10.1099/jmm.0.45576-0.

    PubMed  CAS  Google Scholar 

  97. Marsh PD, Bradshaw DJ: Physiological approaches to the control of oral biofilms. Adv Dental Res. 1997, 11: 176-185. 10.1177/08959374970110010901.

    CAS  Google Scholar 

  98. Socransky SS, Haffajee AD, Cugini MA, Smith C, Kent RL: Microbial complexes in subgingival plaque. J Clin Periodontol. 1998, 25: 134-144. 10.1111/j.1600-051X.1998.tb02419.x.

    PubMed  CAS  Google Scholar 

  99. Human Microbiome Project: Tools & Protocols. [http://www.hmpdacc.org/tools_protocols/tools_protocols.php]

  100. NIH Human Microbiome Project - Core Microbiome Sampling Protocol A (HMP-A). [http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000228.v3.p1]

  101. 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: 7537-7541. 10.1128/AEM.01541-09.

    PubMed  CAS  PubMed Central  Google Scholar 

  102. 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: 7188-7196. 10.1093/nar/gkm864.

    PubMed  CAS  PubMed Central  Google Scholar 

  103. Galaxy/Huttenhower Lab. [http://huttenhower.sph.harvard.edu/lefse/]

  104. 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.

    PubMed  PubMed Central  Google Scholar 

  105. HMR16S - Raw 16S Data and Library Metadata. [http://www.hmpdacc.org/HMR16S]

  106. Schloss PD, Gevers D, Westcott SL: Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS ONE. 2011, 6: e27310-10.1371/journal.pone.0027310.

    PubMed  CAS  PubMed Central  Google Scholar 

  107. HMMCP - mothur Community Profiling. [http://www.hmpdacc.org/HMMCP]

  108. HMIWGS/HMASM - Illumina WGS Reads and Assemblies. [http://www.hmpdacc.org/HMASM]

  109. HMMRC - Metabolic Reconstruction. [http://www.hmpdacc.org/HMMRC]

Download references

Acknowledgements

We thank the other members of the Human Microbiome Project consortium for study design and data production with special thanks to Patrick Schloss and the HMP Metabolic Reconstruction Group for providing the tables from which this analysis is derived. The human subjects who participated in this study are gratefully acknowledged. We thank Joshua A Reyes for participating in the production of Additional file 11. This work was supported by the Army Research Office (ARO) under award W911NF-11-1-0473 to CH, by the National Science Foundation (NSF) under award DBI-1053486 to CH, and by the National Institutes of Health (NIH) under awards HG005969 to CH, HG004969 to DG, CA139193 to JI, DE020751 to KPL, DE020298 and DE021574 to SKH.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jacques Izard.

Additional information

Competing interests

The authors declare that they have no competing interests. While the National Institutes of Health were one of the major drivers for the creation of the Human Microbiome Project, the NIH had no role in data analysis, decision to publish, or preparation of the manuscript.

Authors' contributions

CH, DG, JI, KPL, LW, NS, PM, and SKH analyzed the data. CH, DG, NS, and JI contributed analysis tools. CH, JI, KPL, PM and SKH wrote the paper. All authors have read and approved the manuscript for publication.

Susan Kinder Haake, Peter Mannon, Katherine P Lemon contributed equally to this work.

Electronic supplementary material

13059_2012_2863_MOESM1_ESM.XLSX

Additional file 1: Table s1 - average abundance, expressed in percentage of all microbial clades in the four digestive tract groups and among the ten body habitats. Lettering of groups and body habitats are as in Figure 1. AVG, average; STDEV, standard deviation. (XLSX 245 KB)

13059_2012_2863_MOESM2_ESM.XLSX

Additional file 2: Table s2 - surfaces associated with the sampling sites from which the microbiota of the digestive tract was collected. (XLSX 12 KB)

13059_2012_2863_MOESM3_ESM.PDF

Additional file 3: Figure s1 - higher resolution version of Figure 1bshowing significantly enriched taxa from the four groups of digestive tract sites. This circular cladogram reports significant group-enriched taxa. Differential taxa analysis was performed using LEfSe on all the samples. Colored shading highlights which of the four major bacterial phyla was most enriched in which of the four body site groups. Each colored dot indicates a taxon that was differentially abundant among the groups. Small letters denote bacterial families that were enriched in one of the four body site groups. (PDF 2 MB)

13059_2012_2863_MOESM4_ESM.PDF

Additional file 4: Figure s2 - diversity-based multidimensional scaling (MDS) plot of samples. A distance matrix for all pairwise distances between samples was calculated using Bray-Curtis distance, which was used to project samples to MDS coordinates using the stats::cmdscale R function with default options. Each of the four established groups of body sites (G1, G2, G3, G4) is assigned a color, decreasing in opacity as the density of points of that group decreases, and body sites are denoted with different marker types. G2 and G3 contain the most overlap, while maintaining distinct areas of highest density, while G1 and G4, respectively, increase in distinctness. The distribution of samples in specific body sites does not produce sub-clusters, confirming the homogeneity of bacterial community composition within the four groups. (PDF 92 KB)

13059_2012_2863_MOESM5_ESM.XLSX

Additional file 5: Table s3 - inverse Simpson for each habitat of the digestive tract. The minimum, maximum, average and standard deviation values are reported. (XLSX 12 KB)

13059_2012_2863_MOESM6_ESM.XLSX

Additional file 6: Table s4 - percentages of subjects for whom each taxon was detected in both the upper digestive tract and in the stool. The table is ordered based on the absolute differences between the presence in the stool and in at least one oral body site. Only the subjects with samples in all ten digestive tract body habitats were considered (n = 147) and all the taxonomic units with at least 40% of presence in stool or any oral body site are included. (XLSX 15 KB)

13059_2012_2863_MOESM7_ESM.XLSX

Additional file 7: Table s5 - read counts for all digestive tract samples (after quality control) for each microbial clade. (XLSX 7 MB)

13059_2012_2863_MOESM8_ESM.PDF

Additional file 8: Figure s3 - visual and schematic representation of the oral cavity and oropharyngeal sampling sites. The soft tissues, illustrated here in a 20-year-old healthy male, were sampled by swabbing the tongue dorsum, hard palate, right and left buccal mucosa, the anterior keratinized gingiva, the right and left palatine tonsils, and the throat (posterior wall of the oropharynx). The pooled supragingival and pooled subgingival plaque samples were taken with curettes from molars, premolars and incisors (schematic illlustration). Not shown is the sampling of the saliva, which was collected by having the subject drool accumulated saliva into a collection vial. The complete sampling procedure is described in the Manual of Procedures for Human Microbiome Project (see Materials and methods). (PDF 545 KB)

13059_2012_2863_MOESM9_ESM.PDF

Additional file 9: Figure s4 - higher resolution version of Figure 6showing functional characterization of the digestive microbiota. Differentially abundant metabolic pathways from the buccal mucosa, supragingival plaque, tongue dorsum, and stool are depicted based on metabolic profiling performed with HUMAnN [48] from metagenomic shotgun sequencing data. Lettering indicates metabolic modules. Nucleot./amino acid met., nucleotide and amino acid metabolism; Carbohydrate/lipid met., carbohydrate and lipid metabolism; Energy met., energy metabolism; Met., aminoacyl tRNA and nucleotide sugar metabolism; Genetic information proc., genetic information processes; Environmental inf. proc., environmental information processing. (PDF 2 MB)

13059_2012_2863_MOESM10_ESM.XLSX

Additional file 10: Table s6 - percentages of metagenomic reads assigned to Archaea, Bacteria, and non-human Eukaryota (human reads excluded) in the four digestive tract sites with more than 50 shotgun sequencing samples available. (XLSX 13 KB)

13059_2012_2863_MOESM11_ESM.PDF

Additional file 11: Figure s5 - a subset of significant correlations between metagenomic gene family and organismal abundances. Paired shotgun metagenomic and 16S rRNA gene sequencing samples were associated, resulting in 34 buccal mucosa, 35 stool, 33 supragingival plaque, and 30 tongue microbiomes for joint analysis. Within each body site, Spearman correlations were calculated between the 21 KEGG Orthology gene families described in the Results and all phylotypes at any taxonomic level from phylum to OTU. Significant associations reaching a Benjamini-Hochberg false discovery rate <0.05 are shown here; grey ellipses represent clades, white rectangles KO gene families, and edge width is proportional to -log(q-value). Colors are as in Figure 1 (red, buccal mucosa; green, stool; yellow, plaque; blue, tongue). (PDF 292 KB)

13059_2012_2863_MOESM12_ESM.XLSX

Additional file 12: Table s7 - summary of the read statistics for 16S rRNA gene taxonomic abundances and whole genome shotgun sequencing metagenomic data. (XLSX 11 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Segata, N., Haake, S.K., Mannon, P. et al. Composition of the adult digestive tract bacterial microbiome based on seven mouth surfaces, tonsils, throat and stool samples. Genome Biol 13, R42 (2012). https://doi.org/10.1186/gb-2012-13-6-r42

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/gb-2012-13-6-r42

Keywords