Improved analysis of heart rate variability by methods of nonlinear dynamics

CVRanalysis: a free software for analyzing cardiac, vascular and respiratory interactions

Introduction: Simultaneous beat-to-beat R-R intervals, blood pressure and respiration signals are routinely analyzed for the evaluation of autonomic cardiovascular and cardiorespiratory regulations for research or clinical purposes. The more recognized analyses are i) heart rate variability and cardiac coherence, which provides an evaluation of autonomic nervous system activity and more particularly parasympathetic and sympathetic autonomic arms; ii) blood pressure variability which is mainly linked to sympathetic modulation and myogenic vascular function; iii) baroreflex sensitivity; iv) time-frequency analyses to identify fast modifications of autonomic activity; and more recently, v) time and frequency domain Granger causality analyses were introduced for assessing bidirectional causal links between each considered signal, thus allowing the scrutiny of many physiological regulatory mechanisms. Methods: These analyses are commonly applied in various populations and conditions, including mortality and morbidity predictions, cardiac and respiratory rehabilitation, training and overtraining, diabetes, autonomic status of newborns, anesthesia, or neurophysiological studies. Results: We developed CVRanalysis, a free software to analyze cardiac, vascular and respiratory interactions, with a friendly graphical interface designed to meet laboratory requirements. The main strength of CVRanalysis resides in its wide scope of applications: recordings can arise from beat-to-beat preprocessed data (R-R, systolic, diastolic and mean blood pressure, respiration) or raw data (ECG, continuous blood pressure and respiratory waveforms). It has several tools for beat detection and correction, as well as setting of specific areas or events. In addition to the wide possibility of analyses cited above, the interface is also designed for easy study of large cohorts, including batch mode signal processing to avoid running repetitive operations. Results are displayed as figures or saved in text files that are easily employable in statistical softwares. Conclusion: CVRanalysis is freely available at this website: anslabtools.univ-st-etienne.fr. It has been developed using MATLAB® and works on Windows 64-bit operating systems. The software is a standalone application avoiding to have programming skills and to install MATLAB. The aims of this paper area are to describe the physiological, research and clinical contexts of CVRanalysis, to introduce the methodological approach of the different techniques used, and to show an overview of the software with the aid of screenshots.

Multiscale Cumulative Residual Dispersion Entropy with Applications to Cardiovascular Signals

Heart rate variability (HRV) is used as an index reflecting the adaptability of the autonomic nervous system to external stimuli and can be used to detect various heart diseases. Since HRVs are the time series signal with nonlinear property, entropy has been an attractive analysis method. Among the various entropy methods, dispersion entropy (DE) has been preferred due to its ability to quantify the time series' underlying complexity with low computational cost. However, the order between patterns is not considered in the probability distribution of dispersion patterns for computing the DE value. Here, a multiscale cumulative residual dispersion entropy (MCRDE), which employs a cumulative residual entropy and DE estimation in multiple temporal scales, is presented. Thus, a generalized and fast estimation of complexity in temporal structures is inherited in the proposed MCRDE. To verify the performance of the proposed MCRDE, the complexity of inter-beat interval obtained from ECG signals of congestive heart failure (CHF), atrial fibrillation (AF), and the healthy group was compared. The experimental results show that MCRDE is more capable of quantifying physiological conditions than preceding multiscale entropy methods in that MCRDE achieves more statistically significant cases in terms of p-value from the Mann-Whitney test.

Heart Rate Variability for Evaluating Psychological Stress Changes in Healthy Adults: A Scoping Review

The utility of heart rate variability (HRV) for characterizing psychological stress is primarily impacted by methodological considerations such as study populations, experienced versus induced stress, and method of stress assessment. Here, we review studies on the associations between HRV and psychological stress, examining the nature of stress, ways stress was assessed, and HRV metrics used. The review was performed according to the PRISMA guidelines on select databases. Studies that examined the HRV-stress relationship via repeated measurements and validated psychometric instruments were included (n = 15). Participant numbers and ages ranged between 10 and 403 subjects and 18 and 60 years, respectively. Both experimental (n = 9) and real-life stress (n = 6) have been explored. While RMSSD was the most reported HRV metric (n = 10) significantly associated with stress, other metrics, including LF/HF (n = 7) and HF power (n = 6) were also reported. Various linear and nonlinear HRV metrics have been utilized, with nonlinear metrics used less often. The most frequently used psychometric instrument was the State-Trait Anxiety Inventory (n = 10), though various other instruments have been reported. In conclusion, HRV is a valid measure of the psychological stress response. Standard stress induction and assessment protocols combined with validated HRV measures in different domains will improve the validity of findings.

Caffeine increases performance and leads to a cardioprotective effect during intense exercise in cyclists

The present study was designed to investigate the effects of different caffeine dietary strategies to compare the impact on athletic performance and cardiac autonomic response. The order of the supplementation was randomly assigned: placebo(4-day)-placebo(acute)/PP, placebo(4-day)-caffeine(acute)/PC and caffeine(4-day)-caffeine(acute)/CC. Fourteen male recreationally-trained cyclists ingested capsules containing either placebo or caffeine (6 mg kg-1) for 4 days. On day 5 (acute), capsules containing placebo or caffeine (6 mg kg-1) were ingested 60 min before completing a 16 km time-trial (simulated cycling). CC and PC showed improvements in time (CC vs PP, Δ - 39.3 s and PC vs PP, Δ - 43.4 s; P = 0.00; ƞ2 = 0.33) and in output power (CC vs PP, Δ 5.55 w and PC vs PP, Δ 6.17 w; P = 0.00; ƞ2 = 0.30). At the final of the time-trial, CC and PC exhibited greater parasympathetic modulation (vagal tone) when compared to the PP condition (P < 0.00; ƞ2 = 0.92). Our study provided evidence that acute caffeine intake (6 mg∙kg-1) increased performance (time-trial) and demonstrated a relevant cardioprotective effect, through increased vagal tone.

CEPS: An Open Access MATLAB Graphical User Interface (GUI) for the Analysis of Complexity and Entropy in Physiological Signals

BACKGROUND

We developed CEPS as an open access MATLAB^® GUI (graphical user interface) for the analysis of Complexity and Entropy in Physiological Signals (CEPS), and demonstrate its use with an example data set that shows the effects of paced breathing (PB) on variability of heart, pulse and respiration rates. CEPS is also sufficiently adaptable to be used for other time series physiological data such as EEG (electroencephalography), postural sway or temperature measurements.

METHODS

Data were collected from a convenience sample of nine healthy adults in a pilot for a larger study investigating the effects on vagal tone of breathing paced at various different rates, part of a development programme for a home training stress reduction system.

RESULTS

The current version of CEPS focuses on those complexity and entropy measures that appear most frequently in the literature, together with some recently introduced entropy measures which may have advantages over those that are more established. Ten methods of estimating data complexity are currently included, and some 28 entropy measures. The GUI also includes a section for data pre-processing and standard ancillary methods to enable parameter estimation of embedding dimension m and time delay τ ('tau') where required. The software is freely available under version 3 of the GNU Lesser General Public License (LGPLv3) for non-commercial users. CEPS can be downloaded from Bitbucket. In our illustration on PB, most complexity and entropy measures decreased significantly in response to breathing at 7 breaths per minute, differentiating more clearly than conventional linear, time- and frequency-domain measures between breathing states. In contrast, Higuchi fractal dimension increased during paced breathing.

CONCLUSIONS

We have developed CEPS software as a physiological data visualiser able to integrate state of the art techniques. The interface is designed for clinical research and has a structure designed for integrating new tools. The aim is to strengthen collaboration between clinicians and the biomedical community, as demonstrated here by using CEPS to analyse various physiological responses to paced breathing.

Entropy-Based Measures of Hypnopompic Heart Rate Variability Contribute to the Automatic Prediction of Cardiovascular Events

Surges in sympathetic activity should be a major contributor to the frequent occurrence of cardiovascular events towards the end of nocturnal sleep. We aimed to investigate whether the analysis of hypnopompic heart rate variability (HRV) could assist in the prediction of cardiovascular disease (CVD). 2217 baseline CVD-free subjects were identified and divided into CVD group and non-CVD group, according to the presence of CVD during a follow-up visit. HRV measures derived from time domain analysis, frequency domain analysis and nonlinear analysis were employed to characterize cardiac functioning. Machine learning models for both long-term and short-term CVD prediction were then constructed, based on hypnopompic HRV metrics and other typical CVD risk factors. CVD was associated with significant alterations in hypnopompic HRV. An accuracy of 81.4% was achieved in short-term prediction of CVD, demonstrating a 10.7% increase compared with long-term prediction. There was a decline of more than 6% in the predictive performance of short-term CVD outcomes without HRV metrics. The complexity of hypnopompic HRV, measured by entropy-based indices, contributed considerably to the prediction and achieved greater importance in the proposed models than conventional HRV measures. Our findings suggest that Hypnopompic HRV assists the prediction of CVD outcomes, especially the occurrence of CVD event within two years.

Ventilatory pattern variability as a biometric for severity of acute lung injury in rats

We hypothesize that ventilatory pattern variability (VPV) varies with the magnitude of acute lung injury (ALI). In adult male rats, we instilled a low- or high- dose of bleomycin or saline (PBS) intratracheally. While representative samples of pulmonary tissue indicated graded lung injury, coefficient of variation (CV) of TTOT did not differ among the 3 groups. Broncho-alveolar lavage fluid (BALF), respiratory rate (fR), mutual information were greater in ALI than sham rats; but did not differ between bleomycin doses. However, nonlinear complexity index (NLCI), which is the difference in sample entropy between original and surrogate data sets was greater for high- versus low- dose; but did not differ between low-dose and sham groups. Further, NLCI correlated to an injury index based on protein concentration of BALF and failure to gain weight. Finally, Receiver Operator Curves (ROCs) indicated that both mutual information and NLCI had greater sensitivity and specificity than fR and CVTTOT in identifying ALI. Thus, nonlinear analyses of VPV can distinguish ALI and out performs fR as a biometric.

Multiscale Information Decomposition Dissects Control Mechanisms of Heart Rate Variability at Rest and During Physiological Stress

Heart rate variability (HRV; variability of the RR interval of the electrocardiogram) results from the activity of several coexisting control mechanisms, which involve the influence of respiration (RESP) and systolic blood pressure (SBP) oscillations operating across multiple temporal scales and changing in different physiological states. In this study, multiscale information decomposition is used to dissect the physiological mechanisms related to the genesis of HRV in 78 young volunteers monitored at rest and during postural and mental stress evoked by head-up tilt (HUT) and mental arithmetics (MA). After representing RR, RESP and SBP at different time scales through a recently proposed method based on multivariate state space models, the joint information transfer TRESP,SBP→RR is decomposed into unique, redundant and synergistic components, describing the strength of baroreflex modulation independent of respiration (USBP→RR), nonbaroreflex (URESP→RR) and baroreflex-mediated (RRESP,SBP→RR) respiratory influences, and simultaneous presence of baroreflex and nonbaroreflex respiratory influences (SRESP,SBP→RR), respectively. We find that fast (short time scale) HRV oscillations—respiratory sinus arrhythmia—originate from the coexistence of baroreflex and nonbaroreflex (central) mechanisms at rest, with a stronger baroreflex involvement during HUT. Focusing on slower HRV oscillations, the baroreflex origin is dominant and MA leads to its higher involvement. Respiration influences independent on baroreflex are present at long time scales, and are enhanced during HUT.

Heart rate variability: are you using it properly? Standardisation checklist of procedures

BACKGROUND

Heart rate variability is used as an assessment method for cardiac autonomic modulation. Since the Task Force's publication on heart rate variability in 1996, the European Heart Rhythm Association Position Paper in 2015 and a recent publication in 2017, attention has been paid to recommendations on using heart rate variability analysis methods, as well as their applications in different physiological conditions and clinical studies. This analysis has proved to be useful as a complementary tool for clinical evaluation and to assess the effect of non-pharmacological therapeutic interventions, such as physical exercise programmes, on cardiac autonomic modulation.

OBJECTIVE

The aim of this article is to make recommendations and to develop a checklist of normalisation procedures regarding the use of heart rate variability data collection and analysis methodology, focusing on the cardiology area and cardiac rehabilitation.

METHODS

Based on previous heart rate variability publications, this paper provides a description of the most common shortcomings of using the analysis methods and considers recommendations and suggestions on how to minimise these occurrences by using a specific checklist.

CONCLUSIONS

This article includes recommendations and a checklist regarding the use of heart rate variability collection and analysis methods. This work could help improve reporting on clinical evaluation and therapeutic intervention results and consequently, disseminate heart rate variability knowledge.

Amplitude- and Fluctuation-Based Dispersion Entropy

Dispersion entropy (DispEn) is a recently introduced entropy metric to quantify the uncertainty of time series. It is fast and, so far, it has demonstrated very good performance in the characterisation of time series. It includes a mapping step, but the effect of different mappings has not been studied yet. Here, we investigate the effect of linear and nonlinear mapping approaches in DispEn. We also inspect the sensitivity of different parameters of DispEn to noise. Moreover, we develop fluctuation-based DispEn (FDispEn) as a measure to deal with only the fluctuations of time series. Furthermore, the original and fluctuation-based forbidden dispersion patterns are introduced to discriminate deterministic from stochastic time series. Finally, we compare the performance of DispEn, FDispEn, permutation entropy, sample entropy, and Lempel–Ziv complexity on two physiological datasets. The results show that DispEn is the most consistent technique to distinguish various dynamics of the biomedical signals. Due to their advantages over existing entropy methods, DispEn and FDispEn are expected to be broadly used for the characterization of a wide variety of real-world time series. The MATLAB codes used in this paper are freely available at http://dx.doi.org/10.7488/ds/2326.

Temporal Analysis of Cardiovascular and Respiratory Complexity by Multiscale Entropy Based on Symbolic Dynamics

The effect of an orthostatic stress on cardiovascular and respiratory complexity was investigated to detect impaired autonomic regulation in patients with vasovagal syncope (VVS). A total of 16 female patients and 12 age-matched healthy female subjects were enrolled in a passive 70° head-up tilt test. Also, 12 age-matched healthy male subjects were enrolled to study gender differences. Analysis was performed dynamically using various short-term (5 min) windows shifted by 1 min as well as by 20 min of orthostatic phase (OP) to evaluate local and global complexity. Complexity was determined over multiple time scales by the established method of refined composite multiscale entropy (RCMSE) and by a new proposed method of multiscale entropy based on symbolic dynamics (MSE-SD). Concerning heart rate variability (HRV) during OP, both methods revealed the highest complexity for female controls followed by lower complexity in male controls (p < 0.01) and by the lowest complexity in female patients (p < 0.01). For blood pressure variability (BPV), no gender differences in controls were shown by any method. However, MSE-SD demonstrated highly significantly increased BPV complexity in patients during OP (p < 0.01 on 4 time-scales after 7 min, p < 0.001 on 5 time-scales after 11 min) while RCMSE did not reveal considerable differences (p < 0.05 on 2 time scales after 7 min). Respiratory complexity was further increased in patients primary shown by MSE-SD. Findings indicated impaired autonomic regulation in VVS patients characterized by predominantly increased BPV complexity accompanied with decreased HRV complexity. In addition, results suggested extending the concept of complexity loss with disease.

The sensitivity of 38 heart rate variability measures to the addition of artifact in human and artificial 24-hr cardiac recordings

BACKGROUND

Artifact is common in cardiac RR interval data derived from 24-hr recordings and has a significant impact on heart rate variability (HRV) measures. However, the relative impact of progressively added artifact on a large group of commonly used HRV measures has not been assessed. This study compared the relative sensitivity of 38 commonly used HRV measures to artifact to determine which measures show the most change with increasing increments of artifact. A secondary aim was to ascertain whether short-term and long-term HRV measures, as groups, share similarities in their sensitivity to artifact.

METHODS

Up to 10% of artifact was added to 20 artificial RR (ARR) files and 20 human cardiac recordings, which had been assessed for artifact by a cardiac technician. The added artifact simulated deletion of RR intervals and insertion of individual short RR intervals. Thirty-eight HRV measures were calculated for each file. Regression analysis was used to rank the HRV measures according to their sensitivity to artifact as determined by the magnitude of slope.

RESULTS

RMSSD, SDANN, SDNN, RR triangular index and TINN, normalized power and relative power linear measures, and most nonlinear methods examined are most robust to artifact.

CONCLUSION

Short-term time domain HRV measures are more sensitive to added artifact than long-term measures. Absolute power frequency domain measures across all frequency bands are more sensitive than normalized and relative frequency domain measures. Most nonlinear HRV measures assessed were relatively robust to added artifact, with Poincare plot SD1 being most sensitive.

Baroreflex Coupling Assessed by Cross-Compression Entropy

Estimating interactions between physiological systems is an important challenge in modern biomedical research. Here, we explore a new concept for quantifying information common in two time series by cross-compressibility. Cross-compression entropy (CCE) exploits the ZIP data compression algorithm extended to bivariate data analysis. First, time series are transformed into symbol vectors. Symbols of the target time series are coded by the symbols of the source series. Uncoupled and linearly coupled surrogates were derived from cardiovascular recordings of 36 healthy controls obtained during rest to demonstrate suitability of this method for assessing physiological coupling. CCE at rest was compared to that of isometric handgrip exercise. Finally, spontaneous baroreflex interaction assessed by CCE_BRS was compared between 21 patients suffering from acute schizophrenia and 21 matched controls. The CCE_BRS of original time series was significantly higher than in uncoupled surrogates in 89% of the subjects and higher than in linearly coupled surrogates in 47% of the subjects. Handgrip exercise led to sympathetic activation and vagal inhibition accompanied by reduced baroreflex sensitivity. CCE_BRS decreased from 0.553 ± 0.030 at rest to 0.514 ± 0.035 during exercise (p < 0.001). In acute schizophrenia, heart rate, and blood pressure were elevated. Heart rate variability indicated a change of sympathovagal balance. The CCE_BRS of patients with schizophrenia was reduced compared to healthy controls (0.546 ± 0.042 vs. 0.507 ± 0.046, p < 0.01) and revealed a decrease of blood pressure influence on heart rate in patients with schizophrenia. Our results indicate that CCE is suitable for the investigation of linear and non-linear coupling in cardiovascular time series. CCE can quantify causal interactions in short, noisy and non-stationary physiological time series.

From Anxious to Reckless: A Control Systems Approach Unifies Prefrontal-Limbic Regulation Across the Spectrum of Threat Detection

Here we provide an integrative review of basic control circuits, and introduce techniques by which their regulation can be quantitatively measured using human neuroimaging. We illustrate the utility of the control systems approach using four human neuroimaging threat detection studies (N = 226), to which we applied circuit-wide analyses in order to identify the key mechanism underlying individual variation. In so doing, we build upon the canonical prefrontal-limbic control system to integrate circuit-wide influence from the inferior frontal gyrus (IFG). These were incorporated into a computational control systems model constrained by neuroanatomy and designed to replicate our experimental data. In this model, the IFG acts as an informational set point, gating signals between the primary prefrontal-limbic negative feedback loop and its cortical information-gathering loop. Along the cortical route, if the sensory cortex provides sufficient information to make a threat assessment, the signal passes to the ventromedial prefrontal cortex (vmPFC), whose threat-detection threshold subsequently modulates amygdala outputs. However, if signal outputs from the sensory cortex do not provide sufficient information during the first pass, the signal loops back to the sensory cortex, with each cycle providing increasingly fine-grained processing of sensory data. Simulations replicate IFG (chaotic) dynamics experimentally observed at both ends at the threat-detection spectrum. As such, they identify distinct types of IFG disconnection from the circuit, with associated clinical outcomes. If IFG thresholds are too high, the IFG and sensory cortex cycle for too long; in the meantime the coarse-grained (excitatory) pathway will dominate, biasing ambiguous stimuli as false positives. On the other hand, if cortical IFG thresholds are too low, the inhibitory pathway will suppress the amygdala without cycling back to the sensory cortex for much-needed fine-grained sensory cortical data, biasing ambiguous stimuli as false negatives. Thus, the control systems model provides a consistent mechanism for IFG regulation, capable of producing results consistent with our data for the full spectrum of threat-detection: from fearful to optimal to reckless. More generally, it illustrates how quantitative characterization of circuit dynamics can be used to unify a fundamental dimension across psychiatric affective symptoms, with implications for populations that range from anxiety disorders to addiction.

Symbolic time series analysis of electroencephalographic (EEG) epileptic seizure and brain dynamics with eye-open and eye-closed subjects during resting states

Objective

Epilepsy is a neuronal disorder for which the electrical discharge in the brain is synchronized, abnormal and excessive. To detect the epileptic seizures and to analyse brain activities during different mental states, various methods in non-linear dynamics have been proposed. This study is an attempt to quantify the complexity of control and epileptic subject with and without seizure as well as to distinguish eye-open (EO) and eye-closed (EC) conditions using threshold-based symbolic entropy.

Methods

The threshold-dependent symbolic entropy was applied to distinguish the healthy and epileptic subjects with seizure and seizure-free intervals (i.e. interictal and ictal) as well as to distinguish EO and EC conditions. The original time series data was converted into symbol sequences using quantization level, and word series of symbol sequences was generated using a word length of three or more. Then, normalized corrected Shannon entropy (NCSE) was computed to quantify the complexity. The NCSE values were not following the normal distribution, and the non-parametric Mann–Whitney–Wilcoxon (MWW) test was used to find significant differences among various groups at 0.05 significance level. The values of NCSE were presented in a form of topographic maps to show significant brain regions during EC and EO conditions. The results of the study were compared to those of the multiscale entropy (MSE).

Results

The results indicated that the dynamics of healthy subjects are more complex compared to epileptic subjects (during seizure and seizure-free intervals) in both EO and EC conditions. The comparison of the dynamics of epileptic subjects revealed that seizure-free intervals are more complex than seizure intervals. The dynamics of healthy subjects during EO conditions are more complex compared to those during EC conditions. Further, the results clearly demonstrated that threshold-dependent symbolic entropy outperform MSE in distinguishing different physiological and pathological conditions.

Conclusion

The threshold symbolic entropy has provided improved accuracy in quantifying the dynamics of healthy and epileptic subjects during EC an EO conditions for each electrode compared to the MSE.

HRVanalysis: A Free Software for Analyzing Cardiac Autonomic Activity

Since the pioneering studies of the 1960s, heart rate variability (HRV) has become an increasingly used non-invasive tool for examining cardiac autonomic functions and dysfunctions in various populations and conditions. Many calculation methods have been developed to address these issues, each with their strengths and weaknesses. Although, its interpretation may remain difficult, this technique provides, from a non-invasive approach, reliable physiological information that was previously inaccessible, in many fields including death and health prediction, training and overtraining, cardiac and respiratory rehabilitation, sleep-disordered breathing, large cohort follow-ups, children's autonomic status, anesthesia, or neurophysiological studies. In this context, we developed HRVanalysis, a software to analyse HRV, used and improved for over 20 years and, thus, designed to meet laboratory requirements. The main strength of HRVanalysis is its wide application scope. In addition to standard analysis over short and long periods of RR intervals, the software allows time-frequency analysis using wavelet transform as well as analysis of autonomic nervous system status on surrounding scored events and on preselected labeled areas. Moreover, the interface is designed for easy study of large cohorts, including batch mode signal processing to avoid running repetitive operations. Results are displayed as figures or saved in TXT files directly employable in statistical softwares. Recordings can arise from RR or EKG files of different types such as cardiofrequencemeters, holters EKG, polygraphs, and data acquisition systems. HRVanalysis can be downloaded freely from the Web page at: https://anslabtools.univ-st-etienne.frHRVanalysis is meticulously maintained and developed for in-house laboratory use. In this article, after a brief description of the context, we present an overall view of HRV analysis and we describe the methodological approach of the different techniques provided by the software.

Orthostatic stress causes immediately increased blood pressure variability in women with vasovagal syncope

The cardiovascular and respiratory autonomic nervous regulation has been studied mainly by hemodynamic responses during different physical stressors. In this study, dynamics of autonomic response to an orthostatic challenge was investigated by hemodynamic variables and by diverse linear and nonlinear indices calculated from time series of beat-to-beat intervals (BBI), respiratory cycle duration (RESP), systolic (SYS) and diastolic (DIA) blood pressure. This study included 16 young female patients (SYN) with vasovagal syncope and 12 age-matched female controls (CON). The subjects were enrolled in a head-up tilt (HUT) test, breathing normally, including 5min of baseline (BL, supine position) and 18min of 70° orthostatic phase (OP). To increase the time resolution of the analysis the time series were segmented in five-minute overlapping windows with a shift of 1min. Hemodynamic parameters did not show any statistical differences between SYN and CON. Time domain linear analysis revealed increased respiratory frequency and increased blood pressure variability (BPV) in patients during OP meaning increased sympathetic activity and vagal withdrawal. Frequency domain analysis confirmed a predominance of sympathetic tone by steadily increased values of low over high frequency power in BBI and of low frequency power in SYS and DIA in patients during OP. The nonlinear analysis by symbolic dynamics seemed to be highly suitable for differentiation of SYN and CON in the early beginning of OP, i.e., 5min after tilt-up. In particular the index SYS_plvar3 showed less patterns of low variability in patients reflecting a steadily increase in both BPV and sympathetic activity. The proposed dynamical analysis could lead to a better understanding of the temporal underlying mechanisms in healthy subjects and patients under orthostatic stress.

Ordinal symbolic analysis and its application to biomedical recordings

Ordinal symbolic analysis opens an interesting and powerful perspective on time-series analysis. Here, we review this relatively new approach and highlight its relation to symbolic dynamics and representations. Our exposition reaches from the general ideas up to recent developments, with special emphasis on its applications to biomedical recordings. The latter will be illustrated with epilepsy data.

Heart rate variability: a tool to explore the sleeping brain?

Sleep is divided into two main sleep stages: (1) non-rapid eye movement sleep (non-REMS), characterized among others by reduced global brain activity; and (2) rapid eye movement sleep (REMS), characterized by global brain activity similar to that of wakefulness. Results of heart rate variability (HRV) analysis, which is widely used to explore autonomic modulation, have revealed higher parasympathetic tone during normal non-REMS and a shift toward sympathetic predominance during normal REMS. Moreover, HRV analysis combined with brain imaging has identified close connectivity between autonomic cardiac modulation and activity in brain areas such as the amygdala and insular cortex during REMS, but no connectivity between brain and cardiac activity during non-REMS. There is also some evidence for an association between HRV and dream intensity and emotionality. Following some technical considerations, this review addresses how brain activity during sleep contributes to changes in autonomic cardiac activity, organized into three parts: (1) the knowledge on autonomic cardiac control, (2) differences in brain and autonomic activity between non-REMS and REMS, and (3) the potential of HRV analysis to explore the sleeping brain, and the implications for psychiatric disorders.

Network connectivity modulates power spectrum scale invariance

Measures of complexity are sensitive in detecting disease, which has made them attractive candidates for diagnostic biomarkers; one complexity measure that has shown promise in fMRI is power spectrum scale invariance (PSSI). Even if scale-free features of neuroimaging turn out to be diagnostically useful, however, their underlying neurobiological basis is poorly understood. Using modeling and simulations of a schematic prefrontal-limbic meso-circuit, with excitatory and inhibitory networks of nodes, we present here a framework for how network density within a control system can affect the complexity of signal outputs. Our model demonstrates that scale-free behavior, similar to that observed in fMRI PSSI data, can be obtained for sufficiently large networks in a context as simple as a linear stochastic system of differential equations, although the scale-free range improves when introducing more realistic, nonlinear behavior in the system. PSSI values (reflective of complexity) vary as a function of both input type (excitatory, inhibitory) and input density (mean number of long-range connections, or strength), independent of their node-specific geometric distribution. Signals show pink noise (1/f) behavior when excitatory and inhibitory influences are balanced. As excitatory inputs are increased and decreased, signals shift towards white and brown noise, respectively. As inhibitory inputs are increased and decreased, signals shift towards brown and white noise, respectively. The results hold qualitatively at the hemodynamic scale, which we modeled by introducing a neurovascular component. Comparing hemodynamic simulation results to fMRI PSSI results from 96 individuals across a wide spectrum of anxiety-levels, we show how our model can generate concrete and testable hypotheses for understanding how connectivity affects regulation of meso-circuits in the brain.

[Heart rate variability during sleep in children with multiple disabilities]

OBJECTIVES

To study heart rate variability during sleep in children with multiple disabilities in order to observe the behavior of the autonomic nervous system.

METHODS

The R-R interval variability of 4 to 12 years old children was recorded with a heart rate monitor during one night. Children with multiple disabilities (G1) and healthy children (G2) were compared in time, frequency, and non-linear domains.

RESULTS

Temporal (P<0.01), frequency (P<0.05), and non-linear (P<0.01) variables in the G1 were lower than in the G2 group.

DISCUSSION

The time and frequency analysis confirmed the predominance of the sympathetic nervous system during sleep in G1 children. Reduced non-linear variables can explain sleep with less informations (correlation dimension: 2.203 ± 1.239 versus 3.842 ± 0.378; P<0.001), less complex cycles of sleep (approximate entropy: 1.153 ± 0.200 versus 1.365 ± 0.099; P<0.01), and a reduced correlation in short-term variability (SD1: 42.37 ± 19.0 (ms) versus 73.44 ± 25.3; P<0.01). The fractal structure of the recordings was not affected (P>0.05). The diseases encountered are probably the reason for these findings, but the variety of disorders and medications of the children with multiple disabilities needs to be studied with a larger and more varied sample.

CONCLUSION

Sympathetic predominance during sleep in children with multiple disabilities is associated with a decrease in adaptive abilities of these children's autonomic nervous system.

Heart rate variability in normal and pathological sleep

Sleep is a physiological process involving different biological systems, from molecular to organ level; its integrity is essential for maintaining health and homeostasis in human beings. Although in the past sleep has been considered a state of quiet, experimental and clinical evidences suggest a noteworthy activation of different biological systems during sleep. A key role is played by the autonomic nervous system (ANS), whose modulation regulates cardiovascular functions during sleep onset and different sleep stages. Therefore, an interest on the evaluation of autonomic cardiovascular control in health and disease is growing by means of linear and non-linear heart rate variability (HRV) analyses. The application of classical tools for ANS analysis, such as HRV during physiological sleep, showed that the rapid eye movement (REM) stage is characterized by a likely sympathetic predominance associated with a vagal withdrawal, while the opposite trend is observed during non-REM sleep. More recently, the use of non-linear tools, such as entropy-derived indices, have provided new insight on the cardiac autonomic regulation, revealing for instance changes in the cardiovascular complexity during REM sleep, supporting the hypothesis of a reduced capability of the cardiovascular system to deal with stress challenges. Interestingly, different HRV tools have been applied to characterize autonomic cardiac control in different pathological conditions, from neurological sleep disorders to sleep disordered breathing (SDB). In summary, linear and non-linear analysis of HRV are reliable approaches to assess changes of autonomic cardiac modulation during sleep both in health and diseases. The use of these tools could provide important information of clinical and prognostic relevance.

Auditory stimulation and cardiac autonomic regulation

Previous studies have already demonstrated that auditory stimulation with music influences the cardiovascular system. In this study, we described the relationship between musical auditory stimulation and heart rate variability. Searches were performed with the Medline, SciELO, Lilacs and Cochrane databases using the following keywords: "auditory stimulation", "autonomic nervous system", "music" and "heart rate variability". The selected studies indicated that there is a strong correlation between noise intensity and vagal-sympathetic balance. Additionally, it was reported that music therapy improved heart rate variability in anthracycline-treated breast cancer patients. It was hypothesized that dopamine release in the striatal system induced by pleasurable songs is involved in cardiac autonomic regulation. Musical auditory stimulation influences heart rate variability through a neural mechanism that is not well understood. Further studies are necessary to develop new therapies to treat cardiovascular disorders.

Quantification of compensatory processes of postnatal hypoxia in newborn piglets applying short-term nonlinear dynamics analysis

Background

Newborn mammals suffering from moderate hypoxia during or after birth are able to compensate a transitory lack of oxygen by adapting their vital functions. Exposure to hypoxia leads to an increase in the sympathetic tone causing cardio-respiratory response, peripheral vasoconstriction and vasodilatation in privileged organs like the heart and brain. However, there is only limited information available about the time and intensity changes of the underlying complex processes controlled by the autonomic nervous system.

Methods

In this study an animal model involving seven piglets was used to examine an induced state of circulatory redistribution caused by moderate oxygen deficit. In addition to the main focus on the complex dynamics occurring during sustained normocapnic hypoxia, the development of autonomic regulation after induced reoxygenation had been analysed. For this purpose, we first introduced a new algorithm to prove stationary conditions in short-term time series. Then we investigated a multitude of indices from heart rate and blood pressure variability and from bivariate interactions, also analysing respiration signals, to quantify the complexity of vegetative oscillations influenced by hypoxia.

Results

The results demonstrated that normocapnic hypoxia causes an initial increase in cardiovascular complexity and variability, which decreases during moderate hypoxia lasting one hour (p < 0.004). After reoxygenation, cardiovascular complexity parameters returned to pre-hypoxic values (p < 0.003), however not respiratory-related complexity parameters.

Conclusions

In conclusion, indices from linear and nonlinear dynamics reflect considerable temporal changes of complexity in autonomous cardio-respiratory regulation due to normocapnic hypoxia shortly after birth. These findings might be suitable for non-invasive clinical monitoring of hypoxia-induced changes of autonomic regulation in newborn humans.

Cardiovascular regulation during sleep quantified by symbolic coupling traces

Sleep is a complex regulated process with short periods of wakefulness and different sleep stages. These sleep stages modulate autonomous functions such as blood pressure and heart rate. The method of symbolic coupling traces (SCT) is used to analyze and quantify time-delayed coupling of these measurements during different sleep stages. The symbolic coupling traces, defined as the symmetric and diametric traces of the bivariate word distribution matrix, allow the quantification of time-delayed coupling. In this paper, the method is applied to heart rate and systolic blood pressure time series during different sleep stages for healthy controls as well as for normotensive and hypertensive patients with sleep apneas. Using the SCT, significant different cardiovascular mechanisms not only between the deep sleep and the other sleep stages but also between healthy subjects and patients can be revealed. The SCT method is applied to model systems, compared with established methods, such as cross correlation, mutual information, and cross recurrence analysis and demonstrates its advantages especially for nonstationary physiological data. As a result, SCT proves to be more specific in detecting delays of directional interactions than standard coupling analysis methods and yields additional information which cannot be measured by standard parameters of heart rate and blood pressure variability. The proposed method may help to indicate the pathological changes in cardiovascular regulation and also the effects of continuous positive airway pressure therapy on the cardiovascular system.

Spinal tumor necrosis factor alpha neutralization reduces peripheral inflammation and hyperalgesia and suppresses autonomic responses in experimental arthritis: a role for spinal tumor necrosis factor alpha during induction and maintenance of peripheral inflammation

OBJECTIVE

In addition to the sensitization of pain fibers in inflamed tissues, the increased excitability of the spinal cord is an important mechanism of inflammatory pain. Furthermore, spinal neuronal excitability has been suggested to play a role in modulating peripheral inflammation. This study was undertaken to test the hypothesis that spinal actions of the proinflammatory cytokine tumor necrosis factor alpha (TNFalpha) add significantly to both hyperalgesia and maintenance of peripheral inflammation.

METHODS

Rats with antigen-induced arthritis (AIA) were treated intrathecally with the TNFalpha-neutralizing compound etanercept continuously during the complete time course of AIA, which was 3 days for the acute phase and 21 days for the chronic phase. During this time, inflammation and pain-related behavior were monitored. Since a role for autonomic control of inflammation was proposed, measures from heart rate time series were obtained in the acute phase. Findings were compared with those in vehicle-treated animals and in animals receiving etanercept intraperitoneally.

RESULTS

Spinally administered etanercept acutely reduced pain-related behavior, attenuated both the development and the maintenance of inflammation, and was superior to systemic administration. Parameters indicating autonomic modulation showed a shift toward a sympathetically dominated state in vehicle-treated animals, which was prevented by intrathecal etanercept.

CONCLUSION

Our findings indicate that spinal TNFalpha plays an important role in both pain signaling and modulation of peripheral inflammation. Thus, neutralizing this cytokine at the spinal site not only represents a putative therapeutic option for different pain syndromes, but may be directly used to attenuate peripheral inflammation.

Power spectrum scale invariance quantifies limbic dysregulation in trait anxious adults using fMRI: adapting methods optimized for characterizing autonomic dysregulation to neural dynamic time series

In a well-regulated control system, excitatory and inhibitory components work closely together with minimum lag; in response to inputs of finite duration, outputs should show rapid rise and, following the input's termination, immediate return to baseline. The efficiency of this response can be quantified using the power spectrum density's scaling parameter beta, a measure of self-similarity, applied to the first derivative of the raw signal. In this study, we adapted power spectrum density methods, previously used to quantify autonomic dysregulation (heart rate variability), to neural time series obtained via functional MRI. The negative feedback loop we investigated was the limbic system, using affect-valent faces as stimuli. We hypothesized that trait anxiety would be related to efficiency of regulation of limbic responses, as quantified by power-law scaling of fMRI time series. Our results supported this hypothesis, showing moderate to strong correlations of trait anxiety and beta (r=0.45-0.54) for the amygdala, orbitofrontal cortex, hippocampus, superior temporal gyrus, posterior insula, and anterior cingulate. Strong anticorrelations were also found between the amygdala's beta and wake heart rate variability (r=-0.61), suggesting a robust relationship between dysregulated limbic outputs and their autonomic consequences.

Depression of cardiovascular autonomic function is more pronounced after mitral valve surgery: evidence for direct trauma

The analysis of baroreflex sensitivity (BRS) and heart rate variability (HRV) leads to additional insights into patients' prognosis after cardiovascular events. The following study was performed to assess the differences in the post-operative recovery of autonomic regulation after mitral valve (MV) and aortic valve (AV) surgery with a heart-lung machine. Among the 43 consecutive male patients enrolled in a prospective study, 26 underwent isolated AV surgery and 17 isolated MV surgery. Blood pressure as well as ECG signals were recorded the day before, 24 hours after and one week after surgery. BRS was calculated according to the dual sequence method, and HRV was calculated using standard linear as well as nonlinear parameters. There were no major differences between the two groups in the pre-operative values. At 24 hours a comparable depression of HRV and BRS in both groups was observed, while at 7 days there was partial recovery in AV patients, which was absent in MV patients: p(AV versus MV)<0.001. While the response of the autonomic system to surgery is similar in AV and MV patients, there is obviously a decreased ability to recover in MV patients, probably attributed to traumatic lesions of the autonomic nervous system by opening the atria. Ongoing research is required for further clarification of the pathophysiology of this phenomenon and to establish strategies to restore autonomic function.

Limbic dysregulation is associated with lowered heart rate variability and increased trait anxiety in healthy adults

OBJECTIVES

We tested whether dynamic interaction between limbic regions supports a control systems model of excitatory and inhibitory components of a negative feedback loop, and whether dysregulation of those dynamics might correlate with trait differences in anxiety and their cardiac characteristics among healthy adults.

EXPERIMENTAL DESIGN

Sixty-five subjects received fMRI scans while passively viewing angry, fearful, happy, and neutral facial stimuli. Subjects also completed a trait anxiety inventory, and were monitored using ambulatory wake ECG. The ECG data were analyzed for heart rate variability, a measure of autonomic regulation. The fMRI data were analyzed with respect to six limbic regions (bilateral amygdala, bilateral hippocampus, Brodmann Areas 9, 45) using limbic time-series cross-correlations, maximum BOLD amplitude, and BOLD amplitude at each point in the time-series.

PRINCIPAL OBSERVATIONS

Diminished coupling between limbic time-series in response to the neutral, fearful, and happy faces was associated with greater trait anxiety, greater sympathetic activation, and lowered heart rate variability. Individuals with greater levels of trait anxiety showed delayed activation of Brodmann Area 45 in response to the fearful and happy faces, and lowered Brodmann Area 45 activation with prolonged left amygdala activation in response to the neutral faces.

CONCLUSIONS

The dynamics support limbic regulation as a control system, in which dysregulation, as assessed by diminished coupling between limbic time-series, is associated with increased trait anxiety and excitatory autonomic outputs. Trait-anxious individuals showed delayed inhibitory activation in response to overt-affect stimuli and diminished inhibitory activation with delayed extinction of excitatory activation in response to ambiguous-affect stimuli.

Weaning From Mechanical Ventilation: A Retrospective Analysis Leading to a Multimodal Perspective

Practitioners' decision for mechanical aid discontinuation is a challenging task that involves a complete knowledge of a great number of clinical parameters, as well as its evolution in time. Recently, an increasing interest on respiratory pattern variability as an extubation readiness indicator has appeared. Reliable assessment of this variability involves a set of signal processing and pattern recognition techniques. This paper presents a suitability analysis of different methods used for breathing pattern complexity assessment. The contribution of this analysis is threefold: 1) to serve as a review of the state of the art on the so-called weaning problem from a signal processing point of view; 2) to provide insight into the applied processing techniques and how they fit into the problem; 3) to propose additional methods and further processing in order to improve breathing pattern regularity assessment and weaning readiness decision. Results on experimental data show that sample entropy outperforms other complexity assessment methods and that multidimensional classification does improve weaning prediction. However, the obtained performance may be objectionable for real clinical practice, a fact that paves the way for a multimodal signal processing framework, including additional high-quality signals and more reliable statistical methods.

Correlation dimension analysis of heart rate variability in patients with dilated cardiomyopathy

A correlation dimension analysis of heart rate variability (HRV) was applied to a group of 55 patients with dilated cardiomyopathy (DCM) and 55 healthy subjects as controls. The 24-h RR time series for each subject was divided into segments of 10,000 beats to determine the correlation dimension (CD) per segment. A study of the influence of the time delay (lag) in the calculation of CD was performed. Good discrimination between both groups (p<0.005) was obtained with lag values of 5 or greater. CD values of DCM patients (8.4+/-1.9) were significantly lower than CD values for controls (9.5+/-1.9). An analysis of CD values of HRV showed that for healthy people, CD night values (10.6+/-1.8) were significant greater than CD day values (9.2+/-1.9), revealing a circadian rhythm. In DCM patients, this circadian rhythm was lost and there were no differences between CD values in day (8.8+/-2.4) and night (8.9+/-2.1).

Stochastic model for heart-rate fluctuations

A normal human heart rate shows complex fluctuations in time, which is natural, because the heart rate is controlled by a large number of different feedback control loops. These unpredictable fluctuations have been shown to display fractal dynamics, long-term correlations, and 1/f noise. These characterizations are statistical and they have been widely studied and used, but much less is known about the detailed time evolution (dynamics) of the heart-rate control mechanism. Here we show that a simple one-dimensional Langevin-type stochastic difference equation can accurately model the heart-rate fluctuations in a time scale from minutes to hours. The model consists of a deterministic nonlinear part and a stochastic part typical to Gaussian noise, and both parts can be directly determined from the measured heart-rate data. Studies of 27 healthy subjects reveal that in most cases, the deterministic part has a form typically seen in bistable systems: there are two stable fixed points and one unstable one.

Nonlinear methods of biosignal analysis in assessing terbutaline-induced heart rate and blood pressure changes

The aim of this study was to characterize how different nonlinear methods characterize heart rate and blood pressure dynamics in healthy subjects at rest. The randomized, placebo-controlled crossover study with intravenous terbutaline was designed to induce four different stationary states of cardiovascular regulation system. The R-R interval, systolic arterial blood pressure, and heart rate time series were analyzed with a set of methods including approximate entropy, sample entropy, Lempel-Ziv entropy, symbol dynamic entropy, cross-entropy, correlation dimension, fractal dimensions, and stationarity test. Results indicate that R-R interval and systolic arterial pressure subsystems are mutually connected but have different dynamic properties. In the drug-free state the subsystems share many common features. When the strength of the baroreflex feedback loop is modified with terbutaline, R-R interval and systolic blood pressure lose mutual synchrony and drift toward their inherent state of operation. In this state the R-R interval system is rather complex and irregular, but the blood pressure system is much simpler than in the drug-free state.

Entropies of short binary sequences in heart period dynamics

Dynamic aspects of R-R intervals have often been analyzed by means of linear and nonlinear measures. The goal of this study was to analyze binary sequences, in which only the dynamic information is retained, by means of two different aspects of regularity. R-R interval sequences derived from 24-h electrocardiogram (ECG) recordings of 118 healthy subjects were converted to symbolic binary sequences that coded the beat-to-beat increase or decrease in the R-R interval. Shannon entropy was used to quantify the occurrence of short binary patterns (length N = 5) in binary sequences derived from 10-min intervals. The regularity of the short binary patterns was analyzed on the basis of approximate entropy (ApEn). ApEn had a linear dependence on mean R-R interval length, with increasing irregularity occurring at longer R-R interval length. Shannon entropy of the same sequences showed that the increase in irregularity is accompanied by a decrease in occurrence of some patterns. Taken together, these data indicate that irregular binary patterns are more probable when the mean R-R interval increases. The use of surrogate data confirmed a nonlinear component in the binary sequence. Analysis of two consecutive 24-h ECG recordings for each subject demonstrated good intraindividual reproducibility of the results. In conclusion, quantification of binary sequences derived from ECG recordings reveals properties that cannot be found using the full information of R-R interval sequences.

The frequency of occurrence of anti-cardiac receptor autoantibodies and their correlation with clinical manifestation in patients with hypertrophic cardiomyopathy

The purpose of this study was to investigate the frequency of occurrence of autoantibodies against G-protein coupled cardiovascular receptors and their relation to the clinical manifestation of hypertrophic cardiomyopathy (HCM). Autoantibodies against beta1-receptors, Muscarin-2-receptors, Angiotensin-II-receptor subtype 1 and alpha1-receptors were determined with ELISA in 52 patients with HCM (37 male, 15 female, mean age 55 +/- 15 years) and 40 healthy, age and sex matched controls. The clinical characterization of the HCM-patients included ECG, 24-h Holter, and echocardiography. The results showed that there is no significant difference in the frequency of a single autoantibody between HCM-patients and controls. However, if the number of patients who have autoantibodies against beta1-receptors and/or Muscarin-2-receptors were counted together, there are significantly more autoantibodies in HCM compared to controls (11 vs. 2, p = 0.035). Analysis of clinical data from this pooled group of patients showed that in patients with autoantibodies, heart rate variability (HRV), ultra low frequency (ULF) and very low frequency (VLF) were decreased (HRV by 20%, ULF by 50%, and VLF by 46%, p < 0.008) whereas the QTc-interval was increased by 8% (p < 0.02 each). The ratio of septal to posterior wall thickness was increased by 23% (p = 0.05), and the preejection period was prolonged by 46% in patients with autoantibodies (p < 0.001). These results suggest that the existence of these autoantibodies could be associated with an advanced stage or a severe manifestation of HCM.

Familial and genetic influences on heart rate variability

The authors tested the hypothesis of a genetic influence on heart rate variability (HRV). This genetic influence was assessed in 62, twin pairs (30 monozygotic, 32 dizygotic). From all twins, long-term electrocardiographic records were obtained, edited, and analyzed. Heart rate variability analysis was performed on the basis of parameters from time-domain, frequency-domain, and nonlinear dynamics. First, the parameter distances between the two twins of a pair and between one of the two and a third randomly selected person of another age-matched twin pair (ST1) were compared. Second, the parameter distances between the two twins and the averaged parameter distances of these two twins to all other age-matched persons (ST2) were compared. Finally, the averaged differences in parameter values between monozygotic and dizygotic age-matched twin pairs were compared. For statistical analysis, the nonparametric Wilcoxon's matched-pair signed rank test and parametric t-test for paired samples were used. Twin pairs show a significant lower difference in parameter values than other randomly selected and age-matched couples (P < .001 in ST1 and ST2). This reflects a considerable familial influence. Most parameters of the time-domain, none of the frequency-domain, and half of the nonlinear dynamics show significant differences between twin pairs and nontwin pairs. As a result of the comparison between monozygotic and dizygotic twin pairs, a significant lower parameter difference in the monozygotic pairs (P < .05) is found. These results suggest that there is a genetic component in heart rate generation and HRV, in addition to family environmental influences. Analysis of HRV might become a useful method in phenotyping severe genetic changes in cardiovascular diseases.