Skip to content

Advertisement

Journal of Occupational Medicine and Toxicology

Open Access

Associations of physical activity, fitness, and body composition with heart rate variability–based indicators of stress and recovery on workdays: a cross-sectional study

  • Tiina Teisala1Email author,
  • Sara Mutikainen1,
  • Asko Tolvanen2,
  • Mirva Rottensteiner1,
  • Tuija Leskinen3,
  • Jaakko Kaprio4, 7, 8,
  • Marjukka Kolehmainen5,
  • Heikki Rusko6 and
  • Urho M Kujala1
Journal of Occupational Medicine and Toxicology20149:16

https://doi.org/10.1186/1745-6673-9-16

Received: 21 February 2014

Accepted: 9 April 2014

Published: 18 April 2014

Abstract

Background

The purpose of this study was to investigate how physical activity (PA), cardiorespiratory fitness (CRF), and body composition are associated with heart rate variability (HRV)-based indicators of stress and recovery on workdays. Additionally, we evaluated the association of objectively measured stress with self-reported burnout symptoms.

Methods

Participants of this cross-sectional study were 81 healthy males (age range 26–40 y). Stress and recovery on workdays were measured objectively based on HRV recordings. CRF and anthropometry were assessed in laboratory conditions. The level of PA was based on a detailed PA interview (MET index [MET-h/d]) and self-reported activity class.

Results

PA, CRF, and body composition were significantly associated with levels of stress and recovery on workdays. MET index (P < 0.001), activity class (P = 0.001), and CRF (P = 0.019) were negatively associated with stress during working hours whereas body fat percentage (P = 0.005) was positively associated. Overall, 27.5% of the variance of total stress on workdays (P = 0.001) was accounted for by PA, CRF, and body composition. Body fat percentage and body mass index were negatively associated with night-time recovery whereas CRF was positively associated. Objective work stress was associated (P = 0.003) with subjective burnout symptoms.

Conclusions

PA, CRF, and body composition are associated with HRV-based stress and recovery levels, which needs to be taken into account in the measurement, prevention, and treatment of work-related stress. The HRV-based method used to determine work-related stress and recovery was associated with self-reported burnout symptoms, but more research on the clinical importance of the methodology is needed.

Keywords

Body compositionBody fat percentageBMICardiorespiratory fitnessHRVPhysical activityRecoveryWorking hoursWork stress

Background

Physical activity (PA) is one of the factors that protects against stress [1]. Exercise and psychological stress have similar acute physiological effects because both result in potent increases in cardiovascular, sympathetic, and hypothalamic pituitary adrenocortical responses and decreases in parasympathetic responses [2]. Repeated bouts of exercise lead to physiological adaptations, including decreased resting heart rate (HR) and blood pressure and increased parasympathetic activity [3]. The theory of “cross-stressor adaptation” suggests that regular PA and good fitness lead to adaptations in response to both exercise and psychological stressors [2, 4].

The assessment of heart rate variability (HRV) has gained importance as a technique to explore the function of the autonomic nervous system (ANS) and has been widely used to diagnose both psychological and physiological disorders [5]. HRV is the beat-to-beat variation in time of consecutive heartbeats expressed in normal sinus rhythm on an electrocardiogram [6, 7]. The separate rhythmic contributions from sympathetic and parasympathetic autonomic activity modulate the heart’s R wave to R wave intervals (RR intervals) of the QRS complex at distinct frequencies [8]. Both sympathetic and parasympathetic activity have been suggested to be associated with the low frequency (LF) range (0.04–0.15 Hz) because parasympathetic activity is a major contributor to the high frequency (HF) range (0.15–0.4 Hz) component in modulation of frequencies of the HR [9, 10]. Higher HRV characterizes a healthy person with efficient autonomic mechanisms and good adaptation ability while lower HRV is an indicator of abnormal and insufficient adaptation of the ANS [7, 11]. Reduced HRV is suggested to be associated with harmful events in health [12].

A greater amount of reported stress is suggested to be associated with lower HRV [13]. Also, the effects of PA, cardiorespiratory fitness (CRF), and adiposity on HRV profile have been of interest. Chronic exercise increases HRV [1416], and favorable HRV profiles are associated with higher levels of PA [17] and better CRF [18]. Additionally, obesity has been associated with altered ANS activity [19]. Some studies have found HRV profiles to be relatively poor among obese subjects [20, 21] and to improve with weight loss [22].

Work-related stress and physical inactivity are major concerns in society today. The current interest is in reducing both work-related stress and inactivity, and objective methods, such as the HRV-based approach, are needed to diagnose stress symptoms. To make recording of HRV a practical tool for field measurements, new approaches for HRV analysis are needed. For instance, age is an important determinant of HRV, which also has a large inter-individual variability [9, 23]. Firstbeat Technologies Ltd. has been developing this kind of diagnostic tool to measure and define stress and recovery based on individual values for HR and HRV (for more information, see the company’s web page [24]). Some stress and recovery variables, calculated with this tool, are suggested to be related to chronic work stress and emotions at work [25]. Using this tool, we investigated how different indicators of PA, CRF, and body composition, which are related to each other, are associated with HRV-based indicators of stress and recovery on workdays among healthy young men. The secondary aim was to evaluate the association between objectively measured stress and self-reported burnout symptoms.

Methods

Subjects

The population of this cross-sectional study consisted of groups from two separate studies, the FitFatTwin study (FFT study) and the Body & Future Health study (BFH study). The identical measurements of the cross-sectional FFT study and the baseline measurements of the BHF were pooled for this analysis. In total, 46 men (23 monozygotic twin pairs) participated in the FFT study. The BFH study included 37 men, two of whom were excluded from this study because their HRV recordings included only days off from work. Altogether, 81 healthy males (age range 26–40 y) who were not on regular medication were included in the present study. The twin study was enriched with pairs having within-pair differences in their PA habits. Because the BFH study was originally planned to investigate overweight and physically inactive individuals, it involved the following inclusion criteria: BMI 25.0–35.0 kg/m2, waist circumference ≥94 cm, no vigorous exercise (>20 min/session) more than twice a week, and no smoking; for more information, see the registered controlled trial protocol [26]. Levels of HRV-based stress and recovery did not differ between the participants of these two studies. However, the participants from the FFT study had more favorable body composition, better CRF, and a higher self-rated activity class. The characteristics of the participants of the present study are presented in Table 1.
Table 1

Characteristics of the study population

 

Mean

Range

Age (y)

34

26–40

Weight (kg)

83.5

51.3–123.2

Height (cm)

179.1

156.5–198.0

BMI (kg/[m]2)

26.0

19.8–35.0

Body fat%

24.5

7.6–41.8

VO2maxa (ml/kg/min)

39.6

23.3–72.6

MET indexb

3.5

0.1–27.7

Activity class (0–10)

4.4

0.0–9.5

Working time (h)

8.4

2.8–16.0

Sleeping time (h)

7.7

4.7–10.3

Bergen burnout inventory total scoresc

36

16–74

Firstbeat variables

  

Work stress (%)

72.8

14.5–99.2

Total stress (%)

49.9

12.6–79.1

Stress index

122.3

80.5–193.3

RMSSD (sleep)

56.78

22.6–150.6

Stress balance (−1 to 1)

0.69

−0.43–1.00

Total recovery (%)

29.3

4.6–67.0

Recovery index (24 h)

92.2

74.0–112.2

Recovery index (sleep)

111.4

45.8–246.9

(n = 81).

an = 78 bMET-h/day, based on retrospective PA interview cn = 80.

Both studies were conducted according to the ethical rules stated in the Declaration of Helsinki. All participants were informed about the study, and they signed written informed consent prior to any measurements. The study protocols were approved by the ethics committee of the Central Finland Health Care District.

Measurements

The stress and recovery variables were calculated objectively from the HR and HRV recordings, which were executed over three days using Bodyguard, a measurement device developed by Firstbeat Technologies Ltd. (Jyväskylä, Finland). The HRV recordings from these three days were downloaded with Firstbeat Health Software (version 3.1.1.0.). The Firstbeat Health Software first scans the recorded ambulatory RR interval data through an artifact detection filter to perform an initial correction of falsely detected, missed, and premature heart beats [27]. The data were transferred to the Matlab environment (R2007B), where the analyses of physiological variables describing stress and recovery were performed with the Firstbeat Analysis Server program (version 5.3.0.4). This program calculates HRV indices second-by-second using the short-time Fourier transform method and HR- and HRV-derived variables of respiration rate and oxygen consumption using neural network modeling of data [2729] (for more details, see the white paper by Firstbeat Technologies Ltd. [30]). The program also calculates second-by-second indices of recovery and stress, reflecting activities of the sympathetic (absolute stress vector, ASV) and parasympathetic (absolute recovery vector, ARV) nervous systems. The ASV is calculated from the HR, HF, and LF components and respiratory variables. The ASV is high when HR is elevated, HRV is reduced, and the frequency distribution of HRV is inconsistent because of changes in respiratory period. The ARV, which is calculated from the HR and HF components, is high when HR is close to the basic resting level and HRV is high and regular. The program takes into account individual basic resting HR and HRV values in the determination of the physiological states of the body, including PA, stress state, recovery state, or unrecognized state, and detected based on the above-mentioned variables. The stress and recovery variables used in the present study are based on Firstbeat Analysis Server analyses and are described in Table 2. These variables included the traditional HRV variable RMSSD (root mean square of successive differences), which is a time domain measure of HRV, and plots HRV as the change in normal RR intervals over time [31].
Table 2

The description of variables of stress and recovery based on firstbeat analysis server

Stress variables For non-exercise data segments, continuous indices of stress are used to identify the time when the body is in a stress state. Stress state is defined as an increased activation in the body when sympathetic nervous system activity is dominating and parasympathetic activation is decreased. Stress can be induced by external and internal stress factors, and the definition of the stress state does not take into account the nature of the stress response, i.e., whether it is positive or negative.

Work stress

The percentage of stress time during working hours

Total stress

The percentage of stress time in a workday (including leisure time and night)

Stress index

Absolute value characterizing the magnitude of stress processes on a workday (including leisure time and night)

Recovery variables, 24 hours Continuous indices of recovery are used to identify the time when the body is in a recovery state. Recovery state is defined as a decreased activation in the body during recovery, rest, and/or peaceful working. This state is related to the lack of external and internal stress factors, and parasympathetic activation is dominating.

Total recovery

The percentage of recovery time in a workday (including leisure time and night)

Recovery index

Absolute value characterizing the magnitude of recovery processes in a workday (including leisure time and night)

Recovery variables, sleep time

RMSSD

Average of the RMSSD (root mean square of successive differences) vector values during sleep periods. High RMSSD values are associated with increased parasympathetic activity and good recovery while low values during rest indicate poor recovery. The RMSSD value should be over 20 during sleep in a normal situation.

Stress balance

Indicates proportion of time of stress and recovery reactions during sleep periods in the measurement period; the used values are from −1 to 1. Values from 0.5 to 1 indicate good recovery; values from 0 to 0.5 indicate moderate recovery; and values from 0 to −1 indicate weak recovery.

Recovery index

The recovery index gives an estimate of a person’s recovery during sleep time. The 4-h window for determining the recovery index is set to start 30 min after going to bed. The index summarizes several factors, including RMSSD and other HRV variables, and takes into account the other functions of ANS while determining the value for the index.

The HRV data consisted of successfully recorded (RR intervals corrected max 25 percentages) workdays. The stress and recovery values were the mean values of the workdays with the number of days being one (n = 10), two (n = 70), or three (n = 1). The stress and recovery variables were determined from the whole recording day (24 h) or only from working hours or sleeping hours. Working hours, sleeping time, and alcohol consumption were determined from the participants’ measurement diaries, which they were advised to keep over the measurement period. Alcohol consumption was reported daily in standard units of approximately 12 g of ethanol (one unit: 33 centiliter [cl] beer or 12 cl red or white wine or 8 cl fortified wine or 4 cl liquor).

The Bergen Burnout Inventory (BBI) was used in the assessment of occupational burnout. This inventory is employed in some Finnish occupational health services to monitor and screen stress levels [32] and includes 15 questions concerning three subcategories of burnout: exhaustive fatigue, cynicism, and impaired occupational self-respect. Answers to the 15 questions were given using a Likert-type response scale, which was scored from 1 (totally disagree) to 6 (totally agree), and the total scores were between 15 and 90 [33].

Body weight and height were measured in the morning. The body mass index (BMI) is computed as weight/height2 (kg/m2). The whole body fat percentage was evaluated after fasting 10–12 h using dual-energy X-ray absorptiometry (GE Lunar Prodigy Advance, GE Healthcare). The software used for picture handling and analysis was enCORE™ 2009, version 13.20 (GE Healthcare).

The volume of PA (MET index) was retrospectively inferred using a modified version [34] of the Kuopio Ischemic Heart Disease Risk Factor Study Questionnaire [35], which included questions on leisure-time PA and PA during journeys to and from work, as well as daily activities such as gardening and berry picking. Monthly frequency, mean duration, and mean intensity of each form of activity were evaluated. A multiple of the resting metabolic rate (MET) was assigned for each activity to describe the intensity of the form of PA. The 3-month (BFH study) and 12-month (FFT study) MET indexes for each form of PA were calculated by multiplying the intensity (MET), duration (h), and monthly frequency of the activity, and the MET index was expressed as the sum score of different activities (MET-h/d).

The other method used to estimate PA level was self-reported activity class using a 0–10 scale to represent the activity level of the previous 2–3 months. The values from 0 to 7 were modified from Ross and Jackson’s [36] scale. The values from 7.5 to 10 were added by Firstbeat to include more seriously training individuals and athletes in the scale (for more details, see the white paper by Firstbeat Technologies Ltd. [37]).

CRF was assessed using the cycle ergometer test and a slightly modified World Health Organization protocol [38] with 2-min stages and 25 W/stage increases in workload. The test was submaximal in the BHF study and maximal with breath-by-breath respiratory gas-exchange analysis in the FFT study. The submaximal cycle ergometer test was ended after the workload during which the tested person achieved the submaximal HR level (85–88% from the maximal HR), which was defined from maximal HR based on the participant’s age (210 – (0.65 × age)). In the FFT study, three people did not participate in the maximal cycle ergometer test, and in the case of seven participants, a submaximal test was performed instead of the maximal test for health or motivational reasons. In the analyses, the values used for the maximal oxygen consumption (VO2max) are the results of submaximal workload-based calculations for all participants. VO2max was calculated from the submaximal test values of the maximal cycle ergometer test (n = 36). The correlation between this calculation from submaximal workloads and the value from the maximal oxygen consumption test was significant (P < 0.001).

Statistical analysis

Statistical analyses were performed using Mplus version 7 [39]. The statistical significance was set at P ≤ 0.05. The estimator used was maximum likelihood with robust standard errors that are robust against non-normality. The effect of clustered data (twins) was controlled using the type complex definition in Mplus. A few missing values (VO2max n = 3 and BBI n = 1) were inferred to be missing at random. Because the variables under interest are known to be related to each other, structural equation modeling was used to investigate how PA, fitness, and body composition are associated with HRV-based indicators of stress and recovery. In the first stage, estimates of each independent variable (body fat percentage, BMI, VO2max, MET index, and activity class) were determined to explain the dependent stress and recovery variables. In the second stage, to investigate how much of the variance of the stress and recovery variables is explained by PA, fitness, and body composition together, a latent factor variable was formed from the MET index, activity class, VO2max, body fat percentage, and BMI. In addition, structural equation modeling was used to investigate the association of the HRV-based stress and recovery variables with the BBI total scores.

The model fit was evaluated using the χ2 test, comparative fit index (CFI), Tucker Lewis Index (TLI), root mean square error of approximation (RMSEA), and standardized root mean square residual (SRMR). For a good-fitting model, the χ2 test is non-significant; CFI and TLI are at least 0.95; RMSEA is no more than 0.06; and SRMR no more than 0.08 [40].

Results

Body fat percentage was positively associated (β = 0.306, P = 0.005) whereas VO2max (β = − 0.311, P = 0.019), MET index (β = −0.304, P < 0.001), and activity class (β = − 0.326, P = 0.001) were negatively associated with HRV-based work stress. Similarly to works stress, the other HRV-based stress variables were associated with body composition, CRF, and PA. However, BMI was positively associated with total stress (β = 0.272, P = 0.003) and stress index (β = 0.478, P < 0.001). MET index was not associated with stress index (P = 0.440).

The indicators of recovery during night-time sleep (RMSSD, stress balance, recovery index) were negatively associated with body fat percentage and BMI and positively associated with VO2max. Activity class was positively associated with stress balance value (β = 0.237, P = 0.025).

Body fat percentage (β = −0.445, P < 0.001) and BMI (β = −0.355, P < 0.001) were negatively associated whereas VO2max (β = 0.465, P < 0.001), MET index (β = 0.356, P < 0.001), and activity class (β = 0.418, P < 0.001) were positively associated with total recovery. Body fat percentage (β = −0.343, P = 0.002) and BMI (β = −0.309, P = 0.001) were negatively associated whereas VO2max was positively associated with 24-h recovery index (β = 0.331, P = 0.010). MET index (P = 0.414) and activity class (P = 0.151) were not associated with 24-h recovery index.

Because the age range of the participants was relatively small and adjustment for age changed the results very minimally, the age-adjusted results are not shown. Adjustments for alcohol consumption, working time, or sleeping time influenced the results only minimally or modestly (Table 3). Adjustments for alcohol most influenced the associations of VO2max with the stress and recovery variables. The results are presented in greater detail in Table 3.
Table 3

The association of body composition, cardiorespiratory fitness, and physical activity with stress and recovery

 

No adjustments

Alcohol adjusted

Working time adjusted

Sleeping time adjusted

 

β

S.E.

P

β

S.E.

P

β

S.E.

P

β

S.E.

P

Body fat%

            

WS

0.306

0.108

0.005

0.288

0.114

0.012

0.315

0.114

0.006

0.313

0.115

0.006

TS

0.442

0.101

<0.001

0.405

0.107

<0.001

0.242

0.070

0.001

0.429

0.111

<0.001

SI

0.422

0.115

<0.001

0.395

0.108

<0.001

0.426

0.107

<0.001

0.425

0.107

<0.001

TR

−0.445

0.088

<0.001

−0.409

0.087

<0.001

−0.443

0.090

<0.001

−0.419

0.086

<0.001

RI

−0.343

0.113

0.002

−0.325

0.111

0.003

−0.338

0.103

0.001

−0.330

0.111

0.003

RMSSD

−0.293

0.123

0.017

−0.267

0.115

0.020

−0.285

0.109

0.009

−0.283

0.116

0.015

SB

−0.400

0.083

<0.001

−0.377

0.084

<0.001

−0.408

0.076

<0.001

−0.403

0.077

<0.001

RIS

−0.359

0.130

0.006

−0.328

0.127

0.010

−0.348

0.114

0.002

−0.351

0.125

<0.001

BMI

            

WS

0.162

0.086

0.060

0.139

0.085

0.101

0.176

0.083

0.033

0.171

0.086

0.047

TS

0.272

0.093

0.003

0.210

0.090

0.020

0.268

0.086

0.002

0.251

0.095

0.009

SI

0.478

0.102

<0.001

0.447

0.105

<0.001

0.458

0.091

<0.001

0.482

0.097

<0.001

TR

−0.355

0.084

<0.001

−0.297

0.078

<0.001

−0.351

0.084

<0.001

−0.309

0.084

<0.001

RI

−0.309

0.093

0.001

−0.287

0.094

0.002

−0.298

0.087

0.001

−0.289

0.092

0.002

RMSSD

−0.304

0.106

0.004

−0.274

0.106

0.009

−0.290

0.097

0.003

−0.289

0.107

0.007

SB

−0.379

0.095

<0.001

−0.325

0.101

0.001

−0.395

0.087

<0.001

−0.384

0.092

<0.001

RIS

−0.395

0.102

<0.001

−0.358

0.105

0.001

−0.373

0.095

<0.001

−0.382

0.104

<0.001

VO 2max

            

WS

−0.311

0.132

0.019

−0.282

0.147

0.055

−0.316

0.147

0.031

−0.317

0.145

0.029

TS

−0.416

0.126

0.001

−0.349

0.155

0.025

−0.411

0.145

0.004

−0.396

0.150

0.008

SI

−0.379

0.096

<0.001

−0.334

0.098

0.001

−0.377

0.094

<0.001

−0.375

0.095

<0.001

TR

0.465

0.111

<0.001

0.405

0.122

0.001

0.469

0.110

<0.001

0.429

0.107

<0.001

RI

0.331

0.128

0.01

0.302

0.139

0.030

0.319

0.133

0.016

0.307

0.143

0.032

RMSSD

0.262

0.125

0.037

0.222

0.125

0.076

0.247

0.123

0.045

0.242

0.130

0.062

SB

0.270

0.073

<0.001

0.145

0.091

0.113

0.283

0.073

<0.001

0.274

0.072

<0.001

RIS

0.294

0.125

0.019

0.243

0.131

0.063

0.275

0.121

0.023

0.276

0.136

0.042

MET index

            

WS

−0.304

0.113

<0.001

−0.287

0.087

0.001

−0.308

0.087

<0.001

−0.311

0.086

<0.001

TS

−0.385

0.145

0.008

−0.349

0.084

<0.001

−0.384

0.083

<0.001

−0.368

0.086

<0.001

SI

−0.096

0.124

0.440

−0.065

0.105

0.535

−0.098

0.107

0.360

−0.099

0.105

0.343

TR

0.356

0.131

<0.001

0.316

0.102

0.002

0.355

0.095

<0.001

0.321

0.097

0.001

RI

0.138

0.173

0.414

0.120

0.140

0.390

0.135

0.136

0.319

0.122

0.141

0.389

RMSSD

0.058

0.165

0.723

0.032

0.152

0.834

0.054

0.147

0.713

0.046

0.151

0.761

SB

0.097

0.089

0.255

−0.017

0.081

0.835

0.102

0.074

0.172

0.101

0.069

0.141

RIS

0.074

0.173

0.670

0.038

0.159

0.810

0.067

0.150

0.654

0.063

0.158

0.690

Activity class

            

WS

−0.326

0.102

0.001

−0.298

0.104

0.004

−0.330

0.107

0.002

−0.331

0.109

0.002

TS

−0.434

0.096

<0.001

−0.362

0.102

<0.001

−0.433

0.105

<0.001

−0.424

0.108

<0.001

SI

−0.222

0.113

0.048

−0.175

0.115

0.128

−0.224

0.115

0.051

−0.225

0.114

0.049

TR

0.418

0.097

<0.001

0.336

0.090

<0.001

0.417

0.098

<0.001

0.400

0.101

<0.001

RI

0.175

0.122

0.151

0.146

0.119

0.221

0.172

0.121

0.154

0.164

0.126

0.191

RMSSD

0.090

0.120

0.453

0.047

0.116

0.684

0.086

0.116

0.459

0.082

0.112

0.505

SB

0.237

0.106

0.025

0.082

0.095

0.390

0.241

0.107

0.024

0.239

0.106

0.024

RIS

0.132

0.129

0.305

0.076

0.125

0.544

0.127

0.125

0.311

0.125

0.132

0.343

Standardized estimate (β), standard error (S.E.), and P value for body fat%, BMI, VO2max, MET index, and activity class in relation to Firstbeat stress and recovery variables (work stress [WS], total stress [TS], stress index [SI], total recovery [TR], recovery index 24 h [RI], RMSSD, stress balance [SB], and recovery index sleep [RIS]).

A latent factor variable (Figure 1) was formed from body fat percentage, BMI, VO2max, MET index, and activity class because these variables correlated significantly. The latent variable was a significant explanatory variable for all of the stress and recovery variables except RMSSD. It explained total recovery best at 30.1% of its variance. Among stress variables, the variance of total stress was best explained (27.5%) by the latent factor, which also explained 14.1% of the variance of work stress. The estimated model fit to the data well (χ2 (35) = 41.09, P = 0.22, CFI = 0.99, TLI = 0.98, RMSEA = 0.046, SRMR = 0.06), and all the modification indices were lower than 4, which indicated that were no additional parameters that would have increased the fit of the model for the analysis.
Figure 1

The variance of stress and recovery accounted for by physical activity, cardiorespiratory fitness, and body composition. The factor loadings (standardized estimate and S.E.) of body fat%, BMI, VO2max, MET index, activity class, and the factor loadings (including R2) for Firstbeat variables: work stress (WS), total stress (TS), stress index (SI), total recovery (TR), recovery index 24 h (RI), RMSSD, stress balance (SB), and recovery index sleep (RIS).

In addition, objectively and subjectively measured stress levels were significantly associated (Table 4). HRV-based work stress (β = 0.272, P = 0.003) and total stress (β = 0.304, P = 0.001) explained the variation in total score of the BBI questionnaire. Adjustment for age changed the results only minimally, but after adjustment for body fat percentage, associations between work stress and the BBI total score did not persist.
Table 4

Stress and recovery in relation to the Bergen Burnout Inventory total scores

 

β

S.E.

P

Work stress

0.272

0.091

0.003

Total stress

0.304

0.095

0.001

Stress index

−0.080

0.118

0.500

Total recovery

−0.170

0.092

0.065

Recovery index (24 h)

−0.036

0.109

0.741

RMSSD

−0.009

0.101

0.925

Stress balance

−0.006

0.092

0.947

Recovery index (sleep)

−0.014

0.112

0.898

Standardized estimate (β), standard error (S.E.) and P-value.

Discussion

The results of this study showed that there were significant associations among PA, CRF, body composition, and novel HRV-based levels of stress and recovery in real life on workdays. Greater PA level, better CRF, and more favorable body composition were associated with lower stress levels during working hours as well as with lower stress levels and higher recovery levels throughout the whole day. Better recovery at night was associated with better CRF and a more favorable body composition. Additionally, the relationship between objectively measured stress and self-reported burnout symptoms was significant, and further analysis showed that body composition explained it in part.

The results of the present study are mainly in line with those of earlier studies that used traditional HRV measures. Favorable HRV profiles may be related to PA and CRF [15, 16] while stress appears to be linked with reduced leisure-time PA [41] and lower physical fitness [2, 42]. Earlier studies have proposed that fit individuals would show attenuated physiological reactivity to psychological stressors compared with sedentary individuals [43, 44]. A more recent meta-analytic review by Forcier et al. [2] also provided support for the cross-adaptation hypothesis and suggested that fitness training may increase the ability of cardiovascular systems to control responses to acute stressors and also speed cardiovascular recovery from stress.

The present study found that among men, unfavorable body composition was associated with higher objective stress and lower recovery. According to previous studies that have examined stress as a subjective experience, the association between mental stress and body composition has been inconsistent because evidence both supports [4549] and refutes [5053] it. Additionally, the results suggested that body fat percentage explained part of the association between objectively and subjectively measured levels of stress, indicating that body fat percentage affects both physiological stress and subjective experience of stress among men. However, a previous study by Nyberg et al. [46] found that the relationship between stress and unfavorable body composition may not be linear. In an analysis of pooled European data (n = 161,746), both weight gain and weight loss were associated with the onset of job strain so that the cross-sectional association between job strain and BMI formed a ‘U’-shaped curve [46].

The results showed that PA level was not associated with night-time recovery, although regular PA has consistently been associated with better sleep [54]. However, previous findings indicate that PA in the evening may delay the beginning of recovery. Physical fitness may have a role in recovery of HR and HRV after exercise, as could be seen even between trained and highly trained subjects, so that the recovery takes 60–90 min longer in trained subjects [55]. Recovery of HRV also depends on type, intensity, and duration of exercise [56]. In the study of Myllymäki et al. [57], a late-night exercise caused higher HR during the first three sleeping hours of an exercise day compared to a control day. Mischler et al. [58] reported that prolonged exercise had similar effects on HR during sleep. In addition to these effects of exercise on HR, HRV might decrease during the activity and, for instance, high-intensity exercises reduce HRV for 30–60 min afterward in trained and highly trained subjects [55, 56]. Hynynen et al. [59] found that both moderate and heavy endurance exercise even in the daytime has an increasing effect on HR and decreasing effect on HRV during sleep in a dose–response manner.

Subjective methods were used to assess the level of PA. The MET index was based on a detailed retrospective interview. This index correlated significantly with measured VO2max values, supporting its validity to assess the level of PA. Subjective assessment may include some reporting bias; however, although the method undertaken was not objective, it had the advantage of taking several months of subjective data into account. The other method was self-reported activity class, which is regularly reported before the HRV measurement by Bodyguard.

One limitation of the present study is that the effect of different intensities of PA was not investigated; however, both light and moderate to vigorous leisure-time activity are related to lower likelihood of burnout [41]. Another limitation is that fat distribution was not taken into account, even though, according to earlier studies, it may have a role in relation to the HRV profile [19, 60]. The nature of work in today’s society is rather complex, including inexact and variable working hours and a lack of precise division between work time and leisure time, which causes difficulties in determining the level of HRV-based stress during exact working hours.

The current results cannot be generalized to women because the study population consisted of both physically inactive and active young male employees, who were not heavy drinkers or on regular medication. This population was ideal for this study because alcohol intake, diseases, and medications influence ANS [8]. The effect of possible alcohol intake on HRV was taken into account in the analyses. The fact that our study population consisted of young employees who are still establishing their careers highlights the importance of the findings.

Stress can be explained through several biological processes, such as the inflammatory, dopaminergic, and neuroendocrine systems [61]. However, in the present study, the measurement method for stress and recovery was based on HR and HRV, which are affected by the ANS. Stressful situations result in accentuated sympathovagal antagonism, which may be explained by the interaction of acetylcholine and norepinephrine; consequently, HR may become remarkably unstable [62]. The findings of Tulppo et al. [62] suggest that traditional measures of HRV are not specific for measurement of accentuated sympathovagal interaction. The Firstbeat Analysis Server program connects information from these traditional measures of HRV, HR, and respiratory variables in the determination of stress and recovery states from individual HR and HRV values, but it does not separate eustress from distress.

Conclusions

This study used special HR- and HRV-based indicators of stress and recovery and found that PA, CRF, and body composition were associated with these indicators on workdays. On one hand, detailed interview-based MET index and self-reported activity class, and on the other hand, body fat percentage and BMI, associated rather similarly with the indicators of stress and recovery. This suggests that easily collected self-reported activity class and BMI can be used as indicators of PA and body composition in clinical work. The current findings support the usability of the objective indicators of stress (work stress and total stress) as they were associated with self-reported occupational burnout symptoms. Overall, the results support the usability of this HRV-based method in the evaluation of stress and recovery in line with some previous findings supporting its validity and reliability [25, 63].

Declarations

Acknowledgements

The Body & Future Health study was supported by the SalWe Research Program for Mind and Body (Tekes - the Finnish Funding Agency for Technology and Innovation grant 1104/10). The FitFatTwin study was supported by the Finnish Ministry of Education and Culture. Kaprio J. and the Finnish twin studies have been supported by the Academy of Finland (grants # 265240, 263278, 264146). T. Teisala’s work was supported by a personal grant from the Juho Vainio Foundation.

Authors’ Affiliations

(1)
Department of Health Sciences, University of Jyväskylä, Jyväskylä, Finland
(2)
Methodology Centre for Human Sciences, Faculty of Social Sciences, University of Jyväskylä, Jyväskylä, Finland
(3)
Turku PET Centre, University of Turku, Turku, Finland
(4)
Hjelt Institute, Department of Public Health, University of Helsinki, Helsinki, Finland
(5)
Institute of Public Health and Clinical Nutrition, Clinical Nutrition, University of Eastern Finland, Joensuu, Finland
(6)
Department of Biology of Physical Activity, University of Jyväskylä, Jyväskylä, Finland
(7)
Department of Mental Health and Substance Abuse Services, National Institute for Health and Welfare, Helsinki, Finland
(8)
Institute for Molecular Medicine (FIMM), University of Helsinki, Helsinki, Finland

References

  1. Gerber M, Pühse U: Do exercise and fitness protect against stress-induced health complaints? A review of the literature. Scand J Public Health 2009,37(Suppl 8):801–819.View ArticlePubMedGoogle Scholar
  2. Forcier K, Stroud LR, Papandonatos GD, Hitsman B, Reiches M, Krishnamoorthy J, Niaura R: Links between physical fitness and cardiovascular reactivity and recovery to psychological stressors: a meta-analysis. Health Psychol 2006,25(Suppl 6):723–739.View ArticlePubMedGoogle Scholar
  3. Tulppo MP, Hautala AJ, Mäkikallio TH, Laukkanen RT, Nissilä S, Hughson RL, Huikuri HV: Effects of aerobic training on heart rate dynamics in sedentary subjects. J Appl Physiol 2003,95(Suppl 1):364–372.View ArticlePubMedGoogle Scholar
  4. Sothmann MS, Buckworth J, Claytor RP, Cox RH, White-Welkley JE, Dishman RK: Exercise training and the cross-stressor adaptation hypothesis. Exerc Sport Sci Rev 1996,24(Suppl 1):267–288.PubMedGoogle Scholar
  5. Dishman RK, Nakamura Y, Garcia ME, Thompson RW, Dunn AL, Blair SN: Heart rate variability, trait anxiety, and perceived stress among physically fit men and women. Int J Psychophysiol 2000,37(Suppl 2):121–133.View ArticlePubMedGoogle Scholar
  6. Reed MJ, Robertson CE, Addison PS: Heart rate variability measurements and the prediction of ventricular arrhythmias. QJM 2005,98(Suppl 2):87–95.View ArticlePubMedGoogle Scholar
  7. McMillan DE: Interpreting heart rate variability sleep/wake patterns in cardiac patients. J Cardiovasc Nurs 2002,17(Suppl 1):69–81.View ArticlePubMedGoogle Scholar
  8. Acharya UR, Joseph KP, Kannathal N, Lim CM, Suri JS: Heart rate variability: a review. Med Bio Eng Comput 2006, 44: 1031–1051.View ArticleGoogle Scholar
  9. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology: Heart rate variability: standards of measurement, physiological interpretation and clinical use. Circulation 1996,93(Suppl 5):1043–1065.Google Scholar
  10. Malliani A, Pagani M, Lombardi F, Gerutti S: Cardiovascular neural regulation explored in the frequency domain. Circulation 1991,84(Suppl 2):482–492.View ArticlePubMedGoogle Scholar
  11. Vanderlei LCM, Pastre CM, Hoshi RA, Carvalho TD, Godoy MF: Basic notions of heart rate variability and its clinical applicability. Rev Bras Cir Cardiovasc 2009,24(Suppl 2):205–217.View ArticlePubMedGoogle Scholar
  12. Pumprla J, Howorka K, Groves D, Chester M, Nolan J: Functional assessment of heart rate variability: physiological basis and practical applications. Int J Cardiol 2002,84(Suppl 1):1–14.View ArticlePubMedGoogle Scholar
  13. Chandola T, Heraclides A, Kumari M: Psychophysiological biomarkers of workplace stressors. Neurosci Biobehav Rev 2010,35(Suppl 1):51–57.PubMed CentralView ArticlePubMedGoogle Scholar
  14. Pichot V, Roche F, Denis C, Garet M, Duverney D, Costes F, Barthélémy J-C: Interval training in elderly men increases both heart rate variability and baroreflex activity. Clin Auton Res 2005,15(Suppl 2):107–115.View ArticlePubMedGoogle Scholar
  15. Jurca R, Church TS, Morss GM, Jordan AN, Earnest CP: Eight weeks of moderate-intensity exercise training increases heart rate variability in sedentary postmenopausal women. Am Heart J 2004,147(Suppl 5):8–15.View ArticleGoogle Scholar
  16. Pober DM, Braun B, Freedson PS: Effects of a single bout of exercise on resting heart rate variability. Med Sci Sports Exerc 2004,36(Suppl 7):1140–1148.View ArticlePubMedGoogle Scholar
  17. Rennie KL, Hemingway H, Kumari M, Brunner E, Malik M, Marmot M: Effects of moderate and vigorous physical activity on heart rate variability in a British study of civil servants. Am J Epidemiol 2003,158(Suppl 2):135–143.View ArticlePubMedGoogle Scholar
  18. Buchheit M, Gindre C: Cardiac parasympathetic regulation: respective associations with cardiorespiratory fitness and training load. Am J Physiol Heart Circ Physiol 2006, 291: 451–458.View ArticleGoogle Scholar
  19. Gao YY, Lovejoy JC, Sparti A, Bray GA, Keys LK, Partington C: Autonomic activity assessed by heart rate spectral analysis varies with fat distribution in obese women. Obes Res 2012,4(Suppl 1):55–63.Google Scholar
  20. Jeong AK, Yong-Gyu P, Kyung-Hwan C, Myung-Ho H, Hee-Chul H, Youn-Seon C, Dokyung Y: Heart rate variability and obesity indices: emphasis on the response to noise and standing. J Am Board Fam Pract 2005,18(Suppl 2):97–103.Google Scholar
  21. Zahorska-Markiewicz B, Kuagowska E, Kucio C, Klin M: Heart rate variability in obesity. Int J Obes 1993,17(Suppl 1):21–23.Google Scholar
  22. Rissanen P, Franssila‒Kallunki A, Rissanen A: Cardiac parasympathetic activity is increased by weight loss in healthy obese women. Obes Res 2001,9(Suppl 10):637–643.View ArticlePubMedGoogle Scholar
  23. Martinmäki K, Rusko H, Kooistra L, Kettunen J, Saalasti S: Intraindividual validation of heart rate variability indexes to measure vagal effects on hearts. Am J Physiol Heart Circ Physiol 2006, 290: 640–647.View ArticleGoogle Scholar
  24. Firstbeat Technologies Ltd. http://www.firstbeat.com
  25. Uusitalo A, Mets T, Martinmäki K, Mauno S, Kinnunen U, Rusko H: Heart rate variability related to effort at work. Appl Ergonomics 2011,42(Suppl 6):830–838.View ArticleGoogle Scholar
  26. International standard randomized controlled trial number register http://www.controlled-trials.com/ISRCTN92130721/kujala
  27. Saalasti S: Neural networks for heart rate time series analysis. In PhD Thesis. Finland, Jyväskylä Studies in Computing 33, Jyväskylä: Department of Mathematical Information Technology; University of Jyväskylä; 2003.Google Scholar
  28. Kettunen J, Saalasti S, inventors; Firstbeat Technologies Ltd, assignee: Procedure for Deriving Reliable Information on Respiratory Activity from Heart Period Measurement. United States patent US 20050209521A1; 2005.Google Scholar
  29. Kettunen J, Saalasti S, inventors; Firstbeat Technologies Ltd, assignee: Procedure for Detection of Stress by Segmentation and Analyzing Heart Beat Signal. United States patent US 007330752B2; 2008.Google Scholar
  30. White papers http://www.firstbeat.com/physiology/white-papers
  31. Achten J, Jeukendrup AE: Heart rate monitoring. Sports Med 2003,33(Suppl 7):517–538.View ArticlePubMedGoogle Scholar
  32. Salmela-Aro K, Näätänen P, Nurmi J-E: The role of work-related personal projects during two burnout interventions: a longitudinal study. Work Stress 2004,18(Suppl 3):208–230.View ArticleGoogle Scholar
  33. Näätänen P, Aro A, Matthiesen S, Salmela-Aro K: Bergen Burnout Indicator 15. Helsinki: Edita; 2003.Google Scholar
  34. Waller K, Kaprio J, Kujala UM: Associations between long-term physical activity, waist circumference and weight gain: a 30-year longitudinal twin study. Int J Obes 2007,32(Suppl 2):353–361.Google Scholar
  35. Lakka TA, Salonen JT, The physical activity questionnaires of the Kuopio Ischemic Heart Disease Study (KIHD): A collection of physical activity questionnaires for health related research. Med Sci Sports Exerc 1997, 29: 46–58.Google Scholar
  36. Ross RM, Jackson AS: Exercise Concepts, Calculations, and Computer Applications. Carmel, Indiana: Benchmark press; 1990.Google Scholar
  37. EPOC based training effect assessment http://www.firstbeat.fi/userData/firstbeat/download/white_paper_training_effect.pdf
  38. Andersen KL, Shephard RJ, Denolin H, Varnauskas E, Masironi R: Fundamentals of Exercise Testing. Geneva: World Health Organization; 1971.Google Scholar
  39. Muthén LK, Muthén BO: Mplus User’s Guide. 7th edition. Los Angeles: CA: Muthén & Muthén; 1998–2012.Google Scholar
  40. Hu LT, Bentler PM: Cutoff criteria for fit indexes in covariance structure analysis: conventional criteria versus new alternatives. Struct Equation Modeling 1999,6(Suppl 1):1–55.View ArticleGoogle Scholar
  41. Jonsdottir IH, Rödjer L, Hadzibajramovic E, Börjesson M, Ahlborg G: A prospective study of leisure-time physical activity and mental health in Swedish health care workers and social insurance officers. Prev Med 2010,51(Suppl 5):373–377.View ArticlePubMedGoogle Scholar
  42. Holtermann A, Mortensen OS, Burr H, Søgaard K, Gyntelberg F, Suadicani P: Physical fitness and perceived psychological pressure at work: 30-year ischemic heart disease and all-cause mortality in the copenhagen male study. J Occup Env Med 2011,53(Suppl 7):743–750.View ArticleGoogle Scholar
  43. Anshel MH: Effect of chronic aerobic exercise and progressive relaxation on motor performance and affect following acute stress. Behav Med 1996,21(Suppl 4):186–196.View ArticlePubMedGoogle Scholar
  44. van Doornen LJ, de Geus EJ: Aerobic fitness and the cardiovascular response to stress. Psychophysiol 1989,26(Suppl 1):17–28.View ArticleGoogle Scholar
  45. Kitaoka-Higashiguchi K, Morikawa Y, Miura K, Sakurai M, Ishizaki M, Kido T, Naruse Y, Nakagawa H: Burnout and risk factors for arteriosclerotic disease: follow-up study. J Occup Health 2009,51(Suppl 2):123–131.View ArticlePubMedGoogle Scholar
  46. Nyberg ST, Heikkilä K, Fransson EI, Alfredsson L, De Bacquer D, Bjorner JB, Bonenfant S, Borritz M, Burr H, Casini A, Clays E, Dragano N, Erbel R, Geuskens GA, Goldberg M, Hooftman WE, Houtman IL, Jöckel KH, Kittel F, Knutsson A, Koskenvuo M, Leineweber C, Lunau T, Madsen IEH, Magnusson Hanson LL, Marmot MG, Nielsen ML, Nordin M, Oksanen T, Pentti J, et al.: Job strain in relation to body mass index: pooled analysis of 160 000 adults from 13 cohort studies. J Intern Med 2012, 272: 65–73.PubMed CentralView ArticlePubMedGoogle Scholar
  47. Block JP, He Y, Zaslavsky AM, Ding L: Psychosocial stress and change in weight among US adults. Am J Epidemiol 2009,170(Suppl 2):181–192.PubMed CentralView ArticlePubMedGoogle Scholar
  48. Brunner EJ, Chandola T, Marmot MG: Prospective effect of job strain on general and central obesity in the Whitehall II Study. Am J Epidemiol 2007,165(Suppl 7):828–837.View ArticlePubMedGoogle Scholar
  49. Hannerz H, Albertsen K, Nielsen ML, Tüchsen F, Burr H: Occupational factors and 5-year weight change among men in a danish national cohort. Health Psychol 2004,23(Suppl 3):283–288.View ArticlePubMedGoogle Scholar
  50. Armon G, Shirom A, Berliner S, Shapira I, Melamed S: A prospective study of the association between obesity and burnout among apparently healthy men and women. J Occup Health Psychol 2008,13(1):43–57.View ArticlePubMedGoogle Scholar
  51. Ostry AS, Radi S, Louie AM, LaMontagne AD: Psychosocial and other working conditions in relation to body mass index in a representative sample of Australian workers. BMC Public Health 2006,6(Suppl 1):53–61.PubMed CentralView ArticlePubMedGoogle Scholar
  52. Overgaard D, Gyntelberg F, Heitmann BL: Psychological workload and body weight: is there an association? A review of the literature. Occup Med 2004,54(Suppl 1):35–41.View ArticleGoogle Scholar
  53. Niedhammer I, Goldberg M, Leclerc A, David S, Bugel I, Landre M-F: Psychosocial work environment and cardiovascular risk factors in an occupational cohort in France. J Epidemiol Community Health 1998,52(Suppl 2):93–100.PubMed CentralView ArticlePubMedGoogle Scholar
  54. Driver HS, Taylor SR: Exercise and sleep. Sleep Med Rev 2000,4(Suppl 4):387–402.View ArticlePubMedGoogle Scholar
  55. Seiler S, Haugen O, Kuffel E: Autonomic recovery after exercise in trained athletes: intensity and duration effects. Med Sci Sports Exerc 2007,39(Suppl 8):1366–1373.View ArticlePubMedGoogle Scholar
  56. Kaikkonen P, Rusko H, Martinmäki K: Post‒exercise heart rate variability of endurance athletes after different high‒intensity exercise interventions. Scand J Med Sci Sports 2008,18(Suppl 4):511–519.PubMedGoogle Scholar
  57. Myllymäki T, Kyröläinen H, Savolainen K, Hokka L, Jakonen R, Juuti T, Martinmäki K, Kaartinen J, Kinnunen M‒L, Rusko H: Effects of vigorous late‒night exercise on sleep quality and cardiac autonomic activity. J Sleep Res 2011, 20: 146–53.View ArticlePubMedGoogle Scholar
  58. Mischler I, Vermorel M, Montaurier C, Mounier R, Pialoux V, Péquignot J-M, Cottet-Emard J-M, Coudert J, Fellmann N: Prolonged daytime exercise repeated over 4 days increases sleeping heart rate and metabolic rate. Can J Appl Physiol 2003,28(Suppl 4):616–629.View ArticlePubMedGoogle Scholar
  59. Hynynen E, Vesterinen V, Rusko H, Nummela A: Effects of moderate and heavy endurance exercise on nocturnal HRV. Int J Sports Med 2010,31(Suppl 6):428–432.View ArticlePubMedGoogle Scholar
  60. Grassi G, Dell’Oro R, Facchini A, Trevano FQ, Bolla GB, Mancia G: Effect of central and peripheral body fat distribution on sympathetic and baroreflex function in obese normotensives. J Hypertens 2004,22(Suppl 12):2363–2369.View ArticlePubMedGoogle Scholar
  61. Hamer M: Psychosocial stress and cardiovascular disease risk: the role of physical activity. Psychosom Med 2012,74(Suppl 9):896–903.View ArticlePubMedGoogle Scholar
  62. Tulppo MP, Mäkikallio TH, Seppänen T, Airaksinen JK, Huikuri HV: Heart rate dynamics during accentuated sympathovagal interaction. AJP 1998,274(Suppl 3):H810-H816.Google Scholar
  63. Hynynen E, Konttinen N, Kinnunen U, Kyröläinen H, Rusko H: The incidence of stress symptoms and heart rate variability during sleep and orthostatic test. Eup J Appl Physiol 2011,111(Suppl 5):733–741.View ArticleGoogle Scholar

Copyright

© Teisala et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Advertisement