Multimodal medullary neurons and correlational linkages of the respiratory network

Deglutition and the Regulation of the Swallow Motor Pattern

Despite centuries of investigation, questions and controversies remain regarding the fundamental genesis and motor pattern of swallow. Two significant topics include inspiratory muscle activity during swallow (Schluckatmung, i.e., "swallow-breath") and anatomical boundaries of the swallow pattern generator. We discuss the long history of reports regarding the presence or absence of Schluckatmung and the possible advantages of and neural basis for such activity, leading to current theories and novel experimental directions.

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.

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.

Carotid chemoreceptors tune breathing via multipath routing: reticular chain and loop operations supported by parallel spike train correlations

We tested the hypothesis that carotid chemoreceptors tune breathing through parallel circuit paths that target distinct elements of an inspiratory neuron chain in the ventral respiratory column (VRC). Microelectrode arrays were used to monitor neuronal spike trains simultaneously in the VRC, peri-nucleus tractus solitarius (p-NTS)-medial medulla, the dorsal parafacial region of the lateral tegmental field (FTL-pF), and medullary raphe nuclei together with phrenic nerve activity during selective stimulation of carotid chemoreceptors or transient hypoxia in 19 decerebrate, neuromuscularly blocked, and artificially ventilated cats. Of 994 neurons tested, 56% had a significant change in firing rate. A total of 33,422 cell pairs were evaluated for signs of functional interaction; 63% of chemoresponsive neurons were elements of at least one pair with correlational signatures indicative of paucisynaptic relationships. We detected evidence for postinspiratory neuron inhibition of rostral VRC I-Driver (pre-Bötzinger) neurons, an interaction predicted to modulate breathing frequency, and for reciprocal excitation between chemoresponsive p-NTS neurons and more downstream VRC inspiratory neurons for control of breathing depth. Chemoresponsive pericolumnar tonic expiratory neurons, proposed to amplify inspiratory drive by disinhibition, were correlationally linked to afferent and efferent "chains" of chemoresponsive neurons extending to all monitored regions. The chains included coordinated clusters of chemoresponsive FTL-pF neurons with functional links to widespread medullary sites involved in the control of breathing. The results support long-standing concepts on brain stem network architecture and a circuit model for peripheral chemoreceptor modulation of breathing with multiple circuit loops and chains tuned by tegmental field neurons with quasi-periodic discharge patterns. NEW & NOTEWORTHY We tested the long-standing hypothesis that carotid chemoreceptors tune the frequency and depth of breathing through parallel circuit operations targeting the ventral respiratory column. Responses to stimulation of the chemoreceptors and identified functional connectivity support differential tuning of inspiratory neuron burst duration and firing rate and a model of brain stem network architecture incorporating tonic expiratory "hub" neurons regulated by convergent neuronal chains and loops through rostral lateral tegmental field neurons with quasi-periodic discharge patterns.

Microinjection of kynurenic acid in the rostral nucleus of the tractus solitarius disrupts spatiotemporal aspects of mechanically induced tracheobronchial cough

The importance of neurons in the nucleus of the solitary tract (NTS) in the production of coughing was tested by microinjections of the nonspecific glutamate receptor antagonist kynurenic acid (kyn; 100 mM in artificial cerebrospinal fluid) in 15 adult spontaneously breathing anesthetized cats. Repetitive coughing was elicited by mechanical stimulation of the intrathoracic airway. Electromyograms (EMG) were recorded from inspiratory parasternal and expiratory transversus abdominis (ABD) muscles. Bilateral microinjections of kyn into the NTS rostral to obex [55 ± 4 nl total in 2 locations (n = 6) or 110 ± 4 nl total in 4 locations (n = 5)], primarily the ventrolateral subnucleus, reduced cough number and expiratory cough efforts (amplitudes of ABD EMG and maxima of esophageal pressure) compared with control. These microinjections also markedly prolonged the inspiratory phase, all cough-related EMG activation, and the total cough cycle duration as well as some other cough-related time intervals. In response to microinjections of kyn into the NTS rostral to the obex respiratory rate decreased, and there were increases in the durations of the inspiratory and postinspiratory phases and mean blood pressure. However, bilateral microinjections of kyn into the NTS caudal to obex as well as control vehicle microinjections in the NTS location rostral to obex had no effect on coughing or cardiorespiratory variables. These results are consistent with the existence of a critical component of the cough rhythmogenic circuit located in the rostral ventral and lateral NTS. Neuronal structures of the rostral NTS are significantly involved specifically in the regulation of cough magnitude and phase timing.NEW & NOTEWORTHY The nucleus of the solitary tract contains significant neuronal structures responsible for control of 1) cough excitability, 2) motor drive during cough, 3) cough phase timing, and 4) cough rhythmicity. Significant elimination of neurons in the solitary tract nucleus results in cough apraxia (incomplete and/or disordered cough pattern). The mechanism of the cough impairment is different from that for the concomitant changes in breathing.

Spinal metaplasticity in respiratory motor control

A hallmark feature of the neural system controlling breathing is its ability to exhibit plasticity. Less appreciated is the ability to exhibit metaplasticity, a change in the capacity to express plasticity (i.e., “plastic plasticity”). Recent advances in our understanding of cellular mechanisms giving rise to respiratory motor plasticity lay the groundwork for (ongoing) investigations of metaplasticity. This detailed understanding of respiratory metaplasticity will be essential as we harness metaplasticity to restore breathing capacity in clinical disorders that compromise breathing, such as cervical spinal injury, motor neuron disease and other neuromuscular diseases. In this brief review, we discuss key examples of metaplasticity in respiratory motor control, and our current understanding of mechanisms giving rise to spinal plasticity and metaplasticity in phrenic motor output; particularly after pre-conditioning with intermittent hypoxia. Progress in this area has led to the realization that similar mechanisms are operative in other spinal motor networks, including those governing limb movement. Further, these mechanisms can be harnessed to restore respiratory and non-respiratory motor function after spinal injury.

Physiological and pathophysiological interactions between the respiratory central pattern generator and the sympathetic nervous system

Respiratory modulation seen in the sympathetic nerve activity (SNA) implies that the respiratory and sympathetic networks interact. During hypertension elicited by chronic intermittent hypoxia (CIH), the SNA displays an enhanced respiratory modulation reflecting strengthened interactions between the networks. In this chapter, we review a series of experimental and modeling studies that help elucidate possible mechanisms of sympatho-respiratory coupling. We conclude that this coupling significantly contributes to both the sympathetic baroreflex and the augmented sympathetic activity after exposure to CIH. This conclusion is based on the following findings. (1) Baroreceptor activation results in perturbation of the respiratory pattern via transient activation of postinspiratory neurons in the Bötzinger complex (BötC). The same BötC neurons are involved in the respiratory modulation of SNA, and hence provide an additional pathway for the sympathetic baroreflex. (2) Under hypercapnia, phasic activation of abdominal motor nerves (AbN) is accompanied by synchronous discharges in SNA due to the common source of this rhythmic activity in the retrotrapezoid nucleus (RTN). CIH conditioning increases the CO2 sensitivity of central chemoreceptors in the RTN which results in the emergence of AbN and SNA discharges under normocapnic conditions similar to those observed during hypercapnia in naïve animals. Thus, respiratory-sympathetic interactions play an important role in defining sympathetic output and significantly contribute to the sympathetic activity and hypertension under certain physiological or pathophysiological conditions, and the theoretical framework presented may be instrumental in understanding of malfunctioning control of sympathetic activity in a variety of disease states.

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.

Blood pressure changes alter tracheobronchial cough: computational model of the respiratory-cough network and in vivo experiments in anesthetized cats

We tested the hypothesis, motivated in part by a coordinated computational cough network model, that alterations of mean systemic arterial blood pressure (BP) influence the excitability and motor pattern of cough. Model simulations predicted suppression of coughing by stimulation of arterial baroreceptors. In vivo experiments were conducted on anesthetized spontaneously breathing cats. Cough was elicited by mechanical stimulation of the intrathoracic airways. Electromyograms (EMG) of inspiratory parasternal, expiratory abdominal, laryngeal posterior cricoarytenoid (PCA), and thyroarytenoid muscles along with esophageal pressure (EP) and BP were recorded. Transiently elevated BP significantly reduced cough number, cough-related inspiratory, and expiratory amplitudes of EP, peak parasternal and abdominal EMG, and maximum of PCA EMG during the expulsive phase of cough, and prolonged the cough inspiratory and expiratory phases as well as cough cycle duration compared with control coughs. Latencies from the beginning of stimulation to the onset of cough-related diaphragm and abdominal activities were increased. Increases in BP also elicited bradycardia and isocapnic bradypnea. Reductions in BP increased cough number; elevated inspiratory EP amplitude and parasternal, abdominal, and inspiratory PCA EMG amplitudes; decreased total cough cycle duration; shortened the durations of the cough expiratory phase and cough-related abdominal discharge; and shortened cough latency compared with control coughs. Reduced BP also produced tachycardia, tachypnea, and hypocapnic hyperventilation. These effects of BP on coughing likely originate from interactions between barosensitive and respiratory brainstem neuronal networks, particularly by modulation of respiratory neurons within multiple respiration/cough-related brainstem areas by baroreceptor input.

Effect of baroreceptor stimulation on the respiratory pattern: insights into respiratory-sympathetic interactions

Sympathetic nerve activity (SNA) is modulated by respiratory activity which indicates the existence of direct interactions between the respiratory and sympathetic networks within the brainstem. Our experimental studies reveal that T(E) prolongation evoked by baroreceptor stimulation varies with respiratory phase and depends on the pons. We speculate that the sympathetic baroreceptor reflex, providing negative feedback from baroreceptors to the rostral ventrolateral medulla and SNA, has two pathways: one direct and independent of the respiratory-sympathetic interactions and the other operating via the respiratory pattern generator and is hence dependent on the respiratory modulation of SNA. Our experimental studies in the perfused in situ rat preparation and complementary computational modelling studies support the hypothesis that baroreceptor activation during expiration prolongs the T(E) via transient activation of post-inspiratory and inhibition of augmenting expiratory neurones of the Bötzinger Complex (BötC). We propose that these BötC neurones are also involved in the respiratory modulation of SNA, and contribute to the respiratory modulation of the sympathetic baroreceptor reflex.

Central and peripheral chemoreceptors evoke distinct responses in simultaneously recorded neurons of the raphé-pontomedullary respiratory network

The brainstem network for generating and modulating the respiratory motor pattern includes neurons of the medullary ventrolateral respiratory column (VRC), dorsolateral pons (PRG) and raphé nuclei. Midline raphé neurons are proposed to be elements of a distributed brainstem system of central chemoreceptors, as well as modulators of central chemoreceptors at other sites, including the retrotrapezoid nucleus. Stimulation of the raphé system or peripheral chemoreceptors can induce a long-term facilitation of phrenic nerve activity; central chemoreceptor stimulation does not. The network mechanisms through which each class of chemoreceptor differentially influences breathing are poorly understood. Microelectrode arrays were used to monitor sets of spike trains from 114 PRG, 198 VRC and 166 midline neurons in six decerebrate vagotomized cats; 356 were recorded during sequential stimulation of both receptor classes via brief CO(2)-saturated saline injections in vertebral (central) and carotid arteries (peripheral). Seventy neurons responded to both stimuli. More neurons were responsive only to peripheral challenges than those responsive only to central chemoreceptor stimulation (PRG, 20 : 4; VRC, 41 : 10; midline, 25 : 13). Of 16 474 pairs of neurons evaluated for short-time scale correlations, similar percentages of reference neurons in each brain region had correlation features indicative of a specific interaction with at least one target neuron: PRG (59.6%), VRC (51.0%) and raphé nuclei (45.8%). The results suggest a brainstem network architecture with connectivity that shapes the respiratory motor pattern via overlapping circuits that modulate central and peripheral chemoreceptor-mediated influences on breathing.

Pontine-ventral respiratory column interactions through raphe circuits detected using multi-array spike train recordings

Recently, Segers et al. identified functional connectivity between the ventrolateral respiratory column (VRC) and the pontine respiratory group (PRG). The apparent sparseness of detected paucisynaptic interactions motivated consideration of other potential functional pathways between these two regions. We report here evidence for "indirect" serial functional linkages between the PRG and VRC via intermediary brain stem midline raphé neurons. Arrays of microelectrodes were used to record sets of spike trains from a total of 145 PRG, 282 VRC, and 340 midline neurons in 11 decerebrate, vagotomized, neuromuscularly blocked, ventilated cats. Spike trains of 13,843 pairs of neurons that included at least one raphé cell were screened for respiratory modulation and short-time scale correlations. Significant correlogram features were detected in 7.2% of raphé-raphé (291/4,021), 4.3% of VRC-raphé (292/6,755), and 4.0% of the PRG-raphé (124/3,067) neuron pairs. Central peaks indicative of shared influences were the most common feature in correlations between pairs of raphé neurons, whereas correlated raphé-PRG and raphé-VRC neuron pairs displayed predominantly offset peaks and troughs, features suggesting a paucisynaptic influence of one neuron on the other. Overall, offset correlogram features provided evidence for 33 VRC-to-raphé-to-PRG and 45 PRG-to-raphé-to-VRC correlational linkage chains with one or two intermediate raphé neurons. The results support a respiratory network architecture with parallel VRC-to-PRG and PRG-to-VRC links operating through intervening midline circuits, and suggest that raphé neurons contribute to the respiratory modulation of PRG neurons and shape the respiratory motor pattern through coordinated divergent actions on both the PRG and VRC.

Reconfiguration of the pontomedullary respiratory network: a computational modeling study with coordinated in vivo experiments

A large body of data suggests that the pontine respiratory group (PRG) is involved in respiratory phase-switching and the reconfiguration of the brain stem respiratory network. However, connectivity between the PRG and ventral respiratory column (VRC) in computational models has been largely ad hoc. We developed a network model with PRG-VRC connectivity inferred from coordinated in vivo experiments. Neurons were modeled in the "integrate-and-fire" style; some neurons had pacemaker properties derived from the model of Breen et al. We recapitulated earlier modeling results, including reproduction of activity profiles of different respiratory neurons and motor outputs, and their changes under different conditions (vagotomy, pontine lesions, etc.). The model also reproduced characteristic changes in neuronal and motor patterns observed in vivo during fictive cough and during hypoxia in non-rapid eye movement sleep. Our simulations suggested possible mechanisms for respiratory pattern reorganization during these behaviors. The model predicted that network- and pacemaker-generated rhythms could be co-expressed during the transition from gasping to eupnea, producing a combined "burst-ramp" pattern of phrenic discharges. To test this prediction, phrenic activity and multiple single neuron spike trains were monitored in vagotomized, decerebrate, immobilized, thoracotomized, and artificially ventilated cats during hypoxia and recovery. In most experiments, phrenic discharge patterns during recovery from hypoxia were similar to those predicted by the model. We conclude that under certain conditions, e.g., during recovery from severe brain hypoxia, components of a distributed network activity present during eupnea can be co-expressed with gasp patterns generated by a distinct, functionally "simplified" mechanism.

Functional connectivity in the pontomedullary respiratory network

Current models propose that a neuronal network in the ventrolateral medulla generates the basic respiratory rhythm and that this ventrolateral respiratory column (VRC) is profoundly influenced by the neurons of the pontine respiratory group (PRG). However, functional connectivity among PRG and VRC neurons is poorly understood. This study addressed four model-based hypotheses: 1) the respiratory modulation of PRG neuron populations reflects paucisynaptic actions of multiple VRC populations; 2) functional connections among PRG neurons shape and coordinate their respiratory-modulated activities; 3) the PRG acts on multiple VRC populations, contributing to phase-switching; and 4) neurons with no respiratory modulation located in close proximity to the VRC and PRG have widely distributed actions on respiratory-modulated cells. Two arrays of microelectrodes with individual depth adjustment were used to record sets of spike trains from a total of 145 PRG and 282 VRC neurons in 10 decerebrate, vagotomized, neuromuscularly blocked, ventilated cats. Data were evaluated for respiratory modulation with respect to efferent phrenic motoneuron activity and short-timescale correlations indicative of paucisynaptic functional connectivity using cross-correlation analysis and the "gravity" method. Correlogram features were found for 109 (3%) of the 3,218 pairs composed of a PRG and a VRC neuron, 126 (12%) of the 1,043 PRG-PRG pairs, and 319 (7%) of the 4,340 VRC-VRC neuron pairs evaluated. Correlation linkage maps generated for the data support our four motivating hypotheses and suggest network mechanisms for proposed modulatory functions of the PRG.

Hypercapnia modulates synaptic interaction of cultured brainstem neurons

CO(2) is an important metabolic product whose concentrations are constantly monitored by CO(2) chemoreceptors. However, the high systemic CO(2) sensitivity may not be achieved by the CO(2) chemoreceptors without neuronal network processes. To show modulation of network properties during hypercapnia, we studied brainstem neurons dissociated from embryonic rats (P17-19) in multielectrode arrays (MEA) after initial period (3 weeks) of culture. Spike trains of 33,622 pairs of units were analyzed using peri-event histograms (PEH). The amplitude of peri-central peaks between two CO(2)-stimulated units increased and the peak latency decreased during hypercapnia. Similar enhancement of synaptic strength was observed in those sharing a common input. These phenomena were not seen in CO(2)-unresponsive neurons. The amplitude of peri-central peaks between two CO(2) inhibited units also increased without changing latency. Over 60% CO(2)-stimulated neurons studied received mono-/oligosynaptic inputs from other CO(2)-stimulated cells, whereas only approximately 10% CO(2)-unresponsive neurons had such synaptic inputs. A small number of brainstem neurons showed electrical couplings. The coupling efficiency of CO(2)-stimulated but not CO(2)-unresponsive units was suppressed by approximately 50% with high PCO(2). Inhibitory synaptic projections were also found, which was barely affected by hypercapnia. Consistent with the strengthening of excitatory synaptic connections, CO(2) sensitivity of post-synaptic neurons was significantly higher than presynaptic neurons. The difference was eliminated with blockade of presynaptic input. Based on these indirect assessments of synaptic interaction, our PEH analysis suggests that hypercapnia appears to modulate excitatory synaptic transmissions, especially those between CO(2)-stimulated neurons.

Mechanism of cardioventilatory coupling: insights from cardiac pacing, vagotomy, and sinoaortic denervation in the anesthetized rat

Cardioventilatory coupling (CVC), a temporal alignment between the heartbeat and inspiratory activity, is a major determinant of breath-to-breath variation in observed respiratory rate (f(o)). The cardiac-trigger hypothesis attributes this to adjustments of respiratory timing by baroreceptor afferent impulses to the central respiratory pattern generator. A mathematical model of this hypothesis indicates that apparent CVC in graphical plots of ECG R wave vs. inspiratory time is dependent on the heart rate (HR), the rate of the intrinsic respiratory oscillator (f(i)), and the strength of the hypothetical cardiovascular afferent impulse. Failure to account for HR and f(i) may explain the inconsistent results from previous attempts to identify the neural pathways involved in CVC. Cognizant of these interactions, we factored in the HR-to-f(i) ratio in our examination of the role of the vagus nerve and arterial baroreceptors in CVC by cardiac pacing 29 anesthetized Sprague-Dawley rats and incrementally changing the HR. With the assumption of a relatively constant f(i), CVC could be examined across a range of HR-to-f(o) ratios before and after vagotomy, sinoaortic denervation, and vagotomy + sinoaortic denervation. We confirmed the relation between CVC, HR-to-f(o) ratio, and breath-to-breath respiratory period variability and demonstrated the loss of these relations after baroreceptor elimination. Sham experiments (n = 8) showed that these changes were not due to surgical stress. Our data support the notion that inspiratory timing can be influenced by cardiac afferent activity. We conclude that the putative cardiovascular input arises from the arterial baroreceptors and that the vagus nerve is not critical for CVC.

Neurogenesis of cough, other airway defensive behaviors and breathing: A holarchical system?

Cough and breathing are generated by a common muscular system. However, these two behaviors differ significantly in their mechanical features and regulation. The current conceptualization of the neurogenic mechanism for these behaviors holds that the multifunctional respiratory pattern generator undergoes reconfiguration to produce cough. Our previous results indicate the presence of a functional cough gate mechanism that controls the excitability of this airway defensive behavior, but is not involved in the regulation of breathing. We propose that the neurogenesis of cough, breathing, and other nonbreathing behaviors is controlled by a larger network, of which the respiratory pattern generator is part. This network we term a holarchical system. This system is governed by functional control elements known as holons, which confer unique regulatory features to each behavior. The cough gate is an example of such a holon. Neurons that participate in a cough holon may include behavior selective elements. That is, neurons that are either specifically recruited during cough and/or tonically-active neurons with little or no modulation during breathing but with significant alterations in discharge during coughing. We also propose that the holarchical system is responsible for the orderly expression of different airway defensive behaviors such that each motor task is executed in a temporally and mechanically discrete manner. We further propose that a holon controlling one airway defensive behavior can regulate the excitability of, and cooperate with, holons unique to other behaviors. As such, co-expression of multiple rhythmic behaviors such as cough and swallow can occur without compromising airway defense.

Analyzing functional connectivity using a network likelihood model of ensemble neural spiking activity

Analyzing the dependencies between spike trains is an important step in understanding how neurons work in concert to represent biological signals. Usually this is done for pairs of neurons at a time using correlation-based techniques. Chornoboy, Schramm, and Karr (1988) proposed maximum likelihood methods for the simultaneous analysis of multiple pair-wise interactions among an ensemble of neurons. One of these methods is an iterative, continuous-time estimation algorithm for a network likelihood model formulated in terms of multiplicative conditional intensity functions. We devised a discrete-time version of this algorithm that includes a new, efficient computational strategy, a principled method to compute starting values, and a principled stopping criterion. In an analysis of simulated neural spike trains from ensembles of interacting neurons, the algorithm recovered the correct connectivity matrices and interaction parameters. In the analysis of spike trains from an ensemble of rat hippocampal place cells, the algorithm identified a connectivity matrix and interaction parameters consistent with the pattern of conjoined firing predicted by the overlap of the neurons' spatial receptive fields. These results suggest that the network likelihood model can be an efficient tool for the analysis of ensemble spiking activity.

Two enhancements of the gravity algorithm for multiple spike train analysis

The gravity method for neuronal assembly analysis represents each neuron as a particle in N-space with a time varying charge that is a filtered version of the corresponding spike train, with appropriate rules for forces between and movements of the charged particles. Resulting trajectories reflect neuronal timing relationships. The usual short time constants in the filter restrict aggregation to highly synchronized neurons and reduce the sensitivity for delayed correlations; long time constants in the filter reduce selectivity. Here we describe an enhancement that modifies rules for assigning charge increment times to allow mixtures of short and long lag correlations. Charge increments for each pair are offset from the actual spike times by time lags defined by features in corresponding cross-correlograms; no such charge offsets are invoked if the correlogram is flat. Tuning increases charge products and aggregation of long lag correlated pairs. A second enhancement uses a new three-dimensional display of particle pair trajectories to parse the type of neuronal relationship. For each pair, we record and display the inter-particle distance and the distance each particle moves from its original location in the N-space. The resulting trajectories cluster according to the type of interaction between the represented neurons. Results from simulated networks and in vivo multi-site recordings show that these modifications detect assembly properties not identified by the standard methods.

Increased noise level of purkinje cell activities minimizes impact of their modulation during sensorimotor control

While firing rate is well established as a relevant parameter for encoding information exchanged between neurons, the significance of other parameters is more conjectural. Here, we show that regularity of neuronal spike activities affects sensorimotor processing in tottering mutants, which suffer from a mutation in P/Q-type voltage-gated calcium channels. While the modulation amplitude of the simple spike firing rate of their floccular Purkinje cells during optokinetic stimulation is indistinguishable from that of wild-types, the regularity of their firing is markedly disrupted. The gain and phase values of tottering's compensatory eye movements are indistinguishable from those of flocculectomized wild-types or from totterings with the flocculus treated with P/Q-type calcium channel blockers. Moreover, normal eye movements can be evoked in tottering when the flocculus is electrically stimulated with regular spike trains mimicking the firing pattern of normal simple spikes. This study demonstrates the importance of regularity of firing in Purkinje cells for neuronal information processing.

Arterial pulse modulated activity is expressed in respiratory neural output

Although it is well-established that sympathetic activity is modulated with respiration, it is unknown whether neural control of respiration is reciprocally influenced by cardiovascular function. Even though previous studies have suggested the existence of pulse modulation in respiratory neurons, they could not exclude the possibility that such cells were involved in cardiovascular rather than respiratory motor control, owing to neuroanatomic and functional overlaps between brain stem neurons involved in respiratory and cardiovascular control. The aim of this study was to test the hypothesis that respiratory motoneurons and putative premotoneurons are modulated by arterial pulse. An existing data set composed of 72 well-characterized, respiratory-modulated brain stem motoneurons and putative premotoneurons was analyzed using delta(2), a recently described statistic that quantifies the magnitude of arterial pulse-modulated spike activity [Dick TE and Morris KF. J Physiol 556: 959-970, 2004]. Neuronal activity was recorded in the rostral and caudal ventral respiratory groups of 19 decerebrate, neuromuscular-blocked, ventilated cats. Axonal projections were identified by rectified and unrectified spike-triggered averages of recurrent laryngeal nerve activity or by antidromic activation from spinal stimulation electrodes. The firing rates of approximately 30% of these neurons were modulated in phase with both the respiratory and cardiac cycles. Furthermore, arterial pulse modulation occurred preferentially in the expiratory phase in that only expiratory neurons had high delta(2) values and only expiratory activity had significant delta(2) values after partitioning tonic activity into the inspiratory and expiratory phases. The results demonstrate that both respiratory motoneurons and putative premotoneuronal activity can be pulse modulated. We conclude that a cardiac cycle-related modulation is expressed in respiratory motor activity, complementing the long-recognized respiratory modulation of sympathetic nerve activity.

Quantitative analysis of cardiovascular modulation in respiratory neural activity

We propose the 'delta(2)-statistic' for assessing the magnitude and statistical significance of arterial pulse-modulated activity of single neurones and present the results of applying this tool to medullary respiratory-modulated units. This analytical tool is a modification of the eta(2)-statistic and, consequently, based on the analysis of variance. The eta(2)-statistic reflects the consistency of respiratory-modulated activity on a cycle-by-cycle basis. However, directly applying this test to activity during the cardiac cycle proved ineffective because subjects-by-treatments matrices did not contain enough 'information'. We increased information by dividing the cardiac cycle into fewer bins, excluding cycles without activity and summing activity over multiple cycles. The analysed neuronal activity was an existing data set examining the neural control of respiration and cough. Neurones were recorded in the nuclei of the solitary tracts, and in the rostral and caudal ventral respiratory groups of decerebrate, neuromuscularly blocked, ventilated cats (n= 19). Two hundred of 246 spike trains were respiratory modulated; of these 53% were inspiratory (I), 36.5% expiratory (E), 6% IE phase spanning and 4.5% EI phase spanning and responsive to airway stimulation. Nearly half (n= 96/200) of the respiratory-modulated units were significantly pulse modulated and 13 were highly modulated with delta(2) values exceeding 0.3. In 10 of these highly modulated units, eta(2) values were greater than 0.3 and all 13 had, at least, a portion of their activity during expiration. We conclude that cardiorespiratory interaction is reciprocal; in addition to respiratory-modulated activity in a subset of neuronal activity patterns controlling the cardiovascular system, pulse-modulated activity exists in a subset of neuronal activity patterns controlling the respiratory system. Thus, cardio-ventilatory coupling apparent in respiratory motor output is evident and, perhaps, derived from the neural substrate driving that output.

Repeated sequences of interspike intervals in baroresponsive respiratory related neuronal assemblies of the cat brain stem

Many neurons exhibit spontaneous activity in the absence of any specific experimental perturbation. Patterns of distributed synchrony embedded in such activity have been detected in the brain stem, suggesting that it represents more than "baseline" firing rates subject only to being regulated up or down. This work tested the hypothesis that nonrandom sequences of impulses recur in baroresponsive respiratory-related brain stem neurons that are elements of correlational neuronal assemblies. In 15 Dial-urethan anesthetized vagotomized adult cats, neuronal impulses were monitored with microelectrode arrays in the ventral respiratory group, nucleus tractus solitarius, and medullary raphe nuclei. Efferent phrenic nerve activity was recorded. Spike trains were analyzed with cycle-triggered histograms and tested for respiratory-modulated firing rates. Baroreceptors were stimulated by unilateral pressure changes in the carotid sinus or occlusion of the descending aorta; changes in firing rates were assessed with peristimulus time and cumulative sum histograms. Cross-correlation analysis was used to test for nonrandom temporal relationships between spike trains. Favored patterns of interspike interval sequences were detected in 31 of 58 single spike trains; 18 of the neurons with significant sequences also had short-time scale correlations with other simultaneously recorded cells. The number of distributed patterns exceeded that expected under the null hypothesis in 12 of 14 data sets composed of 4-11 simultaneously recorded spike trains. The data support the hypothesis that baroresponsive brain stem neurons operate in transiently configured coordinated assemblies and suggest that single neuron patterns may be fragments of distributed impulse sequences. The results further encourage the search for coding functions of spike patterns in the respiratory network.

Respiratory neuronal assemblies

This review describes results from in vivo experiments on brain stem network mechanisms that control breathing. Multi-array recording technology and computational methods were used to test predictions derived from simulations of respiratory network models. This highly efficient approach has the advantage that many simultaneously recorded neurons are subject to shared stimulus, history, and state-dependent conditions. Our results have provided evidence for concurrent or parallel network interactions in the generation and modulation of the respiratory motor pattern. Recent data suggest that baroreceptors, chemoreceptors, nociceptors, and airway cough receptors shape the respiratory motor pattern, at least in part, through a system of shared coordinated 'multifunctional' neurons distributed in the brain stem. The 'gravity method' for the analysis and representation of multi-neuron data has demonstrated respiratory phase-dependent impulse synchrony among neurons with no respiratory modulation of their individual firing rates. The detection of this emergent property motivated the development of pattern detection methods that subsequently identified repeated transient configurations of these 'correlational assemblies'. These results support the view that information can be 'coded' in the nervous system by spike timing relationships, in addition to firing rate changes that traditionally have been measured by neurophysiologists.

Respiratory neural activity during long-term facilitation

Intermittent hypoxia results in a long-term facilitation (LTF) of respiratory efferent activity. The studies reviewed here presented data from both anesthetized and decerebrate, paralyzed, vagotomized, artificially ventilated adult cats. Multiple arrays of tungsten microelectrodes were used to record the concurrent responses of brain stem neurons that contribute to respiratory motor pattern generation. Spike trains were analyzed with firing rate histograms, peristimulus time histograms, cycle triggered histograms, spike triggered averages with multiunit phrenic efferent activity, cross correlation histograms, joint peristimulus time histograms and the gravity method. These studies addressed several hypotheses. (1) There is parallel processing of input from carotid chemoreceptors to the brain stem. (2) Respiratory related midline neurons are involved in the induction and maintenance of LTF. (3) There is a change in effective connectivity of brain stem neurons with LTF. (4) Neural networks involved in the induction and maintenance of LTF have patterns of synchrony that recur with a frequency greater than expected by chance.

Transient configurations of baroresponsive respiratory-related brainstem neuronal assemblies in the cat

The regulation of gas exchange requires coordination of the respiratory and cardiovascular systems. Previous work suggested that medullary raphe neurones transform and transmit information from baroreceptors to neurones in the ventral respiratory group. This study tested the hypothesis that distributed brainstem neuronal assemblies are transiently reconfigured during the respiratory cycle and baroreceptor stimulation. Blood pressure was perturbed by intravenous injection of an alpha1-adrenergic receptor agonist, unilateral pressure changes in the carotid sinus, or occlusion of the descending aorta in 14 Dial-urethane anaesthetized, vagotomized, paralysed, artificially ventilated cats. Neurones were monitored simultaneously with microelectrode arrays in two or more of the following sites: n. raphe obscurus, n. raphe magnus, rostral and caudal ventrolateral medulla, and the nucleus tractus solitarii. Transient configurations of baroresponsive assemblies were detected with joint pericycle-triggered histograms, the gravitational representation, and related pattern detection methods. Data were also analysed with cycle-triggered histograms, peristimulus-time and cumulative sum histograms, cross-correlograms, spike-triggered averages of efferent phrenic activity, and joint impulse configuration scatter diagrams (snowflakes). Five to nine simultaneously recorded spike trains from control expiratory phases were compared with data from interleaved equal-duration time blocks from control inspiratory phases. In each of seven animals, significant impulse synchrony detected by gravity analysis was confined to one phase of the respiratory cycle. Repeated patterns of distributed synchrony confined to periods of altered baroreceptor activity were detected and involved neurones that individually did not change firing rates during stimulation. Snowflakes and logical cross-correlation analysis provided evidence for the cooperative actions of impulses in concurrently active parallel channels. In 12 of 17 pairs of neurones with at least one baroresponsive cell, joint pericycle-triggered histograms detected synchrony indicative of shared inputs or functional excitatory interactions that varied as a function of time in the respiratory cycle. Neurones in four of the pairs had no respiratory modulation of their individual firing rates. Data from eight other pairs were indicative of fluctuations in inhibition during the respiratory cycle. The results demonstrate repeated transient configurations of baroresponsive neuronal assemblies during the respiratory cycle, without concomitant firing rate changes in the constituent neurones, and suggest distributed network mechanisms for the modulation of baroreceptor-mediated reflexes.

Responses of simultaneously recorded respiratory-related medullary neurons to stimulation of multiple sensory modalities

This study addresses the hypothesis that multiple afferent systems share elements of a distributed brain stem network that modulates the respiratory motor pattern. Data were collected from 18 decerebrate, bilaterally vagotomized, paralyzed, artificially ventilated cats. Up to 28 neurons distributed in the rostral and caudal ventral respiratory group, nucleus tractus solitarius, and raphe obscurus were recorded simultaneously with microelectrode arrays. Phases of the respiratory cycle and inspiratory drive were assessed from integrated efferent phrenic nerve activity. Carotid chemoreceptors were stimulated by injection of CO2-saturated saline solution via the external carotid artery. Baroreceptors were stimulated by increased blood pressure secondary to inflation of an embolectomy catheter in the descending aorta. Cutaneous nociceptors were stimulated by pinching a footpad. Four hundred seventy-four neurons were tested for respiratory modulated firing rates and responses; 403 neurons were tested with stimulation of all 3 modalities. Chemoreceptor stimulation and pinch, perturbations that tend to increase respiratory drive, caused similar responses in 52 neurons; 28 responded oppositely. Chemoreceptor and baroreceptor stimulation resulted in similar primary responses in 45 neurons; 48 responded oppositely. Similar responses to baroreceptor stimulation and pinch were recorded for 38 neurons; opposite effects were measured in 26 neurons. Among simultaneously recorded neurons, distinct combinations of firing rate changes were evoked in response to stimulation of the different modalities. The results show a functional convergence of information from carotid chemoreceptors, baroreceptors, and cutaneous nociceptors on respiratory-modulated neurons distributed in the medulla. The data are consistent with the hypothesis that brain stem neurons have overlapping memberships in multifunctional groups that influence the respiratory motor pattern.