- Research article
- Open Access
- Open Peer Review
Proarrhythmic risk and determinants of cardiac autonomic dysfunction in collagen-induced arthritis rats
BMC Musculoskeletal Disordersvolume 17, Article number: 491 (2016)
Patients with rheumatoid arthritis (RA) have increased risk of sudden cardiac death (SCD), which is two-fold higher than general population. The driving cause of SCD was considered due to lift-threatening arrhythmia where systemic inflammation acts as the pathophysiological basis linking RA to autonomicdysfunction.
To assess the sympathetic over-activity of “inflammatory reflex”, we measured heart rate variability (HRV) in a rat collagen-induced arthritis (CIA) model, whose arthritis is induced in Lewis rats by intradermal injection of emulsion of type II collagen. Single-lead electrocardiogram (ECG) was recorded for 30 min every two days. Time and frequency-domain parameters, detrended fluctuation analysis (DFA), deceleration (DC) and acceleration capacity (AC) were analyzed.
Compared with 9 control rats, many of HRV parameters of 9 CIA rats revealed significant different. At the beginning of arthritis, LF/HF was significant higher than controls (1st week: 2.41 ± 0.7 vs. 1.76 ± 0.6, p < 0.05; 2nd week: 2.24 ± 0.5 vs. 1.58 ± 0.5, p < 0.05) indicating intensive inflammatory reflex at the initial phase of inflammation but no significant difference was observed in the following recover phase. The similar trend of DFA parameters was noted. However, the DC appeared progressive lower despite of no significant increase of the LF/HF compared with controls since 4th week.
We observed sympathetic over-activation of inflammatory reflex during early stage of arthritis in CIA rats. The ongoing decline of DC indicated advanced cardiac autonomic dysfunction regardless of remission of acute arthritis.
Rheumatoid arthritis (RA) is a common autoimmune disease that is associated with progressive disability, systemic complications, early cardiac death and socioeconomic costs. RA is characterized by loss of tolerance to self-proteins that contain a citrulline residue, where anti-citrulline response can be initiated with T-cell, B-cell and tumor necrosis factor α (TNF-α). Synovitis is launched and perpetuated by positive feedback loops and in turn promotes systemic disorders, including myocardial infarction, stroke and heart failure . The cardiac mortality is 1.5-fold higher and sudden cardiac death (SCD) is twice as likely among RA patients than general population. Although the main focus of clinical setting is on accelerated atherosclerosis and myocardial injury, chronic systemic inflammation may increase arrhythmogenicity via other mechanisms and contribute to higher risk of SCD . Synovitis derived inflammatory cytokines (interleukin-1, interleukin-6 and TNF-α) could induce non-structural myocardium change via direct electrophysiological effect on cardiomyocyte and sympathetic overactivation . Many studies demonstrated inflammatory cytokines could prolong action potential duration via changing expression and function of potassium and calcium channel [4–6]. On the other hand, signs of autonomic dysfunction in RA population have been documented decades ago and involved cardiovascular nervous system . Heart rate variability (HRV) is a non-invasive, practical and reproducible method to assess the effect of sympathovagal balance on the heart. In general population, reduced HRV indicates increased sympathetic tone and decreased parasympathetic system being associated with higher risk of cardiovascular events . Several studies investigating HRV in RA patients, reported a depression of time- and/or frequency-domain measures of HRV, implying an increase in the sympathetic drive of the heart rate in this population. We intend to test the hypothesis that “inflammatory reflex” could be a mechanism leading to heart impairment in RA .
Collagen-induced arthritis (CIA) in rat is an animal model for rheumatoid arthritis (RA) and can be induced in Lewis rat. The CIA model can be used to unravel mechanisms involved in the development of arthritis and frequently used to study the effect of new therapeutics.
Eighteen adult Lewis rats weighing 160–180 g were used. The animals were housed under special pathogen-free conditions at the animal facility. These rats were divided into 2 groups: CIA and control group equally. The healthy condition and physical activity of the rats were monitored daily, and the body weight and paw thickness were measured and recorded every week.
Induction of arthritis
Chick Type II collagen (CII) (Chondrex cat. #20012) was dissolved in 0.05 M acetic acid by gently stirring overnight at 4 °C. Equal amount collagen and complete Freud’s Adjuvant (CFA) (Chondrex cat. #7024) were mixed in an ice-water bath, adding the collagen drop-wise to the CFA while mixing. The same procedure was conducted when mixing Bovine Type II collagen (Chondrex cat. #20022) with Incomplete Freud;s Adjuvant (IFA) (Chondrex cat.#7002) On day 0, 0.5 ml of the emulsion (containing 0.5 mg Chick CII + CFA) was injected at the base of the tail of each rat. On day 7, a second injection (0.5 mg Bovine CII + IFA) was administered in the same way. For arthritis assessment, all rats were monitored three times a week by the same person blinded to the treatment group, and the clinical scores were evaluated. Zero score represents no edema or swelling in the joint; 1 score represents slight edema and erythema limited to the foot and/or ankle; 2 scores represent edema and erythema toes and most joints of the ankles; 3 scores represent severe edema and erythema paw below ankle joint; 4 scores represent edema and erythema of all paws including ankle joint . The cumulative score for all four paws of each rat (maximum possible score of 16) was used to represent the overall disease severity and progression. The model of arthritis was considered successful when the scores were greater than 5.
NN Interval recordings
The ECG of CIA rats was recorded for 30 min every two days (CIA rats and controls). Prior to ECG recordings, animals were conditioned for 7 consecutive days, 30 min in each morning (0700–0900 h) on a self-made platform (Additional file 1: Figure S1). The self-made platform was designed to prevent rats to strip off recording electrodes in non-anesthetized condition. A day before ECG was firstly recorded, the ventral thoracic region of each animal was carefully shaved. We placed the negative electrode near the right shoulder and the positive electrode to the left of the xyphoid space to obtain robust P, R, S, and T waves, as human lead II configuration. We kept the rat awaken and recorded ECG signal using a bioamplifier. The ECG signal was analyzed using custom-made LabVIEW software which automatically detected the R-waves. The QRS complexes were automatically and manually classified as normal sinus rhythm, atrial and ventricular premature beats. The term “NN interval”, or the so-called normal-to-normal intervals isthe intervals between adjacent QRS complexes resulting from sinus node depolarizations. ECGs were bandpass filtered (2–300 Hz) and, after R wave peak detection, 180-s tachograms were generated. The RR intervals were deduced from adjacent normal sinus beats and transferred to a personal computer to be post-processed by a Matlab program (Kubios HRV) .
Time and frequency domain parameters
The mean heart rate, standard deviation of N-N intervals (SDNN) and root mean square of successive differences of N-N interval (RMSSD) were used as time-domain measures of HRV. The power spectrum densities were estimated by Welch’s averaged periodogram method .
The LF power (0.04 to 0.15 Hz) and HF power (0.15 to 0.4 Hz) were derived from the sum of area within specific frequency range under the power spectrum density curve of the entire 30-min segment. The normalized unit (NU) of HF (HF NU) or LF (LF NU) was calculated as the ratio of the absolute powers of HF or LF to TP and multiplied by 100.
Generally, SDNN reflects overall heart rate fluctuation. RMSSD could represent parasympathetic tone and also be sensitive to uneven beat detection. In terms of frequency domain, LF is considered to reflect the combined modulation of efferent vagal and efferent sympathetic nervous system activity and HF is related to modulation of efferent vagal activity by ventilation. LF/HF ratio is often claimed to characterize “sympathovagal balance” or “relative sympathetic activity” because the LF reflects modulation by both the sympathetic and parasympathetic arms of the autonomic nervous system and the HF band reflects parasympathetic activity .
Detrended fluctuation analysis
Detrended fluctuation analysis (DFA) quantified the correlation properties of fractal-like dynamics caused by complex interplay between vagal and sympathetic heart rate regulation . The DFA technique was used to quantify the fractal scaling properties of short- and intermediate-term R-R interval time series. Clinically, reduced DFA might be more precise in predicting fatal arrhythmic events than that based on traditional methods . In this study, both the short-term (DFAα1, 4 to 11 beats) and long-term (DFAα2, >11 beats) scaling exponents were calculated. All the analyses were performed by using software developed in-house provided by Matlab 7.8 (Mathworks, Inc., Natrick, Ma, USA).
Deceleration and acceleration capacity
We used phase-rectified signal averaging method (PRSA) to process sequences of R-R interval obtained from our 30-min recordings . The PRSA software was benchmarked against the freely available download from the Technical University of Munich. Deceleration capacity (DC) quantified spontaneous increases in NN intervals and acceleration capacity (AC) quantified spontaneous decrease in NN intervals. PRSA method could extract periodicities from complex time series, including noise, artefact, non-stationarities as well as periodic component. Non-periodic components were eliminated thereafter. Consequently, the PRSA computation is relatively robust against artefact and ectopic beats with extensive editing of ECG recording . DC quantifies the increases in NN intervals by the phase-rectified signal averaging method, reflecting the capability to slow heart rate and AC characterized the ability to speed heart rate .
The data were expressed as mean ± standard deviation. A two-way, repeated-measures analysis of variance (ANOVA) was used to evaluate statistical differences between groups (control versus CIA group) and among follow-up durations (from 1st to 6th week). A p < 0.05 was considered statistically significant and the alpha level was adjusted during multiple comparisons so as to maintain the rate of Type I error at 5% during Holm-Bonferroni post hoc analysis. When differences were found in main effect, an independent-samples t test was used to analyze differences between two groups at each measurement time. All analyses were performed with SPSS 19.0 (SPSS Inc. Chicago, IL).
The rats in control group were healthy and without inflammation in joints (Fig. 1a) while the rats treated with collagen were unhealthy and had progressive swelling and redness appearance in the paws (Fig. 1b). The swelling of joints in the CIA rats, especially hind paws, was gradually aggravated and approached the peak in the day 18 (Fig. 1c and d). After that, the swelling was alleviated gradually but the joints became stiff and the movements of the limbs were limited.
The HRV parameters are shown in Table 1. The two-way repeated measure ANOVA was conducted to find out the difference HRV between groups and measurement times. It showed a significant difference among follow-up duration (p < 0.001), and between two groups (p < 0.001). Significant interaction between the groups and measurement times was shown (p < 0.05). In terms of liner HRV, heart rate (458 ± 35 vs. 386 ± 25, p = 0.002) and SDNN (8.23 ± 1.6 vs. 4.53 ± 2.5 ms, p = 0.004) of CIA rats were significant higher than control group in the first week. These changes persisted till 4th week (2nd week, heart rate: 438 ± 31 vs. 411 ± 57, p = 0.02; SDNN: 12.28 ± 4.3 vs. 6.08 ± 3.4 ms, p = 0.002 and 3rd week, heart rate: 471 ± 32 vs. 411 ± 6, p = 0.003; SDNN: 5.82 ± 2.7 vs. 3.9 ± 0.9 ms, p = 0.03) when the inflammation of joins in CIA rats was relieved clinically. In frequency-domain parameters, LF (1st week: 69.48 ± 5.7 vs. 61.71 ± 9.3, p = 0.03; 2nd week: 68.41 ± 5.1 vs. 60.00 ± 7.2, p = 0.02 and 3rd week: 68.78 ± 10.1 vs. 62.73 ± 11.1, p = 0.002) was significant higher in CIA rats during the first three weeks. On the other hand, HF (1st week: 30.34 ± 5.6 vs. 38.01 ± 9.1, p = 0.02 and 2nd week: 31.46 ± 5.0 vs. 39.85 ± 7.2, p = 0.03) was significant higher in control group during the first two weeks. In the same way, LF/HF (1st week: 2.41 ± 0.7 vs. 1.76 ± 0.6, p = 0.01 and 2nd week: 2.24 ± 0.5 vs. 1.58 ± 0.5, p = 0.02) of CIA rats was significant higher than control group. During 4th to 6thweek, when arthritis was relieved gradually in CIA rats, the frequency-domain parameters had no significant difference compared with control group (Fig. 2). In the DFA parameters, DFAα1 (1st week: 0.88 ± 0.2 vs. 0.62 ± 0.2, p = 0.001; 2nd week: 0.82 ± 0.2 vs. 0.74 ± 0.3, p = 0.02) and DFAα2 (1st week: 1.12 ± 0.1 vs. 0.71 ± 0.4, p = 0.002; 2nd week: 1.24 ± 0.1 vs. 1.09 ± 0.3, p = 0.03) of CIA rats were significant higher in the first two weeks. Similarly, there was no further change in the following weeks. With respect to non-linear parameters, we observed lower DC (1st week: 9.5 ± 2.4 vs. 13.7 ± 2.7, p = 0.03; 2nd week: 9.8 ± 2.6 vs. 11.1 ± 3.8, p = 0.03; 3rd week : 7.5 ± 2.8 vs. 7.6 ± 0.7, p = 0.51; 4th week: 6.9 ± 1.1 vs. 9.6 ± 2.9, p = 0.02; 5th week: 8.2 ± 1.4 vs. 6.9 ± 1.6, p = 0.02 and 6th week: 7.3 ± 2.2 vs. 8.5 ± 3.5, p = 0.03) in CIA rats compared with control group and there was progressive decline of DC values during whole follow up period (Fig. 3, upper panel). Among AC parameters, there was opposite trend with higher values in CIA rats (1st week: -7.1 ± 2.6 vs. -10.6 ± 2.4, p = 0.02; 2nd week: -7.7 ± 2.0 vs. -8.1 ± 6.3, p = 0.13; 3rdweek : -4.1 ± 3.1 vs. -5.2 ± 0.4, p = 0.42; 4th week: -4.8 ± 0.9 vs. -7.3 ± 2.7, p = 0.03; 5th week: -5.7 ± 1.7 vs. -6.1 ± 1.4, p = 0.04 and 6th week: -5.2 ± 1.8 vs. 6.3 ± 3.1, p = 0.33) than control (Fig. 3, lower panel)
In this study, we examined various HRV parameters in a CIA rat model, which mimicking the joint inflammation of rheumatoid arthritis. We demonstrated higher time-domain parameters, LF and LF/HF ratio during acute arthritis stage (1st to 3rd week). There was no further difference in linear HRV parameters when the arthritis was relieved among 4th to 6th week. On the other hand, we observed progressive incapability of slowing and speeding heart rate in CIA rats, even after relief of joints inflammation. This is the first study to elucidate the temporal change of linear and non-linear HRV parameters in CIA rats and indicate that DC and AC could be potential surrogates of autonomic system dysregulation in patients with rheumatoid arthritis.
CIA is an induced model of rheumatoid arthritis that bears several similarities to human disease, including female vulnerability, symmetrical synovitis and inflammatory-mediated damage to the cartilage and bone of joints [18, 19]. In human RA, autonomic dysfunction were frequently observed and associated with an increased mortality risk, [7, 20] characterized by reduced HRV, an expression of sympathetic over-activity and para-sympathetic withdraw [9, 21–23]. Most of these cross-sectional studies found reduced time and frequency domain HRV parameters in RA patients when compared with control . Because patients with rheumatoid arthritis had long-term disease course and several times of remission and relapsing, it was unlikely to record the temporal change of autonomic system. On the other hand, studies evaluating HRV in RA animal models are very rare. Utilizing continuous ECG recording in CIA rats and HRV analyses, our study observed the high sympathetic outflow in the acute stage of arthritis and ongoing damage of cardiac autonomic function despite of remission of arthritis.
Inflammatory cytokines and local heat and pain causes sympathetic activation by targeting the autonomic centers of the brain, consequently dampening further cytokine production and immune-inflammatory activation [25, 26]. The so-called inflammatory reflex is a self-controlling loop and an adaptive response to inflammation. However, sympathetic over-activation also affect the heart, possibly favoring onset of cardiac arrhythmia . In our study, during the first three weeks after the induction, we observed increased heart rates, indicating cardiac sympathetic activation due to arthritis. There was higher SDNN in CIA rats, reflecting larger heart rate fluctuation [13, 28]. Furthermore, higher LF and LF/HF ratio in CIA rats may also represent activation of sympathetic tone. These results were opposite to most clinical trials evaluating HRV in RA patients, which revealed reduced HRV parameters. It may be because RA patients in most clinical trials had arthritis many times with chronic disease course. However, arthritis was firstly induced in our rats. In the following 4th to 6th week, when arthritis was relieved, the time and frequency domain parameters of CIA rats declined gradually, indicating returning to baseline sympathovagal activity.
Detrended fluctuation analysis method measures the qualitative characteristics and correlation features of heart rate behavior. Heart rate time series are fractal because they display self-similar (scale-invariant) fluctuations over a wide range of time scales . These fluctuations are determined by a delicate interplay between sympathetic and vagal outflow . The significant increased DFAα1 and DFAα2 in our CIA rats during first two weeks reflected increased sympathetic activity accompanied with withdraw of vagal activity . Likewise, after gradual relief of inflammation in joints, both parameters had no difference compared with the control group.
As previously noted in method section, we used PRSA technique to calculate DC and AC, meaning the capability of slowing and speeding heart rates. Autonomic heart-rate modulation due to specific regulation, such as sympathovagal activity, respiratory, baroreflex mediated and circadian, occur on different timescale. PRSA technique could synchronize all periodic components of the signal irrespective of their frequency and timescale. It thus integrated all contributions of above mentioned mechanisms by accumulation the amplitude at the center of PRSA signal. Briefly, DC indicates the overall deceleration capacity of sinus rhythm, without being linked necessarily to one particular physiological process . In our study, we found the progressive decline of DC despite of the recovery of arthritis. One possible explanation was that other pathophysiologic mechanisms persistently affected the SA node function in CIA rats regardless of diminished impact of inflammation reflex or sympathetic activity indicated by HRV. Noteworthy, by using short-duration ECG recording rather than long-term recording, the DC measurement would not have to account for various level of activity of CIA rats . There was similar trend of progressive increase of AC but not achieving significant difference. On the other hand, our findings support that DC might be a more sensitive parameter than HRV and DFA in revealing a residual activation of the inflammatory process . The residual elevated mean clinical score and paw thickness were still noted 6 weeks after the induction (Fig. 1, panels c and d, day 42), suggesting the persistence to a certain degree of inflammatory activation.
Previous studies suggested inflammatory cytokines, particularly TNF-α, could promote cardiac sympathetic activation and arrhythmic events and treatment with TNF-α antagonist could alleviate sympatho-vagal imbalance [33, 34]. HRV impairment was significantly associated with disease duration, activity and inflammation markers [35, 36]. However, few study evaluated the association of HRV impairment and cardiovascular outcomes in RA population , although growing evidences indicate increased cardiovascular mortality [38, 39]. In addition to conventional risk predictors of HRV proposed by the guideline , DC is a good predictor of morality of survivors after myocardial infarction . We also demonstrated that DC is a sensitive markers of cardiac autonomic dysfunction in CIA rats.
There were some limitations of this study. Firstly, our rats were not monitored everyday, leaving some data gaps although our study investigated the temporal change of HRV in CIA rats. Secondly, the conclusions derived from the relatively small numbers of control and CIA rats may require careful consideration of the assumption. Thirdly, the severity of arthritis was based on clinical score measurement, lacking histopathology of joints and inflammatory cytokines from CIA rats. Similarly, we evaluate the activity of sympathetic system by using conventional HRV parameters, lacking baroreflex sensitivity or skeletal muscle sympathetic nerve activity.
Our study demonstrated for the first time the dynamic change of sympathetic activity by using HRV analyses. After relief of arthritis and decline of sympathetic activity, the damage to cardiac autonomic function persisted which reflected by the reduced DC. Using these observations, more in-depth studies could be carried on to evaluate the association of autonomic dysfunction and cardiovascular events in human RA population.
Detrended fluctuation analysis
Heart rate variability
Phase-rectified signal averaging method
Sudden cardiac death
McInnes IB, Schett G: The pathogenesis of rheumatoid arthritis. N Engl J Med. 2011, 365(23):2205-2219. doi: 2210.1056/NEJMra1004965.
Sattar N, McCarey DW, Capell H, McInnes IB. Explaining how “high-grade” systemic inflammation accelerates vascular risk in rheumatoid arthritis. Circulation. 2003;108(24):2957–63.
Lazzerini PE, Capecchi PL, Laghi-Pasini F. Systemic inflammation and arrhythmic risk: lessons from rheumatoid arthritis. Eur Heart J. 2016;1-13. [Epub ahead of print].
Li YH, Rozanski GJ. Effects of human recombinant interleukin-1 on electrical properties of guinea pig ventricular cells. Cardiovasc Res. 1993;27(3):525–30.
Wang J, Wang H, Zhang Y, Gao H, Nattel S, Wang Z. Impairment of HERG K(+) channel function by tumor necrosis factor-alpha: role of reactive oxygen species as a mediator. J Biol Chem. 2004;279(14):13289–92. Epub 12004 Feb 13217.
Hagiwara Y, Miyoshi S, Fukuda K, Nishiyama N, Ikegami Y, Tanimoto K, Murata M, Takahashi E, Shimoda K, Hirano T, et al. SHP2-mediated signaling cascade through gp130 is essential for LIF-dependent I CaL, [Ca2+]i transient, and APD increase in cardiomyocytes. J Mol Cell Cardiol. 2007;43(6):710–6. Epub 2007 Sep 2019.
Stojanovich L: Autonomic dysfunction in autoimmune rheumatic disease. Autoimmun Rev 2009, 8(7):569-572. doi: 510.1016/j.autrev.2009.1001.1018. Epub 2009 Feb 1013.
La Rovere MT, Pinna GD, Maestri R, Mortara A, Capomolla S, Febo O, Ferrari R, Franchini M, Gnemmi M, Opasich C, et al. Short-term heart rate variability strongly predicts sudden cardiac death in chronic heart failure patients. Circulation. 2003;107(4):565–70.
Janse van Rensburg DC, Ker JA, Grant CC, Fletcher L: Autonomic impairment in rheumatoid arthritis. Int J Rheum Dis. 2012, 15(4):419-426. doi: 410.1111/j.1756-1185X.2012.01730.x. Epub 02012 Apr 01727.
Bevaart L, Vervoordeldonk MJ, Tak PP. Collagen-induced arthritis in mice. Methods Mol Biol. 2010;602:181–92. doi:10.1007/1978-1001-60761-60058-60768_60711.
Tarvainen MP, Niskanen J-P, Lipponen JA, Ranta-aho PO, Karjalainen PA. Kubios HRV – Heart rate variability analysis software. Comput Methods Prog Biomed. 2014;113(1):210–20.
Welch PD. The Use of Fast Fourier Transform for the Estimation of Power Spectra: A Method Based on Time Averaging Over Short, Modified Periodograms. IEEE Trans Audio Electroacoust. 1967;AU-15:70–3.
Marek Malik J, Bigger AT, Camm J, Kleiger RE, Malliani A, Moss AJ, Schwartz PJ. 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-381.
Huikuri HV, Perkiomaki JS, Maestri R, Pinna GD: Clinical impact of evaluation of cardiovascular control by novel methods of heart rate dynamics. Philos Trans A Math Phys Eng Sci. 2009, 367(1892):1223-1238. doi:1210.1098/rsta.2008.0294.
Huikuri HV, Makikallio TH, Peng CK, Goldberger AL, Hintze U, Moller M. Fractal correlation properties of R-R interval dynamics and mortality in patients with depressed left ventricular function after an acute myocardial infarction. Circulation. 2000;101(1):47–53.
Kantelhardt JW, Bauer A, Schumann AY, Barthel P, Schneider R, Malik M, Schmidt G. Phase-rectified signal averaging for the detection of quasi-periodicities and the prediction of cardiovascular risk. Chaos. 2007;17(1):015112.
Axel Bauera JWK, Bundec A, Barthela P, Schneidera R, Malikd M, Schmidta G. Phase-rectified signal averaging detects quasi-periodicities in non-stationary data. Physica A. 2006;364(15):423–34.
Trentham DE. Collagen arthritis as a relevant model for rheumatoid arthritis. Arthritis Rheum. 1982;25(8):911–6.
Hegen M, Keith Jr JC, Collins M, Nickerson-Nutter CL. Utility of animal models for identification of potential therapeutics for rheumatoid arthritis. Ann Rheum Dis. 2008;67(11):1505–15. Epub 2007 Nov 1529.
Schwemmer S, Beer P, Scholmerich J, Fleck M, Straub RH. Cardiovascular and pupillary autonomic nervous dysfunction in patients with rheumatoid arthritis - a cross-sectional and longitudinal study. Clin Exp Rheumatol. 2006;24(6):683–9.
Evrengul H, Dursunoglu D, Cobankara V, Polat B, Seleci D, Kabukcu S, Kaftan A, Semiz E, Kilic M. Heart rate variability in patients with rheumatoid arthritis. Rheumatol Int. 2004;24(4):198–202. Epub 2003 Sep 2011.
Aydemir M, Yazisiz V, Basarici I, Avci AB, Erbasan F, Belgi A, Terzioglu E: Cardiac autonomic profile in rheumatoid arthritis and systemic lupus erythematosus. Lupus. 2010, 19(3):255-261. doi: 210.1177/0961203309351540. Epub 0961203309352009 Dec 0961203309351516.
Vlcek M, Rovensky J, Blazicek P, Radikova Z, Penesova A, Kerlik J, Kvetnansky R, Imrich R. Sympathetic nervous system response to orthostatic stress in female patients with rheumatoid arthritis. Ann N Y Acad Sci. 2008;1148:556–61. doi:10.1196/annals.1410.1026.
Lazzerini PE, Capecchi PL, Acampa M, Galeazzi M, Laghi-Pasini F: Arrhythmic risk in rheumatoid arthritis: the driving role of systemic inflammation. Autoimmun Rev 2014, 13(9):936-944. doi: 910.1016/j.autrev.2014.1005.1007. Epub 2014 May 1027.
Elenkov IJ, Wilder RL, Chrousos GP, Vizi ES. The sympathetic nerve--an integrative interface between two supersystems: the brain and the immune system. Pharmacol Rev. 2000;52(4):595–638.
Tracey KJ. The inflammatory reflex. Nature. 2002;420(6917):853–9.
Piccirillo G, Magri D, Ogawa M, Song J, Chong VJ, Han S, Joung B, Choi EK, Hwang S, Chen LS et al: Autonomic nervous system activity measured directly and QT interval variability in normal and pacing-induced tachycardia heart failure dogs. J Am Coll Cardiol 2009, 54(9):840-850. doi: 810.1016/j.jacc.2009.1006.1008.
Kleiger RE, Stein PK, Bigger Jr JT. Heart rate variability: measurement and clinical utility. Ann Noninvasive Electrocardiol. 2005;10(1):88–101.
Goldberger AL, Amaral LA, Hausdorff JM, Ivanov P, Peng CK, Stanley HE. Fractal dynamics in physiology: alterations with disease and aging. Proc Natl Acad Sci U S A. 2002;99 Suppl 1:2466–72.
Tulppo MP, Kiviniemi AM, Hautala AJ, Kallio M, Seppanen T, Makikallio TH, Huikuri HV. Physiological background of the loss of fractal heart rate dynamics. Circulation. 2005;112(3):314–9. Epub 2005 Jul 2011.
Bauer A, Kantelhardt JW, Barthel P, Schneider R, Makikallio T, Ulm K, Hnatkova K, Schomig A, Huikuri H, Bunde A, et al. Deceleration capacity of heart rate as a predictor of mortality after myocardial infarction: cohort study. Lancet. 2006;367(9523):1674–81.
Campana LM, Owens RL, Clifford GD, Pittman SD, Malhotra A: Phase-rectified signal averaging as a sensitive index of autonomic changes with aging. 1985) 2010, 108(6):1668-1673. doi: 1610.1152/japplphysiol.00013.02010. Epub 02010 Mar 00025.
Lazzerini PE, Acampa M, Hammoud M, Maffei S, Capecchi PL, Selvi E, Bisogno S, Guideri F, Galeazzi M, Pasini FL. Arrhythmic risk during acute infusion of infliximab: a prospective, single-blind, placebo-controlled, crossover study in patients with chronic arthritis. J Rheumatol. 2008;35(10):1958–65. Epub 2008 Aug 1915.
Yu X, Patterson E, Huang S, Garrett MW, Kem DC. Tumor necrosis factor alpha, rapid ventricular tachyarrhythmias, and infarct size in canine models of myocardial infarction. J Cardiovasc Pharmacol. 2005;45(2):153–9.
Yadav RK, Gupta R, Deepak KK. A pilot study on short term heart rate variability & its correlation with disease activity in Indian patients with rheumatoid arthritis. Indian J Med Res. 2012;136(4):593–8.
Lazzerini PE, Acampa M, Capecchi PL, Hammoud M, Maffei S, Bisogno S, Barreca C, Galeazzi M, Laghi-Pasini F: Association between high sensitivity C-reactive protein, heart rate variability and corrected QT interval in patients with chronic inflammatory arthritis. Eur J Intern Med. 2013, 24(4):368-374. doi: 310.1016/j.ejim.2013.1002.1009. Epub 2013 Mar 1019.
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. doi: 210.1016/j.semarthrit.2014.1006.1003. Epub 2014 Jul 1022.
Solomon DH, Karlson EW, Rimm EB, Cannuscio CC, Mandl LA, Manson JE, Stampfer MJ, Curhan GC. Cardiovascular morbidity and mortality in women diagnosed with rheumatoid arthritis. Circulation. 2003;107(9):1303–7.
Maradit-Kremers H, Nicola PJ, Crowson CS, Ballman KV, Gabriel SE. Cardiovascular death in rheumatoid arthritis: a population-based study. Arthritis Rheum. 2005;52(3):722–32.
This work was supported by Grant 105WFA0151174 from National Taiwan University College of Medicine and Ministry of Science and Technology, Grant HCH104-073 from the National Taiwan University Hospital Hsinchu Branch, Grant 104 W970 and 105 W970 from the National Chiao-Tung University, and National Natural Science Foundation of China (No. 81271661).
Availability of data and materials
The dataset supporting the conclusions of this article is included within the article.
TTL proposed the study design and wrote the manuscript. YLS performed the induction of arthritis and acquisition of data. CEW calculated the HRV and statistical analysis. HZ and YBL participated in drafting the article or revising it critically for important intellectual content. SFL made substantial contributions to conception and design and gave final approval of the version to be submitted and any revised version. All of the authors read and approved the final version of the manuscript.
Authors’s information were listed in the title page.
The authors declare that they have no competing interests.
Consent for publication
The animal research protocol was approved by the Institutional Animal Care and Use Committee of the National Chiao-Tung University.
Our study methods were carried out in accordance with ARRIVE guideline.
The self-made platform for ECG recording. (JPG 230 kb)