Abstract
Nearby cortical neurons often have correlated trial-to-trial response variability, and a significant fraction of their spikes occur synchronously. These two forms of correlation are both believed to arise from common synaptic input, but the origin of this input is unclear. We investigated the source of correlated responsivity by recording from pairs of single neurons in primary visual cortex of anesthetized macaque monkeys and comparing correlated variability and synchrony for spontaneous activity and activity evoked by stimuli of different orientations and contrasts. These two stimulus manipulations would be expected to have different effects on the cortical pool providing input to the recorded pair: changing stimulus orientation should recruit different populations of cells, whereas changing stimulus contrast affects primarily the relative strength of sensory drive and ongoing cortical activity. Consistent with this predicted difference, we found that correlation was affected by these stimulus manipulations in different ways. Synchrony was significantly stronger for orientations that drove both neurons well than for those that did not, but correlation on longer time scales was orientation independent. Reducing stimulus contrast resulted in a decrease in the temporal precision of synchronous firing and an enhancement of correlated response variability on longer time scales. Our results thus suggest that correlated responsivity arises from mechanisms operating at two distinct timescales: one that is orientation tuned and that determines the strength of temporally precise synchrony, and a second that is contrast sensitive, of low temporal frequency, and present in ongoing cortical activity.
- noise correlation
- visual cortex
- orientation
- contrast
- spontaneous activity
- cross-correlogram
- signal correlation
Introduction
When a stimulus is shown repeatedly, trial-to-trial fluctuations in response strength are correlated between neighboring cortical neurons. This correlation, typically characterized by the spike count correlation rsc, interferes with pooling as a strategy to overcome response variability in a population rate code (Zohary et al., 1994; Shadlen and Newsome, 1998). The timing of action potentials between nearby neurons is also often correlated or synchronized, as shown by a peak in the spike train cross-correlogram (CCG). Temporal coding theories propose that synchronous spikes bind the distributed representation of objects (Gray, 1999; Shadlen and Movshon, 1999). Although synchrony and rsc traditionally have been linked to differing views of how information is encoded by neuronal populations, there is a growing appreciation that they are functionally and mathematically related (Bair et al., 2001; Mazurek et al., 2002).
Both rsc and synchrony arise from shared input, the source of which is unclear. In the visual system, rsc is either assumed (Gawne and Richmond, 1993; Gawne et al., 1996; Reich et al., 2001) or known to be stimulus independent (Zohary et al., 1994; Bair et al., 2001), suggesting that it arises from variations in the state of the cortical network unrelated to the input provided. Synchrony, on the other hand, often is reported to depend on whether a single stimulus is presented or its constituent parts (Gray, 1999). Indirect measures of synchrony also have been shown to depend on basic stimulus parameters such as orientation and speed (Gray and Viana Di Prisco, 1997; Friedman-Hill et al., 2000; Frien et al., 2000). It is unclear, however, how the synchrony between single neurons in primary visual cortex (V1) is affected by such basic stimulus manipulations. Because synchrony and rsc are related, it is also unclear how the former can be stimulus dependent while the latter is not.
To resolve this discrepancy and to determine how each form of correlation arises, we have measured both rsc and synchrony between single V1 neurons during spontaneous activity and for activity evoked by stimuli of different orientations and contrasts. Changing stimulus orientation will likely alter the pool of cortical neurons that provides input to the recorded pair, with the predicted effect on correlation depending on the architecture of V1: if the dominant input to V1 cells is feedforward, correlation should be orientation independent; if recurrent connections predominate, correlation should be sensitive to orientation (Ben-Yishai et al., 1995; Series et al., 2004). Changing stimulus contrast, on the other hand, primarily alters the magnitude of the response of the pool, thus allowing us to compare the relative influence of sensory drive and ongoing cortical activity on correlation.
Characterizing correlation and its stimulus dependence is also important for understanding the computations that can be performed by a neuronal population. For instance, understanding the stimulus dependence of synchrony is necessary for evaluating its role in binding and for determining how it affects the maximal stimulus frequency that can be encoded in a population firing rate (Mazurek and Shadlen, 2002). In addition, theoretical studies suggest that capacity of a population rate code to encode information depends on the magnitude of correlated variability, its relationship to the tuning similarity of two neurons, and its dependence on stimulus drive (Abbott and Dayan, 1999; Sompolinsky et al., 2001; Shamir and Sompolinsky, 2004). Finally, measuring correlation provides a way to test current models for reading out noisy population codes (Pouget et al., 1998), because these make testable predictions about the stimulus dependence of correlation.
We find that stimulus orientation and contrast affect rsc and CCGs in distinct ways. High-contrast stimuli result in sharper synchronization; low-contrast stimuli have stronger overall correlation but on a broad time scale. This suggests that orientation controls the amount of shared stimulus drive that nearby cells receive, whereas contrast may control internal cortical dynamics.
Materials and Methods
Recordings were made in 10 cynomolgus (Macaca fascicularis), one bonnet (Macaca radiata), and one pig-tailed (Macaca nemestrina) adult male monkeys. All experimental procedures were approved by the New York University Animal Welfare Committee.
The procedures used in our laboratory for single-unit recording in anesthetized, paralyzed macaque monkeys have been described previously in detail (Cavanaugh et al., 2002). Briefly, animals were premedicated with atropine (0.05 mg/kg) and diazepam (1.5 mg/kg) and anesthetized initially with ketamine HCl (10 mg/kg). Anesthesia during recording was maintained by intravenous infusion of the opiate anesthetic sufentanil citrate (Sufenta; 4-8 μg · kg-1 · h-1). To minimize eye movements, the animal was paralyzed with intravenous infusion of vecuronium bromide (Norcuron; 0.1 mg · kg-1 · h-1). Vital signs (EEG, electrocardiogram, end-tidal PCO2, temperature, and lung pressure) were monitored continuously. The pupils were dilated with topical atropine and the corneas protected with gas-permeable contact lenses. Refraction was provided by supplementary lenses.
Recordings were made with an array of seven concentrically arranged platinum/tungsten microelectrodes (Thomas Recording, Giessen, Germany). We made a 10-mm-diameter craniotomy, centered 10 mm lateral to the midline and 5-10 mm posterior to the lunate sulcus, and made vertical electrode penetrations or penetrations normal to the cortical surface. Signals from the microelectrode were amplified, bandpass filtered (300 Hz-10 kHz), and fed into a hardware discriminator (BAK Electronics, Germantown, MD) and audio monitor. Spike times were saved with a temporal resolution of 0.25 ms. Most pairs of cells (90.5%) were recorded on separate electrodes; pairs were recorded on one electrode only when the waveforms of the cells were clearly distinct.
To allow histological confirmation of the recording site, we made electrolytic lesions at the end of each recording track by passing DC current (typically 2 μA for 5 s) through the recording electrodes. At the end of the experiment, the monkeys were killed with an overdose of sodium pentobarbital and perfused through the heart with 0.1 m PBS followed by 4% paraformaldehyde in 0.1 m PBS. Sections (40 μm) were stained for Nissl substance with cresyl violet. The close spacing of the electrodes limited our ability to assign lesions to particular electrode tracks. Instead, we used the scatter is lesion location normal to the tracks, together with depth readings from the microdrive system, to calculate the approximate distance between recording sites. Most cells were separated by 300-500 μm (maximum tangential separation of 600 μm and radial separation of 400 μm).
Visual stimuli. Stimuli were luminance modulated, drifting sine-wave gratings presented at a frame rate of 100 Hz using an eight-bit Silicon Graphics (Mountain View, CA) board operating at a resolution of 1024 × 731 pixels. The monitor (Eizo T550) subtended ∼10° of visual angle and had a mean luminance of ∼33 cd/m2. For most cells, we measured, in order, the direction, spatial and temporal frequency, and size tuning for drifting sine-wave gratings (133 of 147 pairs; optimal stimulus parameters in the remaining pairs were determined by hand). After characterizing the stimulus preference of each cell independently, we measured correlation for spontaneous activity and for responses evoked by stimuli of different orientations and contrasts. The spatial and temporal frequency of these test stimuli were set between the preferred values of the two cells or at the value of the cell that responded less vigorously. Stimuli were presented between the receptive field (RF) centers of the cells and covered both RFs. All measurements were made using stimuli presented in a circular aperture to the dominant eye of the less responsive cell. Stimuli were surrounded by a gray field of average luminance.
In the orientation experiments, we presented full-contrast gratings at five orientations that spanned a range between driving the two cells strongly and evoking a weak response; the range spanned 25-122.5°, depending on the tuning bandwidth of the neurons. In the contrast experiments, we used the orientation that was most effective at driving the pair and presented it at four contrasts, typically 1.56, 6.25, 25, and 100%, with contrast defined as the difference between the maximum and minimum luminance divided by the sum. For cells with low contrast sensitivity, we adjusted the range of contrasts to extend from 12.5 to 100% in octave steps. Stimuli were presented for 2.56 s and separated by 3 s intervals during which we presented a gray, isoluminant screen. Each stimulus was presented 30-200 times (mean, 80 ± 3 trials). Contrast and orientation data were collected in separate blocks of trials. In most experiments, the presentation order of each orientation or contrast was block randomized, but in some experiments (∼10% of the data), we presented each orientation or contrast repeatedly before testing the next stimulus. We observed no difference between the data from the randomized and nonrandomized stimulus presentation so the data were pooled.
Spike count and timing correlation. We characterized correlation using two metrics, the spike count correlation rsc and the spike train cross-correlogram (Perkel et al., 1967), calculated using the methods of Bair et al. (2001).
The rsc, the Pearson correlation coefficient of the evoked spike counts of two cells to the repeated presentation of a particular stimulus, was calculated as follows: (1)
where E is the expected value, σ is the SD of the responses, and N1 and N2 are the spike counts of cell 1 and 2, respectively. To avoid contamination of our estimates of rsc by outlier responses, we removed trials on which the response of either neuron was >3 SDs different from its mean response (Zohary et al., 1994). Statistical evaluation was performed after converting rsc to Z-scores using the Fisher transformation (Anderson, 1984) as follows: (2)
To compute the spike train CCG, we represented the spike train of each cell as a binary time series with 1 ms resolution such that: (3)
The CCG was then computed as follows: (4)
where M is the number of trials, N is the number of bins in the trial, xi1 and xi2 are the spike trains of neurons 1 and 2 on trial i, τ is the time lag, and λ1 and λ2 are the mean firing rates of the two cells. Θ(τ) is the following triangular function: (5)
where T is the trial duration in seconds. This function corrects for the degree of overlap of the two spike trains for each time lag (i.e., that there are T opportunities for simultaneous events in a trial of length T but only T - 1 opportunities for coincidences at time lags of 1 ms, etc.).
We chose to normalize the CCG by the geometric mean spike rate, as opposed to the product of spike rates or some other factor, for two reasons. First, it is the most commonly used normalization, so it facilitates comparison with previous studies (Mastronarde, 1983; Eggermont, 1994; Eggermont and Smith, 1996; Bair et al., 2001) [similar to the normalization of the joint peristimulus time histogram (Aertsen et al., 1989)]. Second, this normalization provides results that are most comparable with measures of spike count correlation. We explore the reason for this similarity, and how our conclusions about synchrony depend on the normalization we use, in the appendix [available at www.jneurosci.org as supplemental material (references therein: Tolhurst et al., 1983; Vogels et al., 1989)].
All CCGs were corrected for correlation induced by the stimulus by subtracting a shift predictor calculated from trials 1 to n - 1, with an offset of one trial. Specifically, the shift predictor was calculated as in Equation 4, but with the numerator as the following: (6)
We chose to use a shift predictor rather than an all-way shuffle correction, because the latter is affected more strongly by slow fluctuations in neuronal responsivity. In cells in which the firing rate displayed strong temporal modulation (i.e., simple cells), shuffle correction occasionally resulted in artifactual peaks in the CCG caused by bleed through of this modulation. These peaks were clearly distinct from those we report here: they had widths related to the stimulus modulation frequency (several hundred milliseconds wide), and additional peaks were evident at a range of time lags. Shift correction effectively reduced the presence of these relatively rare artifactual peaks but otherwise resulted in CCGs that were essentially identical to those obtained with shuffle correction.
We smoothed the shift-corrected CCGs with a 5 ms kernel ([0.05 0.25 0.40 0.25 0.05]) before characterizing CCG shape. We computed CCG peak heights and widths for peaks that exceeded 5 SDs of the noise level at time lags from -1000 to -950 ms and from 950 to 1000 ms. Our results were not sensitive to the value of this criterion; we chose a conservative value of Z = 5 so that only CCGs with clearly visible peaks would be considered. In addition, we present data both above and below this significance criterion. Peak height was defined as the difference between the maximal bin size within 10 ms of zero time lag and the mean CCG value at large time lags (τ =±950-1000 ms), which was usually close to zero. Peak widths were calculated as the difference between the first positive and negative time lag at which the CCG peak reached less than half maximum height (full width at half maximum). For pairs recorded from one electrode, we ignored the -1 to +1 ms bins when calculating CCG peak widths because of our inability to detect simultaneous or nearsimultaneous spiking.
We analyzed every neuronal pair for which we recorded at least 30 repetitions of each stimulus (five pairs for the orientation series and one pair for the contrast series failed this criterion), because our estimates of rsc would otherwise be unreliable, and the number of spikes recorded would be insufficient for calculating CCGs. We also removed five pairs in the orientation data set and three pairs in the contrast series, because one of the cells in each pair had an abnormally unreliable response, as indicated by a variance-to-mean response ratio >20.
To calculate the time scale of correlation, we used the following metric introduced by Bair et al. (2001): (7)
where t is the window of integration. In this case, the CCG in the numerator was calculated as the raw CCG (the numerator of Eq. 4) and was corrected by the all-way shuffle predictor. Autocorrelogram 1 (ACG1) and ACG2, the autocorrelograms of the two neurons, were calculated similarly, but the spike train on each trial was compared with itself. That is, the raw ACG for cell 1 was calculated as follows: (8)
and this was corrected by an all-way shuffle predictor. ACG2 was calculated similarly.
Signal correlation analysis. Our measurements of the tuning of the cells for stimulus orientation, spatial and temporal frequency, and size allowed us to explore the relationship between the similarity in receptive field properties of two cells and the magnitude of rsc or CCG peak height. We quantified this relationship in two ways. First, we calculated the tuning similarity as the signal correlation rsignal, which is the Pearson correlation (Eq. 1) between the mean response of each cell to each stimulus. The relationship between rsignal and rsc or CCG peak height was then characterized by the correlation between these values across all cell pairs. In the second method, we characterized the similarity in tuning as the difference in the preferred orientation or direction, the preferred spatial and temporal frequency, and the degree of RF overlap. These differences were then used as independent variables in a multivariate regression model to predict rsc or CCG peak height.
The preferred direction for each cell was defined by the vector sum of the responses to 16 test directions (22.5° steps) as follows: (9)
where θi and Ri are the direction of motion and response strength, respectively, for stimuli i = 1-16. We characterized selectivity as the following: (10)
which is the length of the summed response vector, normalized by the summed magnitude of all response vectors. The selectivity index varies between 0, for a cell that responds equally to all directions, and 1, for a cell that responds to a single direction of motion. Preferred orientation and orientation selectivity were calculated in a similar manner, except that θn was replaced by 2θn in the equations above and the resultant angle halved. We used the preferred direction for cells for which the direction selectivity was greater than the orientation selectivity and the preferred orientation otherwise.
The preferred temporal and spatial frequency were calculated by fitting a difference-of-Gaussians function to the data, using the STEPIT algorithm (Chandler 1969) to minimize the combined χ2 error between the model predictions and the data. From the fitted equation, we defined the optimal stimulus parameter as that which would evoke the strongest predicted response. The difference between the preferred spatial and temporal frequencies of the two cells was then defined in octaves as follows: (11)
To measure the receptive field overlap of the two cells, we determined the receptive field center by hand, using a small (∼0.3°) patch of optimal grating. We then measured responses to gratings of increasing size and fit the data with a difference-of-error function, using the STEPIT algorithm. The receptive field size of each cell was defined as the maximum of the function in the range tested, and the overlap was defined as the percentage of the smaller RF that was included in the larger RF. The mean receptive field size provided by this approach is approximately twice that provided by hand maps using small bars of light (Cavanaugh et al., 2002). As a result, our estimate of RF overlap is substantially higher than that which would result from mapping with small stimuli.
In our regression analysis, we only used data from pairs for which the fits for both cells accounted for at least 50% of the variance (123 of 133 pairs for the spatial and temporal frequency data; 114 of 133 pairs for the area data). The variance accounted for by the fits in these cells was on average 90-92% for each parameter.
Results
We recorded from 147 pairs of single units in 12 anesthetized, paralyzed macaque monkeys. Because our primary objective was to measure the stimulus dependence of neuronal correlation, we recorded from nearby neurons (typically <500 μm apart) that had similar receptive field properties, because distant or dissimilar neurons tend to fire independently (Nelson et al., 1992; Lee et al., 1998; DeAngelis et al., 1999; Nowak et al., 1999; Bair et al., 2001) and have weak correlation in response variability (Zohary et al., 1994, Lee et al., 1998; Bair et al., 2001; Averbeck and Lee, 2003). We recorded in all cortical layers but biased our population toward complex neurons (76% of the population) (Skottun et al., 1991). The receptive field properties of the neurons comprising each pair were similar, with a mean difference of 37° in orientation preference, 0.37 octaves in spatial frequency preference, 0.36 octaves in temporal frequency preference, and a mean receptive field overlap of 75%. The ocular dominance of the two cells was also similar, with a mean difference of 0.83 on the seven-point scale of Hubel and Wiesel (1962).
Orientation dependence of spike count correlation
We evaluated the orientation dependence of rsc, the correlation of evoked spike counts (Eq. 1 in Materials and Methods), by measuring responses to 2.56 s presentations of full-contrast gratings of five orientations. Figure 1A shows the orientation tuning and range of orientations (thick line) used to measure correlation for an example pair. We chose orientations that spanned a range from driving the pair strongly [geometric mean response of 33 impulses per second (ips)] to evoking a relatively weak response (8.6 ips). Scatter plots of the response of the two cells to multiple presentations of each stimulus are shown in Figure 1B-E as Z-scores relative to the mean response for each stimulus. The value of rsc (text in scatter plots) varied among stimulus conditions but did not depend in an obvious way on stimulus orientation or the evoked firing rate. For instance, the correlation for the stimulus that drove both cells strongly (0.24) (Fig. 1D) was similar to that for a stimulus that drove one cell but not the other (0.30) (Fig. 1F).
The data presented in Figure 2 show frequency histograms for rsc in our population of pairs (n = 100), arranged for each pair from the orientation that was most effective at driving the two cells to that which was least effective. We found little relationship between the efficacy of the stimulus and the magnitude of spike count correlation (ANOVA; p = 0.45). Stimuli that drove the pair most strongly (42 ± 2 ips) had an average correlation of 0.18 ± 0.03 (Fig. 2A), similar to the average rsc value of 0.19 ± 0.02 for stimuli that evoked the weakest response (12 ± 1 ips) (Fig. 2E). The mean rsc collapsing across all conditions and pairs was 0.20, a value consistent with previous measurements in the visual system, including those in V1 [0.22 in Gawne et al. (1996) and ∼0.25 in Reich et al. (2001)], middle temporal visual area (MT) [0.19 in Zohary et al. (1994) and 0.20 in Bair et al. (2001)], and inferior temporal cortex [0.23 in Gawne and Richmond (1993)]. Because strong trends between stimulus efficacy and rsc in individual pairs may go undetected in a population analysis, we also calculated the relationship between the evoked firing rate and rsc for each pair individually. We found a significant correlation (p < 0.05) in only 7 of 100 pairs, three of which were positively correlated and four of which were negatively correlated.
We conclude that there is little relationship between the efficacy of an oriented stimulus and the correlation in trial-to-trial variability of evoked spike count, suggesting that this variability arises from orientation-independent variations in trial-to-trial cortical excitability.
Orientation dependence of spike timing correlation
Whereas the orientation independence of rsc agrees well with previous studies (Zohary et al., 1994; Bair et al., 2001), a related form of correlation, spike timing synchrony, has been shown to depend on stimulus drive. However, most studies investigating synchrony between single V1 cortical neurons have focused either on the effect of altering the “gestalt” characteristics of the stimulus (Livingstone, 1996) or have used indirect measurements such as the synchrony of multiunit activity (MUA) (Lamme and Spekreijse, 1998) or the strength of oscillations in single-unit activity, MUA, or the local field potential (LFP) (Gray et al., 1989; Gray and Viana Di Prisco, 1997; Friedman-Hill et al., 2000; Frien et al., 2000). The relationship between the synchronous firing of single neurons and these measurements is unclear. For MUA synchrony, this is because different stimuli might recruit different subpopulations of neurons in the recorded signal. For measures of oscillatory firing, this is because the stimulus dependence of oscillations is controversial (Ghose and Freeman, 1992; Young et al., 1992) and because some studies have failed to find a strong relationship between oscillatory firing and synchronous firing of single neurons (Samonds and Bonds, 2004) (but see Maldonado et al., 2000). To compare the sensitivity of spike count and spike timing correlation between single V1 neurons, we therefore investigated the orientation dependence of synchrony in our population of pairs.
Figure 3, A and B, shows CCGs for two pairs of neurons (left and center column), arranged in each case from the most effective to least effective orientation. The CCGs for both pairs had a clear peak centered on zero time lag (inset in top row of traces), as did all but two of the pairs we recorded, indicating that the spike timing correlation was attributable to common synaptic input rather than a direct connection between cells (Moore et al., 1970). For both pairs, the CCG was clearly dependent on stimulus orientation: the peak of the CCG was smaller for stimuli that evoked a weaker mean firing rate and disappeared entirely for orientations that evoked little response. Figure 3C shows the average CCG for the population of V1 pairs with significant peaks for at least one stimulus condition (n = 56 of 100 pairs; see Materials and Methods for criterion), calculated after normalizing data from each pair by the CCG peak height of the most effective orientation. The peak height of the average CCG was smaller for stimuli that evoked weaker responses, but CCG shape was similar across stimulus conditions. Note that for each pair and stimulus condition, the CCG was normalized by the geometric mean firing rate, so that the reduced peak height is not a trivial consequence of the reduced firing rate. We investigate the effect of firing rate on CCG peak height further below, when we consider the effect of stimulus contrast. Other forms of normalization for spike rate are considered in the appendix [available at www.jneurosci.org as supplemental material (references therein: Tolhurst et al., 1983; Vogels et al., 1989)].
We quantified the orientation dependence of spike timing correlation by measuring CCG peak height (maximum value within 10 ms of 0 ms time lag) and full width at half height (in milliseconds). The scatter plot in Figure 4A shows the change in CCG peak height for all pairs, comparing the stimulus that evoked the highest firing rate with that which evoked the weakest response. Almost all of the points lie below the identity line, indicating that the peak height was consistently higher for the best-compromise stimulus; the mean peak height decreased from 0.015 to 0.007 coincidences per spike (coinc/spk) (57%). CCG peak width, on the other hand, was essentially unaffected by changes in stimulus orientation. To ensure meaningful measurements, we computed CCG peak widths only for those CCGs with a Z-score of 3 or higher (see Materials and Methods; 55 of 56 pairs had significant CCG peaks for at least two stimulus conditions). The mean peak width was 9.0 ms for the best-compromise stimulus and 9.7 ms for stimuli evoking the lowest mean rate (p = 0.63) (Fig. 4B).
We conclude that spike timing correlation is sensitive to stimulus orientation: for orientations ineffective at driving the pair, the CCG peak height is substantially smaller but of similar shape to that for activity evoked by the orientation that evokes the strongest mean response in the pair.
Time scale of correlation
The dependence of the CCG on stimulus orientation stands in stark contrast to the orientation independence of spike count correlation. This is surprising because both are believed to arise from common synaptic input (Shadlen and Newsome, 1998). Stronger still, the CCG and rsc are mathematically related (Eq. 7 in Materials and Methods): rsc is equal to the integral of the CCG, divided by the geometric mean area of the respective autocorrelograms (Bair et al., 2001). Integrating the CCG and ACGs for time lags equal to the trial duration provides spike count correlation values identical to rsc calculated directly from the total number of spikes evoked during the trial (Figs. 1, 2). Integrating over shorter windows can reveal the time scale over which spike count correlation arises and can clarify why these two forms of correlation have a different sensitivity to stimulus orientation. Importantly, the time scale over which V1 responses are integrated in natural viewing conditions is likely to be on the order of 10-100 ms, because stimulus drive changes frequently (e.g., because of frequent eye movements) and because the neural computations underlying perception occur quickly (within 150-300 ms) (Thorpe et al., 1996; Johnson and Olshausen, 2003). For pairs with significant spike timing correlation, the orientation dependence of CCG peak height would suggest that rsc on this more relevant time scale would be largest for the stimulus that evokes the highest mean rate in the pair.
We evaluated the time scale of correlation by integrating the CCG and ACGs at a time resolution of 1 ms, and refer to rsc calculated in this manner as rCCG. We averaged the rCCG curves for each stimulus condition across our population of pairs and scaled the final value to one to facilitate comparison of the time course. The value of rCCG, as a function of the size of the integration window, for the most effective stimulus is shown as a dotted line in Figure 5A. The value of rCCG increased with the size of the integration window and reached its final value for windows of 50-100 ms, a value approximately consistent with previous measurements in MT (32 ms) (Bair et al., 2001) but shorter than the values reported for cat V1 (500-1000 ms) (Reich et al., 2001). For less effective orientations (dashed line for the stimulus that evoked the next highest firing rate, followed by solid lines of increasing thickness), the respective rCCG curves had lower values for small integration windows (e.g., <10 ms), reflecting the loss of area in the central peak of the CCG. We calculated the time at which each average rCCG curve reached 90% of its final value, indicated by the dots on the curves in Figure 5A. As the stimulus orientation became more effective at driving the pair, the 90% rise time decreased from 231 ms for the least effective stimulus to 53 ms for the most effective stimulus.
A consequence of the loss in CCG peak height for stimuli evoking weak mean firing in the pair is that rCCG for brief response epochs will be orientation dependent. To illustrate this, we replotted the relationship between rCCG and the time scale of integration for the most (gray trace) and the least (black trace) effective stimulus in Figure 5B. These curves are not normalized; therefore, rCCG can be read directly off the ordinate. For an integration window of 1 ms, rCCG for the best-compromise stimulus was approximately three times larger than for the least effective stimulus. The difference between the rCCG curves is maintained for integration windows less than ∼200 ms; over this range, the rCCG curve for the least effective stimulus resembles the rCCG curve for the best-compromise stimulus but with a negative offset as a result of the loss of area on short time scales. Consistent with the orientation independence of rsc, the rCCG curves converge for integration windows of >200 ms, and the value for the largest integration window (2600 ms) was not significantly different for the two conditions. These data suggest that orientation-sensitive common input determines the strength of correlation on time scales relevant for cortical computation but that an orientation-independent mechanism dominates on the time scale of hundreds of milliseconds to seconds.
Magnitude and time scale of correlation for unsynchronized pairs
The analysis in the preceding sections focused on the stimulus dependence of spike timing correlation and its effect on rCCG in pairs with significant CCG peaks (56 of 100 pairs). In the remaining pairs, synchronization was weak or absent for all stimulus conditions, as illustrated by the average CCGs for the most effective (33 ± 3 ips; gray trace) and least effective (9 ± 1 ips; black trace) stimulus conditions shown in Figure 6A. Because of the relationship between rsc and the CCG, pairs without significant spike timing correlation should have weak rCCG values at all time scales, and rCCG should be stimulus independent. The rCCG plots shown in Figure 6B show that this was indeed the case. First, there was no difference between the rCCG curves for effective (thick gray line) and ineffective (thick black line) orientations for pairs without CCG peaks, and the correlation for both conditions was nearly zero at short time scales (rCCG < 0.03 for integration windows <30 ms). Second, the correlation for the best-compromise orientation in pairs with significant CCG peaks (thin black line; replotted from Fig. 5B) was substantially larger than for pairs without CCG peaks (thick gray and black traces) at all time scales. For the largest integration window, the rCCG values for pairs with significant CCG peaks was 0.22 ± 0.03 compared with 0.11 ± 0.02 for pairs without (p = 0.007).
In summary, pairs with significant CCG peaks, presumably those that share substantial common input, display a significant orientation dependence of spike timing correlation. As a result, the response variability of neurons is significantly more correlated for effective than ineffective orientations for response epochs of <100 ms. At long time scales, there is little stimulus dependence evident. In pairs without spike timing synchronization, presumably those that share less common input, the correlation of trial-to-trial response variability is significantly lower, and no stimulus dependence is evident. In a later section, we will show that the strength of synchronization and rsc are related to the similarity in receptive field properties of the two neurons.
Effect of stimulus contrast on correlation
Whereas synchrony and correlated response variability on short time scales are stronger for orientations that evoke a robust response in both neurons, correlation is primarily orientation independent on longer time scales and thus determined by some other source. To determine whether correlation on long time scales is at all stimulus sensitive or simply reflects ongoing cortical activity, we measured the effect of reducing stimulus contrast on correlation, a manipulation that would be expected to alter the strength of the evoked cortical response rather than to alter strongly which neurons are recruited by the stimulus. Comparing the effect of stimulus orientation and contrast on correlation is also important for technical reasons: the measured strength of synchrony is notoriously sensitive to firing rate, even when the strength of common input is fixed (Melssen and Epping, 1987; Binder and Powers, 2001; Turker and Powers, 2002; Shu et al., 2003). This makes it difficult to know whether effects on synchrony reflect orientation-dependent changes in the amount of common input or simply a change in our ability to detect this synchrony because of the rate dependence of the CCG. By comparing the effects of orientation and contrast on synchrony, we can distinguish between these possibilities.
We begin by computing rsc for the best-compromise grating (i.e., the orientation that evokes the highest mean firing rate in the pair) presented at a range of contrasts, typically 1.25 to 100% in double-octave steps. The dependence of rsc calculated from the entire response epoch on stimulus contrast for our population of V1 pairs (n = 57 pairs) is shown in Figure 7, with the data arranged from the strongest to the weakest stimulus (almost always the highest to lowest contrast). The mean rsc value was significantly larger for intermediate- and low-contrast stimuli than high-contrast stimuli (ANOVA; p < 0.001). The dependence of rsc on contrast was not linear; rather, rsc was similarly high over a wide range of low to moderate contrasts but dropped strongly at maximum contrast. That the correlation in trial-to-trial response variability was weaker for full-contrast stimuli suggests that strong stimulus drive can reduce the correlation caused by ongoing cortical activity (for similar findings in area MT, see Cardoso de Oliveira et al., 1997; Bair et al., 2001).
The elevated rsc values at low contrast suggest that the area under the CCG increases as stimulus contrast is lowered, as opposed to the sharp loss in peak height observed as the orientation of a high-contrast stimulus is changed. The two example pairs in Figure 8, A and B (left and center column), indicate that this was indeed the case. At low stimulus contrast (bottom rows), CCG peak height was reduced in amplitude, but there was a compensatory increase in peak width. The population average CCG (n = 57) (Fig. 8C) reveals both a substantial increase in peak width at lower contrasts and a loss in the amplitude of the narrow central peak. On a pair-by-pair basis, there was a significant drop in CCG peak amplitude: the mean amplitude decreased by 46% as stimulus efficacy decreased. Distinct from the effect of altering stimulus orientation, however, there was a substantial increase in the CCG peak width for low-contrast stimuli, with the mean width at half maximum more than doubling from 11.5 to 24.0 ms (p = 0.006).
In principle, the increase in rsc and the appearance of a broad CCG at low contrast could arise because of adaptation. For instance, the response of both cells to a low-contrast stimulus could be suppressed by a preceding high-contrast stimulus as a result of rapid changes in contrast gain (Bonds, 1991). The adaptation-induced fluctuations in response strength would inflate rsc and could lead to artifactual CCG peaks (Brody, 1999). To control for the potential effects of adaptation, we compared the correlation for a low-contrast stimulus presented in isolation with that measured when the same stimulus was interleaved with high-contrast stimuli (n = 12 pairs). The rsc values were nearly identical for the two cases, with a mean value of 0.52 and 0.53 for isolated and interleaved stimuli, respectively, and the CCGs were also similar (width of 41 vs 33 ms; height of 0.016 vs 0.015 coinc/spk). In an additional seven pairs of cells, we measured responses to noninterleaved high- and low-contrast stimuli and found that consistently larger rsc values (0.28 for high contrast vs 0.47 for low contrast) and broader CCGs (9.2 vs 24.8 ms) for the low-contrast stimuli. Finally, the correlation we observe at low contrast agrees well with that of spontaneous activity, as discussed below. Together, these results suggest strongly that adaptation does not cause the elevated correlation at low contrast.
The increase in CCG peak width at low contrast suggests a significant increase in the time scale of correlation for these stimuli. Figure 9A shows the population-normalized rCCG curves for high-contrast stimuli (dotted line) and for stimuli of decreasing contrast (dashed line and lines of increasing thickness). The time scale of correlation decreased monotonically as the stimuli became more effective: the time at which rCCG reached 90% of its final value decreased from 499 ms (thick solid line) for the lowest contrast to 70 ms (dotted line) for full contrast. In Figure 9B, unnormalized rCCG plots for high- and low-contrast stimuli allow direct comparison of values on different time scales. The correlation on very brief time scales (<10 ms) was higher for full-contrast stimuli (gray trace) because of the relatively large sharp peak in the CCG for this condition. For intermediate time scales (10-100 ms), the increased CCG peak width at low contrast compensated for the loss of area in the narrow central peak, and the rCCG values are essentially identical. Finally, for integration windows greater than ∼100 ms, rCCG was significantly larger for the low-contrast data, consistent with the results shown in Figure 7. Thus, although the 90% rise times increase for both low-contrast and ineffective orientations, they do so for different reasons. For changes in stimulus orientation, the effect is attributable primarily to a subtractive shift in the rCCG curves caused by reduced sharp synchrony. For contrast, the change in rise time is caused by an increase in the value of rCCG for large integration windows.
Clearly, stimulus contrast and orientation affect correlation differently. However, the preceding analysis involved overlapping but not identical sets of neuronal pairs, and, more importantly, the range of firing rates explored in the contrast experiments was larger than that for orientation (primarily because many V1 neurons respond vigorously to nonpreferred orientations) (Ringach et al., 2002). To make a more direct comparison, we evaluated the stimulus dependence of correlation for 37 pairs of neurons for which we obtained both types of data and for which the CCG peak was significant. We compared the correlation observed at the least effective orientation (mean response of 15 ips) with that of the contrast that evoked the most similar firing rate (mean response of 16 ips). The average CCG for the low-contrast stimulus (Fig. 10A, gray trace) was larger and broader than the average CCG evoked by a high-contrast stimulus, the orientation of which evoked a similar firing rate (black trace). The scatter plots in Figure 10, B and C, show the difference in CCG peak height and width, respectively, on a pair-by-pair basis for the two conditions. Peak height was significantly larger (0.011 vs 0.007 coinc/spk; p = 0.003), and the width at half maximum was substantially broader (27.3 vs 10.1; p = 0.024) for the low-contrast stimulus. The better maintained CCG peak height and increased peak width resulted in significantly higher rsc values for the low-contrast stimulus (0.48 vs 0.28; p < 0.001) (Fig. 10D), a difference that was apparent at all time scales (Fig. 10E, gray rCCG curve for low contrast vs black for ineffective orientation).
In summary, our results suggest that correlated response variability and spike timing correlation are affected differently by changes in stimulus orientation and contrast. Stimulus orientation controls the source of input that causes sharp synchrony, but correlation on longer time scales is orientation independent. This latter source of correlation is, however, affected by stimulus drive: high-contrast stimuli suppress it.
Correlation of spontaneous activity
The strong correlation we observe for activity evoked by low-contrast stimuli suggests that cortical neurons receive substantial common input in the absence of effective stimulus drive. This proposal is consistent with recent optical imaging studies of the spatial distribution of spontaneous activity in V1 (Arieli et al., 1996; Tsodyks et al., 1999), which found that population activity preceding spontaneous spikes was similar to that evoked by an optimal stimulus. Thus, it appears that spontaneous spikes occur because the neurons that provide input to the recorded cell are active, a description that applies equally well to evoked activity. More recently, Fizer et al. (2004) reported that the spatial and temporal structure of ongoing cortical activity was relatively insensitive to stimulus drive. To provide a direct comparison of the correlation of evoked and spontaneous activity, we recorded spontaneous activity for periods of 15-30 min in neuronal pairs that had an appreciable spontaneous firing rate (mean rate of 7.2 ips; n = 42 pairs).
Consistent with the imaging study of Tsodyks et al. (1999), we typically found a significant peak in both the spontaneous and stimulus-evoked CCG (48% of pairs) or in neither (26%). In 10% of pairs, we observed a significant peak in the spontaneous CCG but not the evoked data, and in 17% we observed the opposite, although in many of these cases, the CCGs just failed to meet our significance criterion. Figure 11, A and B, illustrates for two pairs of V1 neurons the differences between CCGs calculated from spontaneous activity (gray traces at bottom) with those for activity evoked by full-contrast gratings, the orientation of which evoked either strong (top traces) or weak (middle traces) mean responses. First, the peak height of CCGs for activity evoked by the best-compromise stimulus was usually either the same as (left column) or larger (right column) than for spontaneous activity. Second, the peak height of CCGs for activity evoked by ineffectively oriented stimuli could be either smaller (left column) or similar (right column) in amplitude to those for spontaneous activity. Finally, the width of spontaneous CCGs was consistently broader than for either high-contrast stimulus.
As suggested by the broad width and substantial amplitude of the spontaneous CCGs, the mean rsc for spontaneous activity in all pairs was relatively strong (0.29 ± 0.04; calculated from 2.56 s epochs). When we limited our analysis to pairs for which the spontaneous CCGs were significant (n = 24), the mean rsc was 0.38, similar to the value at low contrast for pairs with significant CCG peaks and almost twice as high as the correlation of activity evoked by high-contrast stimuli. The mean peak height in this subset of cells was 0.012 coinc/spk, between the mean height for most effective and least effective orientations (0.015 and 0.007 coinc/spk, respectively). Finally, the mean spontaneous CCG peak width was 69 ms, slightly larger than that observed during the presentation of low-contrast stimuli. A direct comparison between the average CCG for the least effective high-contrast stimuli and spontaneous activity recorded in the same pairs is shown in Figure 11C. The greater CCG peak height and width results in a larger rCCG for spontaneous activity (gray trace) at all time scales (Fig. 11D).
Thus, we find that ongoing cortical activity is strongly correlated and that this correlation is not attributable to unspecific changes in cortical excitability. First, correlation of spontaneous activity arises from loose synchrony of spike timing on the time scale of hundreds of milliseconds. Second, synchronized spontaneous activity typically occurs in pairs of neurons whose evoked activity is also synchronized, consistent with previous reports that the spatial structure of activity in cortex is similar for spontaneous and evoked activity (Arieli et al., 1996; Tsodyks et al., 1999). Finally, the signal that gives rise to this correlation can be suppressed by high-contrast stimuli, even when such stimuli evoke little response in the recorded pair.
Relationship between response correlation and receptive field similarity
If the correlation documented in the previous sections is related to common input to the pair of neurons, there should be some relationship between receptive field properties and the overall magnitude of correlation. Specifically, because the RF properties of a neuron are determined primarily by its inputs, cells with similar RF properties are likely to receive more common input than those with dissimilar preferences and should thus be more strongly correlated. Our measurements of orientation, spatial and temporal frequency, and size tuning, used to determine the optimal stimulus parameters for each pair, provided an opportunity to investigate how rsc and spike timing correlation depend on the similarity in tuning between the two cells. A useful measure of tuning similarity is rsignal, which is the correlation between the mean responses of two neurons to the same stimulus set. Because it depends completely on the stimulus set with which it is measured, the magnitude of rsignal is by itself ambiguous (i.e., two neurons can be viewed as having similar or dissimilar tuning depending on the richness of the stimulus set presented), but the relationship between rsignal and rsc remains informative.
We evaluated the relationship between the rsc and rsignal by calculating the signal correlation for each stimulus parameter and comparing it to rsc across our population of pairs. We used the average rsc measured across all orientations or, for pairs for which we only measured responses to stimuli of varying contrasts, at high contrast. We found a positive but weak relationship between rsignal for any particular stimulus parameter and rsc (orientation: r = 0.10, p = 0.21; spatial frequency: r = 0.13, p = 0.12; temporal frequency: r = 0.17, p = 0.05; receptive field size: r = 0.27, p = 0.001), but when we characterized the similarity in tuning with a single rsignal value for all stimulus parameters, the correlation between rsignal and rsc was stronger and statistically significant (r = 0.27, p = 0.001) (Fig. 12). We also found a significant relationship between rsignal and CCG peak height (r = 0.32; p < 0.001). The relationship between rsc and rsignal in V1 (r = 0.27) is similar to that documented for nearby neurons in motor and parietal cortex (∼0.15-0.5) (Lee et al., 1998) and in macaque supplementary motor area (0.35) (Lee, 2003). In MT, Bair et al. (2001) found a stronger correlation of 0.61, but this may be because the pairs were recorded on a single electrode and thus particularly close together (Lee et al., 1998). Our data, on the other hand, may underestimate the relationship between rsignal and rsc in V1, because we recorded few cells with negative signal correlation (Fig. 12, points in the left half) and because we did not explore all stimulus parameters for which V1 neurons are selective (e.g., color and disparity tuning).
In a complementary analysis, we investigated the relationship between tuning similarity and response correlation using a multivariate regression model. The independent parameters in the regression analysis were the differences in preferred orientation, spatial and temporal frequency, and the RF overlap of the two cells (see Materials and Methods). The most informative model included all of the variables and predicted 10% of the variance in rsc (r = 0.31; p = 0.03). The coefficients for each variable were negative so that larger differences in preference predicted lower rsc values, but only the term for orientation preference was significant at the p = 0.05 level (Table 1). Previous studies have suggested that neurons with similar ocular dominance are more likely to fire synchronously than those with different preferences (Kruger and Aiple, 1988; Hata et al., 1991). For a subset of neurons (n = 55), we measured ocular dominance quantitatively and included the difference in the response ratio in the two eyes (in octaves) as an additional variable in the model. This model explained a higher percentage of the variance of rsc (29%; p = 0.02), and orientation and ocular dominance were the most important factors.
We conclude that there is a relationship between the tuning similarity of two cells and both the correlation in their response variability and the strength of spike timing correlation. A simple explanation for these data are that common input to a pair of neurons gives rise to similar tuning but also correlated response variability and a tendency to fire synchronously. It is the modulation of the strength and temporal characteristics of this common input by stimulus orientation and contrast that presumably gives rise to the stimulus-dependent correlation we observe.
Discussion
We have shown that stimulus contrast and orientation have distinct effects on the correlation of single V1 neurons. During spontaneous activity and for weak to moderate stimulus drive, trial-to-trial variability of similarly tuned neurons is often strongly correlated and synchronized with loose temporal precision. High-contrast stimuli reduce correlation on long time scales, make synchronous firing more precise, and either enhance or suppress synchrony strength depending on whether the stimulus drives the pair strongly. Together, our results suggest that correlation arises from two distinct mechanisms: a low-frequency component attributable to ongoing activity, which dominates correlated variability on long time scales and which is suppressed by strong stimulus drive; and a temporally precise component, which is orientation tuned and causes sharp synchronization.
Comparison with previous studies
Previous studies have assumed that spike count correlation in V1 is stimulus independent (Gawne et al., 1996; Reich et al., 2001). We find instead that for response epochs of several seconds, rsc is orientation independent but increases for low-contrast stimuli; on the time scale of tens of milliseconds, rsc is strongest for orientations that drive both cells strongly but primarily independent of stimulus contrast.
Several studies have investigated the dependence of timing correlation between single neurons in V1 on basic stimulus parameters. Snider et al. (1998) found that monosynaptic connections (CCGs with peaks offset from zero) between neurons in cat V1 were more sensitive to stimulus orientation than contrast. Samonds and Bonds (2004) showed that common-input CCGs are larger for effectively oriented stimuli. Llampl et al. (1999) showed that membrane potential fluctuations between pairs of V1 neurons are correlated and stimulus dependent, also consistent with our data. Our study extends these findings by comparing the sensitivity of common input CCGs to stimulus orientation and contrast, by comparing spontaneous and evoked correlation, and by relating the stimulus dependence of synchrony to changes in the spike count correlation over a range of time scales.
Remaining studies of V1 synchrony have used indirect measures or have explored the role of synchrony in binding a distributed representation. Because we used single gratings, our data does not directly evaluate the role of synchrony in binding, but it does suggest that synchrony can be partly predicted by whether a stimulus provides potent common drive to two neurons. We also found substantial synchrony, albeit on a broad time scale, in the absence of any stimulus drive. Both findings suggest that activity will often be synchronized even when there is no need for binding.
Functional implications
Previous measurements of rsc (∼0.2) were based on stimuli presented for several seconds and thus overestimate its value on time scales relevant for sensory information processing (10-100 ms). In a simple pooling model, lower rsc values mean that the population signal-to-noise ratio will improve when a larger number of neurons are averaged (Zohary et al., 1994; Bair et al., 2001). In addition, approximately one-half of our pairs lacked significant synchronization (Fig. 6), so that rCCG on short time scales was nearly zero. Although the preferences of these neurons was more different on average than for neurons that were strongly synchronized, the tuning of the cells was usually sufficiently similar so that a single grating could drive both cells. Any stimulus should thus excite a neuronal population whose variability is essentially uncorrelated during short epochs; pooling over such cells could significantly reduce response variability.
In the pooling model, the fidelity of the neuronal population, the ability to encode time-varying stimuli in the population firing rate, depends critically on the time scale of correlation. Mazurek and Shadlen (2002) showed that for an rsc of 0.2 that arises solely from synchronous spikes, the population is unable to follow temporal modulations over 50 Hz, but this cut-off frequency increases if correlation arises on a slower time scale. We find that the correlation for high-contrast stimuli arises over tens of milliseconds, so that the population could encode stimuli that vary at >200 Hz. At low stimulus contrast, both the time scale and the magnitude of rsc increase, so the fidelity should remain limited.
More sophisticated analysis of population coding has shown that correlation can either improve or impair stimulus encoding, relative to the performance of an uncorrelated population of cells (Vogels, 1990; Abbott and Dayan, 1999; Sompolinsky et al., 2001). The accuracy of the code decreases if similarly tuned neurons are correlated more strongly than dissimilar cells, as we find. For some algorithms for reading population codes, however, stimulus-dependent correlation can improve markedly the accuracy of the population code (Shamir and Sompolinsky, 2004). The net effect of the V1 correlation we report will thus depend on the degree to which subsequent stages of cortical processing can read out stimulus dependent correlation.
Our V1 findings allow us to evaluate the physiological plausibility of current models for reading population codes. Pouget et al. (1998) and Deneve et al. (1999) proposed a recurrent neural network to recover the information available in a population of noisy input neurons. The model predicts that the response co-variance between neurons should be negative for orientations that fall where the slope of the tuning curves are different and positive for stimuli that fall on a common flank. Our results are inconsistent with this suggestion, because correlation is either orientation independent (over seconds) or maximal for stimuli that drive both cells well (on time scales <100 ms), typically at an orientation between the peak preference of the cells, where the slopes of the tuning curves are different. We checked this directly by calculating the average correlation for stimuli that fell where the product of the slopes was either negative or positive but found no difference between the two (0.20 vs 0.22; p = 0.8).
Finally, the stimulus dependence of correlation can be used to constrain models of V1 circuitry. A strictly feedforward model predicts CCGs that are orientation independent (Ben-Yishai et al., 1995), clearly inconsistent with our findings. Models with recurrent cortical connections predict that the tails of the CCG will be negative for stimuli between the peaks of the tuning curves and positive where the slopes of the tuning curves are similar. We did observe a tendency for correlation at long time scales to decrease for the most effective stimuli (Fig. 5), but this effect was weaker than that predicted by Ben-Yishai et al. (1995) and unrelated to the slopes of the tuning curves. Our results are thus not consistent with either implementation of the Ben-Yishai et al. model. More complex models of V1 circuitry, presumably those with realistic patterns of recurrent connectivity and spontaneous cortical activity, may provide a better match.
Mechanisms
Our data suggest that at least two sources of input contribute to correlated responsivity: one of low frequency that is present during spontaneous activity and is sensitive to stimulus contrast and a second that depends on stimulus orientation. The mechanism responsible for the low-frequency component present in spontaneous activity is not clear. In the visual system, it may be partly caused by slow correlation in the activity of retinal ganglion cells (60-100 ms) (Mastronarde, 1989; Brivanlou et al., 1998), which is reduced by increasing light intensity, just as increasing stimulus contrast reduces CCG peak width in V1. Low frequency correlation may also involve correlated transitions to “UP” (depolarized) states in cortical neurons (Llampl et al., 1999). In any case, the low frequency correlation of spontaneous activity is unlikely to be caused by anesthesia because (1) broad CCGs have been reported in the visual cortex of awake monkeys (Gochin et al., 1991); (2) broad correlation of spontaneous subthreshold activity is related to the simultaneously recorded LFP but not to more global measures such as the EEG (Llampl et al., 1999); and (3) the LFP is dominated by low frequencies in both anesthetized and awake animals (Amzica and Steriade, 1995; Llampl et al., 1999; Leopold et al., 2003).
As contrast increases, the CCG peak narrows, and its amplitude either increases or decreases depending on stimulus efficacy. At the network level, tighter synchronization may be related to shifts in the LFP power spectra toward high frequencies (Llampl et al., 1999; Lee, 2003; Henrie, 2004). Cellular mechanisms may also contribute. Azouz and Gray (2003) reported that during depolarization, action potentials are preferentially triggered by high-frequency inputs. Stimulus drive also increases membrane conductance (Borg-Graham et al., 1998; Anderson et al., 2000), thus making the cell less sensitive to low-frequency input. Such changes in conductance also occur for ineffective stimuli, which could explain the disruption of spontaneous synchrony by high-contrast stimuli that fail to evoke a strong response (Hirsch et al., 1998; Anderson et al., 2000). Interestingly, rsc changed primarily between full contrast and the next lowest contrast (typically 0.25), reminiscent of the signal that causes response saturation at high contrasts (Carandini et al., 1997). Finally, the magnitude of sharp synchrony depends in a straightforward way on stimulus orientation: gratings that evoke strong responses result in substantial synchrony; orientations that evoke weak responses result in weak synchrony. Although it is unclear how CCG peak height relates to the amount of common input [see Appendix, available at www.jneurosci.org as supplemental material (references therein: Tolhurst et al., 1983; Vogels et al., 1989)], a simple explanation for the orientation dependence of synchrony is that stimuli that drive two neurons strongly provide more common input that those that evoke weak responses.
Footnotes
This work was supported by a Howard Hughes Medical Institute Investigatorship to J. Anthony Movshon and by a grant to Movshon from the National Institutes of Health (EY2017). Additional support was provided by a grant to Peter Lennie (EY04440) and a Core Grant for Vision Research to New York University (EY13079). We thank Wyeth Bair for helpful discussions and, together with Movshon, Simon Schultz, and Brian Lau, for comments on this manuscript.
Correspondence should be addressed to Adam Kohn, Center for Neural Science, New York University, 4 Washington Place, Room 809, New York, NY 10003. E-mail: adamk{at}cns.nyu.edu.
M. A. Smith's present address: Center for the Neural Basis of Cognition, Carnegie Mellon University, 115 Mellon Institute, 4400 Fifth Avenue, Pittsburgh, PA 15213.
Copyright © 2005 Society for Neuroscience 0270-6474/05/253661-13$15.00/0