
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39237540
52097
10.1038/s41467-024-52097-9
Article
Gut metagenomes of Asian octogenarians reveal metabolic potential expansion and distinct microbial species associated with aging phenotypes
http://orcid.org/0000-0003-2391-0202
Ravikrishnan Aarthi 1
Wijaya Indrik 1
Png Eileen 1
http://orcid.org/0000-0003-2746-9880
Chng Kern Rei 1
Ho Eliza Xin Pei 1
Ng Amanda Hui Qi 1
Mohamed Naim Ahmad Nazri 1
http://orcid.org/0000-0003-2368-8779
Gounot Jean-Sebastien 1
Guan Shou Ping 2
http://orcid.org/0000-0001-7779-8135
Hanqing Jasinda Lee 2
http://orcid.org/0000-0003-0497-8519
Guan Lihuan 2
http://orcid.org/0000-0003-1182-6804
Li Chenhao 1
Koh Jia Yu 1
de Sessions Paola Florez 1
http://orcid.org/0000-0002-5674-6341
Koh Woon-Puay 23
Feng Lei 2
Ng Tze Pin 2
Larbi Anis 4
http://orcid.org/0000-0001-7206-1724
Maier Andrea B. 25
http://orcid.org/0000-0002-5754-1874
Kennedy Brian K. 2
http://orcid.org/0000-0003-0850-5604
Nagarajan Niranjan nagarajann@gis.a-star.edu.sg

12
1 https://ror.org/05k8wg936 grid.418377.e 0000 0004 0620 715X Genome Institute of Singapore (GIS), Agency for Science, Technology and Research (A*STAR), 60 Biopolis Street, Genome, Singapore, 138672 Republic of Singapore
2 https://ror.org/01tgyzw49 grid.4280.e 0000 0001 2180 6431 Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Republic of Singapore
3 grid.185448.4 0000 0004 0637 0221 Singapore Institute for Clinical Sciences (SICS), Agency for Science Technology and Research (A*STAR), 30 Medical Drive, Brenner Centre for Molecular Medicine, Singapore, 117609 Republic of Singapore
4 grid.185448.4 0000 0004 0637 0221 Singapore Immunology Network (SigN), Agency for Science Technology and Research (A*STAR), 8A Biomedical Grove, Immunos, Singapore, 138648 Republic of Singapore
5 https://ror.org/008xxew50 grid.12380.38 0000 0004 1754 9227 Department of Human Movement Sciences, @AgeAmsterdam, Faculty of Behavioural and Movement Sciences, Vrije Universiteit Amsterdam, Amsterdam Movement Sciences, Amsterdam, The Netherlands
5 9 2024
5 9 2024
2024
15 775126 11 2023
23 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
While rapid demographic changes in Asia are driving the incidence of chronic aging-related diseases, the limited availability of high-quality in vivo data hampers our ability to understand complex multi-factorial contributions, including gut microbial, to healthy aging. Leveraging a well-phenotyped cohort of community-living octogenarians in Singapore, we used deep shotgun-metagenomic sequencing for high-resolution taxonomic and functional characterization of their gut microbiomes (n = 234). Joint species-level analysis with other Asian cohorts identified distinct age-associated shifts characterized by reduction in microbial richness, and specific Alistipes and Bacteroides species enrichment (e.g., Alistipes shahii and Bacteroides xylanisolvens). Functional analysis confirmed these changes correspond to metabolic potential expansion in aging towards alternate pathways synthesizing and utilizing amino-acid precursors, vis-à-vis dominant microbial guilds producing butyrate in gut from pyruvate (e.g., Faecalibacterium prausnitzii, Roseburia inulinivorans). Extending these observations to key clinical markers helped identify >10 robust microbial associations to inflammation, cardiometabolic and liver health, including potential probiotic species (e.g., Parabacteroides goldsteinii) and pathobionts (e.g., Klebsiella pneumoniae), highlighting the microbiome’s role as biomarkers and potential targets for promoting healthy aging.

Here, the authors employ shotgun-metagenomic sequencing to profile the microbiome of Singaporean octogenarians, finding key microbial shifts and metabolic changes in the gut, in particular expanded butyrate production from alternate precursors and identifying microbes linked to health in aging populations.

Subject terms

Computational biology and bioinformatics
Computational models
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Over the last few decades, economic changes and advances in healthcare systems have significantly improved life expectancies in Asia1,2, and are expected to lead to a rapid shift in demographics by doubling the population of 60+ individuals by 20503 (from 19.7% in 2015 to 40%). This has been accompanied by a rising incidence of chronic diseases linked to aging4,5 and corresponding socio-economic stress on healthcare systems across Asia6,7. There is therefore an urgent need to identify lifestyle, dietary, and pharmaceutical interventions that promote healthy aging in Asian populations8–10.

Aging is believed to be a complex, multi-factorial phenomenon with progressive decline in several physiological functions including in the gastrointestinal and immune system11,12. Not surprisingly, many studies have therefore identified correlations between gut microbiome composition and age13,14. For example, in recent work, Zhang et al.15 have shown that there are distinct age and gender-associated trajectories in the gut microbiome that are consistent across populations. Several studies have also highlighted that frailty indices and healthy aging are linked to the state of the gut microbiome16,17, along with known gut-microbiome associations with aging-related diseases18,19. In particular, shifts in several families and genera of bacteria including Bacteroidetes, Proteobacteria, Roseburia and Escherichia have been associated via 16S rRNA sequencing to healthy aging20–26. Work by Wilmanski et al. also suggests that compositional uniqueness of the gut microbiome may serve as a marker for healthy aging27. To compensate for the potential low resolution of 16S studies, recent work has sought to leverage shotgun metagenomic sequencing to characterize the gut microbiome in elderly populations, but have been limited by cohort size (n < 50) and the spectrum of age ranges that they study, in their ability to account for potential confounders while identifying microbial mechanisms associated with healthy aging28,29.

To address this, we used deep shotgun metagenomic analysis (n = 234, >15 million reads on average) to study gut microbiomes in a cohort of community-living octogenarians (primarily, age range = [71–100], Table 1) in Singapore30,31. The gut microbiomes of healthy octogenarians exhibited a defined shift in diversity and overall taxonomic composition, as a function of age and independent of other covariates. In addition, taxonomic analysis identified several species-level changes (e.g., enrichment of Alistipes shahii and Bacteroides xylanisolvens), consistent with functional pathway analysis of the data, that points to a metabolic potential expansion in aging towards alternate pathways that synthesize and utilize amino-acids (e.g., L-lysine) as precursor substrates, relative to the classical pyruvate-to-butyrate production pathway. As butyrate derived from gut bacteria has diverse roles in host health (e.g., as energy for colonocytes, and reducing gut inflammation32–34), we next associated gut microbiome composition with key markers for inflammation (e.g., CRP), cardiometabolic health (e.g., fasting blood glucose) and liver health (e.g., AST), to find additional microbial functions associated with healthy aging. Leveraging extensive clinical data and the size of the metagenomic datasets, we identified >10 robust associations (accounting for demographic and clinical covariates that can be confounders) that highlight the role of the microbiome as biomarkers and potential intervention targets for promoting healthy aging.Table 1 Clinical characteristics of the SG90 cohort

	All subjects	Male	Female	
Demographics	
 Number of subjects	234	95	139	
 Age (Mean ± SD)	86 ± 4	86 ± 3	86 ± 4	
Serum measurements (Mean ± SD)	
 High-sensitivity CRP (mg/L)	3.7 ± 8.4	4.6 ± 11.9	3.1 ± 4.4	
 Fasting blood glucose (mmol/L)	5.2 ± 0.4	5.3 ± 0.5	5.2 ± 0.4	
 Triglycerides (mmol/L)	1.4 ± 0.6	1.2 ± 0.6	1.5 ± 0.7	
 Total cholesterol (mmol/L)	4.6 ± 1.1	4.4 ± 1.0	4.8 ± 1.1	
 High-density lipoprotein (mmol/L)	1.4 ± 0.3	1.3 ± 0.3	1.4 ± 0.4	
 Low-density lipoprotein (mmol/L)	2.6 ± 0.9	2.5 ± 0.8	2.7 ± 0.9	
 Alanine aminotransferase (U/L)	18.4 ± 8.0	19.3 ± 8.1	17.8 ± 8.0	
 Aspartate aminotransferase (U/L)	25.0 ± 6.8	25.5 ± 7.6	24.6 ± 6.2	
Disease incidence (%)	
 Type 2 diabetes	18	12	22	
 Heart attack	6	9	5	
 Heart failure	4	4	4	
 Cataracts/glaucoma	36	31	38	
 Kidney failure	0.9	1	0.7	
 Asthma	3	1	4	
 COPD	2	4	0.7	
 Arthritis	9	5	12	
 Osteoporosis	4	1	6	
 Hip fracture	4	1	6	
 Parkinson’s disease	0.9	2	0	
 Other neurodegenerative disorders	0.4	1	0	
 Gastrointestinal problems	6	5	7	
 Thyroid problems	2	3	1	
 Cancer	3	2	3	

Results

Shotgun metagenomics reveals a defined aging-associated shift in microbial richness and functional guilds

To study the role of the gut microbiome in healthy aging, stool samples from a deeply phenotyped cohort of community-living Asian octogenarians (SG90, Table 1, Supplementary Data File 1) were analyzed with deep shotgun metagenomic sequencing (n = 234, 17 million reads on average, Methods). The resulting high-resolution species-level taxonomic profiles were compared to reference shotgun metagenomic data from healthy, younger Singaporeans35,36 from two cohorts (SPMP, n = 109, age range = [53–74]; CPE, n = 96, age range = [21–80]), as well as other Asian populations37 (T2D, n = 171, age range = [21–70]; Supplementary Data File 2). Systematic joint analysis of these cohorts after batch-correction and matching other demographic characteristics (gender, ethnicity; Methods), revealed a progressive shift in taxonomic profiles across different age groups (n = 516, age range = [21–100], Fig. 1A), particularly along the y = x axis relative to the first and second principal components of variation (ordinal logistic test p-value < 0.001, Fig. 1B, Supplementary Fig. 1A). Diversity analysis indicated that this was accompanied by an aging-associated shift (ordinal logistic test p-value < 0.001, Fig. 1C) that is driven by microbial richness and evenness, both of which exhibit significant reduction with age (Supplementary Fig. 1B–D). PERMANOVA analysis showed that Age was the major source of variation compared to other attributes (Supplementary Fig. 2A) and that beta diversity peaked first in the 41–60 age group and then again in the 91–100 age group (Supplementary Fig. 2B), though uniqueness did not consistently increase with age as reported previously in a western cohort27. These results highlight an aging-associated shift in microbial composition that is defined by the loss of gut species and increase in beta diversity in extreme age groups.Fig. 1 Aging-associated shifts in gut microbiome richness and relative abundance of key species.

A Principal coordinates analysis (PCoA) plot based on species-level Bray–Curtis dissimilarity of gut microbiome profiles across age groups after batch correction. Dashed lines indicate the y = x and y = –x axes. B, C Violin plots showing the distribution across age groups. B PCoA1 and PCoA2 values projected on the y = x axes (n = 516 independent samples; The symbol ‘***’ represents p-value < 0.005, obtained from two-sided test with Ordinal logistic regression across all age groups, p-value 8.09 × 10–33). In all the inset boxplots, the center line represents the median, box limits represent upper and lower quartiles, and whiskers represent minimum and maximum values. C Shannon diversity indices; The symbol ‘***’ represents p-value < 0.005, obtained from two-sided test with Ordinal logistic test using n = 516 independent samples across all age groups, p-value 1.58 × 10–12). In all the inset boxplots, the center line represents the median, box limits represent upper and lower quartiles, and whiskers represent minimum and maximum values. Data points that are either less than Q1–1.5 × IQR or more than Q3 + 1.5 × IQR, where Q1, Q3 and IQR refers to the first quartile, third quartile and interquartile range, respectively, were removed for this analysis to reduce the impact of outliers. D Volcano plot showing the β coefficient on the x-axis and FDR-adjusted –log10 p-values (FDR-adjusted using Benjamini–Hochberg test, p-values obtained with two-sided tests from Generalized Linear Model (GLM) for taxa association with age) on the y-axis. Points for statistically significant taxa from “high stringency list” and those with high beta coefficients are highlighted in the volcano plot as red colored points for median abundance in the SG90 cohort greater than zero and green otherwise, with non-significant taxa shown with grey dots. E Relative abundance boxplots for the Bacteroides genus, and corresponding species that was identified to be associated significantly and not significantly with age. The symbols “n.s.”, “*” and “***” represent p-value > 0.05, p-value ≤ 0.05 and p-value < 0.005, respectively (GLM for taxa association with age using n = 516 independent samples; FDR-adjusted p-value using Benjamini–Hochberg test; Bacteroides: 6.4 × 10–3, B. dorei: 7.61 × 10–5). In all boxplots, the center line represents the median, box limits represent upper and lower quartiles, and whiskers represent minimum and maximum values (outlier points are not included in the visualization). Source data are provided as a Source Data file.

In order to identify taxa associated with age, a generalized linear model (GLM) was used to account for demographic and clinical covariates that can be confounders (e.g., gender, body mass index, fasting blood glucose, triglyceride, total cholesterol and high-density lipoprotein; Methods). Overall, 4 phyla, 22 genera, and 40 species were found to be associated significantly with age (False Discovery Rate (FDR)-adjusted p-value < 0.05, Supplementary Table 1, Supplementary Data File 3). The phylum and genus level results were broadly consistent with prior observations (e.g., in Chinese38, Japanese25 and Italian centenarians28,39), with enrichment in Bacteroidetes (β = 0.33, FDR-adjusted p-value < 0.001), and depletion in Firmicutes (β = –0.49, FDR-adjusted p-value < 0.001) with age (Supplementary Table 1). At the genus level, a strong enrichment in unclassified Acidaminococcaceae species (β = 6.69, FDR-adjusted p-value < 0.05) and depletion in Parasutterella species (β = –20.19, FDR-adjusted p-value < 0.001) was observed that has not been noted before, while at the species level a majority of our observations were novel owing to the higher resolution of shotgun-metagenomics20–26 (Supplementary Data File 3, Table 2). In addition, we tested the robustness of this analysis by varying the normalization approach, statistical model and cohorts used to find that most taxa associations were detected in a majority of the conditions tested (Methods, Supplementary Fig. 3). We also noted that while several species associations were detected in individual cohorts (n = 10), the inclusion of the SG90 cohort notably boosted the number of age-associations identified (n = 38), suggesting that having a wider age range and more subjects could benefit such analysis (Supplementary Fig. 4). To mitigate the potential effect of incomplete batch-correction, a high-stringency list of species associations was further derived based primarily on methods that do not depend on this procedure (n = 17; Table 2, Fig. 1D, Supplementary Fig. 3; Methods). Reanalysis of taxonomic profiles and inclusion of results from other publicly available cohorts17,40–44 (n = 6; Methods) highlighted the consistency of these associations with a high percentage being replicated or validated (27/40 = 68%, odds ratio = 21.4, p-value < 6 × 10–15). In addition, the overlap with Asian cohorts (18/40, odds ratio = 8.83, p-value < 10–6) was more statistically significant than with western cohorts (16/40, odds ratio = 3.10, p-value < 7.1 × 10–3), despite having more western cohorts in this comparison (4/6), highlighting the value of having more Asian datasets focused on healthy aging (Supplementary Fig. 5).Table 2 Species-level gut microbial taxa associated with aging

Enriched with age	FDR-adjusted p-value	β coefficient for age [95% C.I.]	Remarks	
Bacteroides xylanisolvens	6.71 × 10–3	2.28 [1.05, 3.51]		
Alistipes shahii	6.81 × 10–4	1.53 [0.99, 2.07]		
Bacteroidales bacterium ph8	1.13 × 10–4	1.51 [0.83, 2.19]		
Klebsiella pneumoniae	8.92 × 10–3	1.12 [0.39, 1.86]		
Escherichia unclassified	3.73 × 10–5	1.02 [0.59, 1.45]	16S OTU enriched in Indian centenarians154	
Parabacteroides merdae	3.42 × 10–2	0.95 [–0.0087, 1.90]		
Bacteroides dorei	7.61 × 10–5	0.68 [0.42, 0.94]		
Alistipes onderdonkii	1.60 × 10–3	0.66 [0.30, 1.02]	Enriched in Italian centenarians28,29	
Escherichia coli	3.10 × 10–2	0.27 [0.054, 0.48]		
Depleted with age	
Bilophila unclassified	6.30 × 10–5	–6.82 [–9.90, –3.74]		
Ruminococcus torques	6.14 × 10–3	–1.22 [–1.99, –0.45]		
Faecalibacterium prausnitzii	5.33 × 10–6	–1.01 [–1.29, –0.74]	Depleted in Italian centenarians28,154	
The table reports statistically significant taxa based on FDR-adjusted p-value ≤ 0.05 using Benjamini–Hochberg test, p-values obtained with two-sided tests from GLM for species association with age, after covariate adjustment (n = 516 independent samples). Species with median relative abundance greater than zero in the SG90 cohort are reported here and the full list is provided in Supplementary Data File 3. The β coefficient estimates strength of association as change in age (in years) for every 1% increment of taxa relative abundance. Taxa in the high-stringency list are shown in bold face.

Among the core gut microbiome taxa identified here, only some species in the Bacteroides genus showed significant positive associations with age27,45 (Table 2). This included Bacteroides dorei and B. xylanisolvens, with B. massiliensis even exhibiting a slight negative trend, though this variation is masked at the genus level highlighting the utility of species-level resolution offered by shotgun metagenomics (Fig. 1E). We also noted an enrichment with age for multiple Alistipes species, including A. shahii, A. onderdonkii and A. senegalensis, that are bile tolerant, have a unique way of fermenting amino acids such as L-lysine46,47, and can produce neurotransmitter precursors such as indole48 (Supplementary Fig. 6A). In contrast, several species that are known to be important for butyrate production in the gut46,47, including Faecalibacterium prausnitzii, Roseburia inulinivorans and Eubacterium hallii, were significantly depleted with age (Fig. 1D, Supplementary Fig. 6B). Several of these associations were seen in other Asian cohorts but not in the western cohorts analyzed here (n = 11, e.g., Dorea formicigenerans, R. inulinivorans, A. shahii and E. hallii) indicating that there may be consistent key differences in ethnic Asian populations. Also, in a few cases the directionality of association differed between Asian and western cohorts (e.g., B. xylanisolvens), highlighting the importance of population-specific studies. Integrating this information through taxonomic co-occurrence analysis highlighted a distinct network structure in the elderly (32 edges and 3 hubs) compared to younger individuals (80 edges and 9 hubs), with a shared pattern of cliques formed by the classical butyrate producers and Alistipes species, indicating that these may represent distinct functional guilds49 that switch roles with age (Supplementary Fig. 7).

Gut microbiomes in elderly Asians exhibit a metabolic potential expansion in butyrate synthesis pathways

To further analyze metabolic changes associated with aging, shotgun metagenomic data from all subjects from the four cohorts (SG90, SPMP, CPE and T2D) was used to quantify gene and pathway abundances (Supplementary Data File 4, Methods). In total, 413 metabolic pathways were identified and quantified of which 70 were found to be differentially abundant across various age groups (Supplementary Data File 5), broadly grouping into the four categories of sugar metabolism (17 pathways), vitamin, energy metabolism, and cell wall biosynthesis (28 pathways), lipid metabolism (9 pathways) and amino acid metabolism (16 pathways, Fig. 2). Specifically, among sugar metabolism pathways we observed distinct age-associated biases, where e.g., pyruvate fermentation to acetate and lactate was enriched in the younger age group (21–40), while degradation of mono-, di- and polysaccharides was enriched in older age groups (Fig. 2A). As many as 8 different pathways for various substrates were enriched in the 61–70 age group, including simple pentose sugars (Pentose phosphate pathways), and complex polysaccharides such as stachyose and starch (Starch degradation V). Notably, microbial glycolytic pathways for metabolism of simple sugars (glucose) were enriched in healthy individuals over 90, while prior work had highlighted the role of corresponding pathways in aging in model organisms50.Fig. 2 Associations of microbial metabolic pathways across age groups.

Bar plot showing pathways that are significantly associated and enriched specifically in different age groups (p-value ≤ 0.05, two-sided Wilcoxon rank-sum test) based on LEfSe analysis. Pathways were grouped into the broad categories of (A) Sugar metabolism, B Vitamin, energy metabolism and cell wall biosynthesis, C Lipid metabolism and (D) Amino acid metabolism. The effect size is depicted in the form of an LDA Score (log 10). Source data are provided as a Source Data file.

In line with the vital role of vitamins as essential micronutrients (particularly in energy metabolism and immune function51,52) with limited synthesis in humans51, several microbial pathways were associated with aging (especially B vitamins, Fig. 2B). In particular, thiamine diphosphate and flavin biosynthesis were enriched in the most elderly (91–100) and could have anti-inflammatory and anti-aging roles52. On the other hand, three alternate biosynthetic pathways to produce thiamin and pantothenate were enriched in younger individuals (21–40).

While microbial lipid metabolism pathways were predominant in younger subjects (21–40 age group, Fig. 2C), pathways related to amino acid metabolism were more often enriched in the elderly (particularly 61–70 age group), including the synthesis of various essential amino acids such as histidine, valine, threonine, isoleucine, tryptophan, and methionine and aromatic amino acids (Fig. 2D). In addition, lysine biosynthesis was over-represented in the gut microbiomes of two older age groups53 (61–70 and 81–90), linking to our earlier observation of age-related enrichment in specific Bacteroides and Alistipes species, which are known to synthesize lysine and produce short-chain fatty acids (SCFA) such as butyrate with lysine as substrate46, respectively.

Analysis of the four major known pathways for butyrate production46,47 (from pyruvate, glutarate, 4-amino butyrate and lysine) confirmed that while several lysine-associated genes show an age-related increasing trend (e.g., kamA, FDR-adjusted p-value < 0.01 and kdd, FDR-adjusted p-value < 0.05), no significant associations were seen in the pyruvate related genes (Fig. 3). Furthermore, the overall gene abundance in the pathway for butyrate production from lysine also showed a significant increase with age (FDR-adjusted p-value < 0.001, Supplementary Fig. 8). Metabolic pathways analysis revealed that the abundance of gut microbial species that provide metabolic support to butyrate producers peaked in older age groups (71–80; Supplementary Fig. 9, Methods). To more directly link taxonomic and functional analysis, we first analysed the contribution of taxa to L-lysine biosynthesis III pathway and observed that Bacteroides sp. were the dominant contributors, with the contribution of B. xylanisolvens increasing with age (Supplementary Fig. 10A, 10B). Next, we systematically generated metagenome-assembled genomes, annotated genes and pathways, and conducted metabolic network analysis (Methods). This analysis showed that the metagenome-assembled genomes represented a diversity of Bacteroides and Alistipes species, with many genomes harboring corresponding complete pathways for lysine biosynthesis (Bacteroides) and lysine-to-butyrate conversion (Alistipes) that are reachable based on metabolic network analysis (Supplementary Fig. 10C–E). These observations, along with those in the previous section, indicate a metabolic potential expansion in butyrate synthesis pathways from precursors other than pyruvate in the gut microbiomes of elderly Asians.Fig. 3 Differential representation of butyrate synthesis pathways in gut metagenomes across age groups.

Pathway diagram where nodes represent major metabolites involved in various steps of butyrate production (but) starting from 4 different precursors (pyruvate, glutarate, lysine and 4-aminobutyrate), and edges are labelled with the genes involved in the conversion. Boxplots present normalized counts for the corresponding gene across different age groups (n = 516 independent samples). The labels “n.s.”, “*”, “**” and “***” represent p-value > 0.05, p-value ≤ 0.05, p-value < 0.01 and p-value < 0.001, respectively (two-sided tests using GLM for gene association with age with FDR-adjustment using Benjamini–Hochberg test per pathway; gcdAB – 2.68 × 10–16, kamA – 9.97 × 10–5, kamDE – 1.65 × 10–2, kdd – 1.51 × 10–3, kce – 1.65 × 10–2, abfH – 3.93 × 10–6). In all the boxplots, the center line represents the median, box limits represent upper and lower quartiles, and whiskers represent minimum and maximum values. Expansions of the abbreviated metabolites are as follows: pep - Phosphoenolpyruvate, pyr - Pyruvate, accoa - Acetyl-CoA, aacoa - Acetoacetyl-CoA, 3hbcoa - (S)-3-Hydroxybutanoyl-CoA, akg - 2-Oxoglutarate, 4h2oglt - 4-Hhydroxy-2-oxoglutarate, 2hydogc - 2-Hydroxy-glutaryl-CoA, 3mgcoa - 3-Methylyglutaconyl-CoA, lys – L-lysine, 36dahx - (3S)-3,6-Diaminohexanoate, 35dahx - (3S)-3,5-Diaminohexanoate, 5a3ohex - 5-Amino-3-oxohexaonate, 3abutcoa - L-3-Aminobutyryl-CoA, 4abut - 4-Aminobutyrate, succsal - Succinic semialdehyde, 4hdxbutn - 4-Hydroxybutyrate, 4-hbutcoa - 4-Hydroxybutyrate-CoA, 3btcoa - Vinylacetyl CoA, b2coa - Crotonyl CoA, but - Butyrate. Source data are provided as a Source Data file.

To further explore these findings, we leveraged the recent development of a healthy aging model in mice, where dietary supplementation with alpha ketoglutarate (AKG) to aged mice (18 months) not only increased lifespan but also improved inflammation and frailty markers54. Shotgun metagenomics of fecal samples from AKG-supplemented and control mice (n = 20 per group, Supplementary Fig. 11A; Methods) showed that even though mice harbor distinct gut microbial species, at the pathway level a similar expansion in butyrate production pathways is seen in healthy aging mice (Supplementary Fig. 11B). In particular, the lysine to butyrate conversion pathways shows the strongest enrichment, as seen in our human data (1.65-fold increase in median abundance, FDR-adjusted p-value < 0.01). At the gene level, several key genes for lysine to butyrate production were strongly enriched including kdd, kamDE, kce and kamA (FDR-adjusted p-value < 0.05; Supplementary Fig. 11C). Such a trend was not observed for genes involved in pyruvate to butyrate conversion, highlighting the role of expansion in alternate butyrate synthesis pathways to support healthy aging.

Robust species-level associations between gut microbial taxa, clinical biomarkers and aging-associated phenotypes

Aging is a significant risk factor for chronic diseases impacting multiple organ functions, including metabolic, immune, and musculoskeletal systems, which could be mediated via gut microbiome function26,55. To further relate the taxonomic and functional changes observed in Asian octogenarians with their clinical phenotypes we harnessed the higher resolution of shotgun metagenomic data, broader age range, and size of our cohorts, to robustly identify associations in a GLM framework while accounting for demographic and clinical confounders (e.g., age, gender, other metabolic markers; Methods, Supplementary Data File 6). Despite the multiplicity of taxa and clinical features, we noted that this analysis typically resulted in a few significant associations (Fig. 4). For example, for fasting blood glucose as a well-studied biomarker for glucose metabolism56, insulin resistance57, and type 2 diabetes58, we identified two associated species (Bifidobacterium adolescentis [β = 0.08, FDR-adjusted p-value = 2.61 × 10−2] and Parabacteroides goldsteinii [β = 0.48, FDR-adjusted p-value = 2.63 × 10−2], Fig. 4A) with prior in vitro and in vivo observations supporting this link59–61, as well as a new association with a Lactobacillus species (Lactobacillus mucosae: β = 0.88, FDR-adjusted p-value = 4.30 × 10−2) that has never been described before despite numerous studies for this genus62–65. Note that β values indicate the change in response variable (e.g., fasting blood glucose levels) for every unit change in relative abundance values, after accounting for other confounding variables. Given the known normal range of fasting blood glucose levels (<6.1 mmol/L66; Supplementary Table 2) the observed β values (0.48 and 0.88) therefore could suggest a large effect size.Fig. 4 Microbiome associations with key clinical markers.

Volcano plots showing the β coefficient on the x-axis and the FDR-adjusted –log10 p-values on the y-axis (FDR-adjusted using Benjamini–Hochberg test, p-values obtained with two-sided tests from GLM) for gut microbial associations with clinical markers for (A) Glucose control as measured by fasting blood glucose levels (Bifidobacterium adolescentis: 2.61 × 10–2, Parabacteroides goldsteinii: 2.63 × 10–2, Lactobacillus mucosae: 4.30 × 10–2), B Lipid control as measured by total cholesterol in plasma (Haemophilus parainfluenzae: 7.95 × 10–4, Lachnospiraceae bacterium 1 4 56FAA: 3.92 × 10–3), C Liver health, as measured by plasma levels of aspartate aminotransferase (Klebsiella pneumoniae: 3.99 × 10–3), and (D) Systemic inflammation, as measured by high-sensitivity CRP levels (Eggerthella unclassified: 4.45 × 10–2, Escherichia coli: 4.45 × 10–2, Streptococcus salivarius: 4.45 × 10–2), E Vitamin B12 (Streptococcus parasanguinis: 4.54 × 10–9, Bacteroides coprocola: 9.76 × 10–5). Points in all plots are colored with red for FDR-adjusted p-value ≤ 0.05 and grey otherwise. Dotted lines represent the thresholds of p-value = 0.05 and p-value = 0.1. Source data are provided as a Source Data file.

We next analyzed cholesterol biomarkers (e.g., high-density lipoproteins—HDL, low-density lipoproteins—LDL and triglycerides) for hypercholesterolemia and cardiovascular diseases, due to their impact on healthy aging67. Association analysis with total cholesterol levels revealed two significant candidates, positive association with abundance of Lachnospiraceae bacterium 1 4 56FAA (β = 1.02, FDR-adjusted p-value = 3.92 × 10–4) and negative association with Haemophilus parainfluenzae (β = –0.06, FDR-adjusted p-value = 7.95 × 10–5, Fig. 4B), neither of which has been reported previously68,69. Interestingly, analysis with LDL levels revealed the opposite trend, with the Lachnospiraceae species (β = –0.88, FDR-adjusted p-value = 4.27 × 10–2) being negatively correlated and H. parainfluenzae (β = 0.06, FDR-adjusted p-value = 1.72 × 10–3) being positively correlated, and similar trends being observed with HDL but not reaching statistical significance (Supplementary Fig. 12A, B). These observations indicate potential beneficial roles of Lachnospiraceae species in maintaining lower levels of LDL. Interestingly, microbial association with triglyceride levels was dominated again by positive association with H. parainfluenzae (Supplementary Fig. 12C), a species that has been described to be strongly enriched in patients with chronic cholecystitis patients70,71 and individuals with increased risk of cardiovascular disease72,73.

As the biological processes of aging are in part believed to be accelerated by inflammation74,75 (Inflammaging), including risk for cardiovascular disease, frailty and disability76–78, we sought to identify gut microbial associations with inflammation markers. Specifically, serum levels of aspartate aminotransferase (AST) were associated with the abundance of Klebsiella pneumoniae (β = 0.65, FDR-adjusted p-value = 3.99 × 10–3, Fig. 4C), a sporadic colonizer of the gut microbiome, and a known pathogen for liver abscesses that are increasingly common in Asian populations79. A similar association was seen with alanine aminotransferase (ALT), a more specific biomarker for liver injury, but it did not reach statistical significance (Supplementary Fig. 12D). Analysis with high-sensitivity C-reactive protein (hs-CRP) levels, a commonly used marker for systemic inflammation80, identified a known association with Escherichia coli81 (β = 0.0020, FDR-adjusted p-value = 4.45 × 10−2), as well as an association with a Streptococcus species (S. salivarius: β = 0.02, FDR-adjusted p-value = 4.45 × 10−2) which has been linked previously with rheumatoid arthritis82 (Fig. 4D). Furthermore, a significant association was observed with a Eggerthella species (β = 0.15, FDR-adjusted p-value = 4.45 × 10−2) distinct from Eggerthella lenta83, whose role is less well-understood68, and the association with host inflammation has not been described before. These results highlight the various bacterial species and potential pathogenic functions by which altered gut microbiomes can contribute to inflammaging.

Aging is a common risk factor that often predisposes elderly to the development of neurological disorders such as Parkinson’s and Alzheimer’s disease84. Many studies have shown that Vitamin B12 is involved in several processes of the nervous system, including the synthesis of myelin and post-injury nerve regeneration85,86, the deficiency of which causes peripheral neuropathy87, gait ataxia and physical frailty that are often observed in the elderly. A strong association was observed between serum levels of Vitamin B12 and the bacterial species Streptococcus parasanguinis (β = 312.50, FDR-adjusted p-value < 5 × 10–9, Fig. 4E), strains of which are known to produce Vitamin B1288, highlighting the metabolic potential for gut bacteria to support healthy aging. Intriguingly, we also observed a significant association for blood levels of Vitamin B12 and the abundance of Bacteroides coprocola (β = 42.89, FDR-adjusted p-value < 10–4), even though Bacteroidetes are not common B12 producers89 and instead compete with the host to capture B12 in the gut90, indicating that dietary intake might be an explanation for the observed positive association. Together, these observations highlight the two-way communication between gut microbial species and their environment that can impact host phenotypes as well as serve as meaningful biomarkers of healthy aging.

Discussion

The complexity of the human gut microbiome, both in terms of its genetic diversity and the myriad host-microbial interactions that shape it, necessitates that association studies be done with sufficiently high-resolution and power to help infer meaningful insights about function and mechanisms91. While prior studies were primarily limited by their resolution (16S rRNA surveys20–26), recent studies have underscored the utility of shotgun metagenomics for species-level and functional pathway analysis, albeit in limited-sized cohorts28,29. By generating deep shotgun metagenomic data for >200 elderly subjects and with a combined dataset of >500 subjects from all cohorts, this study significantly enhances our ability to explore the role of the microbiome in healthy aging. In particular, we provide the first well-powered species-level view of microbiome shifts as a function of age, with a steady decrease in species richness (Fig. 1A, Supplementary Fig. 1B). This is accompanied by an increase in beta diversity in extreme age groups, but further defines that this is due to an overall shift of microbial community composition with age, with specific species serving as markers (e.g., A. shahii, B. xylanisolvens, R. obeum and E. hallii). Our analysis also highlights that the age-associated shift is on an axis of variation (Fig. 1A, B) that is distinct from the axis separating different gut microbiome clusters reminiscent of enterotype structures92 (Supplementary Fig. 1A). In addition, while the overall reduction in microbial diversity observed here (Fig. 1C) is broadly consistent with some studies24,93, our analysis clarifies that this is accompanied by an increase in beta diversity of microbiome composition and a loss of species richness (Supplementary Fig. 2B, 1B–D), defined in part by several classical butyrate producers in the gut (e.g., F. prausnitzii, R. inulinivorans and E. rectale). We did not observe a general loss in core Bacteroides species as reported in Wilmanski et al.27, which could be a function of our Asian cohorts, but may also be due to methodological differences (e.g., covariate adjustment). Also, our results suggest that the trends reported for richness22,94 and uniqueness27 may be variable across cohorts95. Further studies are needed across different global populations to explore the factors impacting this variability using consistently generated datasets and potentially using long-read sequencing to enable more reliable detection/assembly of rare taxa. This could also impact the detection of associations with viruses and fungi (none reported here), which are only sporadically abundant in the human gut microbiome.

In order to have data from a wide range of ages and a large number of subjects to sufficiently power our age-association analysis, we combined data from multiple Asian cohorts (SG90, T2D, SPMP, CPE) with metagenomic sequencing data, and attempted to account for differences in how the data was generated using a batch-correction approach (Methods, Supplementary Table 4). In addition, our association analysis accounted for several demographic and clinical factors that could be confounders, including ethnicity, gender, body mass index (BMI), fasting blood glucose, triglyceride, total cholesterol and high density lipoprotein levels. Nevertheless, it is important to note that batch-correction methods cannot be expected to fully account for all technical sources of variation, and similarly, other potential confounders (e.g., diet, medications) that could not be accounted for as this data was not consistently available across cohorts. Multi-cohort association analysis thus necessarily has to evaluate the robustness of the results obtained as a function of various data subsets used, technical choices made (e.g., normalization) and statistical approaches used. Our experiments suggest that the batch-correction and GLM-based association analysis reported here provide robust results that are reproduced across a majority of conditions tested (Methods, Supplementary Fig. 3). This was further validated by the strong enrichment of age associations detected in other cohorts, particularly in Asian cohorts (Supplementary Fig. 5). In addition, our analysis highlights the importance of having sufficient data from older age groups (Supplementary Figs. 4 and 5). Overall, despite all analytical care, the interpretation of microbiome associations with aging needs care as it is a multi-factorial phenomenon, where e.g., increased lifetime use of antibiotics could explain loss in gut microbial species richness with age, while at the same time such observations could arise from other factors such as lifestyle choices (e.g., smoking) and their independent impact on longevity and the gut microbiome. Our study is structured to identify microbiome associations with age and considers BMI, total cholesterol, fasting blood glucose etc as confounders, though in principle, the influence of age on these variables could be through the microbiome. In that sense, our study is conservative and should identify only the strongest independent associations. These associations can then be further explored in mouse models where we hypothesize that the specific Bacteroides and Alistipes species identified here could promote healthy aging independently or in combination through their metabolic outputs. In addition, the associations identified here can serve as useful biomarkers that in combination with other known clinical biomarkers can help track individuals who are aging in a healthier manner96,97 as well as develop better frailty metrics to identify individuals at risk for aging-associated diseases98.

The importance of incorporating age and other covariates in microbiome association analyses has been emphasized in several recent studies17,99. While this also necessitates larger cohort sizes and extensive metadata collection, a potential advantage is that the few associations that survive covariate adjustments are likely to be more meaningful (e.g., no associations were detected for sleep duration, frailty markers or free thyroxine levels; Supplementary Data File 6). In combination with shotgun metagenomics, where specific species and functional genes/pathways can be directly quantified, and related through metagenome-assembled genomes, this provides for a powerful avenue for developing meaningful mechanistic hypotheses informed by in vivo human data. In particular, a striking feature of our results is the depletion of a ‘microbial guild’ of butyrate-producing species in elderly subjects and its replacement with specific Alistipes species that can produce butyrate through an alternate pathway using amino acids as precursors. At the metabolic potential level, we do not see a significant reduction for the pyruvate-to-butyrate conversion pathway, and it is not clear if dietary or microbially produced amino acids are the key source for the alternate pathways100. Prior studies have shown reduction in SCFA levels, particularly butyrate, in stool samples from the elderly101, association between butyrate levels and frailty14, and hypothesized a shift towards bacterial protein fermentation due to increased intestinal transit time102. We observed a relative reduction of microbial to host reads in elderly subjects (Supplementary Fig. 13) that could reflect a reduced capability to support gut microbial fermentation, though direct measurements are needed to confirm this indirect evidence. Changes in dietary patterns (e.g., low fiber) as well as intestinal physiology (e.g., muscle loss, hypochlorhydria) could thus partly explain the observed shift in species composition and deserves further investigation. Our data from a mouse model of healthy aging further supports the idea that alternate butyrate production pathways in the gut microbiome could play a causal role in promoting healthy aging. The model is established by feeding mice α-ketoglutarate which is rapidly absorbed in the small intestine103, but may still explain the enrichment seen in the glutarate pathway in the gut microbiome. The strong enrichment seen in the lysine-to-butyrate pathways is consistent with our observations in human cohorts and could reflect the systemic metabolic and physiological effects of α-ketoglutarate in promoting healthy aging104, though it is not clear if this effect could be mediated through the gut microbiome. Follow-up experiments in mice with defined microbiome compositions, with and without α-ketoglutarate and L-lysine supplementation, could help further test the hypothesis that alternative butyrate production pathways have a causal role in healthy aging.

As the microbial production of SCFAs in the gut is known to play an important role in many aspects of human health105,106 (e.g., gut barrier function, immune regulation, gut-brain communication), it is not surprising therefore that reduction in butyrate-producing species has been associated with frailty107. Based on the observations in this study we hypothesize that alternate pathways for butyrate production, including using glutarate, 4-amino butyrate, and lysine as precursors, may partially augment the pyruvate-dependent pathway in the gut microbiomes of elderly individuals to promote healthy aging. These observations may be due to a “selection on survival” bias27, but results from the mouse model suggest that this may not provide a full explanation. Lysine is an essential amino acid that is found in lower plasma concentrations in frail elderly subjects108, and lysine supplementation is being investigated in multiple clinical trials for its health benefits. Its use as a supplement to promote healthy aging could thus be another avenue to explore via the pathway of boosting butyrate production in the gut. Among other pathways that were relatively depleted in the elderly, reduction in pyruvate fermentation to acetate and lactate and lipid metabolism pathways (Fig. 2) may have detrimental effects for host health and can be targets for compensation through supplements and probiotics109,110. On the other hand, the relative enrichment of vitamin synthesis pathways in the elderly, including thiamine diphosphate (vitamin B1), flavin (vitamin B2), and phospho-pantothenate (vitamin B5), as well as the association of specific species with serum vitamin B12 levels (Fig. 4E), could contribute to healthy aging and thus deserves further study.

Several aging-related chronic diseases, including type 2 diabetes37, cardiovascular diseases111, and neurological conditions112 have been linked with the gut microbiome and thus it is essential to account for age as a confounding factor when studying gut microbial mechanisms that contribute to disease. We hypothesize that having individuals from as wide a range of ages as possible, particularly those who have aged well, will likely help strengthen and refine such association analysis. We detected a positive association of B. adolescentis with fasting blood glucose, which is well-known to impact glycemic control113. In addition, we also noted a significant association for P. goldsteinii, a potential new probiotic species with emerging evidence for a diet-dependent role in obesity and type 2 diabetes59,114, though the mechanisms for this effect in mouse models remain to be elucidated.

Systemic low-levels of inflammation are believed to contribute to various aging-related phenotypes74,75 (Inflammaging). Correspondingly, associations between gut microbial species and markers of inflammation (e.g., hs-CRP and AST) are of particular interest. While the association of a common-source of liver abscesses in Asia (K. pneumoniae) with a liver disease marker (AST), and a common oral bacterium (S. salivarius) with a systemic inflammation marker (hs-CRP) are consistent with disease biology, the distinctness of these associations is striking i.e., while K. pneumoniae is the only species associated with AST, a more diverse group of bacteria are associated with hs-CRP levels. Furthermore, relatively little is known about the Lachnospiraceae species with strong associations to cholesterol/LDL levels. Together with the strong association of microbial biomarkers such as S. parasanguinis and B. coprocola with serum Vitamin B12 levels, these findings could help develop a non-invasive frailty test based on at-home sample collection115. Overall, these results highlight the value of species-level shotgun metagenomic analysis in large well-characterized cohorts, where further studies in other populations would help determine if some of these associations are specific to Asian environments and lifestyles.

Methods

Cohorts and datasets

SG90 cohort

The SG90 cohort is based on a longitudinal population health study that was set up in the 1990s which involved routine measurement of metabolic and other health variables30,31. The current dataset is based on a subset of 234 elderly individuals (77–97 years old) who are community-living participants (not living in a nursing home, no diagnosis of dementia and not physically unfit) and consented to providing their stool and blood samples. No statistical method was used to predetermine sample size. The participants in the SG90 cohort were recruited under the SLAS-3116 protocol approved by the Institutional Review Board (IRB) at National University of Singapore (reference number: B-15-081). This study was also approved by an IRB for the Singapore Chinese Health Study30 at National University of Singapore (reference number: H-17-027, sample collection: 2019/00439). A written informed consent was obtained from all the participants. Fasting blood glucose (mmol/L), triglyceride (mmol/L), total cholesterol (mmol/L), HDL (mmol/L), LDL (mmol/L), hs-CRP (mg/L), AST (U/L), ALT (U/L) and Vitamin B12 levels (pmol/L) were measured based on a blood draw at NUH Referral Laboratories, Singapore. Stool samples were collected either at the time of blood collection or within a week and were transported in refrigerated conditions to the lab where it was frozen at –80 °C. In addition, physical assessments were performed to calculate BMI (kg/m2) from their weight (kg) and height (m), as well as gait speed (m/s) and handgrip strength (kg) of subjects were measured.

SPMP dataset

The SPMP dataset is based on the recall of a subset of 109 healthy Singaporean subjects (53–74 years old) from a multi-omics study in Singapore117. Stool samples were collected for gut microbiome analysis using shotgun metagenomic sequencing35.

T2D dataset

Shotgun metagenomic datasets were obtained for 171 healthy Chinese individuals from a previously published type 2 diabetes (T2D) study37,118 using the curatedMetagenomicData package118 (v1.16.0). Briefly, the subjects chosen for our study were 21–70 years old and were non-diabetic controls in the study. Clinical data such as fasting blood glucose (mmol/L), triglyceride (mmol/L), total cholesterol (mmol/L), HDL (mmol/L) and LDL (mmol/L) levels were also obtained from this study37.

CPE dataset

The CPE dataset36 is based on a prospective cohort study consisting of CPE-colonized subjects and their healthy family members. For our comparisons, we used shotgun metagenomic data of 82 healthy family members (21–80 years old) with Chinese ethnicity.

Demographic matching

Analysis for all cohorts was restricted to ethnic Chinese individuals and the gender balance across cohorts was found to be comparable (SG90: 59% female, T2D: 52% female, SPMP: 60% female, CPE: 60% female). Age distributions for all cohorts can be found in Supplementary Fig. 14.

DNA library construction and sequencing

PowerSoil DNA Isolation Kit (12888-100, MO Bio Laboratories) was used for the extraction of DNA from frozen stool samples. Minor modifications to the manufacturer’s protocol were made (double the volume of C2, C3, and C4 buffers was added, and the duration of the centrifugation step was extended to twice the original duration). Purified DNA was eluted in 80 µL of C6 solution. DNA libraries were prepared using 50 ng of extracted DNA re-suspended in a volume of 50 µl. This was subjected to shearing using Adaptive Focused AcousticsTM (Covaris) with the following parameters—Duty Factor: 30%, Peak Incident Power (PIP): 450, 200 cycles per burst, Treatment Time: 240 s. Sheared DNA was cleaned up with 1.5× Agencourt AMPure XP beads (A63882, Beckman Coulter) followed by end-repair, A-addition and adapter ligation using the Gene Read DNA Library I Core Kit (180434, Qiagen) according to the manufacturer’s protocol. Custom barcode adapters were used instead of GeneRead Adapter I Set for adapter ligation (Supplementary Table 3). Before enrichment, DNA libraries were cleaned twice using 1.5× Agencourt AMPure XP beads (A63882, Beckman Coulter) using the protocol from Multiplexing Sample Preparation Oligonucleotide kit (PE-400-1001, Illumina). Enrichment PCR was carried out with PE 1.0 and custom index-primers for 12 cycles. DNA Libraries were prepared with Agilent DNA1000 Kit (5067-1504, Agilent Technologies) by pooling equimolar concentrations and quantified using Agilent Bioanalyzer (G2939BA). A pilot set of 50 DNA libraries were sequenced on an Illumina HiSeq 2500 and the rest were sequenced on an Illumina HiSeq X instrument (Supplementary Data File 1) generating >17 million 2 × 101 bp reads on average per library. The sequencing protocol and the median sequencing depth of other cohorts can be found in Supplementary Table 4.

Read pre-processing and profiling

Illumina shotgun metagenomic sequencing reads were processed using a Nextflow pipeline (https://github.com/CSB5/shotgunmetagenomics-nf). Briefly, raw reads were filtered to remove low-quality bases and adapter sequences were removed using fastp119 (v0.20.0) with default parameters. Human reads were removed by mapping to the hg19 reference using BWA-MEM120 (v0.7.17-r1188, default parameters) and samtools121 (v1.7). The remaining reads were used for taxonomic profiling using MetaPhlAn2122 (v2.7.7, default parameters). Functional profiles for the metagenomes were obtained using HUMAnN2123 (v2.8.1). For all statistical tests, Benjamini–Hochberg’s false discovery rate method was used to correct for multiple testing at a significance threshold of 5%.

Statistical analysis with taxonomic profiles

Taxonomic profiles were corrected for batch effects with MMUPHin124 (v1.8.1) using age group as a covariate. MMUPHin is a batch correction method that employs an empirical Bayes approach to model read counts with respect to batch variables and biologically relevant covariates (age group in our case). It then gives as output batch-corrected count data which aims to retain the effects of biologically relevant covariates. For batch correction, we considered samples from the SG90 and CPE cohorts as belonging to the same batch due to their similarity in DNA extraction and library preparation methods i.e., the major sources of variation in metagenomic profiles125 (Supplementary Table 4). To assess if differences in sequencing platform and depth could affect this assumption, we performed differential abundance analyses using MaAsLin2126 (v1.10.0) (correcting for age) between (i) subsets of the SG90 cohort sequenced on the two different platforms (Illumina HiSeq 2500 and Illumina HiSeq X), and (ii) down-sampled version of the CPE dataset with median sequencing depth matching the SG90 cohort compared to the original CPE dataset. In both analyses we did not identify any differentially abundant taxa between the respective groups, supporting our treatment of the two cohorts as the same batch. We conducted a sanity-check on batch-corrected profiles to confirm that cohort effects were reduced (Supplementary Fig. 15). Species-level relative abundances (>0.1%) were then used to compute alpha diversity indices (Shannon and Simpson), Pielou’s evenness and the beta diversity index (Bray–Curtis distance) using the R package vegan127 (v2.6-4) and visualized using ggplot2128 (v3.5.1). Robustness of computed metrics to variations in sequencing depth was confirmed based on sub-sampling and correlation analysis (Supplementary Fig. 16). For sub-sampling analysis, we selected a subset of 30 samples each with at least 15 million reads and randomly subsampled 7, 14, and 21 million reads from each of these samples. We then generated the taxonomic and pathway profiles for each of these sequencing depths using the same Nextflow pipeline and computed Spearman correlation values for alpha diversity metrics and pathway counts with the values obtained originally (Supplementary Fig. 16). Ordinal logistic regression in R was used to test for statistical significance in relation to different age groups. A GLM approach was used to test associations between the relative abundance for each species (SA) and age as a continuous variable accounting for covariates.

For this analysis, all demographic and clinical covariates available in multiple cohorts (n = 7 covariates) were first screened for relevance as potential confounders based on previously reported studies69,81,129 (all retained), and then through Spearman correlation within covariates to exclude collinear variables. Specifically, we selected and included only variables that were not highly correlated with one another, retaining only one of the covariates if their correlation was greater than 0.9 (one covariate removed, i.e., LDL levels). The final set of covariates used in this analysis (n = 6) included gender (G), BMI, fasting blood glucose (FBG), triglycerides (TGL), total cholesterol (TC), and HDL levels (i.e., with the formula: Age ~ SA + G + BMI + FBG + TGL + TC + HDL). Only taxa with >50 non-zero values were tested to avoid spurious associations. Missing clinical covariates were imputed through an Expectation Maximization algorithm130.

To assess reproducibility of taxonomic associations, three distinct approaches were used: (i) GLM analysis using all four cohorts with batch correction as described above, (ii) GLM analysis using the original data for the SG90 and CPE cohorts as they were similarly processed, and (iii) Trend analysis using all four cohorts with the Cochran-Armitage test131 after conversion of relative abundance data into presence-absence values (cutoff of 0.1%). Microbial associations in same direction found in two out of the three methods were considered reproducible with the reported p-values being the median across all methods. In order to further test the robustness of associations, parameters and options were varied for all three approaches including normalization technique (Total Sum Scaling—TSS and Cumulative Sum Scaling—CSS), association analysis technique (ANCOM-BC132, MaAsLin2126) and relative abundance cutoffs (0.05%–1%; Supplementary Fig. 3). To obtain a “high stringency” list of microbial associations we restricted to taxa within the SG90 + CPE cohort that are either (i) reported by all four conditions or (ii) reported by at least three conditions, in which case the presence-absence tests are required to report an association at all tested cutoff values (Supplementary Fig. 3). This list avoids the possibility of a false positive that emerges solely due to batch correction artefacts. The final taxa in this list were visualized using ggplot2 in R. CCREPE133 (v1.0) was used to compute Spearman correlation values and identify bacterial species with strong co-occurrence patterns in different age groups (|ρ| > 0.2, p-value < 0.05).

To further assess the robustness and reproducibility of our associations, we identified microbial associations across six additional cohorts of Asian and Western origin (Supplementary Fig. 5). Metagenomic taxonomic profiles from stool samples for four of these cohorts (labelled as YachidaS_201940, AsnicarF_201741, BritoIL_201642 and LLD_201643) were obtained from curatedMetagenomicData118 (v3.8.0). The subjects in these cohorts were healthy and were selected for not having antibiotic usage. Microbial associations and taxonomic profiles for two other cohorts (labelled as XuWuZhu_202244 and ElderMet_202017) were obtained from their respective publications. All of these cohorts were analyzed independently. We employed a similar GLM framework (with the formula: Age ~ SA + G) as outlined previously to independently identify microbial associations with age. Microbial associations from our discovery cohort (SG90 + SPMP + T2D + CPE) were considered ‘replicated’ or ‘validated’ in other cohorts based on their statistical significance (FDR-adjusted p-value < 0.05 or p-value < 0.05, respectively).

Pathway and gene-level analysis

The HMP Unified Metabolic Analysis Network (HUMAnN2)134 (v 2.8.1) pipeline was used to determine the relative abundance of microbial pathways in different gut metagenomes. The default Kyoto encyclopedia of genes and genomes (KEGG)135 catalog was used as the pathway reference. Unstratified relative abundance values for SG90, CPE and SPMP were integrated with HUMAnN2 results for T2D from curatedMetagenomicData118 (v1.16.0) for shared pathways, and significant differentially abundant pathways across age groups were determined based on linear discriminant analysis with LEfSe136 (p-value < 0.05 and LDA score >3).

For gene-level analysis of metabolic pathways involved in butyrate production, EC numbers of enzymes were mapped to corresponding UniRef90137 identifiers from the HUMAnN2 output. The list of butyrate producing species, genes, and enzymes (EC numbers) in butyrate production pathways were obtained from previously published studies46,47. Using this list, the butyrate production genes of interest from our dataset were obtained using a three step process: (i) UniRef90 identifiers obtained from HUMAnN2, were first cross-referenced with the KEGG database to obtain a list of KEGG orthology (KO) groups, (ii) EC numbers obtained previously were matched with the KO groups (Supplementary Data File 7), and (iii) the KO groups were cross-referenced with UniRef90 to obtain the specific butyrate producing genes of interest. For every gene, read counts were obtained by taking the sum for all corresponding UniRef90 identifier-associated values and statistical analyses were performed using R with visualizations generated using ggplot2128.

Metagenomic assembly

Microbial reads from the SG90 cohort were assembled and binned following an in-house snakemake pipeline. Briefly, the reads were assembled using MEGAHIT138 (v1.2.9). Contigs were binned using MetaBAT2139 (v2.15) based on contig coverage information generated with minimap2140 (v2.26). Assemblies were evaluated based on MIMAG141 definitions, with contamination, completeness, and N50 values obtained from CheckM2142 (v1.0.2), and non-coding RNA identified using barrnap (https://github.com/tseemann/barrnap) (v0.9) and tRNAScan-SE143 (v2.0.12). Potential cases of chimerism were detected using GUNC144 (v1.0.5). All software were used with default arguments. Low-quality MAGs were removed from downstream analysis and taxonomic classification of MAGs at the species level was obtained using GTDB-Tk145 (v2.3.2) ‘classify_wf’ method (GTDB version r214). Genomes were annotated using Prokka146 (v1.14.6) and the EC numbers associated with each gene were extracted. Metabolic pathways were obtained using MinPath147 (v1.6). Reachability analysis was carried out using MetQuest148 (v0.1.33) using the bipartite graphs constructed from these metabolic pathways. The initial seed list for the reachability analysis included source metabolites L-aspartate for L-lysine biosynthesis III and L-lysine for L-lysine to butyrate pathway.

Mouse experiments

All the mouse experiments were carried out with the approval of the Institutional Animal Care and Use Committee (IACUC) of National University of Singapore (NUS) under the R21-0135 protocol. Two independent groups of 16-month-old male and female C57BL6/J mice were purchased from Jackson Laboratories (RRID: IMSR JAX:000664) and housed at 25 °C with 50% humidity on a 12-h light/dark cycle. The animals were housed separately (maximum of four mice per cage), with n = 20 mice per group and were maintained on a regular chow diet until they were 18 months old. At 18 months of age, mice in the control group continued the regular chow diet while those in the healthy aging group were switched to a diet containing 4% calcium alpha-ketoglutarate (Ponce de Leon Health), with stool samples collected at baseline and 3 months after diet switch for metagenomic analysis. Genomic DNA was extracted from mouse stools using QIAamp PowerFecal Pro DNA Kit (51804, Qiagen), according to manufacturer’s instructions. DNA was quantified on a Qubit Fluorometer using the Qubit dsDNA BR Assay Kit (Q32850, ThermoFisher Scientific). Purified genomic DNA (50 ng) was used for library construction steps using NEBNext® Ultra™ II FS DNA Library Prep Kit (E7805L) according to manufacturer’s instructions. Finally, each library sample was assessed for quality based on fragment size and concentration using the Agilent D1000 ScreenTape (G2991BA) system, with samples adjusted to identical concentrations by means of dilution and volume-adjusted pooling. The multiplexed sample pool was paired-end (2 × 151 bp) sequenced on an Illumina HiSeq X Ten system. Mouse gut metagenomes were annotated using eggNOG-mapper149 (v2.1.9). Briefly, reads were mapped against eggNOG protein database (eggnog 5.0150) using DIAMOND151 (v2.0.15.153) in blastx mode. The default KEGG135 catalog was used as the pathway reference. For gene-level analysis of metabolic pathways involved in butyrate production, the corresponding KO IDs of each gene were used from the eggNOG-mapper’s output. For every gene, read counts were obtained by taking the sum for all corresponding KEGG values. Gene normalized abundances were measured as logarithm of counts per millions [log2(CPM + 1)] to adjust for differences in sequencing depths. Pathway normalized abundances were calculated by taking the sum of all of gene counts for each pathway and log-transformed and normalized to sequencing depth as log2(CPM + 1). Statistical tests for groups comparisons were done using two-sided Wilcoxon rank-sum test (FDR adjusted p-value < 0.05 considered to be significant) and were performed using R with visualizations created using ggplot2.

Metabolic network analysis

Metabolic support index (MSI) values for the species identified in Table 2 in relation to other gut microbial species were obtained based on metabolic network analysis, as described previously152. Briefly, MSI uses network flow analysis to quantify the extent of microbial metabolism that is enabled by the presence of other species in the community. The metabolic support network was visualized as a directed graph using Cytoscape153 (v3.8.0). Species that are known butyrate producers46,47 and receive the most metabolic support (based on in-degree of nodes) were highlighted in the network. The combined relative abundances of organisms supporting these butyrate producers were compared across age groups using the Wilcoxon test.

Phenotypic association analysis

Associations of microbial species with clinical phenotypic markers were identified based on linear regression using ‘glm’ function in R, with age and other markers serving as covariates for adjustment. Only species present in at least ten samples were considered for this analysis. Phenotypic markers available in two datasets (T2D and SG90) included BMI, fasting blood glucose, triglycerides, total cholesterol, HDL and LDL levels. For hs-CRP, AST, and ALT levels, gait speed, right handgrip strength, left handgrip strength, sleep duration, MMSE and free thyroxine levels, the associations were only tested for SG90 samples (without batch-correction) where this data was available.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Data file 1

Supplementary Data file 2

Supplementary Data file 3

Supplementary Data file 4

Supplementary Data file 5

Supplementary Data file 6

Supplementary Data file 7

Reporting summary

STORMS checklist

Source data

Source data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52097-9.

Acknowledgements

We thank Siew-Hong Low of the National University of Singapore for supervising the fieldwork in the Singapore Chinese Health Study-SG90 sub-cohort. The SG90 study was jointly supported by the Agency for Science, Technology and Research (A*STAR) under the A*STAR-AMED Joint Grant Call in “Biological and molecular determinants of human ageing (AMED2016-SG002)”, and the NUHS-A*STAR Biology of Human Ageing Program seed funding. This work was also supported by core funding from A*STAR. Authors also acknowledge the funding from National Research Foundation, Singapore (NRFI09-0015, N.N.) and Singapore National Medical Research Clinician Scientist Individual Research Grant (CIRG22jul-0023, J.S.G). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The computational work for this article was partially performed on resources of the National Supercomputing Centre (NSCC), Singapore (https://www.nscc.sg).

Author contributions

N.N., A.L., and P.F.S. planned and designed the project. G.L., W.P.K., L.F., T.P.N., A.B.M., B.K. contributed the SG90 clinical cohort. G.S.P., H.J.L., performed the mouse experiments. E.X.P.N., A.H.Q.N., A.M.N.N. conducted the wet-lab experiments. A.R., I.W., compiled all the data. A.R., I.W., E.P., K.R.C., J.S.G., C.L., K.J.Y. performed the computational data analysis with the guidance of N.N., A.B.M., and B.K. A.R., I.W., and N.N. wrote the manuscript with input from all the authors.

Peer review

Peer review information

Nature Communications thanks Marco Candela, Yan Shao and the other, anonymous, reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

Shotgun metagenomic sequencing data are available from the European Nucleotide Archive (ENA – https://www.ebi.ac.uk/ena/browser/home) under project accession number PRJEB49124. Source data are provided with this paper. Clinical data used in the study is available on request. All the requests should be directed to the corresponding author of this study, N.N. Source data are provided with this paper.

Code availability

Source code for scripts used to analyze the data are available in a GitHub project at https://github.com/CSB5/SG90.

Competing interests

The authors declare no competing interests.

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

These authors contributed equally: Aarthi Ravikrishnan, Indrik Wijaya.
==== Refs
References

1. Goldman DP Consequences of health trends and medical innovation for the future elderly Health Aff. 2005 24 W5R5 W5R17 10.1377/hlthaff.W5.R5
Goldman, D. P. et al. Consequences of health trends and medical innovation for the future elderly. Health Aff. 24, W5R5–W5R17 (2005).10.1377/hlthaff.W5.R5
2. Rowe, J. W., Fulmer, T. & Fried, L. Preparing for better health and health care for an aging population. JAMA 316, 1643–1644 (2016).
3. Malhotra R The aging of a young nation: population aging in Singapore Gerontologist 2019 59 401 410 30517628
Malhotra, R. et al. The aging of a young nation: population aging in Singapore. Gerontologist 59, 401–410 (2019).30517628
4. Jaul E Barron J Age-related diseases and clinical and public health implications for the 85 years old and over population Front. Public Health 2017 5 335 10.3389/fpubh.2017.00335 29312916
Jaul, E. & Barron, J. Age-related diseases and clinical and public health implications for the 85 years old and over population. Front. Public Health 5, 335 (2017).29312916 10.3389/fpubh.2017.00335
5. Brody JA Grant MD Age-associated diseases and conditions: implications for decreasing late life morbidity Aging 2001 13 64 67 11405387
Brody, J. A. & Grant, M. D. Age-associated diseases and conditions: implications for decreasing late life morbidity. Aging 13, 64–67 (2001).11405387
6. Cheng X Population ageing and mortality during 1990-2017: a global decomposition analysis PLoS Med. 2020 17 e1003138 10.1371/journal.pmed.1003138 32511229
Cheng, X. et al. Population ageing and mortality during 1990-2017: a global decomposition analysis. PLoS Med. 17, e1003138 (2020).32511229 10.1371/journal.pmed.1003138
7. Divo MJ Martinez CH Mannino DM Ageing and the epidemiology of multimorbidity Eur. Respir. J. 2014 44 1055 1068 10.1183/09031936.00059814 25142482
Divo, M. J., Martinez, C. H. & Mannino, D. M. Ageing and the epidemiology of multimorbidity. Eur. Respir. J. 44, 1055–1068 (2014).25142482 10.1183/09031936.00059814
8. Kiefte-de Jong JC Mathers JC Franco OH Nutrition and healthy ageing: the key ingredients Proc. Nutr. Soc. 2014 73 249 259 10.1017/S0029665113003881 24503212
Kiefte-de Jong, J. C., Mathers, J. C. & Franco, O. H. Nutrition and healthy ageing: the key ingredients. Proc. Nutr. Soc. 73, 249–259 (2014).24503212 10.1017/S0029665113003881
9. Santos, F. H., Zurek, J. & Janicki, M. P. Efficacy of healthy aging interventions for adults with intellectual and developmental disabilities: a systematic review. Gerontologist 62, e235–52 (2020).
10. Mackowiak PA Recycling metchnikoff: probiotics, the intestinal microbiome and the quest for long life Front. Public Health 2013 1 52 10.3389/fpubh.2013.00052 24350221
Mackowiak, P. A. Recycling metchnikoff: probiotics, the intestinal microbiome and the quest for long life. Front. Public Health 1, 52 (2013).24350221 10.3389/fpubh.2013.00052
11. Ponnappan S Ponnappan U Aging and immune function: molecular mechanisms to interventions Antioxid. Redox Signal 2011 14 1551 1585 10.1089/ars.2010.3228 20812785
Ponnappan, S. & Ponnappan, U. Aging and immune function: molecular mechanisms to interventions. Antioxid. Redox Signal 14, 1551–1585 (2011).20812785 10.1089/ars.2010.3228
12. Haynes L Aging of the immune system: research challenges to enhance the health span of older adults Front. Aging 2020 1 602108 10.3389/fragi.2020.602108 35822168
Haynes, L. Aging of the immune system: research challenges to enhance the health span of older adults. Front. Aging 1, 602108 (2020).35822168 10.3389/fragi.2020.602108
13. O’Toole PW Claesson MJ Gut microbiota: changes throughout the lifespan from infancy to elderly Int. Dairy J. 2010 20 281 291 10.1016/j.idairyj.2009.11.010
O’Toole, P. W. & Claesson, M. J. Gut microbiota: changes throughout the lifespan from infancy to elderly. Int. Dairy J. 20, 281–291 (2010).10.1016/j.idairyj.2009.11.010
14. Claesson MJ Gut microbiota composition correlates with diet and health in the elderly Nature 2012 488 178 184 10.1038/nature11319 22797518
Claesson, M. J. et al. Gut microbiota composition correlates with diet and health in the elderly. Nature 488, 178–184 (2012).22797518 10.1038/nature11319
15. Zhang X Sex- and age-related trajectories of the adult human gut microbiota shared across populations of different ethnicities Nat. Aging 2021 1 87 100 10.1038/s43587-020-00014-2 37118004
Zhang, X. et al. Sex- and age-related trajectories of the adult human gut microbiota shared across populations of different ethnicities. Nat. Aging 1, 87–100 (2021).37118004 10.1038/s43587-020-00014-2
16. Ragonnaud E Biragyn A Gut microbiota as the key controllers of “healthy” aging of elderly people Immun. Ageing 2021 18 2 10.1186/s12979-020-00213-w 33397404
Ragonnaud, E. & Biragyn, A. Gut microbiota as the key controllers of “healthy” aging of elderly people. Immun. Ageing 18, 2 (2021).33397404 10.1186/s12979-020-00213-w
17. Ghosh TS Das M Jeffery IB O’Toole PW Adjusting for age improves identification of gut microbiome alterations in multiple diseases eLife 2020 9 e50240 10.7554/eLife.50240 32159510
Ghosh, T. S., Das, M., Jeffery, I. B. & O’Toole, P. W. Adjusting for age improves identification of gut microbiome alterations in multiple diseases. eLife 9, e50240 (2020).32159510 10.7554/eLife.50240
18. Askarova S The links between the gut microbiome, aging, modern lifestyle and Alzheimer’s disease Front. Cell Infect. Microbiol. 2020 10 104 10.3389/fcimb.2020.00104 32257964
Askarova, S. et al. The links between the gut microbiome, aging, modern lifestyle and Alzheimer’s disease. Front. Cell Infect. Microbiol. 10, 104 (2020).32257964 10.3389/fcimb.2020.00104
19. Nagpal R Gut microbiome and aging: physiological and mechanistic insights Nutr. Healthy Aging 2018 4 267 285 10.3233/NHA-170030 29951588
Nagpal, R. et al. Gut microbiome and aging: physiological and mechanistic insights. Nutr. Healthy Aging 4, 267–285 (2018).29951588 10.3233/NHA-170030
20. Biagi E Through ageing, and beyond: gut microbiota and inflammatory status in seniors and centenarians PLoS ONE 2010 5 e10667 10.1371/journal.pone.0010667 20498852
Biagi, E. et al. Through ageing, and beyond: gut microbiota and inflammatory status in seniors and centenarians. PLoS ONE 5, e10667 (2010).20498852 10.1371/journal.pone.0010667
21. Biagi, E. et al. The gut microbiota of centenarians: signatures of longevity in the gut microbiota profile. Mech. Ageing Dev. 165, 180–184 (2017).
22. Kong F Gut microbiota signatures of longevity Curr. Biol. 2016 26 R832 R833 10.1016/j.cub.2016.08.015 27676296
Kong, F. et al. Gut microbiota signatures of longevity. Curr. Biol. 26, R832–R833 (2016).27676296 10.1016/j.cub.2016.08.015
23. Wang F Gut microbiota community and its assembly associated with age and diet in Chinese centenarians J. Microbiol. Biotechnol. 2015 25 1195 1204 10.4014/jmb.1410.10014 25839332
Wang, F. et al. Gut microbiota community and its assembly associated with age and diet in Chinese centenarians. J. Microbiol. Biotechnol. 25, 1195–1204 (2015).25839332 10.4014/jmb.1410.10014
24. Claesson MJ Composition, variability, and temporal stability of the intestinal microbiota of the elderly Proc. Natl Acad. Sci. USA 2011 108 4586 4591 10.1073/pnas.1000097107 20571116
Claesson, M. J. et al. Composition, variability, and temporal stability of the intestinal microbiota of the elderly. Proc. Natl Acad. Sci. USA 108, 4586–4591 (2011).20571116 10.1073/pnas.1000097107
25. Odamaki T Age-related changes in gut microbiota composition from newborn to centenarian: a cross-sectional study BMC Microbiol. 2016 16 90 10.1186/s12866-016-0708-5 27220822
Odamaki, T. et al. Age-related changes in gut microbiota composition from newborn to centenarian: a cross-sectional study. BMC Microbiol. 16, 90 (2016).27220822 10.1186/s12866-016-0708-5
26. O’Toole PW Jeffery IB Gut microbiota and aging Science 2015 350 1214 1215 10.1126/science.aac8469 26785481
O’Toole, P. W. & Jeffery, I. B. Gut microbiota and aging. Science 350, 1214–1215 (2015).26785481 10.1126/science.aac8469
27. Wilmanski T Gut microbiome pattern reflects healthy ageing and predicts survival in humans Nat. Metab. 2021 3 274 286 10.1038/s42255-021-00348-0 33619379
Wilmanski, T. et al. Gut microbiome pattern reflects healthy ageing and predicts survival in humans. Nat. Metab. 3, 274–286 (2021).33619379 10.1038/s42255-021-00348-0
28. Wu L A cross-sectional study of compositional and functional profiles of gut microbiota in Sardinian centenarians mSystems 2019 4 e00325 19 10.1128/mSystems.00325-19 31289141
Wu, L. et al. A cross-sectional study of compositional and functional profiles of gut microbiota in Sardinian centenarians. mSystems 4, e00325–19 (2019).31289141 10.1128/mSystems.00325-19
29. Rampelli S Shotgun metagenomics of gut microbiota in humans with up to extreme longevity and the increasing role of xenobiotic degradation mSystems 2020 5(2) e00124 20 32209716
Rampelli, S. et al. Shotgun metagenomics of gut microbiota in humans with up to extreme longevity and the increasing role of xenobiotic degradation. mSystems 5, e00124–20 (2020).32209716
30. Hankin JH Singapore Chinese Health Study: development, validation, and calibration of the quantitative food frequency questionnaire Nutr. Cancer 2001 39 187 195 10.1207/S15327914nc392_5 11759279
Hankin, J. H. et al. Singapore Chinese Health Study: development, validation, and calibration of the quantitative food frequency questionnaire. Nutr. Cancer 39, 187–195 (2001).11759279 10.1207/S15327914nc392_5
31. Chua KY Visceral fat area is the measure of obesity best associated with mobility disability in community dwelling oldest-old Chinese adults BMC Geriatr. 2021 21 282 10.1186/s12877-021-02226-6 33910516
Chua, K. Y. et al. Visceral fat area is the measure of obesity best associated with mobility disability in community dwelling oldest-old Chinese adults. BMC Geriatr. 21, 282 (2021).33910516 10.1186/s12877-021-02226-6
32. Louis P Hold GL Flint HJ The gut microbiota, bacterial metabolites and colorectal cancer Nat. Rev. Microbiol. 2014 12 661 672 10.1038/nrmicro3344 25198138
Louis, P., Hold, G. L. & Flint, H. J. The gut microbiota, bacterial metabolites and colorectal cancer. Nat. Rev. Microbiol. 12, 661–672 (2014).25198138 10.1038/nrmicro3344
33. Canani RB Potential beneficial effects of butyrate in intestinal and extraintestinal diseases World J. Gastroenterol. 2011 17 1519 1528 10.3748/wjg.v17.i12.1519 21472114
Canani, R. B. et al. Potential beneficial effects of butyrate in intestinal and extraintestinal diseases. World J. Gastroenterol. 17, 1519–1528 (2011).21472114 10.3748/wjg.v17.i12.1519
34. Liu H Butyrate: a double-edged sword for health? Adv. Nutr. 2018 9 21 29 10.1093/advances/nmx009 29438462
Liu, H. et al. Butyrate: a double-edged sword for health? Adv. Nutr. 9, 21–29 (2018).29438462 10.1093/advances/nmx009
35. Gounot JS Genome-centric analysis of short and long read metagenomes reveals uncharacterized microbiome diversity in Southeast Asians Nat. Commun. 2022 13 6044 10.1038/s41467-022-33782-z 36229545
Gounot, J. S. et al. Genome-centric analysis of short and long read metagenomes reveals uncharacterized microbiome diversity in Southeast Asians. Nat. Commun. 13, 6044 (2022).36229545 10.1038/s41467-022-33782-z
36. Kang JTL Long-term ecological and evolutionary dynamics in the gut microbiomes of carbapenemase-producing Enterobacteriaceae colonized subjects Nat. Microbiol. 2022 7 1516 1524 10.1038/s41564-022-01221-w 36109646
Kang, J. T. L. et al. Long-term ecological and evolutionary dynamics in the gut microbiomes of carbapenemase-producing Enterobacteriaceae colonized subjects. Nat. Microbiol. 7, 1516–1524 (2022).36109646 10.1038/s41564-022-01221-w
37. Qin J A metagenome-wide association study of gut microbiota in type 2 diabetes Nature 2012 490 55 60 10.1038/nature11450 23023125
Qin, J. et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 490, 55–60 (2012).23023125 10.1038/nature11450
38. Pang S Longevity of centenarians is reflected by the gut microbiome with youth-associated signatures Nat. Aging 2023 3 436 449 10.1038/s43587-023-00389-y 37117794
Pang, S. et al. Longevity of centenarians is reflected by the gut microbiome with youth-associated signatures. Nat. Aging 3, 436–449 (2023).37117794 10.1038/s43587-023-00389-y
39. Rampelli S Functional metagenomic profiling of intestinal microbiome in extreme ageing Aging 2013 5 902 912 10.18632/aging.100623 24334635
Rampelli, S. et al. Functional metagenomic profiling of intestinal microbiome in extreme ageing. Aging 5, 902–912 (2013).24334635 10.18632/aging.100623
40. Yachida S Metagenomic and metabolomic analyses reveal distinct stage-specific phenotypes of the gut microbiota in colorectal cancer Nat. Med. 2019 25 968 976 10.1038/s41591-019-0458-7 31171880
Yachida, S. et al. Metagenomic and metabolomic analyses reveal distinct stage-specific phenotypes of the gut microbiota in colorectal cancer. Nat. Med. 25, 968–976 (2019).31171880 10.1038/s41591-019-0458-7
41. Asnicar F Studying vertical microbiome transmission from mothers to infants by strain-level metagenomic profiling mSystems 2017 2 e00164 16 10.1128/mSystems.00164-16 28144631
Asnicar, F. et al. Studying vertical microbiome transmission from mothers to infants by strain-level metagenomic profiling. mSystems 2, e00164–16 (2017).28144631 10.1128/mSystems.00164-16
42. Brito IL Mobile genes in the human microbiome are structured from global to individual scales Nature 2016 535 435 439 10.1038/nature18927 27409808
Brito, I. L. et al. Mobile genes in the human microbiome are structured from global to individual scales. Nature 535, 435–439 (2016).27409808 10.1038/nature18927
43. Tigchelaar EF Cohort profile: lifeLines DEEP, a prospective, general population cohort study in the northern Netherlands: study design and baseline characteristics BMJ Open 2015 5 e006772 10.1136/bmjopen-2014-006772 26319774
Tigchelaar, E. F. et al. Cohort profile: lifeLines DEEP, a prospective, general population cohort study in the northern Netherlands: study design and baseline characteristics. BMJ Open 5, e006772 (2015).26319774 10.1136/bmjopen-2014-006772
44. Xu Q Metagenomic and metabolomic remodeling in nonagenarians and centenarians and its association with genetic and socioeconomic factors Nat. Aging 2022 2 438 452 10.1038/s43587-022-00193-0 37118062
Xu, Q. et al. Metagenomic and metabolomic remodeling in nonagenarians and centenarians and its association with genetic and socioeconomic factors. Nat. Aging 2, 438–452 (2022).37118062 10.1038/s43587-022-00193-0
45. Biagi E Gut microbiota and extreme longevity Curr. Biol. 2016 26 1480 1485 10.1016/j.cub.2016.04.016 27185560
Biagi, E. et al. Gut microbiota and extreme longevity. Curr. Biol. 26, 1480–1485 (2016).27185560 10.1016/j.cub.2016.04.016
46. Vital M Howe AC Tiedje JM Revealing the bacterial butyrate synthesis pathways by analyzing (meta)genomic data mBio 2014 5 e00889 10.1128/mBio.00889-14 24757212
Vital, M., Howe, A. C. & Tiedje, J. M. Revealing the bacterial butyrate synthesis pathways by analyzing (meta)genomic data. mBio 5, e00889 (2014).24757212 10.1128/mBio.00889-14
47. Vital M Karch A Pieper DH Colonic butyrate-producing communities in humans: an overview using omics data mSystems 2017 2 e00130 17 10.1128/mSystems.00130-17 29238752
Vital, M., Karch, A. & Pieper, D. H. Colonic butyrate-producing communities in humans: an overview using omics data. mSystems 2, e00130–17 (2017).29238752 10.1128/mSystems.00130-17
48. Kaur H Bose C Mande SS Tryptophan metabolism by gut microbiome and gut-brain-axis: an in silico analysis Front. Neurosci. 2019 13 1365 10.3389/fnins.2019.01365 31920519
Kaur, H., Bose, C. & Mande, S. S. Tryptophan metabolism by gut microbiome and gut-brain-axis: an in silico analysis. Front. Neurosci. 13, 1365 (2019).31920519 10.3389/fnins.2019.01365
49. Hawkins CP MacMahon JA Guilds: the multiple meanings of a concept Annu. Rev. Entomol. 1989 34 423 451 10.1146/annurev.en.34.010189.002231
Hawkins, C. P. & MacMahon, J. A. Guilds: the multiple meanings of a concept. Annu. Rev. Entomol. 34, 423–451 (1989).10.1146/annurev.en.34.010189.002231
50. Ma Z Epigenetic drift of H3K27me3 in aging links glycolysis to healthy longevity in Drosophila eLife 2018 7 e35368 10.7554/eLife.35368 29809154
Ma, Z. et al. Epigenetic drift of H3K27me3 in aging links glycolysis to healthy longevity in Drosophila. eLife 7, e35368 (2018).29809154 10.7554/eLife.35368
51. Kennedy DO B vitamins and the brain: mechanisms, dose and efficacy-a review Nutrients 2016 8 68 10.3390/nu8020068 26828517
Kennedy, D. O. B vitamins and the brain: mechanisms, dose and efficacy-a review. Nutrients 8, 68 (2016).26828517 10.3390/nu8020068
52. Suwannasom, N., Kao, I., Pruss, A., Georgieva, R. & Baumler, H. Riboflavin: the health benefits of a forgotten natural vitamin. Int. J. Mol. Sci. 21, 950 (2020).
53. Cheng S Distinct metabolomic signatures are associated with longevity in humans Nat. Commun. 2015 6 6791 10.1038/ncomms7791 25864806
Cheng, S. et al. Distinct metabolomic signatures are associated with longevity in humans. Nat. Commun. 6, 6791 (2015).25864806 10.1038/ncomms7791
54. Asadi Shahmirzadi A Alpha-ketoglutarate, an endogenous metabolite, extends lifespan and compresses morbidity in aging mice Cell Metab. 2020 32 447 456.e446 10.1016/j.cmet.2020.08.004 32877690
Asadi Shahmirzadi, A. et al. Alpha-ketoglutarate, an endogenous metabolite, extends lifespan and compresses morbidity in aging mice. Cell Metab. 32, 447–456.e446 (2020).32877690 10.1016/j.cmet.2020.08.004
55. Prasad S Sung B Aggarwal BB Age-associated chronic diseases require age-old medicine: role of chronic inflammation Prev. Med. 2012 54 S29 S37 10.1016/j.ypmed.2011.11.011 22178471
Prasad, S., Sung, B. & Aggarwal, B. B. Age-associated chronic diseases require age-old medicine: role of chronic inflammation. Prev. Med. 54, S29–S37 (2012).22178471 10.1016/j.ypmed.2011.11.011
56. Kalyani RR Egan JM Diabetes and altered glucose metabolism with aging Endocrinol. Metab. Clin. N. Am. 2013 42 333 347 10.1016/j.ecl.2013.02.010
Kalyani, R. R. & Egan, J. M. Diabetes and altered glucose metabolism with aging. Endocrinol. Metab. Clin. N. Am. 42, 333–347 (2013).10.1016/j.ecl.2013.02.010
57. Chia CW Egan JM Ferrucci L Age-related changes in glucose metabolism, hyperglycemia, and cardiovascular risk Circ. Res. 2018 123 886 904 10.1161/CIRCRESAHA.118.312806 30355075
Chia, C. W., Egan, J. M. & Ferrucci, L. Age-related changes in glucose metabolism, hyperglycemia, and cardiovascular risk. Circ. Res. 123, 886–904 (2018).30355075 10.1161/CIRCRESAHA.118.312806
58. Kirkman MS Diabetes in older adults Diabetes Care 2012 35 2650 2664 10.2337/dc12-1801 23100048
Kirkman, M. S. et al. Diabetes in older adults. Diabetes Care 35, 2650–2664 (2012).23100048 10.2337/dc12-1801
59. Wu TR Gut commensal Parabacteroides goldsteinii plays a predominant role in the anti-obesity effects of polysaccharides isolated from Hirsutella sinensis Gut 2019 68 248 262 10.1136/gutjnl-2017-315458 30007918
Wu, T. R. et al. Gut commensal Parabacteroides goldsteinii plays a predominant role in the anti-obesity effects of polysaccharides isolated from Hirsutella sinensis. Gut 68, 248–262 (2019).30007918 10.1136/gutjnl-2017-315458
60. Ko, Yun-Fei, et al. Method to reduce insulin resistance and improve glucose tolerance using Parabacteroides goldsteinii. Chang Gung Biotechnology Corp. (Taipei, Tw) US-10086019-B1 (2017).
61. Qian X Bifidobacterium adolescentis is effective in relieving type 2 diabetes and may be related to its dominant core genome and gut microbiota modulation capacity Nutrients 2022 14 2479 10.3390/nu14122479 35745208
Qian, X. et al. Bifidobacterium adolescentis is effective in relieving type 2 diabetes and may be related to its dominant core genome and gut microbiota modulation capacity. Nutrients 14, 2479 (2022).35745208 10.3390/nu14122479
62. Hsieh MC The beneficial effects of Lactobacillus reuteri ADR-1 or ADR-3 consumption on type 2 diabetes mellitus: a randomized, double-blinded, placebo-controlled trial Sci. Rep. 2018 8 16791 10.1038/s41598-018-35014-1 30429496
Hsieh, M. C. et al. The beneficial effects of Lactobacillus reuteri ADR-1 or ADR-3 consumption on type 2 diabetes mellitus: a randomized, double-blinded, placebo-controlled trial. Sci. Rep. 8, 16791 (2018).30429496 10.1038/s41598-018-35014-1
63. Li G The effects of probiotics supplementation on glycaemic control among adults with type 2 diabetes mellitus: a systematic review and meta-analysis of randomised clinical trials J. Transl. Med. 2023 21 442 10.1186/s12967-023-04306-0 37415167
Li, G. et al. The effects of probiotics supplementation on glycaemic control among adults with type 2 diabetes mellitus: a systematic review and meta-analysis of randomised clinical trials. J. Transl. Med. 21, 442 (2023).37415167 10.1186/s12967-023-04306-0
64. Razmpoosh E The effect of probiotic supplementation on glycemic control and lipid profile in patients with type 2 diabetes: a randomized placebo controlled trial Diabetes Metab. Syndr. 2019 13 175 182 10.1016/j.dsx.2018.08.008 30641692
Razmpoosh, E. et al. The effect of probiotic supplementation on glycemic control and lipid profile in patients with type 2 diabetes: a randomized placebo controlled trial. Diabetes Metab. Syndr. 13, 175–182 (2019).30641692 10.1016/j.dsx.2018.08.008
65. Moroti C Souza Magri LF de Rezende Costa M Cavallini DC Sivieri K Effect of the consumption of a new symbiotic shake on glycemia and cholesterol levels in elderly people with type 2 diabetes mellitus Lipids Health Dis. 2012 11 29 10.1186/1476-511X-11-29 22356933
Moroti, C., Souza Magri, L. F., de Rezende Costa, M., Cavallini, D. C. & Sivieri, K. Effect of the consumption of a new symbiotic shake on glycemia and cholesterol levels in elderly people with type 2 diabetes mellitus. Lipids Health Dis. 11, 29 (2012).22356933 10.1186/1476-511X-11-29
66. Hantzidiamantis, P. J. & Lappin, S. L. In StatPearls (StatPearls Publishing, Treasure Island, FL, 2023).
67. Felix-Redondo FJ Grau M Fernandez-Berges D Cholesterol and cardiovascular disease in the elderly. Facts and gaps Aging Dis. 2013 4 154 169 23730531
Felix-Redondo, F. J., Grau, M. & Fernandez-Berges, D. Cholesterol and cardiovascular disease in the elderly. Facts and gaps. Aging Dis. 4, 154–169 (2013).23730531
68. Vacca M The controversial role of human gut lachnospiraceae Microorganisms 2020 8 573 10.3390/microorganisms8040573 32326636
Vacca, M. et al. The controversial role of human gut lachnospiraceae. Microorganisms 8, 573 (2020).32326636 10.3390/microorganisms8040573
69. Fu J The gut microbiome contributes to a substantial proportion of the variation in blood lipids Circ. Res. 2015 117 817 824 10.1161/CIRCRESAHA.115.306807 26358192
Fu, J. et al. The gut microbiome contributes to a substantial proportion of the variation in blood lipids. Circ. Res. 117, 817–824 (2015).26358192 10.1161/CIRCRESAHA.115.306807
70. Frankard J Rodriguez-Villalobos H Struelens MJ Jacobs F Haemophilus parainfluenzae: an underdiagnosed pathogen of biliary tract infections? Eur. J. Clin. Microbiol. Infect. Dis. 2004 23 46 48 10.1007/s10096-003-1050-z 14669072
Frankard, J., Rodriguez-Villalobos, H., Struelens, M. J. & Jacobs, F. Haemophilus parainfluenzae: an underdiagnosed pathogen of biliary tract infections? Eur. J. Clin. Microbiol. Infect. Dis. 23, 46–48 (2004).14669072 10.1007/s10096-003-1050-z
71. Ren X Bacterial alterations in post-cholecystectomy patients are associated with colorectal cancer Front. Oncol. 2020 10 1418 10.3389/fonc.2020.01418 32903396
Ren, X. et al. Bacterial alterations in post-cholecystectomy patients are associated with colorectal cancer. Front. Oncol. 10, 1418 (2020).32903396 10.3389/fonc.2020.01418
72. Latyshev Y Mathew A Jacobson JM Sturm E Purulent pericarditis caused by Haemophilus parainfluenzae Tex. Heart Inst. J. 2013 40 608 611 24391338
Latyshev, Y., Mathew, A., Jacobson, J. M. & Sturm, E. Purulent pericarditis caused by Haemophilus parainfluenzae. Tex. Heart Inst. J. 40, 608–611 (2013).24391338
73. Kelly TN Gut microbiome associates with lifetime cardiovascular disease risk profile among bogalusa heart study participants Circ. Res. 2016 119 956 964 10.1161/CIRCRESAHA.116.309219 27507222
Kelly, T. N. et al. Gut microbiome associates with lifetime cardiovascular disease risk profile among bogalusa heart study participants. Circ. Res. 119, 956–964 (2016).27507222 10.1161/CIRCRESAHA.116.309219
74. Vasto S Inflammatory networks in ageing, age-related diseases and longevity Mech. Ageing Dev. 2007 128 83 91 10.1016/j.mad.2006.11.015 17118425
Vasto, S. et al. Inflammatory networks in ageing, age-related diseases and longevity. Mech. Ageing Dev. 128, 83–91 (2007).17118425 10.1016/j.mad.2006.11.015
75. Leonardi GC Accardi G Monastero R Nicoletti F Libra M Ageing: from inflammation to cancer Immun. Ageing 2018 15 1 10.1186/s12979-017-0112-5 29387133
Leonardi, G. C., Accardi, G., Monastero, R., Nicoletti, F. & Libra, M. Ageing: from inflammation to cancer. Immun. Ageing 15, 1 (2018).29387133 10.1186/s12979-017-0112-5
76. Franceschi C Inflammaging and anti-inflammaging: a systemic perspective on aging and longevity emerged from studies in humans Mech. Ageing Dev. 2007 128 92 105 10.1016/j.mad.2006.11.016 17116321
Franceschi, C. et al. Inflammaging and anti-inflammaging: a systemic perspective on aging and longevity emerged from studies in humans. Mech. Ageing Dev. 128, 92–105 (2007).17116321 10.1016/j.mad.2006.11.016
77. Liberale L Montecucco F Tardif JC Libby P Camici GG Inflamm-ageing: the role of inflammation in age-dependent cardiovascular disease Eur. Heart J. 2020 41 2974 2982 10.1093/eurheartj/ehz961 32006431
Liberale, L., Montecucco, F., Tardif, J. C., Libby, P. & Camici, G. G. Inflamm-ageing: the role of inflammation in age-dependent cardiovascular disease. Eur. Heart J. 41, 2974–2982 (2020).32006431 10.1093/eurheartj/ehz961
78. Ferrucci L Fabbri E Inflammageing: chronic inflammation in ageing, cardiovascular disease, and frailty Nat. Rev. Cardiol. 2018 15 505 522 10.1038/s41569-018-0064-2 30065258
Ferrucci, L. & Fabbri, E. Inflammageing: chronic inflammation in ageing, cardiovascular disease, and frailty. Nat. Rev. Cardiol. 15, 505–522 (2018).30065258 10.1038/s41569-018-0064-2
79. Fung C-P A global emerging disease of Klebsiella pneumoniae liver abscess: is serotype K1 an important factor for complicated endophthalmitis? Gut 2002 50 420 424 10.1136/gut.50.3.420 11839725
Fung, C.-P. et al. A global emerging disease of Klebsiella pneumoniae liver abscess: is serotype K1 an important factor for complicated endophthalmitis? Gut 50, 420–424 (2002).11839725 10.1136/gut.50.3.420
80. Tang Y Fung E Xu A Lan HY C-reactive protein and ageing Clin. Exp. Pharm. Physiol. 2017 44 9 14 10.1111/1440-1681.12758
Tang, Y., Fung, E., Xu, A. & Lan, H. Y. C-reactive protein and ageing. Clin. Exp. Pharm. Physiol. 44, 9–14 (2017).10.1111/1440-1681.12758
81. Rajkumar H Effect of probiotic (VSL#3) and omega-3 on lipid profile, insulin sensitivity, inflammatory markers, and gut colonization in overweight adults: a randomized, controlled trial Mediat. Inflamm. 2014 2014 348959 10.1155/2014/348959
Rajkumar, H. et al. Effect of probiotic (VSL#3) and omega-3 on lipid profile, insulin sensitivity, inflammatory markers, and gut colonization in overweight adults: a randomized, controlled trial. Mediat. Inflamm. 2014, 348959 (2014).10.1155/2014/348959
82. Zhang X The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment Nat. Med. 2015 21 895 905 10.1038/nm.3914 26214836
Zhang, X. et al. The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nat. Med. 21, 895–905 (2015).26214836 10.1038/nm.3914
83. Alexander M Human gut bacterial metabolism drives Th17 activation and colitis Cell Host Microbe 2022 30 17 30.e19 10.1016/j.chom.2021.11.001 34822777
Alexander, M. et al. Human gut bacterial metabolism drives Th17 activation and colitis. Cell Host Microbe 30, 17–30.e19 (2022).34822777 10.1016/j.chom.2021.11.001
84. Hou Y Ageing as a risk factor for neurodegenerative disease Nat. Rev. Neurol. 2019 15 565 581 10.1038/s41582-019-0244-7 31501588
Hou, Y. et al. Ageing as a risk factor for neurodegenerative disease. Nat. Rev. Neurol. 15, 565–581 (2019).31501588 10.1038/s41582-019-0244-7
85. Adamo AM Nutritional factors and aging in demyelinating diseases Genes Nutr. 2014 9 360 10.1007/s12263-013-0360-8 24311441
Adamo, A. M. Nutritional factors and aging in demyelinating diseases. Genes Nutr. 9, 360 (2014).24311441 10.1007/s12263-013-0360-8
86. Healton EB Savage DG Brust JC Garrett TJ Lindenbaum J Neurologic aspects of cobalamin deficiency Medicine 1991 70 229 245 10.1097/00005792-199107000-00001 1648656
Healton, E. B., Savage, D. G., Brust, J. C., Garrett, T. J. & Lindenbaum, J. Neurologic aspects of cobalamin deficiency. Medicine 70, 229–245 (1991).1648656 10.1097/00005792-199107000-00001
87. Hurvitz JKRAEA Peripheral neuropathy: a true risk factor for falls J. Gerontol. Med. Sci. 1995 50A M211 M215 10.1093/gerona/50A.4.M211
Hurvitz, J. K. R. A. E. A. Peripheral neuropathy: a true risk factor for falls. J. Gerontol. Med. Sci. 50A, M211–M215 (1995).10.1093/gerona/50A.4.M211
88. Rodionov DA Micronutrient requirements and sharing capabilities of the human gut microbiome Front. Microbiol. 2019 10 1316 10.3389/fmicb.2019.01316 31275260
Rodionov, D. A. et al. Micronutrient requirements and sharing capabilities of the human gut microbiome. Front. Microbiol. 10, 1316 (2019).31275260 10.3389/fmicb.2019.01316
89. Magnusdottir S Ravcheev D de Crecy-Lagard V Thiele I Systematic genome assessment of B-vitamin biosynthesis suggests co-operation among gut microbes Front. Genet. 2015 6 148 10.3389/fgene.2015.00148 25941533
Magnusdottir, S., Ravcheev, D., de Crecy-Lagard, V. & Thiele, I. Systematic genome assessment of B-vitamin biosynthesis suggests co-operation among gut microbes. Front. Genet. 6, 148 (2015).25941533 10.3389/fgene.2015.00148
90. Wexler AG Human gut Bacteroides capture vitamin B12 via cell surface-exposed lipoproteins eLife 2018 7 e37138 10.7554/eLife.37138 30226189
Wexler, A. G. et al. Human gut Bacteroides capture vitamin B12 via cell surface-exposed lipoproteins. eLife 7, e37138 (2018).30226189 10.7554/eLife.37138
91. Debelius J Tiny microbes, enormous impacts: what matters in gut microbiome studies? Genome Biol. 2016 17 217 10.1186/s13059-016-1086-x 27760558
Debelius, J. et al. Tiny microbes, enormous impacts: what matters in gut microbiome studies? Genome Biol. 17, 217 (2016).27760558 10.1186/s13059-016-1086-x
92. Arumugam M Enterotypes of the human gut microbiome Nature 2011 473 174 180 10.1038/nature09944 21508958
Arumugam, M. et al. Enterotypes of the human gut microbiome. Nature 473, 174–180 (2011).21508958 10.1038/nature09944
93. Xu C Zhu H Qiu P Aging progression of human gut microbiota BMC Microbiol. 2019 19 236 10.1186/s12866-019-1616-2 31660868
Xu, C., Zhu, H. & Qiu, P. Aging progression of human gut microbiota. BMC Microbiol. 19, 236 (2019).31660868 10.1186/s12866-019-1616-2
94. Kong F Deng F Li Y Zhao J Identification of gut microbiome signatures associated with longevity provides a promising modulation target for healthy aging Gut Microbes 2019 10 210 215 10.1080/19490976.2018.1494102 30142010
Kong, F., Deng, F., Li, Y. & Zhao, J. Identification of gut microbiome signatures associated with longevity provides a promising modulation target for healthy aging. Gut Microbes 10, 210–215 (2019).30142010 10.1080/19490976.2018.1494102
95. Litichevskiy, L. et al. Interactions between the gut microbiome, dietary restriction, and aging in genetically diverse mice. Preprint at bioRxiv; 10.1101/2023.11.28.568137 (2023).
96. Jung HW Frailty as a clinically relevant measure of human aging Ann. Geriatr. Med. Res. 2021 25 139 140 10.4235/agmr.21.0106 34610664
Jung, H. W. Frailty as a clinically relevant measure of human aging. Ann. Geriatr. Med. Res. 25, 139–140 (2021).34610664 10.4235/agmr.21.0106
97. Lara J A proposed panel of biomarkers of healthy ageing BMC Med. 2015 13 222 10.1186/s12916-015-0470-9 26373927
Lara, J. et al. A proposed panel of biomarkers of healthy ageing. BMC Med. 13, 222 (2015).26373927 10.1186/s12916-015-0470-9
98. Howlett SE Rutenberg AD Rockwood K The degree of frailty as a translational measure of health in aging Nat. Aging 2021 1 651 665 10.1038/s43587-021-00099-3 37117769
Howlett, S. E., Rutenberg, A. D. & Rockwood, K. The degree of frailty as a translational measure of health in aging. Nat. Aging 1, 651–665 (2021).37117769 10.1038/s43587-021-00099-3
99. Falahati F The effect of age correction on multivariate classification in Alzheimer’s disease, with a focus on the characteristics of incorrectly and correctly classified subjects Brain Topogr. 2016 29 296 307 10.1007/s10548-015-0455-1 26440606
Falahati, F. et al. The effect of age correction on multivariate classification in Alzheimer’s disease, with a focus on the characteristics of incorrectly and correctly classified subjects. Brain Topogr. 29, 296–307 (2016).26440606 10.1007/s10548-015-0455-1
100. David LA Diet rapidly and reproducibly alters the human gut microbiome Nature 2014 505 559 563 10.1038/nature12820 24336217
David, L. A. et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature 505, 559–563 (2014).24336217 10.1038/nature12820
101. Salazar N Valdes-Varela L Gonzalez S Gueimonde M de Los Reyes-Gavilan CG Nutrition and the gut microbiome in the elderly Gut Microbes 2017 8 82 97 10.1080/19490976.2016.1256525 27808595
Salazar, N., Valdes-Varela, L., Gonzalez, S., Gueimonde, M. & de Los Reyes-Gavilan, C. G. Nutrition and the gut microbiome in the elderly. Gut Microbes 8, 82–97 (2017).27808595 10.1080/19490976.2016.1256525
102. Woodmansey EJ McMurdo ME Macfarlane GT Macfarlane S Comparison of compositions and metabolic activities of fecal microbiotas in young adults and in antibiotic-treated and non-antibiotic-treated elderly subjects Appl. Environ. Microbiol. 2004 70 6113 6122 10.1128/AEM.70.10.6113-6122.2004 15466557
Woodmansey, E. J., McMurdo, M. E., Macfarlane, G. T. & Macfarlane, S. Comparison of compositions and metabolic activities of fecal microbiotas in young adults and in antibiotic-treated and non-antibiotic-treated elderly subjects. Appl. Environ. Microbiol. 70, 6113–6122 (2004).15466557 10.1128/AEM.70.10.6113-6122.2004
103. Dabek M alpha-Ketoglutarate (AKG) absorption from pig intestine and plasma pharmacokinetics J. Anim. Physiol. Anim. Nutr. 2005 89 419 426 10.1111/j.1439-0396.2005.00566.x
Dabek, M. et al. alpha-Ketoglutarate (AKG) absorption from pig intestine and plasma pharmacokinetics. J. Anim. Physiol. Anim. Nutr. 89, 419–426 (2005).10.1111/j.1439-0396.2005.00566.x
104. Naeini SH Mavaddatiyan L Kalkhoran ZR Taherkhani S Talkhabi M Alpha-ketoglutarate as a potent regulator for lifespan and healthspan: evidences and perspectives Exp. Gerontol. 2023 175 112154 10.1016/j.exger.2023.112154 36934991
Naeini, S. H., Mavaddatiyan, L., Kalkhoran, Z. R., Taherkhani, S. & Talkhabi, M. Alpha-ketoglutarate as a potent regulator for lifespan and healthspan: evidences and perspectives. Exp. Gerontol. 175, 112154 (2023).36934991 10.1016/j.exger.2023.112154
105. Parada Venegas D Short Chain Fatty Acids (SCFAs)-mediated gut epithelial and immune regulation and its relevance for inflammatory bowel diseases Front. Immunol. 2019 10 277 10.3389/fimmu.2019.00277 30915065
Parada Venegas, D. et al. Short Chain Fatty Acids (SCFAs)-mediated gut epithelial and immune regulation and its relevance for inflammatory bowel diseases. Front. Immunol. 10, 277 (2019).30915065 10.3389/fimmu.2019.00277
106. Dalile B Van Oudenhove L Vervliet B Verbeke K The role of short-chain fatty acids in microbiota-gut-brain communication Nat. Rev. Gastroenterol. Hepatol. 2019 16 461 478 10.1038/s41575-019-0157-3 31123355
Dalile, B., Van Oudenhove, L., Vervliet, B. & Verbeke, K. The role of short-chain fatty acids in microbiota-gut-brain communication. Nat. Rev. Gastroenterol. Hepatol. 16, 461–478 (2019).31123355 10.1038/s41575-019-0157-3
107. Jackson MA Signatures of early frailty in the gut microbiota Genome Med. 2016 8 8 10.1186/s13073-016-0262-7 26822992
Jackson, M. A. et al. Signatures of early frailty in the gut microbiota. Genome Med. 8, 8 (2016).26822992 10.1186/s13073-016-0262-7
108. Nagao K Kimura T Use of plasma-free amino acids as biomarkers for detecting and predicting disease risk Nutr. Rev. 2020 78 79 85 10.1093/nutrit/nuaa086 33259618
Nagao, K. & Kimura, T. Use of plasma-free amino acids as biomarkers for detecting and predicting disease risk. Nutr. Rev. 78, 79–85 (2020).33259618 10.1093/nutrit/nuaa086
109. Ganesh BP Prebiotics, probiotics, and acetate supplementation prevent hypertension in a model of obstructive sleep apnea Hypertension 2018 72 1141 1150 10.1161/HYPERTENSIONAHA.118.11695 30354816
Ganesh, B. P. et al. Prebiotics, probiotics, and acetate supplementation prevent hypertension in a model of obstructive sleep apnea. Hypertension 72, 1141–1150 (2018).30354816 10.1161/HYPERTENSIONAHA.118.11695
110. Richards JL Yap YA McLeod KH Mackay CR Marino E Dietary metabolites and the gut microbiota: an alternative approach to control inflammatory and autoimmune diseases Clin. Transl. Immunol. 2016 5 e82 10.1038/cti.2016.29
Richards, J. L., Yap, Y. A., McLeod, K. H., Mackay, C. R. & Marino, E. Dietary metabolites and the gut microbiota: an alternative approach to control inflammatory and autoimmune diseases. Clin. Transl. Immunol. 5, e82 (2016).10.1038/cti.2016.29
111. Kazemian N Mahmoudi M Halperin F Wu JC Pakpour S Gut microbiota and cardiovascular disease: opportunities and challenges Microbiome 2020 8 36 10.1186/s40168-020-00821-0 32169105
Kazemian, N., Mahmoudi, M., Halperin, F., Wu, J. C. & Pakpour, S. Gut microbiota and cardiovascular disease: opportunities and challenges. Microbiome 8, 36 (2020).32169105 10.1186/s40168-020-00821-0
112. Cryan JF O’Riordan KJ Sandhu K Peterson V Dinan TG The gut microbiome in neurological disorders Lancet Neurol. 2020 19 179 194 10.1016/S1474-4422(19)30356-4 31753762
Cryan, J. F., O’Riordan, K. J., Sandhu, K., Peterson, V. & Dinan, T. G. The gut microbiome in neurological disorders. Lancet Neurol. 19, 179–194 (2020).31753762 10.1016/S1474-4422(19)30356-4
113. Ruan Y Effect of probiotics on glycemic control: a systematic review and meta-analysis of randomized, controlled trials PLoS ONE 2015 10 e0132121 10.1371/journal.pone.0132121 26161741
Ruan, Y. et al. Effect of probiotics on glycemic control: a systematic review and meta-analysis of randomized, controlled trials. PLoS ONE 10, e0132121 (2015).26161741 10.1371/journal.pone.0132121
114. Logan IE Xanthohumol requires the intestinal microbiota to improve glucose metabolism in diet-induced obese mice Mol. Nutr. Food Res. 2021 65 e2100389 10.1002/mnfr.202100389 34496124
Logan, I. E. et al. Xanthohumol requires the intestinal microbiota to improve glucose metabolism in diet-induced obese mice. Mol. Nutr. Food Res. 65, e2100389 (2021).34496124 10.1002/mnfr.202100389
115. Salles N Herrmann F Sakbani K Rapin CH Sieber C High vitamin B12 level: a strong predictor of mortality in elderly inpatients J. Am. Geriatr. Soc. 2005 53 917 918 10.1111/j.1532-5415.2005.53278_7.x 15877584
Salles, N., Herrmann, F., Sakbani, K., Rapin, C. H. & Sieber, C. High vitamin B12 level: a strong predictor of mortality in elderly inpatients. J. Am. Geriatr. Soc. 53, 917–918 (2005).15877584 10.1111/j.1532-5415.2005.53278_7.x
116. Ng TP Socio-environmental, lifestyle, behavioural, and psychological determinants of biological ageing: the singapore longitudinal ageing study Gerontology 2020 66 603 613 10.1159/000511211 33197920
Ng, T. P. et al. Socio-environmental, lifestyle, behavioural, and psychological determinants of biological ageing: the singapore longitudinal ageing study. Gerontology 66, 603–613 (2020).33197920 10.1159/000511211
117. Saw WY Establishing multiple omics baselines for three Southeast Asian populations in the Singapore Integrative Omics Study Nat. Commun. 2017 8 653 10.1038/s41467-017-00413-x 28935855
Saw, W. Y. et al. Establishing multiple omics baselines for three Southeast Asian populations in the Singapore Integrative Omics Study. Nat. Commun. 8, 653 (2017).28935855 10.1038/s41467-017-00413-x
118. Pasolli E Accessible, curated metagenomic data through ExperimentHub Nat. Methods 2017 14 1023 1024 10.1038/nmeth.4468 29088129
Pasolli, E. et al. Accessible, curated metagenomic data through ExperimentHub. Nat. Methods 14, 1023–1024 (2017).29088129 10.1038/nmeth.4468
119. Chen S Zhou Y Chen Y Gu J fastp: an ultra-fast all-in-one FASTQ preprocessor Bioinformatics 2018 34 i884 i890 10.1093/bioinformatics/bty560 30423086
Chen, S., Zhou, Y., Chen, Y. & Gu, J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890 (2018).30423086 10.1093/bioinformatics/bty560
120. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at https://arxiv.org/abs/1303.3997 (2013).
121. Li H The sequence alignment/map format and SAMtools Bioinformatics 2009 25 2078 2079 10.1093/bioinformatics/btp352 19505943
Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).19505943 10.1093/bioinformatics/btp352
122. Segata N Metagenomic microbial community profiling using unique clade-specific marker genes Nat. Methods 2012 9 811 814 10.1038/nmeth.2066 22688413
Segata, N. et al. Metagenomic microbial community profiling using unique clade-specific marker genes. Nat. Methods 9, 811–814 (2012).22688413 10.1038/nmeth.2066
123. Franzosa EA Species-level functional profiling of metagenomes and metatranscriptomes Nat. Methods 2018 15 962 968 10.1038/s41592-018-0176-y 30377376
Franzosa, E. A. et al. Species-level functional profiling of metagenomes and metatranscriptomes. Nat. Methods 15, 962–968 (2018).30377376 10.1038/s41592-018-0176-y
124. Ma, S. et al. Population structure discovery in meta-analyzed microbial communities and inflammatory bowel disease using MMUPHin. Genome Biol. 23, 208 (2022).
125. Wang Y LeCao KA Managing batch effects in microbiome data Brief. Bioinform. 2020 21 1954 1970 10.1093/bib/bbz105 31776547
Wang, Y. & LeCao, K. A. Managing batch effects in microbiome data. Brief. Bioinform. 21, 1954–1970 (2020).31776547 10.1093/bib/bbz105
126. Mallick H Multivariable association discovery in population-scale meta-omics studies PLoS Comput Biol. 2021 17 e1009442 10.1371/journal.pcbi.1009442 34784344
Mallick, H. et al. Multivariable association discovery in population-scale meta-omics studies. PLoS Comput Biol. 17, e1009442 (2021).34784344 10.1371/journal.pcbi.1009442
127. Dixon P VEGAN, a package of R functions for community ecology J. Veg. Sci. 2003 14 927 930 10.1111/j.1654-1103.2003.tb02228.x
Dixon, P. VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930 (2003).10.1111/j.1654-1103.2003.tb02228.x
128. Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag New York, 2016).
129. Jastroch M Ussar S Keipert S Gut microbes controlling blood sugar: no fire required Cell Metab. 2020 31 443 444 10.1016/j.cmet.2020.02.007 32130877
Jastroch, M., Ussar, S. & Keipert, S. Gut microbes controlling blood sugar: no fire required. Cell Metab. 31, 443–444 (2020).32130877 10.1016/j.cmet.2020.02.007
130. Little, R. J. A. & Rubin, D. B. Statistical Analysis with Missing Data (Wiley, 2002).
131. Margolin, B. H. in Encyclopedia of Statistical Sciences (eds S. Kotz, Johnson, N. L., & Read, C. B.) (John Wiley & Sons, New York; 1988).
132. Lin H Peddada SD Analysis of compositions of microbiomes with bias correction Nat. Commun. 2020 11 3514 10.1038/s41467-020-17041-7 32665548
Lin, H. & Peddada, S. D. Analysis of compositions of microbiomes with bias correction. Nat. Commun. 11, 3514 (2020).32665548 10.1038/s41467-020-17041-7
133. Faust K Microbial co-occurrence relationships in the human microbiome PLoS Comput. Biol. 2012 8 e1002606 10.1371/journal.pcbi.1002606 22807668
Faust, K. et al. Microbial co-occurrence relationships in the human microbiome. PLoS Comput. Biol. 8, e1002606 (2012).22807668 10.1371/journal.pcbi.1002606
134. Abubucker S Metabolic reconstruction for metagenomic data and its application to the human microbiome PLoS Comput. Biol. 2012 8 e1002358 10.1371/journal.pcbi.1002358 22719234
Abubucker, S. et al. Metabolic reconstruction for metagenomic data and its application to the human microbiome. PLoS Comput. Biol. 8, e1002358 (2012).22719234 10.1371/journal.pcbi.1002358
135. Kanehisa M Goto S KEGG: kyoto encyclopedia of genes and genomes Nucleic Acids Res. 2000 28 27 30 10.1093/nar/28.1.27 10592173
Kanehisa, M. & Goto, S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30 (2000).10592173 10.1093/nar/28.1.27
136. Segata N Metagenomic biomarker discovery and explanation Genome Biol. 2011 12 R60 10.1186/gb-2011-12-6-r60 21702898
Segata, N. et al. Metagenomic biomarker discovery and explanation. Genome Biol. 12, R60 (2011).21702898 10.1186/gb-2011-12-6-r60
137. Suzek BE UniRef clusters: a comprehensive and scalable alternative for improving sequence similarity searches Bioinformatics 2015 31 926 932 10.1093/bioinformatics/btu739 25398609
Suzek, B. E. et al. UniRef clusters: a comprehensive and scalable alternative for improving sequence similarity searches. Bioinformatics 31, 926–932 (2015).25398609 10.1093/bioinformatics/btu739
138. Li D MEGAHIT v1.0: a fast and scalable metagenome assembler driven by advanced methodologies and community practices Methods 2016 102 3 11 10.1016/j.ymeth.2016.02.020 27012178
Li, D. et al. MEGAHIT v1.0: a fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods 102, 3–11 (2016).27012178 10.1016/j.ymeth.2016.02.020
139. Kang DD MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies PeerJ 2019 7 e7359 10.7717/peerj.7359 31388474
Kang, D. D. et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ 7, e7359 (2019).31388474 10.7717/peerj.7359
140. Li H Minimap2: pairwise alignment for nucleotide sequences Bioinformatics 2018 34 3094 3100 10.1093/bioinformatics/bty191 29750242
Li, H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094–3100 (2018).29750242 10.1093/bioinformatics/bty191
141. Bowers RM Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea Nat. Biotechnol. 2017 35 725 731 10.1038/nbt.3893 28787424
Bowers, R. M. et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat. Biotechnol. 35, 725–731 (2017).28787424 10.1038/nbt.3893
142. Chklovski A Parks DH Woodcroft BJ Tyson GW CheckM2: a rapid, scalable and accurate tool for assessing microbial genome quality using machine learning Nat. Methods 2023 20 1203 1212 10.1038/s41592-023-01940-w 37500759
Chklovski, A., Parks, D. H., Woodcroft, B. J. & Tyson, G. W. CheckM2: a rapid, scalable and accurate tool for assessing microbial genome quality using machine learning. Nat. Methods 20, 1203–1212 (2023).37500759 10.1038/s41592-023-01940-w
143. Lowe TM Chan PP tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes Nucleic Acids Res. 2016 44 W54 W57 10.1093/nar/gkw413 27174935
Lowe, T. M. & Chan, P. P. tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 44, W54–W57 (2016).27174935 10.1093/nar/gkw413
144. Orakov A GUNC: detection of chimerism and contamination in prokaryotic genomes Genome Biol. 2021 22 178 10.1186/s13059-021-02393-0 34120611
Orakov, A. et al. GUNC: detection of chimerism and contamination in prokaryotic genomes. Genome Biol. 22, 178 (2021).34120611 10.1186/s13059-021-02393-0
145. Chaumeil PA Mussig AJ Hugenholtz P Parks DH GTDB-Tk v2: memory friendly classification with the genome taxonomy database Bioinformatics 2022 38 5315 5316 10.1093/bioinformatics/btac672 36218463
Chaumeil, P. A., Mussig, A. J., Hugenholtz, P. & Parks, D. H. GTDB-Tk v2: memory friendly classification with the genome taxonomy database. Bioinformatics 38, 5315–5316 (2022).36218463 10.1093/bioinformatics/btac672
146. Seemann T Prokka: rapid prokaryotic genome annotation Bioinformatics 2014 30 2068 2069 10.1093/bioinformatics/btu153 24642063
Seemann, T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014).24642063 10.1093/bioinformatics/btu153
147. Ye Y Doak TG A parsimony approach to biological pathway reconstruction/inference for genomes and metagenomes PLoS Comput. Biol. 2009 5 e1000465 10.1371/journal.pcbi.1000465 19680427
Ye, Y. & Doak, T. G. A parsimony approach to biological pathway reconstruction/inference for genomes and metagenomes. PLoS Comput. Biol. 5, e1000465 (2009).19680427 10.1371/journal.pcbi.1000465
148. Ravikrishnan A Nasre M Raman K Enumerating all possible biosynthetic pathways in metabolic networks Sci. Rep. 2018 8 9932 10.1038/s41598-018-28007-7 29967471
Ravikrishnan, A., Nasre, M. & Raman, K. Enumerating all possible biosynthetic pathways in metabolic networks. Sci. Rep. 8, 9932 (2018).29967471 10.1038/s41598-018-28007-7
149. Cantalapiedra CP Hernandez-Plaza A Letunic I Bork P Huerta-Cepas J eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale Mol. Biol. Evol. 2021 38 5825 5829 10.1093/molbev/msab293 34597405
Cantalapiedra, C. P., Hernandez-Plaza, A., Letunic, I., Bork, P. & Huerta-Cepas, J. eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol. Biol. Evol. 38, 5825–5829 (2021).34597405 10.1093/molbev/msab293
150. Huerta-Cepas J eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses Nucleic Acids Res. 2019 47 D309 D314 10.1093/nar/gky1085 30418610
Huerta-Cepas, J. et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314 (2019).30418610 10.1093/nar/gky1085
151. Buchfink B Reuter K Drost HG Sensitive protein alignments at tree-of-life scale using DIAMOND Nat. Methods 2021 18 366 368 10.1038/s41592-021-01101-x 33828273
Buchfink, B., Reuter, K. & Drost, H. G. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat. Methods 18, 366–368 (2021).33828273 10.1038/s41592-021-01101-x
152. Ravikrishnan, A. & Raman, K. Unraveling microbial interactions in the gut microbiome. Preprint at bioRxiv https://www.biorxiv.org/content/10.1101/2021.05.17.444446v2.full (2021).
153. Shannon P Cytoscape: a software environment for integrated models of biomolecular interaction networks Genome Res. 2003 13 2498 2504 10.1101/gr.1239303 14597658
Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504 (2003).14597658 10.1101/gr.1239303
154. Tuikhar N Comparative analysis of the gut microbiota in centenarians and young adults shows a common signature across genotypically non-related populations Mech. Ageing Dev. 2019 179 23 35 10.1016/j.mad.2019.02.001 30738080
Tuikhar, N. et al. Comparative analysis of the gut microbiota in centenarians and young adults shows a common signature across genotypically non-related populations. Mech. Ageing Dev. 179, 23–35 (2019).30738080 10.1016/j.mad.2019.02.001
