Aging in SPF mice is associated with a longitudinal alterations in the structure of the gut microbiota
We monitored six SPF mice and six GF mice from three weeks to until 64 weeks (GF) and 72 weeks (SPF) in separate vinyl isolators. Six mice were divided into two group and housed in a cage by group. 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, and 12 weeks after birth, and every 4 weeks thereafter until 64 weeks (GF) and 72 weeks (SPF) (Fig. 1a). The increase of body weight in SPF mice during aging (Supplementary Fig. 1) is mostly due to an increase in fat mass, even though food intake is slightly decreased[21, 25, 26]. Weight gain in aged GF mice compared to aged SPF mice is probably due to the gradual cecum enlargement characteristic of GF mice because of the accumulation of hydrated dietary fiber components[34].
To assess the overall structure of the gut microbiota, we amplicon-sequenced the V1-V2 region of the 16S rRNA gene. Alpha diversities were negatively correlated with age in weeks (Fig. 1b), indicating that alpha diversities were decreased with age. Mean relative abundance values of multiple bacterial OTUs exhibited shifts with aging in SPF mice (Fig. 1c and Supplementary Fig. 2). Bray-Curtis distance to the data at 3 weeks of age is significantly changed according to age (Fig. 1d). Another longitudinal analysis, liner mixed effects modeling by the q2-longutidutinal package [35] showed that taxa abundance is significantly correlated with age in weeks, as evaluated by the analysis of the predictive accuracy of the regression model (R2 = 0.64, P = 2.5E− 06) (Fig. 1e). To identify taxa that are differentially abundant during aging, we applied the data of taxa to ANCOM[36] and found that the microbial genera Allobaculum, Turicibacter, Bifidobacterium, and an unclassified genus of the Coriobaceriaceae family were significantly increased, whereas the genus of Candidatus Arthromitus and Ruminococcus were significantly decreased with aging (Fig. 1f and 1g). Notably, the gut microbiota drastically changed between 3 and 6 weeks (wean-young periods); after that, it was stable from 8 to 72 weeks. Additionally, mice were transported to a new facility at 3 weeks of age, which might have had an influence on the gut microbiome. Since we expected the statistical tests and liner mixed effects modeling were probably affected by the data of wean-young periods, we reanalyzed using the 8 to 72 weeks data, excluding 3–6 weeks, and found that alpha and beta diversities were significantly changed with age (Supplementary Fig. 3a-c). We concluded from this result that fecal microbiota changed with age even after the weaning period. The data analysis from 8 to 72 weeks also showed that the ANCOM score was high only in the genera of Turicibacter, indicating that Turicibacter might be important bacteria in the young to aged phase (Supplementary Fig. 3e and f). In addition, considering the influence of experimental conditions on aging data, we evaluated the cage effect of SPF mice. We compared the Bray-Curtis distance between cages and found that there were no significant differences between cage 1 and cage 2 (Supplementary Fig. 4).
Aging-associated alterations in 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) of the fecal metabolome data showed that the profiles were significantly separated in GF and SPF mice (Adonis, R2 = 0.6, P < 0.001, Fig. 2a). PCA of GF and SPF mice revealed age-related significant changes in fecal metabolites in GF mice (Adonis, R2 = 0.54, P < 0.001, Fig. 2b) and SPF mice (Adonis, R2 = 0.55, P < 0.001, Fig. 2c).
To focus on metabolites correlated with age, we calculated Spearman’s rank correlations between fecal metabolites and age. In SPF mice, 58 metabolites (29.7% in total metabolites in feces) were positively correlated with age, whereas 65 metabolites (33.3% in total metabolites in feces) were negatively correlated with age. In GF mice, 46 metabolites (26.7% in total metabolites in feces) were positively correlated with age, whereas 62 metabolites (36.0% in total metabolites in feces) were negatively correlated with age. We showed metabolites significantly correlated with age on a heatmap with category of KEGG in SPF mice (Fig. 2d) and GF mice (Supplemental Fig. 5). Comparing the age-related change of fecal metabolites between GF and SPF mice (Fig. 2d, the right yellow- and green- colored heatmap), the correlation coefficients between GF and SPF mice were low and there were few metabolites with significant correlation, which corresponded to the PCA results. These results suggested that the fecal metabolites changed with age in both GF mice and SP mice, however the fluctuating metabolites differed between them.
To investigate which pathway associated with fecal metabolites correlated with age in SPF mice, we performed pathway enrichment analysis. Although there was no significant pathway for fecal metabolites that had a significant positive correlation with age, metabolites associated with degradation of amino acids, such as ammonia recycling and the urea cycle, were decreased with age (Fig. 2e). Indeed, except for aspartate, arginine and histidine, fecal amino acids were negatively correlated or showed no change with age (Supplemental Fig. 6).
Another feature of the fecal metabolites correlated with age in SPF mice is that monosaccharides, except for galactose, were decreased with age (Fig. 2f). In contrast, di- and tri- saccharides, except for raffinose, were increased with age (Fig. 2g). In relation to sugar and energy metabolism, monosaccharides and glucose were decreased with age, and citrate, fumarate, and malate in the downstream TCA cycle were increased with age (Supplemental Fig. 7), suggesting that energy metabolism was promoted with age. Focusing on the fermentative pathway, short-chain fatty acids were unchanged, however lactate was significantly increased with age. These results suggested that monosaccharides might be metabolized to the TCA cycle and lactate. Notably, maltose and maltotriose were detected mainly in SPF, but not in GF mice, suggesting that these sugars might be derived from microbial digestion of dietary starch. However, the concentration of soluble starch in feces was comparable across 8 to 72 weeks in SPF mice (Supplemental Fig. 7). It is possible that the intestinal bacteria involved in degradation of polysaccharides into monosaccharides were reduced with 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 157 metabolites in total across all the murine plasma samples, out of which 148 and 151 metabolites were identified in GF and SPF mice, respectively. PCA showed that the plasma metabolome profiles were significantly clustered into 2 groups of GF and SPF (Adonis, R2 = 0.42, P < 0.001) and also, aging-associated changes of GF and SPF plasma metabolome profiles were observed (Adonis, R2 = 0.44, P < 0.001) (Fig. 3a). Furthermore, PCA calculated for each type of mouse showed aging-associated changes in plasma metabolome profiles in both GF mice (Adonis, R2 = 0.66, P < 0.001, Fig. 3b) and SPF mice (Adonis, R2 = 0.73, P < 0.001, Fig. 3c). SPF mouse number 4 was an outlier in the PCA score at 44 weeks so it was excluded from the plasma analysis (Supplemental Fig. 8).
To focus on metabolites correlated with age, 69 metabolites (46.0% in total metabolites in plasma) had a significant positive correlation with age and 33 metabolites (22.0% in total metabolites in plasma) had a significant negative correlation with age. In GF mice, 35 metabolites (23.6% in total metabolites in plasma) were positively correlated with age, whereas 47 metabolites (31.8% in total metabolites in plasma) were negatively correlated with age. We showed metabolites significantly correlated with age on a heatmap with category of KEGG in SPF mice (Fig. 3d) and GF mice (Supplemental Fig. 9). In the comparison between GF and SPF mice, overall metabolites showed a positive correlation between SPF and GF (Fig. 3d, the right yellow- and green- colored heatmap), suggesting that relatively the same metabolites were changed with age in both GF and SPF mice.
A heatmap of plasma metabolites in SPF mice showed that metabolites associated with amino acid metabolism were significantly accumulated with aging (Fig. 3d). Indeed, total amino acids in plasma were gradually increased with age (Fig. 3e) and 15 amino acids were significantly correlated with age (Supplemental Fig. 10). The pathway enrichment analysis of the plasma metabolites showed that metabolites associated with amino acid metabolism, degradation of amino acids, were increased with age, in contrast, the citric acid cycle was decreased with age (Fig. 3e). These changes in plasma metabolites were opposite to that of fecal metabolites. The pathway analysis suggests that the urea cycle, which decomposes ammonia generated by amino acid catabolism, might be activated to metabolize the increase of amino acids in plasma.
Correlation among 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. Turicibacter, Allobaculum, Bifidobacterium, and an unclassified genus of the Muribaculaceae family were positively correlated with age and disaccharides (sucrose and trehalose), TCA cycle metabolites (citrate, fumarate, and malate) in feces (Fig. 4a). In contrast, an unclassified genus of the Ruminococcaceae family, Ruminococcus, an unclassified genus of the Lachnospiraceae family, Oscillospira, Desulfovibrio, and an unclassified genus of the Clostrididales family were negatively correlated with age. These bacteria also were positively and highly correlated with butyrate and negatively correlated with disaccharides and TCA cycle metabolites. In plasma-microbe correlation, few correlation pairs were detected (Fig. 4b). Turicibacter and Bifidobacterium showed a positive correlation and Ruminococcus showed a negative correlation with age.
Moreover, we also performed Procrustes analysis between Bray-Curtis dissimilarity and fecal metabolome, Bray-Curtis dissimilarity and plasma metabolome, and fecal metabolome and plasma metabolome, respectively (Supplementary Fig. 11). All comparisons had a significant correlation (P < 0.001, 999 permutations). The sum of square distances between paired points in the ordination space were 0.3891 (feces-microbiome), 0.8815 (plasma-microbiome), and 0.8754 (feces-plasma), respectively. This result indicated that the pair between feces-microbiome was highly related compared to other pairs. Therefore, the aging-related phenotype in the fecal metabolome might be related to taxonomic and structural changes in the gut microbiota.
We further conducted an estimation of the cluster of fecal profiles matching with age using change point analysis (CPA) based on a entropy based optimization[37]. We analyzed CPA for the fecal microbiome using scores of principal coordinate (PCo) 1 to 3 of Bray Curtis distance and fecal and plasma metabolome using scores of principal component (PC) 1 to 3. Sum of squared error (SSE) estimated the optimal number of clusters (Supplementary Fig. 12a-c). Interestingly, the clusters of CPA of the fecal microbiome coincided with the fecal metabolome at the three breaking points including 3 to 4 weeks, 6 to 36 weeks and 40 to 72 weeks (Supplementary Fig. 12d and e). The plasma metabolome profiles were categorized into 3 clusters, including 3 to 12 weeks, 16 to 56 weeks and 60 to 72 weeks (Supplementary Fig. 12f), which were little different from fecal microbiome and metabolome, implying that the transition toward host aging might be different from the intestinal environment. Next, we analyzed random forest to identify which metabolites and bacteria contribute to this grouping (Supplementary Fig. 13). Bacteria that have a high contribution to random forest were close to the result of ANCOM (Fig. 1g). Fecal metabolites with a large contribution in the random forest were also highly correlated to fecal microbiome (Fig. 4a), which might construct the common gaps of CPA. Looking at the estimation error of the random forest, the microbiome is 25.2% and the fecal metabolome is 3.4%, suggesting that the fecal metabolome could be used better for estimating the breaking points. Plasma metabolites with a large contribution in the random forest were also highly correlated with fecal bacteria (Fig. 4b). Although the breaking points are slightly different from the fecal microbiota, Allobaculum, Turicibacter, and Ruminococcus have a high contribution of random forest and highly correlated with plasma metabolites, which might influence host aging.
Aged gut microbiota appeared less resistant to HFD-induces obesity.
Our analysis highlighted the features of an aged murine intestinal environment, such as a decrease in monosaccharides and an increase in the metabolites of energy metabolism. In addition, plasma amino acids were increased in aged mice. These features are somewhat similar to the intestinal environment of mice with an obese phenotype[19, 38, 39]. To investigate whether the aged gut microbiota is related to the obese phenotype, we orally transferred fecal microbiota from young or aged SPF mice into young 8-week-old GF mice and maintained them on a high-fat diet (HFD) to develop HFD-induced obesity (Fig. 5a). The GF mice gavaged with aged 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. 5b-d), although these parameters were comparable on a normal diet (Supplementary Fig. 14). Blood insulin levels and the weight of epididymal white adipose tissue (WAT) were relatively higher but not significant, in recipients of aged compared to young fecal microbiota (Supplementary Fig. 15). We also 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-α (TNF-α) mRNA in the small intestine was relatively higher but did not reach significance (P = 0.0611) in recipients of aged compared to young fecal microbiota, however there was no significant difference in the expression of cytokines in the liver and WAT (Supplementary Fig. 16). We additionally examined the expression of inflammatory genes in the small and large intestines, however there were no significant differences between recipients of aged and young fecal microbiota (data not shown). Next, we compared the composition of the fecal microbiota (Fig. 5e). Alpha diversity and beta diversity were comparable in recipients of aged and young fecal microbiota after 4 weeks on a HFD diet (Fig. 5f and 5g), however, the comparison of pairwise distance at 0 and 4 weeks after HFD feeding showed that recipients of aged fecal microbiota were altered compared to those of young fecal microbiota (Fig. 5h), suggesting that young fecal microbiota were more stable than aged fecal microbiota during HFD feeding. In addition, Parabacteroides, an unclassified genus of the Erysipelotrichaceae family, Sutterella, an unclassified genus of the Muribaculaceae family, Allobaculum, Candidatus Arthromitus, and an unclassified genus of the Peptostreptococcaceae family were significantly different at 4 weeks after HFD feeding (Fig. 5i). The increase of Allobaculum and the decrease of Candidatus Arthromitus were also observed in aging of SPF mice (Fig. 1g). From these results, some of the age-related changes in fecal bacteria were detected as features even in recipients of aged fecal microbiota under the HFD. Taken together, these data indicate that the aged SPF gut microbiota has the effect of promoting HFD-induced obesity.