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


     


J Appl Physiol 99: 2036-2044, 2005. First published July 14, 2005; doi:10.1152/japplphysiol.01040.2004
8750-7587/05 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
99/5/2036    most recent
01040.2004v1
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 ISI Web of Science
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 ISI Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Thomson, S.
Right arrow Articles by Semple, S. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Thomson, S.
Right arrow Articles by Semple, S. J.

HIGHLIGHTED TOPICS
Physiology and Pathophysiology of Sleep Apnea

Ventilation is unstable during drowsiness before sleep onset

Stuart Thomson,1 Mary J. Morrell,2 Jeremy J. Cordingley,2 and Stephen J. Semple1

1Department of Respiratory Medicine, Charing Cross Hospital Campus, and 2Clinical and Academic Unit of Sleep and Breathing, Royal Brompton Hospital, National Heart and Lung Institute, Imperial College, London, United Kingdom

Submitted 24 September 2004 ; accepted in final form 20 June 2005


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Ventilation is unstable during drowsiness before sleep onset. We have studied the effects of transitory changes in cerebral state during drowsiness on breath duration and lung volume in eight healthy subjects in the absence of changes in airway resistance and fluctuations of ventilation and CO2 tension, characteristic of the onset of non-rapid eye movement sleep. A volume-cycled ventilator in the assist control mode was used to maintain CO2 tension close to that when awake. Changes in cerebral state were determined by the EEG on a breath-by-breath basis and classified as alpha or theta breaths. Breath duration and the pause in gas flow between the end of expiratory airflow and the next breath were computed for two alpha breaths which preceded a theta breath and for the theta breath itself. The group mean (SD) results for this alpha-to-theta transition was associated with a prolongation in breath duration from 5.2 (SD 1.3) to 13.0 s (SD 2.1) and expiratory pause from 0.7 (SD 0.4) to 7.5 s (SD 2.2). Because the changes in arterial CO2 tension (PaCO2) are unknown during the theta breaths, we made in two subjects a continuous record of PaCO2 in the radial artery. PaCO2 remained constant from the alpha breaths through to the expiratory period of the theta breath by which time the duration of breath was already prolonged, representing an immediate and altered ventilatory response to the prevailing PaCO2. In the eight subjects, the CO2 tension awake was 39.6 Torr (SD 2.3) and on assisted ventilation 38.0 Torr (1.4). We conclude that the ventilatory instability recorded in the present experiments is due to the apneic threshold for CO2 being at or just below that when awake.

apneic threshold for CO2; sleep-related rise in PCO2; arterial to end-tidal PCO2 gradient


IN THE STUDY OF RESPIRATORY variation during the onset of sleep, Bulow (3) recorded that "in some cases periods of total apnea for 10 to 20 s or more occurred occasionally for more than a minute." The end-tidal CO2 tension (PETCO2) at the start of the apneas was described as being relatively low, nearer that when awake than asleep. Subsequently, Skatrud and Dempsey (30) showed, in two normal subjects asleep and on assisted ventilation, that apnea occurred when the PETCO2 was lowered by 1–2 Torr below that when awake. This instability of pulmonary ventilation at sleep onset has been variously ascribed to a fluctuating cerebral state (asleep vs. awake) (24, 25, 32), changes in chemical feedback (17, 21, 25), changes in upper airway resistance (15, 17), and changes in the effectiveness of ventilation in altering arterial blood-gas tensions (plant gain) (4, 20). The fluctuation in cerebral state (24, 32) has been proposed to be the dominant cause of the instability, with other factors playing a contributory role, which may vary in magnitude between individuals (15, 22).

Trinder et al. (32) studied the variation of ventilation at the onset of sleep on a breath-by-breath basis. They concluded that the instability of ventilation at sleep onset was primarily due to the transitions in cerebral state from alpha to theta activity and that the magnitude of the changes was due to the prevailing metabolic drive [as indicated by PETCO2 and end-tidal PO2 (PETO2) levels].

The interpretation of the effect of a change in cerebral state on the control of ventilation is complicated by the fluctuations of ventilation, PCO2, and upper airway resistance, characteristic of sleep onset. Thus it has been shown that the magnitude of a change in ventilation at a transition between cerebral states is critically dependent on the duration of the cerebral state (number of breaths) that precedes that transition (32). Furthermore, it is not always clear whether fluctuations in PCO2 are a determinant of, or determined by, changes in ventilation at a transition between cerebral states. This is particularly relevant at arousal, where the increase in ventilation is in excess of ventilatory demand, and the fall in PCO2 is due to the hyperventilation associated with arousal (11). Finally, there is the sleep-related fall in alveolar ventilation and rise in PCO2 at sleep onset (25), which may prevent the occurrence of apneas and periodic breathing with changes in cerebral state. It is our hypothesis that the sleep-related rise in PCO2 does so by raising the PCO2 above the apneic threshold for CO2.

The purpose and design of our experiments was to determine the change in ventilatory motor output at the alpha-to-theta transition, independent of fluctuations in ventilation, PCO2, and upper airway resistance. We used assisted ventilation to clamp PETCO2 and PETO2 constant during alpha activity and so determine the effect of a change in cerebral activity, independent of fluctuating chemical feedback or feedforward. To avoid changes in tidal volume (VT) from alterations in upper airway resistance, a volume-cycled ventilation was used, set at a predetermined volume. The ventilator was triggered by the subject (assist/control mode), and the assessment of ventilatory drive after the alpha-theta transition was made from the presence or absence of a breath triggered by the subject.

The development of a rapid-responding intra-arterial pH electrode in our laboratory enabled us, for the first time, to record changes in arterial PCO2 (PaCO2) in two subjects during periods of ventilatory instability associated with an apnea, when PETCO2 cannot be used as a surrogate for PaCO2.

We observed, at constant PETCO2, that the transition from alpha to theta activity produced a significant prolongation of breath duration (TT), often with apneas lasting 10–20 s. A preliminary report on these findings has been presented (28).


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Subjects.   Volunteers were recruited from the staff of the medical school. All were men with no history of snoring or ill health. The mean (SD) body mass index was 22.37 kg/m2 (SD 3.04). None of the subjects was aware of the purpose of the experiments. Ethical approval for the experiments was obtained from the local Research Ethics Committee, and full consent was obtained from each subject. The two subjects in whom an arterial record was obtained had normal spirometric lung function studies. Arterial cannulation was carried out under local anesthetic with full aseptic precautions.

Measurements.   The EEG was monitored using two EEGs (C3-A2, C4-A1), two electroocculograms (left and right eye), and electromyogram of submental muscle. These were recorded and calibrated via analog amplifiers (model 18c-P-23, Grass Instrument, Warwick, MA). Airflow was measured by using a Fleisch no. 2 pneumotachometer with a differential pressure transducer (MP 45 ± 2 cmH2O, Validyne, Northridge, CA) attached to a nasal mask. Expired air was sampled through a flexible tube placed within the nostril. From this, PETCO2 and PETO2 were measured by using a mass spectrometer (QP9000, P. K. Morgan, Kent, UK). During apneas and breathing at slow respiratory rates, there was occasional contamination of expired airflow from air within the mask. In these circumstances, PETCO2 was determined by extrapolation of the mid- to terminal part of the expired PCO2 to the end of expiration. The dead space of the mask plus pneumotachygraph was ~60 ml.

Arterial pH and PaCO2.   PaCO2 was continuously measured by using a fast-responding intra-arterial pH electrode, which has been described in detail elsewhere (5). Briefly, the sensor of the system consisted of a pH-sensitive plastic membrane adherent to the tip of a catheter, which was threaded down a radial artery catheter (gauge 18) protruding 2–3 mm into the arterial lumen. The internal reference electrode was Ag-AgCl, housed in a 2-ml syringe containing 0.9% NaCl, attached to the radial artery catheter. To avoid thrombus formation and deposition of proteinaceous material on the pH sensor, a slow infusion of dilute heparin (2,500 units in 500 ml of 0.9% NaCl) was maintained by a flush device at a rate of 3 ml/h. The 10–90% response time of the electrode system has been determined with phosphate buffers and was 0.17 s (0.01) (n = 43) (6). The voltage output of the electrode system was recorded by using an optical isolated electrometer connected to a locally constructed high-impedance amplifier. The electrode system was calibrated in vivo in absolute pH units (pHa) from measurement of arterial blood samples taken during the experiments. Changes in PaCO2 were derived from the relationship between millivolts and pH units from the Nernst equation.

Derivation of PaCO2 from the record of pHa.   Over short time periods at rest, pHa is determined solely by changes in PaCO2 (1). The output of the electrode system was, therefore, calibrated in terms of PCO2 using pHa, the base excess of the blood sample, and the current hemoglobin level of the subject. A blood-gas calculator (29) was used to derive the PaCO2 at intervals of 0.05 units over the full range of pHa covered during an experiment with linear extrapolation between these intervals. The pH, PCO2, and PO2 of the blood samples were measured on a blood-gas machine (CIBA-Corning 248), which also provided the base excess of the samples.

Assisted ventilation.   This was provided by a volume-cycled ventilator (Life Care, PLV-100, Lafayette, IN), triggered by the subject at the start of an inspiratory effort (assist/control mode). The VT, inspiratory-to-expiratory ratio, and flow rate were adjusted to those which the subject found most comfortable.

All signals were recorded on a digital computer via an analog-to-digital interface (model 1401 Plus, Cambridge Electronic Design, Cambridge, UK). Digital signals were then analyzed by using commercially written software (Spike 2, Cambridge Electronic Design).

Protocol.   Subjects attended the laboratory on a day previous to the experiment to familiarize them with the equipment and assisted ventilation. The subjects were asked to restrict their sleep on the night before the experiment (3–5 h). No caffeine was taken for >4 h before the experiment and a light meal in the evening >2 h before the study. Subjects arrived in the laboratory at 10:00 PM. The EEG electrodes were attached, and the catheter for sampling end-tidal gases was placed just inside the nostrils. Subjects then lay on a bed and were asked to relax but not encouraged to sleep. The waking PCO2 was taken over 15–20 breaths at the end of a 10- to 20-min period of relaxed wakefulness. In the two subjects in whom a record of pHa was obtained, the radial artery catheter was then inserted under local anesthesia. The nasal mask was then attached and tested for air leaks. Following this, the assisted ventilation was started. Lights were turned off, and the subjects were asked to sleep. In four subjects, it was necessary to give 5 mg of zopiclone because they experienced no drowsiness, and no theta activity was recorded on the EEG. During the experiment, the presence of air leaks was monitored; if present, the magnitude was assessed from the ratio of expiratory VT (VTE) to inspiratory VT (VTI).

Analysis.   All variables were analyzed during the period of drowsiness before sleep onset. Drowsiness was defined as occurring during intermittent theta activity on the EEG in the absence of stable stage I/II sleep, as defined by the criteria of Rechtschaffen and Kales (27). Theta activity was defined as a low-frequency activity (<10 Hz) in the absence of rapid eye movements and alpha activity as a frequency content of >10 Hz. During drowsiness, the breaths where theta activity occurred for all or part of the breath were defined as theta breaths. All theta breaths were analyzed if they were preceded and followed by breaths associated with continuous alpha activity on the EEG (alpha breaths) (Fig. 1).



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 1. Protocol of 5-breath analysis based on the EEG. Central is the theta breath, with theta activity during part of that breath. This breath is preceded by 2 breaths (Pre-1 and Pre-2) and followed by 2 breaths (Post-1 and Post-2), with continuous alpha activity on the EEG, except there may be some carryover of theta activity into inspiration only of the Post-1 breath. Measurements are made of inspiratory and expiratory tidal volume (VT), inspiratory time (TI), breath duration, and expiratory pause (TP; time between cessation of expiratory gas flow and the next inspiration). TE, expiratory time. End-tidal PCO2 (PETCO2) and, in 2 subjects, arterial pH (pHa) and arterial PCO2 (PaCO2) are recorded at the end of expiration (see arrows).

 
VT, inspiratory time (TI), expiratory time (TE), TT, and expiratory pause (TP) were derived from the airflow signal; TP was the period between cessation of expiratory airflow and the onset of the next inspiration. The end-tidal gas tensions, together with pHa and PaCO2 in two subjects, were recorded at the end of expiration. The maximum pressure generated by the ventilator was recorded during inspiration. All measurements were made on the theta breaths and the two associated breaths before the theta breath (Pre-1 and Pre-2) and two post-alpha breaths (Post-1 and Post-2) (Fig. 1). The group mean (SD) for the five breaths was computed, and the differences between breaths were compared by using a one-way analysis of variance with repeated measures. When the analysis of variance was statistically significant, a post hoc analysis on the averages was carried out using Fisher’s least squares difference. P < 0.05 was considered statistically significant.


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The volume of air delivered by the ventilator and the ratio of TI to TT were set where the subject felt most comfortable. The resulting PETCO2 on ventilatory support was significantly less compared with no support [when fully vigilant at the start of the experiment, no ventilatory support = 39.5 Torr (SD 2.3); ventilatory support = 38.0 Torr (SD 1.4); P <0.05].

Figure 2 shows the effect of the transition from alpha to theta activity on ventilation. Two alpha breaths were followed by an apnea that was associated with theta activity. The record of mask pressure (Pmask) showed the ventilator was triggered by a drop of 1 cmH2O or less; this was the trigger sensitivity used in our experiments.



View larger version (25K):
[in this window]
[in a new window]
 
Fig. 2. Record of the effect of a change in EEG activity on breath duration due to a prolongation of TP. From above downwards, the record shows gas flow, EEG, VT, mask pressure (Pmask), and PETCO2.

 
The onset of theta activity usually occurred during TE (n = 33, 42%) or TP (n = 29, 37%), the remainder being during TI (n = 16, 21%). In the majority of the 78 theta breaths analyzed, alpha activity resumed during the inspiration of the first post-theta breath (Post-1); however, in 21 of the theta breaths (27%), alpha activity clearly resumed before the Post-1 breath, i.e., at the end of the theta breath.

Effect of changes in cerebral activity on components of TT, end-tidal gases, and maximum inspiratory pressure.   The mean results for all eight subjects are included in Table 1. The number of five-breath analyses periods for each subject varied (mean 9.8, range 9–11).


View this table:
[in this window]
[in a new window]
 
Table 1. Effect of changes in EEG activity on breath duration, end-tidal gases, and peak mask pressure

 
There was no statistically significant difference in the mean values for all variables between Pre-1 and Pre-2 breaths. Throughout the five breaths, mean (SD) TI remained constant [1.6 s (SD 0.5)], as did VTI [0.86 liter (SD 0.11)]. There was a significant prolongation of TE and TT during the theta breaths, with a 10-fold increase in TP. If an apnea is defined as a cessation of expiratory airflow for ≥10 s, there were 18 apneas, the mean per subject was 2.3 (SD 1.3); 23% of theta breaths analyzed. There was a small but nonsignificant decrease in TE, TT, and TP during the Post-1 breaths compared with Pre-2 and Post-2 breaths.

The mean PETCO2 of the theta and Post-1 breaths was significantly greater than that of the Pre-2 breaths but not the Post-2 breaths. However, the mean PETCO2 of the Post-2 breaths was not significantly different from that of the Pre-2 breaths. The mean PETO2 of the Post-1 and Post-2 breaths was lower than the preceding three breaths, and this was significant. Maximum pressure was significantly greater in the Post-1 breath only (Table 1).

Change in pHa, PaCO2 before, during, and after theta breaths.   Figure 3 shows a record of pHa, EEG, VT, and PETCO2 in one of the subjects. The theta breath is associated with a change in EEG activity and an apnea. The arrow in the figure joins the VT to the pHa oscillation, which it generates in the lung. The delay is due to the circulatory time between lung and radial artery. The relationship between pHa and VT is maintained throughout the record, such that the fourth pHa oscillation, during the apnea, is generated by the theta breath. Measurements of pHa were made at the end of expiration shown by the short vertical lines on the record of pHa. It can be seen that there is little change in pHa at the beginning of the apnea. Following the apnea, pHa rises rapidly, returning to the preapnea level by the third post-theta breath.



View larger version (25K):
[in this window]
[in a new window]
 
Fig. 3. Record of the effect of a theta breath and the associated alpha breaths on, from above downwards, pHa, EEG, VT, and PETCO2. Each respiratory oscillation in pHa is generated by the immediately preceding VT; thus the second oscillation is generated by the first VT shown on the extreme left of the figure (see arrow). The next breath has a slightly longer TP than the other alpha breaths so that the third oscillation has a slightly longer and deeper downstroke of the pHa. During the theta breath, the fall in pHa does not start until about halfway through TP. Thereafter, pHa rises following the return of ventilation to the same level that preceded the theta breaths. The vertical lines on the pHa record are the points at which pHa is measured, and PaCO2 is computed at the end of expiration. The values of pHa and PaCO2 from left to right are 7.415 and 41.0 Torr, 7.414 and 41.1 Torr, 7.413 and 41.3 Torr, 7.399 and 43.1 Torr, 7.393 and 44.2 Torr, 7.409 and 41.8 Torr, and 7.415 and 41.0 Torr, respectively. The PETCO2 of the Post-1 breath was less than the Pre-1 and the theta breath. However, the pHa fell from 7.399 to 7.393 and the PaCO2 rose from 43.1 to 44.2 between the end of the theta breath and the end of the Post-1 breath. The reason for this discrepancy between PETCO2 and PaO2 is addressed in the DISCUSSION under Changes in lung volume and arterial to end-tidal PCO2 gradient. Due to the slow respiratory frequency on assisted ventilation and the low or absence of gas flow at the end of expiration, the expired CO2 is contaminated from gas within the mask. The PETCO2 was derived from extrapolation of the expired CO2, as described in METHODS.

 
Twelve apneas were measured in the subject shown in Fig. 3, and the mean results are presented in a stylized format in Fig. 4 and in Table 2 (subject A), together with a statistical analysis in the table. At the start of the apnea, there was no change in pHa due to the delayed arrival in the radial artery of the oscillations in pHa from breaths Pre-2 and theta. Thus there was no significant difference in the means of pHa and PaCO2 between breaths Pre-1, Pre-2, and the first missed breath during the apnea (MB-1; Table 2). For the second missed breath during the apnea (MB-2), only four data points were collected, because the apneas terminated before the second missed breath. However, using unpaired t-tests, there was no significant difference between the mean pHa of MB-1 and MB-2, but there was a significant difference between the mean pHa of MB-2 and that of the theta breaths (P < 0.05). The mean pHa of the theta breaths were significantly lower, and the PaCO2 was higher than Pre-2 breaths, being 0.008 units and 1.1 Torr, respectively. After the theta breaths, pHa fell further. The mean pHa of Post-1 breaths was significantly lower and PaCO2 higher than all other means. However, by the Post-3 breaths, none of the means was significantly different from those of the Pre-2 breaths (the control value).



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 4. Cartoon illustrating the relationship between VT and pHa before, during, and after a theta breath. The pHa and PaCO2 values shown are the means of 12 theta breaths in subject A, shown previously in Fig. 3. Details of the means with statistical analysis are included in Table 2. Measurements shown by the arrows were made at the end of expiration, but, in addition, pHa and PaCO2 are recorded for "2 missed breaths," where it was anticipated in the theta breath that inspiration would have occurred based on the respiratory frequency of the preceding alpha breaths. The first oscillation in pHa (interrupted line) is generated by the VT of the Pre-1 breath not included in the figure. There is no significant change in pHa or PaCO2 until the end of the theta breath, when the mean pHa was significantly lower and the PaCO2 higher than the Pre-2 breaths. Thereafter, there was a further fall in pHa and rise in PaCO2, which was statistically significant. By the Post-2 breaths, pHa had begun to rise and PaCO2 fall, and by Post-3 breaths (not shown in figure) the pHa and PaCO2 were not significantly different from Pre-1 and Pre-2 breaths (Table 2).

 

View this table:
[in this window]
[in a new window]
 
Table 2. Arterial pH, PaCO2, PETCO2, and arterial-to-end-tidal PCO2 gradient for 2 subjects, before, during, and after a theta breath

 
The mean results for the second subject in whom pHa was recorded are shown in Table 2 (subject B). In this subject, there was no significant difference in the means of Pre-1, Pre-2, MB-1, MB-2, and Post-3. Thus by the third post-theta breath, the mean pHa and PaCO2 had returned to the Pre-2 breath values. The mean pHa of the theta breaths was significantly lower and the PaCO2 higher than all other breaths, except the Post-1 breaths.

The constant pHa up to and including the two missed breaths in both subjects confirm that the prolongation of TT of the theta breaths was solely due to a change in cerebral state, resulting in an altered ventilatory response to an unchanged chemical drive.

Relationship of pHa of the theta breath to that of the last preceding missed breath.   If the termination of an apnea is due to a chemical stimulus, then the pHa of the theta breath (end expiratory) should be less than the last missed breath. Therefore, a comparison of the mean pHa and PaCO2 of all theta breaths in both subjects with the corresponding means of the last missed breath was made. The mean pHa (SD) of the last missed breath was 7.392 (SD 0.023), and PaCO2 was 41.2 Torr (SD 1.4), whereas that of the theta breaths was 7.388 (SD 0.028) and 41.9 Torr (SD 1.3), respectively. Although the mean differences are small (0.004 units and 0.7 Torr), the mean pHa of the theta breaths was significantly lower and the mean PaCO2 higher than the mean of the last missed breaths.

PaCO2-to-PETCO2 gradient.   The breath-by-breath record of PaCO2 at the end of expiration allowed us to estimate this gradient (PaCO2-PETCO2 gradient) before, during, and after an apnea (Table 2). The gradient during the apnea could only be determined in one of the two subjects, because of contamination of expired CO2 by air from the mask during the long apneas in subject B. There was a small decrease in the mean gradient in the subject in whom it was determined (Table 2). In contrast, in both subjects, there was a significant increase in the PaCO2-PETCO2 gradient following a theta breath for Post-1 and Post-2 breaths. By the Post-3 breaths, the gradient was close to the control means, and the difference was not significantly different (Table 2). The cause of the widening of the PaCO2-PETCO2 gradient was that the rise in PaCO2 following an apnea was not reflected in the PETCO2.

Changes in expired volume and lung volume with changes in cerebral activity.   The mean VTE of the Post-1 breaths was significantly smaller than that of all other breaths (Table 1). Changes in VTE at constant VTI imply changes in lung volume. To determine from our data whether lung volumes had changed between breaths, we calculated the ratio of VTE/VTI for all subjects for all breaths. An increase in lung volume will lead to a reduction in VTE in relation to VTI and hence a decrease in the ratio. A decrease in lung volume will increase the ratio. The results are shown in Table 3. In seven of the eight subjects, the mean VTE/VTI of the Post-1 breaths for each of the seven subjects was less than that of all other breaths. The fall in the group mean from 0.85 to 0.74 was statistically significant. Although there was an increase in the VTE/VTI for the group mean of the theta breaths, this was not statistically significant, as an increase in the ratio occurred in only four subjects.


View this table:
[in this window]
[in a new window]
 
Table 3. Changes in lung volume from the ratio VTE/VTI

 
Air leaks around the mask.   VTE was usually less than VTI, pointing to an air leak between the exit of the Fleisch head and from around the nasal mask. To assess whether the leaks were constant, we used the ratio of VTE to VTI (Table 3). Because there was evidence of a change in lung volume in the Post-1 breath and possibly theta breath, Pre-1, Pre-2, and Post-2 breaths were used to assess the constancy or otherwise of air leaks through the five-breath analysis. It can be seen that the average ratio (0.85) for Pre-1, Pre-2, and Post-2 is constant, with the range of ratios between subjects varying from 0.70 to 1.0. In addition, the variation in the ratio within each subject for the three breaths ranged from 0.06 to 0.02, with a mean of 0.04 (SD 0.01).


    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Our results show that, at the transition between awake and asleep, there is an abrupt reduction of ventilatory motor output, which returned equally abruptly at the sleep-to-awake transition. These changes in ventilation at the alpha-to-theta transition were solely due to a change in cerebral state, independent of fluctuations in ventilation, PCO2, and upper airway resistance. Our observations were made when the subjects were drowsy, before the attainment of stable stage I/II sleep.

The continuous record of PaCO2 was in good agreement with that of the PETCO2 in showing that the PCO2 was stable before the apneas but also during the initial period of the theta breath, by which time the TT was already prolonged. Thus there was no change in PaCO2 at the alpha-theta transition. However, during the rapid changes in PCO2 at the theta-alpha transition, the PETCO2 failed to follow the PaCO2, although, by the third postapnea breath, the relationship between the two records was reestablished, and both records showed that the PCO2 had returned to preapnea levels.

Critique of the methods.   The mean PETCO2 of our subjects was 1.5 Torr below that when awake, and it is possible that the hypocapnia due to assisted ventilation was responsible for the apneas observed. It could be argued that, if the subjects had been breathing spontaneously, then the PCO2 would have been at or above that when awake, and no apneas would have occurred. This interpretation would be correct if our subjects had been in stable slow-wave sleep. Indeed, in this stable state, periodic breathing and apnea can only be produced by lowering the PCO2 with assisted or mechanical ventilation (8, 18). Our subjects were not studied in stable sleep but at the start of the onset of non-rapid eye movement (NREM) sleep. Periodic breathing and apneas have been consistently observed during the onset of NREM sleep (3, 25). Bulow (3) studied 75 normal subjects during drowsiness with alternating alpha and theta activity. He observed in more than two-thirds of subjects periods of periodic breathing, often lasting 10–20 min. Apneas were recorded lasting 10–20 s and occurred when the PETCO2 was "relatively low," closer to that when awake than when asleep. Zhou et al. (35) induced central apneas in stable sleep in men by mechanical ventilation via a nasal mask. A mean reduction in PETCO2 of 3.54 Torr produced apnea at a PETCO2 equivalent to that when awake. The observations of Bulow (3) and of Zhou et al. (35) are consistent with the PCO2 awake being at or very close to the CO2 apneic threshold when asleep [see also the review by Dempsey (7)]. Thus, during drowsiness at the beginning of sleep onset, apneas can be anticipated until the fall in ventilation associated with sleep raises the PCO2 above the apneic threshold.

We conclude that the apneas observed in our studies were not solely a consequence of assisted ventilation and the associated hypocapnia. The reason that apneas are seen both in spontaneously breathing subjects and those on assisted ventilation at sleep onset is that, initially, in both conditions, the PETCO2 is at or below the apneic threshold.

We have chosen to describe the prolongation of TP as an apnea, but it is appreciated that, if the definition of an apnea is cessation of breathing for ≥10 s, then we are including theta breaths, which do not meet that criterion. For the purpose of this paper, the term apnea has been retained.

Derivation of PaCO2.   This was determined continuously from the pHa record with the conversion of units to Torr using the hemoglobin and base excess of the blood samples (29). This relationship between pH and PCO2 in blood has been determined in vitro under steady-state conditions, and this might not hold under non-steady-state conditions in vivo as in our experiments, where there is a continuously changing pH. However, it has been observed that breath-by-breath oscillations in PaCO2 in humans and animals, derived from continuous records of pHa (non-steady state), are equivalent to the corresponding fluctuations in alveolar PCO2 (1, 2). On this basis, it is considered that the derivation of PaCO2 is justified.

Air leaks around the mask.   From the analysis of VTE/VTI (Table 3), we conclude that air leaks were virtually constant between breaths and that the reported changes in VTE are valid. This constancy of air leak between breaths is important when considering the mean fall in VTE of the Post-1 breaths, where there was an increase in mean Pmask and hence a likelihood of a larger mask leak. This is unlikely in that three of the subjects had a rise in Pmask of <0.4 cmH2O, but there was still a fall in VTE. In addition, in subject 2 with the highest increase in Pmask (7 cmH2O), there was a small increase, not reduction, in VTE/VTI.

Constancy of PCO2.   The design of the experiments was to hold PaCO2 constant so as to observe the effect on TT of a change in cerebral state. To determine whether the purpose of this design was achieved, it is necessary to show that the rise in PCO2 following an apnea had been eliminated by the control breaths of the next theta breath. As regards the peripheral chemoreceptors, we consider that the purpose of the design was achieved. There was no significant difference between mean PETCO2 of the Pre-1 and Pre-2 breaths and that of the Post-2 breaths for the group data. There was no significant difference in PaCO2 between Pre-1 and Pre-2 breaths and Post-3 breaths in the two subjects in whom it was measured. Where there were only four breaths between theta breaths, then Post-3 breaths would be the equivalent of the Pre-1 breaths of the following group of five breaths. Assuming the PaCO2 in the radial artery approximates that at the peripheral chemoreceptors in time (14), then the constancy of the PaCO2 at the chemoreceptors between theta breaths can reasonably be claimed to have been achieved.

While PCO2 at the peripheral chemoreceptors can be determined from that in the arterial blood, this is not so for the chemosensitive regions of the brain due to the relatively slow "wash-in" and "washout" of CO2 in these areas. The magnitude of the ventilatory effects of a brief disturbance of PaCO2 at the central chemoreceptor sites has been estimated in experiments in which it was assumed that the contribution to ventilation from the peripheral chemoreceptors was small or nonexistent (13, 16, 19). For example, the ventilatory effect of the inhalation of CO2 for one to three breaths has been studied in humans in hyperoxia awake and asleep (20). The overall cumulative ventilatory response was the same between awake and asleep but was slower in sleep. The peak ventilatory response was 0.08 l/min for a 1-Torr rise in PETCO2; ~2% of expired volume (VE) was due to a change in VT alone. Special techniques not used in the present experiments are required to detect such small changes in ventilation, such as pseudorandom binary CO2 stimulation (16, 19, 31) and ensemble averaging (13). The ventilatory response to a transient disturbance of CO2 has been shown to be predominantly, if not exclusively, due to the peripheral chemoreceptors in the dog (23) and human (31). We conclude that the effect of the transient rise in PCO2 in our experiments on the theta and subsequent breaths was mediated through the peripheral chemoreceptors, and for the reasons quoted above this effect was cleared by the control breaths of the next theta breath.

Nonchemical inhibition of respiratory motor output.   There is evidence of a significant nonchemical inhibition of respiratory motor output in sleeping humans on assist-control mechanical ventilation (33). This inhibition may well have been present in our experiments during the theta breaths. The effect on TT will have been to prolong TE by 13–32% at the transition from alpha to theta activity on the EEG. However, this nonchemical inhibition of respiratory motor output has been shown not to lead to apnea (33) and, therefore, will have had little or no effect on TP.

Changes in lung volume and PaCO2-to-PETCO2 gradient.   These two changes have been discussed together because they may be causally related. Although we did not measure changes in lung volume directly, there was evidence of rapid changes in lung volume over one to two breaths between the alpha and theta transitions. Hudgel and Devadatta (12) measured lung volumes directly during wakefulness and stable sleep and observed a reduction in lung volume in sleep in 8 out of 10 subjects. They attributed this reduction to a loss of muscle tone of the chest muscles with an implied reduction in thoracic cage recoil, leaving the reduction in lung volume due to the unchanged elastic recoil of the lung. Our mean results for the changes in VTE at the alpha and theta transitions are compatible with these proposed mechanisms (Table 1). The reduction in VTE at the theta-to-alpha transition was consistent in all but one of the subjects and was presumably due to a return of muscle tone with the onset of alpha activity and hence increase of lung volume. However, the expected increase of VTE at the alpha-to-theta transition occurred in only four of the eight subjects. We cannot explain this discrepancy between the two transitions, but it could be due to the effect of continuous gas absorption during the theta breaths, an effect on VE opposite to that of the elastic recoil of the lung, thereby leaving VTE unchanged in some subjects.

In both subjects in whom we had a record of PaCO2, there was a significant increase in the PaCO2-PETCO2 gradient in the two breaths following an apnea (Post-1 and Post-2, Table 2). During the apnea in one of the subjects (subject A), there was a small decrease in gradient, but, in subject B, it was not possible to determine the gradient because of contamination of expired air from within the mask. The cause of the increase in gradient after the apneas was not determined, and further research would be required to elicit the cause. However, the changes in lung volume and PaCO2-PETCO2 gradient were closely related in time, and we put forward the possibility that the contraction of lung volume during the apnea altered the ventilation-to-perfusion relationships within the lung, possibly due to airway closure. The effect of this change was manifest in the widening of the PaCO2-PETCO2 gradient in the two subsequent breaths. In addition, the mean Pmask of the first breath after the theta breath (Post-1; Table 1) was significantly higher than all other breaths for the same VTI, indicating a change in compliance and/or resistance during the theta breath.

The contribution of the continuous measurement of pHa to the present experimental results.   PETCO2, a discontinuous measurement, has been used in the present experiment as an indicator of changes in PaCO2. However, this correlate may fail when there is absent expiratory airflow, rapid changes in PaCO2, and a changing arterial-to-alveolar PaCO2 gradient. All of these factors operated in our experiments and hence the need for a continuous measurement of PaCO2 to determine when PETCO2 could be used as a surrogate for PaCO2. While the changes in PaCO2 during and after an apnea could be qualitatively predicted, on known physiological mechanism, we are unaware that they have ever been recorded before.

There was good agreement between record of PaCO2 with that of PETCO2, in showing PaCO2 was stable before the theta breaths. More importantly, PaCO2 remained unchanged during TI, TE, and the beginning of the apnea (TP) of the theta breaths. This is illustrated in Fig. 4 (missed breaths) and confirmed for both subjects in Table 2. Following an apnea, PETCO2 cannot be used to determine changes in PaCO2 because of a changing PaCO2-to-PETCO2 gradient. However, by the third post-theta breath, the relationship between PETCO2 and PaCO2 had been reestablished, and the rise in PaCO2 following the theta breath had been eliminated, thus preserving the constancy of PaCO2 before the next theta breath.

Resumption of breathing after a theta breath.   We consider that the resumption of breathing after a theta breath was due to the increase of PCO2, either acting directly as a chemical stimulus to ventilation or indirectly by causing a change in cerebral state. In two-thirds of the theta breaths in this study, the onset of the Post-1 breaths preceded a return in the EEG to an alpha rhythm. In addition, in all subjects, there were occasions when theta activity persisted beyond the Post-1 breaths, leading to irregular breathing with repetitive apneas (S. J. G. Semple, unpublished results). Following the theta breaths, the mean TT and TP (Table 2) were shorter than all other breaths, although this failed to reach statistical significance. If these mean changes are real, then it can be reasonably attributed to the rise in PaCO2, which continued after the theta breath (Table 2) and, for reasons given earlier, would be mediated primarily through the peripheral chemoreceptors.

Comparison of present findings with previous research.   The results of the present experiments confirm the results of Trinder et al. (32) that ventilatory instability at sleep onset is predominantly due to alternating changes in cerebral state, with state being determined, as in our experiments, by alpha and theta activity in the EEG record. There were important differences in experimental design between the two studies. In the experiments of Trinder et al. (32), their subjects were breathing spontaneously, and there were 2–10 breaths or more between the changes in cerebral state. They observed a cyclical variation in VE and PETCO2, predominantly determined by hyperventilation, and fall in PETCO2 associated with the arousal at the theta-to-alpha transition, an increase in ventilation in excess of metabolic demand (11). The magnitude and duration of the cyclical variation in ventilation depended on a complex interaction between changes in state, their duration, and a fluctuating PCO2. In contrast, in our own experiments, VE and PCO2 were stable at the alpha-to-theta transition, whereas, at arousal, there was a transitory rise in PCO2, with no statistically significant change in TT. Thus, in our experiments, at the alpha-to-theta transition, TT was determined by a change in cerebral state alone.

The most important difference in experimental design and results between our experiments and those of Trinder et al. (32) was that, in their experiments, the subjects at sleep onset were breathing spontaneously, and no apneas were recorded. In the latter experiments, the anticipated rise in PCO2 between awake and asleep had occurred (a mean rise of 2.0 Torr), and it is reasonable to assume that the PCO2 was now above the apneic threshold and hence the absence of apneas. In contrast, in our own experiments, a rise in PCO2 at sleep onset was prevented by the assisted ventilation at fixed volume. These findings illustrate the critical role of the rise in PCO2 at sleep onset in reducing or eliminating central apneas and thereby lessening ventilatory instability (28).

Apneic threshold and control of ventilation.   We did not determine the apneic CO2 threshold directly in our subjects, but plainly the PCO2 at the alpha-theta transition was below the threshold. The mean PCO2 before the transition was 1.5 Torr (SD 1.8) below that awake, with a range of 0–3.6 Torr. In four subjects, the difference was under 2.0 Torr, being 0, 0.4, 0.5, and 1.8 Torr. If, at the start of drowsiness, the PCO2 is that when awake, then a small, augmented VT from whatever cause may precipitate an apnea (8).

The PETCO2 in the present experiments was below the eupneic PCO2 when awake, and this will have the effect of decreasing plant gain (the ventilatory increase required for a given reduction in PCO2) (9). This will render ventilation more stable during the theta breaths than in subjects breathing spontaneously at a higher PCO2. However, the PETCO2 of our subjects was already below the apnea/hypopnea CO2 threshold, and we predict that this would also have been the case if the subjects were breathing spontaneously because, for reasons given earlier, the PCO2 when awake is at or below the threshold. We do not think, therefore, that the predicted increase in ventilatory stability in our experiments affected the conclusions drawn from the results.

Conclusions amd predictions.   At the start of sleep onset when drowsy, and when the PCO2 is close to or at that when awake, apneas and periodic breathing will occur with changes in cerebral state. The most ready explanation for these findings is that the PCO2 when awake is close to or at the apneic/hypopneic CO2 threshold (7). This explanation would also account for the apneas and periodic breathing frequently observed at the onset of NREM sleep and is likely to be most prevalent in subjects with a high-ventilatory sensitivity to CO2 (10). We have also shown that, when the sleep-related rise in PCO2 is prevented, this unstable ventilatory condition is likely to persist. Any condition that maintains PCO2 close to or below that when awake will render ventilation unstable. This phenomenon has been well documented in patients with heart failure and central sleep apnea where there is no rise in PCO2 from wakefulness to sleep and no proportional reduction in their apnea or hypopnea threshold. In contrast, patients with heart failure and no central sleep apnea do show a rise in PCO2 in stable sleep (34).


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The Wellcome Trust, The Royal Society, and The Breathlessness Trust supported this study.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank Abe Guz for continuing support, advice, and encouragement.


    FOOTNOTES
 

Address for reprint requests and other correspondence: S.J.G. Semple, Dept. of Respiratory Medicine, Charing Cross Hospital, Fulham Palace Rd., London W6 8RF, UK (e-mail: s.semple{at}imperial.ac.uk)

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.


    REFERENCES
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 

  1. Band DM, Cameron IR, and Semple SJG. The effect on respiration of abrupt changes in carotid artery pH and PCO2 in the cat. J Physiol 211: 4579–4594, 1970.
  2. Black AMS and Torrance RW. Respiratory oscillations in chemoreceptor discharge in control of breathing. Respir Physiol 13: 221–237, 1971.
  3. Bulow K. Respiration and wakefulness in man. Acta Physiol Scand 59, Suppl 209: 46–49, 1963.
  4. Chapman KR, Bruce EN, Gothe B, and Cherniack NS. Possible mechanisms of periodic breathing during sleep. J Appl Physiol 64: 1000–1008, 1988.
  5. Cordingley JJ, Palazzo MJA, Thomson S, Hynd J, Cross BA, and Semple SJG. A fast responding intra-arterial pH electrode for use in the peripheral artery of adult humans and large mammals: a technical development for use in research. J Med Eng Technol 22: 233–240, 1998.
  6. Cross BA, Stidwill RP, Hughes KR, Peppin R, and Semple SJG. A comparative study of aortic pH oscillations in conscious humans and anaesthetised cats and rabbits. Respir Physiol 102: 51–62, 1995.
  7. Dempsey JA. Crossing the apnoeic threshold: causes and consequences. Exp Physiol 89: 192–199, 2004.
  8. Dempsey JA and Skatrud JB. Sleep-induced apneic threshold and its consequence. Am Rev Respir Dis 133: 1163–1170, 1988.
  9. Dempsey JA, Smith CA, Przybylowski T, Chenuel B, Xie A, Nakayama H, and Skatrud JB. The ventilatory responsiveness to CO2 below eupnoea as a determinant of ventilatory stability in sleep. J Physiol 560: 1–11, 2004.
  10. Dunai J, Kleiman J, and Trinder J. Ventilatory instability during sleep onset in individuals with high peripheral chemosensitivity. J Appl Physiol 87: 661–672, 1999.
  11. Horner RL, Rivera MP, Kozar LF, and Phillipson EA. The ventilatory response to arousal from sleep is not fully explained by differences in CO2 levels between sleep and wakefulness. J Physiol 534: 881–890, 2001.
  12. Hudgel DW and Devadatta P. Decrease in functional residual capacity during sleep in normal humans. J Appl Physiol 57: 1219–1322, 1984.
  13. Jacobi MS, Paul CP, and Saunders KB. The transient response to carbon dioxide at rest and in exercise in man. Respir Physiol 77: 225–238, 1989.
  14. Jain SK, Subramanian S, Julka DB, and Guz A. Search for evidence of lung chemoreflexes in man: study of respiratory and circulatory effects of phenyldiguanide and lobeline. Clin Sci 42: 163–177, 1972.
  15. Kay A, Trinder J, Bowes G, and Kim Y. Changes in airway resistance during sleep onset. J Appl Physiol 76: 1600–1607, 1994.
  16. Lai J and Bruce EN. Ventilatory stability to transient CO2 disturbances in hyperoxia and normoxia in awake humans. J Appl Physiol 83: 466–476, 1997.
  17. Longobardo GS, Gothe B, Goldman MD, and Cherniack NS. Sleep apnoea considered as a control system instability. Respir Physiol 50: 311–333, 1982.
  18. Meza S, Mendez M, Ostroski M, and Yonnes M. Susceptibility to periodic breathing with assisted ventilation during sleep in normal subjects. J Appl Physiol 5: 1929–1940, 1998.
  19. Modarreszadeh M and Bruce EN. Long lasting ventilatory response of humans to a single breath of hypercapnia in hyperoxia. J Appl Physiol 72: 242–250, 1992.
  20. Modarreszadeh M, Bruce EN, Hamilton H, and Hudgel DW. Ventilatory stability to CO2 disturbances in wakefulness and quiet sleep. J Appl Physiol 79: 1071–1081, 1995.
  21. Morrell MJ, Shea SA, Adams L, and Guz A. Effect of inspiratory positive airway pressure upon breathing in man during wakefulness and sleep. Respir Physiol 93: 57–70, 1993.
  22. Morrell MJ, Harty HR, Adams L, and Guz A. Breathing during wakefulness and NREM sleep in humans without an upper airway. J Appl Physiol 81: 274–281, 1996.
  23. Nakayama H, Smith CA, Rodman JR, Skatrud JB, and Dempsey JA. Carotid body denervation eliminates apnea in response to transient hypocapnia. J Appl Physiol 94: 155–164, 2003.
  24. Pack AI, Cola MF, Goldszmidt A, Ogilvie MD, and Gottschalk A. Correlation between oscillations in ventilation and frequency content of the encephalogram. J Appl Physiol 72: 985–992, 1992.
  25. Phillipson EA and Bowes G. Control of breathing during sleep. In: Handbook of Physiology. The Respiratory System. Control of Breathing. Bethesda, MD: Am. Physiol. Soc., 1986, sect. 3, vol. II, pt. 2, chapt. 19, p. 649–690.
  26. Rechtschaffen A and Kales KA. Manual of Standardized Terminology, Techniques and Scoring System for Sleep Stages of Human Subjects. Bethesda, MD: NIH, 1968. (Publ no. 204)
  27. Semple SJG, Cordingley JJ, Thomson S, and Morrell MJ. Prevention of the rise in PCO2 at sleep onset may cause ventilatory instability (Abstract). Am J Respir Crit Care Med 159: A461, 1999.
  28. Severinghaus JW. Blood gas calculator. J Appl Physiol 21: 1108–1116, 1966.
  29. Skatrud JB and Dempsey JA. Interaction of sleep state and chemical stimuli in sustaining rhythmic ventilation. J Appl Physiol 55: 813–822, 1983.
  30. Sohrab S and Yamashiro SM. Pseudorandom testing of ventilatory response to inspired carbon dioxide in man. J Appl Physiol 46: 1000–1009, 1980.
  31. Trinder J, Whitworth F, Kay A, and Wilkin P. Respiratory instability during sleep onset. J Appl Physiol 73: 2462–2469, 1992.
  32. Wilson CR, Satoh M, Skatrud JB, and Dempsey JA. Non-chemical inhibition of respiratory motor output during mechanical ventilation in sleeping humans. J Physiol 518: 605–618, 1999.
  33. Xie A, Skatrud JB, Pules DS, Rahko PS, and Dempsey JA. Apnoea-hypopnea threshold for CO2 in patients with congestive heart failure. Am J Respir Crit Care Med 165: 1245–1250, 2002.
  34. Zhou XS, Shahabuddin S, Zahn BR, Babcock MA, and Badr MS. Effect of gender on the development of hypocapnic apnea/hypopnea during NREM sleep. J Appl Physiol 89: 192–199, 2000.



This article has been cited by other articles:


Home page
ChestHome page
L. H. Ibrahim, S. R. Patel, M. Modarres, N. L. Johnson, R. Mehra, H. L. Kirchner, and S. Redline
A Measure of Ventilatory Variability at Wake-Sleep Transition Predicts Sleep Apnea Severity
Chest, July 1, 2008; 134(1): 73 - 78.
[Abstract] [Full Text] [PDF]


Home page
ThoraxHome page
Y.-L. Lo, A. S Jordan, A. Malhotra, A. Wellman, R. A Heinzer, M. Eikermann, K. Schory, L. Dover, and D. P White
Influence of wakefulness on pharyngeal airway muscle activity
Thorax, September 1, 2007; 62(9): 799 - 805.
[Abstract] [Full Text] [PDF]


Home page
ChestHome page
K. R. Casey, K. O. Cantillo, and L. K. Brown
Sleep-Related Hypoventilation/Hypoxemic Syndromes
Chest, June 1, 2007; 131(6): 1936 - 1948.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
99/5/2036    most recent
01040.2004v1
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 ISI Web of Science
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 ISI Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Thomson, S.
Right arrow Articles by Semple, S. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Thomson, S.
Right arrow Articles by Semple, S. J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE