help button home button
AJRCCM
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by JUBRAN, A.
Right arrow Articles by TOBIN, M. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by JUBRAN, A.
Right arrow Articles by TOBIN, M. J.
Am. J. Respir. Crit. Care Med., Volume 156, Number 4, October 1997, 1129-1139

Effect of Hyperoxic Hypercapnia on Variational Activity of Breathing

AMAL JUBRAN, BRYDON J. B. GRANT, and MARTIN J. TOBIN

The Division of Pulmonary and Critical Care Medicine, Edward Hines Jr. Veterans Administration Hospital and Loyola University of Chicago Stritch School of Medicine, Hines, Illinois, and State University of New York at Buffalo and VAMC, Buffalo, New York

    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Dysrhythmias of breathing occur in several clinical disorders, but their mechanistic basis is obscure. To understand their pathophysiology, factors responsible for the variability of breathing need to be defined. We studied the effect of hyperoxic hypercapnia (CO2) on the variational activity of breathing in 14 volunteers before and after delivering CO2 nonobstrusively via a plastic hood. Compared with air, CO2 increased the gross variability of minute ventilation (V I) and tidal volume (VT), and decreased that of inspiratory time (TI) and expiratory time (TE) (all p < 0.03). CO2 increased the autocorrelation coefficient at a lag of one breath for V I (p < 0.05), the number of consecutive breath lags having significant autocorrelation coefficients for V I and VT (both p < 0.01), and the cycle time of oscillations in V I (p = 0.03) and VT (p = 0.04). Uncorrelated random behavior constituted >=  80% of the variance of each breath component, correlated behavior represented 9 to 20%, and oscillatory behavior represented < 1% during both air and CO2. CO2 increased the correlated behavior of volume components, which was accompanied by development of low-frequency oscillations with a cycle time consistent with central chemoreceptor activation.

    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

We and others have shown that breath components display considerable breath-to-breath variability in healthy human subjects (1). This variability in breathing can be considered to be composed of a fixed part, namely, the mean of the entire breath series, and a variable part, which is the deviation of the magnitude of each breath from the mean. This variable deviation from the mean, in turn, can be considered to have a non-random (correlated) fraction and an uncorrelated random fraction. In a study of steady-state ventilation of eucapnic healthy subjects, we (8) performed autocorrelation analysis to determine the fraction of variational activity that was correlated on a breath-to-breath basis. The group mean autocorrelation coefficient at a lag of one breath for each of the primary breath components, tidal volume (VT), inspiratory time (TI), and expiratory time (TE), were significantly different from zero, indicating that the magnitude of each component in a given breath was positively correlated with the magnitude of the same component in the immediately preceding breath. The autocorrelation coefficients for each component remained significant for approximately three consecutive breaths, suggesting the presence of short-term memory (9).

The factors that determine the variational activity in breathing are largely unknown. Breathing is controlled by multiple feedback loops, and factors such as central neural networks and chemical stimulation are likely to be important in determining the extent of breath-to-breath variability of breathing, which may be important in the development of dyspnea (10). Although an enormous literature pertains to the influence of hypercapnia on mean changes in respiratory motor output, research on the effect of carbon dioxide (CO2) on the variability of breathing in humans is scant to non-existent. On the basis of studies in vagotomized animals, it has been suggested that hyperoxic hypercapnia may strengthen the interrelationship between breath components in neighboring breaths because of the dominant activation of the central chemoreceptors under these conditions (11). The central chemoreceptors have a prolonged time constant, and, thus, induced deviations from the mean level of ventilation should return more slowly to the mean level, leading to an increase in correlated behavior. Several investigators (2, 12) have postulated that CO2 stabilizes the respiratory controller and decreases the occurrence of ventilatory oscillations. This reasoning is partly based on the linear ventilatory response to an increase in CO2, indicating that the gain of the chemoreceptors is constant, and, therefore, a more stable breathing pattern is expected during hypercapnia (2). Indeed, irregularities of breathing have been reported to disappear with CO2 inhalation (2). In addition, CO2 may act as a stabilizing influence by raising the apneic threshold (13).

Ventilatory oscillations are thought to be mediated largely by the peripheral chemoreceptors (12). However, Waggener and colleagues (14) observed ventilatory oscillations in premature infants that had a cycle time considerably longer than that expected for oscillations originating in the peripheral chemoreceptors. They attributed these low-frequency ventilatory oscillations to central chemoreceptor activity, although the role of the latter was not systematically investigated. Ventilatory oscillations are particularly common in the setting of combined hypoxia and sleep (13). In a study of sleeping subjects exposed to hypoxia, Chapman and colleagues (15) found that CO2 gain, quantitated as the ventilatory response to hypercapnia (16), was higher in subjects who developed ventilatory oscillations than in those who remained free of oscillations. The role of CO2 gain as a determinant of this and other fractions of variational activity in breathing has not been systematically examined in awake subjects.

To investigate the role of central chemoreceptor stimulation on variational activity of breathing, we recorded ventilation in healthy volunteers before and after the induction of a steady-state increase in end-tidal CO2 tension (PETCO2). The study was conducted in the presence of hyperoxia since the latter diminishes peripheral chemoreceptor activity in humans, and, in this setting, step changes in PETCO2 are considered to represent a relatively pure central chemoreceptor stimulus (17). The subjects inhaled CO2 through an open-ended plastic hood, and ventilation was measured nonobtrusively with an inductive plethysmograph to avoid the confounding effects resulting from instrumentation attached to the face (18). We tested the following hypotheses: (1) hyperoxic hypercapnia increases the nonrandom fraction of variational activity in breathing; (2) the change in the nonrandom fraction with hyperoxic hypercapnia is a function of CO2 gain.

    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Subjects

Fourteen nonsmoking volunteers (12 men and two women) with a mean age of 35 yr (range, 26 to 44 yr) participated in this study. All had normal pulmonary function. Informed consent was obtained from all subjects, and the study was approved by the Human Studies Subcommittee of Edward Hines Jr. Veterans Administration Hospital.

Respiratory Inductive Plethysmograph

Ventilation was measured nonobtrusively with a respiratory inductive plethysmograph (Non-Invasive Monitoring Systems, Miami Beach, FL). The inductive plethysmograph was calibrated against spirometry by the least squares method (7, 18). Validation was checked during tidal breathing in the horizontal and semirecumbent positions. Inductive plethysmographic data were considered acceptable if they displayed a =< 10% difference from spirometry at both the beginning and the end of the study. Validation of the inductive plethysmograph against spirometry in the horizontal and semirecumbent positions revealed mean ± SD arithmetic differences of 3.5 ± 2.8% and 4.9 ± 3.1%, respectively, before the experiments and 4.2 ± 3.3 and 5.0 ± 3.9%, respectively, after the experiments. The signals from the inductive plethysmograph were recorded on a microprocessor system (Non-Invasive Monitoring Systems, Miami Beach, FL), which sampled the data at 20 Hz. On a breath-by-breath basis, the microprocessor continuously calculated the values of minute ventilation (VI), VT, TI, and TE. The data were processed by taking the duration of the respiratory cycle (Ttot) as the unit of time.

End-tidal CO2

PETCO2 was recorded from nasal prongs using an infrared CO2 analyzer (No. 4700; Ohmeda, Louisville, CO). The infrared analyzer was calibrated before each study. Continuous breath-by-breath recordings of PETCO2 were obtained and stored on a computer disk.

Inspired Gas Administration

To avoid the significant changes in breathing pattern induced by instrumentation attached to the face (19), an open-ended hood was employed for the delivery of CO2. The hood consisted of transparent material that was not in direct contact with the subject's head and did not produce uncomfortable feelings such as claustrophobia. The hood was loosely positioned over the subject's head and it was ventilated with the desired gas composition. Compressed gas cylinders of pure oxygen (O2) and pure CO2 were connected to a flowmeter (Fischer and Porter Co., Warminister, PA), and the outlet of the latter was connected to the hood.

Protocol

The study was conducted while the subject lay semirecumbent on a bed in a quiet room. After successful validation of the inductive plethysmograph, the plastic hood was positioned over the subject's head while he or she rested quietly for 15 min before the initiation of data recording. Subjects were instructed to remain motionless, keeping arms by the side, and to remain awake for the duration of the study; electroencephalography was not monitored. The subjects watched a single recording of a video tape player. Subjects first breathed room air through the hood for 60 min. Then, the mixture of O2 and CO2 flowing into the hood was adjusted to achieve an increase in PETCO2 of ~ 10 mm Hg from baseline; in reality, PETCO2 increased from 33.7 ± 3.4 to 46.6 ± 4.0 mm Hg. This level of PETCO2 was maintained for at least 60 min. Breathing pattern and PETCO2 were continuously recorded on a computer and stored on disk for later analysis (Figure 1).


View larger version (26K):
[in this window]
[in a new window]
 
Figure 1.   Time-series plot of end-tidal PCO2 (PETCO2) and minute ventilation in a subject breathing air and hyperoxic hypercapnia. Compared with air, hyperoxic hypercapnia (section between dashed vertical lines) produced an increase in PETCO2 (40 to 60 mm Hg) and in minute ventilation (5.3 to 20.1 L/min). Two sets of 700 consecutive breaths were selected during air breathing (15 to 66 min) and hyperoxic hypercapnia (111 to 144 min) for analysis.

Ventilatory Response Progressive Hypercapnia

On a separate day, the CO2 gain was assessed by measuring the ventilatory response to hypercapnia using a modification of the Read rebreathing technique (16). The subjects rebreathed a 7% CO2-93% O2 mixture from a 6-L anesthesia bag, and ventilation was measured by a pneumotachograph (Hans Rudolf; Kansas City, MO) attached to the inspiratory side of a mouthpiece. Expired gas was continuously sampled by an infrared CO2 analyzer (Datex, Helsinki, Finland). Signals from the recording amplifiers were digitized continuously at 50 Hz using a 16-bit analog-to-digital converter (DATAQ Instruments, Inc., Akron, OH) connected to a computer and stored on disk. To reduce random variability, three CO2 response curves (separated by a 15-min rest period) were obtained and averaged. In each subject, the ventilatory response to CO2 was linear with a high correlation coefficient (mean, 0.97 ± 0.01).

Data Analysis

Breath components in a breath series display breath-to-breath variability, which may consist of correlated, oscillatory, and random components. Autocorrelation and spectral analysis enable the determination of the relative magnitudes of these components.

Gross variability. The standard deviation, viz., the square root of the variance, for each breath component was calculated in each subject during air and hyperoxic hypercapnia. To ensure that the standard deviations approximated a Gaussian distribution, they were logarithmically transformed. Then, the log-transformed data for each breath component during air were compared with those during hyperoxic hypercapnia using paired t tests.

Autocorrelation analysis. Autocorrelation analysis was employed to determine what fraction of variational activity is correlated on a breath-to-breath basis (see APPENDIX). By its ability to extract correlated activity from data obscured by random noise, autocorrelation analysis can determine if there is a strong relationship between one breath and another at some interval (or lag) away. It also determines the relative strength of "short-term memory" for each breath component, i.e., a given breath is followed by several breaths that persist in similar magnitude, as if by "memory" (19). Although memory is the term used to describe this statistical relationship, it does not necessarily signify that the consecutive serial autocorrelation coefficients have a neural origin (8, 9).

Spectral analysis. Power spectral analysis can also be used to quantitate breath-to-breath variability in a breath component (20). The power spectrum expresses the variance of a signal as a function of frequency. The presence of a significant peak (see APPENDIX) in the spectrum indicates that some of the variance in the data is due to a periodic oscillation with a period equal to the inverse of the frequency of the peak. The area inscribed by the peak (amount of power) reflects the degree of variability in the signal resulting from fluctuations at that frequency (19).

Although spectral analysis and autocorrelation analysis are mathematically related, periodic oscillations as detected by spectral analysis can represent physiologic mechanisms other than autoregressive behavior (4); in particular, spectral analysis can reveal low-frequency (slow) oscillations of breath components that might be missed by autocorrelation analysis; it has also been shown (4) that both types of analysis should be done in order to correctly partition the total variability of breathing without corrupting the autocorrelation coefficients.

Fractionation of variational activity. The variational activity of breathing was partitioned into autoregressive, periodic, and white noise fractions employing the approach described and validated by Modarreszadeh and colleagues (4) (see APPENDIX). This model enables the quantification of each fraction and its contribution to the entire variational activity of breathing. Fractionation of variational activity allows the influence of hyperoxic hypercapnia on the composition of the total variational activity to be investigated.

    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Mean Change in Breathing Component

Mean values of the breath components in each subject during air and hyperoxic hypercapnia are shown in Figure 2. VI increased from 6.47 ± 1.20 (SD) L/min during air breathing to 20.66 ± 6.33 L/min during hyperoxic hypercapnia (p < 0.001); the respective values for VT were 0.47 ± 0.13 and 1.14 ± 0.27 L (p < 0.001); the respective values for TI were 1.7 ± 0.4 and 1.5 ± 0.3 s (p = 0.14); the respective values for TE were 2.6 ± 0.5 and 1.9 ± 0.4 (p < 0.001).


View larger version (16K):
[in this window]
[in a new window]
 
Figure 2.   Minute ventilation (V I), tidal volume (VT), inspiratory time (TI), and expiratory time (TE) during air and hyperoxic hypercapnia (CO2) in 14 subjects (mean of 700 breaths for each condition in each subject). Hyperoxic hypercapnia increased V I and VT, and decreased TE (p < 0.001 in all instances); TI tended to decrease (p = 0.14). Circles and bars represent group mean ± SD.

Gross Variability of Breath Components

Frequency histograms of VI, VT, TI, and TE during air and hyperoxic hypercapnia in a representative subject are shown in Figure 3. Standard deviations of the breath components in each subject during air and hyperoxic hypercapnia are shown in Figure 4. The standard deviation of VI for the group increased from 1.45 ± 0.83 L/min during air breathing to 3.07 ± 1.19 L/min during hyperoxic hypercapnia (p < 0.001); the respective values for VT were 0.143 ± 0.118 and 0.241 ± 0.146 L (p < 0.005); the respective values for TI were 0.46 ± 0.53 and 0.26 ± 0.14 (p < 0.025); the respective values for TE were 0.66 ± 0.38 and 0.35 ± 0.14 (p < 0.001).


View larger version (18K):
[in this window]
[in a new window]
 
Figure 3.   Frequency histogram of minute ventilation (V I), tidal volume (VT), inspiratory time (TI), and expiratory time (TE) during air (open symbols) and hyperoxic hypercapnia (closed symbols) in the same subject as in Figure 1. Hyperoxic hypercapnia broadened the histogram of V I and narrowed the histograms of TI and TE; it had no effect on the histogram of VT. The standard deviation for V I increased from 0.85 L/min during air breathing to 1.73 L/min during hyperoxic hypercapnia; respective values for VT were 0.10 and 0.09 L; respective values for TI were 0.43 and 0.17 s; respective values for TE were 0.65 and 0.22 s.


View larger version (18K):
[in this window]
[in a new window]
 
Figure 4.   Standard deviations (SDs) of minute ventilation (V I), tidal volume (VT), inspiratory time (TI), and expiratory time (TE) during air and hyperoxic hypercapnia (CO2) in 14 subjects (mean of 700 breaths for each component in each subject). Hyperoxic hypercapnia increased the SDs of V I (p < 0.001) and VT (p < 0.005) and decreased SDs of TI (p < 0.025) and TE (p < 0.001). Circles and bars represent group mean ± SD.

Nonrandom Variability Analysis

Autocorrelograms of VI during air and hyperoxic hypercapnia in our representative subject are shown in Figure 5. Hyperoxic hypercapnia increased the autocorrelation coefficient at a lag of one breath and the number of breath lags with significant (p < 0.01) serial correlations. Power spectra of VI during air and hyperoxic hypercapnia in our representative subject are shown in Figure 6. Although no oscillations were observed during air breathing, a low frequency oscillation was detected during hyperoxic hypercapnia.


View larger version (14K):
[in this window]
[in a new window]
 
Figure 5.   Autocorrelogram of minute ventilation in the same subject as in Figure 1 during air and hyperoxic hypercapnia (CO2). The autocorrelation coefficient at a lag of one breath increased from 0.34 during air breathing to 0.44 during hyperoxic hypercapnia. The number of breath lags with significant serial correlations (p < 0.01) increased from four during air breathing to 49 during hyperoxic hypercapnia. On the autocorrelograms, points lying outside the inner pair of isopleths are statistically different than zero, at p < 0.05, and outside the outer pair of isopleths, at p < 0.01.


View larger version (16K):
[in this window]
[in a new window]
 
Figure 6.   Power spectra of minute ventilation during air (dashed tracing) and hyperoxic hypercapnia (solid tracing) in the same subject as in Figure 1. A significant peak was not observed during air breathing, whereas one occurred at 0.02 cycle/breath during hyperoxic hypercapnia. The centroid frequency decreased from 0.11 cycle/breath for V I during air breathing to 0.05 cycle/breath during hyperoxic hypercapnia.

Autocorrelation analysis. The autocorrelation coefficients at a lag of one breath for each breath component during air and hyperoxic hypercapnia in the 14 subjects are shown in Figure 7. The autocorrelation coefficient for VI increased from 0.194 ± 0.163 during air breathing to 0.329 ± 0.212 during hyperoxic hypercapnia (p < 0.05); the respective values for VT were 0.302 ± 0.171 and 0.342 ± 0.157 (p = 0.41); the respective values for TI were 0.232 ± 0.167 and 0.278 ± 0.135 (p = 0.29); the respective values for TE were 0.286 ± 0.137 and 0.286 ± 0.144 (p = 0.99).


View larger version (14K):
[in this window]
[in a new window]
 
Figure 7.   Autocorrelation coefficients at lag of one breath for minute ventilation (V I), tidal volume (VT), inspiratory time (TI), and expiratory time (TE) during air and hyperoxic hypercapnia (CO2) in 14 subjects. The autocorrelation coefficient for V I increased during hyperoxic hypercapnia (p < 0.05), whereas those for VT, TI, and TE did not change. Circles and bars represent group mean ± SD.

The number of consecutive breath lags having significant (non-zero, p < 0.01) autocorrelation coefficients for VI increased from 2.7 ± 4.8 during air breathing to 22.5 ± 25.1 during hyperoxic hypercapnia (p < 0.01); the respective values for VT were 3.2 ± 4.2 and 15.1 ± 16.3 lags (p < 0.01); the respective values for TI were 3.7 ± 4.6 and 6.3 ± 6.3 lags (p = 0.17); the respective values for TE were 4.1 ± 4.9 and 4.4 ± 2.7 lags (p = 0.82).

Spectral analysis. The centroid frequency, i.e., the mathematically weighted median frequency of the entire spectrum (0.0 to 0.5 cycle/breath), decreased from 0.16 ± 0.09 cycle/ breath for VI during air breathing to 0.11 ± 0.08 cycle/breath during hyperoxic hypercapnia (p < 0.05); the respective values for VT were 0.14 ± 0.06 and 0.09 ± 0.07 cycle/breath (p < 0.05); the respective values for TI were 0.17 ± 0.07 and 0.13 ± 0.05 cycle/breath (p = 0.11); the respective values for TE were 0.14 ± 0.06 and 0.12 ± 0.05 cycle/breath (p = 0.17).

The number of subjects who displayed a significant peak, and the corresponding frequency and power of the significant oscillations for each breath component in the combined low- and high-frequency bands, are listed in Table 1. An increase in the number of subjects with significant peaks during hyperoxic hypercapnia was observed for VT only. Hyperoxic hypercapnia caused a significant shift in the oscillations for VI (p = 0.03) and VT (p = 0.04) to lower frequencies, whereas it had no effect on those of TI and TE. Hyperoxic hypercapnia tended to decrease the power of the oscillations for TI (p = 0.11) and TE (p = 0.09), whereas it had no effect on the power of the oscillations for VI and VT.

                              
View this table:
[in this window]
[in a new window]
 

TABLE 1

SPECTRAL CHARACTERISTICS OF SIGNIFICANT OSCILLATIONS DURING AIR AND CO2 BREATHING*

Fractionation of variational activity of breathing. The fractions of variational activity of breathing secondary to oscillatory behavior, correlated behavior, and uncorrelated random [w(n)] behavior for each breath component during air and hyperoxic hypercapnia are listed in Table 2. During both air and hyperoxic hypercapnia, uncorrelated random [w(n)] behavior constituted > 80% of the variance of each breath component, correlated behavior represented 9 to 20%, and oscillatory behavior represented < 1%.

                              
View this table:
[in this window]
[in a new window]
 

TABLE 2

FRACTIONATION OF VARIATIONAL ACTIVITY OF BREATH  COMPONENTS  DURING  AIR AND HYPEROXIC HYPERCAPNIA*

The effect of hyperoxic hypercapnia resulted in an increase in the total variational activity in VI (p < 0.01) and a decrease in the total variational activity in TE (p < 0.04). Hyperoxic hypercapnia decreased the uncorrelated random [w(n)] fraction of VI (p = 0.05), and it tended to increase the correlated fraction of VI (p = 0.06) (Table 2 and Figure 8). Hyperoxic hypercapnia produced a very small, although statistically significant, increase in the oscillatory fraction for VT (p = 0.03).


View larger version (53K):
[in this window]
[in a new window]
 
Figure 8.   The fractions of variational activity in minute ventilation (V I) caused by uncorrelated random [w(n)], correlated, and oscillatory behavior during air and hyperoxic hypercapnia (CO2) in 14 subjects. Hyperoxic hypercapnia caused a marked increase in the variance of V I, represented by an increase in the area of the pie, and it decreased the fraction of w(n) (p = 0.05) and tended to increase the correlated fraction (p = 0.06). (Note, we took minor license in altering the fractions listed in this figure so that their sum was 100%, whereas summation in Table 2 slightly exceeded 100% because of averaging.)

Ventilatory Response to Progressive Hypercapnia

CO2 gain, quantitated as the slope of the ventilatory response to CO2 rebreathing, was 4.29 ± 2.09 (range, 1.29 to 8.00) L/ min/mm Hg in the 14 subjects. The autocorrelation coefficient at a lag of one breath for any breath component during air breathing or hyperoxic hypercapnia was not related to CO2 gain (Table 3). The number of consecutive breath lags having significant (non-zero) autocorrelation coefficients for the breath components was not related to CO2 gain, with the exception of a negative correlation with TI (r = -0.52, p = 0.05) during hyperoxic hypercapnia---that is, the number of significant lags were increased in subjects with a low CO2 gain. Likewise, the fractions of variational activity in each breath component during air and hyperoxic hypercapnia were not related to CO2 gain (Table 4).

                              
View this table:
[in this window]
[in a new window]
 

TABLE 3

RELATIONSHIP BETWEEN CO2 GAIN AND AUTOCORRELATION  COEFFICIENT AND NUMBER OF SIGNIFICANT  LAGS OF BREATH COMPONENTS*

                              
View this table:
[in this window]
[in a new window]
 

TABLE 4

RELATIONSHIP BETWEEN CO2 GAIN AND FRACTIONATION  OF VARIATIONAL ACTIVITY OF BREATH COMPONENTS  DURING AIR AND CO2 BREATHING

We were surprised by the lack of a relationship between measurements of the variational activity in breathing and CO2 gain. We wondered if this occurrence might be due to the fact that the variational activity of breathing was assessed under steady-state conditions, whereas CO2 rebreathing was not. Also, the CO2 rebreathing maneuver and the steady-state measurements were not performed on the same day. Accordingly, we derived a quasi-steady-state measurement of CO2 gain in each individual subject by taking the difference in mean VI between air and hyperoxic hypercapnia and dividing it by the mean change in PETCO2 between air and hyperoxic hypercapnia. This quasi-steady-state measurement of CO2 gain was negatively correlated with the autocorrelation coefficient at a lag of one breath for VI during hyperoxic hypercapnia (r = -0.52, p = 0.05). None of the other measurements of variational activity of any breath component was related to this measurement of CO2 gain.

    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Hyperoxic hypercapnia increased the gross variability of VI and VT, and decreased that of TI and TE. It increased the autocorrelation coefficient at a lag of one breath for VI and the number of consecutive breath lags having sufficient autocorrelation coefficients for VI and VT. Hyperoxic hypercapnia increased the number of subjects displaying oscillations in VT and it produced about a three-fold increase in the cycle time of oscillations in VI and VT. Fractionation revealed that at least 80% of the total variance of each breath component was due to uncorrelated random [w(n)] behavior, correlated behavior represented 9 to 20%, and oscillatory behavior constituted substantially less than 1% during both air and hyperoxic hypercapnia. Hyperoxic hypercapnia decreased the fraction of variational activity in VI because of noise, largely by increasing its fraction of correlated behavior.

Significance of CO2 Stimulus

In human studies, investigators have concluded that hypercapnic stimulation in the presence of hyperoxia causes central chemoreceptor stimulation with little or no contribution of the peripheral chemoreceptors (17, 21, 22); indeed, hyperoxic hypercapnia has been regarded as "an almost pure intracranial chemoreceptor stimulus" (17). Ledlie and colleagues (23) examined the effects of hyperoxic hypercapnia in an "open loop" system in dogs, in which respiratory motor output was prevented by paralysis from affecting ventilation and, thus, influencing the input chemical signals. They found that phrenic nerve activity and TE responses to hypercapnia were the same whether or not the carotid sinus nerves were severed. This further supports the notion that the response to hyperoxic hypercapnia stimulation is mediated by the central chemoreceptors.

Effect of Hypercapnia on Mean Change and Gross Variability of Breathing

CO2 produced increases in VI and VT, a decrease in TE, and no change in TI. Other investigators (24) have reported findings qualitatively similar to ours, although they used instrumentation requiring the use of a mouthpiece and noseclips--- factors that produce marked alterations in the pattern of breathing (18).

In a study that focused primarily on the interrelationships between VT, TI, and TE during air breathing, Newsom-Davis and Stagg (5) studied six healthy subjects who inhaled mixtures of 1.5 and 3% CO2 in air. Unfortunately, these investigators provided few details on this aspect of their investigation, except to state that the coefficients of variation of VT, TI, and TE tended to decrease in one subject, although numerical values were not given (see Figure 5 of their report). However, coefficients of variation are influenced by a change in the mean value (as occurs with CO2) and thus imperfectly reflect the true variability of breath components. We quantitated gross variability in terms of standard deviations, which increased for the volume components and decreased for the time components during hyperoxic hypercapnia (Figure 4). The increase in gross variability of VI with hyperoxic hypercapnia was also supported by inspection of the frequency histogram plots, which exhibited a broader dispersion of values with hyperoxic hypercapnia in 13 of the 14 study subjects. Moreover, the absolute variance of VI because of correlated behavior increased from 0.20 L/min during air breathing to 1.57 L/min during hyperoxic hypercapnia (p < 0.03), and the absolute variance of VI because of random behavior increased from 2.52 L/min during air breathing to 7.92 L/min during hyperoxic hypercapnia (p < 0.009). We are not aware of other systematic studies of gross variability of breath components in awake human subjects during steady-state hypercapnic stimulation with which to compare our results.

Interrelationships between Breath Components

The nonrandom fraction of breathing variability was first examined by Priban (6). Using the analysis of "run lengths" as a measure of randomness of data, Priban found that the average run length was 1.79 breaths compared with a value of 1.50 breaths per run that is expected for a random series. In four rebreathing experiments, he reported that CO2-induced hyperpnea (level of PCO2 or VI not stated) did not alter the number of breaths per run. A similar experimental approach was employed by Bolton and Marsh (25) in three healthy men. They also reported that hyperoxic hypercapnia (magnitude not stated) did not alter the run length. These results suggest that hypercapnia does not influence the overall variability of breathing; however, the investigators did not attempt to determine if correlated or oscillatory behavior was altered by CO2.

Khatib and colleagues (11) fitted a first-order autoregressive model to consecutive values of VT in anesthetized, vagotomized spontaneously breathing rats. While breathing air, five of seven rats had negative autoregressive coefficients, signifying that when VT of a given breath had a value above the mean, the value of the next breath was likely to be below the mean. In contrast, none of our subjects---nor in our previous study of 33 subjects (8)---had negative autocorrelation coefficients for VT while breathing air. When the rats breathed 4% CO2 in O2, however, three of three had positive autoregressive coefficients with an estimated average of ~ 0.29 (extrapolated from Figure 3) (11). We also found a positive autocorrelation coefficient at a lag of one breath for VT in all subjects during hyperoxic hypercapnia with an average autocorrelation coefficient of 0.34---comparable to that of Khatib and colleagues (11).

Hyperoxic hypercapnia increased the strength of the relationship of VI in immediately adjoining breaths, and it also increased the number of breath lags having significant autocorrelation coefficients for VI and VT. Of interest, this conditioning influence was confined to the volume components, which is consistent with the predominant effect of hypercapnia on the mean values (Figure 2). The mechanism underlying the dependence of breath components on the characteristics of the preceding respiratory cycle is not known. We (8) and others (26) have speculated that it might be related to after-discharge (facilitatory memory), that is, the augmentation of ventilatory output that persists after cessation of a primary stimulus (9). Indirect evidence suggests that the level of preceding hyperpnea contributes to the manifestation and magnitude of ventilatory afterdischarge (1, 27). Thus, the considerable hyperpnea in our subjects---a greater than threefold increase in VI---might have enhanced afterdischarge and so contributed to the increase in correlated variational activity. Afterdischarge is extremely sensitive to the prevailing level of PCO2, and even a slight decrease in PCO2, accompanying the increase in VI, greatly attenuates the magnitude and duration of afterdischarge (28). Conversely, Engwall and colleagues (29) showed in awake goats that quite mild hypercapnia (~ 2 mm Hg increase in arterial PCO2) substantially enhanced the magnitude of afterdischarge and prolonged its duration 2.4-fold. The investigators attributed the enhancement of afterdischarge to central CO2 chemical drive rather than to the level of ventilation.

Another contributor to the stronger relationship between volume components of neighboring breaths during hyperoxic hypercapnia may be the blunting of the peripheral chemoreceptors (17). The fast response time of the peripheral chemoreceptors makes them especially important in achieving rapid adjustments around the mean level of ventilation from one breath to the next. Blunting their activity with hyperoxia could have reduced this fine-adjustment mechanism, causing deviations to return more slowly to the mean level, with consequent increase in correlated behavior.

Employing the same mathematical techniques, we recently examined the effect of increasing inspiratory elastic loads on variational activity of breathing in healthy subjects (30). Loading decreased the total variational activity of VT and TE without altering the correlated behavior. The differing responses suggest that the alteration in variational activity of breathing is specific to the nature of the stimulus, probably reflecting behavioral (cortical) influences on the respiratory controller during elastic loading, and automatic (subcortical) influences during hypercapnia.

Oscillatory Behavior

Hyperoxic hypercapnia increased the number of subjects displaying significant ventilatory oscillations and altered their cycle time (Table 1). The development of ventilatory oscillations is thought to result from unstable dynamics in chemical feedback of the respiratory control system (2). In a mathematical model, Khoo and colleagues (12) demonstrated that an increase in loop gain of the respiratory control system resulted in oscillatory behavior. Hypercapnia is considered to decrease loop gain and it eliminates the periodic breathing induced by hypoxia (12, 13). However, the action of CO2 as a stabilizing influence has generally been demonstrated in a system rendered unstable through use of hypoxia, whereas we compared hyperoxic hypercapnia with a stable state, viz., normal resting breathing. Another factor that can lead to ventilatory oscillations is increased circulation time; however, this is an unlikely source of the oscillations with hypercapnia since the latter increases cardiac output (12). Decreased dampening of the respiratory system secondary to a decrease in lung volume can also cause oscillations (9), but this is a less likely mechanism since hypercapnia does not consistently decrease functional residual capacity.

The most striking action of hyperoxic hypercapnia on the oscillatory behavior of breathing was the threefold reduction in the frequency of significant oscillations in VI and VT. For VI, the frequency of oscillations decreased from 0.12 cycle/ breath during air breathing to 0.04 cycle/breath during hyperoxic hypercapnia; this signifies an increase in the number of breaths per cycle from eight during air breathing to 25 during hyperoxic hypercapnia---alternatively, this can be considered in terms of an increase in cycle time from ~ 28 to 88 s (assuming an average breath duration of 3.5 s) (7). This prolongation of cycle time during hyperoxic hypercapnia is compatible with the slower response of the central chemoreceptors compared with the peripheral chemoreceptors (14). Swanson and Bellville (31) quantitated the response time to step changes in PETCO2 while keeping end-tidal O2 concentration constant. The average time constant of the slow component, considered to reflect central chemoreceptor activity, was 75 s. This response time is compatible with the cycle length of the low-frequency oscillations during hyperoxic hypercapnia in our subjects, suggesting that activation of the central chemoreceptors plays an important role in their development. Employing a system to buffer breath-to-breath variations in PETCO2, Modarreszadeh and Bruce (32) achieved a marked reduction in spontaneous variability in VI in healthy subjects. Although the experimental conditions in their study were necessarily constrained compared with the free breathing permitted in our subjects, the results of the two studies complement each other. Modarreszadeh and Bruce (32) found that a decrease in spontaneous fluctuations of PCO2 had its greatest effect on the low-frequency band of the power spectrum for VI, and we found that an increase in mean PCO2 also had its predominant influence on the low-frequency band of the spectra for VI and VT (Figures 6 and 8 and Table 1).

Activation of the peripheral chemoreceptors are considered to play an important role in the development of ventilatory oscillations (2, 12, 13). In awake ponies, however, Brown and colleagues (33) observed ventilatory oscillations when the carotid bodies were denervated, indicating that the peripheral chemoreceptors are not essential for the development of oscillatory behavior. In elderly subjects displaying ventilatory oscillations during sleep, Pack and colleagues (34) observed virtual synchrony of electroencephalographic and VI time series on cross-correlation analysis. This result led them to conclude that the ventilatory oscillations were due to changes occurring in the state-dependent input to the respiratory control system at sleep onset, rather than resulting from instability in the peripheral chemoreceptors. In addition to chemical feedback loops, nonchemical factors almost certainly contribute to the variational activity in breathing (32). Multiple central pattern generators are believed to exist in the neuronal networks of the brainstem (35), which generate a coordinated output to different groups of respiratory motor neurons resulting in multiple breathing patterns (33). It is conceivable that hyperoxic hypercapnia may have a direct effect, or an indirect effect via release of neurotransmitters, on the central pattern generators resulting in altered oscillatory behavior of the type observed in our subjects.

Ventilatory Response to Hypercapnia (CO2 Gain)

Virtually none of the fractions of the variational activity in breathing was related to CO2 gain (Tables 3 and 4). This may be partly due to limitations of the technique for measuring CO2 gain since the ventilatory response to CO2 has limited reproducibility. We sought to minimize this factor by selecting responses that showed a high correlation between VI and PCO2, and taking the average of three measurements. Another reason for the lack of a relationship between the fractions of variational activity and CO2 gain might be the relatively small range of measured CO2 gains; however, the observed four-fold increase in CO2 gain makes this a less likely possibility. In addition to this technical consideration, the lack of a significant relationship between the variability fractions and CO2 gain suggest that primary neuronal processes rather than direct chemoreceptor activity might be responsible for the observed changes in variational activity of breathing.

The lack of a significant relationship between variational activity in breathing and CO2 gain might also be related to the fact that the measurements were obtained on separate days. To overcome this problem, we examined the relationship between the various components of variational activity and a measure of quasi-steady-state CO2 gain, based on the ratio of change in VI between breathing air and hyperoxic hypercapnia to the change in PETCO2 between air and hyperoxic hypercapnia. A negative correlation was observed between the autocorrelation coefficient at a lag of one breath for VI during hyperoxic hypercapnia and this measurement of CO2 gain (r = -0.52, p = 0.05). That is, subjects with a high CO2 gain exhibited a weak correlation between values of VI in immediately adjoining breaths during hyperoxic hypercapnia. This negative correlation suggests that the higher the CO2 gain, the more random or "noisy" the breath-to-breath variability in VI during hyperoxic hypercapnia.

In summary, hyperoxic hypercapnia produced an increase in the gross variational activity of VI and VT and a reduction in that of TI and TE. Greater than 80% of the variational activity for each breath component consisted of uncorrelated random fluctuations (white noise), correlated behavior represented 9 to 20%, and oscillatory behavior represented much less than 1% of the variational activity during both air and hyperoxic hypercapnia. Hyperoxic hypercapnia increased the relationship of VI in immediately adjoining breaths, and it increased the number of consecutive significant autocorrelation coefficients (i.e., short-term memory) of VI and VT. Hyperoxic hypercapnia induced low-frequency oscillations in VI and VT. Virtually none of the fractions of variational activity in breathing correlated with CO2 gain, although the autocorrelation coefficient at a lag of one breath for VI during hyperoxic hypercapnia was negatively correlated with a quasi-steady-state measure of CO2 gain. In conclusion, although the increase in variability of volume components with hyperoxic hypercapnia was predominantly due to uncorrelated random fluctuations, a significant portion resulted from increased correlated behavior accompanied by the development of low-frequency oscillations with a cycle time consistent with central chemoreceptor activation.

    Footnotes

Correspondence and requests for reprints should be addressed to Amal Jubran, M.D., division of Pulmonary and Critical Care Medicine, Edward Hines Jr. Veterans Administration Hospital, Route 111N, Hines, IL 60141.

(Received in original form January 24, 1997 and in revised form June 4, 1997).

Acknowledgments: The writers gratefully thank Rebecca Young, Malinda Mazur, Christine Mullner, and Wilbert Armstrong for their technical assistance, and Eugene N. Bruce, Ph.D., for repeated advice regarding data analysis.

Supported by a Merit Review grant from the Veterans Administration Research Service.

    References
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

1. Badr, M. S., J. Skatrud, and D. A. Dempsey. 1992. Determinants of post-stimulus potentiation in humans during NREM sleep. J. Appl. Physiol. 73: 1958-1971 [Abstract/Free Full Text].

2. Cherniack, N. S., and G. S. Longobardo. 1986. Abnormalities in respiratory rhythm. In N. S. Cherniack and J. G. Widdicombe, editors. Handbook of Physiology, Vol. II: The Respiratory System, Control of Breathing. American Physiological Society, Bethesda, MD. 729-749.

3. Mador, M. J., and M. J. Tobin. 1991. Effect of alterations in mental activity on the breathing pattern in healthy subjects. Am. Rev. Respir. Dis. 144: 481-487 [Medline].

4. Modarreszadeh, M., E. N. Bruce, and B. Gothe. 1990. Non random variability in respiratory cycle parameters of humans during stage 2 sleep. J. Appl. Physiol. 69: 630-639 [Abstract/Free Full Text].

5. Newsom-Davis, J., and D. Stagg. 1975. Interrelationships of the volume and time components of individual breaths in resting man. J. Physiol. (Lond.) 245: 481-498 [Abstract/Free Full Text].

6. Priban, I. P.. 1963. An analysis of some short-term patterns of breathing in man at rest. J. Physiol. (Lond.) 166: 425-434 .

7. Tobin, M. J., M. J. Mador, S. M. Guenther, R. F. Lodato, and M. A. Sackner. 1988. Variability of resting respiratory center drive and timing in healthy subjects. J. Appl. Physiol. 65: 309-317 [Abstract/Free Full Text].

8. Tobin, M. J., K. L. Yang, A. Jubran, and R. F. Lodato. 1995. Interrelationship of breath components in neighboring breaths of normal eupneic subjects. Am. J. Respir. Crit. Care Med. 152: 1967-1976 [Abstract].

9. Eldridge, F. L., and D. E. Millhorn. 1986. Oscillation, gating, and memory in the respiratory control system. In N. S. Cherniack and J. G. Widdicombe, editors. Handbook of Physiology, Vol. II: The Respiratory System, Control of Breathing. American Physiological Society, Bethesda, MD. 93-114.

10. Killian, K. J., C. K. Mahutte, and E. J. M. Campbell. 1980. Resistive load detection during passive ventilation. Clin. Sci. 59: 493-495 [Medline].

11. Khatib, M. F., Y. Oku, and E. N. Bruce. 1991. Contribution of chemical feedback loops to breath-to-breath variability of tidal volume. Respir. Physiol. 83: 115-128 [Medline].

12. Khoo, M. C. K., R. E. Kronauer, K. P. Strohl, and A. S. Slutsky. 1982. Factors inducing periodic breathing in humans: a general model. J. Appl. Physiol. 53: 644-659 [Abstract/Free Full Text].

13. Berssenbrugge, A., J. A. Dempsey, C. Iber, J. Skatrud, and P. Wilson. 1983. Mechanisms of hypoxia-induced periodic breathing during sleep in humans. J. Physiol. (Lond.) 343: 507-534 [Abstract/Free Full Text].

14. Waggener, T. B., I. D. Frantz, A. R. Stark, and R. E. Kronauer. 1982. Oscillatory breathing patterns leading to apneic spells in infants. J. Appl. Physiol. 52: 1288-1295 [Free Full Text].

15. Chapman, K. R., E. N. Bruce, B. Gothe, and N. S. Cherniack. 1988. Possible mechanisms of periodic breathing during sleep in humans. J. Appl. Physiol. 64: 1000-1008 [Abstract/Free Full Text].

16. Read, D. J. C.. 1967. A clinical method for assessing the ventilatory response to carbon dioxide. Australas. Ann. Med. 16: 20-32 [Medline].

17. Gardner, W. N.. 1980. The pattern of breathing following step changes of alveolar partial pressures of carbon dioxide and oxygen in man. J. Physiol. (Lond.) 300: 55-73 [Abstract/Free Full Text].

18. Perez, W., and M. J. Tobin. 1985. Separation of factors responsible for change in breathing pattern induced by instrumentation. J. Appl. Physiol. 59: 1515-1520 [Abstract/Free Full Text].

19. Chatfield, C. 1989. The Analysis of Time Series: An Introduction. Chapman and Hall, London.

20. Pack, A. I., D. A. Silage, R. P. Millman, H. Knight, E. T. Shore, and D. C. C. Chung. 1988. Spectral analysis of ventilation in elderly subjects awake and asleep. J. Appl. Physiol. 64: 1257-1267 [Abstract/Free Full Text].

21. Cunningham, D. J. C., and W. N. Gardner. 1977. A quantitative description of breathing during steady-state CO2 inhalation in man, with special emphasis on expiration. J. Physiol. 272: 613-632 [Abstract/Free Full Text].

22. Miller, J. P., D. J. C. Cunningham, B. B. Lloyd, and J. M. Young. 1974. The transient respiratory effects in man of sudden changes in alveolar CO2 in hypoxia and in high oxygen. Respir. Physiol. 20: 17-31 [Medline].

23. Ledlie, J. F., S. G. Kelsen, N. S. Cherniack, and A. P. Fishman. 1981. Effects of hypercapnia and hypoxia on phrenic nerve activity and respiratory timing. J. Appl. Physiol. 51: 732-738 [Abstract/Free Full Text].

24. Gardner, W. N.. 1977. The relation between tidal volume and inspiratory and expiratory times during steady-state carbon dioxide inhalation in man. J. Physiol. (Lond.) 272: 591-611 [Abstract/Free Full Text].

25. Bolton, D. P. G., and J. Marsh. 1984. Analysis and interpretation of turning points and run lengths in breath-by-breath ventilatory variables. J. Physiol. (Lond.) 351: 451-459 [Abstract/Free Full Text].

26. Benchetrit, G., and F. Bertrand. 1975. A short-term memory in the respiratory centers: statistical analysis. Respir. Physiol. 23: 147-158 [Medline].

27. Ahmed, M., G. G. Giesbrecht, C. Serrette, D. Georgopoulos, and N. R. Anthonisen. 1993. Respiratory short-term potentiation (after-discharge) in elderly humans. Respir. Physiol. 93: 165-173 [Medline].

28. Engwall, M. J. A., L. Daristotle, W. Z. Niu, J. A. Dempsey, and G. E. Bisgard. 1991. Ventilatory afterdischarge in the awake goat. J. Appl. Physiol. 71: 1511-1517 [Abstract/Free Full Text].

29. Engwall, M. J. A., C. A. Smith, J. A. Dempsey, and G. E. Bisgard. 1994. Ventilatory afterdischarge and central respiratory drive interactions in the awake goat. J. Appl. Physiol. 76: 416-423 [Abstract/Free Full Text].

30. Brack, T., A. Jubran, and M. J. Tobin. 1997. Effect of elastic loading on variational activity of breathing. Am. J. Respir. Crit. Care Med. 155: 1341-1348 [Abstract].

31. Swanson, G. D., and J. W. Bellville. 1975. Step changes in end-tidal CO2: methods and implications. J. Appl. Physiol. 39: 377-385 [Abstract/Free Full Text].

32. Modarreszadeh, M., and E. N. Bruce. 1994. Ventilatory variability induced by spontaneous variations of PaCO2in humans. J. Appl. Physiol. 76: 2765-2775 [Abstract/Free Full Text].

33. Brown, D. R., H. V. Forster, A. S. Greene, and T. F. Lowry. 1993. Breathing periodicity in intact and carotid body denervated ponies during normoxia and chronic hypoxia. J. Appl. Physiol. 74: 1073-1082 [Abstract/Free Full Text].

34. Pack, A. I., M. F. Cola, A. Goldszmidt, M. D. Ogilvie, and A. Gottschalk. 1992. Correlation between oscillations in ventilation and frequency content of the electroencephalogram. J. Appl. Physiol. 72: 985-992 [Abstract/Free Full Text].

35. Feldman, J. L. 1986. Neurophysiology of breathing in mammals. In F. E. Bloom, editor. Handbook of Physiology, Vol. I: The Nervous System, Intrinsic Regulatory Systems of the Brain. American Physiological Society, Bethesda, MD. 463-524.
    APPENDIX

Mathematical Methods

Stationary. Time-series analysis such as autocorrelation requires stationary data, i.e., the statistical properties of the breath series should be time-invariant. To ensure data stationarity, 700 consecutive breaths during air and hyperoxic hypercapnia in which the values of VI, VT, TI, and TE did not display large deviations from the mean on visual inspection were chosen for data analysis. Because a small number of outliers can substantially affect a correlogram (19), we employed the clamping technique described by Pack and colleagues (20). Data points that were more than 2 SD above or below the mean value were identified and replaced by values that were equal to the mean plus 2 SD or the mean minus 2 SD, respectively. Next, polynomial equations of increasing order were employed to detrend the chosen data strings. The polynomial equation, with the breath number as the independent variable and the corresponding breath component as dependent variable, that best fitted the clamped data was found by multiple regression (8, 19). This entails calculation of the residual sum of squares as each successive term was added, then calculating the F-value as follows:
F=[(RSS<SUB>n</SUB>−RSS<SUB>n+1</SUB>)/(P<SUB>n+1</SUB>−P<SUB>n</SUB>)]/[RSS<SUB>n+1</SUB>/(N−P<SUB>n+1</SUB>)]

where RSS is the residual sum of squares, subscript n refers to the lower order polynomial, subscript n + 1 is the next higher order polynomial, P is the number of adjustable parameters (order of polynomial plus 1 for intercept), and N is the number of points in the data string. The F-value was calculated for successive orders of polynomials, and the highest order with a significant F-value signified the polynomial that yielded the best fit to the data string. The order of the polynomials were similar for air and hyperoxic hypercapnia for all breath components. The data were detrended by subtracting the polynomial fit from the breath series on a breath-by-breath basis. Autocorrelation analysis was then performed on the detrended data.

Autocorrelation analysis. Autocorrelation analysis was performed to calculate the fraction of variational activity that is correlated among serial breaths, as opposed to random, uncorrelated fluctuations (white noise) in breathing (8, 19). Autocorrelation coefficients at a lag of one breath for VI, VT, TI, and TE were calculated for each data string. The relative strength of short-term memory of each breath component was computed as the number of consecutive lags, starting at a lag of one breath, that displayed autocorrelation coefficients that were statistically different from zero at p < 0.01. Thus, memory extended to the point at which the autocorrelation coefficient first lost statistical significance at the p < 0.01 level.

We rigorously investigated the adequacy of this detrending procedure in achieving stationarity. The data strings of 640 breaths during air and hyperoxic hypercapnia were each divided into two segments: an "early" segment consisting of breaths 0 to 320, and a "late" segment consisting of breaths 321 to 640. Then, polynomial equations of increasing order were fitted to each segment of the clamped data. The data were detrended by subtracting the polynomial fit from the breath series on a breath-by-breath basis. Next, the autocorrelation coefficients and the number of lags with significant serial correlations during air breathing for each breath component were calculated for the early and late segments; a similar analysis was performed on the data obtained during hyperoxic hypercapnia. the autocorrelation coefficients at a lag of one breath for VI, VT, TI, and TE were equivalent in the early and late segments during both air breathing (p > 0.49 in all instances) and hyperoxic hypercapnia (p > 0.22 in all instances). Likewise, the number of breath lags displaying significant correlations for VI, VT, TI, and TE were not different between the early and late segments during either air breathing (p > 0.21 in all instances) or hyperoxic hypercapnia (p > 0.20 in all instances).

The distribution of sample values of a correlation coefficient (r) is restricted to r values between -1.0 and +1.0, and is, in general, asymmetrical (i.e., not Gaussian); thus, a Fisher's z transformation was performed on each r value. After each r value was transformed to a function z, standard parametric statistics were performed. Specifically, the autocorrelation coefficients and number of lags with significant serial correlations for breath components during air breathing were compared with those during hyperoxic hypercapnia using Fisher's z transformation with paired t tests.

Spectral analysis. Power spectra of the data strings during air and hyperoxic hypercapnia were obtained using a fast Fourier transform (FFT) algorithm to detect significant oscillations and to quantify their magnitude. After removing any linear trend, the data were multiplied by a cosine-tapered rectangular (Hanning) window, and a FFT was applied to obtain the power spectrum. Spectral analysis was performed on 640 consecutive breaths, subdivided into 10 segments of 64 breaths each; final spectral estimates were calculated by averaging the spectra of the 10 subsegments (4). Each power spectrum was first characterized by calculation of the centroid frequency, which is the mathematically weighted median frequency of the spectrum. To allow for the possibility of nonuniform baseline power, we divided the spectrum into a low-frequency (0.0 to 0.2 cycle/breath) and high-frequency band (0.2 to 0.5 cycle/ breath); otherwise, the number of significant oscillations in the low-frequency band are likely to be overestimated, whereas the number of significant oscillations in the high-frequency band are likely to be underestimated (4). Significant peaks were detected in each of the frequency bands by comparing the power at each frequency with the mean power in that band. Peaks were considered significant when at least two consecutive data points lay above the 95% confidence interval of each spectral estimate (15). The presence of a significant peak in the power spectrum indicates that part of the variance in the data is due to an oscillation with a period equal to the inverse of the frequency of the peak. The area inscribed by the peak (amount of power) reflects the degree of variability in the signal resulting from oscillations at that frequency. Differences in centroid frequencies, number of significant peaks, frequencies of significant peaks, and the power of significant peaks were compared during air and hyperoxic hypercapnia by paired t tests.

Fractionation of variational activity. The total variational activity of the data was modeled as a mixture of autoregressive variations using a first-order autoregressive model (AR1), oscilatory variations, and uncorrelated white noise [w(n)] using the technique described by Modarreszadeh and colleagues (4):
x(n)=a<SUB>1</SUB>x(n−1)+ 
<LIM><OP>∑</OP><LL>j=1</LL><UL>ρ</UL></LIM>[b<SUB>j1</SUB>sin (2πf<SUB>j</SUB>n)+b<SUB>j2</SUB>cos (2πf<SUB>j</SUB>n)]+w(n)

where x = ventilatory variable, n = breath number, a1 = autoregressive coefficient, rho  = the maximal number of significant peaks, fj = frequency of significant peaks, and bj1 and bj2 = amplitudes of the sinusoids. The expression within the brackets is the contribution of each discrete oscillation, having a frequency of fj, to the variance of the data. The coefficients a1, bj1, and bj2 were calculated by the least-squares method. The variance secondary to autoregressive components and w(n), i.e., sigma 2AR, was calculated using the standard error of the estimate (SEE): (SEE)2/1 - a12; the variance caused by w(n) was calculated as (SEE)2, whereas the variance caused by the autoregressive fraction was calculated as sigma 2AR - (SEE)2. If the a1 coefficient of the data string for a breath component in a subject was statistically different from zero, a second-order autoregressive model (AR2) was then fitted to the data and the entire fractionating process was repeated (19). For the AR2 model, the variance secondary to the autoregressive fraction and w(n), i.e., sigma 2AR2, was calculated as:
σ<SUP>2</SUP><SUB>AR2</SUB>=(SEE)<SUP>2</SUP>(1−a<SUB>2</SUB>)/(1+a<SUB>2</SUB>)<SUP>*</SUP>(1+a<SUB>1</SUB>−a<SUB>2</SUB>)<SUP>*</SUP>(1−a<SUB>1</SUB>−a<SUB>2</SUB>)

The variance secondary to the oscillatory fraction (sigma 2osc) was calculated as: 0.5[bj12 + bj22]. Some caution is needed in interpreting the separation of the nonrandom fractions into correlated and oscillatory behavior. The autoregressive coefficients are the characteristic features of the time domain transformed into the z domain. The Fourier coefficients are the characteristic features of the time domain transformed into the frequency domain. The z domain and frequency domain are closely related to each other. Although the partitioning between the two fractions has implications for inferring different mechanisms, it could result from the differences in statistical efficiency of the two fractions. Analysis of variance (ANOVA) was used to compare the fractions of the total variability caused by correlated, oscillatory, and uncorrelated random behavior in each subject during air and hyperoxic hypercapnia. If ANOVA was statistically significant, then a Newman-Keuls test, with appropriate correction for multiple comparisons, was performed to determine which of the fractions of the variational activity of breathing was different during air and hyperoxic hypercapnia.





This article has been cited by other articles:


Home page
Am. J. Physiol. Regul. Integr. Comp. Physiol.Home page
M.-N. Fiamma, C. Straus, S. Thibault, M. Wysocki, P. Baconnier, and T. Similowski
Effects of hypercapnia and hypocapnia on ventilatory variability and the chaotic dynamics of ventilatory flow in humans
Am J Physiol Regulatory Integrative Comp Physiol, May 1, 2007; 292(5): R1985 - R1993.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
T. Brack, A. Jubran, F. Laghi, and M. J. Tobin
Fluctuations in End-Expiratory Lung Volume during Cheyne-Stokes Respiration
Am. J. Respir. Crit. Care Med., June 15, 2005; 171(12): 1408 - 1413.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
T. Brack, A. Jubran, and M. J. Tobin
Dyspnea and Decreased Variability of Breathing in Patients with Restrictive Lung Disease
Am. J. Respir. Crit. Care Med., May 1, 2002; 165(9): 1260 - 1264.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
E. N. Bruce
Assessing Respiratory Control during Spontaneous Breathing . Practice May Be More Difficult than Theory
Am. J. Respir. Crit. Care Med., April 15, 2002; 165(8): 1033 - 1034.
[Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
J. G. Van den Aardweg and J. M. Karemaker
Influence of Chemoreflexes on Respiratory Variability in Healthy Subjects
Am. J. Respir. Crit. Care Med., April 15, 2002; 165(8): 1041 - 1047.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
H. L. PREAS II, A. JUBRAN, R. W. VANDIVIER, D. REDA, P. J. GODIN, S. M. BANKS, M. J. TOBIN, and A. F. SUFFREDINI
Effect of Endotoxin on Ventilation and Breath Variability . Role of Cyclooxygenase Pathway
Am. J. Respir. Crit. Care Med., August 15, 2001; 164(4): 620 - 626.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
A. JUBRAN and M. J. TOBIN
Effect of Isocapnic Hypoxia on Variational Activity of Breathing
Am. J. Respir. Crit. Care Med., October 1, 2000; 162(4): 1202 - 1209.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
E. N. Bruce
Invited Editorial on "Irregularities and power law distributions in the breathing pattern in preterm and term infants"
J Appl Physiol, September 1, 1998; 85(3): 787 - 788.
[Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
T. BRACK, A. JUBRAN, and M. J. TOBIN
Effect of Resistive Loading on Variational Activity of Breathing
Am. J. Respir. Crit. Care Med., June 1, 1998; 157(6): 1756 - 1763.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by JUBRAN, A.
Right arrow Articles by TOBIN, M. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by JUBRAN, A.
Right arrow Articles by TOBIN, M. J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Proc. Am. Thorac. Soc. Am. J. Respir. Cell Mol. Biol.
Copyright © 1997 American Thoracic Society