Transient interactions between blood pressure, respiration and heart rate in man

Categorizing the Role of Respiration in Cardiovascular and Cerebrovascular Variability Interactions

OBJECTIVE

Respiration disturbs cardiovascular and cerebrovascular controls but its role is not fully elucidated.

METHODS

Respiration can be classified as a confounder if its observation reduces the strength of the causal relationship from source to target. Respiration is a suppressor if the opposite situation holds. We prove that a confounding/suppression (C/S) test can be accomplished by evaluating the sign of net redundancy/synergy balance in the predictability framework based on multivariate autoregressive modelling. In addition, we suggest that, under the hypothesis of Gaussian processes, the C/S test can be given in the transfer entropy decomposition framework as well. Experimental protocols: We applied the C/S test to variability series of respiratory movements, heart period, systolic arterial pressure, mean arterial pressure, and mean cerebral blood flow recorded in 17 pathological individuals (age: 648 yrs; 17 males) before and after induction of propofol-based general anesthesia prior to coronary artery bypass grafting, and in 13 healthy subjects (age: 278 yrs; 5 males) at rest in supine position and during head-up tilt with a table inclination of 60.

RESULTS

Respiration behaved systematically as a confounder for cardiovascular and cerebrovascular controls. In addition, its role was affected by propofol-based general anesthesia but not by a postural stimulus of limited intensity.

CONCLUSION

The C/S test can be fruitfully exploited to categorize the role of respiration over causal variability interactions.

SIGNIFICANCE

The application of the C/S test could favor the comprehension of the role of respiration in cardiovascular and cerebrovascular regulations.

Mathematical modeling of the cardiovascular autonomic control in healthy subjects during a passive head-up tilt test

A mathematical model is proposed for the autonomic control of cardiovascular system, which takes into account two separated self-exciting sympathetic control loops of heart rate and peripheral vascular tone. The control loops are represented by self-exciting time-delay systems and their tone depends on activity of the aortic, carotid, and lower-body baroreceptors. The model is used to study the dynamics of the adaptive processes that manifest in a healthy cardiovascular system during the passive head-up tilt test. Computer simulation provides continuous observation of the dynamics of the indexes and variables that cannot be measured in the direct experiment, including the noradrenaline concentration in vessel wall and heart muscle, tone of the sympathetic and parasympathetic control, peripheral vascular resistance, and blood pressure. In the supine and upright positions, we estimated the spectral characteristics of the model variables, especially in the low-frequency band, and the original index of total percent of phase synchronization between the low-frequency oscillations in heart rate and blood pressure signals. The model demonstrates good quantitative agreement with the dynamics of the experimentally observed indexes of cardiovascular system that were averaged for 50 healthy subjects.

Heart rate variability and critical flicker fusion frequency changes during and after parachute jumping in experienced skydivers

PURPOSE

The purpose of this study was (1) to further explore the heart rate dynamics and assess a potential cardiovascular risk in response to 4000 m jumps in experienced skydivers; (2) to assess whether there is an impact of such jumps on skydivers' cortical arousal or not, which may impact their decision making processes.

METHOD

18 experienced skydivers performed successive jumps from a plane at 4000 m of height. Heart rate dynamics and cortical arousal were assessed by the use of heart rate variability and Critical Flicker Fusion Frequency (CFFF), respectively.

RESULTS

CFFF did not differ between the three measurement time points (p > 0.05). Mean heart rate increased during the jump (p < 0.001) and came back to pre-jump values after the jump (p < 0.001). Percentage of the differences of successive NN intervals greater than 50 ms (pNN50) decreased during the jump (p < 0.001) and kept lower values after the jump compared to pre-jump (p < 0.05). High-frequency power (HF) did not differ during the jump (p > 0.05) but decreased after the jump compared to both pre-jump (p < 0.01) and jump (p < 0.05). Sample entropy decreased during the jump (p < 0.001) and came back to pre-jump values after the jump (p > 0.05).

CONCLUSION

These results confirm a vagal input reduction associated with a rise of the sympathetic tone during the jump and suggests that the experienced skydiver is not exposed to a high cardiovascular risk. This study also shows that environmental stresses induced by free fall could not hamper the perceptual vigilance of experienced skydivers.

Utility of a nonlinear joint dynamical framework to model a pair of coupled cardiovascular signals

We have recently proposed a correlated model to provide a Gaussian mixture representation of the cardiovascular signals, with promising results in identifying rhythm disturbances. The approach provides a transformation of the data into a set of integrable Gaussians distributed over time. Looking into the model from a new joint modeling perspective, it is capable of assembling a filtered estimation, and can be used to derive temporal information of the waveforms. In this paper, we present a step-by-step derivation of the joint model putting correlation assumptions together to conclude a minimal joint description for a pair of ECG-ABP signals. We then probe novel applications of this model, including Kalman filter based denoising and fiducial point detection. In particular, we use the joint model for denoising and employ the denoised signals for pulse transit time (PTT) estimation. We analyzed more than 70 h of data from 76 patients from the MIMIC database to illustrate the accuracy of the algorithm. We have found that this method can be effectively used for robust joint ECG-ABP noise suppression, with mean signal-to-noise ratio (SNR) improvement up to 23.2 (12.0) dB and weighted diagnostic distortion measures as low as 2.1 (3.3)% for artificial (real) noises, respectively. In addition, we have estimated the error distributions for QT interval, systolic and diastolic blood pressure before and after filtering to demonstrate the maximal preservation of morphological features (ΔQT: mean ± std = 2.2 ± 6.1 ms; ΔSBP: mean ± std = 2.3 ± 1.9 mmHg; ΔDBP: mean ± std = 1.9 ± 1.4 mmHg). Finally, we have been able to present a systematic approach for robust PTT estimation (r = 0.98, p <; 0.001, mean ± std of error = -0.26 ± 2.93 ms). These findings may have important implications for reliable monitoring and estimation of clinically important features in clinical settings. In conclusion, the proposed framework opens the door to the possibility of deploying a hybrid system that integrates these algorithmic approaches for index estimation and filtering scenarios with high output SNRs and low distortion.

Welfare implication of measuring heart rate and heart rate variability in dairy cattle: literature review and conclusions for future research

Heart rate (HR) measurements have been used to determine stress in livestock species since the beginning of the 1970s. However, according to the latest studies in veterinary and behaviour-physiological sciences, heart rate variability (HRV) proved to be more precise for studying the activity of the autonomic nervous system. In dairy cattle, HR and HRV indices have been used to detect stress caused by routine management practices, pain or milking. This review provides the significance of HR and HRV measurements in dairy cattle by summarising current knowledge and research results in this area. First, the biological background and the interrelation of the autonomic regulation of cardiovascular function, stress, HR and HRV are discussed. Equipment and methodological approaches developed to measure interbeat intervals and estimate HRV in dairy cattle are described. The methods of HRV analysis in time, frequency and non-linear domains are also explained in detail emphasising their physiological background. Finally, the most important scientific results and potential possibilities for future research are presented.

Life-threatening arrhythmia verification in ICU patients using the joint cardiovascular dynamical model and a Bayesian filter

In this paper, a novel nonlinear joint dynamical model is presented, which is based on a set of coupled ordinary differential equations of motion and a Gaussian mixture model representation of pulsatile cardiovascular (CV) signals. In the proposed framework, the joint interdependences of CV signals are incorporated by assuming a unique angular frequency that controls the limit cycle of the heart rate. Moreover, the time consequence of CV signals is controlled by the same phase parameter that results in the space dimensionality reduction. These joint equations together with linear assignments to observation are further used in the Kalman filter structure for estimation and tracking. Moreover, we propose a measure of signal fidelity by monitoring the covariance matrix of the innovation signals throughout the filtering procedure. Five categories of life-threatening arrhythmias were verified by simultaneously tracking the signal fidelity and the polar representation of the CV signal estimations. We analyzed data from Physiobank multiparameter databases (MIMIC I and II). Performance evaluation results demonstrated that the sensitivity of the detection ranges over 93.50% and 100.00%. In particular, the addition of more CV signals improved the positive predictivity of the proposed method to 99.27% for the total arrhythmic types. The method was also used for false arrhythmia suppression issued by ICU monitors, with an overall false suppression rate reduced from 42.3% to 9.9%. In addition, false critical ECG arrhythmia alarm rates were found to be, on average, 42.3%, with individual rates varying between 16.7% and 86.5%. The results illustrate that the method can contribute to, and enhance the performance of clinical life-threatening arrhythmia detection.

Heart rate and blood pressure interactions during attempts to consciously raise or lower heart rate and blood pressure in normotensive subjects

This study investigated the interaction between heart rate (HR) and blood pressure (BP) during conscious control under visual biofeedback and background noise conditions. Normotensive volunteers were instructed to (i) decrease and (ii) increase HR (group A, n = 16) or BP (group B, n = 16). After instructions to lower HR or BP there was no significant change in HR or BP for either group. After instructions to raise HR, HR increased significantly (13.8 ± 1.3 beats min⁻¹, P < 0.0001) and BP did not change. However, following instructions to raise BP, both HR and BP increased significantly: systolic BP (5.2 ± 1.5 mmHg, P < 0.001), diastolic BP (3.5 ± 0.9 mmHg, P < 0.001) and HR (8.6 ± 1.3 beats min⁻¹, P < 0.0001). Biofeedback and background noise did not alter the relative change in HR or BP. In conclusion, normotensive subjects were unable to reduce BP or HR under conscious control. Subjects were able to increase both HR and BP, and voluntary increases in HR did not alter BP, while voluntary increases in BP also increased HR indicating distinct HR/BP interactions during conscious control.

Nonlinear dynamics of cardiovascular ageing

The application of methods drawn from nonlinear and stochastic dynamics to the analysis of cardiovascular time series is reviewed, with particular reference to the identification of changes associated with ageing. The natural variability of the heart rate (HRV) is considered in detail, including the respiratory sinus arrhythmia (RSA) corresponding to modulation of the instantaneous cardiac frequency by the rhythm of respiration. HRV has been intensively studied using traditional spectral analyses, e.g. by Fourier transform or autoregressive methods, and, because of its complexity, has been used as a paradigm for testing several proposed new methods of complexity analysis. These methods are reviewed. The application of time-frequency methods to HRV is considered, including in particular the wavelet transform which can resolve the time-dependent spectral content of HRV. Attention is focused on the cardio-respiratory interaction by introduction of the respiratory frequency variability signal (RFV), which can be acquired simultaneously with HRV by use of a respiratory effort transducer. Current methods for the analysis of interacting oscillators are reviewed and applied to cardio-respiratory data, including those for the quantification of synchronization and direction of coupling. These reveal the effect of ageing on the cardio-respiratory interaction through changes in the mutual modulation of the instantaneous cardiac and respiratory frequencies. Analyses of blood flow signals recorded with laser Doppler flowmetry are reviewed and related to the current understanding of how endothelial-dependent oscillations evolve with age: the inner lining of the vessels (the endothelium) is shown to be of crucial importance to the emerging picture. It is concluded that analyses of the complex and nonlinear dynamics of the cardiovascular system can illuminate the mechanisms of blood circulation, and that the heart, the lungs and the vascular system function as a single entity in dynamical terms. Clear evidence is found for dynamical ageing.

Multimodal signal processing for the analysis of cardiovascular variability

Cardiovascular (CV) variability as a primary vital sign carrying information about CV regulation systems is reviewed by pointing out the role of the main rhythms and the various control and functional systems involved. The high complexity of the addressed phenomena fosters a multimodal approach that relies on data analysis models and deals with the ongoing interactions of many signals at a time. The importance of closed-loop identification and causal analysis is remarked upon and basic properties, application conditions and methods are recalled. The need of further integration of CV signals relevant to peripheral and systemic haemodynamics, respiratory mechanics, neural afferent and efferent pathways is also stressed.

Low-frequency blood flow oscillations in congestive heart failure and after beta1-blockade treatment

Laser Doppler flowmetry (LDF) of forearm skin blood flow, combined with iontophoretically-administered acetylcholine and sodium nitroprusside and wavelet spectral analysis, was used for noninvasive evaluation of endothelial function in 17 patients newly diagnosed with New York Heart Association class II-III congestive heart failure (CHF). After 20+/-10 weeks' treatment with a beta(1)-blocker (Bisoprolol), the measurements were repeated. Measurements were also made on an age- and sex-matched group of healthy controls (HC). In each case data were recorded for 30 min. In HC, the difference in absolute spectral amplitude of LDF oscillations between the two vasodilators manifests in the frequency interval 0.005-0.0095 Hz (p<0.01); this difference is initially absent in patients with CHF, but appears following the beta(1)-blocker treatment (p<0.01). For HC, the difference between the two vasodilators also manifests in normalised spectral amplitude in 0.0095-0.021 Hz (p<0.05). This latter difference is absent in CHF patients and is unchanged by treatment with beta(1)-blockers. It is concluded that there are two oscillatory skin blood flow components associated with endothelial function. Both are reduced in CHF. Activity in the lower frequency interval is restored by beta(1)-blocker treatment, confirming the association between CHF and endothelial dysfunction but suggesting the involvement of two distinct mechanisms.

Complexity and spectral analysis of the heart rate variability dynamics for distant prediction of paroxysmal atrial fibrillation with artificial intelligence methods

OBJECTIVE

Paroxysmal atrial fibrillation (PAF) is a serious arrhythmia associated with morbidity and mortality. We explore the possibility of distant prediction of PAF by analyzing changes in heart rate variability (HRV) dynamics of non-PAF rhythms immediately before PAF event. We use that model for distant prognosis of PAF onset with artificial intelligence methods.

METHODS AND MATERIALS

We analyzed 30-min non-PAF HRV records from 51 subjects immediately before PAF onset and at least 45min distant from any PAF event. We used spectral and complexity analysis with sample (SmEn) and approximate (ApEn) entropies and their multiscale versions on extracted HRV data. We used that features to train the artificial neural networks (ANNs) and support vector machine (SVM) classifiers to differentiate the subjects. The trained classifiers were further tested for distant PAF event prognosis on 16 subjects from independent database on non-PAF rhythm lasting from 60 to 320 min before PAF onset classifying the 30-min segments as distant or leading to PAF.

RESULTS

We found statistically significant increase in 30-min non-PAF HRV recordings from 51 subjects in the VLF, LF, HF bands and total power (p<0.0001) before PAF event compared to PAF distant ones. The SmEn and ApEn analysis provided significant decrease in complexity (p<0.0001 and p<0.001) before PAF onset. For training ANN and SVM classifiers the data from 51 subjects were randomly split to training, validation and testing. ANN provided better results in terms of sensitivity (Se), specificity (Sp) and positive predictivity (Pp) compared to SVM which became biased towards positive case. The validation results of the ANN classifier we achieved: Se 76%, Sp 93%, Pp 94%. Testing ANN and SVM classifiers on 16 subjects with non-PAF HRV data preceding PAF events we obtained distant prediction of PAF onset with SVM classifier in 10 subjects (58+/-18 min in advance). ANN classifier provided distant prediction of PAF event in 13 subjects (62+/-21 min in advance).

CONCLUSION

From the results of distant PAF prediction we conclude that ANN and SVM classifiers learned the changes in the HRV dynamics immediately before PAF event and successfully identified them during distant PAF prognosis on independent database. This confirms the reported in the literature results that corresponding changes in the HRV data occur about 60 min before PAF onset and proves the possibility of distant PAF prediction with ANN and SVM methods.

Altered Phase Interactions between Spontaneous Blood Pressure and Flow Fluctuations in Type 2 Diabetes Mellitus: Nonlinear Assessment of Cerebral Autoregulation

Cerebral autoregulation (CA) is an important mechanism that involves dilation and constriction in arterioles to maintain relatively s cerebral blood flow in response to changes of systemic blood pressure. Traditional assessments of CA focus on the changes of cerebral blood flow velocity in response to large blood pressure fluctuations induced by interventions. This approach is not feasible for patients with impaired autoregulation or cardiovascular regulation. Here we propose a newly developed technique-the multimodal pressure-flow (MMPF) analysis, which assesses CA by quantifying nonlinear phase interactions between spontaneous oscillations in blood pressure and flow velocity during resting conditions. We show that CA in healthy subjects can be characterized by specific phase shifts between spontaneous blood pressure and flow velocity oscillations, and the phase shifts are significantly reduced in diabetic subjects. Smaller phase shifts between oscillations in the two variables indicate more passive dependence of blood flow velocity on blood pressure, thus suggesting impaired cerebral autoregulation. Moreover, the reduction of the phase shifts in diabetes is observed not only in previously-recognized effective region of CA (<0.1Hz), but also over the higher frequency range from ~0.1 to 0.4Hz. These findings indicate that Type 2 diabetes alters cerebral blood flow regulation over a wide frequency range and that this alteration can be reliably assessed from spontaneous oscillations in blood pressure and blood flow velocity during resting conditions. We also show that the MMPF method has better performance than traditional approaches based on Fourier transform, and is more sui for the quantification of nonlinear phase interactions between nonstationary biological signals such as blood pressure and blood flow.

Multimodal Pressure Flow Analysis: Application of Hilbert Huang Transform in Cerebral Blood Flow Regulation

Quantification of nonlinear interactions between two nonstationary signals presents a computational challenge in different research fields, especially for assessments of physiological systems. Traditional approaches that are based on theories of stationary signals cannot resolve nonstationarity-related issues and, thus, cannot reliably assess nonlinear interactions in physiological systems. In this review we discuss a new technique "Multi-Modal Pressure Flow method (MMPF)" that utilizes Hilbert-Huang transformation to quantify dynamic cerebral autoregulation (CA) by studying interaction between nonstationary cerebral blood flow velocity (BFV) and blood pressure (BP). CA is an important mechanism responsible for controlling cerebral blood flow in responses to fluctuations in systemic BP within a few heart-beats. The influence of CA is traditionally assessed from the relationship between the well-pronounced systemic BP and BFV oscillations induced by clinical tests. Reliable noninvasive assessment of dynamic CA, however, remains a challenge in clinical and diagnostic medicine.In this brief review we: 1) present an overview of transfer function analysis (TFA) that is traditionally used to quantify CA; 2) describe the a MMPF method and its modifications; 3) introduce a newly developed automatic algorithm and engineering aspects of the improved MMPF method; and 4) review clinical applications of MMPF and its sensitivity for detection of CA abnormalities in clinical studies. The MMPF analysis decomposes complex nonstationary BP and BFV signals into multiple empirical modes adaptively so that the fluctuations caused by a specific physiologic process can be represented in a corresponding empirical mode. Using this technique, we recently showed that dynamic CA can be characterized by specific phase delays between the decomposed BP and BFV oscillations, and that the phase shifts are significantly reduced in hypertensive, diabetics and stroke subjects with impaired CA. In addition, the new technique enables reliable assessment of CA using both data collected during clinical test and spontaneous BP/BFV fluctuations during baseline resting conditions.

Detection of Nonlinearity in Cardiovascular Variability Signals using Cyclostationary Analysis

A novel approach for detection of polynomial nonlinearity in the neuro-cardiovascular system based on cyclostationary analysis is presented. Metronome breathing is employed to provide a sinusoidal input to the neuro-cardiovascular system in which Heart Rate Variability (HRV) and Blood Pressure Variation (BPV) are considered as its outputs. The presence of new harmonics of the main respiratory rate in the HRV and BPV is investigated by using the concept of (self) phase and (self) frequency coupling. It is shown that a second order polynomial nonlinear system is actually involved in producing the HRV and BPV. The strength of this nonlinearity decreases with increasing the breathing rate.

Influence of different respiratory maneuvers on exercise-induced cardiac vagal inhibition

Physical exercise inhibits cardiac vagal activity. To study the relationship between heart rate (HR) and respiratory pattern, we applied the 4-s exercise test (4sET) and measured cardiac vagal index (CVI) in 30 healthy subjects who served as their own controls, using the standard plus three additional variations, essentially respiratory, of the original protocol: (a) a maximum inspiratory apnea of 16 s, of which 8 s were in the pre-exercise phase (4sET(insp)); (b) free respiratory pattern (4sET(unc)); and (c) maximum expiratory apnea of 12 s (4sET(exp)). The respective results were expressed by the following CVIs: CVI(insp), CVI(unc) and CVI(exp). CVI was determined in a continuous digital ECG recording through a specific ratio of two RR interval durations. The results [(mean +/- SEM)] for the four different maneuvers were as follows: CVI (1.56 +/- 0.05), CVI(insp) (1.55 +/- 0.05), CVI(unc) (1.63 +/- 0.05) and CVI(exp) (1.37 +/- 0.02). ANOVA-Bonferroni significant differences were only found between CVI(exp) and CVI(insp) (P = 0.009), CVI(unc) (P < 0.001) and CVI (P = 0.003). Dividing our sample in terciles according to CVI values, those with lower CVI, showed an attenuation of biphasic HR response after a 15 s maximum inspiratory apnea. We conclude that cardiac vagal reflex seems to be influencing the biphasic HR response modulation after a 12 s inspiratory apnea as described in the original protocol of 4sET, and this appears to be the option that best discriminates the cardiac vagal reflex, with less variability in the maneuvers when subjects are divided in terciles.

Comparison of the effects of propofol and midazolam on the cardiovascular autonomic nervous system during combined spinal and epidural anesthesia

STUDY OBJECTIVE

To investigate the effect of propofol and midazolam on cardiac autonomic nervous system (CANS) activity during combined spinal-epidural anesthesia.

DESIGN

Prospective, clinical study.

SETTING

Operating room of a university hospital.

PATIENTS

Forty ASA physical status I and II patients scheduled for knee surgery.

INTERVENTION

Patients were randomized to receive sedation with either propofol or midazolam.

MEASUREMENTS

Heart rate (HR), HR variability (HRV), systolic arterial pressure (SAP), and SAP variability (SAPV) were used for the analysis. These values were measured at the preanesthetic period, after intrathecal injection for spinal anesthesia, after sedation with propofol or midazolam, and just before the end of surgery with sedation. Cross-spectral analyses of the HR and SAP data were assessed to quantify the frequency-related coherence spectra and phase spectra.

MAIN RESULTS

Spinal anesthesia itself had no effect on power spectral changes in both groups. After sedation, as for HRV, high-frequency (HF) power (HF, 0.15-0.40 Hz) did not change, whereas low-frequency (LF) power (LF, 0.04-0.15 Hz) and LF/HF, an indicator of CANS balance, significantly decreased with propofol. Further, coherence in cross-spectra presented depression in the LF band area after sedation with propofol. Before the end of surgery with sedation, LF and LF/HF in both HRV and SAPV were correlated with age in those with propofol; however, scarce relation was observed in those who received midazolam.

CONCLUSIONS

Propofol was more potent than midazolam in causing CANS activity to be sympatholytic during combined spinal and epidural anesthesia and which was correlated with age only with propofol.

Effects of RR segment duration on HRV spectrum estimation

Although patterns of heart rate variability (HRV) hold considerable promise for clarifying issues in clinical applications, the inappropriate quantification and interpretation of these patterns may obscure critical issues or relationships and may impede rather than foster the development of clinical applications. The duration of the RR interval series is not a matter of convenience but a fine balance between two important issues: acceptable variance and stationarity of the time series on one hand, and acceptable resolution of the spectral estimate and reduced spectral leakage on the other. Further, in the standard short-term HRV analysis, it has been observed that the previous studies in HRV spectral analysis use a wide range of RR interval segment duration for spectral estimation by Welch's algorithm. The standardization of RR interval segment duration is also important for comparisons among studies and is essential for within-study experimental contrasts. In the present study, a comparative analysis for RR interval segment durations has been made to propose an optimal RR interval segment duration. Firstly a simulated signal was analyzed with Hann window and zero padding for the segment lengths of 1024, 512, 256 and 128 samples resampled at 4 Hz with 50% overlapping. Again, the above procedure was applied to RR interval series and it was concluded that segment length of 256 samples with 50% overlapping provides a smoothed spectral estimate with clearly outlined peaks in low- and high-frequency bands. This easily understandable and interpretable spectral estimate leads to a better visual and automated analysis, which is not only desirable in basic physiology studies, but also a prerequisite for a widespread utilization of frequency domain techniques in clinical studies, where simplicity and effectiveness of information are of primary importance.

Involvement of sympathetic nerve activity in skin blood flow oscillations in humans

We have used the wavelet transform to evaluate the time-frequency content of laser-Doppler flowmetry (LDF) signals measured simultaneously on the surfaces of free microvascular flaps deprived of sympathetic nerve activity (SNA), and on adjacent intact skin, in humans. It was thereby possible to determine the frequency interval within which SNA manifests itself in peripheral blood flow oscillations. The frequency interval from 0.0095 to 2 Hz was examined and was divided into five subintervals: I, approximately 0.01 Hz; II, approximately 0.04 Hz; III, approximately 0.1 Hz; IV, approximately 0.3 Hz; and V, approximately 1 Hz. The average value of the LDF signal in the time domain as well as the mean amplitude and total power in the interval from 0.0095 to 2 Hz and amplitude and power within each of the five subintervals were significantly lower for signals measured on the free flap (P < 0.002). The normalized spectral amplitude and power in the free flap were significantly lower in only two intervals: I, from 0.0095 to 0.021 Hz; and II, from 0.021 to 0.052 Hz (P < 0.05); thus indicating that SNA is manifested in at least one of these frequency intervals. Because interval I has recently been shown to be the result of vascular endothelial activity, we conclude that we have identified SNA as influencing blood flow oscillations in normal tissues with repetition times of 20-50 s or frequencies of 0.02-0.05 Hz.

Effect of immersion, submersion, and scuba diving on heart rate variability

BACKGROUND

Heart rate variability (HRV) describes the cyclic variations in heart rate and offers a non-invasive tool for investigating the modulatory effects of neural mechanisms elicited by the autonomic nervous system on intrinsic heart rate.

OBJECTIVE

To introduce the HRV concept to healthy volunteers under control conditions and during scuba diving. In contrast with more established manoeuvres, diving probably activates both the sympathetic and parasympathetic nervous system through various stimuli-for example, through cardiac stretch receptors, respiration pattern, psychological stress, and diving reflex. A further aim of the study was to introduce a measure for determining a candidate's ability to scuba dive by providing (a) standard values for HRV measures (three from the time domain and three from the frequency domain) and (b) physiological responses to a strenuous manoeuvre such as scuba diving.

METHODS

Twenty five trained scuba divers were investigated while diving under pool conditions (27 degrees C) after the effects of head out immersion and submersion on HRV had been studied.

RESULTS AND CONCLUSIONS

(a) Immersion under pool conditions is a powerful stimulus for both the sympathetic and parasympathetic nervous system. (b) As neither the heart rate nor the HRV changed on going from immersion to submersion, the parasympathetic activation was probably due to haemodynamic alterations. (c) All HRV measures showed an increase in the parasympathetic activity. (d) If a physiological HRV is a mechanism for providing adaptability and flexibility, diving should not provoke circulatory problems in healthy subjects. (e) Either a lower than normal HRV under control conditions or a reduction in HRV induced by diving would be unphysiological, and a scuba diving candidate showing such characteristics should be further investigated.

Modelling couplings among the oscillators of the cardiovascular system

A mathematical model of the cardiovascular system is simulated numerically. The basic unit in the model is an oscillator that possesses a structural stability and robustness motivated by physiological understanding and by the analysis of measured time series. Oscillators with linear couplings are found to reproduce the main characteristic features of the experimentally obtained spectra. To explain the variability of cardiac and respiratory frequencies, however, it is essential to take into account the rest of the system, i.e. to consider the effect of noise. It is found that the addition of noise also results in epochs of synchronization, as observed experimentally. Preliminary analysis suggests that there is a mixture of linear and parametric couplings, but that the linear coupling seems to dominate.

The cardiovascular system as coupled oscillators?

Based on physiological knowledge, and on an analysis of signals related to its dynamics, we propose a model of the cardiovascular system. It consists of coupled oscillators. Each of them describes one of the subsystems involved in the regulation of one passage of blood through the circulatory system. The flow of blood through the system of closed tubes-the blood vessels-is described by wave equations.

RR-arterial pressure variability relationships

Methodological aspects of a causal black-box model of heart period/arterial pressure interaction, arterial pressure closed-loop regulation and respiration effects on both heart period and arterial pressure are revisited in the "time" (more exactly heart beat count) domain. Parameters are estimated from experimental data (model identification) by means of multiple linear regressions of actual samples over the past ones. The elements composing either heart period or systolic arterial pressure variability are visualised as beat-by-beat series. Indexes describing the signal interactions, the loop properties and the spectral components of the variability series are consequently summarised. In 17 normal young volunteers, the analysis was carried out during active standing, rest, mild clinostatic pedalling exercise at 10%, 20%, and 30% of the maximum effort, and recovery. A negative effect of heart period changes on systolic arterial pressure of - 13.3 mm Hg/s was found at rest. This effect, though augmented by exercise, appeared insignificant in explaining arterial pressure variability. Arterial baroreflex was assessed by alphaart index which had a value of 5.18 mm Hg/ms at rest, 3.78 mm Hg/ms during active standing, and decreased progressively with exercise down to 0.55 mmHg/ms. The pressure regulation loop displayed a tendency to amplify disturbances at low frequency (around 0.1 Hz) 5.94 times at rest, augmented to 8.88 times during standing, 7.55 at 30% exercise. The first parameter of the pressure auto-regression was slightly higher than 1 at rest and even more augmented during standing, thus, indicating a tendency of arterial pressure perturbations to persist from one beat to the next. These mechanisms appear important in the genesis of low-frequency pressure waves. Nonetheless, the trace of different sources was evident in the regression residuals. Noticeably, during exercise it explained 10.16% of total heart period variability compared to 12.49% related to the low-frequency oscillations of closed-loops. The origin of high-frequency waves synchronous with respiration appeared miscellaneous as well. Arterial pressure appeared negligibly affected by heart period changes. Conversely, a limited effect of arterial pressure waves was found on heart period superimposed to a large effect of cardiopulmonary reflexes directly modulating the sinus node. In conclusion, both high-frequency and low-frequency waves are composite phenomena and a multi-channel analysis comparing heart period and arterial pressure variability yields a variety of figures assessing cardiovascular regulation and cardiorespiratory coupling.

Quantification of haemodynamic response to auditory stimulus in intensive care

Measuring effects of sensory stimuli on haemodynamics could provide information about the interplay between central and autonomous nervous system (ANS). However, ANS response to sensory stimulus has received little attention. In this paper we present a signal processing scheme to extract the responses of heart rate and systemic arterial pressure on auditory stimulus in intensive care patients (N=5). In short, the effect of mechanical ventilation is rejected by optimal linear modelling. Other disturbances are attenuated by filtering and efficient rejection of outlying sweeps of data. The results show identifiable responses on three out of five cases. The response characteristics may be explained by synchronisation of spontaneous variability in systemic arterial pressure to auditory stimulus.

Dynamic cardiocirculatory control during propofol anesthesia in mechanically ventilated patients

We evaluated dynamic cardiovascular control by spectral analytical methods in 20 young adults anesthetized with propofol (2.5 mg/kg, followed by continuous infusion of 0.1 mg/kg/min) and in an awake control group during cyclic stimulation of the carotid baroreceptors via sinusoidal neck suction at 0.2 Hz (baroreflex response mediated mainly by vagal activity) and at 0.1 Hz (baroreflex response mediated by vagal and sympathetic activity). During anesthesia and mechanical ventilation at 0.25 Hz, major underdampened hemodynamic oscillations occurred at 0.055 +/- 0.012 Hz. The response of RR intervals to baroreceptor stimulation at 0.2 Hz was markedly decreased during anesthesia (median of transfer function magnitude between neck suction and RR intervals 3% of the awake control). Blood pressure response to baroreceptor stimulation at 0.1 Hz was significantly decreased during anesthesia to 26% (systolic blood pressure), and 44% (diastolic blood pressure) of the awake control. There was a significant delay in baroreflex effector responses during anesthesia. Our results demonstrate a markedly depressed vagally mediated heart rate response and an impaired blood pressure response to cyclic baroreceptor stimulation during propofol anesthesia in mechanically ventilated patients. The disturbed baroreflex control is accompanied by an irregular dynamic behavior of cardiovascular regulation, indicating a decreased stability of the control system.

IMPLICATIONS

An irregular dynamic behavior of the cardiovascular control system, associated with an impaired baroreflex control of heart rate and blood pressure, can be observed during propofol anesthesia in mechanically ventilated subjects.

Spectral components of heart rate variability determined by wavelet analysis

Spectral components of heart rate variability (HRV) are determined in the time-frequency domain using a wavelet transform. Based on the finer estimation of low-frequency content enabled by the logarithmic resolution of the wavelet transform, corrections of spectral intervals, already defined by Fourier and model based methods, are proposed. The characteristic peaks between 0.0095 and 0.6 Hz are traced in time and four spectral intervals are defined, I (0.0095-0.021 Hz), II (0.021-0.052 Hz), III (0.052-0.145 Hz) and IV (0.145-0.6 Hz), within which peaks are located for all subjects included. These intervals are shown to be invariant regardless of the age and the state of the system. We also show that the frequency and power of the spectral components are related to age, AMI and particularly to type II diabetes mellitus.

Is fatigue in patients with multiple sclerosis related to autonomic dysfunction?

Time-dependent frequency decomposition of fluctuations in cardiovascular signals (heart rate [HR], blood pressure, and blood flow) provides noninvasive and quantitative evaluation of autonomic activity during transient and steady-state conditions. This method was applied during a change of position from supine to standing in patients with multiple sclerosis (MS) who experienced unexplained fatigue and in age-matched control subjects. No difference in response to standing, as reflected in the time domain parameters (mean HR, mean blood pressure, and mean blood flow), was observed between patients with MS and control subjects. Moreover, no difference was observed in very-low-frequency and low-frequency (related to sympathetic activity) content of HR, blood pressure, blood flow, or high-frequency content of HR (related to parasympathetic activity). The only spectral estimates that showed a significant difference between groups were the ratio of low-frequency to high-frequency content of HR and low-frequency content of HR normalized to total power. Both these parameters provide an estimate of the sympathovagal balance. A significant increase in these two estimates on standing was observed in control subjects only, indicating possible impairment of the sympathovagal balance response to standing in patients with MS who experienced fatigue. The authors observed a significant age dependence between close age subgroups, which occurred in the MS group only and was observed in some of the investigated spectral estimates that reflect vagal activity. Therefore, the authors assumed that age-related reduction in vagal activity occurred earlier in patients with MS who experienced fatigue. This reduction could also explain the lack of increase in the sympathovagal balance on standing. To validate this enhanced age dependence, further investigation should be performed in a larger group of subjects with a wider age range.

Wavelet analysis of oscillations in the peripheral blood circulation measured by laser Doppler technique

The wavelet transform technique, a time-frequency method with logarithmic frequency resolution, was used to analyze oscillations in human peripheral blood flow measured by laser Doppler flowmetry. The oscillations extended over a wide frequency scale and their periods varied in time. Within the frequency range studied, 0.0095-1.6 Hz, five characteristic oscillations were revealed, arising from both local and central regulatory mechanisms. After the insertion of endothelium-dependent and endothelium-independent vasodilators the spectra of blood flow markedly differed in the frequency interval 0.0095-0.02 Hz. In this way it was demonstrated that endothelial activity is a rhythmic process that contributes to oscillations in blood flow with a characteristic frequency of around 0.01 Hz. The study illustrates the potential of laser Doppler flowmetry combined with dynamical systems analysis for studies of both the micro- and macroscopic mechanisms of blood flow regulation in vivo.

Computer-controlled cardiopulmonary bypass increases jugular venous oxygen saturation during rewarming

BACKGROUND

Conventional roller pump apulsatile cardiopulmonary bypass (CPB) was compared with computer-controlled pulsatile bypass, which was designed to recreate biological variability (return of beat-to-beat variability in rate and pressure with superimposed respiratory rhythms). The degree of jugular venous oxygen saturation (SjvO2) less than 50% during rewarming from hypothermic CPB was compared for the two bypass techniques. An SjvO2 less than 50% during rewarming from hypothermic CPB is correlated with cognitive dysfunction in humans.

METHODS

Pigs were placed on CPB for 3 hours using a membrane oxygenator with alpha-stat acid-base management and arterial filtration. After baseline measurements and normothermic CPB, the animals were randomized to apulsatile CPB (n = 12) or computer-controlled pulsatile CPB (roller pump speed adjusted by an average of 2.9 voltage output modulations/s; n = 12). The animals were then cooled to a nasopharyngeal temperature of 28 degrees C. During rewarming to stable normothermic temperatures, SjvO2 was measured at 5-minute intervals. The mean and cumulative areas for an SjvO2 less than 50% were determined for all animals.

RESULTS

No between-group differences in temperature were noted during hypothermic CPB or during rewarming. The rate of rewarming was not different between groups. Mean arterial pressure, partial pressure of oxygen in arterial blood, and partial pressure of carbon dioxide in arterial blood also did not differ between groups. The hemoglobin concentration was within 0.4 g/dL between groups at all time periods. Mean pulse pressure was 10.0 +/- 4.8 mm Hg in the apulsatile CPB group and 20.7 +/- 5.2 mm Hg in the pulsatile CPB group (p = 0.0002; unpaired t test). Markedly greater mean and cumulative areas under the curve for SjvO2 less than 50% were seen with apulsatile CPB (164 +/- 209 versus 1.9 +/- 3.6% x min, p = 0.021; and 1,796 +/- 2,263 versus 23 +/- 45% x min, p = 0.020, respectively).

CONCLUSIONS

Computer-controlled pulsatile CPB was associated with significantly greater SjvO2 during rewarming from hypothermic CPB. Both the mean and cumulative areas under the curve for SjvO2 less than 50% exceeded a ratio of 75:1 for apulsatile versus computer-controlled pulsatile CPB. These experiments suggest that cerebral oxygenation was better preserved during rewarming from moderate hypothermia with computer-controlled pulsatile CPB, which returned biologic variability to the flow pattern.

Transcephalic electrical impedance in the study of cerebral circulation in a juvenile pig model

Transcephalic electrical impedance offers a technique for non-invasive, cot-side monitoring of neonatal cerebral circulation but the exact nature of the signal is somewhat ambiguous. The impedance signal is examined in an animal project where the ventilator settings are adjusted (20 min-1-10 min-1-40 min-1 for 10 min periods each) to produce circulatory changes. Six juvenile pigs are intubated, and ECG, arterial blood pressure, carotid flow (CF) by electromagnetic flowmeter and impedance are continuously monitored and stored on analogue tape. Cardiac output by thermodilution, blood oxygen (pO2) and carbon dioxide (pCO2) tensions are measured. ECG is converted to heart rate, mean blood pressure is integrated, and the high-frequency (1.50-4.00 Hz) component of the impedance signal delta Z is computed using autoregressive spectral estimation. Stroke volume, peripheral vascular resistance (PVR) and cerebral vascular resistance (CVR) are calculated. pCO2 and CF increase and pO2 decreases during hypoventilation. CF correlates positively with cardiac output, stroke volume, delta Z and pCO2, and negatively with pO2 and CVR. delta Z correlates positively with heart rate and cardiac output, and negatively with PVR and CVR. It is concluded that the impedance signal is related to the amount of blood transmitted to the brain by every beat of the heart, depending on the changes in both the systemic circulation and the cerebral vascular compliance.

Multivariate time-variant identification of cardiovascular variability signals: a beat-to-beat spectral parameter estimation in vasovagal syncope

In this paper a bivariate, time-variant model able to continuously measure the mutual interactions between heart rate and systolic blood pressure variability signals is presented. A recursive identification of the model parameters makes it possible to estimate, on a beat-to-beat basis, spectral low-frequency (LF) and high-frequency (HF) power, (LF/HF ratio) and cross-spectral (coherence and phase relationships between spectral peaks) indexes during nonstationary events. These indexes can be helpful in: 1) physiological study of autonomic nervous system mechanisms of cardiovascular control and 2) quantification and clinical evaluation of the neural and mechanical links between the two signals. In addition, an estimate of baroreceptive activation (alpha-gain) is continuously extracted. Before applying the model to cardiovascular signals, the reliability of the estimated parameters was tested on simulated signals. Subsequently, the model was applied to investigating vasovagal syncope episodes, aiming at the assessment of autonomic nervous system status and autonomic role in the dynamic phenomena which lead to syncope. The proposed model, which provides noninvasive beat-to-beat evaluation of the autonomic events, may be useful in the description of the syncopal episodes and in the comprehension of the complex physiological mechanisms of syncope.

Vasomotor instability preceding tilt-induced syncope: does respiration play a role?

This study aimed to determine whether alterations in cardiovascular dynamics before syncope are related to changes in spontaneous respiration. Fifty-two healthy subjects underwent continuous heart rate (HR), arterial blood pressure (BP), and respiratory measurements during 10-min periods of spontaneous and paced breathing (0.25 Hz) in the supine and 60 degrees head-up tilt positions. Data were evaluated by power spectrum and transfer function analyses. During tilt, 27 subjects developed syncope or presyncope and 25 remained asymptomatic. Subjects with tilt-induced syncope had significantly greater increases in low-frequency (0.04-0. 15 Hz) systolic BP, diastolic BP, and HR power during tilt than the asymptomatic subjects (P </= 0.01). This difference was present during spontaneous but not paced breathing. However, average tidal volume, respiratory rate, minute ventilation, proportion of breaths below 0.15 Hz, and low-frequency respiratory power during tilt did not differ between syncopal and nonsyncopal subjects. Transfer magnitudes between low-frequency respiration and BP, and between BP and interbeat interval, were also similar between groups. Thus vasomotor instability before syncope is not related to alterations in respiration or the cardiovagal baroreflex but may reflect oscillating central sympathetic outflow to the vasculature.

Autonomic response to change of posture among normal and mild-hypertensive adults: investigation by time-dependent spectral analysis

In this study, we applied the time-dependent spectral analysis approach (SDA) to investigate the autonomic changes occurring during a transition from supine to standing position (CP), in normal and unmedicated mild hypertensive subjects. The SDA method enables an accurate follow-up of the instantaneous changes in autonomic activity, even during the unsteady phase of the transition, where sudden changes in heart rate (HR) and arterial blood pressure (ABP) are observed. We were able to quantify the vagal withdrawal (reflected in the high frequency component of the time-dependent spectrum of HR fluctuations) in the immediate response to CP and the more slowly following sympathetic increase (reflected in the low frequency component of ABP). This general pattern was observed in both groups. In addition, our results identified an altered sympathetic response to CP in mild-hypertensives, as compared to normal adults. Their basal sympathetic activity is enhanced (higher mean HR and increased low frequency fluctuations in ABP) and their response to CP is reduced, as reflected only in the LF content of ABP fluctuations, relative to normals. No difference was observed in HR fluctuations, showing that there is no parasympathetically mediated alteration of the baroreflex control of HR in mild-hypertension.

Respiration response curve analysis of heart rate variability

A noninvasive study was conducted on intact awake humans to characterize the dynamic response of the heart to the vagus during slow-, comfortable-, and fast-paced respiration (8, 12, and 18 breaths/min), under both sitting and standing conditions. The respiration response curve (RRC) of respiration-associated vagal effects on the heart was estimated, and characteristics of entrainment and frequency dependence on respiration were demonstrated. It was shown that the degree of entrainment and magnitude of phase resetting decrease with increase of pacing rate from 8 to 18 breaths/min. Further, the RRC was examined by overlapping equivalent phase shifts in different respiration cycles. This examination of the RRC can help us not only to find the common pattern underlying the RRC during different respiration cycles but also to perceive its variation related to degree of entrainment.

Transcephalic electrical impedance provides a means for quantifying pulsatile cerebral blood volume changes following head-up tilt

Transcephalic electrical impedance (delta Z) was used to assess pulsatile cerebral blood volume changes while tilting nine premature (30-34 weeks) infants 20 degrees head up. High-frequency (1.50-4.00 Hz) delta Z variability decreased 27% while heart rate did not show any change. We would like to suggest that the variability of transcephalic electrical impedance analysed with spectral estimation seems to provide a means for quantifying alterations in cerebral circulation.

Linear multivariate models for physiological signal analysis: theory

The general linear parametric multivariate modelling concept is presented. This model combines a variety of different kinds of multivariate linear models. The concept of partial spectral analysis is derived from the general model. Some emphasis is laid on the causality demands of the model, and it is shown that the classic strictly-causal structure must be abandoned in order to utilise the modelling in many practical situations. Two special sub-class models are described in detail: the multivariate autoregressive model and the multivariate dynamic adjustment model. Furthermore, time-varying modelling is considered. The modelling of the real system is presented on a general level as a system identification cycle. The application of the methods to real physiological data is presented in the companion paper.

Heart rate variability as an assessment of cardiovascular status

In summary, HRV is a useful tool for cardiovascular risk stratification and assessment of sympathovagal balance. It has been shown to be a useful investigative tool in anesthesia to study autonomic balance and dysfunction but has limitations with regard to a depth-of-anesthesia monitor because of the confounding effects of multiple drugs and surgical stimuli. When assessing HRV, it is important to understand the limitations of the methodology. In particular, absolute power is not a surrogate of autonomic tone, but the relationship between the power in the different frequency bands may better reflect autonomic balance.

Application of time series spectral analysis theory: analysis of cardiovascular variability signals

The paper focuses on the most important application problems commonly encountered in spectral analysis of short-term (less than 10 min) recordings of cardiovascular variability signals (CVSs), critically analysing the different approaches to these problems presented in the literature and suggesting practical solutions based on sound theoretical and empirical considerations. The Blackman-Tukey (BT) and Burg methods have been selected as the most representative of classical and AR spectral estimators, respectively. For realistic simulations, 'synthetic' CVSs are generated as AR processes whose parameters are estimated on corresponding time series of normal, post-myocardial infarction and congestive heart failure subjects. The problem of resolution of spectral estimates is addressed, and an empirical method is proposed for model order selection in AR estimation. The issue of the understandability and interpretability of spectral shapes is discussed. The problem of non-stationarity and removing trends is dealt with. The important issue of identification and estimation of spectral components is discussed, and the main advantages and drawbacks of spectral decomposition algorithms are critically evaluated.

Spectral analysis of systolic blood pressure and heart rate oscillations related to respiration

1. Non-invasive assessment of short-term systolic blood pressure (SBP) and heart rate (HR) variability was obtained with the plethysmographic finger blood pressure measurement device. Respiration was measured with a respiratory inductive plethysmograph, which was calibrated prior to each study. 2. The effects of breathing pattern on the respiratory (high frequency, HF) component of the SBP or HR spectrum were analysed by Fourier transform. 3. Our quantification of the changes in the HF (respiratory) peak of the HR or SBP spectrum with changes in tidal volume (VT) and breathing frequency (BF) indicates that the modulus of this component may be predicted for any combination of depth and frequency of breathing. 4. The modulus of this HF component for the HR or SBP spectrum was linearly related to the respiratory sinus arrhythmia (RSA) or to the SBP oscillation related to respiration.

The effect of severe brainstem injury on heart rate and blood pressure oscillations

To determine whether an intact brainstem is essential for the generation of neurogenically mediated fluctuations of R-R intervals and blood pressure, three patients with cerebellar lesions causing severe brainstem compression or death, one patient with a large pontine infarct and one patient with a pontine haemorrhage, were studied. Time-frequency maps (based on a modified Wigner distribution) were constructed from blood pressure and R-R interval signals in these patients with brainstem injury and were compared with maps of normal control subjects. Low frequency sympathetically mediated rhythms (0.01-0.12 Hz) in systolic and diastolic pressure remained detectable but attenuated in patients with brainstem injury whereas there was an almost complete loss of normal R-R intervals rhythmicity over 0.01 to 0.5 Hz range. These data suggest that fluctuations in R-R intervals require an intact brainstem, whereas low frequency approximately 0.06 +/- 0.02 Hz blood pressure rhythms may be preserved by spinal sympathetic circuitry.

Model for the assessment of heart period and arterial pressure variability interactions and of respiration influences

A model which assesses the closed-loop interaction between heart period (HP) and arterial pressure (AP) variabilities and the influence of respiration on both is applied to evaluate the sources of low frequency (LF approximately 0.1 Hz) and high frequency (HF, respiratory rate approximately 0.25 Hz) in conscious dogs (n = 18) and humans (n = 5). A resonance of AP closed-loop regulation is found to amplify LF oscillations. In dogs, the resonance gain increases slightly during baroreceptor unloading (mild hypotension obtained with nitroglycerine (NTG) i.v. infusion, n = 8) and coronary artery occlusion ((CAO), n = 6), and it is abolished by ganglionic transmission blockade ((ARF), Arfonad i.v. infusion, n = 3). In humans, this gain is considerably increased by passive tilt. Different, possibly central, sources of LF oscillations are also evaluated, finding a strong rhythmic modulation of HP during CAO. At HF, a direct respiratory arrhythmia is dominant in dogs at control, while it is considerably reduced during CAO. On the contrary, in humans, a strong influence of respiration on AP is shown which induces a reflex respiratory arrhythmia. An index of the gain of baroreceptive response, alpha cl, was decreased by NTG and CAO, and virtually abolished by chronic arterial baroreceptive denervation (TABD, n = 4) and ARF.

Short-term variability of blood pressure and heart rate in borderline and mildly hypertensive subjects

Electrocardiogram and intra-arterial blood pressure were recorded in 96 men (aged 35 to 45 years) by the Oxford method over a 30-hour period. The study involved 33 normotensive, 29 borderline hypertensive, and 34 mildly hypertensive individuals, as assessed by the cuff method. Five-minute periods during sleep and with subjects in supine, sitting, and standing positions were extracted from the recordings for frequency domain analysis. Power spectrum density estimates of systolic blood pressure, diastolic blood pressure, and heart rate were calculated by an autoregressive method over the bandwidths of 0.02 to 0.075 (low-frequency), 0.075 to 0.15 (midfrequency), and 0.15 to 0.35 Hz (high-frequency), attributable to thermoregulatory, baroreceptor, and respiratory activity. No significant intergroup differences were observed at nighttime, but in different body positions the borderline hypertensive subjects frequently had either greater low-frequency variability or smaller midfrequency variability than the other groups. In this respect, the power spectra for systolic and diastolic blood pressures provided better statistical differentiation between the groups than those for heart rate. Furthermore, the borderline hypertensive subjects exhibited attenuated night-day changes in the low-frequency band for all time series. The results suggest that in borderline hypertension the baroreceptor oscillations are shifted to lower frequencies, presumably reflecting altered function of the sympathetic nervous system. In conclusion, spectral analysis of blood pressure variability for controlled test situations made it possible to detect differences in the cardiovascular regulatory systems between normotensive, borderline hypertensive, and mildly hypertensive individuals.

Respiration-synchronous fluctuations in stroke volume, heart rate and arterial pressure in humans

1. Simultaneous recordings of beat-to-beat left cardiac stroke volume (SV, pulsed ultrasound Doppler), mean arterial pressure (MAP) and heart rate (HR) were obtained in ten healthy young adults during spontaneous respiration at supine rest, before and after cholinergic blockade by atropine (0.035 mg kg-1). 2. Respiration-synchronous fluctuations in SV, HR, cardiac output (CO) and MAP were quantified by spectral analysis of the recordings of each of these variables. 3. Before atropine administration, respiration-synchronous fluctuations in HR and SV were prominent. The changes in HR and SV were inversely related and variation in SV was the main source of respiratory variability in CO. Respiration-synchronous fluctuations in MAP were mainly caused by variations in CO. 4. After cholinergic blockade, respiratory HR variations were eliminated, whereas the respiratory fluctuations in SV persisted. The fluctuations in CO and MAP increased. In this situation, mechanically induced variations in SV were not counteracted by inverse HR fluctuations and the influence on CO thus increased. 5. The main source of respiratory fluctuations in MAP in supine humans is thus variation in SV, while inverse, vagally mediated HR variations tend to reduce the fluctuations in CO and MAP.

Use of the frequency-tracking locus in estimating the degree of respiratory entrainment in preterm infants

In order to define the complex interactions between external stimuli and non-linear physiological systems, a technique (the frequency-tracking locus, FTL) was devised that describes the cycle-by-cycle changes in phase angle and amplitude between two signals. Qualitative assessment of the nature of interactions between the signals can be made by examining the FTL. Quantitation of the extent of entrainment of the spontaneous physiological rhythm is possible after deriving a numerical index (the path-length index, PLI) describing the departure of the system from a fully entrained state. The FTL was applied to the study of interactions between spontaneous respiratory effort and mechanical inflation in preterm newborn babies undergoing mechanical ventilation. Stable and unstable states of 1:1 interaction were noted while integer-ratio relationships were seen at low rates of mechanical ventilation. Stable states of entrainment corresponded to a PLI value near unity, and the value of PLI increased rapidly as interactions became unstable. The FTL may be used to describe complex interactions in physiological systems, and may be used as a guide to baby-ventilator matching during mechanical ventilation of the newborn.

Effect of breathing pattern on blood pressure and heart rate oscillations in humans

1. The relationships of respiratory sinus arrhythmia (RSA) and respiratory changes in systolic blood pressure (SBP) to tidal volume (VT) and breathing frequency (BF), were quantified during voluntary control of VT and BF in healthy subjects. 2. Respiration was measured non-invasively with a respiratory inductive plethysmograph, which was calibrated prior to each study while breathing through a pneumotachygraph. Finger arterial blood pressure was measured non-invasively by the Finapres. 3. Heart rate (HR) increased during inspiration, with a nearly fixed time delay for most VT and BF approximating 0.9 s. The magnitude of RSA increased with increases in VT and with decreases in BF. SBP decreased during inspiration, with a time delay which increased as BF decreased, resulting in a phase delay approximating 160 degrees. The magnitude of the inspiratory fall in SBP increased with increases in VT. Increased amplitudes of RSA and SBP variation occurred at the lowest BF, consistent with the possibility of interactions between respiratory-related influences and those due to 'slow waves' of vasomotor tone. 4. The present results are consistent with the conclusion that respiratory effects on SBP are caused by a mechanism other that simply changes in HR.

A time domain approach for the fluctuation analysis of heart rate related to instantaneous lung volume

This paper presents a time domain technique for estimating transfer characteristics from fluctuations of instantaneous lung volume (ILV) to heart rate (HR). An effective procedure for estimating the impulse response of HR to ILV is proposed. Pre- and post-processing procedures, including prefiltering of the HR signal, preenhancement of the high frequency content of the ILV signal, and post-filtering of the estimated impulse response, together with a random breathing technique, are shown to effectively reduce spurious transfer gain so as to get a stable estimate of the impulse response. Analysis of the data collected from fourteen healthy male subjects in various conditions revealed that there are three components in the impulse response: fast positive, delayed slow negative, and oscillatory. The effects of the autonomic blocking agents propranolol and atropine on these transfer characteristics are also described.

Periodic posture stimulation of the baroreceptors and the local vasomotor reflexes

In a supine subject, lowering of the foot from heart level to 50 cm below is known to stimulate the local reflex response and the baroreceptor outflow. We lowered and raised the leg of a supine subject periodically, with the leg stationary between movements (square wave). The Traube-Hering-Mayer wave (THM congruent to 0.1 Hz) was captured by or locked on to the leg movement over a certain frequency range, this is usually called the entrainment range. Square wave periodic leg movement in this manner on 10 male subjects, mean age 22 years, demonstrated that the THM frequency can be entrained. The lower limit of the entrainment bandwidth is 0.0841 (SD 0.0030) Hz and the upper limit is 0.1176 (SD 0.0013) Hz. Further examination showed that this phenomenon is independent of the breathing input. Comparison with the Traube-Hering-Mayer entrainment techniques of breathing and periodic neck suction using the Eckberg collar which stimulates the baroreceptors showed similar results. This work supports the hypothesis that the local reflex response and the baroreceptor outflow entrain the THM frequency.

Short-term variability of systolic blood pressure and heart rate in normotensive subjects

Short-term fluctuations in systolic blood pressure (SBP) and heart rate (HR) and their inter-relationship were analysed in a group of normotensive middle-aged men (n = 16) using a multivariate autoregressive modelling technique. This study is the first to evaluate the beat-to-beat variability of SBP and HR in a group of real normotensive subjects. Direct intra-arterial blood pressure was registered together with ECG using an ambulatory tape recording technique (the Oxford method). Power spectrum density estimated (PSD) were used as a measure of the variability. PSDs were calculated over 3-min periods for four basic physiological conditions: during sleep and in the supine, sitting and standing positions. The inter-relationship between the blood pressure and heart rate variabilities was analysed using a closed-loop model. In agreement with results presented earlier in the literature, the beat-to-beat variation in SBP and HR was concentrated in three typical power spectrum regions: the high-frequency (HF = 0.15-0.35 Hz) region (respiration), the mid-frequency (MF = 0.075-0.15 Hz) region (vasomotor oscillation) and the low-frequency (LF = 0.02-0.075 Hz) region (thermoregulation). The variability changes considerably between different situations, especially that of the MF region. The variability was most prominent in the MF region and in the standing position. The variability was generally smallest in the HF region and in sleep. The results also demonstrate that the beat-to-beat variability in SBP and HR can considerably affect one another.