Skip to main content

Multi-trajectories of lipid indices with incident cardiovascular disease, heart failure, and all-cause mortality: 23 years follow-up of two US cohort studies

Abstract

Background

Understanding the distinct patterns (trajectories) of variation in blood lipid levels before diagnosing cardiovascular disease (CVD) might carry important implications for improving disease prevention or treatment.

Methods

We investigated 14,373 participants (45.5% men) aged 45–84 from two large US prospective cohort studies with a median of 23 years follow-up. First, we jointly estimated developmental trajectories of lipid indices, including low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and triglyceride (TG) concentrations using group-based multi-trajectory modeling. Then, the association of identified multi-trajectories with incident CVD, heart failure, and all-cause mortality were examined using Cox proportional hazard model.

Results

Seven distinct multi-trajectories were identified. The majority of participants (approximately 80%) exhibited decreasing LDL-C but rising TG levels and relatively stable HDL-C levels. Compared to the individuals with healthy and stable LDL-C, HDL-C, and TG levels, those in other groups were at significant risk of incident CVD after adjusting for other conventional risk factors. Individuals with the highest but decreasing LDL-C and borderline high and rising TG levels over time were at the highest risk than those in other groups with a 2.22-fold risk of CVD. Also, those with the highest and increased triglyceride levels over time, over optimal and decreasing LDL-C levels, and the lowest HDL-C profile had a nearly 1.84 times CVD risk. Even individuals in the multi-trajectory group with the highest HDL-C, optimal LDL-C, and optimal TG levels had a significant risk (HR, 1.45; 95% CI 1.02–2.08). Furthermore, only those with the highest HDL-C profile increased the risk of heart failure by 1.5-fold (95% CI 1.07–2.06).

Conclusions

The trajectories and risk of CVD identified in this study demonstrated that despite a decline in LDL-C over time, a significant amount of residual risk for CVD remains. These findings suggest the impact of the increasing trend of TG on CVD risk and emphasize the importance of assessing the lipid levels at each visit and undertaking potential interventions that lower triglyceride concentrations to reduce the residual risk of CVD, even among those with the optimal LDL-C level.

Introduction

Cardiovascular diseases (CVDs) remain the leading cause of death and a significant cause of disability worldwide, accounting for 17.9 million deaths per year or 31% of all global deaths [1]. Unfavorable lipid indices, represented by increased serum concentrations of total cholesterol (TC), triglyceride (TG), low‐density lipoprotein cholesterol (LDL‐C), and decreased high-density lipoprotein cholesterol (HDL‐C), are conventionally considered to play an essential role in the development and progression of cardiovascular disease [2].

Some existing evidence suggests age-related changes in the TC, LDL-C, and TG so that they increase up to middle age and then decrease [3,4,5,6]. However, some epidemiological evidence also suggests that long-term exposure to even moderately raised cholesterol levels is associated with CVD in the future [7, 8]. Therefore, understanding the distinct patterns of variation in blood lipid levels before the diagnosis of CVD might carry important implications for improving disease prevention or treatment.

Most of the evidence on the trends in lipid indices comes from comparing observed average lipid levels in sequential cross-sectional surveys, and only a few of these studies have investigated trends in the lipid indices within the same population [3, 9,10,11,12]. However, most of these investigations did not consider the correlation between lipids and assessed the trajectory of lipid indices separately [3, 10,11,12]. Thus, evidence exists on lipid indices’ independent effect, but little is known about their combined impact on CVD risk. Moreover, identifying distinct longitudinal patterns of different lipids might help better understand the variation of lipid indices over time and facilitate targeted cardiovascular prevention programs [9].

Using pooled data from two US cohort studies with a median of 23 years follow-up and repeated measurements of lipid indices, we sought to (1) identify longitudinal multi-trajectories of LDL‐C, HDL‐C, and TG over the adult life course and (2) examine their associations with subsequent risks of incident CVD, heart failure, and all-cause mortality later in life. To capture the overlap between the developments of LDL‐C, HDL‐C, and TG, we estimated the joint developmental trajectories of these lipids. This approach allowed us to identify subgroups that share common patterns of change over time.

Materials and methods

Study population

The present study was based on data from two large, community-based, prospective cohort studies sponsored by the National Heart, Lung, and Blood Institute (NHLBI), the Atherosclerosis Risk in Communities (ARIC) study, and the Multi-Ethnic Study of Atherosclerosis (MESA) study. Details of the method and design of each study have been previously published [13, 14]. In brief, the ARIC Study is a prospective cohort study of 15,792 individuals 45 to 64 years of age recruited from 1987 to 1989 from 4 US communities [13]. The MESA Study recruited 6814 individuals 45 to 84 years of age free of clinical CVD at baseline during 2000–2002, from 6 US communities [14]. So far, both studies have five subsequent examination cycles (a sixth follow-up visit is currently underway), and all participants provided written informed consent. The studies’ website contains details of all the available data through a fully searchable data dictionary [15, 16]. The process of pooling data was carried out using the guidelines developed by Maelstrom Research for rigorous retrospective data harmonization [17].

The current study consists of two samples according to the analyses performed. The flowchart of the participant selection process is shown in Fig. 1. The first sample (sample A) was used to create the multi-trajectory groups for the lipid indices with the objective of a later link to CVD events and all-cause mortality. Since the trajectory approach requires at least three unique time points and becomes more precise with additional time points, this sample was restricted to participants attending all 4-examination cycles. Besides, participants with known CVD and renal failure and missing values for each lipid at baseline were excluded from this sample (Fig. 1, sample A). For incident CVD events from the second exam through the fourth exam, all data until the date of the first incident CVD have been included. Data from examination cycle five was not used for prospective analyses described as follows to ensure sufficient follow-up time for these analyses.

Fig. 1
figure 1

The flowchart of the participant selection process in the current study. Sample a for trajectory analysis; sample b for prospective (survival) analysis

The second sample (sample B) was used in prospective (survival) analyses to link multi-trajectory groups of lipid indices (defined using sample A) with incident CVD, heart failure, and all-cause mortality later in life. This sample (sample B) included participants in sample A who attended the fourth examination cycle and did not have prevalent CVD, renal failure, and missing data on covariates on examination cycle 4 (Fig. 1, sample B). These participants were followed from examination cycle four until December 2014.

Assessment of the lipid indices and covariates

In both studies, demographic characteristics and CVD risk factors were measured using standardized protocols and similar standard and validated methods at each examination cycle. In the current analysis, we used the repeated blood lipid measurements, including LDL-C, HDL-C, and triglycerides (TG).

At each examination cycle, blood samples were collected after a 12-h fast using a standardized venipuncture procedure. EDTA plasma samples were aliquoted on ice and stored at − 70 °C until analysis. Total cholesterol and triglycerides were measured using standard enzymatic processes (Roche Diagnostics). After precipitation of non-HDL-cholesterol with magnesium/dextran, the HDL-C level was measured using the cholesterol oxidase method (Roche Diagnostics). The concentration of LDL-C was calculated from the concentrations of total cholesterol, HDL-C, and triglyceride values < 400 mg/dL by the Friedewald formula. Arterial blood pressure was measured three times on the right arm and in a sitting position after a 5-min rest; the average value of the second and third measurements was used in the analysis. Information on smoking status, alcohol consumption, medications for lowering lipid levels, blood pressure, and diabetes mellitus was obtained from standard questionnaires and bringing in drugs by participants in each examination cycle.

Outcomes

In this study, the primary outcomes of interest for the prospective analysis were incident cardiovascular disease (All-CVD defined as MI, Resuscitated Cardiac Arrest, CHD Death, Stroke, or Stroke Death), heart failure, and all-cause mortality. For these analyses, participants were followed from examination cycle 4 (1996–98 in the ARIC and 2005–2007 in the MESA study) through 2014. In both studies, events were ascertained and adjudicated using each cohort’s specific protocol. All events were ascertained by following each participant at intervals of 9–12 months using telephone calls. Then, a review committee through medical records and death certificates in both studies adjudicated them for end-point classification and assignment of incidence dates.

Statistical analysis

Descriptive statistics, as mean (SD) for continuous variables or as frequencies (%) for categorical variables, were reported for the trajectory sample (sample A) using characteristics measured at the first examination cycle and for the prospective sample (sample B) using data at the fourth examination cycle.

We used the group-based multi-trajectory modeling (GBMTM) with age as the time scale to explore the jointly longitudinal changes of lipid indices. We implemented this technique using Proc Traj [18] in Stata software version 14 (STATA Corp., TX, US). Briefly, this model is a new application of group-based trajectory modeling (GBTM). GBMTM is a semiparametric mixture model, which allows the joint modeling of the trajectories of multiple outcomes. This model identifies latent clusters of individuals who follow similar patterns through multiple outcomes using a maximum likelihood method [18,19,20]. In the GBTM models, each individual is assumed to belong to only one group, where each group has a distinct trajectory. We applied a censored normal model [21] to identify distinct trajectories of lipid concentrations.

This study jointly estimated developmental trajectories of HDL-C, LDL-C, and TG, given that all these lipids completely determine total cholesterol. Values of TG were log-transformed because of its skewed distribution. Varied GBMT models were run before selecting the best model regarding the number of groups and trajectory shapes (e.g., constant, linear, quadratic, cubic) [20, 22]. First, to identify the optimal number of distinct groups to describe heterogeneity in the longitudinal development of lipid indices, various models using 3–9 distinct groups with fixed slope variance within groups were fitted. Then, quadratic slopes were added to the model allowing for curved developmental patterns. The improvement in model fit gained by adding additional groups or shape parameters was assessed based on the Bayesian information criteria (BICs) [22]. When comparing two models with different groups or trajectory shapes, the Bayes factor was also estimated by exp(BIC1−BIC2), where BIC1 and BIC2 represent the BIC values for models 1 and 2, respectively, to assess significant change in BIC value. A tenfold difference in the Bayes factor is considered a significant difference [22]. A model with the least BICs and sufficient sample size in each multi-trajectory group (> 5% of the sample) was chosen as the best model. Finally, to ensure that our chosen model fits the data well, we assessed four model's fit diagnostic criteria as suggested by Nagin [21]: (1) an average posterior probability of assignment for each group j (AvePPj) equal to 0.7 or greater for all groups that are considered as good discrimination in classifying individuals into distinctive groups; (2) the odds of correct classification (OCCj) equal to 5 or higher for all groups; (3) reasonable similarity between the proportion of a sample assigned to a specific group and the group probabilities estimated from the model; and (4) narrow CIs of the estimated proportion.

Since the interaction term between sex and multi-trajectory groups was not statistically significant, and the overall pattern of trajectories was similar in men and women, we used a GBMTM including all participants without sex stratification. Using this algorithm, we identified seven distinct multi-trajectories for HDL-C, LDL-C, and TG. Demographic and health characteristics of the final lipid components multi-trajectory groupings were compared using Pearson’s chi-squared tests (categorical variables) or analysis of variance (ANOVA) (continuous variables).

Finally, we assessed the associations of multi-trajectory group membership (modeled as a categorical independent variable) and incidence of CVD, heart failure, and all-cause mortality on follow-up after exam four. To do so, we conducted separate Cox proportional hazards regression models for each outcome, adjusting for age, sex, race, lipid-lowering medications, smoking status, alcohol consumption, diabetes, hypertension, and BMI measured at exam four.

Finally, as a sensitivity analysis, we repeated survival analyses adjusting for the covariates mentioned above and total caloric intake and physical activity at baseline since we have data on both covariates just at the baseline exam. Collinearity between lipid indices and lifestyle factors, including physical activity, total caloric intake, smoking, and alcohol consumption, was tested using Pearson’s correlation and a variance inflation factor (VIF), which revealed no significant collinearity (VIF < 2 for all variables, Additional file 1: Table S1).

Results

Sample characteristics

The baseline characteristics of the participants are shown in Table 1. A total of 14,373 participants aged 45–84 were included in the trajectory sample to identify the developmental multi-trajectory patterns of LDL-C, HDL-C, and TG (up to 4 visits); 6534 (45.5%) were men, the mean (SD) age at baseline was 56.3 (8.1) years, and 9593 (67.7%) were white. The association between lipid indices’ multi-trajectory groups and incident CVD and all-cause mortality were examined among 12,129 participants; 5379 (44.4%) were men, the mean (SD) age at examination cycle 4 was 63.3 (7.3) years, and 8140 (67.1%) were white.

Table 1 Baseline characteristics of the participants

Characterization of multi-trajectories of lipid indices

Using the procedure and criteria mentioned above, we chose the 7-group multi-trajectory model from all investigated models. Table 2 indicates the estimates of the diagnostic criteria for judging the adequacy of the final model. Average posterior probabilities were high for all seven groups (range, 0.85–0.95) and the odds of correct classification were all well above 5. In all seven multi-trajectory groups, the average posterior probability (AvePP) was greater than 0.85, far greater than the recommended value of 0.7, indicating that the model assigned individuals to different multi-trajectory groups with little ambiguity. Further, the value for the OCC was greater than 20 for all seven groups, which is also greater than the recommendation of 5 as a general guideline for GBTM [22].

Table 2 Diagnostic criteria for judging the adequacy of the final model

Figure 2 shows the plot of the multi-trajectory groups of lipid indices and expected group percentages for each of the groups. The majority of participants (approximately 80%) exhibited decreasing LDL-C but rising TG levels and relatively stable HDL-C levels. Groups one to four and seven have worse values than the rest for most considered lipids, with groups two and three being uniformly worse and the other three groups (one, four, and seven) each having one or two for which values are high. Indeed, group two has worse but decreasing concentrations in LDL-C and borderline high but increasing in TG, representing 18.7% of the sample; group three has worse but increasing concentrations in TG and the lowest trajectory for HDL-C than the rest and represented 16.8% of the sample.

Fig. 2
figure 2

Multi-trajectory groups of LDL-C, HDL-C, and TG using group-based multi-trajectory modeling. Dots show group-specific mean observed levels while solid lines represent fitted trajectories. Lipids were modeled as a function of age, with lipid-lowering medication usage included as a time-varying covariate

Group six has the highest trajectory for HDL-C concentrations but optimal in LDL and TG, representing 6.0% of the sample. Group five has optimal values for all considered lipids and represented 13.1% of the sample. Descriptions of each of the multi-trajectory groups are presented in Table 3.

Table 3 Description of the multi-trajectory groups of lipid indices

For descriptive purposes, specific characteristics of each of the seven groups at baseline and follow-up (exam 4) are given in Table 4. There were significant differences (all p < 0.05) between all participants' demographic and clinical characteristics in the multi-trajectory groups. Briefly, More than 60 percent of individuals in groups four, five, six, and seven were women, whereas most of the individuals in groups one and three were men. The proportion of current drinking was higher in group 6 that has higher levels of HDL-C. In contrast, the proportions of lipid-lowering medications, hypertension, diabetes mellitus, and BMI were lower in this group than the rest. As expected, individuals in group 5 had the lowest values of TC and LDL-C. Furthermore, the prevalence of family history of CVD at baseline and the incidence rate of CVD at the end of follow-up were lower in this group than the rest.

Table 4 Baseline characteristics of participants by multi-trajectory groups of LDL-C, HDL-C, and TG at baseline and exam 4 in the trajectory sample (n = 14,373)

Multi-trajectory groups of lipid indices and incident cardiovascular disease, heart failure, and all-cause mortality

During the follow-up after examination cycle 4 (median = 13.5 years), there were 1133 incident CVD events, 1075 heart failure, and 2315 deaths. Table 4 presents hazard ratios (HRs) and 95% confidence intervals (CIs) of multi-trajectory groups on incident CVD, heart failure, and all-cause mortality. As expected, groups with worse trajectories for the considered lipid indices had a higher risk of developing CVD (Table 5). The proportion of incident events was higher among individuals whose LDL-C levels are the highest and decreasing, TG levels are borderline high and increasing, and HDL-C levels are low (Fig. 2, group 2). Compared with the individuals with an optimal and stable level of the considered lipid indices (Fig. 2, group 5), those in the other multi-trajectory groups were at a statistically significant increased risk of incident CVD after adjusting for sex, age, race, lipid-lowering medication use, diabetes mellitus, hypertension, body mass index, smoking status, and alcohol consumption (Table 5). However, this risk was not similar across multi-trajectory groups. Individuals with the highest but decreasing LDL-C and borderline high and rising TG levels over time (Fig. 2, group 2) were at the highest risk than those in other groups with a 2.22-fold risk of CVD. Also, those with the highest and increased triglyceride levels over time, over optimal and decreasing LDL-C levels, and the lowest HDL-C profile (Fig. 2, group 3) had a nearly 1.84 times CVD risk. Furthermore, individuals with the highest HDL-C level (Fig. 1, group 6) had a 1.45-fold risk of CVD compared to those with the optimal levels of lipid indices (Fig. 2, group 5). Furthermore, only the highest HDL-C level (Fig. 2, group 6) was significantly associated with incident heart failure with a 1.5 times increased risk (95% CI 1.07–2.06).

Table 5 Associations between multi-trajectory groups of lipid indices and incident CVD, heart failure, and all-cause mortality

We did not observe a statistically significant difference in the risk of death between individuals in the referent group and those in other multi-trajectory groups (Table 5). Nevertheless, the multi-trajectory group with the highest HDL-C level increased the risk of death by 20 percent with a marginal significance.

Results from sensitivity analyses additionally adjusted for total caloric intake and physical activity at the baseline exam produced similar findings as our primary analyses, except for multi-trajectory groups 2 and 3 that demonstrated to be significantly associated with incident heart failure with 1.31 and 1.42 times increased risk (95% CI 1.03–1.68 and 1.09–1.84), respectively.

Discussion

This large, pooled cohort study investigated heterogeneity in lipid profile multi-trajectories among individuals aged 45 to 84 years. Using group-based multi-trajectory analyses of the longitudinal data across the four examination cycles, we identified seven distinct multi-trajectory groups of LDL-C, HDL-C, and TG concentrations. Approximately 80% of participants exhibited decreasing LDL-C levels but rising TG levels and relatively stable levels of HDL-C.

We revealed that compared with the individuals presenting with optimal and stable levels of these lipids across the life course, those in the other multi-trajectory groups were at an increased risk of incident CVD after adjusting for confounders. However, individuals in different groups showed a different risk of CVD. We also observed that individuals presenting with very high HDL-C levels throughout the adult life course were at a significantly increased risk of incident CVD and heart failure, and a marginally significant risk of death than those with optimal lipid profile, after adjusting for confounders. Besides, In contrast with very high levels of HDL-C, the risk of death and incidence of heart failure did not differ across other multi-trajectories compared to the optimal lipid indices levels.

Although numerous studies have been suggested the association of lipids with the risk of CVD, the trajectories we identified in this study provide new insights for the common progression of lipid indices that could be expected to be observed during the age of 45 to 84 years in relatively healthy adults. We used a person-centered, multi-trajectory approach that modeled the common progressions of LDL-C, HDL-C, and TG levels. This model is considerably different because for identifying and monitoring the various lipid indices progressions herein, the usage of lipid-lowering medications simultaneously incorporated and accounted for correlation of lipids within the same participant and over time. Such an approach defines a trajectory group in terms of trajectories for multiple indicators, not just one indicator. In so doing, the model efficiently represents the interrelationship of numerous clinically relevant indicators. Such approaches visually represent distinct groups of individuals who display unique lipid indices patterns over time [20].

Previous studies using trajectory modeling only investigated the development of each lipid component separately [3, 10,11,12, 23]. To date, only one previous study by Dayimu et al. has considered the trajectory of these lipids jointly. It has shown three distinct trajectory classes (U-shape class, progressing, and inverse U-shape) in a Chinese population aged 20 to 60 years [9]. Our findings widen previous evidence in this field by showing that lipids can be jointly categorized into seven different multi-trajectory groups over the age of 45–84.

Interestingly, we found a decreasing trend in LDL-C over time, although only 40 percent of participants were on lipid-lowering agents. Moreover, when we excluded those, who had been receiving lipid-lowering medications from our analysis, the trend remained relatively unchanged (Additional file 2: Fig. S1). Thus, this favorite trend might also be attributable to the FDA’s regulations and national guidelines regarding Americans’ dietary fat and cholesterol intake [24,25,26]. Yet, even after falls in LDL-C, a considerable amount of CVD risk remains. So, these findings suggest the impact of the increasing trend of TG on residual cardiovascular risk [27, 28].

Although the exact role of TG in inducing ASCVD has been controversial, the evidence supporting the association of elevated concentrations of triglyceride-rich lipoproteins or remnant cholesterol, reliably marked by raised triglycerides, and cardiovascular disease and all-cause mortality, is increasing. Taken together, numerous observational and genetic studies strongly support the association of mild-to-moderately raised triglyceride concentrations and CVD risk [29,30,31,32,33]. Previous studies have confirmed that elevated TG levels even below 150 mg/dL, previously considered “optimal” levels, were associated with increased CVD risk [34,35,36]. This evidence suggesting that a biologically “optimal” level may be even lower for TG as American Heart Association also indicated that an “optimal” fasting TG level is less than 100 mg/dL [37].

We also revealed that individuals with a higher TG level had a lower HDL-C, depicted in the trajectory plot (Fig. 2, groups 2 and 3), and were at a greater risk of incident CVD than other individuals. Our finding is consistent with evidence from some genetic studies and randomized trials suggesting that low HDL cholesterol might merely be a long-term marker of raised triglycerides and remnant cholesterol, not a cause of CVD [38,39,40,41,42]. A Mendelian randomization study using genetic variants that affect the remnant cholesterol concentrations, HDL cholesterol, or both demonstrated that a 1 mmol/L (39 mg/dL) increase in remnant cholesterol is associated with a 2.8-fold risk for ischemic heart disease, independent of reduced HDL cholesterol [31]. Another study also showed that increasing TGs were associated with more significant increases in CVD risk among individuals with higher HDL-C levels [34].

Our findings support the associations of elevated TG with the risk of future CVD, and the use of trajectory of TGs to identify high-risk individuals for CVD events that are in accordance with some studies showing that using several TG measurements improved prediction of CVD risk more than a single TG measurement [34]. However, whether TGs are an independent cause for incident CVD or serve as a marker for other risk factors remains unclear. Evidence from some genetic studies suggested that the risk associated with elevated TGs might be because of their association with elevations in non–HDL-C or apolipoprotein B (apo B)-containing lipoproteins [43].

Finally, some randomized trials also showed an increased risk for CVD in statin-treated patients with elevated TGs [44,45,46,47]. A randomized trial known as Reduction of Cardiovascular Events with EPA-Intervention Trial (REDUCE-IT) assessed the impact of TG-lowering using Icosapent ethyl among statin-treated patients with TGs ≥ 135 and < 500 mg/dL, and with a history of CVD, diabetes, or other risk factors. Its finding showed that lowering TG was significantly associated with a lower risk of CVD [48]. Recently, a systematic review and meta-regression analysis of randomized controlled trials also showed a significant association between TG-lowering and CV risk reduction, even after adjusting for LDL-C lowering, although the effect attenuated when REDUCE-IT was excluded from the analysis [49].

Another contribution of this study is that we found a group of individuals (approximately 6%) with very high and stable levels of HDL-C but optimal LDL-C and TG levels that were at risk of CVD, heart failure, and even a marginally significant risk of death. Our finding is consistent with evidence from clinical studies showing that both serum HDL-C concentration/quantity and its qualities/properties can play a critical role in determining its overall effects and, hence, its association with clinical outcomes [50,51,52]. So, not only HDL serum cholesterol concentrations but also a range of other properties, including its particle size and composition, its Apo lipoprotein content, its enrichment with proinflammatory properties, and its functional capacity, can play a critical role in determining its overall effects and hence its association with clinical outcomes [50, 53].

Finally, compared to previous studies that reported associations between lipid profile trajectories and all-cause mortality [9, 10], we did not observe a statistically significant association of multi-trajectory groups with death. However, the highest HDL-C but optimal LDL-C and TG profile had a marginally significant risk of death. Several large-scale prospective cohort studies have recently revealed a U-shaped association between HDL-C levels and all-cause mortality, indicating that both very high and low levels are associated with an increased risk of death [54, 55]. Some other studies also reported that high LDL-C and HDL-C levels are inversely associated with mortality, especially in older people [56, 57].

This study has several strengths, including a pooled data set from two large representative cohort studies with multiple examinations across adulthood. Besides, continuous follow-up for collecting data on various variables and CVD and death events at ARIC and MESA studies allowed us for inclusive adjustment for risk factors. A key strength of this study is the use of an innovative multi-trajectory modeling technique to identify subgroups of longitudinal lipid profile trajectories based on multiple lipids. The multi-trajectory analysis used in this study incorporates the inter-correlations among the multiple lipids to improve the accuracy of individual-specific probabilities of group membership, while the conventional group-based trajectory analysis clusters longitudinal trajectories based on one outcome. Despite these strengths, this study has limitations. First, we used repeated measurements of lipids just over the four examination cycles; we could not use the examination cycle 5 to ensure sufficient follow-up time for the prospective analyses, so longer variation in the lipid pattern might be missed. However, some evidence clarifies that at least three unique time points are required for the trajectory approach [58, 59]. Second, despite the large sample size included in this study, the number of individuals with extremely high HDL-C concentrations was relatively small (~ 6% of the total sample), especially in the stratified analyses, limiting statistical power. Third, in common with all modeling approaches, there are limitations attendant to trajectory models. Although several model diagnostic criteria were proposed as guides, the problem of specifying the correct number of groups has not been entirely resolved. Fourth, it should also be noted that GBTM attempts to provide potentially clinically meaningful trajectory groups based on the available data, and they should not be taken to exist literally.

In conclusion, the present observations provide a comprehensive depiction of the joint progression of lipid indices over time. The trajectories and risk of CVD identified in this study demonstrated that despite a decline in LDL-C over time, a significant amount of residual risk for CVD remains. These findings suggest the impact of the increasing trend of TG on CVD risk and emphasize the importance of assessing the lipid levels at each visit and undertaking potential interventions that lower triglyceride concentrations to reduce the residual risk of CVD. Since our findings are exploratory and do not address treatment questions, future research should focus more on the effects of TG-lowering strategies in reducing residual cardiovascular risk, even among those with the optimal LDL-C level.

Availability of data and materials

The datasets used during the current study are available from the corresponding author on reasonable request.

Abbreviations

CVD:

Cardiovascular disease

LDL-C:

Low-density lipoprotein cholesterol

HDL-C:

High-density lipoprotein cholesterol

TG:

Triglyceride

BMI:

Body mass index

MI:

Myocardial infraction

CHD:

Chronic heart disease

ARIC:

Atherosclerosis risk in communities

MESA:

Multi-ethnic study of atherosclerosis

GBTM:

Group-based trajectory model

AvePPj:

Average posterior probability of assignment for each group j

OCCj:

Odds of correct classification

FDA:

Food and Drug Administration

References

  1. World Health Organization (WHO). Cardiovascular diseases. https://www.who.int/health-topics/cardiovascular-diseases/#tab=tab_1. Accessed 18 July 2020.

  2. Ryan A, Heath S, Cook P. Managing dyslipidaemia for the primary prevention of cardiovascular disease. BMJ. 2018. https://doi.org/10.1136/bmj.k946.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Park YM, et al. The effect of cardiorespiratory fitness on age-related lipids and lipoproteins. J Am Coll Cardiol. 2015;65:2091–100. https://doi.org/10.1016/j.jacc.2015.03.517.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Upmeier E, et al. Longitudinal changes in serum lipids in older people the Turku elderly study 1991–2006. Age Ageing. 2011;40:280–3. https://doi.org/10.1093/ageing/afq180.

    Article  PubMed  Google Scholar 

  5. Schubert C, et al. Lipids, lipoproteins, lifestyle, adiposity and fat-free mass during middle age: the Fels longitudinal study. Int J Obes. 2006;30:251–60. https://doi.org/10.1038/sj.ijo.0803129.

    Article  CAS  Google Scholar 

  6. Wilson PW, Anderson KM, Harri T, Kannel WB, Castelli WP. Determinants of change in total cholesterol and HDL-C with age: the Framingham study. J Gerontol. 1994;49:M252–7. https://doi.org/10.1093/geronj/49.6.M252.

    Article  CAS  PubMed  Google Scholar 

  7. Ueda P, Gulayin P, Danaei G. Long-term moderately elevated LDL-cholesterol and blood pressure and risk of coronary heart disease. PLoS ONE. 2018;13: e0200017. https://doi.org/10.1371/journal.pone.0200017.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Navar-Boggan AM, et al. Hyperlipidemia in early adulthood increases long-term risk of coronary heart disease. Circulation. 2015;131:451–8. https://doi.org/10.1161/CIRCULATIONAHA.114.012477.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Dayimu A, et al. Trajectories of lipids profile and incident cardiovascular disease risk: a longitudinal cohort study. J Am Heart Assoc. 2019;8:e013479. https://doi.org/10.1161/JAHA.119.013479.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Duncan MS, Vasan RS, Xanthakis V. Trajectories of blood lipid concentrations over the adult life course and risk of cardiovascular disease and all-cause mortality: observations from the Framingham study over 35 years. J Am Heart Assoc. 2019;8:e011433. https://doi.org/10.1161/JAHA.118.011433.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Elovainio M, et al. Lipid trajectories as predictors of depressive symptoms: the Young Finns study. Health Psychol. 2010;29:237. https://doi.org/10.1037/a0018875.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Lee Y-C, et al. High-density lipoprotein-cholesterol trajectory pattern, associated lifestyle and biochemical factors among Taiwanese. Circ J. 2009. https://doi.org/10.1253/circj.cj-09-0209.

    Article  PubMed  Google Scholar 

  13. Hill C, et al. The atherosclerosis risk in communities (ARIC) study: design and objectives. The ARIC investigators. Am J Epidemiol. 1989;129:687–702. https://doi.org/10.1093/oxfordjournals.aje.a115184.

    Article  Google Scholar 

  14. Bild DE, et al. Multi-ethnic study of atherosclerosis: objectives and design. Am J Epidemiol. 2002;156:871–81. https://doi.org/10.1093/aje/kwf113.

    Article  PubMed  Google Scholar 

  15. Atherosclerosis risk in communities study description. https://sites.cscc.unc.edu/aric/desc_pub.

  16. MESA website. https://www.mesa-nhlbi.org/.

  17. Fortier I, et al. Maelstrom research guidelines for rigorous retrospective data harmonization. Int J Epidemiol. 2017;46:103–5. https://doi.org/10.1093/ije/dyw075.

    Article  PubMed  Google Scholar 

  18. Jones BL, Nagin D. A Stata plugin for estimating group-based trajectory models. Sociol Methods Res. 2018. https://doi.org/10.1184/R1/6470963.v1.

    Article  Google Scholar 

  19. Jones BL, Nagin DS, Roeder K. A SAS procedure based on mixture models for estimating developmental trajectories. Sociol Methods Res. 2001;29:374–93. https://doi.org/10.1177/0049124101029003005.

    Article  Google Scholar 

  20. Nagin DS, Jones BL, Passos VL, Tremblay RE. Group-based multi-trajectory modeling. Stat Methods Med Res. 2018;27:2015–23. https://doi.org/10.1177/0962280216673085.

    Article  PubMed  Google Scholar 

  21. Jones BL, Nagin DS. A note on a Stata plugin for estimating group-based trajectory models. Sociol Methods Res. 2013;42:608–13. https://doi.org/10.1177/0049124113503141.

    Article  Google Scholar 

  22. Nagin DS. Group-based modeling of development. Cambridge: Harvard University Press; 2005. https://doi.org/10.4159/9780674041318.

    Book  Google Scholar 

  23. Tsai C-W, et al. Longitudinal lipid trends and adverse outcomes in patients with CKD: a 13-year observational cohort study. J Lipid Res. 2019;60:648–60. https://doi.org/10.1194/jlr.p084590.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Ford ES, Capewell S. Trends in total and low-density lipoprotein cholesterol among U.S. adults: contributions of changes in dietary fat intake and use of cholesterol-lowering medications. PLoS ONE. 2013;8:e65228. https://doi.org/10.1371/journal.pone.0065228.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Ernst ND, Sempos CT, Briefel RR, Clark MB. Consistency between US dietary fat intake and serum total cholesterol concentrations: the national health and nutrition examination surveys. Am J Clin Nutr. 1997;66:965S-S972. https://doi.org/10.1093/ajcn/66.4.965S.

    Article  CAS  PubMed  Google Scholar 

  26. Johnson CL, et al. Declining serum total cholesterol levels among US adults. The national health and nutrition examination surveys. JAMA. 1993;269:3002–8.

    Article  CAS  Google Scholar 

  27. Ganda OP, Bhatt DL, Mason RP, Miller M, Boden WE. Unmet need for adjunctive dyslipidemia therapy in hypertriglyceridemia management. J Am Coll Cardiol. 2018;72:330–43.

    Article  Google Scholar 

  28. Nejat A, et al. Changes in lipid measures and incident coronary heart disease: Tehran lipid & glucose study. Clin Biochem. 2014;47:1239–44.

    Article  CAS  Google Scholar 

  29. Nordestgaard BG, Benn M, Schnohr P, Tybjærg-Hansen A. Nonfasting triglycerides and risk of myocardial infarction, ischemic heart disease, and death in men and women. JAMA. 2007;298:299–308.

    Article  CAS  Google Scholar 

  30. Sarvar R, Danesh J, Eriksdottir J. Triglycerides and the risk of heart disease. 10 158 incident cases among 262 525 participants in 29 western population studies. Circulation. 2007;115:450–8.

    Article  Google Scholar 

  31. Varbo A, et al. Remnant cholesterol as a causal risk factor for ischemic heart disease. J Am Coll Cardiol. 2013;61:427–36.

    Article  CAS  Google Scholar 

  32. Consortium TCDG, Collaboration ERF. Triglyceride-mediated pathways and coronary disease: collaborative analysis of 101 studies. Lancet. 2010;375:1634–9.

    Article  Google Scholar 

  33. Jørgensen AB, et al. Genetically elevated non-fasting triglycerides and calculated remnant cholesterol as causal risk factors for myocardial infarction. Eur Heart J. 2013;34:1826–33. https://doi.org/10.1093/eurheartj/ehs431.

    Article  CAS  PubMed  Google Scholar 

  34. Aberra T, et al. The association between triglycerides and incident cardiovascular disease: what is “optimal”? J Clin Lipidol. 2020;14:438–47.

    Article  Google Scholar 

  35. Lawler PR, et al. Real-world risk of cardiovascular outcomes associated with hypertriglyceridaemia among individuals with atherosclerotic cardiovascular disease and potential eligibility for emerging therapies. Eur Heart J. 2020;41:86–94.

    Article  Google Scholar 

  36. Kajikawa M, et al. Target of triglycerides as residual risk for cardiovascular events in patients with coronary artery disease—post hoc analysis of the FMD-J study A. Circ J. 2019;83:1064–71.

    Article  CAS  Google Scholar 

  37. Miller M, et al. Triglycerides and cardiovascular disease: a scientific statement from the American Heart Association. Circulation. 2011;123:2292–333.

    Article  Google Scholar 

  38. Nordestgaard B, Langsted A, Freiberg J. Nonfasting hyperlipidemia and cardiovascular disease. Curr Drug Targets. 2009;10:328–35.

    Article  CAS  Google Scholar 

  39. Frikke-Schmidt R, et al. Association of loss-of-function mutations in the ABCA1 gene with high-density lipoprotein cholesterol levels and risk of ischemic heart disease. JAMA. 2008;299:2524–32. https://doi.org/10.1001/jama.299.21.2524.

    Article  CAS  PubMed  Google Scholar 

  40. Haase CL, et al. LCAT, HDL cholesterol and ischemic cardiovascular disease: a Mendelian randomization study of HDL cholesterol in 54,500 individuals. J Clin Endocrinol Metab. 2012;97:E248–56. https://doi.org/10.1210/jc.2011-1846.

    Article  CAS  PubMed  Google Scholar 

  41. Boden WE, et al. Niacin in patients with low HDL cholesterol levels receiving intensive statin therapy. N Engl J Med. 2011;365:2255–67. https://doi.org/10.1056/NEJMoa1107579.

    Article  CAS  PubMed  Google Scholar 

  42. Landray MJ, et al. Effects of extended-release niacin with laropiprant in high-risk patients. N Engl J Med. 2014;371:203–12. https://doi.org/10.1056/NEJMoa1300955.

    Article  CAS  PubMed  Google Scholar 

  43. Ference BA, et al. Association of triglyceride-lowering LPL variants and LDL-C-lowering LDLR variants with risk of coronary heart disease. JAMA. 2019;321:364–73.

    Article  CAS  Google Scholar 

  44. Chapman MJ, et al. Triglyceride-rich lipoproteins and high-density lipoprotein cholesterol in patients at high risk of cardiovascular disease: evidence and guidance for management. Eur Heart J. 2011;32:1345–61.

    Article  CAS  Google Scholar 

  45. Faergeman O, et al. Plasma triglycerides and cardiovascular events in the treating to new targets and incremental decrease in end-points through aggressive lipid lowering trials of statins in patients with coronary artery disease. Am J Cardiol. 2009;104:459–63. https://doi.org/10.1016/j.amjcard.2009.04.008.

    Article  CAS  PubMed  Google Scholar 

  46. Simes RJ, et al. Relationship between lipid levels and clinical outcomes in the long-term intervention with pravastatin in ischemic disease (LIPID) trial: to what extent is the reduction in coronary events with pravastatin explained by on-study lipid levels? Circulation. 2002;105:1162–9. https://doi.org/10.1161/hc1002.105136.

    Article  CAS  PubMed  Google Scholar 

  47. Miller M, et al. Impact of triglyceride levels beyond low-density lipoprotein cholesterol after acute coronary syndrome in the PROVE IT-TIMI 22 trial. J Am Coll Cardiol. 2008;51:724–30.

    Article  CAS  Google Scholar 

  48. Bhatt DL, et al. Cardiovascular risk reduction with icosapent ethyl for hypertriglyceridemia. N Engl J Med. 2019;380:11–22.

    Article  CAS  Google Scholar 

  49. Marston NA, et al. Association between triglyceride lowering and reduction of cardiovascular risk across multiple lipid-lowering therapeutic classes: a systematic review and meta-regression analysis of randomized controlled trials. Circulation. 2019;140:1308–17.

    Article  CAS  Google Scholar 

  50. Chang TI, Streja E, Moradi H. Could high-density lipoprotein cholesterol predict increased cardiovascular risk? Curr Opin Endocrinol Diabetes Obes. 2017;24:140–7. https://doi.org/10.1097/med.0000000000000318.

    Article  CAS  PubMed  Google Scholar 

  51. Kosmas CE, Silverio D, Sourlas A, Montan PD, Guzman E. Dysfunctional high-density lipoprotein and atherogenesis. Vessel Plus. 2019;3:2. https://doi.org/10.20517/2574-1209.2018.79.

    Article  CAS  Google Scholar 

  52. Xiang AS, Kingwell BA. Rethinking good cholesterol: a clinicians’ guide to understanding HDL. Lancet Diabetes Endocrinol. 2019;7:575–82. https://doi.org/10.1016/S2213-8587(19)30003-8.

    Article  PubMed  Google Scholar 

  53. Moradi H, Vaziri ND, Kashyap ML, Said HM, Kalantar-Zadeh K. Role of HDL dysfunction in end-stage renal disease: a double-edged sword. J Ren Nutr. 2013;23:203–6. https://doi.org/10.1053/j.jrn.2013.01.022.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Madsen CM, Varbo A, Nordestgaard BG. Extreme high high-density lipoprotein cholesterol is paradoxically associated with high mortality in men and women: two prospective cohort studies. Eur Heart J. 2017;38:2478–86. https://doi.org/10.1093/eurheartj/ehx163.

    Article  CAS  PubMed  Google Scholar 

  55. Yi S-W, Park S-J, Yi J-J, Ohrr H, Kim H. High-density lipoprotein cholesterol and all-cause mortality by sex and age: a prospective cohort study among 15.8 million adults. Int J Epidemiol. 2020. https://doi.org/10.1093/ije/dyaa243.

    Article  PubMed  Google Scholar 

  56. Ravnskov U, et al. Lack of an association or an inverse association between low-density-lipoprotein cholesterol and mortality in the elderly: a systematic review. BMJ Open. 2016;6: e010401. https://doi.org/10.1136/bmjopen-2015-010401.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Li X, et al. Association between high-density lipoprotein cholesterol and all-cause mortality in the general population of northern China. Sci Rep. 2019;9:1–8. https://doi.org/10.1038/s41598-019-50924-4.

    Article  CAS  Google Scholar 

  58. Bollen KA, Curran PJ. Latent curve models: a structural equation perspective. Hoboken: Wiley; 2006.

    Google Scholar 

  59. Pollock BD, et al. Life course trajectories of cardiovascular risk: impact on atherosclerotic and metabolic indicators. Atherosclerosis. 2019;280:21–7. https://doi.org/10.1016/j.atherosclerosis.2018.11.008.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This article is part of the Ph.D. thesis of the first author supported by the School of Public Health and Safety, Shahid Beheshti University of Medical Sciences, Tehran, Iran. This study was carried out using the ARIC and MESA Research Materials obtained from the NHLBI; BioLINCC Research Materials Distribution Agreement (RMDA) V02 1d20120806. We would like to extend our thanks to the staff and coordinators of BioLINC and the other investigators, the staff, and the participants of the ARIC study and MESA for their valuable contributions.

Funding

No funds, grants, or other support was received.

Author information

Authors and Affiliations

Authors

Contributions

All authors contributed to the study's conception and design. FK, DK, MM, and HS raised the presented idea and designed the study. FK performed statistical analyses under the supervision of senior authors DK, MM, and HS. FK, DK, MM, FH, and HS interpreted the results. FK developed the first draft of the manuscript. DK, MM, FH, and HS critically reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hamid Soori.

Ethics declarations

Ethics approval and consent to participate

The current study was performed following the ethical standards in the Declaration of Helsinki. The institutional review board approved the ARIC and MESA study protocols of each participating center. All methods were carried out following relevant guidelines and regulations. Approval for undertaking the current project was also obtained from the School of Public Health & Neuroscience Research center, Shahid Beheshti University of Medical Sciences, Tehran, Iran (IR.SBMU.PHNS.REC.1398.090).

Consent to participate

All ARIC and MESA participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors have no relevant financial or non-financial interests to disclose.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Table S1.

Correlation coefficients among lipid profile and lifestyle factors.

Additional file 2: Figure S1.

Multi-trajectory groups of LDL-C, HDL-C, and TG among participants not on lipid-lowering treatment. Dots show group-specific mean observed levels while solid lines represent fitted trajectories. Lipids were modeled as a function of age.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Koohi, F., Khalili, D., Mansournia, M.A. et al. Multi-trajectories of lipid indices with incident cardiovascular disease, heart failure, and all-cause mortality: 23 years follow-up of two US cohort studies. J Transl Med 19, 286 (2021). https://doi.org/10.1186/s12967-021-02966-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-021-02966-4

Keywords