| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1Department of Neuroscience, Baylor College of Medicine; 2Computational and Applied Mathematics, Rice University, Houston, Texas; 3Department of Zoology, Cambridge University, Cambridge; and 4Department of Bioengineering, Imperial College London, London, United Kingdom
Submitted 23 January 2006; accepted in final form 5 March 2006
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSACKNOWLEDGMENTSREFERENCES |
|---|
, a membrane time constant of about 8 ms, and exhibited inward rectification and rebound spiking after hyperpolarizing current pulses. Responses to depolarizing current pulses revealed the neuron's intrinsic bursting properties and pronounced spike-frequency adaptation. The characteristics of adaptation, including its time course, the attenuation of the firing rate, the mutual dependency of these two variables, and their dependency on injected current, closely followed the predictions of a model first proposed to describe the adaptation of cat visual cortex pyramidal neurons in vivo. Our results thus validate the model in an entirely different context and suggest that it might be applicable to a wide variety of neurons across species. Spike-frequency adaptation is likely to play an important role in tuning the LGMD and in shaping the variability of its responses to visual looming stimuli. |
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSACKNOWLEDGMENTSREFERENCES |
|---|
; Sah 1996). Spike-frequency adaptation has a variety of important consequences for the processing of sensory stimuli. For example, it allows to emphasize changes in stimulus parameters, but can also protect neurons from firing frequency saturation and tune sensory responses to specific stimulus features, as demonstrated by several in vivo and in vitro studies (Benda et al. 2005; Sobel and Tank 1994; Wang et al. 2003). Recently, a biophysical model of spike-frequency adaptation was proposed (Liu and Wang 2001; Wang 1998) based on intracellular in vivo recordings of cat visual cortical neurons (Ahmed et al. 1998). The model made several specific predictions on how the firing frequency of adapting neurons should evolve over time and on the relation between variables characterizing spike-frequency adaptation that have not yet been tested experimentally. Although this model was based on cortical pyramidal cell data, it is likely to be widely applicable. This is explained by the fact that it is based either on a simple compartmental model endowed with conductances of the Hodgkin–Huxley type (Wang 1998), which replicates the firing properties of a wide class of neurons (Mainen and Sejnowski 1996; Pinsky and Rinzel 1994), or on a leaky integrate-and-fire approximation (Liu and Wang 2001; for related, phenomenological models, see Benda and Herz 2003; La Camera et al. 2004). We set out to experimentally test the model's predictions on an identified interneuron in the locust visual system, the lobula giant movement detector (LGMD).
The LGMD is thought to be involved in collision avoidance and escape behavior (for review, see Gabbiani et al. 2004). It projects onto a postsynaptic target neuron called the descending contralateral movement detector (DCMD) that can easily be recorded from the animal's nerve cord (O'Shea and Williams 1974). The connection between the LGMD and DCMD is mediated by a mixed chemical/electrical synapse (Killmann et al. 1999) and is so strong that spikes in the LGMD and DCMD are in a 1:1 correspondence under visual stimulation (Gabbiani et al. 2005; O'Shea and Williams 1974; Rind 1984). The LGMD can thus be conveniently identified electrophysiologically as the unique neuron in the optic lobe whose spikes precede by nearly 2 ms those of the DCMD. An extensive literature based on DCMD recordings shows that both the LGMD and the DCMD are most sensitive to objects approaching on a collision course (e.g., Hatsopoulos et al. 1995; Rind and Simmons 1992; Schlotterer 1977). Several aspects of the computation the LGMD performs while it is tracking approaching objects were the subject of previous investigations (Gabbiani et al. 1999, 2001; Hatsopoulos et al. 1995). This computation is thought to rely on a multiplicative interaction between two of the inputs impinging on the LGMD's dendritic arborizations (Gabbiani et al. 2002).
The spatial receptive field structure of the LGMD was recently characterized as were the summation properties of local excitatory inputs and the role of feedforward inhibition on the neuron's responses to looming stimuli (Gabbiani et al. 2005; Krapp and Gabbiani 2005). The intrinsic membrane properties of the LGMD have not yet been studied in detail, however. We thus started by characterizing them before comparing the LGMD adaptation characteristics to those of the spike-frequency adaptation models proposed by Wang (1998) and Liu and Wang (2001). Recording from a uniquely identified neuron, such as the LGMD, as opposed to different neurons from a population reduces the number of potential sources of variability. The reproducibility of our LGMD data proved advantageous to test the predictions made by these biophysical adaptation models. Furthermore, one of the models' predictions would have been very difficult to test on a heterogeneous neuronal population.
|
METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSACKNOWLEDGMENTSREFERENCES |
|---|
Dissections and animal preparation were similar to those described in detail in previous publications (Gabbiani et al. 1999, 2001, 2002, 2005). Briefly, experiments were carried out on adult locusts (Schistocerca americana) taken from the laboratory colony 3–4 wk after their final molt. Animals were fixed in a plastic holder and their head was bathed in locust saline. The head capsule was opened and the gut was removed. The connective contralateral to the recorded side was either cut at the level of the subesopageal ganglion and placed in a suction electrode or recorded extracellularly with hook electrodes after cutting a window in the pronotum of the animal. The brain was exposed and desheathed with fine forceps. To minimize brain movement during intracellular recordings, the mandibular muscles were sectioned and a support was placed under the brain.
Electrophysiology and data acquisition
Extracellular recordings were amplified differentially 10,000 times with standard equipment. The DCMD was easily identified as the unit with the largest extracellular action potentials. Intracellular recordings were obtained with sharp electrodes pulled with a horizontal puller (P-87, Sutter Instruments, Novato, CA) using thin-wall borosilicate glass (1.2/0.9-mm outer/inner diameters; WPI, Sarasota, FL). The pull parameters were optimized to minimize electrode resistance and thus facilitate the passage of large currents necessary for this study. The electrodes were filled with potassium acetate (2 M) and their DC resistance, measured at the beginning of each experiment, ranged between 30 and 50 M. The LGMD was identified as the unique neuron whose spikes were in one-to-one correspondence with those of the DCMD (e.g., Gabbiani et al. 2002; Fig. 4A; Gabbiani et al. 2005; Fig. 2A). As explained in RESULTS, we were interested in correlating single spike shape characteristics with the properties of the current versus spike-frequency discharge curve in individual LGMD neurons. We therefore attempted whenever possible to obtain multiple penetrations in each neuron. After a successful penetration, the electrode was carefully retracted and subsequent attempts at impaling the LGMD were directed to dendritic locations in the optic lobe either more proximal or more distal from the spike initiation zone based on anatomical markers. Electrode withdrawal and renewed penetration did not alter the response properties of the LGMD under visual stimulation or affect its input resistance and the stability of subsequent recordings, suggesting that no measurable damage resulted from this procedure. To minimize the effects of electrode impedance, all current-injection protocols were carried out using the discontinuous current-clamp mode of the intracellular amplifier (Axoclamp 2B). Discontinuous current-clamp sampling rates and electrode capacity compensation were adjusted on a cell-by-cell basis to allow settling of the potential to its steady-state value before sampling. Sampling rates ranged between 4 and 8 kHz and were well above the cutoff frequency set by the LGMD membrane time constant (see RESULTS). For a few recordings, a SEL10 amplifier (npi electronic, Tamm, Germany) was used at 15- to 20-kHz sampling rates. Current pulses were programmed using a Master 8 pulse stimulator (AMPI, Jerusalem, Israel) used to drive the external current input of the amplifier. The intra- and extracellular recordings, current traces, and a voice channel were stored on digital audio tape using a professional recorder (DT800, sampling rate: 8 kHz; MicroData Instrument, South Plainfield, NJ). Subsequently, the data were transferred to a personal computer for analysis. Intracellular recordings typically lasted 40–60 min.
|
|
|
|
Each current-injection trial consisted of a square pulse lasting 500 ms. For all but one experiment, positive current pulses were preceded by a 500-ms, –2-nA negative prepulse to minimize possible inactivation of sodium channels. Each trial was repeated eight to ten times before selecting a new current value and two trials were separated by 15 s.
The first protocol consisted of stimulation with negative current pulses of –1, –2, and –5 nA, with eight to ten trials per value. We recorded responses to this protocol in 12 different neurons. In five experiments, we were able to obtain a second or a third penetration in the same neuron. In eight of these penetrations, the protocol could be repeated a second or third time, yielding a total of 23 measurements.
The second protocol consisted of stimulation with negative current pulses of –1 to –12, –15, and –20 nA. We recorded responses to this protocol in two neurons (three penetrations, seven measurements).
The third protocol consisted of stimulation with positive current pulses of +1 to 10, 12, 15, 20, and 25 nA. We recorded responses to this protocol in 13 different neurons, but not all neurons could be recorded at all current values. Three neurons were recorded 
25, five 20, one 15, one 12, two 10, and one 7 nA (Fig. 6A). We obtained multiple penetrations in three neurons for a total of 17 penetrations and 32 measurements. When the currents had been injected in increasing order during a measurement, the subsequent measurement was performed in descending order (i.e., from 25 to 1 nA) to allow detection of the effects of current-injection sequence on the cell's response. No sequence effects (ascending vs. descending) could be observed.
|
|
Multiple measurements during a single penetration were used to monitor the stability and reproducibility of the recordings. Data from different penetrations were treated as independent in pooled analyses. Restricting pooled analyses to different cells did not change any of the results.
RESTING MEMBRANE POTENTIAL. In 11 recordings, we determined the resting membrane potential at the end of the experiment by comparing the intracellular amplifier potential reading within the neuron with that obtained immediately after retracting the electrode to the bath.
INPUT RESISTANCE. An estimate of input resistance was computed for current pulses of –1, –2, and –5 nA by selecting a 300-ms window for each trial, averaging the membrane potential within the window and dividing by the injected current value (Fig. 1A). Usually, the averaging window was centered at the midpoint of the pulse. In some cases, it was shifted earlier or later to avoid clearly visible excitatory postsynaptic potentials. The input resistance was then averaged across the eight to ten trials obtained for each current measurement. For depolarizing pulses of +1, 2, and 3 nA, an estimate of the input resistance was computed similarly, except that each trace was first median-filtered over a 6-ms time window to suppress eventual action potentials (Fig. 1B). The averaging window typically constituted the last 300 ms of the pulse because action potentials were less likely to occur in that period (see RESULTS). In a few trials, one or at most two action potentials were included in the averaging window, but their impact was minimal after median filtering.
MEMBRANE AND EQUALIZATION TIME CONSTANTS.
The membrane and equalization time constants were computed by using the peeling method for current pulses of –1, –2, and –5 nA (Rall 1969). The mean membrane potential averaged across eight to ten trials in response to the pulse was plotted in logarithmic coordinates after subtracting its minimum steady-state value. The data were fitted to a straight line by least squares (Fig. 2B). The absolute value of the inverse line slope yielded the membrane time constant. The fitted line was then subtracted from the data, revealing a second, faster linear decay in logarithmic coordinates that was also fitted to a straight line to obtain the equalization time constant. We verified the accuracy of the double-exponential fit by comparing it to the original mean membrane potential time course and its SD obtained from repeated trials (Fig. 2A). We also verified that a direct, double-exponential fit of the membrane potential time course using least squares yielded identical results (Holmes et al. 1992).
INTERSPIKE INTERVAL DISTRIBUTIONS.
Interspike intervals (ISIs) were computed by subtracting successive spike occurrence times during each positive current pulse. ISIs were pooled across trials and current values. Histograms were obtained by binning the interval distribution in 30 bins of equal size between 0 and 100 ms (Fig. 4B). The coefficient of variation (CV) of the ISI distribution was obtained by dividing the SD of the ISI distribution (
t) by its mean tm (Fig. 8A). The ISI serial correlation coefficient was computed according to the formula
 |
|
|
ADAPTATION PARAMETERS 
ADAP AND FADAP.
The time constant of spike frequency adaptation adap was obtained by fitting an exponentially decaying function with three free parameters (initial firing rate f0 fit, exponential adaptation time constant adap, and steady-state firing rate fss fit) to fm(t) by least squares (Fig. 5, A and D). We verified that double exponential fits did not improve significantly the fit quality, taking into account the SD of fm(t). Following Ahmed et al. (1998) and Wang (1998), the attenuation factor, Fadap, was defined as Fadap = (f0 – fss)/f0 (Fig. 5C).
MEAN SPIKE HEIGHT AND WIDTH. These parameters were measured in 11 neurons (14 different penetrations) using 10 isolated spikes from 10 different trials obtained in response to the lowest positive current above spiking threshold. Spike height was measured as the difference between the value of the membrane potential at the peak of the spike and the value at its inflection point during the initial depolarization leading to the spike. In practice, the inflection point was determined by finding where the first derivative of the membrane potential exceeded 5 mV/ms. Spike width was then measured at half-height. The values obtained for each trial were averaged across trials to obtain the mean spike heights, widths, and their SDs (Fig. 6F, Table 1).
|
FIT OF F0(I) CURVE TO THAT OF A LEAKY INTEGRATE-AND-FIRE (LIF) MODEL.
For positive current pulses, we fitted by least squares the firing rate derived from the first ISI as a function of current, f0(I), to that of an LIF neuron with variable reset. If we denote by vrest the resting membrane potential (relative to an extracellular reference); vth, the spiking threshold; v0, the reset potential after a spike; , the membrane time constant; rin, the input resistance; and tref, the refractory period, then the LIF firing rate, as a function of injected current I, is given by
 | (1) |
th = vth – vrest and 0 = v0 – vrest. In those fits, both vth and v0 were allowed to vary and the other parameters were fixed. The resting membrane potential vrest was set to –70 mV, tref was equal to 1.5 ms, and and rin were determined from the responses to the hyperpolarizing current protocols.
LEAKY INTEGRATE-AND-FIRE (LIF) MODEL OF ADAPTATION.
We compared the adaptation time course of the LGMD firing rate to that of an LIF model similar to that studied by Liu and Wang (2001). The model includes two dynamic variables: 1) the membrane potential v(t) and 2) the intracellular calcium concentration x(t). The subthreshold dynamics of the model is described by the following two differential equations
 | (2) |
 | (3) |
Ca (e.g., Helmchen et al. 1996; Traub 1991). An action potential is generated when v(t) reaches vth and the membrane potential is subsequently reset to v0. The membrane potential stays at this value during the refractory period tref before resuming its dynamical evolution according to Eq. 2. After each action potential, the calcium concentration is incremented by 
(i.e., x 
x + ) and immediately resumes its decay following Eq. 3. The linear dependency of the potassium-dependent current on calcium concentration can be justified by linearizing the model of Wang (1998)
 |
= 0.0267 µM–1. The approximation is very accurate for x < 10 µM, as was the case in our simulations. Note that 
ahp = 
ahp has units of conductance/concentration (in µS/µM). We set vrest = –70 mV, vth = –58 mV, v0 = –62 mV, = 8 ms, rin = 5 M, and tref = 1.5 ms. These parameters were determined from the analyses described above (see RESULTS). The calcium extrusion time constant was set to Ca = 130 ms, a value that will be justified in the RESULTS. We followed Liu and Wang (2001) and set = 0.2 µM and vK = –80 mV. The value of ahp was adjusted to obtain attenuation factors similar to those observed experimentally and was set to 0.12 µS/µM. Note that it is about 10-fold higher than the value used by Liu and Wang (2001; 0.015 µS/µM). This is expected because the input resistance of the LGMD is about a factor 10 lower than that of pyramidal neurons and the degree of adaptation in the LGMD (see RESULTS) is similar to that of pyramidal cells (Ahmed et al. 1998). The model was simulated using a fourth-order Runge–Kutta algorithm (Henrici 1982) with a time step of 0.01 ms.
AHP DECAY IN THE LIF ADAPTATION MODEL.
After a current pulse, Eqs. 2 and 3 imply that both the calcium concentration and the calcium-dependent potassium current decay exponentially toward zero with a time constant Ca. Consequently, the membrane potential v(t) relaxes toward its resting value vrest, with a time course that can be obtained from Eq. 2 by setting dv/dt = 0 (using the approximation << tCa). For typical average values of the peak AHP (–4 mV; see RESULTS), the relaxation is slightly slower than an exponential decay with time constant Ca.
Data analysis was performed with Matlab (The MathWorks, Natick, MA). Data fitting was carried out using the optimization toolbox least-square fitting routines. The LIF integration routine used to implement Eqs. 2 and 3 was coded in C and accessed from Matlab as a Mex file to speed up simulations. In the following, correlation coefficients are indicated by 
. The Smith–Satterthwaite procedure used in Table 1 for comparing means of distributions with unequal variances is described in Milton and Arnold (1995).
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSACKNOWLEDGMENTSREFERENCES |
|---|
and Liu and Wang (2001). Electrotonic properties of the LGMD
The resting membrane potential of the LGMD had a typical mean value of –68 mV (SD: 6.6 mV, n = 11 neurons). Figure 1, A and B illustrates responses to three hyperpolarizing and depolarizing current pulses of varying amplitude, respectively. Figure 1C illustrates the input resistance distributions derived from such experiments by averaging over 300-ms windows (METHODS). Bars of different shades represent input resistance measurements for the six values of injected current. Figure 1D gives the mean and SD of these distributions. There was a slight but nonsignificant trend for the input resistance to decrease with increasing current (P > 0.8, one-way ANOVA; max. difference between –5 and +3 nA: 0.8 M). The overall mean input resistance pooled across all current values was 4.3 M (SD: 1.6; median: 4.0).
Figure 2A shows the response of a single neuron within the first 60 ms of a –5-nA current pulse. Relaxation to the membrane potential steady-state value at the end of the pulse occurred in two phases. The first, rapid phase typically lasted roughly 1.5 ms and was followed by a slower relaxation over the next nearly 40 ms. Accordingly, the membrane potential time course could be fit by a double exponential (Fig. 2A, gray line) whose parameters were determined using the exponential peeling method (Fig. 2B). The first time constant is thought to reflect rapid charge redistribution across the spatial extent of the neuron and will be referred to as the "equalization" time constant, e (Holmes et al. 1992). The second slower time constant reflects local charging of the membrane and thus corresponds to the membrane time constant . In the example depicted in Fig. 2, A and B, = 7.8 ms and e = 0.3 ms. Figure 2, C and D shows histograms of the measured values of and e averaged across three current values (–1, –2, and –5 nA, respectively). The means of and e pooled across currents were 7.3 and 0.34 ms, respectively (SDs: 2.1 and 0.24; medians: 7.8 and 0.26). There was a trend for to increase with the amplitude of the pulse (one-way ANOVA, P < 0.0001). The mean values of at –1, –2, and –5 nA were 6.6, 7.3, and 9.2 ms, respectively (SDs: 1.7, 1.6, and 1.9).
The LGMD is endowed with an inward rectifying current
Figure 1A (arrowheads) suggests the presence of a small sag in the membrane potential toward the end of the –5-nA current pulse and rebound activity immediately after the pulse. The average difference in membrane potential measured over two 50-ms windows centered 100 and 450 ms after pulse onset was equal to 1.1 mV (SD: 1.4). This change corresponded to a significant difference in the mean membrane potentials over the two windows (t-test, P < 0.05). When pooled across 10 neurons (15 penetrations), the average difference (1.2 mV, SD: 1.7) also corresponded to a significant difference in means (t-test, P < 0.05). To further assess the presence of inward rectification, we injected negative currents of larger amplitude in two neurons (up to –20 nA; three penetrations, seven measurements). Figure 3A illustrates the trial-averaged and median filtered responses to such pulses for one measurement obtained in one of the two neurons. The other measurements in the same neuron (two penetrations total) and in the second neuron showed similar results. The magnitude of the sag increased monotonically with the magnitude of the pulse, which is reflected in the increased difference between the peak hyperpolarization (star in Fig. 3, A and B, top) and the hyperpolarization at the end of the pulse (gray circle in Fig. 3, A and B, top). Because the sag was seen at membrane potentials well below the potassium reversal potential it is more likely to be associated with the opening of a mixed sodium/potassium conductance rather than with closure of a potassium conductance (Halliwell and Adams 1982; Kiehn and Harris-Warrick 1992). The peak rebound value is illustrated in the bottom panel of Fig. 3B (triangles; see also Fig. 3A). The abrupt increase in depolarization at –7 nA (arrow) reflects the reliable generation of a rebound spike by the LGMD at this current value. Thereafter, the number of rebound spikes increased with increasing current magnitude, reaching up to nearly 10 spikes at strongly hyperpolarized potentials (–20 nA, about 100 mV below rest). Similar postpulse rebound spiking was also observed in a separate set of experiments during which five LGMD neurons were filled with Lucifer yellow by injection of large negative current pulses (Peron et al. 2003).
The LGMD is an intrinsically bursting neuron
Next, we examined the LGMD responses to depolarizing current pulses. Figure 4A illustrates sample traces in a single neuron stimulated with pulses of 3, 4, 6, and 10 nA (from bottom to top). The current threshold eliciting spikes was typically around 3–4 nA. At threshold, the LGMD usually fired either one spike (Fig. 4A, bottommost trace) or a burst of two spikes (Fig. 4A, second bottommost trace). At currents above threshold, the LGMD always fired a short burst of spikes followed by isolated spikes (Fig. 4A, middle two traces). As current increased, the number of spikes in the burst typically increased and sometimes a second burst riding over a single depolarizing envelope was observed (Fig. 4A, second topmost trace), followed by isolated spikes. A characteristic signature of intrinsically bursting neurons is a bimodal interspike interval (ISI) distribution, with an initial peak corresponding to intervals between spikes occurring in bursts and a second, broader peak corresponding to intervals between bursts and/or isolated spikes (see, e.g., Figs. 1 and 2 of Nowack et al. 2003). Accordingly, the ISI distributions pooled across all positive current-injection values were bimodal in 11 of 13 neurons tested (17 penetrations), similar in shape to that illustrated in Fig. 4B (bottom). In two neurons, the ISI distributions were unimodal when pooled across all current values. This was the case for the neuron illustrated in Fig. 4 (see B, top). A closer examination of the responses showed that at current values >10 nA, the initial burst smoothly merged with subsequent isolated spikes (see Fig. 4A, top trace), thus causing the peak of the ISI distribution associated with the burst to gradually merge with the tail of the distribution. This observation was confirmed by analyzing ISI distributions for each individual current value. These ISI distributions were initially strongly bimodal and gradually became unimodal as current amplitude increased above 10 nA. Thus, when positive currents were particularly effective at driving the LGMD, the large number of action potentials obtained at high current magnitudes sometimes masked the bursting behavior seen at lower current values. This was confirmed by computing the pooled ISI distribution for current values 10 nA (corresponding to steady-state firing frequencies 50 spikes/s), which were bimodal in both neurons (see Fig. 4B, bottom). We were also able to obtain a second penetration in one of the neurons with unimodal ISI distribution. In the second penetration, positive currents were less effective at driving the cell, presumably because the electrode was located farther away from the spike initiation zone (see the next section) and the ISI distribution was bimodal across all currents. Based on these experimental data, we conclude that the LGMD is an intrinsically bursting neuron, with the capacity for firing at high sustained rates.
Adaptation of the LGMD firing frequency to sustained current pulses
Another feature evident from Fig. 4A is that the firing rate of the LGMD adapts over the course of a current pulse. In Fig. 5A the mean instantaneous firing frequency is plotted as a function of time during current pulses of 8, 12, and 15 nA in one LGMD neuron (black lines). These curves were derived from 10 current-injection trials whose spike rasters are shown below the graph in Fig. 5A. The peak instantaneous firing rate f0 was always attained during the first ISI, increased with current magnitude and started to saturate for values >10 nA (Fig. 5B). The shortest intervals reached typical values of about 1.8 ms, corresponding to peak instantaneous firing frequencies of roughly 550 spikes/s. Thus f0 was a nonlinear, saturating function of current and was well fit by the current–firing frequency curve of a leaky integrate-and-fire model with refractory period and a reset potential different from rest (gray line in Fig. 5B; METHODS, Eq. 1). In contrast, the steady-state firing frequency fss was a linear function of injected current (Fig. 5B). Thus in the LGMD as in other neurons, adaptation results in a linear steady-state current–frequency relationship (Ermentrout 1998; Schwindt et al. 1997; Wang 1998). At all current values, the time course of adaptation was well fit by a single-exponential function (Fig. 5A, gray lines). This yielded a time constant of adaptation adap that increased with the magnitude of the current pulse (adap values of 19, 31, and 42 ms for pulses of 8, 12, and 15 nA, respectively). The attenuation of the steady-state firing rate relative to the peak firing rate was defined as Fadap = (f0 – fss)/f0 and decreased with current pulse magnitude (Fadap values of 0.87, 0.78, and 0.70 for pulses of 8, 12, and 15 nA, respectively). The values of adap and Fadap observed in the LGMD were in the same range as those observed in pyramidal cells of cat visual cortex in vivo (0 < adap < 80 ms and Fadap 
0.50; Ahmed et al. 1998; see also Fig. 6 below).
For cat visual cortex pyramidal neurons, a model of adaptation incorporating a variety of electrophysiological data and based on a calcium-dependent mechanism was previously proposed (Liu and Wang 2001; Wang 1998). According to this model, calcium enters the cell through voltage-gated channels during repetitive firing and activates a calcium-dependent potassium conductance that acts as a negative feedback mechanism, decreasing the cell's firing rate. In the model, calcium is extruded from the cytoplasm with a time constant Ca (see also Eqs. 2 and 3, METHODS). The model explains the exponential relaxation of the firing frequency during a current pulse observed in pyramidal neurons and in the LGMD. Analysis of the model using semianalytical and simulation techniques makes three predictions that have not yet been tested experimentally (for details see Liu and Wang 2001; Wang 1998). The first prediction is that the time constant of adaptation should increase linearly with injected current (top panel in Fig. 3B of Wang 1998). The second prediction is that the attenuation factor Fadap should decrease linearly with injected current (bottom panel in Fig. 3B of Wang 1998). In Wang's model, both of these predicted changes originate from a slower rate of intracellular calcium accumulation as the kinetics of the action potential becomes faster with increasing injected current amplitude. The slower rate of calcium buildup results in turn in longer time constants of adaptation and weaker attenuations. These two predictions imply that the attenuation factor Fadap should depend linearly on the time constant of adaptation adap. The third model prediction is that the inverse slope of this linear relation between adap and Fadap is a biophysical characteristic of the investigated neuron that matches the model's time constant of calcium extrusion Ca. In other words
 | (4) |
; their Eq. 47) and holds approximately in the two-compartment model of Wang (Fig. 6D in Wang 1998). In the latter case, Eq. 4 can be explained by noting that Fadap is nearly proportional to adap. The proportionality constant k (units: ms–1) measures the rate of change of calcium concentration influx through voltage-gated Ca2+ channels as a function of intracellular calcium concentration, as determined by the current injection level and instantaneous firing frequency of the cell (k = Gcc in Eq. 18 of Wang 1998; see also his Fig. 3A, bottom). In this model, however, the same constant k also sets the time constant of adaptation through 1/adap = k + 1/Ca (Eq. 10 in Wang 1998). Combining these two results immediately yields Eq. 4 (which corresponds to Eq. 19 of Wang 1998). Note that Eq. 4 implies that adap and Fadap vary in antagonistic manner, i.e., as adap increases, Fadap decreases, as observed for the three example currents considered above. These three predictions are specific to the models considered by Wang (1998) and Liu and Wang (2001). Other biophysical models of spike-frequency adaptation (see, e.g., Koch 1999) do not necessarily fulfill them (Liu and Wang 2001).
To test these theoretical predictions, we plotted the attenuation Fadap and the time constant of adaptation adap as a function of injected current (Fig. 5, C and D). The graphs were remarkably close to linear (correlation coefficients, = 0.996 and –0.997, respectively). The relation between Fadap and adap was also nearly linear ( = –0.989; Ca = 141 ms), in close agreement with theoretical predictions (Liu and Wang 2001; Wang 1998). In Fig. 6A we show the peak instantaneous and steady-state firing frequency for 13 different LGMD neurons (17 penetrations). Both f0 and fss have similar shapes as in Fig. 5A, but there was quite a range of variability across neurons and penetrations. In particular, the slope of the steady-state firing frequency ranged from 2.5 to 13.7 spikes · s–1 · nA–1. Figure 6, B and C shows the time constant of adaptation and the attenuation for the same recordings as a function of injected current. There was again a wide range of variability across neurons and penetrations. Except in one case (Fig. 6B, arrow), the time constant of adaptation increased linearly with current magnitude and the linear relation predicted theoretically was closely followed (mean = 0.864 across 17 penetrations). Excluding the arrow-marked cell (in Fig. 6B), which had a negative correlation coefficient, yielded a mean = 0.971. The same held true for the attenuation (mean = –0.957 across 17 penetrations). Quite remarkably, when the attenuation was plotted as a function of the time constant of adaptation (Fig. 6D) the variability observed in B and C almost entirely disappeared and data points were tightly clustered around a single, common line. The mean between both variables amounted to –0.868 (17 penetrations) and was even higher (–0.957) when excluding the arrow-marked cell in Fig. 6B. In two neurons, the Fadap–adap relation had a similar slope, but was slightly shifted to the right relative to the bulk of the data (arrow in Fig. 6D). A linear fit to individual Fadap–adap relations yielded a mean intercept of 1.02 (SD: 0.05), very close to one, the value predicted by Eq. 4, and a slope of 131 ms (SD: 23 ms). Thus the slope of the Fadap–adap relation is indeed an invariant biophysical property of the LGMD, independent of the particular animal in which the neuron was recorded and independent of the particular penetration in each LGMD neuron.
These results suggest that the variability observed in Fig. 6, A–C may result from differences in the ability of current injected through the intracellular electrode to stimulate firing in the LGMD. This is to be expected if: 1) the LGMD is not electrotonically compact and 2) the distance of the electrode relative to the spike initiation zone varied across recordings. The first assumption is likely true because the LGMD possesses a complex dendritic morphology (e.g., Fig. 1A of Gabbiani et al. 2002) and visual stimulation mediating inputs at different locations across the main dendrite have differential effects on the firing frequency of the LGMD (Krapp and Gabbiani 2005). The second assumption is also likely, given that in the present recordings the LGMD was penetrated blindly, and thus at random locations of its dendritic tree in the optic lobe. Penetrations in the dendrites at different electrotonic distances of the spike initiation zone are expected to yield spike amplitudes of varying height and width (Stuart et al. 1997). To test these assumptions directly, we took advantage of the fact that we were able to obtain in three different neurons two penetrations that led to significantly different spike heights. In all three cases the decrease in spike height was coupled with an increase in spike width, consistent with the assumption that the recording with smaller spike height was farther away from the spike initiation zone. In all three cases, the recording with the smaller spike height led to a shallower slope of the steady-state firing curve (Table 1). To further test these assumptions across our data set, we plotted in Fig. 6E the steady-state firing frequency slope as a function of spike height measured across 11 neurons (14 penetrations). We found a good correlation between spike height and steady-state firing frequency slope ( = 0.81, slope: 0.20 spikes · s–1 · nA–1 · mV–1). As illustrated in Fig. 6F, the spike height was in turn negatively correlated with spike width ( = –0.60; slope: –0.0051 ms/mV). We conclude that differences in electrode location relative to the LGMD spike initiation zone explain to a large extent the variations observed in Fig. 6, A–C across LGMD neurons and penetrations.
Properties of the postpulse afterhyperpolarization
After a depolarizing current pulse, the LGMD membrane potential exhibited an afterhyperpolarization (AHP) before relaxing over time toward its resting value. This is illustrated in Fig. 7A where the inset shows the response to a single pulse and the following AHP time course (shaded gray area). The main panel shows the time course of the membrane potential on an expanded timescale for a range of current values (1–12 nA) immediately after termination of the pulse. The peak AHP is plotted in Fig. 7B as a function of current amplitude. The relation between these variables was close to linear ( = –0.997) in this and all other cases examined (11 neurons, 15 penetrations, 28 measurements). In the model considered in the previous paragraph, the calcium-dependent potassium current causes both spike-frequency adaptation and the postpulse AHP (see also METHODS, AHP DECAY IN THE LIF ADAPTATION MODEL). If this is also predominantly the case for the LGMD, one expects the potassium current to decay after a current pulse with a time constant independent of the recorded neuron and equal to Ca. Consequently, an approximately similar decay is also expected for the AHP (e.g., Sah 1992; Schwindt et al. 1988). To test this hypothesis, we replotted the normalized AHP time course in cases where the absolute peak AHP value exceeded 2 mV (Fig. 7C, gray traces; obtained from those in Fig. 7A). The time course of AHP decay was similar across current-injection amplitudes and the normalized AHP averaged across current values (dotted line in Fig. 7C) could be fitted with a single exponential (Fig. 7C, black line), although these fits were typically less good than those obtained for spike-frequency adaptation. In Fig. 7D we compare the time constant of AHP decay to that obtained from the slope of the Fadap–adap relation in 11 neurons. The mean AHP decay time constant was equal to 114 ms (SD: 37) and was slightly lower than the mean adap–Fadap time constant (131 ms, SD: 23). There was no correlation between these two variables ( = –0.03), indicating no systematic changes across cells. These results are thus consistent with a common mechanism governing spike-frequency adaptation and AHP decay in the LGMD. Other conductances, such as inward rectification, might also contribute to a lesser extent to AHP decay, thus speeding up its dynamics (Lorenzon and Foehring 1992).
Interspike interval characteristics at steady-state
As illustrated by the raster plots shown in Fig. 5A (top), firing of the LGMD was very regular at steady state for high firing frequencies. Correspondingly, the coefficient of variation (CV) of the ISI distribution at steady state was between 0.1 and 0.2 for firing frequencies in excess of 50 spikes/s (Fig. 8A). As the steady-state firing frequencies decreased <50 spikes/s, the variability of ISIs rapidly increased, with CV values between 0.2 and 0.6. When the CV is high at low firing frequencies, the theoretical models of Wang (1998) and Liu and Wang (2001) predict negative correlations between ISIs because randomly occurring short intervals cause sufficient calcium entry to activate the AHP conductance, making the next ISI more likely to be longer and vice versa (Fig. 9C of Wang 1998; Fig. 8 of Liu and Wang 2001). In Fig. 8B we show the ISI correlation coefficient (CC) as a function of the CV at steady state. There was a negative correlation ( = –0.45) between these two variables and the CC was indeed negative at high CVs. In contrast, when strong current pulses were dominant in driving the LGMD responses, low CVs resulted in predominantly positive CCs. This is expected for strongly periodic firing from results obtained in simplified neuronal models in the presence of weak, correlated noise (Lindner 2004).
|
The linear dependency of the adap–I, Fadap–I, and Fadap–adap relations observed in the LGMD were predicted in a conductance-based, two-compartment model (Wang 1998) and in an LIF neuron (Liu and Wang 2001), which both included a calcium-dependent adaptation current and an exponential calcium extrusion mechanism. Because an LIF neuron can be constrained from the electrophysiological data described in the previous paragraphs, we investigated whether this model is able to replicate the adaptation properties of the LGMD. Such a model cannot of course replicate more complex properties such as bursting, but has the advantage of being conceptually simple. It could thus potentially yield insight in the firing properties of the LGMD for stimuli other than current pulses. To be useful, such a model needs to work over a large fraction of the LGMD dynamic range and in particular at high firing frequencies, as observed during stimulation of the LGMD with visual looming stimuli. The inset in Fig. 9 shows a simulation for model parameters matched to the properties of the neuron depicted in Fig. 5 (current pulse amplitude: 8 nA). The firing frequency of the LIF adapts over the course of the current pulse and the model exhibits attenuation and an AHP similar to those observed in the LGMD. An important difference between the model used here and the one studied by Liu and Wang (2001) is that we had to introduce a refractory period (tref = 1.5 ms) to replicate the nonlinear relation between the peak instantaneous firing rate frequency f0 and current over a large fraction of the LGMD firing range (Fig. 5B, gray trace). It follows immediately from Eq. 1 that the firing frequency of an LIF neuron in the presence of a refractory period is given by f(I) = [tref + fnoref(I)–1]–1, where fnoref(I) is the firing frequency with tref set to zero. The refractory period thus plays a negligible role for low currents [tref << fnoref(I)–1], leading to an exponential decay of firing frequency over the time course of the pulse (Fig. 9, bottom black trace and gray line). For higher currents, however, the refractory period is initially predominant and sets its own timescale for firing frequency decay, thus interfering with the negative feedback mechanism mediated by the adaptation current. Accordingly, a single exponential could not fit the decay of the firing frequency observed in the model (Fig. 9, top two black traces and gray lines), in contrast to experimental findings (Fig. 5A). Thus an LIF neuron with refractory period is unable to fit the LGMD adaptation properties over the range of firing frequencies studied experimentally. We conclude that a conductance-based, compartmental model will be required for this purpose.
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONGRANTSACKNOWLEDGMENTSREFERENCES |
|---|
). Adaptation followed remarkably well the predictions made by the two-compartment model of Wang (1998). By using a uniquely identifiable neuron, we were in particular able to verify that the relation between attenuation and the time constant of adaptation is linear and that its slope is a biophysical invariant of the LGMD. This prediction would have been very difficult to test on a heterogeneous neuronal population. Our results thus validate the model's predictions and also rule out a description of the LGMD by a simpler, leaky integrate-and-fire model.
Because we studied the electrotonic properties of the LGMD in vivo, our measurements were presumably affected by more than pure passive membrane properties. Although we did not characterize spontaneous membrane potential fluctuations, intracellular LGMD recordings in vivo typically exhibit a sizable number of spontaneous excitatory and inhibitory postsynaptic potentials. Thus background activity is likely to have affected the input resistance and the membrane time constant of the cell (Bernander et al. 1991). It might also have had an impact on the resting membrane potential. Some active conductances such as the one mediating inward rectification may also be tonically active and contribute to resting membrane properties (Magee 1998). The increase in membrane time constant with hyperpolarizing pulse amplitude would, for example, be consistent with rapid closing of tonically active conductances. Because the input resistance did not change significantly over longer timescales, this might have been compensated for by a slower activation of inward rectification, for instance. The estimation of electrotonic parameters may also depend on the recording technique. In vitro, differences were observed between sharp electrode recordings (used here) and patch-clamp recordings (Spruston and Johnston 1992; Staley et al. 1992), but may be less conspicuous in vivo (Borg-Graham et al. 1996). Currently, it is unknown whether such differences exist in insect visual interneurons because in vivo patch-clamp recordings have proven difficult to obtain from large wide-field cells such as the LGMD. Nonetheless, the values reported here were similar to those observed in fly tangential neurons obtained under similar circumstances (Borst and Haag 1996). These neurons, like the LGMD, integrate on their extended dendritic arborizations a large number of local visual inputs. In vivo electrotonic properties are directly relevant to information processing in the context of natural visual stimuli. In particular, the low input resistance and fast membrane time constant of the LGMD allow for rapid processing of visual information. As is the case for fly tangential neurons involved in visual flight control (e.g., Borst and Haag 2002; Egelhaaf et al. 2002), rapid visual information processing is critical to support successful collision avoidance and escape behaviors.
Inward rectification has been reported in many vertebrate and invertebrate neurons (e.g., Golowasch and Marder 1992; Kiehn and Harris-Warrick 1992; Lüthi and McCormick 1998). Originally, inward rectifying conductances have been linked to bursting (Jahnsen and Llinas 1984). More recently, inward rectification was shown to shorten the window of synaptic integration (Magee 1998; Migliore et al. 2004; Poolos et al. 2002). The role played by inward rectification in the LGMD during synaptic integration of visual inputs is under study. Inward rectification may also play a role in shaping the bursting properties of the LGMD. In our experiments, t
