Dyslipidaemia is associated with elevated cardiovascular risks, with the INTERHEART study observing a tripling of myocardial infarction (MI) risk in patients with dyslipidaemia.1 Most studies focused on mean levels or point estimates of low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C), despite well-known visit-to-visit variability.2 Visit-to-visit cholesterol variability, reflecting fluctuations in cholesterol levels between visits, is prognostic for some adverse cardiovascular outcomes such as cardiac arrhythmias and mortality.3,4 Nonetheless, associations between cholesterol variability and heart failure (HF) remain unclear. This study therefore investigated the associations between LDL-C and HDL-C variabilities and the risk of new-onset HF and major adverse cardiovascular outcomes.

This retrospective cohort study was performed according to the Declaration of Helsinki and approved by the local institutional review board. Patient consent was not required since deidentified data were used. Underlying data are available upon reasonable request to the corresponding authors. Data were retrieved from Clinical Data Analysis and Reporting System, a city-wide medical database in Hong Kong. Diagnoses are coded by International Classification of Diseases Ninth revision codes, which were used for identifying outcomes and comorbidities (see Supplementary material online, Table S1). Mortality data were obtained from the Hong Kong Death registry, a governmental registry of all Hong Kong citizens’ death records. Both have been used in prior publications.5

Adult patients who attended any family medicine clinic in Hong Kong between 1 January 2000–31 December 2003 with more than three sets of fasting lipid profiles were included. Exclusion criteria were known HF, prior MI, using HF medications, pregnancy, known human immunodeficiency virus infection, and less than three sets of lipid profiles. Patients were followed up until 31 December 2019. Primary outcome was new-onset HF. Secondary outcomes were cardiovascular mortality, and MI.

Fasting LDL-C, HDL-C, and total cholesterol during index period were recorded, from which mean, standard deviation (SD), and coefficient of variation (CV) were calculated. As SD is only interpretable together with mean, CV (SDmean) was reported to facilitate interpretation.

Student’s t-test and Fisher’s exact test were used to compare continuous and binary variables between subgroups, respectively. Univariable Cox regression identified covariates (P < 0.010; Supplementary material online, Tables S3–5) for adjustments in multivariable Cox regression which assessed associations between cholesterol measurements and outcomes. A prespecified subgroup analysis by the use of any lipid-lowering medication(s) at baseline was performed. Marginal effects were calculated by identifying partial derivatives of cholesterol measurements in Cox models for each data unit. All p-values were two-sided, with P < 0.05 considered significant. All analyses were performed on Statistical Package for Social Sciences (v25.0), Stata (v13.0), RStudio (v1.1.456), and/or Python (v3.6).

In total, 155 065 patients were identified, of which 5662 were included [2152 (38.0%) males; mean age 63.3 ± 12.4 years; Supplementary material online, Figure S1 and Supplementary material online, Table S2]. Over a mean follow-up of 15.3 ± 4.6 years with 13 ± 9 sets of lipid profile (median of per-patient median between-test duration 0.9 years, interquartile range 0.5–1.1 years), new-onset HF, cardiovascular mortality, and MI occurred in 1196 (21.1%), 548 (9.7%), and 656 (11.6%) patients, respectively. Patients with baseline use of lipid-lowering medication(s) had higher variability of LDL-C, HDL-C, and total cholesterol.

HDL-C measures, but not LDL-C measures, were strongly associated with new-onset HF risk (Table 1). Higher HDL-C variability was associated with significantly higher new-onset HF risk [adjusted hazard ratio (aHR) for CV 48.68 (21.37, 110.90), P < 0.001, Figure 1F]. No significant association was observed between LDL-C variability and new-onset HF risk. Lower mean HDL-C was associated with significantly higher new-onset HF risk [aHR 0.60 (0.50, 0.72), P < 0.001; Figure 1D], but not mean LDL-C. Meanwhile, cholesterol variability was not significantly associated with cardiovascular mortality, whilst higher mean LDL-C [aHR 1.37 (1.20, 1.55), P < 0.001; Supplementary material online, Figure S2A] and lower mean HDL-C [aHR 0.54 (0.41, 0.71), P < 0.001; Supplementary material online, Figure S2D)] were both associated with higher cardiovascular mortality risk. Furthermore, higher variability of LDL-C [aHR for CV 3.64 (1.81, 7.33), P < 0.001; Supplementary material online, Figure S3C] and HDL-C (aHR for CV 39.87 [13.77, 115.47], P < 0.001; Supplementary material online, Figure S3F) were associated with higher MI risk. Lower mean HDL-C was associated with higher MI risk [aHR 0.46 (0.36, 0.60), P < 0.001; Supplementary material online, Figure S3D], but not mean LDL-C. Total cholesterol measures showed congruent results.

Marginal effects of (A) mean low density lipoprotein cholesterol (LDL-C) level, (B) standarddeviation (SD) of LDL-C level, (C) coefficient of variation (CV) of LDL-C level, (D) mean high-density lipoprotein cholesterol (HDL-C) level, (E) SD of HDL-C level, and (F) CV of HDL-C level on the risk of heart failure.
Figure 1

Marginal effects of (A) mean low density lipoprotein cholesterol (LDL-C) level, (B) standarddeviation (SD) of LDL-C level, (C) coefficient of variation (CV) of LDL-C level, (D) mean high-density lipoprotein cholesterol (HDL-C) level, (E) SD of HDL-C level, and (F) CV of HDL-C level on the risk of heart failure.

Table 1

Results of multivariable Cox regression analysis for the overall study cohort

Hazard ratio [95% confidence interval], P value
New-onset heart failureaCardiovascular mortalitybMyocardial infarctionc
Mean LDL-C1.00 [0.91, 1.10], P = 0.9461.37 [1.20, 1.55], P < 0.0010.94 [0.83, 1.07], P = 0.368
SD of LDL-C1.03 [0.85, 1.25], P = 0.7450.94 [0.72, 1.24], P = 0.6771.46 [1.15, 1.85], P = 0.002
CV of LDL-C0.63 [0.36, 1.11], P = 0.1090.19 [0.08, 0.42], P < 0.0013.64 [1.81, 7.33], P < 0.001
Mean HDL-C0.60 [0.50, 0.72], P < 0.0010.54 [0.41, 0.71], P < 0.0010.46 [0.36, 0.60], P < 0.001
SD of HDL-C4.40 [2.72, 7.13], P < 0.0011.05 [0.46, 2.39], P = 0.9012.82 [1.42, 5.59], P = 0.003
CV of HDL-C48.68 [21.37, 110.90], P < 0.0011.63 [0.45, 5.97], P = 0.46039.87 [13.77, 115.47], P < 0.001
Mean total cholesterol0.89 [0.82, 0.96], P = 0.0041.21 [1.08, 1.36], P = 0.0010.87 [0.78, 0.98], P = 0.017
SD of total cholesterol1.38 [1.18, 1.62], P < 0.0011.18 [0.93, 1.51], P = 0.1691.69 [1.39, 2.06], P < 0.001
CV of total cholesterol5.49 [2.29, 13.15], P < 0.0010.58 [0.16, 2.11], P = 0.40629.94 [9.83, 91.19], P < 0.001
Hazard ratio [95% confidence interval], P value
New-onset heart failureaCardiovascular mortalitybMyocardial infarctionc
Mean LDL-C1.00 [0.91, 1.10], P = 0.9461.37 [1.20, 1.55], P < 0.0010.94 [0.83, 1.07], P = 0.368
SD of LDL-C1.03 [0.85, 1.25], P = 0.7450.94 [0.72, 1.24], P = 0.6771.46 [1.15, 1.85], P = 0.002
CV of LDL-C0.63 [0.36, 1.11], P = 0.1090.19 [0.08, 0.42], P < 0.0013.64 [1.81, 7.33], P < 0.001
Mean HDL-C0.60 [0.50, 0.72], P < 0.0010.54 [0.41, 0.71], P < 0.0010.46 [0.36, 0.60], P < 0.001
SD of HDL-C4.40 [2.72, 7.13], P < 0.0011.05 [0.46, 2.39], P = 0.9012.82 [1.42, 5.59], P = 0.003
CV of HDL-C48.68 [21.37, 110.90], P < 0.0011.63 [0.45, 5.97], P = 0.46039.87 [13.77, 115.47], P < 0.001
Mean total cholesterol0.89 [0.82, 0.96], P = 0.0041.21 [1.08, 1.36], P = 0.0010.87 [0.78, 0.98], P = 0.017
SD of total cholesterol1.38 [1.18, 1.62], P < 0.0011.18 [0.93, 1.51], P = 0.1691.69 [1.39, 2.06], P < 0.001
CV of total cholesterol5.49 [2.29, 13.15], P < 0.0010.58 [0.16, 2.11], P = 0.40629.94 [9.83, 91.19], P < 0.001

CV, coefficient of variation; SD, standard deviation.

a

Adjusted for age, hypertension, diabetes mellitus, atrial fibrillation, chronic obstructive pulmonary disease, ischaemic heart disease, peripheral vascular disease, stroke or transient ischaemic attack, dementia or Alzheimer’s disease, and use of lipid-lowering medication(s).

b

Adjusted for age, sex, hypertension, diabetes mellitus, atrial fibrillation, ischaemic heart disease, peripheral vascular disease, stroke or transient ischaemic attack, dementia or Alzheimer’s disease, and use of lipid-lowering medication(s).

c

Adjusted for age, sex, hypertension, diabetes mellitus, chronic obstructive pulmonary disease, ischaemic heart disease, peripheral vascular disease, stroke or transient ischaemic attack, and use of lipid-lowering medication(s).

Table 1

Results of multivariable Cox regression analysis for the overall study cohort

Hazard ratio [95% confidence interval], P value
New-onset heart failureaCardiovascular mortalitybMyocardial infarctionc
Mean LDL-C1.00 [0.91, 1.10], P = 0.9461.37 [1.20, 1.55], P < 0.0010.94 [0.83, 1.07], P = 0.368
SD of LDL-C1.03 [0.85, 1.25], P = 0.7450.94 [0.72, 1.24], P = 0.6771.46 [1.15, 1.85], P = 0.002
CV of LDL-C0.63 [0.36, 1.11], P = 0.1090.19 [0.08, 0.42], P < 0.0013.64 [1.81, 7.33], P < 0.001
Mean HDL-C0.60 [0.50, 0.72], P < 0.0010.54 [0.41, 0.71], P < 0.0010.46 [0.36, 0.60], P < 0.001
SD of HDL-C4.40 [2.72, 7.13], P < 0.0011.05 [0.46, 2.39], P = 0.9012.82 [1.42, 5.59], P = 0.003
CV of HDL-C48.68 [21.37, 110.90], P < 0.0011.63 [0.45, 5.97], P = 0.46039.87 [13.77, 115.47], P < 0.001
Mean total cholesterol0.89 [0.82, 0.96], P = 0.0041.21 [1.08, 1.36], P = 0.0010.87 [0.78, 0.98], P = 0.017
SD of total cholesterol1.38 [1.18, 1.62], P < 0.0011.18 [0.93, 1.51], P = 0.1691.69 [1.39, 2.06], P < 0.001
CV of total cholesterol5.49 [2.29, 13.15], P < 0.0010.58 [0.16, 2.11], P = 0.40629.94 [9.83, 91.19], P < 0.001
Hazard ratio [95% confidence interval], P value
New-onset heart failureaCardiovascular mortalitybMyocardial infarctionc
Mean LDL-C1.00 [0.91, 1.10], P = 0.9461.37 [1.20, 1.55], P < 0.0010.94 [0.83, 1.07], P = 0.368
SD of LDL-C1.03 [0.85, 1.25], P = 0.7450.94 [0.72, 1.24], P = 0.6771.46 [1.15, 1.85], P = 0.002
CV of LDL-C0.63 [0.36, 1.11], P = 0.1090.19 [0.08, 0.42], P < 0.0013.64 [1.81, 7.33], P < 0.001
Mean HDL-C0.60 [0.50, 0.72], P < 0.0010.54 [0.41, 0.71], P < 0.0010.46 [0.36, 0.60], P < 0.001
SD of HDL-C4.40 [2.72, 7.13], P < 0.0011.05 [0.46, 2.39], P = 0.9012.82 [1.42, 5.59], P = 0.003
CV of HDL-C48.68 [21.37, 110.90], P < 0.0011.63 [0.45, 5.97], P = 0.46039.87 [13.77, 115.47], P < 0.001
Mean total cholesterol0.89 [0.82, 0.96], P = 0.0041.21 [1.08, 1.36], P = 0.0010.87 [0.78, 0.98], P = 0.017
SD of total cholesterol1.38 [1.18, 1.62], P < 0.0011.18 [0.93, 1.51], P = 0.1691.69 [1.39, 2.06], P < 0.001
CV of total cholesterol5.49 [2.29, 13.15], P < 0.0010.58 [0.16, 2.11], P = 0.40629.94 [9.83, 91.19], P < 0.001

CV, coefficient of variation; SD, standard deviation.

a

Adjusted for age, hypertension, diabetes mellitus, atrial fibrillation, chronic obstructive pulmonary disease, ischaemic heart disease, peripheral vascular disease, stroke or transient ischaemic attack, dementia or Alzheimer’s disease, and use of lipid-lowering medication(s).

b

Adjusted for age, sex, hypertension, diabetes mellitus, atrial fibrillation, ischaemic heart disease, peripheral vascular disease, stroke or transient ischaemic attack, dementia or Alzheimer’s disease, and use of lipid-lowering medication(s).

c

Adjusted for age, sex, hypertension, diabetes mellitus, chronic obstructive pulmonary disease, ischaemic heart disease, peripheral vascular disease, stroke or transient ischaemic attack, and use of lipid-lowering medication(s).

Subgroup analysis by use of lipid-lowering medication(s) (see Supplementary material online, Tables S6) found that in patients who did not use lipid-lowering medication(s) (N = 4068), associations between cholesterol measures and outcomes remained similar. Contrastingly, in patients who used lipid-lowering medication(s) (N = 1594), LDL-C and HDL-C variability were not prognostic. While lower mean HDL-C was associated with higher new-onset HF risk [aHR 0.43 (0.31, 0.61), P < 0.001], no association was observed for SD of HDL-C [aHR 1.829 (0.659, 5.078), P = 0.247], suggesting that the positive correlation between CV of HDL-C and new-onset HF risk was driven by inverse correlation for mean HDL-C. Meanwhile, HDL-C variability was not associated with MI risk, but higher total cholesterol variability remained associated with higher MI risk [aHR for CV 25.89 (3.74, 179.01), P = 0.001].

In this population-based study, we showed important but varying associations between visit-to-visit cholesterol variability and long-term risks of new-onset HF and major adverse cardiovascular outcomes, which may be partially negated by lipid-lowering medications. As far as we know, this is the first study demonstrating links between cholesterol variability and new-onset HF risk. Using population-based data with long follow-up, our results are widely applicable. Nonetheless, this study was limited by its retrospective nature, and the type of HF was unknown as echocardiographic data were unavailable. Additionally, many were excluded for having less than three sets of fasting lipid profiles during index period, limiting generalizability and possibly biasing for those indicated for lipid testing. Lastly, there may have been residual confounders, such as response to pharmacotherapy, smoking history, and blood pressure. Nonetheless, we adjusted for numerous cardiovascular risk factors which should account for most pertinent confounders.

Although underlying mechanisms remain unclear, cholesterol variability’s links to atherosclerotic progression6,7 and inflammation8,9 may be important, as these play important roles in HF.10 Clinicians should be aware of the importance of controlling cholesterol fluctuations, and further investigations should explore measures which reduce such fluctuations.

Authors’ contributions

J.S.K.C., D.I.S., and G.T. contributed to the conception or design of the work. T.T.L.L., O.H.I.C., A.K.C.W., B.M.Y.C., and G.T. contributed to the acquisition of data for the work. T.T.L.L., O.H.I.C., and J.Z. curated the data for the work. J.S.K.C. and J.Z. contributed to formal analysis and visualization for the work. J.S.K.C. and G.T. contributed to project administration. J.S.K.C. and D.I.S. drafted the manuscript. All authors critically revised the manuscript. All authors gave final approval and agree to be accountable for all aspects of work ensuring integrity and accuracy.

Supplementary material

Supplementary material is available at European Journal of Preventive Cardiology.

Funding

None declared.

Data availability

Underlying data are available upon reasonable request to the corresponding authors.

References

1

Yusuf
PS
,
Hawken
S
,
Ôunpuu
S
,
Dans
T
,
Avezum
A
,
Lanas
F
,
McQueen
M
,
Budaj
A
,
Pais
P
,
Varigos
J
,
Lisheng
L
.
Effect of potentially modifiable risk factors associated with myocardial infarction in 52 countries (the INTERHEART study): case-control study
.
Lancet
2004
;
364
:
937
952
.

2

Nazir
DJ
,
Roberts
RS
,
Hill
SA
,
McQueen
MJ
.
Monthly intra-individual variation in lipids over a 1-year period in 22 normal subjects
.
Clin Biochem
1999
;
32
:
381
389
.

3

Lee
S
,
Jeevaratnam
K
,
Liu
T
,
Chang
D
,
Chang
C
,
Wong
WT
,
Wong
ICK
,
Lip
GYH
,
Tse
G
.
Risk stratification of cardiac arrhythmias and sudden cardiac death in type 2 diabetes mellitus patients receiving insulin therapy: a population-based cohort study
.
Clin Cardiol
2021
;
44
:
1602
1612
.

4

Wang
M-C
,
Li
CI
,
Liu
CS
,
Lin
CH
,
Yang
SY
,
Li
TC
,
Lin
CC
.
Effect of blood lipid variability on mortality in patients with type 2 diabetes: a large single-center cohort study
.
Cardiovasc Diabetol
2021
;
20
:
228
.

5

Chan
JSK
,
Zhou
J
,
Lee
S
,
Li
A
,
Tan
M
,
Leung
KSK
,
Jeevaratnam
K
,
Liu
T
,
Roever
L
,
Liu
Y
,
Tse
G
,
Zhang
Q
.
Fragmented QRS is independently predictive of long-term adverse clinical outcomes in Asian patients hospitalized for heart failure: a retrospective cohort study
.
Front Cardiovasc Med
2021
;
8
:
738417
.

6

Abela
GS
,
Vedre
A
,
Janoudi
A
,
Huang
R
,
Durga
S
,
Tamhane
U
.
Effect of statins on cholesterol crystallization and atherosclerotic plaque stabilization
.
Am J Cardiol
2011
;
107
:
1710
1717
.

7

Navab
M
,
Reddy
ST
,
Van Lenten
BJ
,
Fogelman
AM
.
HDL and cardiovascular disease: atherogenic and atheroprotective mechanisms
.
Nat Rev Cardiol
2011
;
8
:
222
232
.

8

Zhao
L
,
Xu
T
,
Li
Y
,
Luan
Y
,
Lv
Q
,
Fu
G
,
Zhang
W
.
Variability in blood lipids affects the neutrophil to lymphocyte ratio in patients undergoing elective percutaneous coronary intervention: a retrospective study
.
Lipids Health Dis
2020
;
19
:
124
.

9

Lee
S
,
Zhou
J
,
Wong
WT
,
Liu
T
,
Wu
WKK
,
Wong
ICK
,
Zhang
Q
,
Tse
G
, et al.
Glycemic and lipid variability for predicting complications and mortality in diabetes mellitus using machine learning
.
BMC Endocr Disord
2021
;
21
:
94
.

10

Reina-Couto
M
,
Pereira-Terra
P
,
Quelhas-Santos
J
,
Silva-Pereira
C
,
Albino-Teixeira
A
,
Sousa
T
.
Inflammation in human heart failure: major mediators and therapeutic targets
.
Frontiers in Physiology
2021
;
12
:
1753
.

Author notes

J.S.K.C. and D.I.S. are co-first authors.

Part of this study has been presented as an oral presentation at ESC Asia 2021.

Conflict of interest: G.B.-Z. has consulted for Cardionovum, Crannmedical, Innovheart, Meditrial, Opsens Medical, Replycare, and Terumo. Other authors have nothing to disclose.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Supplementary data

Comments

0 Comments
Submit a comment
You have entered an invalid code
Thank you for submitting a comment on this article. Your comment will be reviewed and published at the journal's discretion. Please check for further notifications by email.