Scholarly article on topic 'Association of heart rate variability and inflammatory response in patients with cardiovascular diseases: current strengths and limitations'

Association of heart rate variability and inflammatory response in patients with cardiovascular diseases: current strengths and limitations Academic research paper on "Basic medicine"

Share paper
Academic journal
Frontiers in Physiology
OECD Field of science

Academic research paper on topic "Association of heart rate variability and inflammatory response in patients with cardiovascular diseases: current strengths and limitations"



published: 10 July 2013 doi: 10.3389/fphys.2013.00174

Association of heart rate variability and inflammatory response in patients with cardiovascular diseases: current strengths and limitations

Vasilios Papaioannou1*, loannis Pneumatikos1 and Nikos Maglaveras2

' Intensive Care Unit, Aiexandroupoiis General Hospital, Democritus University of Thrace, Aiexandroupoiis, Greece 2 Laboratory of Medical Informatics, School of Medicine, Aristotle University of Thessaioniki, Thessaioniki, Greece

Edited by:

Karin Trimmei, Medical University of Vienna, Austria Reviewed by:

David R. Van Wagoner, Cleveland Clinic Lerner College of Medicine of Case Western Reserve University, USA

George E. Billman, The Ohio State University, USA


Vasilios Papaioannou, Intensive Care Unit, Alexandroupolis General Hospital, Democritus University of Thrace, Dragana 68100, Alexandroupolis, Greece e-mail:; papabil69@gmail. com

Many experimental and clinical studies have confirmed a continuous cross-talk between both sympathetic and parasympathetic branches of autonomic nervous system and inflammatory response, in different clinical scenarios. In cardiovascular diseases, inflammation has been proven to play a pivotal role in disease progression, pathogenesis and resolution. A few clinical studies have assessed the possible inter-relation between neuro-autonomic output, estimated with heart rate variability analysis, which is the variability of R-R in the electrocardiogram, and different inflammatory biomarkers, in patients suffering from stable or unstable coronary artery disease (CAD) and heart failure. Moreover, different indices derived from heart rate signals' processing, have been proven to correlate strongly with severity of heart disease and predict final outcome. In this review article we will summarize major findings from different investigators, evaluating neuro-immunological interactions through heart rate variability analysis, in different groups of cardiovascular patients. We suggest that markers originating from variability analysis of heart rate signals seem to be related to inflammatory biomarkers. However, a lot of open questions remain to be addressed, regarding the existence of a true association between heart rate variability and autonomic nervous system output or its adoption for risk stratification and therapeutic monitoring at the bedside. Finally, potential therapeutic implications will be discussed, leading to autonomic balance restoration in relation with inflammatory control.

Keywords: heart rate variability, inflammation, autonomic nervous system, coronary artery disease, cardiovascular disease, mortality


Systemic inflammation is a normal response to altered homeostasis and has an important role in several pathophysiological processes, such as infection or trauma. It is characterized by the endocrine release of different cytokines, such as tumor necrosis factor a (TNF-a) and interleukin-1 (IL-1), IL-4, IL-6,

Abbreviations: ACE, angiotensin-converting enzyme; ACTH, adrenocorti-cotropin hormone; ANS, autonomic nervous system; ATRAMI, Autonomic Tone and Reflexes After Myo-cardial Infarction Study; CAD, coronary artery disease; CHF, congestive heart failure; CNS, central nervous system; CRP, C-reactive protein; DHA, docasah-exaenoic acid; DMN, dorsal motor nucleus of the vagus; DVC, dorsal vagal complex; ECG, electrocardiogram; EPA, eicosapentaenoic acid; FFT, Fast Fourier trans-formation; HPA, hypothalamo-pituitary-adrenal axis; HRV, heart rate variability; HF, high frequency; LF, low frequency; LVEF, left ventricular ejection fraction; MODS, multiple organ dysfunction syndrome; NYHA, New York Heart Association; NTS, nucleus tractus solitaries; pNN50, proportion derived from dividing NN50 (number of interval differences of successive intervals greater than 50 ms) by the total NN intervals; PSD, power spectrum density; PUFA, polyunsaturated fatty acids; PVN, paraventricular nucleus; RMSSD, square root of the mean squared differences of successive intervals; RVLM, rostral ventrolateral medulla; TNF-a, tumor necrosis factor a; SAN, sinus atrial node; SDNN, standard deviation of the normal-to-normal intervals which is the square root of the variance; SD, standard deviation; SNS, sympathetic nervous system; THS, Twins Heart Study; ULF, ultra low frequency; VLF, very low frequency; WBC, white blood cell count.

IL-10, and many others, normally confined to paracrine regulation of a local inflammatory response (Koj, 1997; Sporn, 1997). Apart from their involvement in local and systemic inflammation, cytokines may induce activation of brain-derived neuroendocrine immunomodulatory responses. Neuro-endocrine pathways, such as hypothalamo-pituitary-adrenal (HPA) axis and both the sympathetic and parasympathetic divisions of the autonomic nervous system (ANS) are powerful modulators of inflammation, typically through an anti-inflammatory balancing mechanism (Reichilin, 1993; Webster et al., 2002).

Recently, it has been demonstrated that subclinical inflammation and the con-centration of inflammatory markers, such as cytokines, correlate strongly to cardio-vascular mortality and morbidity in both healthy subjects and in those with known coronary artery disease (CAD) (Phillips et al., 1992; Ridker et al., 1997). Furthermore, vascular inflammation plays a critical role in the initiation, evolution, and rupture of atherosclerotic plaque (Ross, 1993).

In the healthy state there is some degree of stochastic variability in physiologic variables, such as heart rate (heart rate variability). This variability is a measure of complexity that accompanies healthy systems and has been suggested to be responsible for their greater adaptability and functionality related

to pathologic systems (Buchman, 2002). Studying physiological signals of patients can easily identify "hidden" information concerning inherent dynamics and overall variability within a time series. Recognition that physiologic time series contain such information defies traditional mechanistic approaches based on conventional biostatistical methodologies and has fueled growing interest in applying techniques from statistical physics for the study of living organisms (Seely and Christou, 2000). Through those techniques different "physiomarkers" can be estimated that fulfill the requirements of contemporary medicine for better and more accurate early warning signs, since they are based on high-frequency measurements and are much easier to measure at the bedside (Seely and Christou, 2000). In this respect, a number of international databases and different processing methods of heart rate signals have been developed with free access from different investigators, such as the Web Site Physionet (

On the contrary, it has been repeatedly demonstrated that various "biomarkers" such as cytokines, exhibit marked interdependence, pleiotropy (multiple effects) and redundancy (multiple cytokines with the same effect) (Friedland et al., 1996). At the same time, their plasma concentrations fluctuate from day to day and correlate poorly with classic physiologic variables in different groups of patients (Friedland et al., 1996; Seely and Christou, 2000; Buchman, 2002). Furthermore, biomarkers are difficult to obtain routinely at the bedside. In addition, the financial cost of various immunoassay techniques for their detection in blood samples tends to become an inhibiting factor for their extensive use, as diagnostic or even prognostic tools, in many Medical Centers.


The pathophysiological link between the communication between the Central Nervous System (CNS) and the immune-regulated inflammation is the capability of the brain to monitor and to affect at the same time the immune status. The first mechanism relies upon activation of vagus nerve afferent fibers that signal the brain that inflammation is occurring. Different kind of mediators such as cytokines can activate visceral vagus afferent fibers which terminate within the dorsal vagal complex (DVC) of the medulla oblongata. The DVC consists of the nucleus tractus solitarius (NTS), the dorsal motor nucleus of the vagus (DMN) and the area postrema (AP) (Berthhoud and Neuhuber, 2000). Ascending projections from the NTS reach hypothalamic paraventricular nucleus (PVN), which is associated with the synthesis and release of corticotropin releasing hormone (CRH). This factor induces the production of adrenocorticotropin hormone (ACTH) from the anterior pituitary, which is the main inducer of the synthesis of immuno-suppressive glucocorticoids from the adrenal cortex. Projections from NTS are connected to the DMN and to rostral ventrolateral medulla (RVLM). This region increases firing of the noradrenergic preganglionic neurons in the spinal cord (Tracey, 2007).

The brain can affect the immunological status through the activation of the HPA axis and increased outflow of sympathetic (SNS) and parasympathetic nervous system. The SNS

activation during the early stages of stress induces local inflammatory response through a2-subtype adrenoreceptor stimulation by norepinephrine (NE), whereas stimulation of ß2-subtype adrenoreceptor-cAMP-protein kinase A pathway is associated with an inhibition of pro-inflammatory cytokines' production (van der Poll et al., 1996; Elenkov et al., 2000; Zhou et al., 2001). It seems that SNS activation protects the organism from the detrimental effects of pro-inflammatory cytokines, while it can increase local inflammatory response (Chrousos, 1995). In addition to the SNS, a link between the parasympathetic part of the ANS and immune-regulatory processes has been suggested (Tracey, 2002). It has been demonstrated that acetylcholine decreases TNF-a production by endotoxin-stimulated human macrophage cultures, through a7-subunit of the nicotinic acetylcholine receptor (Wang et al., 2003; de Jonge et al.,

2005). The vagus nerve cholinergic signaling interacts with the above receptor on immune cells in the spleen and inhibits TNF-a production and release into the circulation (Huston et al.,

2006). Acetylcholine is also effective in suppressing other proinflammatory cytokines such as IL-1ß, IL-6, and high mobility group box 1 (HMGB1) protein (Wang et al., 2004). This "cholinergic anti-inflammatory pathway" is responsible for a "hardwired" connection between the nervous and immune systems and is considered, as the primary component of the "immuno-reflex." A more complete understanding of these reflexes can yield insight into both pathophysiological pathways and therapeutic strategies in many pathological processes, including infections, sepsis and cardiovascular diseases (Tracey, 2002, 2007).

In conclusion, there is strong evidence that CNS controls body's systemic response to inflammation. Recently, different clinical studies investigating a possible association between ANS outflow and various inflammatory indices in patients with heart diseases have appeared in the literature (Aronson et al., 2001; Malave et al., 2003; Janszky et al., 2004; Shehab et al., 2004; Hamaad et al., 2005; Lanza et al., 2006; Madsen et al., 2007; Nolan et al., 2007; Psychari et al., 2007; von Känel et al., 2011). The aim of these studies was to measure ANS activity through a set of different "physio-markers" and correlate them with various biomarkers that can indirectly assess inflammatory response in different clinical scenarios, such as CAD (Janszky et al., 2004; Hamaad et al., 2005; Lanza et al., 2006; Madsen et al., 2007; Nolan et al., 2007; Psychari et al., 2007; von Känel et al., 2011) and heart failure (Aronson et al., 2001; Malave et al., 2003; Shehab et al., 2004). Moreover, the prognostic value of such measurements was tested in different groups of patients with cardiovascular diseases, in terms of mortality and risk of rehospitalization.

The best "physiomarkers" are obtained from analysis of heart rate variability (HRV); that is, the variability of R-R series in the electrocardiogram (ECG), and its frequency components (Lombardi et al., 1987; Task Force, 1996). Beat-to-beat fluctuations reflect the dynamic response of the cardiovascular control systems to a host of naturally occurring physiological perturbations. A variety of animal and human research has established two clear frequency bands in heart rate signals. These bands include high frequency oscillations, between 0.15 and 0.4 Hz that are associated with respiration, and bands with a lower frequency range, below 0.15 Hz (Task Force, 1996). Akselroad et al. (1981)

introduced power spectrum analysis of heart rate fluctuations in order to quantify beat-to-beat cardiovascular control. Power spectrum density (PSD) analysis provides the basic information of how power (variance) distributes as a function of frequency (Akselroad et al., 1981; Malik and Camm, 1993). In 1996, the Task Force of the European Society of Cardiology and the Northern American Society of Pacing and Electrophysiology published guidelines regarding standardization of nomenclature, specification of methods of measurement, definition of physiological and pathophysiological correlates, description of clinical applications and identification of different areas for future research.

The association of higher risk of post-infarction mortality with reduced HRV was first shown by Wolf et al. (1978). The clinical importance of HRV became appreciated in the late 1980s, when it was demonstrated that low HRV was a strong and independent predictor of mortality after an acute myocardial infarction (MI) (Kleiger et al., 1987; Lombardi et al., 1987).


The RR variations may be evaluated by a number of methods: TIME DOMAIN METHODS

Time domain methods determine heart rate or RR intervals in continuous ECG records. Each QRS complex is detected and the normal-to-normal (NN) intervals (all intervals between adjacent QRS complexes) are calculated. Other time domain variables include the mean NN interval, the mean heart rate or the difference between the longest and the shortest NN interval, as well. The simplest of these metrics is the standard deviation of the NN intervals (SDNN), which is the square root of the variance. However, it should be emphasized that SDNN becomes less accurate with shorter monitoring periods. The most commonly used time domain methods are the square root of the mean squared differences of successive NN intervals (RMSSD), the number of interval differences of successive NN intervals greater than 50 ms (NN50) and the proportion derived from dividing NN50 by the total NN intervals (pNN50) (Akselroad et al., 1981; Task Force, 1996; Table 1).


Spectral analysis of heart rate partitions HRV into its frequency components. Most commonly used methods are Fast Fourier Transformation (FFT) and auto-regressive modeling. FFT displays in a plot the relative contribution (amplitude) of each frequency. This plot includes at least three peaks.Fast periodicities in

Table 1 | HRV metrics in time domain.

SDNN Standard deviation of all N-N intervals

SDNN index Average of the standard deviations of N-N intervals for each 5-min period

SDANN Standard deviation of the average N-N intervals for each

5-min period over 24 h NN50 Number of N-N intervals differing by >50 ms from the

preceding interval pNN50 Percentage of adjacent cycles that are >50 ms apart

RMSSD Root mean square of successive differences in ms

the range 0.15-0.4 Hz [high frequency (HF)] are largely due to the influence of the respiratory phase on vagal tone. Low-frequency periodicities (LF), in the region of 0.04-0.15 Hz, are produced by baroreflex feedback loops, affected by both sympathetic and parasympathetic modulation of the heart, whereas very low frequency periodicities (VLF), in the frequency range between 0.003 and 0.04Hz and ultra low frequencies (ULF, <0.003 Hz) have been variously ascribed to modulation by chemoreception, thermoregulation and the influence ofvasomotor activity. The area under the power spectral curve (power) in a particular frequency band is considered to be a measure of HRV at that frequency, whereas the LF/HF ratio has been suggested as an indirect index of sympatho-vagal balance (Task Force, 1996). According to the report of the Task Force, the analyzed ECG signals must satisfy several technical requirements in order to obtain reliable information. The optimal sampling frequency range should be between 250 to 500 Hz. Ectopic beats, arrhythmic events, missing data and noise effects should be properly filtered and omitted. Frequency domain methods must be preferred in cases of short term investigations. The recordings should last for at least 10 times the wavelength of the lower frequency bound, thus recordings of ~1 min can assess the HF component of HRV while 2 min are needed for the LF component. In conclusion, 5-min recordings are preferred, unless the aim of the study dictates a different design (Akselroad et al., 1981; Lombardi et al., 1987; Task Force, 1996; Table 2).



The cyclic variations in intrathoracic pressure perturbate venous return, cardiac out-put and thus, blood pressure. These changes are sensed by baroreceptors and result in changes in autonomic activity to the heart. These perturbations are mediated via the vagus nerve as atropine administration abolishes high frequency oscillations in heart rate (Akselroad et al., 1981). It seems that a major cause of respiratory sinus arrhythmia (RSA) is a central coupling of respiratory drive to cardiac vagal motor neurons. However, the changes in vagal activity are partly induced by baroreceptor sensing of respiratory oscillations in blood pressure and reflect all components of the baroreflex loop (DeBoer et al., 1987). In addition, factors such as reduced respiratory capacity and body position may alter the amplitude of high frequency oscillations in blood pressure and subsequently the HF component of heart rate signals (Malpas, 2002). Thus, heart rate variability analysis cannot be used for comparisons between different patient groups as there is a need for controlling ventilation for both rate and depth. Moreover, and since there is marked inter-individual variation in the relationship between HRV and parasympathetic effect, differences in HRV between individuals may reflect differences in this relationship, as was postulated by Goldberger (Goldberger et al., 2001). In this respect,

Table 2 | HRV metrics in frequency domain.

ULF (ultra low frequency) ms2 24h recordings <0.003Hz VLF (very low frequency) ms2 24 h and 5-min recordings -0.003-0.04 Hz LF (low frequency) ms2 24h and 5-min recordings -0.04-0.15Hz

HF (high frequency) ms2 24h and 5-min recordings -0.15-0.4Hz

this relationship in humans was described by a quadratic function in which there is an initial ascending limb, where HRV increases in parallel with vagal effect until it reaches a plateau level. Beyond this level, HRV decreases with further augmentation of vagal tone, probably due to a saturated HRV response upon intense autonomic stimulation (Goldberger et al., 2001). Finally, age and sex-related differences have also been associated with this variability (Goldberger et al., 2001).


The LF component of HRV is probably the most contentious aspect with respect to cardiovascular variability. There are two opposing theories in the literature proposing different potential origins: (1) the central oscillator theory (Montano et al., 1996) and (2) the baroreflex feedback loop theory (Lanfranchi and Somers, 2002). According to the first theory, it is believed that LF oscillations reflect sympathetic tone and are generated by the brain stem circuits. In cats, Montano (Montano et al., 1996) analyzed the dis-charges of single sympathetic neurons located in the rostral ventrolateral medulla and caudal ventrolateral medulla. He observed activity at 0.12 Hz, which was positively correlated with heart rate and blood pressure variability. As the above oscillations remained after sino-aortic and vagal resection, it was assumed that the central nervous system is able to generate such oscillations.

According to the baroreflex feedback loop theory, a change in blood pressure is sensed by arterial baroreceptors, resulting in heart rate adjustment through the central nervous system and via both the fast vagal and the slower sympathetic actions (Lanfranchi and Somers, 2002). At the same time, baroreceptors induce a slow sympathetic withdrawal from the vessels. The delay in the sympathetic branch of the baroreflex in turn determines a new oscillation, which is sensed by the baroreflex and induces a new oscillation in heart rate. It has been also proposed that the LF oscillation arises from the interaction of slow sympathetic and fast vagal responses, where baroreflex buffering of the slow respiratory induced blood pressure oscillations results in resonant low frequency oscillations, due to the delay in the slow conducting sympathetic loop of the baroreflex (DeBoer et al., 1987).

In conclusion, it must be stressed that the low frequency oscillations of heart rate reflect the ability of the individual components of the baroreflex feedback loop to respond to different inputs that can alter the power of such oscillations and they are not just a measure of sympathetic nerve activity.


The reasons for reduced HRV during cardiovascular diseases have been debated and two theories have been developed. The first theory focuses on reduction of vagal tone and has been introduced by Akselroad et al. (1981). The second theory developed by Goldberger and colleagues (2002) states that normal physiology has fractal-like properties with high levels of complexity that explain phenomena such as HRV. Its reduction during severe disease reflects a "de-complexification," mostly attributed to uncoupling between different restorative mechanisms (Godin and Buchman, 1996). In addition, accumulating evidence from both in vitro and ex vivo experiments support a potential third

mechanism (Griffin et al., 2005), which is associated with an intracardiac origin of HRV. According to this hypothesis, sinus atrial node (SAN) cells can be viewed as an amplifier of various input signals (Zaza and Lombardi, 2001). During cardiovascular diseases, an unfavorable metabolic milieu could affect ion channel gating properties or membrane receptor densities, with significant impact upon level and variability of pacemaker activity. In addition, a possible reduced responsiveness of SAN cells to external stimuli could also negatively affect HRV (Zaza and Lombardi, 2001).

Moreover, different clinical studies in heart transplant recipients have found evidence for heart rate fluctuations originating from the heart itself (Hrusheskyet al., 1984; Bernardi etal., 1990). Bernardi studied intrinsic mechanism regulating HRV in both transplanted and intact heart during exercise (Bernardi et al., 1990). He found that at peak exercise a non-autonomic mechanism, probably intrinsic to the heart muscle, may determine heart rate fluctuations in synchrony with ventilation, in transplanted as well as in intact hearts. Hrushesky and colleagues (1984) quantified respiratory sinus arrhythmia and found that individuals with a transplanted heart had resting RSA values similar to healthy subjects.

In conclusion, there is marked inter-individual variation between HRV response and different levels of autonomic stimulation. Basal autonomic activity, age and sex differences, alterations in expression of ion channel activity or autonomic receptors could be responsible for individualized curves, relating auto-nomic effects to HRV (Eckberg, 1997; Goldberger et al., 2001). In addition, LF/HF ratio has been criticized as an indirect measure of sympathovagal balance, reflecting rather autonomic fluctuations and not absolute measures of autonomic nerve traffic (Eckberg, 1997; Billman, 2013). Thus, interpretation of different studies investigating HRV alterations in different groups of patients should be cautious since variability in time of recordings and methods for HRV analysis, as well as heterogeneity of studying population, limit generalization of their findings.


The first large prospective population study that reported the significant prognostic value of low HRV after an acute myocardial infarction was the Autonomic Tone and Reflexes After Myocardial Infarction Study (ATRAMI) (La Rovera et al., 1998), and included 1284 patients with a recent (<28 days) myocardial infarction. A 24 h Holter recording was done to quantify HRV (using SDNN values) and ventricular arrhythmias. Low values of HRV (SDNN < 70 ms) carried a significant multivariate risk of cardiac mortality. Furthermore, the association of low SDNN with left ventricular ejection fraction (LVEF) <35% carried a relative risk of 6.7, compared with patients with LVEF above 35%. Investigators from the Framingham Heart Study (Tsuji et al., 1994) computed HRV time and frequency domain measures in 736 patients and correlated them with all-cause mortality during 4 years of follow-up. They concluded that HRV offers prognostic information independent of that provided by traditional risk factors.

During the Zutphen study (Dekker et al., 1997), 885 middle-aged (40-60 years old) and elderly Dutch men (aged 65-85)

were followed from 1960 until 1990, whereas SDNN was determined from the resting 12-lead ECG. It was shown that low HRV is predictive of mortality from all causes, indicating that it can be used as an index of compromised health in the general population. It seems that the predictive value of low HRV is independent of other factors, such as depressed left ventricular ejection fraction and presence of late potentials (Kleiger et al., 1987; Lombardi et al., 1987). In addition, retrospective ECG data analysis from 127 patients included in the Veterans Affairs' Survival Trial of Antiarrhythmic Therapy in Congestive Heart Failure (CHF) (Bilchicket al., 2002) demonstrated that CHF patients with SDNN <65.3 ms had a significantly increased risk of sudden death. Moreover, this study demonstrated that every 10 ms increase in SDNN conferred a 20% decrease in risk of mortality.


The relationship between HRV and inflammation has been studied mainly in patients with acute or stable CAD, CHF and metabolic syndrome with impaired glucose tolerance (Brunner et al., 2002). The inflammatory biomarkers that were used included C-reactive protein (CRP), TNF-a, IL-6 and white blood cell count (WBC).


Hamaad et al. (2005) tested the association between time and frequency domain indices of HRV and circulating IL-6, high sensitivity CRP (hs-CRP) and white cell counts, in a sample of 100 patients with proven acute coronary syndrome. In addition, they compared these metrics with healthy controls (n = 49) and estimated possible relationships on repeated measures at 4 months in recovery (n = 51). They found modest negative correlations between all inflammatory biomarkers and mainly SDNN, VLF and LF power. The strongest associations were seen between WBC and SDNN (r = -0.351). However, relationships did not persist on multivariate analyses after a 4-month period. According to the authors, the correlations were observed largely among HRV indices reflecting sympathetic activity, suggesting that the inflammatory response in acute coronary events may be associated with sympathetic activation instead of vagal withdrawal. Furthermore, leukocytosis observed in these patients seems to be a potential source of pro-inflammatory cytokines within the atheromatous plaque and might induce a potential rupture.

In another study, Lanza and colleagues (2006) assessed HRV and measured CRP serum levels within 24 h of admission in 531 patients with unstable angina pectoris. They found a significant negative correlation between CRP levels and all HRV metrics derived from both time and frequency domain, with the highest correlation coefficient with SDNN and VLF. After categorizing patients into 4 subgroups according to CRP quartile levels, significantly lower HRV values were found in the upper CRP quartile. The subsequent multivariate analysis revealed that SDNN and VLF were the most significant predictors of increasing CRP, whereas CRP was a strong predictor of impaired ANS activity as well.

In a study including patients with suspected CAD, Madsen et al. (2007) enrolled 269 subjects referred for elective coronary

angiography. They found that SDNN of heart rate signals was significantly higher in the lower CRP quartile compared to the upper one, whereas associations were stronger for patients with a previous myocardial infarction and with significant coronary stenoses.

In a similar study (Nolan et al., 2007), a negative correlation between CRP and HRV frequency components was reported, whereas a decreased HF power (reflecting vagal tone) in the high CRP quartile, compared to the lowest one, was found. In the study by Janszky and colleagues (2004) that included only female patients who survived hospitalization for acute myocardial infarction and were evaluated 1 year after the event, levels of IL-6 showed an inverse relation with all HRV frequency measures, except for HF. However, this relationship, as well as the association between CRP levels and IL-1 receptor antagonist (IL-1ra) with HRV indices was non-significant. Psychari et al. (2007) also reported a strong inverse association between CRP and several HRV indices (SDNN, HF, and LF) in post-MI patients and after adjustment for left ventricular function.

Recently, von Kanel et al. (2011) investigated the association between HRV measured in the time domain, CRP, IL-6 and fib-rinogen, in a cohort of 862 subjects recruited from the Heart and Soul Study, which assessed health outcomes in 1.024 outpatients with stable CAD. They found that SDNN was inversely and significantly associated with inflammatory indices, after adjustment of all covariates.


In 2001, Aronson et al. (2001) evaluated for the first time the relationship between HRV metrics derived from both time and frequency domains and different biomarkers, such as IL-6, TNF-a, and serum levels of norepinephrine, in 64 patients admitted for decompensated chronic heart failure. TNF-a levels did not correlate with any of the HRV indices. However, IL-6 was inversely correlated with SDNN (r = -0.36), with total power of heart rate signals and ULF (r = -0.37 and r = -0.43, respectively). No correlation was found between IL-6 and time (pNN50 and RMSSD) or frequency domain (HF power) indices of vagal activity.

Malave et al. (2003) examined HRV in relation to circulating levels of TNF-a, TNF-a receptors and norepinephrine in 10 controls, 15 patients with mild CHF and 14 subjects with moderate heart failure. There was a significant inverse linear correlation between increased levels of all biomarkers and SDNN, LF and HF power among CHF patients. In addition, LF power was more closely correlated with circulating levels of TNF-a than was the HF component, whereas multiple linear regression analysis showed that TNF-a was a stronger predictor of reduced HRV than was the circulating levels of norepinephrine. The authors concluded that over-expression of TNF-a and subsequent loss of ^-adrenergic responsiveness contributes to the decrease in HRV, observed in heart failure. According to findings from experimental studies (Chung et al., 1990), TNF-a might inhibit ^-adrenergic signal transduction through either activation of Gi proteins or impairment of activation of Gs proteins, something that could be viewed as an adaptive mechanism in the early stages of CHF, protecting cardiac myocytes from the deleterious

actions of catecholamines. However, in the more advanced stages of the disease, this mechanism could become maladaptive, leading to a reduction in cardiac output (Mann et al., 1992).

Finally, in a small prospective study that included 34 patients with CHF followed for a 2-year period with monthly CRP measurements and 24-h Holter recordings, it was shown that five unexpected deaths that occurred were preceded by progressive increases in both CRP serum levels and autonomic dysfunction (low HRV indices) (Shehab et al., 2004).


As a part of the Copenhagen Holter study, that assessed the value of 24-h Holter recording in the risk assessment of men and women aged 55, 60, 65, 70, and 75 years with no apparent heart disease, Sajadieh et al. (2004) investigated the associations between time domain components of HRV, CRP and WBC in 643 healthy men and women. They found that SDNN was negatively correlated with smoking, inflammatory indices, blood sugar, triglyceride concentration, female gender and diabetes. Moreover, in multivariate regression analysis, increased heart rate and reduced HRV were significantly related to white blood cell count and CRP. The reduction of SDNN was attributed to sympathetic predominance, whereas lack of any association between inflammation and pNN50, which is considered as a marker of vagal activity, indicates that reduced HRV is mainly due to increased sympathetic activity rather than vagal withdrawn.

From the Whitehall II cohort, a multicenter epidemiologic investigation of over 5000 subjects, two studies (Owen and Steptoe, 2003; Sloan et al., 2007) used sub-samples to examine the relation between HRV indices derived from the frequency domain and inflammation, in healthy subjects. In the first study, Sloan found in a sample of 757 people, an inverse correlation between CRP and IL-6 with both LH and HF components of HRV power spectrum. In the second study, Owen and Steptoe did not find any association between IL-6, TNF-a and time domain measures of HRV, in a group of 211 healthy adults.

Other investigators (Albert et al., 2002) reported a strong positive association between CRP levels and the long-term risk of sudden cardiac death, in case-control analysis among healthy individuals, followed for 17 years in the Physician's Health Study. Men in the upper CRP quartile had a 2.8 fold increased risk of sudden cardiac death compared to men in the lower quartile. According to the authors, a low-grade inflammation involved in atherosclerosis shifts ANS balance toward sympathetic activation, making individuals more prone to ventricular arrhythmias and sudden cardiac death.

The Twins Heart Study (THS) (Goldberg et al., 2002) was an investigation of psychological, biological and behavioral risk factors for subclinical cardiovascular diseases in 7.369 middle-aged male-male twin pairs, who served in the United States military during the Vietnam War. From this registry, a cohort of 264 twins free of symptomatic CAD was examined by Lampert et al. (2008), for assessing possible associations between HRV, CRP, and IL-6. They found an inverse relationship between frequency domain HRV metrics (except for HF) and both CRP and IL-6. These associations persisted after adjustment for other traditional CAD risk

factors, such as smoking, hypertension, diabetes, high-density lipoprotein (HDL) and depression.

A significant confounding factor of HRV analysis that has to be considered in these studies includes the presence of depressive symptoms and anxiety. It has been estimated that ~ 12-20% of hospitalized cardiac patients suffer from major depression, whereas 15% of subjects following acute myocardial infarction exhibit a posttraumatic stress disorder (Frasure and Lesperance, 2006; Garder and von Karel, 2006; Pizzi et al., 2008). In a 2-year follow up prospective observational study, Pizzi investigated the relation between time domain HRV indices, IL-6, TNF-a, CRP and depression in a cohort of 415 subjects free of CADs, with at least two CAD risk factors (age, male gender, current smoking, hypertension, dislipidaemia). All HRV and inflammatory indices were significantly associated with depression. Logistic regression further showed that depressive individuals were more likely to have a higher CRP and IL-6 and altered HRV (lower SDNN).

In a recent study of Kop and colleagues (2010) who recruited 908 patients, free of CAD, for a median follow-up period of 13.3 years, it was demonstrated that among depressed participants, HF power of HRV was negatively correlated with CRP (r = -0.205), IL-6 (r = -0.233) and WBC (r = -0.292). Moreover, depression was associated with high IL-6 serum levels and increased cardiovascular mortality risk. In conclusion, in patients without heart disease depression seems to be associated with HRV imbalance and inflammation.

Table 3 summarizes the majority of clinical studies that have evaluated the relationship between different inflammatory biomarkers and HRV metrics, in patients with different cardiovascular diseases.


The presence of sympathetic overactivity, autonomic dysfunction, inappropriately increased heart rate, insulin resistance and in some cases, cardiomyopathy with reduced cardiac contractility, has also been observed during severe sepsis and multiple organ failure (Muller-Werdan et al., 2006). However, in cardiac patients sympathetic activity dominates over vagal tone where in septic patients both branches of ANS are attenuated (Muller-Werdan et al., 2006). For these reasons, it has been hypothesized that during critical illness except for ANS impairment, a defective signal transduction at the level of pacemaker cells could also account for observed differences between cardiac and septic patients (Fairchild et al., 2009).

Alterations in HRV during septic shock and multiple organ dysfunction syndrome (MODS), have been reported from different research groups (Goldstein and Buchman, 1998; Goldstein et al., 1998; Seely and Christou, 2000). In this respect, Goldstein et al. (1998) found that both increased total variability and LF power were associated with recovery and survival, whereas a decrease in total power, LF/HF and LF power correlated with severity of illness and mortality in septic patients, 48 h after being admitted to the Intensive Care Unit. In an animal study of experimental endotoxemia, induced by administration of lipopolysacchraride (LPS, endotoxin derived from the cell wall of Gram-negative bacteria) Fairchild and colleagues demonstrated a

Table 3 | Summary of several clinical studies investigating a possible association between HRV indices and inflammation in patients with CAD, CHF and healthy individuals.

References Study population Duration and HRV measures Inflammatory indices Results

Hamaad etal., 2005 100 patients with acute CAD vs. 29 healthy controls 20 min time, time and frequency domain CRP IL-6 Negative correlation with SDNN, VLF and LF

Lanza etal., 2006 531 patients with unstable CAD 24-h time, time and frequency domain CRP Inverse correlation between CRP with SDNN and VLF

Madsen etal., 2007 269 patients with suspected CAD 24-h time, time domain CRP Upper CRP quartile negatively correlated with SDNN

Nolan etal., 2007 29 patients with CAD 5-min time, frequency domain CRP HF power decreased in high CRP group

Psychari etal., 2007 98 patients with acute CAD (post-MI) 24-h time, time and frequency domain CRP Inverse relation between CRP and SDNN, HF and LF power

von Kanel etal., 2011 862 patients with CAD 24-h time, time domain CRP IL-6, fibrinogen Inverse association between CRP IL-6 and SDNN

Aronson etal., 2001 64 patients with CHF 24-h time, frequency domain TNF-a, IL-6 IL-6 inversely correlated with SDNN and ULF power

Malave etal., 2003 10 healthy controls, 15 patients with mild CHF, 14 patients with moderate CHF 24-h time, frequency domain TNF-a, TNF soluble type 1 and 2 receptors Inverse correlation between inflammatory measures, SDNN, LF and HF

Sajadieh etal., 2004 643 subjects without CHF 24-h time, time domain CRP WBC Inverse correlation between SDNN with CRP and WBC SDNN predictor of CRP

Sloan etal., 2007 757 young healthy adults 10-min time, frequency domain CRP IL-6 CRP and IL-6 inversely correlated with HF and LF

Owen and Steptoe, 2003 211 healthy subjects 20-30 min time, time domain TNF-a, IL-6 No relation between both TNF-a and IL-6 with HRV

Lampert etal., 2008 264 healthy twins individuals 24-h time, frequency domain CRP IL-6 Inverse relation between CRP and IL-6 with all HRV frequency metrics (except for HF)

Abbreviations: CRP, C-reactive protein; CAD, coronary artery disease; CHF, congestive heart failure; MI, myocardial infarction; WBC, white blood cell count.

strong inverse correlation between SDNN and total power of RR time series and peak concentrations of different cytokines, 3-9 h post-LPS (Fairchild et al., 2009). The same results were found after administration of recombinant TNF-a. It was suggested that mechanisms responsible for decrease in HRV could be related with effects of LPS and/or cytokines on various ion channels.

Tateishi et al. (2007) investigated the relationships between HRV and interleukin 6 upon admission in a cohort of 45 septic patients and they found that IL-6 exhibited significant negative correlations with both LF and HF power values. These findings indicate a possible association between low HRV indices and hyper-cytokinemia.

In another study Papaioannou et al. (2009), we investigated possible associations between different HRV indices and various biomarkers of inflammation, in 45 septic patients and during the first 6 days of their stay in the ICU. We daily assessed HRV in

time (SDNN) and frequency domain (LF, HF and LF/HF) and measured C-reactive protein, Interleukin 6 and 10 serum levels in two groups of patients. The first group included subjects suffering from sepsis with mean Sequential Organ Failure Assessment score of severity of illness (SOFA) < 10 (n = 25) and the second group included patients with septic shock (SOFA > 10, n = 20). This study found in the group of patients with SOFA > 10, statistically significant inverse correlations between CRP and LF/HF ratio (r = -0.61), (Figure 1A) and positive correlations with HF (r = 0.80). At the same time, IL-10 proved to be significantly correlated with HF and SOFA score in a positive way and with LF, LF/HF and SDNN in a negative way. Finally, the total variability of heart rate signals (SDNN) was found to be negatively correlated with both CRP (r = -0.79) and SOFA score (r = -0.84) (Figure 1B). IL-6 was not significantly correlated with any HRV parameter. It is possible that, the pleuripotency of this cytokine

A 60 -, j 3,5

1 2 3 4 5 6 days

1 2 3 4 5 6 days

FIGURE 1 | (A) Longitudinal trends over time of mean values of CRP and LF/HF ratio, reflecting sympathovagal balance, for patients with SOFA > 10, during the 6 days of study period. [log transformed data, adapted from Papaioannou et al. (2009)]. It appears that LF/HF changes inversely with CRP (B) Longitudinal trends over time of mean values of CRP and SDNN (secs), for patients with SOFA > 10, during the 6 days of study period. [log transformed data, adapted from Papaioannou et al. (2009)]. There is a progressive increase in SOFA score from day 1 until day 4 (development of septic shock) and a subsequent downward shift in its values. At the same time, the variability of heart rate signals estimated with SDNN seems to be significantly reduced during the development of septic shock.

could be responsible for our results since IL-6 can behave as both a pro-inflammatory activator (induces the production of CRP) and inhibitor (limits Tumor necrosis factor a and IL-1p secretion), at the same time (Chrousos, 1995).

These findings suggest that reduction in HRV and LF/HF is related with an augmented pro- and anti-inflammatory response during sepsis, especially in more severely ill patients. Furthermore, severity of illness is positively associated with HF and IL-10 serum concentrations and changes inversely with variability of heart rate signals. In this respect, elevated levels of IL-10 have also been found in trauma patients who developed sepsis and multiorgan failure (Sherry et al., 1996).

In conclusion, it seems that critical illness and high cytokine levels are associated with reduced HRV, however, existing literature does not elucidate whether loss of HRV is related to an endotoxin effect at the level of ANS output, baroreflex sensitivity or the pacemaker cell itself. However, results from a prospective study in a group of 40 healthy adults who received a single intravenous bolus of 2 ng/kg LPS, suggested that there is no relationship between basal cardiac ANS activity, and the inflammatory response (Kox et al., 2011). Thus, no association was found between frequency components of HRV that were determined hourly and until 8 h after LPS administration and different pro- and anti-inflammatory cytokines, measured at various time

points. According to the authors, vagus nerve innervation of the heart does not reflect outflow to other organs, such as the spleen, one of the major cytokine-producing organs (Tracey, 2002, 2007). As basal vagal input to the spleen may be different from vagal input to the heart, HRV might not be an appropriate method to assess activation of the cholinergic anti-inflammatory pathway (Kox et al., 2011). These findings strengthen the notion that autonomic outflow cannot be regarded as a general response, but appears to be organ-specific. In this respect, results from different studies discussed so far lack generalization and robustness due to different study populations and design, interspecies differences and potential impact of severity of disease, sedation or mechanical ventilation upon HRV (Goldstein and Buchman, 1998). However, it appears that an inverse association between inflammation and total variability of heart rate signals could be found in the most severe cases.


Different clinical trials have shown that fatty acids from fish oil can be considered as powerful disease-modifying nutrients in patients with acute lung injury, sepsis and cardiovascular diseases (Christensen et al., 1999; Abuissa et al., 2005; Pontes-Arruda et al., 2006; Singer et al., 2006). Particularly, feeding with the very-long chain, w-3 polyunsaturated fatty acids (PUFAs) eicosapentaenoic acid (EPA) and docasahexaenoic acid (DHA) has been found to inhibit the activity of the pro-inflammatory transcription factor nuclear factor kB (NF-kB) and subsequently, to attenuate the production of different cytokines, chemokines and other effectors of innate immune response (Singer et al., 2008). In the cardiovascular literature, it has been shown that oral supplementation of w-3 PUFAs increase instantaneous HRV, reduce LF/HF ratio and confer protection against ischemia-induced ventricular tachycardia and sudden cardiac death (Abuissa et al., 2005). Moreover, Christensen (Christensen et al., 1999) demonstrated that fish oil feeding can induce an incorporation of DHA into the membranes of granulocytes, which is associated with a dose-response increase in HRV (SDNN), something that may protect against serious ventricular arrhythmias. Such effects of fish oil reflect an enhanced efferent vagal activity via a central-acting mechanism, due to a possible suppression of pro-inflammatory cytokines that have been found to inhibit central vagal neurons (Singer et al., 2008).

Never-the-less, different interventional studies on w-3 PUFAs and HRV in patients with heart disease have found inconsistent results, with only 8 out of the 20 trials published so far, supporting a beneficial effect on HRV (Christensen, 2011). Indeed, Mozaffarian et al. (2008) reported that individuals with the highest fish consumption (>5 meals/week) only exhibited 1.5 ms greater HRV compared to those with the lowest fish consumption and further that this modest reduction in HRV was associated with only a 1.1% reduction in the relative risk for sudden cardiac death. Reasons for such inconsistency might include heterogeneous populations, limited sample sizes or different study protocols with variable administered doses of w-3 PUFA and length of intervention. Furthermore, different methods of measurement of HRV with variable time of recordings could be an additional confounder. Finally, an animal study with administration of w-3 PUFAs in rabbits showed that a reduction in pacemaker funny

current rather than an alteration in autonomic neural regulation was responsible for heart rate reduction and increase in HRV (Verkerk et al., 2009). However, such experiments were performed in denervated hearts, excluding a potential impact of autonomic tone on HRV. More recently, the HRV response to physiological challenges was not altered by dietary w-3 fatty acids in conscious intact preparations; data that further suggest that these lipids elicited alterations in pacemaker rate rather than cardiac autonomic regulation (Billman and Harris, 2011; Billman, 2012).

Recent evidence suggests that HMG-CoA reductase inhibitors (statins) have pleiotrophic mechanisms in patients with heart failure, such as ANS output modulation (Lefer, 2002). Experiments with animal models of heart failure have found decreased sympathetic activation and autonomic balance restoration with statins, using HRV analysis (Pliquett et al., 2002). In a crossover study of HRV in 30 patients with hyperlipidemia (Welzig et al., 2003), pravastatin administration induced a significant increase in HF power of ECG signals, whereas others (Vrtovec et al., 2005) found that administration of 10 mg of atorvas-tatin for 3 months, was associated with significant increase in SDNN and RMSSD, in 80 patients with CHF and hyperlipi-demia. Moreover, cholesterol lowering was not correlated with HRV changes, suggesting another mechanism than that of lipid-lowering of statins. In this context, Gao et al. (2005) showed that in experimental heart failure states there is intense free radical production and up-regulation of angiotensin receptors, in autonomic areas of the brain. Moreover, simvastatin therapy was proven to inhibit angiotensin II and superoxide pathways in the RVLM of pacing-induced heart failure rabbits, leading to an abolished renal sympathetic nerve activity (Gao et al., 2005).

Different experimental studies have shown that cate-cholamines, except from increasing cardiac contractility and heart rate via interaction with beta adreno-receptors, may induce myocardial damage by calcium overload and subsequent cell necrosis, upon excessive P-adrenoreceptor stimulation (Opie et al., 1985; Mann et al., 1992). Although, toxic cardiac effects of catecholamines have been recognized since 1907, Rona et al. (1959) was the first who observed that isoproterenol injection into rats produced "infarct-like" myocardial necrosis, in the absence of coronary artery lesions. He proposed the theory of "relative hypoxia" as a pathophysiological mechanism, suggesting a possible imbalance between oxygen demand and blood flow, after excessive adrenergic stimulation. Fleckenstein (1971) thought that calcium overload was the result of catecholamine-mediated cell injury, due to extensive activation of Ca-dependent ATPases and subsequent high energy phosphate deficiency, leading to mitochondrial impairment. On the contrary, Opie and co-workers (1985) were the first who demonstrated that catecholamine cell injury was due to calcium overload, mediated by the P-adrenoreceptor. More recent studies have confirmed previous results and have also demonstrated that other mechanisms could be responsible as well, for catecholamine-induced cell injury, such as increased fibrosis of the left ventricle with associated hypertrophy (Briest et al., 2001) or myocardial cellular apoptosis (Commural et al., 1998).

Since heart failure is associated with a sympathetic up-regulation and parasympathetic withdrawal, p blockers have been used to modify the effects of augmented sympathetic tone and restore autonomic imbalance. In this respect, Goldsmith et al. (1997) showed that administration of carvedilol for 4 months was associated with a significant increase in HF power, in patients with CHF under digoxin and angiotensin-converting enzyme (ACE) inhibitors. Similar results were found in post-MI subjects by Lampert et al. (2003), after treatment with propranolol for 6 weeks. An elevation in HF, which reflects restoration of sympatho-vagal balance, was found to increase final outcome (LVEF, exercise capacity, death and CHF development) in both studies.

Recent evidence suggests that a primary site of attenuated vagal control on the heart occurs at the level of the parasympathetic gaglion (Bibevski and Dunlap, 2011). Thus, cervical vagus nerve stimulation (VNS) has been recently assessed as an "add-on" therapy to optimal medical management of CHF. Li and colleagues (2004) were the first who found that VNS performed for 10 s every minute, in rats developed HF after anterior myocar-dial infarction, improved significantly left ventricular function and decreased mortality from 50 to 14%, in comparison with sham treated animals. Zhang et al. (2009) recently demonstrated in a canine model with high rate ventricular pacing induced HF, that VNS over 12 weeks was able to reduce left ventricular end systolic and end-diastolic volumes and increase LVEF significantly. In addition, HRV was significantly improved in VNS dogs whereas plasma norepinephrine and CRP levels were markedly attenuated with VNS treatment. Finally, vagal stimulation has also been found to limit infarct size and inflammatory response to myocardial ischemia and reperfusion in male rats that underwent myocardial ischemia for 30 min and reperfusion for 24 h (Calvillo et al., 2011). According to the authors, the anti-inflammatory and anti-apoptotic properties of the nicotinic pathway were the primary underlying mechanism in the VNS-treated animals.

Based on the results of VNS in animal models of HF, Schwartz et al. (2008) and De Ferrari et al. (2011) were the first who assessed feasibility and safety and tested possible efficacy of chronic VNS in HF patients with New York Heart Association (NYHA) class II-IV symptoms. In a two-staged study, (8-patients feasibility phase plus 24-patients safety and tolerability phase) they used CardioFit, a right cervical VNS implantable system delivering pulses synchronous with heart beats through a multiple contact bipolar cuff electrode. VNS was started 2-4 weeks after implant and patients were followed 1, 3, and 6 months thereafter. VNS was well tolerated whereas, there was a significant improvement in NYHA class and left ventricular end-systolic volume. Moreover, a significant increase in HRV, estimated with pNN50 and slightly but significantly reduced heart rate were found 6 months after VNS onset (Schwartz et al., 2008; De Ferrari et al., 2011).

As a consequence of such preliminary results, a pivotal multicenter international clinical trial, the INOVATE-HF study, has been recently designed to assess safety and efficacy of VNS in 650 CHF patients from 80 sites, with NYHA class III symptoms, sinus rhythm and QRS width less than 120 ms, using the CardioFit system (Hauptman et al., 2012). Thus, in case of

significant decrease in mortality, vagal stimulation will add significant value to current medical therapy in a narrow spectrum of patients with heart failure, through restoration of sympatho-vagal balance. Never-the-less, a lot of questions remain to be addressed, such as optimal stimulation mode (i.e., right vs. left vagus nerve stimulation, continuous vs. pulse-synchronous stimulation etc) or effective and tolerable VNS dose, before adopting this new non-pharmacologic treatment to our therapeutic armamentarium.


Different experimental studies have established an inter-relation between ANS output and inflammatory regulation, whereas the discovery of the "cholinergic anti-inflammatory pathway" has expanded our understanding of how the nervous system modulates the inflammatory response through an immunoreflex. Furthermore, clinical data from large epidemiological studies involving patients with CAD, heart failure and healthy subjects with increased risk factors for heart diseases suggest that there is a rather weak or moderate association between inflammation and ANS activity, estimated through HRV analysis. Finally, investigation of HRV alterations during critical illness, such as sepsis and MODS, has demonstrated loss of variability of heart rate signals that is inversely correlated with immune response, particularly in most severe cases.

Early and more accurate monitoring of cardiovascular patients, particularly in the early stages of life-threatening illnesses through continuous automated detection of abnormal variability of heart rate signals, could alert clinicians to impending clinical deterioration and allow earlier intervention. In this respect, the combination of structural indices, such as the left ventricular ejection fraction, with autonomic function indices derived from heart rate variability analysis has been proposed as the state-of-the-art method for risk assessment among patients with acute myocardial infarction or severe congestive heart failure (Priori et al., 2001). However, in-consistent findings from different studies assessing the relationship between HRV frequency components and inflammation limit adoption of HRV analysis as an indirect estimator of inflammatory response, since there is

a marked heterogeneity in study protocols, time and methods of HRV measurement and patients' characteristics. Moreover, HRV analysis does not simply reflect sympathetic/parasympathetic balance, since HRV data can be influenced by artifacts related to differences in breathing characteristics, genetic factors or basal autonomic tone. Thus, any change in LF/HF ratio may correspond to central, baroreflex or cellular membrane effects of different stimuli during severe stress. Furthermore, sympathetic outflow can either induce or inhibit inflammatory activity, whereas HF component might fail to reflect vagal inputs upon different organs, such as the spleen, which are major cytokine producers. In addition, the explanatory power of HRV analysis is affected by circadian variability as well as by the estimation of inflammatory activity through biomarkers' measurements from a single blood sample (Haensel et al., 2008). Finally, lack of signals' sta-tionarity (stable statistical properties) during measurements and the strong association between HRV data and inflammation that has been demonstrated among the most severely ill patients could reflect a non-linear relationship between ANS and inflammatory response. Thus, it has been suggested that newer methods derived from chaos theory should be implemented to assess ANS output, such as fractal and detrended fluctuation analysis (DFA) of heart rate signals (Goldberger, 1996). In addition, standardization of experimental protocols and methods used for the estimation of HRV is urgently needed, in order to allow comparisons among different studies.

In conclusion, the enormous complexity of neuro-immunological interactions cannot be captured by simple measurements, such as HRV analysis. A rather multi-parameter monitoring of ANS, through different estimators of heart rate variability and complexity has been suggested for assessment of ANS output (Goldberger, 1996). In any case and since HRV metrics are not enough for differentiating between patho-physiological states (poor specificity) or between patients (poor sensitivity), longitu-dinal alterations over time of both HRV and inflammatory markers on an intra-individual basis must be tested for establishing a potential added value of HRV analysis, as an indirect estimator of inflammatory response (Papaioannou et al., 2009).


Abuissa, H., O'Keefe, J. H., Harris, W. Jr., and Lavie, C. J. (2005). Autonomic function, Omega-3, and cardiovascular risk. Chest 127, 1088-1091.

Akselroad, S., Gordon, D., Madwed, J. B., Shannon, D. C., Barger, A. C., and Cohen, R. J. (1981). Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control. Science 213, 220-222.

Albert, C. M., Ma, J., Rifai, N., Stampfer, M. J., and Ridker, P. M. (2002). Prospective study of C-reactive protein, homocystein, and plasma lipid levels as predictors of sudden cardiac death. Circulation

105, 2595-2599. doi: 10.1161/01. CIR.0000017493.03108.1C Aronson, D., Mittleman, M. A., and Burger, A. J. (2001). Interleukin-6 levels are inversely correlated with heart rate variability in patients with decompensated heart failure. J. Cardiovasc. Electrophysiol. 12, 294-300. doi: 10.1046/j.1540-8167.2001.00294.x Bernardi, L., Salvucci, F., Suardi, R., Solda, P. L., Salciati, A., Perlini, S., et al. (1990). Evidence for an intrinsic mechanism regulating heart rate variability in the transplanted and the intact heart during submaximal dynamic exercise? Cardiovasc. Res. 24, 969-981. doi: 10.1093/cvr/24. 12.969

Berthhoud, H. R., and Neuhuber, W. L. (2000). Functional anatomy of afferent vagal system. Auton. Neurosci. 85, 1-17. doi: 10.1016/S1566-0702(00)00215-0 Bibevski, S., and Dunlap, M. E. (2011). Evidence for impaired vagus nerve activity in heart failure. Heart Fail Rev. 16, 129-135. doi: 10.1007/s10741-010-9190-6 Bilchick, K. C., Fetics, B., Djoukeng, R., Fisher, S. G., Fletcher, R. D., Singh, S. N., et al. (2002). Prognostic value of heart rate variability in chronic congestive heart failure (Veterans Affairs' Survival Trial of Anti-arrhythmic Therapy in Congestive Heart Failure). Am. J. Cardiol. 90, 24-28.

Billman, G. E. (2013). The LH/HF ratio does not accurately measure cardiac sympatho-vagal balance. Front. Physiol. 4:26. doi: 10.3389/fphys.2013.00026 Billman, G. E. (2012). Effect of dietary omega-3 polyunsaturated fatty acids on heart rate and heart rate variability in animals susceptible or resistant to ventricular fibrillation. Front. Physiol. 3:71. doi: 10.3389/fphys.2012.00071 Billman, G. E., and Harris, W. S. (2011). Effect of dietary omega-3 fatty acids on heart rate and the heart rate variability responses to myocardial ischemia or exercise. Am. J. Physiol. Heart Circ. Physiol. 300, H2288-H2299. doi: 10.1152/ajpheart.00140.2011

Briest, W., Holzl, A., Raßler, B., Deten, A., Leicht, M., Bada, H. A., et al. (2001). Cardiac remodeling after long term norepinephrine treatment in rats. Cardiovasc. Res. 52, 265-273.

Brunner, E. J., Hemingway, H., Walker, B. R., Page, M., Clarke, P., Juneja, M., et al. (2002). Adreno-cortical, autonomic,

and inflammatory causes of the metabolic syndrome: nested case-control study. Circulation 106, 2659-2665. doi: 10.1161/01.CIR. 0000038364.26310.BD Buchman, T. G. (2002). The community of the self. Nature 420, 246-251. doi: 10.1038/ nature01260 Calvillo, L., Vanoli, E., Andreoli, E., Besana, A., Omodeo, E., Gneechi, M., et al. (2011). Vagal stimulation, through its nicotinic action limits infarct size and the inflammatory response to myocardial ischemia and reperfusion. J. Cardiovasc. Pharmacol. 58, 500-507. doi: 10.1097/FJC.0b013e31822b7204 Christensen, J. H. (2011). Omega-3 polyunsaturated fatty acids and heart rate variability. Front. Physiol. 2:84. doi: 10.3389/fphys.2011.00084 Christensen, J. H., Christensen, M. S., Dyerberg, J., and Schmidt, E. B. (1999). Heart rate variability and fatty acid content of blood cell membranes: a dose-response study with n-3 fatty acids. Am. J. Clin. Nutr. 70,331-337. Chrousos, G. P. (1995). The hypothalamic-pituitary-adrenal axis and immune-mediated inflammation. New Engl. J. Med. 332, 1351-1362. doi: 10.1056/NEJM199505183322008 Chung, M. K., Gulik, T. S., Rotondo, R. E., Schreiner, G. F., and Lange, L. G. (1990). Mechanisms of action of cytokine inhibition of ß-adrenergic agonist stimulation of cyclic AMP in rat cardiac myocytes: impairment of signal transduc-tion. Circ. Res. 67, 753-763. doi: 10.1161/01.RES.67.3.753 Commural, C., Singh, K., Pimentel, D. R., and Colucci, W. S. (1998). Norepinephrine stimulates apopto-sis in adult rat ventricular myocytes by activation of the ß-adrenergic pathway. Circulation 98, 1329-1334. doi: 10.1161/01.CIR.98.13.1329 DeBoer, R. W., Karemaker, J. M., Strackee, J. (1987). Hemodynamic fluctuations and baroreflex sensitivity in humans: a beat-to-beat model. Am. J. Physiol. 253, 680-689. De Ferrari, G. M., Crijns, H. J., Borggrefe, M., Milasinovic, G., Smid, J., Zabel. M., et al. (2011).

Chronic vagus nerve stimulation: a new and promising therapeutic approach for chronic heart failure. Eur. Heart J. 32, 847-855. doi: 10.1093/eurheartj/ehq391

de Jonge, W. J., van der Zanden, E. P., Bijlsma, M. F., van Westerloo, D. J., Bennink, R. J., Berthoud, H. R., et al. (2005). Stimulation of the vagus nerve attenuates macrophage activation by activating the Jak2-STAT3 signaling pathway. Nat. Immun. 6, 844-851. doi: 10.1038/ ni1229

Dekker, J. M., Schouten, E. G., Klootwijk, P., Pool, J, Swenne, C. A., and Kromhout, D. et al. (1997). Heart rate variability from short electro-cardiographic recordings predicts mortality from all causes in middle-aged and elderly men. Am. J. Epidemiol. 145, 899-908. doi: 10.1093/oxfordjournals.aje.a009049

Eckberg, D. L. (1997). Sympathovagal balance. Circulation 96, 3224-3232. doi: 10.1161/01.CIR.96.9.3224

Elenkov, I. J., Wilder, R. L., Chrousos,

G. P., and Vizi, E. S. (2000). The sympathetic nerve-an inte-grative interface between two super-systems: the brain and the immune system. Pharmacol. Rev. 52, 595-638.

Fairchild, K. D., Saucerman, J. J., Raynor, L. L., Sivak, J. A., Xiao, Y., Lake, D. E., et al. (2009). Endotoxin depresses heart rate variability in mice: cytokine and steroid effects. Am. J. Physiol. Regul. Integr. Comp. Physiol. 297, 1019-1027. doi: 10.1152/ajpregu.00132.2009

Fleckenstein, A. (1971). "Specific inhibitors and promoters of calcium action in the excitation-contraction coupling of heart muscle and their role in the prevention or production of myocardial lesions," in Calcium and the Heart, ed H. P. Opie (London, Academic Press), 135-188.

Frasure, S. N., and Lesperance, F. (2006). Depression and coronary artery disease. Herz 31, 64-68.

Friedland, J. S., Porter, J. C., Daryanani, S., Bland, J. M., Screaton, N. J., Vesely, M. J., et al. (1996). Plasma proinflammatory cytokine concentrations, APACHE III scores and survival in patients in an intensive care unit. Crit. Care Med. 24, 1775-1781. doi: 10.1097/00003246199611000-00003

Gao, L., Wang, W., Li, Y-L., Schultz,

H. D., Liu, D., Cornisch, K. G., et al. (2005). Simvastatin therapy normalizes sympathetic neural control in experimental heart failure. Roles of angiotensin II type

1 receptors and NAD(P)H oxidase. Circulation 112, 1763-1770. doi: 10.1161/CIRCULATIONAHA. 105.552174 Garder, M. L., and von Karel, R. (2006). Myocardial infarction and post-traumatic stress disorder: frequency, outcome, and atherosclerotic mechanisms. Eur. J. Cardiovasc. Prev. Rehabil. 13, 165-172. doi: 10.1097/ 01.hjr.0000214606.60995.46 Godin, P. J., and Buchman, T.

G. (1996). Uncoupling of biological oscillators: a complementary hypothesis concerning the pathogenesis of multiple organ dysfunction syndrome. Crit. Care Med. 24, 1107-1116. doi: 10.1097/00003246-19960700000008

Goldberg, J., Curran, B., Vitek, M. E., Henderson, W. G., and Boyko, E. J. (2002). The vietnam era twin registry. Twin Res. 5, 476-481. Goldberger, A. L. (1996). Non-linear dynamics for clinicians: chaos theory, fractals, and complexity at the bedside. Lancet 347, 1312-1314. doi: 10.1016/S0140-6736(96)90948-4 Goldberger, A. L., Amaral, L. A. N., Hausdorff, J. M., Ivanov, P. C., Peng, C. K., and Stanley, H. E. (2002). Fractal dynamics in physiology: alterations with disease and aging. Proc. Natl. Acad. Sci. U.S.A. 99, 2466-2472. doi: 10.1073/pnas.012579499 Goldberger, J. J., Challapalli, S., Tung, R., Parker, M. A., and Kadish, A.

H. (2001). Relationship of heart rate variability to parasympathetic effect. Circulation 103, 1977-1983. doi: 10.1161/01.CIR.103.15.1977

Goldsmith, R. L., Bigger, J. T., Bloomfield, D. M., Krum, H., Steinman, R. C., Sackner-Bernstein, J., et al. (1997). Long-term carvedilol therapy increases parasympathetic nervous system activity in chronic heart failure. Am. J. Cardiol. 80, 1101-1104. doi: 10.1016/S0002-9149(97)00616-4 Goldstein, B., and Buchman, T. G. (1998). Heart rate variability in intensive care. Intensive Care Med. 13.252-265 Goldstein, B., Fiser, D. H., Kelly, M. M., Mickelsen, D., Ruttiman, U., and Pollack, M. M. (1998). Decomplexification in critical illness and injury: relationship between heart rate variability, severity of illness, and outcome. Crit. Care Med. 26, 352-557. doi: 10.1097/00003246-19980200000040

Griffin, M. P., Lake, D. E., Bissonette, E. A., Harrell, F. E., Micheal, O'Shea,

T., et al. (2005). Heart rate characteristics: novel physiomarkers to predict neonatal infection and death. Pediatrics 116, 1070-1074. Haensel, A., Mills, P. J., Nelesen, R. A., Ziegler, M. G., and Dimsdale, J. E. (2008). The relationship between heart rate variability and inflammatory markers in cardio-vascular diseases. Psychoneuroendocrinology 33, 1305-1312. doi: 10.1016/ j.psyneuen.2008.08.007 Hamaad, A., Sosin, M., Blann, A. D., Patel, J., Lip, G. Y. H., and MacFadyen, R. G. (2005). Markers of inflammation in acute coronary syndromes: association with increased heart rate and reductions in heart rate variability. Clin. Cardiol. 28, 570-576. doi: 10.1002/clc.4960281207 Hauptman, P. J., Schwartz, P. J., Gold, M. R., Borggrefe, M., Van Veldhuisen, D. J., Starling, R. C., et al. (2012). Rational and study design of the INcrease Of Vagal TonE in Heart Failure Study: INOVATE-HF. Am. Heart J. 163, 954-962. doi: 10.1016/j.ahj.2012.03.021 Hrushesky, W. J., Fader, D., Schmitt, O., and Gilbertsen, V. (1984). The respiratory sinus arrhythmia: a measure of cardiac age. Science 224, 1001-1004. doi: 10.1126/science.6372092 Huston, J. M., Ochani, M., RosasBallina, M., Liao, H., Ochani, K., Pavlov, V. A., et al. (2006). Splenectomy inactivates the cholinergic anti-inflammatory pathway during lethal endotox-emia and polymicrobial sepsis. J. Exp. Med. 203, 1623-1628. doi: 10.1084/jem.20052362 Janszky, I., Ericson, M., Lekander, M., Blom, M., Buhlin, K., Georgiades, A., et al. (2004). Inflammatory markers and heart rate variabilityin women with coronary heart disease. J. Intern. Med. 256, 421-428. doi: 10.1111/j.1365-2796.2004.01403.x Kleiger, R. E., Miller, J. P., and Bigger, J. T. (1987). Decreased heart rate variability and its association with increased mortality after acute myocardial infarction. Am. J. Cardiol. 59, 256-262. doi: 10.1016/0002-9149(87)90795-8 Koj, A. (1997). Initiation of acute phase response and synthesis of cytokines. Biochim. Biophys. Acta 1317, 84-94.

Kop, W. J., Stein, P. K., Tracy, R. P., Barzilay, J. I., Schulz, R., and Gottdiener, J. S. (2010). Autonomic nervous system dysfunction and inflammation contribute to the increased cardiovascular mortality

risk associated with depression. Psychosom. Med. 72, 626-635. doi: 10.1097/PSY.0b013e3181eadd2b

Kox, M., Ramakers, B. P., Pompe, J. C., van der Hoeven, J. C., Hoedemaekers, C. W., and Pickkers, P (2011). Interplay between the acute inflammatory response and heart rate variability in healthy human volunteers. Shock 36, 115-120. doi: 10.1097/SHK.0b013e31821c2330

La Rovera, M. T., Bigger, J. T., Marcus, F. I., Mortara, A., Maestri, R., and Schwartz, P. J. (1998). Baroreflex sensitivity and heart rate variability in prediction of total cardiac mortality after myocardial infarction. ATRAMI (Autonomic Tone and Reflexes After Myocardial Infarction) Investigators.

Lancet 351, 478-484. doi: 10.1016/S0140-6736(97)11144-8

Lampert, J. R., Ickovics, R., Viscoli,

C. J., Horwitz, R. I., and Lee, F. A. (2003). Effects of propranolol on recovery of heart rate variability following acute myocardial infarction and relation to outcome in the beta-blocker heart attack trial. Am. J. Cardiol. 91, 137-142. doi: 10.1016/S0002-9149(02)03098-9

Lampert, R., Bremmer, J. D., Su, S., Miller, A., Lee, F., Cheema, F., et al. (2008). Decreased heart rate variability is associated with higher levels of inflammation in middle-aged men. Am. Heart J. 156, 759e1-759e7.

Lanfranchi, P. A., and Somers, V. K. (2002). Arterial baroreflex function and cardio-vascular variability: interactions and implications. Am. J. Physiol. Regul. Integr. Comp. Physiol. 283, 815-826.

Lanza, G. A., Sgueglia, G. A., Cianflone,

D., Rebuzzi, A. G., Angeloni, G., Sestito, A., et al. (2006). Relation of heart rate variability to serum levels of C-reactive protein in patients with unstable angina pectoris. Am. J. Cardiol. 97, 1702-1706. doi: 10.1016/j.amjcard.2006.01.029

Lefer, D. J. (2002). Statins as potent anti-inflammatory

drugs. Circulation 106,

2041-2042. doi: 10.1161/01. CIR.0000033635.42612.88

Li, M., Zheng, C., Sato, T., Kawada, T., Sugimachi, M., and Sunagawa, K. (2004). Vagal nerve stimulation markedly improves long-term survival after chronic heart failure in rats. Circulation 109, 120-124. doi: 10.1161/01.CIR. 0000105721.71640.DA

Lombardi, F., Sandrone, G., Pernptuner, S., Sala, R., Garimoldi,

M., Cerutti, S., et al. (1987). Heart rate variability as an index of sympathovagal interaction after myocardial infarction. Am. J. Cardiol. 60, 1239-1245. doi: 10.1016/0002-9149(87)90601-1 Madsen, T., Christensen, J. H., Toft, E., and Schmidt, E. B. (2007). C-reactive protein is associated with heart rate variability. Ann. Noninvasive Electrocardiol. 12, 216-222. doi: 10.1111/j.1542-474X.2007.00164.x Malave, H. A., Taylor, A. A., Nattama, J., Deswal, A., and Mann, D. L. (2003). Circulating levels of tumor necrosis factor correlate with indexes of depressed heart rate variability. Chest 123, 716-724. doi: 10.1378/chest. 123.3.716

Malik, M., and Camm, A. J. (1993). Components of heart rate variability: what they really mean and what we really measure. Am. J. Cardiol. 72, 821-822 doi: 10.1016/0002-9149(93)91070-X Malpas, S. C. (2002). Neural influences on cardiovascular variability: possibilities and pitfalls. Am. J. Physiol. Heart Circ. Physiol. 282, H6-H20. Mann, D. L., Kent, R. L., Parsons, B., and Cooper, G. (1992). Adrenergic effects on the biology of the adult mammalian cardiocyte. Circulation 85, 790-804. doi: 10.1161/01.CIR.85.2.790 Montano, N., Gnecchi-Ruscone, T., Porta, A., Lombardi, F., Malliani, A., and Barman, S. M. (1996). Presence of vasomotor and respiratory rhythms in the discharge of single medullary neurons involved in the regulation of cardiovascular system. J. Auton. Nerv. Syst. 57, 116-122. doi: 10.1016/0165-1838(95)00113-1 Mozaffarian, D., Stein, P. K., Prineas, R. J., and Siscovick, D. S. (2008). Dietary Fish co-3 fatty acid consumption and heart rate variability in US adults. Circulation 117, 1130-1137. doi: 10.1161/ CIRCULATI0NAHA.107.732826 Muller-Werdan, U., Buerke, M., Ebelt, H., Heinroth, K. M., Herklotz, A., Loppnow, H., et al. (2006). Septic cardiomyopathy-a not yet discovered cardiomyopathy? Exp. Clin. Cardiol. 11,226-236. Nolan, R. P., Reid, G. J., Seidelin, P. H., and Lau, H. K. (2007). C-reactive protein modulates vagal heart rate control in patients with coronary artery disease. Clin. Sci. 112, 449-456. doi: 10.1042/CS20060132 0pie, L. H., Walpoth, B., and Barsacchi, R. (1985). Calcium and catecholamines: Rele-vance

to cardiomyopathies and significance to therapeutic strategies. J. Mol. Cell Cardiol. 17, 21-34. doi: 10.1016/0022-2828(85)90005-7 Owen, N., and Steptoe, A. (2003). Natural killer cell and proinflam-matory cytokine responses to mental stress: associations with heart rate and heart rate variability. Biol. Psychol. 63, 101-115. doi: 10.1016/S0301-0511(03)00023-1 Papaioannou, V., Dragoumanis, C., Theodorou, V., Gargaretas, C., and Pneumatikos, I. (2009). Relation of heart rate variability to serum levels of C-reactive protein, interleukin 6 and 10 in patients with sepsis and septic shock. J. Crit. Care 24, 625e1-625e7.

Phillips, A. N., Neaton, J. D., Cook, D. G., Grimm, R. H., and Shaper, A. G. (1992). Leukocyte count and risk of major coronary heart disease events. Am. J. Epidemiol. 136, 59-70. Pizzi, C., Manzoli, L., Mancini, S., and Maria-Costa, G. (2008). Analysis of potential predictors of depression among coronary heart disease risk factors including heart rate variability, markers of inflammation and endothelial function. Eur. Heart J. 29, 1110-1117. doi: 10.1093/eur-heartj/ehn137 Pliquett, R. U., Cornish, K. G., and Zucker, I. H. (2002). Statin therapy restores sympathovagal balance in experimental heart failure. J. Appl. Physiol. 95, 700-704. Pontes-Arruda, A., Aragao, A. M., and Albuquerque, J. D. (2006). Effects of enteral feeding with eicosa-pentaenoic acid, gamma-linolenic acid, and anti-oxidants in mechanically ventilated patients with severe sepsis and septic shock. Crit. Care Med. 34, 2325-2333. doi: 10.1097/ 01.CCM.0000234033.65657.B6 Priori, S. G., Aliot, E., Blomstrom-Lundqvist, C., Bossaert, L., Breithardt, G., Brugada, P., et al. (2001). Task force on sudden cardiac death of the european society of cardiology. Eur. Heart J. 16, 1374-1450. doi: 10.1053/euhj.2001.2824 Psychari, S. N., Apostolou, T. S., Iliodromitis, E. K., Kourakos, P., Liakos, G., and Kremastinos, D. T. (2007). Inverse relation of C-reactive protein levels to heart rate variability in patients after acute myocardial infarction. Hellenic. J. Cardiol. 48, 64-71. Reichilin, S. (1993). Neuroendocrine-immune interactions. New Engl. J. Med. 329, 1246-1253. doi: 10.1056/NEJM199310213291708 Ridker, P. M., Cushman, M., Stampfer, M. J., Tracy, R. P., and Hennekens,

C. H. (1997). Inflammation, aspirin and the risk of cardiovascular disease in apparently healthy men. New Engl. J. Med. 336, 973-979. doi: 10.1056/NEJM199704033361401 Rona, G., Chappel, G. I., Balazs, T., and Gaudry, R. (1959). An infarct-like myocardial lesion and other toxic manifestations produced by isoproterenol in the rat. AMAArch. Pathol 67, 443-455. Ross, T. (1993). The pathogenesis of atherosclerosis: a perspective for the 1990s. Nature 362, 801-809. doi: 10.1038/362801a0 Sajadieh, A., Nielsen, O. W., Rasmussen, V., Hein, H. O., Abedini, S., and Hansen, J. F. (2004). Increased heart rate and reduced heart-rate variability are associated with subclinical inflammation in middle-aged and elderly subjects with no apparent heart disease. Eur. Heart J. 25, 363-370. doi: 10.1016/j.ehj.2003.12.003 Schwartz, P. J., De Ferrari, G. M., Sanzo, A., Landolina, M., Rordorf, R., Raineri, C., et al. (2008). Long term vagal stimulation in patients with advanced heart failure. First experience in man. Eur. J. Heart Fail 10, 884-891. doi: 10.1016/j.ejheart.2008.07.016 Seely, A. J. E., and Christou, N. V. (2000). Multiple organ dysfunction syndrome: exploring the paradigm of complex nonlinear systems. Crit. Care Med. 28, 2193-2200. doi: 10.1097/00003246-20000700000003

Shehab, A. M., MacFadyen, R. J., McLaren, M., Tavendale, R., Belch, J. J., and Struthers, A. D. (2004). Sudden unexpected death in heart failure may be preceded by short term, intra-individual increases in inflammation and in autonomic dys-function: a pilot study. Heart 90, 1263-1268. doi: 10.1136/hrt.2003.028399 Sherry, R. M., Cue, J. I., Gobbard, J. K., Parramore, J. B., and DiPiro, J. T. (1996). Interleukin-10 is associated with the development of sepsis in trauma patients. J. Trauma. 40, 613-616. doi: 10.1097/00005373199604000-00016 Singer, P., Shapiro, H., Theilla, M., Anbar, R., Singer, J., and Cohen, J. (2008). Anti-inflammatory properties of omega-3 fatty acids in critical illness: novel mechanisms and an integrative perspective. Intensive Care Med. 34, 1580-1592. doi: 10.1007/s00134-008-1142-4 Singer, P., Theilla, M., Fisher, H., Gibstein, L., Grozovski, E., and Cohen, J. (2006). Benefit of an enteral diet enriched with

eico-sapentaenoic acid and gamma-linolenic acid in ventilated patients with acute lung injury. Crit. Care Med. 34, 1033-1038. doi: 10.1097/ 01.CCM.0000206111.23629.0A Sloan, R. P., McCreath, H., Tracey, K. J., Sidney, S., Liu, K., and Seeman, T. (2007). RR interval variability is inversely related to inflammatory markers: the CARDIA study. Mol. Med. 13, 178-184. Sporn, M. B. L. (1997). The importance of context in cytokine action. Kidney Int. 51, 1352-1354. doi: 10.1038/ki.1997.184 Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. (1996). Standards of measurement, physiological interpretation and clinical use Circulation 93, 1043-1065. doi: 10.1161/01.CIR.93.5.1043 Tateishi, Y., Oda, S., Nakamura, M., Watanabe, K., Kuwaki, T., Moriguchi, T., et al. (2007). Depressed heart rate variability is associated with high IL-6 blood level and decline in blood pressure in septic patients. Shock 28, 549-553. doi: 10.1097/shk.0b013e 3180638d1 Tracey, K. J. (2002). The inflammatory reflex. Nature 420, 853-859. doi: 10.1038/nature01321 Tracey, K. J. (2007). Physiology and immunology of the choliner-gic anti-inflammatory pathway. J. Clin. Invest. 117, 289-296. doi: 10.1172/JCI30555 Tsuji, H., Venditti, F. J., Manders, E. S., Evans, J. C., Larson, M. G., Feldman, C. L., et al. (1994). Reduced heart

rate variability and mortality risk in an elderly cohort: the Framingham Heart Study. Circulation 90, 878-883. doi: 10.1161/01.CIR. 90.2.878

van der Poll, T., Coyle, S. M., Barbosa, K., Braxton, C. C., and Lowry, S. F. (1996). Epinephrine inhibits tumor necrosis factor-a and potentiates interleukin-10 production during human endotoxemia. J. Clin. Invest. 97, 713-719. doi: 10.1172/JCI118469

Verkerk, A. O., den Ruijiter, H. M., Bourier, J., Boukens, B. J., Brouwer, I. A., Wilders, R., et al. (2009). Dietary fish oil reduces pacemaker current and heart rate in rabbit. Heart Rhythm 6, 1485-1492. doi: 10.1016/j.hrthm.2009.07.024

von Känel, R., Carney, R. M., Zhao, S., and Whooley, M. A. (2011). Heart rate variability and biomark-ers of systemic inflammation in patients with stable coronary heart diseases: finding from the Heart and Soul Study. Clin. Res. Cardiol. 100, 241-247. doi: 10.1007/s00392-010-0236-5

Vrtovec, B., Okrajsek, R., Golisnik, A., Ferjan, M., Starc, V., and Radovancevic, B. (2005). Atorvastatin therapy increases heart rate variability, decreases QT variability and shortens QTc interval duration in patients with advanced chronic heart failure. J. Card. Fail. 11,684-690.

Wang, H., Liao, H., Ochani, M., Justiniani, M., Lin, X., Yang, L., et al. (2004). Cholinergic agonists inhibit HMGB1 release and improve survival in experimental

sepsis. Nat. Med. 10, 1216-1221. doi: 10.1038/nm1124 Wang, H., Yu, M., Ochani, M., Amelia, CA., Tanovic, M., Susarla, S., et al. (2003). Nicotinic acetylcholine receptor a7 subunit is an essential regulator of inflammation. Nature 421, 384-388. doi: 10.1038/nature01339 Webster, J. I., Tonelli, L., and Sternberg, E. M. (2002). Neuroendocrine regulation of immunity. Annu. Rev. Immunol. 20, 125-163. doi: 10.1146/annurev.immunol.20.0824 01.104914 Welzig, C. M., Shin, D.-G., Park, H.-J., Kim, Y. J., Saul, J. P., and Galper, J. B. (2003). Lipid lowering by pravastatin increases parasympa-thetic modulation of heart rate: G(alpha) i2, a possible molecular marker for parasympathetic responsiveness. Circulation 108, 2743-2746. doi: 10.1161/01.CIR. 0000103680.61390.16 Wolf, M. M., Varigos, G. A., Hunt, D., and Sloman, J. G. (1978). Sinus arrhythmia in acute myocardial infarction. Med. J. Aust. 2, 52-53.

Zaza, A., and Lombardi, F (2001). Autonomic indexes based on the analysis of heart rate variability: a view from the sinus node. Cardiovsc. Res. 50, 434-442. doi: 10.1016/S0008-6363(01)00240-1 Zhang, Y., Popovic, Z. B., Bibevski, S., Fakhry, I., Sica, D. A., Van Wagoner, D. R., et al. (2009). Chronic vagus nerve stimulation improves autonomic control and attenuates systemic inflammation and heart failure progression in a

canine high-rate pacing model. Circ. Heart Fail 2, 692-699. Zhou, M., Yang, S., Koo, D. J., Ornan, D. A., Chaudry, I. H., and Wang, P. (2001). The role of Kupffer cell a2 adrenoceptors in norepinephrine-induced TNF-a production. Biochim. Biophys. Acta 1537, 49-57. doi: 10.1016/S0925-4439(01)00055-2

Conflict of Interest Statement: The

authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Received: 09 January 2013; accepted: 20 June 2013; published online: 10 July 2013.

Citation: Papaioannou V, Pneumatikos I and Maglaveras N (2013) Association of heart rate variability and inflammatory response in patients with cardiovascular diseases: current strengths and limitations. Front. Physiol. 4:174. doi: 10.3389/ fphys.2013.00174

This article was submitted to Frontiers in Cardiac Electrophysiology, a specialty of Frontiers in Physiology. Copyright © 2013 Papaioannou, Pneumatikos and Maglaveras. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.

Copyright of Frontiers in Physiology is the property of Frontiers Media S.A. and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.