| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Circulation. 2007;115:468-474.)
© 2007 American Heart Association, Inc.
Exercise Physiology |
From Stanford University Medical Center, Stanford University Medical School, Veterans Affairs Palo Alto Health Care System, Palo Alto, Calif, and Cardiac Science, Bothell, Wash.
Correspondence to Nicholas J. Leeper, MD, Falk Cardiovascular Research Center, 300 Pasteur Dr, Stanford, CA 94305–5406. E-mail nleeper{at}cvmed.stanford.edu
Received September 22, 2006; accepted November 15, 2006.
| Abstract |
|---|
|
|
|---|
Methods and Results— The heart rate rise at the onset of exercise was measured in 1959 patients referred for clinical treadmill testing at the Palo Alto (Calif) Veterans Affairs Medical Center from 1997 to 2004. Multivariable Cox survival analysis was performed for 197 all-cause and 74 cardiovascular deaths that accrued during a mean follow-up of 5.4±2.1 years. Decreased heart rate changes at all initial relative exercise workloads were associated with significantly increased all-cause mortality. The heart rate rise at one-third total exercise capacity, however, was the only early heart rate variable that significantly predicted both all-cause and cardiovascular risk after adjustment for confounders. Failing to reach 1 SD in the heart rate rise at one-third total exercise capacity was associated with a 28% increased all-cause mortality rate (hazard ratio, 0.72; 95% CI, 0.61 to 0.85; P<0.001) and a 35% cardiovascular mortality rate (hazard ratio, 0.65; 95% CI, 0.49 to 0.86; P=0.003). Of all heart rate measurements considered (initial and recovery), the heart rate increase at peak exercise was the most powerful predictor of cardiovascular prognosis after adjustment for potential confounders. The Duke treadmill score, however, was superior to all heart rate measurements in the prediction of cardiovascular mortality.
Conclusions— In the present study population, a rapid initial heart rate rise was associated with improved survival, but the heart rate increase at peak exercise and other conventional measurements such as exercise capacity and the Duke treadmill score were more powerful predictors of prognosis.
Key Words: exercise heart rate mortality nervous system, autonomic
| Introduction |
|---|
|
|
|---|
The exercise test provides an opportunity to examine the interaction of the autonomic nervous system and the CV system at various phases of rest, exercise, and recovery. Initial investigations considered the maximal HR response, revealing the association between chronotropic incompetence, coronary disease, and CV mortality.5,6 Later investigations evaluated the recovery phase after exercise and discovered that impaired HR recovery, which is thought to partially reflect vagal tone, predicted all-cause mortality.7 More recently, investigators have focused on the initial HR response to exercise. Using a selective pharmacological blockade approach, Nobrega et al8 found that a relative tachycardia at the first 4 seconds of exercise onset was due exclusively to vagal withdrawal and was independent of the sympathetic nervous system. The HR response to exercise beyond this initial period of exercise testing is strongly correlated with resting HR and the relative exercise intensity for a given patient.9,10 Falcone et al11 recently evaluated the HR response to a uniform absolute workload (25 W or 2 metabolic equivalents [METs] in nonobese patients) in the first minute of semisupine cycle exercise test in patients with angiographically documented coronary artery disease. The authors reported a hazard ratio of 15 for CV death for an HR rise exceeding 12 bpm in the first minute of exercise. The prognostic value of early HR measurements during treadmill testing remains to be investigated. We aimed to comprehensively evaluate the prognostic value of early HR measurements in 1959 patients referred for treadmill testing for routine clinical indications.
Editorial p 430
Clinical Perspective p 474
| Methods |
|---|
|
|
|---|
Exercise Testing
Subjects underwent symptom-limited treadmill testing using an individualized ramp treadmill protocol and exercised to maximum exertion.12 All tests began at a uniform speed of 2 mph at 0% grade. A pretest questionnaire was used to predict a target maximal MET level that would be reached within 10 minutes.13 HR targets were not used as an end point or to judge the adequacy of the test. Subjects were placed in the supine position immediately after exercise.14 No medications were changed or stopped before testing, and no test was classified as indeterminate.15 The senior authors (V.F., J.M.) read all exercise studies. Blood pressure was taken manually every 2 minutes, and exercise capacity (in METs) was calculated from treadmill speed and grade. The Duke treadmill score was calculated by converting METs estimated from the individualized ramp protocol to minutes in the Bruce protocol.
ECG Study
Twelve-lead ECG data were recorded at 500 samples per second during the exercise test. Visual ST-segment depression was measured at the J junction; ST slope was measured over the following 60 ms and classified as upsloping, horizontal, or downsloping. The ST response considered was the greatest horizontal or downsloping ST-segment depression in any lead except aVR during exercise or recovery. An abnormal response was defined as
1 mm of horizontal or downsloping ST-segment depression.
Median average HRs were computed for 15-second intervals by sorting the R-R intervals longest to shortest and averaging the middle half. Because all treadmill tests began at an identical speed and grade before increasing according to individualized protocols, multiple early HR parameters were measured to study both uniform absolute workload end points (very early time points and at fixed MET levels) and uniform relative workload end points (at peak exertion and by ratio to peak exercise). HRs were computed in the pretest standing state, at 15 seconds of exercise, at 2 METs, at 1 minute of exercise, at one-third the total exercise time, at peak exercise, and at each minute of recovery. The absolute change in HR at each time point was defined as follows:
HR15 seconds is HR at 15 seconds minus standing resting HR;
HR2METs is HR at 2 METs minus standing resting HR;
HR1 minute is HR at 1 minute minus standing resting HR;
HR1/3 exercise is HR at one-third total exercise capacity minus standing resting HR; HR increase is peak HR minus standing resting HR; and HR recovery is peak HR minus HR at 2 minutes of recovery.
Follow-Up
Cardiovascular mortality was the main outcome considered. Mortality data were gathered from the Social Security Death Index and California Death Registry, and cause of death was determined from the registry classification and confirmed by Veterans Affairs medical records. Death status was determined as of March 2006. Cardiovascular mortality was defined as death from a clearly identifiable CV cause or death of subjects with a history of CV disease and no identifiable non-CV cause for death. Classification was made with the observers blinded to the test results and resolved by consensus of 2 observers; conflicts were resolved by the senior author (V.F.).
Statistical Analysis
Clinical and exercise variables for survivors and nonsurvivors were compared through the use of
2 tests (categorical variables) and unpaired t tests (continuous variables). HR measurements were nonnormally distributed and thus were compared by use of the Mann-Whitney U test. To evaluate the relationship between HR measurements and exercise test responses, a correlation matrix was constructed using the Spearman-Rank correlation.
The relationship between tertiles of
HR1 minute and CV mortality was evaluated according to the Kaplan-Meier method. Separate multivariable Cox proportional hazards analyses adjusted for age, ß-blocker use, history of congestive heart failure, smoking history, and Duke treadmill score were used to evaluate associations between HR measurements and all-cause and CV mortality. All continuous variables were standardized and included as such in Cox survival analysis to avoid artifacts of dichotomization, and hazard ratios are reported for a 1-SD increase in each variable. Given 197 total all-cause deaths, 74 total CV deaths, and a 2-sided error rate of 5%, we estimated that the study had a statistical power of 80% to detect a 22% change in the relative risk of all-cause mortality and a 38% change in the relative risk of CV mortality per SD of each HR parameter. Discriminative accuracy was assessed through the right-sided concordance index (C index), and predictive power was evaluated using
2 analysis of deviance. The proportional hazards assumption was evaluated with the scaled Schoenfeld residual, and linearity of associations was assessed with the Martingale residual. No violations were observed.
Survival analyses were repeated in the subgroup of subjects with documented coronary artery disease (prior coronary intervention, prior myocardial infarction, abnormal coronary catheterization, typical angina, or abnormal ST-segment response to exercise). In this subgroup analysis, we estimated that the present study had a statistical power of 80% to detect a 35% change in relative risk of all-cause mortality and a 49% change in the relative risk of CV mortality per SD of each HR parameter. No adjustments were made for multiple comparisons. Subjects were not removed from observation at the time of interventions or infarction because this information was not available. Statistical analyses were performed with NCSS (NCSS, Inc, Salt Lake City, Utah) and the Design and Hmisc libraries in S-Plus 6.0 (Insightful Corp, Seattle, Wash). A 2-sided value of P<0.05 was considered statistically significant.
The authors had full access to the data and take responsibility for the integrity of the data. All authors have read and agree to the manuscript as written.
| Results |
|---|
|
|
|---|
|
At the time of testing, 222 subjects (11.3%) reported typical angina pectoris, and 515 subjects (26.3%) had atypical chest pain symptoms. There were 278 subjects who had suffered a prior myocardial infarction (14.2%), and a similar number (14.1%) had been revascularized. There were 298 diabetics (15%); 497 subjects (25.4%) were current smokers at the time of testing, and 1276 subjects (65.1%) had ever smoked. There were 465 subjects (23.7%) who were taking ß-blockers at the time of enrollment.
Subjects were followed up for a mean period of 5.4±2.1 years, during which time there were 197 total deaths and 74 CV deaths. Nonsurvivors were significantly older; had higher resting systolic blood pressures, prevalence of congestive heart failure, ß-blocker use, angina, and stroke; and were more likely to have ever smoked than survivors (P<0.05).
In the subgroup with a documented history of coronary artery disease, nonsurvivors were significantly older and had a significantly higher prevalence of congestive heart failure, angina, and ß-blocker use than survivors (P<0.05).
Exercise Test Responses
Exercise test responses according to outcome status and coronary artery disease status are described in Table 2. Nonsurvivors achieved lower estimated MET values and Borg perceived exertion scores and exhibited lower Duke treadmill scores than survivors; subjects who died of CV causes also had a higher prevalence of abnormal ST-segment deviation than survivors and subjects who died of non-CV causes (P<0.05). Similar differences in exercise responses between subjects who died of CV causes were observed in the coronary artery disease subgroup (Table 2).
|
Heart Rate Responses and Outcome
The HR responses to exercise according to selected outcome and coronary artery disease status are described in Table 3. Survivors had significantly greater
HR at 1 minute and one third of exercise, HR increase at peak exercise, and HR recovery at 2 minutes than nonsurvivors, regardless of coronary artery disease status (P<0.05). There was no significant difference in median values of
HR15 seconds or
HR2METs.
|
The Spearman-Rank correlation between HR variables and exercise responses is shown in Table 4. Early HR variables were modestly interrelated (r=–0.33 to 0.71, P<0.001) and correlated with HR increase at peak exercise (r=0.31 to 0.58, P<0.05) but did not significantly correlate with other exercise test variables considered. HR increase at peak exercise was strongly correlated with HR recovery (r=0.72, P<0.001) and peak estimated METs (r=0.63, P<0.001) and weakly correlated with perceived exertion score (r=0.16, P<0.001).
|
Results of Kaplan-Meier analysis of the relationship between
HR1 minute and CV mortality are shown in Figure 1. There was a significant association between tertiles of
HR1 minute and CV mortality in the population as a whole, with a protective effect observed for greater
HR1 minute; there was a similar nonsignificant trend in the subgroup of subjects with coronary artery disease. Results of multivariable Cox survival analysis are shown in Figures 2 and 3
. Lower
HR15 seconds and
HR1 minute were significantly associated with increased all-cause mortality (P<0.05) and showed similar nonsignificant trends toward increased CV mortality in the population as a whole. A decrease of 1 SD in
HR1/3 exercise was associated with a 28% increase in all-cause mortality (hazard ratio, 0.72; 95% CI, 0.61 to 0.85; P<0.001) and a 35% increase in CV death (hazard ratio, 0.65; 95% CI, 0.49 to 0.86; P=0.003).
|
|
|
In subjects with coronary artery disease,
HR1/3 exercise was the only early HR parameter that significantly predicted prognosis (P<0.05). Reduced
HR15 seconds and
HR1 minute, however, were associated with nonsignificant trends toward increased all-cause and CV mortality. The HR increase at the uniform absolute workload of 2 METs (
HR2METs) was not significantly associated with all-cause or CV mortality in the population as a whole or in subjects with coronary artery disease.
After adjustment for potential confounders, the HR increase at peak exercise and HR recovery at 2 minutes significantly predicted all-cause and CV mortality in both the population as a whole and the subgroup of subjects with coronary artery disease. Of all HR change variables considered, the HR increase at peak exercise was the most powerful and accurate predictor of all-cause (hazard ratio, 0.62; 95% CI, 0.52 to 0.73; P<0.001) and CV (hazard ratio, 0.68; 95% CI, 0.52 to 0.89; P=0.005) mortality in multivariable Cox survival analysis. The Duke treadmill score was a marginally more powerful and accurate predictor of CV mortality than the HR increase at peak exercise (
2, 12.7 versus 8.0; C index, 0.63 versus 0.61 for Duke treadmill score and HR increase at peak exercise, respectively). In an alternative multivariable Cox regression substituting peak estimated METs for the Duke treadmill score, peak METs performed similarly.
| Discussion |
|---|
|
|
|---|
Although we discovered several HR parameters that predict survival, we found that they all are associated with the degree of HR increase with exercise that ultimately proved to be the most powerful HR prognostic indicator. This absolute change from baseline likely reflects the level of basal vagal tone and the capacity of the heart to respond to maximal sympathetic drive after parasympathetic withdrawal. Heart rate increase was directly related to HR recovery, further implicating the role of both components of the autonomic nervous system in the slowing of HR after exercise.
Our results contrast with a recently published report that found that a rapid HR rise in the first minute of a cycle ergometry exercise protocol was associated with worse outcomes.11 Differences in exercise protocols between the 2 studies might have contributed to these divergent findings. In their study, Falcone et al11 used a nonindividualized semisupine cycle protocol and measured the HR change at 1 minute at a uniform absolute workload (25 W, roughly 2 METs in a 70-kg person). Patients with a rapid HR rise at this fixed absolute workload ultimately achieved lower total METs, and the HR change in this protocol may simply be a surrogate marker for exercise capacity. We used an individualized ramp treadmill protocol, which, in addition to its demonstrated superiority in the determination of the relationship between work rate and oxygen uptake,16 requires a similar relative workload for each subject at each time point. Because the relative workload has been shown to be strongly associated with the physiological HR response to exercise,9,10 this protocol minimizes possible differences in HR that arise solely from different relative exercise workloads. We believe that our individualized ramp protocol adequately controls for differences in exercise capacity. Nevertheless, to control for possible force function differences between these exercise protocols, we measured HR changes at uniform absolute exercise workloads, at uniform relative exercise workloads, and at uniform time intervals in the individualized protocol. We found that a greater early rise in HR in response to uniform relative exercise workloads and uniform time intervals in the individualized protocol was associated with better prognosis. We found no association between the HR rise at the uniform absolute workload of 2 METs and prognosis.
Falcone et al11 recently reported a significantly increased risk for CV mortality in subjects with coronary artery disease who had a
HR1 minute of
12 during nonindividualized semisupine cycle testing. They postulated that this early HR change was associated with an increased risk of death related to sympathetic hyperactivity or premature vagal withdrawal. Furthermore, they postulated that a blunted HR rise in exercise is a marker for persistent protective parasympathetic tone that antagonizes the deleterious effects of the adrenergic system.
On the other hand, others believe that a rapid HR rise is a marker of rapid vagal withdrawal and high resting vagal tone. In an interesting approach, Almeida et al17 studied the response to a 4-second exercise test. They defined the cardiac vagal tone index as the ratio of the R-R intervals at rest and at 4 seconds of resistance-free maximal effort stationary cycling. The cardiac vagal tone reflects vagotonia and CV fitness and importantly was found to correlate with HR recovery, another prognostic indicator thought to partially reflect vagal activity.18,19 The authors postulated that a brisk HR rise reflects the rapid removal of vagal tone, not the early application of sympathetic activity. This hypothesis is supported by pharmacological experiments in which blockade of the parasympathetic system with atropine blunts the cardiac vagal tone, whereas inhibition of the sympathetic system with propranolol fails to suppress the HR rise in this exercise interval.8 Taken together, these findings suggest that the HR profile early in exercise is dominated by the parasympathetic nervous system. The relationship between autonomic modulation, early HR changes in response to exercise, and prognosis remains to be demonstrated in humans, however.
Study Limitations
Our present study population was predominantly male and older, and these results may not be generalizable to other populations. Although our coronary artery disease subgroup had signs and symptoms of coronary artery disease, only
20% had angiographically confirmed coronary artery disease. In addition to the differences in exercise protocols described above, other methodological issues between the present study and the report by Falcone et al11 may make direct comparisons between the 2 studies difficult. Unlike the Falcone et al investigators, we did not provide a pharmacological washout before enrollment, nor did we have access to revascularization data; thus, we were unable to evaluate the predictive power of the early chronotropic parameters after adjusting for these factors. The present study was a retrospective, hypothesis-generating study designed primarily to evaluate associations between chronotropic parameters and outcome. As is the case in all exploratory retrospective cohort studies, the possibility remains that unmeasured clinical factors were unaccounted for. Prospective studies are needed to confirm these results and to elucidate the pathophysiological basis for these findings.
Clinical Implications
The present study shows that the HR changes early in exercise are weakly associated with CV mortality and do not add to conventional measurements or the Duke treadmill score. Contrary to prior evidence, it appears that a faster rise predicts improved survival and may represent overall resting vagal tone and its natural release. Several early time points for HR measurement were prognostic, but the total HR rise ultimately was the most powerful HR variable for predicting CV mortality. Future studies should continue to investigate additional markers for autonomic tone, and particular attention should be paid to the very onset of exercise. For the time being, the initial HR response to the standard exercise test should not be used for clinical decision making.
| Acknowledgments |
|---|
Dr Hadley is an employee of Cardiac Science. The other authors report no conflicts.
| References |
|---|
|
|
|---|
2. La Rovere MT, Bigger JT Jr, Marcus FI, Mortara A, Schwartz PJ. Baroreflex sensitivity and heart-rate variability in prediction of total cardiac mortality after myocardial infarction. Lancet. 1998; 351: 478–484.[CrossRef][Medline] [Order article via Infotrieve]
3. Vivekananthan DP, Blackstone EH, Pothier CE, Lauer MS. Heart rate recovery after exercise is a predictor of mortality, independent of the angiographic severity of coronary disease. J Am Coll Cardiol. 2003; 42: 831–838.
4. Shelter K, Marcus R, Froelicher VF, Vora S, Kalisetti D, Prakash M, Do D, Myers J. Heart rate recovery: validation and methodologic issues. J Am Coll Cardiol. 2001; 38: 1980–1987.
5. Ellestad MH, Wan MKC. Predictive implications of stress testing: follow-up of 2700 subjects after maximal treadmill stress testing. Circulation. 1975; 51: 363–369.
6. Lauer M, Okin P, Martin G, Evans J, Levy D. Impaired heart rate response to graded exercise: prognostic implications of chronotropic incompetence in the Framingham Heart Study. Circulation. 1996; 93: 1520–1526.
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: 552–555.
8. Nobrega AC, Castro CL, Araujo CG. Relative roles of the sympathetic and parasympathetic systems in the 4-s exercise test. Braz J Med Biol Res. 1990; 23: 1259–1262.[Medline] [Order article via Infotrieve]
9. Lewis SF, Taylor WF, Graham RM, Pettinger WA, Schutte JE, Blomqvist CG. Cardiovascular responses to exercise as functions of absolute and relative work load. J Appl Physiol. 1983; 54: 1314–1323.
10. Skinner JS, Gaskill SE, Rankinen T, Leon AS, Rao DC, Wilmore JH, Bouchard C. Heart rate versus %VO2max: age, sex, race, initial fitness, and training response: HERITAGE. Med Sci Sports Exerc. 2003; 35: 1908–1913.
11. Falcone C, Buzzi MP, Klersy C, Schwartz PJ. Rapid heart rate increase at onset of exercise predicts adverse cardiac events in subjects with coronary artery disease. Circulation. 2005; 112: 1959–1964.
12. Myers J, Buchanan N, Walsh D, Kraemer M, McAuley P, Hamilton-Wessler M, Froelicher VF. A comparison of the ramp versus standard exercise protocols. J Am Coll Cardiol. 1991; 17: 1334–1342.[Abstract]
13. Myers J, Do D, Herbert W, Ribisl P, Froelicher VF. A nomogram to predict exercise capacity from a specific activity questionnaire and clinical data. Am J Cardiol. 1994; 73: 591–596.[CrossRef][Medline] [Order article via Infotrieve]
14. Lachterman B, Lehmann KG, Abrahamson D, Froelicher VF. "Recovery only" ST-segment depression and the predictive accuracy of the exercise test. Ann Intern Med. 1990; 112: 11–16.
15. Reid M, Lachs M, Feinstein A. Use of methodological standards in diagnostic test research. JAMA. 1995; 274: 645–651.
16. Tamesis B, Stelken A, Byers S, Shaw L, Younis L, Miller DD, Chaitman BR. Comparison of the Asymptomatic Cardiac Ischemic Pilot and modified Asymptomatic Cardiac Ischemia Pilot versus Bruce and Cornell exercise protocols, Am J Cardiol. 1993; 72: 715–720.[CrossRef][Medline] [Order article via Infotrieve]
17. Almeida MB, Ricardo DR, Araujo CG. Validation of the 4-second exercise test in the orthostatic position. Ar Bras Cardiol. 2004; 83: 160–164.
18. Ricardo DR, de Almeida MB, Franklin BA, Araujo CG. Initial and final exercise heart rate transients: influence of gender, aerobic fitness, and clinical status. Chest. 2005; 127: 318–327.[CrossRef][Medline] [Order article via Infotrieve]
19. Lazzoli JK, Soares PP, da Nobrega AC, de Araujo CG. Electrocardiographic criteria for vagotonia-validation with pharmacological parasympathetic blockade in healthy subjects. Int J Cardiol. 2003; 87: 231–236.[CrossRef][Medline] [Order article via Infotrieve]
![]() |
K. P. Savonen, V. Kiviniemi, J. A. Laukkanen, T. A. Lakka, T. H. Rauramaa, J. T. Salonen, and R. Rauramaa Chronotropic incompetence and mortality in middle-aged men with known or suspected coronary heart disease Eur. Heart J., June 13, 2008; (2008) ehn269v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
K P Savonen, T A Lakka, J A Laukkanen, T H Rauramaa, J T Salonen, and R Rauramaa Workload at the heart rate of 100 beats/min and mortality in middle-aged men with known or suspected coronary heart disease Heart, April 1, 2008; 94(4): e14 - e14. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2007 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |