Cardiorespiratory Information Dynamics during Mental Arithmetic and Sustained Attention

Testing dynamic correlations and nonlinearity in bivariate time series through information measures and surrogate data analysis

The increasing availability of time series data depicting the evolution of physical system properties has prompted the development of methods focused on extracting insights into the system behavior over time, discerning whether it stems from deterministic or stochastic dynamical systems. Surrogate data testing plays a crucial role in this process by facilitating robust statistical assessments. This ensures that the observed results are not mere occurrences by chance, but genuinely reflect the inherent characteristics of the underlying system. The initial process involves formulating a null hypothesis, which is tested using surrogate data in cases where assumptions about the underlying distributions are absent. A discriminating statistic is then computed for both the original data and each surrogate data set. Significantly deviating values between the original data and the surrogate data ensemble lead to the rejection of the null hypothesis. In this work, we present various surrogate methods designed to assess specific statistical properties in random processes. Specifically, we introduce methods for evaluating the presence of autodependencies and nonlinear dynamics within individual processes, using Information Storage as a discriminating statistic. Additionally, methods are introduced for detecting coupling and nonlinearities in bivariate processes, employing the Mutual Information Rate for this purpose. The surrogate methods introduced are first tested through simulations involving univariate and bivariate processes exhibiting both linear and nonlinear dynamics. Then, they are applied to physiological time series of Heart Period (RR intervals) and respiratory flow (RESP) variability measured during spontaneous and paced breathing. Simulations demonstrated that the proposed methods effectively identify essential dynamical features of stochastic systems. The real data application showed that paced breathing, at low breathing rate, increases the predictability of the individual dynamics of RR and RESP and dampens nonlinearity in their coupled dynamics.

Your blush gives you away: detecting hidden mental states with remote photoplethysmography and thermal imaging

Multimodal emotion recognition techniques are increasingly essential for assessing mental states. Image-based methods, however, tend to focus predominantly on overt visual cues and often overlook subtler mental state changes. Psychophysiological research has demonstrated that heart rate (HR) and skin temperature are effective in detecting autonomic nervous system (ANS) activities, thereby revealing these subtle changes. However, traditional HR tools are generally more costly and less portable, while skin temperature analysis usually necessitates extensive manual processing. Advances in remote photoplethysmography (r-PPG) and automatic thermal region of interest (ROI) detection algorithms have been developed to address these issues, yet their accuracy in practical applications remains limited. This study aims to bridge this gap by integrating r-PPG with thermal imaging to enhance prediction performance. Ninety participants completed a 20-min questionnaire to induce cognitive stress, followed by watching a film aimed at eliciting moral elevation. The results demonstrate that the combination of r-PPG and thermal imaging effectively detects emotional shifts. Using r-PPG alone, the prediction accuracy was 77% for cognitive stress and 61% for moral elevation, as determined by a support vector machine (SVM). Thermal imaging alone achieved 79% accuracy for cognitive stress and 78% for moral elevation, utilizing a random forest (RF) algorithm. An early fusion strategy of these modalities significantly improved accuracies, achieving 87% for cognitive stress and 83% for moral elevation using RF. Further analysis, which utilized statistical metrics and explainable machine learning methods including SHapley Additive exPlanations (SHAP), highlighted key features and clarified the relationship between cardiac responses and facial temperature variations. Notably, it was observed that cardiovascular features derived from r-PPG models had a more pronounced influence in data fusion, despite thermal imaging's higher predictive accuracy in unimodal analysis.

Time-varying information measures: an adaptive estimation of information storage with application to brain-heart interactions

Network Physiology is a rapidly growing field of study that aims to understand how physiological systems interact to maintain health. Within the information theory framework the information storage (IS) allows to measure the regularity and predictability of a dynamic process under stationarity assumption. However, this assumption does not allow to track over time the transient pathways occurring in the dynamical activity of a physiological system. To address this limitation, we propose a time-varying approach based on the recursive least squares algorithm (RLS) for estimating IS at each time instant, in non-stationary conditions. We tested this approach in simulated time-varying dynamics and in the analysis of electroencephalographic (EEG) signals recorded from healthy volunteers and timed with the heartbeat to investigate brain-heart interactions. In simulations, we show that the proposed approach allows to track both abrupt and slow changes in the information stored in a physiological system. These changes are reflected in its evolution and variability over time. The analysis of brain-heart interactions reveals marked differences across the cardiac cycle phases of the variability of the time-varying IS. On the other hand, the average IS values exhibit a weak modulation over parieto-occiptal areas of the scalp. Our study highlights the importance of developing more advanced methods for measuring IS that account for non-stationarity in physiological systems. The proposed time-varying approach based on RLS represents a useful tool for identifying spatio-temporal dynamics within the neurocardiac system and can contribute to the understanding of brain-heart interactions.

Changes in the autonomic cardiorespiratory activity in parturient women with severe and moderate features of preeclampsia

Background

Cardiorespiratory coupling (CRC) is a physiological phenomenon that reflects the mutual interaction between the cardiac and respiratory control systems. It is mainly associated with efferent vagal activity from the central autonomic network. Few studies have explored the autonomic changes of CRC in preeclampsia, a critical obstetric complication related to possible autonomic dysfunctions and inflammatory disturbances. This study examined the autonomic mechanisms of CRC in women with severe and moderate preeclampsia and healthy controls by applying nonlinear methods based on information theory, such as mutual information (MI) and Renyi's mutual information (RMI) and the linear and nonlinear analysis of the Pulse-Respiration Quotient (PRQ).

Methods

We studied three groups of parturient women in the third trimester of pregnancy with a clinical diagnosis of preeclampsia without severe symptoms (P, 38.5 ± 1.4 weeks of pregnancy, n=19), preeclampsia with severe symptoms (SP, 37.5 ± 0.9 weeks of pregnancy, n=22), and normotensive control women (C, 39.1 ± 1.3 weeks of pregnancy, n=20). 10-minutes of abdominal electrocardiograms (ECG) and respiratory signals (RESP) were recorded in all the participants. Subsequently, we obtained the maternal beat-to-beat (RR) and breath-to-breath (BB) time series from ECG and RESP, respectively. The CRC between RR and BB was quantified by nonlinear methods based on information theory, such as MI and RMI, along with the analysis of the novel index of PRQ. Subsequently, we computed the mean PRQ (mPRQ) and the normalized permutation entropy (nPermEn_PRQ) from the PRQ time series generated from BB and RR. In addition, we examined the vagal activity in the three groups by the logarithm of the median of the distribution of the absolute values of successive RR differences (logRSA).

Results

The MI and RMI values were significantly lower (p<0.05) in the preeclamptic groups compared to the control group. However, no significant differences were found between the preeclamptic groups. The logRSA and nPermEn_PRQ indices were significantly lower (p<0.05) in SP compared to C and P.

Conclusion

Our data suggest that parturient women with severe and mild preeclampsia may manifest an altered cardiorespiratory coupling compared with normotensive control women. Disrupted CRC in severe preeclampsia could be associated with vagal withdrawal and less complex cardiorespiratory dynamics. The difference in vagal activity between the preeclamptic groups may suggest a further reduction in vagal activity associated with the severity of the disease.

Information-Theoretic Analysis of Cardio-Respiratory Interactions in Heart Failure Patients: Effects of Arrhythmias and Cardiac Resynchronization Therapy

The properties of cardio-respiratory coupling (CRC) are affected by various pathological conditions related to the cardiovascular and/or respiratory systems. In heart failure, one of the most common cardiac pathological conditions, the degree of CRC changes primarily depend on the type of heart-rhythm alterations. In this work, we investigated CRC in heart-failure patients, applying measures from information theory, i.e., Granger Causality (GC), Transfer Entropy (TE) and Cross Entropy (CE), to quantify the directed coupling and causality between cardiac (RR interval) and respiratory (Resp) time series. Patients were divided into three groups depending on their heart rhythm (sinus rhythm and presence of low/high number of ventricular extrasystoles) and were studied also after cardiac resynchronization therapy (CRT), distinguishing responders and non-responders to the therapy. The information-theoretic analysis of bidirectional cardio-respiratory interactions in HF patients revealed the strong effect of nonlinear components in the RR (high number of ventricular extrasystoles) and in the Resp time series (respiratory sinus arrhythmia) as well as in their causal interactions. We showed that GC as a linear model measure is not sensitive to both nonlinear components and only model free measures as TE and CE may quantify them. CRT responders mainly exhibit unchanged asymmetry in the TE values, with statistically significant dominance of the information flow from Resp to RR over the opposite flow from RR to Resp, before and after CRT. In non-responders this asymmetry was statistically significant only after CRT. Our results indicate that the success of CRT is related to corresponding information transfer between the cardiac and respiratory signal quantified at baseline measurements, which could contribute to a better selection of patients for this type of therapy.

Comparison of discretization strategies for the model-free information-theoretic assessment of short-term physiological interactions

This work presents a comparison between different approaches for the model-free estimation of information-theoretic measures of the dynamic coupling between short realizations of random processes. The measures considered are the mutual information rate (MIR) between two random processes X and Y and the terms of its decomposition evidencing either the individual entropy rates of X and Y and their joint entropy rate, or the transfer entropies from X to Y and from Y to X and the instantaneous information shared by X and Y. All measures are estimated through discretization of the random variables forming the processes, performed either via uniform quantization (binning approach) or rank ordering (permutation approach). The binning and permutation approaches are compared on simulations of two coupled non-identical Hènon systems and on three datasets, including short realizations of cardiorespiratory (CR, heart period and respiration flow), cardiovascular (CV, heart period and systolic arterial pressure), and cerebrovascular (CB, mean arterial pressure and cerebral blood flow velocity) measured in different physiological conditions, i.e., spontaneous vs paced breathing or supine vs upright positions. Our results show that, with careful selection of the estimation parameters (i.e., the embedding dimension and the number of quantization levels for the binning approach), meaningful patterns of the MIR and of its components can be achieved in the analyzed systems. On physiological time series, we found that paced breathing at slow breathing rates induces less complex and more coupled CR dynamics, while postural stress leads to unbalancing of CV interactions with prevalent baroreflex coupling and to less complex pressure dynamics with preserved CB interactions. These results are better highlighted by the permutation approach, thanks to its more parsimonious representation of the discretized dynamic patterns, which allows one to explore interactions with longer memory while limiting the curse of dimensionality.

Assessment of Cardiorespiratory Interactions during Apneic Events in Sleep via Fuzzy Kernel Measures of Information Dynamics

Apnea and other breathing-related disorders have been linked to the development of hypertension or impairments of the cardiovascular, cognitive or metabolic systems. The combined assessment of multiple physiological signals acquired during sleep is of fundamental importance for providing additional insights about breathing disorder events and the associated impairments. In this work, we apply information-theoretic measures to describe the joint dynamics of cardiorespiratory physiological processes in a large group of patients reporting repeated episodes of hypopneas, apneas (central, obstructive, mixed) and respiratory effort related arousals (RERAs). We analyze the heart period as the target process and the airflow amplitude as the driver, computing the predictive information, the information storage, the information transfer, the internal information and the cross information, using a fuzzy kernel entropy estimator. The analyses were performed comparing the information measures among segments during, immediately before and after the respiratory event and with control segments. Results highlight a general tendency to decrease of predictive information and information storage of heart period, as well as of cross information and information transfer from respiration to heart period, during the breathing disordered events. The information-theoretic measures also vary according to the breathing disorder, and significant changes of information transfer can be detected during RERAs, suggesting that the latter could represent a risk factor for developing cardiovascular diseases. These findings reflect the impact of different sleep breathing disorders on respiratory sinus arrhythmia, suggesting overall higher complexity of the cardiac dynamics and weaker cardiorespiratory interactions which may have physiological and clinical relevance.

Information Transfer in Linear Multivariate Processes Assessed through Penalized Regression Techniques: Validation and Application to Physiological Networks

The framework of information dynamics allows the dissection of the information processed in a network of multiple interacting dynamical systems into meaningful elements of computation that quantify the information generated in a target system, stored in it, transferred to it from one or more source systems, and modified in a synergistic or redundant way. The concepts of information transfer and modification have been recently formulated in the context of linear parametric modeling of vector stochastic processes, linking them to the notion of Granger causality and providing efficient tools for their computation based on the state-space (SS) representation of vector autoregressive (VAR) models. Despite their high computational reliability these tools still suffer from estimation problems which emerge, in the case of low ratio between data points available and the number of time series, when VAR identification is performed via the standard ordinary least squares (OLS). In this work we propose to replace the OLS with penalized regression performed through the Least Absolute Shrinkage and Selection Operator (LASSO), prior to computation of the measures of information transfer and information modification. First, simulating networks of several coupled Gaussian systems with complex interactions, we show that the LASSO regression allows, also in conditions of data paucity, to accurately reconstruct both the underlying network topology and the expected patterns of information transfer. Then we apply the proposed VAR-SS-LASSO approach to a challenging application context, i.e., the study of the physiological network of brain and peripheral interactions probed in humans under different conditions of rest and mental stress. Our results, which document the possibility to extract physiologically plausible patterns of interaction between the cardiovascular, respiratory and brain wave amplitudes, open the way to the use of our new analysis tools to explore the emerging field of Network Physiology in several practical applications.

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.

A Transfer Entropy Approach for the Assessment of the Impact of Inspiratory Muscle Training on the Cardiorespiratory Coupling of Amateur Cyclists

The strength of cardiorespiratory interactions diminishes with age. Physical exercise can reduce the rate of this trend. Inspiratory muscle training (IMT) is a technique capable of improving cardiorespiratory interactions. This study evaluates the effect of IMT on cardiorespiratory coupling in amateur cyclists. Thirty male young healthy cyclists underwent a sham IMT of very low intensity (SHAM, n = 9), an IMT of moderate intensity at 60% of the maximal inspiratory pressure (MIP60, n = 10) and an IMT of high intensity at the critical inspiratory pressure (CIP, n = 11). Electrocardiogram, non-invasive arterial pressure, and thoracic respiratory movement (RM) were recorded before (PRE) and after (POST) training at rest in supine position (REST) and during active standing (STAND). The beat-to-beat series of heart period (HP) and systolic arterial pressure (SAP) were analyzed with the RM signal via a traditional non-causal approach, such as squared coherence function, and via a causal model-based transfer entropy (TE) approach. Cardiorespiratory coupling was quantified via the HP-RM squared coherence at the respiratory rate (K 2 HP-R M), the unconditioned TE from RM to HP (TER M → HP) and the TE from RM to HP conditioned on SAP (TER M → HP| SAP). In PRE condition we found that STAND led to a decrease of TER M → HP| SAP. After SHAM and CIP training this tendency was confirmed, while MIP60 inverted it by empowering cardiorespiratory coupling. This behavior was observed in presence of unvaried SAP mean and with usual responses of the baroreflex control and HP mean to STAND. TER M → HP and K 2 HP- RM were not able to detect the post-training increase of cardiorespiratory coupling strength during STAND, thus suggesting that conditioning out SAP is important for the assessment of cardiorespiratory interactions. Since the usual response of HP mean, SAP mean and baroreflex sensitivity to postural stressor were observed after MIP60 training, we conclude that the post-training increase of cardiorespiratory coupling during STAND in MIP60 group might be the genuine effect of some rearrangements at the level of central respiratory network and its interactions with sympathetic drive and vagal activity.

Dynamical disentanglement in an analysis of oscillatory systems: an application to respiratory sinus arrhythmia

We develop a technique for the multivariate data analysis of perturbed self-sustained oscillators. The approach is based on the reconstruction of the phase dynamics model from observations and on a subsequent exploration of this model. For the system, driven by several inputs, we suggest a dynamical disentanglement procedure, allowing us to reconstruct the variability of the system's output that is due to a particular observed input, or, alternatively, to reconstruct the variability which is caused by all the inputs except for the observed one. We focus on the application of the method to the vagal component of the heart rate variability caused by a respiratory influence. We develop an algorithm that extracts purely respiratory-related variability, using a respiratory trace and times of R-peaks in the electrocardiogram. The algorithm can be applied to other systems where the observed bivariate data can be represented as a point process and a slow continuous signal, e.g. for the analysis of neuronal spiking. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

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.

Cardiorespiratory Coupling Analysis Based on Entropy and Cross-Entropy in Distinguishing Different Depression Stages

Aims

This study used entropy- and cross entropy-based methods to explore the cardiorespiratory coupling of depressive patients, and thus to assess the values of those entropy methods for identifying depression patients with different disease severities.

Methods

Electrocardiogram (ECG) and respiration signals from 69 depression patients were recorded simultaneously for 5 min. Patients were classified into three groups according to the Hamilton Depression Rating Scale (HDRS) scores: group Non-De (HDRS 0–7), Mid-De (HDRS 8–17), and Con-De (HDRS >17). Sample entropy (SEn), fuzzy measure entropy (FMEn) and high-frequency power (HF) were computed on the original RR interval time series and breath-to-breath interval time series. Cross sample entropy (CSEn) and cross fuzzy measure entropy (CFMEn) were computed on interval time series resampled at 2 Hz and 4 Hz, respectively. The difference among three patient groups and correlation between entropy values and HDRS scores were analyzed by statistical analysis. Surrogate data were also employed to confirm the validation of entropy measures in this study.

Results

A consistent increasing trend has been found among most entropy measures from Non-De, to Mid-De, and to Con-De groups, and a significant (p < 0.05) difference in FMEn of RR intervals exists between Non-De and Mid-De or Con-De groups. Significant differences have been also found in all cross entropies, between Non-De and Con-De groups and between Mid-De and Con-De groups. Furthermore, significant correlations also exist between HDRS scores and FMEn of RR intervals (R = 0.24, p < 0.05), CSEn at 4 Hz (R = 0.26, p < 0.05) or 2 Hz (R = 0.28, p < 0.05) resampling, and CFMEn at 4 Hz (R = 0.31, p < 0.01) or 2 Hz (R = 0.30, p < 0.05) resampling. A significant difference of cardiorespiratory coupling parameters between different depression stages and significant correlations between entropy measures and depression severity both indicate central autonomic dysregulation in depression patients and reflect varying degrees of vagal modulation reduction among different depression levels. Analysis based on surrogate data confirms that the non-linear properties of the physiological signals played a major role in depression recognition.

Conclusion

The current study demonstrates the potential of cardiorespiratory coupling in the auxiliary diagnosis of depression based on the entropy method.

Information Dynamics of the Brain, Cardiovascular and Respiratory Network during Different Levels of Mental Stress

In this study, an analysis of brain, cardiovascular and respiratory dynamics was conducted combining information-theoretic measures with the Network Physiology paradigm during different levels of mental stress. Starting from low invasive recordings of electroencephalographic, electrocardiographic, respiratory, and blood volume pulse signals, the dynamical activity of seven physiological systems was probed with one-second time resolution measuring the time series of the δ, θ, α and β brain wave amplitudes, the cardiac period (RR interval), the respiratory amplitude, and the duration of blood pressure wave propagation (pulse arrival time, PAT). Synchronous 5-min windows of these time series, obtained from 18 subjects during resting wakefulness (REST), mental stress induced by mental arithmetic (MA) and sustained attention induced by serious game (SG), were taken to describe the dynamics of the nodes composing the observed physiological network. Network activity and connectivity were then assessed in the framework of information dynamics computing the new information generated by each node, the information dynamically stored in it, and the information transferred to it from the other network nodes. Moreover, the network topology was investigated using directed measures of conditional information transfer and assessing their statistical significance. We found that all network nodes dynamically produce and store significant amounts of information, with the new information being prevalent in the brain systems and the information storage being prevalent in the peripheral systems. The transition from REST to MA was associated with an increase of the new information produced by the respiratory signal time series (RESP), and that from MA to SG with a decrease of the new information produced by PAT. Each network node received a significant amount of information from the other nodes, with the highest amount transferred to RR and the lowest transferred to δ, θ, α and β. The topology of the physiological network underlying such information transfer was node- and state-dependent, with the peripheral subnetwork showing interactions from RR to PAT and between RESP and RR, PAT consistently across states, the brain subnetwork resulting more connected during MA, and the subnetwork of brain–peripheral interactions involving different brain rhythms in the three states and resulting primarily activated during MA. These results have both physiological relevance as regards the interpretation of central and autonomic effects on cardiovascular and respiratory variability, and practical relevance as regards the identification of features useful for the automatic distinction of different mental states.

Multiscale information storage of linear long-range correlated stochastic processes

Information storage, reflecting the capability of a dynamical system to keep predictable information during its evolution over time, is a key element of intrinsic distributed computation, useful for the description of the dynamical complexity of several physical and biological processes. Here we introduce a parametric approach which allows one to compute information storage across multiple timescales in stochastic processes displaying both short-term dynamics and long-range correlations (LRC). Our analysis is performed in the popular framework of multiscale entropy, whereby a time series is first "coarse grained" at the chosen timescale through low-pass filtering and downsampling, and then its complexity is evaluated in terms of conditional entropy. Within this framework, our approach makes use of linear fractionally integrated autoregressive (ARFI) models to derive analytical expressions for the information storage computed at multiple timescales. Specifically, we exploit state space models to provide the representation of lowpass filtered and downsampled ARFI processes, from which information storage is computed at any given timescale relating the process variance to the prediction error variance. This enhances the practical usability of multiscale information storage, as it enables a computationally reliable quantification of a complexity measure which incorporates the effects of LRC together with that of short-term dynamics. The proposed measure is first assessed in simulated ARFI processes reproducing different types of autoregressive dynamics and different degrees of LRC, studying both the theoretical values and the finite sample performance. We find that LRC alter substantially the complexity of ARFI processes even at short timescales, and that reliable estimation of complexity can be achieved at longer timescales only when LRC are properly modeled. Then, we assess multiscale information storage in physiological time series measured in humans during resting state and postural stress, revealing unprecedented responses to stress of the complexity of heart period and systolic arterial pressure variability, which are related to the different role played by LRC in the two conditions.

Cardiorespiratory Temporal Causal Links and the Differences by Sport or Lack Thereof

Fitness level, fatigue and adaptation are important factors for determining the optimal training schedule and predicting future performance. We think that adding analysis of the mutual relationships between cardiac and respiratory activity enables better athlete profiling and feedback for improving training. Therefore, the main objectives were (1) to apply several methods for temporal causality analysis to cardiorespiratory data; (2) to establish causal links between the signals; and (3) to determine how parameterized connections differed across various subgroups. One hundred elite athletes (31 female) and a control group of 20 healthy students (6 female) took part in the study. All were asked to follow a protocol comprising two 5-min sessions of free breathing - once supine, once standing. The data were collected using Pneumonitor 2. Respiratory-related curves were obtained through impedance pneumography, along with a single-lead ECG. Several signals (e.g., tidal volume, instantaneous respiratory rate, and instantaneous heart rate) were derived and stored as: (1) raw data down-sampled to 25Hz; (2) further down-sampled to 2.5Hz; and (3) beat-by-beat sequences. Granger causality frameworks (pairwise-conditional, spectral or extended), along with Time Series Models with Independent Noise (TiMINo), were studied. The connections enabling the best distinctions were found using recursive feature elimination with a random forest kernel. Temporal causal links are the most evident between tidal volume and instantaneous heart rate signals. Predictions of the “effect” variable were improved by adding preceding “cause” samples, by medians of 20.3% for supine and 14.2% for standing body positions. Parameterized causal link structures and directions distinguish athletes from non-athletes with 83.3% accuracy on average. They may also be used to supplement standard analysis and enable classification into groups exhibiting different static and dynamic components during performance. Physiological markers of training may be extended to include cardiorespiratory data, and causality analysis may improve the resolution of training profiling and the precision of outcome prediction.

Information Domain Analysis of Respiratory Sinus Arrhythmia Mechanisms

Ventilation related heart rate oscillations – respiratory sinus arrhythmia (RSA) – originate in human from several mechanisms. Two most important of them – the central mechanism (direct communication between respiratory and cardiomotor centers), and the peripheral mechanism (ventilation-associated blood pressure changes transferred to heart rate via baroreflex) have been described in previous studies. The major aim of this study was to compare the importance of these mechanisms in the generation of RSA non-invasively during various states by quantifying the strength of the directed interactions between heart rate, systolic blood pressure and respiratory volume signals. Seventy-eight healthy volunteers (32 male, age range: 16.02-25.77 years, median age: 18.57 years) participated in this study. The strength of mutual interconnections among the spontaneous beat-to-beat oscillations of systolic blood pressure (SBP), R-R interval (RR signal) and respiration (volume changes – RESP signal) was quantified during supine rest, orthostatic challenge (head-up tilt, HUT) and cognitive load (mental arithmetics, MA) using bivariate and trivariate measures of cardio-respiratory information transfer to separate baroreflex and nonbaroreflex (central) mechanisms. Our results indicate that both basic mechanisms take part in RSA generation in the intact cardiorespiratory control of human subjects. During orthostatic and mental challenges baroreflex based peripheral mechanism becomes more important.

Paced Breathing Increases the Redundancy of Cardiorespiratory Control in Healthy Individuals and Chronic Heart Failure Patients

Synergy and redundancy are concepts that suggest, respectively, adaptability and fault tolerance of systems with complex behavior. This study computes redundancy/synergy in bivariate systems formed by a target X and a driver Y according to the predictive information decomposition approach and partial information decomposition framework based on the minimal mutual information principle. The two approaches assess the redundancy/synergy of past of X and Y in reducing the uncertainty of the current state of X. The methods were applied to evaluate the interactions between heart and respiration in healthy young subjects (n = 19) during controlled breathing at 10, 15 and 20 breaths/minute and in two groups of chronic heart failure patients during paced respiration at 6 (n = 9) and 15 (n = 20) breaths/minutes from spontaneous beat-to-beat fluctuations of heart period and respiratory signal. Both methods suggested that slowing respiratory rate below the spontaneous frequency increases redundancy of cardiorespiratory control in both healthy and pathological groups, thus possibly improving fault tolerance of the cardiorespiratory control. The two methods provide markers complementary to respiratory sinus arrhythmia and the strength of the linear coupling between heart period variability and respiration in describing the physiology of the cardiorespiratory reflex suitable to be exploited in various pathophysiological settings.

Assessment of Quadratic Nonlinear Cardiorespiratory Couplings During Tilt-Table Test by Means of Real Wavelet Biphase

OBJECTIVE

In this paper, a method for assessment of quadratic phase coupling (QPC) between respiration and heart rate variability (HRV) is presented.

METHODS

First, a method for QPC detection is proposed named real wavelet biphase (RWB). Then, a method for QPC quantification is proposed based on the normalized wavelet biamplitude (NWB). A simulation study has been conducted to test the reliability of RWB to identify QPC, even in the presence of constant delays between interacting oscillations, and to discriminate it from quadratic phase uncoupling. Significant QPC was assessed based on surrogate data analysis. Then, quadratic cardiorespiratory couplings were studied during a tilt-table test protocol of 17 young healthy subjects.

RESULTS

Simulation study showed that RWB is able to detect even weak QPC with delays in the range of [Formula: see text] s, which are usual in the autonomic nervous system (ANS) control of heart rate. Results from the database revealed a significant reduction ([Formula: see text]) of NWB between respiration and both low and high frequencies of HRV in head-up tilt position compared to early supine.

CONCLUSION

The proposed technique detects and quantifies robustly QPC and is able to track the coupling between respiration and various HRV components during ANS changes.

SIGNIFICANCE

The proposed method can help to assess alternations of nonlinear cardiorespiratory interactions related to ANS dysfunction and physiological regulation of HRV in cardiovascular diseases.

Task-induced deactivation in diverse brain systems correlates with interindividual differences in distinct autonomic indices

Neuroimaging research has shown that different cognitive tasks induce relatively specific activation patterns, as well as less task-specific deactivation patterns. Here we examined whether individual differences in Autonomic Nervous System (ANS) activity during task performance correlate with the magnitude of task-induced deactivation. In an fMRI study, participants performed a continuous mental arithmetic task in a task/rest block design, while undergoing combined fMRI and heart/respiration rate acquisitions using photoplethysmograph and respiration belt. As expected, task performance increased heart-rate and reduced the RMSSD, a cardiac index related to vagal tone. Across participants, higher heart rate during task was linked to increased activation in fronto-parietal regions, as well as to stronger deactivation in ventromedial prefrontal regions. Respiration frequency during task was associated with similar patterns, but in different regions than those identified for heart-rate. Finally, in a large set of regions, almost exclusively limited to the Default Mode Network, lower RMSSD was associated with greater deactivation, and furthermore, the vast majority of these regions were task-deactivated at the group level. Together, our findings show that inter-individual differences in ANS activity are strongly linked to task-induced deactivation. Importantly, our findings suggest that deactivation is a multifaceted construct potentially linked to ANS control, because distinct ANS measures correlate with deactivation in different regions. We discuss the implications for current theories of cortical control of the ANS and for accounts of deactivation, with particular reference to studies documenting a "failure to deactivate" in multiple clinical states.

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.

Estimation of confidence limits for descriptive indexes derived from autoregressive analysis of time series: Methods and application to heart rate variability

The growing interest in personalized medicine requires making inferences from descriptive indexes estimated from individual recordings of physiological signals, with statistical analyses focused on individual differences between/within subjects, rather than comparing supposedly homogeneous cohorts. To this end, methods to compute confidence limits of individual estimates of descriptive indexes are needed. This study introduces numerical methods to compute such confidence limits and perform statistical comparisons between indexes derived from autoregressive (AR) modeling of individual time series. Analytical approaches are generally not viable, because the indexes are usually nonlinear functions of the AR parameters. We exploit Monte Carlo (MC) and Bootstrap (BS) methods to reproduce the sampling distribution of the AR parameters and indexes computed from them. Here, these methods are implemented for spectral and information-theoretic indexes of heart-rate variability (HRV) estimated from AR models of heart-period time series. First, the MS and BC methods are tested in a wide range of synthetic HRV time series, showing good agreement with a gold-standard approach (i.e. multiple realizations of the "true" process driving the simulation). Then, real HRV time series measured from volunteers performing cognitive tasks are considered, documenting (i) the strong variability of confidence limits' width across recordings, (ii) the diversity of individual responses to the same task, and (iii) frequent disagreement between the cohort-average response and that of many individuals. We conclude that MC and BS methods are robust in estimating confidence limits of these AR-based indexes and thus recommended for short-term HRV analysis. Moreover, the strong inter-individual differences in the response to tasks shown by AR-based indexes evidence the need of individual-by-individual assessments of HRV features. Given their generality, MC and BS methods are promising for applications in biomedical signal processing and beyond, providing a powerful new tool for assessing the confidence limits of indexes estimated from individual recordings.

Preejection Period as a Sympathetic Activity Index: a Role of Confounding Factors

In previous studies, one of the systolic time intervals – preejection period (PEP) – was used as an index of sympathetic activity reflecting the cardiac contractility. However, PEP could be also influenced by several other cardiovascular variables including preload, afterload and diastolic blood pressure (DBP). The aim of this study was to assess the behavior of the PEP together with other potentially confounding cardiovascular system characteristics in healthy humans during mental and orthostatic stress (head-up tilt test – HUT). Forty-nine healthy volunteers (28 females, 21 males, mean age 18.6 years (SD=1.8 years)) participated in the study. We recorded finger arterial blood pressure by volume-clamp method (Finometer Pro, FMS, Netherlands), PEP, thoracic fluid content (TFC) – a measure of preload, and cardiac output (CO) by impedance cardiography (CardioScreen® 2000, Medis, Germany). Systemic vascular resistance (SVR) – a measure of afterload – was calculated as a ratio of mean arterial pressure and CO. We observed that during HUT, an expected decrease in TFC was accompanied by an increase of PEP, an increase of SVR and no significant change in DBP. During mental stress, we observed a decrease of PEP and an increase of TFC, SVR and DBP. Correlating a change in assessed measures (delta values) between mental stress and previous supine rest, we found that ΔPEP correlated negatively with ΔCO and positively with ΔSVR. In orthostasis, no significant correlation between ΔPEP and ΔDBP, ΔTFC, ΔCO, ΔMBP or ΔSVR was found. We conclude that despite an expected increase of sympathetic activity during both challenges, PEP behaved differently indicating an effect of other confounding factors. To interpret PEP values properly, we recommend simultaneously to measure other variables influencing this cardiovascular measure.

Cardiorespiratory Coordination in Repeated Maximal Exercise

Increases in cardiorespiratory coordination (CRC) after training with no differences in performance and physiological variables have recently been reported using a principal component analysis approach. However, no research has yet evaluated the short-term effects of exercise on CRC. The aim of this study was to delineate the behavior of CRC under different physiological initial conditions produced by repeated maximal exercises. Fifteen participants performed 2 consecutive graded and maximal cycling tests. Test 1 was performed without any previous exercise, and Test 2 6 min after Test 1. Both tests started at 0 W and the workload was increased by 25 W/min in males and 20 W/min in females, until they were not able to maintain the prescribed cycling frequency of 70 rpm for more than 5 consecutive seconds. A principal component (PC) analysis of selected cardiovascular and cardiorespiratory variables (expired fraction of O₂, expired fraction of CO₂, ventilation, systolic blood pressure, diastolic blood pressure, and heart rate) was performed to evaluate the CRC defined by the number of PCs in both tests. In order to quantify the degree of coordination, the information entropy was calculated and the eigenvalues of the first PC (PC1) were compared between tests. Although no significant differences were found between the tests with respect to the performed maximal workload (Wmax), maximal oxygen consumption (VO₂ max), or ventilatory threshold (VT), an increase in the number of PCs and/or a decrease of eigenvalues of PC₁ (t = 2.95; p = 0.01; d = 1.08) was found in Test 2 compared to Test 1. Moreover, entropy was significantly higher (Z = 2.33; p = 0.02; d = 1.43) in the last test. In conclusion, despite the fact that no significant differences were observed in the conventionally explored maximal performance and physiological variables (Wmax, VO₂ max, and VT) between tests, a reduction of CRC was observed in Test 2. These results emphasize the interest of CRC evaluation in the assessment and interpretation of cardiorespiratory exercise testing.

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.