- Open Access
Heart rate and heart rate variability in patients with chronic inflammatory joint disease: the role of pain duration and the insular cortex
BMC Musculoskeletal Disorders volume 23, Article number: 75 (2022)
Chronic inflammatory joint diseases (CIJD) have been linked to increased cardiovascular morbidity and mortality. A decisive reason could be a dysregulation of the autonomic nervous system, which is responsible for the control of cardiovascular function. So far, the cause of changes in autonomic nervous system functions remains elusive. In this study, we investigate the role of chronic pain and the insular cortex in autonomic control of cardiac functioning in patients with CIJD.
We studied the autonomic nervous system through the assessment of heart rate and heart rate variability (HRV) at rest and under cognitive stimulation. Furthermore, we investigated insular cortex volume by performing surface-based brain morphometry with FreeSurfer. For this study, 47 participants were recruited, 22 individual age- and sex-matched pairs for the magnetic resonance imaging analyses and 14 for the HRV analyses. All available patients’ data were used for analysis.
Pain duration was negatively correlated with the resting heart rate in patients with chronic inflammatory joint diseases (n = 20). In a multiple linear regression model including only CIJD patients with heart rate at rest as a dependent variable, we found a significant positive relationship between heart rate at rest and the volume of the left insular cortex and a significant negative relationship between heart rate at rest and the volume of the right insular cortex. However, we found no significant differences in HRV parameters or insular cortex volumes between both groups.
In this study we provide evidence to suggest insular cortex involvement in the process of ANS changes due to chronic pain in CIJD patients.
The study was preregistered with the German Clinical Trials Register (https://www.drks.de; DRKS00012791; date of registration: 28 July 2017).
Chronic inflammatory joint diseases (CIJD), such as rheumatoid arthritis (RA) and ankylosing spondylitis, are commonly associated with disability and pain . CIJDs have also been linked to increased cardiovascular morbidity and mortality, which cannot be attributed to typical risk factors, such as smoking or hypertension [2,3,4,5,6]. Instead, a dysregulation of the autonomic nervous system (ANS) has been suggested as a potential underlying mechanism. This notion is based on an estimate of approximately 60% of patients with RA showing signs of ANS dysfunction .
The function of the ANS can be assessed by measuring heart rate and calculating heart rate variability (HRV), i.e. the variation of the length of the R-to-R-wave time interval [8, 9]. In patients with rheumatic diseases, a reduction of HRV has been observed, indicating reduced cardiac adaptability [7, 10, 11].
However, the reasons for these changes in ANS function in CIJD remain elusive. In 2006, Dunn and Croft for instance showed the importance of pain duration in patients suffering from chronic low back pain. According to their results, the prognosis of chronic back pain was dependent on pain duration. It therefore seems feasible to hypothesize that pain duration might also impact bodily systems like the ANS, represented by changes in HRV .
The ANS is controlled by centers in the brainstem, the hypothalamus, and, ultimately, by several cortical regions, including the insula, anterior cingulate gyrus, and medial frontal cortex [13, 14]. The insula is considered a crucial region for polymodal sensory, cognitive, and affective processing . In addition to the processing function of somatosensory [16, 17], peripheral inflammatory , and nociceptive stimuli , the insular cortex participates in the control of arbitrary and involuntary movements, such as speech production  and swallowing . Oppenheimer et al. were first to demonstrate changes in cardiac function after electrical stimulation of the left and right insula in humans . Corroborating these results, several studies reported cardiac dysfunction following insula lesions due to stroke, tumors, or surgical resections [23,24,25,26].
Wartolowska and colleagues have investigated morphological changes in the brains of patients with rheumatic diseases. Although they did not investigate cortical regions, they described morphological changes in subcortical areas, including the caudate nucleus, putamen, and nucleus accumbens . Similarly to what has been described for other chronic pain conditions such as chronic back and limb pain [28, 29], complex regional pain syndrome , migraine , chronic tension type-headache , and fibromyalgia , patients with CIJD might also display structural changes in cortical regions that are involved in the regulation of the ANS, especially in the insular cortex.
The aims of the present study hence were a) to evaluate differences in HRV at rest and under stimulation between the CIJD and control groups, b) to test whether pain duration might be associated with changes in heart rate and heart rate variability in patients with chronic inflammatory joint diseases, and c) whether altered HRV is associated with structural changes in ANS-controlling cortical areas such as the insular cortex.
The study was approved by the Medical Research Ethics Board of Carl von Ossietzky University of Oldenburg, Germany (#2017–059) and was preregistered with the German Clinical Trials Register (https://www.drks.de; DRKS00012791). All participants gave written informed consent prior to participating in this study. Participants (healthy controls: C, patients with chronic inflammatory joint diseases: CIJD) were recruited between July 2017 and March 2019. All measurements were performed at the Neuroimaging Unit, School of Medicine and Health Sciences, Carl von Ossietzky University of Oldenburg. All data were collected as a part of a larger study, which, in addition to the data evaluated here, included functional magnetic resonance imaging (MRI) and behavioral tests [34, 35].
Patients aged between 30 and 68 years were recruited at a specialized rheumatological outpatient clinic and with the help of a support group for ankylosing spondylitis patients in Oldenburg. Controls were recruited by advertisements in a local newspaper and announcements on the university’s website.
Participants had to be fluent in German and right-handed. Patients needed to suffer from a chronic inflammatory joint disease that caused chronic pain for at least 1 year while controls had to be free of CIJD and chronic pain in general. Controls and CIJD were matched for sex and age (± 5 years).
Pregnant women and participants with contraindications for MRI were excluded, in addition to participants with conditions that could have had an influence on HRV or cerebral functions: previous heart surgeries; heart failure stage > 2 of the New York Heart Association’s classification ; peripheral artery diseases type 2–4 according to Fontaine ; neurological diseases (e.g. stroke, multiple sclerosis, myasthenia gravis, known cerebrovascular disease); active cancer; chronic liver diseases (hepatitis or liver failure); renal failure requiring dialysis; severe asthma and severe chronic obstructive pulmonary disease; poorly adjusted thyroid diseases; infections (tuberculosis, human immunodeficiency virus); psychiatric disorders (schizophrenia, severe depression, borderline personality disorder) and dyscalculia. Finally, patients with medications that might have influenced the conduction system of the heart or/and active or previous alcohol or drug abuse were also excluded.
First, each participant completed a questionnaire including information about demographic and clinical characteristics, such as pain, co-morbidities, and suitability for MRI-examinations. Next, structural and functional MRI was performed while the measurements of HRV at rest and under cognitive stimulation were carried out.
Information about the diagnosis, pain intensity, duration of pain, and medications were collected in the CIJD group. Pain intensity was assessed using an 11-point numerical rating scale (NRS, 0 = no pain and 10 = worst pain imaginable) . Photoplethysmographic pulse oximetry was used to determine HRV. Data recording was performed using the integrated Siemens Physiological Monitoring Unit of the body scanner (MAGNETOM Prisma, Siemens, Erlangen, Germany). In short, an infrared emitter was placed on the left index finger. The pulse curve was recorded at 50 Hz and heart rates, as well as the length of the intervals between heartbeats, were determined. All signals were transmitted wirelessly to the scanner and saved for later analysis. MR images were acquired by a research-dedicated 3 T whole-body scanner (MAGNETOM Prisma, Siemens, Erlangen, Germany) with a 64-channel head/neck coil. T1-weighted images were obtained with a magnetization prepared rapid gradient echo sequence (MP-RAGE) with the following parameters: voxel size: 0.75 × 0.75 × 0.75 mm, 224 sagittal slices, repetition time (TR): 2000 ms, echo time (TE): 2.07 ms, inversion time (TI): 952 ms, anterior to posterior phase encoding direction, flip-angle: 9°, in-plane acceleration with an acceleration factor of 2 (GRAPPA), and an acquisition time of 6:16 min.
Three functional MRI experiments were performed: an arithmetic experiment with two subtraction tasks, a number line identification experiment, and a resting-state measurement. Here, the HRV signals recorded during the arithmetic experiment and the resting-state measurement were analyzed as a marker of autonomic function . The analysis of task-based and resting-state fMRI data will be presented in another publication. As in the Trier social stress test, serial subtraction was employed to stimulate a cardiac response . The arithmetic experiment consisted of four conditions: (1) simple serial subtraction (e.g., starting with 99–7), (2) difficult serial subtraction (e.g., starting with 173–13), (3) continuously reading aloud a number presented on the screen, and (4) a rest condition, fixating a cross on the screen. Every condition lasted 20 s and was repeated with different numbers five times during the experiment. For the subtraction conditions, the participants were instructed to perform serial subtractions as quickly as possible and say out loud every result. Finally, after a short break of less than a minute, a resting state measurement took place. Short breaks between the tasks were sufficient, as the HRV should return rapidly to a normal level after temporary stimulation .
Coupling MRI tasks and pulse oximetry measurements
The Physiological Artifact Removal Tool (PART, https://www.mccauslandcenter.sc.edu/crnl/tools/part) was used to link the starting points of the MRI tasks and pulse oximetry measurements. The MRI measurements included 205 volumes each in the arithmetic task and resting state. For the analysis of the pulse oximetry measurements, the first 4 volumes were deleted. This led to a duration of 6 min and 42 s for each HRV measurement.
Analysis of heart rate variability
Data processing and calculation of the HRV via pulse oximetry measurements were performed with Kubios HRV standard (Version 3.2.0, http://www.kubios.com/), which is an established software for assessing ANS function . Artifacts such as extra, missing, or ectopic beats can cause considerable distortion of HRV results . Therefore, ‘threshold-based correction’ was employed and preset to ‘low’. The correction algorithm detected and corrected every beat-to-beat interval that was 0.35 s longer or shorter compared to the local average calculated using an internal software algorithm. Corrected data were divided into time and frequency domain parameters. Time-domain parameters included the heart rate (HR) in beats per minute (bpm), and the standard deviation of all normal-to-normal beat intervals (SDNN) in ms, which reflects the variation within the R-waves intervals. The frequency data were divided into low frequency (LF) and high frequency (HF) bands. The absolute power of the different frequency bands was specified in ms2/Hz. The frequency ranges for the absolute power measurements were as follows: LF = 0.04–0.15 Hz and HF = 0.15–0.4 Hz . LF and HF were used to characterize the response of the ANS. The HF band is thought to reflect the activity of the cardiac parasympathetic function [8, 44]. The physiological meaning of the LF is controversial. Partly representations of the sympathetic and parasympathetic nervous system, separate sections of the parasympathetic branch of the ANS, or the oscillations of the baroreflex, are discussed [8, 45, 46]. HR, SDNN, LF, and HF are some of the most used parameters for HRV analysis [8, 47].
Analysis of MR-images
The image analysis, performed with FreeSurfer (Version 6.0.0, released January 2017, http://surfer.nmr.mgh.harvard.edu) on Mac OS high sierra (Version 10.13.6), started with the extraction of the cortical (pial) surface. Due to the hybrid watershed/surface deformation procedure, non-brain tissues were removed . Additionally, segmentation of gray and white matter and intensity normalization was performed [49, 50]. The surfaces were constructed as a grid of triangles, placed in areas with the greatest intensity shift. The intensity shift was caused by the transition between different tissues . Next, the software modeled the surface between gray and white matter and the outer-pial surface, which represented the boundary of gray matter to cerebrospinal fluid. The final step was the automated topological correction to fix segmentation errors and produce a closed sphere [52, 53]. Regardless of the MRI scanner and field strength, the morphometric test-retest reliability is high . The quality of the created surfaces and subcortical structures was visually checked for each participant.
Cortical volume analysis
FreeSurfer performed an automatic reconstruction of the brain’s surface. The segmentation procedure and the parcellation of the cerebral cortex were one of the steps of the command,recon-all’, which included the extraction of the white matter and other structures, a vertex-wise reconstruction of the surfaces, smoothing, inflation, and localization of topological defects. In the end, possible tracked defects were automatically fixed. The system used the probability of the localizations of subcortical structures [50, 55]. In the process of cortical parcellation, the units were assigned according to gyral and sulcal structures. Data evaluation was based on the atlas by Desikan et al. [56, 57]. The regions of interest (ROI) for this study were the insulae as processing centers of the ANS in the left and right hemispheres. The cortical brain data were normalized using the estimated total intracranial volume (eTIV) generated by FreeSurfer [58, 59].
All statistical analyses were performed using SPSS 25 (IBM, Ehningen, Germany). Data were analyzed for normal distribution employing Shapiro-Wilk tests. SDNN, LF, and HF were logarithmically transformed to approximate normal distribution. Demographic data were first analyzed descriptively. To determine possible differences between the demographic data of the groups, an independent samples t-test (age) and Chi-Square tests (smoking, education) were used. Mixed-effects analyses of variance (ANOVA) was employed to evaluate possible differences between the HRV parameters in the different measurements and between the matched groups (C = 14 vs. CIJD = 14). Results were considered statistically significant if p < 0.05. The effect size η2 was calculated and interpreted according to Cohen (< 0.01 = small, 0.06 = medium and > 0.14 = large effect) . Next, partial correlation analyses between HRV at rest and under stimulation, together with the pain duration in years were performed with all included CIJD patients independent of age and sex . After Bonferroni correction for multiple testing, p < 0.00625 was considered significant. The obtained correlation coefficients were interpreted according to Cohen as small (r = 0.1), medium (r = 0.3), or large (r = 0.5) . Independent samples t-tests of the left and right insular cortex volume were performed between groups (C = 22, CIJD = 22). Furthermore, four linear regression models were used to test the effect of independent influences. SDNN and HR at rest were included in the first model as dependent variables while group (C = 18, CIJD = 19), age, sex, and volume of the normalized insulae in both hemispheres were defined as independent variables. HR and SDNN at rest were selected for initial analysis. Finally, a second regression with only CIJD (n = 20) was computed with SDNN and HR at rest as dependent variables. The model corresponded to the regression model described above, only that the variable age was exchanged for pain duration and the variable group was deleted. Results were considered statistically significant if p < 0.05.
Characteristics of participants
Figure 1 summarizes the recruitment process; 47 participants were initially recruited (23 controls, 24 CIJD patients). One CIJD patient was excluded due to limited compliance during the MRI acquisition and one control participant due to beta-blocker intake. Five control and four CIJD participants were excluded from HRV analysis because of incomplete pulse recordings or poor data quality. Participants were matched on an individual basis: For every CIJD patient, a control participant of the same sex and similar age (± 5 years) was recruited. For 22 CIJD patients and 22 matched controls, MRI data were available. For 19 CIJD patients and 17 controls, complete pulse recordings at rest and under cognitive stimulation of sufficient quality were available, resulting in 14 matched pairs.
Table 1 summarizes the demographic and clinical characteristics of the patient and control groups. There were no significant differences between groups regarding age (independent samples t-test; t (42) = − 0.131, p = 0.896), smoking status (Chi-Square test; χ2(2) = 0.990, p = 0.605), or formal education (Chi-Square test; χ2(3) = 2.109, p = 0.550). No participants were taking opioids, benzodiazepines, or anticonvulsants at the time of the study. Seven CIJD patients were taking glucocorticoids at the time of the study; four CIJD patients had taken glucocorticoids previously. No control participant had ever taken glucocorticoids. Antidepressants had been taken previously by one person in each group. One CIJD patient was taking a low-dose tetracyclic antidepressant at the time of the study, and none control participant. No included participant was on antihypertensive medication that could have influenced the autonomic nervous system.
Heart rate and HRV
Mean heart rate and heart rate variability at rest and under cognitive stimulation are demonstrated in Table 2. Table 3 displays the results of the mixed-effects ANOVAs for HR, SDNN, LF, and HF of the patient and control groups at rest and during cognitive stimulation. In both, controls and patients with CIJD, cognitive stimulation resulted in a significant increase in HR (C: + 10.68%, CIJD: + 10.73%) with a large effect size and SDNN (C: + 6.11%, CIJD: + 8.08%) with a medium effect size. Importantly, there were no significant effects between groups (Fig. 2, The illustration is based on Allen et al. (2019), using code for R’s ggplot2 package available at https://github.com/RainCloudPlots/RainCloudPlots) . Furthermore, the LF increase in CIJD patients (C: + 2.23%, CIJD: + 12.26%) was significant with a medium effect size and no differences in the HF parameters between groups.
CIJD patients (n = 20) demonstrated a significant negative correlation between heart rate at rest and pain duration independent of age and sex (r = − 0.743 p = < 0.001, Table 4, Fig. 3A). During cognitive stimulation (n = 22), there was no correlation between the heart rate and pain duration (r = − 0.383, p = 0.096, Table 4, Fig. 3B). There were no correlations between SDNN, LF, or HF and pain duration, neither at rest nor during cognitive stimulation independent of age and sex (Table 4). However, age and pain duration were not significantly correlated (r = 0.403, p = 0.057).
Employing independent t-tests of insula volume with eTIV normalization (C = 22, CIJD = 22), we found no significant differences between the matched groups in the left (t (42) = 0.214, p = 0.832) and the right insula (t (42) = 0.165, p = 0.870). The general linear regression model including both groups (C = 18, CIJD = 19) with HR as dependent variable showed no significant differences (adjusted R2 = 0.022, p = 0.349). The SDNN model including group, age, sex, and insula volume in both hemispheres showed significant results and explained 23% of the variability (R2 = 0.227, p = 0.022). However, the only relevant variable was sex. The second subgroup regression model (n = 20) was used to predict HR with pain duration, sex, and volume of the insular cortices in both hemispheres. The subgroup regression explained 66% of the HR and 39% of the SDNN variability (HR: R2 = 0.660, p < 0.001; SDNN: R2 = 0.389, p = 0.021). Tables 5 and 6 shows the significant results of both regression analyses.
Our results provide evidence for an inverse correlation between HR and pain duration in patients with chronic inflammatory joint disease at rest. Moreover, both healthy volunteers and patients with CIJD showed a similarly adequate cardiovascular response during cognitive stimulation (serial subtraction). Finally, our data suggest that structural changes in the insular cortices are related to changes in heart rate in CIJD patients.
Heart rate and heart rate variability
The present data provide evidence for a negative association between HR and pain duration with a large effect size (r = − 0.743) at rest, while controlling for the influence of age and sex. To the best of our knowledge, this is the first report of an association between resting HR and pain duration in patients with CIJD. The dependence of age and HR decrease was already reported in previous studies , but could not be verified here. The present association of HR and pain duration may be explained by an autonomic imbalance, for example, an increase in the parasympathetic or a decrease in the sympathetic outflow, or a generalized downregulation of both ANS branches . In an animal model of chronic stress, a decrease of HR was found after stressing mice on a shaking platform for 7 days. This study demonstrated an initial increase of HR in response to the stressor, followed by a transition to lower heart rates during the experiment . One possible explanation could be an adaptation mechanism under longer-lasting stressful or nociceptive stimulations. For instance, an increased activity of the parasympathetic branch of the ANS could act as a protective factor by lowering the heart rate, e.g., to provide more time for blood flow within coronary arteries in diastole and thus achieve better oxygen supply. This could result in maintaining cardiac adaptability for as long as possible . This notion would fit with our SDNN results in CIJD patients, which are still in the normal range and comparable to the control group. However, no differentiation of the precise mechanism was possible within the setting of our study. Therefore, more research is necessary to elucidate this issue further.
In humans, studies on HR in chronic pain caused by CIJD provided mixed results. For example, Rensburg and colleagues found a significantly higher basal HR in women with RA , whereas Louthrenoo et al. used a non-invasive cardiovascular reflex test (deep breathing) and found no significant difference in maximum and minimum HR in patients with RA compared with healthy controls . The difference between the maximum and minimum HR, however, was significantly smaller in patients with RA compared to controls. Both studies included RA patients with a relatively short mean disease duration of 4.26 ± 1.2 years  and 5.1 ± 3.6 years . In contrast, disease duration in our study was between 1 and 44 years (median: 10 years). It is important to note that there was no significant association between HR under stimulation and pain duration (Fig. 3B), suggesting that autonomic control remains intact even in later stages of CIJD.
Similarly, studies on heart rate variability in patients with chronic pain are inconclusive and the selected HRV parameters differ . For example, in 2006, Anichkov and colleagues detected significantly lower HRV in 23 female RA patients compared to matched healthy controls . In contrast, Vlcek and colleagues did not find changes in HRV at rest and under light physical stress in a comparable group (22 female RA patients vs. 15 matched healthy controls) .
In our study, we not only measured HR and HRV at rest but also during cognitive stimulation. As expected, HR increased in healthy controls during serial subtraction . Interestingly, our CIJD patients demonstrated a similar increase in HR. So far, the most frequently described response to acute stimulation is an HR and LF band increase and a HRV and HF band decrease. This is thought to reflect either a reduction in the parasympathetic or a rise in the sympathetic outflow [63, 70]. In our study, both groups reacted with an increase in HR and SDNN. Since previous studies suggested reduced cardiac function in patients with CIJD, we expected to find a decreased SDNN in our patient group. Nevertheless, CJID patients in our study seemed to have good cardiac adaptability and reacted appropriately to the stressor. However, only the CIJD reacted with a significant increase in the LF band to cognitive stimulation. The increased LF value in CIJD patients might represent a modified baroreflex . Autonomic imbalance is associated with increased cardiovascular morbidity and mortality . However, to assess the baroreflex more thoroughly was beyond the scope of this study.
Volumes of the left and right insular cortices
Processing of nociceptive input is a highly complex function of the brain, involving multiple brain regions as well as excitatory and inhibitory mechanisms. Structural MRI studies in chronic pain patients have been inconclusive so far with some showing increased and others, decreased cerebral gray matter density [72, 73]. For this study, we performed an automated analysis of the volumes of the left and right insular cortex using FreeSurfer. In a multiple linear regression model including only CIJD patients with HR at rest as a dependent variable, we found a significant positive relationship with the volume of the left insular cortex. Moreover, we found a significant negative relationship between HR at rest and the volume of the right insula (Table 5). There is still much debate about the lateralization of insular functions, including autonomic control [14, 74]. A widely acknowledged model of insular lateralization proposes that the right insula predominantly activates sympathetic outflow while the left insula mainly initiates parasympathetic outflow [75, 76]. This model is supported by studies performing unilateral inactivation of the cerebral hemispheres by intracarotid amobarbital injection  and using electrical stimulation of the insulae during epilepsy surgery . Of major clinical importance is the fact that right insular lesions are particularly associated with the development of cardiac arrhythmias [23, 78, 79]. However, this model does not explain all available data and may thus be oversimplified. A recent intraoperative stimulation study did not find differences in HR response between the left and right insula, but between the anterior and posterior parts of both insulae . Left or right posterior insular cortical stimulation predominantly induced tachycardia . Given these controversial findings, the pathomechanistic implications of our results are difficult to assess. Further studies on larger patient populations are needed to replicate these results and to determine if changes in insular volumes might be associated with the increased cardiovascular morbidity and mortality observed in CIJD patients.
Limitations and future research
One of the major strengths of this study is the matching of patients and controls on an individual basis. Further, this study examined diseases that are based on similar pathophysiology. Nevertheless, further research and longitudinal studies are needed to investigate the influence of chronic pain in patients with rheumatic diseases. There is some evidence that changes in brain areas caused by chronic pain were reversible after pain relief . It is possible that some of the morphological brain changes were diminished due to the plasticity of the brain under consistent pain therapy and pain relief or even occurred in other regions that were not examined in this study. It thus would be interesting to follow newly diagnosed patients with rheumatic diseases before and after drug treatment in a long-term study to compare heart and brain modifications. Using only cognitive instead of physical stimulation can be critical. This study chose mental stimulation due to its relevance to everyday life and to maintain the evaluability of pulse-oxymetrical measured data. Further, a focus on physical stimulations would have required movements, e.g. handgrip exercise, which in turn would be a risk for creating MRI and pulse data artifacts [82,83,84]. Therefore, interpretation of the frequency domain parameter measured by photoplethysmographic pulse oximetry, especially HF, would only have been possible to a limited extent. In a future study, it would be interesting to examine whether the CIJD group would reach a limit of their heart adaptability earlier than the control group under physical stimulation using ECG. Finally, our group size corresponds to other studies that dealt with similar MRI questions or methods [27, 28, 30, 33]. Due to the lack of data including both research fields in chronic pain patients (MRI and HRV), no a priori sample size calculation could be performed. Possibly due to the large number of partial correlation analyses together with conservative correction for multiple comparisons, we only found a significant association of pain duration with HR at rest. We assume that other variables, for example HF, are less sensitive than HR to changes in CIJD patients, thus requiring larger sample sizes. Therefore, future studies with larger sample sizes are needed for further evaluation and detailed characterization of the underlying mechanisms of ANS changes in CIJD.
In summary, our results suggest that HR at rest is negatively associated with the pain duration of CIJD patients, independent of the participants’ age and sex. Further, our data established the possibility of insular cortex involvement in the process of ANS changes.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Autonomic nervous system
Beats per minute
Chronic inflammatory joint diseases
Estimated intracranial volume
Heart rate variability
Magnetization prepared rapid gradient echo
Magnetic resonance imaging
Numerical rating scale
Physiological artifact removal tool
Standard deviation of all normal-to-normal beat intervals
Kazis LE, Meenan RF, Anderson JJ. Pain in the rheumatic diseases. Arthritis Rheum. 1983;26(8):1017–22.
Cavalli G, Favalli EG. Cardiovascular disease in patients with rheumatoid arthritis: impact of classic and disease-specific risk factors. Ann Transl Med. 2018;6(Suppl 1):S82.
Solomon DH, Karlson EW, Rimm EB, Cannuscio CC, Mandl LA, Manson JE, et al. Cardiovascular morbidity and mortality in women diagnosed with rheumatoid arthritis. Circulation. 2003;107(9):1303–7.
England BR, Thiele GM, Anderson DR, Mikuls TR. Increased cardiovascular risk in rheumatoid arthritis: mechanisms and implications. BMJ. 2018;361:k1036.
Wolfe F, Mitchell DM, Sibley JT, Fries JF, Bloch DA, Williams CA, et al. The mortality of rheumatoid arthritis. Arthritis Rheum. 1994;37(4):481–94.
Masoud S, Lim PB, Kitas GD, Panoulas V. Sudden cardiac death in patients with rheumatoid arthritis. World J Cardiol. 2017;9(7):562–73.
Adlan AM, Lip GY, Paton JF, Kitas GD, Fisher JP. Autonomic function and rheumatoid arthritis: a systematic review. Semin Arthritis Rheum. 2014;44(3):283–304.
Shaffer F, Ginsberg JP. An overview of heart rate variability metrics and norms. Front Public Health. 2017;5:258.
Hilz MJ, Dütsch M. Quantitative studies of autonomic function. Muscle Nerve. 2006;33(1):6–20.
Janse van Rensburg DC, Ker JA, Grant CC, Fletcher L. Autonomic impairment in rheumatoid arthritis. Int J Rheum Dis. 2012;15(4):419–26.
Evrengul H, Dursunoglu D, Cobankara V, Polat B, Seleci D, Kabukcu S, et al. Heart rate variability in patients with rheumatoid arthritis. Rheumatol Int. 2004;24(4):198–202.
Dunn KM, Croft PR. The importance of symptom duration in determining prognosis. Pain. 2006;121(1–2):126–32.
Bantel C, Trapp S. The role of the autonomic nervous system in acute surgical pain processing - what do we know? Anaesthesia. 2011;66(7):541–4.
Ruiz Vargas E, Soros P, Shoemaker JK, Hachinski V. Human cerebral circuitry related to cardiac control: a neuroimaging meta-analysis. Ann Neurol. 2016;79(5):709–16.
Evrard HC. The organization of the primate insular cortex. Front Neuroanat. 2019;13:43.
Sörös P, Lalone E, Smith R, Stevens T, Theurer J, Menon RS, et al. Functional MRI of oropharyngeal air-pulse stimulation. Neuroscience. 2008;153(4):1300–8.
Sörös P, Marmurek J, Tam F, Baker N, Staines WR, Graham SJ. Functional MRI of working memory and selective attention in vibrotactile frequency discrimination. BMC Neurosci. 2007;8:48.
Kraynak TE, Marsland AL, Wager TD, Gianaros PJ. Functional neuroanatomy of peripheral inflammatory physiology: a meta-analysis of human neuroimaging studies. Neurosci Biobehav Rev. 2018;94:76–92.
Frot M, Faillenot I, Mauguière F. Processing of nociceptive input from posterior to anterior insula in humans. Hum Brain Mapp. 2014;35(11):5486–99.
Ackermann H, Riecker A. The contribution(s) of the insula to speech production: a review of the clinical and functional imaging literature. Brain Struct Funct. 2010;214(5–6):419–33.
Sörös P, Inamoto Y, Martin RE. Functional brain imaging of swallowing: an activation likelihood estimation meta-analysis. Hum Brain Mapp. 2009;30(8):2426–39.
Oppenheimer SM, Gelb A, Girvin JP, Hachinski VC. Cardiovascular effects of human insular cortex stimulation. Neurology. 1992;42(9):1727–32.
Christensen H, Boysen G, Christensen AF, Johannesen HH. Insular lesions, ECG abnormalities, and outcome in acute stroke. J Neurol Neurosurg Psychiatry. 2005;76(2):269–71.
Nagai M, Dote K, Kato M, Sasaki S. The right insular cortex: a critical region for modulating the parasympathetic activity? J Neurosurg Anesthesiol. 2016;28(2):176.
Shu L, Wang Y. Highlights in basic autonomic neuroscience: insular cortex injury leads to cardiovascular dysfunction. Auton Neurosci. 2014;185:1–4.
Sörös P, Hachinski V. Cardiovascular and neurological causes of sudden death after ischaemic stroke. Lancet Neurol. 2012;11(2):179–88.
Wartolowska K, Hough MG, Jenkinson M, Andersson J, Wordsworth BP, Tracey I. Structural changes of the brain in rheumatoid arthritis. Arthritis Rheum. 2012;64(2):371–9.
Apkarian AV, Sosa Y, Sonty S, Levy RM, Harden RN, Parrish TB, et al. Chronic back pain is associated with decreased prefrontal and thalamic gray matter density. J Neurosci. 2004;24(46):10410–5.
Ruscheweyh R, Deppe M, Lohmann H, Stehling C, Floel A, Ringelstein EB, et al. Pain is associated with regional grey matter reduction in the general population. Pain. 2011;152(4):904–11.
Geha PY, Baliki MN, Harden RN, Bauer WR, Parrish TB, Apkarian AV. The brain in chronic CRPS pain: abnormal gray-white matter interactions in emotional and autonomic regions. Neuron. 2008;60(4):570–81.
Schmidt-Wilcke T, Gänssbauer S, Neuner T, Bogdahn U, May A. Subtle grey matter changes between migraine patients and healthy controls. Cephalalgia. 2008;28(1):1–4.
Schmidt-Wilcke T, Leinisch E, Straube A, Kampfe N, Draganski B, Diener HC, et al. Gray matter decrease in patients with chronic tension type headache. Neurology. 2005;65(9):1483–6.
Kuchinad A, Schweinhardt P, Seminowicz DA, Wood PB, Chizh BA, Bushnell MC. Accelerated brain gray matter loss in fibromyalgia patients: premature aging of the brain? J Neurosci. 2007;27(15):4004–7.
Spindler M, Koch K, Borisov E, Ozyurt J, Soros P, Thiel C, et al. The influence of chronic pain and cognitive function on spatial-numerical processing. Front Behav Neurosci. 2018;12:165.
Sörös P, Bantel C. Chronic noncancer pain is not associated with accelerated brain aging as assessed by structural magnetic resonance imaging in patients treated in specialized outpatient clinics. Pain. 2020;161(3):641–50.
Ponikowski P, Voors AA, Anker SD, Bueno H, Cleland JG, Coats AJ, et al. 2016 ESC guidelines for the diagnosis and treatment of acute and chronic heart failure: the task force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC) developed with the special contribution of the heart failure association (HFA) of the ESC. Eur Heart J. 2016;37(27):2129–200.
Aboyans V, Ricco J-B, Bartelink M-LE, Björck M, Brodmann M, Cohnert T, Collet J-P, Czerny M, De Carlo M, Debus S. 2017 ESC Guidelines on the Diagnosis and Treatment of Peripheral Arterial Diseases, in collaboration with the European Society for Vascular Surgery (ESVS) Document covering atherosclerotic disease of extracranial carotid and vertebral, mesenteric, renal, upper and lower extremity arteries Endorsed by: the European Stroke Organization (ESO) The Task Force for the Diagnosis and Treatment of Peripheral Arterial Diseases of the European Society of Cardiology (ESC) and of the European Society for Vascular Surgery (ESVS). Eur Heart J. 2018;39(9):763–816.
Breivik H, Borchgrevink P, Allen S, Rosseland L, Romundstad L, Hals EB, et al. Assessment of pain. Br J Anaesth. 2008;101(1):17–24.
Thayer JF, Ahs F, Fredrikson M, Sollers JJ 3rd, Wager TD. A meta-analysis of heart rate variability and neuroimaging studies: implications for heart rate variability as a marker of stress and health. Neurosci Biobehav Rev. 2012;36(2):747–56.
Kirschbaum C, Pirke K-M, Hellhammer DH. The ‘Trier social stress test’–a tool for investigating psychobiological stress responses in a laboratory setting. Neuropsychobiology. 1993;28(1–2):76–81.
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. Eur Heart J. 1996;17(3):354–81.
Tarvainen MP, Niskanen JP, Lipponen JA, Ranta-Aho PO, Karjalainen PA. Kubios HRV-heart rate variability analysis software. Comput Methods Prog Biomed. 2014;113(1):210–20.
Berntson GG, Stowell JR. ECG artifacts and heart period variability: don't miss a beat! Psychophysiology. 1998;35(1):127–32.
Montano N, Ruscone TG, Porta A, Lombardi F, Pagani M, Malliani A. Power spectrum analysis of heart rate variability to assess the changes in sympathovagal balance during graded orthostatic tilt. Circulation. 1994;90(4):1826–31.
Goldstein DS, Bentho O, Park MY, Sharabi Y. Low-frequency power of heart rate variability is not a measure of cardiac sympathetic tone but may be a measure of modulation of cardiac autonomic outflows by baroreflexes. Exp Physiol. 2011;96(12):1255–61.
Reyes del Paso GA, Langewitz W, Mulder LJ, van Roon A, Duschek S. The utility of low frequency heart rate variability as an index of sympathetic cardiac tone: a review with emphasis on a reanalysis of previous studies. Psychophysiology. 2013;50(5):477–87.
Niskanen JP, Tarvainen MP, Ranta-Aho PO, Karjalainen PA. Software for advanced HRV analysis. Comput Methods Prog Biomed. 2004;76(1):73–81.
Segonne F, Dale AM, Busa E, Glessner M, Salat D, Hahn HK, et al. A hybrid approach to the skull stripping problem in MRI. Neuroimage. 2004;22(3):1060–75.
Sled JG, Zijdenbos AP, Evans AC. A nonparametric method for automatic correction of intensity nonuniformity in MRI data. IEEE Trans Med Imaging. 1998;17(1):87–97.
Fischl B, Salat DH, Busa E, Albert M, Dieterich M, Haselgrove C, et al. Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain. Neuron. 2002;33(3):341–55.
Dale AM, Fischl B, Sereno MI. Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage. 1999;9(2):179–94.
Fischl B, Liu A, Dale AM. Automated manifold surgery: constructing geometrically accurate and topologically correct models of the human cerebral cortex. IEEE Trans Med Imaging. 2001;20(1):70–80.
Segonne F, Pacheco J, Fischl B. Geometrically accurate topology-correction of cortical surfaces using nonseparating loops. IEEE Trans Med Imaging. 2007;26(4):518–29.
Han X, Jovicich J, Salat D, van der Kouwe A, Quinn B, Czanner S, et al. Reliability of MRI-derived measurements of human cerebral cortical thickness: the effects of field strength, scanner upgrade and manufacturer. Neuroimage. 2006;32(1):180–94.
Fischl B, Salat DH, van der Kouwe AJ, Makris N, Segonne F, Quinn BT, et al. Sequence-independent segmentation of magnetic resonance images. Neuroimage. 2004;23(Suppl 1):S69–84.
Fischl B, Van Der Kouwe A, Destrieux C, Halgren E, Ségonne F, Salat DH, et al. Automatically parcellating the human cerebral cortex. Cereb Cortex. 2004;14(1):11–22.
Desikan RS, Ségonne F, Fischl B, Quinn BT, Dickerson BC, Blacker D, et al. An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage. 2006;31(3):968–80.
Buckner RL, Head D, Parker J, Fotenos AF, Marcus D, Morris JC, et al. A unified approach for morphometric and functional data analysis in young, old, and demented adults using automated atlas-based head size normalization: reliability and validation against manual measurement of total intracranial volume. Neuroimage. 2004;23(2):724–38.
Mathalon DH, Sullivan EV, Rawles JM, Pfefferbaum A. Correction for head size in brain-imaging measurements. Psychiatry Res. 1993;50(2):121–39.
Cohen J. Statistical power analysis for the behavioural sciences. Hillsdale: Erlbaum; 1988.
Antelmi I, de Paula RS, Shinzato AR, Peres CA, Mansur AJ, Grupi CJ. Influence of age, gender, body mass index, and functional capacity on heart rate variability in a cohort of subjects without heart disease. Am J Cardiol. 2004;93(3):381–5.
Allen M, Poggiali D, Whitaker K, Marshall TR, Kievit RA. Raincloud plots: a multi-platform tool for robust data visualization. Wellcome Open Res. 2019;4:63.
Tracy LM, Ioannou L, Baker KS, Gibson SJ, Georgiou-Karistianis N, Giummarra MJJP. Meta-analytic evidence for decreased heart rate variability in chronic pain implicating parasympathetic nervous system dysregulation. Pain. 2016;157(1):7–29.
Bernatova I, Key MP, Lucot JB, Morris M. Circadian differences in stress-induced pressor reactivity in mice. Hypertension. 2002;40(5):768–73.
Selye H. Stress without distress. In: Serban G, editor. Psychopathology of human adaptation. Boston: Springer; 1976. p. 137–46.
Louthrenoo W, Ruttanaumpawan P, Aramrattana A, Sukitawut W. Cardiovascular autonomic nervous system dysfunction in patients with rheumatoid arthritis and systemic lupus erythematosus. QJM. 1999;92(2):97–102.
Anichkov DA, Shostak NA, Ivanov DS. Heart rate variability is related to disease activity and smoking in rheumatoid arthritis patients. Int J Clin Pract. 2007;61(5):777–83.
Vlcek M, Rovensky J, Eisenhofer G, Radikova Z, Penesova A, Kerlik J, et al. Autonomic nervous system function in rheumatoid arthritis. Cell Mol Neurobiol. 2012;32(5):897–901.
Allen TM, Struemph KL, Toledo-Tamula MA, Wolters PL, Baldwin A, Widemann B, et al. The relationship between heart rate variability, psychological flexibility, and pain in Neurofibromatosis type 1. Pain Pract. 2018;18(8):969–78.
Kim G, Cheon EJ, Bai DS, Lee YH, Koo BH. Stress and heart rate variability: a Meta-analysis and review of the literature. Psychiatry Investig. 2018;15(3):235–45.
Thayer JF, Yamamoto SS, Brosschot JF. The relationship of autonomic imbalance, heart rate variability and cardiovascular disease risk factors. Int J Cardiol. 2010;141(2):122–31.
Smallwood RF, Laird AR, Ramage AE, Parkinson AL, Lewis J, Clauw DJ, et al. Structural brain anomalies and chronic pain: a quantitative meta-analysis of gray matter volume. J Pain. 2013;14(7):663–75.
Shi H, Yuan C, Dai Z, Ma H, Sheng L. Gray matter abnormalities associated with fibromyalgia: a meta-analysis of voxel-based morphometric studies. Semin Arthritis Rheum. 2016;46:330–7.
Oppenheimer S. Cerebrogenic cardiac arrhythmias: cortical lateralization and clinical significance. Clin Auton Res. 2006;16(1):6.
Craig A. Forebrain emotional asymmetry: a neuroanatomical basis? Trends Cogn Sci. 2005;9(12):566–71.
Nagai M, Hoshide S, Kario K. The insular cortex and cardiovascular system: a new insight into the brain-heart axis. J Am Soc Hypertens. 2010;4(4):174–82.
Hilz MJ, Dütsch M, Perrine K, Nelson PK, Rauhut U, Devinsky O. Hemispheric influence on autonomic modulation and baroreflex sensitivity. Ann Neurol. 2001;49(5):575–84.
Colivicchi F, Bassi A, Santini M, Caltagirone C. Cardiac autonomic derangement and arrhythmias in right-sided stroke with insular involvement. Stroke. 2004;35(9):2094–8.
Abboud H, Berroir S, Labreuche J, Orjuela K, Amarenco P. Insular involvement in brain infarction increases risk for cardiac arrhythmia and death. Ann Neurol. 2006;59(4):691–9.
Chouchou F, Mauguière F, Vallayer O, Catenoix H, Isnard J, Montavont A, et al. How the insula speaks to the heart: cardiac responses to insular stimulation in humans. Hum Brain Mapp. 2019;40(9):2611–22.
Gwilym SE, Filippini N, Douaud G, Carr AJ, Tracey I. Thalamic atrophy associated with painful osteoarthritis of the hip is reversible after arthroplasty: a longitudinal voxel-based morphometric study. Arthritis Rheum. 2010;62(10):2930–40.
Lalande S, Sawicki CP, Baker JR, Shoemaker JK. Effect of age on the hemodynamic and sympathetic responses at the onset of isometric handgrip exercise. J Appl Physiol (1985). 2014;116(2):222–7.
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.
Schäfer A, Vagedes J. How accurate is pulse rate variability as an estimate of heart rate variability? A review on studies comparing photoplethysmographic technology with an electrocardiogram. Int J Cardiol. 2013;166(1):15–29.
The authors want to acknowledge the rheumatological clinic of Dr. Markus Voglau for the great cooperation in patient recruitment. We thank Katharina Grote and Gülsen Yanç for assisting with MRI data acquisition. We also thank Vivian Baker for editing the final manuscript.
Open Access funding enabled and organized by Projekt DEAL. This work was supported by the Neuroimaging Unit, University of Oldenburg, funded by grants from the German Research Foundation (DFG; 3 T MRI INST 184/152–1 FUGG and MEG INST 184/148–1 FUGG).
Ethics approval and consent to participate
The study was approved by the Medical Research Ethics Board of Carl von Ossietzky University of Oldenburg, Germany (#2017–059). All methods were carried out in accordance with relevant guidelines and regulations. All participants gave written informed consent prior to participating in this study.
Consent for publication
Not applicable (the manuscript does not include details, images, or videos relating to an individual person).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Danilin, L.K., Spindler, M., Sörös, P. et al. Heart rate and heart rate variability in patients with chronic inflammatory joint disease: the role of pain duration and the insular cortex. BMC Musculoskelet Disord 23, 75 (2022). https://doi.org/10.1186/s12891-022-05009-1
- Chronic inflammatory joint diseases
- Chronic pain
- Pain duration
- Insular cortex
- Heart rate
- Heart rate variability