Skip to main content

Hibernation reduces GABA signaling in the brainstem to enhance motor activity of breathing at cool temperatures

Abstract

Background

Neural circuits produce reliable activity patterns despite disturbances in the environment. For this to occur, neurons elicit synaptic plasticity during perturbations. However, recent work suggests that plasticity not only regulates circuit activity during disturbances, but these modifications may also linger to stabilize circuits during future perturbations. The implementation of such a regulation scheme for real-life environmental challenges of animals remains unclear. Amphibians provide insight into this problem in a rather extreme way, as circuits that generate breathing are inactive for several months during underwater hibernation and use compensatory plasticity to promote ventilation upon emergence.

Results

Using ex vivo brainstem preparations and electrophysiology, we find that hibernation in American bullfrogs reduces GABAA receptor (GABAAR) inhibition in respiratory rhythm generating circuits and motor neurons, consistent with a compensatory response to chronic inactivity. Although GABAARs are normally critical for breathing, baseline network output at warm temperatures was not affected. However, when assessed across a range of temperatures, hibernators with reduced GABAAR signaling had greater activity at cooler temperatures, enhancing respiratory motor output under conditions that otherwise strongly depress breathing.

Conclusions

Hibernation reduces GABAAR signaling to promote robust respiratory output only at cooler temperatures. Although frogs do not ventilate lungs during underwater hibernation, we suggest this would be beneficial for stabilizing breathing when the animal passes through a large temperature range during emergence in the spring. More broadly, these results demonstrate that compensatory synaptic plasticity can increase the operating range of circuits in harsh environments, thereby promoting adaptive behavior in conditions that suppress activity.

Background

Animals have the mysterious ability to produce reliable behaviors while navigating environments that should otherwise disrupt activity of the nervous system [1]. This is thought to occur, in part, through a set of compensatory mechanisms that sense activity perturbations, and in turn, adjust neuronal properties to counteract the disturbance. For example, if activity falls, neurons increase synaptic excitation, decrease inhibition, and alter ion channel densities to enhance neuronal excitability [2,3,4,5]. This framework, termed “homeostatic plasticity,” has been foundational for the past 30 years [6], but it has been challenging to link ethologically relevant behaviors in animals, as neurons often need to be perturbed artificially or pathologically to reveal homeostatic responses [7,8,9,10,11,12].

Certain organisms inhabit natural environments with potent abiotic stressors, such as pH, temperature, and hypoxia, which can cause activity challenges within the nervous system [13]. An extreme example of this problem is encountered by hibernating frogs. Like most vertebrate species, rhythmic neural circuits in the brainstem generate breathing to meet metabolic demands of the organism. However, in the cold hibernation environment frogs may spend several months underwater using only skin for gas exchange [14], while respiratory motor behavior is completely suspended [15]. To successfully emerge following months underwater, frogs employ various compensatory mechanisms to maintain these networks so that they can work effectively when needed. This includes a classic activity-dependent mechanism of compensatory plasticity called “synaptic scaling” at excitatory synapses [16, 17], and also metabolic adjustments that improve performance in hypoxia [18, 19]. As amphibians use synaptic plasticity to overcome this large environmental challenge, they provide an intriguing opportunity to understand how plasticity mechanisms are integrated within animals to promote adaptive behavior.

Compensatory forms of synaptic plasticity are often interpreted as a conceptually simple regulatory regime used to counteract activity disturbances (e.g., when activity goes down, excitatory synaptic strength goes up). However, experimental and modeling work indicates that when circuits are disrupted, compensatory adjustments can also “prime” circuits to perform better during future perturbations and influence the capacity for subsequent plasticity [20,21,22]. Thus, we hypothesized that the hibernation environment not only triggers synaptic plasticity, but also that it induces network modifications to expand the operating range in challenging environments encountered following the initial activity challenge. To test this in an ethological context, we built upon previous work demonstrating that hibernation drives apparent homeostatic plasticity at excitatory synapses [16, 17] and addressed whether inhibitory synapses contribute as well. We then asked how plasticity through inhibition influences motor function of breathing at colder temperatures, where neuronal activity is otherwise markedly suppressed. Here, we report that hibernation strongly reduced GABAergic inhibition, consistent with a classic response to chronic network inactivity [3, 23]. Despite the critical role of inhibition in this network [24], circuit function appeared surprisingly normal at warm temperatures. However, when we assessed activity across a range of temperatures, GABAAR plasticity served mainly to offset the depressive effects of acute cooling, boosting the strength and frequency of motor output at lower temperatures. Given that these animals must pass through a range of temperatures as they reestablish life on land after hibernation [25, 26], network modifications that improve activity over a temperature range seem to play an adaptive role. Therefore, although GABAergic plasticity has no obvious impact on baseline motor function, it expands the operating range of the respiratory circuit in an otherwise suppressive environment to promote adaptive behavior.

Results

Respiratory motor output in adult frogs is thought to arise via central pattern-generating mechanisms consistent with reciprocal inhibition involving GABA [24]. Rhythmic activity in these premotor neurons is then carried to motor pools through interneuronal pathways that control motor outflow to respiratory muscles [27]. Activity of this network can be assessed in isolated brainstems in vitro which produce rhythmic motor output on a variety of cranial nerves including the trigeminal, vagal, and hypoglossal nerve root. Although the activity pattern of these nerves may differ, the output produced is largely synchronous providing a representation of the motor output that resembles lung ventilation in the whole animal [28, 29]. Here, we measured activity of the vagus nerve to monitor respiratory-related motor output from the intact brainstem preparation (cranial nerve X; CNX; schematic in Fig. 1).

Fig. 1
figure 1

GABAergic tone that controls respiratory burst frequency strongly decreases following hibernation. Integrated vagus nerve activity (CN X) was rhythmic in fully intact brainstem preparations (top right box) from control and hibernated animals. A Top, continuous recording of CNX during perfusion of increasing doses of bicuculline in a control preparation (black). Bottom, zoomed in view of epochs (orange boxes) indicated in continuous trace directly above. In control preparations, bath application of bicuculline led to a dose dependent decrease in respiratory burst frequency and dose-dependent increase in non-respiratory burst frequency. B Top, continuous recording of CNX during perfusion of increasing doses of bicuculline in a preparation from a hibernated animal (blue). Bottom, zoomed in view of epochs (orange boxes) indicated in continuous trace above. Respiratory output from hibernated preparations was insensitive to systemic application of bicuculline, however non-respiratory motor bursting increased in a dose-dependent manner like controls. C Summary of bicuculline-mediated changes in respiratory burst frequency. There was a significant interaction between drug and group by two-way ANOVA (F(4,32) = 6.741; p = 0.0005). Additionally, at 10 µM bicuculline, respiratory burst frequency from controls was significantly lower than hibernated preparations by Holm-Sidak’s post-hoc test (p = 0.0111; n = 5/group). D Summary of bicuculline-mediated changes in non-respiratory burst frequency. Bicuculline mediated similar increases in non-respiratory bursting in control and hibernated groups. There was no significant interaction between drug and group by two-way ANOVA (F(4,32) = 0.7948; p = 0.5373; n = 5/group). a.u., arbitrary units. In summary plots, thick lines represent the group average and thin lines are responses of individual preparations

To investigate plasticity in GABAergic control of the respiratory network, we first bath-applied bicuculline, a GABAA receptor antagonist, to intact brainstem preparations. Consistent with an important role for GABAARs in rhythmic motor activity, bicuculline caused a dose-dependent decrease in burst frequency in control brainstems at room temperature (Baseline, 42.60 ± 16.32 bursts/min; 0.5 µM Bicuculline, 26.32 ± 19.68 bursts/min; 1 µM Bicuculline, 16.86 ± 17.03 bursts/min; 5 µM Bicuculline, 5.66 ± 3.61 bursts/min; n = 5, Fig. 1A,C), with most having little-to-no activity at the highest dose (10 µM Bicuculline, 1.74 ± 1.10 bursts/min; n = 5, Fig. 1C). Under the same experimental conditions, hibernators strongly resisted bicuculline (Baseline, 21.12 ± 6.53 bursts/min; 0.5 µM Bicuculline, 23.80 ± 8.72 bursts/min; 1 µM Bicuculline, 21.52 ± 7.71 bursts/min; 5 µM Bicuculline, 14.50 ± 5.54 bursts/min; n = 5, Fig. 1B,C), with respiratory activity persisting at the highest dose (10 µM Bicuculline, 13.50 ± 4.12 bursts/min; n = 5, Fig. 1B,C). In regard to motor amplitude, two-way ANOVA revealed a significant main effect of bicuculline on normalized motor amplitude (F(2.055, 16.44) = 4.499, P = 0.0269) but not a significant effect of group or an interaction between bicuculline and group as we saw for frequency. However, this bicuculine effect is likely to be minor, as post hoc analysis showed that significantly larger bursts compared to baseline occurred at 5 μM in controls (p = 0.010; Holm-Sidak Multiple Comparison’s test), while all other comparisons to baseline were not significantly different from baseline.

Loss of GABAAR signaling appeared to be localized to breathing circuits. In both controls and hibernators, non-respiratory motor activity (long duration bursts with qualitatively different shape than respiratory bursts [30, 31]) were rare at baseline with only 1 control preparation containing these bursts. Interestingly, and despite a clear group x bicuculine interaction for respiratory burst frequency, we did not observe this same response for non-respiratory bursts (group x bicuculline interaction in two-way ANOVA; p = 0.5373). Non-respiratory bursts emerged similarly in a manner that depended on the dose of bicuculline (Fig. 1A-B,D) in both controls (Baseline, 0.04 ± 0.09 bursts/min; 0.5 µM Bicuculline, 0.44 ± 0.39 bursts/min; 1 µM Bicuculline, 1.44 ± 1.05 bursts/min; 5 µM Bicuculline, 3.18 ± 0.81 bursts/min; 10 µM Bicuculline, 3.08 ± 1.01 bursts/min, n = 5) and hibernators (Baseline, 0.00 ± 0.00 bursts/min; 0.5 µM Bicuculline, 0.28 ± 0.63 bursts/min; 1 µM Bicuculline, 0.56 ± 0.98 bursts/min; 5 µM Bicuculline, 2.62 ± 1.50 bursts/min; 10 µM Bicuculline, 3.60 ± 1.86 bursts/min, n = 5). Given the difference in how respiratory and non-respiratory bursts alter sensitivity to bicuculline, hibernation strongly reduces the role of GABAARs in a way that appears specific to circuits that generate breathing but not non-respiratory behaviors.

The previous experiment demonstrated that GABAARs have a strongly reduced role in the production of the respiratory motor pattern after hibernation. Next, we assessed whether plasticity influenced the capacity of the respiratory network to be modulated by GABAAR activation. Therefore, we bath-applied muscimol, a GABAAR agonist, to globally activate all GABAARs that influence activity. In controls, the lowest dose of muscimol (500 nM) did not alter respiratory burst frequency (Baseline, 8.27 ± 3.7 bursts/min vs. 500 nM muscimol, 7.20 ± 4.0 bursts/min; n = 5). Raising the concentration to 1 µM caused a significant decline (1 µM muscimol, 3.47 ± 1.9 bursts/min; p = 0.0282; n = 5; Fig. 2A,C), and exposure to 3 μM silenced all preparations (5 out of 5). In contrast, preparations from hibernators did not significantly alter respiratory burst frequency between baseline and any dose of muscimol (Baseline, 11.00 ± 2.1 bursts/min; 500 nM muscimol, 11.93 ± 1.8 bursts/min; 1 µM muscimol, 14.47 ± 2.8 bursts/min; 3 µM muscimol, 6.47 ± 7.3 bursts/min; n = 5; Fig. 2B-D), with 5 out of 5 preparations maintaining respiratory bursting at the highest dose. While suppression of network activity with both the antagonist and agonist may seem paradoxical, these results are consistent with the necessity of phasic inhibition for generation of respiratory activity and neuronal silencing through further activation of GABAergic inhibition, respectively. Indeed, in other circuits that use inhibition for rhythmogenesis, antagonism of GABAARs disrupts network dynamics required for network activity [32], while applying agonists of GABAARs also suppress activity through membrane hyperpolarization [33]. With both agonist and antagonist experiments, these results show that hibernation leads to a large reduction in the ability of GABAARs to generate and modulate breathing.

Fig. 2
figure 2

Hibernation decreases sensitivity to GABAAR agonist (muscimol)-mediated respiratory frequency decline. Integrated vagus nerve activity (CN X) was rhythmic in fully intact brainstem preparations (top left box) from control and hibernated animals. A Continuous recording of CNX during perfusion of increasing doses of muscimol in a control preparation (black). In control preparations, respiratory output was unchanged by bath application of 500 nM muscimol but decreased following elevation to 1 µM muscimol and ultimately stopped following further elevation to 3 µM muscimol. Preparations were recovered by exposure to bicuculine, demonstrating specificity of the GABAA receptor agonist. B Continuous recording of CNX during perfusion of increasing doses of muscimol in a preparation from a hibernated animal (blue). Preparations from hibernated animals were relatively insensitive to muscimol mediated changes in respiratory burst frequency compared to controls. C Summary of muscimol-mediated changes in respiratory burst frequency. There was a significant interaction between drug and group by two-way ANOVA(F(3,24) = 4.045; p = 0.0184). Additionally, at 1 µM muscimol, respiratory burst frequency from controls was significantly lower than hibernated preparations by Holm-Šidák’s post-hoc test (p = 0.0007, n = 5/group). D Summary of muscimol-mediated changes in respiratory burst frequency relative to baseline. There was a significant interaction between drug and group by two-way ANOVA(F(3,24) = 7.184; p = 0.0013). Additionally, at 1 µM muscimol, respiratory burst frequency from controls was significantly lower than hibernated preparations relative to baseline by Holm-Sidak’s post-hoc test (p = 0.0030, n = 5 per group). In summary plots, thick lines represent the group average and thin lines are individual preparations

The previous experiments identify that hibernation reduces the role of GABAARs in generating and modulating the respiratory rhythm. At the level of the motoneuron, an activity-dependent mechanism known as “synaptic upscaling” strengthens excitatory synapses in response to hibernation [16, 17], aligning with classic homeostatic responses to chronic network inactivity [4]. Thus, we tested if hibernation also “downscales” postsynaptic GABAARs to reduce GABAAR inhibition. We used patch clamp electrophysiology to record miniature inhibitory postsynaptic currents (mIPSCs) carried by GABAARs in identified vagal motoneurons that cause breathing in amphibians. The GABAAR currents were isolated in tetrodotoxin (TTX), strychnine (strych), and DNQX to block presynaptic action potentials, glycine receptors, and AMPA-glutamate receptors, respectively. GABAAR minis were recorded from a holding voltage of -60 mV. The Nernst potential for Cl under our experimental conditions was ~ 5 mV; thus, inward currents in the presence of strychnine and DNQX represent minis that arise from GABAARs. This was confirmed in each experiment by applying bicuculline (bic) follow each recording. We did not observe changes in the mIPSC amplitude (p = 0.5218; unpaired t test), charge transfer (p = 0.4391; unpaired t test), rise time (p = 0.1029; unpaired t test), frequency (p = 0.5576; unpaired t test), as well as neuronal input resistance after hibernation (p = p = 0.3647; unpaired t test) (Fig. 3A-C). Therefore, hibernation does not influence postsynaptic GABAARs and is unlikely to account for the loss of GABAAR signaling at the network level.

Fig. 3
figure 3

Postsynaptic GABAergic currents in motoneurons are unchanged following hibernation. A Example voltage clamp traces of GABAergic miniature inhibitory postsynaptic currents (mIPSCs; holding potential = -60 mV; 250 nM TTX, 10 µM DNQX, 5 µM strychnine) from control (black) and hibernated (blue) identified vagal motor neurons in brain slice (middle right box). Addition of 50 µM bicuculline (Bic) to the bath eliminated mIPSCs in both groups, confirming GABAergic identity. B Cumulative probability histograms of mIPSC amplitudes from control neurons (black line), hibernated neurons (blue line) appeared nearly identical. Accordingly, there was no statistical difference between distributions from control and hibernated cells by Kolomogorov-Smirnov test (p = 0.5041). C Summary data. There was no statistical difference in mIPSC frequency, average amplitude, rise time, charge transfer, and neuronal input resistance by unpaired t-test (p > 0.05; n = 16 control cells; n = 16 hibernated cells). Error bars are standard deviation (SD)

As the quantal GABAAR currents were unchanged on the postsynaptic motoneuron, we assessed if loss of the GABAergic tone is associated with reduced network-driven presynaptic input onto motoneurons, like that seen in sensory deprivation models of homeostatic plasticity [10, 11]. Specific respiratory synapses are not yet tractable for stimulation experiments in this system; therefore, we developed an approach to circumvent this issue. We used a novel semi-intact preparation that permits recording of individual motoneurons receiving excitation and inhibition from rhythmic premotor circuits that shape motoneuron firing of breathing (Fig. 4A) [34]. As shown in Fig. 4A, vagal motor neurons fire bursts of action potentials during respiratory population motor activity, reflecting activity of breathing at a single-cell resolution. Since the postsynaptic receptor density of GABAARs does not appear to change within these neurons after hibernation (Fig. 3), any difference in the role of GABAergic synapses for the control of motoneuron firing rate would reflect changes in network-driven presynaptic GABA release. To isolate GABAARs that control respiratory-related firing of motoneurons, we focally applied bicuculline to the motoneuron cell body and surrounding region via local pressure injection while recording its activity shaped by excitation and inhibition from rhythmic premotor inputs in the intact network. We simultaneously recorded rhythmic vagal motoneuron activity and the trigeminal nerve (a nerve that innerves the buccal floor to drive air into the lungs and activates near-synchronously with the vagal outflow) to monitor respiratory circuit output. We used the trigeminal nerve (cranial nerve V; CN V) in these experiments because physical constraints in the bath and the focal drug delivery system prevented simultaneous recording of the vagus nerve.

Fig. 4
figure 4

Hibernation decreases the influence of presynaptic GABA on motoneuron firing. A Integrated network output recorded on the trigeminal nerve (CN V) and simultaneous whole-cell current clamp recordings of vagal motor neurons (CN X MN) driven by the respiratory circuit in the semi-intact brainstem preparation (top left box). A pipette was used to focally apply bicuculline onto the recorded vagal motor neuron to isolate the GABAergic tone on individual neurons. An example respiratory neuron (CN X MN, gray, top right) fired rhythmic bursts, characterized by transient increases in firing frequency (Firing Freq, light gray, middle top right) during respiratory motor output (Black, bottom top right, CN V). B Example respiratory motor activity in a control preparation (black). Example vagal motor neuron (CN X MN) fires bursts of action potentials during the respiratory burst (CN V). Focal application of bicuculline (gray box) to the motor neuron increases neuronal burst firing frequency (Firing Freq). C Example respiratory motor activity in a preparation from a hibernated animal (blue). Example vagal motor neuron (CN X MN) also fired bursts of action potentials during the respiratory burst (CN V). Focal application of bicuculline (gray box) to the motor neuron increased neuronal burst firing frequency (Firing Freq, Δ ~ 47 Hz), but not to the same extent as the control neuron (Firing Freq, Δ ~ 74 Hz). D Summary of bicuculline-mediated increases in motor neuron firing frequency. There was a significant interaction between drug and group by two-way ANOVA (F(2,38) = 6.391; p = 0.0041). Furthermore, there was a significant increase in burst firing frequency following focal application of bicuculline in control but not hibernated cells by Holm-Sidak’s post-hoc test (***p < 0.0001; n = 12 control cells; n = 9 hibernated cells). E Summary of bicuculline-mediated fold change in burst firing frequency from baseline. Control cells had a significantly larger bicuculline-mediated fold change in burst firing frequency than hibernated cells by unpaired t-test (p = 0.036; n = 12 control cells; n = 9 hibernated cells). Orange cells were recorded in loose-patch configuration. Box plots represent interquartile range and whiskers represent the minimum to maximum values. a.u., arbitrary units. Error bars are standard deviation (SD)

Motoneurons from both groups had similar average burst firing rates during breathing and resting membrane potentials in between respiratory bursts (firing rate: control, 45.26 ± 34.8 Hz vs. hibernation, 24.19 ± 23.0 Hz; p = 0.133, unpaired t-test; membrane potential: control, -61.52 ± 4.6 mV vs. hibernation, -59.24 ± 6.5 mV, p = 0.4094 unpaired t-test). In controls, local application of bicuculline led to a reversible increase in average firing frequency during the respiratory burst (Baseline, 45.26 ± 34.8 Hz; Bicuculline, 89.36 ± 61.4 Hz; p < 0.0001; Holm-Sidak multiple comparisons test; n = 12; Fig. 4B,D). This was presumably not due to changes in the chloride gradient, because two neurons in the control group that were recorded in the “loose patch” configuration, which does not disturb intracellular milieu (Fig. 4E shown in orange), more than doubled their firing rates during bicuculline exposure, aligning with data from the whole-cell mode. These results demonstrate that phasic inhibition onto motoneurons dampens firing rate during the respiratory burst. In contrast, local bicuculline had no significant effect on the firing rate of motoneurons from hibernators (Baseline, 24.19 ± 23.0 Hz; Bicuculline, 38.75 ± 44.9 Hz; p = 0.2629; n = 9; Holm-Sidak multiple comparisons test; Fig. 4C-E). Additionally, there was no change in the interburst (resting) membrane potential during focal application of bicuculline in either control (Baseline: -61.52 ± 4.6 mV; Bicuculline: -60.92 ± 4.9 mV; n = 10) or hibernation groups (Baseline: -59.24 ± 6.5 mV; Bicuculline: -59.53 ± 6.5 mV; n = 7), further supporting that phasic respiratory-related inhibition influences motoneuron firing rate rather than a tonic GABAAR current. Taken together with the mIPSC results, network-driven GABA release appears to be reduced in a way that no longer dampens the firing rate of motoneurons during activity associated with breathing.

These results demonstrate a critical role for GABAARs in generating and modulating the respiratory rhythm, as well as controlling motoneuron firing rate in control frogs. Therefore, we were surprised that baseline network frequency on average (Fig. 1C, p = 0.2021) and motor burst morphology (Fig. 5) in hibernators did not change despite a rather dramatic loss of GABAAR signaling in the network. These results mirror ventilation data in vivo, whereby breathing frequency and breath volume were the same after hibernation [35]. Despite these similarities at warm temperatures, frogs likely need to restart breathing to some degree at temperatures around 8 °C to maintain aerobic metabolism [35, 36] and then produce reliable rhythmic output at temperatures above 13 °C to maintain metabolic homeostasis as they reestablish a warmer life on land [25, 26]. Yet, cold temperatures ≤ 15 °C strongly depress the network, making it difficult to generate respiratory activity [15, 37]. Interestingly, cooling does not suppress breathing exclusively through “passive” effects on cellular properties (e.g., slower action potentials, decreased rates of synaptic transmission), but rather, through noradrenergic signaling. Indeed, blocking α adrenergic receptors blunts decreases in activity during cooling, and the actions of norepinephrine on respiratory activity in the frog brainstem are known to act through GABAergic pathways [38, 39]. This led us to test whether decreased GABAAR signaling enhances respiratory motor output at cold temperatures faced by frogs when they need to resuscitate breathing while emerging from hibernation.

Fig. 5
figure 5

Respiratory burst morphology after hibernation closely matches controls at room temperature. Integrated vagus motor bursts (CN X) in preparations from 5 controls (red) and 5 hibernators (blue). Overlayed bursts (n = 50) from individual preparations are represented in different shades from top to bottom. Overlaying all bursts from all animals from both groups demonstrates a close matching between controls and hibernators, indicating the loss of GABAA receptors has no obvious impact on motor burst morphology after hibernation

To address this, we use three groups of ex vivo brainstem preparations: controls, hibernators (reduced GABA signaling), and controls with GABAARs experimentally reduced with a subsaturating dose of bicuculline (2 μM). We hypothesized that control preparations with impaired GABAAR signaling would have similar responses to cooling that were like hibernators if GABAAR signaling is part of the process by which cold temperatures depress activity of the respiratory network. That is, less GABAAR signaling would cause enhanced activity at colder temperatures. First, we addressed the respiratory frequency. Cooling control brainstems from 20 °C to 8 °C expectedly reduced the respiratory burst frequency (Fig. 6A). Interestingly, the decline in hibernators was less pronounced, with a temperature that produces 50% of baseline burst rate shifting to colder temperatures (Fig. 6B). Decreasing GABAAR signaling in controls via subsaturating block of GABAARs also produced greater burst frequency at cold temperatures; however, this was more dramatic than hibernators, resulting in some networks with activity at 8 °C (Fig. 6C). Summary statistics are shown in Fig. 6D. Two-way ANOVA revealed a significant temperature, group, and temperature by group interaction effect. These results were driven by several pairwise differences caused by hibernators and controls with bicuculline. Hibernators had statistically greater normalized burst frequency than controls at cooler temperatures (16 °C: p = 0.0267; Holm-Sidak Multiple Comparisons test, Fig. 6D blue). In addition, bicuculline application to controls had significantly greater normalized burst frequency at lower temperatures (Fig. 6D gray). These results are also summarized as the temperature at which burst frequency decreased by 50% (T50) (Fig. 6E). Indeed, hibernators had lower T50 than controls, and controls with bicuculline had the lowest T50 values. Taken together, the network is more active at cooler temperatures when GABAAR signaling is reduced naturally by hibernation or by experimentally blocking GABAARs.

Fig. 6
figure 6

Lowering GABAAR signaling enhances motor performance of breathing at cool temperatures below baseline. Integrated vagus nerve activity (CN X) from intact brainstem preparations (bottom right box) showing the response of control (A, black), hibernators (B, blue), and controls treated with 2 μM bicuculline (C, gray) during cooling ramps from 20°C to 8°C. Zoomed-in traces below the compressed recordings illustrate the effects of cooling on burst amplitude. Arrows indicate the temperature at 50% of baseline burst frequency. D Mean data showing average temperature responses of normalized burst frequency among groups (two-way ANOVA interaction of group x temperature, p<0.0001). E Mean data showing a significant difference in the temperature at 50% of baseline burst rate between control, hibernators, and controls+bicuculline (one-way ANOVA; p<0.0001) (F) Mean data showing change in burst amplitude from baseline in control and hibernators during temperature changes. Controls undergo a decrease in burst amplitude during cooling, while hibernators and controls+bicuculline maintain burst amplitude throughout the cooling ramp (two-way ANOVA, p=0.0001 group x temperature interaction); * represents p<0.05, **<0.01, ***<0.001 in Holm-Sidak multiple comparisons test following one or two-way ANOVAs. n=5 for control and bicuculline, and n=6 for hibernators

In addition to respiratory frequency, we addressed the amplitude of the motor output, as this variable relates to activation of the respiratory muscles. Like frequency, cooling also depressed motor amplitude in controls, reflecting lowered motoneuron firing or recruitment during the fictive breath (Fig. 6A). In hibernators and controls in the presence of bicuculline, burst amplitude appeared to be maintained across the full temperature range (Fig. 6B-C). Summary statistics are shown in Fig. 6F. Two-way ANOVA revealed a significant temperature, group, and temperature by group interaction effect. These results were driven by several pairwise differences in hibernators and controls with impaired GABAA signaling. Hibernators and controls with bicuculline had statistically greater burst amplitude than controls at 12 °C and 10 °C (12 °C: hibernator vs. control: p = 0.0276; bic vs. control: p = 0.0039, 10 °C: hibernator vs. control: p = 0.0254, bic vs. control: p = 0.0254) Holm-Sidak Multiple Comparisons test). Overall, while network output at warm temperatures were similar between controls and hibernators with reduced GABA signaling, these results demonstrate that the downregulation of GABAAR signaling by hibernation drives more robust respiratory motor activity (frequency and motor amplitude) at cooler temperatures.

Discussion

Neural circuits face constant disturbances from the environment, placing a strong pressure on plasticity mechanisms to regulate activity. Here, we corroborate previous work suggesting GABAARs play a role in generating the respiratory motor pattern in adult frogs and identified that they also play a distinct modulatory role to depress activity during cooling. We found that the state of hibernation triggers a downregulation of synaptic inhibition which would have been expected to cause general alterations in network activity. Instead, activity of the respiratory network became more resistant to decreases in activity caused by cooling in a way that would help stabilize respiratory output as the brainstem passes through a large temperature range during emergence from the hibernation environment. Therefore, compensatory forms of synaptic plasticity can expand the operating range of circuits in challenging environments without obvious impacts on activity in unstressed conditions.

Reduced GABAergic inhibition at rhythm generating and premotor loci

Given that the respiratory motor system undergoes a massive activity challenge (inactivity during winter [15]), mechanisms consistent with activity-dependent compensation have been predicted to play a role in maintaining network integrity to promote motor behavior in the spring [16, 17, 35]. In addition to upregulation of excitatory synaptic strength on motoneurons, here we found that hibernation strongly reduced GABAAR signaling throughout the respiratory motor network.

The present data lead us to suggest GABAAR plasticity at two loci within the network (Fig. 7; top panel labeled “warm”): rhythm generating circuits and motoneurons. First, expression of the respiratory rhythm in this species involves GABAARs [24], likely through mechanisms involving reciprocal inhibition. GABAARs in premotor rhythm-generating circuits must have a diminished contribution after hibernation because motor frequency in hibernators was strongly resistant to depression caused by GABAAR antagonists or agonists. Although counterintuitive, these opposing pharmacological manipulations (activating and inhibiting the receptors) both lead to suppression of the rhythm. This is likely to occur because too much inhibition depresses activity by hyperpolarizing neurons, and too little prevents phasic inhibition required for rhythmogenesis. This is reminiscent of what occurs in other species where rhythm generation depends on synaptic inhibition [32, 33]. Second, activity of the respiratory rhythm generator is transmitted through excitatory and inhibitory interneuronal pathways to recruit motoneurons that drive breathing. In contrast to inhibition’s presumed role in facilitating rhythm generation, inhibition onto motoneurons has a conventional role of constraining the firing rate during respiratory-related bursting, which is reduced after hibernation (Fig. 4). As mIPSCs carried by GABAARs on motoneurons were unchanged (Fig. 3), differences in presynaptic properties of inhibitory interneurons within the respiratory network likely account for weakened GABAergic control of firing (Fig. 4). Mechanistically, it is difficult at present to parse out the variables that account for the apparent reduction in network-driven presynaptic GABA (e.g., reduced vesicle release probability, decreased excitability of GABAergic interneurons, etc.), as respiratory-related synapses and GABAergic interneurons have not yet been identified in this network. Of note, rhythm and pattern generation are thought to arise, in part, from distinct populations within the respiratory central pattern generator of mammals [40]. As such, the differential effects we observed involving loss of inhibition in the respiratory circuit may stem from distinct populations within the central pattern generator itself. Overall, the specific mechanisms and network interactions await further experimentation, but the sum of our data points to the idea that rhythm-generating circuits and premotor neurons rely significantly less on GABAARs after hibernation.

Fig. 7
figure 7

Schematic illustrating how reducing inhibition may promote breathing in the cold. In controls at warm temperatures (orange left), inhibition (red, minus symbol) plays an important role in generating the respiratory rhythm (circular arrows) and controlling the excitability to motoneurons. Following hibernation (orange right), the respiratory network decreases reliance on inhibition and potentially increases reliance on mechanisms that involve excitation (dark blue, plus symbol) for respiratory rhythmogenesis (circular arrows). Nonetheless, there is no discernable change in network motor output at warm temperatures following hibernation compared to controls. During acute cooling, inhibitory neurotransmission also plays a modulatory role and is increased. Whether inhibitory neurons that drive the depression in activity during the cold are part of the rhythm generating circuits (as drawn) or distinct modulatory neurons is not yet known. As a result, controls (blue left) have marked reduction in respiratory frequency and motor output at colder temperatures due to elevated GABAergic inhibition. In contrast, output following hibernation is resistant to depressive effects of the cold, likely with a contribution resulting from the decreased reliance on inhibition that normally slows activity and depresses motor output. Thus, output following hibernation is more robust at cooler temperatures, in part, because inhibitory neurotransmission is decreased within rhythm-generating and interneuronal pathways

Given the profound loss of GABAAR signaling, which is normally required to express the respiratory rhythm, one implication of these results is that some other mechanism must take over to produce breathing after hibernation. These results align well with recent work demonstrating circuits that recently recovered from an activity challenge use a different profile of ion channels to generate similar network output under baseline conditions [20, 21]. This leads us to speculate that plasticity resulting from hibernation may switch the mechanism of rhythmogenesis from one which relies heavily on inhibition to others that use excitatory synapses or intrinsic pacemaker mechanisms [24, 41] (Fig. 7; warm panel highlight shift in rhythm generating mechanisms, circular arrows). Consistent with this idea, the amphibian respiratory network appears to use voltage-dependent pacemaker mechanisms early in development, which then is replaced by network-dependent mechanisms via inhibition following metamorphosis [24]. Therefore, hibernation may revert the mechanisms of rhythm generation to that seen in an earlier life stage. Indeed, other aspects of breathing such as modulation of frequency by respiratory gases also revert to a juvenile-like state following hibernation [42], consistent with the idea that emerging from aquatic hibernation may share similar environmental pressures as developing air-breathing during the shift from aquatic to terrestrial habitats [43].

Reduced GABAAR signaling increases the network operating range in the cold

The most intriguing aspect of these results was that networks operating with less GABAAR signaling after hibernation seemed to perform better at cooler temperatures (Fig. 6), despite the fact that activity was largely unchanged at warm temperatures (Fig. 1&5). Our data indicate that GABAAR signaling is part of the process by which cold normally depresses the frequency and amplitude of the motor output in control animals (Fig. 7; bottom panel labeled “acute cooling”). Although temperature effects on neural circuits are multifaceted, we draw this conclusion because subsaturating block of GABAARs in control animals opposed the typical depressive response to cooling. The mechanisms by which GABAAR signaling normally decreases activity in the cold seems to occur through two different processes. First, at the motoneuron level, cold may facilitate inhibitory GABAergic input that we showed to exist on these neurons (Fig. 4) in a way that depresses the population motoneuron output (Fig. 6). Thus, losing this inhibition during hibernation likely explains maintenance of the burst amplitudes across the full temperature range during acute cooling (Fig. 7, bottom right, less inhibitory input in the cold). Second, our pharmacological data suggest that GABAAR signaling plays a role in slowing breathing frequency during cooling. Therefore, the loss of this inhibition in hibernators accounts for the ability to produce stable rhythmic activity at cooler temperatures (Fig. 6E). The network-level mechanisms remain to be uncovered; however, cold temperatures activate noradrenergic neurons of the locus coeruleus [44, 45], which then depress breathing frequency through α-adrenoreceptors [39]. Since noradrenergic signaling depresses respiratory activity mainly through GABAergic pathways [38], modulatory mechanisms involving norepinephrine may be activated by cooling to recruit GABAergic interneurons that depress the breathing frequency at cold temperatures.

Regardless of the specific mechanism by which cold temperature and GABAARs interact, reduced GABAAR signaling in hibernators, therefore, serves to boost respiratory frequency and motor drive at cooler temperatures, with likely relevance for restarting breathing upon exit from hibernation as the animal passes through a range of cool temperatures (Fig. 7; bottom panel labeled “cold”). As such, we interpret plasticity of GABAAR signaling here as “homeostatic,” as it is a change in network function that opposes the activity disturbance and enhances network excitability in an environment that otherwise depresses neuronal output. The loss of GABAAR signaling after hibernation seems to explain most of the maintenance of motor amplitude across temperatures, since hibernators and controls with experimentally reduced GABAAR signaling had similar profiles across temperatures. However, despite hibernators having higher activity at lower temperatures than naïve controls, activity still dropped to a lower extent than we observed in controls with GABAAR signaling experimentally reduced (Fig. 6A-D). We speculate the quantitative differences between these two conditions might point to additional plasticity mechanisms that tamp down enhanced excitability at the coldest temperatures to prevent ectopic breaths while submerged in cold water. While the mechanisms that serve this purpose are not fully known, we have shown that sensory input from CO2/pH chemoreceptors that normally stimulate breathing is reduced in hibernators [35, 42], suggesting that hibernation may also reduce excitability in some parts of the network to help keep breathing off until it is needed. The exact “switch” that restarts breathing also remains unclear, but nevertheless, our results suggest that the loss of GABAAR signaling prepares the network to produce strong motor bursting at cool temperatures when it does restart, which would likely contribute to restoring breathing as the animal passes through a range of colder temperature before returning to a warmer body temperature after emergence.

Conclusions

What does it mean for plasticity to be “homeostatic?” Since its inception [6], compensatory mechanisms within the theme of homeostatic regulation are most often portrayed as a conceptually simple (but mechanistically complex) counterbalancing of altered excitability. Our results provide evidence that mechanisms consistent with homeostatic plasticity can increase the operating range of circuits in harsh environments that otherwise disrupt activity. We acknowledge that a strict definition of homeostatic plasticity might not encompass the plasticity we observed. Through a homeostatic lens, one might expect network activity that drives lung ventilation to restart as compensation enhances network excitability. However, breathing does not restart during hibernation [15], which in itself seems important for animal survival as restarting breathing prematurely could be detrimental underwater. In addition, the integrated output of the entire respiratory network does not appear obviously hyperexcitable after emergence, as breathing in vivo is comparable between controls and hibernators, which could be interpreted as inconsistent with homeostatic theory [35]. While our observations deviate from traditional definitions of negative feedback homeostasis on which the field was founded [4, 6, 46], reduced GABAAR signaling indeed increases the strength and frequency of respiratory motor activity in cool conditions that must be passed through to reestablish life on land after hibernation, thereby widening the circuit’s operating range across temperatures. Intriguingly, hibernation also leads to enhanced desensitization and reduced Ca2+ permeability of NMDA receptors that serves to constrain excitability of this network during severe hypoxia that occurs as the animal emerges, with no impact on baseline activity of the network [47]. Therefore, multiple forms of synaptic plasticity appear to be masked at baseline, but ultimately increase the operating range for motor activity under the variety of stressors that occur during emergence from the hibernation environment. Overall, these results support the view that circuits may implement compensatory mechanisms in ways that enhance activity under specific sets of environmental constraints. These findings may have important implications for diverse neural systems, as the nervous system of most animals operate homeostatically over a range of abiotic conditions, such as changes in local tissue temperature [48], pH [49, 50], oxygen levels [47], hormonal state, and other factors that have the potential to disrupt network activity on acute and chronic timescales. By linking plasticity to a challenging environment in vivo, our results lead us to emphasize that frameworks seeking to link compensation to behavior must ultimately incorporate how plasticity is integrated over the range of environments that neural circuits may encounter. Understanding how the nervous system is tuned to survive extreme conditions in intact animals will continue to provide new insights along this path.

Methods

Animals

All experiments performed were approved by the Animal Care and Use Committee (ACUC) at the University of Missouri (protocol #39,264). Adult Female American Bullfrogs (~ 100 g weight, n = 62), Lithobates catesbeianus, were purchased from Rana Ranch (Twin Falls, Idaho) and were housed in 20-gallon plastic tanks containing dechlorinated water at room temperature bubbled with room air. Control frogs were acclimated for 1 week following arrival before experiments. Frogs had access to wet and dry areas. Frogs were maintained on a 12-h light/dark cycle and fed once per week. Water was cleaned daily and changed weekly. Hibernated frogs were kept in plastic tanks under the same conditions (exception of food being withheld) as control frogs for > 1 week before temperature was gradually lowered to 4 °C over 10 days in a walk-in temperature-controlled environmental chamber. Once water temperature reached 4 °C, air access was blocked using a plastic screen placed in the tank. After 4 weeks of submergence, experiments commenced.

Drugs

Strychnine hydrochloride was from Sigma-Aldrich (St Louis, MO, USA). (-)-Bicuculline methiodide, tetrodotoxin citrate, DNQX disodium salt, and Muscimol were from Hello Bio (Princeton, NJ, USA).

Brainstem-spinal cord preparation

Brainstem-spinal cord preparations were generated as previously described [16]. Briefly, frogs were deeply anesthetized with isoflurane until visible respirations had ceased and response to foot pinch was absent. Frogs were then decapitated with a guillotine. The head was submerged in ice-cold bullfrog artificial cerebrospinal fluid (aCSF; concentrations in [mM]: 104 NaCl, 4 KCl, 1.4 MgCl2, 7.5 glucose, 40 NaHCO3, 2.5 CaCl2 and 1 NaH2PO4, and bubbled with 98.5% O2, 1.5% CO2; pH = 7.85) and the forebrain was pithed. The brainstem-spinal cord was then carefully removed, keeping nerve roots intact. Following the dissection, the brainstem was transferred to a chamber superfused with oxygenated aCSF. In a subset of experiments where bicuculline was bath applied to the in vitro preparation (described below), brainstems were transected at the spinomedullary junction before recording to avoid the contribution of spinal motor populations to non-respiratory motor activity induced through disinhibition [51]. Cranial nerve X (vagus) activity was recorded using glass suction electrodes immediately following dissection. Recordings were AC amplified (1000x, A–M Systems Model 1700, A–M Systems, Carlsborg, WA, USA), filtered (10 Hz – 5 kHz), and digitized (Powerlab 8/35 ADInstruments, Colorado Springs, CO, USA). Nerve activity was rectified and integrated (time constant = 100 ms) online to allow for visualization of nerve output pattern.

Vagus nerve output was then allowed to stabilize for ≥ 1 h before drug application. In a subset of experiments (n = 10 bullfrogs, 5 per group), various does of bicuculline (500 nM, 1 µM, 5 µM, 10 µM), a GABAA receptor antagonist, were superperfused for 10 min in order of increasing concentration and then washed out. In a different subset of experiments (n = 10 bullfrogs, 5 per group), various doses (500 nM, 1 µM, 3 µM) of muscimol, a GABAA receptor agonist, were added to the perfusate. Each dose was superperfused for 10 min in order of increasing concentration. Bicuculline (10 µM) was added to the perfusate following muscimol to reverse muscimol’s effects serving as a control that output depression was not due to some kind of drift over time, but rather activation of GABAA receptors.

For temperature experiments, preparations were allowed to recover for ~ 1 h at 22 °C before exposure to the temperature ramp from 20° to 8 °C. Preparations from controls (n = 5) and hibernators (n = 6) were cooled to each temperature for 10 min and then the temperature was lowered by 2 °C. Following cooling, preparations were rewarmed on a continuous ramp that lasted approximately 5 min. Separate control preparations (n = 5) were treated with 2 μM bicuculine for 30 min prior to starting the cooling ramp. 2 μM was selected as it was the highest dose where respiratory-related bursts were prominent in controls. The temperature was controlled by a Warner Instruments bipolar in-line temperature controller (model CL-100; Hamden, CT), and solution temperature within the chamber was monitored directly next to the brainstem.

Dye labeling of vagal motor neurons

Brainstems were isolated as described above. In a different group of experiments, following isolation, brainstems were transferred to a chamber superfused with oxygenated aCSF. The 4th branch of the vagal nerve root, which primarily innervates the glottal dilator [52] was then isolated and suctioned into a glass pipette (bilaterally). Fluorescent dextran dye dissolved in PBS (tetramethylrhodamine 3000 MW lysine fixable dye; Invitrogen – Thermo Fisher, Waltham, MA, USA) was then backfilled into both pipettes and left to diffuse into the nerve root for 2 h. This generated robust labeling of vagal motor neurons as previously described [16].

Semi-intact brainstem-spinal cord preparation

Following dye labeling, in a subset of experiments (n = 15 bullfrogs) the brainstem was embedded into agarose and sliced, similar to previously described [33] with some modifications. Briefly, partially cooled, but not solidified, agarose was pipetted into a thin layer onto a scored (with fine forceps) agar block. The brainstem was then dragged onto the agarose-covered agar block, laying horizontally dorsal side up. Agarose (~ 57˚ C) was pipetted onto the top of the brainstem caudal to obex. The agar block was then mounted on a vibrating microtome plate (Campden Vibrating Microtome 7000smz, Campden Instruments; Lafayette, IN, USA) and covered in ice-cold bubbled aCSF. The blade was zeroed at the top of the brainstem, then slices were taken caudal to rostral, stopping shortly rostral to obex which was approximately 1/3 of the distance from the hypoglossal root to the vagal root and spanned a fraction of the vagal motor pool extent. The blade was removed after each slice and the resulting attached slice was trimmed off the brainstem with fine spring scissors. Approximately 650 μm were removed from the dorsal surface in total, providing access to labeled motor neurons for recording. The brainstem was removed from the remaining agar and transferred to a 35 mm glass bottom sylgard-coated dish with a stainless-steel mesh insert. The brainstem was then pinned dorsal side up. The dish was placed in a QE-1 platform (Warner Instruments, Holliston, MA, USA) that was mounted on a fixed stage microscope (FN1, Nikon Instruments Inc., Melville, NY, USA). Once on the rig, the preparation was superperfused with bubbled aCSF either using gravity or with a peristaltic pump (Rainin Rabbit). The trigeminal nerve (CN V) was recorded using a glass suction electrode to monitor preparation output over time. Nerve recordings were AC amplified (1000x, A–M Systems Model 1700, A–M Systems, Carlsborg, WA, USA), filtered (10 Hz – 5 kHz), and digitized (Powerlab 8/35 ADInstruments, Colorado Springs, CO, USA). The trigeminal nerve was recorded instead of the vagal nerve root because of physical constraints in the bath with the patch pipette and picospritzer pipette used in experiments described below. The trigeminal provides a sufficient surrogate for efferent respiratory motor output [29, 53], as respiratory activity of the trigeminal and vagus nerves activate near-synchronously in vitro. Labeled motor neurons were then located and visualized using an imaging camera (Hamamatsu ORCA Flash 4.0LT sCMOS, Hamamatsu Photonics, Hamamatsu City, Japan) coupled to Nikon imaging software (NIS elements).

Whole-cell recordings were made from labeled vagal motor neurons with an Axopatch 200B amplifier (Molecular Devices) in current-clamp mode (I = 0). Cells were selected based on firing pattern, that is only cells that fired spontaneous bursts of action potentials or increased firing frequency (respiratory modulated) during respiratory motor bursts were sampled, as these neurons likely innervate the glottal dilator of bullfrogs [54], which is critical to permit airflow into the lungs of bullfrogs during ventilation and as such remains inactive during hibernation. Accordingly, neurons with subthreshold respiratory input or neurons that were active only during the interburst interval were excluded from testing. Glass pipettes (~ 4–7 MΩ) were filled with a solution containing (in mmol): 110 potassium gluconate, 2 MgCl2, 10 Hepes, 1 Na2-ATP, 0.1 Na2-GTP and 2.5 EGTA, pH 7.2 with KOH. Data were acquired in pClamp 11 software using an Axopatch 200B amplifier and Axon Digidata 1550B digitizer (Molecular Devices, San Jose, CA, USA). To determine inhibitory tone onto respiratory motor neurons, bicuculline (50 µM) was focally applied using a Picospritzer II (General Valve Corporation, Fairfield, NJ, USA) to the recorded cell following a baseline period (≥ 10 respiratory bursts). Bicuculline was applied using a 5 s pulse, every 10 secs for 5–10 min with the picospritzer pipette slightly larger (~ 2 μm diameter) than typical for focal application and placed approximately 15 μm away from the cell body. This approach was taken to saturate the cell and surrounding area with the GABAAR antagonist. Importantly, this approach did not recapitulate systemic effects of bicuculline on the vagal motor pool as the population spans nearly 2.5 mm in the rostrocaudal axis in ranid frogs [55] and only a subset of the vagal motor pool was exposed in the preparation as described above. A washout period was used to verify changes in cell firing frequency were due to drug application and not cell drift. Most patched cells were near the surface of the tissue as the thickness of the preparation hindered cell visualization. As a result, changes in firing frequency during bicuculline application and wash were rapid, suggesting the area was saturated with drug upon administration and then quickly washed. In a subset of cells where a gigaohm seal was not formed, recordings of actional potentials were performed in the loose-patch configuration (Vhold = 0 mV; n = 2 in control, and n = 2 in hibernation). Regardless of recording configuration, firing rate data from recorded neurons were combined in group analysis of firing frequency, as data from loose-patch aligned with data from whole cell configuration. All voltages from current-clamp experiments were corrected for a liquid junction potential of 12 mV.

Brain slice electrophysiology

Following dye labeling, in a subset of experiments brainstem slices (300 µm) containing labeled motor neurons were obtained (n = 11 bullfrogs). Whole-cell recordings were made from labeled motor neurons with an Axopatch 200B amplifier (Molecular Devices) in voltage clamp (Vhold = -60 mV) mode. Glass pipettes (2.6–4 MΩ) were filled with a solution containing (in mM): 95 CsCl, 2 MgCl2, 10 Hepes, 1 Na2-ATP, 0.1 Na2-GTP, 10 EGTA, 1 CaCl2 and 10 tetraethylammonium-Cl (TEA), pH ~ 7.2 using CsOH. Data were acquired in pClamp 11 software using an Axopatch 200B amplifier and Axon Digidata 1550B digitizer (Molecular Devices, San Jose, CA, USA). To examine miniature inhibitory postsynaptic currents (mIPSCs), tetrodotoxin (TTX; 250 nM) and DNQX (10 µM) were bath applied for 3 min to block voltage-gated Na+ channels and AMPA receptors. Strychnine (5 µM) was then bath applied in combination with TTX and DNQX for 2 min to block glycinergic mIPSCs and isolate GABAergic mIPSCs which were recorded for 1 min. All GABAergic mIPSCs were blocked following perfusion of bicuculline (50 µM) confirming their isolation. A 10 mV step was performed at the beginning of each recording and before each solution change to monitor series resistance (Rs; estimated from the peak of the capacitive transient) and input resistance (RIN; estimated from the steady-state current). Rs was not compensated but remained < 25 MΩ for inclusion.

Data analysis

For experiments using the intact brainstem preparation, respiratory burst frequency was determined from integrated vagus nerve signals using the peak analysis function in LabChart 8 (ADInstruments Inc., Colorado Springs, CO, USA). Respiratory bursts were identified based on standard metrics used in the field such as ~ 1 s duration [56, 57]. Burst start and stop time points were defined as 5% of the height from baseline. In experiments where muscimol was systemically applied, data were sampled from the last 3 min of each condition (baseline, 500 nM, 1 µM, and 3 µM). In experiments where bicuculline was systemically applied, data were sampled for 10 min during each condition (baseline, 500 nM, 1 µM, 5 µM, 10 µM). In systemic bicuculline experiments where non-respiratory bursting was quantified, non-respiratory were identified by their long duration and qualitatively different shape than respiratory bursts and were counted manually. Non-respiratory bursts appeared in a variety of output patterns as previously characterized in vitro (Reid and Milsom, 1998). Non-respiratory bursts were not quantified during systemic muscimol exposure as they were largely absent from baseline recordings and only emerged with notable frequency following disinhibition with bicuculline. For temperature experiments, burst frequency was averaged in the last 5 min at each temperature, and burst amplitudes were sampled in the last minute or in the case of the failure temperature, bursts were sampled directly before failure. To compare respiratory burst amplitudes in the cold, we analyzed the change from baseline at the coldest temperature the preparation could produce activity. Across control and hibernation groups this typically occurred at 8–10 °C.

For experiments using the semi-intact preparation, action potentials were detected with the peak analysis function in LabChart 8 (ADInstruments Inc., Colorado Springs, CO, USA). Following detection of action potentials, a channel was created that plotted firing frequency over time. Average firing frequency per respiratory burst was determined via the mean firing frequency from neuronal burst start to stop. Average firing frequency was sampled and averaged from 10 respiratory bursts in each condition (baseline, bicuculline, and wash).

For voltage-clamp experiments in brainstem slices, average amplitude (current measurement from baseline to peak), charge transfer (integral of the mIPSC), rise time (time from baseline to peak), and frequency of GABAergic mIPSCs (mIPSC/sec) were analyzed from one minute of gap-free recording following two minutes of exposure to TTX/DNQX/Strychnine (as described above) using the peak analysis function in LabChart 8 (ADInstruments Inc., Colorado Springs, CO, USA). Events below 7.5 pA were excluded. Events were inspected manually to ensure the accurate detection of mIPSCs. For construction of the cumulative probability histogram, the first 50 mIPSCs from the minute of data obtained from each neuron were sampled.

Statistics

Data are presented as mean ± s.d. unless otherwise stated or shown as individual points and means to highlight individual responses. When two groups of dependent samples were compared (“before-after” experiments), a two-tailed paired t-test was used. When two groups of independent samples were compared, a two-tailed unpaired t-test was used. Non-parametric versions of these tests were used if data sets were not normally distributed. In experiments with one main effect, a one-way ANOVA was used to test for the main effect and the interaction effect. One-way ANOVA was followed up with Holm-Sidak multiple comparisons test. In experiments with two main effects, a two-way ANOVA was used to test for the two main effects and the interaction effect. Two-way ANOVA was followed up with Holm-Sidak multiple comparisons test unless otherwise specified. Cumulative distributions were compared with the Kolmogorov–Smirnov test. Significance was accepted when P < 0.05. All analyses were performed using GraphPad Prism (v9.4.1, San Diego, CA, USA).

Data availability

All data generated or analysed during this study are included in this published article and its supplementary information files.

Abbreviations

GABAAR:

GABAAR Receptor

CNX:

Cranial nerve X

strych:

Strychnine

mIPSC:

Miniature inhibitory postsynaptic currents

TTX:

Tetrodotoxin

CNV:

Cranial nerve V

References

  1. Marom S, Marder E. A biophysical perspective on the resilience of neuronal excitability across timescales. Nat Rev Neurosci. 2023;24(10):640–52.

    Article  PubMed  CAS  Google Scholar 

  2. Desai NS, Rutherford LC, Turrigiano GG. Plasticity in the intrinsic excitability of cortical pyramidal neurons. Nat Neurosci. 1999;2(6):515–20.

    Article  PubMed  CAS  Google Scholar 

  3. Kilman V, Van Rossum MC, Turrigiano GG. Activity deprivation reduces miniature IPSC amplitude by decreasing the number of postsynaptic GABAA receptors clustered at neocortical synapses. J Neurosci. 2002;22(4):1328–37.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Turrigiano GG, Leslie KR, Desai NS, Rutherford LC, Nelson SB. Activity-dependent scaling of quantal amplitude in neocortical neurons. Nature. 1998;391(6670):892–6.

    Article  PubMed  CAS  Google Scholar 

  5. Davis GW, DiAntonio A, Petersen SA, Goodman CS. Postsynaptic PKA controls quantal size and reveals a retrograde signal that regulates presynaptic transmitter release in Drosophila. Neuron. 1998;20(2):305–15.

    Article  PubMed  CAS  Google Scholar 

  6. LeMasson G, Marder E, Abbott L. Activity-dependent regulation of conductances in model neurons. Science. 1993;259(5103):1915–7.

    Article  PubMed  CAS  Google Scholar 

  7. Barnes SJ, Franzoni E, Jacobsen RI, Erdelyi F, Szabo G, Clopath C, Keller GB, Keck T. Deprivation-induced homeostatic spine scaling in vivo is localized to dendritic branches that have undergone recent spine loss. Neuron. 2017;96(4):871-882. e875.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Hengen KB, Lambo ME, Van Hooser SD, Katz DB, Turrigiano GG. Firing rate homeostasis in visual cortex of freely behaving rodents. Neuron. 2013;80(2):335–42.

    Article  PubMed  CAS  Google Scholar 

  9. Ransdell JL, Nair SS, Schulz DJ. Rapid homeostatic plasticity of intrinsic excitability in a central pattern generator network stabilizes functional neural network output. J Neurosci. 2012;32(28):9649–58.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Gainey MA, Aman JW, Feldman DE. Rapid disinhibition by adjustment of PV intrinsic excitability during whisker map plasticity in mouse S1. J Neurosci. 2018;38(20):4749–61.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Kuhlman SJ, Olivas ND, Tring E, Ikrar T, Xu X, Trachtenberg JT. A disinhibitory microcircuit initiates critical-period plasticity in the visual cortex. Nature. 2013;501(7468):543–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Wu C-H, Ramos R, Katz DB, Turrigiano GG. Homeostatic synaptic scaling establishes the specificity of an associative memory. Current biology. 2021;31(11):2274-2285. e2275.

    Article  PubMed  CAS  Google Scholar 

  13. Marder E, Rue MC. From the neuroscience of individual variability to climate change. J Neurosci. 2021;41(50):10213–21.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Tattersall GJ, Ultsch GR. Physiological ecology of aquatic overwintering in ranid frogs. Biol Rev. 2008;83(2):119–40.

    Article  PubMed  Google Scholar 

  15. Santin JM, Hartzler LK. Activation of respiratory muscles does not occur during cold-submergence in bullfrogs, Lithobates catesbeianus. J Exp Biol. 2017;220(7):1181–6.

    PubMed  Google Scholar 

  16. Santin JM, Vallejo M, Hartzler LK. Synaptic up-scaling preserves motor circuit output after chronic, natural inactivity. eLife. 2017;6:e30005.

  17. Zubov T. do Amaral-Silva L, Santin JM: Inactivity and Ca2+ signaling regulate synaptic compensation in motoneurons following hibernation in American bullfrogs. Sci Rep. 2022;12(1):1–13.

    Article  Google Scholar 

  18. Amaral-Silva L, Santin JM. Synaptic modifications transform neural networks to function without oxygen. BMC Biol. 2023;21;54. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12915-023-01518-0.

  19. Bueschke N. do Amaral-Silva L, Adams S, Santin JM: Transforming a neural circuit to function without oxygen and glucose delivery. Curr Biol. 2021;31(24):R1564–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Rue MC, Alonso LM, Marder E: Repeated applications of high potassium elicit long-term changes in a motor circuit from the crab, Cancer borealis. Iscience 2022, 25(9).

  21. Alonso LM, Rue MC, Marder E. Gating of homeostatic regulation of intrinsic excitability produces cryptic long-term storage of prior perturbations. Proc Natl Acad Sci. 2023;120(26): e2222016120.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Arendt KL, Sarti F, Chen L. Chronic inactivation of a neural circuit enhances LTP by inducing silent synapse formation. J Neurosci. 2013;33(5):2087–96.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Mitra A, Mitra SS, Tsien RW. Heterogeneous reallocation of presynaptic efficacy in recurrent excitatory circuits adapting to inactivity. Nat Neurosci. 2012;15(2):250–7.

    Article  CAS  Google Scholar 

  24. Broch L, Morales RD, Sandoval AV, Hedrick MS. Regulation of the respiratory central pattern generator by chloride-dependent inhibition during development in the bullfrog (Rana catesbeiana). J Exp Biol. 2002;205(8):1161–9.

    Article  PubMed  Google Scholar 

  25. Wright AH. North American anura: life-histories of the anura of Ithaca. New York: Carnegie institution of Washington; 1914.

    Google Scholar 

  26. Willis YL, Moyle DL, Baskett TS: Emergence, breeding, hibernation, movements and transformation of the bullfrog, Rana catesbeiana, in Missouri. Copeia. 1956;1956(1):30–41.

  27. Kottick A, Baghdadwala MI, Ferguson EV, Wilson RJ. Transmission of the respiratory rhythm to trigeminal and hypoglossal motor neurons in the American Bullfrog (Lithobates catesbeiana). Respir Physiol Neurobiol. 2013;188(2):180–91.

    Article  PubMed  Google Scholar 

  28. Kinkead R. Phylogenetic trends in respiratory rhythmogenesis: insights from ectothermic vertebrates. Respir Physiol Neurobiol. 2009;168(1):39–48.

    Article  PubMed  Google Scholar 

  29. Milsom W, Reid S, Meier J, Kinkead R. Central respiratory pattern generation in the bullfrog, Rana catesbeiana. Comp Biochem Physiol A: Mol Integr Physiol. 1999;124(3):253–64.

    Article  PubMed  CAS  Google Scholar 

  30. Reid SG, Milsom WK. Respiratory pattern formation in the isolated bullfrog (Rana catesbeiana) brainstem-spinal cord. Respir Physiol. 1998;114(3):239–55.

    Article  PubMed  CAS  Google Scholar 

  31. Saunders SE, Santin JM: Activation of respiratory-related bursting in an isolated medullary section from adult bullfrogs. J Exp Biol. 2023;226(18).

  32. Martinez D, Santin JM, Schulz D, Nadim F. The differential contribution of pacemaker neurons to synaptic transmission in the pyloric network of the Jonah crab, Cancer borealis. J Neurophysiol. 2019;122(4):1623–33.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Swensen AM, Golowasch J, Christie AE, Coleman MJ, Nusbaum MP, Marder E. GABA and responses to GABA in the stomatogastric ganglion of the crab Cancer borealis. J Exp Biol. 2000;203(14):2075–92.

    Article  PubMed  CAS  Google Scholar 

  34. Amaral-Silva L, Santin JM. A brainstem preparation allowing simultaneous access to respiratory motor output and cellular properties of motoneurons in American bullfrog. J Exp Biol. 2022.

  35. Santin JM, Hartzler LK: Control of lung ventilation following overwintering conditions in bullfrogs, Lithobates catesbeianus. Journal of Experimental Biology 2016, In press, https://doiorg.publicaciones.saludcastillayleon.es/10.1242/jeb.136259.

  36. Gottlieb G, Jackson DC. Importance of pulmonary ventilation in respiratory control in the bullfrog. Am J Physiol. 1976;230(3):608–13.

    Article  PubMed  CAS  Google Scholar 

  37. Morales RD, Hedrick MS. Temperature and pH/CO2 modulate respiratory activity in the isolated brainstem of the bullfrog Rana catesbeiana. Comp Biochem Physiol- A Mol Int Physiol. 2002;132(2):477–87.

    Article  Google Scholar 

  38. Fournier S, Allard M, Roussin S, Kinkead R. Developmental changes in central O2 chemoreflex in Rana catesbeiana: the role of noradrenergic modulation. J Exp Biol. 2007;210(17):3015–26.

    Article  PubMed  Google Scholar 

  39. Vallejo M, Santin JM, Hartzler LK. Noradrenergic modulation determines respiratory network activity during temperature changes in the in vitro brainstem of bullfrogs. Respir Physiol Neurobiol. 2018;258:25–31.

    Article  PubMed  CAS  Google Scholar 

  40. Kam K, Worrell JW, Janczewski WA, Cui Y, Feldman JL. Distinct inspiratory rhythm and pattern generating mechanisms in the preBötzinger complex. J Neurosci. 2013;33(22):9235–45.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Del Negro CA, Funk GD, Feldman JL: Breathing matters. Nat Rev Neurosci. 2018.

  42. Santin JM, Hartzler LK. Environmentally‐induced return to juvenile‐like chemosensitivity in the respiratory control system of adult bullfrogs, Lithobates catesbeianus. J Physiol. 2016.

  43. Burggren WW. Phenotypic switching resulting from developmental plasticity: fixed or reversible? Front Physiol. 2020;10:1634.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Santin J, Hartzler L. Temperature influences CO2/pH-sensitivity in locus coeruleus neurons of the bullfrog, Lithobates catesbeianus. FASEB J. 2013;27(714.12).

  45. Amaral-Silva L, Santin JM. Molecular profiling of CO2/pH-sensitive neurons in the locus coeruleus of bullfrogs reveals overlapping noradrenergic and glutamatergic cell identity. Comp Biochem Physiol A: Mol Integr Physiol. 2023;283: 111453.

    Article  PubMed  CAS  Google Scholar 

  46. Turrigiano G, LeMasson G, Marder E. Selective regulation of current densities underlies spontaneous changes in the activity of cultured neurons. J Neurosci. 1995;15(5):3640–52.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Bueschke N, Amaral-Silva L, Hu M, Alvarez A, Santin JM. Plasticity in the functional properties of NMDA receptors improves network stability during severe energy stress. J Neuroscie. 2024;44.

  48. Marder E: Individual variability, statistics, and the resilience of nervous systems of crabs and humans to temperature and other perturbations. Eneuro. 2023;10(11).

  49. Chesler M. Regulation and modulation of pH in the brain. Physiol Rev. 2003;83(4):1183–221.

    Article  PubMed  CAS  Google Scholar 

  50. Haley JA, Hampton D, Marder E. Two central pattern generators from the crab, Cancer borealis, respond robustly and differentially to extreme extracellular pH. Life. 2018;7:e41877.

    Google Scholar 

  51. Brockhaus J, Ballanyi K. Anticonvulsant A1 receptor-mediated adenosine action on neuronal networks in the brainstem–spinal cord of newborn rats. Neuroscience. 2000;96(2):359–71.

    Article  PubMed  CAS  Google Scholar 

  52. Matesz C, Szekely G. Organization of the ambiguus nucleus in the frog (Rana esculenta). Journal of Comparative Neurology. 1996;371(2):258–69.

    Article  PubMed  CAS  Google Scholar 

  53. McLean H, Kimura N, Kogo N, Perry SF, Remmers JE. Fictive respiratory rhythm in the isolated brainstem of frogs. J Comp Physiol A. 1995;176(5):703–13.

    Article  PubMed  CAS  Google Scholar 

  54. Kogo N, Remmers JE. Neural organization of the ventilatory activity in the frog. Rana catesbeiana II J neurobiol. 1994;25(9):1080–94.

    Article  PubMed  CAS  Google Scholar 

  55. Stuesse SL, Cruce WL, Powell KS. Organization within the cranial IX–X complex in ranid frogs: a horseradish peroxidase transport study. Journal of Comparative Neurology. 1984;222(3):358–65.

    Article  PubMed  CAS  Google Scholar 

  56. Winmill RE, Chen AK, Hedrick MS. Development of the respiratory response to hypoxia in the isolated brainstem of the bullfrog Rana catesbeiana. J Exp Biol. 2005;208(2):213–22.

    Article  PubMed  Google Scholar 

  57. Kogo N, Perry SF, Remmers JE. Neural organization of the ventilatory activity in the frog. Rana catesbeiana I J neurobiol. 1994;25(9):1067–79.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr. Joe Viteri for comments on an earlier version of the manuscript, as well as several anonymous reviewers that improved the quality of the manuscript.

Funding

National Institutes of Health R01NS114514 to JS.

Author information

Authors and Affiliations

Authors

Contributions

Conceived Research; JS, SS. Designed research; JS, SS. Performed experiments; SS. Analyzed data; SS, JS. Drafted original manuscript; SS, JS. Edited, revised, and finalized manuscript; SS, JS.

Corresponding author

Correspondence to Joseph M. Santin.

Ethics declarations

Ethics approval and consent to participate

Animal experiments were approved by University of Missouri Animal Care and Use Committee protocol #39264.

Consent for publication

The authors each consent to the publication of this work.

Competing interests

Santin is a Guest Editor of the Special Issue “Synaptic plasticity, learning, and memory” at BMC Biology. The authors have no other competing interests to declare.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Saunders, S.E., Santin, J.M. Hibernation reduces GABA signaling in the brainstem to enhance motor activity of breathing at cool temperatures. BMC Biol 22, 251 (2024). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12915-024-02050-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12915-024-02050-5

Keywords