- Original investigation
- Open Access
Glucose variability and the risks of stroke, myocardial infarction, and all-cause mortality in individuals with diabetes: retrospective cohort study
Cardiovascular Diabetology volume 19, Article number: 144 (2020)
Previous research regarding long-term glucose variability over several years which is an emerging indicator of glycemic control in diabetes showed several limitations. We investigated whether variability in long-term fasting plasma glucose (FG) can predict the development of stroke, myocardial infarction (MI), and all-cause mortality in patients with diabetes.
This is a retrospective cohort study using the data provided by the Korean National Health Insurance Corporation. A total of 624,237 Koreans ≥ 20 years old with diabetes who had undergone health examinations at least twice from 2005 to 2008 and simultaneously more than once from 2009 to 2010 (baseline) without previous histories of stroke or MI. As a parameter of variability of FG, variability independent of mean (VIM) was calculated using FG levels measured at least three times during the 5 years until the baseline. Study endpoints were incident stroke, MI, and all-cause mortality through December 31, 2017.
During follow-up, 25,038 cases of stroke, 15,832 cases of MI, and 44,716 deaths were identified. As the quartile of FG VIM increased, the risk of clinical outcomes serially increased after adjustment for confounding factors including duration and medications of diabetes and the mean FG. Adjusted hazard ratios (95% confidence intervals) of FG VIM quartile 4 compared with quartile 1 were 1.20 (1.16–1.24), 1.20 (1.15–1.25), and 1.32 (1.29–1.36) for stroke, MI and all-cause mortality, respectively. The impact of FG variability was higher in the elderly and those with a longer duration of diabetes and lower FG levels.
In diabetes, long-term glucose variability showed a dose–response relationship with the risk of stroke, MI, and all-cause mortality in this nationwide observational study.
For individuals with diabetes, maintaining optimal blood glucose level is required to prevent complications and death . However, as incident cardiovascular disease (CVD) and deaths associated with diabetes cannot be fully explained by increased HbA1c itself [2,3,4], more attention has been focused on non-traditional risk factors such as glucose variability (GV). Patients with diabetes with similar HbA1c values can have different daily glycemic profiles with variations in the duration and frequency of glycemic excursions . Short-term GV usually means same-day or between-day glycemic oscillations and might be influenced by the individuals’ diet, exercise, and treatment modality. Long-term GV represents the visit-to-visit variability in fasting glucose (FG) estimated over months to years; compliance with medication; and deterioration of insulin secretion and resistance, which can be important issues .
During the past decade, detrimental effects of GV on patients with diabetes have been proposed for various medical conditions [7,8,9,10,11,12,13]. GV could not only predict diabetic vascular complications [7, 14, 15], heart failure , and postoperative complications of aortic valve implantation , but also indicate poor prognosis for in-patients with acute lung diseases  and acute coronary syndrome . Moreover, because GV could modify the correlation between time in range and estimated HbA1c, GV consideration was recommended when setting individualized goals for glycemic control .
Although there have been a few attempts to verify if GV can be a contributor to CVD or mortality [6, 7, 15,16,17,18,19,20,21], there is no consensus for patients with diabetes. Previous studies have yielded partial significance according to age, study group, or status of glycemic control [6, 7, 15,16,17,18,19,20,21]. In addition, they have several limitations such as small sample size and short-term follow-up periods.
In the present study, we used a nationwide population-based database in Korea to investigate whether long-term variability in FG can predict the risk of stroke, myocardial infarction (MI), and all-cause mortality during a median 8 years of follow-up, to specify the population prone to the risk of higher GV, and to compare it with single measurements of FG.
Study design and subjects
As shown in Fig. 1, we analyzed the data of patients with diabetes who had undergone health examinations provided by the National Health Insurance Corporation (NHIC) at least twice from 2005 to 2008 and simultaneously more than once between January 1, 2009, and December 31, 2010 (baseline exam). Among them, we excluded individuals aged < 20 years and those with histories of stroke or MI, along with those missing data in the inclusion criteria. After those exclusions, 624,237 individuals were included in this study.
The NHIC is the single health insurance system operated by the Korean government and covers about 97% of the Korean population. It has an eligibility profile, health examination database (general health examinations and questionnaires on lifestyle), medical treatment information identified by the medical bills submitted by healthcare providers, and a medical care institution database [22, 23]. Enrollees are advised to undergo a standardized medical examination annually or biannually. The entire database of the NHIC is now open to all researchers.
The presence of diabetes was defined as at least one claim per year for the prescription of anti-diabetic medications under the International Classification of Diseases, Tenth Revision (ICD-10) codes E10–14 or a fasting plasma glucose level ≥ 7 mm ol/L. Type 1 diabetes (T1DM) was defined as ICD-10 code E10 with at least one prescription history of insulin. The remained subjects were referred to as type 2 diabetes (T2DM).
This study was approved by the official review committee in the NHIC and the Institutional Review Board of the Korea University Ansan Hospital (Institutional Review Board number 2019AS0158) and was carried out in accordance with the Helsinki Declaration of 1975.
Definition of GV
The variability independent of the mean (VIM) of FG was used as a primary variability parameter calculated from FG levels measured at least three times during the 5 years prior to the baseline (Fig. 2). The FG level at baseline was included in calculation of FG variability. Additionally, standard deviation (SD), coefficient of variation (CV, SD/mean), and average real variability (ARV) of FG during the same period were calculated using the following equation:
where β is derived by nonlinear regression analysis, based on the natural logarithm of SD over the natural logarithm of the mean .
where k ranges from 1 to n-1, and n is the number of FG measurements.
Incidences of stroke, MI, or death until December 31, 2017, were the endpoints of this study. Stroke was defined as the recording of ICD-10 codes I63 or I64 during admission with simultaneous claims for brain CT or MRI. MI was diagnosed by the recording of ICD-10 codes I21 or I22 with hospitalization. Deceased cases were identified based on nationwide death certificate data from the Korea National Statistical Office regardless of a previous diagnosis of stroke or MI. The time interval between the baseline examination and the date of study outcome or December 31, 2017 is defined as the follow-up period.
Anthropometric and laboratory measurements
We used data from questionnaires regarding demographic characteristics, lifestyle, medical history, and medications during the medical examination. Alcohol consumption habit was divided into near abstinence, moderate (< 30 g/day), or severe (≥ 30 g/day). Smoking history was categorized as never, ex-, and current smokers. Doing regular exercise was defined > 20 min of vigorous-intensity or > 30 min of moderate-intensity exercise at least once per week . Body mass index (BMI) was estimated by weight in kilograms divided by the square of height in meters. Waist circumference was checked at the middle point between the iliac crest and the rib cage.
The presence of hypertension was defined as systolic blood pressure (BP) ≥ 140 mmHg, diastolic BP ≥ 90 mmHg, or the presence of at least one prescription of anti-hypertensive medications under ICD-10 codes I10–I15 per year. The presence of malignancy was determined by registration in the Korea Central Cancer Registry under the International Classification of Diseases, ICD-10 C00-C96 before baseline examination. Low-income status was defined by the lowest 20%.
Venous blood samples were drawn in the morning after an overnight fast of at least 8 hours to measure the levels of plasma glucose, total cholesterol, triglycerides, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol (LDL-C), and creatinine.
Dyslipidemia was defined as total cholesterol level ≥ 6.21 mmol/L or the presence of at least one claim per year for the prescription of anti-hyperlipidemic drugs under ICD-10 code E78. The estimated glomerular filtration rate (eGFR) was estimated by the Modification of Diet in Renal Disease formula , and eGFR < 60 mL/min/1.73 m2 was classified as chronic kidney disease (CKD) . The number of oral anti-diabetic medication among metformin, sulfonylurea, meglitinide, thiazolidinedione, inhibitors of dipeptidyl peptidase 4 (DPP-4 inhibitors), and acarbose taken in the 12 months prior to baseline was identified. The history of heart disease was identified by self-report.
Quality control of the laboratory tests was conducted in accordance with the Korean Association of Laboratory Quality Control.
Data are presented as mean ± SD, geometric mean (95% confidence intervals [CIs]), or number (%). The baseline characteristics were compared using Chi square tests for categorical variables and analysis of variance for continuous variables after dividing the subjects according to the FG VIM quartile. Triglyceride levels were log-transformed for analysis.
To assess the risks of stroke, MI, and all-cause mortality, we conducted multivariate-adjusted Cox proportional hazards analyses according to FG VIM quartiles and deciles, using quartile 1 or decile 1 as the reference group.
We adjusted for confounders at baseline using two models. Model 1 was adjusted for age, sex, BMI, alcohol drinking, smoking, regular exercise, presence of hypertension and dyslipidemia, CKD, and low-income status. Model 2 is the same as model 1, plus further adjustment for duration of diabetes over 5 years, the number of classes of oral anti-diabetic medication, insulin prescription history taken in the 12 months prior to baseline, and mean FG during the 5 years preceding the baseline exam.
We conducted several subgroup analyses to evaluate the effects of age; sex; BMI; smoking; income status; presence of hypertension, dyslipidemia, and malignancy; duration of diabetes; baseline FG level; subtype of diabetes; and anti-diabetic medication. The hazard ratios (HRs) and 95% CIs for stroke, MI, and all-cause death in FG VIM quartile 4 versus quartiles 1–3, with adjustments for confounders, were estimated.
To compare with the predictive value of a single FG level, we performed the above-mentioned Cox analysis according to baseline FG level, with FG levels of 100-119 mg/dL being the reference group. The mean of FG was excluded as a confounder in this analysis. Additionally, we repeated this analysis after stratifying the subjects according to prescription history of anti-diabetic medication.
We checked a variable inflation factor for all covariates of less than 2.0 and found no relevant multicollinearity in covariates. SAS version 9.3 (SAS Institute Inc., Cary, NC, USA) was used for statistical analysis. A p value of < 0.05 was considered to be statistically significant.
FG VIM quartile 4 had more males and current smokers and showed higher FG and triglyceride levels compared with quartile 1 (Table 1). On the other hand, this group was younger, had a lower BP and LDL-C level, and had lower proportions of hypertension, dyslipidemia, and heart disease versus quartile 1. About 46% of the highest VIM quartile was not treated with any anti-hyperglycemic agents. However, a higher proportion of them was treated with insulin for the preceding 12 months than other tertiles.
During the median (interquartile range) follow-up period of 8.0 (7.3–8.4) years, 25,038 cases of stroke, 15,832 cases of MI, and 44,716 cases of death were identified. As shown in Table 2, the HRs for stroke, MI, and all-cause mortality serially increased as the FG VIM quartile increased. Subjects in FG VIM quartile 4 had a 20% higher risk for stroke and MI, and a 32% higher risk for all-cause mortality, versus those in quartile 1, even after adjustment for several risk factors and mean FG. We also tested the following GV calculation methods: SD, CV, and ARV. A similar relationship was obtained in FG SD, CV, and ARV quartiles instead of VIM (Additional file 1: Table S1).
When we divided individuals into FG VIM deciles, the risks of stroke, MI, and all-cause death showed a positive dose–response association with the FG VIM decile (Fig. 3 and Additional file 1: Table S2).
In the subgroup analyses, FG VIM quartile 4 showed a consistently increased risk for all-cause mortality (Tables 3 and 4) except in young individuals aged 20-39. A similar finding was observed for stroke and MI, with exceptions in younger individuals and those with malignancy, shorter diabetes duration, baseline FG ≥ 126 mg/dL, T1DM, and treatment with thiazolidinedione and insulin.
On the other hand, the relationship between FG status and the risk of clinical outcomes showed a U-shaped association (Fig. 4 and Additional file 1: Table S3). Compared to individuals whose FG levels were in the range 100-119 mg/dL, those with FG < 100 mg/dL or ≥ 160 mg/dL had significantly higher HRs. This relationship was consistent after dividing the subjects by prescription of anti-diabetic medications (Additional file 1: Table S4).
This nationwide population-based study in diabetes revealed that long-term GV for 5 years predicted future development of stroke, MI, and all-cause death, independent of anti-diabetic medication, metabolic risk factors, and mean fasting glucose. This association showed a linear fashion, compared to the U-shaped association between FG and outcomes. In the subgroup analysis, the impact of GV was higher in the elderly and those with a longer duration of diabetes and lower FG levels.
Long-term glucose variability and clinical outcomes
Several medium-sized cohort studies have reported the impact of long-term GV [6, 15, 17,18,19,20,21]. Variability in FG is an independent predictor of all-cause mortality, and the highest tertile group in 1400 type 2 diabetes patients included in the VERONA study had a 67% higher risk [17, 21]. In the sub-analysis of the Action in Diabetes and Vascular Disease: Preterax and Diamicron MR Controlled Evaluation (ADVANCE) trial, variability in FG was associated with an increased risk of death from cardiovascular causes, nonfatal MI, or nonfatal stroke . Furthermore, GV might reflect better the impact of glycemic control on diabetic complications than HbA1c, or at least, GV can complement the power of HbA1c in risk prediction [7, 14, 15]. In the ADVANCE trial, HRs (95% CIs) estimated for each 10-percentile point increase in SD of FG and HbA1c were 1.12 (1.08–1.16) and 1.05 (1.01–1.09), respectively .
However, these studies have limitations such as short follow-up period  and a relatively small population size, and the impact of GV was only shown in some subgroups, but not the entire population . In the present study, we computed long-term GV over 5 years and assessed risk for a median of 8 years of follow-up. In addition, this dataset includes almost all Korean diabetic patients. In this respect, the worst outcome in the highest GV quartile even after adjustment for mean FG provides a clear indication that GV can be considered as a new risk factor for stroke, MI, and all-cause mortality.
Interpretation of the impact of glucose variability
A few mechanisms might be involved in the relationship between FG variability and clinical outcomes. Transient high glucose spikes have been shown to impair endothelial function [28, 29], increase oxidative stress more than sustained chronic hyperglycemia [30, 31], and induce β-cell dysfunction . Furthermore, higher GV was strongly associated with increased lipid level and decreased fibrous content with a larger plaque burden .
In subgroup analysis of this study, the impact of GV was greater in the elderly and those with a longer duration of diabetes and FG level < 126 mg/dL. These findings were partially in accordance with previous publications [18, 20]. In the Verona Diabetes Study, the positive association between GV and mortality risk was confined to patients older than 65 years . The elderly are known to be more susceptible to oxidative stress than younger individuals because their defense mechanisms are less efficient . The limited significance seen in current smokers could be explained by the higher proportion of young subjects in the current smoker group. In the same notion, patients with a longer diabetes duration and lower FG might reflect a population that is more vulnerable to GV. It is interesting that contrary to other studies [15, 19, 20], GV had a greater impact among lower FG patients than in higher FG patients, which may be due to an ethnic difference, but the reason is difficult to ascertain from this study. As we also showed that lower FG level was significantly associated with more CVD and death events, the higher GV in the lower FG group augmented the risk for future CVD and death, indicating the target population for reducing these outcomes.
In addition to this pathophysiologic interpretation, elevated GV might be an indicator of irregular compliance to therapy, comorbidity, poor health, or diabetic complications resulting in the increase of mortality . However, contrary to other studies [14, 17, 20], our study population with higher GV had a shorter duration of diabetes, had fewer comorbidities, and about half of them were not treated with anti-diabetic medication. Although there was a greater proportion of insulin users, the contribution of insulin users to the outcome was small since less than 10% of the study population was treated with insulin. We speculate that, compared to other hospital-based studies, this study included more low-risk diabetic patients since diagnosis of diabetes was based on not only health insurance claims data, but also general health screening exams, which resulted in detection of undiagnosed mild diabetic cases. Therefore, we can evaluate the impact of GV itself, rather than combined poor health status and comorbidities as in other studies.
In subgroup analysis according to anti-diabetic medications, interaction analysis did not find any significance, contrary to previous studies which showed sulfonylurea increases glucose fluctuation and risk of hypoglycemia , and DPP-4 inhibitors and the novel insulin analogue degludec reduced GV [34,35,36]. Given that actual drug exposure could not be determined because of the retrospective study design, a future well-designed, randomized controlled trial is expected to resolve this concern.
Assessment of glycemic variability
There are no standard indices for quantification of long-term glycemic variability , so each value has distinct characteristics. The SD reflects dispersion of measurements around the mean and is sensitive to low sampling frequency, while CV is a standardized variation providing direct comparison among study groups. The ARV index, which averages the absolute differences between successive measurements, might be a reliable index for time series variability [38, 39]. However, SD, CV, and ASV are partially dependent on the mean and its changes over time, and this may not be resolved even if adjusted for mean value . On the other hand, VIM is a measure of variability designed not to correlate with mean level [24, 41] but is sample-specific [42, 43]. Therefore, we used VIM as the main measurement of FG variability, and similar results to other estimates of glycemic variation support the robustness of our study (Additional file 1: Table S1).
Limitations of this study
This is the first large-scale epidemiologic study demonstrating the impact of long-term variability in FG on CVD outcomes separately and all-cause death with a long-term follow-up of 8.0 (7.3–8.4) years. However, we are aware of several limitations of this study. First, given that HbA1c or postprandial glucose level was not measured in this study, incident diabetes might be underestimated. To enhance the accuracy of diabetes diagnosis, we combined ICD-10 codes and patient prescription histories in addition to fasting glucose level. Second, possible selection bias due to extracting patients based on the number of health check-ups is one of the limitations. To overcome it, we conducted various subgroup analyses. Finally, classification of T1DM based on the presence of ICD-10 E10 with at least one prescription of insulin might not be accurate. Because ICD-10 code is based on claims data, but autoantibody, insulin, or c-peptide level was not measured, the prevalence of T1DM could be overestimated. However, we performed subgroup analysis after stratifying the subjects according to subtype of diabetes. There was a significance in all-cause mortality in T1DM. The absence of significance in stroke and MI was thought to be related to the small number of population.
This nationwide population-based study with long-term follow-up showed that GV had a dose–response relationship with the risk of stroke, MI, and all-cause mortality in diabetes, especially in the elderly and those with a longer duration of diabetes and lower FG levels. This important health issue will play a role in reducing future development of CVD and death attributed to the increasing population of diabetes.
Availability of data and materials
The data that support the findings of this study are available from the National Health Insurance Corporation but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of the National Health Insurance Corporation.
Average real variability
Body mass index
Chronic kidney disease
Coefficient of variation
- DPP-4 inhibitors:
Inhibitors of dipeptidyl peptidase 4
Estimated glomerular filtration rate
International Classification of Diseases, tenth revision
low-density lipoprotein cholesterol
National Health Insurance System
Type 1 diabetes
Type 2 diabetes
Variability independent of the mean
6 Glycemic Targets: standards of medical care in diabetes-2018. Diabetes Care. 2018;41:S55-S64.
Wan EYF, Yu EYT, Fung CSC, Chin WY, Fong DYT, Chan AKC, et al. Relation between HbA1c and incident cardiovascular disease over a period of 6 years in the Hong Kong population. Diabetes Metab. 2018;44:415–23.
Arem H, Moore SC, Patel A, Hartge P, de Gonzalez BA, Visvanathan K, et al. Leisure time physical activity and mortality: a detailed pooled analysis of the dose-response relationship. JAMA Intern Med. 2015;175:959–67.
Škrha J, Šoupal J, Škrha J Jr, Prázný M. Glucose variability, HbA1c and microvascular complications. Rev Endocrine Metab Disorders. 2016;17:103–10.
Siegelaar SE, Holleman F, Hoekstra JBL, DeVries JH. Glucose variability; does it matter? Endocr Rev. 2010;31:171–82.
Sun B, He F, Gao Y, Zhou J, Sun L, Liu R, et al. Prognostic impact of visit-to-visit glycemic variability on the risks of major adverse cardiovascular outcomes and hypoglycemia in patients with different glycemic control and type 2 diabetes. Endocrine. 2019;64:536–43.
Zhou JJ, Schwenke DC, Bahn G, Reaven P. Glycemic variation and cardiovascular risk in the Veterans Affairs Diabetes trial. Diabetes Care. 2018;41:2187–94.
Ferreira L, Moniz AC, Carneiro AS, Miranda AS, Fangueiro C, Fernandes D, et al. The impact of glycemic variability on length of stay and mortality in diabetic patients admitted with community-acquired pneumonia or chronic obstructive pulmonary disease. Diabetes Metab Syndr. 2019;13:149–53.
Yokota S, Tanaka H, Mochizuki Y, Soga F, Yamashita K, Tanaka Y, et al. Association of glycemic variability with left ventricular diastolic function in type 2 diabetes mellitus. Cardiovasc Diabetol. 2019;18:166.
Lu J, Ma X, Zhang L, Mo Y, Lu W, Zhu W, et al. Glycemic variability modifies the relationship between time in range and hemoglobin A1c estimated from continuous glucose monitoring: a preliminary study. Diabetes Res Clin Pract. 2020;161:108032.
Besch G, Pili-Floury S, Morel C, Gilard M, Flicoteaux G, du Mont SL, et al. Impact of post-procedural glycemic variability on cardiovascular morbidity and mortality after transcatheter aortic valve implantation: a post hoc cohort analysis. Cardiovasc Diabetol. 2019;18:27.
Zhang Y, Dai J, Han X, Zhao Y, Zhang H, Liu X, et al. Glycemic variability indices determined by self-monitoring of blood glucose are associated with β-cell function in Chinese patients with type 2 diabetes. Diabetes Res Clin Pract. 2020;164:108152.
Takahashi H, Iwahashi N, Kirigaya J, Kataoka S, Minamimoto Y, Gohbara M, et al. Glycemic variability determined with a continuous glucose monitoring system can predict prognosis after acute coronary syndrome. Cardiovasc Diabetol. 2018;17:116.
Cardoso CRL, Leite NC, Moram CBM, Salles GF. Long-term visit-to-visit glycemic variability as predictor of micro- and macrovascular complications in patients with type 2 diabetes: the Rio de Janeiro Type 2 Diabetes Cohort Study. Cardiovasc Diabetol. 2018;17:33.
Hirakawa Y, Arima H, Zoungas S, Ninomiya T, Cooper M, Hamet P, et al. Impact of visit-to-visit glycemic variability on the risks of macrovascular and microvascular events and all-cause mortality in type 2 diabetes: the ADVANCE trial. Diabetes Care. 2014;37:2359–65.
Tang X, Li S, Wang Y, Wang M, Yin Q, Mu P, et al. Glycemic variability evaluated by continuous glucose monitoring system is associated with the 10-y cardiovascular risk of diabetic patients with well-controlled HbA1c. Clin Chim Acta. 2016;461:146–50.
Muggeo M, Zoppini G, Bonora E, Brun E, Bonadonna RC, Moghetti P, et al. Fasting plasma glucose variability predicts 10-year survival of type 2 diabetic patients: the Verona Diabetes Study. Diabetes Care. 2000;23:45–50.
Zoppini G, Verlato G, Targher G, Bonora E, Trombetta M, Muggeo M. Variability of body weight, pulse pressure and glycaemia strongly predict total mortality in elderly type 2 diabetic patients. The Verona Diabetes Study. Diabetes Metab Res Rev. 2008;24:624–8.
Rizvi SI, Maurya PK. Markers of oxidative stress in erythrocytes during aging in humans. Ann N Y Acad Sci. 2007;1100:373–82.
Xu D, Fang H, Xu W, Yan Y, Liu Y, Yao B. Fasting plasma glucose variability and all-cause mortality among type 2 diabetes patients: a dynamic cohort study in Shanghai, China. Sci Rep. 2016;6:39633.
Muggeo M, Verlato G, Bonora E, Zoppini G, Corbellini M, de Marco R. Long-term instability of fasting plasma glucose, a novel predictor of cardiovascular mortality in elderly patients with non-insulin-dependent diabetes mellitus: the Verona Diabetes Study. Circulation. 1997;96:1750–4.
Song SO, Jung CH, Song YD, Park CY, Kwon HS, Cha BS, et al. Background and data configuration process of a nationwide population-based study using the korean national health insurance system. Diabetes Metab J. 2014;38:395–403.
Lee YH, Han K, Ko SH, Ko KS, Lee KU. Data analytic process of a nationwide population-based study using National Health Information Database established by National Health Insurance Service. Diabetes Metab J. 2016;40:79–82.
Rothwell PM, Howard SC, Dolan E, O’Brien E, Dobson JE, Dahlof B, et al. Prognostic significance of visit-to-visit variability, maximum systolic blood pressure, and episodic hypertension. Lancet. 2010;375:895–905.
Oh JY, Yang YJ, Kim BS, Kang JH. Validity and reliability of Korean version of International Physical Activity Questionnaire (IPAQ) short form. J Korean Acad Fam Med. 2007;28:532–41.
Levey AS, Coresh J, Greene T, Stevens LA, Zhang YL, Hendriksen S, et al. Using standardized serum creatinine values in the modification of diet in renal disease study equation for estimating glomerular filtration rate. Ann Intern Med. 2006;145:247–54.
Group KDIGOCW. KDIGO 2012 clinical practice guideline for the evaluation and management of chronic kidney disease: Chapter 1: definition and classification of CKD. Kidney Int Suppl. 2013;2013(3):19–62.
Zhang XG, Zhang YQ, Zhao DK, Wu JX, Zhao J, Jiao XM, et al. Relationship between blood glucose fluctuation and macrovascular endothelial dysfunction in type 2 diabetic patients with coronary heart disease. Eur Rev Med Pharm Sci. 2014;18:3593–600.
Ceriello A, Esposito K, Piconi L, Ihnat MA, Thorpe JE, Testa R, et al. Oscillating glucose is more deleterious to endothelial function and oxidative stress than mean glucose in normal and type 2 diabetic patients. Diabetes. 2008;57:1349–54.
Salisbury D, Bronas U. Reactive oxygen and nitrogen species: impact on endothelial dysfunction. Nurs Res. 2015;64:53–66.
El-Osta A, Brasacchio D, Yao D, Pocai A, Jones PL, Roeder RG, et al. Transient high glucose causes persistent epigenetic changes and altered gene expression during subsequent normoglycemia. J Exp Med. 2008;205:2409–17.
Okada K, Hibi K, Gohbara M, Kataoka S, Takano K, Akiyama E, et al. Association between blood glucose variability and coronary plaque instability in patients with acute coronary syndromes. Cardiovasc Diabetol. 2015;14:111.
Yamazaki M, Hasegawa G, Majima S, Mitsuhashi K, Fukuda T, Iwase H, et al. Effect of repaglinide versus glimepiride on daily blood glucose variability and changes in blood inflammatory and oxidative stress markers. Diabetol Metab Syndr. 2014;6:54.
Suh S, Kim H, Dang-Vu TT, Joo E, Shin C. Cortical thinning and altered cortico-cortical structural covariance of the default mode network in patients with persistent insomnia symptoms. Sleep. 2016;39:161–71.
Heinemann L, Fleming GA, Petrie JR, Holl RW, Bergenstal RM, Peters AL. Insulin pump risks and benefits: a clinical appraisal of pump safety standards, adverse event reporting, and research needs. Diabetes Care. 2015;38:716–22.
Vora J, Cariou B, Evans M, Gross JL, Harris S, Landstedt-Hallin L, et al. Clinical use of insulin degludec. Diabetes Res Clin Pract. 2015;109:19–31.
Inzucchi SE, Umpierrez G, DiGenio A, Zhou R, Kovatchev B. How well do glucose variability measures predict patient glycaemic outcomes during treatment intensification in type 2 diabetes? Diabetes Res Clin Pract. 2015;108:179–86.
Hansen TW, Thijs L, Li Y, Boggia J, Kikuya M, Björklund-Bodegård K, et al. Prognostic value of reading-to-reading blood pressure variability over 24 hours in 8938 subjects from 11 populations. Hypertension. 2010;55:1049–57.
Mena L, Pintos S, Queipo NV, Aizpúrua JA, Maestre G, Sulbarán T. A reliable index for the prognostic significance of blood pressure variability. J Hypertens. 2005;23:505–11.
Lee HJ, Choi EK, Han KD, Lee E, Moon I, Lee SR, et al. Bodyweight fluctuation is associated with increased risk of incident atrial fibrillation. Heart Rhythm. 2020;17:365–71.
Yano Y. Visit-to-visit blood pressure variability-what is the current challenge? Am J Hypertens. 2017;30:112–4.
Bangalore S, Breazna A, DeMicco DA, Wun CC, Messerli FH. Visit-to-visit low-density lipoprotein cholesterol variability and risk of cardiovascular outcomes: insights from the TNT trial. J Am Coll Cardiol. 2015;65:1539–48.
Dolan E, O’Brien E. Blood pressure variability: clarity for clinical practice. Hypertension. 2010;56:179–81.
This study was supported by the Bio & Medical Technology Development Program of the National Research Foundation (NRF) funded by the Korean government (MSIT) (NRF-2019M3E5D3073102 and NRF-2019R1H1A2039682) and supported by Korea University Ansan Hospital Grant (O1903661). However, funders did not participate in the study’s design or reporting.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Hazard ratios (HRs) and 95% confidence intervals (CIs) for the incidence of stroke, myocardial infarction, and all-cause mortality by quartiles of fasting glucose variability, assessed by standard deviation, coefficient of variation, and average real variability. Table S2. Hazard ratios (HRs) and 95% confidence intervals (CIs) for the incidence of stroke, myocardial infarction, and all-cause mortality by deciles of fasting glucose variability. Table S3. Hazard ratios (HRs) and 95% confidence intervals (CIs) for the incidence of stroke, myocardial infarction, and all-cause mortality according to baseline fasting glucose level. Table S4. Hazard ratios (HRs) and 95% confidence intervals (CIs) according to baseline fasting glucose level and antidiabetic medication (ADM).
About this article
Cite this article
Lee, D.Y., Han, K., Park, S. et al. Glucose variability and the risks of stroke, myocardial infarction, and all-cause mortality in individuals with diabetes: retrospective cohort study. Cardiovasc Diabetol 19, 144 (2020). https://doi.org/10.1186/s12933-020-01134-0
- Diabetes mellitus Glucose variability
- Cardiovascular disease
- All-cause mortality
- The Korean National Health Insurance Corporation