Decreased heart rate variability (HRV) is associated with adverse outcomes in cardiovascular diseases and has been observed in patients with systemic lupus erythematosus (SLE). We examined the relationship of HRV with SLE disease activity and selected cytokine pathways.
Fifty-three patients from the Oklahoma Lupus Cohort were evaluated at two visits each. Clinical assessments included the Systemic Lupus Erythematosus Disease Activity Index (SLEDAI), British Isles Lupus Assessment Group (BILAG) index, physician global assessment (PGA), and Safety of Estrogens in Lupus Erythematosus National Assessment-SLEDAI Flare Index. HRV was assessed with a 5-minute electrocardiogram, and the following HRV parameters were calculated: square root of the mean of the squares of differences between adjacent NN intervals (RMSSD), percentage of pairs of adjacent NN intervals differing by more than 50 milliseconds (pNN50), high-frequency power (HF power), and low frequency to high frequency (LF/HF) ratio, which reflects sympathetic/vagal balance. Plasma cytokine levels were measured with a multiplex, bead-based immunoassay. Serum B lymphocyte stimulator (BLyS) and a proliferation-inducing ligand (APRIL) were measured with an enzyme-linked immunosorbent assay. Linear regression analysis was applied.
Baseline HRV (pNN50, HF power, LF/HF ratio) was inversely related to disease activity (BILAG, PGA) and flare. Changes in RMSSD between visits were inversely related to changes in SLEDAI (p = 0.007). Age, caffeine, tobacco and medication use had no impact on HRV. Plasma soluble tumor necrosis factor receptor II (sTNFRII) and monokine induced by interferon gamma (MIG) were inversely related with all baseline measures of HRV (p = 0.039 to <0.001). Plasma stem cell factor (SCF), interleukin (IL)-1 receptor antagonist (IL-1RA), and IL-15 showed similar inverse relationships with baseline HRV, and weaker trends were observed for interferon (IFN)-α, interferon gamma-induced protein (IP)-10, and serum BLyS. Changes in the LF/HF ratio between visits were also associated with changes in sTNFRII (p = 0.021), MIG (p = 0.003), IFN-α (p = 0.012), SCF (p = 0.001), IL-1RA (p = 0.023), and IL-15 (p = 0.010). On the basis of multivariate linear regression, MIG was an independent predictor of baseline HRV after adjusting for plasma IL-1RA, SCF, IFN-α, IP-10, and serum BLyS. In a similar model, the sTNFRII impact remained significant after adjusting for the same variables.
Impaired HRV, particularly the LF/HF ratio, is associated with lupus disease activity and several cytokines related to IFN type II and TNF pathways. The strongest association was with MIG and sTNFRII, expanding previous immune connections of vagal signaling.
Neural influence on immune function can been profound, raising the question of neuroimmune modulation in chronic inflammatory conditions such as systemic lupus erythematosus (SLE) . Autonomic dysfunction is well documented in SLE, with increased sympathetic and decreased parasympathetic activity reported by different studies [2–9]. Such autonomic aberrations can be the result of inflammatory cytokines acting on the central nervous system and may reciprocally modulate inflammatory responses in the periphery, exerting a pathogenic role in chronic inflammation .
Both sympathetic and parasympathetic pathways are relevant in acute and chronic inflammation, based on in vitro studies and animal models. Upon initiation of an acute inflammatory process, the body adopts a nonspecific “inflammatory configuration” with increased systemic sympathetic activity and hypothalamic-pituitary-adrenal axis activation . At the same time, sympathetic neurotransmitters act directly on immune cells, promoting primarily anti-inflammatory mechanisms at the tissue level. As inflammation becomes chronic, tissue regulation is uncoupled from central influence by repulsion of sympathetic nerve fibers that promotes an inflammatory phenotype . Signaling through the vagus nerve has strong anti-inflammatory effects in tissue-specific and systemic inflammation . Efferent potentials from the brainstem inhibit cytokine production in the spleen through a complex pathway that also includes the splanchnic sympathetic nerves . Recent understanding of this pathway suggests a non-neural link between the vagus nerve and the spleen, possibly through vagus-induced mobilization of peripheral lymphocyte pools that subsequently accumulate in the spleen .
Autonomic activity can be measured through a range of techniques that include clinical cardiovascular reflex testing and heart rate variability (HRV) . HRV is a noninvasive tool that measures cyclical fluctuations in resting heart rate to indirectly assess cardiac parasympathetic and sympathetic influences. Although vagal tone is considered organ-specific, HRV is often used as a surrogate for systemic vagal output. Decreased HRV is linked to increased cardiovascular morbidity and mortality, as well as an increased prevalence of heart failure and arrhythmias after myocardial infarction [14, 15]. Measures of HRV have also been inversely correlated with serum inflammatory biomarkers (C-reactive protein, interleukin [IL]-6) in healthy individuals as well as in those with cardiovascular diseases .
We hypothesized that aberrant autonomic function reflected by HRV is related to disease activity in SLE, and we examined this hypothesis cross-sectionally and longitudinally in a cohort of well-characterized patients with active disease who were receiving standard-of-care treatments. We further explored associations of HRV with a large panel of cytokines, including several pertinent to SLE.
We conducted a prospective cohort study at the Department of Arthritis and Clinical Immunology of the Oklahoma Medical Research Foundation (OMRF) as a project of the Oklahoma Lupus Cohort, a longitudinal study that includes 610 patients who meet the 1997 SLE classification criteria . The protocol was approved by the OMRF Institutional Review Board (IRB). All patients underwent informed consent procedures consistent with the governance of the OMRF IRB and the Declaration of Helsinki. Fifty-three consecutive patients with SLE from the Oklahoma Lupus Cohort with active disease at baseline were examined at two time points while standard of care for lupus was rendered. Active disease at baseline was defined as Systemic Lupus Erythematosus Disease Activity Index (SLEDAI) of at least 4 points or at least one B score on the British Isles Lupus Assessment Group (BILAG) index . Patients were evaluated using the hybrid SLEDAI, which is identical to the Safety of Estrogens in Lupus National Assessment (SELENA)-SLEDAI [18, 19] except for the proteinuria definition from the SLEDAI 2000 (SLEDAI-2K) [20, 21]. Presence of flare at each visit was determined by the SELENA-SLEDAI Flare Index (SFI) with the physician global assessment (PGA) [18, 19]. Laboratory tests at each visit included complete blood counts, a comprehensive metabolic panel, urinalysis and C3 and C4 complement levels. In addition to these routine laboratory assessments, serum B lymphocyte stimulator (BLyS) and a proliferation-inducing ligand (APRIL), as well as a panel of plasma cytokines were measured at each visit by enzyme-linked immunosorbent assay and by xMAP® multiplex bead-based immunoassay (eBioscience/Affymetrix, San Diego, California, USA) respectively. The following plasma cytokines were included in the panel: interferon (IFN)-γ, IL-12p70, IL-2, IL-1α, IL-1β, interleukin (IL)-1 receptor antagonist (IL-1RA), IL-15, IL-17A, IL-21, tumor necrosis factor (TNF)-α, soluble tumor necrosis factor receptors I and II (sTNFRI/II), soluble CD40L (sCD40L), BLyS, APRIL, transforming growth factor (TGF)-β, IL-10, IFN-α, interferon gamma-induced protein 10 (IP-10), macrophage inflammatory protein 1 alpha (MIP-1α), monokine induced by interferon gamma (MIG), IL-8, plasminogen activator inhibitor 1 (PAI-1), stem cell factor (SCF), and resistin (see Table S1 in Additional file 1).
Medical records were reviewed and a questionnaire about variables that might affect HRV was completed at each clinic visit. Potentially confounding variables that were evaluated included any heart rhythm other than sinus, cardiomyopathy (left ventricular ejection fraction <40 %), recent (<1 year) myocardial infarction or unstable angina, heart failure (New York Heart Association class III or IV), uncontrolled hypertension (systolic blood pressure >160 mmHg on at least three antihypertensive medications), mitral valve prolapse (other than mild), recurrent vasovagal syncopal episodes, unilateral or bilateral vagotomy, concomitant use of medications potentially affecting autonomic function (e.g., anticholinergic medications and β-blockers), significant neurologic disorder (e.g., Parkinson’s disease, multiple sclerosis, Guillain-Barré syndrome), major untreated depression or psychosis, obstructive sleep apnea, asthma or chronic obstructive pulmonary disease not controlled by medications, any other disease causing clinically significant dyspnea at the time of assessment, and uncontrolled insulin- or non-insulin-dependent diabetes. Patients were also questioned about recent exercise and tobacco use (preceding 12 hours), any caffeine consumed in the 6 hours prior to testing, and alcohol consumed within 12 hours of testing.
A 5 minute electrocardiogram (ECG) tracing was recorded at each visit using a modification of the commercially available AliveCor iPhone ECG device (AliveCor, San Francisco, CA, USA). ECG electrode patches were attached to the right and left forearms and connected to the device with a custom-made attachment to enable stable ECG recording (lead I) over 5 minutes. Recordings were transmitted online to a Health Insurance Portability and Accountability Act (HIPAA)-compliant server with validated encryption technology. ECG results were processed blindly, and the following HRV parameters were obtained by time domain analysis: (1) square root of the mean of the squares of differences between adjacent NN intervals (RMSSD) and (2) percentage of pairs of adjacent NN intervals differing by more than 50 milliseconds (pNN50). The following HRV parameters were obtained by frequency domain analysis: (1) low-frequency power (LF power), (2) high-frequency power (HF power), and (3) the LF/HF ratio . Both short-duration (5 minutes) and long-duration (24 hours) ECG recordings are acceptable options for measuring HRV . Nonetheless, frequency domain analysis is the preferred method when short-duration recordings are examined . Normative values of these HRV measures have yet to be established, and HRV testing in isolation cannot differentiate increased sympathetic from decreased parasympathetic signaling. However, HF power overall reflects activity of the efferent vagus nerve to the heart, and LF power is considered a measure of sympathetic activity, whereas the ratio of low to high frequency (LF/HF) reflects sympathetic to parasympathetic balance . Reduced RMSSD and pNN50 are also indicative of parasympathetic insufficiency.
Categorical variables were expressed as percentages and continuous variables as median and interquartile range (IQR). Paired t tests or Wilcoxon signed-rank tests were used to evaluate changes in disease activity and HRV between visits. A chi-square test with the Yates correction was used to compare immune-suppressive medications and parameters affecting HRV between visits. We used linear regression to evaluate cross-sectional associations of HRV and SLE disease activity, cytokines, medications, and tobacco, caffeine, and alcohol exposure. The natural logarithm of each cytokine was employed to achieve normality. Factors significant in univariate analysis (p < 0.1) were entered into multivariate models. To avoid overfitting, no more than one parameter for each ten subjects was examined in the multivariate regression models, which allows for a nonbiased calculation of the regression coefficients . Changes in HRV between baseline and follow-up were also examined as predictors of changes in disease activity by linear regression. The same approach was used to evaluate whether changes in HRV between baseline and follow-up could be predicted by changes in cytokines. In these models, only cytokines found to be significantly associated with HRV in the baseline regression analysis were included. Statistical significance was declared at p < 0.05. SigmaPlot 12.5 software (Systat Software, San Jose, CA, USA) was used for all statistical analyses.
Characteristics of the study population
Fifty-two female and one male patient were studied. Their median age at baseline was 46 (IQR 36 – 51 years). The sample comprised 18 Caucasian, 18 African American, 7 Native American, 7 Hispanic, and 3 Asian patients. The median time between visits was 1.8 months (IQR 1.1–3.3 months). Disease activity data are shown in Table 1. Patients with active disease in this routine lupus clinic setting most often had arthritis or mucocutaneous manifestations. Arthritis and rash as defined by the SLEDAI were present in 41 (77.4 %) and 14 (26.4 %) patients, respectively, at baseline. Thirty (56.6 %) and thirteen (24.5 %) patients, respectively, had A or B scores on the BILAG musculoskeletal and/or mucocutaneous domains at baseline. At baseline, 21 (39.6 %) patients had low C3 and/or C4 complement levels, and 10 (18.8 %) had positive results for anti-double-stranded DNA (anti-dsDNA) antibodies by Crithidia luciliae assay. Most patients were receiving hydroxychloroquine (Table 1), and about one-fourth were receiving chronic corticosteroids. Azathioprine, methotrexate, or mycophenolate mofetil was prescribed in 15–28 % of the patients, alone or in combination with hydroxychloroquine and prednisone.
Parameters of HRV did not significantly change between baseline and follow-up (see Table S2 in Additional file 1). Exercise and exposure to tobacco, caffeine, and alcohol were similar between visits. Few patients had comorbidities known to affect HRV, including two patients with heart failure, two with recent seizures, one with untreated depression, one with diabetes, and one with asthma not controlled while on medications. Similar results were obtained when these seven individuals were omitted from the analysis; therefore, they were included in the final results. Exposure to medications known to increase HRV (β-blockers, calcium channel blockers, and clonidine) or decrease HRV (tricyclic antidepressants, antihistamines, β-agonists) did not differ between baseline and follow-up (Table 1). At baseline, seven patients were taking at least one medication known to increase HRV, eight were taking at least one known to decrease HRV, and five were taking combinations from both groups.
HRV is inversely associated with disease activity and flare at baseline
Relationships of HRV with disease activity were examined at baseline by univariate linear regression (Table 2). Total BILAG disease activity scores were inversely related to HRV parameters [pNN50 (p = 0.019), HF power (p = 0.020), LF/HF (p = 0.024)], and similar associations were observed when we focused on the most commonly involved organ, the musculoskeletal system [pNN50 (p = 0.043), HF power (p = 0.010)]. Total SLEDAI scores and arthritis by SLEDAI had a similar trend of association with HRV [LF/HF (p = 0.073 for total SLEDAI), HF power (p = 0.088 for SLEDAI arthritis)]. HRV was not associated with the presence of classic serological markers such as anti-dsDNA antibodies and low complement levels. The LF/HF ratio was associated with disease activity by BILAG (p = 0.024) and SFI (p = 0.008), although a trend was also seen in SLEDAI (p = 0.073) and PGA (p = 0.062) (Fig. 1). No association of HRV with tobacco or caffeine exposure or medications affecting autonomic nervous system function was demonstrable in this analysis. At the follow-up visits, no association of HRV with disease activity was observed.
HRV is inversely related to soluble TNFRII and to IFN-induced cytokines
Exploration of associations between HRV and cytokine levels at baseline was performed using linear regression without adjustment for multiple comparisons (Table 3). Elevated plasma sTNFRII and MIG were consistently associated with decreased HRV across all parameters measured (Table 3 and Fig. 2). Inverse associations with one or more parameters of HRV were also observed for plasma IL-1RA, SCF, and IL-15 (Table 3). Trends of similar associations were noted for IFN-α, IP-10 and serum BLyS (p < 0.10). There was no evidence to suggest that plasma IFN-γ, IL-10, IL-12, IL-1β, IL-1α, IL-2, IL-8, IL-17A, IL-21, IL-23, MIP-1α, sCD40L, TNFα, sTNFRI, PAI-1, resistin, TGF-β, or serum APRIL were related to HRV.
Considering the potentially cyclical interplay between HRV and cytokines , we further explored which inflammatory pathway is most closely linked as a potential predictor of HRV in a multivariate model including plasma IL-1RA, SCF, IFN-α, IP-10, and serum BLyS. Associations of HRV with plasma MIG or sTNFRII remained significant (MIG: RMSSD p = 0.036, HF power p < 0.001, LF/HF p = 0.003; sTNFRII: RMSSD p = 0.011, HF power p = 0.040) (Table 4). In a correlation matrix of all parameters included in the models, there was a strong relationship between sTNFRII and MIG (Spearman’s rank correlation R = 0.751, p < 0.001); thus, these were not combined in the same model. Collinearity was also evident between sTNFRII and SCF (R = 0.508, p < 0.001), between MIG and IP-10 (R = 0.585, p < 0.001), and between MIG and SCF (R = 0.462, p < 0.001), as well as between IP-10 and IL-1RA (R = 0.564, p = 0.001), whereas all other R values were <0.4 (Table S3 in Additional file 1). Plasma sTNFRII was weakly correlated with plasma TNF-α (R = 0.357, p = 0.009) and sTNFRI levels (R = 0.273, p = 0.048).
Some cytokines associated with HRV are also associated with disease activity
We further examined associations of baseline disease activity with plasma and serum cytokines by univariate linear regression (see Table S4 in Additional file 1). Plasma MIG was related to BILAG (p = 0.017) and PGA (p = 0.024), but no association was observed for sTNFRII. Several other cytokines appeared to be potentially associated with one or more measures of disease activity, stressing the different pathways that might be related to disease activity in different subsets of patients with SLE, only some of which might have an impact on HRV.
Changes in HRV between visits are associated with changes in SLEDAI
Changes in SLEDAI between visits showed a significant inverse association with changes in RMSSD (p = 0.007), and there was also a trend of association with changes in pNN50 (p = 0.094) (Table 5). Trends of association (p > 0.05) of BILAG and PGA with RMSSD and HF power were also observed.
Changes in the LF/HF ratio are related to changes in cytokines
Using linear regression, we examined changes in HRV in relation to changes in those cytokines that had been inversely related to HRV at baseline (Table 5). Changes of the LF/HF ratio were correlated with similar changes across all of those cytokines, including sTNFRII (p = 0.021) and MIG (p = 0.003). Changes in HF power were inversely related to changes only in serum BLyS (p = 0.036), plasma IFN-α (p = 0.004), and IL-15 (p = 0.012).
Although autonomic dysregulation is well documented in SLE, much remains to be learned about its direct or indirect relationships to disease activity. Although abnormalities in cardiovascular reflex testing and HRV have been described in SLE previously, these did not correlate with disease activity measured by BILAG  or SLEDAI [4, 6, 8]. In the present study, we found potential relationships between HRV and disease activity measured by SLEDAI, BILAG, and PGA, as well as the presence or absence of flare (Table 2). Furthermore, HRV changes between visits (time domain analysis and HF power) were correlated with changes in disease activity across all indices (Table 5). At follow-up, disease activity was generally improved, and the range of disease was far narrower in the population (Table 1), deterring the determination of cross-sectional linear relationships at those visits. These results are consistent with the known anti-inflammatory properties of vagus signaling  and raise the possibility that HRV might be a noninvasive marker for SLE disease activity and improvement, at least in a relevant subset of patients.
To begin a process by which such patients could be identified, we explored the association of HRV with a large panel of cytokines that was previously found at our institution to be relevant to SLE  (arrayed in Table S1 in Additional file 1). HRV was examined not only as a predictor (Tables 3 and 5) but also as an outcome (Table 4) of peripheral inflammation. Although the exact relationship between these variables is not yet known, evidence in both directions supports the hypothesis of a cyclical relationship . This was an exploratory study, which limits the conclusions that can be drawn. Given the exploratory nature of our study and our small sample size, we elected not to adjust for multiple comparisons, in order to be able to detect associations of modest strength . Notably, several strong preliminary associations were observed that provide useful hypotheses for further testing. In analysis that was not adjusted for multiple comparisons or other covariates, all HRV measurements were predictive of plasma MIG, an IFN type II-related cytokine, and sTNFRII, whose shedding has been associated with lupus flares  (Tables 3 and 5). HRV was also predictive of plasma SCF, a major macrophage growth and differentiation factor , while consistent inverse relationships with HRV were observed for two other IFN-related cytokines, IP-10 and IFNα, as well as for IL-1RA and serum BLyS (Tables 3 and 5). Although no relationship of HRV with plasma IFN-γ or TNF-α was noted in this analysis, our results support a potential association of vagal signaling with a Th1-mediated myeloid cell pathway, known to be important in subsets of patients with lupus . The inverse association of HRV with plasma IL-15, a cytokine implicated in intestinal autoimmunity , is also intriguing because vagal signaling is known to reach myenteric neurons in close contact with intestinal macrophages . Certain cytokines previously associated with HRV in experimental models of acute inflammation, such as TNF-α and IL-1β , were not impacted by HRV in this study of patients with lupus. Larger studies, however, might illuminate whether there is a subset of patients where these relationships could be observed or whether there is some reason why these are not linked to HRV or myeloid pathways in lupus. Cytokines can also differ between the systemic circulation and the tissue microenvironment, and autonomic regulation at the tissue level can be discordant with systemic autonomic effects. These considerations further limit the conclusions that can be drawn from this preliminary study.
Clinically evident arthritis was the most prevalent SLE manifestation in our cohort, and arthritis alone was associated with decreased parasympathetic activity. A direct effect of vagal signaling in the joints is unlikely; nonetheless, this finding might be clinically meaningful, considering the promising role of intermittent electrical vagal nerve stimulation in animal models of arthritis  and in patients with rheumatoid arthritis .
Multiple potentially confounding factors can affect HRV, including age, cardiovascular disease, diabetes, renal failure, obstructive sleep apnea, tobacco, caffeine, and drugs [34–36]. Although few patients in this study had such obvious features, cardiovascular disease and prediabetes are often subclinical in SLE [37, 38], so further evaluation of these potential confounders is in order. No association with smoking or caffeine consumption was evident, but more sensitive measures to record and quantify the degree of these exposures might have captured an effect. It is also possible that HRV is more sensitive to psychosocial factors than lupus disease activity itself, especially in patients with mild to moderate disease. We did not formally assess depression, pain, sleep disturbances, and fatigue, all of which are common in SLE and known to be associated with decreased HRV [39, 40]. A more comprehensive assessment of autonomic nervous system function by additional clinical cardiovascular reflex testing, as well as assessment of baroreflex sensitivity and respiratory sinus arrhythmia, might have also allowed better discernment of these effects.
Frequency domain parameters correlated better than time domain indices with disease activity and/or cytokines in this study. This is somewhat expected, as we used only 5-minute ECG recordings, and time domain analysis has higher variability with shorter durations . In addition to the limitations of HRV analysis , those pertinent to the clinical measures of lupus activity should be acknowledged to prevent overinterpretation of these findings .
Our results suggest that impaired HRV, particularly the LF/HF ratio, is associated with lupus disease activity and several cytokines related to IFN type II and TNF pathways. The cytokines with the strongest association were MIG and sTNFRII, confirming and expanding previous connections between vagal signaling and immunity. These data support a potentially important interplay between autonomic nervous system function and lupus disease activity, as well as a possible role of vagal modulation as adjunctive therapy in SLE. This can be approached by pharmacological or electrical treatments or through cognitive behavioral therapy approaches . Such therapies are rapidly entering clinical research programs and may be applicable to SLE and other clinical inflammatory disorders.
A Proliferation-Inducing Ligand
British Isles Lupus Assessment Group
B Lymphocyte Stimulator
Health Insurance Portability and Accountability Act
Heart Rate Variability
Interleukin 1 Receptor Antagonist
Interferon Gamma-Induced Protein 10
Institutional Review Board
Low Frequency Power
Low Frequency to High Frequency ratio
Monokine Induced by Interferon Gamma
Macrophage Inflammatory Protein 1 Alpha
Normal to Normal
Oklahoma Medical Research Foundation
Plasminogen Activator Inhibitor 1
Physician Global Assessment
Percentage of Pairs of Adjacent NN Intervals Differing by more than 50 Milliseconds
Square Root of the Mean of the Squares of Differences between Adjacent NN Intervals
Stem Cell Factor
Safety of Estrogens in Lupus Erythematosus National Assessment
SELENA-SLEDAI Flare index
Systemic Lupus Erythematosus
Systemic Lupus Erythematosus Disease Activity Index
soluble Tumor Necrosis Factor Receptor
Transforming Growth Factor Beta
Tumor Necrosis Factor Alpha
Andersson U, Tracey KJ. Neural reflexes in inflammation and immunity. J Exp Med. 2012;209:1057–68.
Maule S, Quadri R, Mirante D, Pellerito RA, Marucco E, Marinone C, et al. Autonomic nervous dysfunction in systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA): possible pathogenic role of autoantibodies to autonomic nervous structures. Clin Exp Immunol. 1997;110:423–7.
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:97–102.
Harle P, Straub RH, Wiest R, Mayer A, Scholmerich J, Atzeni F, et al. Increase of sympathetic outflow measured by neuropeptide Y and decrease of the hypothalamic-pituitary-adrenal axis tone in patients with systemic lupus erythematosus and rheumatoid arthritis: another example of uncoupling of response systems. Ann Rheum Dis. 2006;65:51–6.
Milovanović B, Stojanović L, Milićevik N, Vasić K, Bjelaković B, Krotin M. Cardiac autonomic dysfunction in patients with systemic lupus, rheumatoid arthritis and sudden death risk. Srp Arh Celok Lek. 2010;138:26–32.
Haensel A, Mills PJ, Nelesen RA, Ziegler MG, Dimsdale JE. The relationship between heart rate variability and inflammatory markers in cardiovascular diseases. Psychoneuroendocrinology. 2008;33:1305–12.
Buyon JP, Petri MA, Kim MY, Kalunian KC, Grossman J, Hahn BH, et al. The effect of combined estrogen and progesterone hormone replacement therapy on disease activity in systemic lupus erythematosus: a randomized trial. Ann Intern Med. 2005;142:953–62.
Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart rate variability: standards of measurement, physiological interpretation and clinical use. Circulation. 1996;93:1043–65.
Austin PC, Steyerberg EW. The number of subjects per variable required in linear regression analyses. J Clin Epidemiol. 2015;68:627–36.
Kaur K, Zheng NY, Smith K, Huang M, Li L, Pauli NT, et al. High Affinity antibodies against influenza characterize the plasmablast response in SLE patients after vaccination. PLoS One. 2015;10:e0125618.
Karonitsch T, Feierl E, Steiner CW, Dalwigk K, Korb A, Binder N, et al. Activation of the interferon-gamma signaling pathway in systemic lupus erythematosus peripheral blood mononuclear cells. Arthritis Rheum. 2009;60:1463–71.
Matteoli G, Gomez-Pinilla PJ, Nemethova A, Di Giovangiulio M, Cailotto C, van Bree SH, et al. A distinct vagal anti-inflammatory pathway modulates intestinal muscularis resident macrophages independent of the spleen. Gut. 2014;63:938–48.
Levine YA, Koopman FA, Faltys M, Caravaca A, Bendele A, Zitnik R, et al. Neurostimulation of the cholinergic anti-inflammatory pathway ameliorates disease in rat collagen-induced arthritis. PLoS One. 2014;9:e104530.
Doria A, Shoenfeld Y, Wu R, Gambari PF, Puato M, Ghirardello A, et al. Risk factors for subclinical atherosclerosis in a prospective cohort of patients with systemic lupus erythematosus. Ann Rheum Dis. 2003;62:1071–7.
We thank Dr. David Albert for technical assistance with acquisition of ECG tracings for HRV analysis.
Patient procedures and cytokine analyses were funded by general funds of the Arthritis and Clinical Immunology Research Program at the OMRF. The authors were not compensated for their contributions to this study.
AT led the acquisition and analysis of data and drafted the manuscript. SS, JWD, MEM, JAJ, and JTM assisted in acquisition of data. All authors participated in data analysis and reviewing, editing, and approval of the final manuscript.
The authors declare that they have no competing interests.
Authors and Affiliations
Arthritis and Clinical Immunology Research Program, Oklahoma Medical Research Foundation, 825 N.E. 13th Street, Oklahoma City, OK, 73104, USA
Aikaterini Thanou, Melissa E. Munroe, Judith A. James & Joan T. Merrill
Heart Rhythm Institute, University of Oklahoma Health Sciences Center, 1200 Mark Everett Drive, TCH 6E103, Oklahoma City, OK, 73104, USA
Stavros Stavrakis & John W. Dyer
Departments of Medicine and Pathology, University of Oklahoma Health Sciences Center, Oklahoma City, OK, USA
Table S1. Plasma cytokines measured at each visit by multiplex bead-based immunoassay. Table S2. Parameters of HRV (RMSSD, pNN50, HF power, and the LF/HF ratio) at baseline and follow-up. HRV parameters are compared between visits by paired analyses (Wilcoxon signed-rank test or paired t test); p values are listed. Table S3. Matrix of Spearman’s rank correlations (R) between cytokines included in the multivariate models. p values <0.1 are included, and R values that are statistically significant (p < 0.05) are shown in bold type. Table S4. Associations of disease activity and flare with plasma cytokines at baseline by univariate linear regression (p values and beta regression coefficients). (DOCX 23 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.
Thanou, A., Stavrakis, S., Dyer, J.W. et al. Impact of heart rate variability, a marker for cardiac health, on lupus disease activity.
Arthritis Res Ther18, 197 (2016). https://doi.org/10.1186/s13075-016-1087-x