OBJECTIVE To examine the correlation between heart rate variability and left ventricular mass in cyclists with an athlete’s heart.
METHODS Left ventricular mass and diastolic function were determined at rest and myocardial high energy phosphates were quantified at rest and during atropine–dobutamine stress in 12 male cyclists and 10 control subjects, using magnetic resonance techniques. Ambulatory 24 hour ECG recordings were obtained, and time and frequency domain heart rate variability indices were computed.
RESULTS In the cyclists, the mean of all RR intervals between normal beats (meanNN), the SD of the RR intervals, and their coefficient of variation were significantly greater than in control subjects (p < 0.01, p < 0.01, and p < 0.05, respectively). For cyclists and control subjects, only meanNN correlated with left ventricular mass (r = 0.48, p = 0.038). The heart rate variability indices that correlated with functional or metabolic variables were: meanNN v E/A peak (the ratio of peak early and peak atrial filling rate) (r = 0.48, p = 0.039); the root mean square of successive differences in RR intervals among successive normal beats v E/A area (ratio of peak early and peak atrial filling volume) (r = 0.48, p = 0.040); percentage of successive RR intervals differing by more than 50 ms v the phosphocreatine to ATP ratio at rest (r = 0.54, p = 0.017); and the SD of the average RR intervals during all five minute periodsv the phosphocreatine to ATP ratio during stress (r = 0.60, p = 0.007).
CONCLUSIONS Highly trained cyclists have increased heart rate variability indices, reflecting increased cardiac vagal control compared with control subjects. Left ventricular mass has no major influence on heart rate variability, but heart rate variability is significantly correlated with high energy phosphate metabolism and diastolic function.
- heart rate variability
- left ventricular mass
- athlete’s heart
Statistics from Altmetric.com
The term heart rate variability refers to variations in both instantaneous heart rate and RR interval length.1 Changes in heart rate primarily result from alterations in autonomic nervous activity in the sinus node. This variable can be measured non-invasively, can be analysed in the time or frequency domain, and can provide general information about the health status of the heart and its autonomic control.1 2
The relevance and significance of heart rate variability in various populations has been well established. Attenuation of heart rate variability, indicating either increased sympathetic activity or reduced baroreflex sensitivity,3 has been shown to be strongly predictive for sudden cardiac death and total cardiovascular mortality in the general population,4 5 in patients with coronary artery disease,6-9 and in patients with mitral valve disease.10 In patients with myocardial infarction, heart rate variability is a better predictor of cardiac events and rhythm or conduction disturbances than left ventricular function and exercise ECG or ambulatory ECG variables.2 7 11Increased sympathetic activity can induce cardiac arrhythmias, whereas increased vagal tone is considered to be protective.2 12
A significant relation between reduced heart rate variability and the severity of left ventricular hypertrophy secondary to hypertension, diabetes mellitus, or aortic valve disease has recently been demonstrated,13-16 suggesting that impaired cardiac autonomic function in pathological left ventricular hypertrophy may contribute to the mechanism of sudden cardiac death. Intensive participation in sports may lead to the development of left ventricular hypertrophy, generally called the athlete’s heart.17 To date, there is continuing debate over whether intense participation in competitive sport is associated with an increased risk of sudden cardiac death. Although exercise leads to increased heart rate variability,18-24 left ventricular hypertrophy may be associated with reduced heart rate variability.13-16Accordingly, it would be appropriate to study the impact of exercise induced left ventricular hypertrophy on the modulation of heart rate variability in athletes. With this in mind, we investigated indices of heart rate variability and its physiological correlates in cyclists with exercise induced left ventricular hypertrophy, and these data were compared with those obtained from non-athletic healthy individuals.
Twelve highly trained male cyclists (mean (SD) age 41 (10) years) and 10 male control subjects (46 (9) years), matched for height and weight, volunteered for the study. The athletes cycled at least 12 000 km/year and had been engaged in competitive cycling for 22 (8) years. None of the 10 control subjects was engaged in regular physical activity or competitive sports. All individuals were free from known cardiovascular disease and they were passed as normal on physical examination. The subjects were non-smoking, normotensive, and none was taking any drugs at the time of the study. A standard 12 lead ECG was recorded at rest. All ECGs were normal in the control subjects. In six cyclists the resting ECG met the Sokolow–Lyon voltage criteria for left ventricular hypertrophy (sum of amplitudes in V1 and R wave in V5 or V6 > 3.5 mV). The study was approved by the human research committee at our institution, and all individuals gave informed consent for their participation.
MAGNETIC RESONANCE IMAGING ACQUISITION
Magnetic resonance imaging was performed using a Philips Gyroscan ACS-NT system (Philips Medical Systems, Eindhoven, the Netherlands) operating at 1.5 Tesla, and ECG gating. Magnetic resonance images were acquired with breath hold multi-shot echo planar imaging as described previously.25 Volume flow was assessed by magnetic resonance velocity mapping of flow across the mitral orifice. Velocity maps were acquired using retrospective gating as described previously.26
MAGNETIC RESONANCE IMAGING ANALYSIS
Multislice multiphase short axis image analysis was performed using the magnetic resonance analytical software system (MASS27) and a Sun IPX workstation (Sun Microsystems Computer Corporation, Palo Alto, California, USA). End diastolic and end systolic epicardial and endocardial contours of the stack of short axis image sections were traced manually. Papillary muscles were outlined and included in the left ventricular wall. Left ventricular mass was calculated as reported previously by our institution.28 Volume flow was calculated by tracing a region of interest along the borders of the mitral valve in all time frames of a velocity map series. For each time frame, instantaneous volume flow was calculated using a computer algorithm by multiplying spatial average flow velocity and the area of the region of interest. Summation of all instantaneous volume flow data yielded total flow per cardiac cycle. The following indices of diastolic function were measured or derived: peak early filling rate, peak early filling volume, peak atrial filling rate, peak atrial filling volume, ratio of peak early and peak atrial filling rate (E/A peak), ratio of peak early and peak atrial filling volume (E/A area).29
31P MAGNETIC RESONANCE SPECTROSCOPY
31P Magnetic resonance spectra of the anterior wall of the left ventricle were acquired at rest and during atropine–dobutamine stress. A 1.5 Tesla Gyroscan S15 (Philips Medical Systems) was used, and a 10 minute three dimensional image selected in vivo spectroscopy (ISIS) protocol was performed as described previously.30 31 Magnetic resonance spectra were transferred to a remote Sun-SPARC workstation to be quantified automatically by model function analysis in the time domain, using a priori spectroscopic knowledge to improve the accuracy of the spectral parameters. The ATP level in the spectra was corrected for the relative contribution of the ATP present in the blood, and the phosphocreatine to adenosine triphosphate (PCr/ATP) ratios were determined. Signal modelling, prior knowledge, and correction for partial saturation effects were applied, according to previously published procedures.30
ATROPINE–DOBUTAMINE INFUSION PROTOCOL
After the baseline spectrum was recorded, atropine sulphate (0.03 μg/kg) was given to achieve complete cholinergic blockade. Thereafter, myocardial stress was induced by administration of incremental intravenous doses of dobutamine. Dobutamine infusion was started at a dose of 10 μg/kg/min and was increased every two minutes by 5 μg/kg/min until a steady target heart rate was reached. The maximum infusion rate allowed was 40 μg/kg/min. The target heart rate in beats/min was 85% of the predicted maximum heart rate.31 Blood pressure was recorded automatically every two minutes at rest, and each minute during dobutamine stress with an automated sphygmomanometer. Acquisition of a spectrum was then carried out during stress.
Ambulatory 24 hour ECG monitoring was performed, using a Marquette 8500 three channel Holter recorder (Marquette Electronics, Milwaukee, Wisconsin, USA). The ECGs were analysed according to the guidelines set out by the European Society of Cardiology1 using a Marquette series 8000 Holter analyser with manual override. Episodes containing beats with normal morphological characteristics were measured. The computed processed Holter ECGs were strategically reviewed by an editor and if necessary, the computer labels applied were modified. The exact point of onset of the QRS complexes was also interactively determined. In particular, episodes with the highest and lowest RR intervals were reviewed. The RR interval series was further analysed using a personal computer.32 Apart from incidental ectopic beats, all subjects were in sinus rhythm throughout the recording period. RR intervals containing ectopic beats were interpolated by dividing the sum of the intervals of the preceding and successive beat by 2. Data segments in which the Marquette Holter analyser detected dual channel noise were not analysed.
Time domain and frequency domain indices of heart rate variability were computed during six hour sleep periods at night. Sleeping hours were defined as the six consecutive hours immediately before rising. Rising was determined either from an entry in the diary or from the rise in heart rate, or, wherever possible, from both. The analysis was limited to the sleeping period in order to avoid influences of daytime heart rate variability owing to the higher activity level of the cyclists.33 The analysed part of the Holter ECG was divided into five minute episodes. Heart rate variability analysis was performed on all separate episodes. Spectral heart rate variability indices were computed by integrating the fast Fourier transform generated, interval series power spectrum in the 0.05–0.15 Hz low frequency (LF) band and in the 0.15–0.40 Hz high frequency (HF) band. Preceding the fast Fourier transform, the interval series was detrended by subtracting the linear trend obtained by linear regression, tapered to zero with a cosine function, and padded with zeros to the nearest power of two data points.
The following time domain indices were computed:
Mean of all RR intervals between normal beats (meanNN)
SD of the RR intervals (SDNN)
Coefficient of variation (CVNN) of all RR intervals
SD of the average RR intervals during all five minute periods (SDANN)
The root mean square of successive differences in RR intervals among successive normal beats (RMSSD)
Percentage of successive RR intervals differing more than 50 ms (pNN50).
The following frequency domain indices were computed34:
LF power, that is, the power between 0.05 and 0.15 Hz
HF power, that is, the power between 0.15 and 0.40 Hz
The sympathovagal balance, defined by the ratio of low frequency power to high frequency power (LF/HF).
Data are presented as mean (SD). Cyclists and control subjects were compared using Student’s t test and the Mann–Whitney test. Significance testing of frequency domain parameters was performed after log transformation of the data. The association between variables was quantified with the Pearson correlation coefficient and, where necessary, this association was adjusted for age using the partial correlation coefficient. A p value of 0.05 or less was considered statistically significant.
The anthropometric characteristics and magnetic resonance data of the 12 cyclists and 10 control subjects are presented in table 1. There were no significant differences between the following subject characteristics of the two groups: age, height, weight, body surface area, E/A peak, E/A area, and myocardial PCr/ATP concentration ratio measured at rest and during atropine–dobutamine stress (table 1). Heart rate was lower in cyclists than in control subjects (p < 0.005). Left ventricular mass and left ventricular mass index were higher in cyclists than in control subjects (both p < 0.001), which is in agreement with the higher training state of the cyclists.
HEART RATE VARIABILITY IN THE TIME DOMAIN
Time domain analysis of heart rate variability revealed relatively high values for all cyclists compared with control subjects (table 2). In cyclists, the meanNN was higher than in control subjects (p < 0.01), reflecting the lower resting heart rate in cyclists. In addition, SDNN (101 (28) v 70 (18) ms, p = 0.007) and CVNN (8.1 (2.1)% v 6.5 (1.5)%, p = 0.042) were higher in athletes than in control subjects. SDNN reflects all the cyclic components responsible for variability during the recording period, including short term and long term variations.1 By analysis of a six hour sleeping period instead of a 24 hour period, variability introduced by different activity levels and differences between day and night time heart rates were eliminated. The higher SDNN in cyclists is therefore a realistic finding, indicative of an intrinsically higher heart rate variability in athletes. This is substantiated by a significantly higher CVNN in cyclists, SD normalised by mean heart rate, which is a highly reproducible and heart rate independent variable.35
HEART RATE VARIABILITY IN THE FREQUENCY DOMAIN
There were no significant differences between cyclists and control subjects after log transformation of the data in the following frequency domain indices: LF (7.42 (0.51) v7.08 (0.54) ms2, NS); HF (6.67 (0.93)v 6.41 (0.76) ms2, NS); or the ln LF/HF ratio (0.75 (0.63) v 0.68 (0.49), NS).
RELATION BETWEEN HEART RATE VARIABILITY AND SUBJECT CHARACTERISTICS
All time and frequency domain parameters of heart rate variability were negatively and significantly correlated with age. Therefore it was necessary to perform a correction for age when calculating the correlation coefficients between heart rate variability and other subject characteristics. A significant positive relation between meanNN and left ventricular mass was found (r = 0.48, p = 0.038), reflecting the low resting heart rate of the cyclists. The other indices of heart rate variability were not correlated with left ventricular mass or left ventricular mass index (table 3).
The heart rate variability indices which correlated with left ventricular diastolic function indices were meanNNv E/A peak (r = 0.48, p = 0.039), and RMSSDv E/A peak (r = 0.42, p = 0.075) and E/A area (r = 0.48, p = 0.040).
The heart rate variability indices that correlated with myocardial high energy phosphate metabolism were SDANN vPCr/ATP ratio at rest (r = 0.44, p = 0.059) and PCr/ATP ratio during stress (r = 0.60, p = 0.007, fig 1), RMSSDv PCr/ATP ratio at rest (r = 0.41, p = 0.079), and pNN50v PCr/ATP ratio at rest (r = 0.54, p = 0.017).
Our aim in this study was to examine the correlation between heart rate variability and cardiovascular functional and metabolic indices in cyclists with training induced left ventricular hypertrophy. Our study yielded no statistically significant correlation between heart rate variability indices and the extent of left ventricular hypertrophy. This finding suggests that exercise induced increase in left ventricular mass in athletes has no influence on heart rate variability. This is in contrast with findings in left ventricular hypertrophy caused by aortic valve disease, hypertension, or diabetes mellitus, where there is a negative correlation between indices of heart rate variability and left ventricular mass.13-16Second, our study shows that meanNN and RMSSD correlated with indices of left ventricular diastolic function. Third, we found a notable correlation between the myocardial PCr/ATP ratio during stress and SDANN, a measure of ultra low frequency heart rate variability. These findings on heart rate variability underscore our previous findings that exercise induced left ventricular hypertrophy in cyclists is a physiological rather than a pathophysiological phenomenon.26 To the best of our knowledge, this is the first study to evaluate the relation between left ventricular hypertrophy and heart rate variability in highly trained cyclists.
INCREASED HEART RATE VARIABILITY IN CYCLISTS
Time domain analysis of heart rate variability revealed relatively high values for cyclists compared with control subjects. This reached statistical significance with respect to SDNN and CVNN, possibly indicating an increased vagal control in cyclists. Reduced SDNN has been shown to carry an increased risk of sudden cardiac death in a heterogeneous patient population independent of other risk factors6; it is a significant predictor of mortality in patients with myocardial infarction7 and is significantly associated with all cause mortality in the general population.4 5 Our results therefore might suggest that these cyclists have a decreased risk of sudden cardiac death.
LEFT VENTRICULAR DIASTOLIC FUNCTION
The correlations of meanNN and RMSSD with indices of left ventricular diastolic function are partly explained by a lower resting heart rate in cyclists. An increasing heart rate favours the atrial contribution versus the early phase of diastolic filling. This results in a lower ratio of the early to late peak filling rate and volume.36 37 This explains the relation between meanNN and E/A ratio, since meanNN is equivalent to heart rate. Part of the relation between RMSSD and E/A ratio can be explained by the fact that RMSSD is an index of vagal tone, and an increased vagal tone and subsequent bradycardia may be held responsible for an increased E/A ratio.38 39 The correlation between indices of heart rate variability and indices of left ventricular diastolic function is not very strong. This is to be expected, because both sets of variables are determined by a wide variety of factors. No previous studies appear to have addressed the relation between heart rate variability and left ventricular diastolic function.
MYOCARDIAL HIGH ENERGY PHOSPHATE METABOLISM
Our study showed similar PCr/ATP ratios at rest and during atropine–dobutamine stress in the hypertrophied hearts of the cyclists compared with control subjects at rest and during atropine–dobutamine stress. Second, a moderate, but significant correlation (r = 0.60, p = 0.007) was found between the myocardial PCr/ATP ratio during stress and SDANN, indicating that high energy phosphate metabolism is associated with heart rate variability. A possible explanation for this correlation might be that both a reduced SDANN and a low stress induced PCr/ATP ratio indicate a decreased energy reserve of the heart. Reduced SDANN may identify individuals at increased risk of dying. This ultra low frequency variable has been found to be a powerful predictor of a poor prognosis in patients with coronary heart disease with8 9 and without myocardial infarction,40 mitral valve disease,10 and chronic heart failure.41However, it is unclear whether the reduced heart rate variability plays a causal role or is merely a risk factor.
Myocardial high energy phosphates can be measured non-invasively in man using spatially localised 31P nuclear magnetic resonance spectroscopy, providing data on the energy status of the heart.42-44 Our study clearly shows that the relatively high myocardial PCr/ATP ratio during stress, as occurs in well trained athletes,26 contributes to prediction of high levels of heart rate variability indices.
CONFOUNDING FACTORS IN THE USE OF HEART RATE VARIABILITY
It has been shown by Bernardi et althat the amount of RR variability and its slower fluctuations largely depend on physical activity, with a very low frequency fluctuation that increases with higher activity levels.33 Highly trained elite cyclists have, by definition, higher daytime activity levels than sedentary control subjects. To compare ultra low frequency power between these two groups, the higher physical activity level of the cyclists should be corrected. This was accomplished by studying a physiologically comparable period in both groups without any physical activity—that is, sleep.
Studies of heart rate variability in athletes have yielded conflicting results. Most investigations have shown that aerobic training may increase cardiac vagal tone at rest, resulting in increased heart rate variability,18-24 but other reports have failed to show such an increase.45-49 Factors that may contribute to the inconsistency of these results are the evaluation of young individuals who already possess high vagal tone and insufficient duration of the training programme to develop left ventricular hypertrophy.47 We therefore chose to study an older population of elite cyclists with a training history of about 20 years to ensure adequate training stimulus of sufficient duration to develop increased left ventricular mass.
Only half of the ECGs of the cyclists met the Sokolow–Lyon voltage criteria for left ventricular hypertrophy. However, the sensitivity of ECG left ventricular hypertrophy is very low and we therefore compared the measured left ventricular mass. The increase in left ventricular mass found in cyclists compared with control subjects is substantial (44%) and in agreement with previous magnetic resonance and echocardiographic studies on the athlete’s heart.50
Our study shows that cyclists with training induced left ventricular hypertrophy have an increased vagal control that is independent of the extent of left ventricular hypertrophy. Heart rate variability was also significantly correlated with high energy phosphate metabolism and diastolic function. If we take into account that hearts of elite cyclists show normal diastolic function and normal high energy phosphate metabolism, we should classify training induced left ventricular hypertrophy as physiological hypertrophy, without a proven augmented risk of sudden cardiac death. The strong correlation between heart rate variability and myocardial PCr/ATP levels found in our study suggests that a low myocardial PCr/ATP ratio in particular individuals might represent a cardiovascular risk factor.
This study was supported by the Netherlands Heart Foundation (grant No 94.107), The Hague, The Netherlands. We wish to thank Janine Voogd and Dick van de Weerd for their assistance in the processing and editing of the data.