Model-Based Evaluation of Methods for Respiratory Sinus Arrhythmia Estimation

One-week test-retest recordings of resting cardiorespiratory data for reliability analysis

Heart rate variability (HRV) is a key indicator of cardiac autonomic function, making reliable assessment crucial. To examine the test-retest stability of resting HRV in healthy individuals, fifty participants attended two lab sessions within a week, at the same time of day. After a 5-minute acclimatization period, electrocardiogram and respiration were recorded at rest. For validation, average heart rate and RMSSD were assessed over 15 minutes using a validated open-source toolbox. Test-retest agreement was evaluated using intra-class correlation (ICC), and coefficients of variation (CV). Mean heart rate showed high stability (ICC = 0.81, CV = 6%), while RMSSD had lower concordance (ICC = 0.75) and greater variation (CV = 30%). These findings indicate good test-retest agreement for standard HRV features. However, a wide range of methodologies exists for assessing various properties of heart rate dynamics. This database is intended to support other researchers in testing additional HRV metrics to evaluate their reliability in healthy individuals.

Characterization of cardiorespiratory coupling via a variability-based multi-method approach: Application to postural orthostatic tachycardia syndrome

There are several mechanisms responsible for the dynamical link between heart period (HP) and respiration (R), usually referred to as cardiorespiratory coupling (CRC). Historically, diverse signal processing techniques have been employed to study CRC from the spontaneous fluctuations of HP and respiration (R). The proposed tools differ in terms of rationale and implementation, capturing diverse aspects of CRC. In this review, we classify the existing methods and stress differences with the aim of proposing a variability-based multi-method approach to CRC evaluation. Ten methodologies for CRC estimation, namely, power spectral decomposition, traditional and causal squared coherence,\;information transfer, cross-conditional entropy, mixed prediction, Shannon entropy of the latency between heartbeat and inspiratory/expiratory onset, conditional entropy of the phase dynamics, synchrogram-based analysis, pulse-respiration quotient, and joint symbolic dynamics, are considered. The ability of these techniques was exemplified over recordings acquired from patients suffering from postural orthostatic tachycardia syndrome (POTS) and healthy controls. Analyses were performed at rest in the supine position (REST) and during head-up tilt (HUT). Although most of the methods indicated that at REST, the CRC was lower in POTS patients and decreased more evidently during HUT in POTS, peculiar differences stressed the complementary value of the approaches. The multiple perspectives provided by the variability-based multi-method approach to CRC evaluation help the characterization of a pathological state and/or the quantification of the effect of a postural challenge. The present work stresses the need for the application of multiple methods to derive a more complete evaluation of the CRC in humans.

Unanticipated evolution of cardio-respiratory interactions with cognitive load during a Go-NoGo shooting task in virtual reality

The cardiovascular system interacts continuously with the respiratory system to maintain the vital balance of oxygen and carbon dioxide in our body. The interplay between the sympathetic and parasympathetic branches of the autonomic nervous system regulates the aforesaid involuntary functions. This study analyzes the dynamics of the cardio-respiratory (CR) interactions using RR Intervals (RRI), Systolic Blood Pressure (SBP), and Respiration signals after first-order differencing to make them stationary. It investigates their variation with cognitive load induced by a virtual reality (VR) based Go-NoGo shooting task with low and high levels of task difficulty. We use Pearson's correlation-based linear and mutual information-based nonlinear measures of association to indicate the reduction in RRI-SBP and RRI-Respiration interactions with cognitive load. However, no linear correlation difference was observed in SBP-Respiration interactions with cognitive load, but their mutual information increased. A couple of open-loop autoregressive models with exogenous input (ARX) are estimated using RRI and SBP, and one closed-loop ARX model is estimated using RRI, SBP, and Respiration. The impulse responses (IRs) are derived for each input-output pair, and a reduction in the positive and negative peak amplitude of all the IRs is observed with cognitive load. Some novel parameters are derived by representing the IR as a double exponential curve with cosine modulation and show significant differences with cognitive load compared to other measures, especially for the IR between SBP and Respiration.

Cardiorespiratory disturbances in focal impaired awareness seizures: Insights from wearable ECG monitoring

PURPOSE

Seizures are characterized by periictal autonomic changes. Wearable devices could help improve our understanding of these phenomena through long-term monitoring. In this study, we used wearable electrocardiogram (ECG) data to evaluate differences between temporal and extratemporal focal impaired awareness (FIA) seizures monitored in the hospital and at home. We assessed periictal heart rate, respiratory rate, heart rate variability (HRV), and respiratory sinus arrhythmia (RSA).

METHODS

We extracted ECG signals across three time points - five minutes baseline and preictal, ten minutes postictal - and the seizure duration. After automatic Rpeak selection, we calculated the heart rate and estimated the respiratory rate using the ECG-derived respiration methodology. HRV was calculated in both time and frequency domains. To evaluate the influence of other modulators on the HRV after removing the respiratory influences, we recalculated the residual power in the high-frequency (HF) and low-frequency (LF) bands using orthogonal subspace projections. Finally, 5-minute and 30-second (ultra-short) ECG segments were used to calculate RSA using three different methods. Seizures from temporal and extratemporal origins were compared using mixed-effects models and estimated marginal means.

RESULTS

The mean preictal heart rate was 69.95 bpm (95 % CI 65.6 - 74.3), and it increased to 82 bpm, 95 % CI (77.51 - 86.47) and 84.11 bpm, 95 % CI (76.9 - 89.5) during the ictal and postictal periods. Preictal, ictal and postictal respiratory rates were 16.1 (95 % CI 15.2 - 17.1), 14.8 (95 % CI 13.4 - 16.2) and 15.1 (95 % CI 14 - 16.2), showing not statistically significant bradypnea. HRV analysis found a higher baseline power in the LF band, which was still significantly higher after removing the respiratory influences. Postictally, we found decreased power in the HF band and the respiratory influences in both frequency bands. The RSA analysis with the new methods confirmed the lower cardiorespiratory interaction during the postictal period. Additionally, using ultra-short ECG segments, we found that RSA decreases before the electroclinical seizure onset. No differences were observed in the studied parameters between temporal and extratemporal seizures.

CONCLUSIONS

We found significant increases in the ictal and postictal heart rates and lower respiratory rates. Isolating the respiratory influences on the HRV showed a postictal reduction of respiratory modulations on both LF and HF bands, suggesting a central role of respiratory influences in the periictal HRV, unlike the baseline measurements. We found a reduced cardiorespiratory interaction during the periictal period using other RSA methods, suggesting a blockade in vagal efferences before the electroclinical onset. These findings highlight the importance of respiratory influences in cardiac dynamics during seizures and emphasize the need to longitudinally assess HRV and RSA to gain insights into long-term autonomic dysregulation.

Spectral Analysis of Heart Rate Variability in Time-Varying Conditions and in the Presence of Confounding Factors

The tools for spectrally analyzing heart rate variability (HRV) has in recent years grown considerably, with emphasis on the handling of time-varying conditions and confounding factors. Time-frequency analysis holds since long an important position in HRV analysis, however, this technique cannot alone handle a mean heart rate or a respiratory frequency which vary over time. Overlapping frequency bands represents another critical condition which needs to be dealt with to produce accurate spectral measurements. The present survey offers a comprehensive account of techniques designed to handle such conditions and factors by providing a brief description of the main principles of the different methods. Several methods derive from a mathematical/statistical model, suggesting that the model can be used to simulate data used for performance evaluation. The inclusion of a respiratory signal, whether measured or derived, is another feature of many recent methods, e.g., used to guide the decomposition of the HRV signal so that signals related as well as unrelated to respiration can be analyzed. It is concluded that the development of new approaches to handling time-varying scenarios are warranted, as is benchmarking of performance evaluated in technical as well as in physiological/clinical terms.

Acute Cardiorespiratory Coupling Impairment in Worsening Sleep Apnea-Related Intermittent Hypoxemia

OBJECTIVE

Hypoxic load is one of the main characteristics of obstructive sleep apnea (OSA) contributing to sympathetic overdrive and weakened cardiorespiratory coupling (CRC). Whether this association changes with increasing hypoxic load has remained obscure. Therefore, we aimed to study our hypothesis that increasing hypoxic load acutely decreases the CRC.

METHODS

We retrospectively analyzed the electrocardiography and nasal pressure signals in 5-min segment pairs (n = 36 926) recorded during clinical polysomnographies of 603 patients with suspected OSA. The segment pairs were pooled into five groups based on the hypoxic load severity described with the the total integrated area under the blood oxygen saturation curve during desaturations. In these severity groups, we determined the frequency-domain heart rate variability (HRV) parameters, the HRV and respiratory high-frequency (HF, 0.15-0.4 Hz) peaks, and the difference between those peaks. We also computed the spectral HF coherence between HRV and respiration in the HF band.

RESULTS

The ratio of low-frequency (LF, 0.04-0.15 Hz) to HF power increased from 1.047 to 1.805 (p < 0.001); the difference between the HRV and respiratory HF peaks increased from 0.001 Hz to 0.039 Hz (p < 0.001); and the spectral coherence between HRV and respiration in the HF band decreased from 0.813 to 0.689 (p < 0.001) as the hypoxic load increased.

CONCLUSION AND SIGNIFICANCE

The vagal modulation decreases and CRC weakens significantly with increasing hypoxic load. Thus, the hypoxic load could be utilized more thoroughly in contemporary OSA diagnostics to better assess the severity of OSA-related cardiac stress.

Technical aspects of cardiorespiratory estimation using subspace projections and cross entropy

BACKGROUND

Respiratory sinus arrhythmia (RSA) is a form of cardiorespiratory coupling. Its quantification has been suggested as a biomarker to diagnose different diseases. Two state-of-the-art methods, based on subspace projections and entropy, are used to estimate the RSA strength and are evaluated in this paper. Their computation requires the selection of a model order, and their performance is strongly related to the temporal and spectral characteristics of the cardiorespiratory signals.

OBJECTIVE

To evaluate the robustness of the RSA estimates to the selection of model order, delays, changes of phase and irregular heartbeats as well as to give recommendations for their interpretation on each case.

APPROACH

Simulations were used to evaluate the model order selection when calculating the RSA estimates explained before, as well as 3 different scenarios that can occur in signals acquired in non-controlled environments and/or from patient populations: the presence of irregular heartbeats; the occurrence of delays between heart rate variability (HRV) and respiratory signals; and the changes over time of the phase between HRV and respiratory signals.

MAIN RESULTS

It was found that using a single model order for all the calculations suffices to characterize the RSA estimates correctly. In addition, the RSA estimation in signals containing more than 5 irregular heartbeats in a period of 5 minutes might be misleading. Regarding the delays between HRV and respiratory signals, both estimates are robust. For the last scenario, the two approaches tolerate phase changes up to 54°, as long as this lasts less than one fifth of the recording duration.

SIGNIFICANCE

Guidelines are given to compute the RSA estimates in non-controlled environments and patient populations.

Assessing stress variations in children during the strange situation procedure: comparison of three widely used respiratory sinus arrhythmia estimation methods

Objective.The respiratory sinus arrhythmia (RSA) is a well-known marker of vagal activity that can be exploited to measure stress changes. RSA is usually estimated from heart rate variability (HRV). This study aims to compare the RSA obtained with three widely adopted methods showing their strengths and potential pitfalls.Approach.The three methods are tested on 69 healthy preschoolers undergoing a stressful protocol, the strange situation procedure (SSP). We compare the RSA estimated by the Porges method, the univariate autoregressive (AR) spectral analysis of the HRV signal, and the bivariate AR spectral analysis of HRV and respirogram signals. We examine RSA differences detected across the SSP episodes and correlation between the estimates provided by each method.Main results.The Porges and the bivariate AR approaches both detected significant differences (i.e. stress variations) in the RSA measured across the SSP. However, the latter method showed higher sensitivity to stress changes induced by the procedure, with the mean RSA variation between baseline and first separation from the mother (the most stressful condition) being significantly different among methods: Porges, -17.5%; univariate AR, -18.3%; bivariate AR, -23.7%. Moreover, the performances of the Porges algorithm were found strictly dependent on the applied preprocessing.Significance.Our findings confirm the bivariate AR analysis of the HRV and respiratory signals as a robust stress assessment tool that does not require any population-specific preprocessing of the signals and warn about using RSA estimates that neglect breath information in more natural experiments, such as those involving children, in which respiratory frequency changes are extremely likely.

Cardiopulmonary coupling indices to assess weaning readiness from mechanical ventilation

The ideal moment to withdraw respiratory supply of patients under Mechanical Ventilation at Intensive Care Units (ICU), is not easy to be determined for clinicians. Although the Spontaneous Breathing Trial (SBT) provides a measure of the patients' readiness, there is still around 15-20% of predictive failure rate. This work is a proof of concept focused on adding new value to the prediction of the weaning outcome. Heart Rate Variability (HRV) and Cardiopulmonary Coupling (CPC) methods are evaluated as new complementary estimates to assess weaning readiness. The CPC is related to how the mechanisms regulating respiration and cardiac pumping are working simultaneously, and it is defined from HRV in combination with respiratory information. Three different techniques are used to estimate the CPC, including Time-Frequency Coherence, Dynamic Mutual Information and Orthogonal Subspace Projections. The cohort study includes 22 patients in pressure support ventilation, ready to undergo the SBT, analysed in the 24 h previous to the SBT. Of these, 13 had a successful weaning and 9 failed the SBT or needed reintubation -being both considered as failed weaning. Results illustrate that traditional variables such as heart rate, respiratory frequency, and the parameters derived from HRV do not differ in patients with successful or failed weaning. Results revealed that HRV parameters can vary considerably depending on the time at which they are measured. This fact could be attributed to circadian rhythms, having a strong influence on HRV values. On the contrary, significant statistical differences are found in the proposed CPC parameters when comparing the values of the two groups, and throughout the whole recordings. In addition, differences are greater at night, probably because patients with failed weaning might be experiencing more respiratory episodes, e.g. apneas during the night, which is directly related to a reduced respiratory sinus arrhythmia. Therefore, results suggest that the traditional measures could be used in combination with the proposed CPC biomarkers to improve weaning readiness.

Benchmarking Transfer Entropy Methods for the Study of Linear and Nonlinear Cardio-Respiratory Interactions

Transfer entropy (TE) has been used to identify and quantify interactions between physiological systems. Different methods exist to estimate TE, but there is no consensus about which one performs best in specific applications. In this study, five methods (linear, k-nearest neighbors, fixed-binning with ranking, kernel density estimation and adaptive partitioning) were compared. The comparison was made on three simulation models (linear, nonlinear and linear + nonlinear dynamics). From the simulations, it was found that the best method to quantify the different interactions was adaptive partitioning. This method was then applied on data from a polysomnography study, specifically on the ECG and the respiratory signals (nasal airflow and respiratory effort around the thorax). The hypothesis that the linear and nonlinear components of cardio-respiratory interactions during light and deep sleep change with the sleep stage, was tested. Significant differences, after performing surrogate analysis, indicate an increased TE during deep sleep. However, these differences were found to be dependent on the type of respiratory signal and sampling frequency. These results highlight the importance of selecting the appropriate signals, estimation method and surrogate analysis for the study of linear and nonlinear cardio-respiratory interactions.

Linear and Non-linear Quantification of the Respiratory Sinus Arrhythmia Using Support Vector Machines

Respiratory sinus arrhythmia (RSA) is a form of cardiorespiratory coupling. It is observed as changes in the heart rate in synchrony with the respiration. RSA has been hypothesized to be due to a combination of linear and nonlinear effects. The quantification of the latter, in turn, has been suggested as a biomarker to improve the assessment of several conditions and diseases. In this study, a framework to quantify RSA using support vector machines is presented. The methods are based on multivariate autoregressive models, in which the present samples of the heart rate variability are predicted as combinations of past samples of the respiration. The selection and tuning of a kernel in these models allows to solve the regression problem taking into account only the linear components, or both the linear and the nonlinear ones. The methods are tested in simulated data as well as in a dataset of polysomnographic studies taken from 110 obstructive sleep apnea patients. In the simulation, the methods were able to capture the nonlinear components when a weak cardiorespiratory coupling occurs. When the coupling increases, the nonlinear part of the coupling is not detected and the interaction is found to be of linear nature. The trends observed in the application in real data show that, in the studied dataset, the proposed methods captured a more prominent linear interaction than the nonlinear one.