Abstract
Synchrony and oscillations in neuronal firing play important roles in information processing in the mammalian brain. Here, we evaluate their role in controlling neuronal output in a well defined motor behavior, breathing, using an in vitro preparation from neonatal rat that generates respiratory-related motor output. In this preparation, phrenic motoneurons (PMNs) receive endogenous rhythmic inspiratory currents with prominent oscillations in the 20-50 Hz range. We recorded these inspiratory currents in individual PMNs and used them as test inputs for the same motoneuron (MN) during the normally silent expiratory periods. The impact of the oscillations on MN output was evaluated by filtering the currents before injection. Responses to unfiltered inspiratory currents were indistinguishable from voltage changes during spontaneous inspiratory periods. More than 90% of action potentials occurred within milliseconds [-2 to +4] of the oscillation peaks. The timing of action potentials was highly reproducible in response to unfiltered currents. Attenuation of the oscillations by low-pass filtering (<50 Hz) decreased the precision in action potential timing and significantly reduced the number of action potentials by ∼35%. The adrenergic agonist phenylephrine increased instantaneous firing frequency in responses evoked by square-wave or low-pass filtered inspiratory currents but had no effect on firing frequency evoked by unfiltered currents. We conclude that oscillations control the precise timing of action potentials, help to maximize synaptic drive efficiency, and constrain MN firing frequencies to those optimal for muscle contraction.
Introduction
Determining how synaptic inputs are transformed into patterns of action potential output is key to understanding information processing in the CNS. The role of intrinsic membrane properties in this transformation is illuminated by approaches that use tonic afferent stimulation or intracellular injection of square-wave or ramp current waveforms to activate neurons. The dynamic characteristics of endogenous synaptic inputs to neurons that are important in transforming input into output are less well characterized. A prominent dynamic component of endogenous inputs is high-frequency (>10 Hz) oscillations. These oscillations can be seen throughout the CNS (Gray, 1994; Fetz et al., 2000), including in networks that provide rhythmic excitatory drive to motoneurons (MNs) in behaviors such as chewing (Smith and Denny, 1990), speaking (Smith and Denny, 1990; Nakazawa et al., 2000), and breathing (Cohen and Gootman, 1969; Cohen et al., 1987; Liu et al., 1990; Smith and Denny, 1990; Christakos et al., 1991; Tarasiuk and Sica, 1997; Funk and Parkis, 2002). Oscillations with peak power in the 20-150 Hz range are ubiquitous in the respiratory network. They are synchronized between different network elements (Liu et al., 1990; Cohen et al., 1997; Tarasiuk and Sica, 1997), yet their function is not known. Here, we explore the importance of endogenous oscillations in synaptic input for MN information processing. We recorded endogenous inspiratory currents in phrenic MNs (PMNs) in a neonatal rat brainstem spinal cord preparation that spontaneously generates rhythmic, respiratory-related motor nerve output (Smith and Feldman, 1987). These endogenous currents were reinjected as command inputs, before and after low-pass filtering, into the same PMN during the normally quiescent expiratory period, and the voltage response was recorded. We also compared the effects of phenylephrine (PE), a potent neuromodulator agonist, on the voltage response of PMNs to these filtered and unfiltered inspiratory currents, as well as to traditional square-wave currents. Analysis of MNs offers the advantage that, as the final common output path from the CNS to skeletal muscles, the physiological significance of changes in their input-output relationships is readily interpretable.
Materials and Methods
Preparation. Brainstem-spinal cords were isolated from neonatal Wistar rats ranging in age from postnatal days 0-3 as described previously (Smith and Feldman, 1987; Parkis et al., 2000). Preparations were pinned down, ventral surface up, on Sylgard resin in a recording chamber (2 ml volume) and perfused at a flow rate of 2-4 ml/min with oxygenated artificial CSF (in mm): 120 NaCl, 3 KCl, 1 CaCl2, 2 MgSO2, 26 NaHCO3, 1.25 NaH2PO4, and 20 d-glucose, pH 7.45 (oxygenated with 95% O2-5% CO2 at 26-27°C).
Recording. Population inspiratory activity was recorded from the severed ends of C1 or C4 nerve roots using suction electrodes, amplified, filtered (0.3-3 kHz), rectified, and integrated (Smith and Feldman, 1987; Parkis et al., 2000). Whole-cell voltage- and current-clamp data were recorded from PMNs using an Axopatch 1D amplifier (Axon Instruments, Foster City, CA) controlled by Axodata version 1.2.2 (Axon Instruments) and custom-written programs developed in LabVIEW (Parkis et al., 2000). Intracellular recording electrodes (3.5-4.5 Ω) contained the following (in mm): 125 K +-gluconate, 5 NaCl, 1 CaCl2, 10 HEPES, 10 BAPTA, and 2 Mg 2+-ATP, pH 7.2-7.3. All signals were recorded on videotape via PCM digital processor (Vetter 3000A; Vetter Instruments, Reberburg, PA).
Protocols. A software program developed in our laboratory using the LabVIEW (National Instruments, Austin, TX) programming environment was used to activate neurons with endogenous waveforms of synaptic current (Parkis et al., 2000). Voltage-clamp recordings of individual inspiratory synaptic currents were triggered from bursts of inspiratory activity on C1 or C4 ventral nerve, and 5-10 different synaptic waveforms were acquired (20 kHz) to a waveform library specific for that cell. The recording configuration was then switched to current clamp, and a waveform was selected to be used as a test input. Because of the negative sign of inward currents measured in voltage clamp, the recorded synaptic waveform was inverted by the LabVIEW program so that injected inspiratory currents would be depolarizing. Each motoneuron was activated with one of its own recorded currents.
When currents were filtered before injection, a low-pass Butterworth filter (order of 3-5) was used. Different waveforms (representing the same initial recorded current filtered using different cutoff frequencies) were presented a minimum of four times in random order, and the number of spikes produced by each stimulus was reported as an average of four responses. Unfiltered waveforms were presented twice as often as other stimuli to ensure consistent control responses. Current injections were triggered to occur in the quiescent (expiratory) period between inspiratory bursts.
To examine the effects of α1 noradrenergic receptor activation on responses of PMNs to injection of smooth and synaptic current waveforms, the α1 receptor agonist PE was locally applied (5-10 min) via a pressure injection pipette over the phrenic nucleus (Parkis et al., 1999). Control and PE responses were elicited from the same membrane potential by injecting DC current to correct for the PE-induced depolarization. To control for time- and activity-dependent changes in firing output that differ between different waveforms, analysis of all responses was limited to the time range over which PMNs fired in response to stimulation with unfiltered waveforms (see Fig. 3a,b, dashed boxes). The magnitude of the square-wave current pulse was selected to produce an average instantaneous firing rate similar to that produced by endogenous unfiltered synaptic currents in control. Maximum effects of PE on C4 nerve inspiratory activity for each animal were taken to be the maximum value in a moving average of five consecutive bursts within 2 min of PE application.
PE induced increases in instantaneous firing frequency responses to smooth, but not endogenous, input waveforms. Membrane potential responses of a PMN to injection of an inspiratory synaptic waveform (A) and a square-wave current pulse (B) in the absence (control) and presence of PE (100 μm). C, D, Instantaneous firing frequencies for the responses evoked in A and B are plotted versus time in C and D, respectively. E, Effect of PE on the response of a PMN to stimulation with a 10 Hz low-pass filtered inspiratory synaptic waveform. F, Effects of PE on instantaneous firing frequency evoked through unfiltered synaptic currents (USC) and smooth (square-wave or 10 Hz filtered) waveforms (n = 6) (dashed boxes in A and B indicate period over which firing frequency was analyzed). Asterisk denotes significant difference in firing frequency relative to control conditions (p = 0.008) and that magnitude of the increase in response to smoothed inputs was greater than that in response to USC (p = 0.01).
Analysis. Offline analysis was performed using AxoGraph (version 4.3) and Microsoft (Seattle, WA) Excel (version 5.0) software. Power spectra of captured currents were generated offline using LabVIEW software. Differences between means were compared using one-way ANOVA (SAS 6.1; SAS Institute, Cary, NC). Linear contrast coefficients were used to partition the sum of squares to determine significant differences between groups. Values of p < 0.05 were considered significant. Data are given as mean ± SE.
Results
Characteristics of inspiratory currents
Experiments were performed in a neonatal rat brainstem-spinal cord preparation (Smith and Feldman, 1987) that generates endogenous, periodic inspiratory-related synaptic inputs to PMNs (Liu et al., 1990; Rekling et al., 2000a). These currents were synchronous with bursts of inspiratory activity on the C4 spinal cord ventral nerve root, had peak amplitudes of -803 ± 55 pA (n = 35), and exhibited a rapidly increasing (140 ± 10 msec to peak) slowly decrementing (530 ± 30 msec) waveform.
Power spectral analysis of the inspiratory drive currents to PMNs (n = 22 neurons; 5-17 consecutive cycles per neuron) revealed two main patterns of synaptic input. Six PMNs received synaptic currents that contained only small-amplitude peaks in the 20-50 Hz range (∼0.04% of total power), centered at 32 ± 1.7 Hz (Fig. 1A, left). Sixteen PMNs received inspiratory synaptic inputs that contained significant peaks at 32 ± 1.4 Hz with mean peak power more than threefold higher (0.13 ± 0.01% of total power) (Fig. 1A, right). These MNs were otherwise indistinguishable on the basis of input resistance, rheobase, membrane capacitance, or magnitude of synaptic current.
Action potentials arise preferentially on peaks of high-frequency oscillations. A, Power spectra of four consecutive inspiratory synaptic currents from two different PMNs, receiving inspiratory drive currents with different amounts of oscillatory activity. B, Stimulation of a PMN with an inspiratory synaptic current waveform (bottom trace) elicited membrane depolarization and eight action potentials. Dashed lines highlight close correlation between peaks of oscillations in the injected current and occurrence of action potentials. Asterisks in B (left) denote action potentials shown in expanded time scale to the right. C, Histogram showing temporal distribution of 689 action potentials recorded in 13 PMNs relative to nearest peak in the injected inspiratory waveform. Time 0 indicates an oscillation peak. Dashed box indicates the 6 msec interval during which 92% of action potentials occurred.
Oscillations control action potential timing
To reproduce as closely as possible the dynamic pattern of membrane depolarization that occurs during actual respiratory behavior, we used our method for repeatedly stimulating each PMN with an inspiratory synaptic current waveform recorded from that MN (Parkis et al., 2000). These endogenous inspiratory current waveforms were injected during the quiescent period between endogenous inspiratory bursts, i.e., during the expiratory period. The membrane depolarization evoked by spontaneous inspiratory inputs, averaged from a representative sample of PMNs that received consistent inspiratory inputs (four consecutive spontaneous inspiratory cycles from four different PMNs (26.5 ± 3.5 mV), was the same as that evoked in the same MNs by reinjected synaptic currents (29.2 ± 4.7 mV).
Consistent with responses of hypoglossal MNs to simulated postsynaptic currents of varying amplitude (Türker and Powers, 1999), the timing of action potentials and the occurrence of oscillations in the endogenous current waveform were closely correlated (Fig. 1B). Most action potentials (92%; 631 of 689; 124 current injections in 13 PMNs) occurred within a 6 msec window around the oscillation peaks (-2 to +4 msec) in which the average interpeak interval was ∼32 msec (Fig. 1C).
We then tested the hypothesis, on the basis of responses of cortical neurons to random white noise inputs (Mainen and Sejnowski, 1995), that synaptically generated oscillations increase the temporal precision, i.e., reproducibility of action potential occurrence during repetitive firing. PMNs were stimulated repeatedly with the same synaptic current waveforms before and after low-pass filtering at 10 Hz. Repeated stimulation with unfiltered currents produced highly consistent spike trains (Fig. 2). In six PMNs tested in this manner, the SDs of the delay between first and second, and first and last action potentials averaged across all cells were 0.67 ± 0.20 and 1.80 ± 0.54 msec (n = 6) (total burst duration of 190 ± 24 msec), respectively. After filtering of the injected current, although significant clustering of action potentials was still apparent (reflecting the contribution of intrinsic membrane properties to this phenomenon), there was a significant reduction in the precision of action potential timing. The second action potential of each train remained tightly locked to the first (SD, 0.96 ± 0.25 msec; n = 6). However, the timing of the last spike was highly variable (SD, 10.08 ± 4.49 msec; n = 6). The increased variability was also evident in the coefficient of variation for the delay between first and last spikes, which was sevenfold greater for responses to filtered (0.063 ± 0.025) compared with unfiltered (0.009 ± 0.002) waveforms. The coefficient of variation for the last interspike interval of responses to filtered currents (0.118 ± 0.037; n = 6) was also significantly increased over that for unfiltered waveforms (0.022 ± 0.009; n = 6). These data suggest that naturally occurring oscillations in synaptic drive play a critical role in determining the timing of action potentials during repetitive firing.
Oscillations on inspiratory synaptic currents increase precision of action potential timing during repetitive firing. Membrane potential responses (top traces) of a PMN to repeated application of an inspiratory synaptic current waveform (Iinj) in its unfiltered (A) and filtered (B; 10 Hz low pass) form. Membrane potential responses are shown superimposed (first 4 trials). C, D, Raster plots of spike times relative to the onset of the first action potential in the train (5 consecutive trials). E, F, Raster plots showing spike times relative to the onset of the first action potential in five consecutive endogenous inspiratory cycles. The plot in E represents activity of a PMN driven by inputs with consistent oscillations. Synaptic inputs to the MN in F lacked oscillations.
One possible function for these oscillations is to ensure optimal motor unit activation by producing a precise pattern of output in each motoneuron. If so, one would predict that the inspiratory discharge patterns of MNs receiving synaptic inputs with strong oscillations would be more consistent than those without. To test this possibility, we measured the variability in spike timing of bursts generated by consecutive spontaneous currents in four PMNs receiving inputs with prominent oscillations and four different PMNs receiving inputs that lacked prominent oscillations. Compared with PMNs receiving inputs that lacked oscillations (Fig. 2F), discharge patterns of PMNs receiving inputs with consistent oscillations showed much less variability for consecutive spontaneous bursts (Fig. 2E). The average coefficient of variation for the interval between the first and seventh spikes was twofold higher in the MNs driven by inputs that lacked oscillations (0.16 ± 0.10) compared with those with strong oscillations (0.08 ± 0.01). Similarly, the coefficients of variation for the last interspike interval for MNs with and without strong oscillations were 0.27 ± 0.11 and 0.41 ± 0.07, respectively. Thus, although variations in the pattern of synaptic input to individual MNs mean that the pattern of action potential output will not be identical between consecutive cycles, oscillations in the input reduce variability in interspike interval as predicted, and play an important role in controlling spike timing.
Oscillations constrain mechanisms by which modulators alter excitability
Neuromodulators such as norepinephrine, substance P, and thyrotropin-releasing hormone increase MN excitability (Rekling et al., 2000a; Powers and Binder, 2001) by reducing a passive resting K+ leak conductance and by inhibiting active conductances, such as the Ca2+-activated K+ conductance. When the impact of these modulators on MN firing behavior is tested using square-wave inputs, these effects produce a leftward shift and an increase in the slope of the relationship between firing frequency and current. The net result for MN output is higher discharge frequencies and, for nonsquare-wave inputs, earlier onset of firing. In other words, one would expect an excitatory modulator to significantly alter the frequency and number of action potentials evoked by a given input waveform.
If, however, oscillations preserve action potential timing under conditions in which excitability is markedly changed, we would expect neuromodulators to have less impact on firing frequency when MNs receive inputs with oscillations. We tested this by comparing the effects of PE on instantaneous firing frequency produced by reinjected, unfiltered, inspiratory waveforms with that produced by smooth waveforms (filtered inspiratory waveforms and square-wave forms). We also assessed the effects of PE on the firing frequency evoked in PMNs by spontaneous inspiratory inputs.
Consistent with its effects on other MNs (Parkis et al., 1995; Binder et al., 1996; Rekling et al., 2000a), PE depolarized PMNs (6.0 ± 1.0 mV; n = 6), decreased Irheo (from 540 ± 96 to 330 ± 80 pA; n = 5), and increased input resistance (Rn) (7-41%). In response to square-wave (Fig. 3B, D) and filtered (Fig. 3E) currents, PE significantly increased PMN instantaneous firing frequency from 46.0 ± 2.1 to 58.0 ± 3.8 Hz (n = 6) (Fig. 3F). In contrast, when the same PMNs were activated with unfiltered inspiratory synaptic currents, PE did not affect instantaneous firing frequency (Fig. 3A,C,F). Instead, its main effect was to double the duration over which firing occurred, from 110 ± 30 to 220 ± 40 msec (n = 6). This can be seen for one PMN in Figure 3A, in which the six action potentials evoked under control conditions occurred at similar time points in the presence of PE, i.e., they arose off the same peaks in the injected current. In PE, however, action potentials seen in control conditions were preceded and followed by additional action potentials that arose from peaks in the input waveform that previously produced subthreshold voltage deflections.
Similar behavior was apparent in response to spontaneous inspiratory inputs. PE had no effect on instantaneous firing frequency evoked by spontaneous inspiratory currents, regardless of whether or not the depolarization induced by PE was offset with DC current injection (control, 43.4 ± 4.8 Hz, n = 6; PE with DC offset, 43.8 ± 5.6 Hz, n = 6; PE without DC offset, 41.2 ± 0.6 Hz, n = 3). The duration of PMN discharge increased from 140 ± 35 msec in control to 250 ± 71 msec in PE (n = 6), which manifested as a significant increase in the duration (57 ± 22%) of inspiratory bursts recorded from the population of PMN axons contained in the C4 nerve root (Fig. 4). Because PE did not increase firing frequency, the PE-induced increase in integrated C4 nerve inspiratory burst amplitude (45 ± 7%) was likely attributable to increased duration of action potential firing and/or recruitment of additional PMNs.
PE increases both the duration and amplitude of inspiratory bursts recorded from the C4 ventral nerve root (which contains the phrenic nerve) but not firing frequency of PMNs during inspiration. A, Inspiratory burst envelopes averaged from 10 consecutive inspiratory cycles (∫C4) before and during application of PE (100 μm). B, Maximum effects of PE on C4 nerve burst duration (∫C4 dur) and amplitude (∫C4 amp) relative to control (control values, 1.0). Effects of PE on the average duration over which individual PMNs produced action potential output (PMN dur) in response to endogenous input is included for comparison. Note that the instantaneous firing frequency (PMN freq) evoked by endogenous inputs was not affected by PE.
Oscillations increase input-output efficiency
To further explore the importance of endogenous synaptically generated current oscillations in controlling repetitive firing behavior, we compared responses to unfiltered synaptic current waveforms and the same waveforms after low-pass filtering with cut-off frequencies between 10 and 100 Hz. Low-pass filtering at 50 Hz and below reduced the number of action potentials (Fig. 5). Filtering at 50, 20-25, and 10 Hz reduced the number of action potentials elicited by injected currents to 91 ± 4, 74 ± 5, and 65 ± 5%, respectively, of the response to the unfiltered current (Fig. 5C) (n = 12). In the example of Figure 5, the number of action potentials produced was reduced from six in control (unfiltered), to five with a 50 Hz cutoff, and to four with 25 and 10 Hz cutoff. This reduction in output corresponded to removal of a large peak in the power spectrum centered at ∼37 Hz (Fig. 5B). Integration of the areas under the injected current waveforms verified that filtering did not reduce the charge transferred by the injected current at any cut-off frequency (total charge transfer was reduced by <0.1% at the lowest cut-off frequency of 10 Hz). Reductions in total power by 2.1 ± 0.7 and 3.9 ± 1.4% with low-pass filtering at 20 and 10 Hz, respectively (n = 10), were also insignificant.
Removal of oscillations decreases efficacy of currents. A, Membrane potential responses of a PMN (top) to the same synaptic current waveform (bottom) injected before (unfiltered) and after low-pass filtering with cut-off frequencies set at 50, 25, and 10 Hz (action potentials are truncated). B, Power spectra corresponding to the waveforms shown in A illustrate the effects of filtering on the power-frequency relationship. C, The number of action potentials evoked by stimulation of PMNs (n = 12) with synaptic current waveforms filtered at a variety of low-pass cut-off frequencies compared relative to the number of action potentials evoked by the unfiltered waveform (dashed line at 100%).
Discussion
Synchronized oscillations in neuronal activity with power in the range of 10-150 Hz appear throughout the brain (Gray, 1994) and are postulated to contribute to a variety of higher functions (Konig and Engel, 1995; Engel et al., 1999), including memory formation, control of arousal state, and “binding” of neuronal arrays for sensory perception. Oscillations in this range are also prominent features of activity in the system controlling breathing. They are present in the inspiratory activities of medullary neurons underlying respiratory rhythm and pattern generation and in respiratory motoneurons, nerves, and muscles in all mammals studied (Cohen et al., 1997; Funk and Parkis, 2002). However, the function of inspiratory oscillations is not known. Here, we provide evidence for two important consequences of these oscillations for MN behavior: they control the timing of action potentials and increase input-output efficiency.
Critique of method
Simulated current waveforms have proven useful for exploring the relationship between current transients and firing probability and highlight the importance for firing behavior and action potential timing of dynamic changes in membrane potential (Mainen and Sejnowski, 1995; Nowak et al., 1997; Tang et al., 1997; Stevens and Zador, 1998; Volgushev et al., 1998; Rekling et al., 2000b; Beierholm et al., 2001; Powers and Binder, 2001), particularly the rise time and amplitude of EPSPs (Türker and Powers, 1999, 2002). To determine whether these effects are seen under more natural conditions of synaptic input, we reinjected endogenous inspiratory currents recorded from the same MNs. These somatic currents are relevant for functional investigations because they result from the membrane integration of behaviorally meaningful endogenous inputs arriving at synapses distributed over the somatodendritic membrane.
An important technical consideration is that currents recorded under voltage clamp may be greater than currents that would normally reach the soma. Under physiological conditions, i.e., no voltage clamp, synaptically evoked membrane depolarization reduces the driving force and magnitude of the non-NMDA receptor-mediated component of the synaptic current. The magnitude of this depolarization-mediated reduction in the non-NMDA component of the inspiratory current, however, will be partially offset by increased NMDA-mediated current (Liu and Feldman, 1992; Cook and Johnston, 1999). Increased membrane conductance associated with activation of glutamatergic inspiratory synapses (Liu et al., 1990) could also attenuate current transfer to the soma under physiological conditions. Reduction of this attenuation under voltage clamp could result in our recording a larger than normal synaptic current. In addition, when currents are reinjected during expiration, membrane conductance may be lower than during inspiration because of the lack of glutamatergic input. If any of these factors played a significant role in our experiments, reinjected currents would have produced larger membrane depolarizations than spontaneous inspiratory inputs. This did not appear to be the case. As discussed previously (Parkis et al., 2000) and in a representative sample of PMNs in the present study, the membrane depolarization evoked by reinjected synaptic currents was not significantly greater (10%) than that evoked by spontaneous synaptic inputs. Thus, inspiratory currents recorded at the soma are similar to those reaching the soma under normal conditions.
Functional significance
Endogenous oscillatory inputs increase the efficiency with which MNs transduce input into output. This is consistent with responses of neocortical neurons to fluctuating inputs simulated by adding Gaussian noise to square-wave pulses (Tang et al., 1997). Thus, by driving PMNs with oscillatory synaptic currents, neuronal synchrony within respiratory motor networks can achieve the same output with less input. The biophysical basis of this efficiency gain extends beyond the benefits of simply summing coincident EPSCs to briefly drive the membrane above threshold. Action potentials coincided with oscillation peaks even when the underlying envelope of synaptic drive maintained membrane potential continuously above the repetitive firing threshold (as determined by square-wave inputs) (Fig. 3, compare A, B). Because action potential threshold decreases (Schlue et al., 1974; Dai et al., 2000) and firing probability increases (Türker and Powers, 1999) with increased rate of depolarization in MNs, coincidence of synaptic inputs functionally amplifies their effect by causing more rapid membrane depolarization (Azouz and Gray, 2000).
From a functional perspective, i.e., motoneuronal discharge producing muscle contraction, oscillatory inputs to MNs mean that generation of a given muscle force will require fewer action potentials from the premotoneuron pool. Similar efficiency gains throughout the respiratory network would reduce the metabolic cost of respiratory rhythm generation and drive transmission (i.e., fewer presynaptic action potentials will be required to evoke a postsynaptic action potential).
Consistent with experiments using simulated current transients (Mainen and Sejnowski, 1995) or oscillatory inputs (Tang et al., 1997; Volgushev et al., 1998), data in cortical neurons (Nowak et al., 1997) and our data in MNs demonstrate that endogenous oscillations in synaptic input strongly influence the timing of action potentials. Action potential timing may be profoundly important within higher neuronal networks by providing an additional coding dimension for information processing and the binding together of spatially distributed neuronal populations. At the level of MN and muscle behavior, the functional consequences of oscillations and their influence on action potential timing are threefold.
Oscillations constrain the mechanisms by which neuromodulators affect MN activity. Potentiation of instantaneous firing frequency is a major effect of neuromodulators on MN responses to steady-state inputs. However, rather than increasing firing frequency of individual PMNs, which would have limited value when MNs are already firing at frequencies close to the fusion frequency of the muscle fibers they innervate (at and above fusion frequency, additional increases in MN firing frequency produce relatively smaller increments in force), the impact of neuromodulators on MNs driven by oscillating inputs will be achieved through recruitment of quiescent MNs and prolongation of firing duration. This conclusion is consistent with the observed PE-mediated potentiation of C4 (Fig. 4) and hypoglossal nerve inspiratory output (Selvaratnam et al., 1998). Modulators could increase discharge frequency during the active period (Lalley, 1986) if subthreshold oscillations interspersed between suprathreshold oscillations became suprathreshold (Tang et al., 1997) or if multiple rather than single spikes arose from each oscillatory peak. However, these mechanisms were not significant contributors under our experimental conditions. Firing frequency did not increase significantly after PE in any PMN driven by inspiratory waveforms or spontaneous inspiratory inputs. That these principles govern input-output processing of PMNs in general is supported by the observation that substance P, like PE, increases firing frequency of PMNs in response to square-wave but not inspiratory synaptic currents (Ptak et al., 2000). Thus, oscillations may serve to maintain MN firing frequency within an optimal range, reducing the variability in discharge frequency that might otherwise occur with the fluctuations in neuromodulator levels, and MN excitability, which accompany changes in activity (transitions from rest to exercise), stress (Stanford, 1995), and arousal state (Aston-Jones and Bloom, 1981).
Implicit in our data is that a change in MN firing frequency requires a change in the oscillation frequency. Indeed, in contrast to neuromodulatory inputs, some mechanosensory and chemo-sensory inputs are associated with small but significant increases in both oscillation frequency and firing frequency (Funk and Parkis, 2002).
Oscillations impact the development of muscle force. The pattern of MN output, not simply the average frequency, determines muscle force (van Lunteren and Sankey, 2000). Thus, although efficient muscle activation is primarily attributed to the matching of MN firing properties and the contractile properties of their innervated muscle fibers, respiratory oscillations must also contribute. By controlling action potential timing, oscillations may increase the efficiency of diaphragmatic contraction by ensuring that muscle fibers are activated within an optimal range of frequencies. Oscillations can also prevent excessive PMN activation and the disastrous consequences of phrenic motoneuronal accommodation or diaphragmatic fatigue causing respiratory failure, by limiting the maximal instantaneous firing frequency. Perhaps for this reason (a) peaks in the power spectra of inspiratory inputs in neonatal rat PMNs (Liu et al., 1990; Parkis et al., 1998) are similar to the fusion frequency of neonatal diaphragm (Martin-Caraballo et al., 2000), (b) the rate of diaphragmatic fatigue development increases significantly with motor unit stimulation frequencies >60 Hz (Martin-Caraballo et al., 2000).
Control of action potential timing by the oscillations also implies that, for MNs receiving consistent oscillatory inputs, increases in muscle force are primarily achieved through recruitment of quiescent MNs. This contrasts with the convention that increases in firing frequency and recruitment both contribute.
PMNs fire synchronously, as evidenced by common peaks in the power spectra of population (Smith and Denny, 1990) and unit activities, as well as peaks in the coherence spectra between phrenic nerve and PMN activities (Christakos et al., 1991). Short-term synchrony, the tendency for some motor units to discharge within a few milliseconds of each other, occurs in many skeletal muscles (Sears and Stagg, 1976; Kirkwood et al., 1982; Datta, 1990; Farmer, 1998; Baker et al., 1999) and has been attributed to common input from synchronization of discharge in presynaptic fibers. Common inputs can also synchronize MN output (Binder and Powers, 2001; Türker and Powers, 2001, 2002), especially those composed of relatively few, large EPSPs (Halliday, 2000; Türker and Powers, 2002). Thus, the degree of synchrony between PMNs is likely to be higher than for many MN pools because PMNs are driven with large-amplitude, oscillatory patterns of synaptic input (Cohen et al., 1997).
Functionally, synchronous MN output may enhance intramuscular force transmission (Murthy and Fetz, 1994; Baker et al., 1999). Whether diaphragmatic muscle fibers are arranged exclusively in parallel, as in rat and rabbit (Gordon et al., 1989), or also in series, as in cat and dog (Gordon et al., 1989), tension is not simply delivered axially along a fiber. Some tension is transmitted through neighboring or serially arranged fibers and the extracellular matrix. Thus, during submaximal contractions typical of the diaphragm, delivery of force by any given motor unit will vary dynamically with the activity patterns of its neighbors or serial partners (Sheard et al., 2002). The problem this presents for generating smooth, consistent gradations of force could be solved by organizing premotor networks to ensure that anatomically coupled motor units are recruited synchronously as “functional units” (Sheard et al., 2002). Oscillations in inspiratory drive could serve this function. At present, however, whether the spatial relationship of endogenously coactive motor units is consistent with existence of such functional units is not known. The degree of synchronization that ultimately occurs between MNs may be task specific and reflect a balance between competing requirements for fine motor control, which decreases with increasing synchronization (Yao et al., 2000), and efficient contraction, which increases with synchronization (Baker et al., 1999). Because breathing is less dependent than other behaviors on the ability to produce fine gradations in force but requires instead that muscles remain rhythmically active throughout life, activation patterns may have evolved to favor efficiency.
Footnotes
This work was supported by grants from the Marsden Fund, Health Research Council of New Zealand, Lotteries Health, Auckland Medical Research Foundation, New Zealand Neurological Foundation, Maurice and Phyllis Paykel Trust, and National Institutes of Health Grant NS 24742. We thank M. Glasson and M. Navakatykyan for technical assistance.
Correspondence should be addressed to Gregory D. Funk, Department of Physiology, Faculty of Medicine and Dentistry, University of Alberta, 7-55 Medical Science Building, 87 Avenue, 114 Street, Edmonton, Alberta, T6G 2H7, Canada. E-mail: gf{at}ualberta.ca.
Copyright © 2003 Society for Neuroscience 0270-6474/03/238152-07$15.00/0