Model for cardiorespiratory synchronization in humans

Sleep-stage dependence and co-existence of cardio-respiratory coordination and phase synchronization

Interactions between the cardiac and respiratory systems play a pivotal role in physiological functioning. Nonetheless, the intricacies of cardio-respiratory couplings, such as cardio-respiratory phase synchronization (CRPS) and cardio-respiratory coordination (CRC), remain elusive, and an automated algorithm for CRC detection is lacking. This paper introduces an automated CRC detection algorithm, which allowed us to conduct a comprehensive comparison of CRPS and CRC during sleep for the first time using an extensive database. We found that CRPS is more sensitive to sleep-stage transitions, and intriguingly, there is a negative correlation between the degree of CRPS and CRC when fluctuations in breathing frequency are high. This comparative analysis holds promise in assisting researchers in gaining deeper insights into the mechanics of and distinctions between these two physiological phenomena. Additionally, the automated algorithms we devised have the potential to offer valuable insights into the clinical applications of CRC and CRPS.

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.

Age-dependent cardiorespiratory directional coupling in wake-resting state

Objective.Cooperation in the cardiorespiratory system helps maintain internal stability. Various types of system interactions have been investigated; however, the characteristics of the interactions have mostly been studied using data collected in well-defined physiological states, such as sleep. Furthermore, most analyses provided general information about the interaction, making it difficult to quantify how the systems influenced one another.Approach.Cardiorespiratory directional coupling was investigated in different age groups (20 young and 19 elderly subjects) in a wake-resting state. The directionality index (DI) was calculated using instantaneous phases from the heartbeat interval and respiratory signal to provide information about the strength and direction of interaction between the systems. Statistical analysis was performed between the groups on the DI and independent measures of directionality (nc_r: influence from cardiac system to respiratory system, and nc_c: influence from the respiratory system to the cardiac system).Main results.The values of DI were -0.52 and -0.17 in the young and elderly groups, respectively (p< 0.001). Furthermore, the values of nc_rand nc_cwere found to be significantly different between the groups (p< 0.001), respectively.Significance.Changes in both directions between the systems influence different aspects of cardiorespiratory coupling between the groups. This observation could be linked to different levels of autonomic modulation associated with ageing. Our approach could aid in quantitatively tracking and comprehending how systems interact in response to physiological and environmental changes. It could also be used to understand how abnormal interaction characteristics influence physiological system dysfunctions and disorders.

Pulse respiration quotient as a measure sensitive to changes in dynamic behavior of cardiorespiratory coupling such as body posture and breathing regime

Objective: In this research we explored the (homeo)dynamic character of cardiorespiratory coupling (CRC) under the influence of different body posture and breathing regimes. Our tool for it was the pulse respiration quotient (PRQ), representing the number of heartbeat intervals per breathing cycle. We obtained non-integer PRQ values using our advanced Matlab^® algorithm and applied it on the signals of 20 healthy subjects in four conditions: supine position with spontaneous breathing (Supin), standing with spontaneous breathing (Stand), supine position with slow (0.1 Hz) breathing (Supin01) and standing with slow (0.1 Hz) breathing (Stand01). Main results: Linear features of CRC (in PRQ signals) were dynamically very sensitive to posture and breathing rhythm perturbations. There are obvious increases in PRQ mean level and variability under the separated and joined influence of orthostasis and slow (0.1 Hz) breathing. This increase was most pronounced in Stand01 as the state of joint influences. Importantly, PRQ dynamic modification showed greater sensitivity to body posture and breathing regime changes than mean value and standard deviation of heart rhythm and breathing rhythm. In addition, as a consequence of prolonged supine position, we noticed the tendency to integer quantization of PRQ (especially after 14 min), in which the most common quantization number was 4:1 (demonstrated in other research reports as well). In orthostasis and slow breathing, quantization can also be observed, but shifted to other values. We postulate that these results manifest resonance effects induced by coupling patterns from sympathetic and parasympathetic adjustments (with the second as dominant factor). Significance: Our research confirms that cardiorespiratory coupling adaptability could be profoundly explored by precisely calculated PRQ parameter since cardiorespiratory regulation in healthy subjects is characterized by a high level of autonomic adaptability (responsiveness) to posture and breathing regime, although comparisons with pathological states has yet to be performed. We found Stand01 to be the most provoking state for the dynamic modification of PRQ (cardiorespiratory inducement). As such, Stand01 has the potential of using for PRQ tuning by conditioning the cardiorespiratory autonomic neural networks, e.g., in the cases where PRQ is disturbed by environmental (i.e., microgravity) or pathologic conditions.

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.

Correlation Analysis of Synchronization Type and Degree in Respiratory Neural Network

Pre-Bötzinger complex (PBC) is a necessary condition for the generation of respiratory rhythm. Due to the existence of synaptic gaps, delay plays a key role in the synchronous operation of coupled neurons. In this study, the relationship between synchronization and correlation degree is established for the first time by using ISI bifurcation and correlation coefficient, and the relationship between synchronization and correlation degree is discussed under the conditions of no delay, symmetric delay, and asymmetric delay. The results show that the phase synchronization of two coupling PBCs is closely related to the weak correlation, that is, the weak phase synchronization may occur under the condition of incomplete synchronization. Moreover, the time delay and coupling strength are controlled in the modified PBC network model, which not only reveals the law of PBC firing transition but also reveals the complex synchronization behavior in the coupled chaotic neurons. Especially, when the two coupled neurons are nonidentical, the complete synchronization will disappear. These results fully reveal the dynamic behavior of the PBC neural system, which is helpful to explore the signal transmission and coding of PBC neurons and provide theoretical value for further understanding respiratory rhythm.

Characteristics of mathematical modeling languages that facilitate model reuse in systems biology: a software engineering perspective

Reuse of mathematical models becomes increasingly important in systems biology as research moves toward large, multi-scale models composed of heterogeneous subcomponents. Currently, many models are not easily reusable due to inflexible or confusing code, inappropriate languages, or insufficient documentation. Best practice suggestions rarely cover such low-level design aspects. This gap could be filled by software engineering, which addresses those same issues for software reuse. We show that languages can facilitate reusability by being modular, human-readable, hybrid (i.e., supporting multiple formalisms), open, declarative, and by supporting the graphical representation of models. Modelers should not only use such a language, but be aware of the features that make it desirable and know how to apply them effectively. For this reason, we compare existing suitable languages in detail and demonstrate their benefits for a modular model of the human cardiac conduction system written in Modelica.

Effect of Stress on Cardiorespiratory Synchronization of Ironman Athletes

The aim of this paper is to investigate the cardiorespiratory synchronization in athletes subjected to extreme physical stress combined with a cognitive stress tasks. ECG and respiration were measured in 14 athletes before and after the Ironman competition. Stroop test was applied between the measurements before and after the Ironman competition to induce cognitive stress. Synchrogram and empirical mode decomposition analysis were used for the first time to investigate the effects of physical stress, induced by the Ironman competition, on the phase synchronization of the cardiac and respiratory systems of Ironman athletes before and after the competition. A cognitive stress task (Stroop test) was performed both pre- and post-Ironman event in order to prevent the athletes from cognitively controlling their breathing rates. Our analysis showed that cardiorespiratory synchronization increased post-Ironman race compared to pre-Ironman. The results suggest that the amount of stress the athletes are recovering from post-competition is greater than the effects of the Stroop test. This indicates that the recovery phase after the competition is more important for restoring and maintaining homeostasis, which could be another reason for stronger synchronization.

The Effect of a Non-Local Fractional Operator in an Asymmetrical Glucose-Insulin Regulatory System: Analysis, Synchronization and Electronic Implementation

For studying biological conditions with higher precision, the memory characteristics defined by the fractional-order versions of living dynamical systems have been pointed out as a meaningful approach. Therefore, we analyze the dynamics of a glucose-insulin regulatory system by applying a non-local fractional operator in order to represent the memory of the underlying system, and whose state-variables define the population densities of insulin, glucose, and β-cells, respectively. We focus mainly on four parameters that are associated with different disorders (type 1 and type 2 diabetes mellitus, hypoglycemia, and hyperinsulinemia) to determine their observation ranges as a relation to the fractional-order. Like many preceding works in biosystems, the resulting analysis showed chaotic behaviors related to the fractional-order and system parameters. Subsequently, we propose an active control scheme for forcing the chaotic regime (an illness) to follow a periodic oscillatory state, i.e., a disorder-free equilibrium. Finally, we also present the electronic realization of the fractional glucose-insulin regulatory model to prove the conceptual findings.

Physical fitness contributes to cardio-respiratory synchronization

Cardio-respiratory synchronization is a phenomenon of particular interest- especially at a 1:1 ratio- and may give greater insight into the underlying mechanisms of cardio-respiratory communication. Synchronization of this ratio is hypothesised to occur when breathing rate exceeds heart rate, which is the premise of this research. A novel experimental design focused on guiding elevated respiration to induce the entrainment of heart rate, and produce an equivalent rise in value. Application of instantaneous phase for identification and analysis of synchronization allowed for a reliable method of measuring the interaction between these stochastic processes. We have identified 1:1 phase synchronization in all volunteers measured. Longer synchronization episodes were observed reliably in athletic individuals, corroborating previous research for spontaneous breathing. This observation suggests that cardio-respiratory synchronization at all respiration rates is associated with a common underlying communication mechanism. Furthermore, it presents cardio-respiratory synchronization as a potential future measurement of fitness and autonomic health.

Control of heart rate through guided high-rate breathing

Understanding the complex dynamics of cardio-respiratory coupling sheds light on the underlying mechanisms governing the communication between these two physiological systems. Previous research has predominantly considered the coupling at respiratory rates slower than the heart rate and shown that respiratory oscillations lead to modulation and/or synchronization of the heart rate. Whereas the mechanisms of cardio-respiratory communication are still under discussion, peripheral nervous regulation is considered to be the predominant factor. This work offers a novel experimental design and applies the concept of instantaneous phase to detect cardio-respiratory entrainment at elevated respiration rates, close to the resting heart rate. If such 1:1 entrainment exists, it would suggest direct neuronal communication between the respiration and heart centres in the brain. We have observed 1:1 entrainment in all volunteers, with consistently longer synchronization episodes seen in physically fitter people, and demonstrated that cardio-respiratory synchronization at both low and high respiration rates is associated with a common underlying communication mechanism.

A quantitative model of relation between respiratory-related blood pressure fluctuations and the respiratory sinus arrhythmia

In order to propose an interpretation of recent experimental findings concerning short-term variability of arterial blood pressure (ABP), heart rate variability (HRV), and their dependence on body posture, we develop a qualitative dynamical model of the short-term cardiovascular variability at respiratory frequency (HF). It shows the respiratory-related blood pressure fluctuations in relation to the respiratory sinus arrhythmia (RSA). Results of the model-based analysis show that the observed phenomena may be interpreted as buffering of the respiratory-related ABP fluctuations by heart rate (HR) fluctuations, i.e., the respiratory sinus arrhythmia. A paradoxical enhancement (PE) of the fluctuations of the ABP in supine position, that was found in experiment, is explained on the ground of the model, as an ineffectiveness of control caused by the prolonged phase shift between the the peak of modulation of the pulmonary flow and the onset of stimulation of the heart. Such phasic changes were indeed observed in some other experimental conditions. Up to now, no other theoretical or physiological explanation of the PE effect exists, whereas further experiments were not performed due to technical problems. Better understanding of the short-term dynamics of blood pressure may improve medical diagnosis in cardiology and diseases which alter the functional state of the autonomous nervous system.

Graphical Abstract

The Baroreflex Mechanism Revisited

We state that the autonomic part of the brain controls the blood pressure (BP) and the heart rate (HR) via the baroreflex mechanism in all situations of human activity (at sleep, at rest, during exercise, fright etc.), in a way which is not, as was hitherto assumed, a mere homeostatic tool or even a resetting device, designed to bring these variables on the road to preset values. The baroreflex is rather a continuous feedback mechanism commanded by the autonomic part of the brain, leading to values appropriate to the situation at hand. Feasibility of this assertion is demonstrated here by using the Seidel-Herzel feedback system outside of its regular practice. Results show indeed that the brain can, and we claim that it does, control the HR and BP throughout life. New responses are demonstrated, e.g., to a sudden fear or apnea. In this event, large BP and HR overshoots are expected before the variables can relax to a new level. Response to abrupt downward change in the controlling parameter shows an undershoot in HR and just a gradual resetting in the BP. The relaxation from sudden external changes to various expected states are calculated and discussed and properties of the Rheos test are explained. Experimental findings for orthostatic tests and for babies under translations and rotations reveal complete qualitative agreement with our model and show no need to invoke the operation of additional body systems. Our method should be the preferred one by the Occam Razor approach. The outcomes may lead to beneficial clinical implication.

Extraction of response waveforms of heartbeat and blood pressure to swallowing. Using mixed signal processing of time domain and respiratory phase domain

BACKGROUND

Evaluating the accurate responses of the cardiovascular system to external stimuli is important for a deeper understanding of cardiovascular homeostasis. However, the responses should be distorted by the conventional time domain analysis when a frequency of the effect of external stimuli matches that of intrinsic fluctuations.

OBJECTIVES

The purpose of this study is to propose a mixed signal processing of time domain and respiratory phase domain to extract the response waveforms of heartbeat and blood pressure (BP) to external stimuli and to clarify the physiological mechanisms of swallowing effects on the cardiovascular system.

METHODS

Measurements were conducted on 12 healthy humans in the sitting and standing positions, with each subject requested to swallow every 30 s between expiration and inspiration. Waveforms of respiratory sinus arrhythmia (RSA) and respiratory-related BP variations were extracted as functions of the respiratory phase. Then, respiratory effects were subtracted from response waveforms with reference to the respiratory phase in the time domain.

RESULTS

As a result, swallowing induced tachycardia, which peaked within 3 s and recovered within 8 s. Tachycardia was greater in the sitting position than during standing. Furthermore, systolic BP and pulse pressure immediately decreased and diastolic BP increased coincident with the occurrence of tachycardia. Subsequently, systolic BP and pulse pressure recovered faster than the R-R interval.

CONCLUSIONS

We conclude that swallowing-induced tachycardia arises largely from the decrease of vagal activity and the baroreflex would yield fast oscillatory responses in recovery.

Central regulation of heart rate and the appearance of respiratory sinus arrhythmia: new insights from mathematical modeling

A minimal model for the neural control of heart rate (HR) has been developed with the aim of better understanding respiratory sinus arrhythmia (RSA)--a modulation of HR at the frequency of breathing. This model consists of two differential equations and is integrated into a previously-published model of gas exchange. The heart period is assumed to be affected primarily by the parasympathetic signal, with the sympathetic signal taken as a parameter in the model. We include the baroreflex, mechanical stretch-receptor feedback from the lungs, and central modulation of the cardiac vagal tone by the respiratory drive. Our model mimics a range of experimental observations and provides several new insights. Most notably, the model mimics the growth in the amplitude of RSA with decreasing respiratory frequency up to 7 breaths per minute (for humans). Our model then mimics the decrease in the amplitude of RSA at frequencies below 7 breaths per minute and predicts that this decrease is due to the baroreflex (we show this both numerically and analytically with a linear baroreflex). Another new prediction of the model is that the gating of the baroreflex leads to the dependency of RSA on mean vagal tone. The new model was also used to test two previously-suggested hypotheses regarding the physiological function of RSA and supports the hypothesis that RSA minimizes the work done by the heart while maintaining physiological levels of arterial CO2. These and other new insights the model provides extend our understanding of the integrative nature of vagal control of the heart.

Cardio-respiratory coordination increases during sleep apnea

Cardiovascular diseases are the main source of morbidity and mortality in the United States with costs of more than $170 billion. Repetitive respiratory disorders during sleep are assumed to be a major cause of these diseases. Therefore, the understanding of the cardio-respiratory regulation during these events is of high public interest. One of the governing mechanisms is the mutual influence of the cardiac and respiratory oscillations on their respective onsets, the cardio-respiratory coordination (CRC). We analyze this mechanism based on nocturnal measurements of 27 males suffering from obstructive sleep apnea syndrome. Here we find, by using an advanced analysis technique, the coordigram, not only that the occurrence of CRC is significantly more frequent during respiratory sleep disturbances than in normal respiration (p-value<10⁻⁵¹) but also more frequent after these events (p-value<10⁻¹⁵). Especially, the latter finding contradicts the common assumption that spontaneous CRC can only be observed in epochs of relaxed conditions, while our newly discovered epochs of CRC after disturbances are characterized by high autonomic stress. Our findings on the connection between CRC and the appearance of sleep-disordered events require a substantial extension of the current understanding of obstructive sleep apneas and hypopneas.

Fine temporal structure of cardiorespiratory synchronization

Cardiac and respiratory rhythms are known to exhibit a modest degree of phase synchronization, which is affected by age, diseases, and other factors. We study the fine temporal structure of this synchrony in healthy young, healthy elderly, and elderly subjects with coronary artery disease. We employ novel time-series analysis to explore how phases of oscillations go in and out of the phase-locked state at each cycle of oscillations. For the first time we show that cardiorespiratory system is engaged in weakly synchronized dynamics with a very specific temporal pattern of synchrony: the oscillations go out of synchrony frequently, but return to the synchronous state very quickly (usually within just 1 cycle of oscillations). Properties of synchrony depended on the age and disease status. Healthy subjects exhibited more synchrony at the higher (1:4) frequency-locking ratio between respiratory and cardiac rhythms, whereas subjects with coronary artery disease exhibited relatively more 1:2 synchrony. However, multiple short desynchronization episodes prevailed regardless of the age and disease status. The same average synchrony level could be alternatively achieved with few long desynchronizations, but this was not observed in the data. This implies functional importance of short desynchronization dynamics. These dynamics suggest that a synchronous state is easy to create if needed but is also easy to break. Short desynchronization dynamics may facilitate the mutual coordination of cardiac and respiratory rhythms by creating intermittent synchronous episodes. It may be an efficient background dynamics to promote adaptation of cardiorespiratory coordination to various external and internal factors.

Evaluating the physiological significance of respiratory sinus arrhythmia: looking beyond ventilation-perfusion efficiency

We conducted a theoretical study of the physiological significance of respiratory sinus arrhythmia (RSA), a phenomenon used as an index of cardiac vagal tone and wellbeing, whereby the heart rate (HR) increases during inspiration and decreases during expiration. We first tested the hypothesis that RSA improves gas exchange efficiency but found that although gas exchange efficiency improved with slow and deep breathing and with increased mean heart rate, this was unrelated to RSA. We then formulated and tested a new hypothesis: that RSA minimizes the work done by the heart while maintaining physiological levels of arterial carbon dioxide. We tested the new hypothesis using two methods. First, the HR for which the work is minimized was calculated using techniques from optimal control theory. This calculation was done on simplified models that we derived from a previously published model of gas exchange in mammals. We found that the calculated HR was remarkably similar to RSA and that this became more profound under slow and deep breathing. Second, the HR was prescribed and the work done by the heart was calculated by conducting a series of numerical experiments on the previously published gas exchange model. We found that cardiac work was minimized for RSA-like HR functions, most profoundly under slow and deep breathing. These findings provide novel insights into potential reasons for and benefits of RSA under physiological conditions.

Automated quantification of the synchrogram by recurrence plot analysis

Recently, the concept of phase synchronization of two weakly coupled oscillators has raised a great research interest and has been applied to characterize synchronization phenomenon in physiological data. Phase synchronization of cardiorespiratory coupling is often studied by a synchrogram analysis, a graphical tool investigating the relationship between instantaneous phases of two signals. Although several techniques have been proposed to automatically quantify the synchrogram, most of them require a preselection of a phase-locking ratio by trial and error. One technique does not require this information; however, it is based on the power spectrum of phase's distribution in the synchrogram, which is vulnerable to noise. This study aims to introduce a new technique to automatically quantify the synchrogram by studying its dynamic structure. Our technique exploits recurrence plot analysis, which is a well-established tool for characterizing recurring patterns and nonstationarities in experiments. We applied our technique to detect synchronization in simulated and measured infants' cardiorespiratory data. Our results suggest that the proposed technique is able to systematically detect synchronization in noisy and chaotic data without preselecting the phase-locking ratio. By embedding phase information of the synchrogram into phase space, the phase-locking ratio is automatically unveiled as the number of attractors.

Computational models for the study of heart-lung interactions in mammals

The operation and regulation of the lungs and the heart are closely related. This is evident when examining the anatomy within the thorax cavity, in the brainstem and in the aortic and carotid arteries where chemoreceptors and baroreceptors, which provide feedback affecting the regulation of both organs, are concentrated. This is also evident in phenomena such as respiratory sinus arrhythmia where the heart rate increases during inspiration and decreases during expiration, in other types of synchronization between the heart and the lungs known as cardioventilatory coupling and in the association between heart failure and sleep apnea where breathing is interrupted periodically by periods of no-breathing. The full implication and physiological significance of the cardiorespiratory coupling under normal, pathological, or extreme physiological conditions are still unknown and are subject to ongoing investigation both experimentally and theoretically using mathematical models. This article reviews mathematical models that take heart-lung interactions into account. The main ideas behind low dimensional, phenomenological models for the study of the heart-lung synchronization and sleep apnea are described first. Higher dimensions, physiology-based models are described next. These models can vary widely in detail and scope and are characterized by the way the heart-lung interaction is taken into account: via gas exchange, via the central nervous system, via the mechanical interactions, and via time delays. The article emphasizes the need for the integration of the different sources of heart-lung coupling as well as the different mathematical approaches.

Time delay correction of the synchrogram for optimized detection of cardiorespiratory coordination

The cardiorespiratory synchrogram, a graphical tool based on the stroboscopic technique, is an established method for evaluating phase-locking between cardiac and respiratory oscillators. In the original method, the phase of the respiratory oscillator is observed at the instants of time when the phase of the cardiac oscillator attains a certain value. In this article, we introduced an additional adaptive delay in the cardiac oscillator based on the maximisation of the cross-correlation or symbolic coupling traces between the phases of respiration and the delayed R-R intervals. We then investigated phase coordination in thirteen normal subjects (five males, eight females; age: 19-24 years) for different body postures. Cardiorespiratory coordination was observed to be significantly reduced in the upright position (supine vs. upright: 11.9 ± 5.1 vs. 6.9 ± 3.6, P < 0.05). Compared to the original algorithm we observed an increase in the detection of average cardiorespiratory coordination (supine original vs. delay: 11.9 vs. 18.9%), together with a decrease in standard deviation of the percentage of coordination in all the subjects, after introducing the heart rate delay (supine original vs. delay: 5.1 vs. 4.4%). In conclusion, the performance of the synchrogram technique was improved by including an adaptive delay in the cardiac oscillator.

HRV strongly depends on breathing. Are we questioning the right suspect?

The fact that the heart rate variability (HRV) depends on breathing is well known. The HRV is an important phenomenon which reflects the functional state of the autonomous nervous system (ANS), although there are some doubts concerning the actual interpretation of spectral components of HRV and their postulated balance. The assessment of the functional state of the ANS is the task of paramount importance in risk stratification of cardiological patients. HRV is considered to depend mainly on the properties of the sinus node (SN), which achieves neurohumoral input from the ANS. Interestingly, there is growing evidence that the relation between the heart rate (HR) and breathing rate (BR) is really strong. The variety of breathing-related effects that are present in HRV is very rich, including respiratory sinus arrhythmia (RSA), cardiorespiratory synchronization and vivid heart rate response to breathing disorders. If the mean frequency of any of rhythms is changed, the other rhythm adjusts itself. This provokes the question on the actual source of the dynamics observed in the HRV. Is it possible that we observe mainly the dynamics of the respiratory rhythm which is just transduced by the heart effector? What might be the role of the intrinsic dynamics of this effector? Is the RSA a product of neural regulation or rather a by-product: what is its teleological role? In consequence: if we concentrate on the sinus node and its properties in order to understand the nature of the HRV - are we questioning the right suspect? The reasoning is supplied by suitable choice of literature and by the analysis of the computational model. Various consequences are discussed.

Development of Deskwork Support System Using Pointing Gesture Interface

The authors have proposed a deskwork support system called an “Attentive Workbench” (AWB), which uses a camera, a projector, and automatically moving trays to support the user both physically and informationally. This research intends to build an interface between the user and the AWB system based on pointing gestures. Considering the specific purposes of the AWB, a simple, reliable, and highly responsive interface is implemented. We demonstrate a deskwork support system with actual automatically moving trays that uses the proposed interface.

Cardiorespiratory phase-coupling is reduced in patients with obstructive sleep apnea

Cardiac and respiratory rhythms reveal transient phases of phase-locking which were proposed to be an important aspect of cardiorespiratory interaction. The aim of this study was to quantify cardio-respiratory phase-locking in obstructive sleep apnea (OSA). We investigated overnight polysomnography data of 248 subjects with suspected OSA. Cardiorespiratory phase-coupling was computed from the R-R intervals of body surface ECG and respiratory rate, calculated from abdominal and thoracic sensors, using Hilbert transform. A significant reduction in phase-coupling was observed in patients with severe OSA compared to patients with no or mild OSA. Cardiorespiratory phase-coupling was also associated with sleep stages and was significantly reduced during rapid-eye-movement (REM) sleep compared to slow-wave (SW) sleep. There was, however, no effect of age and BMI on phase coupling. Our study suggests that the assessment of cardiorespiratory phase coupling may be used as an ECG based screening tool for determining the severity of OSA.

Periodic heart rate decelerations in premature infants

The pacemaking system of the heart is complex; a healthy heart constantly integrates and responds to extracardiac signals, resulting in highly complex heart rate patterns with a great deal of variability. In the laboratory and in some pathological or age-related states, however, dynamics can show reduced complexity that is more readily described and modeled. Reduced heart rate complexity has both clinical and dynamical significance - it may provide warning of impending illness or clues about the dynamics of the heart's pacemaking system. In this paper, we describe simple and interesting heart rate dynamics that we have observed in premature human infants - reversible transitions to large-amplitude periodic oscillations - and we show that the appearance and disappearance of these periodic oscillations can be described by a simple mathematical model, a Hopf bifurcation.

Automated synchrogram analysis applied to heartbeat and reconstructed respiration

Phase synchronization between two weakly coupled oscillators has been studied in chaotic systems for a long time. However, it is difficult to unambiguously detect such synchronization in experimental data from complex physiological systems. In this paper we review our study of phase synchronization between heartbeat and respiration in 150 healthy subjects during sleep using an automated procedure for screening the synchrograms. We found that this synchronization is significantly enhanced during non-rapid-eye-movement (non-REM) sleep (deep sleep and light sleep) and is reduced during REM sleep. In addition, we show that the respiration signal can be reconstructed from the heartbeat recordings in many subjects. Our reconstruction procedure, which works particularly well during non-REM sleep, allows the detection of cardiorespiratory synchronization even if only heartbeat intervals were recorded.

Chaotic dynamics of cardioventilatory coupling in humans: effects of ventilatory modes

Cardioventilatory coupling (CVC), a transient temporal alignment between the heartbeat and inspiratory activity, has been studied in animals and humans mainly during anesthesia. The origin of the coupling remains uncertain, whether or not ventilation is a main determinant in the CVC process and whether the coupling exhibits chaotic behavior. In this frame, we studied sedative-free, mechanically ventilated patients experiencing rapid sequential changes in breathing control during ventilator weaning during a switch from a machine-controlled assistance mode [assist-controlled ventilation (ACV)] to a patient-driven mode [inspiratory pressure support (IPS) and unsupported spontaneous breathing (USB)]. Time series were computed as R to start inspiration (RI) and R to the start of expiration (RE). Chaos was characterized with the noise titration method (noise limit), largest Lyapunov exponent (LLE) and correlation dimension (CD). All the RI and RE time series exhibit chaotic behavior. Specific coupling patterns were displayed in each ventilatory mode, and these patterns exhibited different linear and chaotic dynamics. When switching from ACV to IPS, partial inspiratory loading decreases the noise limit value, the LLE, and the correlation dimension of the RI and RE time series in parallel, whereas decreasing intrathoracic pressure from IPS to USB has the opposite effect. Coupling with expiration exhibits higher complexity than coupling with inspiration during mechanical ventilation either during ACV or IPS, probably due to active expiration. Only 33% of the cardiac time series (RR interval) exhibit complexity either during ACV, IPS, or USB making the contribution of the cardiac signal to the chaotic feature of the coupling minimal. We conclude that 1) CVC in unsedated humans exhibits a complex dynamic that can be chaotic, and 2) ventilatory mode has major effects on the linear and chaotic features of the coupling. Taken together these findings reinforce the role of ventilation in the CVC process.

High-order synchronization, transitions, and competition among Arnold tongues in a rotator under harmonic forcing

We consider a rotator whose equation of motion for the angle theta consists of the zeroth and first Fourier modes. Numerical analysis based on the trailing of saddle-node bifurcations is used to locate the n:1 Arnold tongues where synchronization occurs. Several of them are wide enough for high-order synchronization to be seen in passive observations. By sweeping the system parameters within a certain range, we find that the stronger the dependence of theta[over ] on theta , the wider the regions of synchronization. Use of a synchronization index reveals a vast number of very narrow n:m Arnold tongues. A competition phenomenon among the tongues is observed, in that they "push" and "squeeze" one another: as some tongues widen, others narrow. Two mechanisms for transitions between different n:m synchronization states are considered: slow variation of the driving frequency, and the influence of low-frequency noise on the rotator.

Seidel-Herzel model of human baroreflex in cardiorespiratory system with stochastic delays

The stochastic versus deterministic solution of the Seidel-Herzel model describing the baroreceptor control loop (which regulates the short-time heart rate) are compared with the aim of exploring the heart rate variability. The deterministic model solutions are known to bifurcate from the stable to sustained oscillatory solutions if time delays in transfer of signals by sympathetic nervous system to the heart and vasculature are changed. Oscillations in the heart rate and blood pressure are physiologically crucial since they are recognized as Mayer waves. We test the role of delays of the sympathetic stimulation in reconstruction of the known features of the heart rate. It appears that realistic histograms and return plots are attainable if sympathetic time delays are stochastically perturbed, namely, we consider a perturbation by a white noise. Moreover, in the case of stochastic model the bifurcation points vanish and Mayer oscillations in heart period and blood pressure are observed for whole considered space of sympathetic time delays.

Nonlinear oscillator model reproducing various phenomena in the dynamics of the conduction system of the heart

A dedicated nonlinear oscillator model able to reproduce the pulse shape, refractory time, and phase sensitivity of the action potential of a natural pacemaker of the heart is developed. The phase space of the oscillator contains a stable node, a hyperbolic saddle, and an unstable focus. The model reproduces several phenomena well known in cardiology, such as certain properties of the sinus rhythm and heart block. In particular, the model reproduces the decrease of heart rate variability with an increase in sympathetic activity. A sinus pause occurs in the model due to a single, well-timed, external pulse just as it occurs in the heart, for example due to a single supraventricular ectopy. Several ways by which the oscillations cease in the system are obtained (models of the asystole). The model simulates properly the way vagal activity modulates the heart rate and reproduces the vagal paradox. Two such oscillators, coupled unidirectionally and asymmetrically, allow us to reproduce the properties of heart rate variability obtained from patients with different kinds of heart block including sino-atrial blocks of different degree and a complete AV block (third degree). Finally, we demonstrate the possibility of introducing into the model a spatial dimension that creates exciting possibilities of simulating in the future the SA the AV nodes and the atrium including their true anatomical structure.

Experimental evidence for phase synchronization transitions in the human cardiorespiratory system

Transitions in the dynamics of complex systems can be characterized by changes in the synchronization behavior of their components. Taking the human cardiorespiratory system as an example and using an automated procedure for screening the synchrograms of 112 healthy subjects we study the frequency and the distribution of synchronization episodes under different physiological conditions that occur during sleep. We find that phase synchronization between heartbeat and breathing is significantly enhanced during non-rapid-eye-movement (non-REM) sleep (deep sleep and light sleep) and reduced during REM sleep. Our results suggest that the synchronization is mainly due to a weak influence of the breathing oscillator upon the heartbeat oscillator, which is disturbed in the presence of long-term correlated noise, superimposed by the activity of higher brain regions during REM sleep.

Mechanism of cardioventilatory coupling: insights from cardiac pacing, vagotomy, and sinoaortic denervation in the anesthetized rat

Cardioventilatory coupling (CVC), a temporal alignment between the heartbeat and inspiratory activity, is a major determinant of breath-to-breath variation in observed respiratory rate (f(o)). The cardiac-trigger hypothesis attributes this to adjustments of respiratory timing by baroreceptor afferent impulses to the central respiratory pattern generator. A mathematical model of this hypothesis indicates that apparent CVC in graphical plots of ECG R wave vs. inspiratory time is dependent on the heart rate (HR), the rate of the intrinsic respiratory oscillator (f(i)), and the strength of the hypothetical cardiovascular afferent impulse. Failure to account for HR and f(i) may explain the inconsistent results from previous attempts to identify the neural pathways involved in CVC. Cognizant of these interactions, we factored in the HR-to-f(i) ratio in our examination of the role of the vagus nerve and arterial baroreceptors in CVC by cardiac pacing 29 anesthetized Sprague-Dawley rats and incrementally changing the HR. With the assumption of a relatively constant f(i), CVC could be examined across a range of HR-to-f(o) ratios before and after vagotomy, sinoaortic denervation, and vagotomy + sinoaortic denervation. We confirmed the relation between CVC, HR-to-f(o) ratio, and breath-to-breath respiratory period variability and demonstrated the loss of these relations after baroreceptor elimination. Sham experiments (n = 8) showed that these changes were not due to surgical stress. Our data support the notion that inspiratory timing can be influenced by cardiac afferent activity. We conclude that the putative cardiovascular input arises from the arterial baroreceptors and that the vagus nerve is not critical for CVC.

Empirical mode decomposition and synchrogram approach to cardiorespiratory synchronization

We use the empirical mode decomposition method to decompose experimental respiratory signals into a set of intrinsic mode functions (IMFs), and consider one of these IMFs as a respiratory rhythm. We then use the Hilbert spectral analysis to calculate the instantaneous phase of the IMF. Heartbeat data are finally incorporated to construct the cardiorespiratory synchrogram, which is a visual tool for inspecting synchronization. We perform analysis on 20 data sets collected by the Harvard medical school from ten young (21-34 years old) and ten elderly (68-81 years old) rigorously screened healthy subjects. Our results support the existence of cardiorespiratory synchronization. We also investigate the origin of the cardiorespiratory synchronization by addressing the problem of correlations between regularities of respiratory and cardiac signals. Our analysis shows that regularity of respiratory signals plays a dominant role in the cardiorespiratory synchronization.

Model for complex heart rate dynamics in health and diseases

A physiologically motivated, dynamical model of cardiovascular autonomic regulation is shown to be capable of generating long-range correlated and multifractal heart rate. Virtual disease simulations are carried out systematically to account for the disease-induced relative dysfunction of the parasympathetic and the sympathetic branches of the autonomic control. Statistical agreement of the simulation results with those of real life data is reached, suggesting the possible use of the model as a state-of-the-art basis for further understanding of the physiological correlates of complex heart rate dynamics.

Nonlinear statistical modeling and model discovery for cardiorespiratory data

We present a Bayesian dynamical inference method for characterizing cardiorespiratory (CR) dynamics in humans by inverse modeling from blood pressure time-series data. The technique is applicable to a broad range of stochastic dynamical models and can be implemented without severe computational demands. A simple nonlinear dynamical model is found that describes a measured blood pressure time series in the primary frequency band of the CR dynamics. The accuracy of the method is investigated using model-generated data with parameters close to the parameters inferred in the experiment. The connection of the inferred model to a well-known beat-to-beat model of the baroreflex is discussed.

Estimation of coupling between time-delay systems from time series

We propose a method for estimation of coupling between the systems governed by scalar time-delay differential equations of the Mackey-Glass type from the observed time series data. The method allows one to detect the presence of certain types of linear coupling between two time-delay systems, to define the type, strength, and direction of coupling, and to recover the model equations of coupled time-delay systems from chaotic time series corrupted by noise. We verify our method using both numerical and experimental data.

Inference of a nonlinear stochastic model of the cardiorespiratory interaction

We reconstruct a nonlinear stochastic model of the cardiorespiratory interaction in terms of a set of polynomial basis functions representing the nonlinear force governing system oscillations. The strength and direction of coupling and noise intensity are simultaneously inferred from a univariate blood pressure signal. Our new inference technique does not require extensive global optimization, and it is applicable to a wide range of complex dynamical systems subject to noise.

A quantitative comparison of different methods to detect cardiorespiratory coordination during night-time sleep

Background

The univariate approaches used to analyze heart rate variability have recently been extended by several bivariate approaches with respect to cardiorespiratory coordination. Some approaches are explicitly based on mathematical models which investigate the synchronization between weakly coupled complex systems. Others use an heuristic approach, i.e. characteristic features of both time series, to develop appropriate bivariate methods.

Objective

In this study six different methods used to analyze cardiorespiratory coordination have been quantitatively compared with respect to their performance (no. of sequences with cardiorespiratory coordination, no. of heart beats coordinated with respiration). Five of these approaches have been suggested in the recent literature whereas one method originates from older studies.

Results

The methods were applied to the simultaneous recordings of an electrocardiogram and a respiratory trace of 20 healthy subjects during night-time sleep from 0:00 to 6:00. The best temporal resolution and the highest number of coordinated heart beats were obtained with the analysis of 'Phase Recurrences'. Apart from the oldest method, all methods showed similar qualitative results although the quantities varied between the different approaches. In contrast, the oldest method detected considerably fewer coordinated heart beats since it only used part of the maximum amount of information available in each recording.

Conclusions

The method of 'Phase Recurrences' should be the method of choice for the detection of cardiorespiratory coordination since it offers the best temporal resolution and the highest number of coordinated sequences and heart beats. Excluding the oldest method, the results of the heuristic approaches may also be interpreted in terms of the mathematical models.

Critical scale invariance in a healthy human heart rate

We demonstrate the robust scale-invariance in the probability density function (PDF) of detrended healthy human heart rate increments, which is preserved not only in a quiescent condition, but also in a dynamic state where the mean level of the heart rate is dramatically changing. This scale-independent and fractal structure is markedly different from the scale-dependent PDF evolution observed in a turbulentlike, cascade heart rate model. These results strongly support the view that a healthy human heart rate is controlled to converge continually to a critical state.

Synchronization between main rhythmic processes in the human cardiovascular system

For the cases of spontaneous respiration and paced respiration with a fixed frequency and linearly increasing frequency, we investigate synchronization between three main rhythmic processes governing the cardiovascular dynamics in humans, namely, the main heart rhythm, respiration, and the process whose fundamental frequency is close to 0.1 Hz. The analysis of the experimental records reveals synchronous regimes of different orders n:m between all the three main rhythms. The influence of the regime of breathing and the magnitude of heart rate variability on the degree of synchronization is considered.

Regions of cardiorespiratory synchronization in humans under paced respiration

Cardiorespiratory synchronization under paced respiration is studied systematically as the respiration frequency is changed between 3 and 30 breaths per min. We plot a one-dimensional cut of the classical picture of synchronization regions along the line defining the current breathing amplitude. The existence of n:m synchronization regions of finite width is demonstrated for each of six subjects studied. The statistics of the different types of synchronization and their stability are discussed.