Figures
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.
Citation: Guardino CE, Pan S, Vasan RS, Xanthakis V (2022) Multi-system trajectories and the incidence of heart failure in the Framingham Offspring Study. PLoS ONE 17(5): e0268576. https://doi.org/10.1371/journal.pone.0268576
Editor: Spencer D. Proctor, University of Alberta, CANADA
Received: November 15, 2021; Accepted: May 2, 2022; Published: May 26, 2022
Copyright: © 2022 Guardino et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All data is found at the Biologic Specimen and Data Repository Information Coordinating Center at https://biolincc.nhlbi.nih.gov/studies/framcohort/ (Accession number: HLB00050021a).
Funding: This work was supported by National Heart, Lung, and Blood Institute, Bethesda, MD, contracts NO1‐HC‐25195, HHSN268201500001I, and 75N92019D00031. Dr. Ramachandran was also funded by the Evans Foundation and Jay and Louise Coffman Endowment at Boston University, Boston, MA. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
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 [1–3]. 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 [4–24]. Each organ system follows a different trajectory of aging and decline, which may have a variable impact on HF risk [25–27]. 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.
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.
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.
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)
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.
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.
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).
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).
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.
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 [5–22]. 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 [25–27]. Most prior reports have focused on the relations between traits measured on a single occasion and HF risk [4–24]. 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 [5–16]; 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 [25–27, 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. Kannel WB, Ho K, Thom T. Changing epidemiological features of cardiac failure. Br Heart J. 1994;72(2 Suppl):S3–9. pmid:7946754
- 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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.
Nagin D. Group-based modeling of development. Cambridge, Mass.; London: Harvard University Press; 2005. 201 p. p.
- 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. 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