Skip to main content

Aging progression of human gut microbiota



Human gut microbiota are important for human health and have been regarded as a “forgotten organ”, whose variation is closely linked with various factors, such as host genetics, diet, pathological conditions and external environment. The diversity of human gut microbiota has been correlated with aging, which was characterized by different abundance of bacteria in various age groups. In the literature, most of the previous studies of age-related gut microbiota changes focused on individual species in the gut community with supervised methods. Here, we aimed to examine the underlying aging progression of the human gut microbial community from an unsupervised perspective.


We obtained raw 16S rRNA sequencing data of subjects ranging from newborns to centenarians from a previous study, and summarized the data into a relative abundance matrix of genera in all the samples. Without using the age information of samples, we applied an unsupervised algorithm to recapitulate the underlying aging progression of microbial community from hosts in different age groups and identify genera associated to this progression. Literature review of these identified genera indicated that for individuals with advanced ages, some beneficial genera are lost while some genera related with inflammation and cancer increase.


The multivariate unsupervised analysis here revealed the existence of a continuous aging progression of human gut microbiota along with the host aging process. The identified genera associated to this aging process are meaningful for designing probiotics to maintain the gut microbiota to resemble a young age, which hopefully will lead to positive impact on human health, especially for individuals in advanced age groups.


The human gut, as an eco-system embodying more than 100 trillion microbes, plays an important role in human health [1]. The structure and composition of the gut flora are the result of long-term natural selection acting on both the microbes and host, which finally promotes mutual cooperation and functional stability of this complex ecosystem [2]. Factors such as diet, environment, host genetics and pathological conditions are important factors for explaining the variation of gut microbial community in different individuals [37]. Aging process captures many facets of biological variation of the human body, which leads to functional decline and increased incidence of infection in gut of elderly people [8]. Age-related changes of human gut microbiota have been revealed by several studies [917]. Hopkins et al. found higher numbers of Enterobacteria in children’s fecal than adults through culturing microbes [9]. Using 16S rRNA sequencing, Yatsunenko et al. found Bifidobacterium declined with increasing ages [11]. Odamaki et al. revealed that aging was accompanied by increasing proportion of Bacteroides, Eubacterium and Clostridiaceae; Enterobacteriaceae were enriched in infant and elderly; Bifidobacterium were enriched in infants; Lachnospiraceae were enriched in adults [10]. Using whole genome sequencing, Stewart et al. discovered decline of L-lactate dehydrogenase (milk fermentation) and increase of transketolase (metabolism of fiber) over the first year of life [13]. In these studies, various supervised machine learning methods have been applied, including multi-group comparative analysis with permutational analysis of variance (PERMANOVA) [9, 10, 12, 17], Spearman rank correlation and Random Forest [11], as well as frequency-inverse document frequency and minimum-redundancy maximum-relevance [14], which effectively identified taxonomic or functional signatures showed aging-related changes of gut microbiota.

In this study, we proposed to explore an unsupervised machine learning approach for identifying aging-related progression of microbiota community and bacteria genera associated with the progression. The unsupervised algorithm adopted here is called Sample Progression Discovery (SPD), which was developed to identify progressive changing patterns of gene expression that reflect the biological progression in various biological processes and systems [18]. This idea was first applied to microarray gene expression analysis [18], and then extended to flow cytometry [19] and single-cell RNA-seq analysis [20]. Here, we applied SPD on community profiles extracted from 16S rRNA sequencing data of human gut microbiota samples in various age periods ranging from new-born babies to centenarians. SPD recapitulated the underlying aging progression in the data in an unsupervised fashion, and sorted the gut microbiota samples in an order consistent to the host ages. In addition, SPD identified bacteria genera associated with the aging-related progression of gut microbiota. These findings demonstrated the existence of an aging progression of human gut microbial community, and points to important bacteria genera that characterize the aging of gut microbiota.


Data annotation and samples overview

We obtained a total of 3.2 million high-quality 16S rRNA sequences from 368 samples [10], with 8734±2748 (mean ± deviation) reads per sample. The 16S rRNA sequences were binned into 366 genera using the Mothur pipeline [21] with SILVA [22] as reference database (see Methods). We removed 119 genera with extremely low abundance, the total amount of sequences annotated as these genera only accounted for 0.01% of all the sequences. Also, we excluded one sample with abnormally high proportion of Pseudomonas, which is an indication of abnormal sampling or pathological disorder of this individual ‘Japanese 320’. Overall, we derived a relative abundance matrix of the 247 genera across the 367 samples, which served as the basis for further analyses. To reveal age-related progression of gut microbiota, we divided the samples into 14 age groups considering body transformation periods. New-born babies were grouped according to their weaning status and adults were grouped by decade (Table 1). Except the centenarians, there were at least 10 samples in each age group.

Table 1 Samples were grouped into 14 age-segment groups

We performed PCA to visualize the taxonomic patterns of these samples in a low-dimension space based on the relative abundance matrix of the 247 genera across the 367 samples. The top three principle components explained 33.17%, 15.09% and 10.32% of the original data variance, respectively. As shown in Fig. 1, the samples from children younger than three years old scattered loosely, and were quite different from each other. This observation was consistent with previous literature [11], which concluded that intrapersonal variation decreased as a function of age. Nevertheless, the samples did not form distinct groups when visualized by this linear approach.

Fig. 1

Sample overview using PCA. Using the relative abundance of 247 genera across all the 367 samples as input, we linearly transformed and visualized the data in a three-dimensional space. Each sample is represented by one dot, colored according to age. Samples from children younger than three (the dark blue dots) scattered most distantly, while older age groups were mixed together in the PCA space

Age-related variation of gut microbiota revealed by supervised methods

We applied two previous statistical approaches to identify the age-related variation of the gut microbiota in a univariate fashion. First, we applied permutational one-way ANOVA test [23] to the genus relative abundance matrix to identify genera that significantly varied in different age groups. The abundances of 43 genera showed significant difference across the age groups with P<0.001 (1000 randomizations), and the P values were adjusted using Bonferroni correction (see more details in Additional file 1). We also applied Spearman correlation to find genera that co-vary with age. There were 17 genera positively correlated with aging and one genus negatively correlated with aging (Additional file 2). These results were consistent with multiple previous literatures, showing that individual genus in the gut microbial community varied during the host aging [9, 10, 12, 17]. Further question naturally arose as to whether the gut microbial community as a whole shift continuously during aging.

Aging progression of gut microbiota revealed by unsupervised analysis

Different from the previous supervised univariate methods searching for features that co-varied with aging, we applied an unsupervised method SPD to examine the gut microbiota data in a multivariate fashion. The input to SPD was the averages of genus relative abundance of samples in each age group, which is a 247×14 matrix. The relative abundance of each feature was normalized across samples to release the scale effect. Based on each of the genus features, a minimum spanning tree (MST) was constructed according to Euclidean distance, which represented a putative progression ordering among the 14 sample groups. The 247 genera and the 247 resulting MSTs were cross compared to examine whether multiple genera fitted well the same progression ordering among the samples. Results of these comparisons were summarized into a progression similarity matrix, where each element counted the number of progression orderings that two genera both fit well with. As shown in Fig. 2a and magnified in Fig. 2b, the progression similarity matrix revealed a subset of 35 genera (Additional file 3) that fitted well with a common set of putative progression orderings. Using this subset of genera, an overall minimal spanning tree was constructed to represent the common progression ordering, shown in Fig. 2c. Each node of the tree represented one age group. Nodes were labeled and colored according to their age groups to assist the visualization. However, the age information was not used to determine the structure of the tree. This overall minimal spanning tree is what SPD aimed to identify, a progression ordering among the samples, with respect to which multiple features exhibited gradual changes. The overall minimal spanning in Fig. 2c recapitulated the age progression ordering across the 14 sample groups. Especially, when we further classify these sample groups into four larger groups, i.e. Children and teenagers, Adults, Elderly and Centenarians, the order of sample groups on this minimal spanning tree perfectly matched with the ages of sample groups. This is an interesting result, because SPD was able to recover the correct ordering of aging progression based on the genus relative abundance alone, which implied that there existed an aging progression of the human gut microbiota.

Fig. 2

SPD recovered aging progression with taxonomical composition of human gut microbiota. a Progression similarity matrix for all genera, with each element counting the number of progression orderings the two corresponding genera shared. b We manually picked the highlighted area from (a). These selected genera were consistent with a common set of putative progression orderings. c An overall minimal spanning tree of the 14 age groups based on the selected genera. Each node represents one age group

critical genera underlying the aging progression of gut microbiota

We further examined the 35 selected genera that contributed to the aging progression of gut microbiota, and compared to previous methods. 11 of the 35 genera were significant in the permutational one-way ANOVA analysis with adjusted P<0.001. Among the remaining 24 genera detected by SPD only, a few were previously implicated in the literature, such as Oxalobacter, Butyrivibrio, Lactobacillus which have been experimentally demonstrated to be associated with aging [2426], as well as Prevotellaceae which has been highlighted with lower presence in the gut microbiota of centenarians [27]. Among the 35 genera selected according to the progression similarity defined by SPD, only 9 exhibited monotonic changes with respect to aging, while the rest first increased and then descreased in different age periods (Additional file 4: figure 1). This was because SPD was designed to identify features that exhibited gradual changes with respect to a common underlying progression pattern, and the gradual changes were not limited to be monotonic. Therefore, this analysis was able to identify genera that gradually changed without abrupt fluctuations during aging. We performed extensive literature review of these 35 genera, and found a lot of previous reports of the functional relevance of these genera.

Genera shown in Fig. 3 shared one common feature. Their abundances increased with respect to aging, but decreased in the extremely elderly subjects. Among these genera, Lactobacillus species are commonly used as probiotics [28]. Oscillospira species have been frequently reported as enriched in lean subjects compared to the obese subjects [2932], and are central to the human gut microbiota for degrading fibers [33]. Oxalobacter is responsible for degrading oxalate in the gut. It has been experimentally demonstrated appearing in the gut of almost all young individuals, but these bacterium may later be lost during aging [24]. Prevotellaceae is commonly found in the gastric system of people who maintain a diet low in animal fats and high in carbohydrates [34] and is lost in centenarians [27]. Researchers also found that there was an increased abundance of Prevotellaceae in the guts of healthy people compared with people with Parkinson’s disease [35]. Parascardovia is a genus of Bifidobacteriaceae, which has been shown to provide health-promoting benefits to the host [36]. Butyrivibrio species have been experimentally proved as butyrate producing bacteria, and butyrate is a preferred energy source for colonic epithelial cells and is thought to play an important role in maintaining colonic health in humans [37]. Overall, the decrease of these beneficial genera in the elderly age groups, especially centenarians, maybe manifestation of or causial associations to decline of health in those age groups.

Fig. 3

Genera that first increased and then decreased during aging, especially sharply decreased in the 13th or 14th age groups, or both

In contrast, genera in Fig. 4 showed generally monotonically increasing patterns with respect to aging. Parvimonas has been reported to be enriched in colorectal cancer [3842]. Anaerotruncus was relatively enriched in patients with age-related macular degeneration [43]. Corynebacterium was reported as more abundant in the gut of autistic individuals (autism spectrum disorders) [44]. Many Corynebacterium species were also reported as involved in human and animal diseases [45]. GCA-900066225 is one genus in the Lachnospiraceae family, which has been reported to be associated with ulcerative colitis, Crohn’s and celiac disease, as well as the stress of the host [46]. Desulfovibrio species produce hydrogen sulfide using sulfate as the electron acceptor, and these sulfate-reducing bacteria are positively associated with inflammation [47, 48]. A human stool-derived Bilophila wadsworthia strain caused systemic inflammation in specific-pathogen-free mice [49]. Tumor-bearing mice showed enrichment in species of Odoribacter [50]. Butyricimonas was enriched in the subjects suffering from high rectal temperature, systolic blood pressure, and heart rate and a significantly lower physical activity score [51]. Overall, these monotonically increasing genera were often linked to inflammation and diseases.

Fig. 4

Genera that exhibited general increasing patterns during aging

All these prior literature of the identified genera pointed to one interesting observation. Many of the genera that first increased and then descreased were previously implicated as health beneficial, whereas most of the monotonically increasing genera were frequently reported as disease-related. When individuals turn elderly beyond 90s, their guts tend to lose some of the beneficial genera and gain potentially harmful genera.


Since the variation of gut microbiota is closely linked to the health status of the host body, an ideal dataset for studying aging of human gut microbiota should be collected from healthy subjects of various age groups. Unfortunately, the health status of individuals included in this study is unknown, because that the data were obtained from a published paper [10] which did not provide the health information corresponding to the samples. During our literature search on age-related changes of the human gut microbiota, we found in multiple previously published papers [10, 11, 14, 52, 53] that the health status of individuals in such studies are often not reported.

In order to gauge the health status of samples in the dataset used in this paper, we performed additional analysis by referencing to multiple previous datasets on the human gut microbiome of hosts suffering from different diseases [4, 5, 7, 5457]. In each of the previous datasets, we obtained the relative abundance of the human gut microbial genera, and visualized their distributions for both healthy controls and disease samples.. Multiple genera were significantly enriched in the diseases compared to the healthy controls. Most of these genera have been reported as opportunistic pathogens of the human gut [5867]. These disease-enriched genera typically showed higher abundance and higher variance in disease subjects compared to the healthy ones (first and second columns of Additional file 4: figure 2), while all of these genera exhibit low abundance in the dataset used here (third column of Additional file 4: figure 2). This observation indicated that the samples in the current dataset are more similar to the healthy samples in the previous datasets, and dissimilar to the diseases samples. This comparison demonstrated that majority of the samples in this dataset were derived from healthy subjects.

For 16S rRNA sequencing data analysis, OTU (operational taxonomic unit) is another commonly used classification unit, which allows for binning sequences into features at finer resolution compared to the genus level features. To confirm the observations in our genus level analysis, we applied the progression analysis to the OTU level features. 4663 OTUs were defined by clustering sequences with similarity threshold set as 0.97 for species level classification. After filtering out OTUs with extremely low abundances, the averages of the relative abundances of the remaining 1229 OTUs were calculated for each age group. Progression analysis based on OTU features was able to partially recapitulate the correct order of the age groups (Additional file 4: figure 3), but slightly worse compared to the result from genus level analysis shown in Fig. 2c. It is reassuring that the progression analysis at both OTU level and genus level consistently revealed aging related progression of the human gut microbiota.

In the metagenomics literature, the alpha diversity and the beta diversity are popular metrics for providing quantitative summaries of species diversity. We computed the alpha diversity and the beta diversity based on the averages of genus relative abundance of samples in each age group. The alpha diversity was quantified by the Shannon index and the beta diversity was quantified by Bray-Curtis dissimilarity between different age groups. Additional file 4: figure 4 shows the alpha diversity computed for each individual age group, which showed a steady increase of the alpha diversity as a function of aging, except for the steep drop in the extremely elderly age group [99, 110]. This is consistent with the results shown in Fig. 3, where multiple aging-related genera showed significant decrease in the extremely elderly age group. The beta diversity quantified the dissimilarity between different age groups (Additional file 4: figure 5). Focusing on the beta diversity between neighboring age groups, we observed that the dissimilarity between groups [2, 3] and between groups [13, 14] were notably larger than the dissimilarity between other neighboring age groups. The distinction between group 2 (weaning) and group 3 (weaned) is mainly due to the transformation of weaning status, which is accompanied by drastic dietary changes. However, samples of group 13 and group 14 are all elderly individuals with continuous ages, and the large dissimilarity between groups 13 and 14 cannot be explained by changes of dietary habits. Therefore, we conjecture that the large dissimilarity between groups 13 and 14 is due to the aging of gut microbiota, manifested in the sudden decrease of multiple genera in the extremely elderly samples. Overall, according to both alpha and beta diversity, we can see indications consistent to our observation of the sudden decrease of multiple genera in extremely elderly age samples shown in Fig. 3.


We applied an unsupervised machine learning approach SPD on genera abundance profile of human gut microbiota quantified by 16S rRNA sequencing data. Without using the age information of the samples, SPD sorted sample groups on a minimal spanning tree that recapitulated the aging progression. This result indicated the existence of an aging progression reflected in the human gut microbiota. In the meantime, we found 35 genera associated with this age-related progression. Some of these genera were not identified using the commonly-used statistical approaches for metagenomics analysis. Literature review of these 35 genera led to a lot of evidences of the functional relevance of these genera. The evidences collectively indicated an age-related decline of the beneficial functions of gut microbiota, as well as increase of inflammation and diseases, especially for the elderly people older than 90s.


Data and data annotation

Our study includes 371 samples of subjects ranging from new-born babies to centenarians, which have been described in publication [10]. We downloaded the 16S rRNA data from DNA data bank of Japan with accession number DRA004160. Three samples were discarded because of only one end of paired-end reads were released. We performed 16S rRNA data processing using Mothur [21]. Low quality reads with average quality score <25 or read length <150bp were filtered out. We set the minimum length of reads as 150bp because the overlap region of each pair of reads was about 150bp. The number of reads in each sample was Gaussian distributed (8734±2748), which implied that all the 368 samples were sequenced in normal depth. The high-quality reads with both paired ends were merged as sequences. Those low-quality reads or reads with only one end were discarded. Then we aligned the sequences against Silva reference database version 132 [22] to infer taxonomical composition of samples. Threshold for the alignment was set as bootstrap confidence value 80% (80% identity) during 100 iterations. Based on the alignment result, we revealed the taxonomic composition at genus level. There were 368 genera in all the samples.

Feature matrix

We defined the genus abundance matrix N = {n ij}, where n ij is the number of reads of sample i binned into genus j. One hundred nineteen genera were filtered out for their extremely low abundance, and three genera were combined into one genus cluster as “unclassified”, after which 247 features were obtained for further analysis. To normalize out the variation in sequencing depth of different samples, the genus abundance matrix was transferred into a relative abundance matrix F = {f ij}, where f ij=\(n_{ij}/\sum _{k=1}^{247}n_{i,k}\). One sample from subject “Japanese 320” was discarded for its abnormally high proportion of Pseudomonas. Finally, we have a 367×247 relative abundance matrix F for further analysis.

With decent numbers of observations in different age periods, we estimated the genus relative abundance of population in each age group by calculating the mean value of samples in corresponding group, which partially reduced the variations across individual samples and sparsity of the data matrix. Age segments were defined concerning the physiological transition of the host bodies, wherein the new-born babies were grouped according to their weaning status and the adults were grouped by decade. The number of samples in each age group was depicted in Table 1.

Availability of data and materials

The data we used in this paper was downloaded from previously published paper [10]. The data we generated during data analysis was released as additional files.



Minimum spanning tree


Operational taxonomic unit


Permutational analysis of variance


Sample progression discovery


  1. 1

    Bocci V. The neglected organ: bacterial flora has a crucial immunostimulatory role. Perspect Biol Med. 1992; 35(2):251–60.

    Article  CAS  Google Scholar 

  2. 2

    O’Hara AM, Shanahan F. The gut flora as a forgotten organ. EMBO Rep. 2006; 7(7):688–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. 3

    Human Microbiome Project C.Structure, function and diversity of the healthy human microbiome. Nature. 2012; 486(7402):207–14.

    Article  CAS  Google Scholar 

  4. 4

    Karlsson FH, Fåk F, Nookaew I, Tremaroli V, Fagerberg B, Petranovic D, Bäckhed F, Nielsen J. Symptomatic atherosclerosis is associated with an altered gut metagenome. Nature Commun. 2012; 3:1245.

    Article  CAS  Google Scholar 

  5. 5

    Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, Almeida M, Arumugam M, Batto JM, Kennedy S, Leonard P, Li J, Burgdorf K, Grarup N, Jorgensen T, Brandslund I, Nielsen HB, Juncker AS, Bertalan M, Levenez F, Pons N, Rasmussen S, Sunagawa S, Tap J, Tims S, Zoetendal EG, Brunak S, Clement K, Dore J, Kleerebezem M, Kristiansen K, Renault P, Sicheritz-Ponten T, de Vos WM, Zucker JD, Raes J, Hansen T, Meta HITc, Bork P, Wang J, Ehrlich SD, Pedersen O. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013; 500(7464):541–6.

    Article  CAS  Google Scholar 

  6. 6

    Griffin NW, Ahern PP, Cheng J, Heath AC, Ilkayeva O, Newgard CB, Fontana L, Gordon JI. Prior dietary practices and connections to a human gut microbial metacommunity alter responses to diet interventions. Cell Host Microbe. 2017; 21(1):84–96.

    Article  CAS  Google Scholar 

  7. 7

    Lewis JD, Chen EZ, Baldassano RN, Otley AR, Griffiths AM, Lee D, Bittinger K, Bailey A, Friedman ES, Hoffmann C, Albenberg L, Sinha R, Compher C, Gilroy E, Nessel L, Grant A, Chehoud C, Li H, Wu GD, Bushman FD. Inflammation, antibiotics, and diet as environmental stressors of the gut microbiome in pediatric crohn’s disease. Cell Host Microbe. 2015; 18(4):489–500.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. 8

    B LL. Age related changes in gut physiology and nutritional status. Gut. 1996.

  9. 9

    Hopkins MJ, Sharp R, Macfarlane GT. Age and disease related changes in intestinal bacterial populations assessed by cell culture, 16s rrna abundance, and community cellular fatty acid profiles. Gut. 2001.

  10. 10

    Odamaki T, Kato K, Sugahara H, Hashikura N, Takahashi S, Xiao JZ, Abe F, Osawa R. Age-related changes in gut microbiota composition from newborn to centenarian: a cross-sectional study. BMC Microbiol. 2016; 16:90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. 11

    Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP, Heath AC, Warner B, Reeder J, Kuczynski J, Caporaso JG, Lozupone CA, Lauber C, Clemente JC, Knights D, Knight R, Gordon JI. Human gut microbiome viewed across age and geography. Nature. 2012; 486(7402):222–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 12

    Bian G. The gut microbiota of healthy aged chinese is similar to that of the healthy young. mSphere. 2017.

  13. 13

    Stewart CJ, Ajami NJ, O’Brien JL, Hutchinson DS, Smith DP, Wong MC, Ross MC, Lloyd RE, Doddapaneni H, Metcalf GA, Muzny D, Gibbs RA, Vatanen T, Huttenhower C, Xavier RJ, Rewers M, Hagopian W, Toppari J, Ziegler AG, She JX, Akolkar B, Lernmark A, Hyoty H, Vehik K, Krischer JP, Petrosino JF. Temporal development of the gut microbiome in early childhood from the teddy study. Nature. 2018; 562(7728):583–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. 14

    Lan Y, Kriete A, Rosen GL. Selecting age-related functional characteristics in the human gut microbiome. Microbiome. 2013; 1(1):2.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 15

    Buford TW. (dis)trust your gut: the gut microbiome in age-related inflammation, health, and disease. Microbiome. 2017; 5(1):80.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16

    Minelli EB, Benini A, Beghini AM, Cerutti R, Nardo G. Bacterial faecal flora in healthy women of different ages. Microb Ecol Health Dis. 1993; 6(2):43–51.

    Article  Google Scholar 

  17. 17

    Vatanen T, Franzosa EA, Schwager R, Tripathi S, Arthur TD, Vehik K, Lernmark A, Hagopian WA, Rewers MJ, She JX, Toppari J, Ziegler AG, Akolkar B, Krischer JP, Stewart CJ, Ajami NJ, Petrosino JF, Gevers D, Lahdesmaki H, Vlamakis H, Huttenhower C, Xavier RJ. The human gut microbiome in early-onset type 1 diabetes from the teddy study. Nature. 2018; 562(7728):589–594.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. 18

    Peng Q, Gentles AJ, Plevritis SK. Discovering biological progression underlying microarray samples. PLoS Comput Biol. 2011.

  19. 19

    Peng Q, Simonds EF, Plevritis SK. Extracting a cellular hierarchy from high-dimensional cytometry data with spade. Nat Biotechnol. 2011; 29(10):886–91.

    Article  CAS  Google Scholar 

  20. 20

    Trapnell C, Cacchiarelli D, Grimsby J, Pokharel P, Li S, Morse M, Lennon NJ, Livak KJ, Mikkelsen TS, Rinn JL. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol. 2014; 32(4):381.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. 21

    Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. J Appl Environ Microbiol. 2009; 75(23):7537–41.

    Article  CAS  Google Scholar 

  22. 22

    Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO. The silva ribosomal rna gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012; 41(D1):590–6.

    Article  CAS  Google Scholar 

  23. 23

    Anderson MJ. Permutation tests for univariate or multivariate analysis of variance and regression. Can J Fish Aquat Sci. 2001; 58(3):626–39.

    Article  Google Scholar 

  24. 24

    Duncan SH, Richardson AJ, Kaul P, Holmes RP, Allison MJ, Stewart CS. Oxalobacter formigenes and its potential role in human health. J Appl Environ Microbiol. 2002; 68(8):3841–7.

    Article  CAS  Google Scholar 

  25. 25

    Sovran B, Hugenholtz F, Elderman M, Van Beek AA, Graversen K, Huijskes M, Boekschoten MV, Savelkoul HF, De Vos P, Dekker J. Age-associated impairment of the mucus barrier function is associated with profound changes in microbiota and immunity. Sci Rep. 2019; 9(1):1437.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. 26

    Vemuri R, Shinde T, Gundamaraju R, Gondalia S, Karpe A, Beale D, Martoni C, Eri R. Lactobacillus acidophilus dds-1 modulates the gut microbiota and improves metabolic profiles in aging mice. Nutrients. 2018; 10(9):1255.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. 27

    Tuikhar N, Keisam S, Labala RK, Ramakrishnan P, Arunkumar MC, Ahmed G, Biagi E, Jeyaram K. Comparative analysis of the gut microbiota in centenarians and young adults shows a common signature across genotypically non-related populations. Mech Ageing Dev. 2019.

  28. 28

    Walter J. Ecological role of lactobacilli in the gastrointestinal tract: implications for fundamental and biomedical research. J Appl Environ Microbiol. 2008; 74(16):4985–96.

    Article  CAS  Google Scholar 

  29. 29

    Tims S, Derom C, Jonkers DM, Vlietinck R, Saris WH, Kleerebezem M, De Vos WM, Zoetendal EG. Microbiota conservation and bmi signatures in adult monozygotic twins. ISME J. 2013; 7(4):707.

    Article  CAS  Google Scholar 

  30. 30

    Goodrich JK, Waters JL, Poole AC, Sutter JL, Koren O, Blekhman R, Beaumont M, Van Treuren W, Knight R, Bell JT. Human genetics shape the gut microbiome. Cell. 2014; 159(4):789–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. 31

    Verdam FJ, Fuentes S, de Jonge C, Zoetendal EG, Erbil R, Greve JW, Buurman WA, de Vos WM, Rensen SS. Human intestinal microbiota composition is associated with local and systemic inflammation in obesity. Obesity (Silver Spring). 2013; 21(12):607–15.

    Article  CAS  Google Scholar 

  32. 32

    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:311.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33

    Kohl KD, Amaya J, Passement CA, Dearing MD, McCue MD. Unique and shared responses of the gut microbiota to prolonged fasting: a comparative study across five classes of vertebrate hosts. FEMS Immunol Med Microbiol. 2014; 90(3):883–94.

    CAS  Google Scholar 

  34. 34

    Xu Z, Knight R. Dietary effects on human gut microbiome diversity. Br J Nutr. 2015; 113(S1):1–5.

    Article  CAS  Google Scholar 

  35. 35

    Unger MM, Spiegel J, Dillmann K-U, Grundmann D, Philippeit H, Bürmann J, Faßbender K, Schwiertz A, Schäfer K-H. Short chain fatty acids and gut microbiota differ between patients with parkinson’s disease and age-matched controls. Parkinsonism Relat Disord. 2016; 32:66–72.

    Article  Google Scholar 

  36. 36

    O’Callaghan A, van Sinderen D. Bifidobacteria and their role as members of the human gut microbiota. Front Microbial. 2016; 7:925.

    Google Scholar 

  37. 37

    Barcenilla A, Pryde SE, Martin JC, Duncan SH, Stewart CS, Henderson C, Flint HJ. Phylogenetic relationships of butyrate-producing bacteria from the human gut. J Appl Environ Microbiol. 2000; 66(4):1654–61.

    Article  CAS  Google Scholar 

  38. 38

    Zou S, Fang L, Lee MH. Dysbiosis of gut microbiota in promoting the development of colorectal cancer. Gastroenterol Rep (Oxf). 2018; 6(1):1–12.

    Article  Google Scholar 

  39. 39

    Flemer B, Lynch DB, Brown JM, Jeffery IB, Ryan FJ, Claesson MJ, O’Riordain M, Shanahan F, O’Toole PW. Tumour-associated and non-tumour-associated microbiota in colorectal cancer. Gut. 2017; 66(4):633–43.

    Article  CAS  Google Scholar 

  40. 40

    Nakatsu G, Li X, Zhou H, Sheng J, Wong SH, Wu WK, Ng SC, Tsoi H, Dong Y, Zhang N, He Y, Kang Q, Cao L, Wang K, Zhang J, Liang Q, Yu J, Sung JJ. Gut mucosal microbiome across stages of colorectal carcinogenesis. Nat Commun. 2015; 6:8727.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. 41

    Duvallet C, Gibbons SM, Gurry T, Irizarry RA, Alm EJ. Meta-analysis of gut microbiome studies identifies disease-specific and shared responses. Nat Commun. 2017; 8(1):1784.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. 42

    Chen W, Liu F, Ling Z, Tong X, Xiang C. Human intestinal lumen and mucosa-associated microbiota in patients with colorectal cancer. PLoS One. 2012; 7(6):39743.

    Article  CAS  Google Scholar 

  43. 43

    Zinkernagel MS, Zysset-Burri DC, Keller I, Berger LE, Leichtle AB, Largiader CR, Fiedler GM, Wolf S. Association of the intestinal microbiome with the development of neovascular age-related macular degeneration. Sci Rep. 2017; 7:40826.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. 44

    Strati F, Cavalieri D, Albanese D, De Felice C, Donati C, Hayek J, Jousson O, Leoncini S, Renzi D, Calabro A, De Filippo C. New evidences on the altered gut microbiota in autism spectrum disorders. Microbiome. 2017; 5(1):24.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45

    Padmanabhan R, Dubourg G, Lagier JC, Couderc C, Michelle C, Raoult D, Fournier PE. Genome sequence and description of corynebacterium ihumii sp. nov. Stand Genomic Sci. 2014; 9(3):1128–43.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 46

    Li S, Wang Z, Yang Y, Yang S, Yao C, Liu K, Cui S, Zou Q, Sun H, Guo G. Lachnospiraceae shift in the microbial community of mice faecal sample effects on water immersion restraint stress. AMB Express. 2017; 7(1):82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. 47

    Devkota S, Wang Y, Musch MW, Leone V, Fehlner-Peach H, Nadimpalli A, Antonopoulos DA, Jabri B, Chang EB. Dietary-fat-induced taurocholic acid promotes pathobiont expansion and colitis in il10 −/− mice. Nature. 2012; 487(7405):104–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. 48

    Loubinoux J, Bronowicki J-P, Pereira IA, Mougenel J-L, Le Faou AE. Sulfate-reducing bacteria in human feces and their association with inflammatory bowel diseases. FEMS Immunol Med Microbiol. 2002; 40(2):107–12.

    CAS  Google Scholar 

  49. 49

    Feng Z, Long W, Hao B, Ding D, Ma X, Zhao L, Pang X. A human stool-derived bilophila wadsworthia strain caused systemic inflammation in specific-pathogen-free mice. Gut Pathog. 2017; 9(1):59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. 50

    Zackular JP, Baxter NT, Iverson KD, Sadler WD, Petrosino JF, Chen GY, Schloss PD. The gut microbiome modulates colon tumorigenesis. MBio. 2013; 4(6):00692–13.

    Article  CAS  Google Scholar 

  51. 51

    Shin J-H, Sim M, Lee J-Y, Shin D-M. Lifestyle and geographic insights into the distinct gut microbiota in elderly women from two different geographic locations. J Physiol Anthropol. 2016; 35(1):31.

    Article  PubMed  PubMed Central  Google Scholar 

  52. 52

    Maffei VJ, Kim S, Blanchard IV E, Luo M, Jazwinski SM, Taylor CM, Welsh DA. Biological aging and the human gut microbiota. J Gerontol A Biomed Sci Med Sci. 2017; 72(11):1474–82.

    Article  CAS  Google Scholar 

  53. 53

    de la Cuesta-Zuluaga J, Kelley S, Chen Y, Escobar J, Mueller N, Ley R, McDonald D, Huang S, Swafford A, Knight R, et al.Age-and sex-dependent patterns of gut microbial diversity in human adults. mSystems 4: e00261-19. 2019.

  54. 54

    Karlsson FH, Tremaroli V, Nookaew I, Bergström G, Behre CJ, Fagerberg B, Nielsen J, Bäckhed F. Gut metagenome in european women with normal, impaired and diabetic glucose control. Nature. 2013; 498(7452):99–103.

    Article  CAS  Google Scholar 

  55. 55

    Nielsen HB, Almeida M, Juncker AS, Rasmussen S, Li J, Sunagawa S, Plichta DR, Gautier L, Pedersen AG, Le Chatelier E, et al.Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nature Biotechnol. 2014; 32(8):822.

    Article  CAS  Google Scholar 

  56. 56

    Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, Liang S, Zhang W, Guan Y, Shen D, et al.A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012; 490(7418):55.

    Article  CAS  Google Scholar 

  57. 57

    Qin N, Yang F, Li A, Prifti E, Chen Y, Shao L, Guo J, Le Chatelier E, Yao J, Wu L, et al.Alterations of the human gut microbiome in liver cirrhosis. Nature. 2014; 513(7516):59.

    Article  CAS  Google Scholar 

  58. 58

    Lobatón T, Hoffman I, Vermeire S, Ferrante M, Verhaegen J, Van Assche G. Aeromonas species: an opportunistic enteropathogen in patients with inflammatory bowel diseases? a single center cohort study. Inflamm Bowel Dis. 2014; 21(1):71–8.

    Article  Google Scholar 

  59. 59

    Gophna U, Sommerfeld K, Gophna S, Doolittle WF, van Zanten SJV. Differences between tissue-associated intestinal microfloras of patients with crohn’s disease and ulcerative colitis. J Clin Microbiol. 2006; 44(11):4136–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. 60

    Dalal SR, Chang EB. The microbial basis of inflammatory bowel diseases. J Clin Inv. 2014; 124(10):4190–6.

    Article  CAS  Google Scholar 

  61. 61

    Padmanabhan R, Dubourg G, Lagier J-C, Couderc C, Michelle C, Raoult D, Fournier P-E. Genome sequence and description of corynebacterium ihumii sp. nov. Stand Genomic Sci. 2014; 9(3):1128.

    Article  PubMed  PubMed Central  Google Scholar 

  62. 62

    Feeney A, Kropp KA, O’Connor R, Sleator RD. Cronobacter sakazakii: stress survival and virulence potential in an opportunistic foodborne pathogen. Gut Microbes. 2014; 5(6):711–8.

    Article  Google Scholar 

  63. 63

    Zeng M, Inohara N, Nuñez G. Mechanisms of inflammation-driven bacterial dysbiosis in the gut. Mucosal Immunol. 2017; 10(1):18.

    Article  CAS  Google Scholar 

  64. 64

    Guinane CM, Cotter PD. Role of the gut microbiota in health and chronic gastrointestinal disease: understanding a hidden metabolic organ. Ther Adv Gastroenterol. 2013; 6(4):295–308.

    Article  Google Scholar 

  65. 65

    Liu H, Zhu J, Hu Q, Rao X. Morganella morganii, a non-negligent opportunistic pathogen. Int J Infect Dis. 2016; 50:10–17.

    Article  Google Scholar 

  66. 66

    Hamilton AL, Kamm MA, Ng SC, Morrison M. Proteus spp. as putative gastrointestinal pathogens. Clin Microbiol Rev. 2018; 31(3):00085–17.

    Article  Google Scholar 

  67. 67

    Anderson MJ. Permutation tests for univariate or multivariate analysis of variance and regression. Can J Fish Aquat Sci. 2001; 58(3):626–39.

    Article  Google Scholar 

Download references


We thank Dr. Xingyu Yang, Mr. Qinwei Zhuang, Mr. Hong Seo Bruce Lim and Dr. Cheng Zhu at Georgia Institute of Technology for their interest to the project and useful discussions. We also thank Dr. Kuang Chen and Mr. Zhongjie Xie at Peking University for their helpful technique supporting and useful discussions. Dr. Huaiqiu Zhu shared the corresponding authorship with Dr. Peng Qiu.


National Science Foundation (CCF1552784) and the China Scholarship Council supported the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




CX performed the analysis and prepared the manuscript. CX, PQ and HZ designed the project. PQ supervised the project and revised the manuscript. HZ reviewed the manuscript. All authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Peng Qiu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Xu, C., Zhu, H. & Qiu, P. Aging progression of human gut microbiota. BMC Microbiol 19, 236 (2019).

Download citation


  • Human gut microbiota
  • Aging
  • Sample progression discovery
  • 16S rRNA sequencing