OUP user menu

Distal gut microbiota of adolescent children is different from that of adults

Richard Agans , Laura Rigsbee , Harshavardhan Kenche , Sonia Michail , Harry J. Khamis , Oleg Paliy
DOI: http://dx.doi.org/10.1111/j.1574-6941.2011.01120.x 404-412 First published online: 1 August 2011

Abstract

Human intestinal microbiota plays a number of important roles in human health and is also implicated in several gastrointestinal disorders. Although the diversity of human gut microbiota in adults and in young children has been examined, few reports of microbiota composition are available for adolescents. In this work, we used Microbiota Array for high-throughput analysis of distal gut microbiota in adolescent children 11–18 years of age. Samples obtained from healthy adults were used for comparison. Adolescent and adult groups could be separated in the principal components analysis space based on the relative species abundance of their distal gut microbiota. All samples were dominated by class Clostridia. A core microbiome of 46 species that were detected in all examined samples was established; members of genera Ruminococcus, Faecalibacterium, and Roseburia were well represented among core species. Comparison of intestinal microbiota composition between adolescents and adults revealed a statistically significantly higher abundance of genera Bifidobacterium and Clostridium among adolescent samples. The number of detected species was similar between sample groups, indicating that it was the relative abundances of the genera and not the presence or absence of a specific genus that differentiated adolescent and adult samples. In summary, contrary to the current belief, this study suggests that the gut microbiome of adolescent children is different from that of adults.

Keywords
  • adolescent
  • intestinal microbiota
  • microflora
  • microbial communities
  • phylogenetic microarray

Introduction

The human gastrointestinal tract is an ecosystem rich in microbial species, with an estimated 1013–1014 cells present in the large intestine (Ley et al., 2006). The sheer numbers of microorganisms and the enclosed space of the intestinal tract create close relationships between human host and microorganisms. This intestinal microbial community is an essential component of human health and disease. The intestinal microbiota is responsible for polysaccharide hydrolysis and fermentation, vitamin production, immune system stimulation, modulation of gut motility, as well as for protection of human host from pathogen invasion (Sekirov et al., 2010). At the same time, changes in intestinal microbiota are associated with a variety of gastrointestinal disorders such as irritable bowel syndrome, inflammatory bowel disease, obesity, and gastrointestinal cancer (Turnbaugh et al., 2006; Frank et al., 2007; Swidsinski et al., 2008; Davis & Milner, 2009; Fujimura et al., 2010).

Although the original studies of intestinal microbial communities were limited by the difficulties of utilizing traditional microbiological techniques for the analysis of this complex community, new methods based on the interrogation of total community DNA and RNA provide means to analyze and compare the presence and abundance of hundreds of microbial species in a single sample (Langendijk et al., 1995; Suau et al., 1999; Wang et al., 2004; Ley et al., 2005; Palmer et al., 2006; Biasucci et al., 2010). A number of projects performed in the last several years focused on sampling the diversity of human microbiota and providing the first large-scale comparative analyses (Eckburg et al., 2005; Ley et al., 2005; Bik et al., 2006; Gill et al., 2006; Frank et al., 2007). While several different methods have been utilized for this purpose (Sekirov et al., 2010), next-generation sequencing (NGS) and oligo microarray interrogation of the 16S rRNA gene pool are the primary current approaches that enable high-throughput detection and enumeration of community members. While NGS advantages include new sequence discovery and ability to interrogate any microbial community independent of its source, the advantages of phylogenetic oligo microarrays comprise the ability to conduct direct comparisons of each probeset signal among samples, making data fully quantitative, the ability to interrogate single samples independently, and as of now still lower costs, especially when comparing the total sample signal recorded by microarray vs. the total number of reads in NGS per sample. In recent years, several custom microarrays have been designed and successfully utilized for the characterization of complex human-associated microbial communities including PhyloChip and HITChip (Wang et al., 2004; Brodie et al., 2006; Palmer et al., 2006; Rajilic-Stojanovic et al., 2009), and the Microbiota Array developed in our laboratory (Paliy et al., 2009).

While the intestine in a newborn contains no microorganisms (it is essentially sterile), immediately after birth, it becomes colonized by enterobacteria, enterococci, lactobacilli, and bifidobacteria (Hopkins et al., 2005; Bjorkstrom et al., 2009; Enck et al., 2009). The initial microbial gut colonization is dependent on the mode of infant delivery. While traditional delivery sees the early gut microbiota resembling that of the vaginal canal, infants born by cesarean delivery are exposed to a different bacterial milieu closely related to that of the human skin (Biasucci et al., 2008, 2010; Dominguez-Bello et al., 2010). Later, the establishment of microbiota population is dependent on whether the infant is breast or formula fed (Harmsen et al., 2000). The gut biota of breast milk-fed babies is usually dominated by Bifidobacterium and Lactobacillus, whereas the microbiota of formula-fed babies is more complex and is more similar to that of an adult, with increased counts of Clostridia (Hopkins et al., 2005; Penders et al., 2005). Gradual changes in microbiota composition occur during early childhood, with a general reduction in the number of aerobes and facultative anaerobes and an increase in the populations of obligate anaerobic species (Hopkins et al., 2005). Traditionally, it has been thought that between 1 and 2 years of age, the human gut microbiota start to resemble that of an adult (Palmer et al., 2007), which is dominated by species from phyla Firmicutes (predominantly class Clostridia; 50–70% total bacterial numbers), Bacteroidetes (10–30%), Proteobacteria (up to 10%), and Actinobacteria (up to 10%), with 90% believed to be obligate anaerobes (Eckburg et al., 2005; Lay et al., 2005; Ley et al., 2006). Although the intestinal microbiota has been catalogued extensively in adults, and many recent reports were published for children between 0 and 2 years of age (Hopkins et al., 2005; Palmer et al., 2007; Bjorkstrom et al., 2009), few studies were conducted with gut microbiota of older children (Hopkins et al., 2001; Chernukha et al., 2005; Enck et al., 2009; Paliy et al., 2009). When Hopkins et al. (2001) analyzed fecal samples of young children between 1 and 7 years of age, a less diverse microbiota was evident in child fecal samples, and the numbers of enterobacteria were higher than those in adults. A large-scale study by Enck et al. (2009) used conventional colony plating to assess the numbers of several bacterial genera in children between 0 and 18 years of age. Although there were significant shifts in genus abundances during the first 2 years of life, no noticeable changes were evident in children between 2 and 18 years of age including stable levels of Bifidobacterium and Lactobacillus (see Harmsen et al., 2000, for a description of the limitations of colony plating estimates). A pilot study carried out in our laboratory revealed that adult and preadolescent child samples could be easily separated by principal components analysis (PCA), indicating that gut communities were different between adults and older children (Paliy et al., 2009).

Based on the results of these initial experiments, we hypothesized that older children might harbor different intestinal microbial communities than what has been observed either in adults or in younger children. To this end, we utilized Microbiota Array to quantitatively compare fecal microbiota communities of healthy children of preadolescent and adolescent age group with that of healthy adults.

Materials and methods

Fecal sample collection and processing

Healthy children of preadolescent and adolescent age (n=22; age range: 11–18 years, average – 12.6 years; BMI – 19.4 ± 3.7; 10 males, 12 females) and adults (n=10, age range: 22–61 years, average – 34.3 years; BMI – 22.0 ± 2.7; four males, six females) were recruited at Dayton Children's Medical Center and Wright State University. All volunteers consumed a standard western diet. The study was approved by the Wright State University IRB committee. All volunteers had no gastrointestinal symptoms and did not consume antibiotics or probiotics for 3 months before sample collection. Fresh fecal samples were collected, immediately homogenized, and frozen as described (Paliy et al., 2009; Rigsbee et al., 2011).

Isolation of genomic DNA (gDNA) and hybridization to Microbiota Array

Total gDNA was isolated from 150 mg of fecal samples using the ZR Fecal DNA Isolation kit (Zymo Research Corporation) according to the manufacturer's protocol. gDNA was amplified with phylogenetically conserved primers targeting the full-length 16S rRNA gene as described (Paliy et al., 2009; Rigsbee et al., 2011). PCR reaction conditions were 250 ng starting DNA, 25 cycles of PCR amplification, and total reaction volume – 50 μL. At least four separate PCR reactions were pooled together, fragmented, and then hybridized to Microbiota Array. The pooling of four separate PCR reactions was sufficient to eliminate PCR drift (Rigsbee et al., 2011). Microarray hybridization, washing, and scanning were carried out as described previously (Paliy et al., 2009; Rigsbee et al., 2011). Two replicate arrays were run for each fecal sample, each utilizing independently amplified 16S rRNA gene pools. The replicate arrays produced highly correlated signal values for most probesets – an average Pearson's correlation between replicates was 0.94.

Microarray data analysis

Raw microarray data were analyzed through a previously developed analysis pipeline as described (Rigsbee et al., 2011). Data adjustments and signal summations were carried out in the Microsoft excel template file that accounted for signal intensity, probeset presence call, cross-hybridization reduction, and 16S rRNA gene copy number (Rigsbee et al., 2011). To obtain the average relative abundances of phylogenetic groups among all samples within each sample type, the weighted mean around median was calculated in order to reduce the effect of outliers on the mean estimate. The weighted mean was calculated as (Wilcox, 2003) Embedded Image where Embedded Image is the weighted mean value for a given phylogenetic group, xi is the relative abundance of the phylogenetic group in sample i, and wi is the corresponding weight. Weights were calculated generally following the biweight approach of robust statistics as Embedded Image where m is the median of distribution and Embedded Image is the median absolute deviation Embedded Image. This biweight calculation insured that values close to the median had weights of 1, and the value contribution to the weighted mean was reduced as the distance from the median increased. No weight was lower than 0.5.

The SE of the weighted mean Embedded Image was calculated as (Wilcox, 2003) Embedded Image

Statistical tests of significance of observed differences were carried out in spss v16 (SPSS Inc.) and in carmaweb (Rainer et al., 2006) using moderated T-test and rank correlations. Nonparametric rank tests were used for correlation analyses because the datasets were not found to be normally distributed. Principal components analysis (PCA) was performed in matlab (The Mathworks Inc.).

Phylogenetic analyses were carried out on the fast unifrac web server (Hamady et al., 2010). Phylogenetic trees were constructed in bosque (Ramirez-Flandes & Ulloa, 2008). The results of the weighted and unweighted unifrac analyses were displayed through principal coordinates analysis (PCoA) as described (Hamady et al., 2010).

Quantitative real-time PCR (qPCR)

qPCR was carried out on ABI Prism 7000 Sequence Detection System using Perfecta SYBR Green qPCR mastermix (Quanta Biosciences Inc.) essentially as described (Paliy et al., 2009) including mathematical calculations taking into account unequal amplification efficiency for different primer pairs. Group-specific primers were designed as described (Rigsbee et al., 2011).

Results

Profiling of distal gut microbiome of healthy adolescent children and adults

The goal of this study was to assess whether the intestinal microbiota of teenage children resembles that of adults or whether any specific differences can be detected. Fecal samples representing members of human distal gut microbiota (Gill et al., 2006) were collected from 22 healthy preadolescent and adolescent children (designated kHLT for ‘healthy kids’). A smaller group of healthy adult volunteers was used for comparison (designated aHLT). Microbial populations in the fecal samples were analyzed by hybridizing 16S rRNA gene-enriched gDNA to the Microbiota Array, which provided both detection and quantitative measurements of 775 different species of human gut microbiota. The number of detected species varied appreciably from sample to sample, with a minimum of 209 and 276 species detected in adolescent and adult samples, respectively, and a maximum of 407 and 385 species, respectively. Although we detected slightly higher number of species in adult samples on average, the difference was not statistically significant (306 and 330 species for kHLT and aHLT groups, respectively).

PCA of the species abundance data set was used to test whether we could separate samples that belonged to different groups. As is shown in Fig. 1, PCA achieved a relatively good separation between adolescent and adult groups. Among adolescent samples, two significant outliers were observed. Examination of the bacteria detected in those samples revealed that both harbored an unusually high abundance of members of genus Holdemania (class Mollicutes; 5.1% and 6.0% relative abundance in outlier samples; the weighted mean for adolescent group was 1.0%). The main separation between kHLT and aHLT samples was observed in the principal component 1 (PC1) dimension. The PC1 values for each sample correlated negatively with the age of the corresponding participant; the negative correlation was highly significant (Spearman's rank correlation Rs=−0.55, P<0.01). This can be interpreted to imply that participant's age was one of the main factors differentiating different fecal samples in the PCA space. We have also used both unweighted as well as weighted unifrac analyses (Hamady et al., 2010) to assess whether similar separation could also be achieved when phylogenetic information about detected sequences was taken into account. Both types of unifrac analyses did not produce a clear separation of kHLT and aHLT samples in the PCoA plots (data not shown). The observed lack of concordance between PCA and PCoA results might indicate that the differences between kHLT and aHLT samples are caused by bacterial species and genera spread out among different phylogenetic groups rather than all clustered within one family, order, or class.

1

Results of the PCA of the microarray dataset. PC1 values for each sample are plotted on the x-axis and PC2 values are plotted on the y-axis. aHLT, healthy adult samples (n=10); kHLT, healthy adolescent samples (n=22).

Distribution of class and genera abundances among samples

Figure 2 displays a distribution of the relative abundances of different bacterial classes among samples and groups. As expected, Clostridia was by far the most abundant class in all 32 samples profiled, with an average of 72.2% and 74.4% of the total sample abundance among kHLT and aHLT samples, respectively. Actinobacteria (10.3% and 6.9%), Bacteroidetes (8.1% and 8.7%), and Bacilli (4.4% and 3.6%) were also well represented. Different classes of Proteobacteria combined to account for 2.4% and 2.2% total abundances in the kHLT and aHLT samples, respectively. For most classes, similar average numbers of detected species and relative class abundance were observed between kHLT and aHLT groups, with two exceptions. Members of Actinobacteria were more abundant among adolescent samples, whereas adult samples contained more Betaproteobacteria. However, neither of these differences was statistically significant due to the considerable variability of each class' abundance from sample to sample within each group. Indeed, as can be surmised from Fig. 2, a general trend of higher sample-to-sample class abundance variability that corresponded to lower average class abundance was evident in our experiments. This trend was found to be statistically significant as assessed by both the Spearman's rank and Kendall's tau nonparametric tests (Rs=−0.71, P=0.01; τ=−0.54, P=0.01). One possible explanation for this trend is that abundant classes are indispensable for the normal functioning of human intestinal microbiota in most individuals, whereas low abundance classes can vary from person to person based on a particular diet or environment.

2

Distribution of microbial species detection and abundance among adolescent and adult samples. The main figure displays a group scatter plot, with the x-axis listing different bacterial classes and the y-axis representing relative combined abundance of the members of each class in individual samples. Note that the y-axis is shown in a linear scale separated into three segments. aHLT, healthy adult samples (triangle; n=10); kHLT, healthy adolescent samples (diamond; n=22). Horizontal bars represent the weighted means of the abundance of each class among all samples of a particular group. Three sets of values below the scatter plot show (1) the number of samples in which members of each class were detected (‘Samples’); (2) an average number of detected species of each class per sample (‘Avg det’); and (3) an average combined relative abundance (weighted mean) of class members in each sample (‘Avg abnd’).

The distribution of relative abundances among different bacterial genera is shown in Fig. 3. Similar to the spread of class abundances, the fecal microbiota was dominated by relatively few genera. The top 12 abundant genera contributed between 72% and 76% of the total microbial abundance in both kHLT and aHLT samples. Genus Ruminococcus was by far the most abundant, with about 21–22% of the total abundance in all samples. This genus also contained the largest number of detected species (73 species per sample on average). Among the dozen most abundant genera, nine belonged to class Clostridia (see Fig. 3); the remaining three genera were Streptococcus, Bifidobacterium, and Bacteroides. The relative abundance values for these 12 genera were generally similar between adolescent and adult sample groups for all genera, except Clostridium and Bifidobacterium, both of which displayed a considerably higher level in kHLT samples. In addition, statistically significant differences between kHLT and aHLT sample groups were observed for several low-abundance genera as we describe below.

3

Distribution of genus relative abundances among samples. Different experiments are plotted as columns; 115 analyzed genera are plotted as rows. The relative abundances of each genus are plotted using a gradient scale as shown in the legend. aHLT, healthy adult samples; kHLT, healthy adolescent samples. Vertical line separates aHLT and kHLT samples. The 12 most abundant genera are shown on the right side; numbers represent the relative average abundance of each genus (weighted mean) in aHLT and kHLT samples, respectively. Genus assignments to four most abundant phyla are shown on the left side of the image. Sample kHLT07 lacking members of Bacteroidetes is highlighted.

Plotting genus-level relative abundances for all samples together revealed several specific cases of unusual genus abundance distribution. For example, Fig. 3 highlights fecal sample kHLT07 that lacked members of class Bacteroidetes almost entirely (total abundance <0.01%). This sample was taken from a healthy 13-year-old female adolescent with no indication of gastrointestinal disorders or any other unusual gut physiology. No other genus was found to be uniquely increased in its abundance in this sample to take the place of Bacteroidetes in the community. Rather, the members of several other genera such as Haemophilus, Anaerostipes, Veillonella, Holdemania, and Streptococcus were disproportionately more abundant in this sample compared with their average abundance among other samples.

On the other hand, the detection of microbial species was very similar between sample groups at all levels from genus to phylum. Comparison of the average species detection between kHLT and aHLT groups at the genus level produced a highly statistically significant correlation (Spearman's Rs=0.87, P<0.001; Kendall's τ=0.95, P<0.001), indicating that it was not the presence or absence of particular species or genera that differentiated kHLT and aHLT samples, but rather the relative abundances of the species within each sample group.

Core microbiome

Microbiota Array contains probesets for 775 different intestinal bacterial species. Among these, we were able to detect 647 in at least one sample out of combined 32 samples. Most of these species (543) belonged to the ‘shared’ category of those common to multiple, but not all samples. A total of 58 species were uniquely present in only a single sample. We have also identified a ‘core’ of 46 species shared among all analyzed fecal samples. This core microbiome was dominated by genus Ruminococcus represented by 27 species; Faecalibacterium and Roseburia contributed four core species each. In fact, 44 out of 46 core species belonged to class Clostridia (the other two were Streptococcus species from class Bacilli), again showing the dominance of this class in the fecal microbiota samples examined. Two species of genus Bacteroides, Bacteroides ovatus and Bacteroides thetaiotaomicron, were present in 31 out of 32 analyzed samples, with the sole exception being sample kHLT07, which lacked all members of the Bacteroidetes class (see Fig. 3). Core species factored considerably in the overall microbiota composition – on average, they contributed 24.8% to the total microbiota abundance (range: 17.5–34.4%), even though they only constituted 1/7 fraction of detected species in each sample on average.

Differences in genus relative abundance among groups

Because Microbiota Array provides direct quantitative measurements of 16S rRNA gene amounts and thus directly estimates the abundance of each of the 775 microbiota species, we were able to quantitatively compare genus abundances among the samples examined in this study. Among 115 different genera that microarray can detect, 15 showed substantial differences in their abundance levels between kHLT and aHLT groups (Table 1). The most prominent difference detected was a statistically significantly higher abundance of Bifidobacterium members among adolescent samples (9.0% and 5.4% relative abundance in kHLT and aHLT samples, respectively; P=0.05). This difference was due to an increase in the relative abundance of genus members rather than their presence or absence – the number of Bifidobacterium species detected per sample was practically identical for both groups (6.5). Specifically, both Bifidobacterium longum and Bifidobacterium catenulatum species were more abundant among adolescents (2.9% and 1.2% average abundance in adolescents, respectively; 0.9% and 0.7% in adults, respectively), whereas Bifidobacterium adolescentis had similar abundance in kHLT and aHLT samples (2.0% vs. 1.8% average abundance). Other notable differences included a higher abundance of Clostridium and a lower abundance of Prevotella and Sutterella in adolescent samples. Consistent with previous reports (Iizuka et al., 2004), Clostridium difficile (opportunistic human pathogen causing diarrhea) was detected in most samples at a low abundance level (kHLT samples: 18 out of 22 positive, 0.007% average abundance; aHLT samples: 4 out of 10 positive, 0.003% average abundance). Note that because of the significant variability in the relative genus abundance from sample to sample, not all differences were statistically significant as shown in Table 1. Members of genus Escherichia were detected significantly more often in adolescent than in adult samples (15 out of 22 had detectable levels in kHLT samples; 2 out of 10 – in aHLT samples), although the average abundance level was too low in both groups for us to assess relative abundance difference with confidence.

View this table:
1

Genera present at different relative abundances in samples from adults and adolescent children

GenusCorresponding classAdult samplesChild samplesLog ratioP value
No.%No.%
OxalobacterBetaproteobacteria20.2 ± <0.12<0.1 ± <0.14.90.13
SutterellaBetaproteobacteria50.8 ± 0.250.2 ± <0.12.2 0.02
LimnobacterBetaproteobacteria100.2 ± 0.1200.4 ± 0.1−1.10.14
EnterobacterGammaproteobacteria6<0.1 ± <0.15<0.1 ± <0.12.4 0.04
ClostridiumClostridia101.8 ± 0.3222.8 ± 0.2−0.6 0.01
ParasporobacteriumClostridia90.1 ± <0.115<0.1 ± <0.12.7<0.01
ButyrivibrioClostridia80.1 ± 0.16<0.1 ± <0.14.9 0.01
PeptococcusClostridia40.3 ± 0.16<0.1 ± <0.12.9 0.02
HespelliaClostridia60.2 ± <0.14<0.1 ± <0.13.9<0.01
TuricibacterClostridia90.2 ± 0.1220.4 ± 0.1−1.3 0.05
SlackiaActinobacteria60.1 ± <0.12<0.1 ± <0.13.4 0.01
BifidobacteriumActinobacteria95.4 ± 1.3229.0 ± 0.8−0.7 0.05
PrevotellaBacteroidetes50.7 ± 0.3120.2 ± 0.11.60.18
AnaerophagaBacteroidetes80.3 ± 0.1120.1 ± <0.11.10.24
VictivallisLentisphaerae40.2 ± 0.170.1 ± 0.11.60.15
  • * Number of samples in which members of this genus were detected out of 10 adult samples and 22 adolescent samples.

  • Percent contribution of species in each bacterial genus to the total hybridization signal measured by Microbiota Array. Data are shown as weighted mean ± weighted SE.

  • Logarithm with base 2 of the ratio of relative abundances of adult and adolescent samples.

  • ** Statistical significance of observed differences between sample groups as measured using the moderated T-test. P values ≤0.05 are underlined.

Validation of array results with qPCR

We have used qPCR to validate the findings revealed by the Microbiota Array. Primers developed specifically to amplify 16S rRNA gene from members of three separate genera were used in qPCR tests (Rigsbee et al., 2011). As can be assessed in Table 2, the qPCR results matched microarray data well, with only a slight deviation observed for the Faecalibacterium genus abundance estimate of adult sample aHLT05.

View this table:
2

Validation of microarray results with quantitative PCR

Bacterial groupAdult sample aHLT02Adult sample aHLT05
qPCRArrayqPCRArray
Faecalibacterium11.7%± 5.0%11.0%± 0.3%4.9%± 3.5%9.1%± 0.2%
Bifidobacterium5.7%± 2.2%6.4%± 0.3%0.3%± 0.2%0.1%± <0.1%
Collinsella1.0%± 0.3%0.7%± 0.1%0.1%± 0.1%0.3%± <0.1%
  • * Primers were developed to cover most of the 16S rRNA gene sequences in the RDP database belonging to each group (Rigsbee et al., 2011).

  • Numbers represent the relative abundances of gDNA of each genus in the total sample; no 16S rRNA gene copy adjustment was applied to the normalized values.

  • Data are shown as arithmetic mean ± SE.

Discussion

The present study is the first to assess distal gut microbiota of adolescent children using a high-throughput technology. Although it is generally assumed that older children house gut microbial populations similar to those of adults, our results indicate that such gut communities are sufficiently different (see Fig. 1) in the population examined. Few distinctions were found when the relative abundances were compared at the class level; however, many differences were observed at the genus level between adolescent and adult fecal microbiota. Most strikingly, we found that the relative combined abundance of members of genus Bifidobacterium was statistically significantly higher in adolescent children (close to twofold difference). This high prevalence of Bifidobacterium is well known for children of a young age; in fact, bifidobacteria are some of the first species to colonize newborn gut (Bjorkstrom et al., 2009; Enck et al., 2009). However, an increased prevalence of this genus was not recognized previously among preadolescent and adolescent age groups. Importantly, the level of bifidobacteria correlated negatively with age in our sample set (Spearman's rank Rs=−0.29, P=0.03; Kendall's τ=−0.36, P=0.05). Thus, our data are more consistent with the model where levels of bifidobacteria in children decrease gradually between 2 and 18 years of life until reaching stable levels in the early adulthood, rather than with the model that assumes that the levels of these bacteria decline quickly after weaning.

Considering distal gut microbiota of healthy adolescent children and adults on a global level, microbial populations were dominated by class Clostridia in all cases. In particular, Ruminococcus was the most prevalent genus in all the samples examined, which is consistent with an important role of the members of this genus as the primary carbohydrate degraders in the human gut (Flint et al., 2008). Most of the highly abundant genera and classes showed relatively consistent abundances among samples, whereas the levels of low abundance genera and classes varied widely from sample to sample (see Figs 2 and 3).

We have identified a core set of 46 bacterial species that we detected in every sample examined; as expected, this core set was dominated by members of Clostridia and genus Ruminococcus in particular, which is in general agreement with another study that was focused on core microbiome identification (Tap et al., 2009). It is tempting to speculate that Ruminococcus serves as the ‘irreplaceable’ primary degrader of complex carbohydrates in the ecological food chain in the gut, which would explain its very consistent level of high abundance in all samples. At the same time, functions of the secondary and tertiary degraders can be performed (interchangeably) by many different bacterial species and genera, leading to a wider spread of their abundances among different samples (Flint et al., 2008; Karasov & Carey, 2009). Future mechanistic studies using artificial gut fermentor models and engineered microbiota populations in gnotobiotic mice will be required to shed light on this hypothesis.

Supporting Information

Table S1. Abundance of bacterial classes among kHLT and aHLT samples.

Table S2. Abundance of bacterial genera among kHLT and aHLT samples.

Table S3. Detection of bacterial species among kHLT and aHLT samples.

Table S4. Detection of bacterial species among kHLT and aHLT samples: class level distribution.

Figure S1. Distribution of genus relative abundances among samples.

Appendix S1. List of qPCR primers used.

Please note: Wiley-Blackwell is not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article.

Acknowledgements

Our thanks are due to Vijay Shankar for valuable comments on the manuscript and to the members of the WSU Center for Genomics Research for access to the facility. This work was supported by the National Institutes of Health grants AT003423 and HD065575 to O.P.

References

View Abstract