|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Department of Physiology and Biophysics, University of Washington, Seattle, Washington
Submitted 18 September 2007; accepted in final form 26 March 2008
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
When stimulated with time-varying inputs most neurons studied so far are characterized by a range of input amplitudes and frequencies where firing reliability is highest (Fellous et al. 2001
; Hunter and Milton 2003
; Hunter et al. 1998
). This property of neurons has been called firing resonance. It is not known whether HMs show firing resonance and, if so, how this resonance-firing frequency is related to the dominant frequency of the inspiratory-phase synaptic drive.
To address these questions we intracellularly recorded from HMs in an in vitro brain stem slice preparation from the juvenile mouse. We recorded HM firing patterns in response to injected steady and swept (ZAP) sine-wave current functions. Our data showed that HM firing reliability and spike-timing precision depend on the amplitude and/or frequency of the sine-wave input. We found that reliability and precision are highest for intermediate amplitude stimuli and frequencies where the cells phase-locked to the stimulus. The consequences of these findings for maintenance of synchronous firing among HMs during inspiration are discussed.
| METHODS |
|---|
|
|
|---|
In vitro experiments were performed on nonrhythmically active medullary slice preparations from Swiss–Webster mice (P10–P16). Mice were anesthetized with isoflurane and killed by decapitation in accordance with the regulations of the University of Washington Institutional Animal Care and Use Committee. The brain stem was removed and 300-µm transverse slices containing the hypoglossal motor nucleus were cut. Slices were incubated for 1 h at 37°C before recording. For intracellular recordings slices were transferred to a perfusion chamber mounted on the stage of an upright Zeiss microscope (Axoscope) equipped with infrared differential interference contrast (DIC) optics. The perfusion solution (1–2 ml/min) contained (in mM): 119 NaCl, 2.5 KCl, 1.3 MgSO4, 1 NaH2PO4, 26.2 NaHCO3, 11 D-glucose, and 2.5 CaCl2, gassed with 95% O2-5% CO2 to maintain a pH of 7.4. The temperature in the chamber was maintained at 32 ± 1°C. To reduce background synaptic noise, all recordings were done in the presence of blockers of
-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA)–mediated glutamatergic [6,7-dinitroquinoxaline-2,3-dione (DNQX), 10 µM], GABAergic [2-(3-carboxypropyl)-3-amino-6-methoxyphenyl-pyridazinium bromide (SR95531), 0.5 µM], and glycinergic (strychnine, 1 µM) synaptic transmission.
To assess the firing-resonance properties of neurons it was essential that stable recordings could be maintained for long periods of time to permit examination of responses to repeated stimuli under a range of experimental conditions. We used two different patch solutions to maximize recording stability based on either K-gluconate (Kgluc) or K-methylsulfate (Kmeth). The composition of the Kmeth solution was (in mM): 150 K-meth, 2 HEPES, 2 Mg-ATP, 0.1 K-EGTA, and 15 mM creatine (Chung et al. 1998
). The composition of the Kgluc solution was (in mM): 115 K-gluc, 25 KCl, 1 MgCl2, 9 NaCl, 10 HEPES, 3 K2-ATP, 1 Na-ATP, and 0.2 K-EGTA. The pH of both solutions was adjusted to 7.3 with KOH. For whole cell patch recordings electrodes were pulled from thin-walled borosilicate glass to a resistance of 3–5 M
.
The location of the hypoglossal motor nucleus (nucleus XII) was determined as described in previous publications from this laboratory (Bayliss et al. 1992
; Viana et al. 1993
). The densely packed multipolar cells that together form the XII nucleus were identified under high power using infrared DIC optics. These cells could be easily distinguished from neurons in the adjacent dorsal motor nucleus of the vagus by their large cell bodies and their ability to fire repetitively in response to depolarizing current pulses. To confirm this cell-identification process, several cells were stained intracellularly with biocytin (data not shown). Electrodes were advanced onto visually identified HM neurons, using positive pressure; after formation of a gigaohm seal the cell membrane was ruptured by applying brief suction to the patch.
To avoid intracellular dialysis of the cell with the internal patch solution we also recorded a subset of neurons with "sharp" microelectrodes. Intracellular recording microelectrodes were pulled from thick-walled borosilicate glass to a final resistance of 40–80 M
. First, under visual guidance, the electrodes were placed on the surface of the slice over the hypoglossal motor nucleus and gradually advanced into the slice using a micromanipulator (Sutter Instrument). We continuously monitored bridge balance while the electrodes were advanced by injecting a 5-ms hyperpolarizing current pulse into the electrode. On an increase in electrode resistance we attempted to gain access to the cell interior by briefly ringing the electrode capacitance using the "buzz" control of the amplifier (Axoclamp 2A). We did not find any systematic differences between the firing properties or stability of cells recorded with either patch solution or with sharp or whole cell patch-based techniques (data not shown), and results from all three recording techniques were therefore pooled in this study. We checked for changes in firing pattern over time by repeating the same stimulus at the beginning and end of each trial. We also checked the firing rate in response to DC depolarization at regular intervals during the experimental protocol. We aborted or excluded any recordings where significant changes occurred in firing responses over time.
Analysis of spike-firing resonance
Cells were driven by computer-generated sinusoidal current injections using the external command input of an intracellular amplifier (Axoclamp 2A), operating in bridge mode, through the use of the Clampex 8 software program (Axon Instruments). The same electrode was used for membrane voltage recording and current injections. Sine-wave current injections were superimposed on a DC offset current whose magnitude was adjusted such that the membrane potential of the recorded HM cell was near or at spike threshold. At this membrane potential most neurons fired spontaneously at a low rate. Use of the DC current offset allowed us to use relatively small amplitudes for sine-wave current injections (<0.6 nA). In pilot experiments we found that large-amplitude stimuli were detrimental to the long-term stability of the cells. In HMs current injections of this magnitude typically give rise to membrane potential fluctuations between 5 and 30 mV (see Fig. 1). We stimulated the cells with either steady sine-wave currents (10-s duration) whose frequency was varied stepwise between 2 and 20 Hz or sinusoidal currents of slowly and linearly varying frequency (ZAP functions). The ZAP current stimulus is described by: I(t) = I0 sin [2
f(t) x t], where f(t) = f0 + (fm – f0) x (t/T). The time-dependent frequency of the current f(t) was increased from f0 = 0 Hz to fm = 20 Hz (2 Hz/s) or from f0 = 10 Hz to fm = 40 Hz (3 Hz/s) over a total time period of T = 10 s. In some cells we also tested the reversed protocol (reverse ZAP) where frequency was varied from 20 to 0 or 40 to 10 Hz. Each sine wave or ZAP function was repeated at least twice. Firing responses in 15-s-long digitized data files, including the 10-s period of sine-wave stimulation, were digitized at 2 kHz (Clampex) and detected with an amplitude window discriminator, using Neuromatic software (J. Rothman) written for IgoR Pro (WaveMetrics). Detection threshold was set along the upstroke of the action potential at about 20 mV below spike peak and the time at which the spike occurred was stored with a precision of 500 µs. In addition, in all cells we also injected with short (5-ms) and long (1,000-ms) DC current pulses to characterize action potential characteristics and firing properties at resting membrane potential, respectively. These recordings were digitized at 20 kHz and analyzed using Clampex/Clampfit (Axon Instruments).
|
Reliability and spike-timing precision analysis was performed with custom-written functions implemented in IgorPro. To quantify the reliability of firing (i.e., the likelihood that a cell will fire an action potential at a certain time relative to the onset of the stimulus), the cells were stimulated with repeated ZAP current stimuli (10 trials). The times of occurrence of spikes from the start of the stimulus (T = 0) were collected in 3-ms-wide bins and a peristimulus time histogram (PSTH) was constructed from the number of spikes that fell within each bin. Only bins that contained more than one spike were used for analysis. Reliability was expressed as the number of spikes/bin divided by the number of trials. Highest reliability (1 spike/bin/stimulus) was obtained if the cell always fired a spike at the same latency (i.e., fell within the same bin) at each trial. This reliability measure is sensitive to both the presence of a spike (i.e., firing probability) and spike-timing precision. A low reliability might therefore be due to either a low probability of spike firing or/and low precision. To examine the probability of spike firing we calculated the average number of spikes per sinusoidal cycle and the trial-to-trial variability of the number of spikes/cycle, expressed as the coefficient of variability (CV = SD/mean). Spike-timing precision was examined using phase angle, phase error, and spike jitter of single-trial responses to steady sine-wave stimuli. Phase angle was calculated from the time of the nadir of the trough of the input current time to the time of the first spike in the following half-cycle. The phase error is the SD of this spike latency expressed in degrees, whereas jitter is the SD of spike latency expressed in milliseconds. At the onset of steady sine-wave stimuli, during the first depolarizing half-cycle, the cells often fired more action potentials than in the remainder of the sine-wave stimulus. For this reason the first three-quarter cycle until the first trough in the sine wave was ignored in the spike-timing precision analysis. Values for phase angle, phase error, and jitter were averaged over all remaining cycles of a given 10-s steady sine-wave stimulus. The cells were either entrained by the stimulus or fired irregularly. To quantify the regularity of firing we adopted the modified coefficient of variation (CV2) proposed by Holt et al. (1996)
: CV2 = 2 x abs [(Ii+1 – Ii)]/(Ii+1 + Ii). This measure of variability for an interval at t = i seconds compares only two adjacent spike intervals Ii+1 and Ii. A perfectly regular spike train has a CV2 of 0, whereas spikes occurring randomly (Poisson process) have a CV2 of 1.
Statistics
Statistical comparisons were performed using the t-test, with significance set at P < 0.05. Results are expressed as means ± SE.
| RESULTS |
|---|
|
|
|---|
Response of HMs to sine-wave stimulation
After the firing properties of the recorded HMs at rest were established we depolarized the cells to a membrane potential near or at firing threshold by DC current injection. We stimulated the cells with 10-s-long sine-wave current inputs (Fig. 1), superimposed on this steady DC current injection. This protocol was chosen to simulate the synaptic drive that HMs receive during inspiration, consisting of a depolarizing envelope on top of which membrane potential oscillations are superimposed (Parkis et al. 2003
; Smith and Denny 1990
). In addition, by bringing the cells to firing threshold we were able to reduce the amplitude of the sine-wave current, which greatly improved the stability of our recordings (see METHODS). All cells studied fired only during the depolarizing half-cycle of the sine current stimulus. We used two types of sine-wave stimulation protocols: one in which sine frequency was varied in steps (steady sine-wave stimulation) and another in which sine-wave frequency was varied continuously (swept sine wave or "ZAP" stimulation). Most cells in this study were stimulated with both protocols.
In response to steady sine-wave stimulation (varied stepwise between 2 and 20 Hz) HM neurons showed several distinct firing patterns depending on the input frequency (n = 14). At low input frequencies (generally <5 Hz) cells fired multiple spikes/cycle (i.e., "burst" mode), at intermediate frequency cells fired one spike/cycle, whereas at higher frequencies occasional cycle skipping occurred (Fig. 1A). By plotting the number of spikes per sine-wave cycle against frequency it was apparent that there was a range of input frequencies where the cells fired exactly one spike/cycle and this frequency range we called the 1:1 phase-locking range (Fig. 1B1). Instantaneous firing frequency was highest with low input frequencies when the cells fired in "burst" mode (Fig. 1B2). At intermediate sine frequencies the firing frequency of the cells faithfully followed the input frequency (1:1 phase-locking range; Fig. 1B2). Beyond the maximum firing frequency at which the cells were able to phase-lock 1:1 to the stimulus (i.e., the "critical" or maximum frequency, Fmax) cells were unable to fire every cycle and the firing frequency of the cells was less than the input frequency. All 14 cells in this group showed the three distinct firing ranges ("burst" mode, phase-locking, and cycle skipping), but the sine frequency at which these firing patterns occurred varied from cell to cell. To study the firing patterns to sine-wave stimulation in more detail we stimulated a group of HMs (n = 19 cells) with ZAP functions. Firing patterns obtained with ZAP stimulation showed many similarities compared with those obtained with steady sine-wave stimulation, including "burst" firing, 1:1 phase-locking, and cycle skipping as shown by the representative cell in Fig. 2, A and B (same cell as in Fig. 1). When we plotted instantaneous firing frequency (f) against time (t), the resulting f–t relationships clearly showed these distinct regions of firing responses. In "burst" mode (region "a" in Fig. 2A1) the cells fired more than one spike/cycle in a relatively high frequency burst during the depolarizing half-cycle of the sine wave. The length of the silent period between bursts (i.e., the reciprocal of the interburst frequency) reflected the frequency of the sine-wave input. At intermediate ZAP frequencies (between 7 and 21 Hz in the example shown in Fig. 2) the neurons abruptly transitioned from "burst" firing to firing one spike/cycle. Between this minimum frequency (Fmin) and the maximum frequency (Fmax) HMs coded the frequency content by phase-locking 1:1 to the input stimulus (Fig. 2A, region b). Beyond Fmax irregular or chaotic firing occurred due to occasional cycle skipping (Fig. 2A2, region c). At even higher ZAP frequency additional phase-locking regimens were revealed, often separated by regions of irregular firing (Figs. 2A2 and 3 ). Most cells were able to fire every second (1:2), every third (1:3), and sometimes every fourth (1:4) cycle, phase-locked to subharmonics of the fundamental frequency of the ZAP stimulus (Figs. 2B and 3). On average the minimum input frequency (Fmin) for 1:1 phase-locking was 4.5 ± 0.36 Hz (n = 19 cells). The maximum frequency for 1:1 phase-locking (Fmax) in this group of cells was 11.4 ± 0.93 Hz. Although in individual cells Fmin and Fmax were influenced by stimulus amplitude (see Fig. 4 ), average Fmin and Fmax values were not correlated with stimulus amplitude (average stimulus amplitude, 253 ± 22 pA; range 100–400 pA) in this group of HMs (n = 19 cells; Pearson correlation, P > 0.05). In addition, the maximum and minimum frequencies and width of the frequency range where HMs phase-locked 1:1 to the ZAP stimulus were not correlated with the input resistance of the cells (Pearson correlation, P > 0.05).
|
|
|
These results show that the frequency information present in a time-varying input stimulus is encoded by precisely timed spikes when HMs phase-lock to the stimulus, indicating that in this frequency range they function as temporal encoders. The frequency range where this occurs varies from cell to cell.
Phase-locking region depends on DC depolarization and stimulus amplitude
Firing resonance in HMs depended on the level of DC membrane depolarization (Fig. 3A). When the cells were hyperpolarized from spike threshold the phase-locking regions became smaller and at high-input frequencies the cells stopped firing altogether (Fig. 3A1). Both the minimum and maximum frequencies of the 1:1 phase-locking range were influenced by the level of depolarization. On average, the minimum frequency for 1:1 phase-locking (Fmin) increased by 2.1 ± 0.20 Hz and Fmax by 5.8 ± 1.03 Hz for a +10-mV membrane depolarization in the voltage range between rest and spike threshold (n = 4 cells; Fig. 3B). Since Fmax increased more than Fmin the 1:1 phase-locking range was expanded by 68% in this group of cells by this manipulation. DC depolarization also resulted in an expansion of the regions of irregular or chaotic firing that flanked phase-locking regions (Fig. 3A3, arrows).
Stimulus amplitude also influenced the shape of the f–t relationship (Fig. 4A). With increased stimulus amplitude the instantaneous firing frequency and the frequency range over which the cells fired in "burst" mode increased. Peak instantaneous firing frequency in "burst" mode increased from 15.0 ± 1.7 to 22.2 ± 2.6 Hz for an increase in ZAP amplitude from 200 to 300 pA (P < 0.05, paired t-test, n = 7 cells). Increasing stimulus amplitude also shifted the 1:1 phase-locking region to higher frequencies (Fig. 4A). Mean Fmin and mean Fmax increased by 1.32 ± 0.47 and 2.04 ± 0.44 Hz, respectively, for a 100-pA increase in ZAP amplitude from 200 to 300 pA (both increases were significant at P < 0.05, paired t-test, n = 9 cells; Fig. 4C). Although Fmax increased more on average than Fmin, the resulting 11% increase in width of the 1:1 phase-locking range by increasing stimulus amplitude by 100 pA was not significant (P > 0.05). Increasing stimulus amplitude also increased the irregular firing regions at the transition between phase-locking regions (Fig. 4A). In the example shown in Fig. 4, the cell was able to smoothly switch between 1:1 and 1:2 phase-locking with a 100-pA peak-to-peak amplitude stimulus (Fig. 4B, top trace), but fired irregularly between 1:1 and 1:2 phase-locking when the ZAP stimulus amplitude was doubled (Fig. 4B, bottom trace). These results indicate that the spike-firing rate carries information about the stimulus amplitude when the input frequency is low and the cells fire in "burst" mode (so-called rate coding). The timing of individual spikes carries information about the frequency content of the stimulus when the cells are in "phase-locking" mode ("temporal coding"). In this mode a change in spike-firing frequency requires a change in input frequency. The amplitude of the input stimulus influences at what frequency the switch from rate to temporal coding occurs and also the frequency range for phase-locking and cycle skipping.
Although there was hysteresis in the number of spikes/cycle when forward and reverse stimuli were compared in the same cell (forward ZAP stimuli overall elicited more firing; see Fig. 6), we found that the range for 1:1 phase-locking was similar for these types of stimuli. Fmin and Fmax for forward and reverse ZAP stimuli were not significantly different (n = 10 cells; paired t-test, P > 0.05).
|
Synchronicity of firing of assemblies of neurons depends critically in the ability of the cells to fire with high reliability and precision to repeated stimuli. We examined these properties in HMs by repeatedly stimulating HMs with ZAP stimuli. Such an experiment is illustrated in Fig. 5. This cell was repeatedly (10 trials) stimulated with a 0–20 Hz ZAP function. From the rasterplots (Fig. 5, A and B) and superimposed f–t curves (Fig. 5C) it is clear that this cell fired consistently at a specific sine-wave cycle from trial to trial from the beginning of the stimulus to the end of the 1:1 phase-locking range, but less consistently afterward. This is further demonstrated in Fig. 5D, which shows that the variability in the average number of spikes/cycle (as measured by the coefficient of variation) is smallest when the cell phase-locks 3:1 (burst mode) or 1:1 to the stimulus. Repeatability from trial to trial did not have the same frequency profile as reliability, however. Although the cell in Fig. 5 fired three spikes at every trial at the first cycle of the stimulus (i.e., high repeatability), it did so with low temporal precision (Fig. 5B, left). When the cell was entrained 1:2 the opposite was true, it now fired with high precision but low repeatability because it did not always fire in the same cycle from trial to trial (Fig. 5B, right). Low precision at low ZAP frequencies and low repeatability at high frequencies resulted in roughly bell-shaped reliability curves (Fig. 5E). Reliability (which depends on both the presence and timing of a spike and was expressed as the number of spikes/bin/stimulus; see METHODS) peaked when the cell phase-locked 1:1 to the stimulus at intermediate input frequencies. Experiments similar to those shown in Fig. 5 were performed in eight other HMs. Peak reliability, as determined from cubic spline fits to PSTH data, was found at an average input frequency of 8.68 ± 1.38 Hz (n = 9 cells; ZAP current peak-to-peak amplitude between 150 and 400 pA). Peak reliability was 0.82 ± 0.05 spike/bin/stimulus and in all but one of these neurons peak reliability occurred within the 1:1 phase-locking region. We found that spike reliability curves were similar when repeated steady sine-wave stimuli were compared with repeated ZAP stimuli in the same neuron (n = 3 cells; data not shown).
|
We tested the effect of stimulus amplitude on spike reliability in two HMs. In both cells overall spike reliability increased when stimulus amplitude was increased (see example in Fig. 6C) and the frequency range of high spike reliability was expanded.
In summary, the low reliability at low- and high-input frequencies suggests that HMs act as a band-pass filter and not as a low-pass filter, with highest reliability in the intermediate 1:1 phase-locking range. Spike reliability and the frequency where reliability is highest are functions of stimulus amplitude.
The timing of spikes within individual cycles varied with input frequency when the cells were stimulated with steady sine-wave stimuli (Fig. 7). Spiking occurred early, during the rising phase of the sine-wave stimulus, when the input sine frequency was low (phase advance; Fig. 7A2). The phase advance gradually decreased with increasing input frequency until a delay of zero (= "preferred" frequency) was reached. At even higher input frequency spikes were delayed with respect to the sine-wave stimulus (phase delay). The preferred frequency, as determined by linear interpolation at zero phase delay, was 8.7 ± 1.0 Hz (n = 10 cells), a value very close to the frequency where optimum reliability occurred (8.3 Hz; see previous text). Spike jitter—a measure of spike-timing precision (see METHODS)—generally decreased with increasing sine-wave frequency (Fig. 7A3). Thus the variability in the latency of the occurrence of spikes became smaller with increasing input frequency (i.e., spike-timing precision was enhanced). The relationship between phase error, which in contrast to jitter is frequency invariant, and sine-wave frequency was more complex than jitter-–frequency curves but it also tended to become smaller with increasing sine frequency (Fig. 7, A3 and A4). Within the 1:1 phase-locking region jitter and phase error were not constant. Because jitter—and to a lesser extent phase error—decreased with increasing sine frequency, timing precision was greatest when the cells phase-locked at or near Fmax (Fig. 7, A3 and A4). When the cells were driven at sine frequencies higher than Fmax phase error increased dramatically because the phase angle varied from cycle to cycle when the cells were not entrained by the stimulus (Fig. 7B). In other HMs, where several phase-locking regions were present, additional regions of high spiking precision were evident (Fig. 8). For example, interspike interval (ISI) variability (as measured by the local coefficient of variation; see METHODS), phase error, and jitter reached local minima not only when the cells phase-locked 1:1 to the sine-wave stimulus but also when they phase-locked 1:2 or 1:3 (Fig. 8). Thus HMs encode a time-varying input with highest spike-timing precision when the cell phase-locks to the fundamental frequency or subharmonics of this frequency.
|
|
| DISCUSSION |
|---|
|
|
|---|
Spike-timing reliability and precision have been shown to be frequency dependent in other neuronal types in the brain (Fellous et al. 2001
; Leung and Yu 1998
; Pike et al. 2000
; Schreiber et al. 2004b
). The range of frequencies where precision and reliability are highest in HMs is similar to the preferred frequency range in hippocampal (Leung and Yu 1998
; Pike et al. 2000
) and cortical pyramidal cells (Fellous et al. 2001
), but lower than that in hippocampal and cortical interneurons (Fellous et al. 2001
; Pike et al. 2000
). The preferred frequency is a function of the properties of the ionic currents responsible for action potential generation, ISI duration, and spike-frequency adaptation (Fuhrmann et al. 2002
; Schreiber et al. 2004b
). Differences in the composition of ion channels are most likely responsible for differences in preferred frequencies among neuronal cell types. For instance, slow potassium channels have been shown to be important for setting the preferred frequency (Fellous et al. 2001
; Schreiber et al. 2004b
). In modeling studies (Schreiber et al. 2004b
) it was found that a higher peak conductance for slow potassium channels, such as the muscarinic potassium channel (KM) or a calcium-dependent potassium channel (KCa), decreased the preferred frequency of cortical neurons (but increased spike reliability). In contrast, an increase in sodium conductance, such as the fast (Naf) and persistent (Nap) sodium currents, increased the preferred frequency (but decreased reliability). These conductances and their kinetics also influenced the width of the 1:1 phase-locking range in these studies. It is conceivable that the wide ranges of preferred frequencies and widths of phase-locking ranges that we found in our study are a reflection of the heterogeneity in density and/or kinetics of these conductances among hypoglossal motoneurons. Neuromodulation of ion channels can alter the preferred frequencies and increase spike-timing precision and reliability by adapting the intrinsic oscillatory properties of neurons to facilitate entrainment of these cells to certain input frequencies (Lawrence et al. 2006
).
Phase-locking depends on input frequency, stimulus amplitude, and level of depolarization
When HMs were brought close to firing threshold, spike firing was entrained by the sine-wave input, depending on sine-wave frequency. At low input frequency HMs fired highly repeatable, stereotypical "bursts" of action potentials during the depolarizing phase of the sine-wave cycle (i.e., more than one spike/cycle). This burst firing mode is most suitable to reliably evoke a postsynaptic response when the timing of the spikes is not important, since spike-timing precision in burst mode was low. This study and those by others (Carandini et al. 1996
; Hunter and Milton 2003
; Schreiber et al. 2004b
) show that small increases in the input frequency can lead a neuron to shift from coding spike rate (burst mode) to coding spike timing (phase-locking mode). Phase-locking was most prevalent at intermediate and high-input frequencies, but this frequency range varied from cell to cell. Phase-locked responses could result in the firing of exactly one spike/cycle (1:1 mode), or in one spike for each n cycle(s) of the sine-wave stimulus (with n being an integer number
2). This phase-locked firing mode was characterized by high temporal precision. Phase-locking did not always result in high reliability since the cells did not always fire in the same cycle when cycle skipping was present (1:n mode), similar to the results obtained by Hunter and Milton (2003)
. Thus the 1:1 phase-locking region is optimized for encoding input signals with both high reliability and high precision.
The range of frequencies where the cells phase-locked depended on the mean level of depolarization and the amplitude of the sinusoidal inputs. The relationship between reliability, membrane depolarization, stimulus amplitude, and input frequency has been explored in several types of neurons (Beierholm et al. 2001
; Fellous et al. 2001
; Hunter and Milton 2003
; Smith et al. 2000
). In general, these studies show that the most robust areas of high reliability are found when the cell phase-locks 1:1 to the stimulus using intermediate relative stimulus amplitude. Small-amplitude input fluctuations are unreliable because under these conditions noise becomes the most important determinant of spike timing (Hunter and Milton 2003
). The preferred frequency range increased with increasing amplitude of the input stimulus, similar to results obtained by others (Carandini et al. 1996
; Fellous et al. 2001
), but large-amplitude voltage fluctuations eventually lead to cell damage and decreased reliability (Fellous et al. 2001
). Thus intermediate-amplitude signals in the band-pass frequency range delivered when the neurons are depolarized to near spike threshold will be transmitted more reliably than those outside the band-pass frequency range or when the cells are hyperpolarized.
Several studies have suggested that firing resonance emerges from subthreshold resonance properties (Erchova et al. 2004
; Hutcheon et al. 1996
; Schreiber et al. 2004a
), but others have not found such a clear correlation (Brumberg and Gutkin 2007
; Carandini et al. 1996
; Haas and White 2002
; Nowak et al. 1997
). These discrepancies might be explained in part by differences in subthreshold resonance properties between cell types (Erchova et al. 2004
; Haas and White 2002
). In resonant cells, the amplitude of their membrane potential responses to time-varying subthreshold currents peaks at the resonance frequency, whereas nonresonant cells act as low-pass filters. In resonant cells firing occurs preferentially at the subthreshold resonance frequency when they are stimulated with a suprathreshold stimulus. Subthreshold resonance has been described in HMs that peaks at 4 Hz (Nguyen et al. 2004
). Our own results indicate that firings in HMs are either not resonant or only weakly resonant in the subthreshold range (unpublished observations). HMs therefore act predominantly as low-pass filters in the membrane voltage range between rest and spike threshold for input frequencies >3–5 Hz. This property is responsible for the cessation of firing at high frequencies that was observed in many HM neurons when they are moved away from spike threshold by DC current injection (see Fig. 3A1). The frequency range where firing occurs can be widened by increasing the stimulus amplitude or by membrane depolarization. At depolarized membrane potentials small-amplitude currents can more easily generate spikes and the cell fires at higher input frequencies despite low-pass filtering of the input signal (Hutcheon et al. 1996
). The firing resonance range can also be increased by the addition of stochastic noise to the sine-wave input in thalamocortical neurons (Reinker et al. 2004
). Whether this is also the case in HMs remains to be established.
Relationship between spike jitter, phase, and input frequency
Our results show that HMs show phase advance at low-frequency stimulation and phase delay at high frequency, separated by the preferred frequency at which no phase shift occurs. Similar phase-frequency curves have been described for thalamic relay neurons (Smith et al. 2000
), regular spiking cortical neurons (Carandini et al. 1996
), and hippocampal CA neurons (Leung and Yu 1998
). The mechanism behind this property has been examined in modeling studies (Fuhrmann et al. 2002
) and has been explained as an interplay between two opposing mechanisms: spike-frequency adaptation and the time constant of firing rate dynamics. Firing rate dynamics causes a neuron to delay its response to stimulation and dominates at high frequency. Spike-frequency adaptation causes the response to a sustained input to gradually decrease with time. The slow rising phase with low-frequency sinusoidal stimulation leads to accumulation of adaptation and thus the firing probability will be higher before the peak of the sinusoidal stimulus than after, leading to a phase advance. At the preferred frequency (phase = zero) these two mechanisms balance each other (Fuhrmann et al. 2002
).
Our results show that spike jitter generally decreases with increasing sine frequency. Studies in neocortical neurons have shown that spike jitter is less for high-frequency than that for low-frequency input (Mainen and Sejnowski 1995
; Nowak et al. 1997
). The mechanism behind this observation is not known, but is has been suggested that increased timing precision at high stimulus frequencies is related to the time that the neuron spends near spike threshold. The time spent near threshold is brief for high- than that for low-frequency stimulation. Therefore at high-frequency stimulation there is only a small probability that random fluctuations in membrane potential will exceed spike threshold and thus spike jitter is reduced (Nowak et al. 1997
).
Regions of phase-locking, flanked by regions of irregular firing when HMs were driven with a time-varying stimulus, have also been described in supragranular pyramidal cells under similar conditions (Brumberg 2002
; Brumberg and Gutkin 2007
). Together this firing pattern forms a so-called devil's staircase. The regions of irregular firing are influenced by stimulus amplitude and level of depolarization and they are characterized by low precision and reliability. Our studies suggest that there exists an optimal level of depolarization and stimulus amplitudes where in response to sine-wave stimuli irregular firing in HMs can be largely avoided. If the oscillation frequency or the power of the oscillations change (for some examples of this behavior see Sebe et al. 2006
), HMs will have to adjust their encoding to maintain phase-locking to the stimulus. A widest possible phase-locking range with "smooth" transitions between 1:1 and 1:n phase-locking ranges will ensure that HMs will be able to maintain synchrony over a wide range of time-varying inputs while avoiding regions of irregular firing. The 1:1 phase-locking region is most robust in the presence of noise (Hunter and Milton 2003
), but whether in vivo synaptic noise stabilizes firing in HM neurons remains to be investigated.
Importance of HM spike timing and precision for synchronous oscillation in HM output
The motoneurons in the hypoglossal nucleus are part of the brain stem inspiratory network. From studies performed in in vitro preparations in neonatal or juvenile rodents it is known that inspiratory motoneurons are driven by a slow depolarizing envelope on top of which fast (20- to 50-Hz) membrane potential oscillations are superimposed (Funk and Parkis 2002
; Parkis et al. 2003
). Our results indicate that with time-varying synaptic input in this frequency range HMs operate in the phase-locking mode. Short-timescale inspiratory-phase oscillations in this frequency range recorded from the hypoglossal nerve during these inspiratory bursts are characterized by clusters of action potentials that occur at a regular interval separated by periods of little or no spike firing (Sebe et al. 2006
), suggesting that HMs fire synchronously. In phase-locking mode groups of HMs are able to produce precise and reliable spike output, which would transform common oscillatory synaptic input into synchronous volleys of spikes. High reliability of spike timing enhances synchronization of groups of functionally related neurons activated by a common synaptic input (Usrey and Reid 1999
).
The function of inspiratory-phase synchronous oscillations in the brain stem is not known. Inspiratory-phase synchronous oscillations are not required for respiratory rhythm generation, but they are important in shaping the pattern of inspiratory discharge. At the motoneuron level, oscillations increase the input–output efficiency of and the precision of spike timing (Parkis et al. 2003
). Oscillations provide for a mechanism that avoids accumulation of spike jitter, as occurs with DC stimuli (Mainen and Sejnowski 1995
), by alternating spiking and nonspiking periods. Membrane hyperpolarization between half-cycles allows for recovery from inactivation of voltage-activated ion currents involved in spike generation. Phrenic motoneurons tend to fire within a couple of milliseconds from the short-timescale oscillatory peaks in inspiratory synaptic drive (Parkis et al. 2003
). This timescale is similar to the minimum in spike jitter that we found in this study when the cells phase-lock to the stimulus (see Figs. 7 and 8). The DC current injection that we used in our study mimics the membrane depolarization of the slow depolarizing inspiratory envelope, but it did not account for the increased membrane conductance accompanying this excitatory synaptic drive. It is possible that the increased conductance—and associated decrease in membrane time constant—will improve spike-timing precision and shift the optimum frequency for phase-locking to input frequencies higher than those we found in our quiescent, synaptically isolated slice preparation.
Similar to what has been found in spinal interneurons (Beierholm et al. 2001
), it appears that the combination of long- and short-timescale synaptic inputs to HMs during inspiration reliably drive the cells with high timing precision. The depolarizing envelope during inspiration brings the cells close to spike threshold, so that high reliability can be achieved with relatively small amplitude membrane voltage fluctuations (similar to a very low frequency sine-wave input stimulus). The high-frequency oscillation on top of this depolarizing envelope serves to give high spike reliability and precision.
Power spectra of inspiratory oscillations recorded in vivo from motor nerves in adult animals show two peaks in the 20- to 50- and 50- to 150-Hz ranges (Richardson and Mitchell 1982
; Smith and Denny 1990
). Single-fiber recordings from the phrenic nerve during inspiratory activity showed that the discharge frequency of individual motoneurons is lower than that of the whole nerve during short-timescale oscillations in vivo (Richardson and Mitchell 1982
). In these studies it was found that single fibers sometimes skipped cycles, firing every second, third, or fourth cycle. It is possible that individual HMs in vivo phase-lock to subharmonics of the fundamental oscillation frequency, skipping cycles and firing at a lower frequency than the population fundamental. We hypothesize that, even though spikes are not elicited at every cycle in individual HMs, the high-frequency input signal still serves to synchronize neurons by reducing jitter and improving spike timing. It would therefore be of interest to establish the firing resonance properties of individual HMs during inspiration in vivo and to relate these to the motor output of the XII nucleus as a whole.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: H. van Brederode, UW Physiology and Biophysics, 1705 NE Pacific St., Harris Hydraulics Rm 309, Box 357290, Seattle, WA 98195-7290 (E-mail: hansvb{at}washington.edu)
| REFERENCES |
|---|
|
|
|---|
Beierholm U, Nielsen CD, Ryge J, Alstrom P, Kiehn O. Characterization of reliability of spike timing in spinal interneurons during oscillating inputs. J Neurophysiol 86: 1858–1868, 2001.
Berger AJ, Sebe J. Developmental effects of ketamine on inspiratory hypoglossal nerve activity studied in vivo and in vitro. Respir Physiol Neurobiol 157: 206–214, 2007.[CrossRef][Web of Science][Medline]
Brumberg JC. Firing pattern modulation by oscillatory input in supragranular pyramidal neurons. Neuroscience 114: 239–246, 2002.[CrossRef][Web of Science][Medline]
Brumberg JC, Gutkin BS. Cortical pyramidal cells as non-linear oscillators: experiment and spike-generation theory. Brain Res 1171: 122–137, 2007.[CrossRef][Web of Science][Medline]
Carandini M, Mechler F, Leonard CS, Movshon JA. Spike train encoding by regular-spiking cells of the visual cortex. J Neurophysiol 76: 3425–3441, 1996.
Chung I, Zhang Y, Eubanks JH, Zhang L. Attenuation of hypoxic current by intracellular applications of ATP regenerating agents in hippocampal CA1 neurons of rat brain slices. Neuroscience 86: 1101–1107, 1998.[CrossRef][Web of Science][Medline]
Erchova I, Kreck G, Heinemann U, Herz AV. Dynamics of rat entorhinal cortex layer II and III cells: characteristics of membrane potential resonance at rest predict oscillation properties near threshold. J Physiol 560: 89–110, 2004.
Fellous JM, Houweling AR, Modi RH, Rao RP, Tiesinga PH, Sejnowski TJ. Frequency dependence of spike timing reliability in cortical pyramidal cells and interneurons. J Neurophysiol 85: 1782–1787, 2001.
Fuhrmann G, Markram H, Tsodyks M. Spike frequency adaptation and neocortical rhythms. J Neurophysiol 88: 761–770, 2002.
Funk GD, Parkis MA. High frequency oscillations in respiratory networks: functionally significant or phenomenological? Respir Physiol Neurobiol 131: 101–120, 2002.[CrossRef][Web of Science][Medline]
Holt GR, Softky WR, Koch C, Douglas RJ. Comparison of discharge variability in vitro and in vivo in cat visual cortex neurons. J Neurophysiol 75: 1806–1814, 1996.
Hunter JD, Milton JG. Amplitude and frequency dependence of spike timing: implications for dynamic regulation. J Neurophysiol 90: 387–394, 2003.
Hunter JD, Milton JG, Thomas PJ, Cowan JD. Resonance effect for neural spike time reliability. J Neurophysiol 80: 1427–1438, 1998.
Hutcheon B, Miura RM, Puil E. Subthreshold membrane resonance in neocortical neurons. J Neurophysiol 76: 683–697, 1996.
Kaczorowski CC, Disterhoft J, Spruston N. Stability and plasticity of intrinsic membrane properties in hippocampal CA1 pyramidal neurons: effects of internal anions. J Physiol 578: 799–818, 2007.
Lawrence JJ, Statland JM, Grinspan ZM, McBain CJ. Cell type-specific dependence of muscarinic signalling in mouse hippocampal stratum oriens interneurones. J Physiol 570: 595–610, 2006.
Leung LS, Yu HW. Theta-frequency resonance in hippocampal CA1 neurons in vitro demonstrated by sinusoidal current injection. J Neurophysiol 79: 1592–1596, 1998.
Mainen ZF, Sejnowski TJ. Reliability of spike timing in neocortical neurons. Science 268: 1503–1506, 1995.
Nguyen QT, Wessel R, Kleinfeld D. Developmental regulation of active and passive membrane properties in rat vibrissa motoneurones. J Physiol 556: 203–219, 2004.
Nowak LG, Sanchez-Vives MV, McCormick DA. Influence of low and high frequency inputs on spike timing in visual cortical neurons. Cereb Cortex 7: 487–501, 1997.
Parkis MA, Feldman JL, Robinson DM, Funk GD. Oscillations in endogenous inputs to neurons affect excitability and signal processing. J Neurosci 23: 8152–8158, 2003.
Pike FG, Goddard RS, Suckling JM, Ganter P, Kasthuri N, Paulsen O. Distinct frequency preferences of different types of rat hippocampal neurones in response to oscillatory input currents. J Physiol 529: 205–213, 2000.
Reinker S, Puil E, Miura RM. Membrane resonance and stochastic resonance modulate firing patterns of thalamocortical neurons. J Comput Neurosci 16: 15–25, 2004.[CrossRef][Web of Science][Medline]
Richardson CA, Mitchell RA. Power spectral analysis of inspiratory nerve activity in the decerebrate cat. Brain Res 233: 317–336, 1982.[CrossRef][Web of Science][Medline]
Sawczuk A, Powers RK, Binder MD. Spike frequency adaptation studied in hypoglossal motoneurons of the rat. J Neurophysiol 73: 1799–1810, 1995.
Schreiber S, Erchova I, Heinemann U, Herz AV. Subthreshold resonance explains the frequency-dependent integration of periodic as well as random stimuli in the entorhinal cortex. J Neurophysiol 92: 408–415, 2004.
Schreiber S, Fellous JM, Tiesinga P, Sejnowski TJ. Influence of ionic conductances on spike timing reliability of cortical neurons for suprathreshold rhythmic inputs. J Neurophysiol 91: 194–205, 2004.
Sebe JY, van Brederode JF, Berger AJ. Inhibitory synaptic transmission governs inspiratory motoneuron synchronization. J Neurophysiol 96: 391–403, 2006.
Sejnowski TJ, Paulsen O. Network oscillations: emerging computational principles. J Neurosci 26: 1673–1676, 2006.
Smith A, Denny M. High-frequency oscillations as indicators of neural control mechanisms in human respiration, mastication, and speech. J Neurophysiol 63: 745–758, 1990.
Smith GD, Cox CL, Sherman SM, Rinzel J. Fourier analysis of sinusoidally driven thalamocortical relay neurons and a minimal integrate-and-fire-or-burst model. J Neurophysiol 83: 588–610, 2000.
Turker KS, Powers RK. The effects of common input characteristics and discharge rate on synchronization in rat hypoglossal motoneurones. J Physiol 541: 245–260, 2002.
Usrey WM, Reid RC. Synchronous activity in the visual system. Annu Rev Physiol 61: 435–456, 1999.[CrossRef][Web of Science][Medline]
Viana F, Bayliss DA, Berger AJ. Calcium conductances and their role in the firing behavior of neonatal rat hypoglossal motoneurons. J Neurophysiol 69: 2137–2149, 1993.
Viana F, Bayliss DA, Berger AJ. Repetitive firing properties of developing rat brainstem motoneurones. J Physiol 486: 745–761, 1995.
This article has been cited by other articles:
![]() |
V. Marchenko and R. F. Rogers GABAAergic and Glycinergic Inhibition in the Phrenic Nucleus Organizes and Couples Fast Oscillations in Motor Output J Neurophysiol, April 1, 2009; 101(4): 2134 - 2145. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |