Abstract
Introduction
Panic attacks (PAs) are highly prevalent in many psychiatric disorders but are most commonly associated with panic disorder (PD). PAs are sudden surges of intense fear or discomfort peaking within minutes. During this period, four or more specified symptoms such as palpitations, a feeling of choking, or dizziness develop (American Psychiatric Association, 2013).
In humans, PAs can be reliably provoked in the laboratory by inhaling air enriched with carbon dioxide (CO2). Various concentrations of CO2 have been examined over the years. Exposure to a low concentration (5–7%) for up to 20 min elicits symptoms as seen in a state of general anxiety and is therefore nowadays considered to represent a model for generalized anxiety disorder (Ainsworth et al., 2015; Bailey et al., 2005, 2007), while a deep breath of 35% CO2 provokes the sudden symptoms and intense fear associated with PAs in patients with PD (for extensive review see Leibold et al., 2015). Interestingly, a comparably intense emotional response can be triggered in healthy individuals when using a higher dosage of CO2 (i.e. taking two vital capacity breaths of 35% CO2 instead of one as in PD patients) (Griez et al., 2007), suggesting a spectrum of CO2-sensitivity.
To study the molecular basis of CO2-sensitivity in depth, the use of CO2 as an experimental stimulus in animal models has gained increasing attention in recent years (D’Amato et al., 2011; Johnson et al., 2012; Taugher et al., 2014; Vollmer et al., 2016; Winter et al., 2017; Ziemann et al., 2009). As it is not feasible to apply one or two vital capacity breaths of 35% CO2 in rodents as it is done in humans, a lower percentage of up to 20% is commonly used for 10–20 min to allow assessment of behavior in a testing box filled with CO2. Using this methodology, studies showed that CO2 elicits robust fear-related behavior in rodents. A major aim of such rodent research is the translation of findings to humans. However, it is uncertain how well behavioral assessments in rodents correspond to the self-reported experience in humans. To address this challenge we used CO2 as the experimental stimulus in both rodents and humans and went beyond the traditional outcome measurements by recording additional physiological parameters (Leibold et al., 2016). Next to the finding that CO2 induces a strong fear-related behavioral response in mice and humans alike, we showed comparable cardio-respiratory effects. In detail, we calculated and statistically compared effect sizes, taking into account that the frequency per minute varies between both species. This allowed us to perform a quantitative comparison to demonstrate that exposure to a low concentration of CO2 in rodents corresponds to the response to a 35% CO2 inhalation in humans. Therefore, while the paradigms are not identical, the comparable effects of the cardio-respiratory system in addition to the behavioral fear response suggests that this rodent model represents a model of panic and is a promising approach to bring closer research in rodents and humans.
The mechanisms of how CO2 causes fear and exerts effects such as respiratory stimulation are not fully understood. Under physiological conditions the concentration of CO2 and pH level within the body are maintained within a narrow range. Small fluctuations are not consciously perceived, while larger acute changes typically induce adaptive physiological responses, sensations such as breathlessness, feelings of arousal and panic, and avoidance behavior (Brannan et al., 2001; Diaper et al., 2011; Ziemann et al., 2009). Emotions, in turn, can influence breathing (Homma and Masaoka, 2008) and thereby affect the respiratory response to CO2. Chemoreceptors that detect changes in CO2/pH are located in the peripheral carotid bodies and throughout the entire brain. For instance, Mitchell et al. did groundbreaking work on identifying the location of CO2-chemosensitivity in the medulla (1963b) and showed that increased CO2 in this region mediates an increase in breathing (1963a). Later experiments showed that also the medullary retrotrapezoid nucleus (Guyenet and Bayliss, 2015), and subsets of serotonergic neurons in the medullary raphe nucleus and the dorsal raphe nucleus are pH-sensitive (Richerson, 2004). The retrotrapezoid nucleus receives excitatory serotonergic innervation from both of these raphe nuclei, which can shift the CO2 threshold and stimulate respiration (Guyenet and Bayliss, 2015; Mulkey et al., 2007). Particularly serotonin (5-HT) neurons of the medullary raphe nucleus have been associated with driving the respiratory response to CO2. In contrast, midbrain 5-HT neurons have been associated with arousal and anxiety (Richerson, 2004).
In line with these basic research findings, many studies in humans provided support for a role of the 5-HT system in emotional sensitivity to CO2. For instance, tryptophan depletion, a method to reduce brain 5-HT levels, increased the fear response to CO2 inhalation in PD patients and healthy individuals (Miller et al., 2000; Schruers et al., 2000; Schruers and Griez, 2003), whereas administration of a 5-HT precursor, increasing brain 5-HT levels, reduced the response (Schruers et al., 2002). Blocking the 5-HT transporter (5-HTT) and thereby increasing the extracellular availability of 5-HT also reduced the response to CO2 (e.g. Bertani et al., 1997; Perna et al., 2002; Schruers and Griez, 2004). Further, our group showed that the fear reaction to CO2 is moderated by a functional polymorphism in the transcriptional control region of the gene encoding the 5-HTT (polymorphism referred to as 5-HTT gene-linked polymorphic region, 5-HTTLPR). In brief, subjects homozygously carrying the long (L)-allele reported more fear than short (S)-allele carriers (Schruers et al., 2011). In line with this, we also reported a negative association between DNA methylation in the 5-HTT regulatory region and emotional CO2-sensitivity (Leibold et al., 2020).
Making use of our rodent model and as translational follow-up to our previous studies in humans, we here aimed to explore whether partial or complete loss of 5-HTT moderates the response to CO2. For this purpose, wildtype (WT), hetero- (+/–), and homozygous (–/–) 5-HTT knockout mice were exposed to CO2, while the response was examined in terms of fear-related behavior and by monitoring cardio-respiration. Based on the findings of our previous studies, we hypothesized that wildtype mice would show the strongest sensitivity to CO2.
Materials and methods
For overview see Figure 1.

Overview of experimental design.
Fear-related behavior
Animals
Male 5-HTT deficient (Slc6a4tm1Kpl; 5-HTT+/– (
Behavioral testing
Modified open field test (OFT)
All animals were first tested in the modified OFT. The OFT chamber consisted of a transparent Plexiglas square (50×50×40 cm), covered with a clear lid, allowing live monitoring of behavior during gas exposure. Gas was infused using an infusion port on the upper part of the box. Overpressure was prevented by two holes at the opposite site. In addition, two computer fans (5×5 cm, 26 dB), fixed to the lid, were used to ensure a homogenous concentration of the infused gas throughout the entire chamber. The CO2 concentration was constantly controlled using a 30% CO2 Sampling Data Logger (CO2 meter, Ormond Beach, Florida, USA).
Mice were placed in the center of the chamber that was pre-filled with either compressed air or 10% CO2 (premixed gas tanks, SOL Nederland B.V., Landgraaf, The Netherlands) depending on the experimental group. Movements were scored automatically with a computerized system (Ethovision Color Pro, Noldus, The Netherlands) for 20 min. For analysis purposes, the floor was subdivided into a 30×30 cm central zone, 10×10 cm corners, and 30×10 cm walls. Time spent in the different zones and total distance moved were analyzed. At the end of the trial, the boxes were cleaned thoroughly.
Behavioral testing was performed under low light conditions during the animals’ dark phase and by an investigator blind to genotype.
Two-chamber test (TCT)
After completing the modified OFT, mice were tested in the TCT. The apparatus consisted of two Plexiglas chambers (each 50×25 cm with 40 cm high walls) that were connected by an open door (3.5×3.5 cm) to allow free crossing, whilst limiting mixing of the gasses between chambers. Each chamber had a separate gas infusion port, a hole to prevent overpressure, and a computer fan (5×5 cm, 26 dB) to ensure a homogenous gas concentration throughout the chamber.
For mice assigned to air exposure, both chambers were pre-filled with air. For mice subjected to CO2 exposure, a chamber was pre-filled with 10% CO2 and the other one with air. The chamber that was filled with CO2 and the chamber in which mice were placed at the beginning of the test were randomized. The concentration of CO2 was continuously controlled using a 30% CO2 Sampling Data Logger (CO2 meter, Ormond Beach, Florida, USA; a steady state of 9% and 2% was reached in case of CO2 infusion). Movements, number of crossings, and time spent in each chamber were automatically scored with a computerized system (Ethovision Color Pro, Noldus, The Netherlands) for a period of 10 min. At the end of each trial, the box was thoroughly cleaned. As in the modified OFT, testing was done under low light conditions during the animals’ dark phase and by an investigator blind to genotype.
CO2-evoked freezing
In both behavioral tests, freezing was scored by a trained observer blind to exposure and genotype. Freezing was defined as absence of any movements apart from respiration and is considered to reflect fear-related behavior in rodents (Mongeluzi et al., 2003).
CO2-induced plasma corticosterone secretion
To assess the hypothalamo-pituitary-adrenal (HPA) axis responsiveness to CO2, plasma corticosterone levels were determined. Immediately after taking the mice from their home cage, the first blood sample was drawn via saphenous vein puncture (pre-CO2 level) using heparinized blood collection tubes (Microvette CB300, Sarstedt, Germany). Then, mice were returned to their home cage and exposed to 10% CO2 for 20 min, after which a second blood sample was taken (post-CO2 level). Blood samples were kept on ice and centrifuged at 2000 g for 10 min at 4°C. Plasma was collected and stored at −80°C until further processing. Blood corticosterone levels were analyzed using radioimmunoassay. An ImmunuChem Double antibody corticosterone 125I RIA Kit for rodents (MP Biomedicals, Orangeburg, New York, USA) was used according to the manufacturer’s instructions. Per sample 5 µL plasma was diluted in steroid diluents (1:100 for pre-CO2 and 1:200 for post-CO2 samples). Fifty μL of corticosterone 125I and 100 μL anti-corticosterone were added per 50 µL of diluted sample, after which the tubes were incubated for 2 h at room temperature. Then, 250 µL precipitant solution was added to each sample and subsequently centrifuged at 2300 rpm for 15 min. After removing the supernatant the precipitation was counted in a Wizard gamma counter 2470 (PerkinElmer, Waltham, Massachusetts, USA).
Cardio-respiratory response
In a next step, separate groups of WT and 5-HTT +/– knockout mice (background C57BL/6,
Statistical analysis
All analyses were done using the Statistical Package for the Social Sciences (SPSS 25, Inc., Chicago, Illinois, USA). Statistical significance was considered as
Results
Fear-related behavior
Modified OFT
In the modified OFT, significant overall effects of genotype (

Assessment of the behavioral performance in the modified open field test. (a) Under air exposure serotonin transporter (5-HTT) –/– mice covered less distance than mice of other genotypes. Exposure to carbon dioxide (CO2) significantly reduced the total distance moved in all genotypes. (b) During air exposure 5-HTT –/– mice spent less time in the center zone in comparison with wildtype (WT) mice. When exposed to CO2, only WT mice spent less time in the center compared to air exposure, no difference was found in 5-HTT +/– and 5-HTT –/– mice. (c) In all genotypes, CO2 exposure resulted in a robust freezing response compared to air exposure. Bars represent mean + standard error of the mean (SEM). +/–: heterozygous 5-HTT knockout mice; –/–: homozygous 5-HTT knockout mice. *
With regard to time spent in center (Figure 2(b)), there was no overall genotype effect (
Moreover, regarding freezing (Figure 2(c)), overall analysis showed that there was no main effect of genotype (
TCT
In the TCT, with one chamber filled with CO2 and one with air, or both filled with air, overall effects of genotype (

Assessment of the behavioral performance in the two-chamber test. (a) The total distance moved was strongly reduced under carbon dioxide (CO2) exposure in all genotypes. Serotonin transporter (5-HTT) –/– mice covered less distance than wildtype (WT) and 5-HTT +/– mice during air exposure. (b) Under CO2 exposure the number of crossings was significantly lower than during air exposure. (c) When exposed to air, no genotype differences were found. When exposed to CO2, a marked freezing response was observed only in WT and 5-HTT +/– mice compared to air-exposure. 5-HTT –/– mice froze less to CO2 than WT mice. (d) The ratio of time spent and freezing per chamber indicated that WT and 5-HTT +/– mice froze longer in the chamber filled with CO2 than 5-HTT –/– mice. No effect was found within the chamber filled with air. Bars indicate mean + standard error of the mean (SEM). +/–: heterozygous 5-HTT knockout mice; –/–: homozygous 5-HTT knockout mice. ***
Further, genotype did not significantly affect the number of crossings between the chambers (
Evidently, the freezing duration influences the time spent in a particular chamber. To correct for this, for mice exposed to CO2 in one chamber and to air in the other chamber a ratio of time spent freezing relative to the total time spent in that particular chamber was calculated (Figure 3(d)). Within the CO2 chamber, an overall effect of genotype was found (
CO2-induced plasma corticosterone secretion
Repeated measures ANOVA on plasma corticosterone levels showed a significant effect (

Plasma corticosterone secretion before and after exposure to 10% carbon dioxide (CO2) for 20 min. CO2 exposure strongly increased corticosterone levels in mice of all genotypes (exposure effect
Breathing and heart rate recordings
In a next step, the cardio-respiratory response was assessed in WT and 5-HTT +/– mice only, i.e. excluding 5-HTT –/– mice, because these genotypes might represent the best model for the human 5-HTTLPR that is associated with a reduced and not complete deficiency of 5-HTT, and because re-analyzing the behavioral response of WT and 5-HTT +/– mice retained a significant genotype × exposure interaction (
Analyzing breathing frequency, repeated measures analysis showed no main effect of genotype (

Respiratory and cardiovascular monitoring during exposure to air and carbon dioxide (CO2). (a) Schematic representation of a 2 s epoch of pressure-induced changes to assess breathing frequency and tidal volume (amplitude) during inhaling air (top) and CO2 (bottom). CO2 strongly increased breathing frequency in comparison with air exposure. (b) Heart rate decreased during inhaling CO2. No genotype effect was present in any outcome measurement. Bars represent mean + standard error of the mean (SEM). +/–: heterozygous 5-HTT knockout mice; 5-HTT: serotonin transporter; WT: wildtype mice. ***
Carbon dioxide (CO2)-induced changes in respiratory and cardiovascular measurements.
5-HTT: serotonin transporter; 5-HTT +/–: heterozygous 5-HTT knockout mice; WT: wildtype mice.
Values represent mean (standard error of the mean (SEM)).
Discussion
As follow-up to our previous study in humans showing a differential self-reported fear response to CO2 depending on the 5-HTTLPR genotype (Schruers et al., 2011), we further investigated the relationship between chemosensitivity and the 5-HT system by exposing WT, 5-HTT +/–, and 5-HTT –/– mice to CO2. In line with our previous studies in humans (Leibold et al., 2016; Schruers et al., 2011) and mice (Leibold et al., 2016), CO2 induced robust fear-related behaviors and strongly activated the respiratory system. However, in contrast to our hypothesis, there were only minor effects of differential 5-HTT expression on behavior.
In humans, the use of CO2 as experimental panic-stimulus is well established in PD patients and, when using a higher dosage, also in healthy individuals (for review see Leibold et al., 2015). This is suggestive of the involvement of fundamental neurobiological mechanisms. Inhaling air enriched with CO2 is associated with an acute drop in brain pH (Schuchmann et al., 2006; Ziemann et al., 2008, 2009) that can trigger various adaptive responses, as such a pH drop can have life-threatening consequences. In this study, we investigated the responses to CO2 by using two behavioral tasks (modified OFT and TCT) and by assessing cardio-respiration. Generally, the OFT is considered an anxiety test, in which avoidance of the center zone (open zone) is one of the main outcome measurements. It has to be mentioned that anxiety is distinct from fear/panic and that different brain structures and behaviors are involved (Blanchard and Blanchard, 1990). Anxiety predominates when the distance to a threat is large, i.e. when there is a
Examining the effects of CO2 in these tasks showed that CO2 induced a robust fear-related behavioral response in mice, in line with previous reports (Johnson et al., 2011, 2012; Leibold et al., 2016; Ziemann et al., 2009). However, unexpectedly, there were only minor effects of 5-HTT genotype. More specifically, total distance moved was strongly decreased in both tasks, independent of genotype. Likewise, the total duration of freezing, a correlate of fear and panic in mice (Mongeluzi et al., 2003; Ziemann et al., 2009), in the modified OFT and the number of crossings between compartments in the TCT were similarly affected in all genotypes. Exceptions are the duration spent in the center zone of the modified OFT, an outcome measurement that is relatively independent on locomotor activity, and freezing in the TCT. With regard to the duration spent in the center zone of the modified OFT, 5-HTT –/– mice show no decrease in time spent in the center zone during CO2 exposure compared to air exposure. This suggests an altered threat perception and/or altered coping strategy in mice carrying the full knockout. For instance, the fact that 5-HTT –/– mice already spent very little time in the center during air exposure could represent heightened aversion to the OFT itself in this genotype, leading to a ceiling effect. In the TCT, however, 5-HTT –/– mice spent less time freezing over the total time spent in the aversive compartment. This, again, could indicate that mice carrying the full knockout perceive CO2 differently or have an altered coping strategy. Furthermore, it cannot be excluded that the TCT is less aversive in comparison to the OFT, as the compartments are smaller. These observations are in line with previous work showing, e.g., that mice fully deficient of 5-HTT performed worse in the more stressful water maze test compared to WT and 5-HTT +/– littermates, while no differences were observed in Barnes maze (Karabeg et al., 2013).
Next to direct effects on animal behavior, we aimed to measure avoidance of CO2 as an aversive stimulus. In the TCT, animals had the option to leave the compartment filled with CO2. Unexpectedly, we did not observe any avoidance, in contrast to previous studies in mice, in which WT mice spent most of the time in the chamber filled with air (D’Amato et al., 2011; Ziemann et al., 2009). This discrepancy can be explained by the immense freezing response that we observed in the present study when CO2 was present, which seems to have prevented the animals from leaving the CO2 chamber. This notion is supported by the strongly reduced crossing frequency between both chambers in the presence of CO2 Ziemann et al. (2009) used 15% CO2 in this task in their previous study. In our study, a steady state of 9% CO2 was reached, which is lower than the previously used concentration of 15% CO2. In addition, care was taken to use an optimized methodological approach, including fans to homogenize the air, continuous measurement of the CO2 concentration, and a randomization regarding side of CO2 administration and side, in which mice were placed at the beginning of the test. Moreover, the chambers of our set-up were bigger compared to the set-up used by Ziemann et al. The combination of a larger box in our study and CO2 may have been a particularly aversive stimulus. Possibly, this led to a ceiling effect and prevented finding any subtle genotype differences under CO2 exposure. CO2 exposure is stressful in nature as shown by the observed increase in plasma corticosterone, which is in line with previous reports (D’Amato et al., 2011). This effect of CO2 is independent of genotype and suggests that differences in behavioral performance are not primarily dependent on HPA axis reactivity.
In addition, to studying the behavioral effects of CO2 exposure, also the effects on the cardio-respiratory system were explored. CO2 strongly activated the respiratory system with an increase in breathing frequency and tidal volume; heart rate decreased in comparison to air exposure. These results are consistent with those previously reported in, for example, freely moving rats (Annerbrink et al., 2003) and phrenic nerve recordings of anesthetized rats under CO2 exposure (Dumont et al., 2011). In the present study, effects were independent of genotype. To the best of our knowledge, little research has been done regarding examining the respiratory response to CO2 in mice deficient for 5-HTT or in humans related to the 5-HTTLPR genotype. Li and Nattie (2008) observed in their study that 5-HTT –/– males had a smaller increase in ventilation to CO2 than WT mice. 5-HTT +/– mice were not tested. The lack of a difference between WT and 5-HTT +/– mice in our study might be explained by 5-HT neurons sufficiently acting as chemosensors and initiating an adaptive response (Richerson, 2004), independent of 5-HT availability, or by a functionally insufficient difference in extracellular 5-HT availability to exert effects in this regard.
In addition to the potential explanations above that could contribute to the lack of clear differential 5-HTT expression effects in our study, another possibility is that several genetic and environmental vulnerability factors have to interact to lead to an altered functioning in the brain and a phenotype of heightened CO2-sensitivity. For instance, likewise to 5-HT neurons, many other neurons including noradrenergic neurons in the locus coeruleus (LC) are activated by CO2/pH (Biancardi et al., 2008; Filosa et al., 2002; Pineda and Aghajanian, 1997). The LC receives major input from the medullary nucleus paragigantocellularis, a key region in cardio-respiratory regulation (Aston-Jones et al., 1986), and was shown to modulate the respiratory response to CO2 (Biancardi et al., 2008). Furthermore, LC neurons are involved in the behavioral response to stress (Stanford, 1995) and by projecting to nuclei of the limbic system these neurons likely mediate fear. Based on basic research findings like these and other research, it has been suggested that the central noradrenergic system, particularly the LC, is the common site for mediating the increases in fear and physiology to a 35% CO2 inhalation in humans (Bailey et al., 2003). It is noteworthy that the LC and the dorsal raphe nucleus are connected by extensive reciprocal connections (Kim et al., 2004). CO2 was shown to stimulate the release of 5-HT in the LC. By exerting inhibitory effects, 5-HT reduces the stimulatory effects of the LC on breathing in response to CO2 (de Souza Moreno et al., 2010). Contrary to this effect, in general, projections from the LC have excitatory effects on dorsal raphe nucleus neurons (Brown et al., 2002). These connections make it interesting to investigate, e.g. pharmacologically, whether disturbances in both systems might interact to express a phenotype of heightened CO2-sensitivity. While it is known that both 5-HT and noradrenergic neurons are sensitive to pH/CO2, it is unknown which neuronal molecule senses these changes and an alteration in that molecule could also contribute to a higher sensitivity to CO2. An interesting candidate is the acid-sensing ion channel 1a (ASIC1a). Genetic knockout of ASIC1a significantly blunts the CO2-induced behavioral effects (Ziemann et al., 2009) and antagonism normalizes the hyperreactivity to CO2 that is seen in a specific mouse model (Battaglia et al., 2019). Further, in addition to genetic factors, also their interaction with environmental stimuli shape behavioral phenotypes. For example, early life adversities such as maternal separation and repeated cross-fostering have been shown to increase the sensitivity to CO2 in rodents (Battaglia et al., 2019; D’Amato et al., 2011). Therefore, interplays of genetic vulnerabilities and environmental stressors could provide important insights for understanding the pathophysiology in more depth.
Some potential limitations of the present study should be kept in mind: As only male animals were tested no potential sex effects were examined. In humans, a higher prevalence of PD is reported in females (Dick et al., 1994; Gater et al., 1998), while CO2 studies did not consistently find this effect (Kelly et al., 2006; Leibold et al., 2013; Nardi et al., 2007; Nillni et al., 2012; Poonai et al., 2000). Therefore, more research into the effect of sex and menstrual cycle phase is warranted. Further, cardio-respiration was not assessed in 5-HTT –/– mice. As mice of this genotype were affected differentially in the behavioral tasks, physiological alterations associated with CO2 exposure should also be assessed in this genotype in future studies.
Conclusion
Overall, under CO2 exposure there was no effect of differential 5-HTT expression on physiological parameters and only minor interaction effects on behavior. Future research is needed that also addresses the interplays with other systems. These insights could contribute to a better understanding of the pathophysiology underlying CO2-sensitivity and eventually to develop better treatments.
