Significance of extracellular potassium in central respiratory control studied in the isolated brainstem-spinal cord preparation of the neonatal rat

A novel mechanism for ramping bursts based on slow negative feedback in model respiratory neurons

Recordings from pre-Bötzinger complex neurons responsible for the inspiratory phase of the respiratory rhythm reveal a ramping burst pattern, starting around the time that the transition from expiration to inspiration begins, in which the spike rate gradually rises until a transition into a high-frequency burst occurs. The spike rate increase along the burst is accompanied by a gradual depolarization of the plateau potential that underlies the spikes. These effects may be functionally important for inducing the onset of inspiration and hence maintaining effective respiration; however, most mathematical models for inspiratory bursting do not capture this activity pattern. Here, we study how the modulation of spike height and afterhyperpolarization via the slow inactivation of an inward current can support various activity patterns including ramping bursts. We use dynamical systems methods designed for multiple timescale systems, such as bifurcation analysis based on timescale decomposition and averaging over fast oscillations, to generate an understanding of and predictions about the specific dynamic effects that lead to ramping bursts. We also analyze how transitions between ramping and other activity patterns may occur with parameter changes, which could be associated with experimental manipulations, environmental conditions, and/or development.

Interdependence of cellular and network properties in respiratory rhythmogenesis

How breathing is generated by the preBötzinger Complex (preBötC) remains divided between two ideological frameworks, and the persistent sodium current (INaP) lies at the heart of this debate. Although INaP is widely expressed, the pacemaker hypothesis considers it essential because it endows a small subset of neurons with intrinsic bursting or "pacemaker" activity. In contrast, burstlet theory considers INaP dispensable because rhythm emerges from "pre-inspiratory" spiking activity driven by feed-forward network interactions. Using computational modeling, we discover that changes in spike shape can dissociate INaP from intrinsic bursting. Consistent with many experimental benchmarks, conditional effects on spike shape during simulated changes in oxygenation, development, extracellular potassium, and temperature alter the prevalence of intrinsic bursting and pre-inspiratory spiking without altering the role of INaP. Our results support a unifying hypothesis where INaP and excitatory network interactions, but not intrinsic bursting or pre-inspiratory spiking, are critical interdependent features of preBötC rhythmogenesis.

Putting the theory into 'burstlet theory' with a biophysical model of burstlets and bursts in the respiratory preBötzinger complex

Inspiratory breathing rhythms arise from synchronized neuronal activity in a bilaterally distributed brainstem structure known as the preBötzinger complex (preBötC). In in vitro slice preparations containing the preBötC, extracellular potassium must be elevated above physiological levels (to 7–9 mM) to observe regular rhythmic respiratory motor output in the hypoglossal nerve to which the preBötC projects. Reexamination of how extracellular K⁺ affects preBötC neuronal activity has revealed that low-amplitude oscillations persist at physiological levels. These oscillatory events are subthreshold from the standpoint of transmission to motor output and are dubbed burstlets. Burstlets arise from synchronized neural activity in a rhythmogenic neuronal subpopulation within the preBötC that in some instances may fail to recruit the larger network events, or bursts, required to generate motor output. The fraction of subthreshold preBötC oscillatory events (burstlet fraction) decreases sigmoidally with increasing extracellular potassium. These observations underlie the burstlet theory of respiratory rhythm generation. Experimental and computational studies have suggested that recruitment of the non-rhythmogenic component of the preBötC population requires intracellular Ca²⁺ dynamics and activation of a calcium-activated nonselective cationic current. In this computational study, we show how intracellular calcium dynamics driven by synaptically triggered Ca²⁺ influx as well as Ca²⁺ release/uptake by the endoplasmic reticulum in conjunction with a calcium-activated nonselective cationic current can reproduce and offer an explanation for many of the key properties associated with the burstlet theory of respiratory rhythm generation. Altogether, our modeling work provides a mechanistic basis that can unify a wide range of experimental findings on rhythm generation and motor output recruitment in the preBötC.

Dynamics of ramping bursts in a respiratory neuron model

Intensive computational and theoretical work has led to the development of multiple mathematical models for bursting in respiratory neurons in the pre-Bötzinger Complex (pre-BötC) of the mammalian brainstem. Nonetheless, these previous models have not captured the pre-inspiratory ramping aspects of these neurons' activity patterns, in which relatively slow tonic spiking gradually progresses to faster spiking and a full-blown burst, with a corresponding gradual development of an underlying plateau potential. In this work, we show that the incorporation of the dynamics of the extracellular potassium ion concentration into an existing model for pre-BötC neuron bursting, along with some parameter adjustments, suffices to induce this ramping behavior. Using fast-slow decomposition, we show that this activity can be considered as a form of parabolic bursting, but with burst termination at a homoclinic bifurcation rather than as a SNIC bifurcation. We also investigate the parameter-dependence of these solutions and show that the proposed model yields a greater dynamic range of burst frequencies, durations, and duty cycles than those produced by other models in the literature.

Calcium Imaging Analysis of Cellular Responses to Hypercapnia and Hypoxia in the NTS of Newborn Rat Brainstem Preparation

It is supposed that the nucleus of the solitary tract (NTS) in the dorsal medulla includes gas sensor cells responsive to hypercapnia or hypoxia in the central nervous system. In the present study, we analyzed cellular responses to hypercapnia and hypoxia in the NTS region of newborn rat in vitro preparation. The brainstem and spinal cord were isolated from newborn rat (P0-P4) and were transversely cut at the level of the rostral area postrema. To detect cellular responses, calcium indicator Oregon Green was pressure-injected into the NTS just beneath the cut surface of either the caudal or rostral block of the medulla, and the preparation was superfused with artificial cerebrospinal fluid (25–26°C). We examined cellular responses initially to hypercapnic stimulation (to 8% CO₂ from 2% CO₂) and then to hypoxic stimulation (to 0% O₂ from 95% O₂ at 5% CO₂). We tested these responses in standard solution and in two different synapse blockade solutions: (1) cocktail blockers solution including bicuculline, strychnine, NBQX and MK-801 or (2) TTX solution. At the end of the experiments, the superfusate potassium concentration was lowered to 0.2 from 3 mM to classify recorded cells into neurons and astrocytes. Excitation of cells was detected as changes of fluorescence intensity with a confocal calcium imaging system. In the synaptic blockade solutions (cocktail or TTX solution), 7.6 and 8% of the NTS cells responded to hypercapnic and hypoxic stimulation, respectively, and approximately 2% of them responded to both stimulations. Some of these cells responded to low K⁺, and they were classified into astrocytes comprising 43% hypercapnia-sensitive cells, 56% hypoxia-sensitive cells and 54% of both stimulation-sensitive cells. Of note, 49% of the putative astrocytes identified by low K⁺ stimulation were sensitive to hypercapnia, hypoxia or both. In the presence of a glia preferential blocker, 5 mM fluoroacetate (plus 0.5 μM TTX), the percentage of hypoxia-sensitive cells was significantly reduced compared to those of all other conditions. This is the first study to reveal that the NTS includes hypercapnia and hypoxia dual-sensitive cells. These results suggest that astrocytes in the NTS region could act as a central gas sensor.

Structural and functional identification of two distinct inspiratory neuronal populations at the level of the phrenic nucleus in the rat cervical spinal cord

The diaphragm is driven by phrenic motoneurons that are located in the cervical spinal cord. Although the anatomical location of the phrenic nucleus and the function of phrenic motoneurons at a single cellular level have been extensively analyzed, the spatiotemporal dynamics of phrenic motoneuron group activity have not been fully elucidated. In the present study, we analyzed the functional and structural characteristics of respiratory neuron population in the cervical spinal cord at the level of the phrenic nucleus by voltage imaging, together with histological analysis of neuronal and astrocytic distribution in the cervical spinal cord. We found spatially distinct two cellular populations that exhibited synchronized inspiratory activity on the transversely cut plane at C4–C5 levels and on the ventral surface of the mid cervical spinal cord in the isolated brainstem–spinal cord preparation of the neonatal rat. Inspiratory activity of one group emerged in the central portion of the ventral horn that corresponded to the central motor column, and the other appeared in the medial portion of the ventral horn that corresponded to the medial motor column. We identified by retrogradely labeling study that the anatomical distributions of phrenic and scalene motoneurons coincided with optically detected central and medial motor regions, respectively. Furthermore, we anatomically demonstrated closely located features of putative motoneurons, interneurons and astrocytes in these regions. Collectively, we report that phrenic and scalene motoneuron populations show synchronized inspiratory activities with distinct anatomical locations in the mid cervical spinal cord.

Preinspiratory calcium rise in putative pre-Botzinger complex astrocytes

The neural inspiratory activity originates from a ventrolateral medullary region called the pre-Bötzinger complex (preBötC), yet the mechanism underlying respiratory rhythmogenesis is not completely understood. Recently, the role of not only neurons but astrocytes in the central respiratory control has attracted considerable attention. Here we report our discovery that an intracellular calcium rise in a subset of putative astrocytes precedes inspiratory neuronal firing in rhythmically active slices. Functional calcium imaging from hundreds of preBötC cells revealed that a subset of putative astrocytes exhibited rhythmic calcium elevations preceding inspiratory neuronal activity with a time lag of approximately 2 s. These preinspiratory putative astrocytes maintained their rhythmic activities even during the blockade of neuronal activity with tetrodotoxin, whereas the rhythm frequency was lowered and the intercellular phases of these rhythms were decoupled. In addition, optogenetic stimulation of preBötC putative astrocytes induced firing of inspiratory neurons. These findings raise the possibility that astrocytes in the preBötC are actively involved in respiratory rhythm generation in rhythmically active slices.

Developmental nicotine exposure alters neurotransmission and excitability in hypoglossal motoneurons

Hypoglossal motoneurons (XII MNs) control muscles of the mammalian tongue and are rhythmically active during breathing. Acetylcholine (ACh) modulates XII MN activity by promoting the release of glutamate from neurons that express nicotinic ACh receptors (nAChRs). Chronic nicotine exposure alters nAChRs on neurons throughout the brain, including brain stem respiratory neurons. Here we test the hypothesis that developmental nicotine exposure (DNE) reduces excitatory synaptic input to XII MNs. Voltage-clamp experiments in rhythmically active medullary slices showed that the frequency of excitatory postsynaptic currents (EPSCs) onto XII MNs from DNE animals is reduced by 61% (DNE = 1.7 ± 0.4 events/s; control = 4.4 ± 0.6 events/s; P < 0.002). We also examine the intrinsic excitability of XII MNs to test whether cells from DNE animals have altered membrane properties. Current-clamp experiments showed XII MNs from DNE animals had higher intrinsic excitability, as evaluated by measuring their response to injected current. DNE cells had high-input resistances (DNE = 131.9 ± 13.7 MΩ, control = 78.6 ± 9.7 MΩ, P < 0.008), began firing at lower current levels (DNE = 144 ± 22 pA, control = 351 ± 45 pA, P < 0.003), and exhibited higher frequency-current gain values (DNE = 0.087 ± 0.012 Hz/pA, control = 0.050 ± 0.004 Hz/pA, P < 0.02). Taken together, our data show previously unreported effects of DNE on XII MN function and may also help to explain the association between DNE and the incidence of central and obstructive apneas.

Respiratory and metabolic acidosis differentially affect the respiratory neuronal network in the ventral medulla of neonatal rats

Two respiratory-related areas, the para-facial respiratory group/retrotrapezoid nucleus (pFRG/RTN) and the pre-Bötzinger complex/ventral respiratory group (preBötC/VRG), are thought to play key roles in respiratory rhythm. Because respiratory output patterns in response to respiratory and metabolic acidosis differ, we hypothesized that the responses of the medullary respiratory neuronal network to respiratory and metabolic acidosis are different. To test these hypotheses, we analysed respiratory-related activity in the pFRG/RTN and preBötC/VRG of the neonatal rat brainstem-spinal cord in vitro by optical imaging using a voltage-sensitive dye, and compared the effects of respiratory and metabolic acidosis on these two populations. We found that the spatiotemporal responses of respiratory-related regional activities to respiratory and metabolic acidosis are fundamentally different, although both acidosis similarly augmented respiratory output by increasing respiratory frequency. PreBötC/VRG activity, which is mainly inspiratory, was augmented by respiratory acidosis. Respiratory-modulated pixels increased in the preBötC/VRG area in response to respiratory acidosis. Metabolic acidosis shifted the respiratory phase in the pFRG/RTN; the pre-inspiratory dominant pattern shifted to inspiratory dominant. The responses of the pFRG/RTN activity to respiratory and metabolic acidosis are complex, and involve either augmentation or reduction in the size of respiratory-related areas. Furthermore, the activation pattern in the pFRG/RTN switched bi-directionally between pre-inspiratory/inspiratory and post-inspiratory. Electrophysiological study supported the results of our optical imaging study. We conclude that respiratory and metabolic acidosis differentially affect activities of the pFRG/RTN and preBötC/VRG, inducing switching and shifts of the respiratory phase. We suggest that they differently influence the coupling states between the pFRG/RTN and preBötC/VRG.

Two modes of respiratory rhythm generation in the newborn rat brainstem-spinal cord preparation

Two respiration-related rhythm generators, the pre-Bötzinger complex inspiratory and the parafacial pre-inspiratory rhythm generators, have been identified in the medulla of rodents that produce intrinsic periodic bursts. Although both generators can be independently active under specific conditions, they interact as a coupled oscillator system to regulate the respiratory rhythm. Here, we summarize different mechanisms of modulation of the respiratory rhythm in the brainstem-spinal cord preparation of newborn rats and discuss factors determining rhythm generator dominance. We show two different modes of respiratory rhythm generation in the brainstem-spinal cord preparation that depends on the background stimulation level.

Postnatal developmental changes in activation profiles of the respiratory neuronal network in the rat ventral medulla

Two putative respiratory rhythm generators (RRGs), the para-facial respiratory group (pFRG) and the pre-Bötzinger complex (preBötC), have been identified in the neonatal rodent brainstem. To elucidate their functional roles during the neonatal period, we evaluated developmental changes of these RRGs by optical imaging using a voltage-sensitive dye. Optical signals, recorded from the ventral medulla of brainstem-spinal cord preparations of neonatal (P0-P4) rats (n = 44), were analysed by a cross correlation method. With development during the first few postnatal days, the respiratory-related activity in the pFRG reduced and shifted from a preinspiratory (P0-P1) to an inspiratory (P2-P4) pattern, whereas preBötC activity remained unchanged. The mu-opioid agonist [D-Ala(2),N-Me-Phe(4),Gly(5)-ol]-enkephalin (DAMGO) augmented preinspiratory activity in the pFRG, while the mu-opioid antagonist naloxone induced changes in spatiotemporal activation profiles that closely mimicked the developmental changes. These results are consistent with the recently proposed hypothesis by Janczewski and Feldman that the pFRG is activated to compensate for the depression of the preBötC by perinatal opiate surge. We conclude that significant reorganization of the respiratory neuronal network, characterized by a reduction of preinspiratory activity in the pFRG, occurs at P1-P2 in rats. The changes in spatiotemporal activation profiles of the pFRG neurones may reflect changes in the mode of coupling of the two respiratory rhythm generators.

Dependence on extracellular Ca2+/K+ antagonism of inspiratory centre rhythms in slices and en bloc preparations of newborn rat brainstem

The pre-Bötzinger Complex (preBötC) inspiratory centre remains active in isolated brainstem-spinal cords and brainstem slices. The extent to which findings in these models depend on their dimensions or superfusate [K(+)] and [Ca(2+)] (both of which determine neuronal excitability) is not clear. We report here that inspiratory-related rhythms in newborn rat slices and brainstem-spinal cords with defined boundaries were basically similar in physiological Ca(2+) (1.2 mm) and K(+) (3 mm). Hypoglossal nerve rhythm was 1 : 1-coupled to preBötC activity in slices and to cervical nerve bursts in en bloc preparations lacking the facial motonucleus (VII). Hypoglossal rhythm was depressed in brainstems containing (portions of) VII, while pre/postinspiratory lumbar nerve bursting was present only in preparations with > 79% VII. preBötC-related slice rhythms were inhibited in 1.5 mm Ca(2+) solution, whereas their longevity and burst rate were substantially augmented in 1 mm Ca(2+). Ca(2+) depression of slice rhythms was antagonized by raising superfusate K(+) to 8-10 mm. This strong extracellular Ca(2+)/K(+) antagonism of inspiratory (motor) rhythms was also revealed in brainstem-spinal cords without VII, while the inhibition was progressively attenuated with increasing amount of rostral tissue. We hypothesize that depression of hypoglossal rhythm and decreased Ca(2+) sensitivity of preBötC rhythm are probably not related to an increased content of rostral respiratory structures, but rather to larger brainstem dimensions resulting in interstitial gradients for neuromodulator(s) and K(+), respectively. We discuss whether block of pre/postinspiratory activity in preparations with < 79% VII is due to impairment of the pathway from preinspiratory interneurons to abdominal muscles.

Early postnatal changes in respiratory activity in ratin vitroand modulatory effects of substance P

Developmental changes in the respiratory activity and its modulation by substance P (SP) were studied in the neonatal rat brainstem-spinal cord preparation from the day of birth to day 3 (P0-P3). The respiratory network activity in the ventrolateral medulla was represented by two types of bursts: basic regular bursts with typical decrementing shape and biphasic bursts appearing after augmented biphasic discharges in inspiratory neurons. With advancing postnatal age the respiratory output was considerably modified; the basic rhythm became faster by 20%, whereas the biphasic burst rate, which was originally 15 times slower, declined further by 180% and the C4 burst duration significantly decreased by 20% due to reduced decay time without preceding changes in the central inspiratory drive. SP had an age-dependent excitatory effect on respiratory activity. In the basic rhythm, SP could induce transient rhythm cessations on P0-P2 but not on P3. For the biphasic burst frequency, the sensitivity to SP significantly decreased from P0 to P3, whereas the range of SP-induced changes increased. In both types of bursts, SP prolonged C4 burst duration due to increasing decay time. This effect was three times greater on P3 and did not depend on the central inspiratory drive. Our results suggest that the potency of SP to regulate the respiratory activity elevates during the early postnatal period. The developmental changes in the respiratory activity appear to represent the transient stage in the maturation of rhythm and pattern generation mechanisms facilitating adaptive behavior of a quickly growing organism.

Ischemia-induced disturbance of neuronal network function in the rat spinal cord analyzed by voltage-imaging

Using a voltage-imaging technique, we analyzed the acute effect of ischemia, hypoxia and hypoglycemia on the neuronal network function of the rat spinal cord. Ischemic, hypoxic, or hypoglycemic stress was loaded to spinal cord slices with an oxygen- and glucose-free, oxygen-free, or glucose-free mock cerebrospinal fluid, respectively. Depolarizing signals in the dorsal horn, induced by dorsal root stimulation, consisted of fast (pre-synaptic) and slow (post-synaptic) components. The slow component was attenuated much more than the fast component under an ischemic condition (P<0.0002). Post-synaptic neuronal activities in lamina III-IV were suppressed earlier than those in lamina I-II. The nerve fiber was relatively resistant to ischemia. As long as the fast component was preserved in the dorsal horn, the suppression of the fast and slow components was reversible. There was a significant difference (P<0.05) in the recovered slow component sizes between the group in which the fast component was suppressed by more than 20% by ischemia and the group in which the suppression was less than 20%. Further prolonged stress irreversibly eliminated most of the slow component, and attenuated the fast component (to 59+/-8%) accompanied by cellular damage in histology. Suppression of neural activity by hypoxic or hypoglycemic stress was less prominent than that by ischemia. Prolonged ischemic stress suddenly and irreversibly eliminated depolarizing signals in the ventral horn accompanied by morphological damage of motoneurons. Immunohistochemical staining was negative for apoptosis. We have, for the first time, analyzed the processes of spinal cord disturbance induced by ischemia, hypoxia and hypoglycemia at the neuronal network level by directly observing the regional neuronal network activities within the spinal cord. We conclude that synaptic transmission in the dorsal horn, especially in deep regions, is vulnerable and first affected by these stresses. Severe ischemic stress induces irreversible dysfunction of neurons accompanied by eventual cell death in both dorsal and ventral horns.

Spontaneous crossed phrenic activity in the neonatal respiratory network

Hemisection of the cervical spinal cord causes paralysis of the ipsilateral hemidiaphragm in adult rats. Activation of a latent crossed phrenic motor pathway can restore diaphragmatic function, although structural changes take place before the pathway can be activated. Since mechanisms are employed to eliminate non-functional projections during development, we predicted that this latent neural pathway might be active during development. Therefore, we examined the effect of spinal hemisection (C2) on respiratory-like activity bilaterally using the brainstem--spinal cord preparation from neonatal rats (0-4 days). Spontaneous crossed phrenic activity (respiratory-like activity recorded from the ipsilateral C4 or C5 ventral roots following C2 hemisection) was observed in an age-dependent manner; younger preparations exhibited more than older preparations. Increasing drive (increasing [K+] or superfusion of theophylline) either increased or induced crossed phrenic activity. Hemisection caused no change in the frequency, the burst area, duration or peak amplitude contralateral to hemisection. Unlike adult rats, this study shows that crossed phrenic activity is present in the in vitro respiratory network of neonatal rats suggesting that a crossed neural pathway may be functionally active in neonates.