Computational models of the neural control of breathing

COVID-19 and silent hypoxemia in a minimal closed-loop model of the respiratory rhythm generator

Silent hypoxemia, or "happy hypoxia," is a puzzling phenomenon in which patients who have contracted COVID-19 exhibit very low oxygen saturation ( SaO 2 < 80%) but do not experience discomfort in breathing. The mechanism by which this blunted response to hypoxia occurs is unknown. We have previously shown that a computational model of the respiratory neural network (Diekman et al. in J Neurophysiol 118(4):2194-2215, 2017) can be used to test hypotheses focused on changes in chemosensory inputs to the central pattern generator (CPG). We hypothesize that altered chemosensory function at the level of the carotid bodies and/or the nucleus tractus solitarii are responsible for the blunted response to hypoxia. Here, we use our model to explore this hypothesis by altering the properties of the gain function representing oxygen sensing inputs to the CPG. We then vary other parameters in the model and show that oxygen carrying capacity is the most salient factor for producing silent hypoxemia. We call for clinicians to measure hematocrit as a clinical index of altered physiology in response to COVID-19 infection.

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.

Chemogenetic stimulation of phrenic motor output and diaphragm activity

Impaired diaphragm activation contributes to morbidity and mortality in many neurodegenerative diseases and neurologic injuries. We conducted experiments to determine if expression of an excitatory DREADD (designer receptors exclusively activation by designer drugs) in the mid-cervical spinal cord would enable respiratory-related activation of phrenic motoneurons to increase diaphragm activation. Wild type (C57/bl6) and ChAT-Cre mice received bilateral intraspinal (C4) injections of an adeno-associated virus (AAV) encoding the hM3D(Gq) excitatory DREADD. In wild type mice, this produced non-specific DREADD expression throughout the mid-cervical ventral horn. In ChAT-Cre mice, a Cre-dependent viral construct was used to drive DREADD expression in C4 ventral horn motoneurons, targeting the phrenic motoneuron pool. Diaphragm EMG was recorded during spontaneous breathing at 6-8 weeks post-AAV delivery. The selective DREADD ligand JHU37160 (J60) caused a bilateral, sustained (>1 hr) increase in inspiratory EMG bursting in both groups; the relative increase was greater in ChAT-Cre mice. Additional experiments in a ChAT-Cre rat model were conducted to determine if spinal DREADD activation could increase inspiratory tidal volume (VT) during spontaneous breathing without anesthesia. Three to four months after intraspinal (C4) injection of AAV driving Cre-dependent hM3D(Gq) expression, intravenous J60 resulted in a sustained (>30 min) increase in VT assessed using whole-body plethysmography. Subsequently, direct nerve recordings confirmed that J60 evoked a >50% increase in inspiratory phrenic output. The data show that mid-cervical spinal DREADD expression targeting the phrenic motoneuron pool enables ligand-induced, sustained increases in the neural drive to the diaphragm. Further development of this technology may enable application to clinical conditions associated with impaired diaphragm activation and hypoventilation.

Exploring the role of the Kölliker-Fuse nucleus in breathing variability by mathematical modelling

The Kölliker-Fuse nucleus (KF), which is part of the parabrachial complex, participates in the generation of eupnoea under resting conditions and the control of active abdominal expiration when increased ventilation is required. Moreover, dysfunctions in KF neuronal activity are believed to play a role in the emergence of respiratory abnormalities seen in Rett syndrome (RTT), a progressive neurodevelopmental disorder associated with an irregular breathing pattern and frequent apnoeas. Relatively little is known, however, about the intrinsic dynamics of neurons within the KF and how their synaptic connections affect breathing pattern control and contribute to breathing irregularities. In this study, we use a reduced computational model to consider several dynamical regimes of KF activity paired with different input sources to determine which combinations are compatible with known experimental observations. We further build on these findings to identify possible interactions between the KF and other components of the respiratory neural circuitry. Specifically, we present two models that both simulate eupnoeic as well as RTT-like breathing phenotypes. Using nullcline analysis, we identify the types of inhibitory inputs to the KF leading to RTT-like respiratory patterns and suggest possible KF local circuit organizations. When the identified properties are present, the two models also exhibit quantal acceleration of late-expiratory activity, a hallmark of active expiration featuring forced exhalation, with increasing inhibition to KF, as reported experimentally. Hence, these models instantiate plausible hypotheses about possible KF dynamics and forms of local network interactions, thus providing a general framework as well as specific predictions for future experimental testing. KEY POINTS: The Kölliker-Fuse nucleus (KF), a part of the parabrachial complex, is involved in regulating normal breathing and controlling active abdominal expiration during increased ventilation. Dysfunction in KF neuronal activity is thought to contribute to respiratory abnormalities seen in Rett syndrome (RTT). This study utilizes computational modelling to explore different dynamical regimes of KF activity and their compatibility with experimental observations. By analysing different model configurations, the study identifies inhibitory inputs to the KF that lead to RTT-like respiratory patterns and proposes potential KF local circuit organizations. Two models are presented that simulate both normal breathing and RTT-like breathing patterns. These models provide testable hypotheses and specific predictions for future experimental investigations, offering a general framework for understanding KF dynamics and potential network interactions.

The Modulation by Anesthetics and Analgesics of Respiratory Rhythm in the Nervous System

Rhythmic eupneic breathing in mammals depends on the coordinated activities of the neural system that sends cranial and spinal motor outputs to respiratory muscles. These outputs modulate lung ventilation and adjust respiratory airflow, which depends on the upper airway patency and ventilatory musculature. Anesthetics are widely used in clinical practice worldwide. In addition to clinically necessary pharmacological effects, respiratory depression is a critical side effect induced by most general anesthetics. Therefore, understanding how general anesthetics modulate the respiratory system is important for the development of safer general anesthetics. Currently used volatile anesthetics and most intravenous anesthetics induce inhibitory effects on respiratory outputs. Various general anesthetics produce differential effects on respiratory characteristics, including the respiratory rate, tidal volume, airway resistance, and ventilatory response. At the cellular and molecular levels, the mechanisms underlying anesthetic-induced breathing depression mainly include modulation of synaptic transmission of ligand-gated ionotropic receptors (e.g., γ-aminobutyric acid, N-methyl-D-aspartate, and nicotinic acetylcholine receptors) and ion channels (e.g., voltage-gated sodium, calcium, and potassium channels, two-pore domain potassium channels, and sodium leak channels), which affect neuronal firing in brainstem respiratory and peripheral chemoreceptor areas. The present review comprehensively summarizes the modulation of the respiratory system by clinically used general anesthetics, including the effects at the molecular, cellular, anatomic, and behavioral levels. Specifically, analgesics, such as opioids, which cause respiratory depression and the "opioid crisis", are discussed. Finally, underlying strategies of respiratory stimulation that target general anesthetics and/or analgesics are summarized.

Parabrachial tachykinin1-expressing neurons involved in state-dependent breathing control

Breathing is regulated automatically by neural circuits in the medulla to maintain homeostasis, but breathing is also modified by behavior and emotion. Mice have rapid breathing patterns that are unique to the awake state and distinct from those driven by automatic reflexes. Activation of medullary neurons that control automatic breathing does not reproduce these rapid breathing patterns. By manipulating transcriptionally defined neurons in the parabrachial nucleus, we identify a subset of neurons that express the Tac1, but not Calca, gene that exerts potent and precise conditional control of breathing in the awake, but not anesthetized, state via projections to the ventral intermediate reticular zone of the medulla. Activating these neurons drives breathing to frequencies that match the physiological maximum through mechanisms that differ from those that underlie the automatic control of breathing. We postulate that this circuit is important for the integration of breathing with state-dependent behaviors and emotions.

Respiratory rhythm and pattern generation: Brainstem cellular and circuit mechanisms

Breathing movements in mammals are driven by rhythmic neural activity automatically generated within spatially and functionally organized brainstem neural circuits comprising the respiratory central pattern generator (CPG). This chapter reviews up-to-date experimental information and theoretical studies of the cellular and circuit mechanisms of respiratory rhythm and pattern generation operating within critical components of this CPG in the lower brainstem. Over the past several decades, there have been substantial advances in delineating the spatial architecture of essential medullary regions and their regional cellular and circuit properties required to understand rhythm and pattern generation mechanisms. A fundamental concept is that the circuits in these regions have rhythm-generating capabilities at multiple cellular and circuit organization levels. The regional cellular properties, circuit organization, and control mechanisms allow flexible expression of neural activity patterns for a repertoire of respiratory behaviors under various physiologic conditions that are dictated by requirements for homeostatic regulation and behavioral integration. Many mechanistic insights have been provided by computational modeling studies driven by experimental results and have advanced understanding in the field. These conceptual and theoretical developments are discussed.

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.

A Damaged-Informed Lung Ventilator Model for Ventilator Waveforms

Motivated by a desire to understand pulmonary physiology, scientists have developed physiological lung models of varying complexity. However, pathophysiology and interactions between human lungs and ventilators, e.g., ventilator-induced lung injury (VILI), present challenges for modeling efforts. This is because the real-world pressure and volume signals may be too complex for simple models to capture, and while complex models tend not to be estimable with clinical data, limiting clinical utility. To address this gap, in this manuscript we developed a new damaged-informed lung ventilator (DILV) model. This approach relies on mathematizing ventilator pressure and volume waveforms, including lung physiology, mechanical ventilation, and their interaction. The model begins with nominal waveforms and adds limited, clinically relevant, hypothesis-driven features to the waveform corresponding to pulmonary pathophysiology, patient-ventilator interaction, and ventilator settings. The DILV model parameters uniquely and reliably recapitulate these features while having enough flexibility to reproduce commonly observed variability in clinical (human) and laboratory (mouse) waveform data. We evaluate the proof-in-principle capabilities of our modeling approach by estimating 399 breaths collected for differently damaged lungs for tightly controlled measurements in mice and uncontrolled human intensive care unit data in the absence and presence of ventilator dyssynchrony. The cumulative value of mean squares error for the DILV model is, on average, ≈12 times less than the single compartment lung model for all the waveforms considered. Moreover, changes in the estimated parameters correctly correlate with known measures of lung physiology, including lung compliance as a baseline evaluation. Our long-term goal is to use the DILV model for clinical monitoring and research studies by providing high fidelity estimates of lung state and sources of VILI with an end goal of improving management of VILI and acute respiratory distress syndrome.

Multi-Level Regulation of Opioid-Induced Respiratory Depression

Opioids depress minute ventilation primarily by reducing respiratory rate. This results from direct effects on the preBötzinger Complex as well as from depression of the Parabrachial/Kölliker-Fuse Complex, which provides excitatory drive to preBötzinger Complex neurons mediating respiratory phase-switch. Opioids also depress awake drive from the forebrain and chemodrive.

Methods for frequency and correlation analyses of neurograms

Many physiological functions are based on motor rhythmic activities, among them breathing is a vital issue. The method presented here, or 'temporal grid extraction', aims at characterizing the temporal organization of such an activity. Beyond the measurement of the fundamental frequency, defining the successive cycles, some signal processing tools are helpful in order to look for the presence of higher frequency components that potentially structure these cycles. The method is applied to neurograms recorded from frog brainstem preparations, where two cycle types, buccal and lung cycles, may alternate. It relies on:•Continues Wavelet Transform (CWT) for time-frequency maps and frequency profiles•Crosscorrelation analysis for amplitude maps and amplitude profiles•Cycle-by-cycle autocorrelation analysis for autocorrelation maps and autocorrelation profilesUsing this method, the maps and profiles have revealed that a common high frequency clock drives both buccal and lung cycles.

Inspiratory Off-Switch Mediated by Optogenetic Activation of Inhibitory Neurons in the preBötzinger Complex In Vivo

The role of inhibitory neurons in the respiratory network is a matter of ongoing debate. Conflicting and contradicting results are manifold and the question whether inhibitory neurons are essential for the generation of the respiratory rhythm as such is controversial. Inhibitory neurons are required in pulmonary reflexes for adapting the activity of the central respiratory network to the status of the lung and it is hypothesized that glycinergic neurons mediate the inspiratory off-switch. Over the years, optogenetic tools have been developed that allow for cell-specific activation of subsets of neurons in vitro and in vivo. In this study, we aimed to identify the effect of activation of inhibitory neurons in vivo. Here, we used a conditional transgenic mouse line that expresses Channelrhodopsin 2 in inhibitory neurons. A 200 µm multimode optical fiber ferrule was implanted in adult mice using stereotaxic surgery, allowing us to stimulate inhibitory, respiratory neurons within the core excitatory network in the preBötzinger complex of the ventrolateral medulla. We show that, in anesthetized mice, activation of inhibitory neurons by blue light (470 nm) continuously or with stimulation frequencies above 10 Hz results in a significant reduction of the respiratory rate, in some cases leading to complete cessation of breathing. However, a lower stimulation frequency (4–5 Hz) could induce a significant increase in the respiratory rate. This phenomenon can be explained by the resetting of the respiratory cycle, since stimulation during inspiration shortened the associated breath and thereby increased the respiratory rate, while stimulation during the expiratory interval reduced the respiratory rate. Taken together, these results support the concept that activation of inhibitory neurons mediates phase-switching by inhibiting excitatory rhythmogenic neurons in the preBötzinger complex.

Neural Interactions in Developing Rhythmogenic Spinal Networks: Insights From Computational Modeling

The mechanisms involved in generation of rhythmic locomotor activity in the mammalian spinal cord remain poorly understood. These mechanisms supposedly rely on both intrinsic properties of constituting neurons and interactions between them. A subset of Shox2 neurons was suggested to contribute to generation of spinal locomotor activity, but the possible cellular basis for rhythmic bursting in these neurons remains unknown. Ha and Dougherty (2018) recently revealed the presence of bidirectional electrical coupling between Shox2 neurons in neonatal spinal cords, which can be critically involved in neuronal synchronization and generation of populational bursting. Gap junctional connections found between functionally-related Shox2 interneurons decrease with age, possibly being replaced by increasing interactions through chemical synapses. Here, we developed a computational model of a heterogeneous population of neurons sparsely connected by electrical or/and chemical synapses and investigated the dependence of frequency of populational bursting on the type and strength of neuronal interconnections. The model proposes a mechanistic explanation that can account for the emergence of a synchronized rhythmic activity in the neuronal population and provides insights into the possible role of gap junctional coupling between Shox2 neurons in the spinal mechanisms for locomotor rhythm generation.

On the Organization of the Locomotor CPG: Insights From Split-Belt Locomotion and Mathematical Modeling

Rhythmic limb movements during locomotion are controlled by central pattern generator (CPG) circuits located in the spinal cord. It is considered that these circuits are composed of individual rhythm generators (RGs) for each limb interacting with each other through multiple commissural and long propriospinal circuits. The organization and operation of each RG are not fully understood, and different competing theories exist about interactions between its flexor and extensor components, as well as about left-right commissural interactions between the RGs. The central idea of circuit organization proposed in this study is that with an increase of excitatory input to each RG (or an increase in locomotor speed) the rhythmogenic mechanism of the RGs changes from "flexor-driven" rhythmicity to a "classical half-center" mechanism. We test this hypothesis using our experimental data on changes in duration of stance and swing phases in the intact and spinal cats walking on the ground or tied-belt treadmills (symmetric conditions) or split-belt treadmills with different left and right belt speeds (asymmetric conditions). We compare these experimental data with the results of mathematical modeling, in which simulated CPG circuits operate in similar symmetric and asymmetric conditions with matching or differing control drives to the left and right RGs. The obtained results support the proposed concept of state-dependent changes in RG operation and specific commissural interactions between the RGs. The performed simulations and mathematical analysis of model operation under different conditions provide new insights into CPG network organization and limb coordination during locomotion.

Inhibitory control of active expiration by the Bötzinger complex in rats

KEY POINTS

Contraction of abdominal muscles at the end of expiration during metabolic challenges (such as hypercapnia and hypoxia) improves pulmonary ventilation. The emergence of this active expiratory pattern requires the recruitment of the expiratory oscillator located on the ventral surface of the medulla oblongata. Here we show that an inhibitory circuitry located in the Bötzinger complex is an important source of inhibitory drive to the expiratory oscillator. This circuitry, mediated by GABAergic and glycinergic synapses, provides expiratory inhibition that restrains the expiratory oscillator under resting condition and regulates the formation of abdominal expiratory activity during active expiration. By combining experimental and modelling approaches, we propose the organization and connections within the respiratory network that control the changes in the breathing pattern associated with elevated metabolic demand.

ABSTRACT

The expiratory neurons of the Bötzinger complex (BötC) provide inhibitory inputs to the respiratory network, which, during eupnoea, are critically important for respiratory phase transition and duration control. Here, we investigated how the BötC neurons interact with the expiratory oscillator located in the parafacial respiratory group (pFRG) and control the abdominal activity during active expiration. Using the decerebrated, arterially perfused in situ preparations of juvenile rats, we recorded the activity of expiratory neurons and performed pharmacological manipulations of the BötC and pFRG during hypercapnia or after the exposure to short-term sustained hypoxia - conditions that generate active expiration. The experimental data were integrated in a mathematical model to gain new insights into the inhibitory connectome within the respiratory central pattern generator. Our results indicate that the BötC neurons may establish mutual connections with the pFRG, providing expiratory inhibition during the first stage of expiration and receiving excitatory inputs during late expiration. Moreover, we found that application of GABAergic and glycinergic antagonists in the BötC caused opposing effects on abdominal expiratory activity, suggesting complex inhibitory circuitry within the BötC. Using mathematical modelling, we propose that the BötC network organization and its interactions with the pFRG restrain abdominal activity under resting conditions and contribute to abdominal expiratory pattern formation during active expiration observed during hypercapnia or after the exposure to short-term sustained hypoxia.

Effect of voluntary breathing exercises on stable coronary artery disease in heart rate variability and rate-pressure product: a study protocol for a single-blind, prospective, randomized controlled trial

Background

At present, China has more than 11 million patients with stable coronary heart disease and this is becoming a major public health problem. The pathological changes of coronary heart disease can lead to dysfunction of the cardiac autonomic nervous system, which increases the risk of complications such as malignant arrhythmia (ventricular flutter, ventricular fibrillation, etc.), heart rate, systolic blood pressure, and rate-pressure product (RPP), which is highly correlated with myocardial oxygen consumption and indirectly reflects myocardial blood supply and oxygen consumption. Although the guidelines recommend that such patients take drugs to reduce heart rate and myocardial oxygen consumption, the clinical control of heart rate is still not ideal. Thus, in this trial, we will use voluntary breathing exercises as the strategy of exercise rehabilitation for patients with stable coronary artery disease (SCAD), in order to increase the vagus nerve activity and/or reduce the sympathetic nervous activity, help maintain or rebuild the balance of plant nerve system, improve the time-domain index of heart rate variability, reduce the burden on the heart, and relieve patients’ anxiety and other negative emotions.

Methods

This is a 6-month single-blind, randomized controlled clinical trial that will be conducted in the First Affiliated Hospital of Soochow University. A total of 140 patients who fill out the Informed Consent Form are registered and randomized 1:1 into the Voluntary Breathing Exercises (VBE)-based clinical trial monitoring group (n = 70) or the Routine follow-up group (n = 70). The VBE-based clinical trial monitoring group is given VBE training on the basis of conventional treatment and health education, while the control group received conventional health education and follow-up. The primary outcomes will be measured heart rate variability and RPP. Secondary outcomes will include changes in Self-rating Anxiety Scale, total cholesterol, triglyceride, high-density lipoprotein, low-density lipoprotein, weight, and body mass index.

Discussion

This trial will carry out scientific respiratory exercise for patients with SCAD, which belongs to the category of active secondary prevention for patients, and changes from remedial to pre-protective. VBE is easy to operate and is not limited by time and place. It is important and meaningful to carry out VBE for patients with SCAD. This study will provide considerable evidence for further large-scale trials and alternative strategies for the rehabilitation nursing of patients with SCAD.

Trial registration

Chinese Clinical Trials Registry, 1900024043. Registered on 23 June 2019.

Midbrain control of breathing and blood pressure: The role of periaqueductal gray matter and mesencephalic collicular neuronal microcircuit oscillators

Neural circuitry residing within the medullary ventral respiratory column nuclei and dorsal respiratory group interact with the Kölliker-Fuse and medial parabrachial nuclei to generate the core breathing rhythm and pattern during resting conditions. Triphasic eupnea consists of inspiratory [I], post-inspiratory [post-I], and late-expiratory [E2] phases. Mesencephalic zones exert modulatory influences upon respiratory rhythm-generating circuitry, sympathetic oscillators, and parasympathetic nuclei. The earliest evidence supporting the existence of midbrain control of breathing derives from studies conducted by Martin and Booker in 1878. These authors demonstrated electrical stimulation of the deep layers of the mesencephalic colliculi in the rabbit augmented ventilation and sequentially elicited chest wall tremors and tetany. Investigations performed during the past several decades would demonstrate stimlation of distributed zones within the midbrain reticular formation elicits starkly disparate effects upon respiratory phase switching. Schmid, Böhmer, and Fallert demonstrated electrical stimulation of the nucleus rubre and emanating axon bundles alternately elicits or inhibits the activity of medullary expiratory- or inspiratory-related units and phrenic nerve discharge with differential latency. A series of studies would later indicate the red nucleus mediates hypoxic ventilatory depression. Periaqueductal gray matter neurons exhibit extensive afferent and efferent interconnectivity with suprabulbar, brainstem, and spinal cord zones aptly positioning these units to modulate breathing, autonomic outflow, nociception locomotion, micturtion, and sexual behavior. Experimental stimulatory activation of the tectal colliculi and periaqueductal gray matter via electrical current or glutamate, D,L-homocysteinic acid, or bicuculline microinjections coordinately modulates neuromotor inspiratory bursting frequency and amplitude and discharge of pre-Bötzinger complex, ventrolateral medullary late-I and post-I, and ventrolateral nucleus tractus solitarius decrementing early-I and augmenting and decrementing late-I neurons, elicits expiratory outflow and vocalization, and blunt the Hering-Breuer reflex in unanesthetzed decerebrate and anesthetized preprations of the cat and rat. Stimulation of the mesencephalic colliuli or dorsal divisions of the PAG potently amplifes renal sympathetic neural efferent activity, dynamic arterial pressure magnitude, and myocardial contraction frequency and elicits various behavioral defense responses. Elicited physiological effects exhibit extensive locoregional heterogeneity and variably enlist requisite contributions from the dorsomedial hypothalamus and/or lateral parabrachial nuclei. Stimulation of the dorsal mesencephalon occasionally elicits dynamic increases of arterial pressure magnitude exhibiting prominent oscillatory variability coherent with phrenic nerve discharge, perhaps by generating intra-neuraxial hysteresis, serving to intermittently deliver blood to organ vascular beds under high pressure in order to prevent organ edema, microcirculatory dysfunction, and downregulation of vascular smooth muscle alpha adrenergic receptors. Chemosensitive mesencephalic caudal raphé units and projections of hypoxia-sensitive units in the caudal hypothalamus to the periaqueductal gray matter may imply the existence of a diencephalo-smesencephalic chemosensitive network modulating breathing and sympathetic discharge.

Complex bursting dynamics in an embryonic respiratory neuron model

Pre-Bötzinger complex (pre-BötC) network activity within the mammalian brainstem controls the inspiratory phase of the respiratory rhythm. While bursting in pre-BötC neurons during the postnatal period has been extensively studied, less is known regarding inspiratory pacemaker neuron behavior at embryonic stages. Recent data in mouse embryo brainstem slices have revealed the existence of a variety of bursting activity patterns depending on distinct combinations of burst-generating I_NaP and I_CAN conductances. In this work, we consider a model of an isolated embryonic pre-BötC neuron featuring two distinct bursting mechanisms. We use methods of dynamical systems theory, such as phase plane analysis, fast-slow decomposition, and bifurcation analysis, to uncover mechanisms underlying several different types of intrinsic bursting dynamics observed experimentally including several forms of plateau bursts, bursts involving depolarization block, and various combinations of these patterns. Our analysis also yields predictions about how changes in the balance of the two bursting mechanisms contribute to alterations in an inspiratory pacemaker neuron activity during prenatal development.

Endogenous glutamatergic inputs to the Parabrachial Nucleus/Kölliker-Fuse Complex determine respiratory rate

The Kölliker-Fuse Nucleus (KF) has been widely investigated for its contribution to "inspiratory off-switch" while more recent studies showed that activation of the Parabrachial Nucleus (PBN) shortened expiratory duration. This study used an adult, in vivo, decerebrate rabbit model to delineate the contribution of each site to inspiratory and expiratory duration through sequential block of glutamatergic excitation with the receptor antagonists 2,3-dihydroxy-6-nitro-7-sulfamoyl-benzo[f]quinoxaline-2,3-dione (NBQX) and d(-)-2-amino-5-phosphonopentanoic acid (AP5). Glutamatergic disfacilitation caused large increases in inspiratory and expiratory duration and minor decrease in peak phrenic activity (PPA). Hypoxia only partially reversed respiratory rate depression but PPA was increased to >200 % of control. The contribution of PBN activity to inspiratory and expiratory duration was equal while block of the KF affected inspiratory duration more than expiratory. We conclude that in the in vivo preparation respiratory rate greatly depends on PBN/KF activity, which contributes to the "inspiratory on- "and "off-switch", but is of minor importance for the magnitude of phrenic motor output.

Kölliker-Fuse/Parabrachial complex mu opioid receptors contribute to fentanyl-induced apnea and respiratory rate depression

Overdoses caused by the opioid agonist fentanyl have increased exponentially in recent years. Identifying mechanisms to counter progression to fatal respiratory apnea during opioid overdose is desirable, but difficult to study in vivo. The pontine Kölliker-Fuse/Parabrachial complex (KF/PB) provides respiratory drive and contains opioid-sensitive neurons. The contribution of the KF/PB complex to fentanyl-induced apnea was investigated using the in situ arterially perfused preparation of rat. Systemic application of fentanyl resulted in concentration-dependent respiratory disturbances. At low concentrations, respiratory rate slowed and subsequently transitioned to an apneustic-like, 2-phase pattern. Higher concentrations caused prolonged apnea, interrupted by occasional apneustic-like bursts. Application of CTAP, a selective mu opioid receptor antagonist, directly into the KF/PB complex reversed and prevented fentanyl-induced apnea by increasing the frequency of apneustic-like bursting. These results demonstrate that countering opioid effects in the KF/PB complex is sufficient to restore phasic respiratory output at a rate similar to pre-fentanyl conditions, which could be beneficial in opioid overdose.

Differential impact of two critical respiratory centres in opioid-induced respiratory depression in awake mice

KEY POINTS

The main cause of death from opioid overdose is respiratory depression due to the activation of µ-opioid receptors (MORs). We conditionally deleted MORs from neurons in two key areas of the brainstem respiratory circuitry (the Kölliker-Fuse nucleus (KF) and pre-Bötzinger complex (preBötC)) to determine their role in opioid-induced respiratory disturbances in adult, awake mice. Deletion of MORs from KF neurons attenuated respiratory rate depression at all doses of morphine. Deletion of MORs from preBötC neurons attenuated rate depression at the low dose, but had no effect on rate following high doses of morphine. Instead, high doses of morphine increased the occurrence of apnoeas. The results indicate that opioids affect distributed key areas of the respiratory network in a dose-dependent manner and countering the respiratory effects of high dose opioids via the KF may be an effective approach to combat overdose.

ABSTRACT

The primary cause of death from opioid overdose is respiratory failure. High doses of opioids cause severe rate depression and increased risk of fatal apnoea, which correlate with increasing irregularities in breathing pattern. µ-Opioid receptors (MORs) are widely distributed throughout the brainstem respiratory network, but the mechanisms underlying respiratory depression are poorly understood. The medullary pre-Bötzinger complex (preBötC) and the pontine Kölliker-Fuse nucleus (KF) are considered critical for inducing opioid-related respiratory disturbances. We used a conditional knockout approach to investigate the roles and relative contribution of MORs in KF and preBötC neurons in opioid-induced respiratory depression in awake adult mice. The results revealed dose-dependent and region-specific opioid effects on the control of both respiratory rate and pattern. Respiratory depression induced by an anti-nociceptive dose of morphine was significantly attenuated following deletion of MORs from either the KF or the preBötC, suggesting cumulative network effects on respiratory rate control at low opioid doses. Deletion of MORs from KF neurons also relieved rate depression at near-maximal respiratory depressant doses of morphine. Meanwhile, deletion of MORs from the preBötC had no effect on rate following administration of high doses of morphine. Instead, a severe ataxic breathing pattern emerged with many apnoeas. We conclude that opioids affect distributed areas of the respiratory network and opioid-induced respiratory depression cannot be attributed to only one area in isolation. However, countering the effects of near maximal respiratory depressant doses of opioids in the KF may be a powerful approach to combat opioid overdose.

Effects of persistent sodium current blockade in respiratory circuits depend on the pharmacological mechanism of action and network dynamics

The mechanism(s) of action of most commonly used pharmacological blockers of voltage-gated ion channels are well understood; however, this knowledge is rarely considered when interpreting experimental data. Effects of blockade are often assumed to be equivalent, regardless of the mechanism of the blocker involved. Using computer simulations, we demonstrate that this assumption may not always be correct. We simulate the blockade of a persistent sodium current (I_NaP), proposed to underlie rhythm generation in pre-Bötzinger complex (pre-BötC) respiratory neurons, via two distinct pharmacological mechanisms: (1) pore obstruction mediated by tetrodotoxin and (2) altered inactivation dynamics mediated by riluzole. The reported effects of experimental application of tetrodotoxin and riluzole in respiratory circuits are diverse and seemingly contradictory and have led to considerable debate within the field as to the specific role of I_NaP in respiratory circuits. The results of our simulations match a wide array of experimental data spanning from the level of isolated pre-BötC neurons to the level of the intact respiratory network and also generate a series of experimentally testable predictions. Specifically, in this study we: (1) provide a mechanistic explanation for seemingly contradictory experimental results from in vitro studies of I_NaP block, (2) show that the effects of I_NaP block in in vitro preparations are not necessarily equivalent to those in more intact preparations, (3) demonstrate and explain why riluzole application may fail to effectively block I_NaP in the intact respiratory network, and (4) derive the prediction that effective block of I_NaP by low concentration tetrodotoxin will stop respiratory rhythm generation in the intact respiratory network. These simulations support a critical role for I_NaP in respiratory rhythmogenesis in vivo and illustrate the importance of considering mechanism when interpreting and simulating data relating to pharmacological blockade.

The application of pharmacological agents that affect transmembrane ionic currents in neurons is a commonly used experimental technique. A simplistic interpretation of experiments involving these agents suggests that antagonist application removes the impacted current and that subsequently observed changes in activity are attributable to the loss of that current’s effects. The more complex reality, however, is that different drugs may have distinct mechanisms of action, some corresponding not to a removal of a current but rather to a changing of its properties. We use computational modeling to explore the implications of the distinct mechanisms associated with two drugs, riluzole and tetrodotoxin, that are often characterized as sodium channel blockers. Through this approach, we offer potential explanations for disparate findings observed in experiments on neural respiratory circuits and show that the experimental results are consistent with a key role for the persistent sodium current in respiratory rhythm generation.

Role of fast inhibitory synaptic transmission in neonatal respiratory rhythmogenesis and pattern formation

Several studies have investigated the general role of chloride-based neurotransmission (GABA_A and glycinergic signaling) in respiratory rhythmogenesis and pattern formation. In several brain regions, developmental alterations in these signaling pathways have been shown to be mediated by changes in cation-chloride cotransporter (CC) expression. For instance, CC expression changes during the course of neonatal development in medullary respiratory nuclei and other brain/spinal cord regions in a manner which decreases the cellular import, and increases the export, of chloride ions, shifting reversal potentials for chloride to progressively more negative values with maturation. In slice preparations of the same, this is related to an excitatory-to-inhibitory shift of GABA_A- and glycinergic signaling. In medullary slices, GABA_A-/glycinergic signaling in the early neonatal period is excitatory, becoming inhibitory over time. Additionally, blockade of the Na⁺/K⁺/2Cl^- cotransporter, which imports these ions via secondary active transport, converts excitatory response to inhibitory ones. These effects have not yet been demonstrated at the individual respiratory-related neuron level to occur in intact (in vivo or in situ) animal preparations, which in contrast to slices, possess normal network connectivity and natural sources of tonic drive. Developmental changes in respiratory rhythm generating and pattern forming pontomedullary respiratory circuitry may contribute to critical periods, during which there exist increased risk for perinatal respiratory disturbances of central, obstructive, or hypoxia/hypercapnia-induced origin, including the sudden infant death syndrome. Thus, better characterizing the neurochemical maturation of the central respiratory network will enhance our understanding of these conditions, which will facilitate development of targeted therapies for respiratory disturbances in neonates and infants.

The logic behind neural control of breathing pattern

The respiratory rhythm generator is spectacular in its ability to support a wide range of activities and adapt to changing environmental conditions, yet its operating mechanisms remain elusive. We show how selective control of inspiration and expiration times can be achieved in a new representation of the neural system (called a Boolean network). The new framework enables us to predict the behavior of neural networks based on properties of neurons, not their values. Hence, it reveals the logic behind the neural mechanisms that control the breathing pattern. Our network mimics many features seen in the respiratory network such as the transition from a 3-phase to 2-phase to 1-phase rhythm, providing novel insights and new testable predictions.

Neural network model of an amphibian ventilatory central pattern generator

The neuronal multiunit model presented here is a formal model of the central pattern generator (CPG) of the amphibian ventilatory neural network, inspired by experimental data from Pelophylax ridibundus. The kernel of the CPG consists of three pacemakers and two follower neurons (buccal and lung respectively). This kernel is connected to a chain of excitatory and inhibitory neurons organized in loops. Simulations are performed with Izhikevich-type neurons. When driven by the buccal follower, the excitatory neurons transmit and reorganize the follower activity pattern along the chain, and when driven by the lung follower, the excitatory and inhibitory neurons of the chain fire in synchrony. The additive effects of synaptic inputs from the pacemakers on the buccal follower account for (1) the low frequency buccal rhythm, (2) the intra-burst high frequency oscillations, and (3) the episodic lung activity. Chemosensitivity to acidosis is implemented by an increase in the firing frequency of one of the pacemakers. This frequency increase leads to both a decrease in the buccal burst frequency and an increase in the lung episode frequency. The rhythmogenic properties of the model are robust against synaptic noise and pacemaker jitter. To validate the rhythm and pattern genesis of this formal CPG, neurograms were built from simulated motoneuron activity, and compared with experimental neurograms. The basic principles of our model account for several experimental observations, and we suggest that these principles may be generic for amphibian ventilation.

Pre- and post-inspiratory neurons change their firing properties in female rats exposed to chronic intermittent hypoxia

Obstructive sleep apnea patients face episodes of chronic intermittent hypoxia (CIH), which has been suggested as a causative factor for increased sympathetic activity (SNA) and hypertension. Female rats exposed to CIH develop hypertension and exhibit changes in respiratory-sympathetic coupling, marked by an increase in the inspiratory modulation of SNA. We tested the hypothesis that enhanced inspiratory-modulation of SNA is dependent on carotid bodies (CBs) and are associated with changes in respiratory network activity. For this, in CIH-female rats we evaluated the effect of CBs ablation on respiratory-sympathetic coupling, recorded from respiratory neurons in the working heart-brainstem preparation and from NTS neurons in brainstem slices. CIH-female rats had an increase in peripheral chemoreflex response and in spontaneous excitatory neurotransmission in NTS. CBs ablation prevents the increase in inspiratory modulation of SNA in CIH-female rats. Pre-inspiratory/inspiratory (Pre-I/I) neurons of CIH-female rats have a reduced firing frequency. Post-inspiratory neurons are active for a longer period during expiration in CIH-female rats. Further, using the computational model of a brainstem respiratory-sympathetic network, we demonstrate that a reduction in Pre-I/I neuron firing frequency simulates the enhanced inspiratory SNA modulation in CIH-female rats. We conclude that changes in respiratory-sympathetic coupling in CIH-female rats is dependent on CBs and it is associated with changes in firing properties of specific respiratory neurons types.

Reduced computational modelling of Kölliker-Fuse contributions to breathing patterns in Rett syndrome

KEY POINTS

Reduced computational models are used to test effects of loss of inhibition to the Kölliker-Fuse nucleus (KFn). Three reduced computational models that simulate eupnoeic and vagotomized respiratory rhythms are considered. All models exhibit the emergence of respiratory perturbations associated with Rett syndrome as inhibition to the KFn is diminished. Simulations suggest that application of 5-HT_1A agonists can mitigate the respiratory pathology. The three models can be distinguished and tested based on their predictions about connections and dynamics within the respiratory circuit and about effects of perturbations on certain respiratory neuron populations.

ABSTRACT

Rett syndrome (RTT) is a developmental disorder that can lead to respiratory disturbances featuring prolonged apnoeas of variable durations. Determining the mechanisms underlying these effects at the level of respiratory neural circuits would have significant implications for treatment efforts and would also enhance our understanding of respiratory rhythm generation and control. While experimental studies have suggested possible factors contributing to the respiratory patterns of RTT, we take a novel computational approach to the investigation of RTT, which allows for direct manipulation of selected system parameters and testing of specific hypotheses. Specifically, we present three reduced computational models, developed using an established framework, all of which successfully simulate respiratory outputs across eupnoeic and vagotomized conditions. All three models show that loss of inhibition to the Kölliker-Fuse nucleus reproduces the key respiratory alterations associated with RTT and, as suggested experimentally, that effects of 5-HT_1A agonists on the respiratory neural circuit suffice to alleviate this respiratory pathology. Each of the models makes distinct predictions regarding the neuronal populations and interactions underlying these effects, suggesting natural directions for future experimental testing.

A spatially dynamic network underlies the generation of inspiratory behaviors

Breathing is a vital rhythmic behavior that originates from neural networks within the brainstem. It is hypothesized that the breathing rhythm is generated by spatially distinct networks localized to discrete kernels or compartments. Here, we provide evidence that the functional boundaries of these compartments expand and contract dynamically based on behavioral or physiological demands. The ability of these rhythmic networks to change in size may allow the breathing rhythm to be very reliable, yet flexible enough to accommodate the large repertoire of mammalian behaviors that must be integrated with breathing.

The ability of neuronal networks to reconfigure is a key property underlying behavioral flexibility. Networks with recurrent topology are particularly prone to reconfiguration through changes in synaptic and intrinsic properties. Here, we explore spatial reconfiguration in the reticular networks of the medulla that generate breathing. Combined results from in vitro and in vivo approaches demonstrate that the network architecture underlying generation of the inspiratory phase of breathing is not static but can be spatially redistributed by shifts in the balance of excitatory and inhibitory network influences. These shifts in excitation/inhibition allow the size of the active network to expand and contract along a rostrocaudal medullary column during behavioral or metabolic challenges to breathing, such as changes in sensory feedback, sighing, and gasping. We postulate that the ability of this rhythm-generating network to spatially reconfigure contributes to the remarkable robustness and flexibility of breathing.

Carotid Bodies and the Integrated Cardiorespiratory Response to Hypoxia

Advances in our understanding of brain mechanisms for the hypoxic ventilatory response, coordinated changes in blood pressure, and the long-term consequences of chronic intermittent hypoxia as in sleep apnea, such as hypertension and heart failure, are giving impetus to the search for therapies to "erase" dysfunctional memories distributed in the carotid bodies and central nervous system. We review current network models, open questions, sex differences, and implications for translational research.

Firing properties of ventral medullary respiratory neurons in sino-aortic denervated rats

NEW FINDINGS

What is the central question of this study? After sino-aortic denervation (SAD), rats present normal levels of mean arterial pressure (MAP), high MAP variability and changes in breathing. However, mechanisms involved in SAD-induced respiratory changes and their impact on the modulation of sympathetic activity remain unclear. Herein, we characterized the firing frequency of medullary respiratory neurons after SAD. What is the main finding and its importance? Sino-aortic denervation-induced prolonged inspiration was associated with a reduced interburst frequency of pre-inspiratory/inspiratory neurons and an increased long-term variability of late inspiratory neurons, but no changes were observed in the ramp-inspiratory and post-inspiratory neurons. This imbalance in the respiratory network might contribute to the modulation of sympathetic activity after SAD.

ABSTRACT

In previous studies, we documented that after sino-aortic denervation (SAD) in rats there are significant changes in the breathing pattern, but no significant changes in sympathetic activity and mean arterial pressure compared with sham-operated rats. However, the neural mechanisms involved in the respiratory changes after SAD and the extent to which they might contribute to the observed normal sympathetic activity and mean arterial pressure remain unclear. Here, we hypothesized that after SAD, rats present with changes in the firing frequency of the ventral medullary inspiratory and post-inspiratory neurons. To test this hypothesis, male Wistar rats underwent SAD or sham surgery and 3 days later were surgically prepared for an in situ experiment. The duration of inspiration significantly increased in SAD rats. During inspiration, the total firing frequency of ramp-inspiratory, pre-inspiratory/inspiratory and late-inspiratory neurons was not different between groups. During post-inspiration, the total firing frequency of post-inspiratory neurons was also not different between groups. Furthermore, the data demonstrate a reduced interburst frequency of pre-inspiratory/inspiratory neurons and an increased long-term variability of late-inspiratory neurons in SAD compared with sham-operated rats. These findings indicate that the SAD-induced prolongation of inspiration was not accompanied by alterations in the total firing frequency of the ventral medullary respiratory neurons, but it was associated with changes in the long-term variability of late-inspiratory neurons. We suggest that the timing imbalance in the respiratory network in SAD rats might contribute to the modulation of presympathetic neurons after removal of baroreceptor afferents.

Interaction between the retrotrapezoid nucleus and the parafacial respiratory group to regulate active expiration and sympathetic activity in rats

The retrotrapezoid nucleus (RTN) contains chemosensitive cells that distribute CO₂-dependent excitatory drive to the respiratory network. This drive facilitates the function of the respiratory central pattern generator (rCPG) and increases sympathetic activity. It is also evidenced that during hypercapnia, the late-expiratory (late-E) oscillator in the parafacial respiratory group (pFRG) is activated and determines the emergence of active expiration. However, it remains unclear the microcircuitry responsible for the distribution of the excitatory signals to the pFRG and the rCPG in conditions of high CO₂. Herein, we hypothesized that excitatory inputs from chemosensitive neurons in the RTN are necessary for the activation of late-E neurons in the pFRG. Using the decerebrated in situ rat preparation, we found that lesions of neurokinin-1 receptor-expressing neurons in the RTN region with substance P-saporin conjugate suppressed the late-E activity in abdominal nerves (AbNs) and sympathetic nerves (SNs) and attenuated the increase in phrenic nerve (PN) activity induced by hypercapnia. On the other hand, kynurenic acid (100 mM) injections in the pFRG eliminated the late-E activity in AbN and thoracic SN but did not modify PN response during hypercapnia. Iontophoretic injections of retrograde tracer into the pFRG of adult rats revealed labeled phox2b-expressing neurons within the RTN. Our findings are supported by mathematical modeling of chemosensitive and late-E populations within the RTN and pFRG regions as two separate but interacting populations in a way that the activation of the pFRG late-E neurons during hypercapnia require glutamatergic inputs from the RTN neurons that intrinsically detect changes in CO₂/pH.

Walking with perturbations: a guide for biped humans and robots

This paper provides an update on the neural control of bipedal walking in relation to bioinspired models and robots. It is argued that most current models or robots are based on the construct of a symmetrical central pattern generator (CPG). However, new evidence suggests that CPG functioning is basically asymmetrical with its flexor half linked more tightly to the rhythm generator. The stability of bipedal gait, which is an important problem for robots and biological systems, is also addressed. While it is not possible to determine how biological biped systems guarantee stability, robot solutions can be useful to propose new hypotheses for biology. In the second part of this review, the focus is on gait perturbations, which is an important topic in robotics in view of the frequent falls of robots when faced with perturbations. From the human physiology it is known that the initial reaction often consists of a brief interruption followed by an adequate response. For instance, the successful recovery from a trip is achieved using some basic reactions (termed elevating and lowering strategies), that depend on the phase of the step cycle of the trip occurrence. Reactions to stepping unexpectedly in a hole depend on comparing expected and real feedback. Implementation of these ideas in models and robotics starts to emerge, with the most advanced robots being able to learn how to fall safely and how to deal with complicated disturbances such as provided by walking on a split-belt.

Deficient post-ictal cardiorespiratory compensatory mechanisms mediated by the periaqueductal gray may lead to death in a mouse model of SUDEP

Post-ictal cardiorespiratory failure is implicated as a major cause of sudden unexpected death in epilepsy (SUDEP) in patients. The DBA/1 mouse model of SUDEP is abnormally susceptible to fatal seizure-induced cardiorespiratory failure (S-CRF) induced by convulsant drug, hyperthermia, electroshock, and acoustic stimulation. Clinical and pre-clinical studies have implicated periaqueductal gray (PAG) abnormalities in SUDEP. Recent functional neuroimaging studies observed that S-CRF resulted in selective changes in PAG neuronal activity in DBA/1 mice. The PAG plays a critical compensatory role for respiratory distress caused by numerous physiological challenges in non-epileptic individuals. These observations suggest that abnormalities in PAG-mediated cardiorespiratory modulation may contribute to S-CRF in DBA/1 mice. To evaluate this, electrical stimulation (20 Hz, 20-100 μA, 10 s) was presented in the PAG of anesthetized DBA/1 and C57BL/6 (non-epileptic) control mice, and post-stimulus changes in respiration [inter-breath interval (IBI)] and heart rate variability (HRV) were examined. The post-stimulus period was considered analogous to the post-ictal period when S-CRF occurred in previous DBA/1 mouse studies. PAG stimulation caused significant intensity-related decreases in IBI in both mouse strains. However, this effect was significantly reduced in DBA/1 vis-a-vis C57BL/6 mice. These changes began immediately following cessation of stimulation and remained significant for 10 s. This time period is critical for initiating resuscitation to successfully prevent seizure-induced death in previous DBA/1 mouse experiments. Significant post-stimulus increases in HRV were also seen at ≥60 μA in the PAG in C57BL/6 mice, which were absent in DBA/1 mice. These data along with previous neuroimaging findings suggest that compensatory cardiorespiratory modulation mediated by PAG is deficient, which may be important to the susceptibility of DBA/1 mice to S-CRF. These observations suggest that correcting this deficit pharmacologically or by electrical stimulation may help to prevent S-CRF. These findings further support the potential importance of PAG abnormalities to human SUDEP.

Organization of the core respiratory network: Insights from optogenetic and modeling studies

The circuit organization within the mammalian brainstem respiratory network, specifically within and between the pre-Bötzinger (pre-BötC) and Bötzinger (BötC) complexes, and the roles of these circuits in respiratory pattern generation are continuously debated. We address these issues with a combination of optogenetic experiments and modeling studies. We used transgenic mice expressing channelrhodopsin-2 under the VGAT-promoter to investigate perturbations of respiratory circuit activity by site-specific photostimulation of inhibitory neurons within the pre-BötC or BötC. The stimulation effects were dependent on the intensity and phase of the photostimulation. Specifically: (1) Low intensity (≤ 1.0 mW) pulses delivered to the pre-BötC during inspiration did not terminate activity, whereas stronger stimulations (≥ 2.0 mW) terminated inspiration. (2) When the pre-BötC stimulation ended in or was applied during expiration, rebound activation of inspiration occurred after a fixed latency. (3) Relatively weak sustained stimulation (20 Hz, 0.5-2.0 mW) of pre-BötC inhibitory neurons increased respiratory frequency, while a further increase of stimulus intensity (> 3.0 mW) reduced frequency and finally (≥ 5.0 mW) terminated respiratory oscillations. (4) Single pulses (0.2-5.0 s) applied to the BötC inhibited rhythmic activity for the duration of the stimulation. (5) Sustained stimulation (20 Hz, 0.5-3.0 mW) of the BötC reduced respiratory frequency and finally led to apnea. We have revised our computational model of pre-BötC and BötC microcircuits by incorporating an additional population of post-inspiratory inhibitory neurons in the pre-BötC that interacts with other neurons in the network. This model was able to reproduce the above experimental findings as well as previously published results of optogenetic activation of pre-BötC or BötC neurons obtained by other laboratories. The proposed organization of pre-BötC and BötC circuits leads to testable predictions about their specific roles in respiratory pattern generation and provides important insights into key circuit interactions operating within brainstem respiratory networks.

Circuits in the rodent brainstem that control whisking in concert with other orofacial motor actions

The world view of rodents is largely determined by sensation on two length scales. One is within the animal's peri-personal space; sensorimotor control on this scale involves active movements of the nose, tongue, head, and vibrissa, along with sniffing to determine olfactory clues. The second scale involves the detection of more distant space through vision and audition; these detection processes also impact repositioning of the head, eyes, and ears. Here we focus on orofacial motor actions, primarily vibrissa-based touch but including nose twitching, head bobbing, and licking, that control sensation at short, peri-personal distances. The orofacial nuclei for control of the motor plants, as well as primary and secondary sensory nuclei associated with these motor actions, lie within the hindbrain. The current data support three themes: First, the position of the sensors is determined by the summation of two drive signals, i.e., a fast rhythmic component and an evolving orienting component. Second, the rhythmic component is coordinated across all orofacial motor actions and is phase-locked to sniffing as the animal explores. Reverse engineering reveals that the preBötzinger inspiratory complex provides the reset to the relevant premotor oscillators. Third, direct feedback from somatosensory trigeminal nuclei can rapidly alter motion of the sensors. This feedback is disynaptic and can be tuned by high-level inputs. A holistic model for the coordination of orofacial motor actions into behaviors will encompass feedback pathways through the midbrain and forebrain, as well as hindbrain areas.

Timescales and Mechanisms of Sigh-Like Bursting and Spiking in Models of Rhythmic Respiratory Neurons

Neural networks generate a variety of rhythmic activity patterns, often involving different timescales. One example arises in the respiratory network in the pre-Bötzinger complex of the mammalian brainstem, which can generate the eupneic rhythm associated with normal respiration as well as recurrent low-frequency, large-amplitude bursts associated with sighing. Two competing hypotheses have been proposed to explain sigh generation: the recruitment of a neuronal population distinct from the eupneic rhythm-generating subpopulation or the reconfiguration of activity within a single population. Here, we consider two recent computational models, one of which represents each of the hypotheses. We use methods of dynamical systems theory, such as fast-slow decomposition, averaging, and bifurcation analysis, to understand the multiple-timescale mechanisms underlying sigh generation in each model. In the course of our analysis, we discover that a third timescale is required to generate sighs in both models. Furthermore, we identify the similarities of the underlying mechanisms in the two models and the aspects in which they differ.