|
|
||||||||
Departments of Medicine and Physiology, University of Toronto, Toronto, Ontario, Canada M5S 1A8
| |
ABSTRACT |
|---|
|
|
|---|
Sleep-disordered breathing is associated with adverse clinical consequences such as daytime sleepiness and hypertension. The mechanisms behind these associations have been studied in animal models, especially rats, but intermittent stimuli such as hypoxia have been applied without reference to sleep-wake states. To determine mechanisms underlying the adverse physiological consequences of stimuli associated with sleep-disordered breathing requires criteria for detection of sleep-wake states on-line to trigger stimuli only in sleep. This study aimed to develop such a system for freely behaving rats. Twelve rats with implanted electroencephalogram and neck electromyogram electrodes were studied in the light and dark phases. Electroencephalogram frequencies in the high (20-30 Hz) and low (2-4 Hz) frequency bands distinguished non-rapid eye movement (REM) sleep, whereas neck electromyogram distinguished REM. Using these parameters in a simple algorithm led to detection accuracies of 94.5 ± 1.0 (SE) % for wakefulness, 96.2 ± 0.8% for non-REM sleep, and 92.3 ± 1.6% for REM compared with blinded human judgment. The algorithm was then used to trigger hypoxic stimuli only in sleep. Because frequency and amplitude analysis is readily performed using a variety of commercial systems, incorporation of these parameters into such an algorithm will facilitate studies investigating mechanisms underlying the physiological consequences of sleep-related respiratory stimuli in a fashion that more effectively models clinical disorders.
electroencephalogram; electromyogram; obstructive sleep apnea
| |
INTRODUCTION |
|---|
|
|
|---|
DISORDERS OF SLEEP AND BREATHING constitute a major public health problem (32, 36). Obstructive sleep apnea, for example, affects ~4% of adults (47) and causes recurrent asphyxia and arousals from sleep. Sleep apnea syndromes are associated with important clinical consequences such as excessive daytime sleepiness; impaired driving and work performance (15, 19); hypertension (10); increased risk for stroke, angina, and myocardial infarction (42); and suppressed ventilatory and arousal responses to altered blood gases and airway occlusions (9, 24).
There is currently major interest in determining the potential role of stimuli involved with sleep-disordered breathing (e.g., intermittent hypoxia) and the mechanisms underlying the adverse physiological consequences (33). Rats have been the most used animal model for such studies because of their particular suitability for chronic instrumentation and experiments in sleep. For example, the effects of repetitive intermittent or chronic hypoxia applied for several hours during the day on sleep-wake states (28, 34, 40), blood pressure and cardiovascular responses (3, 16, 17, 20, 26), and ventilation have been studied (30, 34). The effects of repetitive arousal stimuli on cardiovascular responses have also been reported in rats (4).
The applicability of such studies, however, to the understanding of the adverse physiological consequences of sleep-related breathing disorders is limited because the intermittent stimuli in all of those studies were applied without reference to sleep-wake states and only for several hours during the day (3, 16, 17, 20, 26, 28, 34, 40). Such protocols are typically employed because nocturnal rodents like rats tend to consolidate their sleep into the light phase, thus increasing the likelihood that stimuli would be applied in sleep. However, the robust ultradian rhythm of the sleep-wake cycle in rats (e.g., Ref. 43 and see RESULTS) ensures that there are prolonged periods of wakefulness even during the light phase, thereby complicating the assumption that the intermittent stimuli are likely to be applied in sleep. Moreover, even if the intermittent hypoxic episodes happened to be applied coincidentally with sleep episodes, it would also be highly unlikely that the offset of hypoxia would also be coincident with arousal as occurs in clinical disorders. Limiting stimuli only to the light phase also restricts the number of stimuli that can potentially be applied in sleep because the same ultradian rhythm of sleep-wake patterns also ensures that there are prolonged periods of sleep even in the dark phase in rats (e.g., Ref. 43 and see RESULTS). Moreover, because hypoxic stimuli also promote wakefulness (28, 34, 40) and can phase shift or depress circadian rhythms (23, 31), it is also conceivable that repetitive hypoxia in the light phase (3, 16, 17, 20, 26) may result in rats simply shifting their sleep into the dark phase when the stimuli are not applied. These caveats, therefore, confound the interpretation that the intermittent hypoxic stimuli used in previous studies were applied in sleep in a fashion that mimics the clinical situation of sleep-disordered breathing. Rather, these observations suggest that criteria to ensure application of stimuli during sleep and removal of stimuli at arousal are a prerequisite to studying the effects of specific respiratory-related stimuli on sleep mechanisms and other physiological processes in a way that most effectively models the clinical sleep-related breathing problems.
Continuous visual scoring of sleep-wake states from polygraphic records and manual application of stimuli in sleep, however, is obviously not practical for chronic animal experiments that can last many hours or days. For these reasons, Bergmann et al. (6, 7) originally pioneered a computer algorithm for on-line detection of sleep-wake states in behaving rats from chronic electroencephalogram (EEG) and neck electromyogram (EMG) recordings. This algorithm has been used extensively in rats to apply stimuli exclusively in sleep in studies of full or partial sleep deprivation (e.g., Refs. 6, 14, 37, 41). This original algorithm uses the amplitudes of the integrated EEG and neck EMG signals to separate sleep-wake states with two EEG signals filtered between 5-15 and 0.5-15 Hz (6, 7). However, original validation of this algorithm was restricted to the light phase of the light-dark cycle and was performed in few rats (7). Validation data were presented only in the light phase in that study, most likely because constant illumination is routinely used in sleep deprivation experiments to dampen or eliminate circadian rhythms (6, 41). However, such criteria for detection of sleep-wake states in constant illumination may not be readily applicable to studies using the normal light-dark cycle. In addition, restricting analysis of EEG frequencies to <15 Hz precludes analysis of a significant portion of EEG frequencies that change across sleep-wake states in rats.
Given these caveats, the aims of the present study were to develop and validate a simple computer algorithm to accurately detect sleep-wake states on-line in freely behaving rats using an algorithm that would be applicable in both the light and dark cycles. Development and validation of this algorithm would also include measures of EEG frequencies across an appropriate range to include those expected to be significantly modulated by sleep-wake states. Triggering of voltage outputs on sleep detection would also fulfill the requirement of a system capable of applying respiratory-related stimuli (e.g., hypoxia) exclusively in sleep and removing those stimuli at arousal. The hypothesis is that parameters detected by frequency and amplitude analysis of the EEG and neck EMG can distinguish sleep-wake states on-line, irrespective of light-dark cycle. Because such analyses are now readily performed using a variety of commercial systems, such a simple sleep-detection algorithm would facilitate studies investigating the impact of respiratory-related stimuli applied in sleep in a fashion that more closely approximates common clinical disorders.
| |
METHODS |
|---|
|
|
|---|
Animals and Surgical Preparation
Studies were performed on 12 adult male Sprague-Dawley rats (mean body wt: 315 g; range: 290-420 g; Charles River Laboratories). Each rat was housed individually, maintained on a 12:12-h light-dark cycle, and had access to food and water ad libitum. Six animals were studied on an 1100-2300 light and 2300-1100 dark cycle, whereas the other six were studied on the reverse schedule. The rats on the 1100-2300 light cycle had been acclimatized to this regimen for an average of 10.0 days (range, 6-14 days), and the rats on the reverse schedule had been acclimatized for an average of 13.5 days (range, 12-21 days).Sterile surgery was performed under general anesthesia induced by intraperitoneal ketamine (85 mg/kg) and xylazine (15 mg/kg). Before surgery, the rats were also given buprenorphine (0.03 mg/kg), atropine (1 mg/kg), and sterile saline (3 ml, 0.9%). To fix body position, the rats were placed in a stereotaxic apparatus (Kopf model 962) with blunt ear bars. An anesthesia mask (18) was placed over the snout, and halothane in air was administered as necessary for the remainder of the surgery. Body temperature was measured with a rectal probe and maintained between 36 and 38°C with a heating pad (BAS, West Lafayette, IN). All surgical procedures were approved by the Animal Care Committee of the University of Toronto.
Two stainless steel screws attached to insulated wire were implanted in the skull over the frontal-parietal cortex to record the EEG. One of these electrodes was placed ~2 mm anterior and 2 mm to the right of bregma, and the other ~3 mm posterior and 2 mm to the left of bregma (22). Two insulated, multistranded stainless steel wires bared at the tips were sutured onto the dorsal cervical neck muscles to record the EMG. Electrode leads were connected to amphenol pins inserted into a miniature plug (STC-89PI-220ABS, Carleton University, Ottawa, Ontario) affixed to the skull with dental acrylic and anchor screws. The rats recovered for an average of 14.1 ± 3.7 days (range, 7-20 days) before the onset of experiments.
Recording Procedures
For electrophysiological recording, a lightweight shielded cable was connected to the plug on the rat's head that was attached to a counterbalanced swivel that permitted free movement. All rats were studied in their home cages in a noise-attenuated environment, free from interruption. For purposes of habituation, they were connected to the recording apparatus at least 1 day before the experiments. The signals were routed to a Grass model 79D polygraph with 7P511 amplifiers. The EEG and neck EMG were filtered between 1 and 100 and between 10 and 1,000 Hz, respectively, and were recorded on chart paper at 5 mm/s.Computerized Analysis of EEG and EMG Signals
A personal computer (IBM-compatible 386, 16 MHz) received the EEG and EMG signals after analog-to-digital conversion at a sampling rate of 300 Hz (Lab Master DMA, Arrow Electronics, Techmar, OH). In 6-s epochs, the computer then analyzed the periods in the EEG signal, then derived the corresponding frequencies, and calculated the percent signal in each of six bandwidths:
2 (0.5-2 Hz),
1 (2-4 Hz),
(4-7.5 Hz),
(7.5-13.5
Hz),
1 (13.5-20 Hz), and
2
(20-30 Hz). The peak-to-peak EEG amplitude and the moving average
of the EMG signal were also determined. EEG analysis was performed
using a modification of the interval histogram method (27). In this analysis, the amplitude of the EEG signal
was divided into 32 equally spaced horizontal slice lines. A period was
measured as the time interval between two points at which the same
slice line crossed consecutive positive-going slopes of the EEG signal
(27). A histogram was then constructed for these
intervals, and from this histogram the percent distribution of
frequencies was calculated. The analysis software used in the present
study is similar to that previously used in dogs (21, 25),
although it was modified slightly in this study for delivery of stimuli
(see below).
The accuracy for detecting frequencies in the EEG was examined using
signals from a function generator (Wavetek, San Diego, CA). Sine waves
equivalent to inputs of 10-500 µV and frequencies of 1-25
Hz were used in accordance with the amplitudes present in the rat EEG
(see RESULTS), and the range of frequency bands was
analyzed (see above). For a 100-µV peak-to-peak sine wave, the
average signal detection accuracy for all frequency bands (
2 through to
2) was 99.3 ± 0.7 (SE) %. The magnitude of the detected EEG signal was also linearly
related to the input voltage (r = 1.000, P < 0.0001, n = 8, range of inputs
10-500 µV). The average coefficient of variation (i.e., 100 × SD/mean) of the measured EEG signal for a given input was <0.6%.
For the EMG, the magnitude of the detected signal was also linearly
related to the input voltage (r = 0.999, P = 0.03) with the average coefficient of variation
<0.2%.
Protocol
The rats were connected to the cable and swivel apparatus at least 1 day before the experiments. On the day of the experiment, the EEG and EMG signals were calibrated on chart and computer, and these signals were then monitored continuously. Except for brief periods to download and backup the data, computerized sampling was continuous for the remainder of the experiment.Analysis
Wakefulness and non-rapid eye movement (non-REM) and rapid eye movement (REM) sleep were determined visually in 10-s epochs from the chart record using standard EEG and EMG criteria (22). Periodic observations of the animal were also performed as an aid to determination of sleep-wake states, and appropriate notations were made on the chart. In summary, in wakefulness, the rats lay quietly, and the EEG displayed low-voltage, high-frequency activity. During non-REM sleep, the EEG was dominated by high-voltage, low-frequency activity. During REM sleep, the EEG was dominated by low-voltage, high-frequency activity, and neck muscle EMG was silent except for occasional muscle twitches. Brief arousals from sleep lasting from 3 to 10 s were also identified from the EEG and EMG signals (1). From the chart records, sleep efficiency (sleep time/recording time), the percentage of wakefulness, non-REM and REM sleep, and the number of arousals per hour were calculated. These values were also calculated from the computerized records, and, in addition, EEG frequencies and EEG and EMG amplitudes were determined from the computer. When the results of the computerized analysis were compared with visual scoring, both records were time aligned, and direct comparisons of the judgments of wakefulness and non-REM and REM sleep were performed. Also, in six rats chosen randomly (3 studied in the dark phase and 3 studied in the light phase), the percentages of wakefulness and non-REM and REM sleep were also determined by a second independent human scorer without reference to the results of the first scorer. For these rats, comparisons between scorers were performed for the data collected over the whole of the 6-h experimental period in each rat.Triggering of Hypoxic Stimuli in Sleep
To determine the suitability of the on-line sleep-detection system to apply intermittent stimuli exclusively in sleep, three chronically instrumented rats were studied after the threshold criteria used to distinguish sleep-wake states were determined on the first day of recording (see RESULTS). The rats were studied in an airtight chamber continually flushed with room air, and, on detection of two consecutive epochs of sleep by the computer, a voltage pulse was generated (4.4 V) that triggered a solenoid valve to switch from room air to 100% nitrogen. The nitrogen caused oxygen levels to fall rapidly, and, when 10% oxygen was reached, a logic circuit switched the inspired gas to continuous 10% oxygen. On detection of two consecutive epochs of wakefulness after arousal from sleep, the voltage returned to 0 V, which switched the inspired gas transiently to 100% oxygen, which caused oxygen levels to rise rapidly. When 21% oxygen was reached, the rat was returned to room air. This cycle of stimuli application and removal was repeated for several hours in the light phase, and the number, duration, and sleep-wake states in which the stimuli were applied and removed were calculated.Statistical Analysis
The analyses performed for each statistical test are included in the text where appropriate. For all comparisons, differences were considered significant if the null hypothesis was rejected at a level of P < 0.05 using a two-tailed test. Where post hoc comparisons were performed after ANOVA with repeated measures, the Bonferroni corrected P value was used to infer statistical significance. For two-way ANOVA, the two factors were light-dark cycle and sleep-wake states (i.e., wakefulness and non-REM and REM sleep). Analyses were performed using Sigmastat (Jandel Scientific, San Rafael, CA). Data are presented as means ± SE, unless otherwise indicated.| |
RESULTS |
|---|
|
|
|---|
The overall sleep architecture was determined from visual analyses of the chart records to facilitate later comparisons with the computerized records. The chart recorder was run continuously for an average of 5 h 48 min ± 13 min for the rats studied in both the light and dark phases. As such, validation data were recorded for ~6 h per rat per day per lighting condition. For the rats studied in the light, the percentage of the recording time spent in wakefulness and non-REM and REM sleep was 26.9 ± 2.7, 66.3 ± 2.5, and 6.7 ± 1.3%, compared with 69.7 ± 2.3, 23.5 ± 1.8, and 6.7 ± 0.7% for the dark, respectively. There was a significant effect of light-dark cycle on sleep-wake states [F(1,10) = 7.68, P = 0.020] with the differences being significant for wakefulness and non-REM sleep [t(10) = 14.84 and 14.04, respectively, both P < 0.001] but not for REM sleep [t(10) = 0.01, P = 0.989].
Determination of Criteria for On-line Sleep Detection
EEG frequencies across sleep-wake states.
Figure 1 shows an example of the raw EEG
and EMG signals across sleep-wake states and the distribution of EEG
frequencies for the group of 12 rats. Note the typical shift to slower
EEG frequencies in non-REM sleep compared with waking and REM. The frequency bands that showed the most obvious changes were the
1 (2-4 Hz) and
(4-7.5 Hz) bands, which
showed clear increases in non-REM sleep compared with waking and REM,
whereas the
2 (20-30 Hz) band showed clear
decreases in non-REM.
|
1 activity was consistently different
across sleep-wake states [F(2,20) = 137.6, P < 0.0001], being significantly increased in
non-REM sleep compared with wakefulness and REM
[t(11) = 11.7 and 16.0, respectively,
both P < 0.001]. Similarly,
activity was
consistently affected by sleep-wake states
[F(2,20) = 90.2, P < 0.0001], being significantly increased in non-REM sleep compared with
waking and REM [t(11) = 11.4 and 11.9, respectively, P < 0.001]. The relative increase in
activity in non-REM sleep, however, was much less than the increase
in
1 activity, which nearly doubled in non-REM sleep
compared with waking or REM (Fig. 1). A decrease in fast EEG
frequencies in non-REM sleep was also reflected by major changes in
2 activity [F(2,20) = 263.5, P < 0.0001], with significant decreases in
non-REM sleep compared with wakefulness and REM
[t(11) = 20.76 and 18.87, respectively, both P < 0.001; see Fig. 1]. Most importantly, these
effects of sleep-wake states on the distribution of EEG frequencies
were not dependent on whether the rats were studied in the light or dark phases (P > 0.309); i.e., the changes in
1,
, and
2 activities in non-REM sleep
were consistent across the light-dark cycle.
These data showed that there were consistent and highly significant
changes in the mean values of
1,
, and
2 frequencies across sleep-wake states. Figure
2 shows, however, that there was some
overlap between populations of individual EEG frequencies within a rat
that somewhat blurred this apparent discrete separation of sleep-wake
states based on average frequency values alone. However, this overlap
between populations of EEG frequencies was reduced by calculating the
ratio of
2 to
1 activity
(
2/
1), which better reflected the overall
shift to increased slower frequencies and decreased faster frequencies
in non-REM sleep compared with when the individual frequency bands were
considered alone (Fig. 2). Compared with wakefulness and REM sleep,
mean
2/
1 in non-REM decreased by 62.6 and
66.9%, respectively.
|
2/
1
activity, EEG amplitude, and neck EMG amplitude across sleep-wake
states. There were consistent and highly significant differences in
2/
1 activity between states
[F(2,20) = 75.66, P < 0.0001], and importantly this effect was independent of whether the
rats were studied in the light or dark phases
[F(1,10) = 1.18, P = 0.302]. The
2/
1 activity was
significantly decreased in non-REM sleep compared with waking and REM
[t(11) = 9.5 and 11.5, respectively, both
P < 0.001; Fig. 3A], whereas values were
not statistically different between REM sleep and wakefulness
[t(11) = 2.0, P > 0.05;
Fig. 3A].
|
EEG amplitudes across sleep-wake states. There was an effect of sleep-wake states on EEG amplitudes [F(2,20) = 123.0, P < 0.0001], with EEG amplitude being significantly increased in non-REM sleep compared with wakefulness and REM [311.2 ± 22.3 vs. 169.7 ± 11.8 and 183.0 ± 12.4 µV, respectively; t(11) = 14.2 and 12.9, respectively, both P < 0.001; Fig. 3B]. However, EEG amplitudes were not statistically different between REM sleep and waking [t(11) = 1.3, P > 0.10; Fig. 3B]. Importantly, there was also an overall effect of light-dark cycle on EEG amplitude [F(1,10) = 6.78, P = 0.026], with overall amplitudes being greater in the dark phase (189.5 ± 17.3 vs. 253.1 ± 17.3 µV; P < 0.05). Nevertheless, there was no significant interaction between sleep-wake state and light-dark cycle on EEG amplitude [F(2,20) = 2.29, P = 0.128], showing that the direction of change in EEG amplitude across sleep-wake states was consistent across light-dark cycle, although the absolute EEG amplitudes differed between lighting conditions.
Neck EMG across sleep-wake states. Neck EMG varied significantly across sleep-wake states [F(2,20) = 153.3, P < 0.0001] and decreased from wakefulness to non-REM sleep [mean decrease = 62.6%; t(11) = 25.8, P < 0.001; Fig. 3C] and was further decreased in REM [mean decrease = 83.7%; t(11) = 34.46, P < 0.001; Fig. 3C]. Importantly, for neck EMG there was no statistically significant interaction between sleep-wake state and light-dark cycle [F(1,10) = 3.44, P = 0.09].
Choice of algorithm and threshold criteria to separate sleep-wake
states.
The above results showed that, irrespective of the light-dark cycle,
the
2/
1 activity distinguished non-REM
sleep from both wakefulness and REM, and neck EMG effectively
distinguished REM sleep from wakefulness. Based on these results,
2/
1 and neck EMG were then incorporated
into a simple two-step algorithm to distinguish among these three
sleep-wake states. This algorithm is shown in Fig.
4.
|
Validation of Criteria for On-line Sleep Detection
Each rat was studied on 2 separate days to determine the accuracy of the algorithm in detecting sleep-wake states. Day 1 was used to select appropriate thresholds for
2/
1 and neck EMG to distinguish
sleep-wake states according to the algorithm in Fig. 4. These
thresholds were chosen after continuous monitoring of sleep-wake states
on chart and the corresponding computer-generated values of
2/
1 and neck EMG.
The
2/
1 thresholds used to distinguish
non-REM sleep from wakefulness and REM were not statistically different
for the rats studied in the dark and light phases [2.32 ± 0.45 (SD) vs. 1.97 ± 0.29; t(10) = 1.62, P = 0.136, unpaired t-test] as were the EMG
thresholds [t(10) = 0.27 and 1.62;
P = 0.80 and 0.40, unpaired t-tests]. After
choosing the appropriate thresholds on day 1, the accuracy
of these thresholds in the detection of sleep-wake states on the
following day was determined via visual scoring of the chart records
for wakefulness, non-REM, and REM. These data were then compared
blindly with the results of the computer judgment (see below).
Accuracy of algorithm to distinguish sleep-wake states and sleep
architecture.
For periods of unequivocal wakefulness, non-REM and REM sleep scored
blindly without reference to the computerized record; the accuracy of
detecting sleep-wake states using the algorithm in Fig. 4 was 94.5 ± 1.0% for wakefulness, 96.2 ± 0.8% for non-REM sleep, and
92.3 ± 1.6% for REM sleep (Table
1). Because these figures refer to
unequivocally identified sleep-wake periods, computer detection
accuracies were also calculated for periods more difficult to identify
definitively by visual inspection, e.g., periods of drowsiness,
transitions to REM sleep, and brief arousals from sleep. Incorporating
these additional, less definitive periods at the transitions between
states, also identified blindly without reference to the computer
record, reduced the detection accuracy of the computer algorithm to
88.0 ± 2.4% for wakefulness, 86.3 ± 2.1% for non-REM
sleep, and 89.6 ± 1.5% for REM sleep. Figure
5 shows that these additional errors in
computerized detection were clustered at values of
2/
1 close to the threshold used to
separate the sleep-wake states, i.e., periods of transitional EEG
activity. Scoring accuracies for periods of wakefulness and non-REM and
REM sleep by computer were the same in the light and dark phases (mean
difference = 0.8 ± 1.7%).
|
|
2/
1 values were occasionally in the range
of non-REM sleep (see Fig. 2).
Overall sleep architecture was also calculated from the computerized
judgments of sleep-wake states and compared with the same periods
matched from the chart record scored by blinded visual inspection. The
percentages of wakefulness determined from the computer and chart were
not statistically different [mean difference = 0.58 ± 4.42 (SD) %; t(11) = 0.457, P = 0.657, paired t-test]. The percentages of non-REM sleep
determined from the computer algorithm and chart records were also not
statistically different [mean difference = 0.58 ± 4.88 (SD)
%; t(11) = 0.414, P = 0.687]. In contrast, there was a small overdetection of REM sleep by
the computer algorithm [mean difference = 2.81 ± 3.20 (SD)
%; t(11) = 3.038, P = 0.01]. Sleep efficiencies determined by the computerized and chart
records were not statistically different [mean difference = 0.58 ± 4.42 (SD) %; t(11) = 0.457, P = 0.657, paired t-test]. Unlike sleep
architecture, however, the numbers of arousals determined by
computerized detection were different from the number determined by the
standard visual criteria. Over the recording period, an average of
113.7 ± 13.6 brief arousals per hour (i.e., those lasting one
computer epoch of 6 s) were detected by the computer, whereas only
34.3 ± 8.0 were detected visually
[t(11) = 10.654, P < 0.001, paired t-test]. For arousals lasting two computer
epochs, an average of 45.3 ± 10.5 arousals per hour were detected
by computer, whereas 20.5 ± 3.0 were detected visually
[t(11) = 2.806, P = 0.017, paired t-test].
For overall judgments of sleep-wake states, the agreement between the
two independent human scorers was 96.2 ± 1.7% for wakefulness, 93.2 ± 1.9% for non-REM sleep, and 98.3 ± 1.6% for REM
sleep. The discrepancies between scorers were due to determination of transitions from quiet wakefulness to non-REM sleep and determination of brief arousals from sleep.
Application to apply stimuli exclusively in sleep.
Figure 6 shows an example of the hypoxic
stimuli applied in sleep and removal of hypoxia at arousal. A total of
270 hypoxic stimuli were applied in the three rats over a 4-h period in
the light phase (1030-1430). Of these stimuli, 93.7% were
correctly applied in sleep (80.4% in non-REM sleep and 13.3% in REM),
2.2% were incorrectly applied in wakefulness, and 4.1% were applied in drowsiness as judged by an EEG pattern intermediate between light
sleep and quiet wakefulness. Furthermore, 94.4% of these hypoxic
stimuli were correctly removed at arousal from sleep, whereas 4.4%
were incorrectly removed in sleep (3.7% in non-REM sleep and 0.7% in
REM), and 0.7% were removed in drowsiness. The time spent in 10%
O2 after the triggering of stimuli averaged 45.6 ± 3.37 s for non-REM sleep and 72.2 ± 8.7 s for REM
sleep.
|
| |
DISCUSSION |
|---|
|
|
|---|
This study describes an algorithm for on-line detection of sleep-wake states in rats and its use in applying hypoxic stimuli exclusively in sleep. In previous studies in rats, application of episodic hypoxia has typically been performed for several hours during the day without reference to sleep-wake states (e.g., Refs. 3, 16, 17, 20, 26). For the reasons outlined in the Introduction, applying stimuli without reference to sleep has major drawbacks if such studies are designed to mimic episodic hypoxia associated with sleep-disordered breathing and determine the mechanisms underlying the adverse physiological consequences.
EEG and EMG Criteria to Recognize Sleep-Wake States On-line in Rats
EEG frequencies in the
1 (2-4 Hz) and
2 (20-30 Hz) bands showed consistent changes in
non-REM sleep compared with wakefulness and REM sleep. This increase in
slow EEG frequencies in non-REM sleep, and decrease in fast
frequencies, was better reflected in the
2/
1, which showed a larger percent change
in non-REM sleep compared with when the other frequency bands were
considered alone (Figs. 1 and 3). The
2/
1
also varied consistently with sleep-wake states, but not light-dark
cycle, making it a robust parameter to distinguish non-REM sleep. The
experimentally derived threshold values of
2/
1 used in the present sleep-detection algorithm were not statistically different between the light and dark
phases. These results highlight two important factors that distinguish
this algorithm from that described by Bergmann et al. (6,
7). In those previous studies, the EEG was analyzed <15 Hz,
such that significant changes in
2 activity could not be
observed, and validation data were only presented in the light phase
(6, 7).
Like
activity,
also increased in non-REM sleep compared
with REM (Fig. 1). This unexpected result may have been because of the
frequency range used to distinguish
activity. The range of
4-7.5 Hz was based on the range of Kuwahara et al.
(27), whose method we also followed for our EEG analysis.
However, unlike for the
band, there appears to be no uniform
frequency range used to define
, e.g., 5-10, 5-15,
6-9, 5.5-8, and 6.25-9 Hz have been used (6, 7,
39, 43). In our case, we did not see an increase in
activity
in REM sleep compared with non-REM in the 4- to 7.5-Hz range, but an
increase in EEG frequencies in REM was observed in the 7.5- to 13.5-Hz
band (Fig. 1), which partially overlaps the
range used by some
investigators (6, 7, 39, 43). A second set of EEG
electrodes placed more medially may have improved our detection of
rhythm (6, 7). Nevertheless, with our electrode
placements, the EEG signals in REM sleep could be clearly distinguished
from wakefulness (Fig. 1) and had the same visual appearance as in
other studies. Our choice of one set of EEG electrodes was simply for
practical reasons to reduce the number of electrodes during surgery.
However, despite the use of only one set of EEG electrodes, the
algorithm was accurate in detecting sleep-wake states and applying
hypoxic stimuli in sleep. A second set of electrodes for long-term
studies, however, has the advantage of redundancy in case one pair
becomes dislodged or malfunctions.
EEG amplitude also showed characteristic and highly significant changes in non-REM sleep (Fig. 3). However, EEG amplitudes in non-REM sleep were affected by light-dark cycle. This observation fits with the demonstrated time-of-day dependence of homeostatic sleep regulation that produces high-intensity, non-REM sleep at the beginning of the sleep phase after prolonged wakefulness, which then diminishes over time (8, 43). Nevertheless, this observation also shows that EEG amplitude would not be a reliable criterion to distinguish sleep-wake states in rats in experiments that cover several hours, because the ability to detect a difference between non-REM sleep and the other states would vary with time of day and the light-dark cycle. This result highlights another important factor that distinguishes the present algorithm from that described by Bergmann et al. (6, 7), where integrated EEG amplitude is a primary criterion to separate sleep-wake states. The choice of EEG amplitude in those previous studies resulted from the original validations being performed in the light phase (7), because constant illumination is routinely used in sleep deprivation experiments to dampen circadian rhythms (6, 41). Accordingly, the accuracy and/or general applicability of EEG amplitude as a robust criterion was not apparent in those studies, as validation data were not presented for the dark phase (7). Based on our analysis in the light and dark phases, EEG amplitude was not chosen as a criterion in the decision tree of our algorithm (Fig. 4). Nevertheless, EEG amplitude does feature in one part of the EEG analysis because it allows the calculation of EEG frequencies from the distribution of periods (see METHODS).
Therefore, our EEG data showed that the
2/
1 EEG activity was a robust criterion
to distinguish non-REM sleep from wakefulness or REM in both the light
and dark phases. However, no other EEG parameter effectively
distinguished between the sleep-wake states across light-dark cycle,
especially wakefulness from REM sleep, which shared similar overall EEG
characteristics (Figs. 1-3). In contrast, neck EMG showed reliable
and highly consistent decreases in activity from wakefulness to REM
sleep (Fig. 3) that were independent of the light-dark cycle. As such,
the
2/
1 EEG activities and the neck EMG
were used in combination in a simple two-step algorithm to separate
sleep-wake states on-line after appropriate thresholds were chosen.
On-line Detection of Sleep-Wake States
There are a variety of techniques to analyze EEG and EMG signals, such as spectral and interval histogram frequency analyses, autoregressive modeling, fuzzy logic, neural network approaches, and detection of EEG and EMG amplitudes after integration (e.g., see Refs. 6, 11, 35, 38, 39). Review of the major differences and assumptions among techniques is outside the scope of the present discussion. However, it is worth mentioning that the interval histogram method was chosen for EEG analysis in this study because it does not assume that the EEG signal is stationary, unlike another commonly used technique such as fast-Fourier transform. However, this distinction may not be overly significant for practical purposes in the range of frequencies that are important in this algorithm. For example, although differences among techniques exist, robust changes in
and
frequencies have been
recorded across the night with both period amplitude analysis and
fast-Fourier transform (44). Because one of the main aims
of this paper was to identify EEG frequency characteristics that can be
used to distinguish sleep-wake states on-line, for the purposes of
applying stimuli exclusively in sleep, our use of these frequency data
led to the present algorithm that fulfilled these criteria. It remains
to be determined whether this same algorithm can be used with other frequency-detection software using different techniques.
Using the algorithm shown in Fig. 4, the accuracy of detecting unequivocally identified sleep-wake states was 94.5% for wakefulness, 96.2% for non-REM sleep, and 92.3% for REM sleep. This degree of accuracy, especially for REM sleep, is similar or improved compared with other computerized systems typically employed in rodents (e.g., Refs. 5, 7, 13, 39, 45, 46) and also compared with an earlier, similar algorithm used in dogs (21, 25).
In the present study, the accuracies of computerized detection were also calculated for periods more difficult to identify definitively by visual criteria, such as drowsiness, transitions to REM sleep, and brief arousals. Incorporating these additional less definitive periods reduced the detection accuracy of the computer algorithm to 88.0% for wakefulness, 86.3% for non-REM sleep, and 89.6% for REM sleep. Although including all available data reduced computerized detection accuracies compared with visual scoring, this was a necessary and required comparison because such continuous data would occur during an on-line experiment. Nevertheless, this level of computer accuracy is still improved or is at least comparable to other systems that typically did not include such transitional states in the analyses and validation (39). Indeed, a significant strength of this study is that the parameters used in the algorithm to separate sleep-wake states were derived from extensive analysis of EEG and EMG signals in a relatively large number of rats studied in both the light and dark phases. This degree of validation is more stringent than that described in most previous studies (7, 39), and its importance is highlighted by the observation that the analysis revealed parameters that were, and were not, robust enough to distinguish sleep-wake states in long-term experiments covering both the light and dark cycles. As such, it was not sufficient that a parameter such as EEG amplitude showed consistent and highly significant changes across sleep-wake states; these changes had to persist across states and be independent of the light-dark cycle to be identified and used as robust criteria for sleep detection.
Although validation in the present study was more systematic than in
most previous experiments, validation in the light and dark phases for
each rat was for only 6 h/day. However, we do not think that this 6-h
continuous validation (i.e., rather than 24 h) significantly
detracts from the main results. First, the major outcome of the study
was that the threshold criteria used to separate sleep-wake states
using the computer algorithm were sufficiently accurate to apply
hypoxic stimuli exclusively during sleep. Because these threshold
criteria also did not change significantly between the light and dark
phases, the ability of the system to detect sleep and deliver stimuli
only in sleep would not be expected to change appreciably over time to
significantly affect most experimental protocols. Moreover, the 6-h
validation was long enough for large changes in
activity to occur
over the sleep phase (8, 43), and yet the threshold
criteria for sleep-detection were still accurate and similar to the
accuracies at other times of the light-dark cycle compared with visual
scoring. As such, although non-REM sleep may have higher
activity
at the onset of the light phase, the threshold criteria used in the
computer algorithm to actually distinguish non-REM sleep from the other
states remained sufficiently stable. By analogy, although deep non-REM
sleep is dominant in the early phase of sleep in humans, the visual
criteria used to actually determine sleep onset and distinguish it from
wakefulness does not vary as the same visual criteria are used to
separate sleep-wake states across the night.
The percentages of wakefulness and non-REM sleep scored by the computer and the blinded human scorer were indistinguishable, although the percentage of REM sleep was overdetected by the computer by an average of 2.8%. In contrast, arousals from sleep were detected significantly more often than with visual scoring. Such overdetection of arousals has also been observed in other automated systems (e.g., Ref. 12). It remains to be determined whether this difference in arousal detection is due to conventional visual scoring methods (e.g., Ref. 1) being insensitive to the subtle EEG changes that may indicate a true transient waking EEG, or whether computerized detection methods are oversensitive to short-term EEG fluctuations.
Errors in computerized detection tended to be clustered at the threshold values used to separate sleep-wake states (Fig. 5), i.e., periods of transitional activity. Errors also tended to occur at times when EEG activity became temporally dissociated from neck EMG. It is normal, however, for EEG and EMG signals to exhibit somewhat different time courses before achieving the typical patterns stereotypical of wakefulness and non-REM and REM sleep (2, 29). In the present study, such normal variations in neck EMG activity were observed, with some low-EMG activity in wakefulness being in the range more typical of sleep, and conversely with some high-EMG activity in sleep being in the range more typical of waking (Fig. 2F). Sudden changes in body position, especially with arousal, were also associated with bursts of neck EMG activity that decayed slowly compared with the development of sleep, as judged by EEG criteria. These observations highlight the difficulty in constraining the dynamic changes exhibited by EEG and EMG signals across the full spectrum of behavioral states to one of only three classes of sleep-wake states and is a problem encountered in all computerized sleep-detection systems (see Ref. 35 for discussion). Nevertheless, computerized sleep detection with the algorithm shown in Fig. 4 was sufficiently high to comply with the primary aim of this study, i.e., to develop a system capable of applying hypoxic stimuli in periods of sleep and removing those stimuli at arousal (Fig. 6).
Summary
The results of the present study show that EEG frequencies in the high (20-30 Hz)- and low (2-4 Hz)-frequency bands effectively distinguished non-REM sleep from wakefulness and REM in freely behaving chronically instrumented rats. In addition, neck EMG effectively distinguished REM sleep from wakefulness. Importantly, these criteria were sufficiently robust to be applicable in both the light and dark cycles. Incorporation of these parameters into a simple two-step algorithm then showed that sleep-wake states could be accurately detected on-line. The algorithm was also coupled to a system for triggering hypoxic stimuli exclusively in sleep. Because the ability to perform frequency and amplitude analysis of the EEG and neck EMG signals is readily performed using a variety of commercial systems, incorporation of these parameters into such an algorithm will facilitate studies investigating the consequences of respiratory stimuli applied exclusively in sleep in such a way that it more effectively mimics sleep-related breathing disorders.| |
ACKNOWLEDGEMENTS |
|---|
The authors respectfully remember Hanan Al-Mijalli for kind assistance. Drs. Richard Stephenson and Eliot Phillipson, Departments of Physiology and Medicine, respectively, are thanked for critical reading of this manuscript.
| |
FOOTNOTES |
|---|
This work was supported by the Medical Research Council (MRC) of Canada (MT-15563), Ontario Thoracic Society, and Canada Foundation for Innovation. R. L. Horner is a recipient of a MRC of Canada Scholarship.
Address for reprint requests and other correspondence: R. L. Horner, Rm. 6368, Medical Sciences Bldg., 1 Kings College Circle, Toronto, Ontario, Canada M5S 1A8 (E-mail: richard.horner{at}utoronto.ca).
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 2 March 2000; accepted in final form 9 January 2001.
| |
REFERENCES |
|---|
|
|
|---|
1.
American Sleep Disorders Association.
EEG arousals: scoring rules and examples.
Sleep
15:
173-184,
1992[Medline].
2.
Aston-Jones, G,
and
Bloom FE.
Activity of norepinephrine-containing locus coeruleus neurons in behaving rats anticipates fluctuations in the sleep-waking cycle.
J Neurosci
1:
876-886,
1981[Abstract].
3.
Bakehe, M,
Miramand JL,
Chambille B,
Gaultier C,
and
Escourrou P.
Cardiovascular changes during acute episodic repetitive hypoxic and hypercapnic breathing in rats.
Eur Respir J
8:
1675-1680,
1995[Abstract].
4.
Bao, G,
Metreveli N,
and
Fletcher EC.
Acute and chronic blood pressure response to recurrent acoustic arousal in rats.
Am J Hypertens
12:
504-510,
1999[Web of Science][Medline].
5.
Benington, JH,
Kodali SK,
and
Heller HC.
Scoring transitions to REM sleep in rats based on the EEG phenomena of pre-REM sleep: an improved analysis of sleep structure.
Sleep
17:
28-36,
1994[Web of Science][Medline].
6.
Bergmann, BM,
Kushida CA,
Everson CA,
Gilliland MA,
Obermeyer W,
and
Rechtschaffen A.
Sleep deprivation in the rat. II. Methodology.
Sleep
12:
5-12,
1989[Web of Science][Medline].
7.
Bergmann, BM,
Winter JB,
Rosenberg RS,
and
Rechtschaffen A.
NREM sleep with low-voltage EEG in the rat.
Sleep
10:
1-11,
1987[Web of Science][Medline].
8.
Borbély, AA.
Sleep homeostasis and models of sleep regulation.
In: Principles and Practice of Sleep Medicine (2nd ed.), edited by Kryger MH,
Roth T,
and Dement WC.. Philadelphia, PA: Saunders, 1994, p. 309-320.
9.
Brooks, D,
Horner RL,
Kimoff RJ,
Kozar LF,
Render-Teixeira CL,
and
Phillipson EA.
Effect of obstructive sleep apnea versus sleep fragmentation on acute responses to airway occlusion in the dog.
Am J Respir Crit Care Med
155:
1609-1617,
1997[Abstract].
10.
Brooks, D,
Horner RL,
Kozar LF,
Render-Teixeira CL,
and
Phillipson EA.
Obstructive sleep apnea as a cause of systemic hypertension: evidence from a canine model.
J Clin Invest
99:
106-109,
1997[Web of Science][Medline].
11.
Carley, DW,
Önal E,
and
Lopata M.
Quantitative analyses of state of vigilance in periodic breathing.
In: Sleep and Respiration in Aging Adults, , edited by Kuna ST,
Suratt PM,
and Remmers JE.. New York: Elsevier, 1991, p. 287-296.
12.
De Carli, F,
Nobili L,
Gelcich P,
and
Ferrillo F.
A method for the automatic detection of arousals during sleep.
Sleep
22:
561-572,
1999[Web of Science][Medline].
13.
Feng, P,
and
Vogel GW.
A new method for continuous long-term polysomnographic recording of neonatal rats.
Sleep
23:
9-14,
2000[Web of Science][Medline].
14.
Feng, P,
Vogel GW,
Obermeyer W,
and
Kinney GW.
An instrumental method for long-term continuous REM sleep deprivation of neonatal rats.
Sleep
23:
175-183,
2000[Web of Science][Medline].
15.
Findley, LJ,
Levinson MP,
and
Bonnie RJ.
Driving performance and automobile accidents in patients with sleep apnea.
Clin Chest Med
13:
427-435,
1992[Web of Science][Medline].
16.
Fletcher, EC,
and
Bao G.
The rat as a model of chronic recurrent episodic hypoxia and effect upon systemic blood pressure.
Sleep
19:
S210-S212,
1996[Web of Science][Medline].
17.
Fletcher, EC,
Lesske J,
Qian W,
Miller CC, III,
and
Unger T.
Repetitive, episodic hypoxia causes diurnal elevation of blood pressure in rats.
Hypertension
19:
555-561,
1992
18.
Freedman, NL.
A sliding gas-delivery system and ventilating chamber for small animals fixed in standard stereotaxic instruments.
Behav Res Methods Instrum Comput
24:
423-425,
1992[Web of Science].
19.
George, CF,
Nickerson PW,
Hanley PJ,
Millar TW,
and
Kryger MH.
Sleep apnoea patients have more automobile accidents.
Lancet
2:
447,
1987[Web of Science][Medline].
20.
Greenberg, HE,
Sica A,
Batson D,
and
Scharf SM.
Chronic intermittent hypoxia increases sympathetic responsiveness to hypoxia and hypercapnia.
J Appl Physiol
86:
298-305,
1999
21.
Horner, RL,
Brooks D,
Kozar LF,
Leung E,
Hamrahi H,
Render-Teixeira CL,
Makino H,
Kimoff RJ,
and
Phillipson EA.
Sleep architecture in a canine model of obstructive sleep apnea.
Sleep
21:
847-858,
1998[Web of Science][Medline].
22.
Horner, RL,
Sanford LD,
Annis D,
Pack AI,
and
Morrison AR.
Serotonin at the laterodorsal tegmental nucleus suppresses rapid-eye-movement sleep in freely behaving rats.
J Neurosci
17:
7541-7552,
1997
23.
Jarsky, T,
and
Stephenson R.
Effects of hypoxia and hypercapnia on circadian rhythms in the golden hamster (Mesocricetus auratus).
J Appl Physiol
89:
2130-2138,
2000
24.
Kimoff, RJ,
Brooks D,
Horner RL,
Kozar LF,
Render-Teixeira CL,
Champagne V,
Mayer P,
and
Phillipson EA.
Ventilatory and arousal responses to hypoxia and hypercapnia in a canine model of obstructive sleep apnea.
Am J Respir Crit Care Med
156:
886-894,
1997
25.
Kimoff, RJ,
Makino H,
Horner RL,
Kozar LF,
Lue F,
Slutsky AS,
and
Phillipson EA.
A canine model of obstructive sleep apnea: model description and preliminary application.
J Appl Physiol
76:
1810-1817,
1994
26.
Kraiczi, H,
Magga J,
Sun XY,
Ruskoaho H,
Zhao X,
and
Hedner J.
Hypoxic pressor response, cardiac size, and natriuretic peptides are modified by long-term intermittent hypoxia.
J Appl Physiol
87:
2025-2031,
1999
27.
Kuwahara, H,
Higashi H,
Mizuki Y,
Matsunari S,
Tanaka M,
and
Inanaga K.
Automatic real-time analysis of human sleep stages by an interval histogram method.
Electroencephalogr Clin Neurophysiol
70:
220-229,
1988[Web of Science][Medline].
28.
Laszy, J,
and
Sarkadi A.
Hypoxia-induced sleep disturbance in rats.
Sleep
13:
205-217,
1990[Web of Science][Medline].
29.
Mahowald, MK,
and
Schenck CH.
REM sleep behavior disorder.
In: Principles and Practice of Sleep Medicine (2nd ed.), edited by Kryger MH,
Roth T,
and Dement WC.. Philadelphia, PA: Saunders, 1994, p. 574-588.
30.
Megirian, D,
Ryan AT,
and
Sherrey JH.
An electrophysiological analysis of sleep and respiration of rats breathing different gas mixtures: diaphragmatic muscle function.
Electroencephalogr Clin Neurophysiol
50:
303-313,
1980[Web of Science][Medline].
31.
Mortola, JP,
and
Seifert EL.
Hypoxic depression of circadian rhythms in adult rats.
J Appl Physiol
88:
365-368,
2000
32.
National Commission on Sleep Disorders Research.
Wake Up America: A National Sleep Alert. Washington, DC: Government Printing Office, 1993, p. 302.
33.
National Institutes of Health. National Heart, Lung, and
Blood Institute. Oxygen Sensing During Intermittent Hypoxia.
Request for Applications
(http://grants.nih.gov/grants/guide/rfa-files/RFA-HL-00-004.html).
Bethesda, MD: NIH, Nov 29, 1999. (RFA: HL-00-004)
34.
Pappenheimer, JR.
Sleep and respiration of rats during hypoxia.
J Physiol (Lond)
266:
191-207,
1977
35.
Pardey, J,
Roberts S,
Tarassenko L,
and
Stradling J.
A new approach to the analysis of the human sleep/wakefulness continuum.
J Sleep Res
5:
201-210,
1996[Web of Science][Medline].
36.
Phillipson, EA.
Sleep apnea. A major public health problem.
N Engl J Med
328:
1271-1273,
1993
37.
Rechtschaffen, A,
Gilliland MA,
Bergmann BM,
and
Winter JB.
Physiological correlates of prolonged sleep deprivation in rats.
Science
221:
182-184,
1983
38.
Robert, C,
Guilpin C,
and
Limoge A.
Review of neural network applications in sleep research.
J Neurosci Methods
79:
187-193,
1998[Web of Science][Medline].
39.
Robert, C,
Guilpin C,
and
Limoge A.
Automated sleep staging systems in rats.
J Neurosci Methods
88:
111-122,
1999[Web of Science][Medline].
40.
Ryan, AT,
and
Megirian D.
Sleep-wake patterns of intact and carotid sinus nerve sectioned rats during hypoxia.
Sleep
5:
1-10,
1982[Web of Science][Medline].
41.
Shaw, PJ,
Bergmann BM,
and
Rechtschaffen A.
Effects of paradoxical sleep deprivation on thermoregulation in the rat.
Sleep
21:
7-17,
1998[Web of Science][Medline].
42.
Shepard, JW, Jr.
Hypertension cardiac arrhythmias myocardial infarction and stroke in relation to obstructive sleep apnea.
Clin Chest Med
13:
437-458,
1992[Web of Science][Medline].
43.
Trachsel, L,
Tobler I,
and
Borbély AA.
Electroencephalogram analysis of non-rapid-eye-movement sleep in rats.
Am J Physiol Regulatory Integrative Comp Physiol
255:
R27-R37,
1988
44.
Uchida, S,
Feinberg I,
March JD,
Atsumi Y,
and
Maloney T.
A comparison of period amplitude analysis and FFT power spectral analysis of all-night human sleep EEG.
Physiol Behav
67:
121-131,
1999[Medline].
45.
Van Gelder, RN,
Edgar DM,
and
Dement WC.
Real-time automated sleep scoring: validation of a microcomputer-based system for mice.
Sleep
14:
48-55,
1991[Web of Science][Medline].
46.
Vivaldi, EA,
Ocampo A,
Wyneken U,
Roncagliolo M,
and
Zapata AM.
Short-term homeostasis of active sleep and the architecture of sleep in the rat.
J Neurophysiol
72:
1745-1755,
1994
47.
Young, T,
Palta M,
Dempsey J,
Skatrud J,
and
Badr S.
The occurrence of sleep-disordered breathing among middle-aged adults.
N Engl J Med
328:
1230-1235,
1993
This article has been cited by other articles:
![]() |
H. W. Steenland, H. Liu, and R. L. Horner Endogenous Glutamatergic Control of Rhythmically Active Mammalian Respiratory Motoneurons In Vivo J. Neurosci., July 2, 2008; 28(27): 6826 - 6835. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Liu, S. Sood, H. Liu, and R. L Horner Opposing muscarinic and nicotinic modulation of hypoglossal motor output to genioglossus muscle in rats in vivo J. Physiol., June 15, 2005; 565(3): 965 - 980. [Abstract] [Full Text] [PDF] |
||||
![]() |
T.C.B. Schutz, M.L. Andersen, and S. Tufik Influence of Temporomandibular Joint Pain on Sleep Patterns: Role of Nitric Oxide Journal of Dental Research, September 1, 2004; 83(9): 693 - 697. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Stephenson, K. S. Liao, H. Hamrahi, and R. L Horner Circadian rhythms and sleep have additive effects on respiration in the rat J. Physiol., October 1, 2001; 536(1): 225 - 235. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |