Research to date strongly suggests the influence of biological factors in panic disorder
(1 –
4) . There is evidence that the disorder is familial and may be genetically transmitted
(5,
6) . Individuals with panic disorder are susceptible to having panic attacks in response to various metabolic challenges, ranging from sodium lactate infusion
(7 –
9), caffeine ingestion
(10), and carbon dioxide (CO
2 ) inhalation
(11,
12) to sustained hyperventilation
(13) . In subjects with panic disorder, sustained hyperventilation commonly elicits slower end-tidal CO
2 (pCO
2 ) recovery after cessation
(3,
14) . From several studies, respiratory dysregulation also persists following treatment
(15,
16), implicating “trait” respiratory components to the disease.
During alkalotic challenges, such as hyperventilation, a number of studies have investigated peripheral blood lactate response. In most investigations, blood plasma lactate abnormally increases in subjects with panic disorder to these challenges, even when there is control for subject fitness
(13) or glucose availability
(17) . However, the time course and magnitude of acid-base response are different for the peripheral (blood) and central (brain) compartments
(18) . In prior work with single-voxel hydrogen (
1 H) magnetic resonance spectroscopy (MRS) to measure brain lactate response during regulated hyperventilation (20 mm Hg)
(19), more rapid and prolonged brain lactate increases were observed in asymptomatic medication-treated subjects with panic, a group specifically studied to evaluate the underlying features of the disorder. This evidence for abnormal compensation during alkalosis and increasing discussion of acid-base systems in panic disorder
(20) led to the current investigation aimed at characterizing brain pH response to hyperventilation directly.
Phosphorous (
31 P) MRS provides a noninvasive measure of intracellular brain pH derived from the frequency shift between metabolites, most commonly inorganic phosphate relative to phosphocreatine
(21) . To improve time-resolution while ensuring a sufficient
31 P signal-to-noise ratio, a larger volume of interest (a 30-mm slab) was employed compared to our single-voxel
1 H MRS work
(19) . The diffuse brain lactate response to hyperventilation
(22) and the homogeneous pH distribution in the brain
(23) suggest MRS voxel localization need not be markedly specific to assess this regulatory system. Two metabolic models in panic disorder were posited to explain the previously observed excess lactate increase to hyperventilation challenge: 1) subjects with panic disorder would demonstrate greater alkalosis to hyperventilation, implicating lactate as directly compensatory to an exaggerated brain pH response; 2) subjects with panic disorder would demonstrate reduced alkalosis, implicating lactate as overly compensatory to a normal pH response. In both models, we hypothesized that delayed pCO
2 recovery in subjects with panic disorder, if present, would be related to slower pH normalization in the posthyperventilation phase.
Results
There were no group differences in gender distribution (p=0.36, Fisher’s exact test) (df=1, p=0.21), although the healthy subjects were somewhat younger than the panic disorder group (t=–1.99, df=18, p=0.06). Symptom differences between groups in the Acute Panic Inventory at baseline and after the hyperventilation procedure were demonstrated (baseline—healthy subjects: mean=0.50, SD=0.92; subjects with panic disorder: mean=4.56, SD=4.77) (t=2.77, df=18, p=0.01) (posthyperventilation—healthy subjects: mean=1.63, SD=2.06; subjects with panic disorder: mean=7.39, SD=7.19) (t=–2.54, df=18, p=0.02) and for anxiety ratings (baseline—healthy subjects: mean=0.36, SD=0.67; subjects with panic disorder: mean=1.33, SD=1.00; Mann-Whitney U=21.5, df=18, p=0.03; posthyperventilation: healthy subjects: mean=0.09, SD=0.30; subjects with panic disorder: mean=1.56, SD=1.33; Mann-Whitney U=13.5, df=18, p=0.002). No differences were found at baseline for the panic ratings (healthy subjects: mean=0.23, SD=0.75; subjects with panic disorder: mean=1.11, SD=1.83; Mann-Whitney U=37.0, df=18, p=0.18), although significant differences were present after the hyperventilation paradigm (healthy subjects: all=0; subjects with panic disorder: mean=2.33, SD=1.94; Mann-Whitney U=16.5, df=18, p=0.002). Despite greater endorsement of anxiety symptoms in the panic disorder group, no subject reported a response similar to a panic attack, nor did any subject meet DSM-IV criteria for a panic attack in response to hyperventilation.
In both groups, pCO
2 decreases were demonstrated across the hyperventilation phase for the effect of time (F=191.62, df=19, 342, p<0.001; Greenhouse-Geisser effect=0.001), with significantly different group response-by-time interaction (F=2.84, df=19, 342, p<0.001; Greenhouse-Geisser effect=0.02). The group main effect for CO
2 differences was less than significant (F=3.84, df=1, 18, p=0.07). Post hoc testing of the group-by-time effect demonstrated greater decreases in pCO
2 by scan nine in the panic disorder group (t=1.96, df=18, p=0.07) and significant decreases on the last scan of the hyperventilation period (scan 10) (t=2.72, df=18, p=0.01). In the posthyperventilation phase, pCO
2 recovery was significantly delayed in the panic disorder group on scans 11–17 (all t>2.17, df=18, p<0.04), with a difference observed on scans 18 and 19 (all t<2.01, df=18, p<0.06) (
Figure 2 ). To ensure that the delayed recovery was not simply a result of different group pCO
2 minima at the end of the hyperventilation period, data were transformed for each subject to a percent recovery value (B–R)/(B–H), where B was the mean baseline pCO
2, R was the pCO
2 at each point in recovery, and H was the pCO
2 at the end of hyperventilation. Correction for the greater pCO
2 minimum in the panic disorder group at the end of hyperventilation with this percent recovery formula did not markedly affect pCO
2 results, with group differences observed on scans 11, 12, and 20 (all t>2.75, df=18, p<0.10) and significant differences in pCO
2 recovery for scans 13–19 (all t>2.37, df=18, p<0.03).
Across the respiratory challenge, brain pH increases were demonstrated for both groups (time F=5.36, df=19, 342, p<0.001; Greenhouse-Geisser effect=0.001). However, no differential effects of group (F=1.14, df=1, 18, p=0.30) or group by time were demonstrated (F=0.75, df=19, 342, p=0.77; Greenhouse-Geisser effect=0.64). Despite the greater pCO
2 decrease during the hyperventilation phase, the panic disorder group did not demonstrate increased brain alkalosis; instead, pH values for subjects with panic disorder were not significantly different from comparison values throughout the hyperventilation protocol (
Figure 2 ). In the posthyperventilation phase, there was slower pCO
2 recovery, but pH recovery was not significantly different between the subjects with panic disorder and the healthy comparison subjects.
The discordance between pCO
2 and pH measurements in the panic disorder group led to exploratory analyses aimed to model the degree of pH blunting. With healthy comparison subjects’ pCO
2 and pH values, linear regression models were created for the hyperventilation (five scans) and posthyperventilation (10 scans) phases. Comparison group pCO
2 and pH values were entered into a linear regression for each phase, and slope and intercept terms were generated. Resultant linear regression models were as follows: hyperventilation: pH=–0.0124 * pCO
2 + 7.4178; posthyperventilation: pH=–0.0023 * pCO
2 + 7.228. pCO
2 data for subjects with panic disorder were then entered into these models for healthy comparison subjects, generating estimated pH response values for the panic disorder group. The predicted panic disorder pH curve from these models is shown in
Figure 2, illustrating the discrepancy between the estimated and measured acid-base response.
No differences were found when we compared phosphomonoesters, phosphodiesters, or inorganic phosphate amplitudes by group (all F<1.88, df=1, 18, p>0.19), time (all F<1.06, df=19, 342, p>0.39; Greenhouse-Geisser effect>0.39), or group by time (all F<0.83, df=19, 342, p>0.67; Greenhouse-Geisser effect>0.39). Similarly, no differences were found for phosphocreatine or ATP by group (all F<1.21, df=1, 18, p>0.28), time (all F<1.43, df=19, 342, p>0.11; Greenhouse-Geisser effect>0.19), or group by time (all F<1.29, df=19, 342, p>0.19; Greenhouse-Geisser effect>0.25).
Discussion
Despite similar training and feedback provided during the hyperventilation phase, the subjects with panic disorder had difficulty maintaining a hyperventilation rate that produced a pCO
2 of 20 mm Hg. This overbreathing led to more robust pCO
2 changes in the panic disorder group during the hyperventilation phase. Across this same interval, however, the magnitude of alkalosis was similar between groups. This physiological response pattern, wherein a substantially greater decrease in pCO
2 was not accompanied by a similar magnitude of alkalotic response (a pH increase), supports a model of greater alkalotic buffering in the panic disorder group. Consistent with this idea, estimating panic disorder pH response with data from the healthy comparison subjects as a model acid-base buffering system illustrates that the induced alkalosis should have been greater in the panic disorder group, reaching a maximum equal to the healthy comparison subjects at an earlier hyperventilation time point. This predicted asymptotic pH response with decreasing pCO
2 levels is similar to the pH curve observed in animal work
(31), adding support for the validity of this modeling approach and data interpretation.
In the posthyperventilation period, there was delayed pCO
2 recovery in the panic disorder group assessed with either raw or percent recovery values, but the time course of pH normalization was not substantially different between groups. Recent evidence suggests that this may reflect the insensitivity of pH as a metric to predict hyperventilation maintenance. From a study of hypoxia in rats, dichloroacetate, administered to block lactate production, markedly reduced hyperventilation
(32) . If the converse is operating in subjects with panic disorder, it is possible that elevated lactate in response to hyperventilation might be driving or initiating the persistent respiratory behavior after voluntary hyperventilation cessation. Heuristically, it remains undetermined whether the exaggerated lactate response reported in prior work reflects increased production in panic disorder or whether lactate, once produced, is cleared more slowly, perhaps because of altered enzymatic activity, such as lactate dehydrogenase. This latter point has some support from the genetic literature investigating subjects sensitive to anxiety provocation by CO
2 inhalation
(33) . Further work measuring lactate in parallel with brain pH at high time resolution will be necessary to clarify a mechanistic description of events.
Similar to past work, medication-treated subjects with panic disorder were studied to minimize state-specific anxiety contributions to physiological measures and to allow assessment of underlying trait features
(19) . Although the panic disorder group showed increased levels of anxiety at baseline and after hyperventilation challenge in relation to comparison subjects, symptom severity was substantially lower than would be expected among untreated patients with panic disorder
(24) . Future work evaluating medication-treated healthy comparison subjects or a panic group that symptomatically improved with cognitive behavior therapy will aid in evaluating what effect, if any, medication had on the measured results.
A further point to consider when we interpret these results is the rapid cerebral blood flow (CBF) reduction that occurs in response to hypocapnea
(34,
35) . Altered blood flow response to hyperventilation in subjects with panic disorder has been implicated in two studies
(36,
37) . On this point, pilot work conducted in 1998 on six of seven subjects with panic disorder who participated in the
1 H MRS investigation by Dager et al.
(19) did not demonstrate differential middle cerebral artery flow velocity changes using transcranial Doppler during an identical hyperventilation challenge (unpublished data of S. Dager and K. Beach). In the present study, high energy phosphates (ATP and phosphocreatine) were assessed as an indirect indication of hypoxia. Despite the increased hypocapnea in the panic disorder group, ATP and phosphocreatine levels were not significantly altered, suggesting that this factor, per se, was not a major contributor to the observed pH results.
Further study is necessary to evaluate whether lactate increase is a direct buffering response, as we have postulated, or instead a byproduct of another process affected by the acid-base perturbation. For example, it is likely that hyperventilation-induced increases in pH alter the redox state toward glycolysis
(25) . In this model, lactate increases would reflect a redox shift instead of a primary buffering process. Quantifying the time course and buffering components of such acid-base changes will be important to delineate in future work.
Whether pH and CO
2 dysregulation are also present in panic disorder subgroups not susceptible to lactate infusion (which also causes an alkalotic state) or for patients susceptible to only acidotic perturbations (e.g., CO
2 inhalation)
(38) remains an interesting research question. It is also possible that subgroups of panic subjects sensitive to only alkalotic or acidotic paradigms have pH regulatory abnormalities that are specific in direction and a further subgroup sensitive to both challenge types has altered buffering in general. Future work exploiting higher field strengths to increase signal-to-noise ratio and improve the sensitivity of measuring small differences in pH response and energy metabolism will be helpful to further test and extend these observations.