Abstract
Cortical synapses exhibit several forms of short-term plasticity, but the contribution of this plasticity to visual response dynamics is unknown. In part, this is because the simple patterns of stimulation used to probe plasticity in vitro do not correspond to patterns of activity that occur in vivo. We have developed a method of quantitatively characterizing short-term plasticity at cortical synapses that permits prediction of responses to arbitrary patterns of stimulation. Synaptic responses were recorded intracellularly as EPSCs and extracellularly as local field potentials in layer 2/3 of rat primary visual cortical slices during stimulation of layer 4 with trains of electrical stimuli containing random mixtures of frequencies. Responses exhibited complex dynamics that were well described by a simple three-component model consisting of facilitation and two forms of depression, a stronger form that decayed exponentially with a time constant of several hundred milliseconds and a weaker, but more persistent, form that decayed with a time constant of several seconds. Parameters obtained from fits to one train were used to predict accurately responses to other random and constant frequency trains. Control experiments revealed that depression was not caused by a decrease in the effectiveness of extracellular stimulation or by a buildup of inhibition. Pharmacological manipulations of transmitter release and postsynaptic sensitivity suggested that both forms of depression are mediated presynaptically. These results indicate that firing evoked by visual stimuli is likely to cause significant depression at cortical synapses. Hence synaptic depression may be an important determinant of the temporal features of visual cortical responses.
- facilitation
- depression
- temporal integration
- temporal filtering
- contrast adaptation
- presynaptic
- neocortex
- nonlinearity
- EPSC
- field potential
Several decades of concerted effort have yielded important insights into the mechanisms that underlie spatial-filtering properties of primary visual cortical (V1) neurons (for review, see Shapley and Lennie, 1985; Ferster, 1994; DeAngelis et al., 1995) but have revealed much less about the mechanisms that underlie their temporal-filtering properties. Temporal-filtering properties are important because of their role in adaptation and gain control (Bonds, 1991; Nelson, 1991a; Abbott et al., 1997), velocity tuning (Orban et al., 1985), and direction selectivity (Reid et al., 1991; Jagadeesh et al., 1993). Recently, we proposed that some of the specific temporal-filtering properties of visual cortical neurons reflect short-term plasticity at cortical synapses (Abbott et al., 1997; Nelson et al., 1997; for related ideas, see Grossberg, 1984;Thomson and Deuchars, 1994; Finlayson and Cynader, 1995; Fisher et al., 1997; Lisman, 1997; Tsodyks and Markram, 1997). Short-term plasticity at neocortical synapses is typically studied by applying pairs or brief constant frequency trains of electrical impulses (Sutor and Hablitz, 1989; Thomson and West, 1993; Thomson et al., 1993; Metherate and Ashe, 1994; Castro-Alamancos and Connors, 1996; Hess et al., 1996; Stratford et al., 1996). This approach has been invaluable for understanding the mechanisms of short-term plasticity but does not provide a general means of predicting responses during more complex patterns of activity, such as those that occur in vivo. Developing such a method is important for understanding the role these synaptic processes play in the sensory responses of cortical neurons.
Single-unit studies have revealed that the temporal properties of V1 responses can differ substantially from those of neurons in the lateral geniculate nucleus (LGN). In general, V1 neurons respond to lower temporal frequencies and slower velocities than do LGN neurons (Movshon et al., 1978; Orban et al., 1985; Hawken et al., 1996). However, unlike linear low-pass filters, cortical neurons typically respond briskly to transients that contain high-frequency components, although they often have smaller sustained responses (Kulikowski et al., 1979). In addition, the responses of cortical neurons show much more pronounced adaptation during prolonged stimulation than do LGN neurons. Rapid adaptation effects, recovering over the course of a fraction of a second (Bonds, 1991; Nelson, 1991a), and more prolonged adaptation effects, recovering over the course of many seconds (Movshon and Lennie, 1979; Ohzawa et al., 1982), have been observed. Finally, heterogeneity of temporal-filtering properties across the receptive fields of some cortical neurons may play an important role in generating direction selectivity (Reid et al., 1991; Jagadeesh et al., 1993).
Classical studies of synaptic transmission at the neuromuscular junction identified short-term enhancement and depression occurring over several time scales (for review, see Magleby, 1987; Zucker, 1989;Fisher et al., 1997). Enhancement has been particularly well characterized. Four temporally distinct components of response enhancement (facilitation 1, facilitation 2, augmentation, and potentiation) have been identified first at neuromuscular junctions and later at central synapses and have been shown to result from enhanced release of neurotransmitter that persists after presynaptic calcium entry (Swandulla et al., 1991; Regehr et al., 1994).
The time course and mechanisms underlying depression have not been as completely characterized, but depression is believed to result from depletion of a readily releasable pool of neurotransmitter (Betz, 1970;Kusano and Landau, 1975). Depression seems to be a particularly prominent feature of transmission at neocortical synapses (Deisz and Prince, 1989; Nelson and Smetters, 1993; Thomson et al., 1993; Markram and Tsodyks, 1996). Depression and facilitation can coexist, and the balance between the two depends strongly on quantal content. Conditions under which larger numbers of quanta are released per action potential (or probability of release is high) favor depression, whereas conditions under which smaller numbers of quanta are released per action potential (or probability of release is low) favor facilitation (for review, see Zucker, 1989). The objective of the present study was to measure the temporal transfer characteristics of cortical synapses in sufficient detail to permit prediction of responses to arbitrary patterns of stimulation.
MATERIALS AND METHODS
Coronal slices containing primary visual cortex were obtained from Long Evans rats, age postnatal days 14–40. Animals were deeply anesthetized with ketamine (100 mg/kg) and acepromazine (10 mg/kg) or with pentobarbital (35 mg/kg) and decapitated, and their brains were quickly removed and placed in chilled (5°C) artificial CSF (ACSF). Slices of 400 μm thickness were cut on a vibratome. During recording, slices were transilluminated, permitting visualization of the location of primary visual cortex and the boundaries between layers 2/3 and 4 (Domenici et al., 1995).
Slices were maintained at room temperature on semipermeable membranes (Falcon 3090) covered by a thin layer of ACSF continuously oxygenated with 95% O2/5% CO2. They were transferred one at a time to a submerged chamber mounted on a fixed-stage upright microscope (Nikon Optiphot UD) and were slowly warmed to 34–35°C. Slices were equilibrated for 1–2 hr before recording and remained viable for up to 16 hr. Slices were perfused with warmed, oxygenated ACSF at a rate of 2–3 ml/min. For obtaining visually guided whole-cell recordings, slices were illuminated obliquely through an infrared filter and viewed with standard optics using a 40× long working-distance water immersion objective. The resulting image was displayed on a video monitor using a CCD camera.
Solutions. ACSF contained (in mm): 126 NaCl, 3 KCl, 1.25 NaH2PO4, 10 dextrose, 20 NaHCO3, 2 MgSO4, and 2.0 CaCl2. pH was 7.4 when saturated with 95% O2/5% CO2, and osmolarity was 305–310 mOsm.
Pipettes were pulled from 1.0-mm-outer diameter thin-walled capillary tubing (Warner Instruments) on a Flaming–Brown horizontal puller (Sutter) and were filled with (in mm): 130 potassium methylsulfonate, 10 KCl, 10 HEPES, 0.5 EGTA, 2–3 Na2ATP, 0–1 GTP, and 2 MgSO4. Osmolarity of the internal solution was 280–290 mOsm, and pH was 7.2–7.4.
Electrophysiological techniques. Whole-cell recording pipette resistances were 4–8 MΩ in the bath. Voltage-clamp recordings were performed using an Axopatch 1D (Axon Instruments). Seal resistances were 1–10 GΩ. Unless stated otherwise, cells were voltage-clamped at −70 mV. Series resistances (5–20 MΩ) were uncompensated. Voltages were not corrected for liquid junction potentials. EPSCs were judged to be monosynaptic if they occurred with short (1.5–4 msec) and constant (jitter <1 msec) latency that did not change with small changes in stimulus strength. EPSCs were included only if additional inhibitory components were not present after depolarization (40–50 mV above rest).
Field potential recordings were obtained from low resistance (1–2 MΩ) saline-filled pipettes. Biphasic electrical stimuli (80 μsec in duration; 10–150 μA) were delivered via saline-filled pipettes. Electrode placements were performed under low power using transillumination.
Signals were filtered at 1–1000 Hz (field potentials) or DC 2 kHz (EPSCs) and digitized at 5 kHz (Instrutech). Data were collected and analyzed using IGOR (Wavemetrics, Oswego, OR), Pulse Control (Herrington et al., 1995), and additional fitting routines implemented in the C programming language.
Data analysis. Response magnitudes were measured as peak amplitudes within a 1 msec window. Amplitudes, rather than slopes of field potentials, were measured because in many cases the initial slope was contaminated by a nonsynaptic response (Hess et al., 1996). We were concerned initially that truncation of the peak response by overlapping inhibition might significantly influence the kinetics of the measured plasticity. For example, apparent depression of excitatory responses could actually reflect potentiation of inhibition. Several observations indicate that this was not the case. (1) In five slices in which the synaptic and nonsynaptic responses were fully separated, we confirmed that measurements of initial slope and peak amplitude were well correlated (mean Pearson correlation coefficient, r = 0.93l). (2) In four slices, we confirmed that similar synaptic depression was observed in the presence of 20 μmbicuculline (see Fig. 2C). (3) Previous study of monosynaptic IPSCs in neocortex (Deisz and Prince, 1989) and our own preliminary observations (Song et al., 1997) indicate that during repeated stimulation IPSCs also depress, rather than potentiate. (4) Nearly identical depression was observed in recordings of field potentials and of EPSCs. For the intracellular recordings, we used much lower amplitude stimuli and confirmed that no IPSCs were present.
Initial attempts to extract the underlying dynamics of these responses were based on the synaptic decoding method of Sen et al. (1996). This technique is quite generic and makes few assumptions about the temporal profile of short-term changes in synaptic strength. Because our goal was to predict synaptic responses accurately with as concise a description as possible, we subsequently adopted a simpler, but somewhat less generic, model in which changes in response amplitude (A) were assumed to result from the product of an initial amplitude (A0) and dynamic variables representing facilitation and depression. The choice of this type of model was guided by our observation that kernels obtained using the synaptic decoding method were usually well fit by exponentials. A similar class of models has been used previously to describe short-term plasticity at the neuromuscular junction (Liley and North, 1952;Magleby and Zengel, 1976).
Response amplitudes were fit with several related models that differed in the number of facilitation (zero or one) and depression (one to three) factors. The most complex form of the model contained a single facilitation factor (F) and three depression factors (D1, D2, and D3): Equation 1For other less complete variants, one or more of the facilitation and depression factors were omitted (i.e., were constrained to be 1): Equation 2 Equation 3 Equation 4 Equation 5 Equation 6
Dynamic variables representing depression (D) were constrained to be ≤ 1 and depended on the stimulus pattern in the following way: after each stimulus in the train, Dwas multiplied by a constant factor (d) representing the amount of depression per presynaptic action potential: Equation 7Between stimuli, D recovered exponentially back toward 1 with first-order kinetics and time constant τD: Equation 8Accumulation of depression results when the interval between presynaptic action potentials is less than that required for recovery. The different components of depression (D1, D2, and D3) had different constant factors (d1, d2, and d3) and time constants (τd1, τd2, and τd3).
The dynamic variable F representing facilitation was constrained to be ≥ 1. After each stimulus, a constant,f (≥ 0), representing the amount of facilitation per presynaptic action potential was added to F, and between stimuli, F recovered exponentially back toward 1: Equation 9 Equation 10A model of facilitation in which F was updated additively, rather than multiplicatively, was found to provide a superior fit to the data, especially at high stimulus frequencies where multiplicative facilitation can lead to amplitudes that increase without bound.
All of the data were fit with the two-component depression model defined by Equation 4 and the complete form of the model defined by Equation 1. Most of the data were also fit with the model variant containing two depression components and one facilitation component (Eq. 5). This three-component model, which included facilitation, was used for all of the fits illustrated (see Figs. 3, 4, 8, 10, 11, 12, 13), because it was the simplest form of the model that adequately described data obtained under conditions of high and low probability of release (see Fig. 7). For comparing data across different recordings obtained under control conditions (see Fig. 9), the two-component model was used (Eq. 4), because it adequately described that data and had fewer parameters.
Fits were evaluated first by measuring the fractional error [(observed − predicted)/observed] for each response in a train and then by calculating its average (referred to as average error) and its root mean square (referred to as rms error). The average error is a measure of the degree to which the predictions systematically underestimate or overestimate the data, whereas the rms error is a measure of the closeness of fit. The set of model parameters yielding the lowest rms error was found by an automated search algorithm. We also calculated an error index, defined as the ratio of the rms error obtained with the best fit of the model to the rms error obtained when no depression or facilitation factors were included (i.e., whenA = A0).
Typically we fit a set of responses obtained during stimulation with a Poisson-distributed train or with a Poisson-distributed train followed by a period of constant frequency stimulation, and then used the model parameters to predict responses to other random and constant frequency trains. To avoid overlap of individual responses within a train, we clipped the minimum interval at 30 msec. In some experiments, we used smaller minimum intervals. Individual trains were separated by 45–60 sec. We ensured that this was an adequate period of recovery by confirming that the initial amplitudes did not vary systematically on subsequent repetitions (see Fig.1A,C, top).
RESULTS
Synaptic responses in layer 2/3 exhibit prominent synaptic depression
Figure 1 illustrates synaptic responses evoked in layer 2/3 by Poisson-distributed stimulus trains applied to layer 4 immediately below the recording site. Both monosynaptic EPSCs (Fig.1A,B;n = 11 cells), recorded intracellularly, and field potentials (Fig. 1C,D; n = 38 slices), recorded extracellularly, strongly depressed during the course of each stimulus train and recovered during the 45–60 sec period between repetitions of the stimulus train. The amplitude of individual responses within the train typically varied over a twofold to eightfold range (Fig. 1A,C,bottom). Most of this variance was not caused by random fluctuations, particularly for the field potentials, because the trial-to-trial variance of the responses (Fig.1A,C, top; SEM/mean, 8.1% for EPSCs and 4.0% for field potentials) was small compared with the variance of the responses across stimuli within the train. In general, the trial-to-trial variance of the EPSCs evoked by near-minimal stimulation was much greater than the trial-to-trial variance of the field potentials, which were evoked using larger amplitude stimuli. Presumably, EPSCs are more variable because the small number of synapses contributing to the responses are not sufficient to average out random fluctuations in transmitter release. In contrast, field potentials presumably reflect simultaneous activation of large numbers of synapses and hence allow an accurate determination of the mean response to a particular stimulus pattern using a much smaller number of trials.
Although it is widely assumed that field potentials in hippocampus and neocortex reflect primarily activation of excitatory synapses (Mitzdorf and Singer, 1978; Bode-Greuel et al., 1987; Langdon and Sur, 1990;Kirkwood and Bear, 1994), we wished to confirm this fact in our preparation. Two components of the response were distinguished on the basis of latency and ability to follow high-frequency stimulation: (1) a very short latency (0.5–1.5 msec onset; 1–2 msec peak) response that persisted during high-frequency stimulation, and (2) a slightly longer latency response (2–4 msec onset; 4–7 msec peak) that depressed rapidly during high-frequency stimulation. Additional, longer latency negative components were sometimes present but were not studied in detail. Based on these observations and on the effects of pharmacological manipulations described below and in keeping with previous observations in visual cortex (Kirkwood and Bear, 1994) and other cortical areas (Hess et al., 1996), we identified these components as (1) a nonsynaptic fiber volley and antidromic response (referred to as the antidromic response) and (2) a field EPSP (fEPSP). Application of the AMPA glutamate receptor antagonist 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX; 10–20 μm) completely blocked the synaptic component of the response (Fig.2A; n = 6) but spared the antidromic response. Monosynaptic IPSPs and NMDA receptors did not contribute to the potential that remained after application of CNQX, because subsequent addition of bicuculline methiodide (BMI; 20 μm) and aminophosphonovalerate (APV; 50 μm) had no effect (data not shown). Addition of tetrodotoxin (TTX; 1 μm) (Fig. 2A) abolished both the synaptic and antidromic responses. We often observed an upward deflection after the initial negative-going synaptic response. We suspect that this reflects activation of polysynaptic inhibitory pathways because (1) it corresponds temporally to disynaptic inhibitory currents recorded intracellularly; (2) during application of BMI it was reduced or abolished and was replaced by additional negative components that corresponded temporally with polysynaptic excitatory potentials recorded intracellularly; and (3) it tended to fatigue even more rapidly than the short latency response.
During repetitive stimulation, the antidromic response, unlike the synaptic response, did not decline in amplitude. In some cases, like that shown in Figure 2B, there was a small increase in amplitude after the first or second stimulus that then decayed back to the initial levels. In the majority of cases, the amplitude of the antidromic response did not change. This indicates that the depression of the synaptic response is not simply an artifact of decreasing effectiveness of extracellular stimuli in a train.
To avoid large polysynaptic potentials and epileptiform discharges that often occurred after stimulation in the presence of bicuculline, we did not block inhibition in most experiments. During whole-cell recordings, EPSCs were isolated by keeping stimulus strengths below the threshold for evoking an IPSC. In four field potential experiments, we were able to record a stable pharmacologically isolated AMPA receptor component of the fEPSP in the presence of 20 μm BMI to block inhibition and 2 μm MK-801 to block NMDA receptors without evoking epileptiform activity. Synaptic depression observed during blockade of GABAA receptors and NMDA receptors was extremely similar to that observed under control conditions (Fig.2C).
Responses to random stimulus trains can be fit accurately by a simple model of short-term facilitation and depression
Figure 3 illustrates the best fit of the three-component model (one facilitation component and two depression components, see Eq. 5) to EPSCs recorded from a layer 2/3 pyramidal neuron in response to 10 repetitions of a Poisson-distributed train. The model accurately captures the overall time course of depression but is less precise in predicting the response to individual stimuli. Averaged across the responses within the train, the difference between measured and predicted amplitudes was not significantly different from zero (average error, 2.5 ± 15.9%), although predictions of responses to individual stimuli varied appreciably in how well they matched the data (rms error, 16.0%). The error index was 8.6%, indicating that the prediction accounted for >91% of the variation in amplitude across stimuli within the train. The same parameters used to fit the response to the Poisson train were then used to predict responses to constant frequency trains at 5 and 10 Hz (Fig.3C,D). The predicted responses to constant frequency trains agreed well with the measured responses (rms errors of 10.1 and 16.7%, respectively), although, as in the initial fit, responses to individual stimuli deviated in an apparently random way from predicted values. Similar fit parameters were obtained for EPSCs recorded from 8 of 11 cells. In two of the remaining cells, the data were better fit by a single exponentially decaying depression (i.e.,F and D2 both equaled 1), and in the third remaining cell, the data were best fit by a combination of facilitation and a single depression factor (i.e.,D2 was 1).
Much more accurate (i.e., lower rms error) fits were obtained from the field potential data. Presumably the increased accuracy reflects the lower trial-to-trial variability of the field potentials that allow more accurate estimates of the mean response to each stimulus in a train, despite a small number of repetitions (4–5; see Fig. 5 and below). An example of a fit to field potential recordings is shown in Figure 4. Fine details of the fluctuations in response amplitude are accurately captured by the fit. The rms error for this fit was 8.3%, and the error index was only 2.4%, indicating that the prediction accounted for >97% of the variation in amplitude across stimuli within the train. The parameters of the fit (Fig. 4E) were quite similar to those obtained from fits to EPSCs. The same parameters extracted from the fit to the 4 Hz Poisson train (Fig. 4A) accurately predicted responses to constant frequency trains (Fig.4C,D). For these trains, the rms errors were 2.5 and 8.5%, and the error indices were 1.3 and 1.9%. The same parameters also accurately predicted responses to different random trains. The rms errors for additional random trains with mean rates of 1 and 2 Hz were 11.8 and 10.9%, respectively.
To determine whether the fit errors were determined primarily by limitations of the model or by the intrinsic variability of the data, we plotted the precision of the fits against the variability of the data. Figure 5A illustrates the fact that the error of the fits increased with increasing trial-to-trial variability. Each point plots the rms error for a stimulus train against the SEM that, like the rms error, has been expressed as a fraction of the response amplitude and averaged across responses within a train. All data from 2 and 4 Hz Poisson-distributed trains and 1, 5, 10, and 20 Hz constant frequency trains are included. Three to five stimulus trains are shown for each cell and extracellular recording site. In each case, one or two trains were used to generate the fit, and two to four were predicted using the same parameters. The values shown are from fits of the full four-component model (one facilitation and three depression components, Eq. 1), but fits of the three-component model to a subset of the data were quite similar (see Fig. 7). The high degree of correlation between the rms error of the fit and the SEM of the data (r = 0.90) suggest that trial-to-trial variability was an important factor limiting the accuracy of predictions of individual responses within a train. The intercept represents the accuracy of our fits in the limit of zero trial-to-trial variability. An alternative measure of fit accuracy is the residual error remaining after the error attributable to variability is subtracted from the rms error (Fig. 5A,inset). The two measures of fit accuracy independent of intrinsic variability yielded similar values (intercept = 4.7%; residual error mean ± SD, 1.8 ± 5%; median = 3.2%). As seen in the case of the individual set of responses shown in Figure4, errors for response predictions were typically no worse than the errors of the initial fits (Fig. 5B).
Comparison of related models
In most cases, the data were fit better by a model with two depression factors than by a model with only a single depression factor. To test more directly the need for two components of depression with differing time courses, we measured the recovery from depression after a train of 15 stimuli at 20 Hz. As shown in Figure6, the amplitude of field potential responses to a single test pulse given at various times after the end of the train recovered with a time course that was much better fit by the product of two exponentials than by a single exponential. This was true in each of the five experiments in which the recovery time course was measured.
To compare more comprehensively the need for each of the components of the model, we fit a subset of the data with six different related models (Eqs. 1-6). For the sake of comparison, the rms error obtained when assuming no short-term plasticity (A =A0; see error index in Materials and Methods) was 157 ± 28.1%. Under control conditions, the dominant form of short-term plasticity is depression, hence the data were poorly fit by a single facilitation factor but were better fit by a single depression factor. In keeping with the results of the recovery experiments, including a second depression factor produced a substantial improvement in the fit (from rms error of 11.9 ± 1.0% with a single depression factor to 9.8 ± 0.6% with two depression factors). We also assessed the benefit of adding a third depression component having a faster time course (10–100 msec) (Thomson et al., 1993; Stevens and Wang, 1995). Adding this third, very rapidly recovering depression factor produced little improvement for models that lacked facilitation but led to small additional improvements in models that included a facilitation term. It should be noted that the data used for Figure7A may underestimate the magnitude of this faster form of depression because the data were collected using stimulus trains in which the range of frequencies was not truly Poisson but was clipped at 30 Hz. When shorter intervals were included in the stimulus trains (clipping at 60 Hz, n = 4; clipping at 100 Hz, n = 3; no clipping frequency,n = 3), the data were better fit with a model that included this rapid depression component (rms error of 15.3 ± 1.2% for three depression factors vs 16.9 ± 1.9% for two depression factors).
Mixtures of facilitation and depression typically produced the best fits. The improvement over models containing only depression was modest under control conditions (e.g., rms error of 8.7 ± 0.8% for two depression factors and one facilitation factor vs 9.8 ± 0.6% for the two depression factors with no facilitation; Fig. 7A) but became more significant under conditions of reduced transmitter release (Fig. 7B,C). This reflected the reduced depression and enhanced facilitation that was apparent under conditions of reduced release (see below and Figs. 12, 13, 14).
To determine whether or not there was a unique best fit of the model, we varied fit parameters over a wide range and measured the resulting error. For these studies, we used the three-component variant of the model containing a single facilitation component and two depression components (Eq. 5). Varying each parameter individually (Fig.8A,B) produced smoothly changing error functions that had clear minima and were approximately parabolic in shape. The error was more sensitive to variations in some parameters (e.g., d2; Fig. 8A, filled circles) than in others (e.g., τd2; Fig. 8B,filled circles). Varying pairs of parameters simultaneously (Fig. 8C,D) produced error surfaces that were also smoothly varying and that each had a single discrete region of minimum error. The region of minimum error for a given pair of parameters (e.g., f and τd1; Fig. 8C) was typically elongated along a diagonal. For most pairs of parameters examined, the degree of elongation was similar to that shown in Figure 8C. When the degree of facilitation and the faster depression were varied simultaneously (Fig. 8D), the degree of elongation was substantially greater than with any other pair of parameters. Along this diagonal, increases in facilitation could be matched by decreases in depression (or vice versa) without significantly affecting the quality of the fit. This ambiguity is present in all model variants that contain opposing facilitation and depression acting over approximately similar time scales. We performed the same analysis on fits to data from four other slices and obtained very similar error surfaces.
To compare the amount of depression and its rate of recovery across experiments, it was important to minimize the degeneracy of the fits. Data were fit with a variant of the model that included the two depression components described above but in which facilitation was omitted (Eq. 4). This represented a compromise between using a single-component model that in some cases produced rather poor fits but had the fewest free parameters and using the full four-component model that produced better fits but had nine free parameters (A0 plus the amplitudes and time constants of each component). The parameters of the best fits of this model to all of the data recorded under control conditions are plotted as a function of age of the animal in Figure 9. Both forms of depression were evident over the entire range of ages studied and were similar in intracellular and extracellular recordings. There was a tendency for the recovery time constants, especially for the slower of the two forms of depression, to be shorter in the older animals.
Manipulations that alter postsynaptic responsiveness do not alter short-term plasticity
In principle, the prominent short-term synaptic depression we observed could arise as a result of a change in the amount of neurotransmitter released from presynaptic axons (Betz, 1970;Kusano and Landau, 1975; Larkman et al., 1991; Thomson et al., 1993;Stevens and Wang, 1995; Debanne et al., 1996), a change in the sensitivity of postsynaptic cells to that neurotransmitter (Hestrin, 1992; Trussell et al., 1993), or a change in other conductances, such as voltage-dependent conductances (Sutor and Hablitz, 1989; Deisz et al., 1991; Hirsch and Gilbert, 1991; Thomson et al., 1993) activated as a result of the depolarization caused by synaptic current. In an attempt to distinguish between these possibilities, we applied a number of pharmacological manipulations known to affect presynaptic release and postsynaptic responsiveness. We reasoned that if the relevant change were occurring in voltage-dependent conductances activated by postsynaptic depolarization, short-term plasticity should be altered by manipulations that decrease the postsynaptic response.
We compared the best fit of the model for data obtained under control conditions with that obtained at the same recording site after application of a low concentration (0.25 μm) of the AMPA-type glutamate receptor antagonist CNQX. As shown in Figure10, CNQX dramatically decreased the overall amplitude of the responses but had little effect on the observed short-term plasticity. This suggests that changes in conductances that depend strongly on postsynaptic voltage are unlikely to make a major contribution to the observed short-term plasticity. Similar results were obtained in four other experiments. These results were also consistent with the observation that the magnitude and time course of depression were similar during whole-cell voltage-clamp experiments when cells were kept hyperpolarized and much smaller synaptic currents were evoked.
A second possible postsynaptic mechanism that might contribute to short-term synaptic depression is desensitization of AMPA receptors, which occurs at a variety of central synapses including those in rat visual cortex (Hestrin, 1992). To investigate this possibility, we applied cyclothiazide (CTZ), a drug that reduces AMPA receptor desensitization (Trussell et al., 1993), and measured its effect on short-term synaptic plasticity. CTZ (60 μm;n = 5 experiments) markedly increased the duration and slightly increased the amplitude of individual field potential responses (Fig. 11, inset) but had no effect on short-term plasticity (Fig.11A–C).
Manipulations that alter transmitter release alter short-term plasticity
At neuromuscular junctions (del Castillo and Katz, 1954;Betz, 1970), the squid giant synapse (Kusano and Landau, 1975;Swandulla et al., 1991), and synapses in spinal cord (Pinco and Lev-Tov, 1993), hippocampus (Debanne et al., 1996), and neocortex (Thomson et al., 1993; Markram and Tsodyks, 1996), presynaptic changes in quantal content or probability of release have been shown to influence profoundly short-term facilitation and depression. To determine whether the forms of short-term plasticity studied here also depended on presynaptic release, we applied manipulations that have been shown previously to affect release. Presynaptic release was decreased by lowering external calcium from 2.0 to 0.5 mm(with a compensatory increase in external Mg2+ ion concentration; n = 4; Fig.12) and by applying neuromodulatory substances such as adenosine (Prince and Stevens, 1992) (n = 10; Fig. 13) or baclofen (Thompson et al., 1993) (n = 7; Fig.14) that have been shown previously to decrease release. Although adenosine and baclofen can also have postsynaptic effects, we have observed previously (Varela et al., 1995; J. A. Varela and S. B. Nelson, unpublished observations) that at visual cortical synapses the presynaptic effects are far more pronounced and occur at lower concentrations. As shown in Figures 12and 13, reducing release dramatically altered short-term plasticity. Responses recorded during application of lowered calcium were smaller and, in addition, showed much less overall depression. During high-frequency bursts within the random stimulus train, net depression was converted to net facilitation (Fig. 12, arrow). Responses recorded during application of adenosine or baclofen, in the presence of normal extracellular calcium, showed a similar reduction in the degree of depression and enhancement of net facilitation (Fig.13).
In a smaller number of experiments, we increased release by elevating external calcium (n = 4) or by applying the potassium channel blocker 4-aminopyridine (4-AP; n = 4). 4-AP enhances neurotransmitter release at many types of synapses, presumably by causing action potential broadening and hence increased calcium entry (Rutecki et al., 1987). These manipulations produced modest increases in the amplitude of responses to individual stimuli and enhanced depression.
The effects of the various pharmacological manipulations on short-term plasticity are summarized in Figure 14, which shows responses to 20 Hz trains that have been normalized to the initial amplitude and averaged across experiments. Reducing postsynaptic responsiveness with CNQX (Fig. 14A) or blocking postsynaptic receptor desensitization with CTZ (Fig. 14B) did not alter the accumulation of depression that occurred during the train. Manipulations that increase presynaptic transmitter release (Fig.14C,D) increased the rate and relative magnitude of the depression. Conversely, manipulations that decrease transmitter release (Fig. 14E,F) greatly decreased the depression and revealed facilitation during the early portion of the train. Similar results were obtained during application of adenosine (n = 9; data not shown).
To compare fit parameters across cells while minimizing the impact of the ambiguity caused by simultaneous facilitation and depression (see Fig. 8), we fit each set of responses with a model consisting of two components of depression without accompanying facilitation. Data from the various manipulations were pooled into three groups: manipulations of postsynaptic responsiveness (CNQX and CTZ), manipulations that increase presynaptic release (elevated external calcium and 4-AP), and manipulations that decrease presynaptic release (reduced external calcium, adenosine, and baclofen). For each group, we assessed the effect of the manipulations by comparing mean fit parameters obtained from data collected during drug application with those obtained before drug application. The more rapid form of depression (d1) was significantly weaker (i.e., closer to 1) during manipulations that decreased release probability (from 0.88 ± 0.017 to 0.97 ± 0.010;p = 5.7 × 10−5, pairedt test) and significantly stronger (i.e., farther from 1) during manipulations that increased release (from 0.86 ± 0.044 to 0.66 ± 0.084; p = 0.004). The value ofd1 was not significantly altered by manipulations of postsynaptic responsiveness (from 0.87 ± 0.013 to 0.88 ± 0.022; p = 0.64). The amplitude of the slower form of depression (d2) was also significantly decreased during decreased release (p = 0.016) but was not significantly increased during increased release (p = 0.74) or altered postsynaptic responsiveness (p = 0.89). The recovery time constants were not significantly altered during any of the treatments.
DISCUSSION
The present results show that synaptic responses in layer 2/3 evoked by temporally patterned stimulation of layer 4 exhibit a mixture of facilitation and depression with depression predominating under normal conditions. More importantly, they show that most of the apparently complex dynamics of these responses can be well described by a simple three-component model in which facilitation and two forms of depression are updated with each presynaptic stimulus and recover exponentially between stimuli. Describing synaptic dynamics in this way may provide some quantitative constraints that will be useful in delineating underlying mechanisms, but this is not the primary advantage of this approach. Its utility (and that of a similar approach developed independently by Tsodyks and Markram, 1997) is that it provides a concise and portable description that is useful for predicting synaptic responses to more complex patterns of stimulation, for computational studies of circuit dynamics, and for comparing dynamic properties across different synaptic pathways within or between preparations.
The present study focused primarily on understanding the average dynamics of populations of simultaneously activated synapses. As such, the information it provides is complimentary to that obtained from studies in which dynamic properties of individual synaptic connections have been examined (Thomson and Deuchars, 1994; Tsodyks and Markram, 1997). The similarity between the dynamic properties measured extracellularly from large populations of synapses and those measured by activating small numbers of inputs to individual neurons both in this study and in paired recordings of Tsodyks and Markram (1997)suggests that short-term plasticity at individual synapses is the dominant factor determining the dynamics of the population response. An important reason for measuring the dynamics of population responses is that it is large populations of synapses, rather than single connections, that are activated during sensory stimuli (although during sensory stimulation synapses may not be activated as synchronously as during extracellular stimulation). Even for individual cortical neurons, it is likely that sensory stimuli that evoke robust responses activate a substantial fraction of the many thousands of synaptic inputs of the cell. Measurements of the average dynamics of those inputs are therefore crucial to understanding how dynamic properties of synapses contribute to the dynamics of sensory responses.
The depression observed in the present study bears a number of similarities to forms of synaptic depression identified at other synapses. At the neuromuscular junction and squid giant synapse, the amplitudes of responses to successive presynaptic stimuli fall exponentially (Liley and North, 1952; Kusano and Landau, 1975). The mechanism underlying synaptic depression is believed to be depletion of a readily releasable pool of vesicles. The results of the pharmacological manipulations performed here are also consistent with a presynaptic mechanism for depression at cortical synapses. Decreasing the amount of release by lowering external calcium or by applying baclofen or adenosine reduced depression. Conversely, increasing release by increasing extracellular calcium or by applying 4-AP, a potassium channel blocker that augments calcium influx by decreasing action potential repolarization at synaptic terminals, increased depression. We cannot rule out the possibility that some of the manipulations presumed to alter presynaptic release may have also had postsynaptic effects. However, we also observed that reducing postsynaptic responsiveness with CNQX had no effect on the observed depression. In addition, depression was similar during field potential recordings of large synaptic responses and whole-cell voltage-clamp recordings of small synaptic responses. These data suggest that it is not the degree of postsynaptic depolarization that is correlated with the magnitude of depression but rather the amount of presynaptic release. An additional candidate postsynaptic mechanism, AMPA receptor desensitization, also seemed not to contribute significantly to depression over the range of intervals tested (>30 msec), because CTZ, although it had a clear effect on the duration of individual responses, did not alter the kinetics or amplitude of the observed depression. Similar findings have been reported for paired-pulse depression between neurons in hippocampal slice cultures (Debanne et al., 1996).
Another recent quantitative study of depression at neocortical synapses observed a single component of depression comparable to the faster component in this study but did not observe a longer time constant (Tsodyks and Markram, 1997). Several methodological differences could account for this difference. First, it is possible that in the present study, multiple populations of synapses with differing recovery rates are being sampled. Recent studies of cat visual cortex (Stratford et al., 1996) and rodent somatosensory and motor cortex (Thomson and Deuchars, 1994; Castro-Alamancos and Connors, 1996) indicate that different classes of synapses may exhibit different forms of short-term plasticity. Arguing against this interpretation, however, is the fact that the slower form of depression was also observed in EPSCs that, from their amplitudes, were unlikely to represent coactivation of more than a few synapses. As noted above, multiple time constants of recovery are common at single synapses in the case of facilitation and have occasionally been demonstrated in the case of depression (Kusano and Landau, 1975; also see Rosenmund and Stevens, 1996, their Fig. 7). Second, it is possible that the difference is attributable to the different synaptic pathways studied (layer 5 of somatosensory cortex vs layer 2/3 of visual cortex) or to differences in the ages of the animals (postnatal days 13–15 vs postnatal days 14–40). Finally, the longer-lasting component may be more easily recognized when fitting responses to longer duration (20–60 sec) trains. The slower form of depression observed in the present study has an amplitude and time constant similar to that observed for the replenishment of the readily releasable pool of transmitter at synapses between hippocampal neurons grown in culture (Liu and Tsien, 1995; Ryan and Smith, 1995; Rosenmund and Stevens, 1996).
Functional significance
We have suggested previously that depression at cortical synapses is a form of gain control that allows cortical neurons to respond to relative, rather than absolute, changes in firing rate (Abbott et al., 1997; for related ideas, see Grossberg, 1984; Tsodyks and Markram, 1997). Simulations of the responses of depressing synapses to fluctuations in presynaptic firing rate also reveal that such synapses are not performing linear temporal filtering. Instead, postsynaptic responses are enhanced for transients (Abbott et al., 1997; Tsodyks and Markram 1997; i.e., temporal “high-pass” behavior) yet also exhibit temporal “low-pass” behavior in that they respond poorly to high-frequency (>10–20 Hz) sinusoidal fluctuations in presynaptic firing rates (F. Chance, S. B. Nelson, and L. F. Abbott, unpublished observations). To the extent that synapses in the primary visual cortex of higher mammals are similar to those measured in rodent cortex, these properties of synaptic depression may explain some of the temporal nonlinearities noted previously in extracellular responses to visual stimuli. For example, the enhanced response to transients suggests that responses to flashed gratings should be larger than predicted on the basis of temporal frequency tuning for counterphased gratings, as has been observed (Tolhurst et al., 1980).
Perhaps the most dramatic temporal nonlinearity exhibited by cortical neurons is contrast adaptation. Contrast adaptation can occur locally within a portion of the receptive field of a cell (Marlin et al., 1991). It can also decrease the sensitivity of a neuron to particular stimuli without changing the maximal firing rate of the cell (Ohzawa et al., 1982). In single-unit and evoked-potential studies, contrast adaptation recovers with a time constant of 5–10 sec (Ho and Berkley, 1988; Maddess et al., 1988; Giaschi et al., 1993), although a more rapid form of adaptation recovering over the course of hundreds of milliseconds has also been described (Bonds, 1991; Nelson, 1991a). The magnitude and rate of onset of adaptation increase as the temporal frequency of the visual stimulus increases (Maddess et al., 1988;Bonds, 1991; Nelson, 1991a). These properties are consistent with a mechanism that depends on synaptic depression (Nelson, 1991b; Finlayson and Cynader, 1995; Somers et al., 1996; Nelson et al., 1997). Recent pharmacological (McLean and Palmer, 1996) and intracellular studies (Carandini and Ferster, 1997) in vivo support the hypothesis that contrast adaptation is caused by a reduction in excitatory synaptic drive, although the class of synapses affected and the biophysical mechanism involved remain undetermined. Our characterization of visual cortical short-term plasticity and identification of pharmacological means of manipulating that plasticity may permit a more-detailed analysis of its function in shaping the dynamics of cortical responses to visual stimuli.
Footnotes
This work was supported by National Science Foundation Grants 9511094 and 9421388, the Sloan Foundation, National Institutes of Health Grant EY11115, and the W. M. Keck Foundation. We thank Gina Turrigiano for helpful discussions.
Correspondence should be addressed to Dr. Sacha Nelson, Department of Biology MS 008, Brandeis University, Waltham, MA 02254.
Dr. Gibson’s present address: Neuroscience Department, Box 1953, Brown University, Providence, RI 02912.