Aging in SPF mice is associated with an overall alteration in the structure of the gut microbiota
To characterize aging-associated gut microbiota and metabolic phenotype changes, fecal and plasma samples are repetitively collected from GF and SPF mice at 3, 4, 6, 8, 12 weeks after birth, and every 4 weeks thereafter. We divided the experimental time course into 4 categories: 3 to 12 weeks (wean-young), 16 to 32 weeks (young), 36 to 52 weeks (early-middle) and 56 to 72 weeks (middle), respectively (Supplementary Fig. 1), based on the age-related murine frailty index[24] and a report showing that the murine gut microbiota changes significantly during weaning[20]. The increase of body weight in SPF mice during aging (Supplementary Fig. 2) is mostly due to an increase in fat mass, even though food intake is slightly decreased[21, 25, 26]. Weight gain in middle-age GF mice compared to middle-age SPF mice is probably due to the gradual cecum enlargement characteristic of GF mice because of the accumulation of hydrated dietary fiber components[32].
To assess the overall structures of the gut microbiota, we amplicon-sequenced the V1-V2 region of the 16S rRNA gene. Using 97% as the identity threshold from 118 samples with 10,000 reads per sample, 1,377 operational taxonomic units (OTUs) were identified. Principal coordinate analysis (PCoA) of weighted and unweighted UniFrac distances on fecal microbiome data revealed an alteration of the gut microbial structure of individual SPF mice with aging (Fig. 1a and b). Mean relative abundance values of multiple bacterial OTUs exhibited shifts with aging in SPF mice (Fig. 1c and Supplementary Fig. 3). Notably, the gut microbiota in wean-young mice drastically changed between 3 and 8 weeks. Hierarchical clustering analysis (HCA) of the 16S rRNA fecal microbiome datasets revealed that fecal microbes were mainly separated into two clusters consisting of those decreased (Cluster O1) and increased (Cluster O2) with aging (Supplementary Fig. 4). Alpha diversity (Shannon index) in early-middle and middle-age mice declined slightly but significantly compared to wean-young mice (Fig. 1d and Supplementary Fig. 5). The microbial genera, Allobaculum, Turicibacter (now assigned as Erysipelotrichaceae family) and unclassified genus of the S24-7 family (Bacteroidales group) were significantly increased, whereas unclassified genera of the Lachnospiraceae and Ruminococcaceae family were significantly decreased with aging (Fig. 1e and Supplementary Fig. 6).
Aging-associated alterations in the fecal metabolome profiles
We next performed metabolome analysis by capillary electrophoresis-time-of-flight mass spectrometry (CE-TOFMS) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) to identify age-associated fecal metabolite changes in GF and SPF mice. This analysis identified a total of 204 metabolites in all murine feces, of which 183 and 195 metabolites were identified in GF and SPF mice, respectively. Principal component analysis (PCA) on the fecal metabolome data showed that the profiles were clustered into GF and SPF mice groups in Principal component 1 (PC1) (Fig 2a). PCA of GF and SPF mice revealed age-related change in fecal metabolites in both groups (Fig 2b and c). To categorize the metabolite patterns in aging, we performed k-means clustering using the z-score values calculated by concentrations of metabolites in each type of mice (Fig 2d). Fecal metabolites of both SPF and GF mice separated by the k-means cluster showed change with aging; however, most of metabolites had different patterns with aging between SPF and GF mice. Comparison of the concentration of fecal metabolites from GF and SPF mice (Fig. 2d, the rightmost yellow- and green- colored heatmap) showed the large difference in the concentration of various metabolites between GF and SPF mice. Glucosamine, argininosuccinate, glyoxylate and mucate were found only in GF (Fig 2e). In contrast, 20 metabolites including butyrate were found only in SPF and gradually changed with aging (Fig 2f). The concentration of short-chain fatty acids (SCFAs: propionate, butyrate and hexanoate), vitamin B6, aspartate, alanine, lysine, cholate and sugars (maltose and maltotriose) was higher in SPF mice than in GF mice, whereas the other 14 amino acids, derivatives of amino acids and monosaccharides were lower in SPF mice (Supplementary Data 1). These results were mostly consistent with previous reports comparing fecal metabolites in GF and SPF mice[2, 33, 34].
The various fecal metabolites derived from SPF and GF mice changed gradually in association with aging (Fig 2d). Cluster 1, which contains 12 amino acids and their related metabolites involved in amino acid metabolism accumulated at 3 to 4 weeks and then decreased toward young age in SPF mice, whereas those in GF mice were comparable throughout the observed period. Cluster 2 included sugars and amino acids that accumulated at 4 to 12 weeks followed by a gradual reduction with aging in SPF and GF mice. Cluster 3 included nucleotides that accumulated at 4 to 12 weeks in SPF mice, whereas those in GF mice were tended to increase with aging. Clusters 4 and 5 gradually increased with aging in SPF mice. In contrast, cluster 4 was comparable with aging and cluster 5 was slightly accumulated in 3 to 12 weeks in GF mice. Cluster 4 contained amino acids and derivatives of amino acids and vitamins, whereas Cluster 5 included sugars and citrate, fumarate, malate, lactate, belonging to the citric acid cycle, that were especially increased in middle-aged SPF mice. In summary, profiles in SPF mice were clearly changed with aging compared to those in GF mice, and the pattern of alteration in fecal metabolites with aging was different between SPF mice and GF mice, implying that alterations in the intestinal environment associated with aging are likely caused by changes in the gut microbiota with age.
In the gut microbial metabolism, the most prominent microbial activity is the fermentation of dietary or host-derived components, in particular metabolism of non-digestible carbohydrates and host glycans into SCFAs and organic acids[35, 36]. Butyrate and hexanoate were significantly decreased (4- and 2-fold, respectively) in middle-aged SPF mice compared to young ones, whereas lactate was significantly increased by 2-fold in middle-age SPF mice compared to wean-young and young mice (Fig 2g). Moreover, sugars such as sucrose, maltose and maltotriose, thought to be derived from the diet, were also significantly increased more than 2-fold in middle-age SPF mice compared to younger mice (Fig 2h). Notably, maltose and maltotriose were detected mainly in SPF, but not GF mice, suggesting that these sugars might be derived from microbial digestion of dietary starch. We also found the concentration of soluble starch in feces was significantly increased in middle-age GF mice as compared to young GF mice, but was comparable across all ages in SPF mice (Supplementary Fig. 7). Considering the fecal metabolome profiles of GF mice, these notable fecal metabolome alterations in SPF mice imply that the aging-related phenotype of the intestinal environment is largely associated with modulation of the gut microbiota together with host aging.
Aging-associated alterations of plasma metabolome profiles
We further analyzed the plasma metabolome of these mice to capture age-associated plasma metabolic dynamics. This analysis identified 144 metabolites in total across all the murine plasma samples out of which 134 and 142 metabolites are identified in GF and SPF mice, respectively. PCA showed that the plasma metabolome profiles were roughly clustered into 2 groups of GF and SPF in the PC1 direction and aging-associated changes of GF and SPF plasma metabolome profiles were observed in the PC2 direction (Fig 3a). Furthermore, PCA calculated in each type of mouse showed aging-associated changes in plasma metabolome profiles in both GF and SPF mice (Fig 3b and c). K-means clustering of the plasma metabolites showed similar aging-associated shifts from weaning to middle-age between SPF and GF mice (Fig 3d). Argininosuccinate and 5-methoxy-3-indoleaceate were found in GF only (Fig 3e). In contrast, 10 metabolites including 5-aminovalerate and ectoine, which are thought to be mainly derived from gut microbiota[37, 38], were found only in SPF mice and most of them were gradually accumulated with aging (Fig 3f).
Amino acids and their related metabolites categorized into clusters 1, 4 and 5 were accumulated in both GF and SPF middle-aged mice (Fig 3d). The average total amino acid concentration in middle-age mice was significantly higher than that in both GF and SPF young mice (Fig 3g). Although most of the aging-associated amino acid profiles were similar between GF and SPF mice, the concentrations of glycine, serine and threonine were significantly lower in SPF mice than in GF mice almost throughout life (Supplementary Fig 8), whereas glutamine was higher in SPF mice across all ages (Fig 3h). The concentrations of pyruvate, an intermediate metabolite of glycolysis, were also significantly higher in SPF mice than in GF mice across all ages. Notably, however, the glutamate levels were comparable from young to early-middle between GF and SPF mice, and only higher in the middle-age SPF mice (Fig 3h). Cluster 2 included various types of metabolites decreasing with aging in both GF and SPF mice. Metabolites of cluster 3 were dispersed, with no change with aging.
Correlation between metabolites and microbes and transition time point analyses of SPF mice with aging
To examine aging-associated correlations among the fecal microbes, fecal metabolites and plasma metabolites in SPF mice, we calculated Spearman’s rank correlation coefficients among fecal microbe-fecal metabolite and fecal microbe-plasma metabolite pairs in a total of 119 samples (Supplementary Fig 9 and 10). Turicibacter, Allobaculum, Bifidobacterium, unclassified S24-7 and Sutterella, categorized as cluster O2, increased with aging (Supplementary Fig. 4), whereas unclassified Lachnospiraceae, Oscillospira, unclassified Clostridiales, Ruminococcus and unclassified Ruminococcaceae, categorized cluster O1, decreased and were highly and significantly correlated to fecal metabolites. In the correlations of fecal microbes and plasma metabolites, Turicibacter was highly and significantly correlated. Moreover, Procrustes analysis between the fecal metabolome and the fecal microbiome profiles of weighted UniFrac or unweighted UniFrac showed that the aging profiles in fecal metabolome and fecal microbiome are close (Supplementary Fig 11a and b). Therefore, the aging-related phenotype in the fecal metabolome might be related to taxonomic and structural changes in the gut microbiota. Although the fecal microbiome and fecal metabolome appear to be linked, changes in plasma and fecal metabolites were not similar (Supplementary Fig 11c).
We further conducted an estimation of the cluster of intestinal environment matching with age using change point analysis (CPA) based on the Bayesian Information Criterion (BIC)[39]. CPA of fecal microbiome profiles showed that principal coordinate 1 (PCo1) of unweighted and weighted UniFrac was divided into 3 clusters (unweighted) including 3 to 8 weeks, 12 to 36 weeks and 40 to 72 weeks and 3 clusters (weighted) including 3 to 6 weeks, 8 to 48 weeks and 52 to 72 weeks (Supplementary Fig 12a and b). The fecal metabolome profiles of PC1 were categorized into 2 clusters, including 3 to 36 weeks and 40 to 72 weeks (Supplementary Fig 12c), which contain common gaps of 36 to 40 weeks between fecal metabolome and unweighted fecal microbiome. These data suggest that 36 to 40 weeks might be the critical period for the gut microbial and metabolic transition into the middle age phenotype. Interestingly, CPA of the plasma metabolome showed no PC1 gap (Supplementary Fig 12d), implying that the transition toward host aging might occur gradually and continuously, which is different from the intestinal environment. To identify which metabolites and bacteria contribute to this grouping by random forest, we set two groups, 3 to 36 weeks and 40 to 72 weeks, which were determined by a common gap of 36 to 40 of PC1 scores of the fecal metabolome and PCo1 scores of unweighted UniFrac of the fecal microbiome (Supplementary Fig 13). Since fecal metabolites with a large contribution in the random forest contained metabolites mainly produced by gut microbiota, such as spermidine and ectoine [2, 37], the changes in the fecal metabolome profile are linked to the changes in fecal microbiota, which might construct this common gap. Also, Turicibacter is at the top of the contributors of random forest and highly correlated with plasma metabolites (Supplementary Fig 10), which might influence host aging. Turicibacter has already been reported as one of the bacteria increasing with aging[29]. Another report showed that the Erysipelotrichaceae family, including Turicibacter, increased in Western-type diet-induced obese mice [40-42].
Middle-age gut microbiota induces an obese phenotype in a transplanted host.
Our network analysis highlighted the features of a middle-age murine intestinal environment, such as the increase in fecal amino acids, sugars, Erysipelotrichaceae family, and the decrease in Lachnospiraceae and Ruminococcaceae families. These features are somewhat similar to the intestinal environment of an obese phenotype[19, 40, 41]. To directly ask whether the middle-age gut microbiota contributes to the obese phenotype, we orally transferred fecal microbiota from young or middle-age SPF mice into 8-week-old GF mice (Fig 4a). As expected, when placed on a high-fat diet (HFD), mice gavaged with middle-age SPF fecal microbiota displayed enhanced weight gain and higher blood glucose levels after oral glucose challenge compared to those colonized with a young SPF fecal microbiota (Fig 4b-d), although these parameters were comparable on a normal diet (Supplementary Fig. 14). Alpha diversity of the gut microbiota in the ex-GF mice gavaged with a middle-age fecal microbiota was significantly decreased after 4 weeks on a HFD (Fig 4e) and the composition of the fecal microbiota was also altered (Fig 4f). Blood insulin levels and the weight of epididymal white adipose tissue (WAT) were relatively higher but not significant, in recipients of middle-age compared to young fecal microbiota (Supplementary Fig 15). Besides, we assessed the inflammatory state of the liver, WAT and small intestine, because it was recently reported that aged gut microbiota increases intestinal permeability and promotes age-associated inflammation[7]. The expression of tumor necrosis factor α mRNA in the small intestine was significantly higher in recipients of middle-age compared to young fecal microbiota, however there was no significant difference in the expression of cytokines in the liver and WAT (Supplementary Fig 16). Taken together, these data indicate that the middle-age SPF gut microbiota has the potential to induce an obese phenotype upon high caloric intake.