Respiratory sinus arrhythmia is reversed during positive pressure ventilation

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

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

Negative respiratory sinus arrhythmia (nRSA) in the MRI-scanner - a physiologic phenomenon observed during elevated anxiety in healthy persons

Recently, we reported on a rare manifestation of respiratory sinus arrhythmia (RSA), namely the "switched-off" RSA (Rassler et al., 2018), also called negative RSA (nRSA). It was found in a minority of healthy persons during elevated fMRI-related anxiety characterized by slow spontaneous breathing and synchronous slow beat-to-beat interval (RRI) oscillations. From 23 healthy scanner naïve participants of an fMRI study consisting of 4 resting states, we selected resting states with highest state anxiety (AS) from 10 participants (AS=24.6±2.5) and compared them to those with lowest AS of the same participants (AS=15.1±3.8, p<0.001). During elevated anxiety, the percentage of nRSA (nRSA%) was more than twice of RSA (p=0.045), while RSA prevailed during low anxiety. This indicates that nRSA might be related to elevated anxiety. Interestingly, nRSA was not only associated with slow RRI and breathing oscillations, but also occurred at "normal" breathing rates in the 0.20-0.35 Hz range. We often observed coupled RRI oscillations at 0.1 or 0.15 Hz and respiration at 0.3 Hz (rate ratio 1:3 or 1:2) with respiration-synchronous 0.3 Hz-wavelets in the RRI rhythm (termed "superposition") indicating a reduced dominance of the respiratory rhythm over the RRI rhythm. This novel finding is supported by the work of Perlitz et al., (2004) on a "0.15 Hz rhythm" in brainstem. The concept behind such a 1:n ratio is a pacemaker-like rhythm in the brainstem that "drives" the cardiac RRI signal and secondarily also respiration as reflected in the 1:n rate ratio.

Model-based causal closed-loop approach to the estimate of baroreflex sensitivity during propofol anesthesia in patients undergoing coronary artery bypass graft

Cardiac baroreflex is a fundamental component of the cardiovascular control. The continuous assessment of baroreflex sensitivity (BRS) from spontaneous heart period (HP) and systolic arterial pressure (SAP) variations during general anesthesia provides relevant information about cardiovascular regulation in physiological conditions. Unfortunately, several difficulties including unknown HP-SAP causal relations, negligible SAP changes, small BRS values, and confounding influences due to mechanical ventilation prevent BRS monitoring from HP and SAP variabilities during general anesthesia. We applied a model-based causal closed-loop approach aiming at BRS assessment during propofol anesthesia in 34 patients undergoing coronary artery bypass graft (CABG) surgery. We found the following: 1) traditional time and frequency domain approaches (i.e., baroreflex sequence, cross-correlation, spectral, and transfer function techniques) exhibited irremediable methodological limitations preventing the assessment of the BRS decrease during propofol anesthesia; 2) Granger causality approach proved that the methodological caveats were linked to the decreased presence of bidirectional closed-loop HP-SAP interactions and to the increased incidence of the HP-SAP uncoupling; 3) our model-based closed-loop approach detected the significant BRS decrease during propofol anesthesia as a likely result of accounting for the influences of mechanical ventilation and causal HP-SAP interactions; and 4) the model-based closed-loop approach found also a diminished gain of the relation from HP to SAP linked to vasodilatation and reduced ventricular contractility during propofol anesthesia. The proposed model-based causal closed-loop approach is more effective than traditional approaches in monitoring cardiovascular control during propofol anesthesia and indicates an overall depression of the HP-SAP closed-loop regulation.

Mechanical ventilation during anaesthesia: challenges and opportunities for investigating the respiration-related cardiovascular oscillations

The vast majority of the available literature regarding cardiovascular oscillations refers to spontaneously breathing subjects. Only a few studies investigated cardiovascular oscillations, and especially respiration-related ones (RCVO), during intermittent positive pressure mechanical ventilation (IPPV) under anaesthesia. Only a handful considered assisted IPPV, in which spontaneous breathing activity is supported, rather than replaced as in controlled IPPV. In this paper, we review the current understanding of RCVO physiology during IPPV, from literature retrieved through PubMed website. In particular, we describe how during controlled IPPV under anaesthesia respiratory sinus arrhythmia appears to be generated by non-neural mechano-electric feedback in the heart (indirectly influenced by tonic sympathetic regulation of vascular tone and heart contractility) and not by phasic vagal modulation of central origin and/or baroreflex mechanisms. Furthermore, assisted IPPV differs from controlled IPPV in terms of RCVO, reintroducing significant central respiratory vagal modulation of respiratory sinus arrhythmia. This evidence indicates against applying to IPPV interpretative paradigms of RCVO derived from spontaneously breathing subjects, and against considering together IPPV and spontaneously breathing subjects for RCVO-based risk assessment. Finally, we highlight the opportunities that IPPV offers for future investigations of RCVO genesis and interactions, and we indicate several possibilities for clinical applications of RCVO during IPPV.

Role of brainstem centers in cardiorespiratory phase difference during mechanical ventilation

During mechanical ventilation, large inter-patient and intra-patient variations of the phase of respiratory sinus arrhythmia (RSA) were described. To determine whether these variations were neurally mediated, we compared the RSA phase between: (1) 12 control subjects, (2) 23 mechanically ventilated patients without brain injury (MV group) and (3) 12 brain dead, mechanically ventilated patients, whose central nervous functions were abolished (BD group). ECG and ventilatory flow were recorded during 15 min and the RSA phase was then continuously computed by complex demodulation. Control group exhibited RSA phases between 180° and 250° whereas an opposite pattern, between 0° and 90°, was observed in the BD group. For the two groups, the phase was stable over time. In the MV group, the RSA phases were distributed between 0° and 260°, with a greater variability over time than the other groups. Therefore, during mechanical ventilation, brainstem centers may induce large variations of the RSA phase, not synchronous with the mechanical effect of ventilation.

Positive end-expiratory pressure may alter breathing cardiovascular variability and baroreflex gain in mechanically ventilated patients

Background

Baroreflex allows to reduce sudden rises or falls of arterial pressure through parallel RR interval fluctuations induced by autonomic nervous system. During spontaneous breathing, the application of positive end-expiratory pressure (PEEP) may affect the autonomic nervous system, as suggested by changes in baroreflex efficiency and RR variability. During mechanical ventilation, some patients have stable cardiorespiratory phase difference and high-frequency amplitude of RR variability (HF-RR amplitude) over time and others do not. Our first hypothesis was that a steady pattern could be associated with reduced baroreflex sensitivity and HF-RR amplitude, reflecting a blunted autonomic nervous function. Our second hypothesis was that PEEP, widely used in critical care patients, could affect their autonomic function, promoting both steady pattern and reduced baroreflex sensitivity.

Methods

We tested the effect of increasing PEEP from 5 to 10 cm H2O on the breathing variability of arterial pressure and RR intervals, and on the baroreflex. Invasive arterial pressure, ECG and ventilatory flow were recorded in 23 mechanically ventilated patients during 15 minutes for both PEEP levels. HF amplitude of RR and systolic blood pressure (SBP) time series and HF phase differences between RR, SBP and ventilatory signals were continuously computed by complex demodulation. Cross-spectral analysis was used to assess the coherence and gain functions between RR and SBP, yielding baroreflex-sensitivity indices.

Results

At PEEP 10, the 12 patients with a stable pattern had lower baroreflex gain and HF-RR amplitude of variability than the 11 other patients. Increasing PEEP was generally associated with a decreased baroreflex gain and a greater stability of HF-RR amplitude and cardiorespiratory phase difference. Four patients who exhibited a variable pattern at PEEP 5 became stable at PEEP 10. At PEEP 10, a stable pattern was associated with higher organ failure score and catecholamine dosage.

Conclusions

During mechanical ventilation, stable HF-RR amplitude and cardiorespiratory phase difference over time reflect a blunted autonomic nervous function which might worsen as PEEP increases.

Nicotinic receptors partly mediate brainstem autonomic dysfunction evoked by the inhaled anesthetic isoflurane

BACKGROUND

Isoflurane is one of the most commonly used volatile anesthetics, yet the cardiorespiratory depression that occurs with its use remains poorly understood. In this study, the author examined isoflurane modulation of postsynaptic gamma-aminobutyric acid (GABA) receptors in parasympathetic cardiac vagal neurons (CVNs) and alterations of GABAergic function by targeting nicotinic acetylcholine receptors on GABAergic presynaptic terminals.

METHODS

Rhythmic inspiratory-related activity was recorded from the hypoglossal rootlet of 800 microm medullary sections. CVNs were identified by retrograde fluorescent labeling, and GABAergic neurotransmission to CVNs were examined using patch-clamp electrophysiological techniques.

RESULTS

Isoflurane at concentrations of >50 microM significantly suppressed inspiratory bursting frequency, amplitude, and duration. Isoflurane dose-dependently decreased the frequency and increased the decay time of spontaneous GABAergic inhibitory postsynaptic currents (IPSCs) in CVNs. To test whether the inhibition of GABAergic activity to CVNs was mediated by presynaptic nicotinic receptors, the nicotinic antagonist, dihydro-beta-erythroidine in an alpha(4)beta(2)-selective concentration (3 microM), was used. Dihydro-beta-erythroidine (3 microM) prevented the isoflurane-evoked depression of spontaneous GABAergic IPSC frequency, yet isoflurane still increased the IPSC decay time.

CONCLUSIONS

These results suggest clinically relevant concentrations of isoflurane inhibit brainstem respiratory rhythmogenesis, prolong inhibitory GABAergic postsynaptic currents and reduce GABA activity in CVNs. The decrease of GABAergic IPSCs frequency is dependent upon inhibition of presynaptic alpha(4)beta(2) nicotinic receptors.

Breathing cardiovascular variability and baroreflex in mechanically ventilated patients

Heart rate and blood pressure variations during spontaneous ventilation are related to the negative airway pressure during inspiration. Inspiratory airway pressure is positive during mechanical ventilation, suggesting that reversal of the normal baroreflex-mediated pattern of variability may occur. We investigated heart rate and blood pressure variability and baroreflex sensitivity in 17 mechanically ventilated patients. ECG (RR intervals), invasive systolic blood pressure (SBP), and respiratory flow signals were recorded. High-frequency (HF) amplitude of RR and SBP time series and HF phase differences between RR, SBP, and ventilatory signals were continuously computed by Complex DeModulation (CDM). Cross-spectral analysis was used to assess the coherence and the gain functions between RR and SBP, yielding baroreflex sensitivity indices. The HF phase difference between SBP and ventilatory signals was nearly constant in all patients with inversion of SBP variability during the ventilator cycle compared with cycling with negative inspiratory pressure to replicate spontaneous breathing. In 12 patients (group 1), the phase difference between RR and ventilatory signals changed over time and the HF-RR amplitude varied. In the remaining five patients (group 2), RR-ventilatory signal phase and HF-RR amplitude showed little change; however, only one of these patients exhibited a RR-ventilatory signal phase difference mimicking the normal pattern of respiratory sinus arrhythmia. Spectral coherence between RR and SBP was lower in the group with phase difference changes. Positive pressure ventilation exerts mainly a mechanical effect on SBP, whereas its influence on HR variability seems more complex, suggesting a role for neural influences.

Postural-induced phase shift of respiratory sinus arrhythmia and blood pressure variations: insight from respiratory-phase domain analysis

The purpose of this study is to evaluate the multiple effects of respiration on cardiovascular variability in different postures, by analyzing respiratory sinus arrhythmia (RSA) and respiratory-related blood pressure (BP) variations for systolic BP (SBP), diastolic BP (DBP), and pulse pressure (PP) in the respiratory-phase domain. The measurements were conducted for 420 s on healthy humans in the sitting and standing positions, while the subjects were continuously monitored for heart rate and BP variability and instantaneous lung volume. The waveforms of RSA and respiratory-related BP variations were extracted as a function of the respiratory phase. In the standing position, the waveforms of the BP variations for SBP, DBP, and PP show their maxima at around the end of expiration (π rad) and the minima at around the end of inspiration (2 π rad), while the waveform of RSA is delayed by ∼0.35 π rad compared with the BP waveforms. On the other hand, in the sitting position, the phase of the DBP waveform (1.69 π rad) greatly and significantly ( P < 0.01) differs from that in the standing position (1.20 π rad). Also, the phase of PP is delayed and that of RSA is advanced in the sitting position ( P < 0.01). In particular, the phase shift of the DBP waveform is sufficiently large to alter whole hemodynamic fluctuations, affecting the amplitudes of SBP and PP variations. We conclude that the postural change associated with an altered autonomic balance affects not only the amplitude of RSA, but also the phases of RSA and BP variations in a complicated manner, and the respiratory-phase domain analysis used in this study is useful for elucidating the dynamic mechanisms of RSA.

Effects of hypercapnia and hypoxemia on respiratory sinus arrhythmia in conscious humans during spontaneous respiration

Normally, at rest, the amplitude of respiratory sinus arrhythmia (RSA) appears to correlate with cardiac vagal tone. However, recent studies showed that, under stress, RSA dissociates from vagal tone, indicating that separate mechanisms might regulate phasic and tonic vagal activity. This dissociation has been linked to the hypothesis that RSA improves pulmonary gas exchange through preferential distribution of heartbeats in inspiration. We examined the effects of hypercapnia and mild hypoxemia on RSA-vagal dissociation in relation to heartbeat distribution throughout the respiratory cycle in 12 volunteers. We found that hypercapnia, but not hypoxemia, was associated with significant increases in heart rate (HR), tidal volume, and RSA amplitude. The RSA amplitude increase remained statistically significant after adjustment for respiratory rate, tidal volume, and HR. Moreover, the RSA amplitude increase was associated with a paradoxical rise in HR and decrease in low-frequency-to-high-frequency mean amplitude ratio derived from spectral analysis, which is consistent with RSA-vagal dissociation. Although hypercapnia was associated with a significant increase in the percentage of heartbeats during inspiration, this association was largely secondary to increases in the inspiratory period-to-respiratory period ratio, rather than RSA amplitude. Additional model analyses of RSA were consistent with the experimental data. Heartbeat distribution did not change during hypoxemia. These results support the concept of RSA-vagal dissociation during hypercapnia; however, the putative role of RSA in optimizing pulmonary perfusion matching requires further experimental validation.

Beat-by-beat cardiovascular index to predict unexpected intraoperative movement in anesthetized unparalyzed patients: a retrospective analysis

OBJECTIVE

Unexpected intraoperative movement may be detrimental during delicate surgery. This study tested retrospectively an algorithm based on beat-by-beat circulatory variables (incorporated into a Cardiovascular depth of anesthesia index: CARDEAN in relationship to unexpected movement, and compared its performance to that of the electroencephalogram (EEG)-derived index: BIS-XP 4.0.

METHODS

40 ASA I or II patients presenting for knee surgery had EEG (BIS XP 4.0), beat-by-beat (Finapres) finger non-invasive blood pressure (BP), conventional brachial BP and electrocardiogram (EKG) monitors attached. Anesthesia was induced and maintained with propofol and remifentanil. Before incision, the propofol concentration was set to maintain BIS < 60. From incision to emergence, the anesthesiologist was denied access to BIS or Finapres. Anesthesia adjustment was titrated at the discretion of the anesthesiologist according to conventional signs only: brachial BP, EKG, eyelash reflex, movement. Occurrences of movement and eye signs (divergence of eyeballs, tears, corneal reflex, eyelash reflex) were observed. The CARDEAN algorithm was written retrospectively and tested vs. BIS.

RESULTS

11 movements occurred in 8 patients. CARDEAN > 60 predicted movement in 30% of the cases, 15 to 274 s before movement (sensitivity: 100%, specificity: 95%; relative operating curve ROC = 0.98; prediction probability pk = 0.98). BIS > 60 predicted movement in 19% of cases (sensitivity: 64%; specificity: 94%, ROC: 0.85, pk: 0.85).

CONCLUSION

Retrospectively, a cardiovascular index predicted unexpected intraoperative movements. Prospective validation is needed.

Investigation of altered heart rate variability, nonlinear properties of heart rate signals, and organ dysfunction longitudinally over time in intensive care unit patients

PURPOSE

To investigate longitudinally over time heart rate dynamics and relation with mortality and organ dysfunction alterations in patients admitted to a multidisciplinary intensive care unit.

METHODS

Data from 53 patients were used, with heart rate recorded from monitors and analyzed on a daily basis (every morning) for 600 seconds and sampling rate at 250 Hz, from admission to the intensive care unit until final discharge from the unit. Variance, which is a measure of heart rate variability; exponent alpha2; and approximate entropy (ApEn), which assess long-range correlations and periodicity within a signal, respectively; were measured and compared with every day Sequential Organ Failure Assessment Score (SOFA) and mortality.

RESULTS

Nonsurvivors had lower ApEn mean (greater periodicity in their signals) and minimum values compared to survivors (0.53 +/- 0.25 vs 0.62 +/- 0.23, P = .04; 0.24 +/- 0.23 vs 0.48 +/- 0.23, P = .01, respectively). Patients in better conditions with SOFA of less than 7 (mean value) had higher variance and ApEn (more variable, less periodic signals) than those with SOFA of 7 or higher (0.47 +/- 0.51 vs 0.10 +/- 0.65, P < .001; 0.67 +/- 0.28 vs 0.49 +/- 0.24, P < .001, respectively). The alpha2 exponent and variance were correlated with length of stay (r = 0.55, P = .02, and r = 0.53, P = .02, respectively) and minimum ApEn with mortality (r = 0.41, P = .01).

CONCLUSIONS

Loss of variability and increase in periodicity in heart rate of critically ill patients are linked with parallel deterioration of organ dysfunction and high mortality.

Vagal nerve activity contributes to improve the efficiency of pulmonary gas exchange in hypoxic humans

The aim of this study was to test our hypothesis that both phasic cardiac vagal activity and tonic pulmonary vagal activity, estimated as respiratory sinus arrhythmia (RSA) and anatomical dead space volume, respectively, contribute to improve the efficiency of pulmonary gas exchange in humans. We examined the effect of blocking vagal nerve activity with atropine on pulmonary gas exchange. Ten healthy volunteers inhaled hypoxic gas with constant tidal volume and respiratory frequency through a respiratory circuit with a respiratory analyser. Arterial partial pressure of O(2) (P(aO(2))) and arterial oxygen saturation (S(pO(2))) were measured, and alveolar-to-arterial P(O(2)) difference (D(A-aO(2))) was calculated. Anatomical dead space (V(D,an)), alveolar dead space (V(D,alv)) and the ratio of physiological dead space to tidal volume (V(D,phys)/V(T)) were measured. Electrocardiogram was recorded, and the amplitude of R-R interval variability in the high-frequency component (RRIHF) was utilized as an index of RSA magnitude. These parameters of pulmonary function were measured before and after administration of atropine (0.02 mg kg(-1)). Decreased RRIHF (P < 0.01) was accompanied by decreases in P(aO(2)) and S(pO(2)) (P < 0.05 and P < 0.01, respectively) and an increase in D(A-aO(2)) (P < 0.05). Anatomical dead space, V(D,alv) and V(D,phys)/V(T) increased (P < 0.01, P < 0.05 and P < 0.01, respectively) after atropine administration. The blockade of the vagal nerve with atropine resulted in an increase in V(D,an) and V(D,alv) and a deterioration of pulmonary oxygenation, accompanied by attenuation of RSA. Our findings suggest that both phasic cardiac and tonic pulmonary vagal nerve activity contribute to improve the efficiency of pulmonary gas exchange in hypoxic conscious humans.

Fractal ventilation enhances respiratory sinus arrhythmia

Background

Programming a mechanical ventilator with a biologically variable or fractal breathing pattern (an example of 1/f noise) improves gas exchange and respiratory mechanics. Here we show that fractal ventilation increases respiratory sinus arrhythmia (RSA) – a mechanism known to improve ventilation/perfusion matching.

Methods

Pigs were anaesthetised with propofol/ketamine, paralysed with doxacurium, and ventilated in either control mode (CV) or in fractal mode (FV) at baseline and then following infusion of oleic acid to result in lung injury.

Results

Mean RSA and mean positive RSA were nearly double with FV, both at baseline and following oleic acid. At baseline, mean RSA = 18.6 msec with CV and 36.8 msec with FV (n = 10; p = 0.043); post oleic acid, mean RSA = 11.1 msec with CV and 21.8 msec with FV (n = 9, p = 0.028); at baseline, mean positive RSA = 20.8 msec with CV and 38.1 msec with FV (p = 0.047); post oleic acid, mean positive RSA = 13.2 msec with CV and 24.4 msec with FV (p = 0.026). Heart rate variability was also greater with FV. At baseline the coefficient of variation for heart rate was 2.2% during CV and 4.0% during FV. Following oleic acid the variation was 2.1 vs. 5.6% respectively.

Conclusion

These findings suggest FV enhances physiological entrainment between respiratory, brain stem and cardiac nonlinear oscillators, further supporting the concept that RSA itself reflects cardiorespiratory interaction. In addition, these results provide another mechanism whereby FV may be superior to conventional CV.

Time domain correlation analysis of heart rate variability in preterm neonates

BACKGROUND AND AIM

A fuller understanding of the neural control mechanisms of heart rate during the early stages of human development would be of great value to obstetric and neonatal management. In this paper, we investigate the correlation between heart rate variability (HRV) and other physiological parameters such as blood pressure and respiration in preterm neonates with the aim of developing a numerical model to explain and predict heart rate variability.

STUDY DESIGN AND SUBJECTS

All the required data are readily available for premature babies who are routinely monitored while being nursed in intensive care, and we have collected large data sets for a random group of such neonates. For the quantitative analysis of the data, we have developed a time domain correlation method, which has a number of advantages over the more commonly used power spectral analysis. We have been able to study the dynamics of the different frequency components of HRV by this method.

RESULTS

Highly correlated behaviour of the different HRV components, previously observed in our work on fetal HRV, is also present in the neonate, with similar characteristic time constants. Furthermore, the correlation of high-frequency (HF) oscillations of HRV with respiration and that of low-frequency (LF) oscillations of HRV with blood pressure are demonstrated on timescales of a single oscillation. In neonates receiving artificial ventilation, the correlation between HRV and respiration depends on the type of ventilation involved and assumes opposite polarities for the two main types of equipment currently in use.

CONCLUSION

We demonstrate that it is possible to analyse HRV quantitatively by calculating the relative gains and characteristic time constants for the correlated parameters and components.

Phase-averaged characterization of respiratory sinus arrhythmia pattern

A method for the accurate time-domain characterization of respiratory sinus arrhythmia (RSA) pattern is presented and applied to two groups of healthy subjects to lay the baseline of RSA patterns and to underlay their features: response to standing, stability in successive recordings, and individuality of the shape of RSA pattern. RSA pattern is evaluated by selective averaging of heart rate (HR) changes from multiple respiratory cycles over the respiratory phase and represents the complete modulating function of HR by respiration. The RSA pattern is evaluated with free respiration and even in cases of severe arrhythmia. Estimation error is 6–8% in magnitude, phase resolution is 0.2 rad, and sensitivity margin for respiratory-related HR variability (HRV) components is 1%. RSA magnitude, phase lag, and expiration-to-inspiration time ratio are derived in addition to the entire pattern. In a group of 10 healthy young adults, a phase lag difference of 11.4 ± 8.5% (mean ± SD, P < 0.004) was observed between supine and standing postures, possibly ascribed to breathing mechanics. A second group of 15 healthy young adults at supine rest showed stability of the RSA pattern in successive recordings (several weeks apart) as well as individuality among subjects. This may suggest a nonscalar individual long-term index for cardiorespiratory coupling. The method is complementary to the existing statistical and spectral methods. It allows the complete characterization of the primary RSA components and may provide new insight into the effects of vagal activity and changes in clinical conditions.

Paradoxical respiratory sinus arrhythmia in the anesthetized rat

This study examines respiratory sinus arrhythmia (RSA) in the isoflurane-anesthetized rat. In fifteen female Sprague-Dawley rats, we recorded continuous ECG and respiratory airflow before and after bilateral vagotomy. RSA was assessed using power spectral analysis and by plotting the normalised changes in heart period as a function of the time during the respiratory cycle. Contrary to descriptions of RSA in conscious rats, we observed in all rats in the current study a 'reversed' pattern of RSA in which heart rate decelerated during inspiration. Elimination of vagal efferent fibres to the heart by vagotomy did not abolish the presence of reversed RSA suggesting that the pattern of heart period variation is not neural, and may be largely mechanical. Vagotomy altered breathing by increasing respiratory period, tidal volume, and the time to peak inspiratory flow. These changes did not alter the magnitude of RSA but reduced the latency period between inspiratory onset and the onset of respiratory related prolongation of heart period. Periods of positive pressure ventilation were associated with reversal of the inspiratory cardiac-deceleration pattern of RSA to resemble the more widely described pattern of inspiratory cardiac-acceleration. We conclude that RSA is not a suitable measure of vagal tone during anesthesia in the rat and reiterate the caution that needs to be taken when working with anesthetized experimental models of cardiac control.

Contribution of the respiratory rhythm to sinus arrhythmia in normal unanesthetized subjects during positive-pressure mechanical hyperventilation

The precise contribution of the CO2-dependent respiratory rhythm to sinus arrhythmia in eupnea is unclear. The respiratory rhythm and sinus arrhythmia were measured in 12 normal, unanesthetized subjects in normocapnia and hypocapnia during mechanical hyperventilation with positive pressure. In normocapnia (41 +/- 1 mmHg), the respiratory rhythm was always detectable from airway pressure and inspiratory electromyogram activity. The amplitude of sinus arrhythmia (138 +/- 21 ms) during mechanical hyperventilation with positive pressure was not significantly different from that in eupnea. During the same mechanical hyperventilation pattern but in hypocapnia (24 +/- 1 mmHg), the respiratory rhythm was undetectable and the amplitude of sinus arrhythmia was significantly reduced (to 40 +/- 5 ms). These results show a greater contribution to sinus arrhythmia from the respiratory rhythm during hypocapnia caused by mechanical hyperventilation than previously indicated in normal subjects during hypocapnia caused by voluntary hyperventilation. We discuss whether the respiratory rhythm provides the principal contribution to sinus arrhythmia in eupnea.

Macrocirculation is not the sole determinant of respiratory induced variations in the reflection mode photoplethysmographic signal

Photoplethysmography (PPG) is a non-invasive optical technique sensitive to variations in blood volume and perfusion in the tissue. Reflection mode PPG may have clinical advantages over transmission mode PPG. To improve clinical usefulness and further development of the reflection mode PPG, studies on factors that modify the signal are warranted. We studied the coherence between the respiratory induced intensity variations (RIIV) of the PPG signal and respiratory synchronous pressure variations in central venous pressure (CVP), peripheral venous pressure (PVP) and arterial blood pressure (ABP) during positive pressure ventilation on 12 patients under anaesthesia and on 12 patients with spontaneous breathing. During positive pressure ventilation the coherence between all signals was high. Inspiration was followed first by an increase in CVP, then by increases in ABP and PVP and lastly by RIIV indicating less back-scattered light. In spontaneously breathing patients the coherence was high, but the phases between the signals were changed. During inspiration, ABP decreased slightly before CVP, followed by a decrease in RIIV and PVP. The phase relation between RIIV and respiratory induced variation in macrocirculation changed with ventilatory mode, but not in a uniform way, indicating the influence of mechanisms other than macrocirculation involved in generating the RIIV signal.

Linear and nonlinear analysis of heart rate variability during propofol anesthesia for short-duration procedures in children

OBJECTIVE

To determine whether heart rate variability metrics provide an accurate method of monitoring depth of anesthesia, assessing the response to painful stimuli, and assessing neuroautonomic regulation of cardiac activity in children receiving propofol anesthesia for short-duration procedures.

DESIGN

Prospective, case series.

SETTING

Sixteen-bed pediatric intensive care unit, oncology unit, and endoscopy suite in a tertiary care children's hospital and ophthalmology examination rooms in an associated eye institute.

PATIENTS

Thirty-three pediatric patients undergoing propofol anesthesia for short procedures.

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

Heart rate variability metrics studied included mean, SD, low- and high-frequency power, detrended fluctuation analysis (represented by correlation coefficient, alpha), and approximate entropy. Compared with the initial anesthetized state, we found increased heart rate SD (3.17 +/- 1.31 vs. 7.05 +/- 0.26 bpm, p <.0001), heart rate low-frequency power (3.69 +/- 0.36 vs. 4.48 +/- 0.41 bpm(2)/Hz, p <.0001), heart rate low-/high-frequency ratio (1.47 +/- 0.26 vs. 1.26 +/- 0.24, p =.001), and heart rate alpha (1.12 +/- 0.24 vs. 1.35 +/- 0.21, p <.0001) during painful procedure. Mean heart rate (105.8 +/- 13.4 vs. 101.5 +/- 12.4 bpm, p =.005) and heart rate approximate entropy decreased with painful procedure (0.75 +/- 0.19 vs. 0.53 + 0.16, p <.001), whereas there was no significant change in heart rate high-frequency power (3.04 +/- 0.63 vs. 3.16 +/- 0.71 bpm(2)/Hz, p =.26).

CONCLUSIONS

We conclude that power spectral analysis of heart rate variability may be an accurate and clinically useful measure of depth of propofol anesthesia. We speculate that high-frequency heart rate power during propofol anesthesia correlates with depth of anesthesia, whereas low-frequency power allows for assessment of the patient's sympathetic response to pain.

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.

Linear and nonlinear analysis of hemodynamic signals during sepsis and septic shock

OBJECTIVE

Neuroautonomic modulation of heart rate (HR) and blood pressure were assessed in sepsis or septic shock. We hypothesized that these metrics would be diminished in pediatric patients with sepsis and septic shock, indicating uncoupling of the autonomic and cardiovascular systems.

DESIGN

Prospective case series.

SETTING

Pediatric intensive care unit in a tertiary care children's hospital.

PATIENTS

Thirty pediatric patients with sepsis or septic shock.

INTERVENTIONS

None.

MEASURES AND MAIN RESULTS

Metrics used included power spectral analysis, a linear frequency domain measure, and detrended fluctuation analysis, a nonlinear technique that assesses the degree of long-range correlation in HR or blood pressure. We found decreased low-frequency (2.68 +/- 0.24 vs. 3.37 +/- 0.17 [SEM] bpm2; p = .03) and high-frequency HR power (2.18 +/- 0.14 vs. 2.79 +/- 0.23 bpm2; p = .04) and increased detrended fluctuation analysis scaling exponent (1.22 +/- 0.06 vs. 1.00 +/- 0.07 bpm2; p = .02) in sepsis vs. shock patients, respectively. Compared with sepsis or shock, recovery was associated with increases in low-frequency (3.61 +/- 0.15 vs. 3.05 +/- 0.19 bpm2; p < .0001) and high-frequency HR power (3.11 +/- 0.15 vs. 2.50 +/- 0.22 bpm2; p < .0001).

CONCLUSIONS

We conclude that uncoupling of the autonomic and cardiovascular systems occurs over both short- and long-range time scales during sepsis, and the degree of uncoupling may help differentiate between sepsis, septic shock, and recovery states.

Transfer function analysis of the circulation in patients undergoing sevoflurane anesthesia

PURPOSE

The effects of sevoflurane anesthesia on the interactions between heart rate, blood pressure and respiration were assessed using transfer function analysis.

METHODS

Nine ASA 1 or 2 patients undergoing elective surgery were involved. They were paralysed and their lungs were mechanically ventilated during sevoflurane anesthesia. Instantaneous heart rate (IHR) from electrocardiogram, instantaneous lung volume (ILV) by respiratory inductive plethysmography and mean blood pressure (MBP) by arterial tonometry were obtained during conscious state, and 1MAC and 2MAC of sevoflurane anesthesia. Transfer function analysis for the relationships between ILV and IHR, ILV and MBP, MBP and IHR were made for five minute periods during which the respiratory rate was varied in a standardized fashion.

RESULTS

In awake patients transfer magnitudes for the relationships between ILV and IHR and between MBP and IHR in the 0.04-0.5Hz frequency band were 8.9 +/- 7.7 bpm x l(-1) and 0.95 +/- 0.44 bpm x mmHg(-1) respectively. Sevoflurane 2MAC decreased these values to 1.2 +/- 0.7 (P = 0.014) and 0.26 +/- 0.14 (P < 0.01) respectively, but phases were not affected. Neither transfer magnitudes nor phases between ILV and MBP were affected during sevoflurane anesthesia. Coherence for the relationships between ILV and IHR and between MBP and IHR were decreased during 1MAC sevoflurane anesthesia but not affected during 2MAC sevoflurane anesthesia.

CONCLUSIONS

The interactions between heart rate, blood pressure and respiration were altered by sevoflurane anesthesia. These findings could be explained by the attenuation of autonomic nervous system activity.

Decomplexification in critical illness and injury: relationship between heart rate variability, severity of illness, and outcome

OBJECTIVES

To determine if decomplexification of heart rate dynamics occurs in critically ill and injured pediatric patients. We hypothesized that heart rate power spectra, a measure of heart rate dynamics, would inversely correlate with measures of severity of illness and outcome.

DESIGN

A prospective clinical study.

SETTING

A 12-bed pediatric intensive care unit (ICU) in a tertiary care children's hospital.

PATIENTS

One hundred thirty-five consecutive pediatric ICU admissions.

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

We compared heart rate power spectra with the Pediatric Risk of Mortality (PRISM) score, the Pediatric Cerebral Performance Category (PCPC), and the Pediatric Overall Performance Category (POPC). We found significant negative correlations between minimum low-frequency and high-frequency heart rate power spectral values recorded during ICU stay and the maximum PRISM score (log low-frequency heart rate power vs. PRISM, r2 = .293, p < .001; and log high-frequency heart rate power vs. PRISM, r2 = .243, p < .001) and outcome at ICU discharge (log low-frequency heart rate power vs. POPC or PCPC, r2 = .429, p < .001; and log high-frequency heart rate power vs. POPC or PCPC, r2 = .271, p < .001).

CONCLUSIONS

Our data support the hypothesis that measures of heart rate power spectra are inversely related and negatively correlated to severity of illness and outcome in critically ill and injured children. The phenomenon of decomplexification of physiologic dynamics may have important clinical implications in critical illness and injury.

Contributions of tidal lung inflation to human R-R interval and arterial pressure fluctuations

We studied the effects of mechanical lung inflation on respiratory frequency R-R interval and arterial pressure fluctuations in nine healthy young adults undergoing elective orthopedic surgery. We conducted this research to define the contribution of pulmonary and thoracic stretch receptor input to respiratory sinus arrhythmia. We compared fast Fourier transform spectral power during three modes of ventilation: (1) spontaneous, frequency-controlled (0.25 Hz) breathing, (2) intermittent positive pressure ventilation (0.25 Hz, with a tidal volume of 8 ml/kg) and (3) high frequency jet ventilation (5.0 Hz, 2.5 kg/cm2), after sedation and vecuronium paralysis. Mean R-R intervals, arterial pressures and arterial blood gas levels were comparable during all three breathing conditions. Respiratory frequency systolic pressure spectral power was comparable during spontaneous breathing and conventional mechanical ventilation, but was significantly reduced during high frequency jet ventilation (P < 0.05). Respiratory frequency R-R interval spectral power (used as an index of respiratory sinus arrhythmia) declined dramatically with sedation and muscle paralysis (P < 0.05), but was greater during conventional mechanical, than high frequency jet ventilation (P < 0.05). These results suggest that although phasic inputs from pulmonary and thoracic stretch receptors make a statistically significant contribution to respiratory sinus arrhythmia, that contribution is small.

Respiratory sinus arrhythmia during anaesthesia: assessment of respiration related beat-to-beat heart rate variability analysis methods

Beat-to-beat heart rate variability analysis is a powerful tool for the diagnosis of neuropathy. Respiration-related heart rate variability (respiratory sinus arrhythmia, RSA) reflects the function of parasympathetic nervous system during spontaneous ventilation while awake. RSA is also claimed to monitor the depth of anaesthesia. Power spectrum analysis or various averaging techniques of the heart rate variability are usually applied. The current literature, however, does not usually interpret the ground rules and limitations of the method used, and this may sometimes lead to erroneous conclusions on the data. The aim of our study was to compare and analyse critically the performance of different methods of evaluating RSA during anaesthesia and positive pressure ventilation. Power spectrum analysis, the root mean square of the successive RR-interval difference (RMSSD), and two respiration related methods, RSA index and average phase RSA, were included in the comparison. To test these methods, 11 patients were anaesthetised with isoflurane and their lungs were ventilated mechanically with a frequency of 6 cycles min-1. Each patient received a bolus dose of atropine (20 micrograms kg-1) during the trial. Electrocardiogram, electroencephalogram and tracheal pressure signal from respirator were recorded and analyses were performed off-line. We demonstrated that general indices, such as RMSSD, may be strongly affected by heart rate level and other non-respiration related variations in heart rate. We also showed that the effect of unwanted fluctuations on RSA can be reduced with respiration dependent beat-to-beat methods. Furthermore we confirmed that in addition to the amplitude, also the pattern of respiratory sinus arrhythmia is of interest: the pattern is reversed in phase compared to spontaneous breathing while awake, as we have shown earlier. To analyse RSA during anaesthesia, we recommend the use of an average phase RSA method based on beat-to-beat variability that shows both the amplitude and pattern of RSA. Finally, no measure of RSA should be used without a presentation of the actual beat-to-beat heart rate curve.

Changes in heart rate variability during induction of anesthesia with fentanyl and midazolam

OBJECTIVE

The study was designed to evaluate changes in autonomic nervous system function during induction of anesthesia with fentanyl, midazolam, and pancuronium and to answer the question of dose-dependency of these effects.

DESIGN

Prospective, randomized.

SETTING

A university hospital.

PARTICIPANTS

Forty consecutive cardiac surgical patients.

INTERVENTIONS

Anesthesia was induced with fentanyl, midazolam, and pancuronium. The patients were assigned to four groups differing in dosages of fentanyl plus midazolam and speed of injection. Fentanyl, 7.5 micrograms/kg (group A), 12.5 micrograms/kg (group B), and 20.0 micrograms/kg (group C) plus midazolam, 0.075 mg/kg (group A), 0.125 mg/kg (group B), and 0.200 mg/kg (group C) were administered over 10 minutes; in group D, fentanyl, 7.5 micrograms/kg, and midazolam, 0.075 mg/kg, were administered within 1 minute.

MEASUREMENTS AND MAIN RESULTS

Heart rate variability (HRV) was measured using parameters in the time domain and the frequency domain. The comparison of preinduction HRV with the intra-anesthetic epochs did not show significant differences with respect to heart rate, coefficient of variation, and root mean squared successive differences. Spectral analysis showed significant reductions of power in the vasomotor band (0.01 to 0.05 Hz) and the low-frequency band (0.05 to 0.15 Hz) in all groups. Power in the high-frequency band (0.15 to 0.50 Hz) decreased slightly, but this did not reach the significance level. A dose dependency of these changes was found in the low-frequency band only.

CONCLUSIONS

Parameters of HRV suggest that induction with fentanyl, midazolam, and pancuronium decreases sympathetic but not parasympathetic autonomic system activity. The anesthetic induction technique's modulation of autonomic nervous system balance is better represented by means of spectral analysis than by analysis in the time domain. This modulation was largely independent of the doses administered and independent of the speed of injection.

Modulation of muscle sympathetic activity during spontaneous and artificial ventilation and apnoea in humans

Respiratory modulation of muscle sympathetic activity was compared in relaxed subjects breathing spontaneously and in anaesthetized and non-anaesthetized subjects ventilated artificially with intermittent positive pressure. Muscle sympathetic activity was recorded directly from the peroneal nerve using the microneurographic technique. Arterial pressure was monitored continuously either by finger-pulse photoplethysmography (Finapres) or intraarterially. Respiratory modulation of sympathetic activity, heart rate and arterial pressure was measured by averaging consecutive breaths to the ECG R-wave closest to the onset of inspiration. In relaxed subjects (n = 15) breathing quietly the averaged sympathetic activity was greatest during late expiration and the first half of inspiration and minimal after the peak of inspiration, after correcting for delays within the baroreflex loop. Systolic and diastolic pressures fell during inspiration. In anaesthetized or awake subjects ventilated artificially at normal tidal volumes the pattern of respiratory modulation of sympathetic activity was preserved but the changes in arterial pressure were reversed and respiratory sinus arrhythmia abolished. Ventilation with positive end-expiratory pressure (20 cmH2O) increased the overall level of sympathetic activity and enhanced the breath-to-breath modulation. We conclude that, although baroreceptors provide potent modulation of muscle sympathetic activity in humans, the inspiratory inhibition of sympathetic activity does not depend on an increase in arterial pressure and hence an increase in baroreceptor input.

Atropine abolishes electroencephalogram-associated heart rate changes without an effect on respiratory sinus arrhythmia during anaesthesia in humans

Heart rate fluctuates with the electroencephalogram burst suppression pattern during anaesthesia: increasing at burst onset and decreasing at suppression. Heart rate also oscillates with positive pressure ventilation. The effects of atropine on these heart rate changes were studied in 12 patients during isoflurane anaesthesia and positive pressure ventilation at a frequency of 6 cycles min-1. Four additional patients served as controls. A bolus dose of atropine (20 micrograms kg-1 intravenously) abolished the electroencephalogram-correlated heart rate changes; however, the amplitude of respiratory sinus arrhythmia was not changed after atropine. The control mechanism of the burst suppression pattern in electroencephalogram also affects parasympathetic heart rate control. The control mechanisms of respiratory sinus arrhythmia during anaesthesia with positive pressure ventilation differ from those during spontaneous breathing awake.