- Research article
- Open Access
Age-related changes in gut microbiota composition from newborn to centenarian: a cross-sectional study
BMC Microbiology volume 16, Article number: 90 (2016)
It has been reported that the composition of human gut microbiota changes with age; however, few studies have used molecular techniques to investigate the long-term, sequential changes in gut microbiota composition. In this study, we investigated the sequential changes in gut microbiota composition in newborn to centenarian Japanese subjects.
Fecal samples from 367 healthy Japanese subjects between the ages of 0 and 104 years were analyzed by high-throughput sequencing of amplicons derived from the V3-V4 region of the 16S rRNA gene. Analysis based on bacterial co-abundance groups (CAGs) defined by Kendall correlations between genera revealed that certain transition types of microbiota were enriched in infants, adults, elderly individuals and both infant and elderly subjects. More positive correlations between the relative abundances of genera were observed in the elderly-associated CAGs compared with the infant- and adult-associated CAGs. Hierarchical Ward’s linkage clustering based on the abundance of genera indicated five clusters, with median (interquartile range) ages of 3 (0–35), 33 (24–45), 42 (32–62), 77 (36–84) and 94 (86–98) years. Subjects were predominantly clustered with their matched age; however, some of them fell into mismatched age clusters. Furthermore, clustering based on the proportion of transporters predicted by phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt) showed that subjects were divided into two age-related groups, the adult-enriched and infant/elderly-enriched clusters. Notably, all the drug transporters based on Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology groups were found in the infant/elderly-enriched cluster.
Our results indicate some patterns and transition points in the compositional changes in gut microbiota with age. In addition, the transporter property prediction results suggest that nutrients in the gut might play an important role in changing the gut microbiota composition with age.
The microbiota composition of the human gut changes with age, and alterations in this composition influence human health. In the early 1970s, culture-based methods were used to demonstrate that the gut microbiota composition changes during the aging process . Recent studies using molecular methods have also indicated clear differences in the composition of gut microbiota among infants, toddlers, adults and the elderly . After birth, the initial microbiota composition is affected by the mode of birth [3–5] and the mother’s gut microbiota [6, 7]. Subsequently, a significant shift in the composition of the gut microbe community occurs when the infant transitions to a more solid and varied diet. A recent report suggested that the age-affiliated microbiota population shifts from 3 days to 2 years after birth and that major differences are apparent between 2 years and adulthood . Other reports have indicated that the phylogenetic composition of the bacterial communities evolve towards an adult-like configuration within the 3-  or 4-year  period after birth. It has been recently shown that the gut microbiota is not yet established at 5 years of age . Another broad shift in gut microbe populations occurs later in life. However, almost all studies related to the gut microbiota of the elderly have been performed on subjects classified into segmented age groups based on varying definitions of ‘elderly’, such as ‘over 60’ , ‘over 65’ , ‘over 70’  or ‘centenarian’ . It is unclear when and how the microbiota composition shifts from the adult stage to the elderly stage. Yatsunenko et al.  conducted a large study with subjects aged 0–83 years and revealed the sequential changes that occur with age. Their report provided important insights, such as the period required to form an adult-like gut microbiota, greater between-subject variation among children than adults, differences in the phylogenetic composition of gut microbiota among individuals from different countries and an increase in bacterial diversity with age. Nevertheless, the sequential changes that occur in the elderly remained unclear due to the limited number of subjects older than 60 years. Furthermore, although abundant data on gut microbiota composition are available in some public databases, sequential changes cannot be evaluated with these public data because of biases stemming from differences in study methods, especially in DNA extraction [16, 17], and nationality differences among the subjects [9, 18–22].
We identified the sequential changes in gut microbiota composition in Japanese subjects over a wide age range, 0–104 years. Our results provide new insights into the developmental period for gut microbiota composition and the patterns of change with age.
Overview of gut microbiota composition in each age group
A total of 1,839,703 high-quality paired sequences were obtained from the 371 samples, with 4,959 ± 1,813 (average ± standard deviation) reads per sample, which were clustered into 5,952 OTUs and classified into 187 bacterial groups at the genus level (186 genera and one unidentified group). We first calculated UniFrac distances to determine the extent of similarity between microbial communities. UniFrac PCoA (principal co-ordinate analysis) of 5,952 OTUs indicated that age explained the variation in our data set using both weighted and un-weighted analyses (Fig. 1). No gender differences were observed (Additional file 1).
The phylum and genus compositions of gut microbiota in each age group are shown in Fig. 2, Additional file 2 and Additional file 3. In agreement with previous results, the microbiota composition included four predominant phyla. The relative abundance of Actinobacteria substantially decreased after weaning and continued to decrease with age. Firmicutes was the most predominant phylum after weaning but was less abundant in children younger than 4 years compared with subjects older than 4 years. Increases in the relative abundance of Bacteroidetes and Proteobacteria were observed in subjects over 70 years old. The relative abundance of Bacteroidetes did not change sequentially, but a stepwise increase was observed beyond 70 years of age. The change in the relative abundance of Proteobacteria was opposite that of Firmicutes.
Age-related changes in genera and their correlations
To explore the changes in gut microbiota with age in detail, we calculated the co-abundance associations of genera and then clustered the correlated genera into nine CAGs (Additional file 4), which describe the significant differences in microbiota structure among subject groups (Permutational MANOVA, p < 0.001). The transition from infant to centenarian was accompanied by distinctive CAG dominance, with a significant abundance of Bacteroides, [Eubacterium] and Clostridiaceae CAGs (elderly-associated CAGs); Enterobacteriaceae CAGs (infant and elderly-associated CAGs); Bifidobacterium CAGs (infant/child-associated CAGs); and Lachnospiraceae CAGs (adult-associated CAGs) (Fig. 3). Megamonas and Peptoniphilus CAGs were relatively enriched in the elderly. Dorea CAG abundance appeared unrelated to aging. Sequential changes occurred in the relative abundance of Bacteroides, Lachnospiraceae and Bifidobacterium CAGs in the gut microbiota during childhood and adolescence. Not all CAGs were composed exclusively of related species, as shown in Additional files 5 and 6. Wiggum plots showed the relative abundance of each genus and significant associations between nine CAGs (Fig. 4, Additional file 5). Among the 186 genera, 116 had associations with other genera with an absolute coefficient value >0.3 (Additional file 5 and 6 and Fig. 4). Almost all of the correlations were positive; the only four negative correlations were observed between Enterobacteriaceae and Lachnospiraceae, Enterobacteriaceae and Blautia, Bifidobacterium and Parabacteroides and Veillonella and [Mogibacteriaceae]. A greater number of positive correlations were observed in the elderly-associated CAGs, especially the Bacteroides and [Eubacterium] CAGs. We then performed in vitro assays to investigate some relevant relationships among genera. In accordance with the Wiggum plot results, both of Parabacteroides distasonis JCM 5825 and Bifidobacterium longum JCM 1217 growth were suppressed when co-cultivate with each other (Additional file 7). A similar negative relationship was observed between Escherichia coli JCM 1649, belong to family Enterobacteriaceae and Blautia producta JCM 1471. In contrast, Bacteroides uniformis JCM 5828 growth increased when co-cultivated with Parabacteroides distasonis JCM 5825.
As for gut microbiota diversity, all four alpha diversity scores based on PD whole tree, Chao1, the number of observed species and the Shannon index substantially increased after weaning and continued to increase sequentially until the twenties. These scores were stable during adulthood and then increased again at the elderly stage until the centenarian stage (Fig. 5).
Period of age-related change in community structure of gut microbiota
We performed hierarchical Ward’s linkage clustering based on the abundance of genus-like groups to predict the period of age-related reshaping of gut microbiota. Genus-like groups belonging to the same phylum did not comprise the same cluster, suggesting that the age-related pattern of change was different between the phylum and genus levels. The subjects were divided into three age clusters, for infant, adult and elderly cluster (Additional file 8). The relative abundance of Actinobacteria and Clostridia were significantly higher in infant and adult cluster, respectively (Additional file 9). The elderly cluster showed the significantly higher abundance of Bacteroidetes, Betaproteobacteria and Deltaproteobacteria (Additional file 9). When the subjects were divided into five age clusters, the median (interquartile range) age was 3 (0.5–35) in the infant cluster, 33 (24.75–45.5) in adult cluster I, 42 (32–62) in adult cluster II, 74 (34–81.5) in elderly cluster I and 93.5 (85–98) in elderly cluster II (Additional file 8). The age distributions were similar between the adult I and II clusters and between the elderly I and II clusters. The adult clusters and the elderly clusters showed significant differences in the relative abundance of 32 and 28, respectively, of the 61 genus-like groups, which were found in 50 % of the subjects in any cluster (Additional files 10 and 11). Subjects were predominantly clustered by age; however, some subjects, but not those in the elderly II cluster, fell into a mismatched-age cluster (Additional files 8 and 12).
Functional properties predicted by PICRUSt
We performed PICRUSt analysis to predict the relative abundance of transporter genes because an altered diet  and antibiotic treatment  have been reported to be among the most powerful factors that affect the gut microbiota. Clustering based on the relative abundance of the predicted transporters showed that subjects were divided into two age-related groups, the adult-enriched and infant/elderly-enriched clusters (Additional file 13). For example, there was a lower abundance of a predicted xylose transporter (KEGG module: M00215) in pre-weaned infants, probably reflecting the different dietary habits of subjects in each segmented age group (Additional file 14). Interestingly, all drug transporters based on KEGG Orthology groups were found in the infant/elderly-enriched cluster (Additional files 13 and 15).
The composition of gut microbiota is thought to change during the aging process ; however, few reports have utilized molecular techniques to investigate the long-term, sequential changes in gut microbiota composition. Our results are in agreement with those of recent studies indicating clear differences in gut microbiota composition among infants, adults and the elderly [2, 23]. The present study using LEfSe method indicated that Actinobacteria, Clostridia and Bacteroidetes, Betaproteobacteria and Deltaproteobacteria were representative taxa in infant, adults and the elderly cluster. Additionally, our results revealed the sequential changes that occur with age from newborns to centenarians. Furthermore, our results showed that Japanese adults (21–69 years old) have a greater abundance of the genera Blautia and Bifidobacterium (interquartile ranges (IQR) of 18 (12–24) % and 7 (2–14) %, respectively) and a relatively lower abundance of genera related to Bacteroidetes (IQR 4 (1–10) %), compared with those reported in previous studies in other nations. For example, the estimated abundance of Blautia, Bifidobacterium and Bacteroidetes were < 10 %, < 2 % and > 10 %, respectively, in US and Colombia , Korea, China and US  and Ireland . These patterns might be characteristic of the gut microbiota composition of the Japanese population, although they may also reflect the DNA extraction method [16, 17] and the amplified region of the 16S rRNA gene .
Although there are differences among individuals, our analysis of the phylum composition of gut microbiota in each age group showed a significant shift in the relative abundance of Actinobacteria in infants from before to after weaning. The compositional pattern of gut microbiota during childhood has been thought to impact health later in life [11, 25], but children older than 2 years have not been sufficiently investigated. Our data show that some genera belonging to Bacteroides, Lachnospiraceae and Bifidobacterium CAGs and the alpha diversity of gut microbiota continued to change sequentially with age in subjects younger than twenty, reflecting the human gut microbiota maturation process. However, children younger than 20 years fell into both the infant and adult clusters, regardless of their age, when clustered based on the abundance of genus-like groups, thus illustrating the individual differences in the gut microbiota maturation process.
In the present study, the Wiggum plot showed negative relationships between the relative abundance of Enterobacteriaceae, which creates an greater endotoxin challenge for the weakened intestinal barrier and thus results in increased stimulation of the inflammatory response , and the abundance of Blautia and Lachnospiraceae;g (Fig. 4), which belonged to the adult-associated CAG (Lachnospiraceae;g CAG). Other butyrate-producing bacteria, such as Coprococcus, Roseburia and Faecalibacterium, were also clustered in the same CAG. Furusawa et al. reported that microbial-derived butyrate regulates Treg cell differentiation in vitro and in vivo . Given the age-related reduction in the abundance of the genus Bifidobacterium, which down-regulates pro-inflammatory responses in the gut epithelium [28–30], our results suggested that the aging-related dysbiosis in elderly subjects may be a contributing factor to inflammatory responses that occur with advancing age.
We performed in vitro assays to investigate bacterial interactions, including those between Enterobacteriaceae and Blautia. These results were in accordance with the Wiggum plot results. However, these results might have some biases from the different of environmental condition between in human gut and in vitro assay. In addition, it is uncertain that all bacterial interactions are consistent with the relationships in the Wiggum plot, because genera with positive relationships in the Wiggum plot might grow well under the same environmental conditions without a mutualistic relationship. In contrast, some genera combinations have been reported to exhibit mutualistic relationships,although the absolute values of the correlation coefficients were below 0.3 (no visible relationship in the Wiggum plot). Pande et al. revealed that Acinetobacter baylyi and Escherichia coli reciprocally exchange essential amino acids . It has also been reported that Bifidobacterium populations can be stimulated efficiently with a concomitant decrease in Enterobacteriaceae [32, 33]. Acetate, one of the main fermentation products of Bifidobacterium, was reported to promote the growth of butyrate-producing bacteria and the in vitro production of butyrate [34, 35]. Furthermore, Bifidobacterium longum has been reported to alter gut luminal metabolism via interactions with Bacteroides caccae and Eubacterium rectale . Considering these reports, the correlations between gut microbiota members might be more complicated than shown in the Wiggum plot. Our computational analysis results must be interpreted cautiously because they are based on a limited data set. An advanced culture method is needed to clarify the relationships among gut microbiota.
A wide diversity of microorganisms is needed to utilize the many nutrients in adult diets . In addition, a low gut microbiota diversity has been associated with an increasing number of conditions, such as autism , autoimmune disease  and obesity . Maintaining sufficient bacterial richness and diversity is important for providing gut microbiota with functional redundancy, adaptability and thus systematic robustness against environmental challenges . In this study, we observed an increase in gut microbiota diversity with aging until the centenarian stage. Claesson et al. reported that the alpha diversity of the gut microbiota in community dwellers was significantly higher than that of people in long-stay care . Therefore, in the present study, the observed increase in microbiota diversity with age was likely due to the inclusion of community-dwelling elderly participants. In contrast, centenarians in a Chinese longevous village population had a more diverse gut microbiota than did younger elderly aged 85–99 years . It is uncertain why the trend in diversity differed between Chinese and Japanese centenarians.
Biagi et al.  showed a significantly compromised gut microbiota in centenarians but not in elderly subjects aged approximately 70 years compared with a group of younger adults. The authors therefore suggested that a healthy gut microbiota community might be affected by aging-related physiological and behavioral changes that occur after 70 years of age, which was considered the threshold age for defining an individual as elderly. In agreement with this hypothesis, our clustering data suggested that most subjects over 70 years of age comprised two elderly-type clusters (Additional file 8). Nevertheless, the distribution of subjects was not clear from the classification based on subject age. This result suggested that other factors besides age contribute to the composition of gut microbiota communities.
Our study did not allow us to determine why gut microbiota changed with age because we lacked lifestyle and dietary habit information on the subjects. To begin to address this question, we performed PICRUSt analysis to predict the relative abundance of transporter genes because we hypothesized that bacterial transporters that incorporate nutrients from the gut environment would differ if dietary habit was a predominant contributor to gut microbiota composition. Clustering based on the relative abundance of the predicted transporters showed that subjects were divided into two age-related groups, the adult-enriched and infant/elderly-enriched clusters, implying that nutrients in the gut evoke the different gut microbiota compositions between adult and infant/elderly subjects. Not all transporters can explain the relationship with the dietary habits of subjects in each segmented age group, but, for example, the increase in the relative abundance of predicted xylose transporters in subjects after weaning seems to reflect the change from mother’s milk to an omnivorous diet. Interestingly, all the drug transporters based on KEGG Orthology groups were found in the infant/elderly-enriched cluster, perhaps due to the frequent antibiotic treatment of infants and the elderly compared with adults.
In addition, elderly people are known to have decreased intestinal function relative to younger people, which affects digestion, nutrient absorption and immune activity [45, 46] and may also impact the microbiota composition [47, 48]. We found that certain oral bacteria, such as Porphyromonas, Treponema, Fusobacterium and Pseudoramibacter, which have difficulty reaching the intestinal tract due to barriers such as gastric juice and bile acid, were enriched in the elderly-associated CAGs (Additional file 6). Therefore, the decline in gastrointestinal tract functionality in the elderly may also lead to significant changes in gut microbiota.
In conclusion, we provide a description of the changes in gut microbiota with age, thus illustrating the long trajectory throughout human life. Our results indicate some patterns and transition points in gut microbiota composition with age. The gut microbiota in subjects younger than 20 years changed with age as it matured, and that of subjects older than 70 years changed again into the elderly type. In addition, the transporter property prediction results suggested that nutrients in the gut might play an important role in changing the gut microbiota composition with age.
Our findings help clarify the gut microbiota composition in a healthy population at each age period. Further analyses investigating lifestyle traits or prospective cohorts focused on subjects who appear to have a gut microbiota typical of an age group older than their matched age would be valuable for revealing the relationships between gut microbiota and host health, including the aging process.
Fecal samples were collected from a total of 367 community-dwelling Japanese volunteers (one sample per subject, except for two samples from one boy and one girl collected at preweaning and weaning and three samples from one girl at preweaning, weaning and 5 years of age) between 0 and 104 years of age (157 men and 210 women). Subjects over 80 years of age were directly recruited by the authors to confirm that they were community dwellers. The distribution of subjects according to age and individual data are shown in Table 1 and Additional file 12. No significant differences in gender distribution were observed among the age groups (Fisher’s test, p = 0.997). Fecal samples were collected from subjects that participated in three different studies. Two study protocols were for the collection of feces from subjects aged 21–65 years old or from community-dwelling elderly individuals and were approved by the Local Ethics Committee of the nonprofit organization Japan Health Promotion Supporting Network (Wakayama, Japan). The third study protocol was for the collection of feces from subjects aged 0–100 years old and was approved by the ethics committee of Kensyou-kai Incorporated Medical Institution (Osaka, Japan). Written informed consent was obtained from all subjects or from their legal guardians or relatives.
To identify the sequential changes in gut microbiota composition that occur with age, subjects were divided into 10-year age groups, except for subjects aged less than 10 years, who were divided into four groups: preweaning, weaning, weaned to 3 years old and 4–9 years old.
Fresh fecal samples were transferred by the subjects into tubes and immediately enclosed in plastic bags containing AnaeroPouch (Mitsubishi Gas Chemical, Tokyo, Japan) to create an anaerobic environment.
The fecal samples collected from subjects younger than 80 years were stored at −20 °C for three days at most and were transported to the laboratory by logistics companies. The samples collected from subjects older than 80 years were stored at −20 °C and then delivered to the laboratory by the study authors within one day. Immediately upon receipt, the fecal samples were stored at −80 °C until the day of analysis.
A total of 20 mg of each fecal sample was collected from three regions (upper, middle and lower) and mixed well, and DNA was extracted using the bead-beating method as previously described . After centrifugation at 14,000 × g for 5 min, 400 μl of the supernatant was extracted with phenol-chloroform, and 250 μl of the supernatant was precipitated with isopropanol. Purified DNA was suspended in 2,000 μl of Tris-EDTA buffer (pH 8•0).
Sequencing and data processing
The V3-V4 region of the bacterial 16S rRNA gene was then amplified by PCR with the TaKaRa Ex Taq HS Kit (TaKaRa Bio, Shiga, Japan) and the primer sets Tru357F (5′-CGCTCTTCCGATCTCTGTACGGRAGGCAGCAG-3′) and Tru806R (5′-CGCTCTTCCGATCTGACGGACTACHVGGGTWTCTAAT-3′). Each 1-μl sample of DNA, at a concentration of approximately 10–200 ng/μl as measured using a Nanodrop 2000 (Thermo Fisher Scientific, Waltham, MA, USA), was amplified in triplicate using the following protocol: preheating at 94 °C for 3 min; 20 cycles of denaturation at 94 °C for 30 s, annealing at 50 °C for 30 s and extension at 72 °C for 30 s; and a final terminal extension at 72 °C for 10 min. After verifying the amplified DNA based on PCR product size using the QIAxcel system (Qiagen, Valencia, CA, USA), the triplicate samples were combined. A 1-μl sample of the combined PCR products was amplified using the following barcoded primers adapted for Illumina MiSeq: Fwd 5′-AATGATACGGCGACCACCGAGATCTACACXXXXXXXXACACTCTTTCCCTACACGACGCTCTTCCGATCTCTG-3′ and Rev 5′-CAAGCAGAAGACGGCATACGAGATXXXXXXXXGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTGAC-3′, where X represents a barcode base. The DNA was amplified according to the protocol described above except that only 8 cycles were performed. After validating the 2nd amplified DNA product using the QIAxcel system, the PCR products were purified using a QIAquick 96 PCR Purification Kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s protocol. The purified products were quantified using a Quant-iT PicoGreen dsDNA Assay Kit (Life Technologies, Carlsbad, CA, USA). Subsequently, equal amounts of the amplicons from multiple samples were pooled, and primer dimers were removed using a GeneRead Size Selection Kit (Qiagen, Valencia, CA, USA). The pooled libraries were sequenced using an Illumina MiSeq instrument with a MiSeq v3 Reagent Kit (Illumina, Inc., San Diego, CA, USA).
After removing sequences consistent with data from the Genome Reference Consortium human build 37 (GRCh37) or PhiX 174 from the raw Illumina paired-end reads, the 3′ region of each read with a PHRED quality score of less than 17 was trimmed. Trimmed reads less than 150 bp in length with an average quality score of less than 25 or those lacking paired reads were also removed. The reads that passed the quality filters were combined using the fastq-join script in EA-Utils (version 1.1.2–537) .
Potential chimeric sequences were removed using UCHIME, assigned to operational taxonomic units (OTUs) using Open-reference OTU picking  with a 97 % threshold of pairwise identity, and then classified taxonomically using the Greengenes reference database (http://greengenes.secondgenome.com/downloads/database/13_5) . An advantage of using the open reference method is that it minimizes spurious hits of sequencing reads to taxa that are not present in the gut; as a result, the number of obtained OTUs was lower than that in previous reports.
The microbial diversity within each age-segmented group [alpha diversity including Chao1, number of observed species (the number of OTUs), phylogenetic distance (PD) whole tree, and Shannon diversity index] and the distances between subjects (UniFrac distance as beta diversity) were estimated using QIIME version 1.8.0 software.
Hierarchical analysis was performed using the hclust function in R package 3.2.1. Distances based on the squared Euclidean distance were calculated for input into an agglomerative algorithm through Ward’s method. The population densities (z-scores) of genera scaled by color are displayed together with a dendrogram of bacterial genera in a heat map.
Bacterial co-abundance groups (CAGs)
All genus-like-level groups, except for an unidentified group, were entered into this analysis. CAGs were defined by a heat plot showing Kendall correlations between genera clustered by Pearson’s correlation coefficient and Ward’s linkage hierarchical clustering in R using the Made4 package  as previously described . The transition type of each CAG with aging indicates the sum of z-scores converted from the relative abundance of genera belonging to the same CAG. Network plots highlighting correlative relationships were visualized using Cytoscape version 3.2.1 . These associations were controlled for multiple testing using the q-value method, and only those with a false discovery rate <0.05 were retained. The cut-off for line in the Wiggum plots was set at an absolute coefficient value of greater than 0.3.
Bacterial strains and culture conditions
Cultivable strains were selected based on the Wiggum plot results and were obtained from public culture collections (Additional file 16). All the strains were pre-cultured at 37 °C for 16 h under anaerobic conditions in Gifu Anaerobic Medium (GAM) broth (Nissui Seiyaku Co. Ltd., Tokyo, Japan). Then, approximately 1 × 108 cell in the pre-cultures, which were calculated by a density of McFarland, were added to GAM broth and incubated at 37 °C for 5 h under anaerobic conditions. After this culture period, the microorganisms were collected by centrifugation at 10,000 × g and DNA was extracted by the same method as described in the DNA extraction section. Experiments were performed in triplicate.
Real-time PCR for quantitative determination of cell number
Real-time PCR was performed using an ABI PRISM 7500 Fast Real-time PCR System (Life Technologies, Carlsbad, CA, USA) with SYBR Premix Ex Taq (TaKaRa Bio Inc, Shiga, Japan). The primer sets are shown in Additional file 17 [57–60]. The amplification program consisted of one cycle of 94 °C for 10 s, followed by 40 cycles of 94 °C for 5 s, the appropriate annealing temperature for 30 s and 72 °C for 30 s. Fluorescent products were detected during the last step of each cycle. Melting curves were obtained by heating from 60 to 95 °C in 0.2 °C/s increments with continuous fluorescence collection.
Phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt) analysis
PICRUSt analysis was performed to predict the relative abundance of transporter genes . Independent of the taxonomic analysis, 97 % of the OTUs were picked using a closed-reference OTU picking protocol (QIIME 1.8.0 [51, 52]) and the Greengenes database (database/13_8)  pre-clustered at 97 % identity. The obtained OTU table was normalized by 16S rRNA copy number, and functional genes were predicted from the Kyoto Encyclopedia of Genes and Genomes (KEGG) catalogue .
The gender distribution of subjects and intergroup differences at the genus level in each subcluster were analyzed by Fisher’s test and the Mann-Whitney U test, respectively, using SPSS version 23.0 statistical software (IBM, Armonk, NY, USA). Intergroup differences at the phylum, class, order, family and genus level in each cluster were analyzed by the linear discriminant analysis (LDA) effect size (LEfSe) method  with default settings on website (https://huttenhower.sph.harvard.edu/galaxy/root). Permutational MANOVA  was performed to test for significant differences in CAGs using the vegan package in R. Values of p < 0.05 were considered statistically significant.
DNA sequences corresponding to the 16S rRNA gene data have been deposited in DDBJ under accession number DRA004160.
operational taxonomic units
polymerase chain reaction
quantitative insights into microbial ecology
ribosomal ribonucleic acid
Mitsuoka T, Hayakawa K. The fecal flora in man. I. Composition of the fecal flora of various age groups. Zentralbl Bakteriol Orig A. 1973;223:333–42.
Conlon M a, Bird AR. The Impact of Diet and Lifestyle on Gut Microbiota and Human Health. Nutrients. 2014;7:17–44.
Palmer C, Bik EM, DiGiulio DB, Relman D a, Brown PO. Development of the human infant intestinal microbiota. PLoS Biol. 2007;5:e177.
Dominguez-Bello MG, Costello EK, Contreras M, Magris M, Hidalgo G, Fierer N, et al. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci U S A. 2010;107:11971–5.
Tsuji H, Oozeer R, Matsuda K, Matsuki T, Ohta T, Nomoto K, et al. Molecular monitoring of the development of intestinal microbiota in Japanese infants. Benef Microbes. 2012;3:113–25.
Grönlund M, Grześkowiak Ł, Isolauri E, Salminen S. Influence of mother’s intestinal microbiota on gut colonization in the infant. Gut Microbes. 2011;2(4):227–33.
Vaishampayan PA, Kuehl JV, Froula JL, Morgan JL, Ochman H, Francino MP. Comparative metagenomics and population dynamics of the gut microbiota in mother and infant. Genome Biol Evol. 2010;2:53–66.
Avershina E, Storrø O, Øien T, Johnsen R, Pope P, Rudi K. Major faecal microbiota shifts in composition and diversity with age in a geographically restricted cohort of mothers and their children. FEMS Microbiol Ecol. 2014;87:280–90.
Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, et al. Human gut microbiome viewed across age and geography. Nature. 2012;486(7402):222–7.
Bergström A, Skov TH, Bahl MI, Roager HM, Christensen LB, Ejlerskov KT, et al. Establishment of intestinal microbiota during early life: a longitudinal, explorative study of a large cohort of Danish infants. Appl Environ Microbiol. 2014;80:2889–900.
Cheng J, Ringel-Kulka T, Heikamp-de Jong I, Ringel Y, Carroll I, de Vos WM, et al. Discordant temporal development of bacterial phyla and the emergence of core in the fecal microbiota of young children. ISME J. 2015;10(4):1002–14.
Mueller S, Saunier K, Hanisch C, Norin E, Alm L, Midtvedt T, et al. Differences in fecal microbiota in different European study populations in relation to age, gender, and country: a cross-sectional study. Appl Environ Microbiol. 2006;72:1027–33.
Claesson MJ, Cusack S, O’Sullivan O, Greene-Diniz R, de Weerd H, Flannery E, et al. Composition, variability, and temporal stability of the intestinal microbiota of the elderly. Proc Natl Acad Sci U S A. 2011;108 Suppl 1:4586–91.
Mariat D, Firmesse O, Levenez F, Guimarăes V, Sokol H, Doré J, et al. The Firmicutes/Bacteroidetes ratio of the human microbiota changes with age. BMC Microbiol. 2009;9:123.
Biagi E, Nylund L, Candela M, Ostan R, Bucci L, Pini E, et al. Through ageing, and beyond: gut microbiota and inflammatory status in seniors and centenarians. PloS one. 2010;5(5):e10667.
Wesolowska-Andersen A, Bahl MI, Carvalho V, Kristiansen K, Sicheritz-Pontén T, Gupta R, et al. Choice of bacterial DNA extraction method from fecal material influences community structure as evaluated by metagenomic analysis. Microbiome. 2014;2:19.
Wagner Mackenzie B, Waite DW, Taylor MW. Evaluating variation in human gut microbiota profiles due to DNA extraction method and inter-subject differences. Front Microbiol. 2015;6:130.
Nakayama J, Watanabe K, Jiang J, Matsuda K, Chao S-H, Haryono P, et al. Diversity in gut bacterial community of school-age children in Asia. Sci Rep. 2015;5:8397.
Escobar JS, Klotz B, Valdes BE, Agudelo GM. The gut microbiota of Colombians differs from that of Americans, Europeans and Asians. BMC Microbiol. 2014;14:1–14.
Nam Y-D, Jung M-J, Roh SW, Kim M-S, Bae J-W. Comparative Analysis of Korean Human Gut Microbiota by Barcoded Pyrosequencing. PLoS One. 2011;6:e22109.
Schnorr SL, Candela M, Rampelli S, Centanni M, Consolandi C, Basaglia G, et al. Gut microbiome of the Hadza hunter-gatherers. Nat Commun. 2014;5:1–12.
De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci U S A. 2010;107:14691–6.
Ottman N, Smidt H, de Vos WM, Belzer C. The function of our microbiota: who is out there and what do they do? Front Cell Infect Microbiol. 2012;2:1–11.
Campbell JH, Foster CM, Vishnivetskaya T, Campbell AG, Yang ZK, Wymore A, et al. Host genetic and environmental effects on mouse intestinal microbiota. ISME J. 2012;6:2033–44.
Ganu RS, Harris RA, Collins K, Aagaard KM. Early Origins of Adult Disease: Approaches for Investigating the Programmable Epigenome in Humans, Nonhuman Primates, and Rodents. ILAR J. 2012;53:306–21.
Schiffrin EJ, Morley JE, Donnet-Hughes A, Guigoz Y. The inflammatory status of the elderly: the intestinal contribution. Mutat Res. 2010;690:50–6.
Furusawa Y, Obata Y, Fukuda S, Endo TA, Nakato G, Takahashi D, et al. Commensal microbe-derived butyrate induces the differentiation of colonic regulatory T cells. Nature. 2013;504:446–50.
Bermudez BM, Muñoz QS, Gomez LC, Matencio E, Bernal MJ, Romero F, et al. Cell-free culture supernatant of Bifidobacterium breve CNCM I-4035 decreases pro-inflammatory cytokines in human dendritic cells challenged with Salmonella typhi through TLR activation. PLoS One. 2013;8:e59370.
Groeger D, O’Mahony L, Murphy EF, Bourke JF, Dinan TG, Kiely B, et al. Bifidobacterium infantis 35624 modulates host inflammatory processes beyond the gut. Gut Microbes. 2013;4:325–39.
Sagar S, Vos AP, Morgan ME, Garssen J, Georgiou NA, Boon L, et al. The combination of Bifidobacterium breve with non-digestible oligosaccharides suppresses airway inflammation in a murine model for chronic asthma. Biochim Biophys Acta. 1842;2014:573–83.
Pande S, Shitut S, Freund L, Westermann M, Bertels F, Colesie C, et al. Metabolic cross-feeding via intercellular nanotubes among bacteria. Nat Commun. 2015;6:6238.
Wu BB, Yang Y, Xu X, Wang WP. Effects of Bifidobacterium supplementation on intestinal microbiota composition and the immune response in healthy infants. World J Pediatr. 2015. Epub ahead of print.
Pérez Martínez G, Bäuerl C, Collado MC. Understanding gut microbiota in elderly’s health will enable intervention through probiotics. Benef Microbes. 2014;5:235–46.
Mahowald MA, Rey FE, Seedorf H, Turnbaugh PJ, Fulton RS, Wollam A, et al. Characterizing a model human gut microbiota composed of members of its two dominant bacterial phyla. Proc Natl Acad Sci U S A. 2009;106:5859–64.
Falony G, Vlachou A, Verbrugghe K, De Vuyst L. Cross-feeding between Bifidobacterium longum BB536 and acetate-converting, butyrate-producing colon bacteria during growth on oligofructose. Appl Environ Microbiol. 2006;72:7835–41.
Sugahara H, Odamaki T, Fukuda S, Kato T, Xiao J, Abe F, et al. Probiotic Bifidobacterium longum alters gut luminal metabolism through modification of the gut microbial community. Sci Rep. 2015;5:13548.
Cantarel BL, Lombard V, Henrissat B. Complex carbohydrate utilization by the healthy human microbiome. PLoS One. 2012;7:e28742.
Kang DW, Park JG, Ilhan ZE, Wallstrom G, Labaer J, Adams JB, et al. Reduced incidence of prevotella and other fermenters in intestinal microflora of autistic children. PloS one. 2013;8(7):e68322.
Brown CT, Davis-Richardson AG, Giongo A, Gano KA, Crabb DB, Mukherjee N, et al. Gut microbiome metagenomics analysis suggests a functional model for the development of autoimmunity for type 1 diabetes. PLoS One. 2011;6:e25792.
Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, et al. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013;500:541–6.
Gill S, Pop M, DeBoy R, Eckburg P. Metagenomic analysis of the human distal gut microbiome. Science. 2006;312:1355–9.
Claesson MJ, Jeffery IB, Conde S, Power SE, O’Connor EM, Cusack S, et al. Gut microbiota composition correlates with diet and health in the elderly. Nature. 2012;488(7410):178–84.
Wang F, Yu T, Huang G, Cai D, Liang X, et al. Gut microbiota community and its assembly associated with age and diet in Chinese centenarians. J Microbiol Biotechnol. 2015;25:1195–204.
Biagi E, Candela M, Franceschi C, Brigidi P. The aging gut microbiota: new perspectives. Ageing Res Rev. 2011;10:428–9.
Camilleri M, Lee J, Viramontes B, Bharucha A, Tangalos E. Insights into the pathophysiology and mechanisms of constipation, irritable bowel syndrome, and diverticulosis in older people. J Am Geriatr Soc. 2000;48:1142–50.
Magrone T, Jirillo E. The interaction between gut microbiota and age-related changes in immune function and inflammation. Immun ageing. 2013;10:31.
Sato S, Kiyono H, Fujihashi K. Mucosal Immunosenescence in the Gastrointestinal Tract: A Mini-Review. Gerontology. 2014;61(4):336–42.
Clark RI, Salazar A, Yamada R, Fitz-Gibbon S, Morselli M, Alcaraz J, et al. Distinct Shifts in Microbiota Composition during Drosophila Aging Impair Intestinal Function and Drive Mortality. Cell Rep. 2015;12:1656–67.
Odamaki T, Xiao J-Z, Iwabuchi N, Sakamoto M, Takahashi N, Kondo S, et al. Influence of Bifidobacterium longum BB536 intake on faecal microbiota in individuals with Japanese cedar pollinosis during the pollen season. J Med Microbiol. 2007;56(Pt 10):1301–8.
Aronesty E. Comparison of sequencing utility programs. Open Bioinforma J. 2013;7:1–8.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Kuczynski J, Stombaugh J, Walters WA, González A, Caporaso JG, Knight R. Using QIIME to analyze 16S rRNA gene sequences from microbial communities. Curr Protoc Microbiol. 2012;Chapter 1:Unit 1E.5.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.
McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2011;6:610–8.
Culhane AC, Thioulouse J, Perrière G, Higgins DG. MADE4: an R package for multivariate analysis of gene expression data. Bioinformatics. 2005;21:2789–90.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(Karp 2001):2498–504.
Odamaki T, Xiao JZ, Sakamoto M, Kondo S, Yaeshima T, Iwatsuki K, et al. Distribution of different species of the Bacteroides fragilis group in individuals with Japanese cedar pollinosis. Appl Environ Microbiol. 2008;74:6814–7.
Matsuki T, Watanabe K, Fujimoto J, Kado Y, Takada T, Matsumoto K, et al. Quantitative PCR with 16S rRNA-gene-targeted species-specific primers for analysis of human intestinal bifidobacteria. Appl Environ Microbiol. 2004;70:167–73.
Matsuki T, Watanabe K, Fujimoto J, Takada T, Tanaka R. Use of 16S rRNA gene-targeted group-specific primers for real-time PCR analysis of predominant bacteria in human feces. Appl Environ Microbiol. 2004;70:7220–8.
Matsuda K, Tsuji H, Asahara T, Kado Y, Nomoto K. Sensitive quantitative detection of commensal bacteria by rRNA-targeted reverse transcription-PCR. Appl Environ Microbiol. 2007;73:32–9.
Langille MGI, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814–21.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.
Anderson MJ. A new method for non-parametric multivariate analysis of variance. Aust Ecol. 2001;26:32–46.
Availability of data and materials
DNA sequences corresponding to the 16S rRNA gene data have been deposited in DDBJ under accession number DRA004160. The data sets supporting the results of this study are included within the article and its Additional files.
TO is the corresponding author; he contributed to study conception and design, conducted most of the experiments, analyzed and interpreted the data, and wrote the manuscript. KK performed most of the experiments and data interpretation. HS conducted some of the laboratory analyses and data interpretation. NH and SH performed some of the laboratory analyses. JX contributed to study conception and design and reviewed the manuscript. FA contributed to study conception and design and data interpretation. RO designed and supervised the entire project and data interpretation as well as corrected the manuscript. All the authors read and approved the final manuscript.
All authors except for Ro Osawa are employees of Morinaga Milk Industry Co., Ltd., which has several probiotic products marketed worldwide.
Ethics approval and consent to participate
Two study protocols were for the collection of feces from subjects aged 21–65 years old or from community-dwelling elderly individuals and were approved by the Local Ethics Committee of the nonprofit organization Japan Health Promotion Supporting Network (Wakayama, Japan). The third study protocol was for the collection of feces from subjects aged 0–100 years old and was approved by the ethics committee of Kensyou-kai Incorporated Medical Institution (Osaka, Japan). Written informed consent was obtained from all subjects or from their legal guardians or relatives.
UniFrac clustering for gender distribution. (A) Unweighted and (B) weighted UniFrac PCoA of gut microbiota from 371 samples collected from the infant to the centenarian stage. Male and female subjects are displayed as blue and red, respectively. (PDF 155 kb)
Relative abundance of four predominant phyla. The box-plot indicates the interquartile range (IQR) of the relative abundance of each phylum at each age stage. (PDF 211 kb)
Relative abundance of each microbiota at genus level. (XLSX 873 kb)
Definition of bacterial co-abundance groups (CAGs). CAGs were defined by a heat plot showing Kendall correlations between genera clustered by Pearson’s correlation coefficient and Ward’s linkage hierarchical clustering. The colors within the clustering indicate the nine CAGs as shown in Fig. 3. (PDF 148 kb)
Network plot highlighting relationships between genera in nine CAGs. The colors of each node indicate the nine CAGs as shown in Fig. 3. Circle size indicates genus abundance. Pink and blue lines show significant positive and negative correlations between two bacterial genera with an absolute coefficient value greater than 0.3. (PDF 401 kb)
Bacterial groups at genus-like level in each CAG. (XLSX 863 kb)
Cell number of bacteria in vitro assay. (XLSX 855 kb)
Hierarchical Ward’s linkage clustering based on bacterial proportion at the genus level. Age-related groups (Infant, Elderly1, Elderly2, Adult1 and Adult2) were revealed through Ward’s linkage clustering using the squared Euclidean distance. The population densities (z-score) of genera scaled by color are displayed together with a dendrogram of bacterial genera in a heat map. The colors within the horizontal and vertical clustering represent age-segmented groups as shown in Fig. 1 and phylum as shown in Fig. 2. (PDF 245 kb)
LEfSe results on human microbiota in infant, adult and elderly cluster. (A) Histogram of the LDA scores computed for features differentially abundant between three age-related clusters. (B) Taxonomic representation of statistically differences between three age-related clusters. Differences are represented in the color of the most abundant class (yellow non-significant). Each circle’s diameter is proportional to the taxon’s abundance. (PDF 400 kb)
Taxa that are found in more than 50 % of the subjects in any cluster (shown in Additional file 8) with significantly difference between adult 1 and adult 2 clusters. (XLSX 857 kb)
Taxa that are found in more than 50 % of the subjects in any cluster (shown in Additional file 8) with significantly difference between elderly 1 and elderly 2 clusters. (XLSX 859 kb)
Individual data. (XLSX 875 kb)
Hierarchical Ward’s linkage clustering based on the proportion of transporter genes predicted by PICRUSt. Age-related groups (adult-enriched and infant/elderly-enriched clusters) were revealed by Ward’s linkage clustering using the squared Euclidean distance. The population densities (z-score) of the transporters scaled by color are displayed together with a dendrogram of the transporters in a heat map. The colors within the horizontal clustering represent the age-segmented groups as shown in Fig. 1 The color code for the vertical clustering indicates KEGG Orthology (KO) as follows: white, ABC Transporters, Eukaryotic Type; yellow, ABC Transporters, Prokaryotic Type; blue, Solute Carrier Family (SLC); orange, Major Facilitator Superfamily (MFS); red, Phosphotransferase System (PTS); and green, Other Transporters. (PDF 209 kb)
Relative abundance of predicted D-Xylose transporter (KEGG module: M00215). The KEGG module M00215 consists of three KO entries, K10543, K10544 and K10545. Each number indicates a group as shown in Table 1. Box-plots show the interquartile range (IQR) of the relative abundance of the predicted D-Xylose transporter. Open circles indicate outliers from 1.5- to 3.0-fold IQR. (PDF 96 kb)
Bacterial strains used in this study. (XLSX 855 kb)
PCR primers for detection of each species. (XLSX 857 kb)