Mechanisms of causal interaction between short-term RR interval and systolic arterial pressure oscillations during orthostatic challenge

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.

Concomitant evaluation of cardiovascular and cerebrovascular controls via Geweke spectral causality to assess the propensity to postural syncope

The evaluation of propensity to postural syncope necessitates the concomitant characterization of the cardiovascular and cerebrovascular controls and a method capable of disentangling closed loop relationships and decomposing causal links in the frequency domain. We applied Geweke spectral causality (GSC) to assess cardiovascular control from heart period and systolic arterial pressure variability and cerebrovascular regulation from mean arterial pressure and mean cerebral blood velocity variability in 13 control subjects and 13 individuals prone to develop orthostatic syncope. Analysis was made at rest in supine position and during head-up tilt at 60°, well before observing presyncope signs. Two different linear model structures were compared, namely bivariate autoregressive and bivariate dynamic adjustment classes. We found that (i) GSC markers did not depend on the model structure; (ii) the concomitant assessment of cardiovascular and cerebrovascular controls was useful for a deeper comprehension of postural disturbances; (iii) orthostatic syncope appeared to be favored by the loss of a coordinated behavior between the baroreflex feedback and mechanical feedforward pathway in the frequency band typical of the baroreflex functioning during the postural challenge, and by a weak cerebral autoregulation as revealed by the increased strength of the pressure-to-flow link in the respiratory band. GSC applied to spontaneous cardiovascular and cerebrovascular oscillations is a promising tool for describing and monitoring disturbances associated with posture modification.

Differentiating patients with obstructive sleep apnea from healthy controls based on heart rate-blood pressure coupling quantified by entropy-based indices

We introduce an entropy-based classification method for pairs of sequences (ECPS) for quantifying mutual dependencies in heart rate and beat-to-beat blood pressure recordings. The purpose of the method is to build a classifier for data in which each item consists of two intertwined data series taken for each subject. The method is based on ordinal patterns and uses entropy-like indices. Machine learning is used to select a subset of indices most suitable for our classification problem in order to build an optimal yet simple model for distinguishing between patients suffering from obstructive sleep apnea and a control group.

On the Validity of Single Regression Strategy for Granger Causality Assessment in Cardiovascular and Cardiorespiratory Control Studies

Granger causality (GC) analysis is based on the comparison between prediction error variances computed over the full and restricted models after identifying the coefficients of appropriate vector regressions. GC markers can be computed via a double regression (DR) approach identifying two separate, independent models and a single regression (SR) strategy optimizing the description of the dynamics of the target over the full model and, then, reusing some parts of it in the restricted model. The present study compares the SR and DR strategies over heart period (HP), systolic arterial pressure (SAP) and respiration (R) beat-to-beat series collected during a graded orthostatic challenge induced by head-up tilt in 17 healthy individuals (age: 21-36 yrs; median: 29 yrs; 9 females and 8 males). We found that the DR approach was more powerful than the SR one in detecting the expected stronger involvement of the baroreflex during the challenge, while the expected weaker cardiorespiratory coupling was identified by both SR and DR strategies. The less powerful ability of the SR approach was the result of the greater variance of GC markers compared to the DR strategy. We conclude that, contrary to the suggestions present in literature, the SR approach is not necessarily associated with a smaller dispersion of GC markers. Moreover, we suggest that additional factors, such as the strategy utilized to build embedding spaces and metric utilized to compare prediction error variances, might play an important role in differentiating SR and DR approaches.

Spectral decomposition of cerebrovascular and cardiovascular interactions in patients prone to postural syncope and healthy controls

We present a framework for the linear parametric analysis of pairwise interactions in bivariate time series in the time and frequency domains, which allows the evaluation of total, causal and instantaneous interactions and connects time- and frequency-domain measures. The framework is applied to physiological time series to investigate the cerebrovascular regulation from the variability of mean cerebral blood flow velocity (CBFV) and mean arterial pressure (MAP), and the cardiovascular regulation from the variability of heart period (HP) and systolic arterial pressure (SAP). We analyze time series acquired at rest and during the early and late phase of head-up tilt in subjects developing orthostatic syncope in response to prolonged postural stress, and in healthy controls. The spectral measures of total, causal and instantaneous coupling between HP and SAP, and between MAP and CBFV, are averaged in the low-frequency band of the spectrum to focus on specific rhythms, and over all frequencies to get time-domain measures. The analysis of cardiovascular interactions indicates that postural stress induces baroreflex involvement, and its prolongation induces baroreflex dysregulation in syncope subjects. The analysis of cerebrovascular interactions indicates that the postural stress enhances the total coupling between MAP and CBFV, and challenges cerebral autoregulation in syncope subjects, while the strong sympathetic activation elicited by prolonged postural stress in healthy controls may determine an increased coupling from CBFV to MAP during late tilt. These results document that the combination of time-domain and spectral measures allows us to obtain an integrated view of cardiovascular and cerebrovascular regulation in healthy and diseased subjects.

Information decomposition in the frequency domain: a new framework to study cardiovascular and cardiorespiratory oscillations

While cross-spectral and information-theoretic approaches are widely used for the multivariate analysis of physiological time series, their combined utilization is far less developed in the literature. This study introduces a framework for the spectral decomposition of multivariate information measures, which provides frequency-specific quantifications of the information shared between a target and two source time series and of its expansion into amounts related to how the sources contribute to the target dynamics with unique, redundant and synergistic information. The framework is illustrated in simulations of linearly interacting stochastic processes, showing how it allows us to retrieve amounts of information shared by the processes within specific frequency bands which are otherwise not detectable by time-domain information measures, as well as coupling features which are not detectable by spectral measures. Then, it is applied to the time series of heart period, systolic and diastolic arterial pressure and respiration variability measured in healthy subjects monitored in the resting supine position and during head-up tilt. We show that the spectral measures of unique, redundant and synergistic information shared by these variability series, integrated within specific frequency bands of physiological interest and reflect the mechanisms of short-term regulation of cardiovascular and cardiorespiratory oscillations and their alterations induced by the postural stress. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Efficient syncope prediction from resting state clinical data using wavelet bispectrum and multilayer perceptron neural network

Neurally mediated syncope (NMS) is the most common type of syncope, and head up tilt test (HUTT) is, so far, the most appropriate tool to identify NMS. In this work, an effort to predict the NMS before performing the HUTT is attempted. To achieve this, the heart rate variability (HRV) at rest and during the first minutes of tilting position during HUTT was analyzed using both time and frequency domains. Various features from HRV regularity and complexity, along with wavelet higher-order spectrum (WHOS) analysis in low-frequency (LF) and high-frequency (HF) bands were examined. The experimental results from 26 patients with history of NMS have shown that at rest, a time domain entropy measure and WHOS-based features in LF band exhibit significant differences between positive and negative HUTT as well as among 10 healthy subjects and NMS patients. The best performance of multilayer perceptron neural network (MPNN) was achieved by using an input vector consisted of WHOS-based HRV features in the LF zone and systolic blood pressure from the resting period, yielding an accuracy of 89.7%, assessed by 5-fold cross-validation. The promising results presented here pave the way for an early prediction of the HUTT outcome from resting state, contributing to the identification of patients at higher risk NMS. The HRV analysis along with systolic blood pressure at rest predict NMS using a multilayer perceptron neural network.

Multicomponent exercise training in cardiovascular complexity in prefrail older adults: a randomized blinded clinical pilot study

The aim of this study was to investigate the effects of multicomponent training on baroreflex sensitivity (BRS) and heart rate (HR) complexity of prefrail older adults. Twenty-one prefrail community-dwelling older adults were randomized and divided into multicomponent training intervention group (MulTI) and control group (CG). MulTI performed multicomponent exercise training over 16 weeks and CG was oriented to follow their own daily activities. The RR interval (RRi) and blood pressure (BP) series were recorded for 15 min in supine and 15 min in orthostatic positions, and calculation of BRS (phase, coherence, and gain) and HR complexity (sample entropy) were performed. A linear mixed model was applied for group, assessments, and their interaction effects in supine position. The same test was used to assess the active postural maneuver and it was applied separately to each group considering assessments (baseline and post-intervention) and positions (supine and orthostatic). The significance level established was 5%. Cardiovascular control was impaired in prefrail older adults in supine position. Significant interactions were not observed between groups or assessments in terms of cardiovascular parameters. A 16-week multicomponent exercise training did not improve HR complexity or BRS in supine rest or in active postural maneuver in prefrail older adults.

Parasympathetic-Sympathetic Causal Interactions Assessed by Time-Varying Multivariate Autoregressive Modeling of Electrodermal Activity and Heart-Rate-Variability

OBJECTIVE

Most of the bodily functions are regulated by multiple interactions between the parasympathetic (PNS) and sympathetic (SNS) nervous system. In this study, we propose a novel framework to quantify the causal flow of information between PNS and SNS through the analysis of heart rate variability (HRV) and electrodermal activity (EDA) signals.

METHODS

Our method is based on a time-varying (TV) multivariate autoregressive model of EDA and HRV time-series and incorporates physiologically inspired assumptions by estimating the Directed Coherence in a specific frequency range. The statistical significance of the observed interactions is assessed by a bootstrap procedure purposely developed to infer causalities in the presence of both TV model coefficients and TV model residuals (i.e., heteroskedasticity). We tested our method on two different experiments designed to trigger a sympathetic response, i.e., a hand-grip task (HG) and a mental-computation task (MC).

RESULTS

Our results show a parasympathetic driven interaction in the resting state, which is consistent across different studies. The onset of the stressful stimulation triggers a cascade of events characterized by the presence or absence of the PNS-SNS interaction and changes in the directionality. Despite similarities between the results related to the two tasks, we reveal differences in the dynamics of the PNS-SNS interaction, which might reflect different regulatory mechanisms associated with different stressors.

CONCLUSION

We estimate causal coupling between PNS and SNS through MVAR modeling of EDA and HRV time-series.

SIGNIFICANCE

Our results suggest promising future applicability to investigate more complex contexts such as affective and pathological scenarios.

Measuring the Rate of Information Exchange in Point-Process Data With Application to Cardiovascular Variability

The amount of information exchanged per unit of time between two dynamic processes is an important concept for the analysis of complex systems. Theoretical formulations and data-efficient estimators have been recently introduced for this quantity, known as the mutual information rate (MIR), allowing its continuous-time computation for event-based data sets measured as realizations of coupled point processes. This work presents the implementation of MIR for point process applications in Network Physiology and cardiovascular variability, which typically feature short and noisy experimental time series. We assess the bias of MIR estimated for uncoupled point processes in the frame of surrogate data, and we compensate it by introducing a corrected MIR (cMIR) measure designed to return zero values when the two processes do not exchange information. The method is first tested extensively in synthetic point processes including a physiologically-based model of the heartbeat dynamics and the blood pressure propagation times, where we show the ability of cMIR to compensate the negative bias of MIR and return statistically significant values even for weakly coupled processes. The method is then assessed in real point-process data measured from healthy subjects during different physiological conditions, showing that cMIR between heartbeat and pressure propagation times increases significantly during postural stress, though not during mental stress. These results document that cMIR reflects physiological mechanisms of cardiovascular variability related to the joint neural autonomic modulation of heart rate and arterial compliance.

Baroreflex sensitivity with different lags and random forests for staging cardiovascular autonomic neuropathy in subjects with diabetes

Impaired baroreflex sensitivity (BRS) may indicate cardiovascular autonomic neuropathy (CAN), which often remains undiagnosed during the initial course of diabetes mellitus. The baroreflex mechanism can be considered negative feedback because of baroreflex delay, the time delay between a change in blood pressure and the counteracting heart rate response. This work sought to analyze BRS considering lags from 1 to 10 RR intervals. We hypothesized that diabetic patients with subclinical CAN (SCAN) have a detectable delay in autonomic nervous system activity and that this would differ from patients without CAN (NCAN) and with established CAN (ECAN). In the first stage, 30 patients were included in an exploratory analysis using the Principal Component Analysis. Six indexes related to the BRS delay were proposed and considered significant for staging diabetic patients. Three indexes allowed for the differentiating of patients with and without CAN, and three indexes distinguished subjects with SCAN from subjects with NCAN or ECAN. Then, in the second stage, a random forest model was developed with 72 subjects, using the variables selected in the first stage. It was possible to detect SCAN, and to point out those subjects with the potential to change the CAN stage, allowing for the tracking of CAN progression. The model achieved a sensitivity of 96% and specificity of 100% to detect SCAN. Thus, the BRS analysis considering delayed reaction in the dynamics of heart rate variability may contribute to an accurate screening tool to staging CAN, in addition to indicating patients with most insidious disease progress.

Causal analyses to study autonomic regulation during acute head-out water immersion, head-down tilt and supine position

NEW FINDINGS

What is the central question of this study? Can Granger causality analysis of R-R intervals, systolic blood pressure and respiration provide evidence for the different physiological mechanisms induced during thermoneutral water immersion, 6 deg head-down tilt and supine position tests that are not accessible using traditional heart rate variability and baroreflex methods? What is the main finding and its importance? The Granger analysis demonstrated a significant difference in the causal link from R-R intervals to respiration between water immersion and head-down tilt. The underlying physiological mechanism explaining this difference could be the variation in peripheral resistances.

ABSTRACT

Thermoneutral head-out water immersion (WI) and 6 deg head-down tilt (HDT) are used to simulate SCUBA diving, swimming and microgravity, because these models induce an increase in central blood volume. Standard methods to analyse autonomic regulation have demonstrated an increase in parasympathetic activity and baroreflex sensitivity during these experimental conditions. However, such methods are not adapted to quantify all closed-loop interactions involved in respiratory and cardiovascular regulation. To overcome this limitation, we used Granger causality analysis between R-R intervals (RR), systolic blood pressure (SBP) and respiration (RE) in eight young, healthy subjects, recorded during 30 min periods in the supine position, WI and HDT. For all experimental conditions, we found a bidirectional causal relationship between RE and RR and between RR and SBP, with a dominant direction from RR to SBP, and a unidirectional causality from RE to SBP. These causal relationships remained unchanged for the three experimental tests. Interestingly, there was a lower causal relationship from RR to RE during WI compared with HDT. This causal link from RR to RE could be modulated by peripheral resistances. These results highlight differences in cardiovascular regulation during WI and HDT and confirm that Granger causality might reveal physiological mechanisms not accessible with standard methods.

Vascular resistance arm of the baroreflex: methodology and comparison with the cardiac chronotropic arm

Baroreflex response consists of cardiac chronotropic (effect on heart rate), cardiac inotropic (on contractility), venous (on venous return) and vascular (on vascular resistance) arms. Because of the simplicity of its measurement, the cardiac chronotropic arm is most often analyzed. The aim was to introduce a method to assess the vascular baroreflex arm and to characterize its changes during stress. We evaluated the effect of orthostasis and mental arithmetics (MA) in 39 (22 women, 17 men; median age: 18.7 yr) and 36 (21 women, 15 men; 19.2 yr) healthy volunteers, respectively. We recorded systolic (SBP) and mean (MBP) blood pressure by volume-clamp method and R-R interval (RR) by ECG. Cardiac output (CO) was recorded by impedance cardiography. From MBP and CO, peripheral vascular resistance (PVR) was calculated. The directional spectral coupling and gain of cardiac chronotropic (SBP to RR) and vascular (SBP to PVR) arms were quantified. The strength of the causal coupling from SBP to PVR was significantly higher than that of SBP to RR coupling over the whole protocol (P < 0.001). Along both arms, the coupling was higher during orthostasis compared with the supine position (P < 0.001 and P = 0.006); no MA effect was observed. No significant changes in the spectral gain (ratio of RR or PVR change to a unit SBP change) across all phases were found (0.111 ≤ P ≤ 0.907). We conclude that changes in PVR are tightly coupled with SBP oscillations via the baroreflex, providing an approach for baroreflex vascular arm analysis with the potential to reveal new aspects of blood pressure dysregulation.NEW & NOTEWORTHY Baroreflex response consists of several arms, but the cardiac chronotropic arm (blood pressure changes evoking heart rate response) is usually analyzed. This study introduces a method to assess the vascular baroreflex arm with the continuous noninvasive measurement of peripheral vascular resistance as an output considering causality in the interaction between oscillations and slower dynamics of vascular tone changes. We conclude that although vascular baroreflex arm involvement becomes dominant during orthostasis, gain of this interaction is relatively stable.

Respiratory Sinus Arrhythmia Mechanisms in Young Obese Subjects

Autonomic nervous system (ANS) activity and imbalance between its sympathetic and parasympathetic components are important factors contributing to the initiation and progression of many cardiovascular disorders related to obesity. The results on respiratory sinus arrhythmia (RSA) magnitude changes as a parasympathetic index were not straightforward in previous studies on young obese subjects. Considering the potentially unbalanced ANS regulation with impaired parasympathetic control in obese patients, the aim of this study was to compare the relative contribution of baroreflex and non-baroreflex (central) mechanisms to the origin of RSA in obese vs. control subjects. To this end, we applied a recently proposed information-theoretic methodology - partial information decomposition (PID) - to the time series of heart rate variability (HRV, computed from RR intervals in the ECG), systolic blood pressure (SBP) variability, and respiration (RESP) pattern measured in 29 obese and 29 age- and gender-matched non-obese adolescents and young adults monitored in the resting supine position and during postural and cognitive stress evoked by head-up tilt and mental arithmetic. PID was used to quantify the so-called unique information transferred from RESP to HRV and from SBP to HRV, reflecting, respectively, non-baroreflex and RESP-unrelated baroreflex HRV mechanisms, and the redundant information transferred from (RESP, SBP) to HRV, reflecting RESP-related baroreflex RSA mechanisms. Our results suggest that obesity is associated: (i) with blunted involvement of non-baroreflex RSA mechanisms, documented by the lower unique information transferred from RESP to HRV at rest; and (ii) with a reduced response to postural stress (but not to mental stress), documented by the lack of changes in the unique information transferred from RESP and SBP to HRV in obese subjects moving from supine to upright, and by a decreased redundant information transfer in obese compared to controls in the upright position. These findings were observed in the presence of an unchanged RSA magnitude measured as the high frequency (HF) power of HRV, thus suggesting that the changes in ANS imbalance related to obesity in adolescents and young adults are subtle and can be revealed by dissecting RSA mechanisms into its components during various challenges.

Assessment of the Coupling Strength of Cardiovascular Control via Joint Symbolic Analysis during Postural Challenge in Recreational Athletes

Short-term cardiovascular control, comprising cardiac baroreflex and mechanisms governing cardiac contractility and vascular properties, links heart period (HP) and systolic arterial pressure (SAP) fluctuations. It is activated during postural challenge and this activation is traditionally quantified via linear tools such as HP-SAP squared coherence function. In this study the ability of a nonlinear bivariate tool based on joint symbolic analysis (JSA) approach was tested against HP-SAP coherence function during orthostatic challenge in recreational athletes. We studied 9 men healthy cycling practitioners (age: 20-40 yrs) at rest in supine condition (REST) and during active standing (STAND). The JSA method is based on the definition of symbolic HP and SAP patterns and on the evaluation of the rate of their simultaneous occurrence in both HP and SAP series. HP-SAP squared coherence was computed in the low frequency band (LF, from 0.04 to 0.15 Hz). We found the expected response to the postural stimulus, namely the increase of sympathetic modulation and the contemporaneous vagal withdrawal. However, only JSA was able to detect the expected increase of association between HP and SAP variability series over slow time scales. This result suggests that recreational athletes have more relevant nonlinear interactions between HP and SAP that might be missed by traditional linear tools and might require nonlinear tools to be efficiently described.

Strength and Latency of the HP-SAP Closed Loop Variability Interactions in Subjects Prone to Develop Postural Syncope

The coupling and latency between heart period (HP) and systolic arterial pressure (SAP) variability can be investigated along the two arms of the HP-SAP closed loop, namely along the baroreflex feedback from SAP to HP, and along the feedforward pathway from HP to SAP. This study investigates the HP-SAP closed loop variability interactions through cross-correlation function (CCF). Coupling strength and delay between HP and SAP variability series were monitored in 13 subjects prone to develop orthostatic syncope (SYNC, 28±9 yrs, 5 males) and in 13 subjects with no history of postural syncope (noSYNC, age: 27±8 yrs, 5 males). Analysis was carried out at rest in supine position (REST) and during head-up tilt (TILT) at 60 degrees. The null hypothesis of HP-SAP uncoupling was tested through a surrogate analysis associating the HP series of a subject with a SAP sequence of a different one. Results showed that during TILT the coupling strength increased along the baroreflex and latency augmented along the mechanical feedforward pathway exclusively in noSYNC subjects. Finally, closed loop HP-SAP interactions were detected in about one third of subjects and the situation of full uncoupling was rarely found. CCF analysis was found to be a straightforward and easily applicable method to investigate HP-SAP control deserving a direct comparison with more sophisticated signal processing tools assessing causality.

Dynamical interaction between heart rate and blood pressure of end-stage renal disease patients evaluated by cross recurrence plot diagonal analysis

The assessment of spontaneous variability of blood pressure and heart rate is based on specific physiological hypotheses about dynamic features, for example, the baroreflex modulation of heart rate over time in daily life. Usually, arterial baroreflex control of heart rate is explored without delays between blood pressure and heart rate data points, within a narrow range of values, excluding the analysis of saturation regions or low-threshold changes. In this work, we examine the dynamic interactions between systolic blood pressure (SBP) and interbeat interval (IBI), in 15-min length time series and for the first time using the analysis of diagonals derived from a cross-recurrence plots in healthy persons and end-stage renal disease (ESRD) patients. We found that ESRD patients have stronger intermittent dynamical interactions between IBI and SBP, but they lose most of the dynamical interactions. Although healthy subjects exhibit a continuously changing order of precedence between IBI and SBP at different lags, ESRD patients preserve this changing order of precedence only for lags >0 beats.NEW & NOTEWORTHY This study is the first to compare the time-variant pattern of systolic blood pressure (SBP) and interbeat interval (IBI) coupling between ESRD patients and healthy volunteers through the analysis of diagonal in cross-recurrence plots, and in the face of an orthostatic challenge. Our results demonstrated alternant interactions on the order of precedence (IBI → SBP or SBP→ IBI) at different time delays. This pattern is different in resting position and during active standing for the two groups studied, and interestingly, some association patterns are lost in ESRD patients. These patterns of alternant interactions on the order of precedence could be related to autonomic neural activities and cardiovascular synchronization at different scales both in time and space. This could reflect physiological adaptive flexibility of cardiovascular regulation. Losing some association patterns in ESRD may be the result of chronic adjustments of many physiological mechanisms (including chronic sympathetic hyperactivity), which could increase cardiovascular vulnerability to hemodynamic challenges.

Quantifying the Delays Between Multi-Site Photoplethysmography Pulse and Electrocardiogram R-R Interval Changes Under Slow-Paced Breathing

Objective: Objective assessment of autonomic function is important, including the investigation of slow-paced breathing to induce associated periodic changes in the cardiovascular system – such as blood pressure and heart rate. However, pulse changes across a range of peripheral body sites have seldom been explored with this challenge. The primary aim of this pilot study was to utilize multi-site photoplethysmography (MPPG) technology to quantify the phase delays, i.e., correlation lags, between changes in heart rate and changes in key pulse features with slow-paced breathing (0.1 Hz).

Methods: Waveforms were collected simultaneously from the right and left ear lobes, thumbs, and great toes of 18 healthy adult subjects. Cross correlation lags between reference beat-to-beat changes in electrocardiogram (ECG) R-R wave interval and changes in pulse arrival time (foot of pulse; PATf) and also for pulse amplitude (foot-to-peak; AMP) were determined.

Results: Relative to R-R changes, the median ear, thumb, and toe PATf correlation lags were 3.4, 2.9, and 2.1 beats, respectively; contrasting to AMP with 5.7, 6.0, and 6.9 beats, respectively. These PATf correlation lags in beats were significantly lower than for the AMP measure. Segmental differences between sites and timing measure variability have also been quantified.

Conclusion: This pilot study has indicated bilateral similarity plus segmental differences for relative delays in PPG pulse timing and amplitude measures relative to R-R interval changes with paced breathing. These correlation and variability data are now available for comparison with cardiovascular patient groups to support development of autonomic function assessment techniques.

Cardiac baroreflex hysteresis is one of the determinants of the heart period variability asymmetry

In heart period (HP) variability (HPV) recordings the percentage of negative HP variations tends to be greater than that of positive ones and this pattern is referred to as HPV asymmetry (HPVA). HPVA has been studied in several experimental conditions in healthy and pathological populations, but its origin is unclear. The baroreflex (BR) exhibits an asymmetric behavior as well given that it reacts more importantly to positive than negative arterial pressure (AP) variations. We tested the hypothesis that the BR asymmetry (BRA) is a HPVA determinant over spontaneous fluctuations of HP and systolic AP (SAP). We studied 100 healthy subjects (age from 21 to 70 yr, 54 men) comprising 20 subjects in each age decade. Electrocardiogram and noninvasive AP were recorded for 15 min at rest in supine position (REST) and during active standing (STAND). The HPVA was evaluated via Porta's index and Guzik's index, while the BRA was assessed as the difference, and normalized difference, between BR sensitivities computed over positive and negative SAP variations via the sequence method applied to HP and SAP variability. HPVA significantly increased during STAND and decreased progressively with age. BRA was not significantly detected both at REST and during STAND. However, we found a significant positive association between BRA and HPVA markers during STAND persisting even within the age groups. This study supports the use of HPVA indexes as descriptors of BRA and identified a challenge soliciting the BR response like STAND to maximize the association between HPVA and BRA markers.

Entropy Measures as Descriptors to Identify Apneas in Rheoencephalographic Signals

Rheoencephalography (REG) is a simple and inexpensive technique that intends to monitor cerebral blood flow (CBF), but its ability to reflect CBF changes has not been extensively proved. Based on the hypothesis that alterations in CBF during apnea should be reflected in REG signals under the form of increased complexity, several entropy metrics were assessed for REG analysis during apnea and resting periods in 16 healthy subjects: approximate entropy (ApEn), sample entropy (SampEn), fuzzy entropy (FuzzyEn), corrected conditional entropy (CCE) and Shannon entropy (SE). To compute these entropy metrics, a set of parameters must be defined a priori, such as, for example, the embedding dimension m, and the tolerance threshold r. A thorough analysis of the effects of parameter selection in the entropy metrics was performed, looking for the values optimizing differences between apnea and baseline signals. All entropy metrics, except SE, provided higher values for apnea periods (p-values < 0.025). FuzzyEn outperformed all other metrics, providing the lowest p-value (p = 0.0001), allowing to conclude that REG signals during apnea have higher complexity than in resting periods. Those findings suggest that REG signals reflect CBF changes provoked by apneas, even though further studies are needed to confirm this hypothesis.

Phase synchronization of baroreflex oscillations of blood pressure and pulse interval in rats: the effects of cardiac autonomic blockade and gradual blood loss

OBJECTIVE

Phase synchronization of arterial pressure (AP) and pulse interval (PI) oscillations in the low-frequency band (around 0.4 Hz in rats) is governed by baroreflex activity. In long-term stationary data recordings, such synchronization can be estimated by the coherence. The phase synchronization index (PSI) can be used as well. The aim of this study was to correlate PSI and the coherence of AP and PI under stationary conditions and to estimate the informativity of PSI as a measure of baroreflex activity during transient processes.

APPROACH

AP and PI were recorded in conscious Wistar rats using femoral artery catheters. To study the hemodynamics during hemorrhage, blood was gradually withdrawn (20 ml × kg^-1 over 30 min) through a catheter in the carotid artery.

MAIN RESULTS

PSI and coherence spectra calculated from 30-minute AP and PI recordings demonstrated distinct peaks at the frequency of 0.4 Hz; these indicators correlate well with each other (Pearson r  =  0.920, p   <  0.0001). Both PSI and coherence were markedly suppressed by vagal blockade (methylatropine) and tended to reduce after sympathetic blockade (atenolol). Importantly, PSI demonstrated dynamic alterations during gradual hemorrhage. During the initial approx. 10 min of hemorrhage, AP did not change but PI was noticeably shortened, and PSI increased, which indicates the activation of the baroreflex. With further blood loss, baroreflex influences were not enough to prevent blood pressure from falling, and under such conditions PSI decreased.

SIGNIFICANCE

PSI, like coherence, is an informative measure of baroreflex activity under stationary conditions. In addition, PSI permits us to follow the coupling between the baroreflex oscillations of AP and PI during transient processes, which strengthens its informative value.

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.

Mental Workload Alters Heart Rate Variability, Lowering Non-linear Dynamics

Mental workload is known to alter cardiovascular function leading to increased cardiovascular risk. Nevertheless, there is no clear autonomic nervous system unbalance to be quantified during mental stress. We aimed to characterize the mental workload impact on the cardiovascular function with a focus on heart rate variability (HRV) non-linear indexes. A 1-h computerized switching task (letter recognition) was performed by 24 subjects while monitoring their performance (accuracy, response time), electrocardiogram and blood pressure waveform (finger volume clamp method). The HRV was evaluated from the beat-to-beat RR intervals (RRI) in time-, frequency-, and informational- domains, before (Control) and during the task. The task induced a significant mental workload (visual analog scale of fatigue from 27 ± 26 to 50 ± 31 mm, p < 0.001, and NASA-TLX score of 56 ± 17). The heart rate, blood pressure and baroreflex function were unchanged, whereas most of the HRV parameters markedly decreased. The maximum decrease occurred during the first 15 min of the task (P1), before starting to return to the baseline values reached at the end of the task (P4). The RRI dimension correlation (D2) decrease was the most significant (P1 vs. Control: 1.42 ± 0.85 vs. 2.21 ± 0.8, p < 0.001) and only D2 lasted until the task ended (P4 vs. Control: 1.96 ± 0.9 vs. 2.21 ± 0.9, p < 0.05). D2 was identified as the most robust cardiovascular variable impacted by the mental workload as determined by posterior predictive simulations (p = 0.9). The Spearman correlation matrix highlighted that D2 could be a marker of the generated frustration (R = -0.61, p < 0.01) induced by a mental task, as well as the myocardial oxygen consumption changes assessed by the double product (R = -0.53, p < 0.05). In conclusion, we showed that mental workload sharply lowered the non-linear RRI dynamics, particularly the RRI correlation dimension.

Time-Frequency Analysis of Cardiovascular and Cardiorespiratory Interactions During Orthostatic Stress by Extended Partial Directed Coherence

In this study, the linear method of extended partial directed coherence (ePDC) was applied to establish the temporal dynamic behavior of cardiovascular and cardiorespiratory interactions during orthostatic stress at a 70° head-up tilt (HUT) test on young age-matched healthy subjects and patients with orthostatic intolerance (OI), both male and female. Twenty 5-min windows were used to analyze the minute-wise progression of interactions from 5 min in a supine position (baseline, BL) until 18 min of the orthostatic phase (OP) without including pre-syncopal phases. Gender differences in controls were present in cardiorespiratory interactions during OP without compromised autonomic regulation. However in patients, analysis by ePDC revealed considerable dynamic alterations within cardiovascular and cardiorespiratory interactions over the temporal course during the HUT test. Considering the young female patients with OI, the information flow from heart rate to systolic blood pressure (mechanical modulation) was already increased before the tilt-up, the information flow from systolic blood pressure to heart rate (neural baroreflex) increased during OP, while the information flow from respiration to heart rate (respiratory sinus arrhythmia) decreased during the complete HUT test. Findings revealed impaired cardiovascular interactions in patients with orthostatic intolerance and confirmed the usefulness of ePDC for causality analysis.

Baroreflex sensitivity in frailty syndrome

Frailty is related to a decrease in the physiological reserves, which causes difficulties in maintaining homeostasis. An example of physiological mechanisms for cardiovascular homeostasis is the baroreflex. The aim of this study was to compare baroreflex among frail, prefrail, and nonfrail individuals, in supine and orthostatic positions. Community-dwelling older adults were evaluated and categorized into frail, prefrail, or nonfrail groups, according to frailty phenotype. The RR interval (RRi) and systolic blood pressure (SBP) series were recorded for 15 min in the supine and 15 min in the orthostatic positions. Mean and variance of RRi and SBP, and baroreflex evaluated by phase, gain (α), and coherence (K2) were determined. A two-way repeated measures ANOVA, with Tukey's post hoc, was applied for group, position, and their interaction effects. The significance level established was 5%. Prefrail and frail participants did not present a significant decrease in mean values of RRi after postural challenge (893.43 to 834.20 ms and 925.99 to 857.98 ms, respectively). Frail participants showed a reduction in RRi variance in supine to orthostatic (852.04 to 232.37 ms2). Prefrail and frail participants showed a decrease in K2 after postural change (0.69 to 0.52 and 0.54 to 0.34, respectively). Frail participants exhibited lower values of K2 (0.34) compared to nonfrail and prefrail participants (0.61 and 0.52, respectively). Baroreflex indicated the presence of decoupling between heart period and SBP in frail and prefrail. Thus, reduced K2 might be a marker of the frailty process.

Nonlinear Interaction Analysis of Cardiovascular-Respiratory Data by Means of Convergent Cross Mapping

Appropriate analyses of directed complex interactions within the cardiovascular-respiratory system are of growing interest for a better understanding of physiological regulatory mechanisms in healthy subjects and diseased persons. There are various concepts to analyze such interactions. Convergent Cross Mapping (CCM) provides the possibility to define directed interactions in terms of nonlinear stability. A proof-of-principle approach is introduced to apply CCM to cardiovascular-respiratory data of healthy subjects during resting state period. Showing group results of time-invariant as well as single subject results of interval-based CCM, the introduced approach was able to quantify correct directionality and strength of interactions within the cardiovascular-respiratory system and to provide statistical thresholds for significant interactions. These results may serve as a methodological base to compare healthy subjects and diseased persons.

Comparison between probabilistic and Wiener-Granger causality in assessing modifications of the cardiac baroreflex control with age

BACKGROUND

Probabilistic causality (PC) is a framework for checking that the occurrence of a cause raises the probability of the effect by comparing the probability of the effect conditioned and unconditioned to the cause. Even though it is less frequently utilized with respect to the more traditional model-based Wiener-Granger causality (WGC) that is based on the predictability improvement of an effect resulting from the inclusion of the presumed cause in the multivariate linear regression model, PC has the advantage of being model-free.

OBJECTIVE

The aim of the study is to apply the PC framework to assess the evolution of cardiac baroreflex control with age from spontaneous fluctuations of heart period (HP) and systolic arterial pressure (SAP) and to compare it to the more common WGC approach.

APPROACH

We studied 100 healthy humans (54 males, age: 21-70 years, 20 subjects for each 10 years bin). HP and SAP were extracted on a beat-to-beat basis from 5 min recordings of electrocardiogram and plethysmographic arterial pressure at rest in supine position (REST) and during active standing (STAND) under spontaneous breathing. The WGC ratio (WGCR) was computed as the log ratio of the prediction error variance of the autoregressive model on HP to that on HP with exogenous SAP. The PC ratio (PCR) was computed as the probability of observing an HP ramp given an associated parallel SAP variation divided by the probability of observing an HP ramp.

MAIN RESULTS

The WGCR and PCR suggested the gradual impairment of cardiac baroreflex with age, especially during STAND. Moreover, they were significantly associated both at REST and during STAND but the degree of the PCR-WGCR association was weak.

SIGNIFICANCE

PC can be effectively exploited to assess modification of the cardiovascular control during senescence even though a limited agreement was observed with WGC.

Analysis for the Influence of ABR Sensitivity on PTT-Based Cuff-Less Blood Pressure Estimation before and after Exercise

An accurate and continuous measurement of blood pressure (BP) is of great importance for the prognosis of some cardiovascular diseases in out-of-hospital settings. Pulse transit time (PTT) is a well-known cardiovascular parameter which is highly correlated with BP and has been widely applied in the estimation of continuous BP. However, due to the complexity of cardiovascular system, the accuracy of PTT-based BP estimation is still unsatisfactory. Recent studies indicate that, for the subjects before and after exercise, PTT can track the high-frequency BP oscillation (HF-BP) well, but is inadequate to follow the low-frequency BP variance (LF-BP). Unfortunately, the cause for this failure of PTT in LF-BP estimation is still unclear. Based on these previous researches, we investigated the cause behind this failure of PTT in LF-BP estimation. The heart rate- (HR-) related arterial baroreflex (ABR) model was introduced to analyze the failure of PTT in LF-BP estimation. Data from 42 healthy volunteers before and after exercise were collected to evaluate the correlation between the ABR sensitivity and the estimation error of PTT-based BP in LF and HF components. In the correlation plot, an obvious difference was observed between the LF and HF groups. The correlation coefficient r for the ABR sensitivity with the estimation error of systolic BP (SBP) and diastolic BP (DBP) in LF was 0.817 ± 0.038 and 0.757 ± 0.069, respectively. However, those correlation coefficient r for the ABR sensitivity with the estimation error of SBP and DBP in HF was only 0.403 ± 0.145 and 0.274 ± 0.154, respectively. These results indicated that there is an ABR-related complex LF autonomic regulation mechanism on BP, PTT, and HR, which influences the effect of PTT in LF-BP estimation.

Study of impaired cardiovascular and respiratory coupling during orthostatic stress based on joint symbolic dynamics

The present study investigates the instantaneous coupling among the cardiac, vascular, and respiratory systems, using the heart rate, respiration, and systolic and diastolic blood pressure variability in 12 healthy and 16 vasovagal syncope female subjects during a head-up tilt (HUT) testing protocol at 70° This study contributes to the coupling analysis by using a nonlinear joint symbolic dynamics (JSD) in a high-temporal resolution scheme, based on 5 min segments of the time series that are shifted every minute. For each segment, a bivariate JSD matrix was constructed to obtain global and local coupling indices in accordance to Shannon's entropy and the probability of occurrence of various bivariate words, respectively. The novel approach revealed important findings in the coupling dynamics of the systems, thus allowing the detection of group differences during the early orthostatic phase, and during the HUT test, before the occurrence of any pre-syncopal symptoms. In patients, the global indices indicated a significant decrease of cardiovascular coupling, starting at 10 min after the tilt-up, manifested by reduced baroreflex sensitivity and cardiorespiratory coupling that was initiated 8 min after the onset of the orthostatic phase (OP). A decreased autonomic control on cardiovascular-respiratory couplings was further evidenced by increased alterations of the JSD indices during the OP compared to the supine position in patients compared to controls. Furthermore, findings based on local indices demonstrated that female patients showed reductions and disengagements in cardiovascular (p < 0.001) and cardiorespiratory (p < 0.01) couplings, as early as the first 5 min and during the complete OP.

Analysis of diagonals in cross recurrence plots between heart rate and systolic blood pressure during supine position and active standing in healthy adults

The inter beat interval (IBI) duration and systolic blood pressure (SBP) are cardiovascular variables related through several feedback mechanisms. We propose the analysis of diagonal lines in cross recurrence plots (CRPs) from IBI and SBP embedded within the same phase space to identify events where trajectories of both variables concur. The aim of the study was to describe the relationship between IBI and SBP of healthy subjects using CRP and diagonal analysis during baseline condition-supine position (SP)-and how the relationship changes during the physiological stress of active standing (AS). IBI and SBP time series were obtained from continuous blood pressure recordings during SP and AS (15 min each) in 19 young healthy subjects. IBI and SBP time series were embedded within a five-dimensional phase space using an embedding delay estimated from cross correlation between IBI and SBP. During SP, mean CRP showed high determinism (≥85%) and also brief but repeated events where both variables stay within a reduced space. Most quantitative recurrences analysis indexes of CRP increased significantly (p < 0.05) during AS. CRP analysis showed short diagonals indicating a very strong deterministic relationship between IBI and SBP with intermittent unlocking periods. The strength of IBI and SBP relationship increased during the physiological stress of AS. The CRP method allowed a rigorous quantitative description of the deterministic association between these two variables. Diagonal lines were intermittent and not always parallel, showing that there is not a defined and unique rhythm. This suggests the activation of different influences at different times and with different precedence between the heart rate and blood pressure in response to AS.

Comparison of Autonomic Control of Blood Pressure During Standing and Artificial Gravity Induced via Short-Arm Human Centrifuge

Autonomic control of blood pressure is essential toward maintenance of cerebral perfusion during standing, failure of which could lead to fainting. Long-term exposure to microgravity deteriorates autonomic control of blood pressure. Consequently, astronauts experience orthostatic intolerance on their return to gravitational environment. Ground-based studies suggest sporadic training in artificial hypergravity can mitigate spaceflight deconditioning. In this regard, short-arm human centrifuge (SAHC), capable of creating artificial hypergravity of different g-loads, provides an auspicious training tool. Here, we compare autonomic control of blood pressure during centrifugation creating 1-g and 2-g at feet with standing in natural gravity. Continuous blood pressure was acquired simultaneously from 13 healthy participants during supine baseline, standing, supine recovery, centrifugation of 1-g, and 2-g, from which heart rate (RR) and systolic blood pressure (SBP) were derived. The autonomic blood pressure regulation was assessed via spectral analysis of RR and SBP, spontaneous baroreflex sensitivity, and non-linear heart rate and blood pressure causality (RR↔SBP). While majority of these blood pressure regulatory indices were significantly different (p < 0.05) during standing and 2-g centrifugation compared to baseline, no change (p > 0.05) was observed in the same indices during 2-g centrifugation compared to standing. The findings of the study highlight the capability of artificial gravity (2-g at feet) created via SAHC toward evoking blood pressure regulatory controls analogous to standing, therefore, a potential utility toward mitigating deleterious effects of microgravity on cardiovascular performance and minimizing post-flight orthostatic intolerance in astronauts.

Role of respiration in the cardiovascular response to orthostatic and mental stress

The objective of this study was to determine the response of heart rate and blood pressure variability (respiratory sinus arrhythmia, baroreflex sensitivity) to orthostatic and mental stress, focusing on causality and the mediating effect of respiration. Seventy-seven healthy young volunteers (46 women, 31 men) aged 18.4 ± 2.7 yr underwent an experimental protocol comprising supine rest, 45° head-up tilt, recovery, and a mental arithmetic task. Heart rate variability and blood pressure variability were analyzed in the time and frequency domain and modeled as a multivariate autoregressive process where the respiratory volume signal acted as an external driver. During head-up tilt, tidal volume increased while respiratory rate decreased. During mental stress, breathing rate increased and tidal volume was elevated slightly. Respiratory sinus arrhythmia decreased during both interventions. Baroreflex function was preserved during orthostasis but was decreased during mental stress. While sex differences were not observed during baseline conditions, cardiovascular response to orthostatic stress and respiratory response to mental stress was more prominent in men compared with women. The respiratory response to the mental arithmetic tasks was more prominent in men despite a significantly higher subjectively perceived stress level in women. In conclusion, respiration shows a distinct response to orthostatic versus mental stress, mediating cardiovascular variability; it needs to be considered for correct interpretation of heart rate and blood pressure phenomena.

Information transfer and information modification to identify the structure of cardiovascular and cardiorespiratory networks

To fully elucidate the complex physiological mechanisms underlying the short-term autonomic regulation of heart period (H), systolic and diastolic arterial pressure (S, D) and respiratory (R) variability, the joint dynamics of these variables need to be explored using multivariate time series analysis. This study proposes the utilization of information-theoretic measures to measure causal interactions between nodes of the cardiovascular/cardiorespiratory network and to assess the nature (synergistic or redundant) of these directed interactions. Indexes of information transfer and information modification are extracted from the H, S, D and R series measured from healthy subjects in a resting state and during postural stress. Computations are performed in the framework of multivariate linear regression, using bootstrap techniques to assess on a single-subject basis the statistical significance of each measure and of its transitions across conditions. We find patterns of information transfer and modification which are related to specific cardiovascular and cardiorespiratory mechanisms in resting conditions and to their modification induced by the orthostatic stress.

Non-linear Heart Rate and Blood Pressure Interaction in Response to Lower-Body Negative Pressure

Early detection of hemorrhage remains an open problem. In this regard, blood pressure has been an ineffective measure of blood loss due to numerous compensatory mechanisms sustaining arterial blood pressure homeostasis. Here, we investigate the feasibility of causality detection in the heart rate and blood pressure interaction, a closed-loop control system, for early detection of hemorrhage. The hemorrhage was simulated via graded lower-body negative pressure (LBNP) from 0 to -40 mmHg. The research hypothesis was that a significant elevation of causal control in the direction of blood pressure to heart rate (i.e., baroreflex response) is an early indicator of central hypovolemia. Five minutes of continuous blood pressure and electrocardiogram (ECG) signals were acquired simultaneously from young, healthy participants (27 ± 1 years, N = 27) during each LBNP stage, from which heart rate (represented by RR interval), systolic blood pressure (SBP), diastolic blood pressure (DBP), and mean arterial pressure (MAP) were derived. The heart rate and blood pressure causal interaction (RR↔SBP and RR↔MAP) was studied during the last 3 min of each LBNP stage. At supine rest, the non-baroreflex arm (RR→SBP and RR→MAP) showed a significantly (p < 0.001) higher causal drive toward blood pressure regulation compared to the baroreflex arm (SBP→RR and MAP→RR). In response to moderate category hemorrhage (-30 mmHg LBNP), no change was observed in the traditional marker of blood loss i.e., pulse pressure (p = 0.10) along with the RR→SBP (p = 0.76), RR→MAP (p = 0.60), and SBP→RR (p = 0.07) causality compared to the resting stage. Contrarily, a significant elevation in the MAP→RR (p = 0.004) causality was observed. In accordance with our hypothesis, the outcomes of the research underscored the potential of compensatory baroreflex arm (MAP→RR) of the heart rate and blood pressure interaction toward differentiating a simulated moderate category hemorrhage from the resting stage. Therefore, monitoring baroreflex causality can have a clinical utility in making triage decisions to impede hemorrhage progression.

Impaired Baroreflex Function during Orthostatic Challenge in Patients after Spinal Cord Injury

The level of spinal cord injury (SCI) affects baroreflex regulation of blood pressure. While a parasympathetic cardiac chronotropic effect is preserved, baroreflex response could be impaired by sympathetic dysfunction under the SCI level. This study was aimed to evaluate the baroreflex function in SCI patients by the analysis of causal interaction between systolic blood pressure (SBP) and inter-beat intervals (IBI). Blood pressure was continuously recorded in 13 cervical SCI patients (CSCI), nine thoracic SCI (ThSCI) and 13 able-bodied controls (Con) during two phases: sitting (PS) and orthostatic challenge (PO). Beat-to-beat SBP and IBI sequences were obtained from continuous blood pressure recording. Closed loop of SBP-IBI interaction was mathematically opened by bivariate autoregressive model; causal coherence and baroreflex sensitivity (BRS) were calculated in baroreflex direction. Coherence quantifies causal synchronicity between SBP and IBI. The gain of transfer function from SBP to IBI represents BRS. PS (medians of CSCI/ThSCI/Con) coherence was 0.28/0.33/0.25 (no significant difference) and PS BRS was 6.98/7.54/6.66 (no difference). PO coherence was 0.18/0.58/0.45 (CSCI < ThCSI and Con; p < 0.01) and PO BRS was 2.38/5.87/6.22 (CSCI < ThCSI and Con; p < 0.01). For position change effect, there was no change in CSCI coherence; for ThSCI and Con, PS < PO (p < 0.05). For BRS in the CSCI group, PS < PO (p < 0.01); for ThSCI and Con, there was no change. BRS and coherence correlated negatively with SCI level (p < 0.01). In conclusion, baroreflex dysfunction in SCI patients was detected using causal analysis methods during orthostatic challenge only. Baroreflex dysfunction is probably an important mechanism of the more expressed blood pressure decrease associated with CSCI. The severity of autonomic dysfunction was related to SCI level.

Analysis of causal cardio-postural interaction under orthostatic stress using convergent cross mapping

Knowledge of a cause-and-effect relationship between different physiological systems is helpful in predicting their performance under perturbations, such as orthostatic challenge. The causal coupling between representative signals of the cardiovascular and postural systems under orthostatic challenge remains unknown. Understanding the causal relationship between these two systems is critical, as their interplay is vital to maintain stable upright posture of the human body during quiet standing. In this research, convergent cross mapping (CCM) method was applied to study the causal relationship between the cardiovascular and postural systems previously shown to have coherent activity during quiet standing. Causality was studied between Systolic blood pressure (SBP)-EMG (calf muscles), EMG-COPr (resultant center of pressure), and COPr-SBP signal pairs. These signals were simultaneously recorded in a 5-minute sit-to-stand test from five young healthy participants. Strength of causality was obtained between the signal pairs in a 30-second time segments. The results from this study indicate that there exists a bidirectional causal relationship between the cardio-postural signal pairs, indicating a system level interaction to counter perturbation due to orthostatic challenge. Skeletal muscle pump was found to be driving control of SBP and COPr as the value of EMG→SBP (0.54±0.09) and EMG→COPr (0.52±0.07) were higher than the reverse causality of SBP→EMG (0.19±0.16) and COPr→EMG (0.29±0.16).

Multivariate symbolic dynamics for analysis of respiratory-cardiovascular interactions

In this work, a nonlinear method to study multivariate interactions, called multivariate symbolic dynamics (MSD), was introduced. The usefulness of this technique was studied on respiratory-cardiovascular data from young women with vasovagal syncope (VVS) and from healthy subjects. The study included 16 female patients diagnosed with VVS and 24 age-matched healthy subjects (12 women). All subjects were enrolled in a head-up tilt (HUT) test, breathing normally, including 5 min of supine position and 18 to 28 min of 70° orthostatic phase. The MSD parameters were dynamically obtained for 5-min windows shifted by 1 min during HUT test. In supine position there were no considerable differences. During orthostatic phase, parameters from MSD showed a highly significantly (p=0.00005) increased occurrence of impaired respiratory-cardiovascular interactions in female patients susceptible to vasovagal syncope. This study provided promising results for a new multivariate method to investigate respiratory-cardiovascular interactions.

Physiological Mechanisms Mediating the Coupling between Heart Period and Arterial Pressure in Response to Postural Changes in Humans

The upright posture strengthens the coupling between heart period (HP) and systolic arterial pressure (SAP) consistently with a greater contribution of the arterial baroreflex to cardiac control, while paradoxically decreasing cardiac baroreflex sensitivity (cBRS). To investigate the physiological mechanisms that mediate the coupling between HP and SAP in response to different postures, we analyzed the cross-correlation functions between low-frequency HP and SAP fluctuations and estimated cBRS with the sequence technique in healthy male subjects during passive head-up tilt test (HUTT, n = 58), during supine wakefulness, supine slow-wave sleep (SWS), and in the seated and active standing positions (n = 8), and during progressive loss of 1 L blood (n = 8) to decrease central venous pressure in the supine position. HUTT, SWS, the seated, and the standing positions, but not blood loss, entailed significant increases in the positive correlation between HP and the previous SAP values, which is the expected result of arterial baroreflex control, compared with baseline recordings in the supine position during wakefulness. These increases were mirrored by increases in the low-frequency variability of SAP in each condition but SWS. cBRS decreased significantly during HUTT, in the seated and standing positions, and after blood loss compared with baseline during wakefulness. These decreases were mirrored by decreases in the RMSSD index, which reflects cardiac vagal modulation. These results support the view that the cBRS decrease associated with the upright posture is a byproduct of decreased cardiac vagal modulation, triggered by the arterial baroreflex in response to central hypovolemia. Conversely, the greater baroreflex contribution to cardiac control associated with upright posture may be explained, at least in part, by enhanced fluctuations of SAP, which elicit a more effective entrainment of HP fluctuations by the arterial baroreflex. These SAP fluctuations may result from enhanced fluctuations of vascular resistance specific to the upright posture, and not be driven by the accompanying central hypovolemia.

Skeletal Muscle Pump Drives Control of Cardiovascular and Postural Systems

The causal interaction between cardio-postural-musculoskeletal systems is critical in maintaining postural stability under orthostatic challenge. The absence or reduction of such interactions could lead to fainting and falls often experienced by elderly individuals. The causal relationship between systolic blood pressure (SBP), calf electromyography (EMG), and resultant center of pressure (COPr) can quantify the behavior of cardio-postural control loop. Convergent cross mapping (CCM) is a non-linear approach to establish causality, thus, expected to decipher nonlinear causal cardio-postural-musculoskeletal interactions. Data were acquired simultaneously from young participants (25 ± 2 years, n = 18) during a 10-minute sit-to-stand test. In the young population, skeletal muscle pump was found to drive blood pressure control (EMG → SBP) as well as control the postural sway (EMG → COPr) through the significantly higher causal drive in the direction towards SBP and COPr. Furthermore, the effect of aging on muscle pump activation associated with blood pressure regulation was explored. Simultaneous EMG and SBP were acquired from elderly group (69 ± 4 years, n = 14). A significant (p = 0.002) decline in EMG → SBP causality was observed in the elderly group, compared to the young group. The results highlight the potential of causality to detect alteration in blood pressure regulation with age, thus, a potential clinical utility towards detection of fall proneness.

Basic cardiovascular variability signals: mutual directed interactions explored in the information domain

The study of short-term cardiovascular interactions is classically performed through the bivariate analysis of the interactions between the beat-to-beat variability of heart period (RR interval from the ECG) and systolic blood pressure (SBP). Recent progress in the development of multivariate time series analysis methods is making it possible to explore how directed interactions between two signals change in the context of networks including other coupled signals. Exploiting these advances, the present study aims at assessing directional cardiovascular interactions among the basic variability signals of RR, SBP and diastolic blood pressure (DBP), using an approach which allows direct comparison between bivariate and multivariate coupling measures. To this end, we compute information-theoretic measures of the strength and delay of causal interactions between RR, SBP and DBP using both bivariate and trivariate (conditioned) formulations in a group of healthy subjects in a resting state and during stress conditions induced by head-up tilt (HUT) and mental arithmetics (MA). We find that bivariate measures better quantify the overall (direct  +  indirect) information transferred between variables, while trivariate measures better reflect the existence and delay of directed interactions. The main physiological results are: (i) the detection during supine rest of strong interactions along the pathway RR  →  DBP  →  SBP, reflecting marked Windkessel and/or Frank-Starling effects; (ii) the finding of relatively weak baroreflex effects SBP  →  RR at rest; (iii) the invariance of cardiovascular interactions during MA, and the emergence of stronger and faster SBP  →  RR interactions, as well as of weaker RR  →  DBP interactions, during HUT. These findings support the importance of investigating cardiovascular interactions from a network perspective, and suggest the usefulness of directed information measures to assess physiological mechanisms and track their changes across different physiological states.

Development of Causal Interactions Between Systolic Blood Pressure and Inter-Beat Intervals in Adolescents

Systolic blood pressure (SBP) changes control the cardiac inter-beat intervals (IBI) duration via baroreflex. Conversely, SBP is influenced by IBI via non-baroreflex mechanisms. Both causal pathways (feedback – baroreflex and feedforward – non-baroreflex) form a closed loop of the SBP – IBI interaction. The aim of this study was to assess the age-related changes in the IBI – SBP interaction. We have non-invasively recorded resting beat-to-beat SBP and IBI in 335 healthy subjects of different age, ranging from 11 to 23 years. Using a linear autoregressive bivariate model we obtained gain (GainSBP,IBI, used traditionally as baroreflex sensitivity) and coherence (CohSBP,IBI) of the SBP–IBI interaction and causal gain and coherence in baroreflex (GainSBPIBI, CohSBPIBI) and coherence in non-baroreflex (CohIBISBP) directions separately. A non-linear approach was used for causal coupling indices evaluation (CSBPIBI, CIBISBP) quantifying the amount of information transferred between signals. We performed a correlation to age analysis of all measures. CohIBISBP and CIBISBP were higher than CohSBPIBI and CSBPIBI, respectively. GainSBP,IBI increased and CohSBPIBI decreased with age. The coupling indices did not correlate with age. We conclude that the feedforward influence dominated at rest. The increase of GainSBP,IBI with age was not found in the closed loop model. A decrease of CohSBPIBI could be related to a change in the cardiovascular control system complexity during maturation.

Algorithms for the inference of causality in dynamic processes: Application to cardiovascular and cerebrovascular variability

This study faces the problem of causal inference in multivariate dynamic processes, with specific regard to the detection of instantaneous and time-lagged directed interactions. We point out the limitations of the traditional Granger causality analysis, showing that it leads to false detection of causality when instantaneous and time-lagged effects coexist in the process structure. Then, we propose an improved algorithm for causal inference that combines the Granger framework with the approach proposed by Pearl for the study of causality among multiple random variables. This new approach is compared with the traditional one in theoretical and simulated examples of interacting processes, showing its ability to retrieve the correct structure of instantaneous and time-lagged interactions. These approaches for causal inference are then tested on the physiological variability series of heart period, arterial pressure and cerebral blood flow variability obtained in subjects with postural-related syncope during a tilt-test protocol.

Estimating the decomposition of predictive information in multivariate systems

In the study of complex systems from observed multivariate time series, insight into the evolution of one system may be under investigation, which can be explained by the information storage of the system and the information transfer from other interacting systems. We present a framework for the model-free estimation of information storage and information transfer computed as the terms composing the predictive information about the target of a multivariate dynamical process. The approach tackles the curse of dimensionality employing a nonuniform embedding scheme that selects progressively, among the past components of the multivariate process, only those that contribute most, in terms of conditional mutual information, to the present target process. Moreover, it computes all information-theoretic quantities using a nearest-neighbor technique designed to compensate the bias due to the different dimensionality of individual entropy terms. The resulting estimators of prediction entropy, storage entropy, transfer entropy, and partial transfer entropy are tested on simulations of coupled linear stochastic and nonlinear deterministic dynamic processes, demonstrating the superiority of the proposed approach over the traditional estimators based on uniform embedding. The framework is then applied to multivariate physiologic time series, resulting in physiologically well-interpretable information decompositions of cardiovascular and cardiorespiratory interactions during head-up tilt and of joint brain-heart dynamics during sleep.