Abstract
Background and Aims Life’s Essential 8 (LE8) is a comprehensive construct of cardiovascular health. Yet, little is known about LE8 score, its metabolic correlates, and their predictive implications among Black Americans and low-income individuals.
Methods In a nested case-control study of coronary heart disease (CHD) among 598 Black and 596 White low-income Americans, we estimated LE8 score, conducted untargeted plasma metabolites profiling, and used elastic net with leave-one-out cross-validation to identify metabolite signature (MetaSig) of LE8. Associations of LE8 score and MetaSig with incident CHD were examined using conditional logistic regression. Mediation effect of MetaSig on the LE8-CHD association was also examined. The external validity of MetaSig was evaluated in another nested CHD case-control study among 598 Chinese adults.
Results Higher LE8 score was associated with lower CHD risk [standardized OR (95% CI)=0.61 (0.53-0.69)]. The identified MetaSig, consisting of 133 metabolites, showed strong correlation with LE8 score (r=0.61) and inverse association with CHD risk [OR (95% CI)=0.57 (0.49-0.65)], robust to adjustment for LE8 score and across participants with different sociodemographic and health status (ORs: 0.42-0.69; all P<0.05). MetaSig mediated a large portion of the LE8-CHD association: 53% (32%-80%) (P<0.001). Significant associations of MetaSig with LE8 score and CHD risk were found in validation cohort [r=0.49; OR (95% CI)=0.57 (0.46-0.69)].
Conclusions Higher LE8 score and its MetaSig were associated with lower CHD risk among low-income Black and White Americans. Metabolomics may offer an objective and comprehensive measure of LE8 score and its metabolic phenotype relevant to CHD prevention among diverse populations.
Introduction
Coronary heart disease (CHD) is a leading cause of morbidity and mortality in the United States (US) and worldwide, with significant and persistent sociodemographic disparities1, 2. To reduce the burden and life lost due to CHD and other cardiovascular disease (CVD), the American Heart Association (AHA) has recently proposed Life’s Essential 8 (LE8) to assess and promote cardiovascular health (CVH) in individuals and populations3. LE8 includes 4 health behaviors (healthy diet, participation in physical activity, avoidance of nicotine, and healthy sleep [a new component of LE8]) and 4 health factors (weight, blood lipids, glucose, and blood pressure), and has a new scoring system with continuous scale to better reflect inter-individual differences. While a higher LE8 score has been recently associated with lower CVD incidence and mortality4–7, few studies have evaluated LE8 score and its association with incident CVD among Black and White Americans who have low socioeconomic status (SES) and face disproportionate CVD burdens. In addition, although some potential mechanisms have been identified (e.g., reduced inflammation and atherosclerosis)8, 9, beyond those known CVD risk pathways, mechanisms and inter-individual differences underlying the cardioprotective effects of LE8 and its included health behaviors and health factors are not fully understood.
Metabolite profiling (“metabolomics”) comprehensively measures small-molecule metabolites (including substrates, intermediates, and end products) in biological samples and represents a powerful tool for mechanistic investigation, novel biomarker discovery, and precision medicine10, 11. Metabolite profiling of blood samples may improve assessments of individuals’ alignment with LE8, particularly for behavioral factors that are prone to survey and recall biases. In addition, circulating metabolites related to LE8 may capture varied individual metabolic responses to LE8, providing novel mechanistic insights into its cardioprotective effects and informing precision medicine. While previous studies have identified metabolites related to the components of LE8, including diet12–16, physical activity17–19, tobacco exposure20–22, sleep23–26, and body mass index (BMI)27, 28, to our knowledge, no study has applied untargeted plasma metabolomics to identify a comprehensive metabolite signature (MetaSig) for LE8 to enable studies with incident CHD. Given that those health behaviors and factors often correlate and interact with each other, investigating whether plasma metabolomics could provide a good objective assessment of individuals’ alignment with and metabolic responses to overall LE8 and uncovering potential pathways linking LE8 to incident CHD is highly warranted.
Here, leveraging a case-control study of CHD nested within the Southern Community Cohort Study (SCCS) involving 598 Black Americans and 596 White Americans, we assessed the LE8 score and its MetaSig and examined their relations to incident CHD. The results were further replicated in another nested CHD case-control study of racially and geographically different population: 598 Chinese adults from the Shanghai Women’s Health Study and Shanghai Men’s Health Studies (SWMHS). In addition, we identified MetaSigs for the health behaviors and health factors recommended in the LE8 and evaluated their associations with incident CHD.
Methods
Study participants
This study was based on a nested case-control study of CHD within the SCCS29. Briefly, the SCCS enrolled 84,735 primarily low-income, uninsured/underinsured Black and White Americans aged 40-79 years from 12 southeastern US states between 2002-2009, with >50% having household income less than $15,000/y and ∼86% were uninsured or underinsured29. Participants were surveyed for a wide range of information at baseline and followed up regularly for morbidity and mortality outcomes. Venous blood samples were collected at baseline, and plasma samples were aliquoted and stored at -80°C for long-term use. A nested case-control study of CHD within SWMHS was used as the validation cohort, which enrolled 74,940 women, aged 40-70 years, and 61,480 men, aged 40-74 years, from Shanghai, China between 1996–2000 and 2002–2006, respectively30, 31. These cohort studies were approved by the Institutional Review Boards of the Vanderbilt University Medical Center, Meharry Medical College, and/or Shanghai Cancer Institute. Informed consent was obtained from all enrolled participants.
For the nested case-control studies, participant inclusion criteria were 1) no history of CHD, stroke, heart failure, cancer, or end-stage renal disease at baseline; 2) available baseline plasma samples and data on fasting time and in SCCS, time between sample collection and lab processing; 3) no use of antibiotics nor cold/flu in last 7 days before blood collection. In addition, to facilitate CHD case identification and adjudication, in SCCS, participants were eligible for Centers for Medicare & Medicaid Services (CMS) and had ≥2 claims after cohort enrollment; in SWMHS, participants’ medical records were accessible for our study. In SCCS, nonfatal CHD cases were identified through CMS, including acute myocardial infarction, coronary revascularization, and other acute CHD, and CHD deaths were identified through the National Death Index. In SWMHS, CHD cases were first identified by self-reported diagnoses during follow-up visits and then confirmed by medical records. In each race (Black, White, or Chinese) and gender (male or female), 150 incident CHD cases were randomly selected and 1:1 matched with controls who had no CHD, heart failure, stroke, nor cancer at the time of case diagnosis, by enrollment age (±2 years), fasting time (±2 hours), and time between sample collection and lab processing (±4 hours, for SCCS samples; all SWMHS samples were processed within 6 hours after collection). After excluding eight plasma samples that did not pass the metabolite profiling quality control, a total of 1792 participants, including 597 pairs of CHD cases and controls in SCCS (299 pairs of Black Americans and 298 pairs of White Americans), and 299 pairs of Chinese adults in SWMHS were included in the present study (Fig. 1).
Life’s Essential 8
We constructed each LE8 component metric and total LE8 score (ranging from 0 to 100) according to AHA guidelines3, 32, with some modifications based on the characteristics of our data, as shown in Table S1. Briefly, dietary quality was assessed by the Dietary Approaches to Stop Hypertension (DASH) score33–35. Physical activity was measured by total minutes of leisure-time moderate and vigorous physical activity per week, with each minute of vigorous physical activity counted as 2 minutes toward the total minutes. Nicotine exposure was assessed by active tobacco smoking and secondhand smoking exposure. Sleep health was measured by average self-reported sleep hours per day. BMI was calculated as weight (kg)/height (m)2. Blood lipids component was assessed by non-HDL cholesterol (mg/dL) (plasma total cholesterol minus HDL cholesterol). Because concentrations of fasting glucose and HbA1c were not measured in our study, blood glucose component was scored based on diabetes diagnosis and medications and relative abundance of plasma glucose measured simultaneously with other metabolites (see following metabolite profiling methods). Blood pressure was unavailable in SCCS, thus blood pressure component was scored based on hypertension status and medications. Total LE8 score was obtained by calculating the arithmetic mean of individual component scores. We also calculated scores reflecting alignment with health behaviors and health factors.
Metabolite profiling
Baseline plasma samples of selected CHD case-control pairs were retrieved and placed adjacently in the same assay batch. Laboratory persons were blinded to the case-control status of samples. Untargeted metabolite profiling was performed using ultra-high-performance liquid chromatography (UHPLC) coupled with tandem mass spectrometry (MS) by Metabolon Inc. (Morrisville, NC, USA) following a standard assay protocol36. Briefly, plasma samples were extracted with methanol and split into four aliquots for analysis by UHPLC-MS/MS in both positive and negative ion modes using a combination reverse phase and HILIC chromatography methods. Metabolites were identified by automated comparison of mass spectra features to a reference library of >4,000 authenticated standard compounds followed by visual inspection for quality control. Peaks were quantified using area-under-the-curve. A total of 1502 metabolites were detected in our samples. The majority of metabolites (>80%) were annotated based on internal standards. Metabolites that were annotated only by a match to a known MS spectrum or chemical formula were marked by ‘*’ and ‘**’, respectively. We excluded metabolites detected in <10% of participants, resulting in 1322 metabolites. Metabolites with missing values were imputed by half of the minimal value in the non-missing samples. The values of all metabolites were log-transformed and standardized to mean 0 and unit variance.
Statistical analysis
The characteristics of the study participants were presented as mean (standard deviation [SD]) for continuous variables and frequency (percentage) for categorical variables. Spearman correlations between LE8 score and individual component scores were assessed. Elastic net regression was used to identify metabolites associated with LE812, 14, 37, with hyperparameters determined via a ten-repeated tenfold cross-validation framework, using R package caret (version: 6.0-88)38 and glmnet (version: 4.1-3)37, 39. The MetaSig was calculated through a leave-one-out cross-validation approach. In addition, we externally validated the identified signature in SWMHS using weights of the selected metabolites from SCCS. We calculated Spearman correlations between LE8 score and MetaSig among all participants and by race, sex, age (≥60y/<60y), incident CHD status, diabetes status, hypertension status, dyslipidemia status, and fasting status.
We then examined the associations of LE8 and its MetaSig with risk of CHD using conditional logistic regression, adjusting for age, education, income, alcohol intake, and family history of CHD. We also included LE8 score and MetaSig in the same model to assess their independent associations with CHD and potential mediating effect of MetaSig on LE8-CHD association. The potential multicollinearity was assessed by variance inflation factor (VIF) using R package car (version: 3.1-0) with VIF >10 indicating multicollinearity among variables. The ranges of VIF for LE8 score and its MetaSig were 1.08-1.37, confirming no multicollinearity. The mediation analysis was performed using R package mediation (version: 4.5.0)40. Subgroup analyses to evaluate the associations of LE8 and MetaSig with incident CHD were performed by race, age group, sex, education, household income, diabetes status, hypertension status, and dyslipidemia status, with P value for interaction obtained from the corresponding interaction term in the model.
The same methods were used to identify MetaSigs for health behaviors and health factors included in LE8 and evaluate their associations with incident CHD, potential mediation effects, and their external validity. Mutual adjustments of health behaviors and health factors score were performed in corresponding statistical models. All analyses were performed using R (version 4.1.1). Two-sided P < 0.05 was considered statistically significant. An overview of our study design is presented in Fig. 1.
Results
Baseline characteristics of study participants
The mean age at baseline (blood collection) was 55 years in our study participants (Table 1). The mean (SD) of LE8 score was 48.1 (12.3) in Black women, 50.0 (13.6) in Black men, 48.1 (13.4) in White women, and 48.5 (14.4) in White men. The median follow-up time for incident CHD cases was 5 (interquartile range: 3-8) years in SCCS. Incident CHD cases had significantly lower total LE8 score, health behaviors score, and health factors score than controls (all P <0.05). The characteristics of participants in SWMHS (mean age: 61 years; mean LE8 score: 57.2 in women and 50.7 in men) are shown in Table S2. There were moderate correlations between total LE8 score and individual component scores (r ranged from 0.22 with smoking to 0.47 with BMI and blood pressure scores in SCCS; Fig. S1).
Metabolite signature of LE8
We identified 133 metabolites related to LE8 (top 30 are shown in Fig. 2A; the full list can be found in Table S3). The MetaSig was strongly correlated with LE8 score (r = 0.61, P<0.001; Fig. 2B), meanwhile variations in MetaSig were observed among individuals with the same LE8 score, demonstrating interindividual difference in metabolic phenotype of LE8. MetaSig was also correlated with LE8 score in SWMHS (r = 0.49, P<0.001; Fig. 2C). Stratified analyses showed that correlations between LE8 score and MetaSig were consistent regardless of age, sex, race, fasting status, diabetes status, hypertension status, dyslipidemia status, and incident CHD status (r ranged from 0.54 to 0.64; Table S4), suggesting the robustness of our identified LE8 metabolite signature across participants with different sociodemographic backgrounds and metabolic disease status.
Associations with incident CHD
Higher LE8 score and its MetaSig were significantly associated with lower risk of CHD: standardized multivariable-adjusted odds ratio (OR) = 0.61 (95% CI: 0.53-0.69) for LE8 score and 0.57 (0.49-0.65) for MetaSig; both P < 0.001 (Table 2). Sensitivity analysis showed that the MetaSig-CHD associations did not change after excluding any individual metabolites from the signature (Table S5). After further adjusting for LE8 score, the MetaSig-CHD association was only slightly attenuated [OR (95% CI) = 0.66 (0.55-0.78); P<0.001; Table 2], suggesting circulating metabolites may complement LE8 assessment and contribute to CHD risk beyond LE8 score. On the other hand, the LE8-CHD association was moderately attenuated after adjusting for MetaSig [OR (95% CI) = 0.78 (0.66-0.92), P=0.003]. Mediation analysis showed that MetaSig mediated a large portion of the LE8-CHD association [53% (32%-80%); Pmediation<0.001; Fig. 2D].
Both LE8 and its MetaSig were inversely associated with CHD risk in subpopulations by race, age group, education, income, diabetes status, hypertension status, and dyslipidemia status (Pinteraction>0.05; Fig. 3), with stronger associations observed in women than in men [for LE8, OR (95% CI) = 0.53 (0.42-0.66) in women and 0.66 (0.55-0.80) in men, Pinteraction = 0.017; for MetaSig, 0.48 (0.38-0.6) in women and 0.65 (0.54-0.78) in men, Pinteraction = 0.016].
The association of MetaSig with incident CHD was replicated in SWMHS (Table 2), with OR (95% CI) = 0.57 (0.46-0.69) and 0.69 (0.55-0.86) after further adjusting for LE8 score (both P<0.001). The MetaSig also mediated a considerable portion of the LE8-CHD association in SWMHS [27.4% (10%-47%); Pmediation<0.001; Fig. 2E].
Metabolite signatures of health behaviors and health factors and associations with CHD
We further identified MetaSigs for health behaviors (Fig. S2A) and health factors in SCCS (Fig. S3A; full lists of metabolites are shown in Table S3), which showed strong correlations with health behaviors score (r = 0.59, P<0.001; Fig. S2B) and with health factors score (r = 0.76, P<0.001; Fig. S3B). Significant correlations were also found among participant subgroups (Table S4), suggesting the robustness of identified metabolite signatures for health behaviors and health factors. In addition, there were significant inverse associations of health behaviors score, health factors score, and their related signatures with risk of CHD (all P<0.001; Table 2). Specifically, standardized OR (95% CI) was 0.73 (0.63-0.85) for health behaviors MetaSig and 0.57 (0.49-0.66) for health factors MetaSig. Similarly, metabolites mediated large portions of the health behaviors-CHD association [43.9% (13.9%-101%), P=0.004; Fig. S2D] and health factors-CHD association [53.2% (24.8%-89%), Pmediation<0.001; Fig. S3D]. Further, all results on health behaviors MetaSig and health factors MetaSig were replicated in SWMHS (Table 2, Fig. S2, and Fig. S3).
Discussion
Leveraging untargeted plasma metabolites data in a nested case-control study among low-income Black and White Americans, we identified a metabolite signature that could reflect LE8 score and was associated with incident CHD, even after adjusting for LE8 and among participants with varied sociodemographic and metabolic health status, suggesting that circulating metabolite profiling may be used to help assess LE8 alignment across diverse populations and offer additional information on CVH (e.g., inter-individual metabolic phenotypes related to LE8). We also identified MetaSigs for health behaviors and health factors and found consistent results showing that circulating metabolites could reflect the alignment with those recommendations and underlying metabolic phenotypes, which were further linked to incident CHD across diverse populations. All the results were further replicated in another nested case-control study of CHD among Chinese adults. Our findings demonstrate the potential utility of blood metabolomics to improve the assessment of LE8 and its underlying metabolic variations that are linked to incident CHD among sociodemographically diverse populations, towards advancing precision medicine and addressing disparities in CVH.
LE8 is the American Heart Association’s updated and enhanced guideline to measure and promote CVH for individuals and populations3. The beneficial associations of following the LE8 with lower risks of CHD, CVD, and related mortality have been demonstrated in recent studies4–7. However, multi-racial/ethnic populations with low SES remain underrepresented in research studies, even though they have persistently experienced worse CVH and CVD outcomes, as well as systemic disadvantages to improve CVH, than White and middle-class Americans2, 41–43. Leveraging resources from SCCS, a large cohort of predominantly low-income Black and White Americans (in the present study: ∼65% with household income <$15,000/y and ∼95% with household income <$25,000/y), our study assessed CVH based on LE8 and evaluated the association of LE8 score with incident CHD. We found that higher LE8 score (per SD increase) was associated with ∼40-50% lower risk of CHD among Black Americans and individuals with low SES. While LE8 provides a comprehensive approach to quantify CVH, its assessment involves a series of procedures such as questionnaires, anthropometric and blood pressure measures, and blood draw. Particularly, health behaviors (diet, physical activity, smoking, and sleep) are usually assessed by questionnaires, which are time-consuming and prone to measurement errors and low compliance (particularly in the clinical setting). Also, LE8 score cannot capture varied individual metabolic responses to lifestyle exposures. Hence, we incorporated untargeted plasma metabolomics data, and for the first time, identified a robust metabolite signature of LE8 and then examined its association with incident CHD.
Metabolomics has been demonstrated as a powerful tool for improving exposure assessment and identifying potential novel biomarkers and mechanistic pathways in population studies, given its high-throughput characterization of thousands of metabolites in a small amount of biological samples44. Our study provides new evidence that plasma metabolite profiling may provide objective and comprehensive measures of LE8 and CVH among racially and geographically diverse populations. The identified metabolite signature may complement LE8 scores, improve the precision to stratify individuals with different future CHD risk, and potentially facilitate personalized CHD prevention strategies.
Our identified MetaSig consists of metabolites reflecting participants’ alignment with LE8 health behaviors and metabolic health status, majority of which are lipids and amino acids, and many of them have been linked to diet12–16, physical activity17–19, smoking20–22, sleep23–26, obesity27, 28, 45, or composite lifestyles scores46–49 in previous studies. For example, (2,4 or 2,5)-dimethylphenol sulfate, tartronate, and ethyl beta-glucopyranoside are derived from plant-based foods; cotinine is the major metabolite of nicotine from tobacco smoking; cholesterol, sphingomyelin, cortisol, and 1-palmitoleoylglycerolare are related to blood lipids; mannose, metformin, and fructosyllysine are related to prevalent diabetes and blood glucose. Particularly, drug metabolites including hydrochlorothiazide and metformin reflect antihypertensive and antidiabetic medications defined in LE8. Nevertheless, the associations of MetaSig with LE8 score and incident CHD were consistent among participants with or without history of hypertension or diabetes (Fig. 3 and Table S4), and the MetaSig-CHD association did not change after excluding those drug metabolites from the MetaSig [OR (95% CI)=0.55 (0.48-0.64); P<0.001]. Moreover, excluding 29 unknown metabolites (X-) from the MetaSig also did not change the MetaSig-CHD association [OR (95% CI)=0.58 (0.50-0.66); P<0.001]. Notably, the MetaSig also contains microbial metabolites, e.g., maltotetraose, anthranilate, indolebutyrate, and bile acids (taurohyocholate, glycodeoxycholate 3-sulfate, 3b-hydroxy-5-cholenoic acid, and glycohyocholate), suggesting the role of gut microbiome in host’s CVH, which cannot be captured by questionnaires or measurements of glucose, cholesterol, or blood pressure. Moreover, several metabolic pathways related to CVH and CVD development were highlighted. For example, anthranilate, indolebutyrate, picolinate, and serotonin are members of tryptophan metabolism pathway50–52; taurohyocholate, glycodeoxycholate 3-sulfate, 3b-hydroxy-5-cholenoic acid, and glycohyocholate belong to secondary bile acid metabolism pathway53, 54; alpha-tocopherol, delta-tocopherol, and gamma-tocopherol/beta-tocopherol are vitamin E derivatives through tocopherol metabolism pathway55, 56.
Importantly, the MetaSig was related to future CHD risk regardless of participants’ age, sex, race, SES, metabolic disease history, and even after adjustment for LE8. Further analyses indicated that circulating metabolites could play a substantial mediating role linking LE8 and reduced CHD risk. Moreover, our findings were replicated in a racially and geographically different population, suggesting external validity and potential generalizability of our findings. Taken together, our findings demonstrated that circulating metabolites could complement LE8 to improve the precision of CVH assessment and predict CHD risk among sociodemographically and geographically diverse populations.
To our knowledge, this is the first study that assessed the LE8 score, constructed its metabolite signature, and associated LE8 score and its MetaSig with incident CHD in Black and White Americans with low SES. Besides its novelty and inclusion of populations facing socioeconomic challenges and health disparities, other strengths of our study include its prospective design, comprehensive profiling of >1500 blood metabolites for a broad coverage and improved ability to construct metabolite signature for LE8, and robustness of results across populations with different sociodemographic and health status. Meanwhile, several limitations of our current study need to be acknowledged. First, as MetaSig of LE8 was identified using cross-sectional data from baseline blood samples, we cannot be certain as to the directionality of LE8-MetaSig association, and mediation analysis assumed that LE8 score preceded MetaSig. Although several population- or animal-based studies have shown the causal effects of LE8 components on blood metabolites21, 26, 57–60, given that blood metabolites might precede some LE8 components, the longitudinal association between LE8 adherence and circulating metabolites should be investigated. Second, given the observational nature of our study, the causality is unable to be confirmed. However, the prospective design reduces the concern of reverse causation for the LE8/MetaSig-CHD association. Third, we cannot rule out the influence of residual confounding on the LE8/MetaSig-CHD association, although we have adjusted for and stratified by major CHD risk factors. Fourth, the concentrations of fasting glucose and HbA1c and SBP and DBP were not measured in SSCS; thus, glucose score and blood pressure score were defined based on history of diabetes or hypertension and use of medications, which may influence the accuracy of LE8 score. Finally, the nested case-control design may overestimate the predictive ability of the LE8 score and its MetaSig. Therefore, our results should be further validated in other prospective cohort studies.
In summary, our study assessed the LE8 score and identified MetaSig of LE8 among low-income Black and White Americans. We found that both LE8 score and its MetaSig were inversely associated with risk of CHD, consistently among participants with varied sociodemographic and metabolic health status. Our identified metabolite signature may provide an objective and comprehensive measure of LE8 and its metabolic underpinning, which may help improve the precision of CVH assessment and facilitate more effective and personalized CHD prevention strategies in diverse populations. Further examination of our identified metabolites may improve understanding of biological mechanisms as how following LE8 benefits CHD prevention.
Data Availability
All data produced are available online at https://www.southerncommunitystudy.org/ for Southern Community Cohort Study and https://swhs-smhs.app.vumc.org/index.php for Shanghai Women's Health Study & Shanghai Men's Health Study.
Funding
The Southern Community Cohort Study is funded by U01 CA202979, the Shanghai Women’s Health Study is funded by UM1 CA182910, and the Shanghai Men’s Health Study is funded by UM1 CA173640 from the National Cancer Institute (NCI) at the National Institutes of Health (NIH). Biospecimens of these three cohort studies are managed by the Survey and Biospecimen Shared Resource, which is supported in part by the Vanderbilt-Ingram Cancer Center (P30 CA68485). This study is supported by R01 HL149779 from the National Heart, Lung, and Blood Institute (NHLBI) at the NIH. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Conflict of interest
The authors declared no conflict of interest.
Acknowledgments
The authors thank the study participants of the SCCS and the SWMHS.
Footnotes
One author information has been changed.