Spike-Frequency Adaptation

Transcription

1 Spike-Frequency daptation Jan Benda and Joel Tabak Encyclopedia of Computational Neuroscience, Springer 214 Definition When stimulated with a constant stimulus, many neurons initially respond with a high spike frequency that then decays down to a lower steady-state frequency (Fig. 1 ). This dynamics of the spike frequency response is referred to as spike-frequency adaptation. Spike-frequency adaptation is a process that is slower than the dynamics of action potential generation. Spikefrequency adaptation by this definition is an aspect of the neuron s super-threshold firing regime, although the mechanisms causing spike-frequency adaptation could also be at work in the neuron s subthreshold regime. Spike frequency [Hz] f f B Spike frequency [Hz] f (I) f (I) Time [ms] Input I [n] Figure 1: Phenomenon spike-frequency adaptation. Spike train (vertical strokes at the top) of an adapting neuron evoked by the onset of a constant stimulus (green bar at the bottom, I = 15 ns). The spike frequency drops from an initially high onset rate, f, in an approximately exponential way down to a lower steady state rate, f. B Repeating the experiment shown in for various stimulus intensities I (here an injected current) results in the onset f -I curve f (I) and the steady-state f -I curve f (I). Shown are simulations of a single-compartment conductance-based model with an M-type current (the Ermentrout model, Ermentrout, 1998; Benda et al., 21, with ḡ Ca =, ḡ HP =, and ḡ M = 8 µs/cm 2, the input current is scaled such that 1 n equals 1 µ/cm 2 ). 1

2 J. Benda & J. Tabak: Spike-frequency adaptation 2 Detailed Description In the context of spike-frequency adaptation, spike frequency is often measured as the instantaneous rate, i.e., the averaged reciprocal interspike-intervals at each time. This measures the inverse period of super-threshold firing. In contrast, the PSTH (peristimulus time histogram) estimates the probability of a spike at a given time and therefore is also sensitive to variability of the response caused by intrinsic or external noise sources. Important characteristics of a neuron with spike-frequency adaptation are its f -I curves (spike frequency f versus input strength I). Because of adaptation, a certain input intensity will not evoke a constant spike frequency. Instead, the initial response in spike frequency to the onset of a constant input will result in the so-called onset f -I curve (green curve in Fig. 1 B). Later during the constant stimuli, the spike frequency reaches a steady state. Drawn as a function of input intensity, this forms the steady-state f -I curve, which usually lies below the onset f -I curve, because the spike frequency adapts from the onset down to the steady state (red curve in Fig. 1 B). The third type of f -I curves is the adapted f -I curves (Benda and Herz, 23). For these the neuron is first adapted to some adaptation input of fixed intensity and then the onset responses to various test inputs are measured (Fig. 2). The adapted response curves indicate how the neuron being in a certain and fixed adaptation state will initially respond to new stimuli. While each neuron has exactly one onset and one steady-state f -I curve, it has for each adaptation stimulus a different adapted f -I curve. Spike frequency [Hz] Input I f (I ) I f (I ) f (2n,I ) f (1n,I ) Time [ms] I B Spike frequency [Hz] f (I) Input I [n] f (I,I ) Figure 2: dapted f -I curves. In addition to the onset and the steady-state f -I curves, the adapted f -I curves provide important information about the adaptation properties of a neuron. In order to measure an adapted f -I curve for a specific adaptation state, the neuron is first adapted to its steadystate response f (I ) by stimulating it with a constant adapting stimulus I (here 15 n). Then, a test stimulus I is applied and the onset response f (I,I ) to this stimulus is measured. In the figure, the responses (red and orange) to two different test stimuli (I = 2 nin green and I = 1 nin dark cyan) are shown. B Each adapting stimulus I results in one adapted f -I curve f (I,I ). Shown are adapted f -I curves for I = 5, 1, 15, 2 n. In this example the adaptation process mainly shifts the adapted f -I curves to the right. ll curves together characterize the spike-frequency response of an adapting neuron. Same model neuron as in Fig. 1. Simple adaptation currents (see below) will shift the neuron s adapted f -I curves to higher input intensities. Other (less well-understood) mechanisms like presynaptic inhibition may also f (I)

3 J. Benda & J. Tabak: Spike-frequency adaptation 3 tilt the adapted f -I curves (Hildebrandt et al., 211) or distort it in other ways (see also Fig. 5). Spike-frequency adaptation is related to the classification of neurons into phasic, phasictonic, or tonic spiking neurons. Tonic spiking neurons are neurons that keep firing in response to a constant stimulus with an almost constant rate. Many neurons that are classified as tonic spiking, however, still show some weak spike-frequency adaptation. On the other hand, phasic spiking neurons only emit one or several spikes to the onset of the stimulus and then cease spiking. If this onset response consists of several spikes with a decreasing spike frequency, then the phasic response is probably caused by strong spike-frequency adaptation. With less or only a single spike and no clear decrease in spike frequency, this type of response can be that of a type III neuron. The intermediate phasic-tonic responses exactly match the definition of spike-frequency adaptation. Mechanisms There are many different mechanisms that all cause a neuron to adapt its spike-frequency response. First, there are adaptation currents ionic currents that act together with the spikegenerating currents ( output-driven adaptation ). Second, other mechanisms, such as fatigue of receptor currents, adapt the input a neuron receives ( input-driven adaptation ). nd third, there are network effects that also can adapt the neuron s response in both an input-driven or output-driven way. daptation Currents There exists a whole zoo of mechanisms based on ionic currents that are all directly or indirectly activated by the action potentials generated by the neuron the output of the neuron and have at least one time scale involved that is slower than the dynamics of an action potential. M-type currents are a simple example (Fig. 3, Brown and dams, 198). These are voltagegated potassium currents that are mainly activated at high voltages as they occur during an action potential and that deactivate with a slow time constant of about 1 ms. If the neuron initially fires with a frequency of say 2 Hz, then the activation of the M-type currents builds up from one action potential to the next, since there is not enough time for the current to completely deactivate between succeeding action potentials. s a potassium current, the activated M-type currents then counteract the input to the neuron and as a consequence the spike frequency is reduced the neural response adapts. n important class of adaptation currents are HP-type currents, calcium-gated potassium currents, that enhance the afterhyperpolarization following an action potential (Sah, 1996). Here, during each action potential, calcium enters the cell through voltage-gated calcium channels. Depending on the intracellular calcium concentration, the HP-type currents are activated and as potassium currents inhibit the input to the neuron. Calcium is only slowly removed from the cytosol and this time scale is then reflected by the resulting spike-frequency adaptation. Further ionic mechanisms that cause spike-frequency adaptation include sodium-activated potassium channels and slow inactivation of the sodium current (for details see Ion Channels and Pumps in Neuronal daptation ).

4 J. Benda & J. Tabak: Spike-frequency adaptation 4 5 V [mv] B a C IM [n] Time [ms] Figure 3: Dynamics of an adaptation current. The membrane potential V in response to a step input (I = 15 n, indicated by the green bar) simulated with the same model as in Fig. 1. s adaptation progresses the interspike intervals get longer. B The gating variable of the M-type current a is increased by every action potential and decays slowly in between them. This way a slowly builds up during stimulation. C The adaptation current, here an M-type current, I M, gets larger and larger as a increases and inhibits the input. This way it causes the neuron s response to adapt. Input daptation In receptor neurons the whole transduction process from the physical stimulus to the receptor current can adapt in different ways. These can be mechanisms controlling the sensitivity of the receptor organ, like the pupillary light reflex. Within the receptor neuron itself, the transduction machinery may produce adapting receptor currents due to fatigue (e.g., bleaching of rhodopsin) or active adaptation to the intensity of the physical stimulus within a second messenger cascade (e.g., in photoreceptors) or by molecular motors (e.g., in mechanical transduction in hair cells). ll these and many other mechanisms eventually produce an adapting receptor current that then leads to spike-frequency adaptation in the receptor neuron. Gollisch and Herz (24) demonstrated how to unmask input-driven adaptation from output-driven adaptation in auditory receptor neurons of locusts. In all other neurons that receive their input through synapses on their dendrite, the input current reaching the site of action-potential generation can also be adapted. The various forms of synaptic dynamics, in particular short-term depression, are potential sources of spike-frequency adaptation. lso various ionic currents on the dendrite might potentially adapt the input.

5 J. Benda & J. Tabak: Spike-frequency adaptation 5 Network Effects In case the spiking activity of a neuron leads to the activation of recurrent inhibitory synapses, this might cause spike-frequency adaptation in a way quite similar as the adaptation currents (Sutherland et al., 29). Here, the slow time scale would be introduced by the time scale of the postsynaptic current of the inhibitory synapse, by the delay introduced by wiring, and by the integration dynamics of the neurons participating in this negative feedback loop. Similar to input-driven adaptation, some forms of feedforward input may also cause spikefrequency adaptation in the target neuron. In its simplest form, the presynaptic neuron already adapts, and thus the postsynaptic neuron very likely will adapt as well. In addition, the synaptic transmission might adapt through synaptic plasticity or by presynaptic inhibition, or delayed and slower inhibitory input can adapt the response of the target neuron. Network and cell-intrinsic adaptation mechanisms shape in specific ways functionally different response properties of auditory interneurons in locusts (Hildebrandt et al., 29). Models of Spike-Frequency daptation There are many ways to model spike-frequency adaptation. In the following overview, we focus on output-driven adaptation, i.e., adaptation currents that are activated by the neuron s output activity (membrane voltage with action potentials). In contrast, input-driven adaptation is activated by the input signal irrespective of the resulting output. Conductance-Based Models In conductance-based models like the Hodgkin-Huxley model, any ionic current of a neuron is modeled in detail. This requires exact knowledge of the specific type of adaptation current involved. For example, the M-type current I M would be modeled as I M = ḡ M a(v E K ) (1) τ a (V ) da dt = a (V ) a (2) where E K is the potassium reversal potential, ḡ M is the maximum conductance, and a is the gating variable of this voltage-gated current. τ a (V ) and a (V ) are the voltage-dependent time constant and activation variable, respectively, that determine the (slow) temporal evolution of the gating variable. a (V ) is a sigmoidal function depending on the membrane potential V that is close to zero at the resting potential and about one during an action potential (Brown and dams, 198; Benda and Herz, 23). Integrate-and-Fire Models s for conductance-based models, any adaptation current can also be simply added to the input of an integrate-and-fire model. The only difference is that the effect the action potential has on the dynamics of the current has to be modeled explicitly, for example, by increasing a in equation (2) by a certain increment whenever there was a spike. Often, however, a more generalized form of adaptation currents is used in (leaky, quadratic, exponential) integrate-and-fire neurons (Izhikevich, 23; Brette and Gerstner, 25). This has

6 J. Benda & J. Tabak: Spike-frequency adaptation 6 the advantage that the exact nature of the adaptation current in question does not have to be known. First, the electromotoric force E K V is often approximated away, so that the gating variable a then becomes the (appropriately scaled) adaptation current. Second, a (V ) is replaced by incrementing by some appropriate value whenever an action potential occurs at time t i. Between the action potentials decays exponentially with a fixed adaptation timeconstant τ a : τ a d dt I M (3) = + i= δ(t t i ). (4) The adaptation dynamics (4) can be extended by an additional term proportional to V in order to incorporate the fact that some adaptation currents are already activated by subthreshold membrane voltages. Integrate-and-fire models with a dynamic threshold are another class of adapting neuron models, where the firing threshold is incremented by each action potential and follows a slow dynamics between the action potentials. lthough this mechanism also produces spike-frequency adaptation (Liu and Wang, 21), the resulting adapted f -I curves do not only shift to higher inputs, as any adaptation current does, but also reduce their slopes (Benda et al., 21). Firing-Rate Models By simply approximating the neuron s spike generator by its f -I curve f (I) and using Kirchhoff s law that ionic currents over the cell s membrane add up yields the following equation for the spike frequency f (t): f = f (I ) (5) d τ a dt = + ( f ) (6) (Benda and Herz, 23). Because the adaptation current acts inhibitory on the input current, it is subtracted from the input. The dynamics for is the same as for the simplified integrateand-fire neurons, but it is driven via the term ( f ) by the output spike frequency and not by individual spikes any more. This dependence on the spike frequency arises from averaging over the spike train in time (Wang, 1998). This averaging procedure is possible if the spike frequency is higher than the reciprocal adaptation time-constant. Only then the adaptation dynamics can be separated from the one of the membrane potential. Choosing ( f ) to be proportional to f often is a reasonable approximation. Then, the only nonlinearity remaining is the neuron s onset f -I curve. By linearizing the onset f -I curve around the current spike frequency, one can easily compute the neuron s filter properties that are introduced by the adaptation process. Signal Processing The slow adaptation variable (gating variable of M-type current, intracellular calcium concentration for HP-type currents, etc.) is a low-pass filtered version of the output spike train (output-driven adaptation, Eqs. (4) or (6)) or of the input to the neuron (in the case of inputdriven adaptation). Since this low-pass filtered signal is then subtracted from the input current

7 J. Benda & J. Tabak: Spike-frequency adaptation f (1n,1n) Gain [Hz/n] f (1n) Stimulus frequency [Hz] Figure 4: High-pass filter of adapting neurons. daptation processes add high-pass filter characteristics to a neuron s transfer function between the input to the neuron (here, e.g., an injected current) and the resulting spike-frequency. Shown is the gain, i.e., the amplitude of the modulation of the spike-frequency response divided by the amplitude of the injected current input as a function of the frequency of that input, obtained for the same model neuron as in Fig. 1. s a stimulus a low-pass filtered Gaussian white noise with cutoff frequency at 16 Hz, mean 1 n, and standard deviation 2 n was used. The slopes of the steady-state f -I curve (red) at 1 n and the one of the adapted f -I curve (green) for I = 1 n are plotted for comparison. I, this basically (i.e., with linear f -I curves) results in a high-pass filter between the input to the neuron and the resulting output spike frequency. The spike frequency of an adapting neuron thus encodes the high-pass filtered input signal (Fig. 4; Benda and Herz, 23). For example, in electroreceptors, this high-pass filter resulting from adaptation processes enhances the response to fast stimulus components in communication signals (Benda et al., 25). The gain of the filter at low stimulus frequencies that are below the reciprocal adaptation time-constant is given by the slope of the steady-state f -I curve. The gain of stimulus frequencies faster than the reciprocal adaptation time-constant is given by the slope of the adapted f -I curves. In the case of purely subtractive adaptation the latter equals the slope of the onset f -I curve. It follows that neurons with flat steady-state f -I curves (zero slope) implement perfect intensity invariance (zero gain at low stimulus frequencies), i.e., no matter on what mean intensity a certain input waveform is delivered, the output spike-frequency will always be the same, since the mean intensity is completely filtered out by the adaptation process (Benda and Hennig, 28). The high-pass filter effect of spike-frequency adaptation considers the linear aspects of spike-frequency adaptation only. t least the steady-state f -I curve is indeed linearized by adaptation currents (as shown in Ermentrout, 1998, and also illustrated in Fig. 1 B). In addition, the non-linear shape of a neuron s f -I curves (e.g., rectification at zero firing rate and saturation) together with the high-pass filter introduces additional aspects. For example, the subtractive shift of the adapted f -I curves to higher input intensities might be accompanied by a reduction in their maximum firing rate and/or by a reduction of their slope (see Fig. 5 for simple examples) and thereby influence the filter properties of a neuron in a non-linear way.

8 J. Benda & J. Tabak: Spike-frequency adaptation 8 The various forms of this interplay of adaptation dynamics and the resulting high-pass characteristics with nonlinearities lead to diverse functional consequences. For example, in the auditory system of crickets, the slow calcium dynamics leads to forward masking of weaker stimuli (Sobel and Tank, 1994). collision-detecting neuron in the visual system of locusts responds very selectively to looming stimuli because of spike-frequency adaptation (Peron and Gabbiani, 29). Spike frequency f [Hz] f (I) f (I) f (I,5) f (I,9) f (I,16) B Spike frequency f [Hz] f (I) f (I) f (I, 4.6) f (I, 3.2) f (I, 1.1) Input I 1 1 Input I Figure 5: Nonlinearities and adaptation. Simple subtractive adaptation can have unexpected effects when it occurs after a nonlinearity. The examples show two simple scenarios where the firing rate f of a neuron is directly proportional to the input current J, and an adaptation current acts subtractive on the current, as in Eq. (5): f = c (J ), with c >, and f = if J <. The input current, however, depends non-linearly on some (sensory) input I: J(I). The input nonlinearity is a sigmoidal function of the sensory input: J(I) = (1 + exp( (J 6)/4)) 1.2 (c = 32 Hz). This is the case, for example, in auditory or olfactory receptor neurons. lthough the adaptation acts subtractive on the input current, the resulting adapted f -I curves appear to be shifted to the right as well as downwards. B J(I) = cos(.2j)+.2 (c = 18 Hz). This models, for example, an orientation-selective cell in the visual system. Here, the adapted f -I curves appear to be shifted downward. daptation currents can even alter the type of action-potential dynamics. type I neuron where repetitive firing occurs via a saddle-node bifurcation on an invariant cycle can switch into a type II neuron where the resting potential loses stability via a Hopf bifurcation by the activation of an M-type adaptation current (Ermentrout et al., 21; Prescott et al., 26). Multiple daptation Time-Scales Usually, spike-frequency adaptation occurs not only on a single but on several time scales. The shortest time scales of spike-frequency adaptation are in the range of 1 to 1 ms. This at the same time usually also is the strongest adaptation process in the neuron. dditional slower adaptation processes are weaker and have time constants of seconds, 1 seconds and slower. daptation on several time scales could resemble power-law adaptation (Clarke et al., 213) or logarithmic adaptation (Xu et al., 1996) in contrast to the single adaptation process showing exponential adaptation of the firing rate.

9 J. Benda & J. Tabak: Spike-frequency adaptation 9 Noise daptation currents are a potential source of ion channel noise. Since these ionic currents are characterized by their slow adaptation dynamics, the resulting channel noise is a colored noise with correlation time given by the adaptation time-constant. Such a noise results in positive correlations between successive interspike intervals and in interspike-interval histograms with a sharper peak and heavier tail in comparison to the ones obtained with white noise (Fig. 6 D F; Schwalger et al., 21). Output-driven adaptation also interacts with white noise, e.g., generated by ion channel noise from fast spike-generating currents. In this case, adaptation introduces negative serial correlations between successive interspike intervals (Fig. 6 C; Chacron et al., 2; Benda et al., 21) that improve information transmission at low frequencies (Chacron et al., 21). The interspike-interval statistics of auditory receptor neurons in locusts have been shown to be shaped by both white noise sources and colored noise presumably generated by an adaptation current (Fisch et al., 212). p [Hz] D p [Hz] ISI [ms] ISI [ms] B ISIi+1 [ms] E ISIi+1 [ms] ISI i [ms] ISI i [ms] C ρ F ρ lag lag Figure 6: Spike-frequency adaptation and noise. C Interspike -interval (ISI) histogram (), return map B, and ISI correlations C for the conductance-based model with an M-type adaptation current (same model as in Fig. 1) driven with white noise (mean I = 1 n, noise strength D = 1 n 2 /Hz). D F The same measures for colored noise arising from the stochasticity of the adaptation current. Noise strength was scaled to result in the same standard deviation of the ISIs (D = 3 n 2 /Hz). Note the more peaked shape of the ISI distribution as well as the positive ISI correlations. Rhythmogenesis in Neural Networks The slow decrease in firing frequency due to adaptation currents provides delayed negative feedback that can result in the production of oscillations by excitatory neural networks. The population activity of an excitatory neural network can be roughly described by its mean firing

10 J. Benda & J. Tabak: Spike-frequency adaptation 1 rate, which varies according to τ d f dt = f (I) f (7) (Wilson and Cowan, 1972). Here, the firing frequency f is averaged spatially (over the population) and temporally (over the synaptic integration time scale). It is assumed that the neurons fire asynchronously. f is the steady-state firing rate that will be reached with a network integration time-constant τ. It is similar to the function f in Eq. (5), but is not strictly the firing rate of one neuron, as it also incorporates information about synaptic dynamics, heterogeneity in the neuron population, etc. Because of the recurrent excitatory connectivity, the input I to the network is given by I = w f, where f is the output firing rate of the network, and w is a connectivity parameter akin to synaptic weight or the number of connections per neuron. If neurons are adapting, the effective input becomes I = w f, where is the average adaptation current in the population. This description is valid as long as adaptation is much slower than network integration (τ a > τ). 1 a B 1 a f, b c d 1 2 Time [arbitrary unit] f d c 1 b Figure 7: Oscillatory activity in an excitatory network with adaptation. Time courses of f, the average firing rate (black), and, the average adaptation current (blue) in the population. During an episode of activity ( f 1), is increasing, until activity drops to. During the inter-episode interval, decreases until the network is excitable enough to start a new high-activity episode. B The two time courses from describe a closed trajectory in the (, f ) plane. The red curve is the curve of steadystates of f, solutions of Eq. (8) at each value of. There is an interval of over which there are two stable steady states for f, one high and one low (and an intermediate, unstable steady state, dashed). The blue curve represents ( f ); when the phase point (, f ) is above this curve increases, below this curve decreases, according to Eq. (6). During an episode of network activity (a), f is high so increases, moving the trajectory to the right. fter the point where the high steady state meets the unstable steady state, the only available steady state is at low f. The trajectory jumps down (b), starting the inter-episode interval (c), during which decreases. When is low enough, the lower steady state disappears and the trajectory jumps up (d), beginning a new high activity episode. The population firing rate equilibrates to f = f (w f ) (8) and varies according to Eq. (6). Note that Eq. (8) can have two stable solutions over a range of values. This bistability, due to the fast positive feedback provided by excitatory connections, together with the slow dynamics, produce relaxation oscillations as illustrated in Fig. 7 (Tabak et al., 2; van Vreeswijk and Hansel, 21; Giugliano et al., 24; Gigante et al., 27).

11 J. Benda & J. Tabak: Spike-frequency adaptation 11 Note that the exact dependence of adaptation on frequency (here given by Eq. (6) and represented as the blue curve in Fig. 7 B) may affect the type of bifurcation that occurs between a silent and an oscillatory network as a source of excitation is increased (Nesse et al., 28), but it does not change the fundamental mechanism of the oscillations. So the excitatory network can oscillate whether adaptation is input- or output-driven (Nesse et al., 28). Synaptic adaptation (often called short-term synaptic depression) also produces oscillations in excitatory networks, through the same mechanism as the intrinsic adaptation mechanisms described here (Tabak et al., 2; Tsodyks et al., 2; Wiedman and Luthi, 23). Note, however, that the two types of adaptation (intrinsic or synaptic) can be distinguished in models (even in firing rate models) as well as in experiments based on the models predictions (Tabak et al., 2, 21; Wiedman and Luthi, 23). Thus, it is important to pay careful attention to what adaptation model is most appropriate in a particular modeling context. In addition to its rhythmogenic properties in excitatory networks, activity-dependent adaptation facilitates network-based oscillations due to the interaction of excitatory and inhibitory populations (ugustin et al., 213). The rhythmogenic mechanism based on excitatory connectivity and slow activity-dependent adaptation may underlie various types of network activity such as episodic activity in developing networks (Butts et al., 1999; Tabak et al., 2), respiratory rhythms (Kosmidis et al., 24), or slow oscillations between up and down states in cortex (Compte et al., 23). References ugustin M, Ladenbauer J, Obermayer K (213) How adaptation shapes spike rate oscillations in recurrent neuronal networks. Front. Comput. Neurosci. 7: 9. Benda J, Hennig RM (28) Dynamics of intensity invariance in a primary auditory interneuron. J. Comput. Neurosci. 24: Benda J, Herz VM (23) universal model for spike-frequency adaptation. Neural Comput. 15: Benda J, Longtin, Maler L (25) Spike-frequency adaptation separates transient communication signals from background oscillations. J. Neurosci. 25: Benda J, Maler L, Longtin (21) Linear versus nonlinear signal transmission in neuron models with adaptation-currents or dynamic thresholds. J. Neurophysiol. 14: Brette R, Gerstner W (25) daptive exponential integrate-and-fire model as an effective description of neuronal activity. J. Neurophysiol. 94: Brown D, dams PR (198) Muscarinic suppression of a novel voltage-sensitive K + current in a vertebrate neuron. Nature 183: Butts D, Feller M, Shatz C, Rokhsar D (1999) Retinal waves are governed by collective network properties. J. Neurosci. 19: Chacron MJ, Longtin, Maler L (21) Negative interspike interval correlations increase the neuronal capacity for encoding time-dependent stimuli. J. Neurosci. 21: Chacron MJ, Longtin, St-Hilaire M, Maler L (2) Suprathreshold stochastic firing dynamics with memory in P-type electroreceptors. Phys. Rev. Lett. 85: Clarke SE, Naud R, Longtin, Maler L (213) Speed-invariant encoding of looming object

12 J. Benda & J. Tabak: Spike-frequency adaptation 12 distance requires power law spike rate adaptation. Proc. Natl. cad. Sci. US 11: Compte, Sanchez-Vives M, McCormick D, Wang XJ (23) Cellular and network mechanisms of slow oscillatory activity (<1 Hz) and wave propagations in a cortical network model. J. Neurophysiol. 89: Ermentrout B (1998) Linearization of f -I curves by adaptation. Neural Comput. 1: Ermentrout B, Pascal M, Gutkin B (21) The effects of spike frequency adaptation and negative feedback on the synchronization of neural oscillators. Neural Comput. 13: Fisch K, Schwalger T, Lindner B, Herz VM, Benda J (212) Channel noise from both slow adaptation currents and fast currents is required to explain spike-response variability in a sensory neuron. J. Neurosci. 32: Gigante G, Mattia M, Giudice PD (27) Diverse population-bursting modes of adapting spiking neurons. PhysRevLett 98: Giugliano M, Darbon P, rsiero M, Lüscher HR, Streit J (24) Single-neuron discharge properties and network activity in dissociated cultures of neocortex. J. Neurophysiol. 92: Gollisch T, Herz VM (24) Input-driven components of spike-frequency adaptation can be unmasked in vivo. J. Neurosci. 24: Hildebrandt KJ, Benda J, Hennig RM (29) The origin of adaptation in the auditory pathway of locusts is specific to cell type and function. J. Neurosci. 29: Hildebrandt KJ, Benda J, Hennig RM (211) Multiple arithmetic operations in a single neuron: the recruitment of adaptation processes in the cricket auditory pathway depends on sensory context. J. Neurosci. 31: Izhikevich EM (23) Simple model of spiking neurons. IEEE TNN 14: Kosmidis E, Pierrefiche O, Vibert JF (24) Respiratory-like rhythmic activity can be produced by an excitatory network of non-pacemaker neuron models. J. Neurophysiol. 92: Liu YH, Wang XJ (21) Spike-frequency adaptation of a generalized leaky integrate-and-fire model neuron. J. Comput. Neurosci. 1: Nesse W, Borisyuk, Bressloff P (28) Fluctuation-driven rhythmogenesis in an excitatory neuronal network with slow adaptation. J. Comput. Neurosci. 25: Peron S, Gabbiani F (29) Spike frequency adaptation mediates looming stimulus selectivity in a collision-detecting neuron. Nat. Neurosci. 12: Prescott S, Ratté S, Sejnowski TJ (26) Nonlinear interaction between shunting and adaptation controls a switch between integration and coincidence detection in pyramidal neurons. J. Neurosci. 26: Sah P (1996) Ca 2+ -activated K + currents in neurones: types, physiological roles and modulation. Trends Neurosci. 19: Schwalger T, Fisch K, Benda J, Lindner B (21) How noisy adaptation of neurons shapes interspike interval histograms and correlations. PLoS Comput. Biol. 6: e1126. Sobel EC, Tank DW (1994) In vivo Ca 2+ dynamics in a cricket auditory neuron: an example of chemical computation. Science 263: Sutherland C, Doiron B, Longtin (29) Feedback-induced gain control in stochastic spiking networks. Biol. Cybern. 1: Tabak J, Mascagni M, Bertram R (21) Mechanism for the universal pattern of activity in

13 J. Benda & J. Tabak: Spike-frequency adaptation 13 developing neuronal networks. J. Neurophysiol. 13: Tabak J, Senn W, O Donovan M, Rinzel J (2) Modeling of spontaneous activity in the developing spinal cord using activity-dependent depression in an excitatory network. J. Neurosci. 2: Tsodyks M, Uziel, Markram H (2) Synchrony generation in recurrent networks with frequency-dependent synapses. J. Neurosci. 2: RC5. van Vreeswijk C, Hansel D (21) Patterns of synchrony in neural networks with spike adaptation. Neural Comput. 13: Wang XJ (1998) Calcium coding and adaptive temporal computation in cortical pyramidal neurons. J. Neurophysiol. 79: Wiedman U, Luthi (23) Timing of network synchronization by refractory mechanisms. J. Neurophysiol. 9: Wilson H, Cowan J (1972) Excitatory and inhibitory interactions in localized populations of model neurons. Biophys. J. 12: Xu Z, Payne JR, Nelson ME (1996) Logarithmic time course of sensory adaptation in electrosensory afferent nerve fibers in a weakly electric fish. J. Neurophysiol. 76: