Temporal Changes in Resting Heart Rate and Risk of Diabetes Mellitus

Article information

Diabetes Metab J. 2024;48(4):752-762
Publication date (electronic) : 2024 February 2
doi : https://doi.org/10.4093/dmj.2023.0305
1Division of Population Health Research, Department of Precision Medicine, Korea National Institute of Health, Cheongju, Korea
2Korea National Institute of Health, Cheongju, Korea
Corresponding author: Hyun-Young Park https://orcid.org/0000-0002-6698-7368 Korea National Institute of Health, 187 Osongsaengmyeong 2-ro, Heungdeok-gu, Cheongju 28159, Korea E-mail: hypark65@korea.kr
Received 2023 August 29; Accepted 2023 November 13.

Abstract

Background

To investigate the association between the time-varying resting heart rate (RHR) and change in RHR (∆RHR) over time and the risk of diabetes mellitus (DM) by sex.

Methods

We assessed 8,313 participants without DM or atrial fibrillation/flutter from the Korean Genome and Epidemiology Study, a community-based prospective cohort study that was initiated in 2001 to 2002. The participants were followed up until December 31, 2018. Updating RHR with biennial in-study re-examinations, the time-varying ∆RHR was calculated by assessing the ∆RHR at the next follow-up visit.

Results

Over a median follow-up of 12.3 years, 1,345 participants (16.2%) had DM. As compared with RHR of 60 to 69 bpm, for RHR of ≥80 bpm, the incidence of DM was significantly increased for both male and female. A drop of ≥5 bpm in ∆RHR when compared with the stable ∆RHR group (–5< ∆RHR <5 bpm) was associated significantly with lower risk of DM in both male and female. However, an increase of ≥5 bpm in ∆RHR was significantly associated with higher risk of DM only in female, not in male (hazard ratio for male, 1.057 [95% confidence interval, 0.869 to 1.285]; and for female, 1.218 [95% confidence interval, 1.008 to 1.471]).

Conclusion

In this community-based longitudinal cohort study, a reduction in ∆RHR was associated with a decreased risk of DM, while an increase in ∆RHR was associated with an increased risk of DM only in female.

GRAPHICAL ABSTRACT

Highlights

• Time-varying RHR ≥80 bpm was found to be strongly related to the risk of incident DM.

• A decrease in RHR of ≥5 bpm was associated with a lower risk of developing DM.

• A decrease in RHR was significantly associated with a reduced risk of DM only in men.

• An increase in RHR was significantly associated with a increased risk of DM only in women.

• RHR may be a useful and easily measured biomarker for DM risk detection.

INTRODUCTION

Diabetes mellitus (DM) is a leading cause of disability and mortality [1,2]. It is a public health problem worldwide, with prevalence in adults aged 20 to 79 years expected to approach 9.3% (approximately 578 million) by 2030, and is projected to be the seventh leading cause of death by 2030 [3,4]. Thus, identification of groups at high risk for developing DM is urgently needed to better understand and prevent this public health epidemic.

Resting heart rate (RHR) is regulated by the autonomic nervous system (ANS) [5]. It is now possible to continuously and accurately self-measure RHR using a mobile phone or watch bracelet, such that monitoring RHR has become very popular in the general population [6]. Therefore, there is an opportunity for RHR to be used for chronic disease prevention in healthy adults if the clinical significance of RHR can be understood. Many previous studies have reported that a high RHR is associated with an increased risk of DM [7-9]. RHR may change over time due to environmental factors such as physical activity or clinical conditions; however, few studies have assessed serial measurements of RHR over time [10]. Furthermore, RHR is slightly higher in female than in male [11]. However, only a few previous studies have explored the correlation of RHR with DM by sex [12]. Thus, the association between RHR and DM is not fully understood in the general population.

Therefore, the objective of this analysis was to determine whether RHR and temporal changes in RHR (∆RHR) are useful risk markers for the development of DM in individuals who do not have DM. We also assessed by sex the association between DM and time-varying RHR and ∆RHR using repeated RHR measurements.

METHODS

Data source and study population

The Ansan-Ansung Cohort Study within the Korean Genome and Epidemiology Study (KoGES) is an ongoing prospective cohort started in 2001 to 2002. Participants are followed biennially through face-to-face visits, where information on sociodemographics, lifestyle factors, medical history, and anthropometric measures as well as blood and urine tests are collected by trained interviewers and staff using standard methods and protocols [13].

Initially, the cohort comprised 10,030 participants (5,010 from a farming community in the Ansung cohort and 5,020 from an industrial community in the Ansan cohort) aged 40 to 69 years old between 2001 and 2002. Biennial follow-up examinations were conducted in 2018. All the participants provided written informed consent. The study protocol was approved by the Institutional Review Board of Korea Centers for Disease Control and Prevention (approval No. 2021-04-02-2C-A). All research procedures were performed in accordance with relevant guidelines and regulations.

Participants were excluded if they had DM at baseline (n=826), atrial fibrillation or flutter at baseline (n=40), were lost to follow-up (n=840), or were missing heart rate information from the baseline examination (n=11). Finally, this study included 8,313 participants.

RHR and calculation of changes in RHR

At every visit, the RHR was obtained, after participants had rested in the supine position for 10 minutes, by well-trained specialists who measured the radial artery pulse over 1 minute using an automatic blood pressure monitor (FT-500, Jawon Medical Co. Ltd., Gyeongsan, Korea). As the RHR can vary over time, we created a time-varying RHR using each available study visit over the course of the study, so that all longitudinal study visit data were incorporated. That is, we used data from visits one to eight such that the RHR was updated up to seven times after the baseline for each participant. To assess ∆RHR, we calculated ∆RHR to the next follow-up visit by subtracting the time-varying RHR value at each biennial visit from the value at the next follow-up visit (Supplementary Fig. 1). Patients missing RHR information at each visit and those whose RHR measurements were from a visit when atrial fibrillation/flutter was electrocardiographically confirmed were not included in the analysis.

Outcome

The primary outcome was the incidence of DM, defined as a fasting serum glucose concentration of ≥126 mg/dL, post 2-hour glucose after a 75 g oral glucose tolerance test of ≥200 mg/dL, glycosylated hemoglobin ≥6.5%, current antidiabetic treatment, or DM-related death. Mortality data for the participants were obtained from the National Death Records. The participants were followed until the index date—the date of DM incidence, DM-related death, or the last interview before the end of the study period on December 31, 2018, whichever came first.

Statistical analysis

The characteristics of the study population were described as frequency (percentage) for categorical variables and mean±standard deviation (SD) for continuous variables. The incidence rate of DM was reported as the number of patients per 1,000 person-years.

To assess the association of longitudinal RHR and ∆RHR measures with incident DM, the time-varying Cox proportional hazards model, which is an extension of the widely used traditional Cox proportional hazards model, was used (Supplementary Tables 1 and 2, Supplementary Fig. 1). Although this method is usually used to focus on short-term effects, it provides better accuracy for assessing the relationship between RHR and DM than use of a single measurement that is fixed at the baseline [14,15]. Analyses of the association between time-varying RHR and time-varying ∆RHR with the risk of DM were adjusted for all variables in Table 1 and for time-varying ∆RHR and time-varying RHR, respectively. All variables in Table 1 could be treated as time-varying since they were assessed and obtained at every visit. The results are presented as hazard ratios (HRs) with 95% confidence intervals (CIs), using subjects with an RHR of 60 to <70 bpm as a reference. In addition, an adjusted model using a restricted cubic spline with five knots was constructed to display the association between the hazards of developing the outcome and the continuous covariates of the time-varying RHR using a reference value of 60 bpm, as previously reported [16-18]. For ∆RHR, 0 was used as the reference.

Characteristics of participants at visits 1–8

Mixed-effects regression models were used for RHR and ∆RHR trajectories in female and male, respectively, with time before incident DM or last examination used as the common timescale for all analyses. We used repeated RHR measures to fit the mixed-effects linear regression models, with each RHR measure as the outcome, participant identification as random intercepts, and time as a fixed effect expressed using restricted cubic splines with four knots to allow for nonlinear relationships. The mean values were estimated for each RHR measure over the time range of –15 to 0 years for female and male. Mixed-effects regression models were adjusted for the covariates listed in Table 1.

All statistical tests were 2-sided and P<0.05 was accepted as statistical significance. All analyses were performed using SAS version 9.4 (SAS Institute, Cary, NC, USA).

RESULTS

The characteristics of the study participants during visits 1–8 are presented in Table 1. The mean age of the cohort at baseline (visit 1) was 51.9 years, and 4,397 (52.9%) participants were female. During the follow-up period (mean, 12.3 years; range, 1.4 to 16.6 years), there was an increase in prevalence of hypertension (29.3% to 41.6%), cardiovascular disease (CVD; 2.7% to 5.0%), and chronic kidney disease (6.8% to 12.0%). Of the 8,313 subjects without DM, 1,345 (16.2%) developed DM. Their starting mean±SD value of RHR was 63.5±7.3 bpm compared to 62.8±7.4 bpm at the last visit.

Trajectory of RHR and ∆RHR according to the incidence of DM

According to the distribution of RHR from following visit, both male and female showed that RHR values were more frequently observed in individuals with incidence of DM (Supplementary Fig. 2A). The ∆RHR appeared to exhibit a similar distribution between male with and without incident DM, while in female, a slightly greater increase in RHR was observed in those with incident DM (Supplementary Fig. 2B). Fig. 1 shows the estimated trajectories of RHR and ∆RHR before the incidence of DM or last examination. The estimated mean RHR from 15 years before the incidence of DM was approximately 1 bpm higher in male and 1.5 bpm higher in female with incident DM compared to those without DM (Fig. 1A and B). However, the predicted mean ∆RHR in subjects with and without incident DM was similar both in male and female, respectively (Fig. 1C and D).

Fig. 1.

Trajectories of resting heart rate (RHR) before incident diabetes mellitus (DM) or last examination. Fitted lines denote mean RHR (incident DM vs. without incident DM), calculated from restricted cubic spline of mixed-effects linear regression models. Time 0=diagnosis or last examination. The dashed lines delimit the 95% confidence interval. (A) Time-varying RHR and DM in male. (B) Time-varying RHR and DM in female. (C) Time-varying change in RHR (∆RHR) and DM in male. (D) Time-varying ∆RHR and DM in female.

Association of time-varying RHR and ∆RHR with incidence of DM by sex

In fully adjusted time-varying Cox analysis (Table 2, Supplementary Tables 3 and 4), compared to participants with time-varying RHR 60 to 69 bpm, those with time-varying RHR ≥80 bpm had higher risk of DM both in male and female (HR in male and female, 2.195 [95% CI, 1.514 to 3.183] and 2.224 [95% CI, 1.562 to 3.165], respectively), and those with time-varying RHR <60 bpm had lower risk of DM in male, but not in female (HR in male and female, 0.816 [95% CI, 0.667 to 0.997] and 0.819 [95% CI, 0.665 to 1.010], respectively). In subgroup analysis, these patterns were significant in male without antihypertensive drug use and CVD, and in female with and without antihypertensive drug use, without CVD and who had body mass index (BMI) <23 kg/m2 (Fig. 2). Whereas, time-varying RHR <60 bpm also had significantly lower risk of DM in female without antihypertensive drug use and who had BMI ≥23 kg/m2.

Risk of diabetes mellitus by time-varying RHR and ΔRHR, time-varying Cox regression analysis

Fig. 2.

Subgroup analysis for the risk of diabetes mellitus according to the resting heart rate (RHR). Adjusted for sex, area, and time-updated age, body mass index (BMI), physical activity, smoking, drinking, systolic blood pressure, antihypertensive drug use, chronic kidney disease, cardiovascular disease (CVD), glycosylated hemoglobin, total cholesterol, and change in RHR. HR, hazard ratio; CI, confidence interval.

Compared with a stable change in RHR (∆RHR <5 bpm change), a decrease in RHR more than 5 bpm was associated with lower DM risk in both male and female, whereas an increase in RHR more than 5 bpm was associated with higher DM risk only in female (HR for ∆RHR ≥5 bpm drop in male and female, 0.602 [95% CI, 0.471 to 0.770] and 0.806 [95% CI, 0.651 to 0.998], respectively; HR for ∆RHR ≥5 bpm increase in female, 1.218 [95% CI, 1.008 to 1.471]). In subgroup analysis, a ∆RHR ≥5 bpm drop was associated with lower DM risk in male without antihypertensive drug use and CVD and who had BMI <23 kg/m2, and in female without CVD and who had BMI ≥23 kg/m2. A ∆RHR ≥5 bpm increase was associated with higher DM risk in female without antihypertensive drug use and CVD, and who had BMI ≥23 kg/m2 (Fig. 3).

Fig. 3.

Subgroup analysis for the risk of diabetes mellitus according to the change in resting heart rate (∆RHR). Adjusted for sex, area, and time-updated age, body mass index (BMI), physical activity, smoking, drinking, systolic blood pressure, antihypertensive drug use, chronic kidney disease, cardiovascular disease (CVD), glycosylated hemoglobin, total cholesterol, and RHR. HR, hazard ratio; CI, confidence interval.

When stratified by time-varying RHR (Supplementary Table 5), in male with time-varying RHR 70 to 79 bpm, a ∆RHR ≥5 bpm drop was significantly associated with lower risk of DM (HR, 0.569; 95% CI, 0.347 to 0.933) compared to male with a stable ∆RHR. In female with time-varying RHR <60 bpm, a ∆RHR ≥5 increase was significantly associated with higher risk of DM (HR, 1.296; 95% CI, 1.005 to 1.672) compared to those with a stable ∆RHR. A ∆RHR ≥5 drop and increase were also significantly associated with lower (HR, 0.681; 95% CI, 0.499 to 0.928) and higher (HR, 1.324; 95% CI, 1.049 to 1.670) risk of DM in female with time-varying RHR 60 to 69 bpm, respectively.

The restricted cubic spline model showed a relationship between the time-varying RHR/∆RHR and incident DM (Fig. 4). We found that a time-varying RHR of <60 bpm was associated with a lower risk of incident DM than a time-varying RHR of 60 bpm in male, not in female (Fig. 4A and B). Time-varying RHR values above 89 and 70 bpm relative to 60 bpm were associated with an increased risk of incident DM in male and female, respectively. An increase in RHR was significantly associated with an increased risk of incident DM in female but not in male (Fig. 4C and D).

Fig. 4.

Association between resting heart rate (RHR) and the incidence of diabetes mellitus (DM). The adjusted cubic spline model demonstrates the flexible association between time-varying RHR, including all available RHR values before an event or end of the study, and the hazard of DM incidence overall (A), in male (B), and in female (C) with 60 bpm taken as the reference RHR. The adjusted cubic spline model demonstrates the flexible association between time-varying change in RHR (∆RHR) from the follow-up visit and the hazard of DM incidence overall (D) when a RHR change of 0 bpm is taken as the reference. The dashed black curves represent the upper and lower 95% confidence limits. The horizontal red line represents the hazard ratio (HR) of 1.

DISCUSSION

This study demonstrates the sex-specificity of the RHR and DM risk. In this community-based prospective cohort study of middle-aged and older adults without DM at baseline, we found that RHR ≥80 bpm, as a time-varying factor, was found to be strongly related to the risk of DM. A ∆RHR ≥5 bpm drop was associated with lower DM risk. However, when stratified by sex, decrease in RHR was significantly associated with decreased risk of DM only in male, and increase in ∆RHR was significantly associated with increased risk of DM only in female.

A novel aspect of our study includes the use of the time-varying models; to our knowledge, this approach has not been used by any previous study assessing the association between RHR or ∆RHR and DM in general population. RHR and other factors such as lifestyle, physical activity, and clinical conditions may change over time. However, previous studies relied on RHR measurements single points in time and have consistently linked elevated RHR with an increased risk of DM [7-9]. Additionally, we used a traditional Cox proportional hazards model to identify the long-term effects of fixed baseline risk factors and short-term effects of a fixed last measurement (the most recent RHR) before the index date (Supplementary Fig. 3). We also found similar results in baseline model to previous studies (HR per 5 bpm increase, 1.046; 95% CI, 1.006 to 1.086) (Supplementary Table 6). However, considering that the RHR changes over time, the results of our analysis using the most recent RHR before the index date confirmed that the linear relationship between the RHR and DM was no longer significant, indicating that a U- or J-shaped relationship between the RHR and DM is possible (Supplementary Table 6). Few epidemiological studies have evaluated longitudinal changes in RHR with the risk of DM using repeated assessments of RHR [9]. However, these studies used only two time points (start and end time). Exposure and/or confounding factor changes over time generally lead to regression dilution bias, resulting in underestimation of the true association between exposure and outcome [19]. It could be addressed by performing time-varying analysis [20]. Our study population has been biennially evaluated. Thus, we assessed the associations of RHR and ∆RHR with DM during up to 14 years of follow-up by including these longitudinal repeated measurements of RHR and covariates as time-varying variables. Our findings showed that a higher RHR was associated with increased risk. Although the underlying mechanisms are yet to be fully elucidated, several biological mechanisms have been suggested to be associated with RHR and DM. RHR is a crude index of ANS status, and an elevated RHR may reflect an ANS imbalance favoring sympathetic activation [7]. Sympathetic over-activity has been linked to reduced insulin sensitivity, high blood pressure, obesity, subclinical inflammation, and metabolic syndrome [21-24], all of which are known to be associated with the risk of DM. Higher RHRs have also been linked to increased sympathetic nervous system activity which may increase both acute and chronic insulin resistance, thereby inducing DM [25,26]. In contrast, a decrease in RHR may reflect improved cardiac function, physical fitness, and a lower sympathetic tone [27].

A novel finding in our study was that the relationship between time-varying RHR and ∆RHR with incidence of DM was different according to the sex, which was not shown in any previous study. In our study, RHR was slightly higher (approximately 1 bpm) in female than in male, and the ∆RHR showed patterns in which RHR increased with increasing age only in female with incident DM (Supplementary Fig. 4). Few studies have assessed the relationship between RHR and risk of DM by sex [12]. In previous study, the risk of type 2 diabetes mellitus (T2DM) was increased with increasing RHR and decreased with decreasing RHR in both male and female [12]. While the reason for these differing results by sex is not clear, we speculate on a few possibilities. One possibility is related to adipose tissue, which is known as an immune organ releasing various immunomodulatory substances such as tumor necrosis factor-α, interleukin 6, leptin, resistin and adiponectin, serving as a notable origin of inflammatory signals associated with insulin resistance [28]. Especially, leptin is associated with increased heart rate and sympathetic activity [29,30]. High leptin concentrations are directly associated with the obesity subsequent development of metabolic disease sequelae such as insulin resistance, T2DM, and CVD [31]. Previous study showed that leptin was significantly higher in female than male [32]. In addition, female had significantly lower sympathetic tone compared to male [33], and were found to have lower resting sympathetic nerve activity to muscles than male, especially in younger than 50 years of age [34]. As expected, heart rate variability (HRV) indices vary across the menstrual cycle in healthy female, and ANS reactivity is different across phases [35,36]. Hence, sex differences in the ANS may play a role. Another potential explanation for sex-specific differences could be attributed to the characteristics of our study population. In male, the group with ∆RHR ≥5 increase showed a trend of higher DM risk compared to those with stable ∆RHR, although the statistical significance was not achieved. Notably, in male with time-varying RHR <60 bpm, there was a relatively high frequency (48.2%) of observations with an RHR increase of ≥5 bpm and the incidence of DM was comparatively low in this group (3.72% vs. 4.50% in male and female, respectively) (Supplementary Table 7). Hence, with a larger sample size, statistical significance might be achieved. Additionally, we found that male with a time-varying RHR 70 to 79 bpm had a lower risk of DM compared to reference group. However, upon subgroup analysis, the significance was observed in male without CVD and with a BMI <23 kg/m2. Thus, it is hypothesized that relatively healthy male without CVD and with a normal BMI may have additional protective effects against DM, even if their RHR falls within the range of 70 to 79 bpm, which is slightly higher than the reference group. As a result, our findings suggest that monitoring RHR over time may be useful in identifying groups at high risk for DM.

This study had some limitations. First, changes in the RHR were only assessed biennially; therefore, further studies measuring heart rate at shorter intervals should prove useful. Second, the use of heart rate-limiting drugs such as β-blockers, which are the most frequently used drugs for the treatment of hypertension, is a possible explanation for the reduction in RHR from visit to visit; thus, HRV may have more impact than ∆RHR over time. Unfortunately, we could not obtain information regarding HRV and the use of β-blockers from the KoGES data. However, in Korea, the treatment rate of β-blockers was found to be relatively low (18.5%) [37]. Additionally, subgroup analysis in subjects without the use of antihypertensive drugs and CVD was conducted, and similar patterns were found in these groups. Another factor that may play a role in the ∆RHR over time is physical activity [38]. Thus, we also controlled for physical activity using a time-varying analysis.

Nevertheless, our study has several strengths. First, we used a large homogeneous population based on a well-designed 14-year longitudinal database. Second, time-varying models that used serial measurements of risk factors as determinants of subsequent survival were used.

In a prospective longitudinal cohort study, time-varying RHR and ∆RHR were associated with the risk of DM, and higher time-varying RHR and lower time-varying ∆RHR were associated with an increased and decreased risk of DM, respectively. In addition, the relationship between RHR and ∆RHR with incidence of DM were different by gender; lower RHR was associated with decreased risk of DM in male, and higher ∆RHR was associated with increased risk of DM in female. Our findings support the importance of frequent monitoring of RHR and the identification of changes over time to identify subjects who may be at an increased and/or decreased risk of incident DM.

SUPPLEMENTARY MATERIALS

Supplementary materials related to this article can be found online at https://doi.org/10.4093/dmj.2023.0305.

Supplementary Table 1.

Example of the study

dmj-2023-0305-Supplementary-Table-1.pdf
Supplementary Table 2.

Data structures for modeling, analysis for time-varying Cox regression

dmj-2023-0305-Supplementary-Table-2.pdf
Supplementary Table 3.

Risk of diabetes mellitus for time-varying resting heart rate by univariate and multivariate time-varying Cox regression analysis in male

dmj-2023-0305-Supplementary-Table-3.pdf
Supplementary Table 4.

Risk of diabetes mellitus for time-varying resting heart rate by univariate and multivariate time-varying Cox regression analysis in female

dmj-2023-0305-Supplementary-Table-4.pdf
Supplementary Table 5.

Adjusted hazard ratio (95% confidence interval) for diabetes mellitus based on changes in the resting heart rate

dmj-2023-0305-Supplementary-Table-5.pdf
Supplementary Table 6.

Association between resting heart rate incidence of diabetes mellitus

dmj-2023-0305-Supplementary-Table-6.pdf
Supplementary Table 7.

Observation rate of diabetes mellitus incident according to the time-varying RHR and ∆RHR

dmj-2023-0305-Supplementary-Table-7.pdf
Supplementary Fig. 1.

Design for the analysis using time-varying resting heart rate (RHR) and change in RHR (∆RHR).

dmj-2023-0305-Supplementary-Fig-1.pdf
Supplementary Fig. 2.

Distribution of time-varying (A) resting heart rate (RHR) and (B) change in RHR (∆RHR) in subjects with and without incident diabetes mellitus (DM) by sex.

dmj-2023-0305-Supplementary-Fig-2.pdf
Supplementary Fig. 3.

Design for the Cox proportional hazards regression analysis using baseline and the most recent resting heart rate (RHR).

dmj-2023-0305-Supplementary-Fig-3.pdf
Supplementary Fig. 4.

Trajectories of resting heart rate (RHR) according to the age in subjects with and without incident diabetes mellitus (DM). Fitted lines denote mean RHR (incident DM vs. without incident DM), calculated from restricted cubic spline of mixed-effects linear regression models. The dashed lines delimit the 95% confidence interval. (A) Time-varying RHR and DM in male. (B) Time-varying RHR and DM in female. (C) Time-varying change in RHR (∆RHR) and DM in male. (D) Time-varying ∆RHR and DM in female.

dmj-2023-0305-Supplementary-Fig-4.pdf

Notes

CONFLICTS OF INTEREST

No potential conflict of interest relevant to this article was reported.

AUTHOR CONTRIBUTIONS

Conception or design: M.K.S., H.Y.P.

Acquisition, analysis, or interpretation of data: M.K.S.

Drafting the work or revising: K.L., H.Y.P.

Final approval of the manuscript: all authors.

FUNDING

This research was supported by the National Institute of Health research project (2021-NI-015-02).

Acknowledgements

Data used in this study were obtained from the Korean Genome and Epidemiology Study (KoGES 4851-302) of the Korea National Institute of Health, Republic of Korea. The authors wish to thank the participants and staff of the KoGES.

References

1. GBD 2017 Risk Factor Collaborators. Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018;392:1923–94.
2. World Health Organization. Global report on diabetes. Available from: https://www.who.int/publications/i/item/9789241565257 (cited 2023 Nov 22).
3. International Diabetes Federation. IDF Diabetes Atlas 9th edth ed. Brussels: IDF; 2019.
4. Whiting DR, Guariguata L, Weil C, Shaw J. IDF diabetes atlas: global estimates of the prevalence of diabetes for 2011 and 2030. Diabetes Res Clin Pract 2011;94:311–21.
5. Sabbah HN, Ilsar I, Zaretsky A, Rastogi S, Wang M, Gupta RC. Vagus nerve stimulation in experimental heart failure. Heart Fail Rev 2011;16:171–8.
6. Shcherbina A, Mattsson CM, Waggott D, Salisbury H, Christle JW, Hastie T, et al. Accuracy in wrist-worn, sensor-based measurements of heart rate and energy expenditure in a diverse cohort. J Pers Med 2017;7:3.
7. Aune D, Hartaigh BO, Vatten LJ. Resting heart rate and the risk of type 2 diabetes: a systematic review and dose: response meta-analysis of cohort studies. Nutr Metab Cardiovasc Dis 2015;25:526–34.
8. Wang L, Cui L, Wang Y, Vaidya A, Chen S, Zhang C, et al. Resting heart rate and the risk of developing impaired fasting glucose and diabetes: the Kailuan prospective study. Int J Epidemiol 2015;44:689–99.
9. Zhao Y, Zhang M, Liu Y, Yin Z, Li H, Sun H, et al. 6-Year change in resting heart rate is associated with incident type 2 diabetes mellitus. Nutr Metab Cardiovasc Dis 2019;29:236–43.
10. Nauman J, Janszky I, Vatten LJ, Wisloff U. Temporal changes in resting heart rate and deaths from ischemic heart disease. JAMA 2011;306:2579–87.
11. Eppinga RN, Hagemeijer Y, Burgess S, Hinds DA, Stefansson K, Gudbjartsson DF, et al. Identification of genomic loci associated with resting heart rate and shared genetic predictors with all-cause mortality. Nat Genet 2016;48:1557–63.
12. Liu D, Qin P, Liu Y, Sun X, Li H, Wu X, et al. Sex-specific association of resting heart rate with type 2 diabetes mellitus. J Diabetes Complications 2020;34:107754.
13. Kim Y, Han BG, ; KoGES group. Cohort profile: the Korean Genome and Epidemiology Study (KoGES) consortium. Int J Epidemiol 2017;46e20.
14. Allison PD. Survival analysis using SAS: a practical guide 2nd edth ed. Cary: SAS Institute; 2010.
15. Fisher LD, Lin DY. Time-dependent covariates in the Cox proportional-hazards regression model. Annu Rev Public Health 1999;20:145–57.
16. Vazir A, Claggett B, Pitt B, Anand I, Sweitzer N, Fang J, et al. Prognostic importance of temporal changes in resting heart rate in heart failure and preserved ejection fraction: from the TOPCAT Study. JACC Heart Fail 2017;5:782–91.
17. Cui X, Mandalenakis Z, Thunstrom E, Fu M, Svardsudd K, Hansson PO. The impact of time-updated resting heart rate on cause-specific mortality in a random middle-aged male population: a lifetime follow-up. Clin Res Cardiol 2021;110:822–30.
18. Archangelidi O, Pujades-Rodriguez M, Timmis A, Jouven X, Denaxas S, Hemingway H. Clinically recorded heart rate and incidence of 12 coronary, cardiac, cerebrovascular and peripheral arterial diseases in 233,970 men and women: a linked electronic health record study. Eur J Prev Cardiol 2018;25:1485–95.
19. Hutcheon JA, Chiolero A, Hanley JA. Random measurement error and regression dilution bias. BMJ 2010;340:c2289.
20. Smabrekke B, Rinde LB, Hindberg K, Hald EM, Vik A, Wilsgaard T, et al. Atherosclerotic risk factors and risk of myocardial infarction and venous thromboembolism; time-fixed versus time-varying analyses: the Tromso Study. PLoS One 2016;11e0163242.
21. Sajadieh A, Nielsen OW, Rasmussen V, Hein HO, Abedini S, Hansen JF. Increased heart rate and reduced heart-rate variability are associated with subclinical inflammation in middle-aged and elderly subjects with no apparent heart disease. Eur Heart J 2004;25:363–70.
22. Shibao C, Gamboa A, Diedrich A, Ertl AC, Chen KY, Byrne DW, et al. Autonomic contribution to blood pressure and metabolism in obesity. Hypertension 2007;49:27–33.
23. Flanagan DE, Vaile JC, Petley GW, Moore VM, Godsland IF, Cockington RA, et al. The autonomic control of heart rate and insulin resistance in young adults. J Clin Endocrinol Metab 1999;84:1263–7.
24. Mancia G, Bousquet P, Elghozi JL, Esler M, Grassi G, Julius S, et al. The sympathetic nervous system and the metabolic syndrome. J Hypertens 2007;25:909–20.
25. Jamerson KA, Julius S, Gudbrandsson T, Andersson O, Brant DO. Reflex sympathetic activation induces acute insulin resistance in the human forearm. Hypertension 1993;21:618–23.
26. Kim DI, Yang HI, Park JH, Lee MK, Kang DW, Chae JS, et al. The association between resting heart rate and type 2 diabetes and hypertension in Korean adults. Heart 2016;102:1757–62.
27. Vazir A, Claggett B, Cheng S, Skali H, Shah A, Agulair D, et al. Association of resting heart rate and temporal changes in heart rate with outcomes in participants of the atherosclerosis risk in communities study. JAMA Cardiol 2018;3:200–6.
28. Wisse BE. The inflammatory syndrome: the role of adipose tissue cytokines in metabolic disorders linked to obesity. J Am Soc Nephrol 2004;15:2792–800.
29. Shek EW, Brands MW, Hall JE. Chronic leptin infusion increases arterial pressure. Hypertension 1998;31(1 Pt 2):409–14.
30. Haynes WG, Morgan DA, Walsh SA, Mark AL, Sivitz WI. Receptor-mediated regional sympathetic nerve activation by leptin. J Clin Invest 1997;100:270–8.
31. Ghadge AA, Khaire AA. Leptin as a predictive marker for metabolic syndrome. Cytokine 2019;121:154735.
32. Shamsuzzaman AS, Winnicki M, Wolk R, Svatikova A, Phillips BG, Davison DE, et al. Independent association between plasma leptin and C-reactive protein in healthy humans. Circulation 2004;109:2181–5.
33. Cowan MJ, Pike K, Burr RL. Effects of gender and age on heart rate variability in healthy individuals and in persons after sudden cardiac arrest. J Electrocardiol 1994;27 Suppl:1–9.
34. Dart AM, Du XJ, Kingwell BA. Gender, sex hormones and autonomic nervous control of the cardiovascular system. Cardiovasc Res 2002;53:678–87.
35. Sato N, Miyake S, Akatsu J, Kumashiro M. Power spectral analysis of heart rate variability in healthy young women during the normal menstrual cycle. Psychosom Med 1995;57:331–5.
36. Saeki Y, Atogami F, Takahashi K, Yoshizawa T. Reflex control of autonomic function induced by posture change during the menstrual cycle. J Auton Nerv Syst 1997;66:69–74.
37. Kim SH, Shin DW, Kim S, Han K, Park SH, Kim YH, et al. Prescribing patterns of antihypertensives for treatment-naive patients in South Korea: from Korean NHISS claim data. Int J Hypertens 2019;2019:4735876.
38. Black A, Murray L, Cardwell C, Smith GD, McCarron P. Secular trends in heart rate in young adults, 1949 to 2004: analyses of cross sectional studies. Heart 2006;92:468–73.

Article information Continued

Fig. 1.

Trajectories of resting heart rate (RHR) before incident diabetes mellitus (DM) or last examination. Fitted lines denote mean RHR (incident DM vs. without incident DM), calculated from restricted cubic spline of mixed-effects linear regression models. Time 0=diagnosis or last examination. The dashed lines delimit the 95% confidence interval. (A) Time-varying RHR and DM in male. (B) Time-varying RHR and DM in female. (C) Time-varying change in RHR (∆RHR) and DM in male. (D) Time-varying ∆RHR and DM in female.

Fig. 2.

Subgroup analysis for the risk of diabetes mellitus according to the resting heart rate (RHR). Adjusted for sex, area, and time-updated age, body mass index (BMI), physical activity, smoking, drinking, systolic blood pressure, antihypertensive drug use, chronic kidney disease, cardiovascular disease (CVD), glycosylated hemoglobin, total cholesterol, and change in RHR. HR, hazard ratio; CI, confidence interval.

Fig. 3.

Subgroup analysis for the risk of diabetes mellitus according to the change in resting heart rate (∆RHR). Adjusted for sex, area, and time-updated age, body mass index (BMI), physical activity, smoking, drinking, systolic blood pressure, antihypertensive drug use, chronic kidney disease, cardiovascular disease (CVD), glycosylated hemoglobin, total cholesterol, and RHR. HR, hazard ratio; CI, confidence interval.

Fig. 4.

Association between resting heart rate (RHR) and the incidence of diabetes mellitus (DM). The adjusted cubic spline model demonstrates the flexible association between time-varying RHR, including all available RHR values before an event or end of the study, and the hazard of DM incidence overall (A), in male (B), and in female (C) with 60 bpm taken as the reference RHR. The adjusted cubic spline model demonstrates the flexible association between time-varying change in RHR (∆RHR) from the follow-up visit and the hazard of DM incidence overall (D) when a RHR change of 0 bpm is taken as the reference. The dashed black curves represent the upper and lower 95% confidence limits. The horizontal red line represents the hazard ratio (HR) of 1.

Table 1.

Characteristics of participants at visits 1–8

Characteristic Visit 1 (2001–2002) Visit 2 (2003–2004) Visit 3 (2005–2006) Visit 4 (2007–2008) Visit 5 (2009–2010) Visit 6 (2011–2012) Visit 7 (2013–2014) Visit 8 (2015–2016)
No. of participants 8,313 7,118 6,130 5,394 5,198 4,754 4,407 4,389
Demographics
 Age, yr 51.9±8.80 53.7±8.8 55.5±8.6 57.3±8.5 59.1±8.5 60.8±8.4 62.3±8.2 64.4±8.2
 Female sex 4,397 (52.9) 3,754 (52.7) 3,257 (53.1) 2,883 (53.4) 2,819 (54.2) 2,582 (54.3) 2,427 (55.1) 2,393 (54.5)
 Residence, urban 4,026 (48.4) 3,373 (47.4) 2,937 (47.9) 2,673 (49.6) 2,589 (49.8) 2,385 (50.2) 2,268 (51.5) 2,233 (50.9)
Lifestyle variables
 Current smoker 1,835 (22.4) 1,523 (21.7) 1,262 (20.9) 1,093 (20.5) 1,022 (29.9) 902 (19.2) 838 (19.3) 856 (19.7)
 Current drinker 3,944 (47.9) 3,375 (47.8) 2,917 (48.0) 2,554 (47.7) 2,471 (47.9) 2,266 (48.1) 2,097 (48.0) 2,066 (47.5)
 Physical activity, MET-hr/day 30.2±17.6 30.8±17.6 38.3±15.4 36.1±16.1 36.3±15.8 36.1±15.8 36.6±15.7 36.0±16.0
Clinical characteristics
 BMI, kg/m2 24.5±3.1 24.5±3.1 24.3±3.0 24.3±3.0 24.3±3.0 24.3±3.1 24.3±3.1 24.3±3.1
 Hypertension 2,433 (29.3) 2,032 (28.6) 1,808 (29.5) 1,737 (32.2) 1,926 (37.1) 1,898 (39.9) 1,854 (42.1) 1,824 (41.6)
 CVD 227 (2.7) 135 (1.9) 129 (2.1) 161 (3.0) 172 (3.3) 192 (4.0) 207 (4.7) 221 (5.0)
 CKD 565 (6.8) 528 (7.4) 356 (5.8) 248 (4.6) 216 (4.2) 574 (12.1) 588 (13.3) 520 (12.0)
 SBP, mm Hg 120.9±18.1 117.8±17.0 115.9±16.4 117.4±16.1 119.1±16.1 119.3±15.9 120.5±16.7 120.4±16.0
 DBP, mm Hg 80.2±11.4 78.4±10.9 77.8±10.3 77.2±9.7 78.2±9.6 77.4±9.6 76.9±9.9 75.8±9.4
 RHR, bpm 63.5±7.3 62.9±6.8 64.1±6.8 62.8±6.9 63.1±7.4 62.8±7.4 61.7±7.2 62.8±7.4
 Total cholesterol, mg/dL 190.2±34.6 197.4±35.1 191.5±34.5 196.0±34.1 197.0±34.8 193.7±34.1 191.0±33.7 189.2±34.4
 HDL-C, mg/dL 44.9±10.1 46.6±9.9 44.7±10.3 44.8±10.7 44.1±10.8 48.1±12.1 46.8±11.7 47.1±12.1
 Triglyceride, mg/dL 157.2±97.7 136.8±91.6 134.1±93.6 134.6±84.3 138.4±92.2 136.2±86.0 131.8±86.5 130.7±81.0
 Fasting glucose, mg/dL 83.6±10.1 91.2±10.6 90.2±9.2 92.1±9.5 94.5±9.9 92.4±9.3 90.8±9.0 92.4±9.4
 HbA1c, % 5.59±0.4 5.58±0.4 5.42±0.4 5.49±0.4 5.58±0.4 5.53±0.4 5.53±0.3 5.60±0.3
 Fasting insulin, µIU/mL 7.60±4.8 7.69±3.9 7.44±3.4 8.70±5.3 8.65±4.5 7.96±3.6 8.41±3.8 8.94±5.1

Values are presented as mean±standard deviation or number (%).

MET, metabolic equivalent of task; BMI, body mass index; CVD, cardiovascular disease; CKD, chronic kidney disease; SBP, systolic blood pressure; DBP, diastolic blood pressure; RHR, resting heart rate; HDL-C, high-density lipoprotein cholesterol; HbA1c, glycosylated hemoglobin.

Table 2.

Risk of diabetes mellitus by time-varying RHR and ΔRHR, time-varying Cox regression analysis

Variable Male
Female
No. of observations No. of cases (%) HR (95% CI) P value No. of observations No. of cases (%) HR (95% CI) P value
Time-varying RHR, bpm
 <60 5,965 144 (2.41) 0.816 (0.667–0.997) 0.047 5,272 131 (2.48) 0.819 (0.665–1.010) 0.062
 60–69 11,993 630 (3.00) Reference 15,343 442 (2.88) Reference
 70–79 2,613 85 (3.25) 0.700 (0.523–0.937) 0.017 3,217 101 (3.14) 1.126 (0.892–1.422) 0.320
 ≥80 620 35 (5.65) 2.195 (1.514–3.183) <0.001 680 47 (6.91) 2.224 (1.562–3.165) <0.001
∆RHR, bpm
 ≥5 drop 4,552 123 (2.70) 0.602 (0.471–0.770) <0.001 5,282 153 (2.90) 0.806 (0.651–0.998) 0.048
 <5 change 11,884 353 (2.97) Reference 14,175 374 (2.64) Reference
 ≥5 increase 4,755 148 (3.11) 1.057 (0.869–1.285) 0.582 5,055 194 (3.84) 1.218 (1.008–1.471) 0.041

Time-varying RHR was adjusted for sex, area, and time-updated age, body mass index, physical activity, smoking, drinking, systolic blood pressure, antihypertensive drug use, chronic kidney disease, cardiovascular disease, glycosylated hemoglobin (HbA1c), total cholesterol, and change in RHR (ΔRHR). For ΔRHR, adjusted for sex, area and time-updated age, body mass index, physical activity, smoking, drinking, systolic blood pressure, antihypertensive drug use, chronic kidney disease, cardiovascular disease, HbA1c, total cholesterol, and RHR.

RHR, resting heart rate; ΔRHR, change in resting heart rate; HR, hazard ratio; CI, confidence interval.