Phasic Vagal Sensory Feedback Transforms Respiratory Neuron

The Journal of Neuroscience, September 15, 2001, 21(18):7363–7371
Phasic Vagal Sensory Feedback Transforms Respiratory Neuron
Activity In Vitro
Nicholas M. Mellen and Jack L. Feldman
Department of Neurobiology, University of California Los Angeles, Los Angeles, California 90095-1763
The isolated neonatal rat medulla generates respiratory-related
rhythms recorded from cervical spinal cord ventral roots. When
lungs and their vagal innervation are retained, respiratory activity is modulated by lung mechanoreceptor feedback: transient lung inflation triggered off inspiratory onset (phasic inflation) shortens inspiration and increases respiratory frequency.
In this study, the activity of six respiratory neuron classes before
and during phasic inflation was studied. Type 1 and 2 inspiratory neurons, identified in the transverse slice, were distinguished by the presence of a transient outward current or a
hyperpolarization-activated inward current, respectively. Cell
types only identified in the en bloc medulla included type II and
III inspiratory neurons, distinguished by delayed onset and
peri-inspiratory inhibition, respectively, and preinspiratory neurons, active before and after but silent during inspiration. Biphasic neurons, identified in the preparation used here, fired
briskly during lung inflation but are otherwise quiescent. During
phasic inflation, biphasic neurons showed a decrementing expiratory pattern of activity, matched by augmented postinspiratory hyperpolarization in type 1 neurons only, suggesting that
biphasic neurons inhibit type 1 neurons, removing drive to other
inspiratory neurons and terminating the inspiratory burst. This
mechanism could account for a phasic inflation-induced increase in respiratory frequency via resetting effects. Alternatively, the phasic inflation-induced respiratory frequency increase may be attributable to slow facilitation. Slow modulation
consistent with facilitation was apparent in the earlier onset of
pre-I firing before inspiration and loss of postinspiratory firing
and in the earlier onset of depolarization in type 2 neurons. On
the basis of relative onset times and responses to phasic
inflation, connectivity between these cell types is proposed.
In mammals, the hypothesized kernel for respiratory rhythmogenesis is localized to the pre-Bötzinger complex (pre-BötC)
ventral and just caudal to the compact formation of the nucleus
ambiguus in the rostroventrolateral medulla (Smith et al., 1991;
Rekling and Feldman, 1998; Gray et al., 1999). In transverse
brainstem slices containing the pre-BötC, respiratory-related
rhythmic motor output can be obtained (Smith et al., 1991).
Although these preparations establish that pre-BötC networks
generate stable rhythmic activity, how these circuits function as
part of a more complete network is for the most part unknown. In
the intact brainstem, the excitability of respiratory neurons in the
pre-BötC is modulated by medullary circuits rostral and caudal to
the pre-BötC and by sensory feedback. In the more intact in vitro
en bloc brainstem preparation (Suzue, 1984; Smith and Feldman,
1987), neurons within the pre-BötC have patterns of activity
(Onimaru et al., 1997) not seen in the slice, presumably because
of network inputs absent in the slice.
We investigated the effects of afferent feedback in an in vitro
neonate rat brainstem–spinal cord preparation retaining the
lungs and their vagal innervation (Murakoshi and Otsuka, 1985;
Mellen and Feldman, 1997). In this preparation, respiratory frequency and pattern are modulated by lung inflation in a manner
congruent with effects in vivo: midexpiratory lung inflation
lengthens expiration [Breuer-Hering expiratory reflex (BHE);
Breuer, 1868; Murakoshi and Otsuka, 1985; Mellen and Feldman,
1997], and inflation during inspiration terminates inspiration
(Breuer-Hering inspiratory reflex; Breuer, 1868; Murakoshi and
Otsuka, 1985; Mellen and Feldman, 1997, 2000). In addition,
when lung inflation concurrent with inspiratory onset is applied in
consecutive cycles (henceforth referred to as phasic inflation),
respiratory frequency increases (Mellen and Feldman, 2000). The
effect on frequency declines monotonically as the delay between
inspiratory burst onset and inflation onset increases; i.e., the
effect is phase-dependent (Mellen and Feldman, 2000). Because
these responses are obtained with lung pressure changes within
the physiological range (Widdicombe, 1961), slowly adapting pulmonary receptor afferents (SARs) are likely being selectively
activated.
The pathway by which mechanoreceptor feedback modulates
respiratory frequency has been partially characterized in vivo:
SARs excite neurons in the medial nucleus tractus solitarius
(Bonham et al., 1993), which project to the ventrolateral medulla
to excite subsets of expiratory neurons in and rostral to the
pre-BötC. In turn, inspiratory neurons are inhibited (Hayashi et
al., 1996).
In this study, phasic inflation increased respiratory frequency
and shortened inspiratory duration. We observed consistent
changes in the activity of six distinct classes of pre-BötC respiratory neurons. These responses constrain the possible cellular and
synaptic mechanisms underlying the generation and modulation
of respiratory rhythm and pattern and connectivity of identified
respiratory neuron types.
To optimally classify respiratory neurons based on the existing
literature, it was necessary to use two distinct nomenclatures. The
Received Feb. 14, 2000; revised May 30, 2000; accepted July 3, 2001.
This research was funded by National Institutes of Health Grants HL40959 and
HL37941 and by Research Grant RG-105-N from the American Lung Association.
Correspondence should be addressed to Nicholas M. Mellen, Department of
Neurobiology, Box 951763, University of California Los Angeles, Los Angeles, CA
90095-1763. E-mail: [email protected].
Copyright © 2001 Society for Neuroscience 0270-6474/01/217363-09$15.00/0
Key words: mammal; neonate rat; respiration; pulmonary
afferents; vagus; medulla; pre-Bötzinger complex
7364 J. Neurosci., September 15, 2001, 21(18):7363–7371
first classification scheme, developed in the mouse transverse slice
(Rekling et al., 1996), differentiates between two types of inspiratory neurons based on the presence of a transient outward current
(IA; type 1) or a hyperpolarization-activated inward current (Ih;
type 2). The second classification scheme, developed in the en
bloc medulla, identifies three neuron types not reported in the
slice on the basis of their firing patterns (Onimaru et al., 1997):
neurons with delayed, steeply rising onset of inspiratory activity
(type II), neurons with peri-inspiratory inhibition (type III), and
neurons active before and after but inhibited during inspiration
(pre-I).
MATERIALS AND METHODS
Dissection
Neonatal Sprague Dawley rats (postnatal days 0 –3; n ⫽ 36) were used. In
accordance with methods approved by the Institutional Animal C are and
Use Committee, rat pups were cooled to 5°C, decerebrated immediately
rostral to the superior colliculus, and transferred to a bath continuously
perf used with artificial C SF (AC SF) containing (in mM): 128.0 NaC l, 3.0
KC l, 1.5 C aC l2, 1.0 MgSO4, 21.0 NaHC O3, 0.5 NaH2PO4, and 30.0
glucose, equilibrated with 95% O2 and 5% C O2, at 5°C. As described
previously (N. M. Mellen and J. L. Feldman, unpublished procedures),
the brainstem and spinal cord, together with the lungs and the intact
vagus nerve, were isolated and transferred to a chamber warmed to 27°C.
The dura was removed, and the medulla was transected at the level of the
facial nucleus. A cannula filled with AC SF was inserted into the trachea
and held in place with a suture. Pressure changes to the lungs were
applied using a computer-controlled precision syringe pump (C arnegie
Medecin M100). Viability of the vagal afferent pathway was tested with
sustained midexpiratory inflation (pressure, 2–5 mmH2O), which elicited
expiratory lengthening.
Recording methods
Ventral root recordings. Inspiration-related activity from ventral root C2
was recorded using a saline-filled glass suction electrode (100 k⍀).
Signals were amplified 15,000 –30,000 times and bandpass-filtered (0.3–3
kHz) using Grass P5 differential amplifiers (Grass Instruments, Quincy,
M A).
Intracellular recordings. Blind whole-cell patch recordings (Blanton et
al., 1989) were performed. Neurons were recorded at or just rostral to the
most rostral hypoglossal rootlet at depths of 200 – 600 ␮m, thus within the
typical boundaries of the transverse slice used to isolate the pre-BötC
(Smith et al., 1991). Electrodes with resistances ranging from 4 to 12 M⍀,
pulled from filamented glass capillary tubing (1.5 mm outer diameter,
0.86 mm inner diameter; A-M Systems) and filled with (in mM): 126
CH4SO3, 126 KOH, 15 NaC l, 1 MgC l2, 0.01 C aC l2, 10 H EPES, 3
2ATP(Mg 2⫹), and 0.1 bis(2-aminophenoxy)ethane-N,N,N⬘,N⬘-tetraacetic acid, were used.
Recordings were performed in current-clamp bridge mode using the
Axoclamp-2A amplifier (Axon Instruments) with a 0.01 gain head stage.
The electrode was quickly advanced to a depth ⬃200 ␮m below the
surface under positive pressure. Thereafter, while passing a highfrequency 120 pA square wave, the electrode was slowly advanced. When
a cell was contacted, the voltage deflection attributable to the current
pulse would increase sharply. Bias current was then applied to hold the
electrode at ⫺90 mV, and the membrane was ruptured by applying
negative pressure to the electrode tip. Access resistance was 10 –20 M⍀,
based on corrected bridge balance under current clamp. Only neurons
with Vm of less than ⫺40 mV and action potential overshoot were
analyzed. Access resistance was balanced before applying bias currents.
A voltage deflection associated with a 2 sec 20 pA hyperpolarizing pulse,
applied in midexpiration, was used to estimate input resistance. A liquid
junction potential of ⫹8 mV was not compensated. Step currents were
applied using a pulse generator (AM PI Master 8) triggered off the
instrumentation control computer (see below).
Data acquisition, signal processing, and
instrumentation control
C ervical ventral root activity, respiratory neuron membrane potential,
applied bias currents, and syringe pump control signals were digitized at
20 kHz (AT-M IO-64-E3; National Instruments, Austin, TX) and written
to a hard disk using software developed in LabView (National Instru-
Mellen and Feldman • Lung Afferents Transform Rat Respiratory Network In Vitro
ments). In addition, time stamps of inspiratory onsets, lung inflation, and
deflation were stored.
A second computer was used to inflate lungs, with onset of inflation
triggered off inspiratory onset (delay, ⬍100 msec), using an analog-todigital–digital-to-analog board (AT-MIO-16E-10; National Instruments).
Experimental protocol
Lung inflation protocols. Before recording data, the minimum injectate
volume necessary to obtain consistent reflex effects was identified. This
was done by applying lung inflations triggered off inspiratory onset, with
gradually increasing injectate volumes, beginning with 0.2 ml and never
⬎0.4 ml, at ⬃0.1 ml /sec. Because inflation triggered off inspiratory onset
was immediately followed by deflation, pressure changes could not be
monitored, because resistance to flow through the cannula produced
large pressure transients. Under conditions in which inflation is maintained (Mellen et al., 1997), this volume of injectate produced pressure
changes of 2–5 mmH2O, within the physiological range (4 –7 mm H2O;
Widdicombe, 1961). This protocol of phasic inflation was maintained for
at least 10 cycles. Together, the ⬎10 control cycles and the ⬎10 phasic
inflation cycles obtained for each cell in this study are referred to as a
bout.
Measurement of membrane properties. Delayed excitation, i.e., delay in
action potential onset in response to step depolarization consistent with
the presence of IA (Connor and Stevens, 1971), and sag rebound response
consistent with a Ih (Pape and McCormick, 1989), was tested for. By
triggering current pulses 2– 4 sec after inspiratory onset, i.e., in midexpiration, the contribution of phasic peri-inspiratory synaptic inputs to the
membrane potential in these tests was minimized.
To test for IA, a hyperpolarizing bias current sufficient to bring the
membrane potential below ⫺80 mV was applied, the bridge was balanced, and a depolarizing step current was applied sufficient to elicit
tonic impulse activity. To test for Ih, a series of step hyperpolarizations
were applied from the resting membrane potential. In addition, this
protocol tested for the presence of low-voltage-activated C a 2⫹
conductances.
Data analysis
Motor output anal ysis. For both control and test cycles within a bout, the
cycle period was defined as the interval from onset of one inspiratory
burst to the next. Estimates of burst onset time were insensitive to the
threshold set for inspiratory burst detection, because ventral root inspiratory onset was abrupt. Within each bout, mean control (n ⬎ 10) and test
(n ⬎ 10) cycle periods were calculated. Changes in burst duration and
respiratory period accompanying phasic inflation were tested using a
paired t test on bout means within Origin (Microcal).
Burst-trig gered average anal ysis. In a previous study of inspiratory
neurons, onset of inspiratory neuronal depolarization was used to estimate relative onset times (Rekling et al., 1996). By this criterion, activity
onset of some inspiratory neuron classes preceded motor nerve output.
Here, we used the steep rising portion of the averaged inspiratory burst
as the reference point (Fig. 1), because we found that it was less sensitive
to measurement error, recording conditions, and threshold. Thus, in
addition to characterizing membrane properties as described above, we
quantified inspiratory neuron activity from averages triggered off the
cervical ventral root inspiratory burst onset (C2 burst) (Fig. 1, gray
dashed line), with the width of the averaging window varying from ⫺2 to
⫹ 3 s.
For all neurons, calculation of depolarization or hyperpolarization
onset and duration (Fig. 1, top, horizontal arrows) were measured relative
to the half-height of the depolarization (for inspiratory neurons) or
hyperpolarization (for pre-I or expiratory neurons) accompanying ventral root inspiratory burst (Fig. 1, vertical double arrow).
Onset and offset of preinspiratory or postinspiratory modulation, i.e.,
preinspiratory and postinspiratory depolarization in pre-I neurons and
preinspiratory hyperpolarization in type III inspiratory neurons, as
shown in Figure 1, were also referenced to ventral root inspiratory onset
(Fig. 1, bottom, horizontal double arrow), using resting membrane potential as the threshold for onset and offset. Measurements in the amplitude
domain were made relative to deflections from each the resting membrane potential of each cell. For hyperpolarizations, the voltage change
from resting potential to hyperpolarization minimum was measured (Fig.
1, bottom, vertical arrow). For preinspiratory depolarization, amplitude
was estimated at the midpoint between onset of depolarization from
baseline to ventral root inspiratory onset. For each cell, measurements
were made from averages of at least 5 cycles. Timing and amplitude
Mellen and Feldman • Lung Afferents Transform Rat Respiratory Network In Vitro
J. Neurosci., September 15, 2001, 21(18):7363–7371 7365
0.01) (Fig. 2 B). The delayed effect was the advanced onset of the
subsequent inspiratory burst, which followed inflation offset by
seconds. Thus, respiratory periods with (7.5 ⫾ 0.5 sec) and
without (9.5 ⫾ 0.5 sec) phasic inflation were significantly different
( p ⬍ 0.01) (Fig. 2C).
Thirty-three respiratory-modulated neurons were recorded,
sorted into six groups: biphasic neurons, type 1 and 2 inspiratory
neurons (Rekling et al., 1996), type II neurons, type III inspiratory neurons, and pre-I neurons (Onimaru et al., 1997). Periinspiratory firing patterns, relative onset times, and inflationinduced changes in activity for these cell types were as follows.
Biphasic neurons (n ⴝ 7)
Figure 1. Schematic of method for measuring timing and amplitude of
respiratory neuron modulation during the inspiratory burst. Top trace,
Rectified ventral root activity (C2) averaged using ventral root burst onset
as the trigger (vertical dashed line). Bottom trace, Type III inspiratory
neuron activity triggered off ventral root burst. Duration (top horizontal
arrow) and activity onset (second horizontal arrow) of inspiration-related
activity in respiratory neurons are measured from C2 burst onset to
half-maximal depolarization (types 1, 2, II, and III neurons) or hyperpolarization (pre-I neurons); these points fall on the line labeled half-height.
Preinspiratory or postinspiratory deflections in membrane potential are
measured relative to baseline membrane potential indicated by the line
labeled resting Vm. Here, preinspiratory hyperpolarization onset and
offset are measured relative to C2 burst onset (hyp duration, bottom,
double arrow); the magnitude of hyperpolarization is measured relative to
resting membrane potential (hyp amplitude, bottom, vertical arrow).
measurements within each cell type were averaged; all bar or column
graphs display means ⫾ SE. For each cell type, the statistical significance
of inflation-induced changes in activity was tested using paired t tests on
means.
Cycle-triggered histogram generation. To quantify spike frequency over the
course of the respiratory cycle, cycle-triggered histograms (CTHs) were
generated from at least six control and test cycles. Each cycle was divided
into 30 bins, and the number of spikes in each bin was counted. Spike
frequency in each bin was obtained by dividing the number of spikes in each
bin by the bin duration. The first 15 bins of the subsequent cycle were
appended to each cycle, so that the CTHs displayed activity over 1.5 cycles
in 45 bins. Spike frequency was averaged across cycles by bin.
Comparison of onset times. We used the steeply rising (or falling) phase
of the burst-triggered average to estimate relative inspiratory-related
activity onset. Because individual spikes preceding maximal inspiratory
firing contributed little to the burst-triggered averages, estimates of onset
times were skewed late (see Results). Because this measurement error
was consistent across cell types, ventral root inspiratory onset served as
an index to sort cell types by onset delay. We tested for differences in
onset delay using a two-way fixed-effect ANOVA on delay means, in
which one fixed effect was cell type and the other fixed effect was control
versus phasic inflation. The SAS procedure M I X ED (SAS Institute,
C ary, NC) was used for the analysis.
RESULTS
Phasic lung inflation modulated respiratory frequency and burst
pattern (Mellen and Feldman, 2000) (Fig. 2 A,B). For motor nerve
output, modulatory effects can be divided into “immediate” and
“delayed.” The immediate effect of phasic inflation was significant
inspiratory burst shortening in the same cycle, from 690 ⫾ 80
msec in control cycles to 470 ⫾ 50 msec with phasic inflation ( p ⬍
Biphasic neurons had a resting membrane potential of ⫺54 ⫾ 3
mV and an input resistance of 700 ⫾ 190 M⍀. In control cycles,
these neurons were silent or fired tonically at a low frequency
(Fig. 3A,C); during inspiration they were weakly excited and then
inhibited (Fig. 3B). With phasic inflation, these neurons were
strongly excited and fired in a decrementing pattern (Fig. 3B),
reaching a maximal frequency of 29 ⫾ 1 spikes/sec during inflation, and returned to baseline at varying rates (Fig. 3C). Thus,
although firing onset during phasic inflation showed little dispersion across neurons (560 ⫾ 80 msec after inspiratory onset), firing
duration varied considerably (2900 ⫾ 760 msec).
Type 1 inspiratory neurons (n ⴝ 8)
Type 1 inspiratory neurons (Rekling et al., 1996) have a ramplike expiratory depolarization and IA (Connor and Stevens,
1971). These cells had a resting membrane potential of ⫺53 ⫾ 2
mV and an input resistance of 640 ⫾ 110 M⍀. Phasic inflation
resulted in strong, transient postinspiratory inhibition (Fig. 4 A,
arrows, B). Inflation hyperpolarized these cells by ⫺5.9 ⫾ 0.8 mV
(n ⫽ 8) (Fig. 4 B) and was transient, lasting 710 ⫾ 110 msec (n ⫽
8). Inflation significantly shortened inspiratory firing ( p ⬍ 0.01)
from 1000 ⫾ 120 to 460 ⫾ 60 msec, a reduction of 51%. These
neurons began to depolarize before inspiratory onset, and the
steep portion of the rise to firing occurred 11 ⫾ 14 msec after
ventral root inspiratory onset. Onset of firing relative to ventral
root inspiratory onset changed little with phasic inflation (8 ⫾ 5
msec) (Fig. 4 B).
Type 2 inspiratory neurons (n ⴝ 4)
The defining features of type 2 inspiratory neurons (Rekling et
al., 1996) are Ih (Pape and McCormick, 1989) and a flat membrane potential during expiration with low-frequency tonic spiking. Here, EPSPs rather than spikes were observed during expiration, perhaps because recordings were performed without the
elevated excitability (via elevated extracellular [K ⫹], lowered
extracellular [Ca 2⫹], or both) necessary in the transverse slice
(Fig. 5A). These neurons had a resting membrane potential of
⫺58 ⫾ 2 mV and an input resistance of 440 ⫾ 160 M⍀ and began
to depolarize before and fired shortly after ventral root inspiratory onset. In control cycles, this depolarization began 480 ⫾ 250
msec before ventral root inspiratory onset (Fig. 5B). With phasic
inflation, depolarization onset began significantly earlier (720 ⫾
270 msec; p ⬍ 0.05), but preinspiratory depolarization, measured
at the midpoint between onset of preinspiratory depolarization
and inspiratory discharge, was not significantly greater with phasic inflation (2.0 ⫾ 1.5 mV) than control (0.6 ⫾ 0.2 mV). The
steep phase of inspiratory depolarization followed ventral root
inspiratory onset by 30 ⫾ 40 msec in control cycles and by 45 ⫾
50 msec in cycles with phasic inflation. Inflation significantly
7366 J. Neurosci., September 15, 2001, 21(18):7363–7371
Mellen and Feldman • Lung Afferents Transform Rat Respiratory Network In Vitro
Figure 2. A, Top, Rectified integrated C2 population
activity from one experiment (␶ ⫽ 20 msec). Bottom,
Raster plot of periods from five experiments before
(left of gray dashed line) and during (right of gray
dashed line) phasic inflation. B, C, Paired t tests on
means reveal that phasic inflation significantly shortens burst duration (700 ⫾ 80 vs 470 ⫾ 50 msec; p ⬍
0.01; B) and cycle period (9.5 ⫾ 0.5 vs 7.5 ⫾ 0.5 sec;
p ⬍ 0.01; C).
Figure 3. Biphasic neurons fire briskly during phasic inflation but are silent in control cycles. A, Raw traces of biphasic neuron membrane potential in
control cycles (top) and with phasic inflation (bottom). In each panel, rectified integrated C2 ventral root activity (␶ ⫽ 10 msec) is shown in the top trace,
and biphasic neuron membrane potential is shown in the bottom trace. Transient inflations are indicated by gray boxes. The same format is used in Figures
4 – 8. B, Top, Average of biphasic neuron membrane potential over six cycles, triggered off ventral root inspiratory onset (vertical dashed line) with (black)
and without ( gray) phasic inflation in one biphasic neuron. Bottom, Bar graph of average inflation-induced firing duration in biphasic neurons (4 ⫾ 0.8
sec; n ⫽ 7), scaled to the x-axis of the burst-triggered average. The same format is used in Figures 4 – 8. Note excitatory and inhibitory drive during
inspiration in control cycles and brisk firing in the presence of phasic inflation. Both firing onset and firing offset of each neuron were calculated relative
to ventral root inspiratory onset. C, Cycle-triggered histograms (30 bins per cycle) of two biphasic neurons. The first neuron fired tonically at low rates
in control cycles ( gray) and burst at a high rate during phasic inflation followed by slow adaptation afterward (dotted line), whereas the second neuron
showed faster adaptation (solid line).
shortened inspiratory firing ( p ⬍ 0.01) from 1130 ⫾ 100 to 690 ⫾
80 msec, a reduction of 48%.
Type II neurons (n ⴝ 6)
In type II neurons (Onimaru et al., 1997) onset of inspiratory
depolarization was rapid and steep (Fig. 6 A) and followed ventral
root inspiratory onset by 150 ⫾ 50 msec in control cycles (Fig.
6 B). These cells had a resting membrane potential of ⫺57 ⫾ 3
mV and an input resistance of 730 ⫾ 170 M⍀. To test whether
burst onset delay was caused by preinspiratory hyperpolarization,
hyperpolarizing bias currents were applied to hold neurons below
⫺85 mV; no outward current reversal was observed (Fig. 6 B).
With phasic inflation, inspiratory depolarization began 190 ⫾ 80
msec after ventral root inspiratory onset. For each neuron, onset
of firing relative to ventral root inspiratory onset, which coincided
with peak ventral root activity, showed very little variability.
Phasic inflation shortened inspiratory burst duration from 920 ⫾
260 msec in control cycles to 680 ⫾ 140 msec, a reduction of 26%.
No inflation-induced inhibition was apparent (Fig. 6 B).
Type III inspiratory neurons (n ⴝ 5)
The defining feature of type III inspiratory neurons (Onimaru et
al., 1997) is peri-inspiratory inhibition (Fig. 7A). These neurons
had a resting membrane potential of ⫺56 ⫾ 3 mV and an input
resistance of 510 ⫾ 170 M⍀. Phasic inflation significantly advanced (n ⫽ 5; p ⬍ 0.05) onset of preinspiratory inhibition from
440 ⫾ 110 to 580 ⫾ 150 msec before ventral root inspiratory onset
(Fig. 7B) and delayed onset of inspiratory firing from 110 ⫾ 40
Mellen and Feldman • Lung Afferents Transform Rat Respiratory Network In Vitro
J. Neurosci., September 15, 2001, 21(18):7363–7371 7367
Figure 4. Type 1 neurons are transiently inhibited
during phasic inflation but show no change in preinspiratory activity. A, Raw traces of type 1 neuron
membrane potential with and without phasic inflation
are displayed as in Figure 3. In control cycles, type 1
neurons have a ramp-like depolarization during expiration and can begin to spike before ventral root
inspiratory onset. Although preinspiratory trajectories with and without inflation coincide closely, in the
presence of phasic inflation there is pronounced
postinspiratory hyperpolarization coinciding with inflation (arrows). B, Burst-triggered average of seven
cycles from one type 1 neuron, displayed as in Figure
3. Bar and column graphs show mean values for all
eight type 1 neurons and are scaled to the axes of the
averaged traces. Phasic inflation causes postinspiratory inhibition; the black column to the right shows
average postinspiratory inhibition in type 1 neurons
(⫺5.8 ⫾ 0.7 mV), which lasted 710 ⫾ 100 msec (bottom, black bar). Phasic inflation significantly shortens
( p ⬍ 0.01) type 1 inspiratory firing duration, from
1000 ⫾ 120 msec (top, white bar) to 460 ⫾ 60 msec
(top, black bar).
Figure 5. Preinspiratory depolarization in type 2
neurons begins earlier with phasic inflation, but
postinspiratory activity is unchanged. A, Raw traces
of type 2 neuron membrane potential with and without phasic inflation are displayed as in Figure 3. In
control cycles, type 2 neuron membrane potential
varies little during expiration and begins to depolarize
before ventral root inspiratory onset. B, Bursttriggered average of seven cycles from one type 2
neuron, displayed as in Figure 3. Bar and column
graphs show mean values for all three type 2 neurons
and are scaled to the axes of the averaged traces.
Although postinspiratory trajectories with and without inflation coincide closely, in the presence of phasic
inflation, preinspiratory depolarization begins significantly earlier ( p ⬍ 0.05; bottom, black bar) than in
control cycles (bottom, white bar). Mean inflationinduced preinspiratory depolarization amplitude during phasic inflation (left, black column) was not significantly different from control (left, white column).
Phasic inflation significantly shortens ( p ⬍ 0.01) type
2 inspiratory firing duration, from 1160 ⫾ 130 msec
(top, white bar) to 670 ⫾ 110 msec (top, black bar).
msec after ventral root inspiratory onset in control cycles to
170 ⫾ 40 msec after ventral root inspiratory onset (Fig. 7B).
Preinspiratory inhibition amplitude was unchanged by phasic
inflation. Inflation significantly shortened inspiratory firing ( p ⬍
0.01) from 1440 ⫾ 170 msec (Fig. 7B) to 890 ⫾ 80 msec (Fig. 7B,
top), a reduction of 37%. Postinspiratory inhibition was not a
consistent feature of these neurons and was not changed by phasic
inflation.
Pre-I neurons (n ⴝ 6)
In control cycles, pre-I neurons were active before and after but
hyperpolarized during inspiration (Onimaru et al., 1997). These
neurons had a resting membrane potential of ⫺47 ⫾ 2 mV in
midexpiration and an input resistance of 520 ⫾ 160 M⍀. Because
the duration of preinspiratory and postinspiratory activity varied
from cycle to cycle, tests for Ih and IA could not be performed
without synaptic blockade. In control cycles, pre-I neurons began
firing 640 ⫾ 100 msec before and fell silent 130 ⫾ 30 msec after
ventral root inspiratory onset (Fig. 8 B). In the presence of phasic
inflation, preinspiratory firing began significantly earlier (1250 ⫾
220 msec; p ⬍ 0.05) but fell silent at the same delay relative to
ventral root inspiratory onset as in control cycles (130 ⫾ 30 msec)
(Fig. 8 B). To test whether loss of postinspiratory firing during
phasic inflation was caused by postinspiratory synaptic inhibition,
hyperpolarizing bias currents were applied to shift the resting
potential to ⫺75 mV. No outward current reversal was observed
(Fig. 8 B). In control cycles, pre-I neurons resumed firing 1180 ⫾
100 msec after inspiratory offset and fired for 2000 ⫾ 260 msec
(Fig. 8B).
Comparison of relative onset times
Phasic inflation did not significantly alter relative onset times for
any of the neuron classes ( p ⫽ 0.6). Significant differences in
onset times between cell types were found, however (Fig. 9A).
Type 1 neuron onsets were significantly earlier than all other cell
types with the exception of type 2 neurons; type 2 neuron onsets
were significantly earlier than type II and III neurons but not
pre-I neurons; type II neuron onsets were not significantly different from either type III neurons or pre-I neurons; and pre-I
neuron onsets were significantly earlier than type III neurons.
DISCUSSION
In intact animals, central circuits regulating respiratory frequency
and amplitude operate under closed loop conditions: rhythmo-
7368 J. Neurosci., September 15, 2001, 21(18):7363–7371
Mellen and Feldman • Lung Afferents Transform Rat Respiratory Network In Vitro
Figure 6. Type II neuron inspiratory burst duration
is shortened by phasic inflation, but preinspiratory
and postinspiratory activity are unchanged. A, Raw
traces of type II neuron membrane potential with and
without phasic inflation are displayed as in Figure 3.
In both control and phasic inflation cycles, type II
neuron membrane potential is flat during expiration
and rises abruptly during inspiration. B, Bursttriggered average of eight cycles from one type II
neuron, displayed as in Figure 3. Bar graphs show
mean values for all four type II neurons and are scaled
to the axes of the averaged traces. In the averaged
trace, the slope of inspiration-related depolarization
remains nearly vertical, indicating tight coupling between motor output and type II activity. Firing onset
during phasic inflation (bottom, black bar) and in
control cycles (bottom, white bar) is unchanged. Preinspiratory membrane potential remained flat after
application of bias currents to hold the cell at ⫺87
mV; thus delayed onset to firing is not attributable to
preinspiratory inhibition. Phasic inflation significantly
shortens ( p ⬍ 0.05) type II inspiratory firing duration,
from 1010 ⫾ 120 msec (top, white bar) to 640 ⫾ 130
msec (top, black bar).
Figure 7. Preinspiratory hyperpolarization begins
earlier with phasic inflation in type III neurons. A,
Raw traces of type III neuron membrane potential
with and without phasic inflation are displayed as in
Figure 3. In both control and phasic inflation cycles,
type III neuron membrane potential shows consistent
preinspiratory inhibition. B, Burst-triggered average
of seven cycles from one type III neuron, displayed as
in Figure 3. Bar and column graphs show mean values
for all five type II neurons and are scaled to the axes
of the averaged traces. Preinspiratory hyperpolarization onset occurs significantly earlier ( p ⬍ 0.05) during phasic inflation (bottom, black bar) than in control
cycles (bottom, white bar). Preinspiratory hyperpolarization amplitude (right, black column) and inspiratory burst onset are the same with and without phasic
inflation. Phasic inflation significantly shortens ( p ⬍
0.01) type III inspiratory firing duration, from 1450 ⫾
170 msec (top, white bar) to 890 ⫾ 80 msec (top, black
bar).
genic circuits drive phrenic and other respiratory pump motoneurons resulting in lung inflation; this activates pulmonary afferents,
which project centrally to rhythmogenic circuits. Removal of this
sensory feedback by bilateral vagotomy in vivo significantly slows
respiratory frequency (Smith et al., 1990; Bruce, 1997). We have
approximated these closed loop conditions in vitro by transiently
inflating the lungs triggered by ventral root inspiratory onset. We
observed changes congruent with effects in vivo (Feldman and
Gautier, 1976): (1) inspiratory shortening, which immediately
follows lung inflation; and (2) advancement of the onset time of
the next inspiratory burst, which follows seconds after lungs have
relaxed back to their resting volume. Changes in single-unit
activity accompanying these changes in respiratory pattern and
rhythm suggest constraints on the network organization of respiratory neurons.
Biphasic neurons
Biphasic neurons, typically silent in control cycles, show changes
in membrane potential in the peri-inspiratory interval consistent
with both excitatory and inhibitory drives. During lung inflation,
these neurons fired briskly and continued to fire after inflation
offset with a decrementing discharge pattern that varied from cell
to cell. In the presence of phasic inflation, these neurons resembled decrementing expiratory neurons in vivo, which are also
excited by lung inflation or vagal stimulation (Feldman et al.,
1976, Hayashi et al., 1996), or late inspiratory neurons with a
putative inspiratory off-switch function (Cohen et al., 1993). Bath
application of bicuculline blocks the BHE (Murakoshi and Otsuka, 1985), but biphasic neurons continue to fire in response to
lung inflation after bath application of bicuculline (our unpublished results). We hypothesize that these neurons provide
GABAergic inhibition of their postsynaptic targets (Fig. 9B) and
thus are essential constituents in the Breuer-Hering reflex
pathway.
Type 1 neurons
Type 1 neurons are hypothesized to be essential constituents of
the respiratory rhythm generator (Rekling et al., 1996; Rekling
and Feldman, 1998; Gray et al., 1999). Several observations made
here are consistent with this hypothesis: (1) they fired earliest in
relation to ventral root inspiratory onset and significantly earlier
than all other inspiratory neurons except type 2; (2) their burst
onset showed the least variability with respect to ventral root
inspiratory onset, both with and without phasic inflation; and (3)
Mellen and Feldman • Lung Afferents Transform Rat Respiratory Network In Vitro
J. Neurosci., September 15, 2001, 21(18):7363–7371 7369
Figure 8. Preinspiratory depolarization begins earlier with phasic inflation in pre-I neurons, and postinspiratory firing is lost. A, Raw traces of pre-I neuron
membrane potential with and without phasic inflation
are displayed as in Figure 3. In phasic inflation cycles,
the characteristic preinspiratory and postinspiratory
firing patterns of pre-I neurons in control cycles are
transformed into pure preinspiratory firing. B, Bursttriggered average of nine cycles from one pre-I neuron, displayed as in Figure 3. Bar and column graphs
show mean values for all six pre-I neurons and are
scaled to the axes of the averaged traces. Preinspiratory firing begins significantly earlier ( p ⬍ 0.05) during phasic inflation (bottom, black bar) than in control
cycles (bottom, white bar). Because no reversed outward currents are apparent after bias current application to hyperpolarize the cell to ⫺75 mV, loss of
postinspiratory activity is unlikely to be caused by
Cl ⫺-mediated inhibition.
Figure 9. A, Mean onset times of maximal firing in
type 1 (T-1), type 2 (T-2), type II (T-II ), and type III
(T-III ) neurons or offset of firing in pre-I neurons.
Onset times in control cycles (white columns) were not
significantly different from onset times with phasic
inflation ( gray columns). Statistically significant differences in onset times between cell types were found,
however, and are indicated by arrows between pairs of
cell types. These differences support inferences about
connectivity. B, Diagram of inferred connectivity between respiratory neurons. SAR afferent input (SAR)
causes brisk firing of biphasic neurons (biphasic) after
lung inflation. Postinspiratory inhibition seen only in
type 1 neurons (Type 1) is consistent with biphasic
neuron inhibition of type 1 neurons. Because type 1
inspiratory onset coincides closely with cervical ventral root motor output and fires earliest, type 1 neurons are hypothesized to drive inspiratory premotoneurons, such as type III (Type III ) and
motoneurons (MN ). Delay in onset, together with low
cycle-to-cycle variability in burst onset times in type II
neurons (Type II ), suggests an efference copy pathway
from motoneurons to type II neurons. Because onset
of type II neuron firing precedes but is not significantly different from onset of pre-I inhibition, type II neurons are hypothesized to inhibit pre-I neurons. Peri-inspiratory inhibition in type III neurons
(Type III ) is consistent with inhibition from pre-I neurons; the absence of this inhibition in type II neurons suggests that reciprocal inhibition from pre-I
neurons to type II is not present. Because type III firing onset is significantly later than pre-I inhibition, reciprocal inhibition from type III neurons to
pre-I neurons is not inferred. Neuron classes proposed to be causal to respiratory rhythm generation include pre-I neurons (Onimaru et al.; 1997), type
1 neurons (Rekling et al., 1996; Gray et al., 1999), and type 2 neurons (Thoby-Brisson et al.; 2000).
in the presence of phasic inflation, they showed consistent postinspiratory inhibition that coincided with biphasic neuron activity,
absent in all other respiratory neurons. The first observation is
consistent with the idea that if type 1 neurons generate respiratory rhythm, then they should be active before other inspiratory
neurons. Similarly, the observation of low variability in delay
between type 1 neurons and motor output burst onset is consistent with a causal relationship between their activation and the
onset of inspiratory motor activity. Finally, the inhibition
uniquely seen in type 1 neurons during phasic inflation provides
a mechanism for inspiratory burst shortening. When type 1 neurons are hyperpolarized below the firing threshold, excitatory
premotor drive to other inspiratory neurons is removed, leading
to early inspiratory burst termination (Fig. 9B). In addition, the
postinspiratory hyperpolarization seen in type 1 neurons during
phasic inflation could account for the observed increase in respiratory frequency via resetting of endogenous bursters among type
1 neurons (Rekling et al., 1996), as predicted by modeling studies
(Butera et al., 1999).
Type 2 neurons
The feature used here to differentiate between type 1 and 2
neurons is the presence of Ih and absence of IA in type 2 neurons.
The data here, albeit from a small sample, identified two additional features that distinguish type 2 from type 1 neurons: the
absence of immediate postinspiratory inflation-induced inhibi-
7370 J. Neurosci., September 15, 2001, 21(18):7363–7371
Mellen and Feldman • Lung Afferents Transform Rat Respiratory Network In Vitro
tion and delayed preinspiratory depolarization. The absence of
postinspiratory inhibition is not consistent with a role for type 2
neurons in inspiratory shortening. The earlier onset of preinspiratory depolarization of type 2 neurons with phasic inflation
may contribute to the shortened respiratory cycle. Because other
groups have identified Ih-positive inspiratory neurons with endogenous bursting properties in neonatal mice (Thoby-Brisson et al.,
2000), the criteria used here to identify type 2 neurons may select
electrophysiologically and functionally heterogeneous neurons.
Thus, more selective criteria are required to differentiate between functionally distinct classes of neurons sharing Ih.
apparent at either resting or hyperpolarized membrane potentials. In addition, lung inflation applied during midexpiration,
when pre-I neurons are silent, does not elicit hyperpolarization
(our unpublished results). Thus, loss of postinspiratory firing in
pre-I neurons is unlikely to be caused by hyperpolarization. In the
presence of phasic inflation, pre-I neurons in vitro resemble the
firing pattern of augmenting expiratory neurons in vivo (Hayashi
et al., 1996). In addition, the absence of immediate effects in
response to SAR activation in vitro matches the absence of
postsynaptic potentials in augmenting expiratory neurons after
vagal stimulation in vivo (Hayashi et al., 1996).
Type II neurons
Mechanisms for inflation-induced frequency increase
A consistent feature of the type II neuron firing pattern was a
steplike depolarization during firing onset (Fig. 6 B). This suggests that type II firing onset was tightly coupled to ventral root
inspiratory activity. Because neither bias currents below the estimated reversal potential for Cl ⫺-mediated inhibition in the
neonatal rat medulla (⫺75 mV; Shao and Feldman, 1997) (Fig.
6 B, bottom trace) nor Cl ⫺ loading (Onimaru et al., 1997) reveals
outward currents, preinspiratory inhibition cannot account for
the 100 msec delay between ventral root inspiratory onset and
type II firing onset. Although inspiratory burst duration was
shortened during phasic inflation, the burst-triggered averages
with and without phasic inflation otherwise matched. Type II
neuron burst onset occurs before pre-I inhibition onset but is not
significantly different. Thus type II neurons may inhibit pre-I
neurons, but because there is no evidence for peri-inspiratory
inhibition of type II neurons, pre-I inhibition of type II neurons
is unlikely (Fig. 9B).
Type III neurons
Because the preinspiratory hyperpolarization of type III neurons
is concurrent with the preinspiratory firing of pre-I neurons,
reciprocal inhibition between them has been proposed (Arata et
al., 1998). Consistent with this hypothesis, phasic inflation advances both the onset of preinspiratory firing in pre-I neurons and
the preinspiratory hyperpolarization in type III neurons. The
onset delay of depolarization in type III neurons is significantly
later than hyperpolarization onset in pre-I neurons (Fig. 9A), so
although pre-I neurons may inhibit type III neurons, type III
neurons are unlikely to be the source of the inspiratorymodulated inhibition of pre-I neurons (Fig. 9B).
Pre-I neurons
Although it has been proposed that pre-I neurons provide both
excitatory (Onimaru et al., 1992) and inhibitory drive to postsynaptic cells (Onimaru et al., 1997), here only effects consistent with
pre-I inhibitory drive were observed. As mentioned above, relative onset times suggest that type II neurons inhibit pre-I neurons,
which in turn inhibit type III neurons.
The transformation of pre-I activity during phasic inflation
cannot be accounted for by this hypothesized connectivity and
suggests that phasic inflation leads to a reorganization of network
inputs to pre-I neurons, modulation of pre-I intrinsic properties,
or both. Phasic inflation elicited both immediate cessation of
postinspiratory firing and advance in the onset of preinspiratory
firing effects in pre-I neurons, occurring seconds later. This pattern of activity, accompanied by an increase in respiratory frequency has also been obtained in vitro by bath application of
5-HT2A agonists (Onimaru et al., 1998). Although loss of postinspiratory firing is compatible with pre-I neuron inhibition during
and just after inflation, no postinspiratory hyperpolarization is
Inflation-induced inspiratory shortening matches the classic
Breuer-Hering reflex, whose effect is to shorten inspiration and
prolong expiration. Although it has been proposed that inspiratory shortening is required for shortening of subsequent expiration (Knox, 1973), lung inflation applied after inspiratory offset
can also significantly increase respiratory frequency (Mellen and
Feldman, 2000); thus, inspiratory shortening is not required for
shortening of the subsequent expiration. In addition, although
transient hyperpolarization of putatively rhythmogenic type 1
neurons might account for the observed increase in respiratory
frequency via resetting of endogenous burster neurons, this mechanism cannot account for the delayed increase in excitability seen
in type 2 and pre-I neurons, because type 1 neurons begin to fire
after these effects are seen.
In addition to the inspiration-terminating, expirationlengthening effect of vagal afferent feedback on central respiratory circuits, a parallel, more slowly adapting facilitatory pathway,
whose effects persist after offset of the inflation-induced inhibition, has been proposed (Younes and Polacheck, 1985). This
hypothesis is supported by the observation in vitro that midexpiratory inflation lengthens the period of the cycle in which it is
applied but also significantly shortens the subsequent period, well
after inflation offset (Mellen and Feldman, 1997). The observation that the effect of phasic inflation on respiratory frequency
decreases as inflation onset is delayed seems incompatible with
slowly adapting facilitation lasting tens of seconds. If the time
course of biphasic neuron firing reflects the duration of
expiration-lengthening GABAergic inhibition, then the phase
dependence of inflation-induced increases in respiratory frequency may arise out of the interaction between transient
GABAergic inhibition and tonic facilitation. Within this framework, inhibition resulting from transient inflation during inspiration or early in expiration will have decayed early enough in the
respiratory cycle to allow tonic facilitatory drive to advance the
onset of the subsequent inspiratory burst. Inhibitory drive resulting from transient inflation later in the respiratory cycle will mask
facilitatory drive, eliminating the frequency effect.
Afferent modulation of central rhythmogenic circuits
In the study of central pattern generators (CPGs), the dynamic
interplay between sensory feedback and central rhythmogenic
circuits has been mostly neglected because of the highly reduced
preparations used (Wilson, 1961, Cohen and Wallén, 1980).
When afferent feedback to CPGs was restored in invertebrate
systems, timing relationships between CPG constituents were
transformed, indicating that in the intact organism, phasic afferent feedback reorganizes rhythmogenic networks (Wolf and
Pearson, 1987). In the lung-attached preparation used here, it has
been possible to approximate in vivo feedback conditions using
Mellen and Feldman • Lung Afferents Transform Rat Respiratory Network In Vitro
transient lung inflations triggered off inspiratory motor output.
The transformation of firing patterns in subsets of respiratory
neurons suggests that in mammals also, proprioceptive afferent
feedback reorganizes central rhythmogenic circuits. Thus, the
transformation of primarily silent biphasic neurons into briskly
firing early expiratory neurons, complemented by the inhibition
observed in type 1 neurons, suggests that the postulated neural
substrate for an inspiratory off switch, critical to in vivo models of
respiratory rhythm generation (Richter, 1982; von Euler, 1983;
Sammon, 1994), is latent in vitro. In addition, the transformation
of pre-I neurons from peri-inspiratory to late expiratory suggests
mappings from in vivo to in vitro classification schemes. Finally,
these transformations taken together suggest that sensory feedback may play a critical role in generating the putative threephase respiratory cycle in vivo (Richter, 1982).
REFERENCES
Arata A, Onimaru H, Homma I (1998) Possible synaptic connections of
expiratory neurons in the medulla of newborn rat in vitro. NeuroReport 9:743–746.
Blanton MG, Lo Turco JJ, Kriegstein AR (1989) Whole cell recording
from neurons in slices of reptilian and mammalian cerebral cortex.
J Neurosci Methods 30:203–210.
Bonham AC, Coles SK, McCrimmon DR (1993) Pulmonary stretch receptor afferents activate excitatory amino acid receptors in the nucleus
tractus solitarii in rats. J Physiol (Lond) 464:725–745.
Breuer J (1868) Self-steering of respiration through the nervus vagus
[English translation]. KK Akad Wissenschaft Wien Sitzungsberichte
Mathematisch Naturwissenschaft Classe 67:909 –937.
Bruce EN (1997) Chemoreflex and vagal afferent mechanisms enhance
breath to breath variability of breathing. Respir Physiol 110:237–244.
Butera RJ, Rinzel J, Smith JC (1999) Models of respiratory rhythm
generation in the pre-Bötzinger complex. II. Populations of coupled
pacemaker neurons. J Neurophysiol 82:398 – 415.
Cohen AH, Wallén P (1980) The neuronal correlate of locomotion in
fish. “Fictive swimming” induced in an in vitro preparation of the
lamprey spinal cord. Exp Brain Res 41:11–18.
Cohen MI, Huang WX, Barnhardt R, See WR (1993) Timing of medullary late-inspiratory neuron discharges: vagal afferent effects indicate
possible off-switch function. J Neurophysiol 69:1784 –1787.
Connor JA, Stevens CF (1971) Voltage clamp studies of a transient
out-ward membrane current in gastropod neural somata. J Physiol
(Lond) 213:21–30.
Feldman JL, Gautier H (1976) Interaction of pulmonary afferents and
pneumotaxic center in control of respiratory pattern in cats. J Neurophysiol 39:31– 44.
Feldman JL, Cohen MI, Wolotsky P (1976) Powerful inhibition of pontine respiratory neurons by pulmonary afferent activity. Brain Res
104:341– 634.
Gray PA, Rekling JC, Bocchiaro CM, Feldman JL (1999) Modulation of
respiratory frequency by peptidergic input to rhythmogenic neurons in
the preBötzinger complex. Science 286:1566 –1568.
Hayashi F, Coles SK, McCrimmon DR (1996) Respiratory neurons mediating the Breuer-Hering reflex prolongation of expiration in rat.
J Neurosci 16:6526 – 6536.
Knox C (1973) Characteristics of inflation and deflation reflexes during
expiration of the cat. J Neurophysiol 36:284 –295.
J. Neurosci., September 15, 2001, 21(18):7363–7371 7371
Mellen NM, Feldman JL (1997) Vagal stimulation induces expiratory
lengthening in the in vitro neonate rat. J Appl Physiol 83:1607–1611.
Mellen NM, Feldman JL (2000) Phasic lung inflation shortens inspiration and respiratory period in the lung-attached neonate rat brain stem
spinal cord. J Neurophysiol 83:3165–3168.
Murakoshi T, Otsuka M (1985) Respiratory reflexes in an isolated
brainstem-lung preparation of the newborn rat: possible involvement of
gamma-aminobutyric acid and glycine. Neurosci Lett 62:63– 68.
Onimaru H, Homma I, Iwatsuki K (1992) Excitation of inspiratory neurons by preinspiratory neurons in rat medulla in vitro. Brain Res Bull
29:879 – 882.
Onimaru H, Arata A, Homma I (1997) Neuronal mechanisms of respiratory rhythm generation: an approach using in vitro preparation. Jpn
J Physiol 47:385– 403.
Onimaru H, Shamoto A, Homma I (1998) Modulation of respiratory
rhythm by 5-HT in the brainstem-spinal cord preparation from newborn rat. Eur J Physiol 435:485– 494.
Pape HC, McCormick DA (1989) Noradrenaline and serotonin selectively modulate thalamic burst firing by enhancing a hyperpolarizationactivated cation current. Nature 340:715–718.
Rekling JC, Feldman JL (1998) Pre-Bötzinger complex and pacemaker
neurons: hypothesized site and kernel for respiratory rhythm generation. Annu Rev Physiol 60:385– 405.
Rekling JC, Champagnat J, Denavit-Saubié M (1996) Electroresponsive
properties and membrane potential trajectories of three types of inspiratory neurons in the newborn mouse brainstem in vitro. J Neurophysiol 75:795– 810.
Richter DW (1982) Generation and maintenance of the respiratory
rhythm. J Exp Biol 100:93–107.
Sammon M (1994) Geometry of respiratory phase switching. J Appl
Physiol 77:2468 –2480.
Shao XM, Feldman JL (1997) Respiratory rhythm generation and synaptic inhibition of expiratory neurons in pre-Bötzinger complex: differential roles of glycinergic and GABAergic neural transmission.
J Neurophysiol 77:1853–1860.
Smith JC, Feldman JL (1987) In vitro brainstem-spinal cord preparations for study of motor systems for mammalian respiration and locomotion. J Neurosci Methods 21:321–333.
Smith JC, Greer JJ, Liu GS, Feldman JL (1990) Neural mechanisms
generating respiratory pattern in mammalian brain stem-spinal cord in
vitro. I. Spatiotemporal patterns of motor and medullary neuron activity. J Neurophysiol 64:1149 –1169.
Smith JC, Ellenberger HH, Ballanyi K, Richter DW, Feldman JL (1991)
Pre-Bötzinger complex: a brainstem region that may generate respiratory rhythm in mammals. Science 254:726 –729.
Suzue T (1984) Respiratory rhythm generation in the in vitro brain
stem-spinal cord preparation of the neonatal rat. J Physiol (Lond)
354:173–183.
Thoby-Brisson M, Telgkamp P, Ramirez J-M (2000) The role of
hyperpolarization-activated current in modulating rhythmic activity in
the isolated respiratory network of mice. J Neurosci 20:2994 –3005.
von Euler C (1983) On the central pattern generator for the basic breathing rhythmicity. J Appl Physiol 55:1647–1659.
Widdicombe J (1961) Respiratory reflexes in man and other mammalian
species. Clin Sci 21:163–170.
Wilson D (1961) The central nervous control of flight in the locust. J
Exp Biol 38:471– 490.
Wolf H, Pearson KG (1987) Comparison of motor patterns in the intact
and deafferented flight system of the locust. J Comp Physiol [A]
160:269 –279.
Younes M, Polacheck J (1985) Central adaptation to inspiratoryinhibiting expiratory-prolonging vagal input. J Appl Physiol 59:1072–
1084.