Physiological and morphological properties of Dbx1-derived respiratory neurons in the pre-Botzinger complex of neonatal mice

Human IPSC-Derived PreBötC-Like Neurons and Development of an Opiate Overdose and Recovery Model

Opioid overdose is the leading cause of drug overdose lethality, posing an urgent need for investigation. The key brain region for inspiratory rhythm regulation and opioid-induced respiratory depression (OIRD) is the preBötzinger Complex (preBötC) and current knowledge has mainly been obtained from animal systems. This study aims to establish a protocol to generate human preBötC neurons from induced pluripotent cells (iPSCs) and develop an opioid overdose and recovery model utilizing these iPSC-preBötC neurons. A de novo protocol to differentiate preBötC-like neurons from human iPSCs is established. These neurons express essential preBötC markers analyzed by immunocytochemistry and demonstrate expected electrophysiological responses to preBötC modulators analyzed by patch clamp electrophysiology. The correlation of the specific biomarkers and function analysis strongly suggests a preBötC-like phenotype. Moreover, the dose-dependent inhibition of these neurons' activity is demonstrated for four different opioids with identified IC50's comparable to the literature. Inhibition is rescued by naloxone in a concentration-dependent manner. This iPSC-preBötC mimic is crucial for investigating OIRD and combating the overdose crisis and a first step for the integration of a functional overdose model into microphysiological systems.

The hypoxic respiratory response of the pre-Bötzinger complex

Since the discovery of the pre-Bötzinger Complex (preBötC) as a crucial region for generating the main respiratory rhythm, our understanding of its cellular and molecular aspects has rapidly increased within the last few decades. It is now apparent that preBötC is a highly flexible neuronal network that reconfigures state-dependently to produce the most appropriate respiratory output in response to various metabolic challenges, such as hypoxia. However, the responses of the preBötC to hypoxic conditions can be varied based on the intensity, pattern, and duration of the hypoxic challenge. This review discusses the preBötC response to hypoxic challenges at the cellular and network level. Particularly, the involvement of preBötC in the classical biphasic response of the respiratory network to acute hypoxia is illuminated. Furthermore, the article discusses the functional and structural changes of preBötC neurons following intermittent and sustained hypoxic challenges. Accumulating evidence shows that the preBötC neural circuits undergo substantial changes following hypoxia and contribute to several types of the respiratory system's hypoxic ventilatory responses.

Inspiratory and sigh breathing rhythms depend on distinct cellular signalling mechanisms in the preBötzinger complex

Breathing behaviour involves the generation of normal breaths (eupnoea) on a timescale of seconds and sigh breaths on the order of minutes. Both rhythms emerge in tandem from a single brainstem site, but whether and how a single cell population can generate two disparate rhythms remains unclear. We posit that recurrent synaptic excitation in concert with synaptic depression and cellular refractoriness gives rise to the eupnoea rhythm, whereas an intracellular calcium oscillation that is slower by orders of magnitude gives rise to the sigh rhythm. A mathematical model capturing these dynamics simultaneously generates eupnoea and sigh rhythms with disparate frequencies, which can be separately regulated by physiological parameters. We experimentally validated key model predictions regarding intracellular calcium signalling. All vertebrate brains feature a network oscillator that drives the breathing pump for regular respiration. However, in air-breathing mammals with compliant lungs susceptible to collapse, the breathing rhythmogenic network may have refashioned ubiquitous intracellular signalling systems to produce a second slower rhythm (for sighs) that prevents atelectasis without impeding eupnoea. KEY POINTS: A simplified activity-based model of the preBötC generates inspiratory and sigh rhythms from a single neuron population. Inspiration is attributable to a canonical excitatory network oscillator mechanism. Sigh emerges from intracellular calcium signalling. The model predicts that perturbations of calcium uptake and release across the endoplasmic reticulum counterintuitively accelerate and decelerate sigh rhythmicity, respectively, which was experimentally validated. Vertebrate evolution may have adapted existing intracellular signalling mechanisms to produce slow oscillations needed to optimize pulmonary function in mammals.

Inhibitory subpopulations in preBötzinger Complex play distinct roles in modulating inspiratory rhythm and pattern

Inhibitory neurons embedded within mammalian neural circuits shape breathing, walking, chewing, and other rhythmic motor behaviors. At the core of the neural circuit controlling breathing is the preBötzinger Complex (preBötC), a nucleus in the ventrolateral medulla necessary for generation of inspiratory rhythm. In the preBötC, a recurrently connected network of glutamatergic Dbx1-derived (Dbx1 + ) neurons generates rhythmic inspiratory drive. Functionally and anatomically intercalated among Dbx1 + preBötC neurons are GABAergic (GAD1/2 + ) and glycinergic (GlyT2 + ) neurons, whose roles in breathing remain unclear. To elucidate the inhibitory microcircuits within preBötC, we first characterized the spatial distribution of molecularly-defined inhibitory preBötC subpopulations in double reporter mice expressing either the red fluorescent protein tdTomato or EGFP in GlyT2 + , GAD1 + , or GAD2 + neurons. We found that, in postnatal mice, the majority of inhibitory preBötC neurons expressed a combination of GlyT2 and GAD2 while a much smaller subpopulation also expressed GAD1. To determine the functional role of these subpopulations, we used holographic photostimulation, a patterned illumination technique with high spatiotemporal resolution, in rhythmically active medullary slices from neonatal Dbx1 tdTomato ;GlyT2 EGFP and Dbx1 tdTomato ;GAD1 EGFP double reporter mice. Stimulation of 4 or 8 preBötC GlyT2 + neurons during endogenous rhythm prolonged the interburst interval in a phase-dependent manner and increased the latency to burst initiation when bursts were evoked by stimulation of Dbx1 + neurons. In contrast, stimulation of 4 or 8 preBötC GAD1 + neurons did not affect interburst interval or latency to burst initiation. Instead, photoactivation of GAD1 + neurons during the inspiratory burst prolonged endogenous and evoked burst duration and decreased evoked burst amplitude. We conclude that the majority of preBötC inhibitory neurons express both GlyT2 and GAD2 and modulate breathing rhythm by delaying burst initiation while a smaller GAD1 + subpopulation shapes inspiratory patterning by altering burst duration and amplitude.

Role of NaV1.6-mediated persistent sodium current and bursting-pacemaker properties in breathing rhythm generation

Inspiration is the inexorable active phase of breathing. The brainstem pre-Bötzinger complex (preBötC) gives rise to inspiratory neural rhythm, but its underlying cellular and ionic bases remain unclear. The long-standing "pacemaker hypothesis" posits that the persistent Na+ current (INaP) that gives rise to bursting-pacemaker properties in preBötC interneurons is essential for rhythmogenesis. We tested the pacemaker hypothesis by conditionally knocking out and knocking down the Scn8a (Nav1.6 [voltage-gated sodium channel 1.6]) gene in core rhythmogenic preBötC neurons. Deleting Scn8a substantially decreases the INaP and abolishes bursting-pacemaker activity, which slows inspiratory rhythm in vitro and negatively impacts the postnatal development of ventilation. Diminishing Scn8a via genetic interference has no impact on breathing in adult mice. We argue that the Scn8a-mediated INaP is not obligatory but that it influences the development and rhythmic function of the preBötC. The ubiquity of the INaP in respiratory brainstem interneurons could underlie breathing-related behaviors such as neonatal phonation or rhythmogenesis in different physiological conditions.

Microcircuit Synchronization and Heavy-Tailed Synaptic Weight Distribution Augment preBötzinger Complex Bursting Dynamics

The preBötzinger Complex (preBötC) encodes inspiratory time as rhythmic bursts of activity underlying each breath. Spike synchronization throughout a sparsely connected preBötC microcircuit initiates bursts that ultimately drive the inspiratory motor patterns. Using minimal microcircuit models to explore burst initiation dynamics, we examined the variability in probability and latency to burst following exogenous stimulation of a small subset of neurons, mimicking experiments. Among various physiologically plausible graphs of 1000 excitatory neurons constructed using experimentally determined synaptic and connectivity parameters, directed Erdős-Rényi graphs with a broad (lognormal) distribution of synaptic weights best captured the experimentally observed dynamics. preBötC synchronization leading to bursts was regulated by the efferent connectivity of spiking neurons that are optimally tuned to amplify modest preinspiratory activity through input convergence. Using graph-theoretic and machine learning-based analyses, we found that input convergence of efferent connectivity at the next-nearest neighbor order was a strong predictor of incipient synchronization. Our analyses revealed a crucial role of synaptic heterogeneity in imparting exceptionally robust yet flexible preBötC attractor dynamics. Given the pervasiveness of lognormally distributed synaptic strengths throughout the nervous system, we postulate that these mechanisms represent a ubiquitous template for temporal processing and decision-making computational motifs.SIGNIFICANCE STATEMENT Mammalian breathing is robust, virtually continuous throughout life, yet is inherently labile: to adapt to rapid metabolic shifts (e.g., fleeing a predator or chasing prey); for airway reflexes; and to enable nonventilatory behaviors (e.g., vocalization, breathholding, laughing). Canonical theoretical frameworks-based on pacemakers and intrinsic bursting-cannot account for the observed robustness and flexibility of the preBötzinger Complex rhythm. Experiments reveal that network synchronization is the key to initiate inspiratory bursts in each breathing cycle. We investigated preBötC synchronization dynamics using network models constructed with experimentally determined neuronal and synaptic parameters. We discovered that a fat-tailed (non-Gaussian) synaptic weight distribution-a manifestation of synaptic heterogeneity-augments neuronal synchronization and attractor dynamics in this vital rhythmogenic network, contributing to its extraordinary reliability and responsiveness.

µ-Opioid Receptor Activation Reduces Glutamate Release in the PreBötzinger Complex in Organotypic Slice Cultures

The inspiratory rhythm generator, located in the brainstem preBötzinger complex (preBötC), is dependent on glutamatergic signaling and is affected profoundly by opioids. Here, we used organotypic slice cultures of the newborn mouse brainstem of either sex in combination with genetically encoded sensors for Ca2+, glutamate, and GABA to visualize Ca2+, glutamatergic and GABAergic signaling during spontaneous rhythm and in the presence of DAMGO. During spontaneous rhythm, the glutamate sensor SF-iGluSnFR.A184S revealed punctate synapse-like fluorescent signals along dendrites and somas in the preBötC with decay times that were prolonged by the glutamate uptake blocker (TFB-TBOA). The GABA sensor iGABASnFR showed a more diffuse fluorescent signal during spontaneous rhythm. Rhythmic Ca2+- and glutamate transients had an inverse relationship between the spontaneous burst frequency and the burst amplitude of the Ca2+ and glutamate signals. A similar inverse relationship was observed when bath applied DAMGO reduced spontaneous burst frequency and increased the burst amplitude of Ca2+, glutamate, and GABA transient signals. However, a hypoxic challenge reduced both burst frequency and Ca2+ transient amplitude. Using a cocktail that blocked glutamatergic, GABAergic, and glycinergic transmission to indirectly measure the release of glutamate/GABA in response to an electrical stimulus, we found that DAMGO reduces the release of glutamate in the preBötC but has no effect on GABA release. This suggest that the opioid mediated slowing of respiratory rhythm involves presynaptic reduction of glutamate release, which would impact the ability of the network to engage in recurrent excitation, and may result in the opioid-induced slowing of inspiratory rhythm.SIGNIFICANCE STATEMENT Opioids slow down breathing rhythm by affecting neurons in the preBötzinger complex (preBötC) and other brainstem regions. Here, we used cultured slices of the preBötC to better understand this effect by optically recording Ca2+, glutamate, and GABA transients during preBötC activity. Spontaneous rhythm showed an inverse relationship between burst frequency and burst amplitude in the Ca2+ and glutamate signals. Application of the opioid DAMGO slowed the rhythm, with a concomitant increase in Ca2+, glutamate, and GABA signals. When rhythm was blocked pharmacologically, DAMGO reduced the presynaptic release of glutamate, but not GABA. These data suggest the mechanism of action of opioids involves presynaptic reduction of glutamate release, which may play an important role in the opioid-induced slowing of inspiratory rhythm.

Single cell transcriptome sequencing of inspiratory neurons of the preBötzinger complex in neonatal mice

Neurons in the brainstem preBötzinger complex (preBötC) generate the rhythm and rudimentary motor pattern for inspiratory breathing movements. We performed whole-cell patch-clamp recordings from inspiratory neurons in the preBötC of neonatal mouse slices that retain breathing-related rhythmicity in vitro. We classified neurons based on their electrophysiological properties and genetic background, and then aspirated their cellular contents for single-cell RNA sequencing (scRNA-seq). This data set provides the raw nucleotide sequences (FASTQ files) and annotated files of nucleotide sequences mapped to the mouse genome (mm10 from Ensembl), which includes the fragment counts, gene lengths, and fragments per kilobase of transcript per million mapped reads (FPKM). These data reflect the transcriptomes of the neurons that generate the rhythm and pattern for inspiratory breathing movements.

Optogenetic stimulation of pre-Bötzinger complex reveals novel circuit interactions in swallowing-breathing coordination

The coordination of swallowing with breathing, in particular inspiration, is essential for homeostasis in most organisms. While much has been learned about the neuronal network critical for inspiration in mammals, the pre-Bötzinger complex (preBötC), little is known about how this network interacts with swallowing. Here we activate within the preBötC excitatory neurons (defined as Vglut2 and Sst neurons) and inhibitory neurons (defined as Vgat neurons) and inhibit and activate neurons defined by the transcription factor Dbx1 to gain an understanding of the coordination between the preBötC and swallow behavior. We found that stimulating inhibitory preBötC neurons did not mimic the premature shutdown of inspiratory activity caused by water swallows, suggesting that swallow-induced suppression of inspiratory activity is not directly mediated by the inhibitory neurons in the preBötC. By contrast, stimulation of preBötC Dbx1 neurons delayed laryngeal closure of the swallow sequence. Inhibition of Dbx1 neurons increased laryngeal closure duration and stimulation of Sst neurons pushed swallow occurrence to later in the respiratory cycle, suggesting that excitatory neurons from the preBötC connect to the laryngeal motoneurons and contribute to the timing of swallowing. Interestingly, the delayed swallow sequence was also caused by chronic intermittent hypoxia (CIH), a model for sleep apnea, which is 1) known to destabilize inspiratory activity and 2) associated with dysphagia. This delay was not present when inhibiting Dbx1 neurons. We propose that a stable preBötC is essential for normal swallow pattern generation and disruption may contribute to the dysphagia seen in obstructive sleep apnea.

Inspiratory rhythm generation is stabilized by Ih

Cellular and network properties must be capable of generating rhythmic activity that is both flexible and stable. This is particularly important for breathing, a rhythmic behavior that dynamically adapts to environmental, behavioral, and metabolic changes from the first to the last breath. The pre-Bötzinger complex (preBötC), located within the ventral medulla, is responsible for producing rhythmic inspiration. Its cellular properties must be tunable, flexible as well as stabilizing. Here, we explore the role of the hyperpolarization-activated, nonselective cation current (I_h) for stabilizing PreBötC activity during opioid exposure and reduced excitatory synaptic transmission. Introducing I_h into an in silico preBötC network predicts that loss of this depolarizing current should significantly slow the inspiratory rhythm. By contrast, in vitro and in vivo experiments revealed that the loss of I_h minimally affected breathing frequency, but destabilized rhythmogenesis through the generation of incompletely synchronized bursts (burstlets). Associated with the loss of I_h was an increased susceptibility of breathing to opioid-induced respiratory depression or weakened excitatory synaptic interactions, a paradoxical depolarization at the cellular level, and the suppression of tonic spiking. Tonic spiking activity is generated by nonrhythmic excitatory and inhibitory preBötC neurons, of which a large percentage express I_h. Together, our results suggest that I_h is important for maintaining tonic spiking, stabilizing inspiratory rhythmogenesis, and protecting breathing against perturbations or changes in network state.NEW & NOTEWORTHY The I_h current plays multiple roles within the preBötC. This current is important for promoting intrinsic tonic spiking activity in excitatory and inhibitory neurons and for preserving rhythmic function during conditions that dampen network excitability, such as in the context of opioid-induced respiratory depression. We therefore propose that the I_h current expands the dynamic range of rhythmogenesis, buffers the preBötC against network perturbations, and stabilizes rhythmogenesis by preventing the generation of unsynchronized bursts.

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 Ca2+ dynamics and activation of a calcium-activated nonselective cationic current. In this computational study, we show how intracellular calcium dynamics driven by synaptically triggered Ca2+ influx as well as Ca2+ 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.

Transcriptomes of electrophysiologically recorded Dbx1-derived respiratory neurons of the preBötzinger complex in neonatal mice

Breathing depends on interneurons in the preBötzinger complex (preBötC) derived from Dbx1-expressing precursors. Here we investigate whether rhythm- and pattern-generating functions reside in discrete classes of Dbx1 preBötC neurons. In a slice model of breathing with ~ 5 s cycle period, putatively rhythmogenic Type-1 Dbx1 preBötC neurons activate 100-300 ms prior to Type-2 neurons, putatively specialized for output pattern, and 300-500 ms prior to the inspiratory motor output. We sequenced Type-1 and Type-2 transcriptomes and identified differential expression of 123 genes including ionotropic receptors (Gria3, Gabra1) that may explain their preinspiratory activation profiles and Ca2+ signaling (Cracr2a, Sgk1) involved in inspiratory and sigh bursts. Surprisingly, neuropeptide receptors that influence breathing (e.g., µ-opioid and bombesin-like peptide receptors) were only sparsely expressed, which suggests that cognate peptides and opioid drugs exert their profound effects on a small fraction of the preBötC core. These data in the public domain help explain the neural origins of breathing.

Role of DSCAM in the Development of Neural Control of Movement and Locomotion

Locomotion results in an alternance of flexor and extensor muscles between left and right limbs generated by motoneurons that are controlled by the spinal interneuronal circuit. This spinal locomotor circuit is modulated by sensory afferents, which relay proprioceptive and cutaneous inputs that inform the spatial position of limbs in space and potential contacts with our environment respectively, but also by supraspinal descending commands of the brain that allow us to navigate in complex environments, avoid obstacles, chase prey, or flee predators. Although signaling pathways are important in the establishment and maintenance of motor circuits, the role of DSCAM, a cell adherence molecule associated with Down syndrome, has only recently been investigated in the context of motor control and locomotion in the rodent. DSCAM is known to be involved in lamination and delamination, synaptic targeting, axonal guidance, dendritic and cell tiling, axonal fasciculation and branching, programmed cell death, and synaptogenesis, all of which can impact the establishment of motor circuits during development, but also their maintenance through adulthood. We discuss herein how DSCAM is important for proper motor coordination, especially for breathing and locomotion.

Analysis of temperature-dependent abnormal bursting patterns of neurons in Aplysia

Temperature affects the firing pattern and electrical activity of neurons in animals, eliciting diverse responses depending on neuronal cell type. However, the mechanisms underlying such diverse responses are not well understood. In the present study, we performed in vitro recording of abdominal ganglia cells of Aplysia juliana, and analyzed their burst firing patterns. We identified atypical bursting patterns dependent on temperature that were totally different from classical bursting patterns observed in R15 neurons of A. juliana. We classified these abnormal bursting patterns into type 1 and type 2; type 1 abnormal single bursts are composed of two kinds of spikes with a long interspike interval (ISI) followed by short ISI regular firing, while type 2 abnormal single bursts are composed of complex multiplets. To investigate the mechanism underlying the temperature dependence of abnormal bursting, we employed simulations using a modified Plant model and determined that the temperature dependence of type 2 abnormal bursting is related to temperature-dependent scaling factors and activation or inactivation of potassium or sodium channels.

Role of Synaptic Inhibition in the Coupling of the Respiratory Rhythms that Underlie Eupnea and Sigh Behaviors

The preBötzinger complex (preBötC) gives rise to two types of breathing behavior under normal physiological conditions: eupnea and sighing. Here, we examine the neural mechanisms that couple their underlying rhythms. We measured breathing in awake intact adult mice and recorded inspiratory rhythms from the preBötC in neonatal mouse brainstem slice preparations. We show previously undocumented variability in the temporal relationship between sigh breaths or bursts and their preceding eupneic breaths or inspiratory bursts. Investigating the synaptic mechanisms for this variability in vitro, we further show that pharmacological blockade of chloride-mediated synaptic inhibition strengthens inspiratory-to-sigh temporal coupling. These findings contrast with previous literature, which suggested glycinergic inhibition linked sigh bursts to their preceding inspiratory bursts with minimal time intervals. Furthermore, we verify that pharmacological disinhibition did not alter the duration of the prolonged interval that follows a sigh burst before resumption of the inspiratory rhythm. These results demonstrate that synaptic inhibition does not enhance coupling between sighs and preceding inspiratory events or contribute to post-sigh apneas. Instead, we conclude that excitatory synaptic mechanisms coordinate inspiratory (eupnea) and sigh rhythms.

Emergent Elements of Inspiratory Rhythmogenesis: Network Synchronization and Synchrony Propagation

We assessed the mechanism of mammalian breathing rhythmogenesis in the preBötzinger complex (preBötC) in vitro, where experimental tests remain inconsistent with hypotheses of canonical rhythmogenic cellular or synaptic mechanisms, i.e., pacemaker neurons or inhibition. Under rhythmic conditions, in each cycle, an inspiratory burst emerges as (presumptive) preBötC rhythmogenic neurons transition from aperiodic uncorrelated population spike activity to become increasingly synchronized during preinspiration (for ∼50-500 ms), which can trigger inspiratory bursts that propagate to motoneurons. In nonrhythmic conditions, antagonizing GABAA receptors can initiate this synchronization while inducing a higher conductance state in nonrhythmogenic preBötC output neurons. Our analyses uncover salient features of preBötC network dynamics where inspiratory bursts arise when and only when the preBötC rhythmogenic subpopulation strongly synchronizes to drive output neurons. Furthermore, downstream propagation of preBötC network activity, ultimately to motoneurons, is dependent on the strength of input synchrony onto preBötC output neurons exemplifying synchronous propagation of network activity.

Muscarinic Modulation of Morphologically Identified Glycinergic Neurons in the Mouse PreBötzinger Complex

The cholinergic system plays an essential role in central respiratory control, but the underlying mechanisms remain elusive. We used whole-cell recordings in brainstem slices from juvenile mice expressing enhanced green fluorescent protein (EGFP) under the control of the glycine transporter type 2 (GlyT₂) promoter, to examine muscarinic modulation of morphologically identified glycinergic neurons in the preBötzinger complex (preBötC), an area critical for central inspiratory rhythm generation. Biocytin-filled reconstruction of glycinergic neurons revealed that the majority of them had few primary dendrites and had axons arborized within their own dendritic field. Few glycinergic neurons had axon collaterals extended towards the premotor/motor areas or ran towards the contralateral preBötC, and had more primary dendrites and more compact dendritic trees. Spontaneously active glycinergic neurons fired regular spikes, or less frequently in a "burst-like" pattern at physiological potassium concentration. Muscarine suppressed firing in the majority of regular spiking neurons via M₂ receptor activation while enhancing the remaining neurons through M₁ receptors. Interestingly, rhythmic bursting was augmented by muscarine in a small group of glycinergic neurons. In contrast to its heterogeneous modulation of glycinergic neuronal excitability, muscarine generally depressed inhibitory and excitatory synaptic inputs onto both glycinergic and non-glycinergic preBötC neurons, with a stronger effect on inhibitory input. Notably, presynaptic muscarinic attenuation of excitatory synaptic input was dependent on M₁ receptors in glycinergic neurons and on M₂ receptors in non-glycinergic neurons. Additional field potential recordings of excitatory synaptic potentials in the M₂ receptor knockout mice indicate that glycinergic and non-glycinergic neurons contribute equally to the general suppression by muscarine of excitatory activity in preBötC circuits. In conclusion, our data show that preBötC glycinergic neurons are morphologically heterogeneous, and differ in the properties of synaptic transmission and muscarinic modulation in comparison to non-glycinergic neurons. The dominant and cell-type-specific muscarinic inhibition of synaptic neurotransmission and spiking may contribute to central respiratory disturbances in high cholinergic states.

Evaluating the Burstlet Theory of Inspiratory Rhythm and Pattern Generation

The preBötzinger complex (preBötC) generates the rhythm and rudimentary motor pattern for inspiratory breathing movements. Here, we test “burstlet” theory (Kam et al., 2013a), which posits that low amplitude burstlets, subthreshold from the standpoint of inspiratory bursts, reflect the fundamental oscillator of the preBötC. In turn, a discrete suprathreshold process transforms burstlets into full amplitude inspiratory bursts that drive motor output, measurable via hypoglossal nerve (XII) discharge in vitro. We recap observations by Kam and Feldman in neonatal mouse slice preparations: field recordings from preBötC demonstrate bursts and concurrent XII motor output intermingled with lower amplitude burstlets that do not produce XII motor output. Manipulations of excitability affect the relative prevalence of bursts and burstlets and modulate their frequency. Whole-cell and photonic recordings of preBötC neurons suggest that burstlets involve inconstant subsets of rhythmogenic interneurons. We conclude that discrete rhythm- and pattern-generating mechanisms coexist in the preBötC and that burstlets reflect its fundamental rhythmogenic nature.

Insights into the dynamic control of breathing revealed through cell-type-specific responses to substance P

The rhythm generating network for breathing must continuously adjust to changing metabolic and behavioral demands. Here, we examined network-based mechanisms in the mouse preBötzinger complex using substance P, a potent excitatory modulator of breathing frequency and stability, as a tool to dissect network properties that underlie dynamic breathing. We find that substance P does not alter the balance of excitation and inhibition during breaths or the duration of the resulting refractory period. Instead, mechanisms of recurrent excitation between breaths are enhanced such that the rate that excitation percolates through the network is increased. We propose a conceptual framework in which three distinct phases of inspiration, the burst phase, refractory phase, and percolation phase, can be differentially modulated to control breathing dynamics and stability. Unraveling mechanisms that support this dynamic control may improve our understanding of nervous system disorders that destabilize breathing, many of which involve changes in brainstem neuromodulatory systems.

Presynaptic Mechanisms and KCNQ Potassium Channels Modulate Opioid Depression of Respiratory Drive

Opioid-induced respiratory depression (OIRD) is the major cause of death associated with opioid analgesics and drugs of abuse, but the underlying cellular and molecular mechanisms remain poorly understood. We investigated opioid action in vivo in unanesthetized mice and in in vitro medullary slices containing the preBötzinger Complex (preBötC), a locus critical for breathing and inspiratory rhythm generation. Although hypothesized as a primary mechanism, we found that mu-opioid receptor (MOR1)-mediated GIRK activation contributed only modestly to OIRD. Instead, mEPSC recordings from genetically identified Dbx1-derived interneurons, essential for rhythmogenesis, revealed a prevalent presynaptic mode of action for OIRD. Consistent with MOR1-mediated suppression of presynaptic release as a major component of OIRD, Cacna1a KO slices lacking P/Q-type Ca2+ channels enhanced OIRD. Furthermore, OIRD was mimicked and reversed by KCNQ potassium channel activators and blockers, respectively. In vivo whole-body plethysmography combined with systemic delivery of GIRK- and KCNQ-specific potassium channel drugs largely recapitulated these in vitro results, and revealed state-dependent modulation of OIRD. We propose that respiratory failure from OIRD results from a general reduction of synaptic efficacy, leading to a state-dependent collapse of rhythmic network activity.

Defining the Rhythmogenic Elements of Mammalian Breathing

Breathing's remarkable ability to adapt to changes in metabolic, environmental, and behavioral demands stems from a complex integration of its rhythm-generating network within the wider nervous system. Yet, this integration complicates identification of its specific rhythmogenic elements. Based on principles learned from smaller rhythmic networks of invertebrates, we define criteria that identify rhythmogenic elements of the mammalian breathing network and discuss how they interact to produce robust, dynamic breathing.

Trpm4 ion channels in pre-Bötzinger complex interneurons are essential for breathing motor pattern but not rhythm

Inspiratory breathing movements depend on pre-Bötzinger complex (preBötC) interneurons that express calcium (Ca2+)-activated nonselective cationic current (ICAN) to generate robust neural bursts. Hypothesized to be rhythmogenic, reducing ICAN is predicted to slow down or stop breathing; its contributions to motor pattern would be reflected in the magnitude of movements (output). We tested the role(s) of ICAN using reverse genetic techniques to diminish its putative ion channels Trpm4 or Trpc3 in preBötC neurons in vivo. Adult mice transduced with Trpm4-targeted short hairpin RNA (shRNA) progressively decreased the tidal volume of breaths yet surprisingly increased breathing frequency, often followed by gasping and fatal respiratory failure. Mice transduced with Trpc3-targeted shRNA survived with no changes in breathing. Patch-clamp and field recordings from the preBötC in mouse slices also showed an increase in the frequency and a decrease in the magnitude of preBötC neural bursts in the presence of Trpm4 antagonist 9-phenanthrol, whereas the Trpc3 antagonist pyrazole-3 (pyr-3) showed inconsistent effects on magnitude and no effect on frequency. These data suggest that Trpm4 mediates ICAN, whose influence on frequency contradicts a direct role in rhythm generation. We conclude that Trpm4-mediated ICAN is indispensable for motor output but not the rhythmogenic core mechanism of the breathing central pattern generator.

Neural mechanisms for sigh generation during prenatal development

The respiratory network of the preBötzinger complex (preBötC), which controls inspiratory behavior, can in normal conditions simultaneously produce two types of inspiration-related rhythmic activities: the eupneic rhythm composed of monophasic, low-amplitude, and relatively high-frequency bursts, interspersed with sigh rhythmic activity, composed of biphasic, high-amplitude, and lower frequency bursts. By combining electrophysiological recordings from transverse brainstem slices with computational modeling, new advances in the mechanisms underlying sigh production have been obtained during prenatal development. The present review summarizes recent findings that establish when sigh rhythmogenesis starts to be produced during embryonic development as well as the cellular, membrane, and synaptic properties required for its expression. Together, the results demonstrate that although generated by the same network, the eupnea and sigh rhythms have different developmental onset times and rely on distinct network properties. Because sighs (also known as augmented breaths) are important in maintaining lung function (by reopening collapsed alveoli), gaining insight into their underlying neural mechanisms at early developmental stages is likely to help in the treatment of prematurely born babies often suffering from breathing deficiencies.

Fate mapping neurons and glia derived from Dbx1-expressing progenitors in mouse preBötzinger complex

The brainstem preBötzinger complex (preBötC) generates the inspiratory breathing rhythm, and its core rhythmogenic interneurons are derived from Dbx1‐expressing progenitors. To study the neural bases of breathing, tamoxifen‐inducible Cre‐driver mice and Cre‐dependent reporters are used to identify, record, and perturb Dbx1 preBötC neurons. However, the relationship between tamoxifen administration and reporter protein expression in preBötC neurons and glia has not been quantified. To address this problem, we crossed mice that express tamoxifen‐inducible Cre recombinase under the control of the Dbx1 gene (Dbx1^CreERT2) with Cre‐dependent fluorescent reporter mice (Rosa26^tdTomato), administered tamoxifen at different times during development, and analyzed tdTomato expression in the preBötC of their offspring. We also crossed Rosa26^tdTomato reporters with mice that constitutively express Cre driven by Dbx1 (Dbx1^Cre) and analyzed tdTomato expression in the preBötC of their offspring for comparison. We show that Dbx1‐expressing progenitors give rise to preBötC neurons and glia. Peak neuronal tdTomato expression occurs when tamoxifen is administered at embryonic day 9.5 (E9.5), whereas tdTomato expression in glia shows no clear relationship with tamoxifen timing. These results can be used to bias reporter protein expression in neurons (or glia). Tamoxifen administration at E9.5 labels 91% of Dbx1‐derived neurons in the preBötC, yet only 48% of Dbx1‐derived glia. By fate mapping Dbx1‐expressing progenitors, this study illustrates the developmental assemblage of Dbx1‐derived cells in preBötC, which can be used to design intersectional Cre/lox experiments that interrogate its cellular composition, structure, and function.

The Dynamic Basis of Respiratory Rhythm Generation: One Breath at a Time

Rhythmicity is a universal timing mechanism in the brain, and the rhythmogenic mechanisms are generally dynamic. This is illustrated for the neuronal control of breathing, a behavior that occurs as a one-, two-, or three-phase rhythm. Each breath is assembled stochastically, and increasing evidence suggests that each phase can be generated independently by a dedicated excitatory microcircuit. Within each microcircuit, rhythmicity emerges through three entangled mechanisms: ( a) glutamatergic transmission, which is amplified by ( b) intrinsic bursting and opposed by ( c) concurrent inhibition. This rhythmogenic triangle is dynamically tuned by neuromodulators and other network interactions. The ability of coupled oscillators to reconfigure and recombine may allow breathing to remain robust yet plastic enough to conform to nonventilatory behaviors such as vocalization, swallowing, and coughing. Lessons learned from the respiratory network may translate to other highly dynamic and integrated rhythmic systems, if approached one breath at a time.

Breathing matters

Breathing is a well-described, vital and surprisingly complex behaviour, with behavioural and physiological outputs that are easy to directly measure. Key neural elements for generating breathing pattern are distinct, compact and form a network amenable to detailed interrogation, promising the imminent discovery of molecular, cellular, synaptic and network mechanisms that give rise to the behaviour. Coupled oscillatory microcircuits make up the rhythmic core of the breathing network. Primary among these is the preBötzinger Complex (preBötC), which is composed of excitatory rhythmogenic interneurons and excitatory and inhibitory pattern-forming interneurons that together produce the essential periodic drive for inspiration. The preBötC coordinates all phases of the breathing cycle, coordinates breathing with orofacial behaviours and strongly influences, and is influenced by, emotion and cognition. Here, we review progress towards cracking the inner workings of this vital core.

Dbx1 Pre-Bötzinger Complex Interneurons Comprise the Core Inspiratory Oscillator for Breathing in Unanesthetized Adult Mice

The brainstem pre-Bötzinger complex (preBötC) generates inspiratory breathing rhythms, but which neurons comprise its rhythmogenic core? Dbx1-derived neurons may play the preeminent role in rhythm generation, an idea well founded at perinatal stages of development but incompletely evaluated in adulthood. We expressed archaerhodopsin or channelrhodopsin in Dbx1 preBötC neurons in intact adult mice to interrogate their function. Prolonged photoinhibition slowed down or stopped breathing, whereas prolonged photostimulation sped up breathing. Brief inspiratory-phase photoinhibition evoked the next breath earlier than expected, whereas brief expiratory-phase photoinhibition delayed the subsequent breath. Conversely, brief inspiratory-phase photostimulation increased inspiratory duration and delayed the subsequent breath, whereas brief expiratory-phase photostimulation evoked the next breath earlier than expected. Because they govern the frequency and precise timing of breaths in awake adult mice with sensorimotor feedback intact, Dbx1 preBötC neurons constitute an essential core component of the inspiratory oscillator, knowledge directly relevant to human health and physiology.

The interdependence of excitation and inhibition for the control of dynamic breathing rhythms

The preBötzinger Complex (preBötC), a medullary network critical for breathing, relies on excitatory interneurons to generate the inspiratory rhythm. Yet, half of preBötC neurons are inhibitory, and the role of inhibition in rhythmogenesis remains controversial. Using optogenetics and electrophysiology in vitro and in vivo, we demonstrate that the intrinsic excitability of excitatory neurons is reduced following large depolarizing inspiratory bursts. This refractory period limits the preBötC to very slow breathing frequencies. Inhibition integrated within the network is required to prevent overexcitation of preBötC neurons, thereby regulating the refractory period and allowing rapid breathing. In vivo, sensory feedback inhibition also regulates the refractory period, and in slowly breathing mice with sensory feedback removed, activity of inhibitory, but not excitatory, neurons restores breathing to physiological frequencies. We conclude that excitation and inhibition are interdependent for the breathing rhythm, because inhibition permits physiological preBötC bursting by controlling refractory properties of excitatory neurons.

Dendritic A-Current in Rhythmically Active PreBötzinger Complex Neurons in Organotypic Cultures from Newborn Mice

The brainstem preBötzinger complex (preBötC) generates the inspiratory rhythm for breathing. The onset of neural activity that precipitates the inspiratory phase of the respiratory cycle may depend on the activity of type-1 preBötC neurons, which exhibit a transient outward K⁺ current, I_A Inspiratory rhythm generation can be studied ex vivo because the preBötC remains rhythmically active in vitro, both in acute brainstem slices and organotypic cultures. Advantageous optical conditions in organotypic slice cultures from newborn mice of either sex allowed us to investigate how I_A impacts Ca²⁺ transients occurring in the dendrites of rhythmically active type-1 preBötC neurons. The amplitude of dendritic Ca²⁺ transients evoked via voltage increases originating from the soma significantly increased after an I_A antagonist, 4-aminopyridine (4-AP), was applied to the perfusion bath or to local dendritic regions. Similarly, glutamate-evoked postsynaptic depolarizations recorded at the soma increased in amplitude when 4-AP was coapplied with glutamate at distal dendritic locations. We conclude that I_A is expressed on type-1 preBötC neuron dendrites. We propose that I_A filters synaptic input, shunting sparse excitation, while enabling temporally summated events to pass more readily as a result of I_A inactivation. Dendritic I_A in rhythmically active preBötC neurons could thus ensure that inspiratory motor activity does not occur until excitatory synaptic drive is synchronized and well coordinated among cellular constituents of the preBötC during inspiratory rhythmogenesis. The biophysical properties of dendritic I_A might thus promote robustness and regularity of breathing rhythms.SIGNIFICANCE STATEMENT Brainstem neurons in the preBötC generate the oscillatory activity that underlies breathing. PreBötC neurons express voltage-dependent currents that can influence inspiratory activity, among which is a transient potassium current (I_A) previously identified in a rhythmogenic excitatory subset of type-1 preBötC neurons. We sought to determine whether I_A is expressed in the dendrites of preBötC. We found that dendrites of type-1 preBötC neurons indeed express I_A, which may aid in shunting sparse non-summating synaptic inputs, while enabling strong summating excitatory inputs to readily pass and thus influence somatic membrane potential trajectory. The subcellular distribution of I_A in rhythmically active neurons of the preBötC may thus be critical for producing well coordinated ensemble activity during inspiratory burst formation.

Transcriptome of neonatal preBötzinger complex neurones in Dbx1 reporter mice

We sequenced the transcriptome of brainstem interneurons in the specialized respiratory rhythmogenic site dubbed preBötzinger Complex (preBötC) from newborn mice. To distinguish molecular characteristics of the core oscillator we compared preBötC neurons derived from Dbx1-expressing progenitors that are respiratory rhythmogenic to neighbouring non-Dbx1-derived neurons, which support other respiratory and non-respiratory functions. Results in three categories are particularly salient. First, Dbx1 preBötC neurons express κ-opioid receptors in addition to μ-opioid receptors that heretofore have been associated with opiate respiratory depression, which may have clinical applications. Second, Dbx1 preBötC neurons express the hypoxia-inducible transcription factor Hif1a at levels three-times higher than non-Dbx1 neurons, which links core rhythmogenic microcircuits to O₂-related chemosensation for the first time. Third, we detected a suite of transcription factors including Hoxa4 whose expression pattern may define the rostral preBötC border, Pbx3 that may influence ipsilateral connectivity, and Pax8 that may pertain to a ventrally-derived subset of Dbx1 preBötC neurons. These data establish the transcriptomic signature of the core respiratory oscillator at a perinatal stage of development.

Morphology of Dbx1 respiratory neurons in the preBötzinger complex and reticular formation of neonatal mice

The relationship between neuron morphology and function is a perennial issue in neuroscience. Information about synaptic integration, network connectivity, and the specific roles of neuronal subpopulations can be obtained through morphological analysis of key neurons within a microcircuit. Here we present morphologies of two classes of brainstem respiratory neurons. First, interneurons derived from Dbx1-expressing precursors (Dbx1 neurons) in the preBötzinger complex (preBötC) of the ventral medulla that generate the rhythm for inspiratory breathing movements. Second, Dbx1 neurons of the intermediate reticular formation that influence the motor pattern of pharyngeal and lingual movements during the inspiratory phase of the breathing cycle. We describe the image acquisition and subsequent digitization of morphologies of respiratory Dbx1 neurons from the preBötC and the intermediate reticular formation that were first recorded in vitro. These data can be analyzed comparatively to examine how morphology influences the roles of Dbx1 preBötC and Dbx1 reticular interneurons in respiration and can also be utilized to create morphologically accurate compartmental models for simulation and modeling of respiratory circuits.

Functional Interactions between Mammalian Respiratory Rhythmogenic and Premotor Circuitry

Breathing in mammals depends on rhythms that originate from the preBötzinger complex (preBötC) of the ventral medulla and a network of brainstem and spinal premotor neurons. The rhythm-generating core of the preBötC, as well as some premotor circuits, consist of interneurons derived from Dbx1-expressing precursors (Dbx1 neurons), but the structure and function of these networks remain incompletely understood. We previously developed a cell-specific detection and laser ablation system to interrogate respiratory network structure and function in a slice model of breathing that retains the preBötC, the respiratory-related hypoglossal (XII) motor nucleus and XII premotor circuits. In spontaneously rhythmic slices, cumulative ablation of Dbx1 preBötC neurons decreased XII motor output by ∼50% after ∼15 cell deletions, and then decelerated and terminated rhythmic function altogether as the tally increased to ∼85 neurons. In contrast, cumulatively deleting Dbx1 XII premotor neurons decreased motor output monotonically but did not affect frequency nor stop XII output regardless of the ablation tally. Here, we couple an existing preBötC model with a premotor population in several topological configurations to investigate which one may replicate the laser ablation experiments best. If the XII premotor population is a "small-world" network (rich in local connections with sparse long-range connections among constituent premotor neurons) and connected with the preBötC such that the total number of incoming synapses remains fixed, then the in silico system successfully replicates the in vitro laser ablation experiments. This study proposes a feasible configuration for circuits consisting of Dbx1-derived interneurons that generate inspiratory rhythm and motor pattern.

SIGNIFICANCE STATEMENT

To produce a breathing-related motor pattern, a brainstem core oscillator circuit projects to a population of premotor interneurons, but the assemblage of this network remains incompletely understood. Here we applied network modeling and numerical simulation to discover respiratory circuit configurations that successfully replicate photonic cell ablation experiments targeting either the core oscillator or premotor network, respectively. If premotor neurons are interconnected in a so-called "small-world" network with a fixed number of incoming synapses balanced between premotor and rhythmogenic neurons, then our simulations match their experimental benchmarks. These results provide a framework of experimentally testable predictions regarding the rudimentary structure and function of respiratory rhythm- and pattern-generating circuits in the brainstem of mammals.

Respiratory Network Stability and Modulatory Response to Substance P Require Nalcn

Respiration is a rhythmic activity as well as one that requires responsiveness to internal and external circumstances; both the rhythm and neuromodulatory responses of breathing are controlled by brainstem neurons in the preBötzinger complex (preBötC) and the retrotrapezoid nucleus (RTN), but the specific ion channels essential to these activities remain to be identified. Because deficiency of sodium leak channel, non-selective (Nalcn) causes lethal apnea in humans and mice, we investigated Nalcn function in these neuronal groups. We found that one-third of mice lacking Nalcn in excitatory preBötC neurons died soon after birth; surviving mice developed apneas in adulthood. Interestingly, in both preBötC and RTN neurons, the Nalcn current influences the resting membrane potential, contributes to maintenance of stable network activity, and mediates modulatory responses to the neuropeptide substance P. These findings reveal Nalcn's specific role in both rhythmic stability and responsiveness to neuropeptides within the respiratory network.

Breathing control center neurons that promote arousal in mice

Slow, controlled breathing has been used for centuries to promote mental calming, and it is used clinically to suppress excessive arousal such as panic attacks. However, the physiological and neural basis of the relationship between breathing and higher-order brain activity is unknown. We found a neuronal subpopulation in the mouse preBötzinger complex (preBötC), the primary breathing rhythm generator, which regulates the balance between calm and arousal behaviors. Conditional, bilateral genetic ablation of the ~175 Cdh9/Dbx1 double-positive preBötC neurons in adult mice left breathing intact but increased calm behaviors and decreased time in aroused states. These neurons project to, synapse on, and positively regulate noradrenergic neurons in the locus coeruleus, a brain center implicated in attention, arousal, and panic that projects throughout the brain.

Breathing control center neurons that promote arousal in mice

Slow, controlled breathing has been used for centuries to promote mental calming, and it is used clinically to suppress excessive arousal such as panic attacks. However, the physiological and neural basis of the relationship between breathing and higher-order brain activity is unknown. We found a neuronal subpopulation in the mouse preBötzinger complex (preBötC), the primary breathing rhythm generator, which regulates the balance between calm and arousal behaviors. Conditional, bilateral genetic ablation of the ~175 Cdh9/Dbx1 double-positive preBötC neurons in adult mice left breathing intact but increased calm behaviors and decreased time in aroused states. These neurons project to, synapse on, and positively regulate noradrenergic neurons in the locus coeruleus, a brain center implicated in attention, arousal, and panic that projects throughout the brain.

Transient Suppression of Dbx1 PreBötzinger Interneurons Disrupts Breathing in Adult Mice

Interneurons derived from Dbx1-expressing precursors located in the brainstem preBötzinger complex (preBötC) putatively form the core oscillator for inspiratory breathing movements. We tested this Dbx1 core hypothesis by expressing archaerhodopsin in Dbx1-derived interneurons and then transiently hyperpolarizing these neurons while measuring respiratory rhythm in vitro or breathing in vagus-intact adult mice. Transient illumination of the preBötC interrupted inspiratory rhythm in both slice preparations and sedated mice. In awake mice, light application reduced breathing frequency and prolonged the inspiratory duration. Support for the Dbx1 core hypothesis previously came from embryonic and perinatal mouse experiments, but these data suggest that Dbx1-derived preBötC interneurons are rhythmogenic in adult mice too. The neural origins of breathing behavior can be attributed to a localized and genetically well-defined interneuron population.

Microcircuits in respiratory rhythm generation: commonalities with other rhythm generating networks and evolutionary perspectives

Rhythmicity is critical for the generation of rhythmic behaviors and higher brain functions. This review discusses common mechanisms of rhythm generation, including the role of synaptic inhibition and excitation, with a focus on the mammalian respiratory network. This network generates three phases of breathing and is highly integrated with brain regions associated with numerous non-ventilatory behaviors. We hypothesize that during evolution multiple rhythmogenic microcircuits were recruited to accommodate the generation of each breathing phase. While these microcircuits relied primarily on excitatory mechanisms, synaptic inhibition became increasingly important to coordinate the different microcircuits and to integrate breathing into a rich behavioral repertoire that links breathing to sensory processing, arousal, and emotions as well as learning and memory.

Defining preBötzinger Complex Rhythm- and Pattern-Generating Neural Microcircuits In Vivo

Normal breathing in rodents requires activity of glutamatergic Dbx1-derived (Dbx1(+)) preBötzinger Complex (preBötC) neurons expressing somatostatin (SST). We combined in vivo optogenetic and pharmacological perturbations to elucidate the functional roles of these neurons in breathing. In transgenic adult mice expressing channelrhodopsin (ChR2) in Dbx1(+) neurons, photoresponsive preBötC neurons had preinspiratory or inspiratory firing patterns associated with excitatory effects on burst timing and pattern. In transgenic adult mice expressing ChR2 in SST(+) neurons, photoresponsive preBötC neurons had inspiratory or postinspiratory firing patterns associated with excitatory responses on pattern or inhibitory responses that were largely eliminated by blocking synaptic inhibition within preBötC or by local viral infection limiting ChR2 expression to preBötC SST(+) neurons. We conclude that: (1) preinspiratory preBötC Dbx1(+) neurons are rhythmogenic, (2) inspiratory preBötC Dbx1(+) and SST(+) neurons primarily act to pattern respiratory motor output, and (3) SST(+)-neuron-mediated pathways and postsynaptic inhibition within preBötC modulate breathing pattern.

Voltage-Dependent Rhythmogenic Property of Respiratory Pre-Bötzinger Complex Glutamatergic, Dbx1-Derived, and Somatostatin-Expressing Neuron Populations Revealed by Graded Optogenetic Inhibition

The rhythm of breathing in mammals, originating within the brainstem pre-Bötzinger complex (pre-BötC), is presumed to be generated by glutamatergic neurons, but this has not been directly demonstrated. Additionally, developmental expression of the transcription factor Dbx1 or expression of the neuropeptide somatostatin (Sst), has been proposed as a marker for the rhythmogenic pre-BötC glutamatergic neurons, but it is unknown whether these other two phenotypically defined neuronal populations are functionally equivalent to glutamatergic neurons with regard to rhythm generation. To address these problems, we comparatively investigated, by optogenetic approaches, the roles of pre-BötC glutamatergic, Dbx1-derived, and Sst-expressing neurons in respiratory rhythm generation in neonatal transgenic mouse medullary slices in vitro and also more intact adult perfused brainstem-spinal cord preparations in situ. We established three different triple-transgenic mouse lines with Cre-driven Archaerhodopsin-3 (Arch) expression selectively in glutamatergic, Dbx1-derived, or Sst-expressing neurons for targeted photoinhibition. In each line, we identified subpopulations of rhythmically active, Arch-expressing pre-BötC inspiratory neurons by whole-cell recordings in medullary slice preparations in vitro, and established that Arch-mediated hyperpolarization of these inspiratory neurons was laser power dependent with equal efficacy. By site- and population-specific graded photoinhibition, we then demonstrated that inspiratory frequency was reduced by each population with the same neuronal voltage-dependent frequency control mechanism in each state of the respiratory network examined. We infer that enough of the rhythmogenic pre-BötC glutamatergic neurons also have the Dbx1 and Sst expression phenotypes, and thus all three phenotypes share the same voltage-dependent frequency control property.

Mixed-mode oscillations and population bursting in the pre-Bötzinger complex

This study focuses on computational and theoretical investigations of neuronal activity arising in the pre-Bötzinger complex (pre-BötC), a medullary region generating the inspiratory phase of breathing in mammals. A progressive increase of neuronal excitability in medullary slices containing the pre-BötC produces mixed-mode oscillations (MMOs) characterized by large amplitude population bursts alternating with a series of small amplitude bursts. Using two different computational models, we demonstrate that MMOs emerge within a heterogeneous excitatory neural network because of progressive neuronal recruitment and synchronization. The MMO pattern depends on the distributed neuronal excitability, the density and weights of network interconnections, and the cellular properties underlying endogenous bursting. Critically, the latter should provide a reduction of spiking frequency within neuronal bursts with increasing burst frequency and a dependence of the after-burst recovery period on burst amplitude. Our study highlights a novel mechanism by which heterogeneity naturally leads to complex dynamics in rhythmic neuronal populations.

DOI:http://dx.doi.org/10.7554/eLife.13403.001

Each breath we take removes carbon dioxide from the body and exchanges it for oxygen. A structure called the brainstem, which connects the brain with the spinal cord, generates the breathing rhythm and controls its rate. While this process normally occurs automatically, we can also control our breathing voluntarily, such as when singing or speaking.

Within the brainstem, a group of neurons in the area known as the pre-Bötzinger complex is responsible for ensuring that an animal breathes in at regular intervals. Recordings of the electrical activity from slices of brainstem show that pre-Bötzinger neurons display rhythmic activity with characteristic patterns called “mixed-mode oscillations”. These rhythms consist of bursts of strong activity (“large amplitude bursts”), essential for triggering regular breathing, separated by a series of bursts of weak activity (“small amplitude bursts”). However, it is not clear how mixed-mode oscillations arise.

Bacak, Kim et al. now provide insights into this process by developing two computational models of the pre-Bötzinger complex. The first model consists of a diverse population of 100 neurons joined by a relatively small number of weak connections to form a network. The second model is a simplified version of the first, consisting of just three neurons. By manipulating the properties of the simulated networks, and analysing the data mathematically, Bacak, Kim et al. identify the properties of the neurons that allow them to generate mixed-mode oscillations and thus rhythmic breathing.

The models suggest that mixed-mode oscillations result from the synchronization of many neurons with different levels of activity (excitability). Neurons with low excitability have low bursting frequencies, but generate strong activity and recruit other neurons, ultimately producing large amplitude bursts that cause breathing.

Many parts of the nervous system are also made up of networks of neurons with diverse excitability. A challenge for future studies is thus to investigate whether other networks of neurons similar to the pre-Bötzinger complex generate rhythms that control other repetitive actions, such as walking and chewing.

DOI:http://dx.doi.org/10.7554/eLife.13403.002

Dbx1 precursor cells are a source of inspiratory XII premotoneurons

All behaviors require coordinated activation of motoneurons from central command and premotor networks. The genetic identities of premotoneurons providing behaviorally relevant excitation to any pool of respiratory motoneurons remain unknown. Recently, we established in vitro that Dbx1-derived pre-Bötzinger complex neurons are critical for rhythm generation and that a subpopulation serves a premotor function (Wang et al., 2014). Here, we further show that a subpopulation of Dbx1-derived intermediate reticular (IRt) neurons are rhythmically active during inspiration and project to the hypoglossal (XII) nucleus that contains motoneurons important for maintaining airway patency. Laser ablation of Dbx1 IRt neurons, 57% of which are glutamatergic, decreased ipsilateral inspiratory motor output without affecting frequency. We conclude that a subset of Dbx1 IRt neurons is a source of premotor excitatory drive, contributing to the inspiratory behavior of XII motoneurons, as well as a key component of the airway control network whose dysfunction contributes to sleep apnea.

DOI:http://dx.doi.org/10.7554/eLife.12301.001

Organotypic slice cultures containing the preBötzinger complex generate respiratory-like rhythms

Study of acute brain stem slice preparations in vitro has advanced our understanding of the cellular and synaptic mechanisms of respiratory rhythm generation, but their inherent limitations preclude long-term manipulation and recording experiments. In the current study, we have developed an organotypic slice culture preparation containing the preBötzinger complex (preBötC), the core inspiratory rhythm generator of the ventrolateral brain stem. We measured bilateral synchronous network oscillations, using calcium-sensitive fluorescent dyes, in both ventrolateral (presumably the preBötC) and dorsomedial regions of slice cultures at 7-43 days in vitro. These calcium oscillations appear to be driven by periodic bursts of inspiratory neuronal activity, because whole cell recordings from ventrolateral neurons in culture revealed inspiratory-like drive potentials, and no oscillatory activity was detected from glial fibrillary associated protein-expressing astrocytes in cultures. Acute slices showed a burst frequency of 10.9 ± 4.2 bursts/min, which was not different from that of brain stem slice cultures (13.7 ± 10.6 bursts/min). However, slice cocultures that include two cerebellar explants placed along the dorsolateral border of the brainstem displayed up to 193% faster burst frequency (22.4 ± 8.3 bursts/min) and higher signal amplitude (340%) compared with acute slices. We conclude that preBötC-containing slice cultures retain inspiratory-like rhythmic function and therefore may facilitate lines of experimentation that involve extended incubation (e.g., genetic transfection or chronic drug exposure) while simultaneously being amenable to imaging and electrophysiology at cellular, synaptic, and network levels.

Facing the challenge of mammalian neural microcircuits: taking a few breaths may help

Breathing in mammals is a seemingly straightforward behaviour controlled by the brain. A brainstem nucleus called the preBötzinger Complex sits at the core of the neural circuit generating respiratory rhythm. Despite the discovery of this microcircuit almost 25 years ago, the mechanisms controlling breathing remain elusive. Given the apparent simplicity and well-defined nature of regulatory breathing behaviour, the identification of much of the circuitry, and the ability to study breathing in vitro as well as in vivo, many neuroscientists and physiologists are surprised that respiratory rhythm generation is still not well understood. Our view is that conventional rhythmogenic mechanisms involving pacemakers, inhibition or bursting are problematic and that simplifying assumptions commonly made for many vertebrate neural circuits ignore consequential detail. We propose that novel emergent mechanisms govern the generation of respiratory rhythm. That a mammalian function as basic as rhythm generation arises from complex and dynamic molecular, synaptic and neuronal interactions within a diverse neural microcircuit highlights the challenges in understanding neural control of mammalian behaviours, many (considerably) more elaborate than breathing. We suggest that the neural circuit controlling breathing is inimitably tractable and may inspire general strategies for elucidating other neural microcircuits.

Mechanisms Leading to Rhythm Cessation in the Respiratory PreBötzinger Complex Due to Piecewise Cumulative Neuronal Deletions

The mammalian breathing rhythm putatively originates from Dbx1-derived interneurons in the preBötzinger complex (preBötC) of the ventral medulla. Cumulative deletion of ∼15% of Dbx1 preBötC neurons in an in vitro breathing model stops rhythmic bursts of respiratory-related motor output. Here we assemble in silico models of preBötC networks using random graphs for structure, and ordinary differential equations for dynamics, to examine the mechanisms responsible for the loss of spontaneous respiratory rhythm and motor output measured experimentally in vitro. Model networks subjected to cellular ablations similarly discontinue functionality. However, our analyses indicate that model preBötC networks remain topologically intact even after rhythm cessation, suggesting that dynamics coupled with structural properties of the underlying network are responsible for rhythm cessation. Simulations show that cumulative cellular ablations diminish the number of neurons that can be recruited to spike per unit time. When the recruitment rate drops below 1 neuron/ms the network stops spontaneous rhythmic activity. Neurons that play pre-eminent roles in rhythmogenesis include those that commence spiking during the quiescent phase between respiratory bursts and those with a high number of incoming synapses, which both play key roles in recruitment, i.e., recurrent excitation leading to network bursts. Selectively ablating neurons with many incoming synapses impairs recurrent excitation and stops spontaneous rhythmic activity and motor output with lower ablation tallies compared with random deletions. This study provides a theoretical framework for the operating mechanism of mammalian central pattern generator networks and their susceptibility to loss-of-function in the case of disease or neurodegeneration.

Sigh and Eupnea Rhythmogenesis Involve Distinct Interconnected Subpopulations: A Combined Computational and Experimental Study

How a single neural network can generate several rhythmic activities at different time scales remains an open question. Here, in addition to the already described reconfiguring process, we propose a new mechanism by which the respiratory network can generate simultaneously two distinct inspiration-related activities (eupnea and sigh) at different frequencies.

Neural networks control complex motor outputs by generating several rhythmic neuronal activities, often with different time scales. One example of such a network is the pre-Bötzinger complex respiratory network (preBötC) that can simultaneously generate fast, small-amplitude, monophasic eupneic breaths together with slow, high-amplitude, biphasic augmented breaths (sighs). However, the underlying rhythmogenic mechanisms for this bimodal discharge pattern remain unclear, leaving two possible explanations: the existence of either reconfiguring processes within the same network or two distinct subnetworks. Based on recent in vitro data obtained in the mouse embryo, we have built a computational model consisting of two compartments, interconnected through appropriate synapses. One compartment generates sighs and the other produces eupneic bursts. The model reproduces basic features of simultaneous sigh and eupnea generation (two types of bursts differing in terms of shape, amplitude, and frequency of occurrence) and mimics the effect of blocking glycinergic synapses. Furthermore, we used this model to make predictions that were subsequently tested on the isolated preBötC in mouse brainstem slice preparations. Through a combination of in vitro and in silico approaches we find that (1) sigh events are less sensitive to network excitability than eupneic activity, (2) calcium-dependent mechanisms and the I_h current play a prominent role in sigh generation, and (3) specific parameters of I_h activation set the low sensitivity to excitability in the sigh neuronal subset. Altogether, our results strongly support the hypothesis that distinct subpopulations within the preBötC network are responsible for sigh and eupnea rhythmogenesis.

Identification of the pre-Bötzinger complex inspiratory center in calibrated "sandwich" slices from newborn mice with fluorescent Dbx1 interneurons

Inspiratory active pre‐Bötzinger complex (preBötC) networks produce the neural rhythm that initiates and controls breathing movements. We previously identified the preBötC in the newborn rat brainstem and established anatomically defined transverse slices in which the preBötC remains active when exposed at one surface. This follow‐up study uses a neonatal mouse model in which the preBötC as well as a genetically defined class of respiratory interneurons can be identified and selectively targeted for physiological recordings. The population of glutamatergic interneurons whose precursors express the transcription factor Dbx1 putatively comprises the core respiratory rhythmogenic circuit. Here, we used intersectional mouse genetics to identify the brainstem distribution of Dbx1‐derived neurons in the context of observable respiratory marker structures. This reference brainstem atlas enabled online histology for generating calibrated sandwich slices to identify the preBötC location, which was heretofore unspecified for perinatal mice. Sensitivity to opioids ensured that slice rhythms originated from preBötC neurons and not parafacial respiratory group/retrotrapezoid nucleus (pFRG/RTN) cells because opioids depress preBötC, but not pFRG/RTN rhythms. We found that the preBötC is centered ~0.4 mm caudal to the facial motor nucleus in this Cre/lox reporter mouse during postnatal days 0–4. Our findings provide the essential basis for future optically guided electrophysiological and fluorescence imaging‐based studies, as well as the application of other Cre‐dependent tools to record or manipulate respiratory rhythmogenic neurons. These resources will ultimately help elucidate the mechanisms that promote respiratory‐related oscillations of preBötC Dbx1‐derived neurons and thus breathing.

Breathing movements emanate from Dbx1‐derived interneurons of the brainstem pre‐Bötzinger complex (preBötC). We generated a histology atlas of the medulla in newborn Dbx1 Cre/lox reporter mice and performed physiological tests to pinpoint the preBötC location and map the Dbx1 neuron distribution, which will facilitate neurobiological studies of respiratory rhythm generation.

Laser ablation of Dbx1 neurons in the pre-Bötzinger complex stops inspiratory rhythm and impairs output in neonatal mice

To understand the neural origins of rhythmic behavior one must characterize the central pattern generator circuit and quantify the population size needed to sustain functionality. Breathing-related interneurons of the brainstem pre-Bötzinger complex (preBötC) that putatively comprise the core respiratory rhythm generator in mammals are derived from Dbx1-expressing precursors. Here, we show that selective photonic destruction of Dbx1 preBötC neurons in neonatal mouse slices impairs respiratory rhythm but surprisingly also the magnitude of motor output; respiratory hypoglossal nerve discharge decreased and its frequency steadily diminished until rhythm stopped irreversibly after 85±20 (mean ± SEM) cellular ablations, which corresponds to ∼15% of the estimated population. These results demonstrate that a single canonical interneuron class generates respiratory rhythm and contributes in a premotor capacity, whereas these functions are normally attributed to discrete populations. We also establish quantitative cellular parameters that govern network viability, which may have ramifications for respiratory pathology in disease states.

Atoh1-dependent rhombic lip neurons are required for temporal delay between independent respiratory oscillators in embryonic mice

All motor behaviors require precise temporal coordination of different muscle groups. Breathing, for example, involves the sequential activation of numerous muscles hypothesized to be driven by a primary respiratory oscillator, the preBötzinger Complex, and at least one other as-yet unidentified rhythmogenic population. We tested the roles of Atoh1-, Phox2b-, and Dbx1-derived neurons (three groups that have known roles in respiration) in the generation and coordination of respiratory output. We found that Dbx1-derived neurons are necessary for all respiratory behaviors, whereas independent but coupled respiratory rhythms persist from at least three different motor pools after eliminating or silencing Phox2b- or Atoh1-expressing hindbrain neurons. Without Atoh1 neurons, however, the motor pools become temporally disorganized and coupling between independent respiratory oscillators decreases. We propose Atoh1 neurons tune the sequential activation of independent oscillators essential for the fine control of different muscles during breathing.DOI: http://dx.doi.org/10.7554/eLife.02265.001.

Morphological characterization of respiratory neurons in the pre-Bötzinger complex

Although the pre-Bötzinger complex (preBötC) was defined as the inspiratory rhythm generator long ago, the functional-anatomical characterization of its neuronal components is still being achieved. Recent advances have identified the expression of molecular markers in the preBötC neurons that, however, are not exclusive to specific respiratory neuron subtypes and have not always been related to specific cell morphologies. Here, we evaluated the morphology and the axonal projections of electrophysiologically defined respiratory neurons in the preBötC using whole-cell recordings and intracellular biocytin labeling. We found that respiratory pacemaker neurons are larger than expiratory neurons and that inspiratory neurons are smaller than pacemaker and expiratory neurons. Other morphological features such as somata shapes or dendritic branching patterns were not found to be significantly different among the preBötC neurons sampled. We also found that both pacemaker and inspiratory nonpacemaker neurons, but not expiratory neurons, show extensive axonal projections to the contralateral preBötC and show signs of electrical coupling. Overall, our data suggest that there are morphological differences between subtypes of preBötC respiratory neurons. It will be important to take such differences in consideration since morphological differences would influence synaptic responses and action potential propagation.

Flufenamic acid as an ion channel modulator

Flufenamic acid has been known since the 1960s to have anti-inflammatory properties attributable to the reduction of prostaglandin synthesis. Thirty years later, flufenamic acid appeared to be an ion channel modulator. Thus, while its use in medicine diminished, its use in ionic channel research expanded. Flufenamic acid commonly not only affects non-selective cation channels and chloride channels, but also modulates potassium, calcium and sodium channels with effective concentrations ranging from 10(-6)M in TRPM4 channel inhibition to 10(-3)M in two-pore outwardly rectifying potassium channel activation. Because flufenamic acid effects develop and reverse rapidly, it is a convenient and widely used tool. However, given the broad spectrum of its targets, experimental results have to be interpreted cautiously. Here we provide an overview of ion channels targeted by flufenamic acid to aid in interpreting its effects at the molecular, cellular, and system levels. If it is used with good practices, flufenamic acid remains a useful tool for ion channel research. Understanding the targets of FFA may help reevaluate its physiological impacts and revive interest in its therapeutic potential.