|
|
||||||||
1Merikoski Rehabilitation and Research Center, 90100 Oulu; 2Department of Medicine, Division of Cardiology, University of Oulu, 90220 Oulu; 3Polar Electro, 90440 Kempele, Finland; and 4Department of Kinesiology, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1
Submitted 14 August 2002 ; accepted in final form 20 March 2003
| ABSTRACT |
|---|
|
|
|---|
heart rate variability; cardiovascular regulation; detrended fluctuation analysis; R-wave-R-wave interval
There are no controlled studies comparing the effects of moderate- and high-volume training on HR dynamics among sedentary subjects. The present study was designed to examine the effects of aerobic training and training volume on time domain and on spectral and fractal HR dynamics over a 24-h period in a population of healthy sedentary men.
| METHODS |
|---|
|
|
|---|
|
Protocol. Subjects were not allowed to eat or drink coffee for 2 h before the exercise test, and all physical exercise and alcohol were forbidden for 48 h before the day of testing. A 24-h ambulatory R-R interval recording was started 25 h before the exercise test. Before the exercise test, subjects remained in a supine position in a quiet room for 10 min, after which blood pressure was measured in a supine position by manual auscultation. Fat percentage was calculated from four skinfolds (sites: biceps, triceps, subscapular, and suprailiac) according to the equation described by Durnin and Womersley (11). All tests were performed at the same time of day before and after the training period.
Measurement of peak O2 consumption. Subjects performed a graded maximal exercise test on a treadmill (Telineyhtymä, Kotka, Finland), starting at 4.5 km/h followed by a work rate increase at a rate of 0.5 km/h every 2 min until voluntary exhaustion. Minute ventilation and gas exchange (M909 ergospirometer, Medikro, Kuopio, Finland) were monitored continuously during the ramp protocol. Minute ventilation and gas exchange were calculated on a breath-by-breath basis and reported as mean values for 1 min. The highest value of O2 consumption measured during the test was used as the peak O2 consumption.
Measurement of R-R intervals. The R-R intervals were recorded over 24 h with Polar R-R Recorder (Polar Electro, Kempele, Finland) at an accuracy of 1 ms (41) and saved in a computer for further analysis of HR variability with Hearts software (Heart Signal, Kempele, Finland). All R-R intervals were edited by visual inspection based on electrocardiogram (ECG) portions, to exclude all the undesirable beats, which accounted for <2% in every subject. The measures of R-R interval dynamics were calculated from the entire 24-h recording and also separately for the hours representing nighttime (midnight to 6 AM) and daytime (9 AM to 6 PM) hours to detect possible diurnal differences (21, 40). Subjects were asked to go to bed before midnight and to stay in bed until 6 AM on the R-R interval recording days. The R-R intervals were recorded during a nonexercise day before and after training intervention. At the end of training intervention, the R-R intervals were recorded after a 48-h nonexercise period.
Time and frequency domain analysis of HR variability. The mean length of R-R intervals and the standard deviation of all R-R intervals were used as time-domain measures of HR variability. An autoregressive model (order 20) was used to estimate the power spectrum densities of R-R interval variability. Ultralow-frequency power (<0.0033 Hz) and very low-frequency (LF) power (0.00330.04 Hz) were calculated from the entire 24-h segment. LF and high-frequency (HF) power values were calculated from the segments of 512 R-R intervals over 24-h recording (48). LF and HF were also calculated from 1-h segments of the 24-h recording (by using segments of 512 R-R intervals), and the mean values of these segments were used to detect day and night HR variability values (40). The spectral values are expressed as absolute values (ln; ms2/Hz) and in normalized units (nu), which were obtained by dividing the power of each component by total variance from which the very LF component had been subtracted and multiplying this value by 100 (34). The CCV% were calculated as in the following equation: CCV(%) = 100 · (power of component)1/2/(mean R-R interval) (16). The frequency at the HF power peak was analyzed for every hour (from the segments of 512 R-R intervals), and the mean values were calculated separately for day and night to evaluate the effects of training on average breathing frequency.
Fractal analysis of R-R intervals. DFA quantifies the fractal correlation properties of the R-R interval data (38). The root-mean-square fluctuations of the integrated and detrended data are measured in observation windows of different sizes and then plotted against the size of the window on a log-log scale (Fig. 1). The scaling exponent
represents the slope of this line, which relates (log)fluctuation to (log)window size. In this study, the short-term (from 4 to 11 beats) scaling exponent (
1) was used on the basis of previous experiments (32). The scaling exponent was calculated from 1-h segments of the 24-h recording, and the mean value of these segments was used (40). The power-law relationship of R-R interval variability (
) was also calculated from the frequency range of 10-4 to 10-2 by a previously described method to assess the longer-term correlation of HR (4).
|
Training program. The training period was 8 wk, including six 30-min sessions a week for the moderate-volume training group and six 60-min sessions a week for the high-volume training group at an intensity of 7080% of maximum HR. The training mode was walking and jogging. All subjects wore a telemetric HR monitor (Polar Smart Edge, Polar Electro) during all training sessions to reach and stay at the correct training intensity. Subjects also recorded the duration and average training HR of each session. An extra resting day once a week was allowed if they felt exhausted.
Statistical methods. Results are expressed as means ± SD. The normal Gaussian distribution of the data was verified by the Kolmogorov-Smirnov goodness-of-fit test. The spectral values of 24-h HR variability were skewed. Therefore, these data were transformed by taking the natural logarithms of the absolute values. The differences within the groups after training were analyzed by a two-factor ANOVA with time and interventions followed by post hoc analysis (Student's paired t-test). Pearson's correlation analysis was performed on changes in the different HR variability parameters from baseline to the end of intervention. A P value of <0.05 was considered significant.
| RESULTS |
|---|
|
|
|---|
Time and frequency domain analysis over a 24-h period. Representative examples of the effects of training on 24-h spectral HR variability are shown in Fig. 1B, and all mean values are shown in Tables 2 and 3. Mean HR decreased significantly for the moderate- and high-volume groups (P < 0.001; Fig. 2A). The mean value of the standard deviation of all R-R intervals increased for both training groups (P < 0.001). The absolute spectral values of HF, LF, and very LF power showed significant (P < 0.01) and similar increases for both training groups after the interventions. The mean value of ultralow-frequency power did not change after the training.
|
|
|
HF power normalized by the average R-R interval (CCV%) increased in both training groups after the interventions (P < 0.01; Fig. 2B), but LF power (CCV%) only increased significantly in the moderate-volume training group (P < 0.05). The LF/HF ratio decreased in both training groups (P < 0.05; Fig. 2C). The mean value of normalized HF (HFnu) power increased in both training groups (P < 0.05). Similarly, the normalized LF (LFnu) power decreased significantly (P < 0.05).
Scaling analysis of R-R intervals over a 24-h period. A representative example of the effects of training on 24-h short-term fractal HR behavior is shown in Fig. 1C, and all the mean values are shown in Table 3. The mean value of the short-term scaling exponent (
1) decreased for both training groups (P < 0.001; Fig. 2D). The mean value of
, the power law scaling exponent, did not change during the intervention.
Effects of training on the HF spectral peak. The frequency at the HF spectral peak did not change after the intervention (moderate-volume group: from 0.25 ± 0.02 to 0.25 ± 0.02 Hz at night and from 0.32 ± 0.03 to 0.31 ± 0.03 Hz during daytime; high-volume group: from 0.23 ± 0.02 to 0.23 ± 0.03 Hz at night and from 0.30 ± 0.03 to 0.30 ± 0.03 Hz during daytime; P = not significant).
Correlation between the
1 and normalized spectral components. The change in the scaling exponent
1 from the baseline to the end of intervention correlated moderately with the changes in HFnu (r = -0.76, P
< 0.001), LFnu (r = 0.76, P < 0.001), and the LF/HF ratio (r = 0.69, P < 0.001) as analyzed over the 24-h recording. The correlation between the change in
1
and the change in the LF/HF ratio became weaker during the daytime hours (r = 0.57, P < 0.001) compared with the nighttime correlation (r = 0.68, P < 0.001).
| DISCUSSION |
|---|
|
|
|---|
Methodological considerations for HR variability analysis. A decreased mean HR is normally associated with an increased total variance of R-R interval fluctuation and increases the absolute HF and LF spectral components (16). Therefore, the time domain HR variability indexes and the absolute spectral components may not be ideal methods for detecting changes in autonomic regulation after physical training interventions, where a decrease in basal HR occurs. In the present study, too, the time domain and absolute spectral calculations of HR variability are influenced by decreased HR because of the long-term aerobic training. The absolute spectral components are not able to indicate which physiological mechanisms, the enhancement in vagal tone and/or the decreased sympathetic activity, are responsible for the altered HR after training.
The normalized spectral values and LH/HF ratio represent the relative changes in the characteristics of the HR fluctuation pattern rather than the magnitude of HR variability (34). Therefore, the nu were used in the present study to assess the effects of training on the autonomic regulation of HR. Spectral components normalized by the average R-R interval length (CCV%) have also been used to evaluate the altered autonomic regulation in various pharmacological and clinical settings, where differences in the HR levels disturb the comparison of the absolute spectral values (1618). Specifically, the HF power (CCV%) showed a good linear correlation with pharmacologically assessed cardiac vagal tone (16).
The DFA technique is a method to detect qualitative rather than quantitative changes in HR dynamics. From the mathematical point of view, spectral measures are associated with fractal scaling measures when analyzed as nu during strictly controlled external conditions, because both describe relative changes in the characteristics of HR fluctuations rather than the magnitude of HR variability, and, therefore, these indexes are not related to changes in the average HR. Indeed, our laboratory has recently reported (49, 52) a relationship between
1 values and the LF/HF ratio in controlled external situations with a fixed respiratory rate. However, the relationship is weak during "free-running" ambulatory ECG recordings (31, 40, 49) because measurement of scaling exponents by the DFA method provide information on the scaling properties of HR fluctuations over several segmented time windows, whereas conventionally computed spectral ratios describe HR fluctuations only in two predetermined time windows. Therefore, both normalized spectral units of HR variability and fractal scaling exponent were used in this study to describe changes in HR dynamics caused by exercise training.
From the physiological point of view, we have described the altered short-term fractal HR behavior in various pharmacological and physiological interventions (49, 52). Briefly, a change in the R-R interval dynamics from fractal behavior (
1 of
1.0) toward a stronger correlation in short-term R-R interval dynamics (
1 of
1.5) was observed during a withdrawal of vagal activity by graded atropine infusion and dynamic low-intensity exercise (49, 52). Similarly, in the present study, the reduced
1 values correspond to increased vagal activity of the heart after aerobic training.
Effects of training on HR dynamics. Bradycardia is a well-known consequence of aerobic training and has been attributed to changes in the autonomic nervous system, i.e., either an increase in parasympathetic activity or a decrease in sympathetic activity, or else to a reduced intrinsic HR. Good aerobic fitness has been suggested to be associated with increased HR variability, especially with vagally mediated respiratory sinus arrhythmia, which is expressed as various time domain indexes or as absolute spectral HF power (810, 13, 14). However, some studies have failed to show any association between good aerobic fitness and enhanced HR variability (7, 26). The results of prospective, well-controlled interventions are also controversial. All of the studies have indicated decreased HR, and most have revealed increased vagal activity measured by absolute HF power after aerobic training (1, 28, 4244, 46). Contrariwise, some controlled studies have failed to show any association between aerobic training and HR variability (5, 29, 30). The results of these previous contradictory studies are partly limited by either the small sample size (26, 30) and subject characteristics, e.g., a narrow range in aerobic fitness (7). Also, the duration of interventions (from 6 to 36 wk), training frequency (from 3 to 7 sessions/wk), and the intensity of exercise (from
60 to 90% of maximal HR) varies between the previous studies, and thus the results of different studies are difficult to compare. However, long-duration and/or intense training may not necessarily lead to greater enhancement in HR variability as prolonged (12 mo) and intense training (from walking to running) may return these changes in HR variability back to the control level (22). The training sessions performed during the previous days also affect HR variability and autonomic function (2, 12, 15). In our laboratory's recent study (15), the autonomic function recovered back to the baseline level at the second day after strenuous aerobic exercise. Therefore, 2 days of nonexercise were used in the present study before R-R intervals were recorded.
Several other methodological issues may also explain the controversial results. First, in our population-based cross-sectional study, good aerobic fitness was associated with an augmentation of vagally mediated beat-to-beat R-R interval fluctuation during a controlled exercise but not at controlled rest (51). This suggests that short-term HR variability recording at rest (
10 min) is not an ideal method for evaluating changes in autonomic regulation. Although the spectral HR variability methods with short-term recording at rest may not be sensitive enough to detect subtle changes in autonomic regulation after training, standard spectral methods over 24 h also failed to show any association between aerobic training and spectral HR variability indexes (26, 29). Second, sampling frequency is important for an accurate detection of beat-to-beat fluctuations in R-R intervals, as the precision of sampling of 8 vs. 1 ms resulted in a significant difference in the values of HR variability indexes, e.g., the difference was 14% for the LF/HF ratio in a 24-h recording (47). In the present study, the R-R intervals were recorded with a real-time microprocessor QRS detector system with 1-ms timing accuracy (41).
In the present study, the increased HFnu power and the decreased LFnu power as well as the decreased LF/HF ratio showed an alteration in the autonomic regulation of HR toward increased vagal dominance after aerobic training. The LF/HF ratio decreased also during the night hours after moderate-volume training intervention. This is an important finding because the night hours may reflect a more standardized condition and a subject's behavior pattern does not disturb the R-R interval recordings. Concurrently, the present data also showed a decreased short-term scaling exponent
1 in R-R interval dynamics during the free-running daytime hours as well as during the night hours after aerobic training. The normalized spectral components remained unaltered during the night hours after a high-volume training period, suggesting that the scaling indexes detect subtle alterations in HR more effectively than simple spectral ratios.
In the previous pharmacological studies, both HFnu power indexes (CCV% and nu) are associated with vagal activity (16, 33). Also in the present study, both HFnu power indexes increased similarly as evidence of increased vagal activity after aerobic training. The opposite changes occur in the LF power normalized by average R-R interval (CCV%) and LFnu after aerobic training. These changes are concurrent with pharmacological studies where LF (CCV%) power decreased significantly during atropine infusion, revealing that this index is partly related to vagal activity (16). On the contrary, LFnu is related to sympathetic activity studied by muscle sympathetic activity recordings from the peroneal nerve (33, 37), and decreased LFnu in the present study may therefore indicate a decreased sympathetic activity after aerobic training.
Although not a focus of this study, it was observed that both moderate- and high-volume training groups had reductions of 57 beats/min in maximal HR. This could have been a consequence of altered autonomic balance, although vagal activity is elevated only at lower HR. Alternatively, a reduction in intrinsic HR and altered sinoatrial node responses to norepinephrine have been identified as potential contributors to reduced maximal HR after training (19). This might suggest that altered autonomic balance and reduced intrinsic HR contribute to the changes that we observed in HR after training. Another interesting finding is an increased diastolic blood pressure after intervention for the high-volume training group (Fig. 3). This may be caused by an early stage of over-training usually referred to as "overreaching" syndrome. However, specific and sensitive parameters are not available to identify overreaching syndrome, and we did not focus on that issue in the present study (25).
|
Study limitations. Endurance training produces well-known functional improvements in cardiac performance. An enlarged blood volume, together with greater ventricular compliance and distensibility (27), may influence HR variability with changes of blood pressure variability as the input signal of the baroreflex system. In conscious dogs, changes in cardiac mechanics affect blood pressure variability, especially in HF by influencing stroke volume (36). Furthermore, the tidal volume becomes larger, and breathing frequency is considerably reduced after training interventions, at least in submaximal exercise intensity levels (23, 53). The reduced breathing frequency is known to increase the absolute spectral values of HF power of R-R intervals (39). Therefore, changes in the respiratory components of HR variability, reflecting cardiac outflow, may also be due to altered respiratory patterns caused by physical training. Consequently, the HF power of the R-R interval may not be an ideal index to estimate vagal activity if the cardiac mechanical properties and breathing pattern are changed. The effects of altered cardiac mechanics and respiration on fractal scaling properties of HR fluctuation are not well known. Therefore, more studies of the physiological background, particularly of the effects of respiratory patterns, on the fractal HR dynamics are warranted.
Analysis of cardiac autonomic regulation from 24-h recordings by frequency domain as well as by other measures may also have potential biases. Most importantly, many indexes may be influenced by physical activity during recordings (3). However, several studies have shown good short-term and long-term reproducibility of 24-h recordings when measured from surface ECGs (20, 24, 35). Furthermore, 24-h recordings have appeared to be a reliable and practical method that is not dependent on a subject's cooperation and does not involve any placebo effect on the results (45, 48). This may result from a rather equal level of physical activity in individual subjects between the different 24-h periods. Therefore, 24-h HR variability measurement can be considered as an ideal method for assessing temporal changes in the autonomic regulation in individual subjects. In fact, reproducibility of 24-h HR variability measures may be even better than assessment of HR dynamics from short-term recordings during controlled situations (6, 50).
In conclusion, the reduction of HR after aerobic training is associated with increased indexes of cardiac vagal outflow in healthy sedentary males. An 8-wk aerobic training intervention, including six 30-min sessions a week at an intensity of 7080% of maximum HR, is a sufficient intervention to induce these effects on cardiac autonomic regulation.
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
This research was funded by grants from the Ministry of Education (Helsinki, Finland), the Heart and Stroke Foundation of Ontario (no. T-3081, Toronto, Canada), and the Medical Council of the Academy of Finland (Helsinki, Finland).
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. P. Sloan, P. A. Shapiro, R. E. DeMeersman, E. Bagiella, E. N. Brondolo, P. S. McKinley, I. Slavov, Y. Fang, and M. M. Myers The Effect of Aerobic Training and Cardiac Autonomic Regulation in Young Adults Am J Public Health, May 1, 2009; 99(5): 921 - 928. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Tumiati, G. Mazzoni, E. Crisafulli, B. Serri, C. Beneventi, C. M Lorenzi, G. Grazzi, F. Prato, F. Conconi, L. M Fabbri, et al. Home-centred physical fitness programme in morbidly obese individuals: a randomized controlled trial Clinical Rehabilitation, October 1, 2008; 22(10-11): 940 - 950. [Abstract] [PDF] |
||||
![]() |
K. S. Heffernan, J. J. Sosnoff, C. A. Fahs, K. K. Shinsako, S. Y. Jae, and B. Fernhall Fractal scaling properties of heart rate dynamics following resistance exercise training J Appl Physiol, July 1, 2008; 105(1): 109 - 113. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. L. T. Uusitalo, E. Vanninen, E. Levalahti, M. C. Battie, T. Videman, and J. Kaprio Role of genetic and environmental influences on heart rate variability in middle-aged men Am J Physiol Heart Circ Physiol, August 1, 2007; 293(2): H1013 - H1022. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. A. Dewland, A. S. Androne, F. A. Lee, R. J. Lampert, and S. D. Katz Effect of acetylcholinesterase inhibition with pyridostigmine on cardiac parasympathetic function in sedentary adults and trained athletes Am J Physiol Heart Circ Physiol, July 1, 2007; 293(1): H86 - H92. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Felber Dietrich, C. Schindler, J. Schwartz, J.-C. Barthelemy, J.-M. Tschopp, F. Roche, A. von Eckardstein, O. Brandli, P. Leuenberger, D. R. Gold, et al. Heart rate variability in an ageing population and its association with lifestyle and cardiovascular risk factors: results of the SAPALDIA study Europace, July 1, 2006; 8(7): 521 - 529. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Buchheit and C. Gindre Cardiac parasympathetic regulation: respective associations with cardiorespiratory fitness and training load Am J Physiol Heart Circ Physiol, July 1, 2006; 291(1): H451 - H458. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. Hautala, T. Rankinen, A. M. Kiviniemi, T. H. Makikallio, H. V. Huikuri, C. Bouchard, and M. P. Tulppo Heart rate recovery after maximal exercise is associated with acetylcholine receptor M2 (CHRM2) gene polymorphism Am J Physiol Heart Circ Physiol, July 1, 2006; 291(1): H459 - H466. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Tulppo and H. V. Huikuri Origin and significance of heart rate variability J. Am. Coll. Cardiol., June 16, 2004; 43(12): 2278 - 2280. [Full Text] [PDF] |
||||
![]() |
A. L. T. Uusitalo, T. Laitinen, S. B. Vaisanen, E. Lansimies, and R. Rauramaa Physical training and heart rate and blood pressure variability: a 5-yr randomized trial Am J Physiol Heart Circ Physiol, May 1, 2004; 286(5): H1821 - H1826. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. A. Lipsitz Physiological Complexity, Aging, and the Path to Frailty Sci. Aging Knowl. Environ., April 21, 2004; 2004(16): pe16 - pe16. [Abstract] [Full Text] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |