Abstract
The modification of first-spike latencies by low-threshold and inactivating K+ currents (IA) have important implications in neuronal coding and synaptic integration. To date, cells in which first-spike latency characteristics have been analyzed have shown that increased hyperpolarization results in longer first-spike latencies, producing a monotonic relationship between first-spike latency and membrane voltage. Previous work has established that cerebellar stellate cells express members of the Kv4 potassium channel subfamily, which underlie IA in many central neurons. Spike timing in stellate cells could be particularly important to cerebellar output, because the discharge of even single spikes can significantly delay spike discharge in postsynaptic Purkinje cells. In the present work, we studied the first-spike latency characteristics of stellate cells. We show that first-spike latency is nonmonotonic, such that intermediate levels of prehyperpolarization produce the longest spike latencies, whereas greater hyperpolarization or depolarization reduces spike latency. Moreover, the range of first-spike latency values can be substantial in spanning 20-128 ms with preceding membrane shifts of <10 mV. Using patch clamp and modeling, we illustrate that spike latency characteristics are the product of an interplay between IA and low-threshold calcium current (IT) that requires a steady-state difference in the inactivation parameters of the currents. Furthermore, we show that the unique first-spike latency characteristics of stellate cells have important implications for the integration of coincident IPSPs and EPSPs, such that inhibition can shift first-spike latency to differentially modulate the probability of firing.
Introduction
Many neurons actively generate long latencies to first spike that are sensitive to the preceding membrane voltage (Byrne et al., 1979; Storm, 1988; Kanold and Manis, 1999). Traditionally, a low threshold and fast inactivating K+ current (IA) is required. IA is activated by hyperpolarization and subsequently opposes any membrane depolarization, resulting in a delay to reach firing threshold. In cells expressing IA, the relationship between latency and membrane voltage has been shown to be monotonic, such that greater hyperpolarization causes longer delays to first spike (Shibata et al., 2000; Kanold and Manis, 2005). First-spike latency has been studied in many sensory systems in which the magnitude of a sensory input can be more effectively encoded by action potential timing than through changes in spike frequency (Thorpe et al., 1996; Heil, 2004; Johansson and Birznieks, 2004; VanRullen et al., 2005). In addition to spike timing, IA has been shown to play an important role in synaptic integration by influencing the subthreshold membrane properties (Hoffman et al., 1997; Schoppa and Westbrook, 1999). Thus, understanding the intrinsic mechanisms involved in setting first-spike latency is of general relevance to processing in central neurons.
The cerebellum is a motor structure that requires high-temporal resolution to maintain appropriate motor control (Timmann et al., 1999; Kistler et al., 2000). Cerebellar stellate cells are inhibitory interneurons that receive excitatory and inhibitory inputs from parallel fibers and other stellate cells, respectively (Palay and Chan-Palay, 1974; Midtgaard, 1992a,b). By projecting onto Purkinje cell dendrites, stellate cells can modify the output of the cerebellar cortex (Midtgaard, 1992b; Hausser and Clark, 1997; Jorntell and Ekerot, 2002; Mittmann et al., 2005). Indeed, even single spikes in stellate cells have been shown to significantly delay Purkinje cell firing (Hausser and Clark, 1997), revealing the potential impact of stellate cell spike timing. Stellate cells have been shown to express members of the Kv4 family (Hsu et al., 2003), a channel that often underlies IA. Despite the expression of IA, stellate cells have not been reported to generate first-spike latencies that are sensitive to hyperpolarization (Hausser and Clark, 1997; Carter and Regehr, 2002; Suter and Jaeger, 2004). This would suggest that an inward current working in the subthreshold region could be interacting with IA to modify first-spike latencies. Previous work has shown that IA can interact with other subthreshold currents like a transient Ca2+ current (IT) and/or a hyperpolarizing activated cation current (IH) to influence neuronal excitability (Pape et al., 1994; MacLean et al., 2003; Russier et al., 2003).
In the present work, we studied the biophysical mechanisms that regulate first-spike latency in cerebellar stellate cells and its implications on synaptic integration and spike output. We found that stellate cells showed the longest delays to first spike within an intermediate membrane voltage range (-70 to -74 mV). Shorter delays to first spike occurred with more negative or positive membrane voltages. Thus, in stellate cells, the relationship between first-spike latency and membrane voltage is nonmonotonic. We also show that this property is established by the specific kinetic and steady-state properties of IA and IT. Moreover, these interactions have implications for synaptic integration, in that inhibitory inputs can change the probability for EPSPs to evoke spike discharge by shifting first-spike latency.
Materials and Methods
Animal care. Sprague Dawley rats [postnatal day 14 (P14) to P24] were obtained from Charles River (Saint-Constant, Quebec, Canada), and procedures were conducted according to guidelines approved by the local Animal Care Committee and the Canadian Council for Animal Care.
Electrophysiology. All chemicals were obtained from Sigma (St. Louis, MO) unless otherwise noted. Whole-cell current- and voltage-clamp recordings were obtained using an Axoclamp 2A and a 700A amplifier (Molecular Devices, Union City, CA), respectively, and data were collected with pClamp 8.1 software (Molecular Devices). Tissue slices of the cerebellum were prepared as detailed by McKay et al. (2005). Briefly, artificial CSF was composed of the following (in mm): 125 NaCl, 3.25 KCl, 1.5 CaCl2, 1.5 MgCl2, 25 NaHCO3, and 25 d-glucose preoxygenated by carbogen (95% O2, 5% CO2) gas. Rats were anesthetized with sodium pentobarbital (MTC Pharmaceuticals, Cambridge, Ontario, Canada), and the cerebellum was removed. Parasagittal tissue slices of 300 μm thickness were cut and subsequently transferred to the recording chamber of a Zeiss (Thornwood, NY) Axioskop FS-2 microscope and maintained at 35°C as a submerged preparation. Neurons were visualized using differential interference contrast optics and infrared light transmission. All recordings were performed in synaptic blockers that were bath applied after obtaining the initial seal: picrotoxin (50 μm), dl-2-amino-5-phosphonopentanoic acid (25 μm), 6,7-dinitroquinoxolinedione (10 μm; Tocris Cookson, Ellisville, MO), and (2S)-3-([(15)-1-(3,4-dichlorophenyl)ethyl]amino-2-hydroxypropyl) (phenylmethyl) phosphinic acid (1 μm; Tocris Cookson). Pipettes had a resistance of 3-6 MΩ with access resistance of 6-12 MΩ (80-90% compensation in voltage clamp). Our internal solution for current-clamp recordings consisted of the following (in mm): 130 K-gluconate, 0.1 EGTA, 10 HEPES, 7 NaCl, 0.3 MgCl2, 5 di-Tris-creatine phosphate, 2 Tris-ATP, and 0.5 Na-GTP, pH 7.3, with KOH. For all voltage-clamp recordings, TTX (100 nm) was bath applied. For voltage-clamp recordings isolating IA, the internal solution consisted of the following (in mm): 140 KCl, 5 EGTA, 10 HEPES, 2.5 MgCl2, with 1 NiCl2 and 2 CsCl in the external bath. For voltage-clamp recordings isolating IT, the internal solution consisted of the following (in mm): 100 CsCl, 10 EGTA, 10 HEPES, and 3 MgCl2, with 2 CsCl in the external bath.
Voltage- and current-clamp analysis. Voltage-clamp analysis consisted of fitting the time course of current inactivation with a single exponential decay function and the steady-state activation and inactivation curves with Boltzmann functions. Analysis of current-clamp data was accomplished using custom software written in MatLab 6.5 (MathWorks, Natick, MA) or IGOR PRO 4.07 (WaveMetrics, Oswego OR) software. Spike threshold was determined through analysis of the voltage derivative. A calculated junctional potential of 11 mV was subtracted from all current-clamp recordings. In voltage clamp, junction potential was not subtracted because the value was small (∼2 mV) for internal KCl and CsCl. Significant differences were assessed using one- or two-way ANOVA using Tukey's honest significant difference criterion. Averaged data are presented as mean ± SEM.
Simulations. Simulations were constructed in MatLab 6.5 using a fourth-order Runge-Kutta algorithm with a time step (dt) of 0.005 ms using the following equations: (1)
Na+ current activation. (2)
Na+ current inactivation. (3)
steady-state Na+ current inactivation. (4)
Ca2+ current activation. (5)
Ca2+ current inactivation. (6)
steady-state Ca2+ current inactivation. (7)
noninactivating K+ current activation. (8)
steady-state noninactivating K+ current activation. (9)
inactivating K+ current activation. (10)
inactivating K+ current inactivation. (11)
steady-state K+ current inactivation.
The time constant for inactivation of Na+ current in the model was voltage dependent and described by a Lorentzian function: (12)
for h: Vc =-74, w = 28, A = 232, y0 =-0.15.
Voltage was integrated according to the following: (13)
Constants in the model consisted of the following:
C = 1.5 μF/cm2; ENa+ = 45 mV; EK+ =-90 mV; Eleak =-70 mV; ECa2+ = 22 mV; gNa = 30 μS/cm2; gK = 7 μS/cm2; gleak = 0.1 μS/cm2; gA = 16 μS/cm2; gT = 0.55 μS/cm2.
The conductances for synaptic waveforms were modeled as an α function: gPSP(t) = te(-αt), with α = 3 for EPSP and α = 1.25 for IPSP and current equal to IE(t) = gEPSP(V) and II(t) = gIPSP(V + 82).
Results
Stellate cells show a nonmonotonic first-spike latency-voltage relationship
Stellate cells were recorded from the upper half of the cerebellar molecular layer and identified as small cells of 7-10 μm diameter with a multipolar soma that tapered to the primary trunk of small-diameter dendrites (n = 90). The input resistance of cells ranged between 500 and 700 MΩ, requiring bath perfusion of synaptic blockers to reduce background membrane fluctuations. In the presence of synaptic blockers, stellate cells spontaneously discharged Na+ spikes at a rate of 10-30 Hz. This differs from the activity seen in stellate cells in intact networks in vivo, in which spikes can discharge at a rate of only ∼3 Hz and then primarily in response to synaptic activity (Jorntell and Ekerot, 2003). This study then indicates that stellate cells ordinarily rest at a potential just hyperpolarized enough to allow single EPSPs to trigger spikes. To more closely approximate this activity and compare first-spike latencies, cells were set just below spike threshold using small amounts of bias current (10-20 pA) to establish a quiescent resting membrane ordinarily between -60 and -65 mV. Spike threshold for stellate cells was -58 ± 1.4 mV (n = 14). The resting potential was thus established at a level at which either evoked or simulated EPSPs readily evoked single spike discharge.
Previous immunohistochemical work has shown that stellate cells express Kv4.3 channels that are known to underlie fast-inactivating K+ currents (IA) (Hsu et al., 2003). To investigate the effects of IA on first-spike latency, we used hyperpolarizing pulses to vary its availability. Our standard current-clamp protocol consisted of a 500 ms hyperpolarization followed by a 300 ms depolarization (Fig. 1A). We found that stellate cells did not respond with an increase in firing frequency following a strong hyperpolarization to -90 mV (Fig. 1A). Interestingly, the cells did show a dramatic increase in first-spike latency from intermediate hyperpolarized membrane potentials of -70 to -74 mV (Fig. 1A, middle). This latency was reduced by using either a preceding hyperpolarization or depolarization (Fig. 1A, left and right). To gain a better understanding of the relationship between first-spike latency and membrane voltage, we systematically varied the hyperpolarizing step magnitude from -90 to -60 mV over the course of 30 current steps. Because the input resistance of stellate cells was very large, the size of the current steps was adjusted for each cell. We used the mean membrane voltage during the 25 ms preceding the depolarized step to obtain a more accurate estimate of the state of voltage-dependent currents.
Unlike previously described latency plots (Kanold and Manis, 1999; Shibata et al., 2000), our relationship between spike latency and the preceding membrane voltage was nonmonotonic, in which spike latency could be delayed or reduced in a voltage-dependent manner (Fig. 1B,C). To quantify the voltage-latency relationship, we created a scatter plot of spike latencies against preceding potentials for each cell (60-120 points per plot) (Fig. 1B). The voltage associated with the longest delay in first-spike latency was established within a sliding 4 mV bin width, identifying a mean voltage of -72 mV ± 1.01 (n = 12) (Fig. 1C). All other data points were then referenced to this mean value and binned at 4 mV increments. As seen in Figure 1C, first-spike latency was highly voltage dependent. At membrane voltages of approximately -72 mV, the cells had a maximum latency to first spike of 128 ± 15 ms, whereas at -80 mV, latency was 70 ± 10.8 ms, and at -60 mV, spike latency was 20 ± 4.1 ms (Fig. 1C) (n = 12). Additionally, the first-spike latency-voltage relationship was dependent on the basal level of depolarization and excitability of the cell. When the bias current was adjusted such that the cell fired tonically at low frequencies (1-12 Hz), the nonmonotonic first-spike latency relationship was retained (supplemental Fig. 1A, available at www.jneurosci.org as supplemental material). When the basal level of excitability was adjusted to allow for tonic high-frequency firing (>20 Hz), the nonmonotonic first-spike latency relationship was abolished (supplemental Fig. 1B, available at www.jneurosci.org as supplemental material). We note that in vivo recordings from stellate cells suggest that these cells most often discharge spontaneously in the low-frequency range in the presence of synaptic inputs and do not support regular tonic high-frequency firing (Ekerot and Jorntell, 2001; Jorntell and Ekerot, 2003). The higher basal level of excitability that allows for tonic firing of stellate cells may simply be a phenomenon of the in vitro preparation. The nonmonotonic voltage-latency relationship observed here is thus expected to be active under physiological conditions.
To fully assess the range of stellate cell output, we extended the protocol to include progressively larger depolarizing steps, there by assessing the effects of a preceding hyperpolarization to approximately -90 mV on a full frequency-current (F-I) plot of stellate cell output. This established that the F-I relationship was linear, with a firing range between 5 and 50 Hz that was not affected by the presence or absence of a preceding hyperpolarization (Fig. 1D). A reduction in first-spike latency from hyperpolarizing voltages might be expected to be accompanied by an increase in firing frequency after the generation of a rebound burst. This, however, did not occur in stellate cells, because firing frequency did not change significantly as a function of preceding membrane voltage (Fig. 1E) (p = 0.4; n = 12).
Long delays to first-spike latency are consistent with the presence of a fast-inactivating K+ current (Byrne et al., 1979; Storm, 1988). However, a reduced first-spike latency following strong hyperpolarization suggested the interaction of an additional inward subthreshold current, which could include IT or IH. To first test for the presence of IA in stellate cells, we bath applied 4-AP at 6 mm, a concentration shown to block a wide spectrum of A-type K+ currents (Coetzee et al., 1999), and performed the same protocol and analysis as in Figure 1, A-C. The application of 4-AP significantly decreased first-spike latency from controls and eliminated the nonmonotonic spike latency relationship (Fig. 2A,B) (p < 0.05). The latency values were also now uniform, ranging from 26 to 33 ms, and not significantly different from the -72 mV bin that originally distinguished the shift from short to long first-spike latencies (p = 0.78; n = 5) (Fig. 2B). Surprisingly, 4-AP still did not unmask a rebound burst, because firing frequency after the hyperpolarization remained independent of the preceding membrane voltage (data not shown). However, the application of 4-AP did increase baseline firing frequency, suggesting that IA contributed a window current sufficiently depolarized to influence baseline firing frequency. The ability to fire spikes after a hyperpolarization in the presence of 4-AP was also transient, with spike firing terminating before the end of the excitatory current pulse (Fig. 2A). We reasoned that this was because of the very high frequency of spike firing (>100 Hz), which could elicit a large hyperpolarization through Ca2+-activated K+ currents. This indicated that the 4-AP-sensitive current had some baseline contribution to steady-state firing frequency.
We also considered that 4-AP may block other nontransient K+ currents, which could change our interpretation of the above results. To test this, we voltage clamped stellate cells and separated the inactivating and noninactivating K+ currents (Fig. 3) (see below for additional details). We found that 6 mm 4-AP produced a substantial block of a fast-inactivating K+ current (35 ± 4.6% reduction; p < 0.05; n = 5) but did not significantly block noninactivating K+ currents (8 ± 3.5%; p > 0.05). Thus, the effects of 6 mm 4-AP on first-spike latency in stellate cells can be primarily attributed to IA.
Although there are no well established specific blockers for IT, we attempted to block IT under current-clamp conditions using mibefradil at 500 nm, which has been shown to specifically block IT in Purkinje cells (McDonough and Bean, 1998). Mibefradil eliminated the nonmonotonic first-spike latency characteristic of stellate cells, resulting in a graded increase in first-spike latency with the magnitude of a preceding hyperpolarization (Fig. 2C,D). First-spike latency values now ranged from 187 ± 17 ms with a preceding membrane voltage of -92 mV to 30 ± 3 ms with a preceding membrane voltage of -60 mV. Furthermore, the first-spike latencies generated with preceding membrane voltages from -92 to -84 mV were significantly different from those generated in the control (p < 0.05; n = 5). This result is consistent with mibefradil blocking IT, which would isolate IA in the subthreshold region and lead to membrane hyperpolarizations inducing a more conventional monotonic increase in first-spike latency.
It is well established that a hyperpolarization-activated cationic current (IH) can be recruited at negative voltages to increase cell excitability (Pape, 1996). Such a profile of recruitment could then help account for the shorter latencies observed after strong hyperpolarizations. To test whether IH was involved in setting the first-spike latency relationship in stellate cells, we bath applied 4-(N-ethyl-N-phenylamino)-1,2-dimethyl-6-(methylamino)-pyrimidinium chloride (ZD 7288; 10 μm) (Fig. 2E). Application of ZD 7288 uniformly increased the first-spike latency at all voltages by 18.9 ± 3.6 ms (n = 11) and eliminated the characteristic hyperpolarization activated membrane voltage sag associated with IH. It did not, however, alter the nonmonotonic first-spike latency-voltage relationship when compared with control data (p = 0.87). Furthermore, the preceding membrane voltage that produced the longest first-spike latency did not significantly change (control, -72 ± 1.1 mV vs ZD 7288, -74 ± 1.9 mV; p = 0.19) (Fig. 2E). Thus, IH can uniformly modify first-spike latency times but is not involved in setting the nonmonotonic first-spike latency profile of stellate cells.
These data establish that first-spike latency of stellate cells is determined in large part through an interplay between IA and IT. These currents share the property of exhibiting overlapping activation and inactivation curves that establish window currents to precisely regulate their availability over a voltage range near resting potential. The nonmonotonic relationship between membrane voltage and first-spike latency must then reflect the balance between IA and IT activation according to their specific voltage-dependent kinetic and steady-state properties. A key prediction to explain the decrease in spike latency found with prestep voltages more negative than -72 mV would be a relatively greater deinactivation of IT relative to IA. Thus, to further investigate the role of IA and IT in producing first-spike latency characteristics, we performed voltage-clamp analysis and then modeling to test whether differences in channel kinetics could account for the experimental data.
Biophysical characterization of IA and ITB
To confirm the presence of IA and IT and determine their biophysical parameters, we used whole-cell voltage clamp and isolated the currents pharmacologically and electrophysiologically. To isolate IA, we bath applied TTX (100 nm), Ni2+ (1 mm), and Cs+ (2 mm) or ZD 7288 (10 μm). In all voltage-clamp recordings, synaptic blockers were used to increase cell input resistance and reduce spontaneous activity that contaminated current traces. These experiments revealed two classes of K+ currents, a low-threshold and fast-inactivating current consistent with IA and a higher-threshold and noninactivating delayed rectifier. IA was insensitive to 20 mm tetraethylammonium (TEA) but blocked with 4-AP at concentrations >4 mm (35% block at 6 mm) and completely blocked with 20 mm 4-AP. However, we found that long-term recordings with 20 mm 4-AP had negative effects on cell health. Thus, to isolate IA, we used a 1 s depolarizing voltage step to -45 mV, which inactivated >95% of IA, and subtracted this record from control traces that used a prevoltage step of -100 mV. Our voltage-clamp protocol consisted of step depolarizations from -100 to 30 mV in 10 mV increments and were followed by a test pulse to 0 mV to determine the steady-state inactivation curve. IA started activating between -70 and -60 mV with peak current ranging from 500 pA to 5 nA (at 0 mV) and showed rapid inactivation (Fig. 3A). Activation was very fast, with time to peak ranging from 0.8 to 2.5 ms. The conductance-voltage relationship was fit with a Boltzmann function and had a V1/2 of -27 ± 3.0 mV and a slope factor of 10.4 ± 1.1 (n = 7). The steady-state inactivation curve had a V1/2 of -68.2 ± 1.2 mV and a slope factor of 5.1 ± 0.5 mV (n = 9) (Fig. 3B). The combination of these factors produced a small window current for IA between -60 and -40 mV, which, in conjunction with the high density of expression in stellate cells, explains why the application of 6 mm 4-AP increased baseline firing frequency, as shown in Figure 2A.
To isolate IT, we used TTX and Cs+ in the external medium to block INa and IH, respectively, and used intracellular Cs+ to block K+ currents. Our voltage-clamp protocol consisted of step depolarizations from -120 to -30 mV in 10 mV increments. Under these conditions, we observed a low-threshold and fast-inactivating Ca2+ current that activated at approximately -80 mV (Fig. 3A,C) and was insensitive to bath application of the dihydropyridine nifedipine (1 μm). Peak IT in the I-V plot was seen between -50 and -40 mV and ranged from -20 to -80 pA. Like IA, IT showed rapid activation kinetics with time to peak ranging from 1.7 to 3.8 ms. At voltages more positive than -30 mV, a second slower and noninactivating Ca2+ current was observed in most cases. Because of the large discrepancy in the activation voltage between the two Ca2+ currents, we could restrict our voltage-clamp protocol to the lower-voltage range to activate IT in isolation. The activation- and inactivation-voltage relationships were fit with a Boltzmann function. The activation curve had a V1/2 of -57 ± 3.3 mV and a slope factor of 5.3 ± 1.1 (n = 6), whereas the inactivation curve had a V1/2 of -77.6 ± 0.73 mV and a slope factor of 4.7 ± 0.4 (n = 7). Note that the steady-state inactivation curve for IT was negatively shifted compared with that of IA by ∼10 mV (Fig. 3D) (IA, -68.2 ± 1.2 vs IT, -77.6 ± 0.73 mV; p < 0.05) and that the conductance-voltage relationship was more negative and steeper for IT than IA (Fig. 3E). Thus, a hyperpolarizing current pulse causing a membrane voltage shift more negative than -72 mV recruits more IT relative to IA and as a result may decrease first-spike latency.
We also characterized the inactivation time constant-voltage relationship for both IA and IT, because this parameter could play an important role in influencing the latency to first spike in cells (Storm, 1988). With both currents, the inactivation profile was accurately fit with a single exponential decay function. The inactivation time constants for IA showed weak voltage dependence with values ranging from 13.2 ms at -60 mV to 8.8 ms at -10 mV (Fig. 3F). Conversely, the inactivation time constants for IT showed a slightly stronger voltage dependence with values ranging from 12.7 ms at -70 mV to 5.2 ms at -40 mV (Fig. 3F). Note, however, that the inactivation time constants for both currents approached one another (9-16 ms) at hyperpolarized potentials and were not significantly different at -60 mV (Fig. 3F). Thus, at voltages just below the firing threshold for stellate cells (approximately -65 to -70 mV), both currents inactivated with a similar time course.
The kinetics obtained in the voltage-clamp experiments lead to predictions on how IA and IT interact to produce the first-spike latency characteristics of stellate cells. The key factor in this process is expected to be the negatively shifted steady-state inactivation curve of IT relative to IA. Two other contributing factors are expected to be the similar rates of IA and IT inactivation near spike threshold and the relative density of IA versus IT. These properties would establish three key regions of membrane potential that could differentially modulate first-spike latency. The first region is spike firing immediately after a strong hyperpolarization (-90 to -74 mV), in which sufficient IT can be deinactivated to oppose IA. In the second region (-74 to -70 mV), an intermediate hyperpolarization significantly inactivates more IT than IA, thereby causing a long first-spike latency. Finally, after a weak or no hyperpolarization (less than -70 mV), both IA and IT are inactivated, leading once again to a decrease in spike latency in direct relation to the membrane depolarization. By combining both IA and IT, the stellate cell can achieve a nonmonotonic relationship between the first-spike latency and the preceding membrane voltage.
A firing model using IA and IT reproduces stellate cell first-spike latency characteristics
To further test our hypothesis that IA and IT can synergistically determine the first-spike latency characteristics of stellate cells, we built a reduced firing model incorporating their biophysical parameters. To simplify the model, we assumed that activation of both currents equilibrated with voltage instantaneously, which is reasonable considering that both currents activated very rapidly. Our model also incorporated a high ratio of IA relative to IT (IA, 16 μS/cm2; IT, 0.55 μS/cm2), which was within the values observed experimentally (gmax for IA, 1 to 0.17 vs gmax for IT, 0.1 to 0.01; values normalized to IA gmax). Finally, the inactivation time constant was fixed at 15 ms, because the region of interest and the principal function for IA and IT is in the subthreshold region in which both currents showed similar inactivation time constants (Fig. 3F) (both ranged from 9 to 16 ms at -60 mV). In addition to IA and IT, the model incorporated a standard fast inactivating Na+ and noninactivating K+ current, the parameters of which were set to reproduce the spike shape of stellate cells and the firing frequency range over a F-I plot (5-50 Hz).
The model had a current threshold for spiking of 0.83 μA/cm2, which corresponded to a voltage threshold of approximately -64 mV, which compared favorably with the threshold determined in stellate cells of -58 ± 1.4 mV (n = 14; range, -51 to -68 mV). When the model was depolarized to 0.9 μA/cm2 following a negative current injection of -2 μA/cm2 it produced a spike latency of 74 ms (Fig. 4A, left). Using the same protocol but without negative current injection (0 μA/cm2), the model generated a longer latency of 122 ms (Fig. 4A). Finally, following a positive current injection of 0.8 μA/cm2, the model had a short latency with a value of 72 ms (Fig. 4A). The F-I plot of the model was also similar to that of the stellate cell, with a relatively linear F-I relationship that was insensitive to prehyperpolarization (-2 μA/cm2) (Fig. 4B).
We quantified the first-spike latency-voltage relationship in the model by systematically applying a range of current injections that produced preceding membrane voltages between -90 and -64 mV. In this way, we could assess the accuracy of the model with respect to the experimental results of Figure 1C. In the presence of both IA and IT, we found that the model reproduced the experimental data, with a maximum delay of 130 ms following a membrane voltage of -74 mV and delays of 74 and 20 ms following voltages of -85 and -65 mV, respectively. In the absence of IT, the influence of IA in the model produced long first-spike latencies that increased in a monotonic manner with membrane hyperpolarization (Fig. 4C) and was thus similar to the experimental results in the presence of mibefradil (Fig. 2D). Conversely, the removal of IA in the model increased firing frequency and almost entirely eliminated delays to first spike, with values ranging from 22 to 10 ms following membrane voltages from -90 to -68 mV, respectively (Fig. 4C). Finally, the removal of both currents produced a slight increase in first-spike latency compared with IT expression alone.
We note that the removal of IA did not promote a rebound burst discharge in the model cell (data not shown), a result that matched experimental data in the presence of 6 mm 4-AP (Fig. 2A,B). The lack of rebound burst discharge in the absence of IA was because of two reasons. First, the stellate cell and model had relatively small amounts of IT that inactivated very rapidly. Thus, IT is sufficient to influence latency near spike threshold in which small and rapid changes in net current can substantially change the first-spike latency but is insufficient to produce an increase in firing frequency. Second, the removal of IA significantly increased firing frequency both in the model and stellate cells (>100 Hz), such that the system was operating near the maximum possible firing frequency even with the small depolarizing current injection presented after the hyperpolarization. Consequently, the recruitment of IT with a strong hyperpolarization could not increase firing frequency, because the cell was in the saturated region of the F-I relationship.
To test the hypothesis that the difference in steady-state inactivation curves between IA and IT was responsible for the first-spike latency-voltage relationship, we used the model to shift the V1/2 inactivation values. In the first test, we shifted the inactivation V1/2 of IA to match that of IT, whereas in the second test, we shifted the inactivation V1/2 of IT to match that of IA (Fig. 4D). Although there was a small shift in the firing threshold for the second test (because of the larger window current for IT), matching the V1/2 of both currents prevented the longest first-spike latency from occurring at intermediate voltages and turned the latency-voltage relationship into a monotonically increasing function (Fig. 4D). This is because the cell expressed substantially larger amounts of IA than IT. Matching the V1/2 inactivation values thus removed any net increase in inward current that would result from hyperpolarizing the model cell more negative than -75 mV, resulting in a net outward K+ current. This result depended on the ratio of IA to IT, such that reversing the relative density of these currents under conditions of identical V1/2 for inactivation produced a bursting cell (data not shown). Finally, the similar rates of inactivation were necessary to allow IT to offset the influence of IA in delaying spike onset.
To test the above ideas, we measured the difference in peak activation between IA and IT during the conditioning hyperpolarizing step and plotted it against the membrane voltage during the step. For comparison, we have provided the steady-state inactivation curves for IA and IT used in the model and indicated the key voltage range of interest by dashed lines (Fig. 5A). As shown in Figure 5B, the maximum difference between the variables controlling IA (hA) and IT (hT) inactivation occurred at approximately -74 mV, which also corresponded to the longest spike latency in the model (Fig. 4C). Moreover, by measuring the amount of the h variable for IA and IT recruited per membrane voltage change in the conditioning prestep, we could quantify the sensitivity of each h variable to membrane voltage (dh/dV). With this analysis, a greater value indicates a greater sensitivity to membrane voltage, such that hyperpolarization greatly increases the net amount of available current. This test revealed that the dh/dV curve for IA started decreasing at more positive membrane voltages than IT (Fig. 5C). Thus, IT deinactivation continued to grow with more negative membrane voltages, whereas IA began saturating (Fig. 5 A, C). Note that at intermediate voltages, the dh/dV value for IA is significantly greater than IT, which corresponded to the region of longest spike latencies in the model (Fig. 5C).
Our results show that first-spike latency characteristics of cerebellar stellate cells can be accounted for by the specific kinetic and steady-state properties of IA and IT. The most critical factor is a left-shifted inactivation voltage for IT compared with IA, with additional influence of a similar inactivation time constant for both currents near resting potential, and the relative density of expression.
Spike latency characteristics of stellate cells determine responses to coincident synaptic input
We next sought to test the interaction between IA and IT to consider its implications for synaptic integration. Stellate cells receive excitation from parallel fibers and inhibition from other stellate cells (Palay and Chan-Palay, 1974). In vitro dynamic-clamp experiments further indicate that stellate cells are optimized to process fast transient inputs (Suter and Jaeger, 2004), suggesting that first-spike latency effects may contribute significantly to modulating synaptically triggered spike firing in vivo. The effects of a nonmonotonic spike latency-voltage relationship would conceivably be most readily detected when EPSPs are immediately preceded by IPSPs. A study by Jorntell and Ekerot (2003) has shown, through in vivo patch recordings, that the baseline membrane potential of stellate cells can be shaped by EPSP and IPSP sequences with no apparent temporal pattern. In contrast, recent work has proposed that a feedforward inhibitory network can restrict synaptic activation of stellate cells to an EPSP-IPSP temporal sequence (Mittmann et al., 2005). It is unlikely that the mechanism we describe here would function within the context of feedforward inhibition. However, a randomly activated series of synaptic potentials could well act to trigger specific delays under conditions of an IPSP-EPSP sequence. We hypothesized that the nonmonotonic spike latency-voltage relationship could further alter the probability for spike firing after coincident inhibition, because this would make membrane voltage immediately before excitation more negative and thus reduce first-spike latency. To test this, we examined the effects of near coincident IPSP-EPSP responses in the model and in stellate cell recordings on the ability to trigger a spike.
Previous work has determined the kinetics of stellate cell EPSPs and IPSPs in vivo and has shown that there is a lack of heterogeneity of either EPSPs or IPSPs in terms of kinetics and peak amplitudes (Jorntell and Ekerot, 2003). The reversal potential of AMPA-mediated EPSPs is 0 mV, whereas that of GABAA-mediated IPSPs is -82.5 mV (Carter and Regehr, 2002), with EPSPs exhibiting faster kinetics than those of IPSPs (Jorntell and Ekerot, 2003). To reproduce the EPSP and IPSP waveforms, we injected conductance waveforms into the model as α-functions with parameters selected to match the in vivo characteristics reported by Jorntell and Ekerot (2003) (Fig. 6A). We then applied the EPSP in conjunction with different sizes of IPSPs to the model cell to examine the effects of coincident synaptic inputs on first-spike latency.
In the model, the EPSP presented alone attained a level of depolarization that was sufficient to reliably generate a single spike (Fig. 6C). We next applied the IPSP and EPSP in combination with the IPSP (peak voltage, -73 mV) triggered 15 ms before the EPSP. When coupled with a preceding IPSP, the EPSP failed to generate a spike and simply depolarized the membrane voltage by ∼10 mV (Fig. 6C). The inhibition of firing was not simply attributable to an algebraic summation of the IPSP and EPSP reducing overall EPSP amplitude, because the EPSP was reduced by <0.5 mV. Moreover, a more slowly decaying EPSP of the same amplitude was sufficient to generate a spike when presented in conjunction with the IPSP (data not shown). We reasoned that the block of spike firing was attributable to an IA-mediated delay to first spike associated with the hyperpolarization evoked by the single IPSP. In fact, the hA variable was increased by the IPSP from 0.068 to 0.71. When we doubled the conductance underlying the IPSP to simulate two coincident inhibitory events (peak voltage, -80 mV), the ability of the EPSP to generate a spike was restored (Fig. 6C). Note that this occurred despite a small decrease in EPSP peak amplitude but was associated with an increase in the hT variable from 0.205 to 0.29. Thus, the nonmonotonic latency-voltage relationship in stellate cells allows coincident IPSPs to promote spike firing by an EPSP by removing a delay to spike firing that otherwise dominates with a single IPSP.
To investigate this in the stellate cell, we applied the same conductance waveforms shown in Figure 6A as current injection protocols (we reversed the inhibitory conductance to simulate inhibition) in the presence of synaptic blockers and adjusted the peak current injection values to produce similar IPSPs and EPSPs to those in the model. Although shunting could have an additional influence, we reasoned that the relatively negative Cl- reversal potential of -82.5 mV would reduce this problem. Furthermore, Suter and Jaeger (2004) showed that in stellate cells, conductance injection with a dynamic-clamp setup did not produce substantially different results from direct current injection. As shown in Figure 6D, stellate cells were able to reproduce the model results, with single IPSP waveforms inhibiting spike firing by the EPSP, and coincident inhibition restoring spike firing. We tested this in seven cells and found that it could be evoked in all cells, although the reliability varied from 30 to 50%. In the remaining cases, IPSPs were often more effective in inhibiting spike firing. Part of the variability was associated with the difficulties in maintaining a constant membrane voltage throughout the duration of the IPSP-EPSP waveform in a cell with such a high input resistance. Nevertheless, these experiments are sufficient to establish the potential for a nonmonotonic first-spike latency relationship to affect stellate cell output in response to membrane perturbations relevant to afferent input in vivo.
Discussion
The fast-inactivating class of K+ channels is known to be involved in actively generating a delay to first spike by activating in the subthreshold region after hyperpolarization (Byrne et al., 1979; Storm, 1988; Kanold and Manis, 1999). Conversely, IT has been generally implicated in generating rebound burst firing that is also promoted by a preceding hyperpolarization (Llinas and Jahnsen, 1982; Huguenard and Prince, 1992; Huguenard, 1996). We have shown that the coexpression of these two channels in stellate cells gives rise to a novel spike latency-voltage relationship, allowing preceding membrane hyperpolarizations to either reduce or lengthen the latency to first spike. Furthermore, we show that these interactions can be invoked by synaptic inputs inherent to stellate cells, indicating that this interaction will be effective in modulating the latency to discharge of stellate cells in vivo and hence their influence on Purkinje cell output.
Recordings and modeling indicate that first-spike latency in stellate cells is shaped by the coexpression, kinetic properties, and relative densities of IA and IT. The presence of IA allows moderate hyperpolarizations (more positive than -72 mV) to invoke a well recognized increase in the delay to first spike. A nonmonotonic spike latency-voltage relationship arises as the result of IT having a more negative steady-state inactivation curve relative to IA. Thus, spike latencies after strong hyperpolarizing steps (more negative than -72 mV) can disproportionately recruit more IT than IA because of the negatively shifted inactivation curve of IT. This results in a net increase in inward current after depolarization that translates into a significant reduction in first-spike latency. This further illustrates that IT can interact with other currents in the subthreshold region to function in roles outside of generating rebound burst firing. Finally, our model and recordings demonstrate that a nonmonotonic spike latency-voltage relationship can have important implications for synaptic integration. In particular, the data indicate that even subtle increases in inhibitory coincidence can essentially reverse first-spike latencies and change the probability of firing in response to excitatory synaptic depolarizations.
Low density of expression and fast inactivation kinetics of IT preclude burst firing
Stellate cells did not show rebound burst firing under physiological conditions. Even in the absence of IA, neither the cell nor model were able to produce rebound burst firing after strong hyperpolarizations. This result differs from a previous study that showed that the removal of IA in thalamic interneurons, which also coexpress IT, could lead to rebound burst firing (Pape et al., 1994). IT in stellate cells, however, differed from those recorded in thalamic cells in a few important aspects. First, stellate cells express relatively low levels of IT, especially when compared with IA. In thalamic interneurons, this ratio is more balanced (Pape et al., 1994). Second, IT in stellate cells is faster to inactivate near firing threshold compared with thalamic neurons (Pape et al., 1994). Thus, the absence of rebound burst firing in stellate cells is likely a result of the inability of a small and rapidly inactivating IT to evoke an increase in firing frequency. Furthermore, in the absence of IA, the cell had an intrinsic high firing frequency that was operating near the saturated region of the F-I relationship. As a result, a small additional recruitment of IT may not be able to significantly influence firing frequency resulting from the low gain (slope) in the F-I relationship of the neuron in 4-AP.
Potential molecular identity of IA and IT
We found that IA in stellate cells was insensitive to TEA at concentrations of up to 20 mm but was sensitive to 4-AP (>4 mm). The lack of effect of TEA rules out the potential contribution of a number of A-type or inactivating K+ currents of the Kv1 or Kv3 family (Coetzee et al., 1999). Rather, these results are consistent with members of the Kv4 family, with stellate cells reportedly positive for at least Kv4.3 immunolabel (Hsu et al., 2003). The biophysical parameters for IA in stellate cells are also similar to those for Kv4 channels in other central neurons as well as expression systems when Kv4 subunits are coexpressed with the intracellular protein DPPX (dipeptidyl aminopeptidase-like protein) (Nadal et al., 2003). The activation V1/2 of Kv4 subunits are relatively positive (V1/2 more than -30 mV), with shallow slope factors (k > 8) and a relatively negative and steep steady-state inactivation (V1/2 less than -65 mV; k <6) (Nadal et al., 2003). This compares very favorably to an activation V1/2 of -27 mV (k = 10) and inactivation V1/2 of -68 mV (k = 5) for IA recorded in stellate cells. The potential contribution of different Kv4 subunits in producing this current are not currently known.
Previous in situ hybridization studies have not detected measurable levels of mRNA for Cav3.x calcium channel subtypes in cerebellar stellate cells (Craig et al., 1999; Kase et al., 1999; Talley et al., 1999). Nevertheless, our recordings of a nifedipine-insensitive, low-threshold, and fast-inactivating calcium current that is blocked by 500 nm mibefradil provides strong support for its identification as a current generated by channels of the Cav3 family. Indeed, we now have immunocytochemical evidence for the expression of Cav3.2 and Cav3.3 channel subtypes in cerebellar stellate cells (M. Molineux, unpublished observations). Our measurements of IT are also consistent with the biophysical parameters measured in some central neurons, such as thalamic relay neurons (Huguenard and Prince, 1992; McCormick and Huguenard, 1992; Huguenard, 1996). One notable difference is the slightly faster inactivation kinetics measured in stellate cells. It is difficult, however, to compare Cav3 current kinetics recorded in different neuronal cell types or in expression systems given the lack of a full comparison of Cav3 channel distributions and the potential for posttranslational modifications to contribute to the wide variety of biophysical parameters reported for IT (Huguenard, 1996; Perez-Reyes et al., 1998; McRory et al., 2001; Chemin et al., 2002).
A single IPSP inhibits firing, and coincident IPSPs promote firing
The activation of IA is traditionally thought to delay first-spike latency, a result that can be induced with synaptic inhibition (Berman and Maler, 1998; Kanold and Manis, 2005). Because the first-spike latency-voltage relationship in stellate cells is nonmonotonic, our model and data suggest that interactions between EPSPs and IPSPs can be more complicated. We show that EPSPs following weak inhibition by single IPSPs can be effectively filtered to prevent spike firing by an IA-mediated delay in first-spike latency. However, we also show that stronger inhibition, which can result from coincident IPSPs, can overcome the filtering properties of IA by recruiting more IT. This mechanism is particularly suited in stellate cells because of the relatively negative Cl- reversal potential of -82.5 mV (Carter and Regehr, 2002). The decrease in first-spike latency observed here occurs with premembrane voltages more negative than -72 mV, which minimally requires that IPSPs reverse slightly more negative than this value. The control over spike latency by this mechanism proves to be substantial, in providing a range of mean latencies that span from 20 to 128 ms after a hyperpolarization. Previous studies reveal that even a single spike in stellate cells can have an important influence on Purkinje cells by delaying spike onset and reducing firing frequency (Hausser and Clark, 1997). The spike latency-voltage relationships presented here thus have significant potential to control cerebellar output. Given the complexity of cerebellar circuitry, the precise contribution that these interactions could have on cerebellar computation is not yet understood, but the current mechanism is likely to be a factor in establishing the temporal requirements expected for cerebellar output.
A previous study looking at parallel fiber activation and EPSP dynamics in Purkinje and stellate cells found that feedforward inhibition in stellate cells strongly influenced the precise spike timing in Purkinje cells (Mittmann et al., 2005). In that study, the EPSP in stellate cells preceded the IPSP. Hence, our mechanism of synaptic control of first-spike latency is unlikely to work under those conditions. We suggest that the mechanism of synaptic control of first-spike latency outlined in the present study will emerge during the activation of stellate cell receptive fields. In vivo recordings have shown conditions in which EPSPs and IPSPs may arrive in any sequence (Jorntell and Ekerot, 2003). The network activity that generates this form of synaptic response has not been established but could arise from the high degree of interconnectedness of stellate cells combined with parallel fiber activation (Kondo and Marty, 1998). In addition, the first-spike latency mechanism here could potentially modulate a feedforward inhibitory loop by affecting the precise timing of inhibition between stellate cells. The network activity that could take advantage of a nonmonotonic spike latency-voltage relationship thus remains to be determined. Nevertheless, our results identify a novel interaction between two transient currents that serves to modulate both the latency and probability for discharging a spike in stellate cells by excitatory synaptic inputs.
Footnotes
This work was supported by the Canadian Institutes of Health Research (CIHR) (R.W.T.), CIHR Doctoral Studentships (F.R.F., W.H.M.), Province of Alberta Graduate Scholarships (F.R.F., W.H.M.), and an Alberta Heritage Foundation for Medical Research Scientist Award (R.W.T.). We gratefully acknowledge the technical assistance of M. Kruskic and B. E. McKay for comments on this manuscript.
Correspondence should be addressed to Ray W. Turner, Hotchkiss Brain Institute, University of Calgary, 3330 Hospital Drive Northwest, Calgary, Alberta, Canada T2N 4N1. E-mail: rwturner{at}ucalgary.ca.
Copyright © 2005 Society for Neuroscience 0270-6474/05/2510863-11$15.00/0
↵* M.L.M. and F.R.F. contributed equally to this work.