The association between multisite musculoskeletal pain and cardiac autonomic modulation during work, leisure and sleep – a cross-sectional study
BMC Musculoskeletal Disorders volume 19, Article number: 405 (2018)
The prevention and rehabilitation of multisite musculoskeletal pain would benefit from studies aiming to understand its underlying mechanism. Autonomic imbalance is a suggested mechanism for multisite pain, but hardly been studied during normal daily living. Therefore, the aim of the study is to investigate the association between multisite musculoskeletal pain and cardiac autonomic modulation during work, leisure and sleep.
This study is based on data from the “Danish Physical activity cohort with objective measurements” among 568 blue-collar workers. Pain intensity scales were dichotomized according to the median of each scale, and the number of pain sites was calculated. No site was regarded as the pain-free, one site was considered as single-site musculoskeletal pain and pain in two or more sites was regarded as multisite musculoskeletal pain. Heart rate variability (HRV) was measured by an electrocardiogram system (ActiHeart) and physical activity using accelerometers (Actigraph). Crude and adjusted linear mixed models were applied to investigate the association between groups and cardiac autonomic regulation during work, leisure and sleep.
There was no significant difference between groups and no significant interaction between groups and domains in the crude or adjusted models for any HRV index. Significant differences between domains were found in the crude and adjusted model for all indices, except SDNN; sleep time showed higher values than leisure and work time, except for LF and LF/HF, which were higher during work.
This cross-sectional study showed that multisite musculoskeletal pain is not associated with imbalanced cardiac autonomic regulation during work, leisure and sleep time.
Chronic musculoskeletal pain has high prevalence and large consequences for the society [12, 34]. Although most studies focus on pain localized in a particular body region (single-site pain), such as low back pain or neck/shoulder pain , musculoskeletal pain usually occurs concurrently in several anatomical sites [4, 7]. This condition is called multisite pain  and has been shown to be associated with increased healthcare utilization, sick leave, early retirement, sickness and social welfare benefit [11, 13, 17].
In contrast to single-site pain, which is considered to be due to overload or insufficient use of a particular body region , multisite pain may be driven by more generalized mechanisms, such as an imbalance in autonomic cardiac modulation .
Autonomic cardiac modulation can be assessed by heart rate variability (HRV) reflecting parasympathetic and sympathetic regulation of beat-to-beat heart rate. Recent systematic reviews show moderate evidence supporting a decrease in parasympathetic modulation in chronic pain patients [20, 40]. However, we are aware of only one study (the Netherlands Study of Depression and Anxiety - NESDA) investigating the association between multisite musculoskeletal pain and HRV, finding no relationship between HRV and pain onset  or recovery . However, more studies on the association between multisite musculoskeletal pain and HRV are required before any conclusion can be drawn.
For better prevention and rehabilitation of multisite musculoskeletal pain, it is important to understand the underlying mechanism , e.g., if it is related to an imbalanced autonomic regulation. It is also relevant to investigate large populations with a wide variation in the number of pain sites, and the variation of autonomic activity throughout daily living (not only artificial conditions). Moreover, because physical activity and body postures influence HRV in ambulatory recordings [2, 5, 33], it is important to use valid technical information of physical activity and body postures during the measurement of HRV. Also, a previous study of HRV in a working population showed differences between work and leisure time only for the sitting posture . None of the previous studies on multisite musculoskeletal pain and HRV has taken all these factors into account. Thus, the aim of this study was to investigate the association between multisite musculoskeletal pain and cardiac autonomic modulation during work, leisure and sleep, and the interaction between multisite musculoskeletal pain and time domains.
Study population and exclusion criteria
This is a cross-sectional study based on data from “The Danish Physical activity cohort with objective measurements” (DPhacto) cohort, conducted on blue-collar workers recruited in the cleaning in public and private sector (i.e. hospitals, schools, municipalities, and private firms), manufacturing and production companies in metal, plastic, and food industries, and transportation (i.e. mail, and parcel service companies). The recruitment was performed in collaboration with a labor union and the data were collected between 2012 and 2013. The inclusion criteria were: to be allowed to participate during the paid working time, to be employed for more than 20 h per week and being between 18 and 65 years. Exclusion criteria were declining to sign the informed consent, pregnancy, fever on the testing day, and allergy to adhesives.
Among 2107 potentially eligible workers in the DPhacto cohort, objective measurement data from 759 blue-collar workers were available for analysis. Workers having less than 4 h of valid HRV recordings during work, leisure and sleep time (n = 163) and with no information about pain intensity (n = 3) were excluded, resulting in a final sample of 568 blue-collar workers (Fig. 1). The response rate was 63%. Detailed information about the DPhacto cohort can be found elsewhere .
Multisite musculoskeletal pain
Musculoskeletal pain was assessed by modified questions from the validated Nordic Musculoskeletal Questionnaire  with questions on pain intensity in seven anatomical areas (neck/shoulders; elbows; hands/wrists; low back; hips; knees and feet/ankles) during the past three months using a scale from 0 (no pain) to 10 (worst pain). Pain intensity scales were dichotomized from the median, i.e. less than median = ‘no pain’, more than median = ‘pain’ . The cut-off values were 3 points for neck/shoulder and low back, and 0 points for elbows, hands/wrists, hips, knees, feet/ankles. The elbows and hands/wrists were grouped to represent the upper limbs and the hips; knees and feet/ankles represented the lower limbs. The four regions were summed to provide information about the number of pain sites (0 = none to 4 = 4 pain sites). Pain in none of the regions was regarded as ‘pain-free’, one pain site was regarded as ‘single-site pain’, while the pain in two or more sites was defined as ‘multisite musculoskeletal pain’.
Technical measurements of heart rate variability and physical activity
HRV was measured by the ActiHeart system (Camntech Ltd., Cambridge, UK) with electrocardiography sensitivity of 0.250 mV. The sensor was attached by a two-led configuration at the recommended position . The analogue signal was band-pass filtered (10–35 Hz), sampled with a frequency of 128 Hz, and processed by a real-time QRS-detection algorithm to achieve a 1 ms time resolution of the RR intervals. Since data was collected during daily living conditions, respiratory rate was not controlled. Abnormal beats were removed using an automatic algorithm before analyzing HRV .
HRV data were analyzed using a robust method  from 5-min windows with less than 10% erroneous interbeat intervals (IBI). For time domain, the measures obtained were the standard deviation of R-R intervals (SDNN), which is a measure of overall variability; and the root mean square of successive differences of R-R intervals (RMSSD), a measure of beat-to-beat variability, which is related to the vagal modulation. For frequency domain, the low (LF, 0.04–0.15 Hz) and high frequency (HF, 0.15–0.40 Hz) components were analyzed, as well as the sympathovagal balance (LF/HF ratio). HF indicates the parasympathetic modulation of the cardiac rhythm, while LF is an indicator of both sympathetic and parasympathetic cardiac modulations [26, 28].
Physical activity and body posture were objectively measured using multiple accelerometers (ActiGraph GT3X+, Actigraph, Florida, USA). The accelerometers were attached to the thigh and upper back for several days, including work, leisure and sleep time. A diary was also filled out by the workers, including information about times getting up in the morning, starting and finishing work, and going to bed, as well as times of the reference position (upright stance) for calibration of the accelerometer records .
Physical activity and HRV data were processed using the Acti4 software (The National Research Centre for the Working Environment, Copenhagen, Denmark and BAuA, Berlin, Germany). The Acti4 software classifies different physical activities (walking, moving, cycling and running) and body postures (sitting, standing, and lying down) with high sensitivity and specificity . The HRV indices obtained during 5-min non-overlap intervals in sitting posture at work and leisure were analyzed, as well as three periods with the lowest R-R intervals during nocturnal sleep without movement .
Assessment of individual and occupational factors
A self-reported questionnaire was administered to the workers including age, gender, alcohol and tobacco use, medication prescription, job seniority, lifting and carrying during work. Height (cm) was measured using a scale (Seca, model 123) and weight (kg) was measured by a digital scale (Tanita model BC 418 MA). Body mass index (BMI) was calculated according to the formulae BMI = weight (kg)/height2 (m).
Groups were compared for the characterization variables using one-way ANOVA and Tukey’s post hoc tests for continuous variables and Chi-square test for categorical variables. The HRV indices, except IBI, showed a non-normal distribution (Kolmogorov Smirnov test; see Additional file 1: Figures S1, S2 and S3). Thus a natural logarithm (ln) transformation was applied.
Linear mixed models with two fixed factors were applied to verify the association between groups (pain-free, single-site pain, multisite pain), domains (work, leisure, sleep), and interaction of groups and domains. Linear mixed models were chosen because it includes fixed and random effects and increase the study power including in the analysis subjects with missing data. Subject and intercept were included as random effects. The covariance type was unstructured, and the restricted maximum likelihood (REML) estimation method was chosen. When the interaction was significant, the mean difference (MD), standard error (SE) and P value for the pairwise comparison, based on estimated marginal means, were reported.
Crude and adjusted models were tested. In the adjusted model, age, sex, BMI, tobacco use, objectively measured moderate and vigorous physical activity (i.e. fast walking, running, stair climbing and cycling) at work and leisure, and sitting time at work and leisure were included as covariates since these factors may affect both multisite pain and HRV [1, 10, 19, 21, 29, 32]. Sensitivity analyses, excluding workers reporting prescribed medication in the last three months and considering the cut-off point of zero to dichotomize the pain intensity scales, were also performed. Stratified analyses were applied to explore possible effect modification of age (< 50 years; ≥50 years) and sex (male; female). All analyses were performed using SPSS software (version 24.0), and the significance level was set at 1% to control for type I error in multiple comparison tests.
About 44% of the sample was composed of females, the mean age was 45 years and mean BMI was 27 kg/m2. About 29% of the workers smoked daily or occasionally; 19% reported to use analgesic medication, 12% antihypertensive, 3% heart medication and 3% antidepressants. Most of the workers were from the manufacturing sector (72%), and 41% reported to carry and lift for at least half of the work time. The prevalence of single-site pain was 23% (CI 95%, 20–27%), and the prevalence of multisite pain was 63% (CI 95%, 58–66%). In both groups, the most affected body parts were knees (single-site pain: 32%; multisite pain: 67%), lower back (single-site pain: 13%; multisite pain: 65%) and neck/shoulder (single-site pain: 15%; multisite pain: 65%).
The three groups were similar in most of the sociodemographic variables, except for a larger proportion of lifting and carrying almost all the time to ¼ of the time in the multisite pain group. Pain intensity, the number of pain sites and the proportion of workers with prescribed analgesic medication were also higher in the multisite pain group compared with the pain-free and single-site pain groups. There were no differences between the groups for the other types of medication (Table 1).
Table 2 shows the mean and interquartile range for the HRV indices during work, leisure and sleep time for each group. The between groups comparison showed similar values for all indices during work, leisure and sleep time.
The results from the linear mixed models showed no significant differences between groups and no significant interaction between groups and domains in the crude or adjusted models for any HRV index (Table 3). Significant differences between domains were found in the crude and adjusted model for all indices, except SDNN; sleep time showed higher values than leisure and work time, except for LF and LF/HF, which were higher during work.
The sensitivity analysis, excluding workers with prescribed medication in the last three months (n = 234), showed no significant differences between groups and significant differences among domains for all indices in the crude and adjusted model, except for SDNN. The interaction between groups and domains was not significant for all indices. When using a more restrictive definition of a pain-free group, i.e., including in the pain-free group workers reporting no pain in all body parts (cut-off point = 0), the sensitivity analysis also showed no significant differences between groups, significant differences between domains (except for SDNN) and no significant interaction between groups and domains for all indices (Additional file 1: Table S1). The same results were found in the stratified analyses for age (Additional file 1: Table S2) and sex (Additional file 1: Table S3).
This cross-sectional study showed no association between multisite musculoskeletal pain and cardiac autonomic modulation during work, leisure and sleep.
Previous systematic reviews indicated that chronic pain is associated with a decrease in parasympathetic modulation [20, 40]. However, in agreement with our findings, Generaal et al. [14, 15] showed that autonomic cardiac modulation was not impaired with chronic multisite pain. One possible explanation for this divergence between the systematic reviews and the MSP studies could be that the results from the systematic reviews showing an autonomic imbalance in chronic pain are mainly based on studies with fibromyalgia patients. So, the autonomic imbalance involved in the pathophysiology of fibromyalgia  does not necessarily extend to active workers with multisite pain.
Another possible explanation for our findings can be related to the study population and data collection conditions. Most studies are based on patient samples evaluated in artificial conditions which can differ from active workers concerning physical, cognitive and psychosocial characteristics. Thus, it is possible that people with pain who is still at work are healthier, more physically active and have better coping mechanisms with their pain compared to clinical samples with pain [9, 24], which may be reflected in better autonomic function throughout daily living conditions, including work, leisure and sleep time. It is still possible that there are sub-groups of people with other features of multisite pain not captured in our study. Likewise, other factors like sleep disturbance, depressed mood, somatising tendency and psychosocial aspects of work may be involved in chronic multisite pain [6, 38]. These factors were not examined in our study, which is a limitation. Our findings may also suggest publication bias of positive results in previous studies .
The sensitivity analysis for medication yielded the same findings as the main analysis, including workers with prescribed medication. Additionally, the adoption of a more restrictive definition of pain-free resulted in the same findings. Also, stratified analysis on sex and age did not show association between multisite musculoskeletal pain and cardiac autonomic modulation during work, leisure and sleep.
The prevalence of multisite pain was very high in our sample, as it affected 63% of the DPhacto blue-collar workers, with the mean peak pain intensity of 6.9 points on a 0 to 10 scale. This high prevalence may be specific to the study population consisting of blue-collar workers, i.e., it refers to a disadvantaged socioeconomic group. Other studies have also shown a high prevalence of multisite pain in different populations, ranging from 35 to 64% [17, 30]. Although our findings showed that the imbalance of the autonomic modulation was not associated with multisite pain, this issue still deserves attention, since the theoretical framework support this relationship [25, 42]. Future studies should have a longitudinal design to verify if the autonomic imbalance precedes the occurrence of multisite pain.
Strengths and limitations
A strength of our study is the large sample size, which allowed for stratified analyses. A further strength was the relatively homogenous group of blue-collar workers, which minimized potential socioeconomic confounding. The control for lifestyle factors such as smoking, physical activity, sitting time, and individual characteristics such as sex, age, and BMI is also highly relevant as these factors are closely related to HRV and pain. A potential limitation is the lack of control for the respiration rate, circadian variation, sleep quality and psychosocial aspects of work. Additionally, we have only looked at HRV during sleep, work and leisure, while there exist several other ways of evaluating autonomic function, such as assessing autonomic reactivity to functional tests (e.g. Valsalva, cold pressor and handgrip tests). Finally, the cross-sectional design of this study does not allow determining whether an autonomic imbalance may occur before the development of multisite pain.
This cross-sectional study showed that multisite musculoskeletal pain is not associated with imbalanced cardiac autonomic regulation during work, leisure and sleep time.
Analysis of variance
Body mass index
Danish Physical activity cohort with objective measurements
Heart rate variability
Inter beat intervals
Moderate and vigorous physical activity
Restricted maximum likelihood
Root mean square of successive differences of R-R intervals
Standard deviation of R-R intervals
Acharya UR, Paul JK, Kannathal N, Lim CM, Suri JS. Heart rate variability: a review. Med Biol Eng Comput. 2006;44(12):1031–51.
Bernardi L, Valle F, Coco M, Calciati A, Sleight P. Physical activity influences heart rate variability and very-low-frequency components in Holter electrocardiograms. Cardiovasc Res. 1996;32(2):234–7.
Brage S, Brage N, Franks PW, Ekelund U, Wareham NJ. Reliability and validity of the combined heart rate and movement sensor Actiheart. Eur J Clin Nutr. 2005;59(4):561–70. https://doi.org/10.1038/sj.ejcn.1602118.
Carnes D, Parsons S, Ashby D, Breen A, Foster NE, Pincus T, Vogel S, Underwood M. Chronic musculoskeletal pain rarely presents in a single body site: results from a UK population study. Rheumatology. 2007;46(7):1168–70. https://doi.org/10.1093/rheumatology/kem118.
Chan H-L, Lin M-A, Chao P-K, Lin C-H. Correlates of the shift in heart rate variability with postures and walking by time-frequency analysis. Comp Meth Program Biomed. 2007;86(2):124–30. https://doi.org/10.1016/j.cmpb.2007.02.003.
Coggon D, Ntani G. Trajectories of multisite musculoskeletal pain and implications for prevention. Occup Environ Med. 2017;74(7):465–6. https://doi.org/10.1136/oemed-2016-104196.
Coggon D, Ntani G, Palmer KT, Felli VE, Harari R, Barrero LH, Felknor SA, Gimeno D, Cattrell A, Vargas-Prada S, Bonzini M, Solidaki E, Merisalu E, Habib RR, Sadeghian F, Masood Kadir M, Warnakulasuriya SS, Matsudaira K, Nyantumbu B, Sim MR, Harcombe H, Cox K, Marziale MH, Sarquis LM, Harari F, Freire R, Harari N, Monroy MV, Quintana LA, Rojas M, Salazar Vega EJ, Harris EC, Serra C, Martinez JM, Delclos G, Benavides FG, Carugno M, Ferrario MM, Pesatori AC, Chatzi L, Bitsios P, Kogevinas M, Oha K, Sirk T, Sadeghian A, Peiris-John RJ, Sathiakumar N, Wickremasinghe AR, Yoshimura N, Kelsall HL, Hoe VC, Urquhart DM, Derrett S, McBride D, Herbison P, Gray A. Patterns of multisite pain and associations with risk factors. Pain. 2013;154(9):1769–77. https://doi.org/10.1016/j.pain.2013.05.039.
Croft P, Dunn KM, Von Korff M. Chronic pain syndromes: you can't have one without another. Pain. 2007;131(3):237–8. https://doi.org/10.1016/j.pain.2007.07.013.
de Vries HJ, Brouwer S, Groothoff JW, Geertzen JH, Reneman MF. Staying at work with chronic nonspecific musculoskeletal pain: a qualitative study of workers’ experiences. BMC Musculoskelet Disord. 2011;12:126. https://doi.org/10.1186/1471-2474-12-126.
Dinas PC, Koutedakis Y, Flouris AD. Effects of active and passive tobacco cigarette smoking on heart rate variability. Int J Cardiol. 2013;163(2):109–15. https://doi.org/10.1016/j.ijcard.2011.10.140.
Eckhoff C, Straume B, Kvernmo S. Multisite musculoskeletal pain in adolescence as a predictor of medical and social welfare benefits in young adulthood: the Norwegian Arctic adolescent health cohort study. Eur J Pain. 2017;21:1697–706. https://doi.org/10.1002/ejp.1078.
Fayaz A, Croft P, Langford RM, Donaldson LJ, Jones GT. Prevalence of chronic pain in the UK: a systematic review and meta-analysis of population studies. BMJ Open. 2016;6(6):e010364. https://doi.org/10.1136/bmjopen-2015-010364.
Fernandes RCP, Burdorf A. Associations of multisite pain with healthcare utilization, sickness absence and restrictions at work. Int Arch Occup Environ Health. 2016;89(7):1039–46. https://doi.org/10.1007/s00420-016-1141-7.
Generaal E, Vogelzangs N, Macfarlane GJ, Geenen R, Smit JH, de Geus EJ, Penninx BW, Dekker J. Biological stress systems, adverse life events and the onset of chronic multisite musculoskeletal pain: a 6-year cohort study. Ann Rheum Dis. 2016;75(5):847–54. https://doi.org/10.1136/annrheumdis-2014-206741.
Generaal E, Vogelzangs N, Macfarlane GJ, Geenen R, Smit JH, de Geus EJ, Dekker J, Penninx BW. Biological stress systems, adverse life events, and the improvement of chronic multisite musculoskeletal pain across a 6-year follow-up. J Pain. 2017;18(2):155–65. https://doi.org/10.1016/j.jpain.2016.10.010.
Hallman DM, Birk Jørgensen M, Holtermann A. On the health paradox of occupational and leisure-time physical activity using objective measurements: effects on autonomic imbalance. PLoS One. 2017;12(5):e0177042. https://doi.org/10.1371/journal.pone.0177042.
Haukka E, Kaila-Kangas L, Ojajärvi A, Saastamoinen P, Holtermann A, Jørgensen MB, Karppinen J, Heliövaara M, Leino-Arjas P. Multisite musculoskeletal pain predicts medically certified disability retirement among Finns. Eur J Pain. 2015;19(8):1119–28. https://doi.org/10.1002/ejp.635.
Jørgensen MB, Korshøj M, Lagersted-Olsen J, Villumsen M, Mortensen OS, Skotte J, Søgaard K, Madeleine P, Thomsen BL, Holtermann A. Physical activities at work and risk of musculoskeletal pain and its consequences: protocol for a study with objective field measures among blue-collar workers. BMC Musculoskelet Disord. 2013;20(14):213. https://doi.org/10.1186/1471-2474-14-213.
Kamaleri Y, Natvig B, Ihlebaek CM, Benth JS, Bruusgaard D. Number of pain sites is associated with demographic, lifestyle, and health-related factors in the general population. Eur J Pain. 2008;12(6):742–8.
Koenig J, Falvay D, Clamor A, Wagner J, Jarczok MN, Ellis RJ, Weber C, Thayer JF. Pneumogastric (vagus) nerve activity indexed by heart rate variability in chronic pain patients compared to healthy controls: a systematic review and meta-analysis. Pain Physician. 2016;19(1):E55–78.
Koenig J, Jarczok MN, Warth M, Ellis RJ, Bach C, Hillecke TK, Thayer JF. Body mass index is related to autonomic nervous system activity as measured by heart rate variability – a replication using short term measurements. J Nutr Health Aging. 2014;18(3):300–2. https://doi.org/10.1007/s12603-014-0022-6.
Kristiansen J, Korshøj M, Skotte JH, Jespersen T, Søgaard K, Mortensen OS, Holtermann A. Comparison of two systems for long-term heart rate variability monitoring in free-living conditions—a pilot study. Biomed Eng Online. 2011. https://doi.org/10.1186/1475-925X-10-27.
Kuorinka I, Jonsson B, Kilbom Å, Vinterberg H, Biering-Sörensen F, Anderson G, Jörgensen K. Standardised Nordic questionnaires for the analysis of musculoskeletal symptoms. Appl Ergon. 1987;18:233–7.
Linton SJ, Buer N. Working despite pain: factors associated with work attendance versus dysfunction. Int J Behav Med. 1995;2(3):252–62.
Maletic V, Raison CL. Neurobiology of depression, fibromyalgia and neuropathic pain. Front Biosci. 2009;14:5291–338.
Malik M, Bigger JT, Camm AJ, Kleiger RE, Malliani A, Moss AJ, Schwartz PJ. Heart rate variability standards of measurement, physiological interpretation, and clinical use. Eur Heart J. 1996;17(3):354–81.
Martínez-Lavín M. Is fibromyalgia a generalized reflex sympathetic dystrophy? Clin Exp Rheumatol. 2001;19(1):1–3.
Michael S, Graham KS, Oam DGM. Cardiac autonomic responses during exercise and postexercise recovery using heart rate variability and systolic time intervals - a review. Front Physiol. 2017;29(8):301. https://doi.org/10.3389/fphys.2017.00301.
Molfino A, Fiorentini A, Tubani L, Martuscelli M, Rossi Fanelli F, Laviano A. Body mass index is related to autonomic nervous system activity as measured by heart rate variability. Eur J Clin Nutr. 2009;63(10):1263–5. https://doi.org/10.1038/ejcn.2009.35.
Neupane S, Leino-Arjas P, Nygård CH, Oakman J, Virtanen P. Developmental pathways of multisite musculoskeletal pain: what is the influence of physical and psychosocial working conditions? Occup Environ Med. 2017;74(7):468–75. https://doi.org/10.1136/oemed-2016-103892.
Neupane S, Virtanen P, Leino-Arjas P, Miranda H, Siukola A, Nygård CH. Multi-site pain and working conditions as predictors of work ability in a 4-year follow-up among food industry employees. Eur J Pain. 2013;17(3):444–51. https://doi.org/10.1002/j.1532-2149.2012.00198.x.
Pan F, Laslett L, Blizzard L, Cicuttini F, Winzenberg T, Ding C, Jones G. Associations between fat mass and multisite pain: a five-year longitudinal study. Arthritis Care Res. 2017;69(4):509–16. https://doi.org/10.1002/acr.22963.
Perini R, Veicsteinas A. Heart rate variability and autonomic activity at rest and during exercise in various physiological conditions. Eur J Appl Physiol. 2003;90(3–4):317–25. https://doi.org/10.1007/s00421-003-0953-9.
Phillips CJ. The cost and burden of chronic pain. Rev Pain. 2009;3(1):2–5. https://doi.org/10.1177/204946370900300102.
Sato TO, Hallman DM, Kristiansen J, Skotte JH, Holtermann A. Different autonomic responses to occupational and leisure time physical activities among blue-collar workers. Int Arch Occup Environ Health. 2017. https://doi.org/10.1007/s00420-017-1279-y.
Skotte J, Korshøj M, Kristiansen J, Hanisch C, Holtermann A. Detection of physical activity types using triaxial accelerometers. J Phys Act Health. 2014;11(1):76–84. https://doi.org/10.1123/jpah.2011-0347.
Skotte JH, Kristiansen J. Heart rate variability analysis using robust period detection. Biomed Eng Online. 2014. https://doi.org/10.1186/1475-925X-13-138.
Solidaki E, Chatzi L, Bitsios P, Markatzi I, Plana E, Castro F, Palmer K, Coggon D, Kogevinas M. Work-related and psychological determinants of multisite musculoskeletal pain. Scand J Work Environ Health. 2010;36(1):54–61.
Tak LM, Riese H, de Bock GH, Manoharan A, Kok IC, Rosmalen JG. As good as it gets? A meta-analysis and systematic review of methodological quality of heart rate variability studies in functional somatic disorders. Biol Psychol. 2009;82(2):101–10. https://doi.org/10.1016/j.biopsycho.2009.05.002.
Tracy LM, Ioannou L, Baker KS, Gibson SJ, Georgiou-Karistianis N, Giummarra MJ. Meta-analytic evidence for decreased heart rate variability in chronic pain implicating parasympathetic nervous system dysregulation. Pain. 2016;157(1):7–29. https://doi.org/10.1097/j.pain.0000000000000360.
Walker-Bone K, Reading I, Coggon D, Cooper C, Palmer KT. The anatomical pattern and determinants of pain in the neck and upper limbs: an epidemiologic study. Pain. 2004;109(1–2):45–51. https://doi.org/10.1016/j.pain.2004.01.008.
Woda A, Picard P, Dutheil F. Dysfunctional stress responses in chronic pain. Psychoneuroendocrinology. 2016;71:127–35. https://doi.org/10.1016/j.psyneuen.2016.05.017.
We would like to thank the entire DPhacto research team for the project planning, data collection, handling and management.
This study was conducted with the financial support from The Danish Work Environment Research Fund and São Paulo Research Foundation (FAPESP), São Paulo, Brazil (grant#2015/18310–1). The funding body had no participation in the design of the study, data collection, data analysis and interpretation and in writing the manuscript.
Availability of data and materials
Data cannot be made publicly available but may be available upon reasonable request. Any requests concerning the availability of the data should be directed to the project leader Andreas Holtermann, firstname.lastname@example.org, no/+ 45 39 16 53 52 at the National Research Centre for the Working Environment, Denmark.
Ethics approval and consent to participate
The study was approved by the Ethics Committee for the Capital Region of Denmark (journal number H-2-2012-011) and conducted in accordance with the Helsinki declaration. Written informed consent was obtained from all participants included in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Estimates, standard error and P value from the linear mixed models for heart rate variability indices showing the main effect of group, domain and the interaction (group × domain) in the crude and adjusted models using a strict definition of pain-free workers in DPhacto (n = 568). Table S2. Stratified analysis for age (< 50 years; ≥50 years). Estimates, standard error and P values from the linear mixed models for heart rate variability indices showing the main effect of group, domain and the interaction (group × domain) in the crude and adjusted models in DPhacto. Table S3. Stratified analysis for sex (male; female). Estimates, standard error and P values from the linear mixed models for heart rate variability indices showing the main effect of group, domain and the interaction (group × domain) in the crude and adjusted models in DPhacto. Figure S1. Original distribution of the HRV indices during work. Figure S2. Original distribution of the HRV indices during leisure. Figure S3. Original distribution of the HRV indices during sleep. (DOCX 1424 kb)
About this article
Cite this article
de Oliveira Sato, T., Hallman, D.M., Kristiansen, J. et al. The association between multisite musculoskeletal pain and cardiac autonomic modulation during work, leisure and sleep – a cross-sectional study. BMC Musculoskelet Disord 19, 405 (2018). https://doi.org/10.1186/s12891-018-2312-3