Heart rate variability

Heart rate variability (HRV) is the physiological phenomenon of variation in the time interval between heartbeats. It is measured by the variation in the beat-to-beat interval.

Other terms used include: "cycle length variability", "RR variability" (where R is a point corresponding to the peak of the QRS complex of the ECG wave; and RR is the interval between successive Rs), and "heart period variability".

See also Heart rate turbulence, Sinus rhythm.

Methods used to detect beats include: ECG, blood pressure, ballistocardiograms, and the pulse wave signal derived from a photoplethysmograph (PPG). ECG is considered superior because it provides a clear waveform, which makes it easier to exclude heartbeats not originating in the sinoatrial node. The term "NN" is used in place of RR to emphasize the fact that the processed beats are "normal" beats.

Clinical significance
Reduced HRV has been shown to be a predictor of mortality after myocardial infarction although others have shown that the information in HRV relevant to acute myocardial infarction survival is fully contained in the mean heart rate. A range of other outcomes/conditions may also be associated with modified (usually lower) HRV, including congestive heart failure, diabetic neuropathy, depression, post-cardiac transplant, susceptibility to SIDS and poor survival in premature babies.

Mental and social aspects
In the field of psychophysiology, there is interest in HRV. For example, HRV is related to emotional arousal. High-frequency (HF) activity has been found to decrease under conditions of acute time pressure and emotional strain and elevated state anxiety, presumably related to focused attention and motor inhibition. HRV has been shown to be reduced in individuals reporting a greater frequency and duration of daily worry. In individuals with post-traumatic stress disorder (PTSD), HRV and its HF component (see below) is reduced compared to controls whilst the low-frequency (LF) component is elevated. Furthermore, unlike controls, PTSD patients demonstrated no LF or HF reactivity to recalling a traumatic event.

The Polyvagal Theory  derives from a psychophysiologic imputation of importance to HRV. This theory emphasizes the role of heart rate variability in understanding the magnitude and nature of vagal outflow to the heart. This theory decomposes heart rate variability based on frequency domain characteristics with an emphasis on respiratory sinus arrhythmia and its transmission by a neural pathway that is distinct from other components of HRV. There is anatomic and physiological  evidence for a polyvagal control of the heart.

Variation
Variation in the beat-to-beat interval is a physiological phenomenon. The SA node receives several different inputs and the instantaneous heart rate or RR interval and its variation are the results of these inputs.

The main inputs are the sympathetic and the parasympathetic nervous system (PSNS) and humoral factors. Respiration gives rise to waves in heart rate mediated primarily via the PSNS, and it is thought that the lag in the baroceptor feedback loop may give rise to 10 second waves in heart rate (associated with Mayer waves of blood pressure), but this remains controversial.

Factors that affect the input are the baroreflex, thermoregulation, hormones, sleep-wake cycle, meals, physical activity, and stress.

Decreased PSNS activity or increased SNS activity will result in reduced HRV. High frequency (HF) activity (0.15 to 0.40 Hz), especially, has been linked to PSNS activity. Activity in this range is associated with the respiratory sinus arrhythmia (RSA), a vagally mediated modulation of heart rate such that it increases during inspiration and decreases during expiration. Less is known about the physiological inputs of the low frequency (LF) activity (0.04 to 0.15 Hz). Though previously thought to reflect SNS activity, it is now widely accepted that it reflects a mixture of both the SNS and PSNS.

Heart rate variability phenomena
There are two primary fluctuations:
 * Respiratory arrhythmia (or Respiratory sinus arrhythmia). This heart rate variation is associated with respiration and faithfully tracks the respiratory rate across a range of frequencies.
 * Low-frequency oscillations. This heart rate variation is associated with Mayer waves (Traube–Hering–Mayer waves) of blood pressure and is usually at a frequency of 0.1 Hz or a 10-second period.

HRV artifact
Errors in the location of the instantaneous heart beat will result in errors in the calculation of the HRV. HRV is highly sensitive to artifact and errors in as low as even 2% of the data will result in unwanted biases in HRV calculations. To ensure accurate results therefore it is critical to manage artifact and R-R errors appropriately prior to performing any HRV analyses.

HRV analysis
The most widely used methods can be grouped under time-domain and frequency-domain. Other methods have been proposed, such as non-linear methods.

Time-domain methods
These are based on the beat-to-beat or NN intervals, which are analysed to give variables such as:
 * SDNN, the standard deviation of NN intervals. Often calculated over a 24-hour period. SDANN, the standard deviation of the average NN intervals calculated over short periods, usually 5 minutes. SDANN is therefore a measure of changes in heart rate due to cycles longer than 5 minutes. SDNN reflects all the cyclic components responsible for variability in the period of recording, therefore it represents total variability.


 * RMSSD ("root mean square of successive differences"), the square root of the mean of the sum of the squares of the successive differences between adjacent NNs.


 * SDSD ("standard deviation of successive differences"), the standard deviation of the successive differences between adjacent NNs.


 * NN50, the number of pairs of successive NNs that differ by more than 50 ms.


 * pNN50, the proportion of NN50 divided by total number of NNs.


 * NN20, the number of pairs of successive NNs that differ by more than 20 ms.


 * pNN20, the proportion of NN20 divided by total number of NNs.


 * EBC ("estimated breath cycle"), the range (max-min) within a moving window of a given time duration within the study period. The windows can move in a self-overlapping way or be strictly distinct (sequential) windows. EBC is often provided in data acquisition scenarios where HRV feedback in real time is a primary goal. EBC derived from PPG over 10-second and 16-second sequential and overlapping windows has been shown to correlate highly with SDNN.

Geometric methods
The series of NN intervals also can be converted into a geometric pattern such as the sample density distribution of NN interval durations, sample density distribution of differences between adjacent NN intervals, Lorenz plot of NN or RR intervals, and so forth, and a simple formula is used that judges the variability on the basis of the geometric and/or graphics properties of the resulting pattern.

Frequency-domain methods
Frequency domain methods assign bands of frequency and then count the number of NN intervals that match each band. The bands are typically high frequency (HF) from 0.15 to 0.4 Hz, low frequency (LF) from 0.04 to 0.15 Hz, and the very low frequency (VLF) from 0.0033 to 0.04 Hz.

Several methods of analysis are available. Power spectral density (PSD), using parametric or nonparametric methods, provides basic information on the power distribution across frequencies. One of the most commonly used PSD methods is the discrete Fourier transform. Methods for the calculation of PSD may be generally classified as nonparametric and parametric. In most instances, both methods provide comparable results. The advantages of the nonparametric methods are (1) the simplicity of the algorithm used (fast Fourier transform [FFT] in most of the cases) and (2) the high processing speed, while the advantages of parametric methods are (1) smoother spectral components that can be distinguished independent of preselected frequency bands, (2) easy postprocessing of the spectrum with an automatic calculation of low- and high-frequency power components with an easy identification of the central frequency of each component, and (3) an accurate estimation of PSD even on a small number of samples on which the signal is supposed to maintain stationarity. The basic disadvantage of parametric methods is the need of verification of the suitability of the chosen model and of its complexity (that is, the order of the model).

In addition to classical FFT-based methods used for the calculation of frequency parameters, a more appropriate PSD estimation method is the Lomb–Scargle (LS) periodogram. Analysis has shown that the LS periodogram can produce a more accurate estimate of the PSD than FFT methods for typical RR data. Since the RR data is an unevenly sampled data, another advantage of the LS method is that in contrast to FFT-based methods it is able to be used without the need to resample and detrend the RR data.

A newly used HRV index, which depends on the wavelet entropy measures, is an alternative choice. The wavelet entropy measures are calculated using a three-step procedure defined in the literature. First, the wavelet packet algorithm is implemented using the Daubechies 4 (DB4) function as the mother wavelet with a scale of 7. Once the wavelet coefficients are obtained, the energy for each coefficient are calculated as described in the literature. After calculating the normalized values of wavelet energies, which represent the relative wavelet energy (or the probability distribution), the wavelet entropies are obtained using the definition of entropy given by Shannon.

Non-linear methods
Given the complexity of the mechanisms regulating heart rate, it is reasonable to assume that applying HRV analysis based on methods of non-linear dynamics will yield valuable information. Although chaotic behavior has been assumed, more rigorous testing has shown that heart rate variability cannot be described as a chaotic process. The most commonly used non-linear method of analysing heart rate variability is the Poincaré plot. Each data point represents a pair of successive beats, the x-axis is the current RR interval, while the y-axis is the previous RR interval. HRV is quantified by fitting mathematically defined geometric shapes to the data. Other methods used are the correlation dimension, nonlinear predictability, pointwise correlation dimension and approximate entropy.

Long term correlations
Sequences of RR intervals have been found to have long-term correlations. However, one flaw with these analyses is their lack of goodness-of-fit statistics, i.e. values are derived that may or may not have adequate statistical rigor.

Duration and circumstances of ECG recording
Time domain methods are preferred to frequency domain methods when short-term recordings are investigated. This is due to the fact that the recording should be at least 10 times the wavelength of the lowest frequency bound of interest. Thus, recording of approximately 1 minute is needed to assess the HF components of HRV (i.e., a lowest bound of 0.15 Hz is a cycle of 6.6 seconds and so 10 cycles require ~60 seconds), while more than 4 minutes are needed to address the LF component (with a lower bound of 0.04 Hz).

Although time domain methods, especially the SDNN and RMSSD methods, can be used to investigate recordings of long durations, a substantial part of the long-term variability is day-night differences. Thus, long-term recordings analyzed by time domain methods should contain at least 18 hours of analyzable ECG data that include the whole night.

Autonomic influences of heart rate
Although cardiac automaticity is intrinsic to various pacemaker tissues, heart rate and rhythm are largely under the control of the autonomic nervous system. The parasympathetic influence on heart rate is mediated via release of acetylcholine by the vagus nerve. Muscarinic acetylcholine receptors respond to this release mostly by an increase in cell membrane K+ conductance. Acetylcholine also inhibits the hyperpolarization-activated "pacemaker" current If the "Ik decay" hypothesis65 proposes that pacemaker depolarization results from slow deactivation of the delayed rectifier current, Ik, which, due to a time-independent background inward current, causes diastolic depolarization. Conversely, the "If activation" hypothesis suggests that after action potential termination, If provides a slowly activating inward current predominating over decaying Ik, thus initiating slow diastolic depolarization.

The sympathetic influence on heart rate is mediated by release of epinephrine and norepinephrine. Activation of ß-adrenergic receptors results in cAMP-mediated phosphorylation of membrane proteins and increases in ICaL and in If the end result is an acceleration of the slow diastolic depolarization.

Under resting conditions, vagal tone prevails and variations in heart period are largely dependent on vagal modulation. The vagal and sympathetic activity constantly interact. Because the sinus node is rich in acetylcholinesterase, the effect of any vagal impulse is brief because the acetylcholine is rapidly hydrolyzed. Parasympathetic influences exceed sympathetic effects probably through two independent mechanisms: (1) a cholinergically induced reduction of norepinephrine released in response to sympathetic activity and (2) a cholinergic attenuation of the response to an adrenergic stimulus.

Components of HRV
The RR interval variations present during resting conditions represent beat-by-beat variations in cardiac autonomic inputs. However, efferent vagal (parasympathetic) activity is a major contributor to the HF component, as seen in clinical and experimental observations of autonomic maneuvers such as electrical vagal stimulation, muscarinic receptor blockade, and vagotomy. More problematic is the interpretation of the LF component, which was considered by some as a marker of sympathetic modulation (especially when expressed in normalized units) but is now known to include both sympathetic and vagal influences. For example, during sympathetic activation the resulting tachycardia is usually accompanied by a marked reduction in total power, whereas the reverse occurs during vagal activation. Thus the spectral components change in the same direction and do not indicate that LF faithfully reflects sympathetic effects.

It is important to note that HRV measures fluctuations in autonomic inputs to the heart rather than the mean level of autonomic inputs. Thus, both withdrawal and saturatingly high levels of autonomic input to the heart can lead to diminished HRV.

Changes of HRV related to specific pathologies
A reduction of HRV has been reported in several cardiovascular and noncardiovascular diseases.

Myocardial infarction
Depressed HRV after MI may reflect a decrease in vagal activity directed to the heart. HRV in patients surviving an acute MI reveal a reduction in total and in the individual power of spectral components. The presence of an alteration in neural control is also reflected in a blunting of day-night variations of RR interval. In post-MI patients with a very depressed HRV, most of the residual energy is distributed in the VLF frequency range below 0.03 Hz, with only a small respiration-related variations.

Diabetic neuropathy
In neuropathy associated with diabetes mellitus characterized by alteration in small nerve fibers, a reduction in time domain parameters of HRV seems not only to carry negative prognostic value but also to precede the clinical expression of autonomic neuropathy. In diabetic patients without evidence of autonomic neuropathy, reduction of the absolute power of LF and HF during controlled conditions was also reported.

Cardiac transplantation
A very reduced HRV with no definite spectral components has been reported in patients with a recent heart transplant. The appearance of discrete spectral components in a few patients is considered to reflect cardiac reinnervation. This reinnervation may occur as early as 1 to 2 years after transplantation and is assumed to be of sympathetic origin. In addition, a correlation between respiratory rate and the HF component of HRV observed in some transplanted patients also indicates that a nonneural mechanism may generate a respiration-related rhythmic oscillation.

Myocardial dysfunction
A reduced HRV has been observed consistently in patients with cardiac failure. In this condition characterized by signs of sympathetic activation such as faster heart rates and high levels of circulating catecholamines, a relation between changes in HRV and the extent of left ventricular dysfunction was reported. In fact, whereas the reduction in time domain measures of HRV seemed to parallel the severity of the disease, the relationship between spectral components and indices of ventricular dysfunction appears to be more complex. In particular, in most patients with a very advanced phase of the disease and with a drastic reduction in HRV, an LF component could not be detected despite the clinical signs of sympathetic activation. This reflects that, as stated above, the LF may not accurately reflect cardiac sympathetic tone.s

Tetraplegia
Patients with chronic complete high cervical spinal cord lesions have intact efferent vagal neural pathways directed to the sinus node. However, an LF component can be detected in HRV and arterial pressure variabilities of some tetraplegic patients. Thus, the LF component of HRV in those without intact sympathetic inputs to the heart represent vagal modulation.

Modifications of HRV by specific interventions
Interventions that augment HRV may be protective against cardiac mortality and sudden cardiac death. Although the rationale for changing HRV is sound, it also contains the inherent danger of leading to the unwarranted assumption that modification of HRV translates directly into cardiac protection, which may not be the case. Despite the growing consensus that increases in vagal activity can be beneficial, it is not as yet known how much vagal activity (or HRV as a marker) has to increase in order to provide adequate protection.

ß-Adrenergic blockade and HRV
The data on the effect of ß-blockers on HRV in post-MI patients are surprisingly scant. Despite the observation of statistically significant increases, the actual changes are very modest. In conscious post-MI dogs, ß-blockers do not modify HRV. The unexpected observation that before MI, ß-blockade increases HRV only in the animals destined to be at low risk for lethal arrhythmias after MI may suggest novel approaches to post-MI risk stratification.

Antiarrhythmic drugs and HRV
Data exist for several antiarrhythmic drugs. Flecainide and propafenone but not amiodarone were reported to decrease time domain measures of HRV in patients with chronic ventricular arrhythmia. In another study, propafenone reduced HRV and decreased LF much more than HF. A larger study confirmed that flecainide, also encainide and moricizine, decreased HRV in post-MI patients but found no correlation between the change in HRV and mortality during follow-up. Thus, some antiarrhythmic drugs associated with increased mortality can reduce HRV. However, it is not known whether these changes in HRV have any direct prognostic significance.

Scopolamine and HRV
Low-dose muscarinic receptor blockers, such as atropine and scopolamine, may produce a paradoxical increase in vagal effects on the heart, as suggested by a decrease in heart rate. In addition, scopolamine and low dose atropine can markedly increase HRV. However, though the heart rate slowing in proportional to the (low) does of atropine, the increase in HRV varies widely across and within individuals. This suggests that even for vagal activity to the heart, HRV may be a limited marker.

Thrombolysis and HRV
The effect of thrombolysis on HRV (assessed by pNN50) was reported in 95 patients with acute MI. HRV was higher 90 minutes after thrombolysis in the patients with patency of the infarct-related artery. However, this difference was no longer evident when the entire 24 hours were analyzed.

Exercise training and HRV
Exercise training may decrease cardiovascular mortality and sudden cardiac death. Regular exercise training is also thought to modify cardiac autonomic control. Individuals who exercise regularly have a 'training bradycardia' (i.e., low resting heart rate) and generally have higher HRV than sedentary individuals.

Biofeedback
The technique called resonant breathing biofeedback teaches how to recognize and control involuntary heart rate variability. A randomized study by Sutarto et al. assessed the effect of resonant breathing biofeedback among manufacturing operators; depression, anxiety and stress significantly decreased.

Wind Instruments
One study that surveyed the physiological effects of playing Native American flutes found a significant HRV increase when playing both low-pitched and high-pitched flutes.

Normal values of standard measures of HRV
There are no widely accepted standard values for HRV that can be used for clinical purposes.