Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Multi-system trajectories and the incidence of heart failure in the Framingham Offspring Study

  • Cara E. Guardino,

    Roles Investigation, Methodology, Writing – original draft

    Affiliations Division of Cardiology, Warren Alpert Medical School of Brown University and Lifespan Cardiovascular Institute, Providence, Rhode Island, United States of America, Department of Medicine, Boston University School of Medicine and Boston Medical Center, Boston, Massachusetts, United States of America

  • Stephanie Pan,

    Roles Data curation, Formal analysis, Writing – review & editing

    Affiliation Department of Biostatistics, Boston University School of Public Health, Boston, Massachusetts, United States of America

  • Ramachandran S. Vasan,

    Roles Investigation, Writing – review & editing

    Affiliations National Heart, Lung, and Blood Institute’s and Boston University’s Framingham Heart Study, Framingham, Massachusetts, United States of America, Section of Preventive Medicine and Epidemiology, Boston University School of Medicine, Boston, Massachusetts, United States of America, Section of Cardiology, Boston University School of Medicine, Boston, Massachusetts, United States of America, Department of Epidemiology, Boston University School of Public Health, Boston, Massachusetts, United States of America

  • Vanessa Xanthakis

    Roles Conceptualization, Investigation, Methodology, Supervision, Writing – review & editing

    vanessax@bu.edu

    Affiliations Department of Biostatistics, Boston University School of Public Health, Boston, Massachusetts, United States of America, National Heart, Lung, and Blood Institute’s and Boston University’s Framingham Heart Study, Framingham, Massachusetts, United States of America, Section of Preventive Medicine and Epidemiology, Boston University School of Medicine, Boston, Massachusetts, United States of America

Abstract

Background

Heart failure is a multi-system disease, with non-cardiac systems playing a key role in disease pathogenesis.

Objective

Investigate whether longitudinal multi-system trajectories incrementally predict heart failure risk compared to single-occasion traits.

Methods

We evaluated 3,412 participants from the Framingham Heart Study Offspring cohort, free of heart failure, who attended examination cycle 5 and at least one examination between 1995–2008 (mean age 67 years, 54% women). We related trajectories for the following organ systems and metabolic functions to heart failure risk using Cox regression: kidney (estimated glomerular filtration rate), lung (forced vital capacity and the ratio of forced expiratory volume in one second/forced vital capacity), neuromotor (gait time), muscular (grip strength), cardiac (left ventricular mass index and heart rate), vascular function (pulse pressure), cholesterol (ratio of total/high-density lipoprotein), adiposity (body mass index), inflammation (C-reactive protein) and glucose homeostasis (hemoglobin A1c). Using traits selected via forward selection, we derived a trajectory risk score and related it to heart failure risk.

Results

We observed 276 heart failure events during a median follow up of 10 years. Participants with the ‘worst’ multi-system trajectory profile had the highest heart failure risk. A one-unit increase in the trajectory risk score was associated with a 2.72-fold increase in heart failure risk (95% CI 2.21–3.34; p<0.001). The mean c-statistics for models including the trajectory risk score and single-occasion traits were 0.87 (95% CI 0.83–0.91) and 0.83 (95% CI 0.80–0.86), respectively.

Conclusion

Incorporating multi-system trajectories reflective of the aging process may add incremental information to heart failure risk assessment when compared to using single-occasion traits.

Introduction

The prevalence of heart failure (HF) among adults 20 years of age and older in the United States is estimated to be 6.2 million based on data accrued from the National Health and Nutrition Examination Survey between 2013–2016 [1]. From 2012–2030, the prevalence of HF is expected to increase by 46% [1]. Despite advancements in treatment, the diagnosis of HF continues to portend an extremely poor prognosis [13]. Given the high mortality rate once the diagnosis is made, a better understanding of the clinical syndrome before the manifestation of the disease is imperative. The increasing incidence of HF has been attributed to the aging population [3]. However, it remains unclear why the extent of cardiac dysfunction, a potential surrogate of cardiac aging, does not consistently correlate with the clinical severity of the condition (in terms of symptoms); a potential explanation could be the fact that HF is a multi-system syndrome with several non-cardiac predictors [4].

Prior studies have focused on the relation between a set of traits measured at one point in time and the risk of HF [424]. Each organ system follows a different trajectory of aging and decline, which may have a variable impact on HF risk [2527]. However, to our knowledge, data are lacking on whether the decline of the function of multiple organ systems over time is associated with HF risk. Moreover, it is not clear whether the trajectories of the function of various organ systems conjointly could provide better discrimination ability with regards to HF risk compared to single-occasion measurements of the same organ systems.

In this investigation, we hypothesized that the distinct long-term trajectory patterns of multiple organ systems are individually related to HF risk. Moreover, we postulated that collective-trajectory profiles of key organ systems might provide better discrimination information in terms of HF risk compared to single-occasion measurements. We evaluated these hypotheses using data from the Framingham Heart Study (FHS) Offspring cohort. We chose 12 distinct organ systems, previously related to cardiovascular disease (CVD), CVD mortality, and all-cause mortality, as a comprehensive representation of the body’s aging [25].

Methods

Data for the Framingham Offspring Study examinations can be obtained from the National Heart, Lung, and Blood Institute’s data repository at https://biolincc.nhlbi.nih.gov/studies/.

Study sample

Prior investigations have published on the inclusion and exclusion criteria of the FHS cohorts and its study design [28]. Of the 3,799 Offspring cohort participants who attended the fifth examination cycle (1991–1995), we excluded 250 participants who did not attend at least one other examination after the fifth, i.e., the sixth (1995–1998), seventh (1998–2001), or eighth (2005–2008) examination cycles. Further, we excluded participants with prevalent HF at their most recent examination attendance (n = 111), had missing baseline covariates (n = 11), or had no follow-up information (n = 15), resulting in a final sample size of 3,412 participants. We applied the eligibility criteria at the start of the fifth examination cycle since the purpose of this investigation was to examine mid- to late-life biological system trajectories in relation to HF risk, and serial measurements for the selected traits were available starting at the fifth examination cycle.

All participants provided written informed consent. Approval for this investigation was obtained from the Institutional Review Board at the Boston University Medical Center.

Organ system traits

For the present investigation, we used the following traits and their corresponding representative body organ system: estimated glomerular filtration rate, eGFR (kidney function); forced vital capacity, FVC, and the ratio of forced expiratory volume in one second and FVC, FEV1/FVC (pulmonary function); gait time (neuromotor function); grip strength (muscular function); left ventricular mass indexed by height, LVMI and heart rate, HR (cardiac function); pulse pressure, PP (vascular function); the ratio of total cholesterol and high-density lipoprotein cholesterol, TC/HDL (cholesterol metabolism); C-reactive protein, CRP (inflammation); body mass index, BMI (adiposity); and hemoglobin A1c, HbA1c (glucose homeostasis).

We selected these unique traits because they have previously been associated with CVD, CVD-related death, and all-cause mortality [25]. These traits provide a comprehensive representation of multiple organ systems and the metabolism of individuals [25].

Outcome of interest

The primary outcome of this investigation was incident HF. Participants were followed for the development of HF from their most recent examination cycle attendance (sixth, seventh or eighth) through December 31, 2017. Surveillance of the participants for HF outcomes included collecting their health history, reviewing FHS visits and medical records from outpatient visits, and hospitalizations related to HF. An adjudication panel consisting of three FHS investigators evaluated all outcomes. The diagnosis of HF was based on the previously established Framingham criteria for congestive HF requiring the presence of two major criteria or one major and two minor criteria [29].

Statistical analysis

Descriptive statistics for the sample characteristics were reported from the most recent examination attended by the participants.

For each trait, we created two exposure variables: one representing a single occasion measurement and another representing the trajectory-based longitudinal pattern. The most recent measurements collected between 1991 and 2008 were used as the single time point exposures.

We estimated sex-specific age-adjusted Spearman partial correlation coefficients among the single-occasion traits.

Association of single-occasion traits with HF risk

After confirming that the proportional hazards assumption was met, we estimated Cox proportional hazards regression models to relate each single-occasion trait to the time to incident HF (separate model for each of the 12 traits). Models were adjusted for age, sex, smoking, BMI (or weight for traits indexed by height), antihypertensive treatment, diabetes, TC/HDL, and SBP (except when evaluating PP to avoid multicollinearity). The most recent measurements collected between 1991 and 2008 were used as the single time point exposures.

Group-based trajectories of traits

We created group-based model trajectories for each of the 12 traits, using at least two measurements for each trait from each participant (one from the fifth examination cycle and all other available measurements from exams six through eight). Group-based trajectory modeling assumes unique subgroups within the sample, each with its underlying trajectory. We used SAS Proc Traj, a finite mixture model procedure using maximum likelihood, to identify groups of participants who follow a similar progression of traits [30]. In determining the group-based trajectory model that best fits each trait, we followed an algorithm proposed by Jones and Nagin [31, 32]. We decided a priori to create models ranging from one to five groups with each set to a quadratic order for comparison purposes. All traits were continuous and assumed to have a censored normal distribution. Using the Bayesian Information Criterion (BIC), we compared the goodness-of-fit between the models in a stepwise manner to establish our decision on the number of groups with the aim to have the estimated membership in each group be at least 5%. Next, we fitted the model by varying the polynomial degrees (i.e., quartic, cubic, quadratic, and linear) to determine the shape of each group’s trajectory over time. The final fitted model computes an individual’s posterior probability of group membership and assigns the person to the group with the highest posterior probability. To confirm that our model adequately fit the data, we performed a diagnostic assessment by verifying that the average posterior probabilities of the group an individual was assigned to exceeded 0.7. Once we determined the best model fit and confirmed the model diagnostics for each trait, each participant was classified into a trajectory group based on their estimated highest posterior probability of group membership. CRP was winsorized at the 95th percentile by exam and sex, and was then log-transformed before model fitting. In fitting the group-based trajectory models, all traits were adjusted for age at examination cycle five with the assumption that it is a risk factor affecting the probability of group membership. Depending on the trait, we further included a time-varying covariate in estimating the trajectory groups: height (for models including FEV1/FVC and FVC), diabetes medication use (for HbA1c), beta-blocker use (for HR), and lipid-lowering medication use (for TC/HDL). As an exploratory analysis, we excluded participants with prevalent MI and/or previous cardiac surgery when creating the group-based trajectories. Moreover, we also performed sex-specific exploratory analysis of the single-occasion and group trajectories association in relation to HF risk (with corresponding sex-specific trajectory profiles).

Association of group-based trajectories with HF risk

After estimating the trajectory profile group for each trait using group-based modelling, we created binary variables for traits resulting in two distinct trajectory profiles (1 = ‘best’ and 2 = ‘worst’) and 3-level variables for traits resulting in three distinct trajectory profiles (1 = ‘best’, 2 = ‘intermediate’, and 3 = ‘worst’). We related the derived variables (‘best’ serving as the referent group) to HF risk using Cox proportional hazards regression models (separate model for each trait): (1) adjusting for age and sex and (2) adjusting for age, sex, smoking status, BMI (or weight for traits indexed by height), antihypertensive treatment, diabetes status, TC/HDL, and SBP (except when evaluating PP to preclude multicollinearity). The proportional hazards assumption was met for all models. Follow up for the development of HF started after the most recent examination for each participant.

Model discrimination for single-occasion models vs. group-based trajectory models

We compared the discrimination ability of models including single-occasion variables versus those including group-based trajectory variables for each trait, using 500 bootstrapped samples by estimating the mean Harrell’s c-statistic index (with 95% CI). In each bootstrap, a c-statistic index was obtained for three Cox proportional hazards models: (1) covariates-only model (S0), (2) covariates plus the single-occasion variable (S1), and (3) covariates plus the group-based trajectory variable (S2). In addition, we computed the mean change in c-statistic between a model with only covariates and a model including covariates plus the single-occasion variable (ΔS1 = S1 − S0). We also calculated the mean change in c-statistic between a model with only covariates and a model including covariates plus the group-based trajectory variable (ΔS2 = S2 − S0). The set of covariates (i.e., age, sex, smoking, BMI, etc.) in these models were consistent with the prior models. The bootstrap process was performed using the SAS macro %boot.

Derivation of the trajectory risk score

We used a forward selection algorithm in a Cox proportional hazards model that incorporated all traits and selected variables using the likelihood ratio test. Multicollinearity was inspected among the group-based trajectory variables prior to the forward selection process. The proportional hazards assumption was evaluated and confirmed at each step of the variable selection. We then created a composite trajectory risk score using the variables that were selected from the forward algorithm as follows: we calculated the sum of the scores previously assigned to the trajectory variables (1 = ‘best,’ 2 = ‘intermediate,’ 3 = ‘worst’), weighted by the corresponding regression coefficients. Then, we related the derived composite trajectory risk score to HF risk using a Cox proportional hazards model, adjusted for age and sex. Finally, we created tertiles of the composite trajectory risk score and created unadjusted cumulative incidence curves.

Model discrimination for multiple single-occasion traits model vs. trajectory risk score model

To compare the discrimination ability of a model including the trajectory risk score to a model including multiple single-occasion variables, we used a process similar to what was described above for discriminating single-occasion models and group-based trajectory models. The c-statistic index obtained for the three Cox proportional hazards models include: (1) age-sex only model (M0); (2) age, sex, plus trajectory risk score (M1); and (3) age, sex, plus multiple single-occasion variables (M2). The traits selected for the multiple single-occasion model are based on the aforementioned forward selection process. As described above, we calculated the mean change in c-statistic from the age- and sex-adjusted model to the model, including age, sex, and the trajectory risk score (ΔM1 = M1 − M0), as well as the mean change in c-statistic from the age- sex-adjusted model to the model, including age, sex, and the multiple single-occasion variables (ΔM2 = M2 − M0).

Statistical significance was evaluated using a two-sided P<0.05 unless otherwise stated. All analyses were performed using SAS version 9.4 (Cary, NC).

Results

Study sample

The characteristics of our study sample, stratified by sex, are shown in Table 1. Our sample had an average age of 67 years, and slightly over half were women (54%). On average, participants were overweight, and approximately half were on antihypertensive treatment. Details on the average and maximum number of measurements per participant for each trait are included in S3 Table.

thumbnail
Table 1. Baseline characteristics of study sample, stratified by sex.

https://doi.org/10.1371/journal.pone.0268576.t001

Correlation among traits

Age-adjusted and sex-stratified Spearman correlations among the 12 surrogate traits are illustrated in Fig 1. In men, we observed moderate correlations of FVC with grip strength, and of BMI with LVMI (Fig 1). In women, there were also modest correlations of BMI with TC/HDL, LVMI, and CRP.

thumbnail
Fig 1. Sex-specific spearman correlations, adjusted for age.

Correlations among single-occasion traits using Spearman partial correlation coefficients, adjusting for age, are depicted. Men (n = 1179) are shown in the lower triangle, and women (n = 1456) are shown in the upper triangle. Deeper red = stronger positive correlation; deeper blue = stronger negative correlation. BMI = body mass index; CRP = C-reactive protein; eGFR = estimated glomerular filtration rate; FEV1 = forced expiratory volume; FVC = forced vital capacity; HbA1c = hemoglobin A1c; HDL = high-density lipoprotein; HR = heart rate; LVMI = left ventricular mass index; PP = pulse pressure; TC = total cholesterol.

https://doi.org/10.1371/journal.pone.0268576.g001

Association of single-occasion traits with HF risk

We observed 276 new-onset HF events over a median follow-up period of 10.3 years. Adjusting for age, sex, smoking status, BMI (or weight for traits indexed by height), antihypertensive treatment, diabetes status, TC/HDL, and SBP (except when evaluating PP), we observed positive associations of HbA1c, BMI, PP, CRP, HR, TC/HDL, LVMI, and gait with HF risk. We also observed an inverse association of eGFR and FVC with HF risk. (Table 2)

thumbnail
Table 2. Associations of single-occasion traits and group trajectories with HF risk.

https://doi.org/10.1371/journal.pone.0268576.t002

Group-based trajectories of traits

The trajectories are shown in Fig 2. For most traits, the trajectory patterns were parallel over time. For LVMI and gait time, the ‘worst’-risk group had a steeper increase while the ‘worst’-risk group trajectory for TC/HDL started relatively high but had a steeper decrease over time. For HbA1c, the ‘worst’-risk group had a non-linear pattern. In exploratory analysis, when excluding participants with prevalent MI and/or previous cardiac surgery (n = 169, 5%) the group-based trajectories remained consistent (S1 Fig and S2 Table). Moreover, the sex-specific exploratory analysis of the single-occasion and group trajectories association in relation to HF risk is shown in S4 and S5 Tables and the corresponding trajectory profiles are in S2 Fig.

thumbnail
Fig 2. Group-based trajectories for traits.

Group-based trajectories for the twelve surrogate traits are shown. For traits resulting in two distinct trajectory profiles, green = ‘best’ and red = ‘worst,’ and if three different profiles emerged, green = ‘best,’ blue = ‘intermediate,’ and red = ‘worst.’ The x-axis represents the number of years since the fifth examination cycle. BMI = body mass index; CRP = C-reactive protein; eGFR = estimated glomerular filtration rate; LVMI = left ventricular mass index; PP = pulse pressure.

https://doi.org/10.1371/journal.pone.0268576.g002

Association of group-based trajectories with HF risk

Adjusting for the covariates noted above, participants with the ‘worst’ trajectory for PP, LVMI, HR, FVC, HbA1c, BMI, CRP, and gait were at higher risk for HF. In contrast to the single-occasion models, we observed a significant positive relation between the trajectory groups for FEV1/FVC and HF risk (Table 2).

The mean change in c-statistic (from a base model) was higher for most traits when using the model with the group-based trajectories than the model with the single-occasion variables (S1 Table).

Association of trajectory risk score with HF

Based on the forward selection algorithm, the 7 traits selected and used to derive the trajectory risk score were CRP, FVC, gait time, HbA1c, HR, LVMI, and PP. A one-unit increase in the composite trajectory risk score, adjusting for age and sex, was associated with a 2.72-fold increased risk of HF (95% CI [2.21, 3.34]; p<0.001). Moreover, adding the trajectory risk score to a model including age and sex resulted in a higher c-statistic compared to the c-statistic derived from a multiple single-occasion model that included the same 7 single-occasion traits used to derive the risk score, plus age and sex (Table 3).

thumbnail
Table 3. C-statistics for models including the trajectory risk score and single-occasion traits.

https://doi.org/10.1371/journal.pone.0268576.t003

The cumulative incidence of HF among participants with moderate risk scores (2nd tertile) was 2.3% and among those with high-risk scores (highest tertile) was approximately 16% (Fig 3).

thumbnail
Fig 3. Cumulative incidence of HF by tertiles of the trajectory risk score.

The cumulative incidence of heart failure by tertiles of the trajectory risk score is illustrated. Trajectory risk scores are segregated into 3 groups: high = red, moderate = blue, and low = green. The number of participants at risk, every two years since the most recent exam, is included along the x-axis. The median overall follow up for participants with a trajectory risk score was 10.7 years. Log-rank p<0.001.

https://doi.org/10.1371/journal.pone.0268576.g003

Discussion

Principal findings

First, we observed relations of LVMI, PP, CRP, BMI, HbA1c, HR, FVC and gait time with HF risk using either single-occasion measurements or trajectory groups [522]. Notably, gait time was positively related to HF risk in both the single-occasion model and the group trajectory models; limited prior data exist on this relation among people free of HF. Second, we did not observe a relation between grip strength and HF risk using either model. Third, we observed a positive relation between the FEV1/FVC and HF risk using the trajectory model, but not as a single occasion measurement. Fourth, we observed positive relations of TC/HDL and eGFR with HF risk when the single-occasion model was used, but not with trajectory groups. Fifth, using the forward selection algorithm, seven traits were derived to create the trajectory risk score: CRP, FVC, gait time, HbA1c, HR, LVMI, and PP. Sixth, we observed that a one-unit increase in the composite multi-system trajectory risk score was associated with an approximately three-fold increase in HF risk adjusting for standard HF risk factors. Seventh, adding the trajectory risk score to a model including age and sex resulted in a higher increase in the c-statistic compared to the addition of single-occasion variables to an age-, and sex-adjusted model. Finally, participants with values in the highest tertile of the trajectory score had an absolute HF risk that was substantially elevated compared to the lower two tertiles.

Comparison with the literature

Prior studies have demonstrated a link between the aging of the population and an increased incidence of HF [3, 5]. However, it has become apparent that not all aging is the same [2527]. Most prior reports have focused on the relations between traits measured on a single occasion and HF risk [424]. Our investigation demonstrated that most single-occasion traits are weakly correlated with each other, indicating that the function of individual organ systems may follow different trajectories during the aging process.

Previous studies have reported on the relations of single-occasion traits including TC/HDL, PP, HR, CRP, BMI, LVMI, and HbA1c, with HF risk, consistent with our study findings [516]; this supports the established associations of a pro-inflammatory state (CRP), sympathetic hyperactivity (HR), coronary heart disease (TC/HDL, PP, BMI, LVMI, HbA1c), with the risk of incident HF. We also confirmed an inverse relation of FVC with HF risk [19, 20]. The association of lung dysfunction with HF risk has been attributed to altered pre-load and after-load conditions on the heart, an inflammatory state, and increased oxidative stress on the myocardium. Furthermore, our findings supported the previously noted inverse relation of eGFR with HF risk [17, 18]. Kidney dysfunction has been associated with an increased risk of HF due to an altered renin-angiotensin-aldosterone system which is integral to the body’s management of blood volume and blood pressure [18]. Although grip strength has been linked to cardiovascular fitness and risk of HF, we did not observe an association between grip strength and incident HF; prior data for this trait are limited, but two previous investigations reported an inverse relation [23, 24]. Gait time has been associated with mortality, cardiovascular fitness, and recurrent HF in individuals with prevalent HF; we extended these findings by confirming a direct association between gait time and future HF risk in individuals without HF, which, to our knowledge, had only been reported in one previous study [22, 33, 34].

Given that organ systems have been shown to follow different aging patterns, we created trajectories of the selected traits to provide additional insight into future disease risk [2527, 32]. To our knowledge, the association of organ system trajectories to incident HF has not been previously investigated. In our analysis, we leveraged serial trait measurements to create trajectory groups that were associated with a higher risk of HF. In contrast to the single-occasion trait model that did not show an association between FEV1/FVC and incident HF, the ‘worst’ subgroup trajectory for FEV1/FVC was associated with higher HF risk. Based on our results, lung dysfunction may only be clinically relevant to HF risk in the subgroup of individuals with the most drastic decline in lung function over time [21]. Similar to the single occasion trait model, we did not observe an association between grip strength trajectory and incident HF, indicating that the trajectory for grip strength may not be an ideal marker for future HF risk. Additionally, the ‘worst’ TC/HDL trajectory did not show a significant relation to incident HF possibly because there was improvement in the participants’ cholesterol profiles over time. Lastly, we did not observe an association of the eGFR trajectory with HF risk, suggesting that the most recent measurement of eGFR may be a better indicator of future HF risk than its trajectory.

Our investigation demonstrated that individuals with suboptimal trajectory scores across multiple systems were at the highest risk of HF. The incremental discrimination ability of models including trajectory variables (above and beyond that of a model including age and sex) was greater than the ability of a model including single-occasion traits. Furthermore, we observed a positive association between the trajectory score and HF risk. Our results suggest that trajectories may provide a more powerful tool in discriminating the most at-risk individuals for future HF. Ultimately, we hope that trajectory scores may provide clinicians with better insight into ascertaining which individuals are at highest risk of future HF. This may allow for earlier intervention and targeted therapy that would hopefully alter the disease course leading to clinical HF. Trajectory scores may provide clinicians with a more comprehensive understanding of a person’s risk of future HF by incorporating the aging of multiple organ systems over time into one model. Future studies are needed to assess whether incorporating a trajectory score into an individual’s assessment alters disease management and reduces risk of future HF.

Study strengths and limitations

The participants of this investigation were from a large, community-based cohort, followed over a prolonged period. We chose 12 traits to provide a comprehensive representation of the body’s functions and metabolism. There are several limitations of this investigation. First, there was a potential for survival bias given that participants who attended the fifth examination cycle and at least one other examination cycle after 1995 were eligible for inclusion. Second, the participants consisted mainly of white adults of European ancestry, limiting the generalizability of our findings to other races and ethnicities. Third, this was a prospective cohort investigation, and therefore, cause and effect cannot be inferred from the associations we observed.

Conclusion

Our findings support the incremental predictive value of long-term repeated measurements of multiple organ systems in evaluating HF risk, which may ultimately provide better insight into ascertaining which individuals are at greatest risk of future HF. Additional studies are needed to determine whether earlier interventions targeted towards individuals with poor trajectory risk profiles may prevent future HF.

Supporting information

S1 Fig. Group trajectories for the surrogate traits.

https://doi.org/10.1371/journal.pone.0268576.s001

(DOCX)

S2 Fig. Sex-specific group trajectories for surrogate traits.

https://doi.org/10.1371/journal.pone.0268576.s002

(DOCX)

S1 Table. C-statistics for single trait and group-based trajectory models.

Bias-corrected c-statistics are reported based on 500 bootstrap samples using the Cox proportional hazards model. Model 0 includes the covariates: age, sex, smoking status, BMI = body mass index (or weight for traits indexed by height), antihypertensive treatment, diabetes status, TC/HDL = total cholesterol/high-density lipoprotein, and systolic blood pressure (except when evaluating pulse pressure). Model 1 includes the covariates plus the single occasion variable. Model 2 includes the covariates plus the group-based trajectory variable. ΔS1 = mean change in c-statistic between Model 0 and Model 1. ΔS2 = mean change in c-statistic between Model 0 and Model 2. CRP = C-reactive protein; eGFR = estimated glomerular filtration rate; FEV1 = forced expiratory volume; FVC = forced vital capacity; HbA1c = hemoglobin A1c; HR = heart rate; LVMI = left ventricular mass index.

https://doi.org/10.1371/journal.pone.0268576.s003

(DOCX)

S2 Table. Associations of single-occasion traits and group trajectories with HF risk (excluding prevalent MI and/or previous cardiac surgery).

https://doi.org/10.1371/journal.pone.0268576.s004

(DOCX)

S4 Table. Associations of single-occasion traits and group trajectories with HF risk for men.

https://doi.org/10.1371/journal.pone.0268576.s006

(DOCX)

S5 Table. Associations of single-occasion traits and group trajectories with HF risk for women.

https://doi.org/10.1371/journal.pone.0268576.s007

(DOCX)

Acknowledgments

We want to acknowledge the participants of the Framingham Heart Study whose selfless participation in research has allowed us to conduct not only this investigation but also many others in hopes of bringing us closer to a more comprehensive understanding of cardiovascular disease.

References

  1. 1. Virani SS, Alonso A, Benjamin EJ, Bittencourt MS, Callaway CW, Carson AP, et al. Heart Disease and Stroke Statistics-2020 Update: A Report From the American Heart Association. Circulation. 2020;141(9):e139–e596. pmid:31992061
  2. 2. Goldberg RJ, Ciampa J, Lessard D, Meyer TE, Spencer FA. Long-term survival after heart failure: a contemporary population-based perspective. Arch Intern Med. 2007;167(5):490–6. pmid:17353497
  3. 3. Ho KK, Pinsky JL, Kannel WB, Levy D. The epidemiology of heart failure: the Framingham Study. J Am Coll Cardiol. 1993;22(4 Suppl A):6A–13A. pmid:8376698
  4. 4. Lam CS, Lyass A, Kraigher-Krainer E, Massaro JM, Lee DS, Ho JE, et al. Cardiac dysfunction and noncardiac dysfunction as precursors of heart failure with reduced and preserved ejection fraction in the community. Circulation. 2011;124(1):24–30. pmid:21670229
  5. 5. Gottdiener JS, Arnold AM, Aurigemma GP, Polak JF, Tracy RP, Kitzman DW, et al. Predictors of congestive heart failure in the elderly: the Cardiovascular Health Study. J Am Coll Cardiol. 2000;35(6):1628–37. pmid:10807470
  6. 6. de Simone G, Gottdiener JS, Chinali M, Maurer MS. Left ventricular mass predicts heart failure not related to previous myocardial infarction: the Cardiovascular Health Study. Eur Heart J. 2008;29(6):741–7. pmid:18204091
  7. 7. Drazner MH, Rame JE, Marino EK, Gottdiener JS, Kitzman DW, Gardin JM, et al. Increased left ventricular mass is a risk factor for the development of a depressed left ventricular ejection fraction within five years: the Cardiovascular Health Study. J Am Coll Cardiol. 2004;43(12):2207–15. pmid:15193681
  8. 8. Chae CU, Pfeffer MA, Glynn RJ, Mitchell GF, Taylor JO, Hennekens CH. Increased pulse pressure and risk of heart failure in the elderly. JAMA. 1999;281(7):634–9. pmid:10029125
  9. 9. Nanchen D, Stott DJ, Gussekloo J, Mooijaart SP, Westendorp RG, Jukema JW, et al. Resting heart rate and incident heart failure and cardiovascular mortality in older adults: role of inflammation and endothelial dysfunction: the PROSPER study. Eur J Heart Fail. 2013;15(5):581–8. pmid:23250912
  10. 10. Kardys I, Knetsch AM, Bleumink GS, Deckers JW, Hofman A, Stricker BH, et al. C-reactive protein and risk of heart failure. The Rotterdam Study. Am Heart J. 2006;152(3):514–20. pmid:16923423
  11. 11. Vasan RS, Sullivan LM, Roubenoff R, Dinarello CA, Harris T, Benjamin EJ, et al. Inflammatory markers and risk of heart failure in elderly subjects without prior myocardial infarction: the Framingham Heart Study. Circulation. 2003;107(11):1486–91. pmid:12654604
  12. 12. Ndumele CE, Matsushita K, Lazo M, Bello N, Blumenthal RS, Gerstenblith G, et al. Obesity and Subtypes of Incident Cardiovascular Disease. J Am Heart Assoc. 2016;5(8). pmid:27468925
  13. 13. Iribarren C, Karter AJ, Go AS, Ferrara A, Liu JY, Sidney S, et al. Glycemic control and heart failure among adult patients with diabetes. Circulation. 2001;103(22):2668–73. pmid:11390335
  14. 14. Matsushita K, Blecker S, Pazin-Filho A, Bertoni A, Chang PP, Coresh J, et al. The association of hemoglobin a1c with incident heart failure among people without diabetes: the atherosclerosis risk in communities study. Diabetes. 2010;59(8):2020–6. pmid:20484138
  15. 15. Kannel WB, Ho K, Thom T. Changing epidemiological features of cardiac failure. Br Heart J. 1994;72(2 Suppl):S3–9. pmid:7946754
  16. 16. Velagaleti RS, Massaro J, Vasan RS, Robins SJ, Kannel WB, Levy D. Relations of lipid concentrations to heart failure incidence: the Framingham Heart Study. Circulation. 2009;120(23):2345–51. pmid:19933936
  17. 17. Shara NM, Wang H, Mete M, Al-Balha YR, Azalddin N, Lee ET, et al. Estimated GFR and incident cardiovascular disease events in American Indians: the Strong Heart Study. Am J Kidney Dis. 2012;60(5):795–803. pmid:22841159
  18. 18. Kottgen A, Russell SD, Loehr LR, Crainiceanu CM, Rosamond WD, Chang PP, et al. Reduced kidney function as a risk factor for incident heart failure: the atherosclerosis risk in communities (ARIC) study. J Am Soc Nephrol. 2007;18(4):1307–15. pmid:17344421
  19. 19. Kannel WB, Seidman JM, Fercho W, Castelli WP. Vital capacity and congestive heart failure. The Framingham study. Circulation. 1974;49(6):1160–6. pmid:4275471
  20. 20. Jankowich M, Elston B, Liu Q, Abbasi S, Wu WC, Blackshear C, et al. Restrictive Spirometry Pattern, Cardiac Structure and Function, and Incident Heart Failure in African Americans. The Jackson Heart Study. Ann Am Thorac Soc. 2018;15(10):1186–96. pmid:30011374
  21. 21. Georgiopoulou VV, Kalogeropoulos AP, Psaty BM, Rodondi N, Bauer DC, Butler AB, et al. Lung function and risk for heart failure among older adults: the Health ABC Study. Am J Med. 2011;124(4):334–41. pmid:21435424
  22. 22. Chaudhry SI, McAvay G, Chen S, Whitson H, Newman AB, Krumholz HM, et al. Risk factors for hospital admission among older persons with newly diagnosed heart failure: findings from the Cardiovascular Health Study. J Am Coll Cardiol. 2013;61(6):635–42. pmid:23391194
  23. 23. Sillars A, Celis-Morales CA, Ho FK, Petermann F, Welsh P, Iliodromiti S, et al. Association of Fitness and Grip Strength With Heart Failure: Findings From the UK Biobank Population-Based Study. Mayo Clin Proc. 2019;94(11):2230–40. pmid:31685151
  24. 24. Laukkanen JA, Khan H, Lavie CJ, Voutilainen A, Kurl S, Jae SY, et al. Inverse Association of Handgrip Strength With Risk of Heart Failure. Mayo Clin Proc. 2021;96(6):1490–9. pmid:33775422
  25. 25. Niiranen TJ, Enserro DM, Larson MG, Vasan RS. Multisystem Trajectories Over the Adult Life Course and Relations to Cardiovascular Disease and Death. J Gerontol A Biol Sci Med Sci. 2019;74(11):1778–85. pmid:30358808
  26. 26. Belsky DW, Caspi A, Houts R, Cohen HJ, Corcoran DL, Danese A, et al. Quantification of biological aging in young adults. Proc Natl Acad Sci U S A. 2015;112(30):E4104–10. pmid:26150497
  27. 27. Newman AB, Arnold AM, Sachs MC, Ives DG, Cushman M, Strotmeyer ES, et al. Long-term function in an older cohort—the cardiovascular health study. J Am Geriatr Soc. 2009;57(3):432–40. pmid:19187412
  28. 28. Feinleib M, Kannel WB, Garrison RJ, McNamara PM, Castelli WP. The Framingham Offspring Study. Design and preliminary data. Prev Med. 1975;4(4):518–25. pmid:1208363
  29. 29. McKee PA, Castelli WP, McNamara PM, Kannel WB. The natural history of congestive heart failure: the Framingham study. N Engl J Med. 1971;285(26):1441–6. pmid:5122894
  30. 30. Jones BL, Nagin DS. Advances in Group-Based Trajectory Modeling and an SAS Procedure for Estimating Them. Sociological Methods & Research. 2007;35(4):542–71.
  31. 31. Jones BL, Nagin DS, Roeder K. A SAS Procedure Based on Mixture Models for Estimating Developmental Trajectories. Sociological Methods & Research. 2001;29(3):374–93.
  32. 32. Nagin D. Group-based modeling of development. Cambridge, Mass.; London: Harvard University Press; 2005. 201 p. p.
  33. 33. Studenski S, Perera S, Patel K, Rosano C, Faulkner K, Inzitari M, et al. Gait speed and survival in older adults. JAMA. 2011;305(1):50–8. pmid:21205966
  34. 34. Pulignano G, Del Sindaco D, Di Lenarda A, Alunni G, Senni M, Tarantini L, et al. Incremental Value of Gait Speed in Predicting Prognosis of Older Adults With Heart Failure: Insights From the IMAGE-HF Study. JACC Heart Fail. 2016;4(4):289–98. pmid:26970831