Association between heart rate recovery after exercise and renal function in patients referred for treadmill exercise test

Abstract

Introduction

Heart rate recovery (HRR) is a marker of parasympathetic activity recovery after exercise, and it is associated with cardiovascular mortality and total mortality. Impaired renal function is also associated with cardiac mortality. The aim of this study was to investigate the association between HRR after exercise and renal function in patients referred for a treadmill exercise test.

Patients and methods

This cross-sectional study was conducted at a regional hospital in southern Taiwan. Patients who completed a symptom-limited treadmill exercise test from January 2015 to February 2018 were recruited. Before the treadmill exercise test, patients were asked to complete a questionnaire on the past disease history and lifestyle factors. Serum creatinine measurement within two years prior to or after the date of the treadmill exercise test of the patients was also obtained from the medical records for these patients. Estimated glomerular filtration rate (eGFR) was calculated. Simple and multiple linear regression analyses were performed to investigate the association between one-minute HRR and eGFR.

Results

A total of 2,825 patients completed the treadmill exercise test, and serum creatinine measurement was identified from medical records for 2,153 patients (76.2%). Multiple linear regression analysis revealed that a lower eGFR was significantly associated with lower one-minute HRR (P< 0.001), adjusting for other significant independent factors, including age, waist circumference, type 2 diabetes mellitus, and smoking.

Conclusions

In this cross-sectional observational study, a lower eGFR was significantly and independently associated with decreased one-minute HRR, suggesting that parasympathetic activity recovery after exercise could be impaired by a decrease in renal function.

About the Authors:

Rei-Yeuh Chang

Roles Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Writing – original draft

Affiliations Division of Cardiology, Department of Internal Medicine, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi City, Taiwan, Chung Jen Junior College of Nursing, Health Sciences and Management, Chiayi, Taiwan, Min-Hwei Junior College of Health Care Management, Tainan City, Taiwan

Han-Lin Tsai

Roles Data curation, Investigation, Writing – review & editing

Affiliation: Division of Cardiology, Department of Internal Medicine, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi City, Taiwan

Ping-Gune Hsiao

Roles Data curation, Investigation, Writing – review & editing

Affiliation: Division of Cardiology, Department of Internal Medicine, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi City, Taiwan

Chao-Wen Tan

Roles Data curation, Investigation, Writing – review & editing

Affiliation: Division of Cardiology, Department of Internal Medicine, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi City, Taiwan

Chi-Pin Lee

Roles Data curation, Investigation, Writing – review & editing

Affiliation: Division of Cardiology, Department of Internal Medicine, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi City, Taiwan

I-Tseng Chu

Roles Data curation, Investigation, Writing – review & editing

Affiliation: Division of Cardiology, Department of Internal Medicine, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi City, Taiwan

Yung-Ping Chen

Roles Data curation, Investigation, Writing – review & editing

Affiliation: Division of Cardiology, Department of Internal Medicine, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi City, Taiwan

Malcolm Koo

Roles Formal analysis, Writing – review & editing

* E-mail: [email protected]

Affiliations Graduate Institute of Long-term Care, Tzu Chi University of Science and Technology, Hualien City, Hualien, Taiwan, Dalla Lana School of Public Health, University of Toronto, Toronto, Ontario, Canada

ORCID logo http://orcid.org/0000-0002-9242-9167

Introduction

The treadmill exercise test is a widely used tool to diagnose obstructive coronary artery disease. Hemodynamic data during treadmill exercise test can also provide valuable diagnostic information, including exercise capacity, heart rate reserve, blood pressure reserve, rate-pressure product, chronotropic reserve, heart rate recovery, rate-pressure product reserve, and arrhythmia. Heart rate recovery (HRR) is a decrease in heart rate immediately after termination of exercise, which is related to the restoration of the parasympathetic nervous system after exercise [1]. Older age, low physical activity [2], metabolic syndrome [3], type 2 diabetes mellitus [4, 5], and smoking [6] had been reported to be associated with a low HRR.

Previous studies indicated that a low HRR was independently associated with cardiovascular mortality and all-cause mortality [7, 8]. A low HRR was also related to an increased risk of sudden death from myocardial infarction, but not non-sudden death from myocardial infarction [9]. Vagal tone restoration post exercise appeared to play a critical role in the prevention of potentially lethal arrhythmias and sudden death. An attenuated HRR was also associated with advanced coronary artery calcification independent of traditional coronary artery disease risk factors and the development of new-onset atrial fibrillation [10, 11].

Chronic kidney disease (CKD) is condition characterized by a gradual and progressive loss of kidney function over time, and is defined as a reduced glomerular filtration rate (GFR), increased urinary albumin excretion, or both, present for over 3 months. Globally, the burden of CKD continues to increase and the prevalence is estimated to be 8 to 16% [12]. CKD was first defined based on the National Kidney Foundation (NKF) Kidney Disease Outcomes Quality Initiative (KDOQI) introduced in 2002 [13] and endorsed at subsequent Kidney Disease Improving Global Outcomes (KDIGO) Controversies Conferences with minor modifications in 2004 and further in 2018 [14, 15]. The classification is widely used in clinical practice, research, prevention, and treatment of CKD. A large community-based study showed that a low estimated glomerular filtration rate (eGFR) was associated with increased risk for cardiovascular events, hospitalization, and mortality [16–18]. In addition, a recent meta-analysis of 83 studies reporting 30,392 strokes revealed that the risk of stroke increased linearly and additively with declining eGFR [19].

Although both a low HRR and eGFR have been reported to be associated with an increased cardiovascular risk, few studies have investigated the association between HRR and eGFR. A study on 107 patients with CKD had directly investigated the association between HRR and eGFR. Results from multiple linear regression analysis indicated that reduced eGFR was a significant and independent risk factor for a lower HRR [20]. Therefore, the aim of this cross-sectional study with a large sample size was to evaluate the association between eGFR and HRR in Taiwanese patients referred for a treadmill exercise test.

Patients and methods

Study design and study population

A cross-sectional observational study design was used. The study participants included all patients referred for treadmill exercise test at a regional hospital in southern Taiwan between January 2015 and February 2018. Patients were excluded if they had an implanted pacemaker or atrial fibrillation. The study protocol was approved by the Institutional Review Board of Chia-Yi Christian Hospital, Taiwan (CYCH-IRB No. 104008). All patients provided written informed consent.

Treadmill exercise testing

After a supine and standing electrocardiography (ECG), heart rate, and blood pressure were obtained, a symptom-limited treadmill exercise testing using the standard Bruce protocol was performed (GE T-2000 treadmill, CardioSoft diagnostic software version 5.20, Marquette ECG analysis program, GE Medical Systems IT, Inc. Milwaukee, USA). Blood pressure was measured non-invasively using a Suntech 4240 monitor (Suntech Medical Instruments, Raleigh, NC) from left brachial artery. Patients were encouraged to achieve a goal of 85% of maximum age-predicted heart rate (HRmax) (beats/min), which was calculated by 220 minus the age of the participants. Percentage of maximum age predicted heart rate (%HRmax) was calculated by (peak heart rate/maximum age predicted heart rate) × 100. Functional capacity was estimated based on a range of speeds and grades of the treadmill. It was expressed in metabolic equivalent of task (MET), where one MET is equivalent to 3.5 mL/kg/min of oxygen consumption. ECG, heart rate, and blood pressure were recorded at each stage of exercise, peak exercise, and every minute in the recovery phase up to 6 minutes. A positive ischemic ST-segment response was defined as the horizontal or downsloping ST-segment depression of > 1 mm below baseline taken 80 ms after the J-point. Rate-pressure product (RPP) (mmHg × bpm × 10−2) was calculated as the product of systolic blood pressure (mmHg) and heart rate (beats/min) at peak exercise and then divided by 100. The dependent variable in this study was one-minute HRR, which was defined as the decrease in heart rate from peak exercise to 1 minute after exercise. All treadmill exercise testing data were obtained from chart review.

Clinical data and questionnaire

A medical chart review was conducted to obtain measurements on serum creatinine of the patients who had completed the treadmill exercise test. Only patients with serum creatinine measurement within two years prior to or after the date of the treadmill exercise test were included in the study. eGFR was calculated with the Modification of Diet in Renal Disease (MDRD) equation [eGFR (mL/min/1.73m2) = 186 × serum creatinine -1.154 × age-0.203 (× 0.742 if female)] [21]. In addition, the values of eGFR≥ 90 mL/min/1.73 m2 are defined as CKD stage 1 or normal. The values of eGFR between 60‒89 mL/min/1.73 m2 are defined as CKD stage 2. The values of eGFR between 45‒59 mL/min/1.73 m2 are defined as CKD stage 3a. The values of eGFR between 30‒44 mL/min/1.73 m2 are defined as CKD stage 3b. The values of eGFR between 15‒29 mL/min/1.73 m2 are defined as CKD stage 4. The values of eGFR < 15 mL/min/1.73 m2 are defined as CKD stage 5.

Body weight, body height, and waist circumference of the patients were measured prior to the treadmill exercise test. Moreover, the patients were asked to complete a questionnaire on past history of diseases (including hypertension, diabetes mellitus, coronary artery disease, hyperlipidemia, chronic obstructive pulmonary disease, and renal disease), lifestyle factors (smoking, drinking, exercise, and perceived health status), and sleep quality. The Pittsburgh Sleep Quality Index (PSQI) questionnaire was used to assess sleep quality over the previous one-month period. A PSQI global score was calculated with a range from 0 to 21 points. A score of > 5 is indicative of poor sleep quality [22].

Statistical analysis

Continuous variables are expressed as mean and standard deviation (SD). Categorical data are expressed as number and percentage. Simple linear regression and stepwise multiple linear regression analyses were conducted to investigate the association between HRR and renal function. Two multiple linear regression models were evaluated with renal function expressed as either stages of CKD (Model 1) or eGFR (Model 2). A stepwise regression method was used to obtain the final model, adjusting for potential confounding variables. All statistical analyses were performed with PASW Statistics for Windows, version 18.0 (IBM Corp., Chicago, USA). A two-tailed P value of < 0.05 was considered statistically significant.

Results

Of the 2,825 patients completed the treadmill exercise test, serum creatinine measurement could be identified for 2,153 patients (76.2%) from the medical records. The mean duration between the treadmill exercise test and creatinine measurement was 20.4 days (SD 159.8 days) with 40% of the duration were within 30 days. The basic demographic and clinical data of the patients are shown in Table 1. The mean age was 53.4 years and 55.6% were male patients. The mean body mass index was 25.6 kg/m2 and the mean waist circumference was 86.3 cm. The mean eGFR was 106.6 mL/min/1.73m2. Moreover, 72.7% of the patients were CKD stage 1, 23.2% were stage 2, 2.6% were stage 3a, 0.7% were stage 3b, 0.3% were stage 4, and 0.6% were stage 5.

[Figure omitted. See PDF.]

Table 1. Demographic and clinical characteristics of study participants (N = 2153).

https://doi.org/10.1371/journal.pone.0222236.t001

Treadmill exercise test results showed that the mean one-minute HRR was 22.1 beats/min (SD 9.5), RPP reserve was 16.2 mmHg × bpm × 10−2 (SD 4.6), functional capacity was 9.5 METs (SD 2.9), ventricular premature contractions (VPCs) was 11.9%, atrial premature contractions (APCs) was 5.2%, and ischemic ST-T change during exercise was 19.7% (Table 2).

[Figure omitted. See PDF.]

Table 2. Results of treadmill maximal exercise test.

https://doi.org/10.1371/journal.pone.0222236.t002

Table 3 shows the results of simple linear regression analysis, and two models of multiple linear regression analyses, using either stages of CKD (Model 1) or eGFR (Model 2) as a measure of renal function. Simple linear regression analysis revealed that age, sex, body mass index, waist circumference, hypertension, type 2 diabetes mellitus, coronary artery disease, hyperlipidemia, smoking, poor or very poor perceived health status, eGFR, and stages of CKD were significantly associated with one-minute HRR. In addition, multiple linear regression analyses showed that renal function expressed by either CKD stage (stage 2, stage 3a, and stage 3b, 4, 5) or eGFR was significantly and independently associated with one-minute HRR, adjusting for age, waist circumference, type 2 diabetes mellitus, and smoking. The regression coefficients associated with stage 2, stage 3a, and stage 3b, 4, 5 were −1.94, −3.96, and −5.20, respectively, with a dose-response pattern (P for trend < 0.001). The regression coefficient for eGFR was 0.036 (P < 0.001).

[Figure omitted. See PDF.]

Table 3. Simple and multiple linear regression analyses of one-minute heart rate recovery.

https://doi.org/10.1371/journal.pone.0222236.t003

Discussion

This cross-sectional study revealed that renal function, expressed as either CKD stage (stage 2, stage 3a, and stage 3b, 4, 5) or eGFR was significantly and independently associated with one-minute HRR in adults. Késoi et al. reported a significant and independent association between eGFR and HRR in Hungarian patients with biopsy-confirmed chronic IgA nephropathy [20]. A slower one-minute HRR after exercise stress test was also observed in patients with primary nephrotic syndrome [23]. Our findings generalize these results to the Taiwanese population with a larger sample size of 2,153 patients. We also showed that a significant association between impaired renal function and HRR could be observed regardless of whether renal function was expressed either as a continuous variable eGFR or when it was categorized as stages of CKD.

The correlation between renal function and HRR could be related to autonomic dysfunction in patients with CKD [24]. Several possible mechanisms have been proposed to the dysfunction, including impaired reflex control of autonomic activity, activation of the renin-angiotensin-aldosterone system, activation of renal afferents, structural remodeling of the heart and vasculature, decreased nitric oxide bioavailability, and increased mental stress. Nevertheless, the precise mechanism is still unclear [25].

In addition to CKD stages, an older age, a larger waist circumference, type 2 diabetes mellitus, and smoking were also found to be independently associated with a low one-minute HRR in the present study. Sympathetic responses and parasympathetic responses are known to decline with increasing age [26]. The Coronary Artery Risk Development in Young Adults (CARDIA) cohort study showed that HRR declined 2.5 beats/min over 7 years. The study also found that participants with increased physical activity during the study period had the least declined HRR (−1.3 beats/min), whereas those with decreased physical activity had the most declined HRR (−3.6 beats/min) [2]. Many studies aimed to modify HRR by exercise training or rehabilitation program. For example, an exercise training study revealed that high intensity interval training could improve HRR, but there was no improvement after a high volume of low intensity training [27]. Another cardiac rehabilitation study revealed that a cardiac rehabilitation exercise program had a positive effect on HRR, while home-based exercise group was not able to improve HRR [28]. In addition, a study on patients with stable chronic heart failure reported that continuous but not interval exercise training could improve HRR [29]. In the present study, exercise ≥ 3 days/week with a duration ≥ 30 min was not significantly associated with HRR. Since exercise intensity was not well defined by our questionnaire, further study on the association between HRR and exercise should include measurement of intensity in addition to duration and frequency. Further studies using a longitudinal design should be conducted in patients with CKD to clarify whether the restoration of parasympathetic tone by increasing HRR could decrease the progression of CKD.

Waist circumference was found to be significantly associated with one-minute HRR in the present study. Waist circumference, but not body mass index has previously been shown to inversely associate with HRR in obese individuals [3]. Visceral abdominal fat was found to be associated with autonomic function on HRV assessment [30], which might be related to catecholamine production and adipocytokines secretion from fat cells [31, 32]. A study from the National Health and Nutrition Examination Survey in Taiwan also revealed that waist circumference, systolic blood pressure, serum glucose, serum C-reactive protein were inversely associated with HRR in girls, but only waist circumference associated was associated with HRR in boys [33].

In our study, type 2 diabetes was inversely associated with one-minute HRR, indicating that patients with diabetes had declined parasympathetic recovery post exercise. A cohort study of 1,818 patients underwent a routine coronary artery disease screening program also found an increased risk of abnormal HRR in patients with diabetes, even after adjusting for a number of potential confounding factors [4]. A recent study on 123 patients also reported that post-exercise cardiac autonomic recovery as assessed by HRR and heart rate variability was found to be impaired in patients with type 2 diabetes, and the effects were more pronounced in patients with poor glycemic control [34]. Another study aimed to investigate the risk factors of HRR in patients with type 2 diabetes reported that fasting blood glucose, glycosylated hemoglobin, low-density lipoprotein cholesterol, and resting and maximum heart rates were significantly associated with HRR. The study also found that a combined aerobic and resistance training program of moderate intensity was able to improve HRR, possibly due to better glycemic control, resting heart rate, and physical fitness [35].

Smoking was found to be associated with a lower one-minute HRR in the present study. Smoking is one of the major cardiovascular risk factors and its pathophysiology includes platelet aggregation, endothelial dysfunction, and coronary vasoconstriction [36]. Smoking can also induce autonomic dysfunction. Heavy smoking could impair heart rate variability, heart rate turbulence [37], and HRR [38]. The CARDIA cohort study with 2,730 participants followed for a period of 20 years revealed that current smoking, in addition, to a higher body mass index, larger waist, low education, and fasting glucose level were significant factors independently associated with incident slow HRR [6]. Since smoking is a modifiable risk factor for low HRR and autonomic dysfunction, our study highlights the importance of cessation.

Our findings should be interpreted in light of some limitations. First, previous research indicated that exercise intensity is a key factor in improving HRR. However, exercise, dichotomized based on a cut-off of ≥ 3 days/week with a duration of ≥ 30 min, was not found to be significantly associated with HRR in our study. Future studies should include questions that specifically address the intensity of exercise. Second, our patients were recruited from a single regional hospital and those referred for treadmill exercise test, which might limit the generalizability of our results.

Conclusions

In this cross-sectional observational study, renal function expressed as either eGFR or CKD stages were significantly and independently associated with one-minute HRR. Our findings suggested that impaired parasympathetic activity recovery after exercise is associated with a decline in renal function.

Supporting information

[Figure omitted. See PDF.]

S1 Dataset.

https://doi.org/10.1371/journal.pone.0222236.s001

(SAV)

Citation: Chang R-Y, Tsai H-L, Hsiao P-G, Tan C-W, Lee C-P, Chu I-T, et al. (2019) Association between heart rate recovery after exercise and renal function in patients referred for treadmill exercise test. PLoS ONE 14(9): e0222236. https://doi.org/10.1371/journal.pone.0222236

References

1. Pierpont GL, Stolpman DR, Gornick CC. Heart rate recovery post-exercise as an index of parasympathetic activity. J Auton Nerv Syst. 2000;80(3):169–174. pmid:10785283.

2. Carnethon MR, Jacobs JD, Sidney S, Sternfeld B, Gidding SS, Shoushtari C, et al. A longitudinal study of physical activity and heart rate recovery: CARDIA, 1987–1993. Med Sci Sports Exerc. 2005;37(4):606–612. pmid:15809559.

3. Fidan-Yaylali G, Yaylali YT, Erdogan C, Can B, Senol H, Gedik-Topcu B, et al. The association between central adiposity and autonomic dysfunction in obesity. Med Princ Pract. 2016;25(5):442–448. pmid:27194294.

4. Seshadri N, Acharya N, Lauer MS. Association of diabetes mellitus with abnormal heart rate recovery in patients without known coronary artery disease. Am J Cardiol. 2003;91(1):108–111. pmid:12505588.

5. Sydó N, Sydó T, Merkely B, Carta KG, Murphy JG, Lopez-Jimenez F, et al. Impaired heart rate response to exercise in diabetes and its long-term significance. Mayo Clin Proc. 2016; 91(2):157–65. pmid:26769183.

6. Carnethon MR, Sternfeld B, Liu K, Jacobs DR Jr., Schreiner PJ, Williams OD, et al. Correlates of heart rate recovery over 20 years in a healthy population sample. Med Sci Sports Exerc. 2012;44(2):273–279. pmid:21796053.

7. Cole CR, Foody JM, Blackstone EH, Lauer MS. Heart rate recovery after submaximal exercise testing as a predictor of mortality in a cardiovascularly healthy cohort. Ann Intern Med. 2000;132(7):552–555. pmid:10744592.

8. Qiu S, Cai X, Sun Z, Li L, Zuegel M, Steinacker JM, et al. Heart rate recovery and risk of cardiovascular events and all‐cause mortality: a meta‐analysis of prospective cohort studies. J Am Heart Assoc. 2017;6(5):e005505. pmid:28487388.

9. Jouven X, Empana J-P, Schwartz PJ, Desnos M, Courbon D, Ducimetière P. Heart-rate profile during exercise as a predictor of sudden death. N Engl J Med.2005;352(19):1951–1958. pmid:15888695.

10. Sabbag A, Berkovitch A, Sidi Y, Kivity S, Ben Zekry S, Beinart R, et al. Poor heart rate recovery is associated with the development of new-onset atrial fibrillation in middle-aged adults. Mayo Clin Proc. 2016;91(12):1769–1777. pmid:28126153.

11. Jae SY, Kurl S, Laukkanen JA, Yoon ES, Choi YH, Fernhall B, et al. Relation of heart rate recovery after exercise testing to coronary artery calcification. Ann Med. 2017;49(5):404–410. pmid:28166422.

12. Jha V, Garcia-Garcia G, Iseki K, Li Z, Naicker S, Plattner B, et al. Chronic kidney disease: global dimension and perspectives. Lancet. 2013;382(9888):260–272. pmid:23727169.

13. Levey AS, Coresh J, Bolton K, Culleton B, Harvey KS, Ikizler TA, et al. K/DOQI clinical practice guidelines for chronic kidney disease: evaluation, classification, and stratification. Am J Kidney Dis. 2002;39(2 SUPPL. 1). pmid:11904577.

14. Kidney Disease: Improving Global Outcomes (KDIGO) Hepatitis C Work Group. KDIGO 2018 clinical practice guideline for the prevention, diagnosis, evaluation, and treatment of hepatitis c in chronic kidney disease. Kidney Int Suppl. 2018;8(3):91–165. pmid:30675443.

15. Levey AS, Eckardt K-U, Tsukamoto Y, Levin A, Coresh J, Rossert J, et al. Definition and classification of chronic kidney disease: a position statement from Kidney Disease: Improving Global Outcomes (KDIGO). Kidney Int. 2005;67(6):2089–2100. pmid:15882252.

16. Go AS, Chertow GM, Fan D, McCulloch CE, Hsu CY. Chronic kidney disease and the risks of death, cardiovascular events, and hospitalization. N Engl J Med. 2004;351(13):1296–1305. pmid:15385656.

17. Van Der Velde M, Matsushita K, Coresh J, Astor BC, Woodward M, Levey AS, et al. Lower estimated glomerular filtration rate and higher albuminuria are associated with all-cause and cardiovascular mortality. A collaborative meta-analysis of high-risk population cohorts. Kidney Int. 2011;79(12):1341–1352. pmid:21307840.

18. Chronic Kidney Disease Prognosis Consortium, Matsushita K, van der Velde M, Astor BC, Woodward M, et al. Association of estimated glomerular filtration rate and albuminuria with all-cause and cardiovascular mortality in general population cohorts: a collaborative meta-analysis. Lancet. 2010;375(9731):2073–2081. pmid:20483451.

19. Masson P, Webster AC, Hong M, Turner R, Lindley RI, Craig JC. Chronic kidney disease and the risk of stroke: a systematic review and meta-analysis. Nephrol Dial Transplant. 2015;30(7):1162–1169. pmid:25681099.

20. Késői I, Sági B, Vas T, Kovács T, Wittmann I, Nagy J. Heart rate recovery after exercise is associated with renal function in patients with a homogenous chronic renal disease. Nephrol Dial Transplant. 2009;25(2):509–513. pmid:19783602.

21. Crowe E, Halpin D, Stevens P, Guideline Development Group. Early identification and management of chronic kidney disease: summary of NICE guidance. BMJ. 2008;337:a1530. pmid:18824486.

22. Zhu B, Xie M, Park CG, Kapella MC. Adaptation of the Pittsburgh Sleep Quality Index in Chinese adults with type 2 diabetes. J Chin Med Assoc. 2018;81(3):242–247. pmid:29258729.

23. Ozkayar N, Akyel S, Dede F, Akyel F, Turgut D, Bulut M, et al. Evaluation of heart rate recovery in patients with primary nephrotic syndrome. Hippokratia. 2015;19(2):109–113. pmid:27418757.

24. Meyer C, Schueller P, Balzer J, Lauer T, Westenfeld R, Schauerte P, et al. Sympathetic hyperactivity influences chemosensor function in patients with end-stage renal disease. Eur J Med Res. 2009;14 Suppl 4:151–155. pmid:20156747.

25. Salman IM. Cardiovascular autonomic dysfunction in chronic kidney disease: a comprehensive review. Curr Hypertens Rep. 2015;17(8):59. pmid:26071764.

26. Parashar R, Amir M, Pakhare A, Rathi P, Chaudhary L. Age related changes in autonomic functions. J Clin Diagn Res. 2016;10(3):CC11–15. pmid:27134865.

27. Stöggl TL, Björklund G. High intensity interval training leads to greater improvements in acute heart rate recovery and anaerobic power as high volume low intensity training. Front Physiol. 2017;8:562. pmid:28824457.

28. Wu SK, Lin YW, Chen CL, Tsai SW. Cardiac rehabilitation vs. home exercise after coronary artery bypass graft surgery: a comparison of heart rate recovery. Am J Phys Med Rehabil. 2006;85(9):711–717. pmid:16924183.

29. Dimopoulos S, Anastasiou-Nana M, Sakellariou D, Drakos S, Kapsimalakou S, Maroulidis G, et al.Effects of exercise rehabilitation program on heart rate recovery in patients with chronic heart failure. Eur J Cardiovasc Prev Rehabil. 2006;13(1):67–73. pmid:16449866.

30. Christou DD, Jones PP, Pimentel AE, Seals DR. Increased abdominal-to-peripheral fat distribution contributes to altered autonomic-circulatory control with human aging. Am J Physiol Heart Circ Physiol. 2004;287(4):H1530–1537. pmid:15178547.

31. Dart AM, Du X-J, Kingwell BA. Gender, sex hormones and autonomic nervous control of the cardiovascular system. Cardiovasc Res. 2002;53(3):678–687. pmid:11861039.

32. Kreier F, Fliers E, Voshol PJ, Van Eden CG, Havekes LM, Kalsbeek A, et al. Selective parasympathetic innervation of subcutaneous and intra-abdominal fat—functional implications. J Clin Invest. 2002;110(9):1243–1250. pmid:12417562.

33. Lin LY, Kuo HK, Lai LP, Lin JL, Tseng CD, Hwang JJ. Inverse correlation between heart rate recovery and metabolic risks in healthy children and adolescents: insight from the National Health and Nutrition Examination Survey 1999–2002. Diabetes Care. 2008;31(5):1015–1020. pmid:18268066.

34. Verma S, Alam R, Ahmad I, Singla D, Ali K, Hussain ME. Effect of glycemic control and disease duration on cardiac autonomic function and oxidative stress in type 2 diabetes mellitus. J Diabetes Metab Disord. 2018;17(2):149–158. pmid:30918849.

35. Liu Y, Liu SX, Zheng F, Cai Y, Xie KL, Zhang WL. Cardiovascular autonomic neuropathy in patients with type 2 diabetes. J Diabetes Investig. 2016;7(4):615–621. pmid:27181277.

36. Ambrose JA, Barua RS. The pathophysiology of cigarette smoking and cardiovascular disease: an update. J Am Coll Cardiol. 2004;43(10):1731–1737. pmid:15145091.

37. Cagirci G, Cay S, Karakurt O, Eryasar N, Kaya V, Canga A, et al. Influence of heavy cigarette smoking on heart rate variability and heart rate turbulence parameters. Ann Noninvasive Electrocardiol. 2009;14(4):327–332. pmid:19804508.

38. Erat M, Dogan M, Sunman H, Dinc Asarcikli L, Efe TH, Bilgin M, et al. Evaluation of heart rate recovery index in heavy smokers. Anatol J Cardiol. 2016;16(9):667–672. pmid:27488749.