Front. Physio.Frontiers in PhysiologyFront. Physio.1664-042XFrontiers Research Foundation10.3389/fphys.2011.00057PhysiologyOriginal ResearchHeart Rate Dynamics after Exercise in Cardiac Patients with and without Type 2 DiabetesNevesVictor R.12KiviniemiAntti M.1HautalaArto J.1KarjalainenJaana3PiiraOlli-Pekka3CataiAparecida M.2MäkikallioTimo H.3HuikuriHeikki Veli3TulppoMikko P.13*1Department of Exercise and Medical PhysiologyVerve, Oulu, Finland2Department of Physiotherapy, Federal University of São CarlosSão Carlos, Brazil3Department of Internal Medicine, Institute of Clinical Medicine, University of OuluOulu, Finland
Edited by: Jaakko Hartiala, Turku University Hospital, Finland
Reviewed by: Antti Saraste, University of Turku, Finland; Petri Haapalahti, HUSLAB, Finland; Tomi Laitinen, University of Eastern Finland and Kuopio University Hospital, Finland
*Correspondence: Mikko P. Tulppo, Department of Exercise and Medical Physiology, Verve, P.O. Box 404, FI-90101 Oulu, Finland. e-mail: mikko.tulppo@verve.fi
This article was submitted to Frontiers in Clinical and Translational Physiology, a specialty of Frontiers in Physiology.
This is an open-access article subject to a non-exclusive license between the authors and Frontiers Media SA, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and other Frontiers conditions are complied with.
Purpose: The incidence of cardiovascular events is higher in coronary artery disease patients with type 2 diabetes (CAD + T2D) than in CAD patients without T2D. There is increasing evidence that the recovery phase after exercise is a vulnerable phase for various cardiovascular events. We hypothesized that autonomic regulation differs in CAD patients with and without T2D during post-exercise condition. Methods: A symptom-limited maximal exercise test on a bicycle ergometer was performed for 68 CAD + T2D patients (age 61 ± 5 years, 78% males, ejection fraction (EF) 67 ± 8, 100% on β-blockade), and 64 CAD patients (age 62 ± 5 years, 80% males, EF 64 ± 8, 100% on β-blockade). Heart rate (HR) recovery after exercise was calculated as the slope of HR during the first 60 s after cessation of exercise (HRRslope). R–R intervals were measured before (5 min) and after exercise from 3 to 8 min, both in a supine position. R–R intervals were analyzed using time and frequency methods and a detrended fluctuation method (α1). Results: BMI was 30 ± 4 vs. 27 ± 3 kg m2 (p < 0.001); maximal exercise capacity, 6.5 ± 1.7 vs. 7.7 ± 1.9 METs (p < 0.001); maximal HR, 128 ± 19 vs. 132 ± 18 bpm (p = ns); and HRRslope, −0.53 ± 0.17 vs. −0.62 ± 0.15 beats/s (p = 0.004), for CAD patients with and without T2D, respectively. There was no differences between the groups in HRRslope after adjustment for METs, BMI, and medication (ANCOVA, p = 0.228 for T2D and, e.g., p = 0.030 for METs). CAD + T2D patients had a higher HR at rest than non-diabetic patients (57 ± 10 vs. 54 ± 6 bpm, p = 0.030), but no other differences were observed in HR dynamics at rest or in post-exercise condition. Conclusion: HR recovery is delayed in CAD + T2D patients, suggesting impairment of vagal activity and/or augmented sympathetic activity after exercise. Blunted HR recovery after exercise in diabetic patients compared with non-diabetic patients is more closely related to low exercise capacity and obesity than to T2D itself.
The incidence of cardiovascular events is higher in coronary artery disease patients (CAD) with type 2 diabetes (T2D) than in CAD patients without T2D (Haffner et al., 1998; Junttila et al., 2010), but the physiological or pathophysiological mechanisms causing these differences are not well known. Altered autonomic regulation is one potential mechanism resulting in the increased number of cardiovascular events in CAD + T2D patients (Okada et al., 2010; Pop-Busui et al., 2010; Lanza et al., 2011). Autonomic regulation can be studied using various methods at the laboratory, or ambulatory condition tests such as passive head-up tilt and handgrip tests (Montano et al., 1994; Tulppo et al., 2001, 2005; Fu et al., 2002), or ambulatory heart rate (HR) variability and blood pressure measurements (Pagani et al., 1985; Kleiger et al., 1987; Piira et al., 2011). Analysis of autonomic regulation during and after exercise has also been used in various physiological (Yamamoto et al., 1991; Gregoire et al., 1996; Tulppo et al., 1996, 1998b) and clinical settings (Cole et al., 1999, 2000; Jouven and Ducimetiere, 2000; Jouven et al., 2005; Kiviniemi et al., 2011).
There is increasing evidence that the recovery phase after exercise is a vulnerable phase for various cardiovascular events. Case-crossover studies have shown that exercise as a trigger of acute myocardial infarction is not limited to the time of exercise, but extends for a certain time period after cessation of physical activity (Siscovick et al., 1982, 1984; Albert et al., 2000; von Klot et al., 2008). Similarly, the risk of sudden cardiac death is transiently increased in the 30-min after vigorous exercise, and atrial fibrillation episodes occur more commonly after rather than before exercise (Siscovick et al., 1982, 1984; Coumel, 1994; Albert et al., 2000; Huikuri, 2008). Measurement of autonomic function in the very early phase of recovery after exercise has also provided prognostic information. For example, delayed HR recovery 1–2 min after exercise has been shown to predict cardiovascular events in the general population and in various patient groups and animal studies (Cole et al., 1999; Lauer and Froelicher, 2002; Nissinen et al., 2003; Jouven et al., 2005; Smith et al., 2005). Recently, we have shown that co-activation of both autonomic arms may occur during the recovery phase of exercise (Tulppo et al., 2011), which may partly explain the clustering of various cardiovascular events in the recovery phase of exercise.
The present research was designed to study the behavior of HR dynamics during exercise and in the recovery phase of exercise in CAD patients with and without T2D, matched for age, ejection fraction (EF), and sex. We hypothesized that autonomic regulation, measured by HR recovery and HR variability methods, differs in CAD patients with and without T2D, particularly during post-exercise condition.
Materials and MethodsSubjects and study protocol
The present study was conducted in the Department of Exercise and Medical Physiology at Verve (Oulu, Finland) and Oulu University Hospital, and the patients were selected from patients in the ARTEMIS (Innovation to Reduce Cardiovascular Complications of Diabetes at the Intersection) study database who had stable CAD without (n = 64) and with T2D (n = 68). The exclusion criteria were inability to perform the exercise stress test, unstable angina at the time of recruitment, advanced age (>75 years), a recent (<6 months) myocardial infarction, severe nephropathy, heart failure, scheduled cardiac revascularization therapy, T2D diagnosed less than 6 months from the beginning of study, diabetes autonomic neuropathy, dementia, alcoholism, drug abuse, or any condition that could impair the subject’s capacity to give informed consent. CAD was documented by coronary angiography and T2D was verified by oral glucose tolerance test according to current recommendations (WHO, 1999). The study was performed according to the Declaration of Helsinki, the local committees of research ethics in the Northern Ostrobothnia Hospital District (Oulu, Finland) approved the protocol, and all the subjects gave written informed consent.
The laboratory measurements were performed in the Department of Exercise and Medical Physiology at Verve (Oulu, Finland). The patients were not allowed to eat or consume caffeine for 3 h before the tests. Physical exercise and use of alcohol were prohibited for 24 h before testing. Electrocardiography (ECG) and R–R intervals were collected for 10 min at supine rest, during exercise, and 10 min after exercise in supine position. Breathing was spontaneous in all phases. A capillary blood sample was obtained for analysis of HbA1c concentration before testing (Afinion™ AS100, Axis-Shield PoC AS, Oslo, Norway).
Exercise stress test
The patients performed an incremental maximal test on a bicycle ergometer (Monark Ergomedic 839 E; Monark Exercise AB, Vansbro, Sweden), starting at 30 W with the work rate increasing at a rate of 10 and 15 W every 1 min until exhaustion for females and males, respectively. The patients moved to the supine position within 30 s after cessation of exercise. The patients were not allowed to move or talk during the recovery phase.
Ventilation (VE) and gas exchange (M909 Ergospirometer, Medikro, Kuopio, Finland) were measured and reported as the mean value for every minute. The highest 1-min mean value of oxygen consumption was expressed as the peak oxygen consumption (VO2peak). Maximal workload (W) and maximal METs were calculated as average workload and METs during the last 1 min of the test. ECG was monitored and recorded using a standard 12-lead ECG (GE Healthcare, Cam-14, Waukesha, WI, USA) and at the same time R–R intervals were recorded with a Polar R–R recorder with a sampling frequency of 1,000 Hz (Polar Electro, Kempele, Finland). Blood pressure was measured with an electronic sphygmomanometer (Tango, Sun-Tech, Raleigh, NC, USA) at rest and during exercise testing. The patients were encouraged to reach a symptom-limited maximal workload, and exercise was stopped if ST depression was >0.2 mV.
The maximal HR was calculated as a mean value of 5 R–R intervals before cessation of exercise from Polar R–R interval data and converted to maximal HR in beats/min with the following equation: maximal HR (beats/min) = 60·(mean R–R interval in ms/1000)−1. The chronotropic response to exercise was calculated by chronotropic response index (CRI) with the following equation: CRI = 100·(maximal HR − resting HR)·(220 − age − resting HR)−1 (Kiviniemi et al., 2011), and HR reserve as maximal HR − resting HR. HR recovery values were calculated as mean values of 5 R–R intervals around the time points of 15 (HRR15), 30 (HRR30), 60 (HRR60), and 120 (HRR120) s after cessation of exercise from Polar data. Thereafter, HR recovery (beats/min) was calculated as the change in HR from maximal HR to recovery HR at the time points of 15, 30, 60, and 120 s after cessation of exercise. The slope of HR during the 60-s after cessation of exercise (HRRslope) was calculated by linear model using above described HR values at the maximum and 15, 30, and 60 s time points (Figure 1).
Example of heart rate recovery as the slope of heart rate during the first 60 s after cessation of exercise for two patients (with and without type 2 diabetes) with similar age and maximal heart rate. The slope was calculated by linear model using mean values of 5 R–R intervals [converted to the HR values in beats/min as 60·(mean R–R interval in ms/100)−1] at maximum and around the time points of 15, 30, and 60 s after cessation of exercise.
Heart rate variability
R–R intervals were analyzed from the last 5-min period at rest, at every load during exercise (1-min periods), and from 3 to 8 min (5-min) after exercise. Before HR variability analysis (using Polar R–R data), all of the R–R intervals were edited manually to exclude all premature beats and noise, which accounted for <8% in every subject. ECG data was used during the editing to confirm sinus origin of the beats. An autoregressive model was used to estimate the power spectrum densities of HR variability after period level detrending (Huikuri et al., 1996; Tulppo et al., 1996). The power spectra were quantified by measuring the area under two frequency bands: LF power (0.04–0.15 Hz) and HF power (0.15–0.4 Hz). Details of the spectrum analyses have been described previously elsewhere (Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 1996). Detrended fluctuation analysis (DFA) method was used to calculate fractal-like correlation properties of the R–R interval data (Peng et al., 1995; Iyengar et al., 1996). In this study the short-term (4–11 beats) scaling exponent (α-1) was calculated based on previous experiments (Mäkikallio et al., 1996).
Statistics
The data were presented as mean ± SD. The normal Gaussian distribution of the data was verified with the Kolmogorov–Smirnov goodness-of-fit test (z value > 1.0). For repeated measurements of post-exercise HR recovery (HRR15, HRR30, HRR60, HRR120), two-factor ANOVA (time × group) was used to assess main effects followed by post hoc comparison (unpaired t-test). Significant post hoc differences in HR recovery between the groups were adjusted for maximal METs, BMI, and medications (calcium antagonists and nitrates) using ANCOVA. For other variables, the unpaired Student’s t-test was used to assess the differences between CAD and CAD + T2D groups in normally distributed data. The Mann–Whitney U-test was used to assess the difference in HR variability indices between groups. The Chi-square test was used for categorical variables. The data were analyzed using SPSS software (SPSS 19.0, SPSS Inc., Chicago, USA). A p-value < 0.05 was considered statistically significant.
ResultsStudy population
The characteristics of the study population, including comparisons between diabetic and non-diabetic patients, are given in Table 1. The diabetic patients’ body composition, e.g., weight (p < 0.001), waist circumference (p < 0.001), BMI (p < 0.001), and HbA1c (p < 0.001) differed from that of the non-diabetic patients. However, blood pressure, EF, history of infarction, and revascularization or smoking did not differ between the diabetic and non-diabetic patients (Table 1). Medication did not differ between the groups, except that calcium antagonists, nitrates, and antidiabetics were more common among the diabetic than the non-diabetic patients (Table 1).
Characteristics of patients.
CAD + T2D (n = 68)
CAD (n = 64)
p-Level
Sex (F/M)
15/53
13/51
p = 0.488
Age (years)
61 ± 5
62 ± 5
p = 0.514
Height (cm)
171 ± 7
171 ± 8
p = 0.678
Weight (kg)
89 ± 15
80 ± 12
p < 0.001
Waist (cm)
105 ± 12
94 ± 10
p < 0.001
Hip (cm)
105 ± 10
99 ± 6
p < 0.001
Waist/Hip ratio
1.00 ± 0.07
0.95 ± 0.07
p < 0.001
BMI (kg m−2)
30 ± 4
27 ± 3
p < 0.001
HbA1c (%)
6.0 ± 0.7
5.5 ± 0.2
p < 0.001
SBP (mmHg)
138 ± 18
136 ± 16
p = 0.672
DBP (mmHg)
82 ± 10
82 ± 8
p = 0.875
HISTORY OF AMI
NSTEMI
18 (26%)
19 (30%)
P = 0.469
STEMI
13 (19%)
14 (21%)
p = 0.662
REVASCULARIZATION
CAGB
14 (21%)
13 (28%)
p = 0.506
PCI
42 (67%)
27 (59%)
p = 0.426
Current smokers
3 (4%)
8 (12%)
p = 0.147
LVEF (%)
67 ± 8
65 ± 8
p = 0.107
MEDICATION
Anticougulants
66 (97%)
63 (98%)
p = 1.000
Beta blockers
68 (100%)
64 (100%)
p = 1.000
Calcium antagonists
20 (29%)
9 (14%)
p = 0.037
ACEI
31 (46%)
23 (36%)
p = 0.291
AT2
13 (19%)
14 (22%)
p = 0.830
Diuretics
26 (38%)
15 (23%)
p = 0.090
Statin
63 (93%)
57 (89%)
p = 0.553
Insulines
10 (7%)
0 (0%)
p = 0.001
Oral antidiabetics
50 (74%)
0 (0%)
p < 0.001
Nitrates
28 (41%)
8 (13%)
p < 0.001
Arrhythmia
1 (1%)
1 (1%)
p = 1.000
Values are mean ± SD. BMI, body mass index; HbA1c, Glycosylated hemoglobin; SBP, systolic blood pressure; DBP, diastolic blood pressure; HR, heart rate; AMI, acute myocardial infarction; NSTEMI, no-ST segment elevation myocardial infarction; STEMI, ST segment elevation myocardial infarction; revascularized, the patients who had at least one of the procedures (CABG, coronary artery by-pass grafting or PCI, Percutaneus coronary intervention); LVEF, left ventricular ejection fraction; ACEI, angiotensin conversion enzyme inhibitor; AT2, angiotensin II receptor blocker; CAD, coronary artery disease group; CAD + T2D, coronary artery disease with type 2 diabetes.
HR response to maximal exercise
The results of the maximal exercise test, including comparisons between diabetic and non-diabetic patients, are given in Table 2. Maximal exercise capacity was lower (p < 0.001) in the diabetes patients, but maximal HR (p = 0.178) was at the same level compared with the non-diabetic patients (Table 2). HR reserve was smaller (p = 0.027) and CRI tended to be smaller (p = 0.057) in the diabetic patients than in the non-diabetic patients. However, there were not significant differences between groups in HR reserve or CRI after adjustment for maximal METs, BMI, and medication, e.g., for HR reserve ANCOVA p = 0.686, p < 0.001, p = 0.126, and p = 0.169 for diabetes, Mets, BMI, and calcium antagonists, respectively. There was no difference between the groups in the prevalence of ST segment depression >0.1 mV during the exercise test.
Maximal exercise capacity and heart rate recovery.
CAD + T2D (n = 68)
CAD (n = 64)
p-Level
Max load (W)
139 ± 41
155 ± 43
p = 0.033
Max METs
6.5 ± 1.7
7.7 ± 1.9
p < 0.001
MAX VO2PEAK
ml kg−1 min−1
21.6 ± 5.8
25.9 ± 6.5
p < 0.001
l min−1
1.89 ± 0.52
2.05 ± 0.58
p = 0.101
Max ST depression > 1 mm
32 (47%)
57 (49%)
p = 0.228
Max HR (bpm)
128 ± 19
132 ± 18
p = 0.178
HR reserve (bpm)
70 ± 20
78 ± 18
p = 0.027
CRI (bpm)
69 ± 18
75 ± 17
p = 0.057
HR RECOVERY
15 s (bpm)
10.9 ± 4.6
10.8 ± 5.0
p = 0.376
30 s (bpm)
15.8 ± 7.0
17.8 ± 6.9
p = 0.098
60 s (bpm)
32.7 ± 10.5
37.8 ± 9.7
p = 0.006
120 s (bpm)
44.3 ± 13.4
48.9 ± 11.1
p = 0.037
Slope 60 (beats/s)
−0.53 ± 0.17
−0.62 ± 0.14
p = 0.004
Values are mean ± SD. METs, metabolic equivalents; VO2peak, peak uptake oxygen; HR, heart rate; CRI, maximal chronotropic response index.
HR recovery
Type 2 diabetes modified post-exercise HR recovery (main effect for time × group interaction p = 0.037). The diabetic patients had delayed HRR60 (p = 0.006), HRR120 (p = 0.037), and HRRslope (p = 0.004) compared with the non-diabetic patients (Table 2). However, there were not significant differences between groups in HRR60, HRR120, or HRRslope after adjustment for maximal METs, BMI, and medication, e.g., for HRR60 ANCOVA p = 0.223, p = 0.061, p = 0.387, and p = 0.094 for diabetes, Mets, BMI, and calcium antagonists, respectively. The corresponding p values for HRR120 were p = 0.980, p = 0.001, p = 0.451, and p = 0.006 and for HRRslopep = 0.228, p = 0.030, p = 0.404, and p = 0.079 for diabetes, Mets, BMI, and calcium antagonists, respectively. Nitrates did not modify HR behavior between groups at any condition. Examples of HRRslope for non-diabetic and diabetic subjects are shown in Figure 1.
HR and HR variability at rest, during submaximal exercise and after exercise
Heart rate and HR variability before and after exercise, including comparisons between diabetic and non-diabetic patients, are given in Table 3. Five diabetic patients were excluded from the analysis due to the significant number of extra systoles or technical artifacts in the R–R interval data. One non-diabetic patient was excluded due to the extra systoles in both pre- and post-exercise conditions and one in post-exercise condition due to the same reason.
Average values of linear and non-linear heart rate variability before exercise (5 min) and in post-exercise condition from 3 to 8 min after exercise, both in supine position.
CAD + T2D
CAD
p-Level
Pre-exercise 5 min
n = 63
n = 63
HR (bpm)
57 ± 10
54 ± 6
p = 0.030
SDNN (ms)
36 ± 19
36 ± 21
p = 0.969
HF power (ms2)
299 ± 353
479 ± 1347
p = 0.490
LF power (ms2)
385 ± 521
409 ± 563
p = 0.940
LF/HF ratio
1.81 ± 1.71
1.47 ± 1.34
p = 0.210
α-1
1.04 ± 0.27
0.99 ± 0.26
p = 0.280
Post-exercise 5 min
n = 63
n = 62
HR (bpm)
76 ± 10
75 ± 10
p = 0.765
SDNN (ms)
29 ± 12
28 ± 17
p = 0.386
HF power (ms2)
102 ± 172
114 ± 246
p = 0.268
LF power (ms2)
162 ± 212
175 ± 352
p = 0.608
LF/HF ratio
2.8 ± 2.3
3.9 ± 3.9
p = 0.257
α-1
1.18 ± 0.30
1.24 ± 0.28
p = 0.277
Values are mean ± SD. HR, heart rate; SDNN, SD of normal-to-normal R–R intervals; HF, high frequency power of R–R intervals; LF, low frequency power of R–R intervals; α-1, fractal scaling exponent of R–R intervals.
The diabetic patients had a higher HR in resting condition than the non-diabetic patients (p = 0.030), but no other differences were observed at rest or in post-exercise condition (Table 3). HR or HR variability did not differ between the groups during submaximal exercise at the levels of 40, 60, or 80% of maximal oxygen uptake, e.g., HR was 92 ± 12 vs. 95 ± 14 bpm (p = 0.099) and α1 1.02 ± 0.37 vs. 1.02 ± 0.32 (p = 0.93) at the level of 61 ± 5% and 61 ± 5% of maximal oxygen uptake for diabetic and non-diabetic patients, respectively.
Discussion
The present study showed that post-exercise HR recovery was slower in CAD + T2D patients than in CAD patients, suggesting an impairment of vagal modulation and/or augmented sympathetic activity initially after exercise for diabetic patients. The impaired HR recovery in diabetic patients was the most obvious 1 min after exercise, documented by HRR60 and HRRslope indices. However, blunted HR recovery initially after exercise in diabetic patients compared with non-diabetic patients was more closely related to low exercise capacity and obesity than to T2D itself. Taken together, these findings suggest that delayed HR recovery after exercise in T2D patients may be reversible through prevention and treatment, including, e.g., regular exercise and weight management.
Post-exercise HR dynamics
The interplay between sympathetic and vagal regulation of HR during exercise is organized in a reciprocal fashion, i.e., increased sympathetic activity is accompanied by decreased vagal activity in the heart during dynamic exercise (Robinson et al., 1966; Maciel et al., 1986; Orizio et al., 1988; Yamamoto and Hughson, 1991; Tulppo et al., 1996, 1998b). However, this reciprocal behavior is altered in the recovery phase after exercise due to temporal differences in the recovery pattern of the autonomic arms in post-exercise condition (Tulppo et al., 2011). A rapid restoration of vagal activity occurs after cessation of exercise (Imai et al., 1994; Goldberger et al., 2006; Martinmaki and Rusko, 2008; Tulppo et al., 2011). On the contrary, the sympathetic nervous system seems to have a longer latency to return to the baseline after cessation of exercise, resulting in long-lasting hyperactivity of sympathetic activity (Ray, 1993; Tulppo et al., 2011). Taken together, these changes in autonomic regulation may result in dual activation of the sympathetic and vagal arms in post-exercise condition (Tulppo et al., 2011). These changes in autonomic activity in post-exercise condition may partly explain the clinical findings, since particularly the recovery phase of exercise has been shown to be a vulnerable phase for various cardiovascular events (Siscovick et al., 1982, 1984; Albert et al., 2000; von Klot et al., 2008).
It is well known that CAD + T2D patients are at a higher risk for cardiac events than CAD patients without diabetes. Altered autonomic regulation in post-exercise condition is one potential mechanism, since slow HR recovery after exercise has also been associated with cardiovascular events in various clinical and subclinical populations (Cole et al., 1999; Lauer and Froelicher, 2002; Nissinen et al., 2003; Jouven et al., 2005). In the present study, HR recovery from 1 to 2 min after exercise was the only marker of autonomic activity separating diabetic and non-diabetic cardiac patients matched with age, sex, and EF all in optimal medications, including β-blockade. Since a complex interaction of autonomic regulation occurs in the initial phase after exercise, it is difficult to detect the difference in HR recovery between diabetic and non-diabetic patients due to impaired vagal activation or augmented sympathetic activation or both. Also circulating catecholamines have an important contribution for the sympathetic post-exercise hyperactivity (Krock and Hartung, 1992), but unfortunately we did not measure epinephrine or norepinephrine in the present study. Thirdly, as a novel finding of the present study, there were no differences between patients groups in HR recovery after adjustment for METs and BMI. These findings emphasized regular exercise training and weight management as potential treatments to improve post-exercise HR recovery in CAD + T2D patients. Importantly, also calcium antagonists modified HR recovery particularly 2 min after exercise. Patients with calcium antagonists usually have more severe hypertension which is associated with delayed HR recovery (Carnethon et al., 2011).
HR response to maximal exercise
Measurement of the chronotropic response of HR to exercise has been used to assess particularly sympathetic influence on the heart and it has been shown to be a powerful predictor of cardiac mortality in both asymptomatic populations (Azarbal et al., 2004; Gulati et al., 2005; Jouven et al., 2005; Savonen et al., 2006; Kiviniemi et al., 2011) and in patients with various cardiac diseases (Myers et al., 2007; Savonen et al., 2008; Kiviniemi et al., 2011). In the present study HR reserve was lower and the maximal chronotropic response adjusted for age (CRI) tended to be lower in CAD + T2D patients than in CAD patients without diabetes. However, there were no differences between patients groups in these indices after adjustment for fitness and BMI.
HR variability
Type 2 diabetes has been shown to decrease HR variability and baroreflex sensitivity in T2D patients without CAD (Masaoka et al., 1985; Frattola et al., 1997). There are previous studies where T2D has shown no additional decrease in baroreflex sensitivity among patients with CAD (Wykretowicz et al., 2005) or in HR variability among T2D patients with heart failure or CAD (Burger and Aronson, 2001; Kiviniemi et al., 2010). Also in the present study the short-term HR variability indices measured at resting or post-exercise condition were not able to separate CAD patients with and without T2D. There are several potential explanations for these findings. First, the populations are carefully matched according to age, sex, EF, and β-blockade which all are known to effect on HR variability and may partly explain the present findings. Secondly, short-term HR variability measures at laboratory condition are not so well reproducible than 24-h recordings (Huikuri et al., 1990; Tulppo et al., 1998b) and may be influenced by “white coat” effects in some subjects (Grassi et al., 1999). Thirdly, high level of circulating norepinephrine particularly in post-exercise condition may results in abrupt changes in beat-to-beat R–R interval dynamics which are not detected by used HR variability methods (Tulppo et al., 1998a). On the contrary to HR variability measurements at rest or few minutes after exercise, rapid and marked changes occurs in both autonomic arms in tens of seconds from sympathetic dominance to restoration of vagal activity resulting decreased HR after exercise (Tulppo et al., 2011). The unique behavior of autonomic regulation initially after maximal exercise is the major candidate to explain differences in HR recovery between diabetic and non-diabetic patient groups.
Study limitations
The measurements of the present study were performed under continued prescribed medication because of ethical reasons and the well known withdrawal effect of beta-block cessation. However, the present results will have more practical implications when the analyses are performed at a time when the patients are on their normal medication.
Clinical applications
Based on the present study, it is important to emphasize the role of physical exercise in maintaining good cardiovascular function, improving physical performance, and losing weight. Although all the patients were on optimized medication, it is important to underscore that behavioral modification, exercises, and diet should be given priority similar to other preventive medications (Chow et al., 2010).
Conclusion
Heart rate recovery 1 min after exercise is the most powerful index separating CAD patients with and without T2D in the modern treatment era, beyond short-term HR variability measurements at rest, during exercise, or in post-exercise condition. HR recovery is delayed in CAD + T2D patients compared with CAD patients without T2D, suggesting impairment of vagal activity and/or augmented sympathetic activity after exercise. However, blunted HR recovery initially after exercise in diabetic patients compared with non-diabetic patients is more closely related to exercise capacity and obesity than to T2D itself.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by grants from the Finnish Technology Development Centre (TEKES. Helsinki. Finland), Paavo Nurmi Foundation, Turku, Finland, Polar Electro Oy, Kempele, Finland and Hur Oy, Kokkola, Finland, and the National Council for Scientific and Technological Development, Brasília, Brazil. (proc. 200630/2010-5).
ReferencesAlbertC. M.MittlemanM. A.ChaeC. U.LeeI. M.HennekensC. H.MansonJ. E. (2000). Triggering of sudden death from cardiac causes by vigorous exertion. 343, 1355–1361.10.1056/NEJM20001109343190211070099AzarbalB.HayesS. W.LewinH. C.HachamovitchR.CohenI.BermanD. S. (2004). The incremental prognostic value of percentage of heart rate reserve achieved over myocardial perfusion single-photon emission computed tomography in the prediction of cardiac death and all-cause mortality: superiority over 85% of maximal age-predicted heart rate. 44, 423–430.10.1016/j.jacc.2004.02.06015261942BurgerA. J.AronsonD. (2001). Effect of diabetes mellitus on heart rate variability in patients with congestive heart failure. 24, 53–59.10.1046/j.1460-9592.2001.00053.x11227970CarnethonM. R.SternfeldB.LiuK.JacobsD. R.Jr.SchreinerP. J.WilliamsO. D.LewisC. E.SidneyS. (2011). Correlates of heart rate recovery over 20 years in a healthy population sample. . [Epub ahead of print].10.1249/MSS.0b013e31822cb19021796053ChowC. K.JollyS.Rao-MelaciniP.FoxK. A.AnandS. S.YusufS. (2010). Association of diet, exercise, and smoking modification with risk of early cardiovascular events after acute coronary syndromes. 121, 750–758.10.1161/CIRCULATIONAHA.109.89152320124123ColeC. R.BlackstoneE. H.PashkowF. J.SnaderC. E.LauerM. S. (1999). Heart-rate recovery immediately after exercise as a predictor of mortality. 341, 1351–1357.10.1056/NEJM19991028341180410536127ColeC. R.FoodyJ. M.BlackstoneE. H.LauerM. S. (2000). Heart rate recovery after submaximal exercise testing as a predictor of mortality in a cardiovascularly healthy cohort. 132, 552–555.10744592CoumelP. (1994). Paroxysmal atrial fibrillation: a disorder of autonomic tone?15(Suppl. A), 9–16.8070496FrattolaA.ParatiG.GambaP.PaleariF.MauriG.Di RienzoM.CastiglioniP.ManciaG. (1997). Time and frequency domain estimates of spontaneous baroreflex sensitivity provide early detection of autonomic dysfunction in diabetes mellitus. 40, 1470–1475.10.1007/s0012500508519447956FuQ.LevineB. D.PawelczykJ. A.ErtlA. C.DiedrichA.CoxJ. F.ZuckermanJ. H.RayC. A.SmithM. L.IwaseS.SaitoM.SugiyamaY.ManoT.ZhangR.IwasakiK.LaneL. D.BuckeyJ. C.Jr.CookeW. H.RobertsonR. M.BaischF. J.BlomqvistC. G.EckbergD. L.RobertsonD.BiaggioniI. (2002). Cardiovascular and sympathetic neural responses to handgrip and cold pressor stimuli in humans before, during and after spaceflight. 544, 653–664.10.1113/jphysiol.2002.02337412381834GoldbergerJ. J.LeF. K.LahiriM.KannankerilP. J.NgJ.KadishA. H. (2006). Assessment of parasympathetic reactivation after exercise. 290, H2446–H2452.10.1152/ajpheart.01118.200516415073GrassiG.TurriC.VailatiS.Dell’oroR.ManciaG. (1999). Muscle and skin sympathetic nerve traffic during the “white-coat” effect. 100, 222–225.10411843GregoireJ.TuckS.YamamotoY.HughsonR. L. (1996). Heart rate variability at rest and exercise: influence of age, gender, and physical training. 21, 455–470.10.1139/h96-0408959312GulatiM.BlackH. R.ShawL. J.ArnsdorfM. F.MerzC. N.LauerM. S.MarwickT. H.PandeyD. K.WicklundR. H.ThistedR. A. (2005). The prognostic value of a nomogram for exercise capacity in women. 353, 468–475.10.1056/NEJMoa04415416079370HaffnerS. M.LehtoS.RönnemaaT.PyoräläK.LaaksoM. (1998). Mortality from coronary heart disease in subjects with type 2 diabetes and in nondiabetic subjects with and without prior myocardial infarction. 339, 229–234.10.1056/NEJM1998072333904049673301HuikuriH. V. (2008). Heart rate dynamics as a marker of vulnerability to atrial fibrillation. 19, 913–914.10.1111/j.1540-8167.2008.01197.x18479324HuikuriH. V.KesslerK. M.TerracallE.CastellanosA.LinnaluotoM. K.MyerburgR. J. (1990). Reproducibility and circadian rhythm of heart rate variability in healthy subjects. 65, 391–393.10.1016/0002-9149(90)90308-N2301268HuikuriH. V.SeppanenT.KoistinenM. J.AiraksinenJ.IkaheimoM. J.CastellanosA.MyerburgR. J. (1996). Abnormalities in beat-to-beat dynamics of heart rate before the spontaneous onset of life-threatening ventricular tachyarrhythmias in patients with prior myocardial infarction. 93, 1836–1844.8635263ImaiK.SatoH.HoriM.KusuokaH.OzakiH.YokoyamaH.TakedaH.InoueM.KamadaT. (1994). Vagally mediated heart rate recovery after exercise is accelerated in athletes but blunted in patients with chronic heart failure. 24, 1529–1535.10.1016/0735-1097(94)90150-37930286IyengarN.PengC. K.MorinR.GoldbergerA. L.LipsitzL. A. (1996). Age-related alterations in the fractal scaling of cardiac interbeat interval dynamics. 271, R1078–R1084.8898003JouvenX.DucimetiereP. (2000). Recovery of heart rate after exercise. 342, 662–663.10.1056/NEJM20000302342091510702064JouvenX.EmpanaJ. P.SchwartzP. J.DesnosM.CourbonD.DucimetiereP. (2005). Heart-rate profile during exercise as a predictor of sudden death. 352, 1951–1958.10.1056/NEJMoa04301215888695JunttilaM. J.BarthelP.MyerburgR. J.MakikallioT. H.BauerA.UlmK.KiviniemiA.TulppoM.PerkiomakiJ. S.SchmidtG.HuikuriH. V. (2010). Sudden cardiac death after myocardial infarction in patients with type 2 diabetes. 7, 1396–1403.10.1016/j.hrthm.2010.07.03120682359KiviniemiA. M.TiinanenS.HautalaA. J.SeppanenT.NortonK. N.FrancesM. F.NolanR. P.HuikuriH. V.TulppoM. P.ShoemakerJ. K. (2010). Low-frequency oscillations in R-R interval and blood pressure across the continuum of cardiovascular risk. 158, 92–99.10.1016/j.autneu.2010.05.01020573550KiviniemiA. M.TulppoM. P.HautalaA. J.MakikallioT. H.PerkiomakiJ. S.SeppanenT.HuikuriH. V. (2011). Long-term outcome of patients with chronotropic incompetence after an acute myocardial infarction. 43, 33–39.10.3109/07853890.2010.52176420977382KleigerR. E.MillerJ. P.BiggerJ. T.Jr.MossA. J. (1987). Decreased heart rate variability and its association with increased mortality after acute myocardial infarction. 59, 256–262.10.1016/0002-9149(87)90795-83812275KrockL. P.HartungG. H. (1992). Influence of post-exercise activity on plasma catecholamines, blood pressure and heart rate in normal subjects. 2, 89–97.10.1007/BF018196631638110LanzaG. A.BaroneL.ScaloneG.PitoccoD.SguegliaG. A.MolloR.NerlaR.ZaccardiF.GhirlandaG.CreaF. (2011). Inflammation-related effects of adjuvant influenza A vaccination on platelet activation and cardiac autonomic function. 269, 118–125.10.1111/j.1365-2796.2010.02285.x20964738LauerM. S.FroelicherV. (2002). Abnormal heart-rate recovery after exercise. 360, 1176–1177.10.1016/S0140-6736(02)11223-212387983MacielB. C.GalloL.Jr.Marin NetoJ. A.Lima FilhoE. C.MartinsL. E. (1986). Autonomic nervous control of the heart rate during dynamic exercise in normal man. 71, 457–460.3757442MäkikallioT. H.SeppänenT.NiemeläM.AiraksinenK. E.TulppoM.HuikuriH. V. (1996). Abnormalities in beat to beat complexity of heart rate dynamics in patients with a previous myocardial infarction. 28, 1005–1011.10.1016/S0735-1097(96)00243-4MartinmakiK.RuskoH. (2008). Time-frequency analysis of heart rate variability during immediate recovery from low and high intensity exercise. 102, 353–360.10.1007/s00421-007-0594-517943304MasaokaS.Lev-RanA.HillL. R.VakilG.HonE. H. (1985). Heart rate variability in diabetes: relationship to age and duration of the disease. 8, 64–68.10.2337/diacare.8.1.643971850MontanoN.RusconeT. G.PortaA.LombardiF.PaganiM.MallianiA. (1994). Power spectrum analysis of heart rate variability to assess the changes in sympathovagal balance during graded orthostatic tilt. 90, 1826–1831.7923668MyersJ.TanS. Y.AbellaJ.AletiV.FroelicherV. F. (2007). Comparison of the chronotropic response to exercise and heart rate recovery in predicting cardiovascular mortality. 14, 215–221.10.1097/HJR.0b013e328088cb9217446799NissinenS. I.MakikallioT. H.SeppanenT.TapanainenJ. M.SaloM.TulppoM. P.HuikuriH. V. (2003). Heart rate recovery after exercise as a predictor of mortality among survivors of acute myocardial infarction. 91, 711–714.10.1016/S0002-9149(02)03410-012633804OkadaN.TakahashiN.YufuK.MurozonoY.WakisakaO.ShinoharaT.AnanF.NakagawaM.HaraM.SaikawaT.YoshimatsuH. (2010). Baroreflex sensitivity predicts cardiovascular events in patients with type 2 diabetes mellitus without structural heart disease. 74, 1379–1383.10.1253/circj.CJ-10-085020453396OrizioC.PeriniR.ComandeA.CastellanoM.BeschiM.VeicsteinasA. (1988). Plasma catecholamines and heart rate at the beginning of muscular exercise in man. 57, 644–651.10.1007/BF004184773396585PaganiM.FurlanR.LombardiF.PizzinelliP.LanziG.CastelliP.CeruttiS.SantoliC.MallianiA. (1985). Technique for 24 hour recording of continuous high fidelity arterial pressure and electrocardiogram in ambulatory patients. 7, 401–405.10.3109/106419685090735644039993PengC. K.HavlinS.StanleyH. E.GoldbergerA. L. (1995). Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series. 5, 82–87.10.1063/1.16614111538314PiiraO. P.HuikuriH. V.TulppoM. P. (2011). Effects of emotional excitement on heart rate and blood pressure dynamics in patients with coronary artery disease. 160, 107–114.10.1016/j.autneu.2010.11.01121183409Pop-BusuiR.EvansG. W.GersteinH. C.FonsecaV.FlegJ. L.HoogwerfB. J.GenuthS.GrimmR. H.CorsonM. A.PrineasR. (2010). Effects of cardiac autonomic dysfunction on mortality risk in the action to control cardiovascular risk in diabetes (ACCORD) trial. 33, 1578–1584.10.2337/dc10-012520215456RayC. A. (1993). Muscle sympathetic nerve responses to prolonged one-legged exercise. 74, 1719–1722.8514687RobinsonB. F.EpsteinS. E.BeiserG. D.BraunwaldE. (1966). Control of heart rate by the autonomic nervous system. Studies in man on the interrelation between baroreceptor mechanisms and exercise. 19, 400–411.5914852SavonenK. P.KiviniemiV.LaukkanenJ. A.LakkaT. A.RauramaaT. H.SalonenJ. T.RauramaaR. (2008). Chronotropic incompetence and mortality in middle-aged men with known or suspected coronary heart disease. 29, 1896–1902.10.1093/eurheartj/ehn26918556711SavonenK. P.LakkaT. A.LaukkanenJ. A.HalonenP. M.RauramaaT. H.SalonenJ. T.RauramaaR. (2006). Heart rate response during exercise test and cardiovascular mortality in middle-aged men. 27, 582–588.10.1093/eurheartj/ehi70816399774SiscovickD. S.WeissN. S.FletcherR. H.LaskyT. (1984). The incidence of primary cardiac arrest during vigorous exercise. 311, 874–877.10.1056/NEJM1984100431114026472399SiscovickD. S.WeissN. S.HallstromA. P.InuiT. S.PetersonD. R. (1982). Physical activity and primary cardiac arrest. 248, 3113–3117.10.1001/jama.248.23.31137143686SmithL. L.KukielkaM.BillmanG. E. (2005). Heart rate recovery after exercise: a predictor of ventricular fibrillation susceptibility after myocardial infarction. 288, H1763–H1769.10.1152/ajpheart.00919.200315563535Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. (1996). Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. 17, 354–381.8737210TulppoM. P.HughsonR. L.MakikallioT. H.AiraksinenK. E.SeppanenT.HuikuriH. V. (2001). Effects of exercise and passive head-up tilt on fractal and complexity properties of heart rate dynamics. 280, H1081–H1087.11179050TulppoM. P.KiviniemiA.HautalaA. J.KallioM.SeppänenT.TiinanenS.MäkikallioT. H.HuikuriH. (2011). Sympatho-vagal interaction in the recovery phase of exercise. 31, 272–281.TulppoM. P.KiviniemiA. M.HautalaA. J.KallioM.SeppanenT.MakikallioT. H.HuikuriH. V. (2005). Physiological background of the loss of fractal heart rate dynamics. 112, 314–319.10.1161/CIRCULATIONAHA.104.52371216009791TulppoM. P.MakikallioT. H.SeppanenT.AiraksinenJ. K.HuikuriH. V. (1998a). Heart rate dynamics during accentuated sympathovagal interaction. 274, H810–H816.TulppoM. P.MakikallioT. H.SeppanenT.LaukkanenR. T.HuikuriH. V. (1998b). Vagal modulation of heart rate during exercise: effects of age and physical fitness. 274, H424–H429.TulppoM. P.MakikallioT. H.TakalaT. E.SeppanenT.HuikuriH. V. (1996). Quantitative beat-to-beat analysis of heart rate dynamics during exercise. 271, H244–H252.8760181von KlotS.MittlemanM. A.DockeryD. W.HeierM.MeisingerC.HormannA.WichmannH. E.PetersA. (2008). Intensity of physical exertion and triggering of myocardial infarction: a case-crossover study. 29, 1881–1888.10.1093/eurheartj/ehn23518534976WHO. (1999). Definition, diagnosis and classification of diabetes mellitus and its complications, part 1: diagnosis and classification of diabetes mellitus. Geneva: WHO.WykretowiczA.GuzikP.BartkowiakG.KrauzeT.KasinowskiR.DziarmagaM.WesselingK. H.WysockiH. (2005). Endothelial function and baroreflex sensitivity according to the oral glucose tolerance test in patients with coronary artery disease and normal fasting glucose levels. 109, 397–403.10.1042/CS2005009515948715YamamotoY.HughsonR. L. (1991). Coarse-graining spectral analysis: new method for studying heart rate variability. 71, 1143–1150.1757311YamamotoY.HughsonR. L.PetersonJ. C. (1991). Autonomic control of heart rate during exercise studied by heart rate variability spectral analysis. 71, 1136–1142.1757310