This document is a Task Force report for the standardization of heart rate variability (HRV) measurement, written as HRV has emerged as an important metric amid growing recognition of the relationship between the autonomic nervous system and cardiovascular mortality. It presents various HRV measurement methods, physiological correlates, clinical applications, and future research directions, with particular detail on how HRV can be used for risk assessment after acute myocardial infarction and diagnosis of diabetic neuropathy. It also includes technical recommendations for improving the accuracy of HRV measurements and guidelines for testing commercial equipment, contributing to greater reliability and comparability of HRV analysis.


1. Importance of Heart Rate Variability (HRV) and Background of the Task Force

Over the past 20 years, as significant relationships between the autonomic nervous system and cardiovascular mortality -- particularly sudden cardiac death -- have been revealed, heart rate variability (HRV) has drawn attention as one of the most promising metrics for quantitatively assessing autonomic nervous system activity. HRV refers to variations in the intervals between consecutive heartbeats. While it may appear straightforward to measure, the meaning and significance of the various measurement methods are more complex than expected, potentially leading to erroneous conclusions or overinterpretation.

Recognizing these issues, the European Society of Cardiology and the North American Society of Pacing and Electrophysiology formed a Task Force to promote standardization of HRV analysis. The goals of this Task Force were as follows:

  1. Standardization of terminology and development of definitions
  2. Specification of standard measurement methods
  3. Definition of physiological and pathophysiological correlates
  4. Description of currently appropriate clinical applications
  5. Identification of areas for future research

To achieve these goals, experts in mathematics, engineering, physiology, and clinical medicine convened to advance HRV research and establish standards that enable appropriate comparison and careful interpretation.


2. Historical Background and Development of HRV

The clinical importance of HRV was first recognized in 1965 when Hon and Lee observed that changes in fetal interbeat intervals could predict fetal distress before changes in heart rate itself became apparent. This suggested that HRV could be a more sensitive indicator than heart rate alone.

Subsequently, in the 1970s, Sayers and others focused on the physiological rhythms inherent in heart rate signals, and Ewing and colleagues developed simple RR interval difference tests to detect autonomic neuropathy in diabetic patients. In 1977, Wolf et al. first demonstrated the association between reduced post-myocardial infarction mortality and decreased HRV.

In 1981, Akselrod et al. introduced power spectral analysis of heart rate fluctuations, contributing to the quantitative assessment of heart rate control. Such frequency-domain analysis helped understand the autonomic nervous system background of RR interval fluctuations. In the late 1980s, HRV was confirmed as a powerful and independent predictor of mortality after acute myocardial infarction, further highlighting its clinical significance. In the modern era, the advent of new digital high-frequency 24-hour multichannel ECG recording equipment has given HRV the potential to provide valuable insights into physiological and pathological states and to enhance risk stratification.


3. HRV Measurement Methods: Time-Domain and Frequency-Domain Analysis

HRV is assessed by analyzing the time intervals between heart rate or successive heartbeats (NN intervals or RR intervals), and can be broadly divided into time-domain methods and frequency-domain methods.

3.1. Time-Domain Methods

Time-domain methods are the simplest approach, directly measuring changes in heart rate in units of time. Each QRS complex is detected from the ECG recording to determine normal-to-normal (NN) intervals or instantaneous heart rate.

Statistical Methods

Statistical indices are calculated from instantaneous heart rate or cycle interval data recorded over long periods (typically 24 hours).

  • SDNN (Standard Deviation of NN intervals): The standard deviation of all NN intervals, reflecting overall variability during the recording period. In 24-hour recordings, it includes everything from short-term high-frequency fluctuations to the lowest-frequency components, and is therefore used to estimate total HRV. Since the value varies with recording duration, SDNN values from recordings of different lengths are difficult to compare. It is therefore important to standardize the recording duration to 5 minutes or 24 hours.
  • SDANN (Standard Deviation of the Average NN intervals): The standard deviation of mean NN intervals calculated over 5-minute segments, estimating heart rate changes with cycles longer than 5 minutes.
  • SDNN Index: The mean of the standard deviations of all 5-minute NN interval segments calculated over 24 hours, measuring variability with cycles shorter than 5 minutes.
  • RMSSD (Square Root of the Mean of the Sum of the Squares of Differences Between Adjacent NN intervals): The square root of the mean of the squared differences between consecutive NN intervals, estimating short-term variability and high-frequency fluctuations. Due to superior statistical properties, it is preferred over pNN50 and NN50.
  • NN50 count: The number of pairs of consecutive NN intervals differing by more than 50 ms.
  • pNN50 (proportion of NN50): The ratio of NN50 to the total number of NN intervals.

Geometric Methods

Geometric methods convert the NN interval series into geometric patterns such as histograms or Lorenz plots, then evaluate variability based on the geometric properties of those patterns.

  • HRV Triangular Index: The total number of all NN intervals divided by the maximum of the density distribution. This represents overall HRV and is more influenced by low frequencies.
  • TINN (Triangular Interpolation of NN interval Histogram): The baseline width of a triangular approximation of the NN interval distribution. This also represents total HRV measured over 24 hours and is more influenced by low frequencies.

Geometric methods have the advantage of being less sensitive to the analytical quality of NN interval data, but they require a sufficient amount of NN interval data -- at least 20 minutes (ideally 24 hours) -- for accurate pattern construction, and are unsuitable for assessing short-term HRV changes.

Summary and Recommendations

Since the various time-domain measures are closely correlated with one another, the following four indices are recommended for time-domain HRV assessment:

  1. SDNN (total HRV estimate)
  2. HRV Triangular Index (total HRV estimate)
  3. SDANN (long-term component estimate of HRV)
  4. RMSSD (short-term component estimate of HRV)

The choice of measurement method should depend on the research objectives, and it should be kept in mind that measures representing total HRV cannot be compared across different recording durations.


3.2. Frequency-Domain Methods

Frequency-domain analysis calculates the power spectral density (PSD) of heart rate fluctuations to determine the energy (variance) distributed across each frequency band. This helps in understanding the rhythmic characteristics of the various autonomic nervous system activities that contribute to heart rate regulation.

Spectral Components

Three main spectral components are distinguished in short-term recordings (2-5 minutes):

  • VLF (Very Low Frequency, 0.003-0.04 Hz): The physiological explanation is most ambiguous, and caution is needed when interpreting short-term recordings.
  • LF (Low Frequency, 0.04-0.15 Hz): There is controversy over whether it reflects both sympathetic and vagal activity, but when expressed in normalized units (nu), it is sometimes considered an index of sympathetic modulation.
  • HF (High Frequency, 0.15-0.4 Hz): Primarily driven by vagal (parasympathetic) activity, reflecting changes associated with the respiratory cycle.

LF and HF can be measured in absolute values (ms^2) or normalized units (nu). Normalized units emphasize the balanced behavior of the two branches of the autonomic nervous system and tend to minimize the effect of changes in total power.

In long-term recordings (24 hours), a ULF (Ultra Low Frequency, <=0.003 Hz) component is added in addition to VLF, LF, and HF. The 24-hour spectrum can also be evaluated by fitting a slope (alpha) on a log-log scale. Since the physiological mechanisms of LF and HF components cannot be considered stable throughout a 24-hour period, spectral analysis of an entire long-term recording or averaging results from short segments may obscure detailed information about autonomic modulation. It should be remembered that HRV components measure the degree of autonomic modulation, not the average level of the autonomic nervous system.

Technical Requirements and Recommendations

Several important technical requirements must be met to obtain reliable spectral estimates:

  • Clearly distinguish short-term/long-term ECG analysis: Due to the importance of interpretation, this must always be strictly maintained.
  • Signal stationarity: The mechanisms regulating heart rate must not change during the recording.
  • Adequate sampling rate: 250-500 Hz or above is optimal; at rates below 100 Hz, interpolation may be needed to reduce jitter in R-wave fiducial point estimation.
  • Baseline and trend removal: Must be applied carefully, as they can affect low-frequency spectral components.
  • QRS fiducial point selection: A well-tested algorithm must be used to find stable, noise-independent fiducial points.
  • Handling of abnormal beats and noise: Ectopic beats, arrhythmias, missing data, and noise can distort PSD estimates. Short-term recordings free of abnormal beats should be used when possible, or appropriate interpolation should be applied with its effects taken into account.

Algorithm Standards and Recommendations

The data series used for spectral analysis can be obtained in various ways:

  • Interpolation of RR interval tachogram or DES (Discrete Event Series): Use the RR interval tachogram (RR duration vs. sequential beat number) or interpolate the irregularly time-sampled DES to obtain a continuous signal.
  • Non-parametric methods (FFT-based): For 5-minute recordings, a minimum of 512 points, preferably 1024 points, should be used. The DES interpolation formula, sampling frequency, number of samples used for spectral calculation, and spectral window (e.g., Hann, Hamming) must be specified.
  • Parametric methods: Must include the model type used, number of samples, central frequency of each spectral component, and model order (number of parameters). Statistical measures to test model reliability (e.g., prediction error whiteness test (PEWT), optimal order test (OOT)) should also be computed.

Correlation and Differences Between Time-Domain and Frequency-Domain Measures

In short-term recording analysis under stationary conditions, frequency-domain measures have more experience and theoretical knowledge regarding physiological interpretation than time-domain measures. However, many time- and frequency-domain variables measured over 24 hours are strongly correlated with each other. Since physiological interpretation of 24-hour spectral components faces many difficulties, unless a specific investigation dictates otherwise, frequency-domain analysis results may be considered equivalent to the easier-to-perform time-domain analysis results.


3.3. Rhythmic Pattern Analysis and Nonlinear Methods

Rhythmic Pattern Analysis

Rhythmic pattern analysis of HRV has been proposed to overcome the limitations of time and spectral methods caused by irregularities in the RR interval series. These approaches focus on measuring RR interval blocks determined by the characteristics of the rhythm and examining the relationships between these blocks without considering internal variability.

  • Interval spectrum and spectrum of counts methods: Suitable for investigating the relationship between HRV and variability in other physiological measurements. The interval spectrum is useful for linking RR intervals with variables defined in beat units (e.g., blood pressure), while the spectrum of counts is preferred for correlating RR intervals with continuous signals (e.g., respiration) or the occurrence of specific events (e.g., arrhythmias).
  • 'Peak-valley' procedure: Characterizes various oscillations based on heart rate acceleration or deceleration, wavelength, and/or amplitude. In short- to medium-term recordings, it correlates with frequency components of HRV, but the correlation tends to decrease as wavelength and recording duration increase.

Nonlinear Methods

Nonlinear phenomena are certainly involved in HRV generation, determined by complex interactions of hemodynamic, electrophysiological, and humoral variables as well as autonomic and central nervous regulation. HRV analysis based on nonlinear dynamic methods is presumed to provide valuable information for the physiological interpretation of HRV and risk assessment of sudden death. To date, nonlinear dynamic methods applied to biomedical data (including HRV analysis) have not produced major breakthroughs. This may be because integrative complexity measures are not appropriate for analyzing biological systems or are too insensitive to detect nonlinear perturbations in RR intervals. More encouraging results have emerged using differential complexity measures (e.g., scaling exponent methods), but systematic studies investigating large patient populations using these methods have not yet been conducted. Currently, nonlinear methods are considered potential tools for HRV assessment, but standards are still lacking and the full scope of these methods cannot yet be evaluated. Further technical and interpretive advances are needed before they can be applied in physiological and clinical research.

3.4. Stability and Reproducibility of HRV Measurements

Short-term HRV measurements return rapidly to baseline after transient perturbations such as mild exercise or short-acting vasodilator administration. In contrast, potent stimuli such as maximal exercise or long-acting drug administration may require much longer intervals to return to control values.

There is far less data on the stability of long-term HRV measures obtained from 24-hour ambulatory monitoring, but the limited data suggest excellent stability of HRV measures derived from 24-hour ambulatory monitoring in both normal subjects and post-myocardial infarction and ventricular arrhythmia patients. There is also fragmentary data suggesting that the stability of HRV measures may persist for months to years. Since 24-hour indices are stable and free from placebo effects, they may be ideal variables for evaluating interventional therapies.

3.5. Recording Requirements

ECG Signal

The fiducial point of the ECG tracing that identifies the QRS complex can be based on various approaches, including the maximum of the complex, determining the maximum of an interpolation curve, and matching with a template. The voluntary standards of diagnostic ECG equipment (signal-to-noise ratio, common-mode rejection, bandwidth, etc.) are satisfactory for fiducial point localization. However, a lower upper-band frequency cutoff (approximately 200 Hz) compared to diagnostic equipment can introduce jitter in QRS complex fiducial point recognition, causing RR interval measurement errors. Similarly, a limited sampling rate introduces errors in the HRV spectrum, with greater impact on high-frequency components. Using appropriate interpolation, a 100 Hz sampling rate can be adequate.

When using solid-state storage recorders, data compression techniques must be carefully considered in terms of effective sampling rate and the quality of the reconstruction method.

ECG Recording Duration and Environment

The recording duration in HRV studies is determined by the nature of each investigation. Standardization is particularly necessary in studies investigating the physiological and clinical potential of HRV.

  • Short-term recordings: Frequency-domain methods are preferred over time-domain methods. Recordings should be at least 10 times as long as the wavelength of the component being investigated and should not be extended too long to maintain signal stationarity. For example, approximately 1 minute is needed for HF component assessment and approximately 2 minutes for the LF component. For standardization of short-term HRV studies, 5-minute recordings under stationary conditions are preferred unless the nature of the investigation requires a different design.
  • Long-term recordings: Time-domain methods are ideal for long-term recording analysis (since the low stationarity of heart rate regulation during long-term recordings can make interpretation of frequency-method results difficult). Since a substantial portion of long-term HRV values is determined by day-night differences, long-term recordings analyzed by time-domain methods should include a minimum of 18 hours of analyzable ECG data including the entire night.

Little is known about the influence of the environment (type and nature of physical activity, emotional state) during long-term ECG recordings. Some experimental designs require controlling environmental variables, and the nature of the environment should always be described in each study.

RR Interval Sequence Editing

Inaccuracies in the NN interval sequence are known to significantly affect the results of statistical time-domain and all frequency-domain methods. While rough editing of RR interval data may be sufficient for approximate evaluation of overall HRV using geometric methods, it is not known how precisely editing must be performed to ensure accurate results from other methods. Therefore, when using statistical time-domain and/or frequency-domain methods, manual editing of RR data must be performed to a very high standard to ensure correct identification and classification of all QRS complexes. Automatic "filters" should not be relied upon to ensure the quality of RR interval sequences, as they may function unsatisfactorily and potentially introduce errors.

Standardization Recommendations for Commercial Equipment

  • Standard HRV measures: Commercial equipment for short-term HRV analysis should include non-parametric and preferably parametric spectral analysis. To minimize confusion from reporting time-frequency components of heartbeat-based analysis, analysis based on regular sampling of the tachogram should be provided in all cases. Non-parametric spectral analysis should use a minimum of 512 points, preferably 1024 points, for 5-minute recordings.
  • Long-term recording HRV analysis equipment: Should implement time-domain methods including all four standard measures (SDNN, SDANN, RMSSD, HRV Triangular Index). Frequency analysis should be performed on 5-minute segments, and spectral analysis of the entire 24-hour recording should be performed with comparable spectral sampling precision such as 2^18 points to calculate the full range of HF, LF, VLF, and ULF components.
  • Precision and testing: To ensure the quality of the various equipment involved in HRV analysis, testing independent of manufacturers is needed. Testing should include all recording, playback, and analysis stages, including inaccuracies in QRS complex fiducial point identification. Therefore, it is ideal to test equipment with computer-simulated signals with known HRV properties.

4. Physiological Correlates of HRV

4.1. Autonomic Nervous System Influences on HRV Components

Heart rate and rhythm are primarily controlled by the autonomic nervous system.

  • Parasympathetic (vagal) influence: Slows heart rate through acetylcholine release, which increases K+ conductance and suppresses the 'pacemaker' current If. Due to the abundance of acetylcholinesterase, the effects of vagal stimulation are brief.
  • Sympathetic influence: Activates beta-adrenergic receptors through the release of epinephrine and norepinephrine, increasing ICaL and If, accelerating slow diastolic depolarization, and increasing heart rate.
  • Interaction: At rest, vagal tone predominates, and cardiac cycle variation depends primarily on vagal modulation. Parasympathetic activity appears to override sympathetic effects by reducing norepinephrine release and attenuating the response to adrenergic stimulation.

HRV Components

Under resting conditions, RR interval fluctuations represent fine-tuning of beat-to-beat control mechanisms. Vagal afferent stimulation reflexively excites vagal efferent activity and inhibits sympathetic efferent activity, while the opposite reflex action is mediated by sympathetic afferent activity.

  • HF (High Frequency) component: Primarily a major contributor to vagal efferent activity. This has been observed clinically and experimentally through autonomic manipulation including vagal electrical stimulation, muscarinic receptor blockade, and vagotomy.
  • LF (Low Frequency) component: Interpretation is controversial. Some studies consider it a quantitative index of sympathetic modulation (especially when expressed in normalized units), while others view it as reflecting both sympathetic and vagal influences. This is because absolute LF values can decrease during sympathetic activation due to reduced total power.
  • ULF (Ultra Low Frequency) and VLF (Very Low Frequency) components: Account for 95% of total power in 24-hour recordings, but their physiological correlates remain unknown.

LF and HF expressed in normalized units demonstrate autonomic balance changes well, such as increased LF and decreased HF during 90-degree head-up tilt. Spectral analysis of 24-hour recordings shows that LF and HF exhibit day-night patterns and mutual fluctuations -- LF is higher during the day and HF is higher at night. However, these patterns are not detected in a single spectrum of the entire 24-hour recording or when averaging spectra from short segments.

Summary and Recommendations for HRV Component Interpretation

  • HF: A major contributor to vagal activity.
  • LF: There are opposing views -- one holding that it is a quantitative index of sympathetic modulation when expressed in normalized units, and the other that it reflects both sympathetic and vagal activity. Consequently, the LF/HF ratio is also interpreted as reflecting sympathovagal balance or sympathetic modulation.
  • VLF and ULF: The low-frequency components of HRV require further research for physiological interpretation.
  • HRV measures fluctuations in autonomic input to the heart, not the average level. Therefore, both autonomic withdrawal and excessive sympathetic input lead to reduced HRV.

4.2. HRV Changes Associated with Specific Pathological States

Reduced HRV has been reported in several cardiac and non-cardiac conditions.

Myocardial Infarction (MI)

Reduced HRV after MI reflects decreased vagal activity directed to the heart, which can lead to sympathetic dominance and cardiac electrical instability. Decreased 24-hour SDNN during acute MI is significantly associated with left ventricular dysfunction, peak creatine kinase, and Killip class. Although the mechanisms by which HRV transiently decreases after MI and predicts neural responses to acute MI have not yet been defined, they are likely related to abnormalities in neural activity of cardiac origin.

Spectral analysis of HRV in MI survivors showed reduced total power and individual spectral components, but when LF and HF were calculated in normalized units, increased LF and decreased HF were observed. This may indicate a shift in sympathovagal balance toward sympathetic dominance and decreased vagal tone. In post-MI patients with very reduced HRV, most of the residual energy is distributed in the VLF frequency range, and respiratory-related HF appears small.

Diabetic Neuropathy

Diabetes-associated autonomic neuropathy is characterized by early and widespread nerve fiber degeneration of small nerve fibers. A decrease in time-domain parameters of HRV not only has negative prognostic value but appears to precede clinical manifestations of autonomic neuropathy. Reduced power across all spectral bands is the most common finding in diabetic patients, and failure of LF to increase upon standing reflects impaired sympathetic responses or decreased baroreceptor sensitivity.

Cardiac Transplantation

Very reduced HRV has been reported in recent cardiac transplant patients. The appearance of discrete spectral components in a minority of patients is considered to reflect cardiac reinnervation. This reinnervation can occur within 1-2 years after transplantation and is predominantly sympathetic in origin.

Myocardial Dysfunction

Reduced HRV is consistently observed in heart failure patients. In this condition, characterized by signs of sympathetic activation such as increased heart rate and elevated catecholamine levels, a relationship has been reported between HRV changes and the degree of left ventricular dysfunction. In most patients with very advanced disease and drastic HRV reduction, the LF component was not detected despite clinical signs of sympathetic activation.

Tetraplegia

Patients with chronic complete high cervical spinal cord injuries have intact vagal and sympathetic efferent pathways, but spinal sympathetic neurons lose modulatory control (particularly baroreceptor reflex supraspinal inhibitory input). Early studies showed no detectable LF in tetraplegic patients, suggesting an important role of supraspinal mechanisms in determining the 0.1 Hz rhythm, but recent studies have shown that the LF component can be detected in some tetraplegic patients.

4.3. HRV Changes Induced by Specific Interventions

Attempts to alter HRV after MI stem from numerous observations that post-MI patients with reduced HRV have higher cardiac mortality. While there is an inference that interventions that increase HRV may protect against cardiac mortality and sudden cardiac death, the assumption that HRV changes directly lead to cardiac protection can be hazardous. The goal is improvement of cardiac electrical stability, and HRV is merely an index of autonomic activity.

Beta-Adrenergic Blockers and HRV

Surprisingly, data on the effects of beta-blockers on HRV in post-MI patients are scarce. Despite statistically significant increases, the actual changes are very modest. However, it is noteworthy that beta-blockers prevent the increase in the LF component observed during morning hours.

Antiarrhythmic Drugs and HRV

Flecainide and propafenone have been reported to reduce time-domain HRV measures in patients with chronic ventricular arrhythmias. Although some antiarrhythmic drugs can reduce HRV, it is not known whether these HRV changes have direct prognostic significance.

Scopolamine and HRV

Low-dose muscarinic receptor blockers atropine and scopolamine can paradoxically increase vagal efferent activity, as evidenced by heart rate reduction. Scopolamine markedly increases HRV, indicating that pharmacological manipulation can effectively increase vagal activity. However, efficacy during long-term treatment has not been evaluated.

Thrombolysis and HRV

A study of 95 acute MI patients investigating the effect of thrombolysis on HRV (assessed by pNN50) found that HRV was higher in patients with patent lesion-related arteries 90 minutes after thrombolysis. However, when the full 24 hours were analyzed, this difference was no longer apparent.

Exercise Training and HRV

Exercise training can reduce cardiovascular mortality and sudden cardiac death and is thought to modulate autonomic balance. In experimental studies, HRV (SDNN) increased by 74% after exercise training, and all animals survived a new ischemia test. This demonstrates that exercise training can accelerate the restoration of physiological sympathovagal interaction in post-MI patients.


5. Clinical Applications of HRV

HRV has been the subject of numerous clinical studies across various cardiac and non-cardiac conditions, but general consensus on the practical use of HRV in adult medicine has been reached in only two clinical scenarios: risk prediction after acute myocardial infarction and early warning signs of diabetic neuropathy.

5.1. Risk Assessment After Acute Myocardial Infarction

Since the initial observation that the absence of respiratory sinus arrhythmia in acute MI patients was associated with increased 'in-hospital' mortality, numerous reports have followed demonstrating the prognostic value of HRV assessment in identifying high-risk patients.

  • Powerful predictive indicator: Reduced HRV is a powerful predictor of mortality and arrhythmic complications (e.g., symptomatic sustained ventricular tachycardia) in post-acute MI patients. It has predictive value independent of other post-MI risk stratification factors such as reduced left ventricular ejection fraction, increased ventricular ectopic activity, and presence of late potentials.
  • Mortality prediction: The value of HRV in predicting all-cause mortality is comparable to left ventricular ejection fraction. However, for predicting arrhythmic events (sudden cardiac death and ventricular tachycardia), it is superior to left ventricular ejection fraction.
  • Optimal measurement timing: The time point at which HRV has the highest predictive value post-MI has not yet been comprehensively investigated, but measurement just before discharge (i.e., approximately 1 week after acute MI) is generally recommended. HRV decreases immediately after acute MI, begins recovering within weeks, and reaches a maximum at 6-12 months but does not fully recover.
  • Multivariate risk stratification: The predictive value of HRV alone is moderate. However, when combined with other factors such as mean heart rate, left ventricular ejection fraction, frequency of ventricular ectopic activity, high-resolution ECG parameters (presence of late potentials), and clinical assessment, the positive predictive accuracy of HRV improves significantly in the clinically important sensitivity range (25%-75%). However, the optimal combination of risk factors and corresponding dichotomous cutoffs has not yet been established.

"Reduced HRV is a predictor of mortality and arrhythmic complications, independent of other risk factors. The general consensus is that HRV should be measured approximately 1 week after acute myocardial infarction."

5.2. Diabetic Autonomic Neuropathy (DAN) Assessment

Diabetic autonomic neuropathy (DAN) is characterized by widespread nerve degeneration of small nerve fibers in parasympathetic and sympathetic pathways. Since the estimated 5-year mortality rate is approximately 50% once clinical symptoms of DAN appear, early subclinical detection of autonomic dysfunction is important for risk stratification and subsequent management. Short-term and/or long-term HRV analysis has proven useful for diagnosing DAN.

  • Long-term time-domain measures: HRV calculated from 24-hour Holter recordings is more sensitive for detecting DAN than simple bedside tests (Valsalva maneuver, standing test, deep breathing). The NN50 and SDSD methods have been most commonly used.
  • Frequency-domain measures: Abnormal findings in frequency HRV analysis associated with DAN include:
    1. Reduced power across all spectral bands: The most common finding.
    2. Failure of LF to increase upon standing: Reflects impaired sympathetic response or decreased baroreceptor sensitivity.
    3. Abnormal reduction in total power with unchanged LF/HF ratio.
    4. Left shift of LF central frequency: The physiological significance requires further elucidation.

In advanced neuropathy states, resting supine power spectra often have extremely low amplitudes of all spectral components, making it difficult to separate signal from noise. Therefore, including interventions such as standing or tilt is recommended.

5.3. Other Clinical Potential

Beyond myocardial infarction and diabetic neuropathy, HRV has been studied in various cardiac conditions and other clinical settings including hypertension, chronic heart failure (CHF), cardiac transplantation, chronic mitral regurgitation, mitral valve prolapse, cardiomyopathy, cardiac arrest survivors, ventricular arrhythmias, and supraventricular arrhythmias. As presented in Table 4, HRV has the potential to reflect autonomic nervous system dysfunction and indicate specific pathophysiological changes in these conditions. For example, reduced HF power and increased LF/HF ratio have been observed in CHF patients, suggesting withdrawal of parasympathetic tone and sympathetic dominance, with a tendency for HRV to increase with ACE inhibitor treatment.


6. Future Research Possibilities and Challenges

HRV has considerable potential for assessing the role of autonomic nervous system fluctuations in normal subjects and patients with various cardiovascular and non-cardiovascular diseases. HRV research will enhance understanding of physiological phenomena, drug actions, and disease mechanisms.

6.1. Development of HRV Measurement Methods

Current time-domain methods are adequate for long-term HRV profile assessment, with room for improvement in numerical robustness. Non-parametric and parametric spectral methods are adequate for short-term ECG analysis where there are no transient changes in heart rate regulation.

Going forward, attention should be directed to three areas:

  1. Dynamics and transient changes in HRV: Current capabilities for characterizing and quantifying the dynamics and transient changes in RR interval series are insufficient, and mathematical development is underway. Proper assessment of HRV dynamics will greatly enhance understanding of cardiac cycle regulation and its physiological and pathophysiological correlates.
  2. PP and PR intervals: Little is known about the interaction between PP and PR autonomic modulation. While current technology makes it difficult to accurately locate P-wave fiducial points on the surface ECG, technological advances will enable investigation of PP and PR interval variability in future studies.
  3. Multi-signal analysis: Cardiac cycle regulation is not the sole manifestation of autonomic modulatory mechanisms. Commercial equipment capable of simultaneously recording ECG, respiration, blood pressure, and other signals exists, but there are no widely accepted methods for comprehensive multi-signal analysis.

6.2. Research for Physiological Understanding

Efforts to identify the physiological correlates and biological relevance of the various HRV measures currently in use are needed. While this has been achieved in some cases, such as the HF component, in other cases such as the VLF and ULF components, physiological correlates remain largely unknown. This uncertain knowledge limits the interpretation of associations between these variables and the risk of cardiac events.

6.3. Future Clinical Applications

Normal Reference Values

Large-scale prospective population studies with long-term follow-up are needed to establish normal reference values for HRV across various age and gender subgroups. A recent study from the Framingham Heart Study has shown that HRV provides prognostic information beyond traditional risk factors.

Physiological Phenomena

Evaluating HRV under various circadian rhythm patterns -- normal day-night cycles, reversed day-night cycles (night shifts), and temporarily altered day-night cycles (international travel) -- is an interesting area of research. Autonomic nervous system responses to exercise training and rehabilitation exercise programs, as well as deconditioning during prolonged bed rest and weightlessness, can also be understood through HRV data.

Pharmacological Responses

Many drugs act directly or indirectly on the autonomic nervous system, and HRV can be used to explore the effects of various medications on sympathetic and parasympathetic activity. There is currently little data on the effects of calcium channel blockers, sedatives, anxiolytics, analgesics, narcotics, chemotherapeutic agents, and other drugs on HRV.

Risk Stratification

Both time- and frequency-domain HRV measures calculated from 24-hour and short-term ECG recordings have been used to predict time of death after MI, all-cause mortality, and risk of sudden cardiac death in patients with structural heart disease. Diagnostic equipment capable of measuring HRV alongside the frequency and complexity of ventricular arrhythmias, signal-averaged ECG, ST segment variability, and repolarization heterogeneity could significantly enhance the ability to identify patients at risk for sudden cardiac death and arrhythmic events.

Fetal and neonatal HRV is an important area of investigation that can provide early information about fetal and neonatal distress and identify infants at risk for sudden infant death syndrome (SIDS).

Disease Mechanisms

Using HRV techniques to explore the role of autonomic nervous system changes in disease mechanisms is a promising area of research, particularly in conditions where sympathovagal factors are thought to play an important role. The question of the primary or secondary role of enhanced sympathetic activity in essential hypertension could be addressed through longitudinal studies of initially normotensive subjects.


Conclusion

Heart rate variability (HRV) is an important physiological metric that reflects the complex regulatory processes of the autonomic nervous system. Through decades of research, HRV has demonstrated its value in predicting cardiovascular disease mortality -- particularly after acute myocardial infarction -- and in diagnosing diabetic neuropathy. This report provides practical guidelines for standardizing HRV measurement methods, clarifying physiological interpretation, and clinical application, establishing a foundation for more reliable and effective use of HRV analysis in clinical practice.

Of course, many challenges remain to be addressed. In particular, elucidating the physiological significance of low-frequency components (VLF, ULF), advancing HRV dynamics analysis techniques, and optimizing multivariate risk stratification models combining multiple risk factors are important directions for future research. With continued technological advances and in-depth research, HRV is expected to play a decisive role not only in the diagnosis of various diseases, prognosis prediction, and treatment response assessment, but also in deepening our understanding of normal physiological phenomena. Continued efforts to maximize the potential of HRV analysis are essential, and these will ultimately contribute to improved patient care and public health.

Related writing

HarvestStartup · Engineering LeadershipEnglish

Fighting Cancer: How GitLab Founder Sid Sijbrandij Beat Cancer with Founder Mode

The remarkable story of GitLab founder Sid Sijbrandij who, after being diagnosed with a rare bone cancer in 2022, approached treatment in 'Founder Mode' - applying maximized diagnostics, developing 10+ personalized therapies, and running treatments in parallel to overcome the limitations of modern medicine and achieve remission, pointing toward the future of cancer treatment.

Mar 29, 2026Read more
HarvestData & Decision-Making · HealthEnglish

Comparison of Deep Learning Models for Non-Invasive Interstitial Glucose Prediction Using Wearable Sensors in Healthy Cohorts

This pilot study compares deep learning models for predicting blood glucose (interstitial glucose, IG) using only biosignals from wearable devices in healthy individuals. Without requiring food logs or exercise data, the BiLSTM model achieved the best prediction performance (RMSE 13.42 mg/dL) using only wrist-worn sensor data, opening possibilities for convenient non-invasive blood glucose monitoring in daily life.

Feb 2, 2026Read more
HarvestData & Decision-Making · HealthEnglish

A Deep Learning Framework for Virtual Continuous Glucose Monitoring and Prediction Based on Life-Log Data

This study proposes a 'Virtual CGM (V-CGM)' deep learning model that predicts current and future blood glucose levels using only life-log data such as food intake and exercise records, without actual blood glucose measurements. Built on LSTM-based encoder-decoder architecture with food name embeddings and autocorrelation techniques, the model was developed from data of 171 healthy adults and offers a cost-effective alternative to expensive CGM devices.

Feb 2, 2026Read more