Abstract
How rhythms are generated by neuronal networks is fundamental to understand rhythmic behaviors such as respiration, locomotion, and mastication. Respiratory rhythm is generated by the preBötzinger complex (preBötC), an anatomically and functionally discrete population of brainstem neurons, central and necessary for respiratory rhythm. In specific in vitro conditions, preBötC neurons depend on voltage-dependent inward currents to generate respiratory rhythm. In the mature and intact organism, where preBötC neurons are deeply embedded in the respiratory network, the contribution of ionic currents to respiratory rhythm is unclear. We propose that a set of ionic currents plays a key role in generating respiratory rhythm in the mature organism in vivo. By microperfusing ionic current blockers into the preBötC of adult rats, we identify the hyperpolarization-activated cation current as a critical component of the mechanism promoting respiratory rhythm, and that this current, in combination with the persistent sodium current, is essential to respiratory rhythm in vivo. Importantly, both currents contribute to rhythmic activity in states of anesthesia, quiet wakefulness, and sleep, but not when the organism is engaged in active behaviors. These data show that a set of ionic currents at the preBötC imparts the network with rhythmicity in reduced states of arousal, although the network can override their contribution to adjust its activity for nonrhythmic behaviors in active wakefulness.
Introduction
Neuronal networks mediate essential rhythmic behaviors, such as locomotion, mastication, and breathing (Feldman and Del Negro, 2006; Zhong et al., 2007). Breathing emerges from the respiratory rhythmic network which expresses a regular and simple rhythmic motor activity, yet orchestrated by a complex network of brainstem neurons. At its core is the preBötzinger complex (preBötC), an anatomically and functionally discrete population of brainstem neurons, central and necessary for respiratory rhythm (Smith et al., 1991; Gray et al., 2001). In the preBötC, voltage-dependent ionic currents, such as the persistent sodium current INaP or the calcium-activated cation current, play key roles in mediating respiratory rhythm (Pena et al., 2004), and their roles are apparent in specific in vitro (Pena et al., 2004; Paton et al., 2006) and in vivo conditions (Pena and Aguileta, 2007; St-John et al., 2007). INaP often interacts with other inward currents to promote rhythm (Pape, 1996). The hyperpolarization-activated cation current (Ih), for instance, mediates rhythms in various networks when combined with INaP (Biel et al., 2009). Ih was detected in preBötC neurons (Mironov et al., 2000) but blockade of Ih current in immature rodents studied in vitro elicits conflicting results as it had either no effect on rhythm (Mironov et al., 2000) or increased respiratory rhythm (Thoby-Brisson et al., 2000). Ih and INaP expressions, however, substantially increase as the organism matures (Bayliss et al., 1994; Cho et al., 2011), therefore limiting our understanding of their contributions in the mature network. Here, we aim to understand the role of Ih and INaP in the generation of respiratory rhythm in the intact and mature organism in vivo.
Respiratory rhythm is a robust autonomic motor behavior that is continuously generated to maintain homeostasis. In reduced states of brain arousal, such as sleep or anesthesia, respiratory rhythm is generated by mechanisms within the brainstem respiratory network, that do not require voluntary activation of respiratory muscles. The network has, however, the capacity to override automatic rhythm to coordinate its activity with behaviors, such as grooming or feeding. In fact, respiratory rhythm is typically erratic when an organism is behaviorally active and closely follows behaviors, but it becomes more regular when the organism progresses from quiet wakefulness to deep sleep (Phillipson and Bowes, 1986). PreBötC neurons are likely involved in these state-dependent changes because their destruction or inhibition induces unstable rhythm and respiratory depression especially in states of reduced brain arousal (McKay et al., 2005; Montandon et al., 2011). Here, we propose that the robust rhythmicity observed in reduced arousal states is dependent on ionic currents, but that this mechanism is overridden when the organism is engaged in active behaviors. By applying agents to the preBötC with anatomical and functional specificity (Montandon et al., 2011), we show that Ih or INaP blockers reduced respiratory rhythm and abolished breathing only when applied in combination. Importantly, we demonstrate that these blockers also decreased respiratory rhythm in the freely behaving organism, with this effect only present in quiet wakefulness and sleep, but not in active wakefulness. In summary, Ih and INaP contribute to respiratory rhythm in states of reduced brain arousal in the intact rodent, although the respiratory network can override their contributions to coordinate its activity with nonrespiratory behaviors in active wakefulness.
Materials and Methods
All procedures were performed in accordance with the recommendations of the Canadian Council on Animal Care, and were approved by the University of Toronto Animal Care Committee.
Anesthetized preparations.
To determine the contribution of Ih and INaP at the preBötC on respiratory activities, we used reverse-microdialysis to microperfuse selected agents into the preBötC of anesthetized adult male Wistar rats (n = 34) weighing between 280 and 320 g. The experimental procedures were described previously (Montandon et al., 2011). Briefly, we recorded diaphragm and genioglossus muscle activities in isoflurane-anesthetized (2–2.5%), tracheotomized and spontaneously breathing (50% oxygen–gas mixture, balance nitrogen) adult rats while microperfusing agents into the preBötC. Electrical signals were amplified, filtered, and averaged with a moving time window (100 ms). Raw and averaged signals were recorded on a computer with Spike 2 software (v 6 and Micro-1401, Cambridge Electronic Design). Using a dorsal approach, a microdialysis probe (CX-I-12–01, diameter 200 μm, length of diffusing membrane 1 mm, Eicom) was inserted into the brainstem 2.0 mm dorsal to the preBötC using a stereotaxic frame and micromanipulator with a resolution of 50 μm (ASI Instruments). The microdialysis probe was continuously perfused at 3 μl/min with freshly made artificial CSF (aCSF). The composition of the aCSF was (in mm): 125 NaCl, 3 KCl, 1 KH2PO4, 2 CaCl2, 1 MgSO4, 25 NaHCO3, and 30 glucose. The pH was adjusted at 7.4 by bubbling CO2 in the aCSF. The probe was initially placed 12.2 mm posterior, 2 mm lateral, and 8.5 mm ventral to bregma, was then progressively lowered into the brainstem while recording all physiological variables, and was left in place when genioglossus muscle activity decreased by ∼30% which typically occurred with the probe ∼10.5 mm ventral to bregma. We note that this response is a good initial marker of the appropriate depth of the probe before subsequent confirmation by postmortem histology (Montandon et al., 2011).
Baseline levels of the physiological variables were recorded for at least 30 min. Following this control period, the Ih blockers ZD7288 (100 μm) and zatebradine (50 μm; Tocris Bioscience), the INaP blocker riluzole (50 μm, riluzole hydrochloride; Tocris Bioscience), the voltage-gated sodium channel (NaV) activator veratridine (100 nm; Tocris Bioscience), and the NaV1 blocker ranolazine (100 μm, Tocris Bioscience) were added to the aCSF for microperfusion into the preBötC region. The responses to ZD7288 (n = 9), zatebradine (n = 3), riluzole (n = 10), veratridine (n = 5), and ranolazine (n = 4) were recorded for the next 60 min. For the riluzole experiments, riluzole microperfusion was then followed by washout with aCSF for 60 min. In a separate set of experiments, we also bilaterally microperfused either ZD7288 (100 μm) or a combination of 50 μm riluzole and 100 μm ZD7288 into both preBötCs to determine whether it would abolish respiratory activity altogether (n = 3). To evaluate the frequency-dependent effect of ZD7288 and riluzole, we correlated the percentage change in respiratory frequency with the initial (i.e., baseline) conditions. The initial baseline frequency varied between 30 and 50 breath/min due to the inherent variation between preparations due to the variable sensitivity to anesthetized agents and probe placement. To evaluate instantaneous variability of the cyclic period (T), we plotted Tn+1 as a function of the preceding cycle Tn and calculated the index of cycle-by-cycle variability. This index is equivalent to the SD of the successive differences of cyclic periods. By comparing each cycle with its preceding cycle, this index excludes the variability due to slow degradation of rhythm over time.
Changes from aCSF in response to drug treatment were determined with one-way repeated-measure ANOVAs (drug application being the repeated factor) followed by Dunnet's post hoc tests for comparison with a single control (i.e., aCSF). p < 0.05 was considered statistically significant. Data are presented as means ± SEM.
Freely behaving preparations.
To determine the effects of Ih and INaP blockade at the preBötC on respiratory activity across sleep–wake states, we perfused riluzole (50 μm), ZD7288 (100 μm), or a combination of ZD7288 (100 μm) and riluzole (50 μm) bilaterally into both preBötCs of freely behaving male adult Wistar rats (n = 13). One week before the experiments, sterile surgery was performed under isoflurane anesthesia to implant the rats with electroencephalogram (EEG) and postural (neck) muscle electrodes to identify sleep–wake states, and diaphragm and genioglossus muscle electrodes for respiratory muscle recordings as previously described (Montandon et al., 2011). Two microdialysis guide cannulas (AG-8, Eicom) were positioned 5.0 mm above the preBötC by placing them 12.2 mm posterior, 2 mm either side, and 5.5 mm ventral to bregma with the guide cannulas secured in place with dental acrylic. The rat recovered for 1 week before the experiment.
On the day of the experiment, the rat was transiently anesthetized (isoflurane 2–2.5%) for the careful and accurate placement of the microdialysis probes (CX-I-12–01, Eicom) as described above for the anesthetized preparation, i.e., to ensure consistent placements within and between preparations. After placements of the microdialysis probes, the rat was connected to the recording apparatus which allowed electrophysiological signals to be recorded while the rat moved freely in a large open-topped Plexiglas bowl filled with fresh bedding, food, and water. The bowl was placed on a rotating turntable (Raturn, BASi) which automatically adjusts its position when the rat moves to avoid entanglements of the microdialysis tubing and recording cable. After a recovery period from anesthesia of at least 120 min (with the rat apparently behaving normally after ∼45 min), the EEG, neck and diaphragm muscle activities were recorded during bilateral perfusion of aCSF into both preBötC and subsequent drugs. For every 10 s epoch, sleep–wake state was classified as active wakefulness, quiet wakefulness, rapid eye movement (REM), and non-REM sleep, and physiological values were calculated for every epoch. Prevailing sleep–wake states were identified according to standard criteria (Morrison et al., 2003). Briefly, EEG frequencies in the following frequency bands: δ2 (0.5–2 Hz), δ1 (2–4 Hz), θ (4–7.5 Hz), α (7.5–13.5 Hz), β1 (13.5–20 Hz), and β2 (20–30 Hz) were calculated. Active wakefulness was characterized by low δ frequencies and high neck muscle activity, whereas quiet wakefulness low δ frequencies and low neck muscle activity. Non-REM sleep was characterized by high δ frequencies, high EEG amplitude, and low neck muscle activity, whereas REM sleep presented low δ frequencies, high θ frequencies, and low neck muscle activity. Averaged baseline values were calculated over a 30 min period before blockers were microperfused. Riluzole and/or ZD7288 were then added to the perfusing solution for 240 min and average values were calculated over the last 30 min of this period. In the event that drugs did not change rhythm, DAMGO, a μ-opioid receptor agonist known to depress frequency by its action on preBötC neurons (Montandon et al., 2011), was microperfused as positive control for accurate placement of probe. Measurements started at ∼11.00 and ended at ∼17.00. Data were amplified, filtered, moving-time averaged, sampled, and analyzed as described previously (Montandon et al., 2011).
To determine synchronization of postural neck muscle and diaphragm muscle activities, we computed the cross-correlation function using MATLAB Signal Processing Toolbox (function xcorr MATLAB R12, Mathworks). This function is a measure of similarity of two waveforms. Cross-correlation function was calculated for each 10 s segment and associated with a specific state.
Two-way repeated-measure ANOVAs (sleep–wake states and drug application being the repeated factors) followed by Holm–Sidak post hoc tests were used to determine the state-dependent effect of manipulation of preBötC for each physiological variable. p < 0.05 was considered significant.
Construction of correlation maps.
To determine the locations of the intervention (perfusion) sites, we cut 50 μm sections with cryostat and stained them with neutral red. Under light microscope (BX41, Olympus), we followed the track made by the microperfusion probe and identified the end of the track as the tip of the probe. To determine the anterior–posterior (AP) coordinate, we used the caudal end of the facial nucleus as the reference point (located 11.6 mm posterior to bregma) and counted the number of sections caudal to the facial nucleus until the perfusion site was reached. The dorsal-ventral (DV) and medial-lateral (ML) coordinates were then identified using the nucleus ambiguus, inferior olive, and standard brain maps (Paxinos and Watson, 1998). We constructed correlation maps to relate the location of the intervention sites with the resultant effect on respiratory activity as previously described (Montandon et al., 2011). The rationale for the construction of correlation maps is that for a locus of effect of drugs at any particular brainstem site, the latency for the drug to diffuse through the tissue and to progressively change respiratory activity will vary as a function of the distance of the probe from the effective site. We then calculated the distance from an arbitrary set of coordinates such as the preBötC to the perfusion site. We estimated the latency to a 10% change in respiratory frequency in response to ZD7288 or riluzole. Using this approach, one latency and one distance was obtained for each animal/experiment. For a specific set of experiments including many rats, we then calculated the relationship between latency of drug effects and distances from microperfusion sites to the preBötC. Each experiment in each rat contributed a single point to the correlation. This approach can also be applied to every possible set of coordinates within a 3-D grid (resolution 50 μm) spanning from 11.6 to 12.6 mm caudal to bregma. Correlation coefficients were then calculated between distances and latencies to see whether it is likely that the drugs are acting on other neural sites. Overall, this approach identifies the sites in the brainstem that respond the fastest to agent perfusion. Using MATLAB 12 software (Mathworks), correlation coefficients (0<r2<1) were calculated for every set of coordinates and plotted as color pixels (blue to red) in standard brain maps. To determine the relationship between the location of microperfusion site and the latency to induce an effect on respiratory rhythm, this approach requires to use a unilateral approach where only one neural site is progressively affected by drug microperfusion. For this reason, correlation maps were only constructed for the experiments involving unilateral microperfusion.
Immunohistochemistry of neurokinin-1 receptors.
We used neurokinin-1 receptor (NK1R) immunohistochemistry to locate preBötC neurons in brainstem sections (Gray et al., 2001) using techniques previously described (Montandon et al., 2011). Following fixation in 4% paraformaldehyde and cryoprotection in 30% sucrose, the brains (n = 3) were frozen and cut in 50 μm sections. The antibodies used were rabbit anti-NK1R (1:1000, AB-N04; Advanced Targeting Systems) and donkey anti-rabbit immunoglobin G (1:100, Jackson ImmunoResearch Laboratories). Sections were stained with a chromogen diaminobenzidine solution and were mounted on slides, dried, and sealed with Cytoseal 280. Sections spanning from 11.6 to 12.6 mm posterior to bregma were digitized with a CCD camera (Infinity 1, Olympus) and microscope (BX-41, Olympus). We then counted the number of NK1R-positive cells in a 1-mm-diameter circle placed ventral to the nucleus ambiguus on each section of the ventrolateral medulla (every 200 μm, 11.6–12.6 caudal to bregma). We then calculated the density of NK1R-positive cells as number of positive cells per mm2 and identified the preBötC as the region with the highest density.
Results
Ih contributes to respiratory rhythm in vivo
To determine the contribution of Ih to the generation of respiratory rhythm in vivo, we locally microperfused the Ih blocker ZD7288 into the preBötC of anesthetized adult rats while recording respiratory activity. Unilateral microperfusion of ZD7288 (100 μm) into the preBötC (Fig. 1A–C), at concentration known to substantially block Ih in vitro (Gasparini and DiFrancesco, 1997; Thoby-Brisson et al., 2000), decreased respiratory frequency (fR) by 37.7% (aCSF, fR = 44.1 ± 3.3 breath/min; ZD7288, fR = 28.0 ± 2.6 breath/min; p = 0.029, n = 5) (Fig. 1C), but did not significantly change the amplitude of diaphragm muscle activity (p = 0.273, n = 5). We used Poincare plots to quantify cycle-by-cycle variability and showed that ZD7288 increased short-term cycle variability by 180.2% compared with aCSF (p = 0.007) (Fig. 1D). Microperfusion into the preBötC of zatebradine (50 μm), another Ih blocker (Matt et al., 2011), decreased respiratory frequency by 25.8% (aCSF, fR = 42.2 breath/min; zatebradine, fR = 31.3 ± 0.7 breath/min; p = 0.002, n = 3) (Fig. 1E), replicating the effect produced by ZD7288. There also was a significant decreased in diaphragm muscle amplitude (decrease of 13.4 ± 2.1%, p = 0.023, n = 3).
To overcome the spatial limitations inherent to the localized pharmacological manipulation of selected neuronal groups in specific brain regions, we used a functional and anatomical approach to identify the sites in the medulla most sensitive to the applied agents (Montandon et al., 2011). Using the capacity of ZD7288 to diffuse through the tissue and progressively induce a reduction in respiratory frequency depending on physical proximity of the intervention site to the effective site, we related the latency for ZD7288 to decrease respiratory frequency by 10% to the proximity of the microperfusion site to the preBötC (Fig. 2A). Microperfusion of ZD7288 within the preBötC caused a relatively fast decrease in respiratory frequency, whereas perfusion further from the preBötC induced a slower decrease in respiratory frequency (Fig. 2B). There was a significant relationship between latency of response and distance from the preBötC to the perfusion sites (r = 0.836, p = 0.005, n = 9). For all the possible coordinates within the volume of brainstem surrounding the preBötC (see Materials and Methods), we calculated coefficients of correlation relating the latencies for respiratory slowing with distances from these coordinates to the perfusion sites (Fig. 2C). Construction of an anatomical map reveals “hotspot” regions that are statistically highly correlated with the latencies for respiratory slowing. The anatomical region most strongly associated with respiratory frequency slowing (illustrated in red on correlation map) corresponds to the preBötC as identified by high NK1R expression (Fig. 2C–D), suggesting that the sites overlapping with the preBötC are most rapidly sensitive to ZD7288. Importantly, this approach also showed that it is Ih blockade at the preBötC that elicits a significant decrease in respiratory frequency and that effects of ZD7288 were not due to diffusion of the drug to nearby respiratory nuclei.
INaP contributes to respiratory rhythm in vivo
To determine the contribution of INaP to the generation of respiratory rhythm in vivo, we locally microperfused the INaP blocker riluzole into the preBötC of anesthetized adult rats. While recording rhythmic output from the diaphragm, riluzole was unilaterally microperfused at concentrations known to substantially block INaP in vitro (Del Negro et al., 2002; Pena et al., 2004). Riluzole (50 μm) into the preBötC (Fig. 3A) decreased respiratory frequency by 25.1% (p = 0.025, n = 5) (Fig. 3B–D), increased respiratory rate variability by 160.5% (p = 0.035, n = 5) (Fig. 3E), but did not affect diaphragm muscle amplitude (p = 0.086, n = 5) (Fig. 3F). Riluzole increased cycle-by-cycle short-term variability (p = 0.033, n = 5) (Fig. 3G). The riluzole effect on respiratory rhythm was not transient (Fig. 3C) as previously reported (Del Negro et al., 2005) but persisted as long as riluzole was perfused, and this effect was reversed after 45 ± 5 min of aCSF washout (Fig. 3B,D).
Using a similar approach as with ZD7288, we related the latency for riluzole to decrease respiratory frequency by 10% to the proximity of the perfusion sites to the preBötC (Fig. 4A). Again there was a significant relationship between the latencies of responses and the distances of the microperfusion sites to the preBötC (r = 0.76, p = 0.010, n = 10) (Fig. 4B). A correlation map (Fig. 4C) showed that the anatomical region most strongly associated with slowing of respiratory frequency (illustrated in red in the correlation maps) corresponds to the preBötC as identified by NK1R expression. This result shows that the sites overlapping with preBötC are most rapidly sensitive to riluzole and that riluzole effect was not due to diffusion of the drugs to nearby respiratory nuclei, such as the medullary raphe as previously suggested (Del Negro et al., 2010).
To determine whether the observed decrease in respiratory frequency is specific to INaP blockade, we used ranolazine, another INaP blocker known to preferentially block persistent over transient sodium currents of NaV1 channels (Kahlig et al., 2010). These channels are highly expressed in preBötC neurons (Ptak et al., 2005) and may contribute to INaP in these cells. Microperfusion of ranolazine (100 μm) into the preBötC of anesthetized adult rats decreased respiratory frequency by 25.6% (p = 0.001, n = 4) (Fig. 5A–C), without affecting diaphragm muscle amplitude (p = 0.103). Subsequent addition of riluzole (50 μm) to the perfusion medium containing ranolazine did not further decrease respiratory frequency (Fig. 5C) suggesting that riluzole and ranolazine are blocking similar channels.
If reduction of INaP in preBötC neurons decreases respiratory frequency, then INaP activation should increase it. Accordingly, we examined the effect on the INaP activator veratridine (Fekete et al., 2009) applied to the preBötC on respiratory rhythm. Veratridine targets the neurotoxin receptor site 2 of NaV channels and evokes a persistent sodium influx (Fekete et al., 2009). Microperfusion into the preBötC of veratridine at 100 nm, a concentration known to elicit INaP in vitro, increased respiratory frequency by 18.0% (aCSF, fR = 36.7 ± 3.2 breath/min; veratridine, fR = 43.3 ± 3.4 breath/min; p = 0.039, n = 5) (Fig. 5D–E), without affecting diaphragm amplitude (p = 0.207, n = 5).
Frequency-dependency of Ih but not INaP
One of the properties of Ih is that it can regulate oscillatory activity over a dynamic range from 0.5 to 4 Hz (Luthi and McCormick, 1998; Biel et al., 2009). This range corresponds to respiratory frequencies >30 breath/min. On the other hand, INaP can support frequencies over a broad dynamic range including lower frequencies (Koizumi and Smith, 2008). To test the frequency-dependent contribution of Ih and INaP to respiratory rhythm, we correlated the capacity of ZD7288 (100 μm) or riluzole (50 μm) to decrease respiratory frequency with the initial respiratory frequency observed in baseline conditions, i.e., in the absence of blockers. At similar levels of anesthetic, different rats exhibited various initial baseline frequencies likely due to variable sensitivities to anesthetics and/or variable placements of microperfusion probes into the preBötC region. This inherent variability provided us with a wide range of initial baseline frequencies to test the frequency-dependent contribution of Ih or INaP blockers. A significant correlation of the response to ZD7288 and the initial baseline frequency was observed (r = −0.88, p = 0.009) (Fig. 6A,B). It showed that Ih blockade has a substantial effect on reducing rhythm when it is initially high, but less so when it is low and this effect fits with the operating range of Ih. Conversely, there was no significant correlation between frequencies and INaP blockade (r = 0.26, p = 0.615) (Fig. 6C,D), showing that riluzole reduced respiratory rhythm regardless of frequency.
Ih and INaP contribute to respiratory rhythm
To determine whether Ih and INaP are essential to generate respiratory rhythm in vivo, we tested whether concomitant blockade of these currents would abolish respiratory activity. To first test the essential role of Ih, we perfused ZD7288 (100 μm) bilaterally into both preBötCs and observed the consequent respiratory decline in anesthetized rats. After 45 min of microperfusion, Ih blockade decreased significantly respiratory frequency by 44.0% (aCSF, fR = 48.7 ± 1.5 breath/min; bilateral ZD7288, fR = 24.0 ± 7.3 breath/min; p = 0.049, n = 3) (Fig. 7A) without significantly affecting the amplitude of diaphragm activity (p = 0.184), and did not abolished breathing within the 2 h recording period. In a separate set of experiments, we microperfused the Ih and INaP blockers into both preBötC. ZD7288 (100 μm) and riluzole (50 μm) at the preBötCs (Fig. 7B) substantially decreased respiratory frequency by 86.7% within 45 min (aCSF, fR = 43.2 ± 3.9 breath/min; bilateral ZD7288+riluzole, fR = 6.7 ± 6.7 breath/min, p = 0.023, n = 3) and reduced diaphragm amplitude by 75.5% (p = 0.041). Breathing was completely abolished after 43.3 ± 17.4 min (n = 3). This effect strongly differs from the effects observed with ZD7288 applied alone to both preBötCs.
State-dependent contribution of Ih and INaP to respiratory rhythm
Respiratory rhythm is altered by anesthetics and modulated by states of arousal (Phillipson and Bowes, 1986). To determine whether Ih and INaP contribute to respiratory rhythm in the absence of anesthetics and across sleep–wake states (Fig. 8A), we microperfused Ih and INaP blockers into preBötCs of freely behaving and nonanesthetized adult rats as previously described (Montandon et al., 2011). Representative recordings (Fig. 8B,C) and group data (Fig. 8D) showed that concomitant bilateral microperfusion of ZD7288 and riluzole disrupted rhythm by decreasing respiratory frequency (p = 0.031, n = 5) (Fig. 8B,C) in quiet wakefulness (by 22.5%, p = 0.027), non-REM sleep (by 28.4%, p = 0.010), and REM sleep (by 25.7%, p = 0.010), but not in active wakefulness (p = 0.720). The respiratory pattern was also affected by Ih and INaP blockade at the preBötCs. Inspiratory durations were overall increased by the blockers across the sleep–wake cycle (p = 0.045) and expiratory durations were increased only in states of quiet wakefulness, non-REM, and REM sleep, but not in active wakefulness (p < 0.001), therefore closely following the effect observed on frequency.
Using an approach similar to the anesthesia experiments, we related the average distance from perfusion sites to preBötCs with the efficacy of the blockers to depress respiratory frequency. Two microperfusion sites were located for each experiment and the two distances from each microperfusion site to its respective preBötC were calculated and averaged. We then correlated these average distances collected in 5 experiments and the respective efficacy of blockers to decrease respiratory frequency in each experiment. There was a significant correlation (r = 0.80, p = 0.044, n = 6) (Fig. 8F) between average distances from microperfusion sites to preBötCs and the severity of reductions by ZD7288/riluzole suggesting that the closer the probes were to the preBötC, the stronger were the frequency reductions. Incidentally, sleep was also disrupted by perfusion of ZD7288 and riluzole into preBötCs with a reduced time spent in non-REM sleep (p = 0.010) (Fig. 8G). Bilateral microperfusion of ZD7288 (100 μm) or riluzole (50 μm) alone into the preBötC did not significantly change respiratory frequency (p = 0.967, n = 5 and p = 0.351, n = 3, respectively) (Fig. 9) across sleep–wake states.
Synchronization of respiratory and behavioral motor activities
The data described above showed that Ih and INaP at the preBötC contribute to rhythmogenesis in reduced states of brain arousal, such as sleep and anesthesia, but not in active wakefulness when the organism is engaged in complex behaviors. Active behaviors such as grooming or feeding are not automatic behaviors generated by the brainstem respiratory network, but rather arise from distinct cortical and cerebellar areas that impinge on the respiratory network and influence its activity to coordinate respiratory network activity to these nonrespiratory behaviors. To determine whether respiratory network activity may be modulated by other brain areas involved in distinct behaviors, we evaluated synchronization of two separate muscle activities originating from separate brain circuits. The degree of synchronization of respiratory network activity and postural motor activity was calculated using cross-correlation functions (Katz et al., 2002) between respiratory diaphragm and postural neck muscle activities. This measure shows whether two signals are synchronized independently of their intensity, which is essential to determine synchronization in signals of low amplitude such as motor activity during sleep. We showed that diaphragm muscle activity was tightly synchronized to neck muscle activity when the animal was behaviorally active or in quiet wakefulness, but not when the animal was in non-REM sleep (p = 0.002, n = 3) (Fig. 10A) or REM sleep (p = 0.002). When Ih and INaP were blocked, however, postural neck and diaphragm muscles were synchronized at all states (p > 0.737) (Fig. 10B,C) which differ from the state-dependent synchronization of muscle activities observed in control conditions (p = 0.001).
Discussion
Here, we determined the contribution of voltage-dependent inward currents in respiratory rhythmogenesis in vivo. We showed that Ih and INaP, two ionic currents known to mediate rhythms in many networks, mediate respiratory rhythm in states of anesthesia, quiet wakefulness, and sleep, although the network can supersede their contribution to comply its activity with active behaviors that originate from distinct neuronal circuits.
The ionic players
The unique property of Ih to activate upon hyperpolarization beyond resting potential makes it a potential current for providing rhythmicity in the respiratory network (Luthi and McCormick, 1998). When the Ih blocker ZD7288 was applied to rhythmically active brainstem sections containing the preBötC, Ih was reduced, but either no reduction of respiratory rhythm was observed (Mironov et al., 2000) or respiratory rhythm increased (Thoby-Brisson et al., 2000). This shows that Ih in preBötC neurons does not promote rhythmicity when cells are immature and isolated from the network. Ih expression is, however, reduced under these conditions, implying that its role may differ compared with mature and intact preparations (Bayliss et al., 1994; Biel et al., 2009). Here, we microperfused two distinct Ih blockers into the preBötC, at concentrations known to substantially block Ih (Williams et al., 2002; Van Bogaert and Pittoors, 2003), and we reduced respiratory rhythm. A rhythmogenic role for Ih is further evidenced by the frequency range at which it operates. Ih is involved in rhythms from 0.5 to 4 Hz (Biel et al., 2009) corresponding to respiratory frequencies between 30 and 240 breath/min. ZD7288 blocked respiratory rhythm only at frequencies >30 breath/min, consistent with Ih activation kinetics. There was also a relationship between Ih blockade and the initial baseline frequency before drug application. Frequency decline by ZD7288 was more pronounced at relatively high than low frequencies, which may explain why Ih blockers failed to depress frequency and even increased frequency in neonatal in vitro preparations where frequency is slow (Mironov et al., 2000; Thoby-Brisson et al., 2000). Similarly, respiratory rhythm was sensitive to NK1R blockade or to locus ceruleus stimulation only when the initial frequency was elevated (Doi and Ramirez, 2010), suggesting that distinct mechanisms mediate respiratory rhythm in low compared with high frequencies.
Surprisingly, we did not observe a significant effect of Ih blockade alone in the freely behaving experiments. The fact that Ih is activated by hyperpolarization suggests that, in strongly inhibited conditions such as anesthesia, Ih expression is increased and its physiological role is potentiated. In freely behaving animals, excitatory inputs to the preBötC may raise membrane potential above Ih activation potentials. It is, only when INaP is blocked, that Ih activation threshold would be reached. However, Ih role may be underestimated in our study. Although ZD7288 is the most potent Ih blocker currently available, it only partially blocks Ih (Gasparini and DiFrancesco, 1997). ZD7288 also inhibits T-type calcium channels (Sanchez-Alonso et al., 2008) and glutamate receptor-mediated currents (Chen, 2004), and depresses synaptic transmission (Chevaleyre and Castillo, 2002). Because these effects are only substantial when ZD7288 is applied at high concentrations for long periods, it is unlikely that ZD7288 effect was due to its action on other targets, but it cannot be excluded. Also, the fact that the Ih blocker zatebradine, which does not block calcium channels (Satoh and Yamada, 2002), also decreased respiratory rhythm, is consistent with a role of Ih in promoting respiratory rhythm. Zatebradine, however, inhibits potassium currents in some cell types (Satoh and Yamada, 2002), but such blockade would likely increase respiratory rhythm by reducing potassium outward current. Overall, our data reveal the importance of an essential ionic player in promoting respiratory rhythm in mature mammals in vivo, a contribution that substantially differs from its role in neonatal rodents in vitro.
However, Ih is not sufficient to generate respiratory rhythm by itself, because its blockade did not completely abolish rhythm. Ih interacts with other ionic currents. INaP for instance is important for rhythm generation in specific configurations of the respiratory network in hypoxia (Pena et al., 2004; Paton et al., 2006) and in neonatal rodents in vivo (Pena and Aguileta, 2007). Here, we show that blockade, or activation, of INaP at the preBötC, respectively decreased or increased, frequency. Correlation maps strongly suggest that riluzole is acting in the preBötC region and not in adjacent regions, such as the medullary raphe (Del Negro et al., 2010). In addition, low concentration of riluzole avoids unspecific effects observed at higher concentrations (Urbani and Belluzzi, 2000). Voltage-gated sodium channels, which are present in high density in preBötC neurons (Ptak et al., 2005), contribute to INaP in many types of neurons (Kahlig et al., 2010). Blockade of the persistent component of sodium channels with ranolazine, or activation with veratridine, substantially suppressed or increased respiratory rhythm, respectively, supporting the contribution of INaP and its underlying channels to respiratory rhythm in vivo. In the freely behaving animal, riluzole did not affect respiratory rhythm which is consistent with previous data in conscious mature rats (St-John et al., 2007). In immature rodents studied in vivo, however, respiratory rhythm was increased after systemic administration of riluzole (Pena and Aguileta, 2007), again highlighting differences between immature and mature respiratory networks.
Concomitant blockade of Ih and INaP
During hyperpolarization, Ih induces slow membrane depolarization (Pape, 1996) which activates INaP. This mechanism has been described in the cerebellum (Williams et al., 2002), in subthreshold oscillations of the entorhinal cortex (Dickson et al., 2000), and in single-spike activity in subthalamic neurons (Beurrier et al., 2000). Consistent with this mechanism, we demonstrated that when Ih or INaP alone were blocked in the anesthetized preparation, respiratory rhythm was reduced but not abolished, and only the concomitant blockade of both currents led to rhythm cessation. Ih/INaP blockade abolishes breathing only under anesthesia, but not in freely behaving organisms, suggesting either that anesthetics are essential to abolish rhythm or that an alternative mechanism of rhythmogenesis is triggered when these currents are absent in nonanesthetized animals. Similarly, reduced ventilation due to Ih/INaP blockade likely increases carbon dioxide levels, which may stimulate rhythm, and minimize the effects of Ih/INaP blockade. Overall, these currents are necessary but probably not sufficient to generate respiratory rhythm, and the contribution of other currents cannot be excluded (Pena et al., 2004; Koch et al., 2013). In summary, under normal physiological conditions and when ionic currents and the underlying channels are fully mature, respiratory rhythmogenesis depends on a set of currents which includes Ih and INaP.
State dependency
The brain can override automatic respiratory rhythm to synchronize breathing with other behaviors (such as feeding or grooming) which also engage the respiratory system (Ramirez and Garcia, 2007). Here, we showed that Ih and INaP play key roles in respiratory rhythmogenesis, but that their contribution is overridden when the organism is behaviorally active. The locus ceruleus, a pontine structure with noradrenergic projections to the preBötC (Doi and Ramirez, 2010), is active when the organism is engaged in behaviors and is inactive in sleep (Aston-Jones and Bloom, 1981; Takahashi et al., 2010). Combined with the fact that noradrenaline directly inhibits Ih (Wang et al., 2007), high locus ceruleus activity may therefore reduce Ih expression in preBötC cells to supersede rhythm and coordinate diaphragm activity to nonrhythmic motor behaviors. When Ih and INaP are blocked, however, postural and respiratory muscles are synchronized across states, suggesting that a mechanism independent of Ih/INaP maintains respiratory rhythm. Under anesthesia, Ih/INaP blockade abolishes rhythm, and volatile anesthetics may prevent the respiratory network to recruit an alternative mechanism to sustain rhythm as observed in the freely behaving animal. The mechanism mediating synchronized activity of motor outputs is unknown and it cannot be excluded that indirect effects, such as hypercapnia may trigger their synchronization. The idea of multiple mechanisms generating rhythm within the preBötC has been proposed previously (Pena et al., 2004; Smith et al., 2009) and these state-dependent changes identify that there are multiple and redundant mechanisms that functionally cooperate to maintain respiratory rhythm, with one or other predominating as circumstances vary.
Mechanisms of rhythmogenesis
Voltage-dependent inward currents, such as Ih and INaP, are found in cells of various rhythmic networks and are essential to trigger spontaneous bursting (Pena et al., 2004; Biel et al., 2009). It is unknown, however, whether subthreshold voltage-activated inward currents contribute to spontaneous pacemaker-like bursting of mature preBötC neurons. Here, we did not record the bursting of preBötC neurons in vivo, but highlight Ih/INaP importance in generating respiratory rhythm in freely behaving mature organisms. Also, Ih/INaP mediate the overall excitability of the rhythmic network (Biel et al., 2009; Del Negro et al., 2010), and their blockade may decrease neuronal excitability which would decrease rhythm independently of an action on pacemaker properties. Blockade of INaP abolishes respiratory rhythm in vitro (Pena et al., 2004), but rhythm can be reestablished by exciting the respiratory network through AMPA or NK1Rs (Del Negro et al., 2010). Similarly, our study shows that Ih/INaP blockade and the consequent rhythm inhibition may be reactivated by excitatory inputs during active states. This is consistent with the idea that ionic conductances are important for respiratory rhythm in some conditions or states (quiet wakefulness and sleep), but are superseded by excitatory inputs when the organism is behaviorally active.
Footnotes
This work was supported by the Parker B. Francis Fellowship (G.M.), the Ontario Ministry of Research and Innovation Fellowship (G.M.), the Canadian Institutes for Health Research (R.L.H.), and the Tier 1 Canada Research in Sleep and Respiratory Neurobiology (R.L.H.).
The authors declare no competing financial interests.
- Correspondence should be addressed to either Gaspard Montandon or Richard L. Horner, University of Toronto, 1 King's College Circle, Toronto, ON M5S 1A8, Canada, gaspard.montandon{at}utoronto.ca or richard.horner{at}utoronto.ca