Respiratory rhythm generation in neonatal and adult mammals: the hybrid pacemaker-network model

Galaninergic and hypercapnia-activated neuronal projections to the ventral respiratory column

In mammals, the ventral respiratory column (VRC) plays a pivotal role in integrating neurochemically diverse inputs from brainstem and forebrain regions to generate respiratory motor patterns. VRC microinjection of the neuropeptide galanin has been reported to dampen carbon dioxide (CO2)-mediated chemoreflex responses. Additionally, we previously demonstrated that galaninergic neurons in the retrotrapezoid nucleus (RTN) are implicated in the adaptive response to hypercapnic stimuli, suggesting a link between RTN neuroplasticity and increased neuronal drive to the VRC. VRC neurons express galanin receptor 1, suggesting potential regulatory action by galanin, however, the precise galaninergic chemoreceptor-VRC circuitry remains to be determined. This study aimed to identify sources of galaninergic input to the VRC that contribute to central respiratory chemoreception. We employed a combination of retrograde neuronal tracing, in situ hybridisation and immunohistochemistry to investigate VRC-projecting neurons that synthesise galanin mRNA. In an additional series of experiments, we used acute hypercapnia exposure (10% CO2, 1 h) and c-Fos immunohistochemistry to ascertain which galaninergic nuclei projecting to the VRC are activated. Our findings reveal that a total of 30 brain nuclei and 51 subnuclei project to the VRC, with 12 of these containing galaninergic neurons, including the RTN. Among these galaninergic populations, only a subset of the RTN neurons (approximately 55%) exhibited activation in response to acute hypercapnia. Our findings highlight that the RTN is the likely source of galaninergic transmission to the VRC in response to hypercapnic stimuli.

Interdependence of cellular and network properties in respiratory rhythm generation

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

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.

Modeling Effects of Variable preBötzinger Complex Network Topology and Cellular Properties on Opioid-Induced Respiratory Depression and Recovery

The preBötzinger complex (preBötC), located in the medulla, is the essential rhythm-generating neural network for breathing. The actions of opioids on this network impair its ability to generate robust, rhythmic output, contributing to life-threatening opioid-induced respiratory depression (OIRD). The occurrence of OIRD varies across individuals and internal and external states, increasing the risk of opioid use, yet the mechanisms of this variability are largely unknown. In this study, we utilize a computational model of the preBötC to perform several in silico experiments exploring how differences in network topology and the intrinsic properties of preBötC neurons influence the sensitivity of the network rhythm to opioids. We find that rhythms produced by preBötC networks in silico exhibit variable responses to simulated opioids, similar to the preBötC network in vitro. This variability is primarily due to random differences in network topology and can be manipulated by imposed changes in network connectivity and intrinsic neuronal properties. Our results identify features of the preBötC network that may regulate its susceptibility to opioids.

Interdependence of cellular and network properties in respiratory rhythmogenesis

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

Multiple neuron clusters on Micro-Electrode Arrays as an in vitro model of brain network

Understanding the brain functioning is essential for governing brain processes with the aim of managing pathological network dysfunctions. Due to the morphological and biochemical complexity of the central nervous system, the development of general models with predictive power must start from in vitro brain network engineering. In the present work, we realized a micro-electrode array (MEA)-based in vitro brain network and studied its emerging dynamical properties. We obtained four-neuron-clusters (4N) assemblies by plating rat embryo cortical neurons on 60-electrode MEA with cross-shaped polymeric masks and compared the emerging dynamics with those of sister single networks (1N). Both 1N and 4N assemblies exhibited spontaneous electrical activity characterized by spiking and bursting signals up to global activation by means of network bursts. Data revealed distinct patterns of network activity with differences between 1 and 4N. Rhythmic network bursts and dominant initiator clusters suggested pacemaker activities in both assembly types, but the propagation of activation sequences was statistically influenced by the assembly topology. We proved that this rhythmic activity was ivabradine sensitive, suggesting the involvement of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, and propagated across the real clusters of 4N, or corresponding virtual clusters of 1N, with dominant initiator clusters, and nonrandom cluster activation sequences. The occurrence of nonrandom series of identical activation sequences in 4N revealed processes possibly ascribable to neuroplasticity. Hence, our multi-network dissociated cortical assemblies suggest the relevance of pacemaker neurons as essential elements for generating brain network electrophysiological patterns; indeed, such evidence should be considered in the development of computational models for envisaging network behavior both in physiological and pathological conditions.

The protective effect of inhaled corticosteroid on lung inflammation and breathing pattern complexity in a rat model of asthma

Asthma is a heterogeneous disease in which the complexity of the breathing pattern reduces as the severity of the disease increases. Since the pathophysiological basis of reduced breathing pattern complexity in asthma is unclear, in this study, we investigated the effect of reducing inflammation using an inhaled corticosteroid (fluticasone propionate) on the breathing pattern of a rat model of asthma. Detrended fluctuation analysis, sample entropy, and cross-sample entropy analysis of both inter-breath interval and respiratory volume time series showed that early treatment with inhaled corticosteroids not only diminishes lung inflammation and airway hyper-responsiveness, but also has a protective effect against the reduction of breathing pattern complexity due to asthma. However, late treatment had a partial effect on asthma-induced respiratory pattern changes. Since inflammation is a key factor in shifting breathing dynamics away from normal fluctuations, these findings further emphasize the importance of early treatment of asthma with corticosteroids.

Predictions and experimental tests of a new biophysical model of the mammalian respiratory oscillator

Previously our computational modeling studies (Phillips et al., 2019) proposed that neuronal persistent sodium current (INaP) and calcium-activated non-selective cation current (ICAN) are key biophysical factors that, respectively, generate inspiratory rhythm and burst pattern in the mammalian preBötzinger complex (preBötC) respiratory oscillator isolated in vitro. Here, we experimentally tested and confirmed three predictions of the model from new simulations concerning the roles of INaP and ICAN: (1) INaP and ICAN blockade have opposite effects on the relationship between network excitability and preBötC rhythmic activity; (2) INaP is essential for preBötC rhythmogenesis; and (3) ICAN is essential for generating the amplitude of rhythmic output but not rhythm generation. These predictions were confirmed via optogenetic manipulations of preBötC network excitability during graded INaP or ICAN blockade by pharmacological manipulations in slices in vitro containing the rhythmically active preBötC from the medulla oblongata of neonatal mice. Our results support and advance the hypothesis that INaP and ICAN mechanistically underlie rhythm and inspiratory burst pattern generation, respectively, in the isolated preBötC.

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.

Correlation Analysis of Synchronization Type and Degree in Respiratory Neural Network

Pre-Bötzinger complex (PBC) is a necessary condition for the generation of respiratory rhythm. Due to the existence of synaptic gaps, delay plays a key role in the synchronous operation of coupled neurons. In this study, the relationship between synchronization and correlation degree is established for the first time by using ISI bifurcation and correlation coefficient, and the relationship between synchronization and correlation degree is discussed under the conditions of no delay, symmetric delay, and asymmetric delay. The results show that the phase synchronization of two coupling PBCs is closely related to the weak correlation, that is, the weak phase synchronization may occur under the condition of incomplete synchronization. Moreover, the time delay and coupling strength are controlled in the modified PBC network model, which not only reveals the law of PBC firing transition but also reveals the complex synchronization behavior in the coupled chaotic neurons. Especially, when the two coupled neurons are nonidentical, the complete synchronization will disappear. These results fully reveal the dynamic behavior of the PBC neural system, which is helpful to explore the signal transmission and coding of PBC neurons and provide theoretical value for further understanding respiratory rhythm.

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.

Optimizing phase variability threshold for automated synchrogram analysis of cardiorespiratory interactions in amateur cyclists

We propose a procedure suitable for automated synchrogram analysis for setting the threshold below which phase variability between two marker event series is of such a negligible amount that the null hypothesis of phase desynchronization can be rejected. The procedure exploits the principle of maximizing the likelihood of detecting phase synchronization epochs and it is grounded on a surrogate data approach testing the null hypothesis of phase uncoupling. The approach was applied to assess cardiorespiratory phase interactions between heartbeat and inspiratory onset in amateur cyclists before and after 11-week inspiratory muscle training (IMT) at different intensities and compared to a more traditional approach to set phase variability threshold. The proposed procedure was able to detect the decrease in cardiorespiratory phase locking strength during vagal withdrawal induced by the modification of posture from supine to standing. IMT had very limited effects on cardiorespiratory phase synchronization strength and this result held regardless of the training intensity. In amateur athletes training, the inspiratory muscles did not limit the decrease in cardiorespiratory phase synchronization observed in the upright position as a likely consequence of the modest impact of this respiratory exercise, regardless of its intensity, on cardiac vagal control. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Apnoea of Prematurity and Neurodevelopmental Outcomes: Current Understanding and Future Prospects for Research

Infants who are born prematurely are at significant risk of apnoea. In addition to the short-term consequences such as hypoxia, apnoea of prematurity has been associated with long-term morbidity, including poor neurodevelopmental outcomes. Clinical trials have illustrated the importance of methylxanthine drugs, in particular caffeine, in reducing the risk of long term adverse neurodevelopmental outcomes. However, the extent to which apnoea is causative of this secondary neurodevelopmental delay or is just associated in a background of other sequelae of prematurity remains unclear. In this review, we first discuss the pathophysiology of apnoea of prematurity, previous studies investigating the relationship between apnoea and neurodevelopmental delay, and treatment of apnoea with caffeine therapy. We propose a need for better methods of measuring apnoea, along with improved understanding of the neonatal brain's response to consequent hypoxia. Only then can we start to disentangle the effects of apnoea on neurodevelopment in preterm infants. Moreover, by better identifying those infants who are at risk of apnoea, and neurodevelopmental delay, we can work toward a risk stratification system for these infants that is clinically actionable, for example, with doses of caffeine tailored to the individual. Optimising treatment of apnoea for individual infants will improve neonatal care and long-term outcomes for this population.

Volatile anesthetics maintain tidal volume and minute ventilation to a greater degree than propofol under spontaneous respiration

Background

Although general anesthetics depress spontaneous respiration, the comprehensive effect of general anesthetics on respiratory function remains unclear. We aimed to investigate the effects of general anesthetics on spontaneous respiration in non-intubated mice with different types and doses of general anesthetic.

Methods

Adult C57BL/6 J mice were administered intravenous anesthetics, including propofol and etomidate, and inhalational anesthetics, including sevoflurane and isoflurane in vivo at doses of 0.5-, 1.0-, and 2.0-times the minimum alveolar concentration (MAC)/median effective dose (ED₅₀) to induce loss of the righting reflex (LORR). Whole-body plethysmography (WBP) was applied to measure parameters of respiration under unrestricted conditions without endotracheal intubation. The alteration in respiratory sensitivity to carbon dioxide (CO₂) under general anesthesia was also determined. The following respiratory parameters were continuously recorded during anesthesia or CO₂ exposure: respiratory frequency (FR), tidal volume (TV), minute ventilation (MV), expiratory time (TE), inspiratory time (TI), and inspiratory–expiratory time ratio (I/E), and peak inspiratory flow.

Results

Sub-anesthetic concentrations (0.5 MAC) of sevoflurane or isoflurane increased FR, TV, and MV. With isoflurane and sevoflurane exposure, the CO₂-evoked increases in FR, TV, and MV were decreased. Compared with inhalational anesthetics, propofol and etomidate induced respiratory suppression, affecting FR, TV, and MV. In 100% oxygen (O₂), FR in the group that received propofol 1.0-times the ED₅₀ was 69.63 ± 33.44 breaths/min compared with 155.68 ± 64.42 breaths/min in the etomidate-treated group. In the same groups, FR was 88.72 ± 34.51 breaths/min and 225.10 ± 59.82 breaths/min, respectively, in 3% CO₂ and 144.17 ± 63.25 breaths/min and 197.70 ± 41.93 breaths/min, respectively, in 5% CO₂. A higher CO₂ sensitivity was found in etomidate-treated mice compared with propofol-treated mice. In addition, propofol induced a greater decrease in FR, MV, and I/E ratio compared with etomidate, sevoflurane, and isoflurane at equivalent doses (all P < 0.05).

Conclusions

General anesthetics differentially modulate spontaneous breathing in vivo. Volatile anesthetics increase FR, TV, and MV at sub-anesthetic concentrations, while they decrease FR at higher concentrations. Propofol consistently depressed respiratory parameters to a greater degree than etomidate.

The lamprey respiratory network: Some evolutionary aspects

Breathing is a complex behaviour that involves rhythm generating networks. In this review, we examine the main characteristics of respiratory rhythm generation in vertebrates and, in particular, we describe the main results of our studies on the role of neural mechanisms involved in the neuromodulation of the lamprey respiration. The lamprey respiratory rhythm generator is located in the paratrigeminal respiratory group (pTRG) and shows similarities with the mammalian preBötzinger complex. In fact, within the pTRG a major role is played by glutamate, but also GABA and glycine display important contributions. In addition, neuromodulatory influences are exerted by opioids, substance P, acetylcholine and serotonin. Both structures respond to exogenous ATP with a biphasic response and astrocytes there located strongly contribute to the modulation of the respiratory pattern. The results emphasize that some important characteristics of the respiratory rhythm generating network are, to a great extent, maintained throughout evolution.

Lack of association between heart period variability asymmetry and respiratory sinus arrhythmia in healthy and chronic heart failure individuals

Temporal asymmetry is a peculiar aspect of heart period (HP) variability (HPV). HPV asymmetry (HPVA) is reduced with aging and pathology, but its origin is not fully elucidated. Given the impact of respiration on HPV resulting in the respiratory sinus arrhythmia (RSA) and the asymmetric shape of the respiratory pattern, a possible link between HPVA and RSA might be expected. In this study we tested the hypothesis that HPVA is significantly associated with RSA and asymmetry of the respiratory rhythm. We studied 42 middle-aged healthy (H) subjects, and 56 chronic heart failure (CHF) patients of whom 26 assigned to the New York Heart Association (NYHA) class II (CHF-II) and 30 to NYHA class III (CHF-III). Electrocardiogram and lung volume were monitored for 8 minutes during spontaneous breathing (SB) and controlled breathing (CB) at 15 breaths/minute. The ratio of inspiratory (INSP) to expiratory (EXP) phases, namely the I/E ratio, and RSA were calculated. HPVA was estimated as the percentage of negative HP variations, traditionally measured via the Porta’s index (PI). Departures of PI from 50% indicated HPVA and its significance was tested via surrogate data. We found that RSA increased during CB and I/E ratio was smaller than 1 in all groups and experimental conditions. In H subjects the PI was about 50% during SB and it increased significantly during CB. In both CHF-II and CHF-III groups the PI was about 50% during SB and remained unmodified during CB. The PI was uncorrelated with RSA and I/E ratio regardless of the experimental condition and group. Pooling together data of different experimental conditions did not affect conclusions. Therefore, we conclude that the HPVA cannot be explained by RSA and/or I/E ratio, thus representing a peculiar feature of the cardiac control that can be aroused in middle-aged H individuals via CB.

Physiological Basis of Neonatal Aerodigestive Difficulties in Chronic Lung Disease

In the United States, preterm birth rates have steadily increased since 2014. Despite the recent advances in neonatal-perinatal care, more than 40% of very low-birth-weight infants develop chronic lung disease (CLD) and almost 25% have feeding difficulties resulting in delayed achievement of full oral feeds and longer hospital stay. Establishment of full oral feeds, a major challenge for preterm infants, becomes magnified among those on respiratory support and/or with CLD. The strategies to minimize aerodigestive disorders include supporting nonnutritive sucking, developing infant-directed feeding protocols, sensory oromotor stimulation, and early introduction of oral feeds.

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.

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.

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.

Pituitary adenylate cyclase-activating polypeptide: Postnatal development in multiple brain stem respiratory-related nuclei in the rat

The pituitary adenylate cyclase-activating polypeptide (PACAP) plays an important role in anterior pituitary hormone secretion, neurotransmission, and the control of breathing. Mice lacking PACAP die suddenly mainly in the 2^nd postnatal week, coinciding temporally with a critical period of respiratory development uncovered by our laboratory in the rat. The goal of the current study was to test our hypothesis that PACAP expression is reduced during the critical period in normal rats. We undertook immunohistochemistry and optical densitometry of PACAP (specifically PACAP38) in several brain stem respiratory-related nuclei of postnatal days P2-21 rats, and found that PACAP immunoreactivity was significantly reduced at P12 in the pre-Bötzinger complex, nucleus ambiguus, hypoglossal nucleus, and the ventrolateral subnucleus of the nucleus tractus solitarius. No changes were observed in the control, non-respiratory cuneate nucleus at P12. Results imply that the down-regulation of PACAP during normal postnatal development may contribute to the critical period of vulnerability, when the animals' response to hypoxia is at its weakest.

State-dependent rhythmogenesis and frequency control in a half-center locomotor CPG

The spinal locomotor central pattern generator (CPG) generates rhythmic activity with alternating flexion and extension phases. This rhythmic pattern is likely to result from inhibitory interactions between neural populations representing flexor and extensor half-centers. However, it is unclear whether the flexor-extensor CPG has a quasi-symmetric organization with both half-centers critically involved in rhythm generation, features an asymmetric organization with flexor-driven rhythmogenesis, or comprises a pair of intrinsically rhythmic half-centers. There are experimental data that support each of the above concepts but appear to be inconsistent with the others. In this theoretical/modeling study, we present and analyze a CPG model architecture that can operate in different regimes consistent with the above three concepts depending on conditions, which are defined by external excitatory drives to CPG half-centers. We show that control of frequency and phase durations within each regime depends on network dynamics, defined by the regime-dependent expression of the half-centers' intrinsic rhythmic capabilities and the operating phase transition mechanisms (escape vs. release). Our study suggests state dependency in locomotor CPG operation and proposes explanations for seemingly contradictory experimental data. NEW & NOTEWORTHY Our theoretical/modeling study focuses on the analysis of locomotor central pattern generators (CPGs) composed of conditionally bursting half-centers coupled with reciprocal inhibition and receiving independent external drives. We show that this CPG framework can operate in several regimes consistent with seemingly contradictory experimental data. In each regime, we study how intrinsic dynamics and phase-switching mechanisms control oscillation frequency and phase durations. Our results provide insights into the organization of spinal circuits controlling locomotion.

Effect of Rho-kinase inhibition on complexity of breathing pattern in a guinea pig model of asthma

Asthma represents an episodic and fluctuating behavior characterized with decreased complexity of respiratory dynamics. Several evidence indicate that asthma severity or control is associated with alteration in variability of lung function. The pathophysiological basis of alteration in complexity of breathing pattern in asthma has remained poorly understood. Regarding the point that Rho-kinase is involved in pathophysiology of asthma, in present study we investigated the effect of Rho-kinase inhibition on complexity of respiratory dynamics in a guinea pig model of asthma. Male Dunkin Hartley guinea pigs were exposed to 12 series of inhalations with ovalbumin or saline. Animals were treated by the Rho-kinase inhibitor Y-27632 (1mM aerosols) prior to each allergen challenge. We recorded respiration of conscious animals using whole-body plethysmography. Exposure to ovalbumin induced lung inflammation, airway hyperresponsiveness and remodeling including goblet cell hyperplasia, increase in the thickness of airways smooth muscles and subepithelial collagen deposition. Complexity analysis of respiratory dynamics revealed a dramatic decrease in irregularity of respiratory rhythm representing less complexity in asthmatic guinea pigs. Inhibition of Rho-kinase reduced the airway remodeling and hyperreponsiveness, but had no significant effect on lung inflammation and complexity of respiratory dynamics in asthmatic animals. It seems that airway hyperresponsiveness and remodeling do not significantly affect the complexity of respiratory dynamics. Our results suggest that inflammation might be the probable cause of shift in the respiratory dynamics away from the normal fluctuation in asthma.

Eupnea, tachypnea, and autoresuscitation in a closed-loop respiratory control model

How sensory information influences the dynamics of rhythm generation varies across systems, and general principles for understanding this aspect of motor control are lacking. Determining the origin of respiratory rhythm generation is challenging because the mechanisms in a central circuit considered in isolation may be different from those in the intact organism. We analyze a closed-loop respiratory control model incorporating a central pattern generator (CPG), the Butera-Rinzel-Smith (BRS) model, together with lung mechanics, oxygen handling, and chemosensory components. We show that 1) embedding the BRS model neuron in a control loop creates a bistable system; 2) although closed-loop and open-loop (isolated) CPG systems both support eupnea-like bursting activity, they do so via distinct mechanisms; 3) chemosensory feedback in the closed loop improves robustness to variable metabolic demand; 4) the BRS model conductances provide an autoresuscitation mechanism for recovery from transient interruption of chemosensory feedback; and 5) the in vitro brain stem CPG slice responds to hypoxia with transient bursting that is qualitatively similar to in silico autoresuscitation. Bistability of bursting and tonic spiking in the closed-loop system corresponds to coexistence of eupnea-like breathing, with normal minute ventilation and blood oxygen level and a tachypnea-like state, with pathologically reduced minute ventilation and critically low blood oxygen. Disruption of the normal breathing rhythm, through either imposition of hypoxia or interruption of chemosensory feedback, can push the system from the eupneic state into the tachypneic state. We use geometric singular perturbation theory to analyze the system dynamics at the boundary separating eupnea-like and tachypnea-like outcomes.NEW & NOTEWORTHY A common challenge facing rhythmic biological processes is the adaptive regulation of central pattern generator (CPG) activity in response to sensory feedback. We apply dynamical systems tools to understand several properties of a closed-loop respiratory control model, including the coexistence of normal and pathological breathing, robustness to changes in metabolic demand, spontaneous autoresuscitation in response to hypoxia, and the distinct mechanisms that underlie rhythmogenesis in the intact control circuit vs. the isolated, open-loop CPG.

Seizure-associated central apnea in a rat model: Evidence for resetting the respiratory rhythm and activation of the diving reflex

Respiratory derangements, including irregular, tachypnic breathing and central or obstructive apnea can be consequences of seizure activity in epilepsy patients and animal models. Periods of seizure-associated central apnea, defined as periods >1s with rapid onset and offset of no airflow during plethysmography, suggest that seizures spread to brainstem respiratory regions to disrupt breathing. We sought to characterize seizure-associated central apneic episodes as an indicator of seizure impact on the respiratory rhythm in rats anesthetized with urethane and given parenteral kainic acid to induce recurring seizures. We measured central apneic period onsets and offsets to determine if onset-offset relations were a consequence of 1) a reset of the respiratory rhythm, 2) a transient pausing of the respiratory rhythm, resuming from the pause point at the end of the apneic period, 3) a transient suppression of respiratory behavior with apnea offset predicted by a continuation of the breathing pattern preceding apnea, or 4) a random re-entry into the respiratory cycle. Animals were monitored with continuous ECG, EEG, and plethysmography. One hundred ninety central apnea episodes (1.04 to 36.18s, mean: 3.2±3.7s) were recorded during seizure activity from 7 rats with multiple apneic episodes. The majority of apneic period onsets occurred during expiration (125/161 apneic episodes, 78%). In either expiration or inspiration, apneic onsets tended to occur late in the cycle, i.e. between the time of the peak and end of expiration (82/125, 66%) or inspiration (34/36, 94%). Apneic period offsets were more uniformly distributed between early and late expiration (27%, 34%) and inspiration (16%, 23%). Differences between the respiratory phase at the onset of apnea and the corresponding offset phase varied widely, even within individual animals. Each central apneic episode was associated with a high frequency event in EEG or ECG records at onset. High frequency events that were not associated with flatline plethysmographs revealed a constant plethysmograph pattern within each animal, suggesting a clear reset of the respiratory rhythm. The respiratory rhythm became highly variable after about 1s, however, accounting for the unpredictability of the offset phase. The dissociation of respiratory rhythm reset from the cessation of airflow also suggested that central apneic periods involved activation of brainstem regions serving the diving reflex to eliminate the expression of respiratory movements. This conclusion was supported by the decreased heart rate as a function of apnea duration. We conclude that seizure-associated central apnea episodes are associated with 1) a reset of the respiratory rhythm, and 2) activation of brainstem regions serving the diving reflex to suppress respiratory behavior. The significance of these conclusions is that these details of seizure impact on brainstem circuitry represent metrics for assessing seizure spread and potentially subclassifying seizure patterns.

The Onset of the Fetal Respiratory Rhythm: An Emergent Property Triggered by Chemosensory Drive?

The mechanisms responsible for the onset of respiratory activity during fetal life are unknown. The onset of respiratory rhythm may be a consequence of the genetic program of each of the constituents of the respiratory network, so they start to interact and generate respiratory cycles when reaching a certain degree of maturation. Alternatively, generation of cycles might require the contribution of recently formed sensory inputs that will trigger oscillatory activity in the nascent respiratory neural network. If this hypothesis is true, then sensory input to the respiratory generator must be already formed and become functional before the onset of fetal respiration. In this review, we evaluate the timing of the onset of the respiratory rhythm in comparison to the appearance of receptors, neurotransmitter machinery, and afferent projections provided by two central chemoreceptive nuclei, the raphe and locus coeruleus nuclei.

Effects of ion channel noise on neural circuits: an application to the respiratory pattern generator to investigate breathing variability

Neural activity generally displays irregular firing patterns even in circuits with apparently regular outputs, such as motor pattern generators, in which the output frequency fluctuates randomly around a mean value. This "circuit noise" is inherited from the random firing of single neurons, which emerges from stochastic ion channel gating (channel noise), spontaneous neurotransmitter release, and its diffusion and binding to synaptic receptors. Here we demonstrate how to expand conductance-based network models that are originally deterministic to include realistic, physiological noise, focusing on stochastic ion channel gating. We illustrate this procedure with a well-established conductance-based model of the respiratory pattern generator, which allows us to investigate how channel noise affects neural dynamics at the circuit level and, in particular, to understand the relationship between the respiratory pattern and its breath-to-breath variability. We show that as the channel number increases, the duration of inspiration and expiration varies, and so does the coefficient of variation of the breath-to-breath interval, which attains a minimum when the mean duration of expiration slightly exceeds that of inspiration. For small channel numbers, the variability of the expiratory phase dominates over that of the inspiratory phase, and vice versa for large channel numbers. Among the four different cell types in the respiratory pattern generator, pacemaker cells exhibit the highest sensitivity to channel noise. The model shows that suppressing input from the pons leads to longer inspiratory phases, a reduction in breathing frequency, and larger breath-to-breath variability, whereas enhanced input from the raphe nucleus increases breathing frequency without changing its pattern.

NEW & NOTEWORTHY

A major source of noise in neuronal circuits is the "flickering" of ion currents passing through the neurons' membranes (channel noise), which cannot be suppressed experimentally. Computational simulations are therefore the best way to investigate the effects of this physiological noise by manipulating its level at will. We investigate the role of noise in the respiratory pattern generator and show that endogenous, breath-to-breath variability is tightly linked to the respiratory pattern.

Neurons as oscillators

Regularly spiking neurons can be described as oscillators. In this article we review some of the insights gained from this conceptualization and their relevance for systems neuroscience. First, we explain how a regularly spiking neuron can be viewed as an oscillator and how the phase-response curve (PRC) describes the response of the neuron's spike times to small perturbations. We then discuss the meaning of the PRC for a single neuron's spiking behavior and review the PRCs measured from a variety of neurons in a range of spiking regimes. Next, we show how the PRC can be related to a number of common measures used to quantify neuronal firing, such as the spike-triggered average and the peristimulus histogram. We further show that the response of a neuron to correlated inputs depends on the shape of the PRC. We then explain how the PRC of single neurons can be used to predict neural network behavior. Given the PRC, conduction delays, and the waveform and time course of the synaptic potentials, it is possible to predict neural population behavior such as synchronization. The PRC also allows us to quantify the robustness of the synchronization to heterogeneity and noise. We finally ask how to combine the measured PRCs and the predictions based on PRC to further the understanding of systems neuroscience. As an example, we discuss how the change of the PRC by the neuromodulator acetylcholine could lead to a destabilization of cortical network dynamics. Although all of these studies are grounded in mathematical abstractions that do not strictly hold in biology, they provide good estimates for the emergence of the brain's network activity from the properties of individual neurons. The study of neurons as oscillators can provide testable hypotheses and mechanistic explanations for systems neuroscience.

Bio-Inspired Controller on an FPGA Applied to Closed-Loop Diaphragmatic Stimulation

Cervical spinal cord injury can disrupt connections between the brain respiratory network and the respiratory muscles which can lead to partial or complete loss of ventilatory control and require ventilatory assistance. Unlike current open-loop technology, a closed-loop diaphragmatic pacing system could overcome the drawbacks of manual titration as well as respond to changing ventilation requirements. We present an original bio-inspired assistive technology for real-time ventilation assistance, implemented in a digital configurable Field Programmable Gate Array (FPGA). The bio-inspired controller, which is a spiking neural network (SNN) inspired by the medullary respiratory network, is as robust as a classic controller while having a flexible, low-power and low-cost hardware design. The system was simulated in MATLAB with FPGA-specific constraints and tested with a computational model of rat breathing; the model reproduced experimentally collected respiratory data in eupneic animals. The open-loop version of the bio-inspired controller was implemented on the FPGA. Electrical test bench characterizations confirmed the system functionality. Open and closed-loop paradigm simulations were simulated to test the FPGA system real-time behavior using the rat computational model. The closed-loop system monitors breathing and changes in respiratory demands to drive diaphragmatic stimulation. The simulated results inform future acute animal experiments and constitute the first step toward the development of a neuromorphic, adaptive, compact, low-power, implantable device. The bio-inspired hardware design optimizes the FPGA resource and time costs while harnessing the computational power of spike-based neuromorphic hardware. Its real-time feature makes it suitable for in vivo applications.

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.

Organization of left-right coordination of neuronal activity in the mammalian spinal cord: Insights from computational modelling

KEY POINTS

Coordination of neuronal activity between left and right sides of the mammalian spinal cord is provided by several sets of commissural interneurons (CINs) whose axons cross the midline. Genetically identified inhibitory V0D and excitatory V0V CINs and ipsilaterally projecting excitatory V2a interneurons were shown to secure left-right alternation at different locomotor speeds. We have developed computational models of neuronal circuits in the spinal cord that include left and right rhythm-generating centres interacting bilaterally via three parallel pathways mediated by V0D , V2a-V0V and V3 neuron populations. The models reproduce the experimentally observed speed-dependent left-right coordination in normal mice and the changes in coordination seen in mutants lacking specific neuron classes. The models propose an explanation for several experimental results and provide insights into the organization of the spinal locomotor network and parallel CIN pathways involved in gait control at different locomotor speeds.

ABSTRACT

Different locomotor gaits in mammals, such as walking or galloping, are produced by coordinated activity in neuronal circuits in the spinal cord. Coordination of neuronal activity between left and right sides of the cord is provided by commissural interneurons (CINs), whose axons cross the midline. In this study, we construct and analyse two computational models of spinal locomotor circuits consisting of left and right rhythm generators interacting bilaterally via several neuronal pathways mediated by different CINs. The CIN populations incorporated in the models include the genetically identified inhibitory (V0D ) and excitatory (V0V ) subtypes of V0 CINs and excitatory V3 CINs. The model also includes the ipsilaterally projecting excitatory V2a interneurons mediating excitatory drive to the V0V CINs. The proposed network architectures and CIN connectivity allow the models to closely reproduce and suggest mechanistic explanations for several experimental observations. These phenomena include: different speed-dependent contributions of V0D and V0V CINs and V2a interneurons to left-right alternation of neural activity, switching gaits between the left-right alternating walking-like activity and the left-right synchronous hopping-like pattern in mutants lacking specific neuron classes, and speed-dependent asymmetric changes of flexor and extensor phase durations. The models provide insights into the architecture of spinal network and the organization of parallel inhibitory and excitatory CIN pathways and suggest explanations for how these pathways maintain alternating and synchronous gaits at different locomotor speeds. The models propose testable predictions about the neural organization and operation of mammalian locomotor circuits.

Modeling the effects of extracellular potassium on bursting properties in pre-Bötzinger complex neurons

There are many types of neurons that intrinsically generate rhythmic bursting activity, even when isolated, and these neurons underlie several specific motor behaviors. Rhythmic neurons that drive the inspiratory phase of respiration are located in the medullary pre-Bötzinger Complex (pre-BötC). However, it is not known if their rhythmic bursting is the result of intrinsic mechanisms or synaptic interactions. In many cases, for bursting to occur, the excitability of these neurons needs to be elevated. This excitation is provided in vitro (e.g. in slices), by increasing extracellular potassium concentration (K out) well beyond physiologic levels. Elevated K out shifts the reversal potentials for all potassium currents including the potassium component of leakage to higher values. However, how an increase in K out , and the resultant changes in potassium currents, induce bursting activity, have yet to be established. Moreover, it is not known if the endogenous bursting induced in vitro is representative of neural behavior in vivo. Our modeling study examines the interplay between K out, excitability, and selected currents, as they relate to endogenous rhythmic bursting. Starting with a Hodgkin-Huxley formalization of a pre-BötC neuron, a potassium ion component was incorporated into the leakage current, and model behaviors were investigated at varying concentrations of K out. Our simulations show that endogenous bursting activity, evoked in vitro by elevation of K out , is the result of a specific relationship between the leakage and voltage-dependent, delayed rectifier potassium currents, which may not be observed at physiological levels of extracellular potassium.

Somatostatin 2a receptors are not expressed on functionally identified respiratory neurons in the ventral respiratory column of the rat

Microinjection of somatostatin (SST) causes site-specific effects on respiratory phase transition, frequency, and amplitude when microinjected into the ventrolateral medulla (VLM) of the anesthetized rat, suggesting selective expression of SST receptors on different functional classes of respiratory neurons. Of the six subtypes of SST receptor, somatostatin 2a (sst2a ) is the most prevalent in the VLM, and other investigators have suggested that glutamatergic neurons in the preBötzinger Complex (preBötC) that coexpress neurokinin-1 receptor (NK1R), SST, and sst2a are critical for the generation of respiratory rhythm. However, quantitative data describing the distribution of sst2a in respiratory compartments other than preBötC, or on functionally identified respiratory neurons, is absent. Here we examine the medullary expression of sst2a with particular reference to glycinergic/expiratory neurons in the Bötzinger Complex (BötC) and NK1R-immunoreactive/inspiratory neurons in the preBötC. We found robust sst2a expression at all rostrocaudal levels of the VLM, including a large proportion of catecholaminergic neurons, but no colocalization of sst2a and glycine transporter 2 mRNA in the BötC. In the preBötC 54% of sst2a -immunoreactive neurons were also positive for NK1R. sst2a was not observed in any of 52 dye-labeled respiratory interneurons, including seven BötC expiratory-decrementing and 11 preBötC preinspiratory neurons. We conclude that sst2a is not expressed on BötC respiratory neurons and that phasic respiratory activity is a poor predictor of sst2a expression in the preBötC. Therefore, sst2a is unlikely to underlie responses to BötC SST injection, and is sparse or absent on respiratory neurons identified by classical functional criteria. J. Comp. Neurol. 524:1384-1398, 2016. © 2015 Wiley Periodicals, Inc.

Organization of the Mammalian Locomotor CPG: Review of Computational Model and Circuit Architectures Based on Genetically Identified Spinal Interneurons(1,2,3)

The organization of neural circuits that form the locomotor central pattern generator (CPG) and provide flexor–extensor and left–right coordination of neuronal activity remains largely unknown. However, significant progress has been made in the molecular/genetic identification of several types of spinal interneurons, including V0 (V0_D and V0_V subtypes), V1, V2a, V2b, V3, and Shox2, among others. The possible functional roles of these interneurons can be suggested from changes in the locomotor pattern generated in mutant mice lacking particular neuron types. Computational modeling of spinal circuits may complement these studies by bringing together data from different experimental studies and proposing the possible connectivity of these interneurons that may define rhythm generation, flexor–extensor interactions on each side of the cord, and commissural interactions between left and right circuits. This review focuses on the analysis of potential architectures of spinal circuits that can reproduce recent results and suggest common explanations for a series of experimental data on genetically identified spinal interneurons, including the consequences of their genetic ablation, and provides important insights into the organization of the spinal CPG and neural control of locomotion.

Differences in respiratory changes and Fos expression in the ventrolateral medulla of rats exposed to hypoxia, hypercapnia, and hypercapnic hypoxia

Respiratory responses to hypoxia and/or hypercapnia, and their relationship to neural activity in the ventrolateral medulla (VLM), which includes the respiratory center, have not yet been elucidated in detail. We herein examined respiratory responses during exposure of 10% O2 (hypoxia), 10% CO2 (hypercapnia), and 10% O2-10% CO2 (hypercapnic hypoxia) using plethysmography. In addition to recording respiration, Fos expressions were examined in the VLM of the rat exposed to each gas to analyze neural activity. Respiratory frequency was increased in rats exposed to hypoxia, and Fos-positive neurons were observed in the caudal VLM (cVLM) and medial VLM (mVLM). Tidal volume was increased in rats exposed to hypercapnia, and Fos-positive neurons were observed in the rostral VLM (rVLM) includes the retrotrapezoid nucleus (RTN) and mVLM. Tidal volume was enhanced in rats exposed to hypercapnic hypoxia, similar to that in hypercapnia-exposed rats, and Fos-positive neurons were observed in the entire region of the VLM. In the mVLM and cVLM, double immunofluorescence showed Fos-immunoreactive nerve cells were also immunoreactive to dopamine β-hydroxylase, the marker for A1/C1 catecholaminergic neuron. These results suggested that hypoxia and hypercapnia modulated rhythmogenic microcircuits in the mVLM via A1/C1 neurons and the RTN, respectively.

Respiratory neuron characterization reveals intrinsic bursting properties in isolated adult turtle brainstems (Trachemys scripta)

It is not known whether respiratory neurons with intrinsic bursting properties exist within ectothermic vertebrate respiratory control systems. Thus, isolated adult turtle brainstems spontaneously producing respiratory motor output were used to identify and classify respiratory neurons based on their firing pattern relative to hypoglossal (XII) nerve activity. Most respiratory neurons (183/212) had peak activity during the expiratory phase, while inspiratory, post-inspiratory, and novel pre-expiratory neurons were less common. During synaptic blockade conditions, ∼10% of respiratory neurons fired bursts of action potentials, with post-inspiratory cells (6/9) having the highest percentage of intrinsic burst properties. Most intrinsically bursting respiratory neurons were clustered at the level of the vagus (X) nerve root. Synaptic inhibition blockade caused seizure-like activity throughout the turtle brainstem, which shows that the turtle respiratory control system is not transformed into a network driven by intrinsically bursting respiratory neurons. We hypothesize that intrinsically bursting respiratory neurons are evolutionarily conserved and represent a potential rhythmogenic mechanism contributing to respiration in adult turtles.

Respiratory rhythm generation in vivo

The cellular and circuit mechanisms generating the rhythm of breathing in mammals have been under intense investigation for decades. Here, we try to integrate the key discoveries into an updated description of the basic neural processes generating respiratory rhythm under in vivo conditions.

GABAergic and glycinergic inputs modulate rhythmogenic mechanisms in the lamprey respiratory network

We have previously shown that GABA and glycine modulate respiratory activity in the in vitro brainstem preparations of the lamprey and that blockade of GABAA and glycine receptors restores the respiratory rhythm during apnoea caused by blockade of ionotropic glutamate receptors. However, the neural substrates involved in these effects are unknown. To address this issue, the role of GABAA, GABAB and glycine receptors within the paratrigeminal respiratory group (pTRG), the proposed respiratory central pattern generator, and the vagal motoneuron region was investigated both during apnoea induced by blockade of glutamatergic transmission and under basal conditions through microinjections of specific antagonists. The removal of GABAergic, but not glycinergic transmission within the pTRG, causes the resumption of rhythmic respiratory activity during apnoea, and reveals the presence of a modulatory control of the pTRG under basal conditions. A blockade of GABAA and glycine receptors within the vagal region strongly increases the respiratory frequency through disinhibition of neurons projecting to the pTRG from the vagal region. These neurons were retrogradely labelled (neurobiotin) from the pTRG. Intense GABA immunoreactivity is observed both within the pTRG and the vagal area, which corroborates present findings. The results confirm the pTRG as a primary site of respiratory rhythm generation, and suggest that inhibition modulates the activity of rhythm-generating neurons, without any direct role in burst formation and termination mechanisms.

Computational models and emergent properties of respiratory neural networks

Computational models of the neural control system for breathing in mammals provide a theoretical and computational framework bringing together experimental data obtained from different animal preparations under various experimental conditions. Many of these models were developed in parallel and iteratively with experimental studies and provided predictions guiding new experiments. This data-driven modeling approach has advanced our understanding of respiratory network architecture and neural mechanisms underlying generation of the respiratory rhythm and pattern, including their functional reorganization under different physiological conditions. Models reviewed here vary in neurobiological details and computational complexity and span multiple spatiotemporal scales of respiratory control mechanisms. Recent models describe interacting populations of respiratory neurons spatially distributed within the Bötzinger and pre-Bötzinger complexes and rostral ventrolateral medulla that contain core circuits of the respiratory central pattern generator (CPG). Network interactions within these circuits along with intrinsic rhythmogenic properties of neurons form a hierarchy of multiple rhythm generation mechanisms. The functional expression of these mechanisms is controlled by input drives from other brainstem components,including the retrotrapezoid nucleus and pons, which regulate the dynamic behavior of the core circuitry. The emerging view is that the brainstem respiratory network has rhythmogenic capabilities at multiple levels of circuit organization. This allows flexible, state-dependent expression of different neural pattern-generation mechanisms under various physiological conditions,enabling a wide repertoire of respiratory behaviors. Some models consider control of the respiratory CPG by pulmonary feedback and network reconfiguration during defensive behaviors such as cough. Future directions in modeling of the respiratory CPG are considered.

Airflow pattern complexity during resting breathing in patients with COPD: effect of airway obstruction

We investigated the influence of airway obstruction in the complexity of the airflow pattern in COPD and its use as a marker of disease activity. The sample entropy (SampEnV') and the variability (SDV') of the airflow pattern were measured in a group of 88 subjects with various levels of airway obstruction. Airway obstruction resulted in a reduction in the SampEnV' (p<0.0001) that was significantly correlated with spirometric indices of airway obstruction (R=0.50, p<0.001). The early adverse effects in mild airway obstruction were detected by the SampEnV' with an accuracy of 84%. SDV' increased with airway obstruction (p<0.002). We conclude that (1) the airflow patterns in COPD exhibit reduced complexity compared with healthy subjects; (2) this reduction in complexity is proportional to airway obstruction; and (3) the evaluation of SampEnV' may provide novel respiratory biomarkers suitable to facilitate the diagnosis of respiratory abnormalities in COPD.

Postnatal development of glycine receptor subunits α1, α2, α3, and β immunoreactivity in multiple brain stem respiratory-related nuclear groups of the rat

The respiratory system is immature at birth and significant development occurs postnatally. A critical period of respiratory development occurs in rats around postnatal days 12-13, when enhanced inhibition dominates over suppressed excitation. The mechanisms underlying the heightened inhibition are not fully understood. The present study tested our hypothesis that the inhibition is marked by a switch in glycine receptor subunits from neonatal to adult form around the critical period. An in-depth immunohistochemical and single neuron optical densitometric study was undertaken on four respiratory-related nuclear groups (the pre-Bötzinger complex, nucleus ambiguus, hypoglossal nucleus, and ventrolateral subnucleus of solitary tract nucleus), and a non-respiratory cuneate nucleus in P2-21 rats. Our data revealed that in the respiratory-related nuclear groups: (1) the expressions of GlyRα2 and GlyRα3 were relatively high at P2, but declined after 1-1½ weeks to their lowest levels at P21; (2) the expression of GlyRα1 increased with age and reached significance at P12; and (3) the expression of GlyRβ rose from P2 to P12 followed by a slight decline until P21. No distinct increase in GlyRα1 at P12 was noted in the cuneate nucleus. Thus, there is a switch in dominance of expression from neonatal GlyRα2/α3 to the adult GlyRα1 and a heightened expression of GlyRα1 around the critical period in all respiratory-related nuclear groups, thereby supporting enhanced inhibition at that time. The rise in the expression of GlyRβ around P12 indicates that it plays an important role in forming the mature heteropentameric glycine receptors in these brain stem nuclear groups.