Association of health behaviour with heart rate variability: a population-based study

Epidemiological studies provide evidence that a healthy lifestyle (e.g., physical activity, not smoking, healthy diet, and moderate consumption of alcohol) reduces morbidity and mortality due to cardiovascular causes [1]. In addition, reduced heart rate variability (HRV), a marker of autonomic dysfunction, has been shown to be associated with an increased risk of incident myocardial infarction, cardiovascular mortality, and death from other causes in general populations [2–4], and to be associated with a poor prognosis of cardiovascular diseases (CVD) [5, 6]. Reduced HRV has been shown to be related to risk factors for cardiovascular disease [7–13], so autonomic dysfunction could be a mediator of the association of cardiovascular risk factors with CVD. There is some evidence that reduced HRV is amenable to intervention that may improve future health outcomes [14, 15].

Several studies investigated the association of physical activity and HRV, but results are inconsistent and difficult to compare due to differences in study populations, differences in the assessment of physical activity or in training regimes, differences in recording conditions of the electrocardiograms (ECG) and due to methodological problems inherent in the analysis of HRV [9, 11, 15–21]. Furthermore, only few studies have analyzed the association of HRV with other behavioural risk factors such as smoking [10, 11, 22–24], alcohol consumption [10, 11, 23, 25, 26], and diet [10, 25, 27].

The Cardiovascular Disease, Living and Ageing in Halle (CARLA) study is a large population-based study with a comprehensive and highly standardized HRV measurement protocol. The aim of the present analyses was to assess the association of behavioural risk factors with HRV as the first step of the potential pathway to CVD. We hypothesized that physical activity and a favourable dietary pattern are directly associated with HRV, smoking is associated with reduced HRV, and alcohol consumption shows a nonlinear J-shaped association with HRV.

The present analyses are based on data from the baseline examination of the prospective, population-based CARLA study. Details of the study have been described elsewhere [28, 29]. In brief, the CARLA study is a prospective cohort study of a representative sample of the inhabitants of the city of Halle (Eastern Germany) comprising 1779 (812 females, 967 males) participants aged 45-83 years at baseline. The baseline examination took place between December 2002 and January 2006.

The study protocol was approved by the Ethics Committee of the Martin Luther University Halle-Wittenberg (Halle, Germany) and conforms to the tenets of the Declaration of Helsinki. All participants were informed about the study and written consent was obtained from them.

The medical examination involved recording of: blood pressure taken while seated; heart rate (HR) while seated and lying down; circumference of the waist and hip; weight and height; a 20-min 12-lead resting ECG (CardioControl Working Station, Welch Allyn, Delft, the Netherlands); and an echocardiogram. A venous blood sample was also taken.

A standardized, computer-assisted interview was undertaken to collect information on: socio-demographic and socio-economic variables; behavioural, biomedical and psychosocial factors; medical history; and use of medication within the preceding 7 days.

Physical activity was recorded using the Baecke questionnaire, describing the typically physical activity during the previous 12 months [30]. Being physically active at sport was defined as giving an affirmative answer to the question "Do you play sports?". A sport index was calculated based on the intensity, frequency and duration of sports activity [30]. Furthermore, physical activity during leisure time (excluding sport) was summarized in the leisure-time index. The total physical activity index was calculated as the sum of the sports and leisure-time index. Moreover, for the subset of participants who were characterized as being physically active in sports, we calculated the energy expenditure during sports activities in metabolic equivalents (MET) [31]. Time spent per week on each sporting activity was multiplied by the MET value of the activity to give MET-hours per week.

Information on smoking habits involved questions on: past and current smoking status; duration of smoking; and on the quantity of tobacco products smoked per day. We used two continuous measures of smoking for the present analyses, involving the number of: (i) currently smoked tobacco products per day; (ii) packyears of tobacco products ever smoked.

Self-reported usual consumption of alcohol in g/day was calculated from the answer to the questions "How much beer (in units of 0.5 l)...", "How much wine or champagne (in units of 0.2 l)...", and "How many glasses of spirits (2 cl/glass)... do you usually drink during one week?" The underlying concentration of alcohol was 4.8 volume-% for beer, 11.0 volume-% for wine or champagne, and 33.0 volume-% for spirits [32].

The mean dietary pattern during the previous 12 months was recorded by a validated food frequency questionnaire [33]. The use of selected food groups was classified according to the recommendations of the German Society of Nutrition as used in the KORA-MONICA Studies [33]. The classifications per food category were coded in an overall summary score per subject, which ranged from 0 (worst dietary pattern) to 30 (best dietary pattern). A favourable dietary pattern is characterized by a diet rich in carbohydrates (potatoes, pasta, rice and whole-grain bread), vegetables and fruits.

The 20-min ECG was recorded after a resting period (in the supine position) of ≥20 min. Throughout the ECG, subjects were asked to breathe at 15 breaths/min (0.25 Hz) - guided by a visual metronome - to standardize the influence of the respiratory rate on spectral HRV parameters.

All electrocardiograms were processed by the Modular ECG Analysis System (MEANS) [34] to obtain the location and type of the QRS complexes of the 20-min ECG. This information was used to compute standard time and frequency domain parameters of HRV for 5-min segments of the ECG according to the current guidelines for the analysis of HRV [35]. Artefacts and ectopic beats were replaced by interpolated normal sinus beats. We used the standard deviation of normal intervals (SDNN) - a time domain parameter - and the frequency domain parameters low frequency power (LF) (0.04 to <0.15 Hz), high-frequency (HF) power (0.15-0.4 Hz), and the ratio of LF to HF (LF/HF). To calculate frequency domain parameters, tachograms of RR intervals were adjusted for linear trends, tapered and zero-padded, and a Fast Fourier transformation was employed. For further analyses, HRV derived from the first 5-min segment of the 20-min ECG which fulfilled the following quality criteria was used: <10% of abnormal beats; stationarity of the tachogram; absence of atrial fibrillation, atrial flutter, artificially paced beats and other arrhythmias.

The 5-min HRV for 1671 subjects (94% of 1777 participants with a 20-min ECG) could be used for further analyses after 106 electrocardiograms had been excluded due to atrial fibrillation or flutter (n = 45), artificial pacemaker (n = 18), >10% non-sinus beats (n = 27), other abnormal rhythms (n = 4), or technical problems during processing and analysis of electrocardiograms (n = 12). For HR analyses, we used the mean value of the second and third measurement.

All HRV variables were log-transformed before analyses because of their skewed distribution. Potential confounders were selected for adjustment based on directed acyclic graphs (DAG) taking into account prior knowledge regarding their associations with health behaviours and HRV [36]. According to the resulting DAG, only age had to be adjusted for in analyses of the association of health behaviour with HRV (Figure 1). However, especially in cross-sectional studies as our study a clear temporal differentiation of cause and effect is impossible. We cannot separate the effect of health behaviour on biomedical risk factors or diseases from the effect of disease on behavioural and biomedical risk factors. For example, we cannot distinguish whether physical inactivity has caused diabetes or vice versa. We therefore performed several sensitivity analyses. For example, we additionally calculated models with adjustment for further potential confounders such as CVD, diabetes mellitus, heart rate, hypertension, body mass index, education, health behaviour, beta blockers, ACE inhibitors, diuretic, calcium channel blockers and antiarrhythmic agents.

Since categorization of originally continuous variables does not use within-category information efficiently [37], and non-linearities in the association of HRV and health behaviours were anticipated, we fitted age-adjusted generalized additive models (GAMs) [38]. GAMs can be used to check the assumption of the linearity of the relationship between HRV and health behaviours. For each smoothing effect in a model, a χ²-test comparing the deviance between the full model (including age and the smoothed exposure variable) and the model without the exposure variable was performed.

We calculated age-adjusted geometric means (± 95% confidence interval (CI)) of HRV parameters by categories of health-behaviour variables using linear regression models. We calculated quartiles for physical-activity indices and the dietary pattern index. For the consumption of tobacco products and alcohol, zero consumption was defined as the lowest category, and the remaining subjects were categorized in tertiles. The F-test was used to test the difference in adjusted means of HRV between categories of health-behaviour variables.

To assess the influence of prevalent disease status, we carried out sensitivity analyses of the association of HRV with health behaviour in the whole population (893 males and 778 females), as well as in a "healthy" subgroup without prevalent CVD (defined as myocardial infarction (MI), self-reported coronary artery bypass graft, self-reported percutaneous transluminal coronary angioplasty, self-reported physician-diagnosed stroke, or carotid surgery), without diabetes mellitus (defined as self-reported physician-diagnosed diabetes and/or use of anti-diabetic medication) and without HRV-relevant medication (beta-blockers, angiotensin-converting enzyme (ACE) inhibitors, anti-arrhythmic drugs) (411 males and 363 females). The association of HRV with MET-hours per week was analyzed in the subgroup of participants who played sports.

Hypothesis tests were conducted at the significance level of α = 0.01 to account for the problem of multiple testing because each independent variable was tested for its association with four HRV indices and heart rate. All analyses were undertaken using SAS 9.1 (SAS Institute, Cary, NC, USA).

GAM analyses confirmed the association of SDNN with packyears of smoking in both sexes (Figure 2) and the association of HR with currently smoked tobacco products and packyears in males. Furthermore, there was a weak (but statistically significant) positive association of sport index with HR, SDNN and with LF power in males. No further statistically significant associations between health behaviour indices and HRV parameters were observed.

Analyses of the association of HRV with MET-hours per week in participants performing sport showed no significant association of SDNN, HF, and LF/HF with MET-hours per week. However, a weak (but statistically significant) association of MET-hours per week with LF was observed in females (Figure 3).

To exclude the influence of prevalent diseases and potential confounders on the association of behavioural variables and HRV, we calculated age-adjusted models in healthy subjects, and we calculated multivariate adjusted models in the whole group (adjustment for age and additionally for CVD, diabetes mellitus, HR, education, hypertension, body mass index, beta blockers, ACE inhibitors, diuretics, calcium-channel blockers and anti-arrhythmic agents). However, there was no relevant change in the association of behavioural variables with HRV (data not shown).

We aimed to provide population-based data on the association of HRV with behavioural risk factors in an elderly general population. We found only weak and inconsistent associations of HRV with physical activity, smoking, alcohol consumption and dietary pattern. To our knowledge, the CARLA study is the largest population-based study with such a comprehensive protocol to increase reliability and standardization of HRV measurements. Our results thus add to the literature by providing more reliable estimates of true HRV values. This enables study of the associations between behavioural factors and HRV with less bias due to measurement error and variability in examination conditions.

A limitation of our analyses was its cross-sectional nature, which hampers evaluation of causal inferences. HRV was assessed simultaneously with prevalent disease and health behaviour, so we cannot exclude the possibility that prevalent disease may have affected the levels of health behaviour and HRV. We addressed this limitation by conducting subgroup analyses in healthy subjects and additionally by calculating multivariate adjusted models in the whole group. However, we did not find relevant changes in the lack of associations.

Measurement of behavioural factors by questionnaire is susceptible to misclassification [39, 40]. It cannot be ruled out that the lack of association of behavioural risk factors with HRV in the present study occurred due to measurement error in self-reported behaviours (particularly with respect to physical activity and dietary pattern). However, the used questionnaires have been validated [30, 33, 41] and are frequently used.

In general, comparison of the present HRV results with the literature must be done cautiously because of methodological differences. In particular, the wide range of recording conditions of applied ECG as well as the lengths and pre-processing steps used before spectral analyses may result in considerable differences in the HRV parameters themselves between studies. Furthermore, a limitation generally true for studies of HRV is that spectral parameters are sensitive to several physiological and environmental influences as well as to the pre-processing protocol [42]. Moreover, within- and between-subject reliability of HRV measurements can be poor [43]. However, in the present study, environmental conditions during ECG recording were controlled and standardized to improve reliability.

The literature is inconsistent with respect to the effect of physical activity on HRV. Physical activity causes a resting bradycardia that is thought to be partly due to the higher vagal tone of HR, thus increasing in particular HF and total HRV power. Some studies showed a positive association of physical activity with some or all HRV measures examined [7, 9, 15, 16, 19, 21, 44], but others did not [17, 18]. A possible explanation for these conflicting results may be that, in studies showing no association of physical activity with HRV, the frequency, duration and intensity of physical activity or physical training may have been insufficient to increase vagal regulation of HR. This interpretation is supported by studies which showed that a short period of intense aerobic training as well as long-term (but low-intensity) endurance training did not cause any change in resting HRV [18, 45]. In contrast, a 6-month randomized trial of moderate-intensity training (8 kcal/kg per week) revealed increased HRV in post-menopausal sedentary women [44]. A further explanation for conflicting results might be that the effect of physical activity on HRV is modified by age. Older subjects showed smaller effects than younger ones, which might be due to a reduced "trainability" of the heart [19]. However, our results showed no modification of effects by age (data not shown).

A methodological explanation for discrepant results regarding physical activity and HRV might be that positive associations tended to be derived primarily from long-term ECG recordings. This could be due to the variability of current physical activity during long-term ECG recording, which influences HRV measurements particularly low-frequency HRV [46].

Numerous studies with different recording lengths of ECG and different populations showed a negative association of smoking with various measures of HRV [3, 10–12, 22–24], although significant effects were not always observed for all HRV measures determined [8, 11, 12, 24] or for both sexes [20], or were not confirmed in multivariable analyses [8]. Stolarz et al. found conflicting results for different populations [9], and no association was observed in the Rotterdam Study [4]. In the present study, we also did not find a clear pattern of a smoking effect on HRV, and the observed effects of smoking were mostly small. We assumed that current HRV is mostly affected by current smoking, but we did not find an association with the currently smoked number of tobacco products. In contrast, we found significantly higher HRV in males who never smoked compared with smokers with any accumulated packyears. This argues for a long-term effect of smoking on autonomic function.

Findings regarding alcohol consumption and HRV are inconsistent. Some authors reported a positive association in women [22, 26], but others found a negative association [10, 11] or no association [4, 23, 25]. Our results are in accordance with the latter findings. One problem impeding real comparability of the studies is that alcohol consumption was defined differently across studies with varying cutoff points for exposure classification. Perhaps measures of patterns of alcohol consumption (e.g., binge drinking) would be more relevant for the association of HRV and CVD risk.

Only one study found an association of reduced HRV with an "unhealthy eating pattern" - defined as frequent consumption of white bread and full-fat milk (instead of low-fat milk) and little consumption of fruits - in men [10]. Several authors evaluated the effect of unsaturated fatty acids on HRV [25, 27, 47]. Christensen et al. found a positive association of SDNN with increased consumption of fish [25]. However, a randomized trial of the effect of the consumption of industrially produced trans-fatty acids and n-3-unsaturated fatty acids on HRV did not replicate this association [27]. In the CARLA study, we did not find an association of HRV with dietary pattern and with fish consumption.

A possible explanation of our overall results could be that HRV is primarily determined by genetic factors [48–50]. A recent study found that genetic factors accounted for a major portion of the inter-individual differences in HRV, and no single behavioural determinant appeared to have a major influence on HRV [50]. Another possibility is that HRV is not a mediator of the effect of behavioural risk factors on CVD, but becomes effective after the stage of biomedical risk factors. Thus, it could just be a marker of ill health and for prospective adverse health outcomes. HRV may therefore have been found to be more strongly related to those biomedical risk factors which have a high predictive power for future CVD events (e.g., diabetes and obesity) [10, 20].

We hypothesized that HRV is associated with behavioural factors and therefore might be a mediator of the effect of behavioural CVD risk factors on CVD events, but this hypothesis was not confirmed by our results. These findings support the interpretation that there may be no true causal association of behavioural factors with HRV, which is also compatible with the inconsistent literature. However, a final conclusion as to whether HRV is an intermediate on the causal path or a marker of subclinical or impending disease cannot be drawn on the basis of cross-sectional analyses. Fortunately, the ongoing follow-up investigation of the study subjects will give us the opportunity to examine the prospective associations of behavioural risk factors with HRV and incident events, and to assess correct temporal relations.