Skip to main content
  • Research article
  • Open access
  • Published:

Musculoskeletal pains and cardiovascular autonomic function in the general Northern Finnish population

Abstract

Background

Heart rate variability (HRV) and baroreflex sensitivity (BRS) measurements provide means for the objective assessment of cardiovascular autonomic function. As previous studies have associated chronic pain with abnormal autonomic function, we aimed to characterize the relationship between the number of musculoskeletal pain sites (NPS), pain intensity, and cardiovascular autonomic function among the population-based Northern Finland Birth Cohort 1966.

Methods

At the age of 46, cohort members self-reported their musculoskeletal pains (enabling the determination of NPS [0–8] and pain intensity [Numerical Rating Scale, NRS, 0–10]) and underwent clinical assessments of cardiovascular autonomic function in seated and standing positions (HRV variables: heart rate [HR] and root mean square of successive differences in beat-to-beat intervals [rMSSD] for the entire cohort; BRS variables: low-frequency systolic blood pressure variability [SBPV] and cross-spectral baroreflex sensitivity [BRS] for those attending the examination in Oulu, Finland). Extensive confounder data were also collected (body mass index, physical activity, smoking, Hopkins Symptom Checklist-25, comorbidities, and medications). The full samples included 4186 and 2031 individuals (HRV and BRS samples, respectively). Three subanalyses focused on individuals with intense and frequent pain, individuals with symptoms of depression and anxiety, and the relationship between pain intensity and autonomic parameters.

Results

Linear regression models showed varying associations between NPS, pain intensity, and cardiovascular autonomic parameters. However, after all adjustments NPS was only associated with one outcome among women (BRS, standing: beta = − 0.015, p = 0.048) and two among men (HR, seated: beta = − 0.902, p = 0.003; HR, standing: beta = − 0.843, p = 0.014). Pain intensity was not associated with any outcome after full adjustments. Significant sex*pain interactions were found in the data.

Conclusions

Our data suggest that musculoskeletal pain has, at most, a limited independent association with cardiovascular autonomic function. Future studies should carefully account for the potential confounders and sex interactions that this study revealed.

Peer Review reports

Background

The sensation of pain, mediated by the peripheral and central nervous system to the cerebral cortex where it is ultimately perceived, serves as an indispensable protective function in the human body [1]. However, it notoriously commonly converts to disabling chronic conditions such as low back pain, fibromyalgia, or other complex pain syndromes [2, 3]. Although much remains to be studied in terms of the mechanisms of pain [1, 2], the interconnections between the autonomic nervous system (ANS) and pain have long been recognized [4]. Being able to adjust autonomic functions such as heart rate (HR), the ANS contributes to maintaining physiological equilibrium under changing conditions. Importantly, the ANS also regulates and modulates the perception of pain [5, 6].

Non-invasive methods for assessing cardiovascular autonomic function include measurements of vagally-mediated heart rate variability (HRV) and spontaneous baroreflex sensitivity (BRS), in which depressed vagal activity and augmented sympathetic activity suggest dysregulation [6,7,8]. Abnormal HRV and BRS parameters, indicating impaired cardiovascular autonomic function, seem to increase the risk of sudden cardiac death [7,8,9,10,11]. Together with various lifestyle-related factors and medical conditions [12,13,14,15,16,17], chronic pain has been associated with abnormal ANS function [6, 18], although there are also contradictory findings [19]. These considerations establish the need for further data on the relationship between pain characteristics and ANS function, not solely in patients with severe pain conditions but also from a wider perspective, including the general population with multiple pains. Potentially, this knowledge could be exploited to develop feasible tools for objectively quantifying the individual pain experience, to enhance health care resource allocation and ultimately, to improve pain treatment.

In this study we investigated the association between musculoskeletal pains and cardiovascular autonomic function in a large unselected sample of Northern Finns. At the age of 46, the participants self-reported the location, frequency and intensity of their pains with the aid of a questionnaire, and attended recordings of HRV and BRS, allowing the analysis of cardiovascular autonomic function. Based on previous findings regarding individuals suffering from chronic pain [6, 18], we hypothesized that an increasing number of pain sites (NPS) and increasing pain intensity would be associated with cardiovascular autonomic dysregulation in our general population sample.

Methods

Initiation and progression of the cohort study

We used a large middle-aged sample from the prospective, population-based Northern Finland Birth Cohort 1966 (NFBC1966) as our study population. The cohort was initiated in 1965–1966, when pregnant women who lived in the Northern Finnish provinces of Oulu and Lapland and had expected dates of delivery in the year 1966 (Jan 1 to Dec 31) were recruited into the cohort [20]. The cohort base initially covered up to 96% of births in the area (12,068 mothers and 12,231 children). The main data collections have been arranged when the cohort members (i.e., those born into the cohort) have been 1, 14, 31, and 46 years old. Data have been collected through postal questionnaires and clinical examinations. Figure 1 presents the progression of NFBC1966.

Fig. 1
figure 1

Flow chart of the study. BMI Body mass index, BRS Baroreflex sensitivity, HRV Heart rate variability, HSCL-25 Hopkins Symptom Checklist-25, LTPA Leisure-time physical activity

In 2012–2014, when the cohort members were 46 years old, postal questionnaires were sent to all those eligible (i.e., those alive with known contact details; n = 10,321), the response rate being 66% (n = 6825). Subsequently, those who were living in Finland were invited to participate in clinical examinations at laboratory units which were set up across the country by the NFBC project, with the attendance rate of 57% (n = 5861). The postal questionnaires included a pain survey, and the clinical examination included recordings of cardiovascular autonomic function, as specified below. Those who attended the clinical examinations in the Oulu laboratory unit were assessed for both HRV and BRS, whereas those attending the examination at laboratories outside Oulu were assessed for HRV only. From those who had attended the clinical examination at the 46-year follow-up, we excluded 1675 individuals due to missing data in the HRV analyses, and 3830 individuals due to missing data in the BRS analyses. The sample sizes and origins of the data are specified in Fig. 1.

Assessment of musculoskeletal pains at age 46

Musculoskeletal pains were self-reported via a questionnaire that elicited the individual’s previous history of pain symptoms by the question: ‘Have you had any aches or pains in the following areas of your body within the last 12 months (yes/no)? 1) neck, 2) shoulders, 3) arms/elbows 4) wrists/hands 5) low back, 6) hips 7) knees, 8) ankles/feet’. The question was accompanied with an illustration defining these anatomical areas. Each positive response was followed by a subsequent question: ‘How often have you had aches or pains in this area during the last 12 months? 1) On 1–7 days, 2) On 8–30 days, 3) On more than 30 days but not daily, 4) Daily’. In addition, the participant was asked to estimate the overall intensity of their pains using a Numerical Rating Scale (NRS) from 0 (defined as ‘no pain’) to 10 (defined as ‘extremely severe pain’).

Based on the pain location data, we created a variable representing the number of anatomical sites at which the individual had experienced pain during the previous 12 months (NPS, 0–8). NPS is a widely used measure of musculoskeletal pain [21] that has been associated with, e.g., overall health and sleep quality in the general population [22]. In addition to the full study population, we were also interested in two subgroups: those with at least moderately intense and frequent pain (corresponding to an intensity of ≥5 on NRS [23, 24] and frequency of > 30 days/year), i.e., Subsample 1 and corresponding Subanalysis 1; and those with pain-related clinically relevant symptoms of depression and anxiety (corresponding to Hopkins Symptom Checklist-25 (HSCL-25) total mean score of ≥1.55 [25, 26]), i.e., Subsample 2 and corresponding Subanalysis 2 (Fig. 1). Subanalysis 3 evaluated the association between overall pain intensity (according to NRS) and autonomic function among the full study population.

Assessment of cardiovascular autonomic function at age 46

We assessed cardiovascular autonomic function by means of HRV and BRS parameters at the 46-year clinical examinations. The detailed methodology behind the acquisition, inspection and processing of the cardiovascular ANS data is provided in earlier publications [15,16,17] which also introduced the HRV and BRS variables used here. For the assessment of HRV, we used an HR monitor (RS800CX, Polar Electro Oy, Kempele, Finland) to record R-R intervals (RRi) to an accuracy of 1 ms. For the assessment of BRS, we recorded standard lead-II electrocardiograms (ECG; Cardiolife, Nihon Kohden, Tokyo, Japan) together with ventilation rate (VR; nasal temperature probe MLT415/D by ADInstruments, Bella Vista, New South Wales, Australia) and blood pressure (BP; finger plethysmography by Nexfin, BMEYE Medical Systems, Amsterdam, the Netherlands), with a 1000 Hz sampling frequency (PowerLab 8/35, ADInstruments). At the measurement laboratory, we first ensured that the participant was correctly positioned (sitting position) and that the measurement devices were properly adjusted. We then allowed one minute for the participant’s HR to stabilize, after which we obtained a three-minute recording of data in a sitting position. Then, the participant stood up and remained standing for another three minutes while we recorded the data a second time. We analyzed the first 150 s of sitting data and the last 150 s of standing data.

For the analysis of HRV, we processed the RRi data by visual inspection (Hearts 1.2, University of Oulu, Oulu, Finland), replaced short artefacts and ectopic beats with the local average, and deleted longer sequences of noise or ectopy (≥ 10 consecutive defective beats). Recordings containing ≥80% of adequate data were accepted for further processing. Eligible HRV data in both sitting and standing positions were available for 5473 individuals (96% of those 5679 who attended the recording). From the eligible HRV data we calculated mean HR (bpm) and the root mean square of successive differences in RRi (rMSSD, ms), which represent cardiac vagal modulation [27].

For the analysis of BRS, we used custom-made software (Biosignal Processing Team, University of Oulu, Oulu, Finland) to process the ECG, BP and VR data. From the continuous ECG and BP recordings, we extracted systolic blood pressure (SBP) and RRi values, respectively. We replaced artefacts and ectopic beats using linear interpolation and then resampled at 2 Hz. Recordings with < 5% of defective data were accepted for further processing. We used the Savitzky-Golay method to delete very low frequency components (< 0.04 Hz) from the data, and we conducted fast Fourier transform (Welch method, segments of 128 samples with 50% overlap, length 1024 samples) to analyze the low frequency power of RRi and SBP oscillations (ms2, mmHg2) for the subsequent analysis of BRS by the alpha method, after ensuring a coherence of ≥0.5 between the low frequency oscillations in RRi and SBP [28]. Eligible BRS data in both sitting and standing positions were available for 2641 and 2617 individuals (97 and 96% of the 2726 who attended the recording), respectively. Low frequency SBP variability (SBPV, mmHg2), an estimate of peripheral sympathetic regulation [29], and cross-spectral BRS (ms/mmHg) are principal variables representing ANS function obtained by this protocol.

Confounders

We assessed the participants’ body mass index (BMI) [17], leisure-time physical activity (LTPA) [15, 17], smoking [30], and symptoms of depression and anxiety [31] as potential confounders, together with comorbidities and medications influencing ANS function. BMI (kg/m2) was calculated on the basis of systematic height and weight measurements, which were obtained in the 46-year clinical examination. The frequency of LTPA was elicited in the 46-year questionnaires by asking: ‘How often do you participate in brisk physical activity (defined as causing at least some sweating and getting out of breath) during your leisure-time? 1) daily, 2) 4–6 times a week, 3) 2–3 times a week, 4) once a week, 5) 2–3 times a month, 6) once a month or less often’. We combined the first two and last two categories due to small group sizes. Lifetime smoking history was elicited in the 46-year questionnaire by the questions: ‘Have you ever smoked cigarettes (yes/no)?’ and ‘Do you currently smoke (yes/no)?’ As previously described [32], we created three categories according to the responses: non-smoker, former smoker and current smoker. Symptoms of depression and anxiety were assessed in the 46-year questionnaire using the acknowledged HSCL-25 tool [33, 34]. We used the total HSCL-25 mean score in the analysis. The participants were also asked to report all their medications, dietary supplements, and diseases diagnosed by a physician in the 46-year questionnaire. The routine use of antihypertensives and/or painkillers was regarded as confounding medication, while confounding comorbidities included diabetes mellitus, cardiovascular diseases and fibromyalgia. Data on diabetes mellitus were based not only on self-reports but also on a two-hour oral glucose tolerance test and the measurement of glycated hemoglobin in the 46-year follow-up, as described in our previous publication [16]. Age was not assessed as a confounder because the sample was coeval.

Statistical analysis

We explored the characteristics of the sample by calculating the frequencies and percentages of the categorical variables, means and standard deviations (SD) of continuous variables with normal distributions, and the medians and interquartile ranges (IQR) of continuous variables with skewed distributions. NPS distribution was illustrated using histograms. We calculated the correlation between NPS and overall pain intensity using Spearman’s Rho (R).

We analyzed the association between musculoskeletal pain and cardiovascular autonomic parameters using linear regression models. The autonomic parameters served as independent outcomes, for which all models were run. The parameters which had skewed distributions were natural logarithm-transformed to Gaussian, as specified in Additional file 1. NPS was the primary predictor in all models except for Subanalysis 3 in which pain intensity (NRS) was used instead of NPS. Four sex-stratified models were run for each outcome: Model I (crude); Model II (adjusted for lifestyle); Model III (adjusted for comorbidities); and Model IV (adjusted for lifestyle and comorbidities). The variables used in the models are specified in the Results section; we have also summarized the outcomes, explanatory variables, variable coding, and regression models in Additional files 1, 2 and 3. Beta coefficients (β) of the predictors were collected with their 95% confidence intervals (CI). In addition to the Primary Analysis, we ran three subanalyses; the full sample and subsamples have been specified in the earlier part of the Methods section and in Fig. 1, and the subanalyses are further introduced in the Results section.

In addition to continuous modeling, the NPS data were also assessed using various categorizations (0,1,2,3,4,5,6,7,8; 0,1,2−3,4−5,6–8; 0,1−2,3−5,6–8) in order to rule out potential non-linear relationship patterns in the data and thus justify the linear modeling of NPS. The need for sex stratification was explored by running all models with pooled sexes and including the sex*NPS (or sex*NRS) interaction term in the models. The need for sex stratification was evaluated on the basis of the statistical significance of the interaction terms, taking into account the large sample size which would allow stratification without critically reducing the statistical power of the analyses.

The data were analyzed using SPSS version 25 (IBM, Armonk, NY, USA). P values of ≤ 0.05 were considered statistically significant.

Ethical considerations

The study protocol followed the Declaration of Helsinki and was approved by the Ethical Committee of the Northern Ostrobothnia Hospital District in accordance with the Declaration of the World Medical Association. Each participant granted their written informed consent.

Results

The study sample consisted of 1813 men and 2373 women, whose detailed characteristics are presented in Table 1. The median NPS was 3 (interquartile range 2–5) among men and 4 (2–5) among women. Figure 2 presents the distribution of NPS among the sample, together with the mean reported intensity of pain in each NPS category. The overall intensity of pain (NRS) correlated positively with NPS (R = 0.382 among men and 0.430 among women, p < 0.001).

Table 1 Characteristics of the full sample, i.e., those with pain data, covariate data, and either HRV or BRS data (n = 4186)
Fig. 2
figure 2

Distribution of NPS among men and women, and mean intensity of pain in NPS categories. NPS Number of pain sites, NRS Numerical Rating Scale (0–10)

Non-linear relationships between NPS and cardiovascular autonomic parameters were not detected. However, the pooled-sex analyses showed statistically significant sex*NPS and sex*NRS interactions (p < 0.05, data not shown), and all models were consequently stratified by sex. The regression coefficients for the associations between NPS and cardiovascular autonomic parameters (i.e., HR, rMSSD, SBPV, BRS) are presented in Table 2 (Primary Analysis), Table 3 (Subanalysis 1), and Table 4 (Subanalysis 2). Correspondingly, the regression coefficients for the associations between pain intensity (NRS) and autonomic parameters are shown in Table 5 (Subanalysis 3). For clarity, we only provide the regression coefficients of the primary predictors (i.e., NPS/NRS) in the main manuscript and present the full models in Additional files 4, 5, 6, 7, 8, 9, 10 and 11.

Table 2 Primary Analysis. Results from sex-stratified linear regression models regarding the association between musculoskeletal pain sites (NPS) and cardiovascular autonomic function (HR, rMSSD, SBPV, BRS) in the full sample (for HR and rMSSD, n = 1813 men and 2373 women; for SBPV and BRS, n = 919 men and 1112 women). The complete models are presented in Additional files 4 and 5
Table 3 Subanalysis 1. Results from sex-stratified linear regression models regarding the association between musculoskeletal pain sites (NPS) and cardiovascular autonomic function (HR, rMSSD, SBPV, BRS) in Subsample 1, i.e., individuals with intense (NRS ≥ 5) and frequent (> 30 days/year) pain (for HR and rMSSD, n = 450 men and 761 women; for SBPV and BRS, n = 250 men and 375 women). The complete models are presented in Additional files 6 and 7
Table 4 Subanalysis 2. Results from sex-stratified linear regression models regarding the association between musculoskeletal pain sites (NPS) and cardiovascular autonomic function (HR, rMSSD, SBPV, BRS) in Subsample 2, i.e., individuals with clinically relevant symptoms of depression and anxiety (HSCL-25 score ≥ 1.55) (for HR and rMSSD, n = 302 men and 507 women; for SBPV and BRS, n = 163 men and 241 women). The complete models are presented in Additional files 8 and 9
Table 5 Subanalysis 3. Results from sex-stratified linear regression models regarding the association between between musculoskeletal pain intensity (according to NRS) and cardiovascular autonomic function (HR, rMSSD, SBPV, BRS) in the full sample (for HR and rMSSD, n = 1619 men and 2096 women; for SBPV and BRS, n = 832 men and 979 women). The complete models are presented in Additional files 10 and 11

Among women, NPS showed several statistically significant associations with cardiovascular autonomic parameters in the crude and partially adjusted models of the Primary Analysis, Subanalysis 1 and Subanalysis 2, but these associations were mostly attenuated by adjustments (Tables 24). Only in the Primary Analysis was one outcome (i.e., BRS, standing) consistently associated with NPS after full adjustments (β = − 0.015, p = 0.048). In Subanalysis 3, pain intensity showed several statistically significant crude associations with the autonomic parameters, but all of these became non-significant after full adjustments (Table 5).

Among men, some statistically significant associations between NPS and autonomic parameters were obtained from Subanalysis 1 and Subanalysis 2 (Tables 24), but only in Subanalysis 1 were two outcomes (i.e., HR, seated and standing) associated with NPS after full adjustments (HR, seated: β = − 0.902, p = 0.003; HR, standing: β = − 0.843, p = 0.014). In Subanalysis 3, pain intensity showed several statistically significant associations with the autonomic parameters in the crude and lifestyle-adjusted models but not after full adjustments (Table 5).

Generally, lifestyle variables (i.e., BMI, LTPA and smoking) and other covariates (i.e., HSCL-25, comorbidity, medication) showed stronger and more consistent associations with cardiovascular autonomic parameters than NPS or pain intensity among both sexes (Additional files 4, 5, 6, 7, 8, 9, 10 and 11).

Discussion

The present study was a cross-sectional birth cohort study that assessed the association between musculoskeletal pains (in terms of NPS and NRS) and cardiovascular autonomic function (in terms of HR, rMSSD, SBPV, and BRS) in a large unselected sample of Northern Finns. The analyses showed some associations between NPS, pain intensity, and altered cardiovascular autonomic function, but these associations were mostly attenuated by adjustments. Thus, the results of the present study suggest that musculoskeletal pain has, at most, a limited independent association with cardiovascular autonomic function among the general Northern Finnish population.

Previous studies have investigated several groups of patients suffering from, e.g., low back pain [35, 36] and chronic pain syndromes [37], but surprisingly few large studies [18, 19] have addressed the general population with multiple pains. Importantly, results derived from large unselected samples should be of great relevance due to the high prevalence of musculoskeletal pains in the general population [38] and the subsequent economic burden on health care systems and societies [39]. Improved pain quantification methods would enhance resource allocation if they were applicable to a wide population base. As such, the present study aimed to explore the association between musculoskeletal pains and cardiovascular autonomic parameters in a large general population sample of Northern Finns.

In our data, some of the crude and partially adjusted models showed increasing NPS and pain intensity to be significantly associated with higher HR, lower rMSSD, and lower BRS, all of which point towards impaired cardiovascular autonomic function with increasing NPS and pain intensity. These findings were varyingly obtained from both sexes, such that the association between NPS and autonomic function was somewhat more consistent among women and the association between pain intensity and autonomic function among men. However, after full adjustment for lifestyle factors (BMI, LTPA, smoking) and comorbidities (HSCL-25, somatic comorbidity, medication), the associations were mostly attenuated, with only a few exceptions. The pooled-sex analyses showed significant sex*pain interactions to be related to the cardiovascular autonomic parameters.

The present results suggest that NPS is not independently associated with cardiovascular autonomic function among the general middle-aged Northern Finnish population (Primary Analysis), not even when the pain is intense and frequent (intensity ≥5 on NRS and frequency > 30 days/year; Subanalysis 1), or when it is associated with clinically relevant symptoms of depression and anxiety (HSCL-25 score ≥ 1.55; Subanalysis 2). Similarly, the overall intensity of musculoskeletal pain seems to lack an independent association with cardiovascular autonomic function in this study population. Thus, the present findings evoke two principal research questions for future studies. First, given that the confounder variables attenuated the crude associations of NPS and pain intensity with cardiovascular autonomic parameters, and also showed more systematic associations with the outcomes than NPS or pain intensity did, the question arises whether these factors may mediate the association between pain and cardiovascular autonomic function. For example, obesity and physical activity have both been associated with experiencing pain [40, 41] and ANS function [17, 42], and were systematically associated with the autonomic parameters in the present data. Due to the limitations of our study setting, future studies should aim to investigate the potential causal relationships between pain and ANS function, together with an in-depth characterization of potential mediators. Second, our results regarding sex*pain interactions suggest that the association between pain and cardiovascular autonomic function differs according to sex. This finding has also been reported in an earlier study [43], but further investigations of other datasets are needed to confirm it. In the meantime, it seems advisable to carefully assess and account for potential sex*pain interactions in future studies.

Two previous studies have investigated the association between pain and ANS parameters in large population-based (or partially population-based) samples. In a recent Norwegian study [18], 1143 individuals who suffered clinically meaningful chronic pain were compared to 5640 controls without pain in terms of several HRV and BRS parameters. After adjustment for age, sex and BMI, chronic pain was associated with reduced HRV and BRS, indicating impaired cardiovascular autonomic function. In addition to the differences between the pain variables (clinically meaningful chronic pain vs. calculation of NPS) and the modeling technique (case-control vs. continuous modeling of pain sites), also the adjustments of the Norwegian study and the present study varied considerably (age, sex, BMI vs. sex-stratification, BMI, LTPA, smoking, HSCL-25, comorbidity, medication). In the present regression models, most confounders showed clear associations with the autonomic parameters and were thus seen as relevant covariates. Similar findings regarding covariates were also reported in the Norwegian study, although it assessed fewer covariates than the present study. The present data also manifested significant sex*pain interactions which, interestingly, were tested for but not detected in the Norwegian data. It remains unclear whether the differing findings reflect differences between the studies’ definition of pain (presence of clinically meaningful chronic pain vs. NPS), comprehensiveness of the covariates, or other methodological factors.

In contrast to the Norwegian study, a Dutch study [19] found evidence for the association between widespread chronic pain and ANS dysregulation in their partially population-based sample, which consisted of both healthy individuals and patients with depressive and anxiety disorders (731 individuals reporting chronic widespread pain and 843 controls without pain). Notably, the assessed confounders were extensive (age, sex, education, alcohol use, BMI, smoking, medications, somatic and psychiatric comorbidities, physical activity, insomnia). Similar to the present study, the Dutch study detected several significant associations between pain and the ANS variables in the crude models, but not after controlling for the confounders. The authors stated that previous studies have been potentially hampered by insufficient control of potential confounders, small sample sizes and the use of patient samples instead of population-based approaches. Prospective, longitudinal studies utilizing large representative samples are needed to further address these considerations.

The main strength of the present study is its large general population sample with extensive data on pain, lifestyle and health, allowing the assessment of NPS, pain intensity, and cardiovascular autonomic function, while accounting for numerous confounding factors. In comparison with patient samples, the population-based setting of the present study increased the generalizability of its results. Furthermore, the present study was designed to address musculoskeletal pains at several anatomical locations instead of focusing on a certain pain site, indicating that the present approach was comprehensive, and took into account the possible widespread nature of pain. The present study was also able to assess several cardiovascular autonomic parameters, addressing both HRV and BRS.

A major limitation of this study is its cross-sectional setting. The outcomes (i.e., cardiovascular autonomic parameters) and primary predictors (i.e., NPS/NRS) were only recorded at the most recent follow-up of the cohort study, i.e., at the age of 46. Despite the large, well-characterized sample, the lack of longitudinal data on pain may have reduced the strength of the variable and may thus partially explain the lack of consistent association between pain and cardiovascular autonomic function. The pain data may also be limited because the pain survey enquired musculoskeletal pains over a period of 12 months. This may have increased the likelihood of recall bias, and also decreased the temporal accordance between the exposure (i.e., pain) and outcome (i.e., cardiovascular autonomic function). Unfortunately, the study subjects’ musculoskeletal pain status was not enquired on the clinical examination day when the autonomic parameters were obtained due to NFBC1966’s data collection protocol. However, it is notable that many crude models demonstrated clear associations between NPS, pain intensity, and cardiovascular autonomic parameters, and these associations were only attenuated after adjustments. As such, these findings suggest that the lack of association is primarily explained by the presence of confounding factors and not by the limitations of the pain data themselves. In the calculation of NPS, each anatomical location was weighted equally due to the lack of more accurate calculation protocols. However, NPS is a widely used measure of musculoskeletal pain [21] and has been associated with, e.g., overall health and sleep quality in the general population [22]. The linear relationship between NPS and cardiovascular autonomic parameters was only assumed after no other relationship patterns were detected. With regard to the outcomes, we acknowledge that our ANS parameters were indicators of ANS-based cardiovascular function, as opposed to other aspects of ANS function. Thus, the present conclusions primarily concern cardiovascular autonomic function, leaving the more peripheral ANS to be explored in future studies.

Conclusions

This cross-sectional population-based study of Northern Finns found some evidence of a relationship between NPS, pain intensity, and cardiovascular autonomic dysfunction among both sexes. However, these associations were mostly attenuated after adjustment for lifestyle factors and comorbidities. Thus, the present results suggest that musculoskeletal pain has, at most, a limited independent association with cardiovascular autonomic function in the general Northern Finnish population, even when the pain is intense and frequent, or associated with clinically relevant symptoms of depression and anxiety. Future studies should confirm the present findings in other general population samples with longitudinal designs, and take into account the potential confounders and sex interactions that this study revealed. If a causal relationship between musculoskeletal pain and ANS is detected, the potential mediators should also be characterized.

Abbreviations

ANS:

Autonomic nervous system

BMI:

Body mass index

BRS:

Baroreflex sensitivity

CI:

Confidence interval

HR:

Heart rate

HRV:

Heart rate variability

HSCL-25:

Hopkins Symptom Checklist-25

LTPA:

Leisure-time physical activity

NPS:

Number of pain sites

NRS:

Numerical Rating Scale

P:

P value

rMSSD:

Root mean square of the successive differences in R-R intervals

SBPV:

Systolic blood pressure variability

β:

Beta coefficient

References

  1. Kuner R. Central mechanisms of pathological pain. Nat Med. 2010;16:1258–66.

    Article  CAS  Google Scholar 

  2. Diatchenko L, Fillingim R, Smith S, Maixner W. The phenotypic and genetic signatures of common musculoskeletal pain conditions. Nat Rev Rheumatol. 2013;9:340–50.

    Article  CAS  Google Scholar 

  3. Marinus J, Moseley G, Birklein F, Baron R, Maihofner C, Kingery W, et al. Clinical features and pathophysiology of complex regional pain syndrome. Lancet Neurol. 2011;10:637–48.

    Article  Google Scholar 

  4. Haugen F. The autonomic nervous system and pain. Anesthesiology. 1968;29:785–92.

    Article  CAS  Google Scholar 

  5. Bruehl S, Chung O. Interactions between the cardiovascular and pain regulatory systems: an updated review of mechanisms and possible alterations in chronic pain. Neurosci Biobehav Rev. 2004;28:395–414.

    Article  Google Scholar 

  6. Tracy L, Ioannou L, Baker K, Gibson S, Georgiou-Karistianis N, Giummarra M. Meta-analytic evidence for decreased heart rate variability in chronic pain implicating parasympathetic nervous system dysregulation. Pain. 2016;157:7–29.

    Article  Google Scholar 

  7. Vaseghi M, Shivkumar K. The role of the autonomic nervous system in sudden cardiac death. Prog Cardiovasc Dis. 2008;50:404–19.

    Article  CAS  Google Scholar 

  8. Fukuda K, Kanazawa H, Aizawa Y, Ardell J, Shivkumar K. Cardiac innervation and sudden cardiac death. Circ Res. 2015;116:2005–19.

    Article  CAS  Google Scholar 

  9. La Rovere M, Bigger J, Marcus F, Mortara A, Schwartz P. Baroreflex sensitivity and heart-rate variability in prediction of total cardiac mortality after myocardial infarction. ATRAMI (autonomic tone and reflexes after myocardial infarction) investigators. Lancet. 1998;351:478–84.

    Article  Google Scholar 

  10. Kiviniemi A, Tulppo M, Hautala A, Perkiomaki J, Ylitalo A, Kesaniemi Y, et al. Prognostic significance of impaired baroreflex sensitivity assessed from phase IV of the valsalva maneuver in a population-based sample of middle-aged subjects. Am J Cardiol. 2014;114:571–6.

    Article  Google Scholar 

  11. Shen M, Zipes D. Role of the autonomic nervous system in modulating cardiac arrhythmias. Circ Res. 2014;114:1004–21.

    Article  CAS  Google Scholar 

  12. Tsuji H, Venditti F, Manders E, Evans J, Larson M, Feldman C, et al. Determinants of heart rate variability. J Am Coll Cardiol. 1996;28:1539–46.

    Article  CAS  Google Scholar 

  13. Kardos A, Watterich G, de Menezes R, Csanady M, Casadei B, Rudas L. Determinants of spontaneous baroreflex sensitivity in a healthy working population. Hypertension. 2001;37:911–6.

    Article  CAS  Google Scholar 

  14. Kiviniemi A, Hautala A, Kinnunen H, Nissila J, Virtanen P, Karjalainen J, et al. Daily exercise prescription on the basis of HR variability among men and women. Med Sci Sports Exerc. 2010;42:1355–63.

    Article  Google Scholar 

  15. Kiviniemi A, Perkiomaki N, Auvinen J, Herrala S, Hautala A, Ahola R, et al. Lifelong physical activity and cardiovascular autonomic function in midlife. Med Sci Sports Exerc. 2016;48:1506–13.

    Article  Google Scholar 

  16. Perkiomaki N, Auvinen J, Tulppo M, Hautala A, Perkiomaki J, Karhunen V, et al. Association between birth characteristics and cardiovascular autonomic function at mid-life. PLoS One. 2016;11:e0161604.

    Article  Google Scholar 

  17. Kiviniemi A, Perkiomaki N, Auvinen J, Niemela M, Tammelin T, Puukka K, et al. Fitness, fatness, physical activity, and autonomic function in midlife. Med Sci Sports Exerc. 2017;49:2459–68.

    Article  Google Scholar 

  18. Bruehl S, Olsen R, Tronstad C, Sevre K, Burns J, Schirmer H, et al. Chronic pain-related changes in cardiovascular regulation and impact on comorbid hypertension in a general population: the tromso study. Pain. 2018;159:119–27.

    Article  Google Scholar 

  19. Barakat A, Vogelzangs N, Licht C, Geenen R, MacFarlane G, de Geus E, et al. Dysregulation of the autonomic nervous system and its association with the presence and intensity of chronic widespread pain. Arthritis Care Res (Hoboken). 2012;64:1209–16.

    Google Scholar 

  20. Rantakallio P. The longitudinal study of the northern Finland birth cohort of 1966. Paediatr Perinat Epidemiol. 1988;2:59–88.

    Article  CAS  Google Scholar 

  21. Blyth F, Noguchi N. Chronic musculoskeletal pain and its impact on older people. Best Pract Res Clin Rheumatol. 2017;31:160–8.

    Article  Google Scholar 

  22. Kamaleri Y, Natvig B, Ihlebaek C, Benth J, Bruusgaard D. Number of pain sites is associated with demographic, lifestyle, and health-related factors in the general population. Eur J Pain. 2008;12:742–8.

    Article  Google Scholar 

  23. Oldenmenger W, de Raaf P, de Klerk C, van der Rijt C. Cut points on 0-10 numeric rating scales for symptoms included in the Edmonton symptom assessment scale in cancer patients: a systematic review. J Pain Symptom Manag. 2013;45:1083–93.

    Article  Google Scholar 

  24. Boonstra A, Schiphorst Preuper H, Balk G, Stewart R. Cut-off points for mild, moderate, and severe pain on the visual analogue scale for pain in patients with chronic musculoskeletal pain. Pain. 2014;155:2545–50.

    Article  Google Scholar 

  25. Veijola J, Jokelainen J, Laksy K, Kantojarvi L, Kokkonen P, Jarvelin M, et al. The Hopkins symptom checklist-25 in screening DSM-III-R axis-I disorders. Nord J Psychiatry. 2003;57:119–23.

    Article  Google Scholar 

  26. Nettelbladt P, Hansson L, Stefansson C, Borgquist L, Nordstrom G. Test characteristics of the Hopkins symptom check list-25 (HSCL-25) in Sweden, using the present state examination (PSE-9) as a caseness criterion. Soc Psychiatry Psychiatr Epidemiol. 1993;28:130–3.

    Article  CAS  Google Scholar 

  27. Heart rate variability. Standards of measurement, physiological interpretation and clinical use. Task force of the european society of cardiology and the north american society of pacing and electrophysiology. Circulation. 1996;93:1043–65.

    Article  Google Scholar 

  28. Pagani M, Somers V, Furlan R, Dell'Orto S, Conway J, Baselli G, et al. Changes in autonomic regulation induced by physical training in mild hypertension. Hypertension. 1988;12:600–10.

    Article  CAS  Google Scholar 

  29. Furlan R, Porta A, Costa F, Tank J, Baker L, Schiavi R, et al. Oscillatory patterns in sympathetic neural discharge and cardiovascular variables during orthostatic stimulus. Circulation. 2000;101:886–92.

    Article  CAS  Google Scholar 

  30. Middlekauff H, Park J, Moheimani R. Adverse effects of cigarette and noncigarette smoke exposure on the autonomic nervous system: Mechanisms and implications for cardiovascular risk. J Am Coll Cardiol. 2014;64:1740–50.

    Article  CAS  Google Scholar 

  31. Kemp A, Quintana D. The relationship between mental and physical health: insights from the study of heart rate variability. Int J Psychophysiol. 2013;89:288–96.

    Article  Google Scholar 

  32. Oura P, Paananen M, Niinimaki J, Tammelin T, Herrala S, Auvinen J, et al. Effects of leisure-time physical activity on vertebral dimensions in the northern Finland birth cohort 1966. Sci Rep. 2016;6:27844.

    Article  CAS  Google Scholar 

  33. Winokur A, Winokur D, Rickels K, Cox D. Symptoms of emotional distress in a family planning service: stability over a four-week period. Br J Psychiatry. 1984;144:395–9.

    Article  CAS  Google Scholar 

  34. Nabbe P, Le Reste J, Guillou-Landreat M, Munoz Perez M, Argyriadou S, Claveria A, et al. Which DSM validated tools for diagnosing depression are usable in primary care research? A systematic literature review. Eur Psychiatry. 2017;39:99–105.

    Article  CAS  Google Scholar 

  35. Sodervall J, Karppinen J, Puolitaival J, Kyllonen E, Kiviniemi A, Tulppo M, et al. Heart rate variability in sciatica patients referred to spine surgery: a case control study. BMC Musculoskelet Disord. 2013;14:149.

    Article  Google Scholar 

  36. Hautala A, Karppinen J, Seppanen T. Short-term assessment of autonomic nervous system as a potential tool to quantify pain experience. Conf Proc IEEE Eng Med Biol Soc. 2016;2016:2684–7.

    PubMed  Google Scholar 

  37. Meeus M, Goubert D, De Backer F, Struyf F, Hermans L, Coppieters I, et al. Heart rate variability in patients with fibromyalgia and patients with chronic fatigue syndrome: a systematic review. Semin Arthritis Rheum. 2013;43:279–87.

    Article  Google Scholar 

  38. Cimmino M, Ferrone C, Cutolo M. Epidemiology of chronic musculoskeletal pain. Best Pract Res Clin Rheumatol. 2011;25:173–83.

    Article  Google Scholar 

  39. Henschke N, Kamper S, Maher C. The epidemiology and economic consequences of pain. Mayo Clin Proc. 2015;90:139–47.

    Article  Google Scholar 

  40. Naugle K, Riley J. Self-reported physical activity predicts pain inhibitory and facilitatory function. Med Sci Sports Exerc. 2014;46:622–9.

    Article  Google Scholar 

  41. Okifuji A, Hare B. The association between chronic pain and obesity. J Pain Res. 2015;8:399–408.

    Article  Google Scholar 

  42. Guarino D, Nannipieri M, Iervasi G, Taddei S, Bruno R. The role of the autonomic nervous system in the pathophysiology of obesity. Front Physiol. 2017;8:665.

    Article  Google Scholar 

  43. Walker L, Stone A, Smith C, Bruehl S, Garber J, Puzanovova M, et al. Interacting influences of gender and chronic pain status on parasympathetically mediated heart rate variability in adolescents and young adults. Pain. 2017;158:1509–16.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the late professor Paula Rantakallio (launch of NFBC1966), the participants in the study and the NFBC project center.

Funding

The funds of Tekes – the Finnish Funding Agency for Technology and Innovation were received in support of this work. NFBC1966 received financial support from University of Oulu (Grant no. 24000692), Oulu University Hospital (Grant no. 24301140), and ERDF European Regional Development Fund (Grant no. 539/2010 A31592). A.K. received financial support from The Paulo Foundation (Espoo, Finland) and the Finnish Foundation for Cardiovascular Research (Helsinki, Finland). The funding sources had no involvement in study design, data collection, analysis, interpretation, or writing the paper.

Availability of data and materials

The datasets generated and analysed during the current study are not publicly available due to local privacy regulations but are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

AH, AK, JA, KP, MT, HH, TS, and JK collected the data. PO analysed the data, wrote the manuscript draft, and prepared the figures. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Petteri Oura.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Ethical Committee of the Northern Ostrobothnia Hospital District. Each participant granted their written informed consent.

Consent for publication

Not applicable.

Competing interests

Funding has been disclosed above (A.K.). The remaining authors declare no conflict of interest.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Summary of outcome variables. (DOCX 27 kb)

Additional file 2:

Predictor and covariate variables of the regression models. (DOCX 24 kb)

Additional file 3:

Construction of regression models. (DOCX 27 kb)

Additional file 4:

Primary Analysis, men. (DOCX 52 kb)

Additional file 5:

Primary Analysis, women. (DOCX 51 kb)

Additional file 6:

Subanalysis 1, men. (DOCX 50 kb)

Additional file 7:

Subanalysis 1, women. (DOCX 50 kb)

Additional file 8:

Subanalysis 2, men. (DOCX 48 kb)

Additional file 9:

Subanalysis 2, women. (DOCX 48. kb)

Additional file 10:

Subanalysis 3, men. (DOCX 52 kb)

Additional file 11:

Subanalysis 3, women. (DOCX 52 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Oura, P., Hautala, A., Kiviniemi, A. et al. Musculoskeletal pains and cardiovascular autonomic function in the general Northern Finnish population. BMC Musculoskelet Disord 20, 45 (2019). https://doi.org/10.1186/s12891-019-2426-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-019-2426-2

Keywords