The current study was undertaken to attempt to replicate our earlier studies by using 5% and 7% CO2 inhalation and hyperventilation with improved methodology and to thoroughly examine the respiratory response of our largest series of panic disorder patients to date. Specifically, we aimed to determine which one of the two most frequently employed CO2 concentrations should be favored in assessing the respiratory response of patients with panic disorder; to characterize this response by applying novel analytic techniques to respiratory measures identified in our earlier work as discriminating patients and normal comparison subjects; to identify respiratory predictors of panic; and to describe the respiratory pattern accompanying CO2-induced panic. We were particularly interested in the clinical and treatment implications of respiratory abnormalities. We also wanted to address previous methodological shortcomings such as possible order effects of the challenges, absence of self-ratings of panic attacks, possible selection and assessment bias of patients and comparison subjects, and relatively small group size.
Fifty-nine patients, 21 men and 38 women between the ages of 19 and 58 years (mean=34.0, SD=9.8), with DSM-III-R panic disorder with (N=38) or without agoraphobia were evaluated with standard psychiatric interview. The diagnoses were confirmed by the Structured Clinical Interview for DSM-III-R (SCID) (
+17). Patients were excluded if they met criteria for current major depressive disorder, obsessive-compulsive disorder, or substance use disorder or if they had a lifetime history of schizophrenia or bipolar disorder. Patients with an anxiety disorder other than panic disorder and those with dysthymia were included if the disorder was secondary to panic disorder.
The comparison group consisted of 20 men and 19 women between the ages of 20 and 60 years (mean=32.1, SD=9.9). They were assessed by either the Schedule for Affective Disorders and Schizophrenia—Lifetime Version Modified for the Study of Anxiety Disorders (
+18) or the SCID. They were free of any lifetime history of anxiety disorders, major affective disorders, schizophrenia, and current substance use disorders. All subjects underwent medical evaluations and were in good health.
Subjects were required to have been drug free for at least 2 weeks before baseline assessment. However, the use of benzodiazepines up to the equivalent of 5 mg of diazepam/day was permitted until 48 hours before testing. Only one patient took benzodiazepines within 2 weeks of testing. All subjects signed informed consent forms that explained that they would breathe air mixed with CO2 and also would be asked to hyperventilate room air, that the procedures were not dangerous, and that anxiety or panic could occur during the study. The study was approved by the institutional review board.
Subjects were tested on two mornings separated by not more than 10 days. The first morning involved CO
2 response testing that used the Read rebreathing method (
+19), a traditional pulmonary assessment of ventilatory sensitivity. These results have been reported elsewhere (
+20). On the second morning, following an overnight fast, subjects underwent the "canopy" procedure reported here.
First, subjects were randomly assigned to watch one of two videotapes that gave different sets of instructions regarding the effects of CO
2 inhalation. The data showing that the instructions did not meaningfully influence the responses to respiratory challenges (including no interactions) have been reported separately (
+21). Next, the subject was asked to lie down, and his or her head was placed in a clear plastic canopy described in detail elsewhere (
+16). The canopy was sealed, but the subject could see and hear and be seen and heard at all times. The subject knew how to open the canopy quickly by flipping a latch (no subject actually did this). Subjects could also signal to have the procedure terminated at any time. The canopy was vented by a source of external air at 40 liters per minute. End-tidal CO
2 concentration was measured by capnograph. End-tidal CO
2 correlates well with arterial CO
2 concentration (
+22). Respiratory rate and tidal volume were recorded in a spirometer.
The experiment consisted of seven periods: 1) room air breathing for 20 minutes, 2) according to random assignment either 5% CO2 in room air for 20 minutes or room air hyperventilation for 15 minutes (for the hyperventilation, a metronome was placed on top of the canopy with a flashing light, and the subject was told to take a deep breath every time the light flashed [30 breaths per minute]), 3) room air breathing for 15 minutes, 4) the intervention that was not assigned in period 2, 5) room air breathing for 15 minutes, 6) 7% CO2 in room air for 20 minutes, and 7) room air breathing for 15 minutes. Subjects were not informed about the timing and type of each period.
Ten raters were trained to identify panic attacks by watching videotapes of lactate infusions. The occurrence of a panic attack was determined in two ways. The rater, blind to diagnosis, used DSM-III-R criteria to assess the presence of a panic attack (four or more DSM-III-R symptoms plus a clear and sudden increase in anxiety). A few months after the study began a second determination of panic was added. Here the subject was asked to indicate whether he or she had experienced a panic attack by pointing to answers on a card. The rater was not permitted to observe the subject's rating. Because the subject's rating of panic was included after the study began, fewer subject assessments, compared to rater assessments, were available. In addition to the varying extent to which subjects participated in the experiment, this is the other reason why the number varies in certain analyses.
Four rating scales were administered at the end of each test period to determine anxiety and distress: the 27-item Acute Panic Inventory (
+23), 10-point Likert-scored Anxiety and Apprehension Scales, and the Borg Scale of Exertion (
+24), which measures effort required to breathe. Whenever a subject could not tolerate the intervention (indicated by raising his or her hand), room air breathing resumed, and the scales were administered.
Panic rate. Rates of panic among the three interventions, between patients and normal volunteers, and also separately for raters' and subjects' ratings and by gender were compared with chi-square (or Fisher's exact) tests.
Physiology. The respiratory physiology data consisted of continuous (every other breath) recordings of end-tidal CO2, respiratory rate, and tidal volume. Before statistical analysis aberrant values were edited according to preset criteria (i.e., end-tidal CO2 less than 10 torr during room air, less than 33 torr during 5% CO2, less than 48 torr during 7% CO2; respiratory rate greater than 80; tidal volume less than 100 ml). Data points with three standard deviations outside the mean were also eliminated. These corrections resulted in the elimination of less than 3% of the data. Vital capacity was covaried in comparisons involving tidal volume.
In order to correct for multiple dependent variables, whenever appropriate, multivariate analyses were performed first with subject groups, respiratory measures (end-tidal CO
2, tidal volume, respiratory rate), and gender included as factors. Significant differences were followed up by separate comparisons. However, given our a priori directional hypotheses (
+25), the reported prominent gender differences in respiratory response (
+26), and our unequal sex distribution, all analyses were repeated in men and women separately whenever group sizes were meaningful (greater than five), regardless of the results of the multivariate analyses. These "exploratory" analyses are listed separately within each section.
Five-minute mean values (and standard deviations) were calculated for each variable throughout the study. In the first set of analyses the last 5-minute mean of each period was used as a single value. Following up on significant multivariate tests, first we compared patients with normal subjects within each period for each variable by using analysis of variance (ANOVA). These analyses assessed group differences at baseline and during the interventions regardless of panic rates. The three baselines (room air breathing) were also compared with ANOVAs for each variable in order to determine if the interventions caused any differences in "recovery."
Second, we divided the patients into those panicking in response to 5% CO2 and those who did not (also into those who did and did not panic in response to 7% CO2, to hyperventilation, and to both 5% and 7% CO2), separately according to the raters' and the subjects' ratings, and performed three group ANOVAs for each period and each variable.
During CO2 inhalation and hyperventilation it takes approximately 10 minutes to reach steady state respiratory conditions in the canopy. Therefore, the last 5-minute means of these interventions were used as single values. This method of analysis is adequate for subjects who were able to complete the challenges. Subjects who panicked during any of the interventions and stopped the procedure early were exposed to the challenge for a shorter time period than nonpanicking subjects. For these subjects we considered the 5-minute mean immediately before the panic attack. Subjects whose panic attack occurred sooner than 5 minutes into the period (eight for 5% CO2, four for hyperventilation, and 19 for 7% CO2) were not included in this analysis.
Since the overall respiratory response to CO2 is hyperbolic, we also performed a time analysis to examine the first few minutes of linear respiratory response to CO2. This ANOVA included all subjects. One-minute means of respiratory rate, tidal volume, minute ventilation, and end-tidal CO2 were calculated at the following six time points: 2 minutes before CO2 and 1, 2, 3, 4, and 5 minutes of CO2 inhalation. Many patients who panicked had their panic attack by minute 5, so further comparison was meaningless because of diminishing group size. CO2 sensitivity (delta minute ventilation divided by delta end-tidal CO2) was calculated on a minute-by-minute basis by using the difference between the values at each time point and baseline. ANOVAs were used to compare the same time periods across the three groups (panickers, nonpanickers, and comparison subjects). We also performed this comparison by using respiratory rate (delta respiratory rate divided by delta end-tidal CO2) and tidal volume (delta tidal volume divided by delta end-tidal CO2) and self- and raters' ratings separately.
In order to compare the pattern of response in the three subject groups exposed to CO2 for the full 15 minutes, we also plotted tidal volume, respiratory rate, minute ventilation, and end-tidal CO2 curves by using 1-minute means throughout the CO2 challenges for comparison subjects, nonpanicking patients, and self-rated panickers and performed ANOVAs followed up by pairwise comparisons when appropriate.
The recovery analyses compared the pattern of change in tidal volume, respiratory rate, minute ventilation, and end-tidal CO2 from the end of the interventions through the subsequent 10 minutes of room air breathing by using ANOVAs on 1-minute means. The recovery analyses were repeated by using self- and raters' ratings separately.
In order to identify predictors of panic response, logistic regression analyses were performed with the forward stepping procedure by using eight preselected baseline measures (respiratory rate, tidal volume, minute ventilation, end-tidal CO2, and their standard deviations) with self-rated panic as the outcome measure and sex and age as control measures (probability in and out criteria were set at 0.20 and 0.25, respectively).
All p values reported are two-tailed, with significance levels of 0.05. When indicated by unequal variances, the separate variance estimate statistic was used.
+Table 1 shows the panic rates during the interventions according to the rater and according to the subject. Significantly more patients panicked in response to 7% than 5% CO
2, and the panic rates in response to CO
2 were significantly higher than those in response to hyperventilation, determined by using either ratings. The order of the interventions did not affect the rates of panic. During 5% CO
2 inhalation, of the 51 patients for whom both ratings were available, 31 panic patients rated themselves as panicking, but only 14 were noted by the clinician. There was only one case in which the clinician-rated panic attack was not endorsed by the subject. During 7% CO
2 inhalation, of the 47 patients for whom both ratings were available, 34 panic patients rated themselves as panicking, and 26 were so rated by clinicians. Four clinician-rated panic attacks were not endorsed by the patient. Agreement between blind and self-ratings of panic did not reach statistically significant differences. The only sex difference was a significantly higher self-rating of panic in female patients (N=26 of 31) than in male patients (N=8 of 16) during 7% CO
2 inhalation (χ
2=6.05, df=1, p<0.01).
The behavioral rating scales demonstrated more anxiety, apprehension, and breathlessness in patients than in comparison subjects. Compared to our preliminary report (
+16), these differences became more accentuated.
Baseline. The physiology data (means and standard deviations) are summarized in
+table 2. The physiological differences were most prominent between panicking and nonpanicking patients when the self-ratings to 5% CO
2 inhalation were used. Unless stated otherwise, this rating is used in all subsequent analyses and tables.
At first baseline the multivariate test revealed significant sex effects only (F=3.84, df=3,79, p<0.01). Follow-up exploratory comparisons showed that baseline end-tidal CO2 for patients was significantly lower than that for the comparison group (F=7.64, df=1,91, p<0.007). Both panicking and nonpanicking patients had significantly lower end-tidal CO2 values than normal subjects (F=3.12, df=2,81, p<0.05) but did not differ from each other. Respiratory rate, tidal volume, and minute ventilation were all higher in patients than in comparison subjects, but the differences did not reach significance. The overall baseline comparison of standard deviations again revealed significant sex effects (F=4.53, df=1,90, p<0.04). The subsequent comparisons within gender showed significantly higher individual standard deviations in tidal volume (t=2.25, df=51, p<0.03, separate variance estimate) and minute ventilation (t=3.19, df=46, p<0.003) in female patients than in female comparison subjects.
We compared the three baselines (periods 1, 3, and 5) in order to assess possible differential order effects of CO2 or hyperventilation in periods 2 and 4. The overall multivariate test showed no significant order effects.
5% CO2 inhalation. Both nonpanickers (mean=16.65 minutes, SD=5.33) (t=2.05, df=50, p<0.05) and comparison subjects (mean=19.06 minutes, SD=2.10) (t=4.83, df=36,13, p<0.001, separate variance estimate) continued to breathe 5% CO2 significantly longer than panickers (mean=12.95 minutes, SD=6.89) (F=10.36, df=2,86, p<0.001). The multivariate test with three groups, three respiratory measures, and sex, which used pre- and intervention values, showed overall sex effect (F=3.19, df=3,69, p<0.03), time effect (F=101.1, df=3,69, p<0.0001), and sex-by-time interaction (F=3.64, df=3,69, p<0.02). The follow-up ANOVAs showed group effect for end-tidal CO2 (F=4.14, df=2,69, p<0.02). Both panicking and nonpanicking patients had lower end-tidal CO2 than comparison subjects. The group-by-sex interaction for end-tidal CO2 (F=8.98, df=1,81, p<0.004) was due to significantly lower values in both panicking and nonpanicking patients than in comparison subjects; female patients had the lowest values followed by male comparison subjects, male patients, and female comparison subjects.
Through use of the last 5-minute means the multivariate test showed a significant sex effect (F=2.85, df=3,75, p<0.04), and the follow-up exploratory ANOVAs found that the differences were in respiratory rate (patients: mean=20.37, SD=5.09; comparison subjects: mean=18.11, SD=3.92) (F=4.54, df=1,83, p<0.04) and end-tidal CO2 (patients: mean=49.41, SD=3.78; comparison subjects: mean=51.87, SD=3.03) (F=10.66, df=1,81, p<0.002). The group-by-time near-significant interaction for respiratory rate (F=2.65, df=2,43, p<0.08, with Greenhouse-Geisser correction) indicated a larger increase in panicking female patients than in the other two female groups.
The time analysis elaborated on this finding. The multivariate comparison showed significant sex effect (F=3.63, df=3,61, p<0.02), time effect (F=46.16, df=15,49, p<0.001), and significant sex-by-time (F=3.94, df=15,49, p<0.001) and group-by-time (F=1.65, df=30,96, p<0.04) interactions. The follow-up comparison showed a near-significant group difference and significant group-by-time and group-by-time-by-sex interactions in respiratory rate (
+figure 1). Patients who panicked in response to 5% CO
2 increased their respiratory rate much sooner than comparison subjects or nonpanickers. The interactions were due to a more variable course for panickers than for nonpanickers and comparison subjects and a more accentuated increase in women than in men.
Using 1-minute means we performed the multivariate test comparing the three groups (comparison subjects, self-rated panickers, nonpanicking patients) exposed to CO
2 for the full 15 minutes. The significant overall sex effect (F=6.18, df=3,44, p<0.001) and group-by-sex interaction (F=2.93, df=6,86, p<0.01) were due to a steeper, smoother hyperbolic tidal volume rise in comparison subjects than nonpanicking patients or self-rated panickers (F=1.80, df=7,180, p<0.09, with Greenhouse-Geisser correction). Respiratory rate showed an almost opposite pattern; here the significant group effects in women (F=3.85, df=2,26, p<0.03) was the result of a steep and chaotic increase in patients, with relatively smaller changes in comparison subjects. The group-by-sex interaction for end-tidal CO
2 (
+figure 2) was due to significant group differences in women (F=6.26, df=2,27, p<0.006); comparison and nonpanicking women developed steady levels by minute 5, while the end-tidal CO
2 in self-rated panickers continued to rise. The within-subject correlation coefficients between minute ventilation and end-tidal CO
2 confirmed this pattern; the mean correlation in self-rated panickers (r=0.75) was significantly higher than that in either comparison subjects (r=0.59) (t=2.13, df=43, p<0.04;) or nonpanicking patients (r=0.41) (t=2.70, df=28, p<0.02).
Hyperventilation. Panicking patients hyperventilated in room air for a mean of 11.48 minutes (SD=4.62), nonpanickers for a mean of 12.35 minutes (SD=3.44), and the comparison subjects for a mean of 12.99 minutes (SD=2.98) (n.s.). Panicking patients reached a mean end-tidal CO
2 of 31.27 torr (SD=6.23), while nonpanicking patients reached a mean of 30.28 torr (SD=6.13) and comparison subjects reached a mean of 31.87 torr (SD=6.13). All subjects were encouraged to maintain a respiratory rate of 30 breaths per minute. The actual respiratory rates, shown in
+table 2, were in fact very close to 30 breaths per minute for all groups. The multivariate comparison of pre- and intervention means showed sex effect (F=6.12, df=3,65, p<0.001), time effect (F=94.71, df=3,65, p<0.001), and group-by-sex (F=2.27, df=6,128, p<0.04) and sex-by-time (F=3.84, df=3,65, p<0.01) interactions. Follow-up ANOVAs showed no group differences for any of the measures.
Through use of the last 5-minute means the multivariate analysis showed sex effect only (F=3.49, df=3,68, p<0.02). The follow-up exploratory comparisons showed no significant differences between patients and comparison subjects with respect to tidal volume and respiratory rate. For end-tidal CO2 significant group-by-sex interaction (F=4.49, df=2,78, p<0.01) showed that panicking (mean=26.62 torr, SD=3.51) and nonpanicking (mean=27.35 torr, SD=5.82) female patients had significantly lower end-tidal CO2 than comparison women (mean=32.88 torr, SD=7.63). The same analysis in men would have been meaningless, since only two men panicked during this intervention.
7% CO2 inhalation. The mean times of 7% CO2 exposure for each group were as follows: panickers, mean=8.74 minutes, SD=7.01; nonpanickers, mean=10.87 minutes, SD=6.06; normal comparison subjects, mean=15.75 minutes, SD=5.43 (F=9.88, df=2,75, p<0.001). Comparison subjects lasted significantly longer than either panickers (t=4.38, df=61, p<0.001) or nonpanickers (t=2.60, df=40, p<0.02). The multivariate comparison of pre- and intervention values showed sex effect (F=8.41, df=3,48, p<0.001) and group-by-sex (F=2.80, df=6,94, p<0.02) and sex-by-time (F=9.53, df=3,48, p<0.001) interactions. The follow-up ANOVAs showed significant sex effect (F=20.1, df=1,52, p<0.001) and sex-by-time interaction (F=22.82, df=1,52, p<0.001) for tidal volume, with panicking male patients having the lowest tidal volume response followed by nonpanickers and then by comparison subjects. This difference resulted in significantly higher tidal volume during the last 5 minutes of CO2 exposure in comparison subjects (mean=1,956 ml, SD=729) than in patients (mean=1,613 ml, SD=642) (F=4.23, df=1,69, p<0.04). The significant group-by-sex-by-time interaction for respiratory rate (F=3.25, df=2,52, p<0.05) was due to panicking female patients having a significantly larger respiratory rate increase than comparison subjects (group-by-time interaction in women: F=4.65, df=2,30, p<0.02). Panicking female patients started with a lower respiratory rate but ended up with a higher rate than the other two groups of women.
No overall conventional CO2 sensitivity (delta minute ventilation divided by delta end-tidal CO2) differences emerged with use of either the 5% or the 7% values against baseline. However, sensitivity attributable to respiratory rate showed significant overall sex effect (F=4.67, df=1,71, p<0.03). Female patients (mean=0.20, SD=0.19) had significantly higher response than female comparison subjects (mean=0.10, SD=0.11) during 5% CO2 inhalation (t=–2.25, df=41.44, p<0.03, separate variance estimate). Similar but more accentuated respiratory rate response differences were found during 7% CO2. Here, the significant overall sex effect (F=7.42, df=1,55, p<0.01), group effect (F=3.39, df=2,55, p<0.04), and group-by-sex interaction (F=3.31, df=2,55, p<0.05) were due to significantly higher respiratory rate response in female patients (mean=0.30, SD=0.29) than in female comparison subjects (mean=0.10, SD=0.13) (F=6.03, df=2,31, p<0.01), resulting in a significant overall patient/comparison subject difference in respiratory rate response (patients: mean=0.23, SD=0.25; comparison subjects: mean=0.13, SD=0.13). While it was not significant, panicking patients had the highest response (0.26) followed by nonpanicking patients (0.17) and comparison subjects (0.13).
During 5% CO2 female comparison subjects had higher tidal volume response (mean=40.97, SD=21.78) than female patients (mean=28.63, SD=12.66) (F=2.87, df=2,31, p<0.07), which produced an overall patient/comparison subject difference in tidal volume response (F=2.63, df=2,71, p<0.08).
The recovery period showed similar patterns after all three interventions, with the most pronounced differences emerging after 5% CO
2 was switched to room air. Here, the multivariate analysis showed sex effect (F=3.51, df=3,52, p<0.02), time effect (F=53.61, df=30,25, p<0.001), and a near-significant group effect (F=2.07, df=6,102, p<0.06). The follow-up ANOVA for tidal volume demonstrated that patients maintained their increased tidal volume longer than comparison subjects (group-by-time interaction: F=2.03, df=10,285, p<0.03). While patients lowered their respiratory rate sharply, they showed a delay in lowering their minute ventilation during the recovery phase, resulting in a more pronounced drop in end-tidal CO
2 than in comparison subjects (
+figure 3).
The logistic regression analysis identified low end-tidal CO2 (Wald F=4.16, df=1, p<0.05) for the entire group and low end-tidal CO2 (F=6.20, df=1, p<0.02) and high within-subject standard deviation of minute ventilation (F=4.98, df=1, p<0.01) in women only as baseline predictors of panic.
Received Oct. 15, 1996; revisions received April 14 and May 21, 1997; accepted May 29, 1997. From the Biological Studies Unit, New York State Psychiatric Institute, College of Physicians and Surgeons, Columbia University; New York University Medical Center; and the Phobia, Anxiety and Stress Disorders Clinic, Hillside Hospital, New York. Address reprint requests to Dr. Papp, 722 West 168th St., New York, NY 10032. Supported in part by NIMH grants MH-41778 and MH-30906, Scientist Development Award for Clinicians MH-00858 (Dr. Papp), Research Scientist Award MH-00416 (Dr. Gorman), and Scientist Development Award for Clinicians MH-01039 (Dr. Coplan).