Journal of Applied Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 83: 2005-2011, 1997;
8750-7587/97 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
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
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Trinder, J.
Right arrow Articles by Kay, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Trinder, J.
Right arrow Articles by Kay, A.

Vol. 83, Issue 6, 2005-2011, December 1997

Correlation between ventilation and EEG-defined arousal during sleep onset in young subjects

John Trinder, John A. Van Beveren, Philip Smith, Jan Kleiman, and Amanda Kay

School of Behavioural Science, University of Melbourne, Parkville, Victoria 3052, Australia

ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
FOOTNOTES
REFERENCES


ABSTRACT

Trinder, John, John A. Van Beveren, Philip Smith, Jan Kleiman, and Amanda Kay. Correlation between ventilation and EEG-defined arousal during sleep onset in young subjects. J. Appl. Physiol. 83(6): 2005-2011, 1997.---In studies of elderly individuals, ventilation and EEG-defined arousal have been shown to vary periodically and synchronously. Such results have been interpreted as indicating the primacy of sleep/wake state in causing ventilatory instability during sleep onset. However, because the elderly individuals studied were periodic breathers, the results do not unequivocally support this conclusion. In this study the relationship between ventilation and EEG-defined arousal was assessed in a group of 21 young, healthy men in whom ventilatory instability during sleep onset was not periodic. Ventilation and EEG (O1-A2) recordings were collected, and the longest uncontaminated periods from early and late in sleep onset were selected for subsequent analysis. The 84 time series (21 subjects, 2 variables, and 2 occasions in sleep onset) were subjected to spectral analysis to identify periodicity, and the relationship between the two variables was determined by cross-correlational methods. The results indicated that the time series were nonperiodic, yet significant correlations were observed between the two variables. The data support the view that during sleep onset ventilatory instability is driven primarily by variations in sleep/wake arousal level.

electroencephalogram; periodic breathing; respiratory instability; sleep onset


INTRODUCTION

VENTILATION IS UNSTABLE in most individuals during the process of getting to sleep [stage 1/2 non-rapid-eye-movement (NREM) sleep] (3, 14). Under particular conditions such as hypoxia (15, 21), or in particular groups of individuals such as elderly subjects (13, 16, 20), the variations in ventilation are periodic with a cycle time of ~40-60 s over individuals (the term "periodic" will be used to refer to variations in a variable such that the cycle of hyper- and hypoactivity has a constant duration within a subject). In other individuals, for example young and healthy individuals, the variation is not periodic (3, 7). Early hypotheses as to the mechanism underlying instability in the respiratory system during sleep onset focused on the observed periodicity and attributed the instability to the operation of feedback mechanisms in which feedback delay and overall loop gain are the critical concepts (4, 11).

Other explanations for periodic breathing have emphasized the effect of sleep/wake state on ventilation (3, 10, 14, 17). The effect is typically modeled as an abrupt change in the gain of the system in association with changes in state, with higher gain in wakefulness than sleep. The higher gain has been conceptualized as a tonic input to central respiratory drive that is specific to the awake state. It is frequently referred to as the "wakefulness stimulus." The state change, as identified by the transition between dominant alpha (wakefulness) and theta (light sleep) activity in the electroencephalogram (EEG), is associated with abrupt alterations in ventilation. Indeed, changes in state as brief as single breaths are associated with changes in ventilation (3, 5, 18).

In an attempt to demonstrate the critical role of state changes in the occurrence of respiratory instability during sleep onset, Pack et al. (12) showed that, during periodic breathing in elderly individuals, changes in ventilation were correlated with zero lag with changes in arousal state. Hyperventilation was shown to occur in association with relative high-frequency EEG activity, indicative of arousal, whereas hypoventilation or apnea was shown to occur in association with relative low-frequency EEG activity, indicative of light sleep. This relationship was interpreted as being consistent with the view that changes in state play a critical role in producing respiratory instability during sleep onset.

However, studies of the relationship between sleep/wake state and ventilation during periodic breathing may not fully elucidate the role of state changes. This is because the existence of periodicity in the range of ~40-60 s indicates a contribution from respiratory control mechanisms as there is not an intrinsic periodicity in sleep/wake oscillations during sleep onset. Thus, because the periodicity of changes in ventilation in these elderly individuals was due to respiratory mechanisms, it seems hazardous to conclude that the actual changes in ventilation were due to changes in state.

A more definitive test of the hypothesis that changes in state contribute to respiratory instability during sleep onset would be to investigate the relationship between ventilation and arousal state in individuals in whom ventilation was variable, but not periodic. Some evidence suggests respiratory instability during sleep onset in young adults is nonperiodic (3, 7). This is consistent with our informal observations over a number of studies of such subjects (5, 17, 18). The purposes of the present experiment were, first, to assess the periodicity of ventilation during sleep onset in young healthy men, and second, to assess the relationship between ventilation and EEG-defined arousal level in those subjects shown to have nonperiodic breathing.


METHODS

Subjects

The data from 21 young male subjects between the ages of 18 and 24 yr (average = 20.3 yr) are reported in this study. The subjects were nonsmokers with an average height of 179 cm (167-192 cm) and weight of 71.7 kg (60-84 kg) and had no known history of sleep or respiratory pathology.

The subjects were originally run in two studies by using identical laboratory procedures. Their data were selected for inclusion in this analysis if they contained two periods of at least 50 consecutive breaths in which there were no major movements, other disruptions of the recordings, or sustained changes in state. A further requirement was that one of the periods occur early in the sleep onset process, during what will be referred to as "phase 2" of sleep onset, and the other late in the sleep onset process, or phase 3 of sleep onset (8). The studies from which the subjects were selected were approved by the University of Melbourne's Human Ethics Committee, and subjects gave written informed consent before their participation in the studies.

General Laboratory Procedures

Subjects spent two nonconsecutive nights in the sleep laboratory of the School of Behavioural Science, the University of Melbourne. They were requested to refrain from consuming caffeine and alcohol on the day of the experimental sessions. Subjects arrived at ~2100, prepared for bed, and then the recording equipment was attached. Data collection commenced at ~2230. The equipment consisted of surface electrodes for the assessment of sleep/wake state, a full face mask to allow for the measurement of airflow and, in 17 subjects, a Millar catheter, which was introduced transnasally to the level of the epiglottis, to measure upper airway resistance. The variables of interest for the present analysis were sleep/wake state and ventilation, the latter being derived from the rate of airflow measurement.

Subjects were required to maintain a supine position during the experiment. The purpose was to record respiratory and sleep/wake state variables during sleep onset. Accordingly, subjects were instructed to go to sleep after lights out, but they were then awakened as soon as they had attained ~10-15 min of stage 2 NREM sleep. After subjects were awakened, the experimenter ensured that subjects were alert before instructing them to return to sleep. These procedures were repeated until data collection terminated (0300-0400). On average, four to five sleep onsets were obtained each night. The two data sets, one early and one late in a sleep onset, were selected from among these sleep onsets according to rules described below.

Sleep/Wake State-Recording Procedures

The EEG was recorded from central (C3-A2) and occipital (O1-A2) sites, the electrooculogram from two electrodes placed ~1 cm above and slightly lateral to the outer canthus of one eye and the same position but slightly below the other eye, and the electromyogram was recorded from the submental muscles. Grass gold cup electrodes were used, with the C3 and O1 electrodes being secured with collodion, and the remaining electrodes were attached with surgical tape.

Respiratory Measurements

Subjects wore either a modified Vital Signs "Downs Continuous Positive Airway Pressure (CPAP)" or a Commonwealth Industrial Gases (CIG) anesthetic face mask that covered the mouth and nose and was held tightly in place by a head strap. Route of breathing was not controlled. A heated Morgan pneumotachograph was attached to the mask. For subjects using the Downs CPAP mask, the pneumotachograph was attached directly to a port in the mask so that inspiratory and expiratory flow were obtained. For subjects using the CIG mask, the pneumotachograph was attached to the expiratory side of a two-way nonrebreathing valve so that expiratory flow was obtained. The total dead space of the mask and pneumotachograph varied as a function of facial configuration and was ~155 ml for the Downs CPAP mask and 120 ml for the CIG mask. For measurement of airflow the pneumotachograph was connected to a differential pressure transducer (Validyne DP45-14), the output of which was converted to a voltage signal by using a carrier demodulator (Validyne CD15). Airflow was calibrated by using a Shorate flowmeter (1355). For the subjects who originally participated in a study in which upper airway resistance was measured, a transducer-tipped catheter was inserted transnasally and positioned visually ~1 cm below the tongue base, 16-18 cm from the nares. The measurements from this device are not relevant to this study.

All sleep/wake and respiratory measurements were recorded by using a model 7D Grass polygraph. The occipital EEG and rate of airflow were also recorded and stored on an IBM-compatible 486 PC via a 16-bit analog-to-digital converter sampling at 100 Hz. Automated off-line analysis was used to compute expiratory tidal volume (VT) and respiratory-timing information.

Statistical Procedures

The steps in the statistical analyses are summarized in Table 1.

Table  1.   Steps in statistical analysis
Step 1: Initial preparation of data
  Determination of phase of sleep onset
  Identification of time series
  Conversion of analog EEG into a weighted ratio EEG
  Conversion of the irregularly sampled time series to a regularly spaced time series
  Test for stationarity and detrending of each time series
Step 2: Data analysis
  Estimation of variability of ventilation
  Spectral analysis of each series
  Correlations between tidal volume and weighted ratio

EEG, electroencephalogram.

Initial preparation of data. DETERMINATION OF PHASE OF SLEEP ONSET. Each sleep onset was divided into three phases (1, 2, and 3) according to previously published procedures (8). Briefly, phase 1 was defined as the initial period of wakefulness and was identified by sustained EEG alpha activity. The phase began with lights out and was terminated when alpha activity was first replaced by theta activity. This point was defined as that when the EEG associated with three of five consecutive breaths was classified as predominantly theta activity. Phase 2 was defined by alternating periods of alpha and theta EEG activity. This phase was terminated, and phase 3 commenced at the first occurrence of a sleep spindle or K complex. Phase 3 corresponded to late in sleep onset and was identified by periods of theta activity with sleep spindles and K complexes, which were repeatedly broken by arousals. Phase 3 ended when stable stage 2 sleep (absence of arousals) was achieved.

IDENTIFICATION OF TIME SERIES. For each subject, two data sets were identified, one each from phases 2 and 3. A data set consisted of two sequences, or time series, one of VT values and the other of the digitized EEG (O1-A2 recording). VT was selected as the ventilatory variable because it accords better with the assumptions of time series analyses and because previous work has indicated that changes in ventilation at transitions between states are a consequence of changes in VT (18). To identify the time series, the raw data from the two nights were inspected, and the longest segment of continuous artifact-free data in both phases 2 and 3 was selected. An artifact was defined as any disruption of the record that prevented the computation of VT or obscured the EEG tracing. The lengths of the identified time series for each phase and subject are shown in Table 2.

Table  2.   Time series length and periodicity in ventilation and EEG for phases 2 and 3 
Subject No. Phase 2
Phase 3
TS Length, s Ventilation
EEG
TS Length, s Ventilation
EEG
Fisher test Cycle time, s Fisher test Cycle time, s Fisher test Cycle time, s Fisher test Cycle time, s

 1 908 6.7 14.8* 22.6 1,192 12.7* 59.6 18.1* 29.8
 2 680 5.4 16.6* 324 10.3* 28.8 6.6
 3 272 3.4 7.4* 85.1 352 7.0 7.7
 4 1,364 24.1* 27.2 20.1* 960 6.4 8.0
 5 464 5.1 5.0 228 5.9 6.4
 6 1,728 6.7 8.5 920 7.0 24.2* 24.0
 7 772 4.8 14.9* 1,484 8.0 10.9* 74.0
 8 460 7.7 15.8* 40.8 288 3.7 4.2
 9 948 17.0* 47.2 19.5* 23.6 1,708 27.8* 42.6 10.7* 67.8
10 228 3.5 4.7 284 7.0 4.3
11 792 10.6* 31.7 21.2* 1,496 20.2* 7.5
12 744 15.8* 5.2 900 8.1 12.2*
13 1,592 11.4* 31.8 18.7* 1,416 15.3* 63.7 7.9
14 928 4.1 6.4 972 6.4 15.2* 29.0
15 932 10.3* 10.6* 578 7.5 9.5* 28.8
16 912 13.5* 47.9 17.7* 26.6 1,064 7.0 6.2
17 852 5.0 4.5 980 6.9 6.3
18 784 5.9 25.0* 588 5.6 6.7
19 788 6.3 6.5 752 5.5 5.2
20 1,084 8.8 29.8* 812 4.5 8.4
21 800 4.8 20.3* 904 6.1 15.6*

* Significant periodicity (periods between 20 and 90 s). Fisher test values with an asterisk, but without a corresponding cycle time, indicate significant periodicity outside interval 20-90s. Phase 2, sleep period with alternating alpha and theta activity; phase 3, later sleep period with alternating alpha and theta activity but with first appearance of sleep spindle and K complexes and repeated arousals.

CONVERSION OF THE ANALOG EEG INTO A DISCRETE WEIGHTED RATIO EEG. The two time series (VT and EEG) within each of the two data sets (phases 1 and 2) were sampled at different frequencies: VT at the respiratory rate (~0.25 Hz) and the EEG at the analog-to-digital-sampling frequency (100 Hz). Thus it was necessary to transform the EEG in two ways: first, to convert the EEG into a time series with the same frequency as VT, and, second, to represent the frequency and amplitude characteristics of the EEG by a single value.

To achieve the first of these purposes, the digitized EEG was broken into nonoverlapping segments corresponding to the breath intervals identified in the respiratory analysis. Thus a segment of EEG corresponded to each breath.

To achieve the second purpose, each EEG segment was subjected to a Fourier transform algorithm, and the resulting frequency f)-by-power (P) distribution over the frequency range of 3 to 13 Hz was transformed to a single value [the weighted ratio (R)]. This value reflected the extent to which power was concentrated in the high (alpha, 8-13 Hz) or low (theta, 3-7 Hz) end of the frequency distribution. The following procedures were used to calculate R. Two linear functions were constructed. The first of these (a) was such that at 3 Hz, it equaled zero, and at 13 Hz, it equaled one. The second function (b) equaled one at 3 Hz and zero at 13 Hz. The quantities <LIM><OP>∑</OP></LIM>faf )Pf ) and <LIM><OP>∑</OP></LIM>fbf)Pf) were then calculated. The first of these quantities expressed the power of higher frequency activity during a given breath, whereas the latter expressed the power of lower frequency activity during the same breath (over the frequencies 3-13 Hz). R was then defined as
R = <LIM><OP>∑</OP><LL><IT>f</IT></LL></LIM><IT> a</IT>(<IT>f</IT>)<IT>P</IT>(<IT>f</IT>)<FENCE> </FENCE><LIM><OP>∑</OP><LL><IT>f</IT></LL></LIM> <IT>a</IT>(<IT>f</IT>)<IT>P</IT>(<IT>f</IT>*) + <LIM><OP>∑</OP><LL><IT>f</IT></LL></LIM> <IT>b</IT>(<IT>f</IT>)<IT>P</IT>(<IT>f</IT>)
Thus R could vary from ~1, indicating predominantly higher frequency activity, to ~0, indicating predominantly lower frequency activity.

In summary, these procedures allowed the EEG to be represented as a time series that expressed changes in the frequency content of the EEG and that corresponded in time with the time series for VT.

CONVERSION OF THE IRREGULARLY SAMPLED TIME SERIES TO A REGULARLY SPACED TIME SERIES. The statistical procedures subsequently applied to the data required that the time series had constant intervals. Thus the irregularly sampled time series for VT and R were resampled at intervals of 4 s by using the method described by Waggener et al. (19). The term "time series" as used in the remainder of the study refers to these evenly spaced (4-s) values for VT and R.

TEST FOR STATIONARITY AND DETRENDING OF EACH TIME SERIES. An additional requirement of subsequent statistical procedures was that each time series be stationary. That is, its mean and variance do not change with time. Each series was tested for stationarity by testing the autocorrelation coefficients obtained from the series by using the formula r2 < 1/lag (6), where lag refers to the number of adjacent values within the time series for which the series itself is correlated. If a series was not stationary, linear trends were removed, and in all cases this was sufficient to achieve stationarity. Of the 84 time series (21 subjects × 2 phases × 2 variables), 75 were initially stationary.

Data analysis. ESTIMATION OF VARIABILITY OF VENTILATION. A coefficient of variation (CV) was determined for VT for each subject in wakefulness and phases 2 and 3 of sleep onset. For wakefulness (phase 1), all breaths during phase 1 of the onset with the longest phase 1 were entered into the computation. In phases 2 and 3 this was calculated over the series of breaths that made up the time series.

SPECTRAL ANALYSIS OF EACH SERIES. The purpose of this component of the analysis was to determine whether variations in VT and R were periodic. The range of acceptable periodicities was restricted to between 20 and 90 s because periodicities outside this range were unlikely to occur as a function of respiratory mechanisms. The power spectral density of each time series was estimated by using periodogram methods (1). The periods of VT and R were extracted from the relevant spectral estimate by testing for periodicities by using the Fisher test (2). The frequency values of periodicity identified by the Fisher test were determined by examination of the periodogram.

CORRELATIONS BETWEEN VT AND R. The final stage of the analysis was to determine whether VT and R were related. The two series were cross-correlated for each subject in each phase of sleep onset by using three different procedures. First, the two time series were correlated. Second, the two series were correlated after the sequential dependencies in both series had been removed. This was done by fitting appropriate ARMA (autoregressive/moving average) models individually to both the VT and R time series in the manner previously described (2, 6). The ARMA model is a linear filter that removes the autocovariance structures within a series (2). Mathematically, the value of the series at any time is represented as the sum of two components, one predictable from the history of the series before that time and one that is not. The predictable part is represented as a linear combination of previous values of the series and random statistical fluctuations. The residual series obtained after the predictable structure has been estimated and removed is uncorrelated. In the second part of the analysis, the two residual series obtained in this way were cross-correlated.

The third part of the analysis investigated whether ventilation was predictable from the EEG, under the assumption that the VT and R series are, respectively, the input and output of a linear system. To do this, the Box-Jenkins "prewhitening" procedure was used (2, 6). Residuals were obtained for both series by fitting the R ARMA model to both the R and VT series. The latter series is the output that would be obtained from an unknown (black box) linear system if it were given an uncorrelated series as input. To determine whether such a representation was appropriate, the coherency spectrum was calculated between the two series. This function is a measure of the squared correlation between the frequency components of the two series. If a correlation exists, it is appropriate to attempt to estimate the transfer function of the linear system that relates the two series. In such cases, the R residuals were entered into the transfer function modeling program (1) as the input series, and the VT residuals, obtained by fitting the R ARMA model to VT, were entered as the output series. The predicted output series obtained from the transfer model were then cross-correlated with the actual VT series to assess the quality of the linear system representation.


RESULTS

The CVs averaged over subjects for wakefulness and early (phase 2) and late (phase 3) in sleep onset were 25.6 (SD = 15.2), 24.3 (SD = 14.8), and 31.2 (SD = 26.3), respectively. The values for each phase did not differ significantly from each other. These values are comparable with those reported for elderly subjects in two studies (7, 16) and for equivalently aged young subjects in one of these (16) but were greater than those reported for young subjects in one other study (7). Of note were the large differences between subjects in CV values, as indicated by the SDs.

The results of the spectral analysis of the time series are shown in Table 2. Of the 84 series, 34 showed significant periodicity. Fourteen of the 34 were outside the range of periods that might reflect periodicity generated by the respiratory system (20-90 s). Of the 20 time series that showed significant periodicity with periods between 20 and 90 s, 9 were from the ventilation series (5 from phase 2 and 4 from phase 3) and 11 were from the EEG series (5 from phase 2 and 6 from phase 3). However, only one of the 21 subjects (subject 9) had significant ventilatory periodicity with approximately the same period for both phases 2 and 3. This subject also had significant periodicity in the EEG time series but at different periods from those observed for ventilation and at different periods for phases 2 and 3. No subject had significant periodicity with the same period length in both ventilation and the EEG within a phase, although in three cases (S1, S9, and S16) both time series were significant such that the ventilatory period was approximately double the EEG period.

In summary, within subjects neither the presence of periodicity, nor the cycle time of any periodicity, was consistent between VT and R, or between phases 2 and 3. The unsystematic pattern of the observed significant periodicities suggests the test is particularly sensitive and likely to produce type 1 errors. Thus the data are consistent with a previous report (3) and with our informal observations, suggesting that in young healthy individuals instability in both ventilation and sleep/wake state during sleep onset is essentially nonperiodic.

The finding that neither ventilation nor changes in sleep/wake state were periodic, despite substantial variability in ventilation, provided the appropriate preconditions for assessing the hypothesis that variations in ventilation need not be periodic for changes in ventilation and arousal state to be synchronous.

The relationship between ventilation and the EEG measurement of arousal level is shown in Table 3. The data presented are the maximum correlation coefficients between series for lags between ±3 (VT phase delayed/advanced with respect to R by ~12 s). Values were considered only within this restricted range because Pack et al. (12) have indicated that ventilation and arousal level oscillate with an approximately zero lag, and also to reduce the likelihood of type 1 errors given the length of the time series, and thus the number of coefficients computed between any two series.

Table  3.   Cross-correlations between EEG measurements of arousal and ventilation
Subject No. Phase 2
Phase 3
A
B
C
A
B
C
Coefficient Lag Coefficient Lag Coefficient Lag Coefficient Lag Coefficient Lag Coefficient Lag

 1 0.38* 1 0.38* 1 0.35* 1 0.23* 1 0.16* 1 0.18* 1
 2 0.08 3 0.09  -1 0.08 0 0.18 0 0.42* 0 0.29* 0
 3 0.11 3 0.12* 3 0.03 1 0.06 3 0.18  -1 0.12 0
 4 0.17* 1 0.15* 1 0.19* 0 0.19* 1 0.16* 3 0.09 0
 5 0.55* 1 0.37* 1 0.38* 1 0.49* 1 0.25 1 0.41* 1
 6 0.09 2 0.09 0 0.04 2 0.31* 1 0.21* 1 0.21* 1
 7 0.14 1 0.11 1 0.15 1 0.40* 1 0.41* 1 0.40* 1
 8 0.16 2 0.17 0 0.15 0 0.14 0 0.13 3 0.16 0
 9 0.12 1 0.08 3 0.08 1 0.15* 2 0.05 2 0.16* 2
10 0.33* 1 0.36* 1 0.36* 1 0.36* 1 0.25* 1 0.27* 1
11 0.19*  -1 0.11 1 0.09 1 0.27* 1 0.24* 1 0.26* 1
12 0.19 1 0.17 1 0.23* 1 0.49* 1 0.37* 1 0.42* 1
13 0.58* 1 0.35* 1 0.33* 1 0.26* 2 0.26* 1 0.27* 1
14 0.15 1 0.22* 1 0.13 1 0.43* 1 0.34* 1 0.32* 1
15 0.11 2 0.12 1 0.10 1 0.11 2 0.08 1 0.15 2
16 0.12 1 0.15 1 0.17 1 0.45* 2 0.20* 2 0.20* 2
17 0.09 1 0.16 1 0.08 1 0.11 1 0.17 2 0.16 1
18 0.36* 2 0.16 1 0.16 1 0.41* 2 0.22* 2 0.22* 2
19 0.35* 1 0.23* 1 0.24* 1 0.25* 1 0.24* 1 0.25* 1
20 0.13 2 0.06 3 0.05 2 0.14 1 0.16 1 0.16 1
21 0.39* 1 0.22* 1 0.21* 1 0.49* 1 0.40* 1 0.40* 1

A, equal interval raw data; B, residuals after fitting autoregressive/moving average (ARMA) models; C, transfer-function analysis. Lag units, 4-s intervals; negative lags, weighted ratio lagged tidal volume. * Significant, P < 0.01.

Consider, first, columns A and B, for phases 2 and 3, in Table 3. The results indicate that the removal of sequential dependencies (column B) from the original series (column A) slightly reduced the magnitude of the relationship in most, but not all, subjects but did not affect the number of subjects for whom there were statistically significant associations between the two series. The results also indicate that the relationship became stronger as sleep onset developed (phase 3 compared with phase 2). In phase 2, coefficients for 9 subjects were significant, both for the raw series and after removal of sequential dependencies, whereas for phase 3 the subjects with significant coefficients numbered 15 and 14 before and after the removal of sequential dependencies, respectively. The extent to which the R series predicted the VT series (column C) showed a similar pattern, with significant transfer functions for 8 subjects in phase 2 and 15 in phase 3. The nature of the relationship between VT and R in one subject is illustrated in Fig. 1.


Fig. 1. Evenly spaced time series for weighted average EEG (R) and tidal volume (VT) for a single subject (subject 21) late in sleep onset with periods of theta activity with sleep spindles and K complexes, which were repeatedly broken by arousals (phase 3).
[View Larger Version of this Image (20K GIF file)]

The magnitude of the relationship between ventilation and EEG-defined arousal level was lower in these young subjects than in elderly subjects studied previously (12). Thus Pack et al. (12) reported that all of their eight elderly subjects had significant coefficients during stage 1/2 NREM sleep after sequential dependencies had been removed. The average value was 0.39. In contrast, in the present study 14 of 22 younger subjects had significant coefficients in phase 3 (approximately equivalent to stage 1/2 NREM), with an average value of 0.23. In addition, the phase relationship between the two series was slightly different between studies. The previous paper reported a zero phase lag between ventilation and EEG variables for six of eight subjects, whereas this study found that in most subjects ventilation was delayed with respect to the EEG by 4-8 s.

Finally, as noted above, there were marked individual differences in the variability of ventilation. The magnitude of the relationship between VT and R tended to be lower in subjects who exhibited relatively low ventilatory variability. Thus, of the five subjects in whom there was no relationship between the two time series during phase 3, three had very low ventilatory variability.


DISCUSSION

The present data are consistent with an earlier report (12) showing a synchronous and positive relationship between ventilation and level of arousal during periods of unstable ventilation characteristic of sleep onset. Furthermore, the data suggest that the level of arousal, as defined by EEG frequency in the 3- to 13-Hz range, predicts the level of ventilation. However, it should be noted that despite the general similarities between the results of this study and the previous study (12), there were three notable differences between the two.

The first was that although the younger subjects in this study exhibited respiratory variability during sleep onset, neither ventilation nor EEG activity was periodic. This contrasts with the clear periodicity in ventilation characteristic of the elderly subject sample studied previously (12). Second, the magnitude of the relationship between arousal level and ventilation was greater in the older subjects. Finally, the phase relationship between the two series was slightly different between the two studies. In the elderly subjects the phase of the two oscillations was essentially identical, whereas in the younger subjects changes in ventilation tended to be delayed relative to changes in the EEG by 4-8 s (~1-2 breaths).

The positive relationship between arousal level and ventilation provides strong support for the view that in normal young adults the instability in ventilation that characterizes sleep onset is determined primarily by instability in sleep/wake state. Three aspects of the data support this conclusion. First, the absence of periodicity in ventilatory fluctuations suggests that respiratory mechanisms do not make a major contribution to these ventilatory oscillations. This is because fluctuations in the output of the respiratory control system are, by nature of the system, periodic. Second, the transfer function, in which the R residuals were entered as the input function, was significant in 15 of 21 subjects during phase 3 of sleep onset. The interpretation of this finding is that the level of arousal, as reflected in EEG activity, determines the level of ventilation. Third, the change in ventilation lagged behind the change in arousal such that maximum ventilation typically occurred, on average, one breath after the initiation of an increase in EEG frequency.

It may be questioned whether the lag in ventilation was significant and whether it might have occurred as an artifact of data processing and reduction. It is concluded that the effect is significant because it occurred in the vast majority of subjects in both phases of sleep onset. Furthermore, careful inspection of the data-reduction methods did not identify a systematic bias on either variable. Indeed, this seems unlikely because the variables were converted to a common time base early in the analyses and were treated identically from that point onward. Furthermore, both variables were maintained as continuous, avoiding artifacts associated with conversion to categorical data. Although it is not possible to entirely rule out an artifact, the data suggest a lag in changes in ventilation compared with changes in state in these young subjects, an effect that is possibly related to the absence of periodicity in their ventilation.

In the young subjects in this study the relationship between an EEG measurement of arousal level and ventilation was weaker than that reported for older subjects (12). Although this may reflect any number of differences between the two studies, it is also possible that the difference was a consequence of a greater contribution of respiratory mechanisms to respiratory instability in the older subjects. Thus in individuals who are prone to respiratory instability, such as elderly subjects, periodic instability in the respiratory system and instability in level of arousal, may, through mutual entrainment, become synchronous. The combined effect of these two sources of instability in elderly subjects might then result in a stronger relationship between ventilation and level of arousal than that observed in young individuals in whom ventilatory instability is driven solely by instability in level of arousal. Furthermore, such a mechanism could also account for the difference between the two age groups in phase relationship between EEG and ventilation.

The absence of periodic breathing in these young adult men was not due to lower levels of respiratory variability in these compared with older subjects because the average CV in all phases was comparable to that reported previously for older subjects (7, 16). There was, however, some evidence to suggest that in the younger individuals, in whom there was no relationship between arousal level and ventilation, there was low ventilatory variability. Such an effect would be anticipated on purely statistical grounds.

In this study the relationship between fluctuations in sleep/wake state and ventilation became stronger as sleep onset developed. In addition to the fact that more subjects showed statistically significant relationships between EEG and ventilation during phase 3 than during phase 2, 17 of the 21 subjects had higher transform functions in phase 3 than in phase 2. The cause of this is uncertain. There was no consistent change in the periodicity from phase 2 to phase 3. There was, however, a nonsignificant increase in the variability of ventilation during phase 3, and our laboratory has reported previously that changes in ventilation at state transitions between alpha and theta EEG activity are larger during phase 3 than phase 2 (17). Furthermore, by definition, the difference in arousal level between the wake and sleep states is greater in phase 3 because the sleep state is "deeper." Thus the stronger statistical relationship between the two time series during phase 3 may have been due to greater variability in the variables during this phase.

The present study has not addressed the issue of the mechanism by which state changes might mediate changes in ventilation. It is widely accepted that state changes are associated with changes in central respiratory drive (14). However, it is likely that this effect is mediated through respiratory pump muscle activity and/or activity of the dilator muscles of the upper airway. Although it is known that airway resistance increases at state transitions from alpha to theta activity and decreases at transitions from theta to alpha (8, 9), this is not the only mechanism involved because there is only a weak relationship between the magnitude of the change in upper airway resistance and changes in ventilation (9). Furthermore, measurements of muscle electromyogram activity at these transitions indicate changes in both respiratory pump and upper airway muscles (22). Nevertheless, a significant relationship between changes in ventilation and airway resistance has been demonstrated during periodic breathing in healthy elderly subjects (7), suggesting that changes in ventilation in the present study are due, in part, to changes in resistance.

In conclusion, the data from this study support the view that fluctuations in arousal level play the major role in producing ventilatory instability during sleep onset in healthy young adults. In other groups, such as elderly subjects, this mechanism may be augmented by respiratory mechanisms.


FOOTNOTES

Address for reprint requests: J. Trinder, Dept. of Psychology, School of Behavioural Science, Univ. of Melbourne, Parkville, Victoria 3052, Australia.

Received 28 October 1996; accepted in final form 23 July 1997.


REFERENCES

1. Brockwell, P. J., and R. A. Davis. ITSM: An Interactive Time Series Modelling Package for the PC. New York: Springer-Verlag, 1991.
2. Brockwell, P. J., and R. A. Davis. Time Series: Theory and Methods. New York: Springer-Verlag, 1991.
3. Bulow, K. Respiration and wakefulness in man. Acta Physiol. Scand. 209, Suppl. 59: 5-110, 1963.
4. Carley, D. W., and D. C. Shannon. A minimal mathematical model of human periodic breathing. J. Appl. Physiol. 65: 1400-1409, 1988[Abstract/Free Full Text].
5. Colrain, I. M., J. Trinder, G. Fraser, and G. V. Wilson. Ventilation during sleep onset. J. Appl. Physiol. 63: 2067-2074, 1987[Abstract/Free Full Text].
6. Gottman, J. M. Time Series Analysis: A Comprehensive Introduction for Social Scientists. Cambridge, UK: Cambridge Univ. Press, 1992.
7. Hudgel, D. W., P. Devadatta, and H. Hamilton. Pattern of breathing and upper airway mechanics during wakefulness and sleep in healthy elderly humans. J. Appl. Physiol. 74: 2198-2204, 1993[Abstract/Free Full Text].
8. Kay, A., J. Trinder, G. Bowes, and Y. Kim. Changes in airway resistance during sleep onset. J. Appl. Physiol. 76: 1600-1607, 1994[Abstract/Free Full Text].
9. Kay, A., J. Trinder, and Y. Kim. Individual differences in the relationship between upper airway resistance and ventilation during sleep onset. J. Appl. Physiol. 79: 411-419, 1995[Abstract/Free Full Text].
10. Khoo, M. C. K., A. Gottschalk, and A. I. Pack. Sleep induced periodic breathing and apnea: a theoretical study. J. Appl. Physiol. 70: 2014-2024, 1991[Abstract/Free Full Text].
11. Longobardo, G. S., B. Gothe, M. D. Goldman, and N. S. Cherniack. Sleep apnea considered as a control system instability. Respir. Physiol. 50: 311-333, 1982[Medline].
12. Pack, A. I., M. F. Cola, A. Goldzmidt, M. D. Ogilvie, and A. Gottschalk. Correlation between oscillations in ventilation and frequency content of the electroencephalogram. J. Appl. Physiol. 72: 985-992, 1992[Abstract/Free Full Text].
13. Pack, A. I., D. A. Silage, R. P. Millman, H. Knight, E. T. Shore, and D. C. Chung. Spectral analysis of ventilation in elderly subjects awake and asleep. J. Appl. Physiol. 64: 1257-1267, 1988[Abstract/Free Full Text].
14. Phillipson, E. A., and G. Bowes. Control of breathing during sleep. In: Handbook of Physiology. The Respiratory System. Bethesda, MD: Am. Physiol. Soc., 1986, sect. 3, vol. II, pt. 2, chapt. 19, p. 649-690.
15. Reite, M., D. Jackson, R. L. Cahoon, and J. V. Weil. Sleep physiology at high altitude. Electroencephalogr. Clin. Neurophysiol. 38: 463-471, 1975[Medline].
16. Shore, E. T., R. P. Millman, D. A. Silage, D. C. Chung, and A. I. Pack. Ventilatory and arousal patterns during sleep in normal young and elderly subjects. J. Appl. Physiol. 59: 1607-1615, 1985[Abstract/Free Full Text].
17. Trinder, J., J. Dunai, A. Kay, and I. M. Colrain. The contribution of state changes to respiratory instability during sleep onset. In: Sleep-Wakefulness, edited by V. M. Kumar, H. N. Mallick, and U. Nayar. New Delhi: Wiley, 1993, p. 205-211.
18. Trinder, J., W. Whitworth, A. Kay, and P. Wilken. Respiratory instability during sleep onset. J. Appl. Physiol. 73: 2462-2469, 1992[Abstract/Free Full Text].
19. Waggener, T. B., A. R. Stark, B. A. Cohlan, and I. D. Frantz. Apnea duration is related to ventilatory oscillation characteristics in newborn infants. J. Appl. Physiol. 57: 536-544, 1984[Abstract/Free Full Text].
20. Webb, P. Periodic breathing during sleep. J. Appl. Physiol. 37: 899-903, 1974[Free Full Text].
21. Weil, J. V., M. H. Kryger, and C. H. Scoggin. Sleep and breathing at high altitude. In: Sleep Apnea Syndromes, edited by C. Guilleminault, and W. C. Dement. New York: Liss, 1978, p. 119-136.
22. Worsnop, C., A. Kay, J. Trinder, Y. Kim, and R. J. Pierce. The activity of respiratory pump and upper airway muscles at sleep onset (Abstract). Am. J. Respir. Crit. Care Med. 153: A689, 1996.

0161-7567/97 $5.00 Copyright © 1997 the American Physiological Society



This article has been cited by other articles:


Home page
Biol Res NursHome page
B. W. Carlson and V. J. Neelon
Evaluation of Variables to Characterize Respiratory Periodicity during Sleep in Older Adults
Biol Res Nurs, April 1, 2002; 3(4): 176 - 188.
[Abstract] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
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
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Trinder, J.
Right arrow Articles by Kay, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Trinder, J.
Right arrow Articles by Kay, A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online