Heart rate variability is related to self-reported physical activity in a healthy adolescent population.

This study investigated whether there is a relationship between heart rate variability (HRV) versus lifestyle and risk factors for cardiovascular disease in a population of healthy adolescents. HRV is as an index of tonic autonomic activity and in adults HRV is related to lifestyle and risk factors for cardiovascular disease, but it is not known if this is the case in adolescents. HRV was registered for 4 min in sitting position in 99 healthy adolescents (age range 15 years 11 months-17 years 7 months) and repeated after 6 months. On both occasions there were significant correlations (P < 0.05) between physical activity and HRV, with respective r values: high frequency (HF) 0.26, 0.30; low frequency power (LF) 0.35, 0.29 and the standard deviation of inter-beat intervals (SDNN) 0.28, 0.37. There was no significant interaction between first and second measurements. In contrast, there were no correlations to sleeping patterns, eating habits and smoking. Risk factors for cardiovascular disease [body mass index (BMI = weight (kg)/length in m(2)), systolic blood pressure and p-glucose] did not show any repeatable significant correlations to HRV. Multiple regression models showed that physical activity was a predictor for HF, LF and SDNN in both measurements. In conclusion HF, LF and SDNN were reproducible after 6 months and were related to physical activity on both occasions.

Introduction

Heart rate variability (HRV) is related to lifestyle in adults. Regular physical exercise in adults has beneficial effect on autonomic control of the heart, with decreased resting heart rate often accompanied with increase in HRV (Hottenrott et al. 2006; Tuomainen et al. 2005). Disturbed sleeping patterns or lack of sleep effect HRV. Shift workers and patients with sleep apnoea have deranged HRV (Furlan et al. 2000; Gula et al. 2003). Smoking impairs the sympatho-vagal balance and decreases the HRV in healthy adults (Alyan et al. 2008; Karakaya et al. 2007) and in heavy smokers vagal modulation of the heart is blunted (Barutcu et al. 2005).

In adolescent HRV decreases with augmented BMI (Rabbia et al. 2003) as it does with blood glucose levels (Singh et al. 2000) and risk factor influence on autonomic heart rate control is present before manifest cardiovascular disease (Jensen-Urstad et al. 1998). The relationship between HRV and lifestyle is not well studied in children and adolescents and it is not known at what age the relationship is established.

There is a modification of HRV with age, which reflects the evolution of the autonomic nervous system (Silvetti et al. 2001). LF, HF and total power increase from 0–6 years, followed by a decrease to adulthood (Finley and Nugent 1995). Early life individual differences in HRV can be explained by hereditary factors (Singh et al. 2002), diabetic neuropathy (Colhoun et al. 2001) or an underdeveloped ANS in prematurely born children (Longin et al. 2006).

Decreased HRV is associated with risk for cardiac events in adults (Tsuji et al. 1996) and is a predictor of imminent hypertension (Schroeder et al. 2003; Singh et al. 1998). A decreased HRV is related to an increased risk of death and HRV could have a predictive value for life expectancy and health (Tsuji et al. 1996).

Our hypothesis is that HRV changes due to lifestyle factors precedes or parallels the development of cardiovascular disorders, i.e. it should present also in healthy adolescents, where long-term effects of poor lifestyle has not yet generated manifest cardiovascular disease. To this end we investigated whether there is an impact of lifestyle (physical activity, eating habits, sleeping pattern and smoking) on autonomic heart rate control in healthy adolescents, i.e. at an age when lifestyle-associated cardiovascular disorders are not yet frequent. The impacts of cardiovascular risk factors (systolic blood pressure, p-glucose and BMI) on HRV were also studied.

Methods

High school students (66 girls and 33 boys) with mean age 16.5 years at the first measurement (age range 15 years 11 months–17 years 7 months) constituted the original sample. The subjects were recruited from high schools: in a small rural town, in the Stockholm city, in a northern affluent suburb and in a less affluent southern suburb with a large immigrant population. The first author met the students and gave oral and written information about the study and an informed consent form was signed by the student and at least one parent. The measurements had the same order in all subjects. The measurements were repeated after 6 months to evaluate the stability of the intra-individual HRV. When both measurements were completed, the students were offered two cinema tickets each. About 80% of the students informed participated and no one was denied participation. The participation ratio was similar for all schools. The main reasons to decline participation were fear of blood sampling and not wanting to miss school hours. Exclusion criteria were diabetes, thyroid dysfunction and more than 5% missing or distorted data in any registered HRV segment. The study was approved by the Regional Ethics Committee at Karolinska Institutet.

The measurements were carried out at the school nurses offices by four nurses and one physician. The subjects were sitting upright, in silence and no body movements were allowed during the procedure. None of the subjects had clinical signs or symptoms of infectious disease. Use of tobacco or intake of tea, coffee, caffeine containing soft drinks or beta stimulant asthma medication was not allowed 1 h prior to the measurements.

The HRV registration was preceded by 15 min of rest. HRV was measured during 2 min × 2, in between which blood pressure was checked. This was a modified version of a 12-min protocol (von Scheele et al. 2005).

SDNN was used as a time domain measure and HF and LF as frequency domain measures. In spectral analyses variability distributes as a function of frequency (Camm 1996). HF (0.15–0.4 Hz) is related to vagal activity and includes the respiratory sinus arrhythmia when the breathing rate is normal. LF (0.04–0.15 Hz) is a quantitative marker for sympathetic modulations, but in some studies it is viewed as reflecting both sympathetic and vagal activity (Berntson et al. 1997). The equipment was a I-330-C-2 Physiological Monitoring System (J&J Engineering; Poulsbo, WA), and c-Stress customized software (PBM Systems, Stockholm, Sweden). The electrocardiogram (ECG) was recorded from electrodes placed on the left and right wrist with a sampling rate of 1,024 Hz. Inter-beat intervals were calculated online using an R-wave peak detection algorithm and stored on a PC for off-line processing. Fourier analysis was performed on 2-min segments of detrended data passed through a Hamming window. The HRV mean values were logarithmically transformed. Inter-beat intervals were scanned manually for ectopic beats and artefacts, which were replaced with cubic spline interpolation.

Lifestyle was self-reported in a form with five-point scales. Sleeping patterns were described by frequency of going to bed later than midnight (“never”, “seldom”, “1 time/week”, “2 times/week”, “>2 times per week”), food habits by the frequency of missing breakfast and/or lunch (“never”, “seldom”, “1-2 days/week”, “3–4 days/week” or “every day”), physical exercise by the frequency of exercising with hard breathing and sweating (“never”, “seldom”, “once a week”, “twice a week” and “more than twice a week”). Smoking was defined as daily smoking or not. Socio-demographic data was represented by school, number of parents in the household, parent unemployment and parent ethnicity.

Capillary p-glucose and b-haemoglobin were analysed with a portable Heamocue Glucose and Haemoglobin System (Banauch et al. 1975; Vanzetti 1966). Weight and length were measured and BMI calculated. The same equipment was used throughout the whole study for measuring weight and blood pressure.

Statistics

Results are given as means and standard deviations (SD). Non-normally distributed values were logarithmically transformed before analyses. Two-tailed group analyses were done with t test, but when normal distributions were not present it was performed with Mann–Whitney’s U test. Bi- and multi-variate relations among the variables were assessed by Pearson’s product–moment correlation or, in case of variables of an ordinal nature, by Spearman’s rank correlation. Partial correlations were used to remove the effect of confounding variables. Stability of repeated measures was calculated with paired t test or, in case of ordinal variables, the Wilcoxon matched pairs test was used. Possible differences in correlation structure between the two test periods were checked with structural equation modelling. Multivariate relations were further explored by principal component analysis in order to get a stable foundation for the construction of regression models. Such regression models were then further refined by using stepwise models or ‘best subset’ models with Mallow’s Cp as the criterion for best fit. Probability levels of 0.05 or less were considered significant. All analyses were done in Statistica 8.0 (http://www.statsoft.com).

Results

Study population

One subject with diabetes and one with thyroid dysfunction were excluded. In the first measurements, 20 subjects were excluded because of poor quality of registered HRV data (>5% of missing or distorted beats). The final sample included 77 subjects for the first occasion. In the repeated measurements, four subjects failed to appear (change of school or sick leave) and one handed in incomplete data. One exclusion was done due to poor quality of HRV registration. 71 subjects had full data from both measurements, and constituted the basis for all correlations that are reported. In the original sample 68.4% were girls. 64.4% lived with both of their parents, 75.0% had both parents employed and 67.1% had both parents born in Sweden. The excluded subjects did not differ significantly from the original sample, but the exclusions were unevenly distributed among the schools (P < 0.1) (Table 1).

Table 1

Comparison between included and excluded subjects on background variables in the original sample

	Included	Excluded	χ2	P
n	71	28
Sex			0.02	0.87
Males	24	9
Females	47	19
School			10.0	0.02
1. Rural	13	10
2. Northern affluent suburb	18	4
3. City centre	18	1
4. Southern suburb	22	13
Family situation			0.78	0.38
Live with two parents	52	18
Live with single parent	19	10
Parents occupational status			0.63	0.43
Both parents employed	53	23
At least one parent unemployed	18	5
Parents ethnical background			7.63	0.01
Both parents of foreign origin	10	11
At least one parent Swedish	61	17

Repeated measures

Significant correlations (P < 0.05) were found between the first and repeated measurement, with r values for HF (0.54), LF (0.62) and SDNN (0.64). The correlations between measurements increased slightly when controlled for heart rate, with r values for HF (0.66), LF (0.65), SDNN (0.66). Structural equation modelling did not identify any differences in correlation structure between the two test periods. The correlations between the first and repeated measurement were also significant for BMI (0.94), systolic blood pressure (0.69) and p-glucose (0.32). No significant differences were found in socio-demographic and lifestyle parameters between the two measurements (Table 2).

Table 2

Comparison between the first and the repeated measurement after 6 months with paired t test (N = 71)

	First measurement mean (SD)a	Second measurement mean (SD)	t value	P
Log HFb	5.88 (0.93)	5.90 (1.04)	−0.13	0.89
Log LFc	6.03 (0.90)	6.20 (1.04)	−1.72	0.09
SDNNc	4.11 (0.35)	4.18 (0.38)	−1.86	0.07
Heart rate	74.5 (9.63)	71.43 (10.65)	2.62	0.01
p-glucose	5.6 (0.76)	5.4 (0.71)	2.57	0.01
BMId	22.3 (3.39)	22.7 (3.43)	−2.31	0.02
Systolic blood pressure	115 (10.4)	111 (9.77)	3.99	0.001
Single parent status	28.8%	27.8%
Both parents of foreign origin	15.2%	15.2%
At least one parent unemployed	26.0%	20.8%
Physical activity	3.93 (1.17)	3.86 (1.18)	0.63	0.53
Irregular eating habits	2.14 (1.24)	2.02 (1.20)	1.42	0.15
Late sleeping pattern	2.64 (1.52)	2.54 (1.44)	0.94	0.34

aStandard deviation, b Logarithmically transformed high frequency power of HRV, c Logarithmically transformed low frequency power of HRV, d Logarithmically transformed standard deviation of inter-beat intervals, e Body mass index

Gender

In the final sample 64.4% were girls. Boys had a higher frequency of physical activity and girls a higher frequency of skipping meals. No gender differences were found in sleeping patterns. Only eight subjects were daily smokers of which six were girls. BMI and p-glucose showed no gender differences, whereas systolic blood pressure was higher in boys, as expected from gender-specific reference values. Girls had a significantly higher heart rate compared to boys but no significant differences were found in the HF, LF or SDNN. Further analyses were therefore presented as one group (Table 3).

Table 3

Sex differences in the HRV measures, heart rate, clinical characteristics and blood chemistry

	Mean (SD)a/nb (%) girls	Mean (SD)/n (%) boys	t value/χ2	P
n	47	24
Living with single parent	12 (25.5%)	7 (29.2%)	0.11	0.74
At least one parent unemployed	12 (25.5%)	6 (25.0%)	0.00	0.96
Both parents of foreign origin	8 (17.0%)	2 (8.3%)	0.99	0.32
Number of daily smokers	6 (12.8%)	1 (4.212.5%)	1.33	0.25
Physical activity (frequency index)	3.5 (1.14)	4.7 (0.82)	−4.44	<0.001
Irregular eating habits (frequency index)	2.4 (1.25)	1.5 (0.98)	3.03	0.003
Late sleeping habits (frequency index)	2.7 (1.49)	2.7 (1.65)	−0.07	0.94
P−glucose (mmol/L)	5.52 (0.74)	5.90 (0.79)	−1.99	0.05
Haemoglobin (g/L)	129.7 (12.25)	145.5 (10.25)	−5.32	<0.001
Body Mass Index (kg/m2)	22.4 (3.73)	22.3 (2.75)	0.03	0.98
Systolic blood pressure (mmHg)	112.1 (9.53)	121.2 (9.75)	−3.79	<0.001
Heart rate (beats per min)	76.74 (9.51)	70.30 (8.75)	2.77	0.007
Log HFc	5.94 (0.81)	5.80 (1.17)	0.60	0.55
Log LFd	5.99 (0.88)	6.14 (0.97)	−0.64	0.52
Log SDNNe	4.09 (0.32)	4.17 (0.40)	−0.90	0.37

aStandard deviation, b Number, c Logarithmically transformed high frequency power of HRV, d Logarithmically transformed low frequency power of HRV, e Logarithmically transformed standard deviation of inter-beat intervals

Relationship between HRV versus lifestyle and risk factors for cardiovascular disease

Significant correlations (P < 0.05) were found between self-reported physical activity and HF, LF, SDNN and heart rate on both the first and second measurement (Table 4).

Table 4

Correlation coefficients between lifestyle parameters and heart rate variability

	First measurement				Second measurement
	Heart rate	Log HFa	Log LFb	Log SDNNc	Heart rate	Log HF	Log LF	Log SDNN
Physical activity
r (partial r)	−0.37* (−0.35)*	0.26*(0.23)	0.35* (0.36)*	0.28*(0.28)*	−0.42*(−0.44)*	0.30*(0.32)*	0.29*(0.29)*	0.37*(0.37)*
Irregular eating habits
r (partial r)	0.14 (0.11)	−0.06 (−0.03)	−0.16(−0.07)	−0.18 (−0.11)	0.17(0.19)	−0.05(−0.11)	−0.16(−0.15)	−0.09(−0.08)
Late sleeping pattern
r (partial r)	−0.03 (0.01)	−0.00(−0.01)	−0.08(−0.03)	−0.07(−0.05)	−0.02 (−0.02)	−0.00(−0.03)	−0.19(−0.20)	−0.05(−0.05)

Correlations significant at P < 0.05 are marked with* (N = 71), partial correlations controlling for BMI, p-glucose, systolic blood pressure are in parenthesis

aLogarthmically transformed high frequency power of HRV, b Logarithmically transformed low frequency power of HRV, c Logarithmically transformed standard deviation of inter-beat intervals

The correlations between physical activity and LF remained significant when controlled for heart rate, with a correlation coefficient of 0.27 in the first and 0.23 in the second measurement. This was not the case for HF and SDNN. Late sleeping patterns, eating habits and smoking showed no significant correlations to HRV in any of the measurements. BMI, p-glucose and systolic blood pressure did not show any significant correlations to HF, LF, SDDN or heart rate on any of the measurements except p-glucose, which correlated to SDNN in the first measurement (0.24).

Multiple regression models showed that physical activity (but not sleeping pattern, eating habits or smoking) predicted HF, LF and SDNN on both measurements. Physical activity predicted LF (but not HF and SDNN) also when heart rate was controlled for in both measurements. Of the cardiovascular risk factors only p-glucose showed a weak but significant prediction power for LF in the second measurement, but BMI and systolic blood pressure did not predict HRV (Table 5).

Table 5

Regression models showing the contribution of physical activity to HRV parameters and heart rate in the first and second measurement

	First measurement		Second measurement
	Multiple R2	Adjusted R2	Multiple R2	Adjusted R2
Log HF	0.06	0.05	0.13	0.11
Log LF	0.11	0.09	0.17*	0.13*
SDNN	0.07	0.06	0.14	0.12
Heart rate	0.13	0.12	0.21	0.18

These are ‘best subset’ models where physical activity was the sole significant predicting variable except in one case

* When only physical activity and not p-glucose was allowed in the general regression model the contribution to LF was estimated by multiple R2 0.08 and adjusted R2 0.07

Discussion

The main finding of the present study was that self-reported physical activity but not eating habits, sleeping patterns or smoking predicted HF, LF and SDNN and heart rate in a healthy study population of adolescents. Of the cardiovascular risk factors only p-glucose showed a weak but significant prediction power for LF in the second measurement, but BMI and systolic blood pressure did not predict HRV.

No main effects were noted for gender except that heart rate was significantly lower in boys. SDNN showed no gender difference in the present study, but data are being reported as higher in adolescent boys compared to girls (Silvetti et al. 2001). However, the present data was in line with the Silvetti study in that boys had a lower heart rate and that vagus-associated parameters of HRV showed no gender differences.

In adults regular aerobic exercise results in reduction in heart rate, reflecting an increase in autonomic efferent activity with enhanced vagal influence on the cardiac rhythm (Hottenrott et al. 2006). This is in line with the present findings of adolescents. The effect of physical exercise on different HRV parameters is less clear-cut in previous studies. Low to moderate intensity exercise is associated with favourable changes in HRV in middle aged men (Tuomainen et al. 2005), but limited influence of physical training on HRV has also been reported (Buchheit and Gindre 2006). High intensity intermittent training in children during 7 weeks does not affect HRV (Gamelin et al. 2009). Physical activity in the present study is self-reported as a lifestyle parameter and probably reflects a longer duration of regular training than 7 weeks, which may explain the impact on HRV. Sleeping patterns are not related to HRV, which is in contrast to what has been reported in adult samples (Furlan et al. 2000). This may reflect that our sleep measure was less directed towards sleep deprivation and rather measured a shift in diurnal rhythm. The fact that smoking did not correlate to HRV may be a lack of sensitivity, as the subjects were few and had not smoked for extended periods, combined with the low reliability in self-reported tobacco use in this age group (Hennrikus et al. 2005). The lack of correlations between BMI and systolic blood pressure versus HRV are also probably due to lack of sensitivity, as few subjects suffered from obesity or hypertension in the sample. Multiple regression models in the second measurements showed a weak predictive power of p-glucose on LF. Our data support (Faulkner et al. 2003) in that HRV data allows for the early identification of subclinical cardiac autonomic changes in youth.

The inferential power of the study is limited by that the sample encompassed only of 99 subjects, that we had some loss of data and the group did not represent a fully randomised population sample. However, the repeated measures design demonstrated that the reported regressions were robust within the sample.

The present findings support that even a short registration without any standardised interventions can capture HRV differences of possible predictive and clinical value. Five-minute Holter recordings in a non-laboratory setting has been shown to be stable over 2 months and generate individually characteristic HRV data (Sinnreich et al. 1998). Our data show stable results with 4-min registrations over a period of 6 months. The registering of HRV only during 4 min (and not 24 h) has obvious advantages for screening purposes. The HRV measurement device is portable and easily managed.

In conclusion, we have reported that HRV is related to self-reported physical activity in adolescents and thereby life style factors influence HRV at an age preceding the development of manifest lifestyle-related cardiovascular disorders. Within the age range of 15 years and 11 months to 17 years and 7 months, physical activity, but not sleeping pattern, eating habits or smoking was a predictor for HF, LF and SDNN. Of the cardiovascular risk factors, only p-glucose showed a weak predictive power for LF on one occasion. Whether the impact of physical exercise on HRV in adolescents is associated with a prognostic benefit regarding overall mortality, needs to be studied further. It would be of interest to target the impact of increased physical activity in young individuals with high exposure to known risk factors for cardiovascular disease using both HRV and future morbidity as endpoints. The present findings support that an ambulatory short registration of HRV can have possible predictive and clinical value.