Cardiovascular and cardiorespiratory coupling analyses: a review

Quantifying the Suitability of Biosignals Acquired During Surgery for Multimodal Analysis

Goal: Recently, large datasets of biosignals acquired during surgery became available. As they offer multiple physiological signals measured in parallel, multimodal analysis - which involves their joint analysis - can be conducted and could provide deeper insights than unimodal analysis based on a single signal. However, it is unclear what percentage of intraoperatively acquired data is suitable for multimodal analysis. Due to the large amount of data, manual inspection and labelling into suitable and unsuitable segments are not feasible. Nevertheless, multimodal analysis is performed successfully in sleep studies since many years as their signals have proven suitable. Hence, this study evaluates the suitability to perform multimodal analysis on a surgery dataset (VitalDB) using a multi-center sleep dataset (SIESTA) as reference. Methods: We applied widely known algorithms entitled "signal quality indicators" to the common biosignals in both datasets, namely electrocardiography, electroencephalography, and respiratory signals split in segments of 10 s duration. As there are no multimodal methods available, we used only unimodal signal quality indicators. In case, all three signals were determined as being adequate by the indicators, we assumed that the whole signal segment was suitable for multimodal analysis. Results: 82% of SIESTA and 72% of VitalDB are suitable for multimodal analysis. Unsuitable signal segments exhibit constant or physiologically unreasonable values. Histogram examination indicated similar signal quality distributions between the datasets, albeit with potential statistical biases due to different measurement setups. Conclusions: The majority of data within VitalDB is suitable for multimodal analysis.

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.

Long Sudarshan Kriya Yoga enhances cardiovascular and respiratory synchronization: An observational study

BACKGROUND

Sudarshan Kriya Yoga (SKY - a rhythmic cyclic breathing) is known to produce several physiological changes in human body. Earlier it has been reported that SKY improves cardiovascular modulations, namely increase in heart rate variability.

OBJECTIVE

To observe the synchronization in oscillatory modulations in cardiac autonomic tone and ventilatory exchange during Long Sudarshan Kriya Yoga (LSKY). LSKY is a sequential combination of pranayama in ujjayi breath, bhastrika, and cyclic rhythmic breathing followed by yog-nidra.

METHODS

Regular LSKY practitioners from the Art of Living community with more than two years of experience participated in the study (n = 22; age 40.09 ± 12.68). The Electrocardiogram (ECG), respiration, oxygen saturation, and concentrations of oxygen and carbon-di-oxide from exhaled air were recorded before and during LSKY. The time domain parameters of heart rate variability (HRV) were calculated from ECG. All parameters were compared and correlated at each stage of LSKY.

RESULTS

Highly significant reciprocal correlation was found between HRV parameters and respiration rate during LSKY. Both O2 consumption and CO2 production increased significantly during three stages of pranayama and decreased towards the end of cyclic breathing. We also saw increased SPO2 simultaneously.

CONCLUSION

The autonomic parameters exhibited reciprocal response to respiratory rate and correlated well to the ventilatory parameters. Further during LSKY we observed enhanced synchronization. In conclusion the LSKY has a potential to influence cardiorespiratory parameters for improving the performance of both systems. LSKY - enhances oscillations in HRV that resets the autonomic system, indicative of better cardiac health and prepares body for better metabolic response. Such changes are capable of inducing resilience along with physiological, psychological relaxation and emotional well-being.

Bradycardia May Decrease Cardiorespiratory Coupling in Preterm Infants

Bradycardia, frequently observed in preterm infants, presents significant risks due to the immaturity of their autonomic nervous system (ANS) and respiratory systems. These infants may face cardiorespiratory events, leading to severe complications like hypoxemia and neurodevelopmental disorders. Although neonatal care has advanced, the influence of bradycardia on cardiorespiratory coupling (CRC) remains elusive. This exploratory study delves into CRC in preterm infants, emphasizing disparities between events with and without bradycardia. Using the Preterm Infant Cardio-Respiratory Signals (PICS) database, we analyzed interbeat (R-R) and inter-breath intervals (IBI) from 10 preterm infants. The time series were segmented into bradycardic (B) and non-bradycardic (NB) segments. Employing information theory measures, we quantified the irregularity of cardiac and respiratory time series. Notably, B segments had significantly lower entropy values for R-R and IBI than NB segments, while mutual information was higher in NB segments. This could imply a reduction in the complexity of respiratory and cardiac dynamics during bradycardic events, potentially indicating weaker CRC. Building on these insights, this research highlights the distinctive physiological characteristics of preterm infants and underscores the potential of emerging non-invasive diagnostic tools.

Cardiorespiratory dynamics during respiratory maneuver in athletes

Introduction: The modern practice of sports medicine and medical rehabilitation requires the search for subtle criteria for the development of conditions and recovery of the body after diseases, which would have a prognostic value for the prevention of negative effects of training and rehabilitation tools, and also testify to the development and course of mechanisms for counteracting pathogenetic processes in the body. The purpose of this study was to determine the informative directions of the cardiorespiratory system parameters dynamics during the performing a maneuver with a change in breathing rate, which may indicate the body functional state violation.

Methods: The results of the study of 183 healthy men aged 21.2 ± 2.3 years who regularly engaged in various sports were analyzed. The procedure for studying the cardiorespiratory system included conducting combined measurements of indicators of activity of the respiratory and cardiovascular systems in a sitting position using a spiroarteriocardiograph device. The duration of the study was 6 min and involved the sequential registration of three measurements with a change in breathing rate (spontaneous breathing, breathing at 0.1 Hz and 0.25 Hz).

Results: Performing a breathing maneuver at breathing 0.1 Hz and breathing 0.25 Hz in comparison with spontaneous breathing leads to multidirectional significant changes in heart rate variability indicators–TP (ms2), LF (ms2), LFHF (ms2/ms2); of blood pressure variability indicators–TPDBP (mmHg2), LFSBP (mmHg2), LFDBP (mmHg2), HFSBP (mmHg2); of volume respiration variability indicators - LFR, (L×min-1)2; HFR, (L×min-1)2; LFHFR, (L×min-1)2/(L×min-1)2; of arterial baroreflex sensitivity indicators - BRLF (ms×mmHg-1), BRHF (ms×mmHg-1). Differences in indicators of systemic hemodynamics and indicators of cardiovascular and respiratory systems synchronization were also informative.

Conclusion: According to the results of the study, it is shown that during performing a breathing maneuver with a change in the rate of breathing, there are significant changes in cardiorespiratory parameters, the analysis of which the increments made it possible to determine of the changes directions dynamics, their absolute values and informative limits regarding the possible occurrence of the cardiorespiratory interactions dysregulation.

Cardiorespiratory dynamics during respiratory maneuver in athletes

Introduction: The modern practice of sports medicine and medical rehabilitation requires the search for subtle criteria for the development of conditions and recovery of the body after diseases, which would have a prognostic value for the prevention of negative effects of training and rehabilitation tools, and also testify to the development and course of mechanisms for counteracting pathogenetic processes in the body. The purpose of this study was to determine the informative directions of the cardiorespiratory system parameters dynamics during the performing a maneuver with a change in breathing rate, which may indicate the body functional state violation. Methods: The results of the study of 183 healthy men aged 21.2 ± 2.3 years who regularly engaged in various sports were analyzed. The procedure for studying the cardiorespiratory system included conducting combined measurements of indicators of activity of the respiratory and cardiovascular systems in a sitting position using a spiroarteriocardiograph device. The duration of the study was 6 min and involved the sequential registration of three measurements with a change in breathing rate (spontaneous breathing, breathing at 0.1 Hz and 0.25 Hz). Results: Performing a breathing maneuver at breathing 0.1 Hz and breathing 0.25 Hz in comparison with spontaneous breathing leads to multidirectional significant changes in heart rate variability indicators-TP (ms2), LF (ms2), LFHF (ms2/ms2); of blood pressure variability indicators-TPDBP (mmHg2), LFSBP (mmHg2), LFDBP (mmHg2), HFSBP (mmHg2); of volume respiration variability indicators - LFR, (L×min-1)2; HFR, (L×min-1)2; LFHFR, (L×min-1)2/(L×min-1)2; of arterial baroreflex sensitivity indicators - BRLF (ms×mmHg-1), BRHF (ms×mmHg-1). Differences in indicators of systemic hemodynamics and indicators of cardiovascular and respiratory systems synchronization were also informative. Conclusion: According to the results of the study, it is shown that during performing a breathing maneuver with a change in the rate of breathing, there are significant changes in cardiorespiratory parameters, the analysis of which the increments made it possible to determine of the changes directions dynamics, their absolute values and informative limits regarding the possible occurrence of the cardiorespiratory interactions dysregulation.

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 cardiorespiratory coordination assessed by the photoplethysmography imaging technique

Cardiorespiratory coordination (CRC) probes the interaction between cardiac and respiratory oscillators in which cardiac and respiratory activity are synchronized, with individual heartbeats occurring at approximately the same temporal positions during several breathing cycles. An increase of CRC has previously been related to pathological stressful states. We studied CRC employing coordigrams computed from non-contact photoplethysmography imaging (PPGI) and respiratory data using the optical flow method. In a blocked study design, we applied the cold pressure test (CPT), water at ambient temperature (AWT), and intermittent resting conditions. In controls (no intervention), CRC remained on initial low levels throughout measurements. In the experimental group (AWT and CPT intervention), CRC decreased during AWT and CPT. Following both interventions, CRC increased significantly, with a rebound effect following AWT. In controls, HR increased steadily over time. CPT evoked a significant HR increase which correlated with subjective stress/pain ratings. The CRC increase following AWT correlated significantly with subjective pain (r = .79) and stress (r = .63) ratings. Furthermore, we observed a significant correlation (r = - .80) between mean RMSSD and mean duration of CRC, which further supports an association between autonomic state and CRC level. CRC analysis obtained from cutaneous tissue perfusion data therefore appears to be a sensitive and useful method for the study of CRC and ANS activity. Future studies need to investigate the physiological principles and clinical significance of these findings.

Classifier for the functional state of the respiratory system via descriptors determined by using multimodal technology

Currently, intelligent systems built on a multimodal basis are used to study the functional state of living objects. Its essence lies in the fact that a decision is made through several independent information channels with the subsequent aggregation of these decisions. The method of forming descriptors for classifiers of the functional state of the respiratory system includes the study of the spectral range of the respiratory rhythm and the construction of the wavelet plane of the monitoring electrocardiosignal overlapping this range. Then, variations in the breathing rhythm are determined along the corresponding lines of the wavelet plane. Its analysis makes it possible to select slow waves corresponding to the breathing rhythm and systemic waves of the second order. Analysis of the spectral characteristics of these waves makes it possible to form a space of informative features for classifiers of the functional state of the respiratory system. To construct classifiers of the functional state of the respiratory system, hierarchical classifiers were used. As an example, we took a group of patients with pneumonia with a well-defined diagnosis (radiography, X-ray tomography, laboratory data) and a group of volunteers without pulmonary pathology. The diagnostic sensitivity of the obtained classifier was 76% specificity with a diagnostic specificity of 82%, which is comparable to the results of X-ray studies. It is shown that the corresponding lines of the wavelet planes are correlated with the respiratory system and, using their Fourier analysis, descriptors can be obtained for training neural network classifiers of the functional state of the respiratory system.

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.

Assessing cardiorespiratory interactions via lagged joint symbolic dynamics during spontaneous and controlled breathing

Introduction: Joint symbolic analysis (JSA) can be utilized to describe interactions between time series while accounting for time scales and nonlinear features. JSA is based on the computation of the rate of occurrence of joint patterns built after symbolization. Lagged JSA (LJSA) is obtained from the more classical JSA by introducing a delay/lead between patterns built over the two series and combined to form the joint scheme, thus monitoring coordinated patterns at different lags. Methods: In the present study, we applied LJSA for the assessment of cardiorespiratory coupling (CRC) from heart period (HP) variability and respiratory activity (R) in 19 healthy subjects (age: 27-35 years; 8 males, 11 females) during spontaneous breathing (SB) and controlled breathing (CB). The R rate of CB was selected to be indistinguishable from that of SB, namely, 15 breaths·minute-1 (CB15), or slower than SB, namely, 10 breaths·minute-1 (CB10), but in both cases, very rapid interactions between heart rate and R were known to be present. The ability of the LJSA approach to follow variations of the coupling strength was tested over a unidirectionally or bidirectionally coupled stochastic process and using surrogate data to test the null hypothesis of uncoupling. Results: We found that: i) the analysis of surrogate data proved that HP and R were significantly coupled in any experimental condition, and coupling was not more likely to occur at a specific time lag; ii) CB10 reduced CRC strength at the fastest time scales while increasing that at intermediate time scales, thus leaving the overall CRC strength unvaried; iii) despite exhibiting similar R rates and respiratory sinus arrhythmia, SB and CB15 induced different cardiorespiratory interactions; iv) no dominant temporal scheme was observed with relevant contributions of HP patterns either leading or lagging R. Discussion: LJSA is a useful methodology to explore HP-R dynamic interactions while accounting for time shifts and scales.

Decomposing the Mutual Information Rate of Heart Period and Respiration Variability Series to Assess Cardiorespiratory Interactions

Heart rate variability results from the coupled activity of the cardiovascular and cardiorespiratory systems, which have their own internal regulation mechanisms but also interact with each other and with the autonomic nervous system to maintain homeostasis. In this work, the assessment of these physiological mechanisms is carried out decomposing the Mutual Information Rate (MIR), an information-theoretic measure of the interdependence between coupled processes, into terms of entropy rate or conditional mutual information related respectively to complexity and causality measures. These measures are computed using a non-parametric approach based on nearest-neighbors. The proposed framework is first tested on simulated autoregressive processes and then applied to experimental data consisting of heart period and respiratory time series measured in healthy subjects monitored at rest and during head-up tilt. Our results evidence that MIR decomposition is able to highlight the interdependence of short-term physiological mechanisms of cardiorespiratory interactions during postural stress.

Directional coupling detection through cross-distance vectors

Inferring the coupling direction from measured time series of complex systems is challenging. We propose a state-space-based causality measure obtained from cross-distance vectors for quantifying interaction strength. It is a model-free noise-robust approach that requires only a few parameters. The approach is applicable to bivariate time series and is resilient to artefacts and missing values. The result is two coupling indices that quantify coupling strength in each direction more accurately than the already established state-space measures. We test the proposed method on different dynamical systems and analyze numerical stability. As a result, a procedure for optimal parameter selection is proposed, circumventing the challenge of determining the optimal embedding parameters. We show it is robust to noise and reliable in shorter time series. Moreover, we show that it can detect cardiorespiratory interaction in measured data. A numerically efficient implementation is available at https://repo.ijs.si/e2pub/cd-vec.

Analysis of Respiratory Sinus Arrhythmia and Directed Information Flow between Brain and Body Indicate Different Management Strategies of fMRI-Related Anxiety

Background: Respiratory sinus arrhythmia (RSA) denotes decrease of cardiac beat-to-beat intervals (RRI) during inspiration and RRI increase during expiration, but an inverse pattern (termed negative RSA) was also found in healthy humans with elevated anxiety. It was detected using wave-by-wave analysis of cardiorespiratory rhythms and was considered to reflect a strategy of anxiety management involving the activation of a neural pacemaker. Results were consistent with slow breathing, but contained uncertainty at normal breathing rates (0.2–0.4 Hz). Objectives and methods: We combined wave-by-wave analysis and directed information flow analysis to obtain information on anxiety management at higher breathing rates. We analyzed cardiorespiratory rhythms and blood oxygen level-dependent (BOLD) signals from the brainstem and cortex in 10 healthy fMRI participants with elevated anxiety. Results: Three subjects with slow respiratory, RRI, and neural BOLD oscillations showed 57 ± 26% negative RSA and significant anxiety reduction by 54 ± 9%. Six participants with breathing rate of ~0.3 Hz showed 41 ± 16% negative RSA and weaker anxiety reduction. They presented significant information flow from RRI to respiration and from the middle frontal cortex to the brainstem, which may result from respiration-entrained brain oscillations, indicating another anxiety management strategy. Conclusion: The two analytical approaches applied here indicate at least two different anxiety management strategies in healthy subjects.

Two Operational Modes of Cardio-Respiratory Coupling Revealed by Pulse-Respiration Quotient

Due to the fact that respiratory breath-to-breath and cardiac intervals between two successive R peaks (BBI and RRI, respectively) are not temporally concurrent, in a previous paper, we proposed a method to calculate both the integer and non-integer parts of the pulse respiration quotient (PRQ = BBI/RRI = PRQ_int + b1 + b2), b1 and b2 being parts of the border RRIs for each BBI. In this work, we study the correlations between BBI and PRQ, as well as those between BBI and mean RRI within each BBI (mRRI), on a group of twenty subjects in four conditions: in supine and standing positions, in combination with spontaneous and slow breathing. Results show that the BBI vs. PRQ correlations are positive; whereas the breathing regime had little or no effect on the linear regression slopes, body posture did. Two types of scatter plots were obtained with the BBI vs. mRRI correlations: one showed points aggregated around the concurrent PRQ_int lines, while the other showed randomly distributed points. Five out of six of the proposed aggregation measures confirmed the existence of these two cardio-respiratory coupling regimes. We also used b1 to study the positions of R pulses relative to the respiration onsets and showed that they were more synchronous with sympathetic activation. Overall, this method should be used in different pathological states.

On the significance of estimating cardiorespiratory coupling strength in sports medicine

The estimation of cardiorespiratory coupling (CRC) is attracting interest in sports physiology as an important tool to characterize cardiac neural regulation genuinely driven by respiration. When applied in sports medicine, cardiorespiratory coupling measurements can provide information on the effects of training, pre-competition stress, as well as cardiovascular adjustments during stressful stimuli. Furthermore, since the cardiorespiratory coupling is strongly affected by physical activity, the study of the cardiorespiratory coupling can guide the application of specific training methods to optimize the coupling between autonomic activity and heart with possible effects on performance. However, a consensus about the physiological mechanisms, as well as methodological gold standard methods to quantify the cardiorespiratory coupling, has not been reached yet, thus limiting its application in experimental settings. This review supports the relevance of assessing cardiorespiratory coupling in the sports medicine, examines the possible physiological mechanisms involved, and lists a series of methodological approaches. cardiorespiratory coupling strength seems to be increased in athletes when compared to sedentary subjects, in addition to being associated with positive physiological outcomes, such as a possible better interaction of neural subsystems to cope with stressful stimuli. Moreover, cardiorespiratory coupling seems to be influenced by specific training modalities, such as inspiratory muscle training. However, the impact of cardiorespiratory coupling on sports performance still needs to be better explored through ad hoc physical exercise tests and protocols. In addition, this review stresses that several bivariate and multivariate methods have been proposed to assess cardiorespiratory coupling, thus opening new possibilities in estimating cardiorespiratory interactions in athletes.

Confounding effects of heart rate, breathing rate, and frontal fNIRS on interoception

Recent studies have established that cardiac and respiratory phases can modulate perception and related neural dynamics. While heart rate and respiratory sinus arrhythmia possibly affect interoception biomarkers, such as heartbeat-evoked potentials, the relative changes in heart rate and cardiorespiratory dynamics in interoceptive processes have not yet been investigated. In this study, we investigated the variation in heart and breathing rates, as well as higher functional dynamics including cardiorespiratory correlation and frontal hemodynamics measured with fNIRS, during a heartbeat counting task. To further investigate the functional physiology linked to changes in vagal activity caused by specific breathing rates, we performed the heartbeat counting task together with a controlled breathing rate task. The results demonstrate that focusing on heartbeats decreases breathing and heart rates in comparison, which may be part of the physiological mechanisms related to "listening" to the heart, the focus of attention, and self-awareness. Focusing on heartbeats was also observed to increase frontal connectivity, supporting the role of frontal structures in the neural monitoring of visceral inputs. However, cardiorespiratory correlation is affected by both heartbeats counting and controlled breathing tasks. Based on these results, we concluded that variations in heart and breathing rates are confounding factors in the assessment of interoceptive abilities and relative fluctuations in breathing and heart rates should be considered to be a mode of covariate measurement of interoceptive processes.

Feasibility of Ultra-Short-Term Analysis of Heart Rate and Systolic Arterial Pressure Variability at Rest and during Stress via Time-Domain and Entropy-Based Measures

Heart Rate Variability (HRV) and Blood Pressure Variability (BPV) are widely employed tools for characterizing the complex behavior of cardiovascular dynamics. Usually, HRV and BPV analyses are carried out through short-term (ST) measurements, which exploit ~five-minute-long recordings. Recent research efforts are focused on reducing the time series length, assessing whether and to what extent Ultra-Short-Term (UST) analysis is capable of extracting information about cardiovascular variability from very short recordings. In this work, we compare ST and UST measures computed on electrocardiographic R-R intervals and systolic arterial pressure time series obtained at rest and during both postural and mental stress. Standard time-domain indices are computed, together with entropy-based measures able to assess the regularity and complexity of cardiovascular dynamics, on time series lasting down to 60 samples, employing either a faster linear parametric estimator or a more reliable but time-consuming model-free method based on nearest neighbor estimates. Our results are evidence that shorter time series down to 120 samples still exhibit an acceptable agreement with the ST reference and can also be exploited to discriminate between stress and rest. Moreover, despite neglecting nonlinearities inherent to short-term cardiovascular dynamics, the faster linear estimator is still capable of detecting differences among the conditions, thus resulting in its suitability to be implemented on wearable devices.

Cardiorespiratory Coordination in Collegiate Rowing: A Network Approach to Cardiorespiratory Exercise Testing

Collegiate rowing performance is often assessed by a cardiopulmonary exercise test (CPET). Rowers' on-water performance involves non-linear dynamic interactions and synergetic reconfigurations of the cardiorespiratory system. Cardiorespiratory coordination (CRC) method measures the co-variation among cardiorespiratory variables. Novice (n = 9) vs. Intermediate (n = 9) rowers' CRC (H₀: Novice CRC = Intermediate CRC; H_A: Novice CRC < Intermediate CRC) was evaluated through principal components analysis (PCA). A female NCAA Division II team (N = 18) grouped based on their off-water performance on 6000 m time trial. Rowers completed a customized CPET to exhaustion and a variety of cardiorespiratory values were recorded. The number of principal components (PCs) and respective PC eigenvalues per group were computed on SPSS vs28. Intermediate (77%) and Novice (33%) groups showed one PC₁. Novice group formed an added PC₂ due to the shift of expired fraction of oxygen or, alternatively, heart rate/ventilation, from the PC₁ cluster of examined variables. Intermediate rowers presented a higher degree of CRC, possible due to their increased ability to utilize the bicarbonate buffering system during the CPET. CRC may be an alternative measure to assess aerobic fitness providing insights to the complex cardiorespiratory interactions involved in rowing during a CPET.

Cardiorespiratory Phase Synchronization in Elderly Patients with Periodic and non-Periodic Breathing Patterns

Cardiorespiratory Phase Synchronization (CRPS) is the manifestation of the non-linear coupling between the cardiac and the respiratory systems, different to the Respiratory Sinus Arrythmia (RSA). This takes place when the heartbeats occur at the same relative phase of the breathing, during a succession of respiratory cycles. In this study, we investigated the CRPS in 45 elderly patients admitted to the semi-critical unit of a hospital. The patients were classified according to their respiratory state as non-Periodic Breathing (nPB), Periodic Breathing (PB) and Cheyne-Stokes Respiration (CSR). The phase synchrogram between the electrocardiographic and respiratory signals was computed using the Hilbert transform technique. A continuous measure of the CRPS was obtained from the synchrogram, and was characterized by the average duration of synchronized epochs (A vgDurSync), the percentage of synchronized time (%Sync), the number of synchronized epochs (NumSync), and the frequency ratio between the cardiac and respiratory oscillators (FreqRat). These measures were studied using two different thresholds (0.1 and 0.05) for the amplitude of the synchronization and a minimum duration threshold of 10s. According to the results, the AvgDurSync and %Sync had a decreasing trend in patients with breathing periodicity. In addition, CSR patients presented the lowest values A vgDurSync and %Sync. Therefore, the CRPS method could be a useful tool for characterizing periodic respiratory patterns in elderly patients, which might be related to chronic heart failure. Clinical Relevance- This study analyzes the synchronization between cardiac and respiratory systems in elderly patients with a possible progressive decompensation in the cardiac function.

Effects of exercise anticipation on cardiorespiratory coherence

In this study, we explored the role of feedforward mechanisms in triggering cardiorespiratory adjustments before the onset of exercise. To isolate the feedforward aspects, we examined the effect of exercise anticipation on cardiorespiratory coherence. Twenty-nine healthy males (age = 18.8 [0.96] years) were subjected to bicycle (BE) and handgrip exercise (H) at two different intensities, viz., low and high. Bicycle exercise was performed in a unilateral (left- and right-sided) or bilateral mode, whereas handgrip was performed only in a unilateral mode. Single-lead ECG and respiratory rhythm, measured in the 5 min of anticipation phase before the onset of exercise, were used for analysis. Coherence was computed between ECG-derived instantaneous heart rate and respiratory signal. Average coherence in the high-frequency band (0.15-0.4 Hz) was used to estimate respiratory sinus arrhythmia (RSA). We found that coherence decreased with the anticipation of exercise relative to baseline (baseline = 0.54 [0.16], BE = 0.41 [0.12], H = 0.39 [0.12], p < 0.001). The decrease was greater for high intensity exercise (low = 0.42 [0.11], high = 0.37 [0.1], p < 0.001). The fall of coherence with intensity was stronger for bicycle exercise (BE: low = 0.44 [0.12], high = 0.37 [0.12], H: low = 0.4 [0.12], high = 0.37 [0.12], p = 0.00433). The expectation of bilateral exercise resulted in lower coherence compared to unilateral exercise (right-sided = 0.45 [0.16], left-sided = 0.4 [0.16], bilateral = 0.36 [0.15], unilateral vs. bilateral: p < 0.001), and the left-sided exercise had lower coherence compared to that of the right (left-sided vs. right-sided: p = 0.00925). Handgrip exercise showed similar trend (right-sided = 0.4 [0.15], left-sided = 0.37 [0.14], p = 0.0056). In conclusion, feedforward RSA adjustments in anticipation of exercise covaried with subsequent exercise-related features like intensity, muscle mass (unilateral vs. bilateral), and the exercise side (left vs. right). The left versus the right difference in coherence indicates autonomic asymmetry. Feedforward changes in RSA are like those seen during actual exercise and might facilitate the rapid phase transition between rest and exercise.

Comparative Analysis of Cardiorespiratory Parameters of Basketball and Soccer Players Using Principal Component Analysis

Principal component analysis (PCA) is a statistical technique used to identify variations in multivariate data obtained during the performance of the maximum ergospirometry test (MET). To use the PCA to compare the coefficients of change of the principal component (PC1) using the eigenvalue and the maximum values ​​of the cardiorespiratory variables obtained in the athletes' in MET. 10 soccer players and 10 basketball players, all male, were evaluated. The PCA analyzed the values ​​of the variables during the performance of the MET. The PC1 for each variable was calculated, and the eigenvalue was generated, representing the coefficients of variation of the PC1 of all variables. In the quantitative assessment (maximum values), a higher VO2max (3.93±0.62 vs. 3.41±0.37 l·min-1) was observed in basketball players compared to soccer players (p<0.05). The qualitative evaluation using PC1 of cardiorespiratory parameters (heart rate, minute volume, O2 consumption, CO2 production, expired fraction of O2 and expired fraction CO2) was observed as an eigenvalue (6.50±0.27 vs. 6.22±0.19) high for basketball players compared to soccer players (p<0.05). It is concluded that the basketball players showed more significant variability in their cardiorespiratory variables during the performance of the MET and higher VO2max at the end of the MET. These findings indicate that basketball players were less efficient in buffering the ventilatory acidosis observed during the MET. The results of this study highlight the importance of making complex assessments of the cardiorespiratory system, providing qualitative information to complement the quantitative data.

The Reconstruction of Causal Networks in Physiology

We systematically compare strengths and weaknesses of two methods that can be used to quantify causal links between time series: Granger-causality and Bivariate Phase Rectified Signal Averaging (BPRSA). While a statistical test method for Granger-causality has already been established, we show that BPRSA causality can also be probed with existing statistical tests. Our results indicate that more data or stronger interactions are required for the BPRSA method than for the Granger-causality method to detect an existing link. Furthermore, the Granger-causality method can distinguish direct causal links from indirect links as well as links that arise from a common source, while BPRSA cannot. However, in contrast to Granger-causality, BPRSA is suited for the analysis of non-stationary data. We demonstrate the practicability of the Granger-causality method by applying it to polysomnography data from sleep laboratories. An algorithm is presented, which addresses the stationarity condition of Granger-causality by splitting non-stationary data into shorter segments until they pass a stationarity test. We reconstruct causal networks of heart rate, breathing rate, and EEG amplitude from young healthy subjects, elderly healthy subjects, and subjects with obstructive sleep apnea, a condition that leads to disruption of normal respiration during sleep. These networks exhibit differences not only between different sleep stages, but also between young and elderly healthy subjects on the one hand and subjects with sleep apnea on the other hand. Among these differences are 1) weaker interactions in all groups between heart rate, breathing rate and EEG amplitude during deep sleep, compared to light and REM sleep, 2) a stronger causal link from heart rate to breathing rate but disturbances in respiratory sinus arrhythmia (breathing to heart rate coupling) in subjects with sleep apnea, 3) a stronger causal link from EEG amplitude to breathing rate during REM sleep in subjects with sleep apnea. The Granger-causality method, although initially developed for econometric purposes, can provide a quantitative, testable measure for causality in physiological networks.

Directional couplings between the respiration and parasympathetic control of the heart rate during sleep and wakefulness in healthy subjects at different ages

Cardiorespiratory interactions are important, both for understanding the fundamental processes of functioning of the human body and for development of methods for diagnostics of various pathologies. The properties of cardiorespiratory interaction are determined by the processes of autonomic control of blood circulation, which are modulated by the higher nervous activity. We study the directional couplings between the respiration and the process of parasympathetic control of the heart rate in the awake state and different stages of sleep in 96 healthy subjects from different age groups. The detection of directional couplings is carried out using the method of phase dynamics modeling applied to experimental RR-intervals and the signal of respiration. We reveal the presence of bidirectional couplings between the studied processes in all age groups. Our results show that the coupling from respiration to the process of parasympathetic control of the heart rate is stronger than the coupling in the opposite direction. The difference in the strength of bidirectional couplings between the considered processes is most pronounced in deep sleep.

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'.

Assessment of long-range cross-correlations in cardiorespiratory and cardiovascular interactions

We propose higher-order detrending moving-average cross-correlation analysis (DMCA) to assess the long-range cross-correlations in cardiorespiratory and cardiovascular interactions. Although the original (zeroth-order) DMCA employs a simple moving-average detrending filter to remove non-stationary trends embedded in the observed time series, our approach incorporates a Savitzky-Golay filter as a higher-order detrending method. Because the non-stationary trends can adversely affect the long-range correlation assessment, the higher-order detrending serves to improve accuracy. To achieve a more reliable characterization of the long-range cross-correlations, we demonstrate the importance of the following steps: correcting the time scale, confirming the consistency of different order DMCAs, and estimating the time lag between time series. We applied this methodological framework to cardiorespiratory and cardiovascular time series analysis. In the cardiorespiratory interaction, respiratory and heart rate variability (HRV) showed long-range auto-correlations; however, no factor was shared between them. In the cardiovascular interaction, beat-to-beat systolic blood pressure and HRV showed long-range auto-correlations and shared a common long-range, cross-correlated factor. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

The Added Value of Nonlinear Cardiorespiratory Coupling Indices in the Assessment of Depression. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2021;2021:5473-5476. [PMID: 34892364 DOI: 10.1109/embc46164.2021.9631096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The present study investigates the differences in autonomic nervous system (ANS) function and stress response between patients with major depressive disorder (MDD) and healthy subjects by measuring changes in ANS biomarkers. ANS-related parameters are derived from various biosignals during a mental stress protocol consisting of a basal, stress, and recovery phase. The feature set consists of ANS biomarkers such as the heart rate (HR) derived from the electrocardiogram, the respiratory rate derived from the respiration signal, vascular parameters obtained from a model-based photoplethysmographic pulse waveform analysis, and cardiorespiratory coupling indices derived from the joint analysis of the heart rate variability (HRV) and respiratory signals. In particular, linear cardiorespiratory interactions are quantified by means of time-frequency coherence, while interactions of quadratic nonlinear nature between HRV and respiration are quantified by means of real wavelet biphase. The intra-subject difference of a feature value between two phases of the protocol, the so-called autonomic reactivity, is considered as a ANS biomarker as well. The performance of ANS biomarkers on discriminating MDD patients is evaluated using a classification pipeline. The results show that the most discriminative ANS biomarkers are related with differences in HR and autonomic reactivity of both vascular and nonlinear cardiorespiratory coupling indices. Differences in autonomic reactivity imply that MDD and healthy subjects differ in their ability to cope with stress. Considering only HR and vascular characteristics a linear support-vector machine classifier yields to accuracy 72.5% and F1-score 73.2%. However, taking into account the nonlinear cardiorespiratory coupling indices, the classification performance improves, yielding to accuracy 77.5% and F1-score 78.0%.Clinical relevance- Changes in the nonlinear properties of the cardiorespiratory system during stress may yield additional information on the assessment of depression.

Autonomic Function in Parkinson's Disease Subjects Across Repeated Short-Term Dry Immersion: Evidence From Linear and Non-linear HRV Parameters

Several studies have shown that “dry” immersion appears as a promising method of rehabilitation for Parkinson's disease. Still, little is known about the cardiovascular reaction in “dry” immersion (DI), especially in Parkinson's disease (PD). Therefore, this study was aimed to evaluate the effect of repeated 45-min DI sessions on autonomic function in subjects with PD. The study group consisted of 20 subjects with PD [13 men, seven women, aged 51–66 years old, Hoehn & Yahr (H&Y) staged 1–3] were enrolled in the study according to inclusion and non-inclusion criteria. The DI program was comprised of seven 45-min DI sessions, applied within 25–30 days. Blood pressure (BP), heart rate (HR), and electrocardiogram (ECG) in the standard lead II were recorded at 1st, 4th, and 7th DI, before, on the 15, 30, and 40th min of DI session. Autonomic function was assessed with analysis of heart rate variability (HRV) using Kubios Standard version 2 software. Linear (time- and frequency-domain) and non-linear (correlation dimension, entropies, DFA1 and DFA2, percent of determinism, and recurrence) were computed. At baseline condition, time- and frequency-domain HRV parameters showed low variability of HR, which indicates reduced autonomic neurogenic control of HR. Throughout the DI session, systolic and diastolic BP has decreased by 5–7 mm Hg (p < 0.001), and time- and frequency-domain parameters of HRV have significantly increased, what can be regarded as compensatory mechanisms of hemodynamics during DI. The structure of the regulatory input to the heart seen by HRV was characterized by low complexity and reduced autonomic neurogenic control of HR. Across the program of DI sessions, the hypotensive effect was documented, but no notable modification of the HRV-parameters was found. The absence of long-term modification of the studied parameters can be attributed both to deconditioning environmental effect of DI and limited adaptation of the organism due to neurodegeneration in PD. That should be taken into consideration when planning rehabilitation measures in subjects of older age and chronic somatic diseases with modeled microgravity.

Cardiorespiratory Coordination During Exercise in Adults With Down Syndrome

Introduction: Down syndrome (DS) is a chromosomal disorder affecting simultaneously cardiovascular and respiratory systems. There is no research studying the coupling between these systems during cardiorespiratory exercise testing in a population with DS. Cardiorespiratory coordination (CRC), evaluated through principal component analysis (PCA), measures the covariation of cardiorespiratory variables during exercise.

Objective: To investigate and compare CRC in adults with and without DS during maximal cardiorespiratory exercise testing.

Methods: Fifteen adults with DS and 15 adults without disabilities performed a maximal cardiorespiratory exercise test on a treadmill. First, the slope, and afterward the velocity was increased regularly until participants reached exhaustion. The time series of six selected cardiorespiratory variables [ventilation per minute, an expired fraction of O₂, the expired fraction of CO₂, heart rate, systolic blood pressure (SBP), and diastolic blood pressure (DBP)] were extracted for the analysis. The number of principal components (PCs), the first PC eigenvalues (PC₁), and the information entropy were computed for each group (non-DS and DS) and compared using a t-test or a Mann-Whitney U test.

Results: Two PCs in the non-DS group and three PCs in the DS group captured the variance of the studied cardiorespiratory variables. The formation of an additional PC in the DS group was the result of the shift of SBP and DBP from the PC₁ cluster of variables. Eigenvalues of PC₁ were higher in the non-DS (U = 30; p = 0.02; d = 1.47) than in the DS group, and the entropy measure was higher in the DS compared with the non-DS group (U = 37.5; p = 0.008; d = 0.70).

Conclusion: Adults with Down syndrome showed higher CRC dimensionality and a higher entropy measure than participants without disabilities. Both findings point toward a lower efficiency of the cardiorespiratory function during exercise in participants with DS. CRC appears as an alternative measure to investigate the cardiorespiratory function and its response to exercise in the DS population.

Technical aspects of cardiorespiratory estimation using subspace projections and cross entropy

BACKGROUND

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

OBJECTIVE

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

APPROACH

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

MAIN RESULTS

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

SIGNIFICANCE

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

The Value of Heart Rhythm Complexity in Identifying High-Risk Pulmonary Hypertension Patients

Pulmonary hypertension (PH) is a fatal disease—even with state-of-the-art medical treatment. Non-invasive clinical tools for risk stratification are still lacking. The aim of this study was to investigate the clinical utility of heart rhythm complexity in risk stratification for PH patients. We prospectively enrolled 54 PH patients, including 20 high-risk patients (group A; defined as WHO functional class IV or class III with severely compromised hemodynamics), and 34 low-risk patients (group B). Both linear and non-linear heart rate variability (HRV) variables, including detrended fluctuation analysis (DFA) and multiscale entropy (MSE), were analyzed. In linear and non-linear HRV analysis, low frequency and high frequency ratio, DFAα1, MSE slope 5, scale 5, and area 6–20 were significantly lower in group A. Among all HRV variables, MSE scale 5 (AUC: 0.758) had the best predictive power to discriminate the two groups. In multivariable analysis, MSE scale 5 (p = 0.010) was the only significantly predictor of severe PH in all HRV variables. In conclusion, the patients with severe PH had worse heart rhythm complexity. MSE parameters, especially scale 5, can help to identify high-risk PH patients.

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.

Two-Tiered Response of Cardiorespiratory-Cerebrovascular Network to Orthostatic Challenge

Dynamic interdependencies within and between physiological systems and subsystems are key for homeostatic mechanisms to establish an optimal state of the organism. These interactions mediate regulatory responses elicited by various perturbations, such as the high-pressure baroreflex and cerebral autoregulation, alleviating the impact of orthostatic stress on cerebral hemodynamics and oxygenation. The aim of this study was to evaluate the responsiveness of the cardiorespiratory-cerebrovascular networks by capturing linear and nonlinear interdependencies to postural changes. Ten young healthy adults participated in our study. Non-invasive measurements of arterial blood pressure (from that cardiac cycle durations were derived), breath-to-breath interval, cerebral blood flow velocity (BFV, recorded by transcranial Doppler sonography), and cerebral hemodynamics (HbT, total hemoglobin content monitored by near-infrared spectroscopy) were performed for 30-min in resting state, followed by a 1-min stand-up and a 1-min sit-down period. During preprocessing, noise was filtered and the contribution of arterial blood pressure was regressed from BFV and HbT signals. Cardiorespiratory-cerebrovascular networks were reconstructed by computing pair-wise Pearson-correlation or mutual information between the resampled signals to capture their linear and/or nonlinear interdependencies, respectively. The interdependencies between cardiac, respiratory, and cerebrovascular dynamics showed a marked weakening after standing up persisting throughout the sit-down period, which could mainly be attributed to strikingly attenuated nonlinear coupling. To summarize, we found that postural changes induced topological changes in the cardiorespiratory-cerebrovascular network. The dissolution of nonlinear networks suggests that the complexity of key homeostatic mechanisms maintaining cerebral hemodynamics and oxygenation is indeed sensitive to physiological perturbations such as orthostatic stress.

Time-evolving coupling functions for evaluating the interaction between cerebral oxyhemoglobin and arterial blood pressure with hypertension

PURPOSES

This study aimed to investigate the network coupling between arterial blood pressure (ABP) and changes in cerebral oxyhemoglobin concentration (Δ [O₂ Hb]/Δ [HHb]) oscillations based on dynamical Bayesian inference in hypertensive subjects.

METHODS

Two groups of subjects, consisting of 30 healthy (Group Control, 55.1 ± 10.6 y), and 32 hypertensive individuals (Group AH, 58.9 ± 8.7 y), participated in this study. A functional near-infrared spectroscopy system was used to measure the Δ [O₂ Hb] and Δ [HHb] signals in the bilateral prefrontal cortex (LPFC/RPFC), motor cortex (LMC/RMC), and occipital lobe (LOL/ROL) during the resting state (12 min). Based on continuous wavelet analysis and coupling functions, the directed coupling strength (CS) between ABP and cerebral hemoglobin was identified and analyzed in three frequency intervals (I: 0.6-2 Hz, II: 0.145-0.6 Hz, III: 0.01-0.08 Hz). The Pearson correlations between the CS and blood pressure parameters were calculated in the hypertension group.

RESULTS

In interval I, Group AH exhibited a significantly higher CS for the coupling from ABP to Δ [O₂ Hb] than Group Control in LMC, RMC, LOL, and ROL. In interval III, the CS from ABP to Δ [O₂ Hb] in LPFC, RPFC, LMC, RMC, LOL, and ROL was significantly higher in Group AH than in Group Control. For the patients with hypertension, diastolic blood pressure was negatively and pulse pressure was positively related to the CS from ABP to Δ [O₂ Hb] oscillations in interval III.

CONCLUSIONS

The higher CS from ABP to Δ [O₂ Hb] in interval I indicated that the components of cardiac activity in cerebral hemoglobin oscillations were more directly responsive to the changes in systematic ABP in patients with hypertension than in healthy subjects. Meanwhile, the higher CS from ABP to Δ [O₂ Hb] in interval III indicated that the cerebral hemoglobin oscillations were susceptible to changes in blood pressure in hypertensive subjects. The results may serve as evidence of impairment in cerebral autoregulation after hypertension. The Pearson correlation results showed that diastolic blood pressure and pulse pressure might be regarded as predictors of cerebral autoregulation function in patients with hypertension, and may be useful for hypertension stratification. This study provides novel insights into the interaction mechanism between ABP and cerebral hemodynamics and could help in the development of new assessment techniques for cerebral vascular disease.

Scale-Free Coupled Dynamics in Brain Networks Captured by Bivariate Focus-Based Multifractal Analysis

While most connectivity studies investigate functional connectivity (FC) in a scale-dependent manner, coupled neural processes may also exhibit broadband dynamics, manifesting as power-law scaling of their measures of interdependence. Here we introduce the bivariate focus-based multifractal (BFMF) analysis as a robust tool for capturing such scale-free relations and use resting-state electroencephalography (EEG) recordings of 12 subjects to demonstrate its performance in reconstructing physiological networks. BFMF was employed to characterize broadband FC between 62 cortical regions in a pairwise manner, with all investigated connections being tested for true bivariate multifractality. EEG channels were also grouped to represent the activity of six resting-state networks (RSNs) in the brain, thus allowing for the analysis of within- and between- RSNs connectivity, separately. Most connections featured true bivariate multifractality, which could be attributed to the genuine scale-free coupling of neural dynamics. Bivariate multifractality showed a characteristic topology over the cortex that was highly concordant among subjects. Long-term autocorrelation was higher in within-RSNs, while the degree of multifractality was generally found stronger in between-RSNs connections. These results offer statistical evidence of the bivariate multifractal nature of functional coupling in the brain and validate BFMF as a robust method to capture such scale-independent coupled dynamics.

Child-parent cardiac transference is decreased following disrupted/absent early care

Parental input shapes youth self-regulation development, and a lack of sensitive caregiving is a risk factor for youth mental health problems. Parents may shape youth regulation through their influence over biological (including neural) and behavioral development during childhood at both micro (moment-to-moment) and macro (global) levels. Prior studies have shown that micro-level parent-child interactions shape youth's biology contributing to youth mental health. However, it remains unclear whether prior disrupted/absent care affects those moment-to-moment parent-youth interactions in ways that increase risk for youth psychopathology. In the current study, we calculated transfer entropy on cardiac data from parent-youth dyads where the youth had either been exposed to disrupted care prior to adoption or not. Transfer Entropy (TE) tracks information flow between two signals, enhancing quantification of directional coupling, allowing for the examination of parent-child and child-parent influences. Using this novel approach, we identified lower cardiac information transfer from youth-to-parents in dyads where the youth had been exposed to disrupted/absent early care. Moreover we showed that the degree to which the parent's physiology changed in response to youth's physiology was negatively related to the youth's mental health, representing a potential pathway for elevated mental health risk in populations exposed to disrupted/absent early care.

Non-Contact HR Monitoring via Smartphone and Webcam During Different Respiratory Maneuvers and Body Movements

As a reliable indicator for individual's healthiness conditions, heart rate (HR) has been widely considered and used. Imaging photoplethysmography (iPPG) is recently highlighted as a promising HR measurement method, due to its non-contact characteristics, by extracting the HR from facial video recordings. In this study, we propose a camera-based HR monitoring technique that estimates HR information from iPPG signals extracted from a video sequence. Videos were recorded using a smartphone or a laptop camera. We adopted the plane-orthogonal-to-skin (POS) method to compute iPPG. The proposed method is evaluated by applying it to extract HR of 9 subjects at rest and during two motion conditions (lateral and frontal) while they were performing several respiratory maneuvers-spontaneous, metronome, and forced. Automatic face detection algorithms were implemented in the proposed method. Our experimental results show that mean values of HR have 0.56% error and 99.4% accuracy when compared to HR calculated from the gold-standard electrocardiography (ECG) reference in diverse conditions of motions and respiratory maneuvers.

Multivariate Correlation Measures Reveal Structure and Strength of Brain-Body Physiological Networks at Rest and During Mental Stress

In this work, we extend to the multivariate case the classical correlation analysis used in the field of network physiology to probe dynamic interactions between organ systems in the human body. To this end, we define different correlation-based measures of the multivariate interaction (MI) within and between the brain and body subnetworks of the human physiological network, represented, respectively, by the time series of δ, θ, α, and β electroencephalographic (EEG) wave amplitudes, and of heart rate, respiration amplitude, and pulse arrival time (PAT) variability (η, ρ, π). MI is computed: (i) considering all variables in the two subnetworks to evaluate overall brain-body interactions; (ii) focusing on a single target variable and dissecting its global interaction with all other variables into contributions arising from the same subnetwork and from the other subnetwork; and (iii) considering two variables conditioned to all the others to infer the network topology. The framework is applied to the time series measured from the EEG, electrocardiographic (ECG), respiration, and blood volume pulse (BVP) signals recorded synchronously via wearable sensors in a group of healthy subjects monitored at rest and during mental arithmetic and sustained attention tasks. We find that the human physiological network is highly connected, with predominance of the links internal of each subnetwork (mainly η-ρ and δ-θ, θ-α, α-β), but also statistically significant interactions between the two subnetworks (mainly η-β and η-δ). MI values are often spatially heterogeneous across the scalp and are modulated by the physiological state, as indicated by the decrease of cardiorespiratory interactions during sustained attention and by the increase of brain-heart interactions and of brain-brain interactions at the frontal scalp regions during mental arithmetic. These findings illustrate the complex and multi-faceted structure of interactions manifested within and between different physiological systems and subsystems across different levels of mental stress.

Psychophysiological responses to various slow, deep breathing techniques

Deep breathing exercises are commonly used for several health conditions including pain and hypertension. Various techniques are available to practice deep breathing, whereas possible differential psychophysiological effects have not been investigated. We compared four deep breathing techniques and examined outcomes in blood pressure variability, respiratory sinus arrhythmia, baroreflex function, and emotional state. Healthy adult volunteers performed pursed-lips breathing, left and right unilateral nostril breathing, and deep breathing with an inspiratory threshold load (loaded breathing), all at a frequency of 0.1 Hz (i.e., controlled breathing) and for three minutes each. Results showed that blood pressure variability was higher during loaded breathing versus other conditions and higher during pursed-lips breathing versus left and right unilateral nostril breathing. Respiratory sinus arrhythmia was higher during loaded breathing versus other conditions and higher during pursed-lips breathing versus left unilateral nostril breathing. The effect of breathing condition on respiratory sinus arrhythmia was mediated by alterations in blood pressure variability. There was no difference between the breathing conditions in baroreflex sensitivity or effectiveness. Participants rated pursed-lips breathing as more calming and pleasant and with more sense of control (vs. other conditions). Overall, among the four tested deep breathing techniques, loaded breathing was associated with enhanced cardiovascular effects and pursed-lips breathing with better emotional responses, while also enhancing cardiovascular effects (albeit less than loaded breathing). These findings can be informative in applying deep breathing techniques as self-management interventions for health conditions, in which baroreceptors stimulation and autonomic and emotional modulations can be beneficial, such as pain and hypertension.

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.

Mutual Information Analysis of Brain-Body Interactions during different Levels of Mental stress

In this work, we analyze brain-heart interactions during different mental states computing mutual information (MI) between the dynamic activity of different physiological systems. In 18 healthy subjects monitored in a relaxed resting state and during a mental arithmetic and a serious game task, multichannel EEG, one lead ECG, respiration and blood volume pulse were collected via wireless non-invasive biosensors. From these signals, synchronous 300-second time series were extracted measuring brain activity via the δ, θ, α, and β EEG power, and activity of the body district via the ECG R-R interval η, the respiratory amplitude ρ and the pulse arrival time π. MI was computed using a linear estimator: (i) between {η,ρ,π} and {δ,θ,α,β}, to measure overall brain-body interactions; (ii) between each time series and the others of the same district, to measure information shared within a district; and (iii) between each time series of a district and all series of the other district, to evaluate individual contributions to the information shared between brain and body. Results document the existence of statistically significant brain-body interactions, with high MI values involving mainly the η body dynamics and the δ and β brain dynamics. State-dependent variations were mostly relevant to the MI of the brain system involving δ, θ, α during mental arithmetic, and α and β during serious game. Thus, MI can be useful to detect correlated activity within and between brain and body systems monitored simultaneously during different mental states.

The Cardiorespiratory Network in Healthy First-Degree Relatives of Schizophrenic Patients

Impaired heart rate- and respiratory regulatory processes as a sign of an autonomic dysfunction seems to be obviously present in patients suffering from schizophrenia. Since the linear and non-linear couplings within the cardiorespiratory system with respiration as an important homeostatic control mechanism are only partially investigated so far for those subjects, we aimed to characterize instantaneous cardiorespiratory couplings by quantifying the casual interaction between heart rate (HR) and respiration (RESP). Therefore, we investigated causal linear and non-linear cardiorespiratory couplings of 23 patients suffering from schizophrenia (SZO), 20 healthy first-degree relatives (REL) and 23 healthy subjects, who were age-gender matched (CON). From all participants' heart rate (HR) and respirations (respiratory frequency, RESP) were investigated for 30 min under resting conditions. The results revealed highly significant increased HR, reduced HR variability, increased respiration rates and impaired cardiorespiratory couplings in SZO in comparison to CON. SZO were revealed bidirectional couplings, with respiration as the driver (RESP → HR), and with weaker linear and non-linear coupling strengths when RESP influencing HR (RESP → HR) and with stronger linear and non-linear coupling strengths when HR influencing RESP (HR → RESP). For REL we found only significant increased HR and only slightly reduced cardiorespiratory couplings compared to CON. These findings clearly pointing to an underlying disease-inherent genetic component of the cardiac system for SZO and REL, and those respiratory alterations are only clearly present in SZO seem to be connected to their mental emotional states.

Application of the nonlinear methods in pneumocardiogram signals

In this work, the pneumocardiogram signals of nine rats were analysed by scale index, Boltzmann Gibbs entropy and maximum Lyapunov exponents. The scale index method, based on wavelet transform, was proposed for determining the degree of aperiodicity and chaos. It means that the scale index parameter is close to zero when the signal is periodic and has a value between zero and one when the signal is aperiodic. A new entropy calculation method by normalized inner scalogram was suggested very recently. In this work, we also used this method for the first time in an empirical data. We compared the both methods with maximum Lyapunov exponents and observed that using together the scale index and the entropy calculation method by normalized inner scalogram increases the reliability of the pneumocardiogram signal analysis. Thus, the analysis of the pneumocardiogram signals by those methods enables to compare periodical and/or nonlinear aspects for further understanding of dynamics of cardiorespiratory system.

Time Window Determination for Inference of Time-Varying Dynamics: Application to Cardiorespiratory Interaction

Interacting dynamical systems abound in nature, with examples ranging from biology and population dynamics, through physics and chemistry, to communications and climate. Often their states, parameters and functions are time-varying, because such systems interact with other systems and the environment, exchanging information and matter. A common problem when analysing time-series data from dynamical systems is how to determine the length of the time window for the analysis. When one needs to follow the time-variability of the dynamics, or the dynamical parameters and functions, the time window needs to be resolved first. We tackled this problem by introducing a method for adaptive determination of the time window for interacting oscillators, as modeled and scaled for the cardiorespiratory interaction. By investigating a system of coupled phase oscillators and utilizing the Dynamical Bayesian Inference method, we propose a procedure to determine the time window and the propagation parameter of the covariance matrix. The optimal values are determined so that the inferred parameters follow the dynamics of the actual ones and at the same time the error of the inference represented by the covariance matrix is minimal. The effectiveness of the methodology is presented on a system of coupled limit-cycle oscillators and on the cardiorespiratory interaction. Three cases of cardiorespiratory interaction were considered-measurement with spontaneous free breathing, one with periodic sine breathing and one with a-periodic time-varying breathing. The results showed that the cardiorespiratory coupling strength and similarity of form of coupling functions have greater values for slower breathing, and this variability follows continuously the change of the breathing frequency. The method can be applied effectively to other time-varying oscillatory interactions and carries important implications for analysis of general dynamical systems.

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.

Slow 0.1 Hz Breathing and Body Posture Induced Perturbations of RRI and Respiratory Signal Complexity and Cardiorespiratory Coupling

Objective: We explored the physiological background of the non-linear operating mode of cardiorespiratory oscillators as the fundamental question of cardiorespiratory homeodynamics and as a prerequisite for the understanding of neurocardiovascular diseases. We investigated 20 healthy human subjects for changes using electrocardiac RR interval (RRI) and respiratory signal (Resp) Detrended Fluctuation Analysis (DFA, α1RRI, α2RRI, α1Resp, α2Resp), Multiple Scaling Entropy (MSERRI1-4, MSERRI5-10, MSEResp1-4, MSEResp5-10), spectral coherence (CohRRI-Resp), cross DFA (ρ1 and ρ2) and cross MSE (XMSE1-4 and XMSE5-10) indices in four physiological conditions: supine with spontaneous breathing, standing with spontaneous breathing, supine with 0.1 Hz breathing and standing with 0.1 Hz breathing. Main results: Standing is primarily characterized by the change of RRI parameters, insensitivity to change with respiratory parameters, decrease of CohRRI-Resp and insensitivity to change of in ρ1, ρ2, XMSE1-4, and XMSE5-10. Slow breathing in supine position was characterized by the change of the linear and non-linear parameters of both signals, reflecting the dominant vagal RRI modulation and the impact of slow 0.1 Hz breathing on Resp parameters. CohRRI-Resp did not change with respect to supine position, while ρ1 increased. Slow breathing in standing reflected the qualitatively specific state of autonomic regulation with striking impact on both cardiac and respiratory parameters, with specific patterns of cardiorespiratory coupling. Significance: Our results show that cardiac and respiratory short term and long term complexity parameters have different, state dependent patterns. Sympathovagal non-linear interactions are dependent on the pattern of their activation, having different scaling properties when individually activated with respect to the state of their joint activation. All investigated states induced a change of α1 vs. α2 relationship, which can be accurately expressed by the proposed measure-inter-fractal angle θ. Short scale (α1 vs. MSE1-4) and long scale (α2 vs. MSE5-10) complexity measures had reciprocal interrelation in standing with 0.1 Hz breathing, with specific cardiorespiratory coupling pattern (ρ1 vs. XMSE1-4). These results support the hypothesis of hierarchical organization of cardiorespiratory complexity mechanisms and their recruitment in ascendant manner with respect to the increase of behavioral challenge complexity. Specific and comprehensive cardiorespiratory regulation in standing with 0.1 Hz breathing suggests this state as the potentially most beneficial maneuver for cardiorespiratory conditioning.

The Resting-State Pulse-Respiration Quotient of Humans: Lognormally Distributed and Centered Around a Value of Four

The pulse-respiration quotient (heart rate divided by the respiration rate, PRQ = HR/RR) is a parameter capturing the complex state of cardiorespiratory interactions. We analysed 482 single PRQ values obtained from measurement on 134 healthy adult subjects (49 men, 85 women, age: 24.7 ± 3.4, range: 20-46 years) during rest. We found that the distribution of PRQ values (i) has a global maximum at around a value of 4 (median: 4.19) and (ii) follows a lognormal distribution function. A multimodality of the distribution, associated with several PRQ attractor states was not detected by our group-level based analysis. In summary, our analysis shows that in healthy humans the resting-state PRQ is around 4 and lognormally distributed. This finding supports claims about the special role of the 4 to 1 cardiorespiratory coupling in particular and the PRQ in general for physiological and medical views and applications. To the best of our knowledge, our study is the largest conducted so far in healthy adult humans about reference values of the PRQ during a resting-state at day.

An Overview of the Cardiorespiratory Hypothesis and Its Potential Contribution to the Care of Neurodegenerative Disease in Africa

One hypothesis that could explain the beneficial effects of physical exercise on cognitive function is the cardiorespiratory hypothesis. This hypothesis proposes that improved cognitive functioning may be in part a result of the physiological processes that occur after physical exercise such as: Increased cerebral perfusion and regional cerebral blood flow. These processes ensure increased oxygenation and glucose transportation to the brain, which together can improve cognitive function. The objective of this narrative review is to examine the contribution of this hypothesis in the care of African older adults with neurodegenerative conditions (i.e., dementia (Alzheimer’s disease)) or with mild cognitive impairments. Although studies in developed countries have examined people of African descent (i.e., with African Americans), only the limited findings presented in this review reflect how these conditions are also important for the African continent. This review revealed that no studies have examined the effects of cardiorespiratory fitness on neurodegenerative disease in Africa. African nations, like many other developing countries, have an aging population that is growing and will face an increased risk of neurodegenerative declines. It is therefore imperative that new research projects be developed to explore the role of the cardiorespiratory fitness in neurodegenerative disease prevention in African nations.