Neonatal heart rate variability: a contemporary scoping review of analysis methods and clinical applications

Establishing normative values for short-term heart rate variability indices in healthy infants in the emergency department

Background

Heart rate variability (HRV) has been used as a marker of cardiovascular health and a risk factor for mortality in the adult and paediatric populations, and as an indicator of neonatal sepsis. There has been an increasing interest in using short-term (5 minutes) HRV to identify infants ≤90 days of life with serious bacterial infections. However, there has not been any normative data range reported for short-term HRV indices in this infant population. The aim of this study was to evaluate short-term HRV indices in awake, healthy young infants >48 hours and ≤90 days of life and to establish a reference range. We also aimed to produce a clinical calculator that can be used in this population for evaluation of short-term HRV variables in young infants in the emergency department (ED) setting that can be potentially used in future clinical validation and research.

Methods

We conducted a prospective observational study of short-term HRV analysis of awake, well infants ≤90 days of life in the ED setting.

Results

One hundred and eight infants with complete data [51.9% male, median age 9 days (interquartile range, 4-35 days)] were included. We found that heart rate (HR) is correlated with HRV. Thus, normalisation of HRV parameters was done to remove their dependence on HR. We then provided normative reference range of widely used short-term HRV time-domain, frequency-domain, and non-linear HRV metrics in our cohort.

Conclusions

We established normative values and HRV calculator for evaluation of these short-term HRV variables in young infants in ED settings that can be used for further clinical validation and clinical research.

Exploring Computational Techniques in Preprocessing Neonatal Physiological Signals for Detecting Adverse Outcomes: Scoping Review

BACKGROUND

Computational signal preprocessing is a prerequisite for developing data-driven predictive models for clinical decision support. Thus, identifying the best practices that adhere to clinical principles is critical to ensure transparency and reproducibility to drive clinical adoption. It further fosters reproducible, ethical, and reliable conduct of studies. This procedure is also crucial for setting up a software quality management system to ensure regulatory compliance in developing software as a medical device aimed at early preclinical detection of clinical deterioration.

OBJECTIVE

This scoping review focuses on the neonatal intensive care unit setting and summarizes the state-of-the-art computational methods used for preprocessing neonatal clinical physiological signals; these signals are used for the development of machine learning models to predict the risk of adverse outcomes.

METHODS

Five databases (PubMed, Web of Science, Scopus, IEEE, and ACM Digital Library) were searched using a combination of keywords and MeSH (Medical Subject Headings) terms. A total of 3585 papers from 2013 to January 2023 were identified based on the defined search terms and inclusion criteria. After removing duplicates, 2994 (83.51%) papers were screened by title and abstract, and 81 (0.03%) were selected for full-text review. Of these, 52 (64%) were eligible for inclusion in the detailed analysis.

RESULTS

Of the 52 articles reviewed, 24 (46%) studies focused on diagnostic models, while the remainder (n=28, 54%) focused on prognostic models. The analysis conducted in these studies involved various physiological signals, with electrocardiograms being the most prevalent. Different programming languages were used, with MATLAB and Python being notable. The monitoring and capturing of physiological data used diverse systems, impacting data quality and introducing study heterogeneity. Outcomes of interest included sepsis, apnea, bradycardia, mortality, necrotizing enterocolitis, and hypoxic-ischemic encephalopathy, with some studies analyzing combinations of adverse outcomes. We found a partial or complete lack of transparency in reporting the setting and the methods used for signal preprocessing. This includes reporting methods to handle missing data, segment size for considered analysis, and details regarding the modification of the state-of-the-art methods for physiological signal processing to align with the clinical principles for neonates. Only 7 (13%) of the 52 reviewed studies reported all the recommended preprocessing steps, which could have impacts on the downstream analysis.

CONCLUSIONS

The review found heterogeneity in the techniques used and inconsistent reporting of parameters and procedures used for preprocessing neonatal physiological signals, which is necessary to confirm adherence to clinical and software quality management system practices, usefulness, and choice of best practices. Enhancing transparency in reporting and standardizing procedures will boost study interpretation and reproducibility and expedite clinical adoption, instilling confidence in the research findings and streamlining the translation of research outcomes into clinical practice, ultimately contributing to the advancement of neonatal care and patient outcomes.

Machine learning model with output correction: Towards reliable bradycardia detection in neonates

Bradycardia is a commonly occurring condition in premature infants, often causing serious consequences and cardiovascular complications. Reliable and accurate detection of bradycardia events is pivotal for timely intervention and effective treatment. Excessive false alarms pose a critical problem in bradycardia event detection, eroding trust in machine learning (ML)-based clinical decision support tools designed for such detection. This could result in disregarding the algorithm's accurate recommendations and disrupting workflows, potentially compromising the quality of patient care. This article introduces an ML-based approach incorporating an output correction element, designed to minimise false alarms. The approach has been applied to bradycardia detection in preterm infants. We applied five ML-based autoencoder techniques, using recurrent neural network (RNN), long-short-term memory (LSTM), gated recurrent unit (GRU), 1D convolutional neural network (1D CNN), and a combination of 1D CNN and LSTM. The analysis is performed on ∼440 hours of real-time preterm infant data. The proposed approach achieved 0.978, 0.73, 0.992, 0.671 and 0.007 in AUC-ROC, AUC-PRC, recall, F1 score, and false positive rate (FPR) respectively and a false alarms reduction of 36% when compared with methods without the correction approach. This study underscores the imperative of cultivating solutions that alleviate alarm fatigue and encourage active engagement among healthcare professionals.

Full-Term and Preterm Newborns Differ More Significantly in Photoplethysmographic Waveform Variability than Heart Rate Variability

BACKGROUND

Features of cardiovascular autonomic regulation in infants are poorly studied compared with adults. However, the clinical significance of autonomic dysfunction in infants is very high. The goal of our research was to study the temporal and frequency-dependent features, as well as low-frequency synchronization in cardiovascular autonomic regulation in full-term vs. preterm newborns, based on the analysis of their heart rate variability (HRV) and photoplethysmographic waveform variability (PPGV).

METHODS

The study included three groups of newborns: 64 full-term newborns (with a gestational age at birth of 37-40 weeks) with a physiological course of the neonatal adaptation; 23 full-term newborns (with a gestational age at birth of 37-40 weeks) with a pathological course of the neonatal adaptation; and 17 preterm newborns (with a postconceptional age of 34 weeks or more). We conducted spectral analysis of HRV and PPGV, along with an assessment of the synchronization strength between low-frequency oscillations in HRV and in PPGV (synchronization index). We employed several options for the boundaries of the high-frequency (HF) band: 0.15-0.40 Hz, 0.2-2 Hz, 0.15-0.8 Hz, and 0.24-1.04 Hz.

RESULTS

Preterm newborns had higher heart rate, RMSSD, and PNN50 values relative to both groups of full-term newborns. Values of SDNN index and synchronization index (S index) were similar in all groups of newborns. Differences in frequency domain indices of HRV between groups of newborns depended on the considered options of HF band boundaries. Values of frequency domain indices of PPGV demonstrated similar differences between groups, regardless of the boundaries of considered options of HF bands and the location of PPG signal recording (forehead or leg). An increase in sympathetic influences on peripheral blood flow and a decrease in respiratory influences were observed along the following gradient: healthy full-term newborns → preterm newborns → full-term newborns with pathology.

CONCLUSIONS

Differences in frequency domain indices of autonomic regulation between the studied groups of newborns depended on the boundaries of the considered options of the HF band. Frequency domain indices of PPGV revealed significantly more pronounced differences between groups of newborns than analogous HRV indicators. An increase in sympathetic influences on peripheral blood flow and a decrease in respiratory influences were observed along the following gradient: healthy full-term newborns → preterm newborns → full-term newborns with pathology.

Prediction of outcome of hypoxic-ischemic encephalopathy in newborns undergoing therapeutic hypothermia using heart rate variability

OBJECTIVE

To assess the use of continuous heart rate variability (HRV) as a predictor of brain injury severity in newborns with moderate to severe HIE that undergo therapeutic hypothermia.

STUDY DESIGN

Two cohorts of newborns (n1 = 55, n2 = 41) with moderate to severe hypoxic-ischemic encephalopathy previously treated with therapeutic hypothermia. HRV was characterized by root mean square in the short time scales (RMSS) during therapeutic hypothermia and through completion of rewarming. A logistic regression and Naïve Bayes models were developed to predict the MRI outcome of the infants using RMSS. The encephalopathy grade and gender were used as control variables.

RESULTS

For both cohorts, the predicted outcomes were compared with the observed outcomes. Our algorithms were able to predict the outcomes with an area under the receiver operating characteristic curve of about 0.8.

CONCLUSIONS

HRV assessed by RMSS can predict severity of brain injury in newborns with HIE.

Heart rate variability in school-age children born moderate-to-late preterm

BACKGROUND

Prematurity is associated with reduced cardiac autonomic function. This study aimed to investigate the heart rate variability (HRV) in school-age children born moderately to late preterm (MLPT).

METHODS

This cross-sectional study investigated school-age children, aged 5 to 10 years, born moderate-to-late preterm. Electrocardiograms recordings were performed during fifteen-minutes. Time and frequency domain parameters were calculated, corrected for heart rate and compared between the groups.

RESULTS

A total of 123 children were evaluated and 119 were included in this study. HRV measures, studied in the time and frequency domains, were similar in both groups. Corrected values of root mean square of successive differences between normal cycles (RMSSD), percentage of successive cycles with a duration difference >50 ms (pNN50%), and high frequency (HF), indices that predominantly represent the parasympathetic activity of the autonomic nervous system, were 1.6E-7 and 1.8E-7 (p=0.226); 1.6E-13 and 1.6E-13 (p=0.506); 6.9E-12 and 7.4E-12 (p=0.968) in the preterm and control groups, respectively.

CONCLUSION

This study did not find differences in heart rate variability between school-age children born MLPT and those born at term, suggesting that plasticity of cardiac autonomic modulation continues to occur in children up to school age or there is less impairment of the autonomic system in MLPT.

Comprehensive viewpoints on heart rate variability at high altitude

OBJECTIVES

Hypoxia is a physiological state characterized by reduced oxygen levels in organs and tissues. It is a common clinicopathological process and a major cause of health problems in highland areas.  Heart rate variability (HRV) is a measure of the balance in autonomic innervation to the heart. It provides valuable information on the regulation of the cardiovascular system by neurohumoral factors, and changes in HRV reflect the complex interactions between multiple systems. In this review, we provide a comprehensive overview of the relationship between high-altitude hypoxia and HRV. We summarize the different mechanisms of diseases caused by hypoxia and explore the changes in HRV across various systems. Additionally, we discuss relevant pharmaceutical interventions. Overall, this review aims to provide research ideas and assistance for in-depth studies on HRV. By understanding the intricate relationship between high-altitude hypoxia and HRV, we can gain insights into the underlying mechanisms and potential therapeutic approaches to mitigate the effects of hypoxia on cardiovascular and other systems.

METHODS

The relevant literature was collected systematically from scientific database, including PubMed, Web of Science, China National Knowledge Infrastructure (CNKI), Baidu Scholar, as well as other literature sources, such as classic books of hypoxia.

RESULTS

There is a close relationship between heart rate variability and high-altitude hypoxia. Heart rate variability is an indicator that evaluates the impact of hypoxia on the cardiovascular system and other related systems. By improving the observation of HRV, we can estimate the progress of cardiovascular diseases and predict the impact on other systems related to cardiovascular health. At the same time, changes in heart rate variability can be used to observe the efficacy of preventive drugs for altitude related diseases.

CONCLUSIONS

HRV can be used to assess autonomic nervous function under various systemic conditions, and can be used to predict and monitor diseases caused by hypoxia at high altitude. Investigating the correlation between high altitude hypoxia and heart rate variability can help make HRV more rapid, accurate, and effective for the diagnosis of plateau-related diseases.

The Effect of Caffeine on Heart Rate Variability in Newborns: A Pilot Study

Neonatal apnoea can be treated with caffeine, which affects the central nervous and cardiovascular systems. Heart rate variability (HRV) reflects the activity of the autonomic nervous system (ANS) and might be used as a measure of ANS maturation in newborns. We aimed to establish the effect of caffeine on HRV in newborns and investigated the potential correlation between HRV and postmenstrual age (PMA). In 25 haemodynamically stable newborns hospitalized due to apnoea and treated with caffeine (2.5 mg/kg), we assessed breathing frequency, arterial oxygen saturation, body temperature, and the heart rate while they were sleeping. We assessed HRV by spectral analysis using fast Fourier transformation. The same protocol was reapplied 100 h after caffeine withdrawal to assess the control parameters. Caffeine increased breathing frequency (p = 0.023) but did not affect any other parameter assessed including HRV. We established a positive correlation between postmenstrual age and HRV during treatment with caffeine as well as after caffeine had been withdrawn (total power: p = 0.044; low-frequency band: p = 0.039). Apparently, the maintenance dose of caffeine is too low to affect the heart rate and HRV. A positive correlation between PMA and HRV might reflect maturation of the ANS, irrespective of caffeine treatment.

Methods for analyzing infant heart rate variability: A preliminary study

Heart rate (HR) and heart rate variability (HRV) reflect autonomic development in infants. To better understand the autonomic response in infants, reliable HRV recordings are vital, yet no protocol exists. The purpose of this paper is to present reliability of a common procedure for analysis from two different file types. In the procedure, continuous electrocardiograph recordings of 5-10 min are obtained at rest in infants at 1 month of age by using a Hexoskin Shirt-Junior's (Carre Technologies Inc., Montreal, QC, Canada). Electrocardiograph (ECG; .wav) and R-R interval (RRi; .csv) files are extracted. The RRi of the ECG signal is generated by VivoSense (Great Lakes NeuroTechnologies, Independence, OH). Two MATLAB (The MathWorks, Inc., Natick, MA) scripts converted files for analysis with Kubios HRV Premium (Kubios Oy, Kuopio, Finland). A comparison was made between RRi and ECG files for HR and HRV parameters, and then tested with t tests and correlations via SPSS. There are significant differences in root mean squared successive differences between recording types, with only HR and low-frequency measures significantly correlated together. Recording with Hexoskin and analysis with MATLAB and Kubios enable infant HRV analysis. Differences in outcomes exist between procedures, and standard methodology for infant HR analysis is needed.

Comparing Full and Pre-Term Neonates' Heart Rate Variability in Rest Condition and during Spontaneous Interactions with Their Parents at Home

BACKGROUND

Preterm neonates show decreased HRV compared to those at full-term. We compared HRV metrics between preterm and full-term neonates in transfer periods from neonate rest state to neonate-parent interaction, and vice versa.

METHODS

Short-term recordings of the HRV parameters (time and frequency-domain indices and non-linear measurements) of 28 premature healthy neonates were compared with the metrics of 18 full-term neonates. HRV recordings were performed at home at term-equivalent age and HRV metrics were compared between the following transfer periods: from first rest state of the neonate (TI1) to a period in which the neonate interacted with the first parent (TI2), from TI2 to a second neonate rest state (TI3), and from TI3 to a period of neonate interaction with the second parent (TI4).

RESULTS

For the whole HRV recording period, PNN50, NN50 and HF (%) was lower for preterm neonates compared to full-terms. These findings support the reduced parasympathetic activity of preterm compared to full-term neonates. The results of comparisons between transfer period simply a common coactivation of SNS and PNS systems for both full and pre-term neonates.

CONCLUSIONS

Spontaneous interaction with the parent may reinforce both full and pre-term neonates' ANS maturation.

Donkey Heart Rate and Heart Rate Variability: A Scoping Review

Heart rate (HR) and heart rate variability (HRV) are commonly used physiological measures in animals. While several studies exist on horse HRV, less information is available for donkeys. This scoping review aims to understand the extent and type of published evidence on donkey HR and HRV, their clinical and research applications, the devices used, and the analysis performed. Only quantitative primary studies published in English were considered. Four different databases were queried through the Web of Science platform, with additional evidence identified by citation chasing. After a two-stage screening phase, data were extracted considering study and population characteristics, information on HR/HRV analysis, and applications. The majority of the 87 included articles (about 80%) concerned a sample size of up to 20 individuals and were published since 2011 (about 65%). Forty-one articles employed an electronic device for signal acquisition (mainly electrocardiographs and heart rate monitors), yet only two articles reported HRV parameters. The literature on donkey HRV is lacking, and this gap can be filled by gaining knowledge on donkey characteristics and finding useful tools for welfare assessment. Comparison with what is known about the horse allows a discussion of the technical and interpretative difficulties that can be encountered with donkeys.

Heart rate variability analysis for the prediction of EEG grade in infants with hypoxic ischaemic encephalopathy within the first 12 h of birth

Background and aims

Heart rate variability (HRV) has previously been assessed as a biomarker for brain injury and prognosis in neonates. The aim of this cohort study was to use HRV to predict the electroencephalography (EEG) grade in neonatal hypoxic-ischaemic encephalopathy (HIE) within the first 12 h.

Methods

We included 120 infants with HIE recruited as part of two European multi-centre studies, with electrocardiography (ECG) and EEG monitoring performed before 12 h of age. HRV features and EEG background were assessed using the earliest 1 h epoch of ECG-EEG monitoring. HRV was expressed in time, frequency and complexity features. EEG background was graded from 0-normal, 1-mild, 2-moderate, 3-major abnormalities to 4-inactive. Clinical parameters known within 6 h of birth were collected (intrapartum complications, foetal distress, gestational age, mode of delivery, gender, birth weight, Apgar at 1 and 5, assisted ventilation at 10 min). Using logistic regression analysis, prediction models for EEG severity were developed for HRV features and clinical parameters, separately and combined. Multivariable model analysis included 101 infants without missing data.

Results

Of 120 infants included, 54 (45%) had normal-mild and 66 (55%) had moderate-severe EEG grade. The performance of HRV model was AUROC 0.837 (95% CI: 0.759-0.914) and clinical model was AUROC 0.836 (95% CI: 0.759-0.914). The HRV and clinical model combined had an AUROC of 0.895 (95% CI: 0.832-0.958). Therapeutic hypothermia and anti-seizure medication did not affect the model performance.

Conclusions

Early HRV and clinical information accurately predicted EEG grade in HIE within the first 12 h of birth. This might be beneficial when EEG monitoring is not available in the early postnatal period and for referral centres who may want some objective information on HIE severity.