Abstract
Insight into associations between the gut microbiome with metabolism and aging is crucial for tailoring interventions to promote healthy longevity. In a discovery cohort of 10,207 individuals aged 40–93 years, we used 21 metabolic parameters to classify individuals into five clusters, termed metabolic multimorbidity clusters (MCs), that represent different metabolic subphenotypes. Compared to the cluster classified as metabolically healthy (MC1), clusters classified as ‘obesity-related mixed’ (MC4) and ‘hyperglycemia’ (MC5) exhibited an increased 11.1-year cardiovascular disease (CVD) risk by 75% (multivariable-adjusted hazard ratio (HR): 1.75, 95% confidence interval (CI): 1.43–2.14) and by 117% (2.17, 1.72–2.74), respectively. These associations were replicated in a second cohort of 9,061 individuals with a 10.0-year follow-up. Based on analysis of 4,491 shotgun fecal metagenomes from the discovery cohort, we found that gut microbial composition was associated with both MCs and age. Next, using 55 age-specific microbial species to capture biological age, we developed a gut microbial age (MA) metric, which was validated in four external cohorts comprising 4,425 metagenomic samples. Among individuals aged 60 years or older, the increased CVD risk associated with MC4 or MC5, as compared to MC1, MC2 or MC3, was exacerbated in individuals with high MA but diminished in individuals with low MA, independent of age, sex and other lifestyle and dietary factors. This pattern, in which younger MA appears to counteract the CVD risk attributable to metabolic dysfunction, implies a modulating role of MA in cardiovascular health for metabolically unhealthy older people.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$209.00 per year
only $17.42 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
In accordance with the Medical Ethics Committee of Ruijin Hospital and the institutional review board of BGI-Shenzhen related to protecting individual privacy, metagenomic sequencing data of the 4,491 fecal DNA samples from the JD_2014 cohort have been deposited at the China Nucleotide Sequence Archive under accession number CNP0004479 (https://doi.org/10.26036/CNP0004479) for controlled access. The controlled-access sequencing data are available upon reasonable request from the corresponding Data Access Committee (tiange.wang@shsmu.edu.cn). All phenotypic data that support the findings of this study are available from the corresponding author upon reasonable request.
Code availability
Computational analyses for microbiota were conducted using the bioBakery suite of tools. Species-level microbial abundances were determined using MetaPhlAn version 3.0.7 (https://github.com/biobakery/MetaPhlAn), and functional potential profiling was performed with HUMAnN version 3.0.0 (https://github.com/biobakery/humann).
References
Palmer, A. K. & Jensen, M. D. Metabolic changes in aging humans: current evidence and therapeutic strategies. J. Clin. Invest. 132, e158451 (2022).
Vaduganathan, M., Mensah, G. A., Turco, J. V., Fuster, V. & Roth, G. A. The global burden of cardiovascular diseases and risk: a compass for future health. J. Am. Coll. Cardiol. 80, 2361–2371 (2022).
Forman, D. E. et al. Multimorbidity in older adults with cardiovascular disease. J. Am. Coll. Cardiol. 71, 2149–2161 (2018).
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).
Wilmanski, T. et al. Gut microbiome pattern reflects healthy ageing and predicts survival in humans. Nat. Metab. 3, 274–286 (2021).
Ghosh, T. S., Shanahan, F. & O’Toole, P. W. Toward an improved definition of a healthy microbiome for healthy aging. Nat. Aging 2, 1054–1069 (2022).
Wilmanski, T., Gibbons, S. M. & Price, N. D. Healthy aging and the human gut microbiome: why we cannot just turn back the clock. Nat. Aging 2, 869–871 (2022).
Ghosh, T. S., Shanahan, F. & O’Toole, P. W. The gut microbiome as a modulator of healthy ageing. Nat. Rev. Gastroenterol. Hepatol. 19, 565–584 (2022).
Biagi, E. et al. Gut microbiota and extreme longevity. Curr. Biol. 26, 1480–1485 (2016).
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).
Li, C. et al. Deep insights into the gut microbial community of extreme longevity in south Chinese centenarians by ultra-deep metagenomics and large-scale culturomics. NPJ Biofilms Microbiomes 8, 28 (2022).
Pang, S. et al. Longevity of centenarians is reflected by the gut microbiome with youth-associated signatures. Nat. Aging 3, 436–449 (2023).
Fromentin, S. et al. Microbiome and metabolome features of the cardiometabolic disease spectrum. Nat. Med. 28, 303–314 (2022).
Wagner, R. et al. Pathophysiology-based subphenotyping of individuals at elevated risk for type 2 diabetes. Nat. Med. 27, 49–57 (2021).
Fu, J. et al. Ageing trajectory of the gut microbiota is associated with metabolic diseases in a chronological age-dependent manner. Gut 72, 1431–1433 (2023).
Salosensaari, A. et al. Taxonomic signatures of cause-specific mortality risk in human gut microbiome. Nat. Commun. 12, 2671 (2021).
Ruuskanen, M. O. et al. Gut microbiome composition is predictive of incident type 2 diabetes in a population cohort of 5,572 Finnish adults. Diabetes Care 45, 811–818 (2022).
Zeevi, D. et al. Personalized nutrition by prediction of glycemic responses. Cell 163, 1079–1094 (2015).
Zhernakova, A. et al. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science 352, 565–569 (2016).
Forslund, S. K. et al. Combinatorial, additive and dose-dependent drug–microbiome associations. Nature 600, 500–505 (2021).
Asnicar, F. et al. Microbiome connections with host metabolism and habitual diet from 1,098 deeply phenotyped individuals. Nat. Med. 27, 321–332 (2021).
Davies, D. L. & Bouldin, D. W. A cluster separation measure. IEEE Trans. Pattern Anal. Mach. Intell. 1, 224–227 (1979).
Krzanowski, W. J. & Lai, Y. T. A criterion for determining the number of groups in a data set using sum of squares clustering. Biometrics 44, 23–34 (1988).
Hennig, C. Dissolution point and isolation robustness: robustness criteria for general cluster analysis methods. J. Multivar. Anal. 99, 1154–1176 (2008).
Wu, H. et al. The gut microbiota in prediabetes and diabetes: a population-based cross-sectional study. Cell Metab. 32, 379–390 (2020).
Zhong, H. et al. Distinct gut metagenomics and metaproteomics signatures in prediabetics and treatment-naïve type 2 diabetics. EBioMedicine 47, 373–383 (2019).
Liu, R. et al. Gut microbiome and serum metabolome alterations in obesity and after weight-loss intervention. Nat. Med. 23, 859–868 (2017).
Tran, L. & Greenwood-Van Meerveld, B. Age-associated remodeling of the intestinal epithelial barrier. J. Gerontol. A Biol. Sci. Med. Sci. 68, 1045–1056 (2013).
Thevaranjan, N. et al. Age-associated microbial dysbiosis promotes intestinal permeability, systemic inflammation, and macrophage dysfunction. Cell Host Microbe 21, 455–466.e4 (2017).
Arumugam, M. et al. Enterotypes of the human gut microbiome. Nature 473, 174–180 (2011).
Faust, K. et al. Microbial co-occurrence relationships in the human microbiome. PLoS Comput. Biol. 8, e1002606 (2012).
Langenberg, C., Hingorani, A. D. & Whitty, C. J. M. Biological and functional multimorbidity—from mechanisms to management. Nat. Med. 29, 1649–1657 (2023).
Stefan, N. & Schulze, M. B. Metabolic health and cardiometabolic risk clusters: implications for prediction, prevention, and treatment. Lancet Diabetes Endocrinol. 11, 426–440 (2023).
Powell-Wiley, T. M. et al. Obesity and cardiovascular disease: a scientific statement from the American Heart Association. Circulation 143, e984–e1010 (2021).
Czech, M. P. Insulin action and resistance in obesity and type 2 diabetes. Nat. Med. 23, 804–814 (2017).
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).
Bárcena, C. et al. Healthspan and lifespan extension by fecal microbiota transplantation into progeroid mice. Nat. Med. 25, 1234–1242 (2019).
Boehme, M. et al. Microbiota from young mice counteracts selective age-associated behavioral deficits. Nat. Aging 1, 666–676 (2021).
Parker, A. et al. Fecal microbiota transfer between young and aged mice reverses hallmarks of the aging gut, eye, and brain. Microbiome 10, 68 (2022).
Wang, D. D. et al. The gut microbiome modulates the protective association between a Mediterranean diet and cardiometabolic disease risk. Nat. Med. 27, 333–343 (2021).
Pedersen, H. K. et al. Human gut microbes impact host serum metabolome and insulin sensitivity. Nature 535, 376–381 (2016).
Kovatcheva-Datchary, P. et al. Dietary fiber-induced improvement in glucose metabolism is associated with increased abundance of Prevotella. Cell Metab. 22, 971–982 (2015).
Wu, G. D. et al. Linking long-term dietary patterns with gut microbial enterotypes. Science 334, 105–108 (2011).
Tett, A. et al. The Prevotella copri complex comprises four distinct clades underrepresented in Westernized populations. Cell Host Microbe 26, 666–679 (2019).
Tett, A., Pasolli, E., Masetti, G., Ercolini, D. & Segata, N. Prevotella diversity, niches and interactions with the human host. Nat. Rev. Microbiol. 19, 585–599 (2021).
Gacesa, R. et al. Environmental factors shaping the gut microbiome in a Dutch population. Nature 604, 732–739 (2022).
Brown, J. M. & Hazen, S. L. Microbial modulation of cardiovascular disease. Nat. Rev. Microbiol. 16, 171–181 (2018).
Kong, L. et al. Diabesity phenotype and the risks of cardiovascular disease and subclinical atherosclerosis: a prospective cohort study. Obesity 30, 1681–1690 (2022).
Wang, T. et al. Association of insulin resistance and β-cell dysfunction with incident diabetes among adults in China: a nationwide, population-based, prospective cohort study. Lancet Diabetes Endocrinol. 8, 115–124 (2020).
Pasolli, E. et al. Accessible, curated metagenomic data through ExperimentHub. Nat. Methods 14, 1023–1024 (2017).
Craig, C. L. et al. International physical activity questionnaire: 12-country reliability and validity. Med. Sci. Sports Exerc. 35, 1381–1395 (2003).
Matthews, D. R. et al. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 28, 412–419 (1985).
Levey, A. S. et al. A new equation to estimate glomerular filtration rate. Ann. Intern. Med. 150, 604–612 (2009).
Zhou, B. F. & Cooperative Meta-Analysis Group of the Working Group on Obesity in China. Predictive values of body mass index and waist circumference for risk factors of certain related diseases in Chinese adults—study on optimal cut-off points of body mass index and waist circumference in Chinese adults. Biomed. Environ. Sci. 15, 83–96 (2002).
World Health Organization. Regional Office for the Western Pacific. The Asia-Pacific perspective: redefining obesity and its treatment. https://apps.who.int/iris/handle/10665/206936 (Health Communications Australia, 2000).
American Diabetes Association. Diagnosis and classification of diabetes mellitus. Diabetes Care 33, S62–S69 (2010).
Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults. Executive Summary of the Third Report of the National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (Adult Treatment Panel III). JAMA 285, 2486–2497 (2001).
Yang, F. et al. Assessment of fecal DNA extraction protocols for metagenomic studies. Gigascience 9, giaa071 (2020).
Chen, S., Zhou, Y., Chen, Y. & Gu, J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890 (2018).
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).
Beghini, F. et al. Integrating taxonomic, functional, and strain-level profiling of diverse microbial communities with bioBakery 3. eLife 10, e65088 (2021).
Charrad, M., Ghazzali, N., Boiteau, V. & Niknafs, A. NbClust: an R package for determining the relevant number of clusters in a data set. J. Stat. Softw. 61, 1–36 (2014).
World Health Organization. Ageing and health. https://www.who.int/news-room/fact-sheets/detail/ageing-and-health (2022).
The UN Refugee Agency. Emergency handbook. Older persons. https://emergency.unhcr.org/protection/persons-risk/older-persons (2020).
Vieira-Silva, S. et al. Statin therapy is associated with lower prevalence of gut microbiota dysbiosis. Nature 581, 310–315 (2020).
Acknowledgements
This work was supported by the grants from the National Key R&D Program of China (2021YFA1301103), the National Natural Science Foundation of China (82088102, 91857205, 82370820 and 81930021), the ‘Shanghai Municipal Education Commission–Gaofeng Clinical Medicine Grant Support’ from Shanghai Jiao Tong University School of Medicine (20171901 Round 2) and the Innovative Research Team of High-Level Local Universities in Shanghai. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. The authors are grateful to K. Kristiansen from the University of Copenhagen for his valuable suggestions in interpreting the findings of this study.
Author information
Authors and Affiliations
Contributions
T.W., G.N., J. Li, H.Z., Y.B. and W.W. contributed to the conception and design of the study. T.W., Z.S., H.R., F.Y. and H.Z. performed bioinformatic analyses. T.W., H.Z., Z.S. and H.R. drafted the manuscript. T.W., H.Z., G.N., J. Li, Y.B. and W.W. critically revised the manuscript. T.W., G.N., Y.B. and W.W. obtained funding. T.W., M.X., J. Lu, C.Y., K.W., M.C., X.X., D.L., L.K., R.Z., J.Z., M.L., Y.X., Z.Z., Y.C., H.Y. and J.W. contributed to the acquisition or interpretation of data. All authors performed proofreading of the manuscript for important intellectual content and gave final approval of the version to be published. G.N., Y.B. and W.W. are the guarantors of this work and take responsibility for the integrity of the data and the accuracy of the data analysis.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Ethical approval
This study was approved by the Medical Ethics Committee of Ruijin Hospital (no. 2020 Lin-Lun-Shen 359), Shanghai Jiao Tong University School of Medicine. All study participants provided written informed consent before study participation. The metagenomic analysis on the JD_2014 cohort was also approved by the institutional review board of BGI Research (no. BGI-IRB 21144).
Peer review
Peer review information
Nature Medicine thanks Sean Gibbons, Cristina Menni and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editor: Michael Basson, in collaboration with the Nature Medicine team.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Validation of the five metabolic MCs in JD_2014.
a. In JD_2014 (N = 4,491), k-means-based unsupervised clustering was applied to 21 Z-score-transformed metabolic variables (including BMI, WC, SBP, DBP, FPG, PPG, HbA1c, Fins, Pins, HOMA-IR, TC, HDL-C, LDL-C, TG, ApoA-1, ApoB, ALT, AST, GGT, eGFR, and UA), revealing five robust MCs with distinctive subphenotypes. The heatmap displays the transformed values of 21 metabolic variables in all participants. b. Radar charts show the mean standardized Z-scores of metabolic variables of the JD_2014 cohort. Each MC exhibited a unique metabolic subphenotype, with MC1 characterized by a relatively healthy metabolic profile and other MCs by feature variables highlighted in colors. Colored feature variables for MCs2-5 were identified using the criteria of P < 0.05 and Cliff’s delta effect size>0.30 (two-sided Wilcoxon rank-sum test). c. Point plots show the average Jaccard similarity for the five MCs obtained using k-means clustering and random assignments in JD_2010, CM_2010, and JD_2014 (100 bootstrap resampling). d. Feature variables for each MC (MCs2-5 versus MC1) were identified using the criteria of P < 0.05 and Cliff’s delta effect size>0.3 (two-sided Wilcoxon rank-sum test) and are highlighted with an asterisk. e. The prevalence of metabolic disturbances closely aligned with the feature variables in each respective MC in JD_2010, CM_2010, and JD_2014.
Extended Data Fig. 2 Stability assessment of the five metabolic MCs in men and women.
a. Heatmaps depict the distribution of all clustering metabolic variables in JD_2010, CM_2010, and JD_2014, stratified by sex. b. Point plots show the Jaccard similarity index for each MC among men and women in JD_2010, CM_2010, and JD_2014. c. PCoA based on Euclidean distance shows the distributions of metabolic variables across the five MCs between men and women. Data are presented as mean values with standard errors. d. Bar plots show the distribution of MCs between men and women. The percentages of each MC within each sex group are presented. Percentages might not sum to 100% because of rounding.
Extended Data Fig. 3 Associations between metabolic MCs, age, and the risk of incident CVD in JD_2010 and CM_2010.
a. The distributions of five metabolic MCs varied by age stratifications in JD_2010 and CM_2010 cohorts. Percentages might not sum to 100% because of rounding. b. Forest plots display the multivariable-adjusted HRs with 95% CIs for incident CVD associated with older versus younger age groups in the overall, low-risk (MCs123) and high-risk (MCs45) participants in JD_2010 and CM_2010. Cox proportional hazards regression models were adjusted for MCs (not for MCs-stratified models), sex, educational attainment, smoking status, alcohol drinking status, and physical activity. The error bars represent the two-sided 95% CIs determined by multivariable analysis in Cox regression.
Extended Data Fig. 4 Distinct gut microbiome signatures associated with age and metabolism in JD_2014.
a. Box plots display the log 10-transformed relative abundances of the four MCs-related core genera (Fig. 3b), including Faecalibacterium, Alistipes, Ruminococcus, and Megamonas, among five MCs (Top) and between age groups (Bottom). b. Box plots display the levels of microbial richness, two uniqueness indices (Kendall and Aitchison), and inter-individual beta diversity (Bray-Curtis dissimilarity at the genus level) across five MCs and between age groups. The central line represents the median value, the box represents the interquartile range (IQR), and the whiskers extend to the minimum and maximum values within 1.5 times the IQR from the lower and upper quartiles, respectively. Statistical significance was determined using the Kruskal-Wallis test followed by Dunn’s test (two-sided) for MCs groups, and the two-sided Wilcoxon rank-sum test for age groups. P < 0.05 is considered statistical significance. c. Spearman’s rank correlation analysis (two-sided) between uniqueness indices and age or microbial richness. The four uniqueness indices were calculated using different distance metrics to capture different aspects of within-sample microbial variations: Aitchison distance, Bray-Curtis dissimilarity, Jaccard distance, and Kendall’s tau coefficient (details see Methods). Spearman’s rho and P-values are presented. d. The heatmap shows the associations between two uniqueness indices (Aitchison and Kendall) and age-related microbial species, with colors indicating respective phyla. Statistical significance was determined using Spearman’s correlation (two-sided), with a BH-adjusted P < 0.05 denoted by an asterisk. e. The heatmap shows the associations between two uniqueness indices (Aitchison and Kendall) and MCs-related microbial species, with colors indicating respective phyla. Statistical significance was determined using Spearman’s correlation (two-sided), with a BH-adjusted P < 0.05 denoted by an asterisk.
Extended Data Fig. 5 Estimation of the impact of commonly used oral medications on gut microbiota in JD_2014.
a. Bar plot shows the explained variance (R2, one-sided) of ten types of commonly used oral medications (adjusted R2, one-sided) on gut microbiota determined by PERMANOVA with Bray-Curtis dissimilarities at the genus level. Adjusted R2 for each type of medication was estimated with adjustment for other medicines. Five out of ten oral medications showed significant associations with gut microbiota, including AGIs, biguanides, ARBs, TZDs, and calcium antagonists. Asterisks denote a BH-adjusted P < 0.05. b. The bar plot shows the distribution of drug-negative and drug-positive participants (without or with the use of any of the five types of medications or statins) in each MC. The number and percentage of drug-negative and drug-positive participants in each MC are presented. c. Top: Heatmap displays the significantly differential species between drug-negative MC1 (metabolically healthy cluster) and drug-negative MCs2-5 (metabolically unhealthy clusters), and between drug-negative MC1 and drug-positive unhealthy clusters (MCs2-5) (two-sided Wilcoxon rank-sum test, BH-adjusted P < 0.05). Bottom: Heatmap shows the significantly differential species between metabolically unhealthy participants (MCs2-5) with and without a certain type of medication (two-sided Wilcoxon rank-sum test, BH-adjusted P < 0.05). Asterisks denote a BH-adjusted P < 0.05.
Extended Data Fig. 6 Identification of age-related gut microbial pathways in JD_2014.
Left: Bar plot shows the gut microbial pathways significantly associated with age in MC1 using the partial Spearman’s correlation (two-sided) after adjusting for sex, educational attainment, smoking status, alcohol drinking status, physical activity, fruit and vegetable intake, and medication use. Red indicates significant positive associations and blue indicates significant negative associations. Right: Heatmap shows associations between relative abundances of age-related species and pathways in JD_2014 using Spearman’s correlation analysis (two-sided). A BH-adjusted P < 0.05 is considered as statistical significance. Asterisks denote a BH-adjusted P < 0.05.
Extended Data Fig. 7 Validation of the age-related gut microbial species and pathways in four external datasets.
Heatmaps show the associations of age with age-related gut microbial species and pathways in JD_2014 and four external validation cohorts obtained from the curatedMetagenomicData (details see Methods). For JD_2014, models were adjusted for sex, educational attainment, smoking status, alcohol drinking status, physical activity, fruit and vegetable intake, and medication use; for external datasets, models were adjusted for sex, using partial Spearman’s correlation analysis (two-sided). The ‘*’ denotes a BH-adjusted P < 0.05 and the ‘#’ denotes a P < 0.05. Purple indicates significant positive associations, green indicates significant negative associations, and gray indicates non-significant associations.
Extended Data Fig. 8 Associations between MA and the risk of incident CVD in JD_2014.
a. Cox proportional hazards regression models were applied to estimate the associations between MA groups and incident CVD during a 6.8-year follow-up period among the entire JD_2014 cohort (45–91 years), the younger age group (45- < 60 years), and the older age group (≥60 years), as well as the narrower age subgroups stratified by 10-year age spans (for example, 60–70 years, 61–71 years, and 62–72 years). Forest plots display the unadjusted HRs (95% CIs) (Left) and multivariable-adjusted HRs (95% CIs) (Right) for incident CVD associated with the HMA group compared with the LMA group. b. Forest plots display the multivariable-adjusted HRs (95% CIs) for incident CVD associated with the MCs45-HMA group compared with the overall MCs123 group among the older age group and the narrower age subgroups stratified by 10-year age spans. c. Forest plots display the multivariable-adjusted HRs (95% CIs) for incident CVD associated with the MCs45-LMA group compared with the overall MCs123 group among the older age group and the narrower age subgroups stratified by 10-year age spans. For a, b, and c, multivariable-adjusted models were adjusted for age, sex, educational attainment, smoking status, alcohol drinking status, physical activity, fruit and vegetable intake, and medication use. BH-adjusted P < 0.05 is considered statistical significance. The error bars represent the two-sided 95% CIs determined by multivariable analysis in Cox regression.
Extended Data Fig. 9 Distinct association patterns of Prevotella with metabolically healthy enriched microbial species among three enterotypes in JD_2014.
a. Histograms are the mean of relative abundance or prevalence (%) of Prevotella and the annotated species in three enterotypes (ETs) in older (aged≥60 years) MCs45 participants. b. Spearman’s rank correlation analysis (two-sided) reveals distinct association patterns between Prevotella and the other top abundant genera within the three ETs in older MCs45 participants. These include Megamonas, Bacteroides, Faecalibacterium, Eubacterium, Roseburia, Ruminococcus, Blautia, Klebsiella, and Escherichia. Asterisks denote a BH-adjusted P < 0.05. c. Spearman’s rank correlation analysis (two-sided) reveals distinct association patterns between Prevotella, P. copri, and the metabolically healthy enriched species within the three ETs in older MCs45 participants. The ‘*’ denotes a BH-adjusted P < 0.05, and the ‘#’ denotes a P < 0.05.
Extended Data Fig. 10 Interaction between host age and metabolic MCs on gut microbial composition in JD_2014.
a. Bar plot shows the significant impact of host age on gut microbiota in MCs123 and MCs45 groups, separately, determined by PERMANOVA (one-sided) with Bray-Curtis dissimilarities at the genus level. Asterisks denote a P < 0.05. b. The relationships between age and a set of age-related gut microbial variables, including 55 age-related species, MA, and microbial richness, in MCs subgroups, separately (MCs123 and MCs45, MC1 to MC5). Boxplot shows the absolute Spearman’s rho (two-sided) between MCs123 and MCs45. Boxes represent the interquartile ranges (IQRs) between the first and third quartiles, the center line represents the median, and the whiskers extend 1.5 times the IQR from the top and bottom of the box. Asterisks denote a BH-adjusted P < 0.05. c. Detailed information on health metrics used for stratifying healthy and unhealthy groups in JD_2014. d. The point plot shows the partial Spearman’s correlations (two-sided) between age and age-related gut microbial signatures, including MA, richness, four uniqueness indices, and the relative abundances of genera Bacteroides and Prevotella, within healthy and unhealthy groups. All models were adjusted for sex, educational attainment, smoking status, alcohol drinking status, physical activity, fruit and vegetable intake, and medication use (not for medication-stratified groups). P < 0.05 is considered as statistical significance. Red indicates significant associations in unhealthy groups, green indicates significant associations in healthy groups, and gray indicates non-significant associations. e. PERMANOVA analysis (one-sided) indicates the interaction effect between age and MCs (MCs1 and MCs45) on gut microbial composition. The interaction P values from the unadjusted model (P = 0.01), the model including sex and all environmental factors excluding medications (adjusted Model 1, P = 0.04), and the fully adjusted Model (adjusted Model 1 + medications, P = 0.07) are shown.
Supplementary information
Supplementary Tables 1–16
Supplementary Table 1a. Baseline characteristics of participants in the JD_2010, CM_2010 and JD_2014 cohorts according to five metabolic MCs. Supplementary Table 1b. Characteristics of participants in the JD_2010 and JD_2014 cohorts in 2010. Supplementary Table 2. Distribution of the clustering metabolic variables for each metabolic MC. Supplementary Table 3. Prevalence of related metabolic disorders in each metabolic MC. Supplementary Table 4. Associations of age, MC, MA and ET with incident CVD risk. Supplementary Table 5. Associations of gut microbiota with host and environmental factors (PERMANOVA, Bray–Curtis dissimilarity at the genus level, one-sided F-test). Supplementary Table 6. Summary of drug-negative and drug-positive individuals within five MCs. Supplementary Table 7. Partial Spearman’s correlations between species abundances and metabolic variables. Supplementary Table 8. Identification of gut microbial species associated with metabolic health. Supplementary Table 9. Identification of gut microbial features (richness, uniqueness indices and species abundances) associated with host age. Supplementary Table 10. Spearman’s correlations between species abundances and two uniqueness indices. Supplementary Table 11. Identification of gut microbial pathways associated with host age. Supplementary Table 12. List of gut metagenomic datasets used in this study from curatedMetagenomicData. Supplementary Table 13. Partial Spearman’s correlations between abundances of microbial features (species and pathways) and host age in four external cohorts. Supplementary Table 14. Partial Spearman’s correlations between MA and age-related species in the JD_2014 and validation cohorts. Supplementary Table 15. Distribution of metabolic variables between low-MA and high-MA groups within each metabolic cluster. Supplementary Table 16. Spearman’s correlations between Prevotella/P. copri and MC1-enriched species within each ET.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Wang, T., Shi, Z., Ren, H. et al. Divergent age-associated and metabolism-associated gut microbiome signatures modulate cardiovascular disease risk. Nat Med 30, 1722–1731 (2024). https://doi.org/10.1038/s41591-024-03038-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41591-024-03038-y