The Center for Neuroscience, University of California, 1544 Newton Ct, Davis, CA 95616, USA
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Although studies of neuronal synchronization in monkey striate cortex are few in number, they have yielded a wide range of findings. In anesthetized macaques, correlated firing with millisecond precision has been reported to occur in ~50% of neuronal pairs separated by <1 mm (Krüger and Aiple, 1988; Ts'o and Gilbert, 1988
). Paired recordings in V1 and V2 of anesthetized macaques have demonstrated that synchronization also occurs between these cortical areas (Roe and Ts'o, 1997
; Nowak et al., 1999
), but correlations having millisecond precision were rare. The dominant time scale of correlations in several of these studies was in the range of 20100 ms (Krüger and Aiple, 1988
; Roe and Ts'o, 1997
; Nowak et al., 1999
), suggesting some commonality with the cat (Nelson et al., 1992
; Nowak et al., 1995
). In contrast to these findings, one study in the anesthetized squirrel monkey has shown that correlated firing with millisecond precision is abundant, stimulus dependent and dominated by an oscillatory temporal structure in the gamma frequency band (Livingstone, 1996
). These data are similar to some studies of cat striate cortex (Gray et al., 1989
, 1990
; Engel et al., 1990
; Munk et al., 1996
) and suggest a possible species difference in the temporal structure of striate cortical activity.
Disparate findings have also been reported in two studies of temporal correlation in striate cortex of alert macaques. Frien et al. (Frien et al., 1994) found robust stimulus-dependent synchronization between neurons in V1 and V2 that was dominated by oscillations in the gamma frequency band. Lamme and Spekreijse (Lamme and Spekreijse, 1998
), on the other hand, found stimulus-evoked temporal correlations at a slower time scale (i.e. 2050 ms) that showed no evidence of gamma-band oscillations and only weak variations across differing stimuli. While these data rule out a strict separation on the basis of species, they do suggest that the properties of correlated activity may depend on stimulus parameters and possibly on behavioral task as well. Thus, the picture that emerges from the few studies of anesthetized and alert monkeys is one of widespread variation in the incidence and properties of correlated activity in striate cortex. This fact is responsible in large part for the ongoing debate concerning the functional significance of synchronous activity in the visual system (Ghose and Maunsell, 1999
; Gray, 1999
; Shadlen and Movshon, 1999
; Singer, 1999
).
In an effort to resolve these disparate findings, we have performed multineuron recordings in area V1 of alert macaque monkeys to determine the incidence, stimulus-dependence and temporal properties of local and long range synchronous activity occurring on a millisecond time scale. Our goal was to employ stimuli that would lead to coactivation of pairs of recorded cells, and to compare the properties of synchronous activity occurring spontaneously with those evoked by visual stimulation. To accomplish this, we utilized single gratings to coactivate pairs of cells with overlapping receptive fields, and either large single gratings or pairs of independent gratings to coactivate cells with nonoverlapping receptive fields. We varied the direction, spatial frequency or drift speed of the gratings in order to evaluate responses to stimuli that produced the greatest coactivation of the cells. Using these stimuli, and cross-correlation analysis of the spike trains, we find that synchronized firing on a millisecond time scale occurs in a substantial fraction of cell pairs during the response to visual stimuli. Ninety percent of this synchronized firing occurs when one or both cells of a pair display gammaband oscillations in firing probability, and correlation strength is closely correlated with oscillation amplitude. These findings demonstrate that millisecond time scale synchronization is a robust and stimulus-dependent feature of striate cortical activity in the alert macaque, and suggest that neuronal rhythmicity is an important mechanism contributing to its generation (König et al., 1995a; Gray and McCormick, 1996
).
The results of this study have been reported previously in abstract form (Friedman-Hill et al., 1995).
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The methods for surgical preparation, behavioral training and electrophysiological recording are described in detail in the accompanying paper (Friedman-Hill et al., 2000). All surgical and experimental techniques were in accordance with NIH guidelines and were approved by the University of California at Davis animal care and use committee.
Visual Stimuli
Our goal was to evaluate the incidence and properties of temporally correlated activity under conditions where the cells were coactivated by visual stimuli. To accomplish this, we utilized drifting sine-wave or square-wave gratings presented on a dark background (Friedman-Hill et al., 2000). A small portion of the data from the initial phase of the experiments was collected using drifting bars. During each session, our recording configuration (Friedman-Hill et al., 2000
) resulted in two nonoverlapping populations of cellular receptive fields in the lower visual hemifield that were typically separated from each other by 23°. We attempted to coactivate as many cells as possible in the two groups by presenting large grating stimuli that covered both sets of receptive fields and pairs of gratings positioned to activate each set of receptive fields simultaneously. We used different combinations of orientation, direction of movement and spatial frequency in an effort to activate as many cells as possible. Once this procedure was completed, we then presented grating patches over each set of receptive fields separately and varied the direction, spatial frequency and drift speed over a range of values to enable the calculation of tuning curves for these parameters. Because recording stability varied from session to session, only a fraction of the cell pairs in the study were subjected to the complete range of stimulus parameters.
Data Analysis
Spike separation, peristimulustime histograms (PSTH), autoand crosscorrelation histograms (ACH, CCH), including trial-shuffled controls (shift-predictor), associated power spectra and statistical confidence limits, were computed as described elsewhere (Gray et al., 1995; Gray and Viana Di Prisco, 1997
; Friedman-Hill et al., 2000
). All correlation histograms (single trial and cumulative across trials) had time lags of ±128 ms and a bin width of 1 ms, and were computed over two userdefined epochs, one prior to (win1) and another during (win2) the response to the stimulus. The former ranged in duration from 0.2 to 1.0 s and the latter from 1.5 to 2.5 s. The occurrence and properties of oscillatory firing were evaluated according to the methods described in Friedman-Hill et al. (Friedman-Hill et al., 2000
).
Our principal objective in quantifying response synchronization was to establish an unbiased measure of the magnitude and statistical significance of correlated firing that occurred within ±10 ms of zero time-lag. In particular, we sought to avoid assumptions about the underlying temporal structure of the spike trains, and to determine the magnitude of correlated firing above and beyond that which might have been introduced by the stimulus. We therefore explicitly avoided using spectral analysis or curve-fitting methods to detect correlation peaks, and utilized the variance in the shift-predictor control correlograms to establish the confidence limits for significance. Our measure, referred to as the significance ratio (SR), is the ratio of two integral values: A peak value (P), representing the magnitude of spike correlation, computed by taking the sum of the bins in the central 20 ms of the CCH that exceed the 95% confidence limits, and a variance value (V), representing the expected occurrence of coincident spikes, computed from the sum of the central 20 ms in each histogram lying between the 95% confidence limits and the mean value of the correlogram. This ratio measure is illustrated in Figure 1 and was computed as follows:
|
![]() | (1) |
where
![]() | (2) |
and
![]() | (3) |
![]() | (4) |
and where
![]() | (5) |
and
![]() | (6) |
![]() | (7) |
and
![]() | (8) |
![]() | (9) |
where X is the mean value of all the bins in the CCH, is the standard deviation of all the bins in the corresponding shift-predictor histogram, and i and j range from 10 to +10 ms. Consequently, this measure treats negative and positive deviations from the correlation mean with equal weight. By using the standard deviation of the shift-predictor histogram to establish the confidence limits for significance, we were not only able to estimate the correlated firing that occurred beyond that introduced by the stimulus, but also to estimate significance utilizing spike trains whose temporal structure was identical to that of the signals being evaluated. This measure increased the stringency of our significance test, but it may also have contributed to the occurrence of false negatives. However, we considered it preferrable to err in this direction rather than falsely assign correlograms as significant.
Because the SR measure yielded values over a wide range, it became necessary to establish a separate confidence limit for the statistical significance of each SR value. To do this, we utilized the Monte Carlo procedure described in the preceding report (Friedman-Hill et al., 2000). In this analysis, however, we employed two different simulations of the experimental data. For each experimental spike train, we computed an equivalent pseudorandom spike train (i.e. random samples taken from a uniform distribution, with no spike times having the same value) and an interval-shuffled spike train obtained by randomly shuffling the time of occurrence of the interspike intervals in the experimental data. Both simulations yielded control spike trains identical in spike count, mean firing rate and duration to the experimental data. The latter simulation, however, preserved the interspike interval distribution of the original data but resulted in a random shuffling of interval times. All of the results in this study are based on the confidence limits computed from the pseudorandom spike trains. The interval-shuffled spike trains are used to control for the possibility that the temporal structure of the firing patterns might yield statistically significant correlations of a spurious origin.
As with the autocorrelation analysis, both simulations were repeated 500 times for each window of data analyzed on each trial. The CCH was computed from each simulated spike train and corresponding cumulative session correlograms were calculated. The resulting correlograms were edge corrected as described earlier (Friedman-Hill et al., 2000), and the SR value was then computed for each CCH. This yielded two sets of 500 control SR values (pseudorandom and ISI-shuffled) for each experimental CCH [i.e. single trial (win1,win2), session (win1,win2) and shiftpredictor (win1,win2)]. We assigned a confidence limit for statistical significance by choosing the SR values in the control distributions that were greater than 99% of all the values. This resulted in two different cutoff values, one computed from the pseudorandom spike trains [cutoff(rnd)] and the other computed from the interval-shuffled spike trains [cutoff(isi)]. Any experimental SR value equal to or greater than either of these cutoff values was considered statistically significant at a P-value of <0.01. As with the autocorrelation analysis, we calculated the ratio of the SR/cutoff to provide a measure of correlation strength that was independent of the total spike count in the data. Values of this ratio >1.0 were considered significant. Our use of two different Monte Carlo procedures to estimate statistical significance further enabled us to compare the incidence and magnitude of significant correlation in our sample between pseudorandom and interval-shuffled controls.
For those correlograms classified as significant, we calculated the width and time lag of the corresponding peaks utilizing the full CCH (i.e. ±128 ms time lag). Full widths at half height (i.e. midpoint between peak value and correlogram mean) of CCH peaks were measured after subjecting each CCH to a low pass, boxcar filter with a 3 ms bin width. In order to measure the time-lag of each significant CCH, we fit a Gaussian function to the central peak using the maximum likelihood method (Maldonado and Gray, 1996). The time-lag was computed by taking the difference between the peak of the Gaussian function and the 0 ms bin. Power spectra and peak frequencies of the CCHs were computed and tested for significance as described in the preceding report (Friedman-Hill et al., 2000
).
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Incidence of Response Synchronization
After separating the spike trains arising from different cells (Gray et al., 1995), our sample consisted of 608 paired recordings. For each pair, we used from one to 32 different stimuli varying in spatial frequency, direction and speed. We selected for analysis the stimulus condition that elicited the highest combined firing rate of the cells, computed by taking the square root of the product of the mean firing rates of the two spike trains. Using this selection criterion, and the confidence limits derived from the pseudorandom spike trains, the incidence of significant response synchronization is summarized in Table 1
. Several notable findings are apparent from these results. Consistent with previous correlation measurements in the cat and monkey, the incidence of significant correlation decreased as a function of spatial separation between the cells: 38% of all shortrange pairs and 17% of all long-range pairs exhibited significant synchronized firing in response to the stimulus condition that evoked the highest combined firing rates of the cells. Similarly, the strength of correlation also decreased with spatial separation. The mean correlation magnitude was significantly greater (P < 0.01, U-test) for short-range as compared with long-range CCHs (Figure 2A
).
|
|
|
|
In addition to the cumulative measure of cross-correlation, we also found that synchronous activity was frequently of sufficient strength to be detected on single trials, and therefore did not require extensive averaging to be resolved. We applied our analysis to the data collected on single trials and found a pattern of results similar to that seen in the cumulative measurements (Figure 2B). Table 2
shows summary statistics for this analysis on the basis of the maximum combined firing rate criterion. Overall, the incidence of single-trial correlations was approximately one-third of that observed in the session measurements, with short-range correlations being more prevalent and greater in magnitude than long-range correlations. These differences were due in part, however, to the increase in variance relative to the mean that is present in the single-trial CCHs. As a result, our measure of significance is likely to have led to an underestimate of significant correlation when the number of sampled spikes was small. This effect is illustrated in Figure 5
, where the CCHs for 10 consecutive trials (B) are shown for a short-range pair that displayed a strong correlation in the cumulative correlogram (A). Each of the 10 single-trial CCHs displays a clear indication of response synchronization, but the central peaks in four of the correlograms (indicated by *) failed to reach statistical significance.
|
|
|
Dependence of Temporal Correlations on Orientation and Direction Preference
In both the cat and the monkey, it has been reported that correlated firing occurs most often and with greatest magnitude among cells that are tuned to the same orientation and direction (Ts'o et al., 1986; Ts'o and Gilbert, 1988
; Schwarz and Bolz, 1991
). In our own experience, however, we have found that significant correlations can occur among cells with widely different orientation preferences as long as they can be coactivated, particularly with a single stimulus (Gray et al., 1989
; Engel et al., 1990
; Kreiter and Singer, 1996
). To address this issue, we examined the dependence of short-range temporal correlation on orientation and direction preferences by presenting a single drifting grating over a range of 360° (22.5° steps) in pseudorandom order. For each cell in a pair, we fit a Gaussian function, utilizing the maximum likelihood method (Maldonado and Gray, 1996
), to a 180° region of the tuning curve centered over the stimulus that elicited the maximum response and determined the difference in orientation preference for each pair of tuned cells. We then analyzed the responses to all 16 stimuli to determine if at least one of them evoked significant correlated firing (Table 3
). Of 105 pairs examined, 72 (~69%) displayed significant peaks in their CCHs. The majority of these 72 cell pairs differed in their orientation preferences by <45°. This bias in our sample is likely to be due to the physical arrangement of our short-range measurements in which the recorded cells were separated by no more than 400 µm. When cells did differ widely in their preferred orientations, it was often difficult to coactivate them with a single drifting grating. The resulting low firing rates are likely to have further biased our estimate of correlation because of the sensitivity of our measure to spike count. In spite of these factors, we did find a number of significantly correlated cell pairs in which the preferred orientation of the cells differed by 4590°. Interestingly, the percentage of significant correlations did not differ substantially as a function of the difference in orientation preference (Table 3
).
|
|
Stationary versus Drifting Stimuli
In the accompanying paper (Friedman-Hill et al., 2000), we reported that stationary stimuli elicited lower frequency oscillatory responses when compared with those evoked by drifting stimuli. This implies that if response synchronization is to be maintained for both stationary and drifting stimuli then oscillatory responses must covary in frequency among cells that are coactivated by the same stimulus. We addressed this issue by comparing the frequency of oscillatory synchronization evoked by stationary and moving gratings. Of 133 cell pairs investigated, 17 (12.8%) showed significant synchronization for both stationary and drifting stimuli. Another 21 (15.7%) pairs showed correlation only in response to a stationary grating, and 7 (5.2%) pairs showed significant correlation only during the presentation of a drifting grating. For those cell pairs that exhibited synchronous oscillations in response to both stimuli, the mean frequency was significantly higher when the stimulus was drifting (45 ± 7 Hz) than when stationary (33 ± 2 Hz; P << 0.001, paired t-test).
Correlation Peak Width and Time Lag
Previous studies, employing cross-correlation analysis in areas 17 and 18 of the cat and monkey, have emphasized the presence of at least three distinct temporal scales of neuronal synchronization (Krüger and Aiple, 1988; Nelson et al., 1992
; Nowak et al., 1994
, 1995
, 1999
). By computing the width of correlation peaks, these authors have demonstrated that correlated firing tends to occur on time scales of ~120, 20100 and 1001000 ms. To determine if our data could be categorized in a similar manner, we computed the width at half-height of all the significant correlation peaks in our sample. Figure 8A
shows the resulting distribution of peak widths for longand short-range pairs. We found a clustering of values below 12 ms and little or no evidence for peak widths >15 ms (mean values were 6.2 ± 1.7 ms and 5.7 ± 1.6 ms for short- and long-range pairs, respectively). Because the maximum time-lag of our calculations did not exceed 128 ms, we were unable to assess the presence of correlation peaks broader than ~200 ms. The distribution of correlation widths in our data clearly reflects the characteristic temporal structure of the synchronous neuronal oscillations found in our sample.
|
In spite of the clustering of time-lags within 23 ms, we also found a substantial fraction of cell pairs with time-lags in the range of 49 ms. We suspected that these large delays may have resulted from the combined dependence on stimulus orientation and orientation preference reported in a recent study (König et al., 1995b). These authors demonstrated that optimally driven neurons tend to fire slightly before suboptimally activated cells, with the magnitude of the difference being dependent on both the differences in preferred orientation of the cells and the orientation of the stimulus. To examine this relation, we measured the time-lag as a function of stimulus orientation among significantly correlated short-range pairs of cells having overlapping receptive fields. The data were taken from paired recordings during the execution of direction tuning curves. For each CCH having a significant peak, we determined the time-lag of the peak and then computed the slope of the linear fit to the data points comparing time-lag with stimulus orientation. We found that all the slopes were positive and their values were significantly correlated with the difference in preferred orientation of the cell pairs. An example of this result is illustrated in Figure 9A
. Here, the CCHs of a neuronal pair stimulated by a drifting grating presented at four different directions show a continuous shift in time-lag. This transition occurs as the stimulus changes from the nonpreferred direction of the reference cell (Figure 9A
, top) to the preferred direction (Figure 9A
, bottom). The magnitude of the shift for a given change in stimulus orientation differed among cell pairs. The rate of change of the time-lag with respect to stimulus orientation was greater for cell pairs with large differences in orientation preference. Plotting the slope of the change in time-lag as a function of the orientation preference difference for the entire sample revealed a positive correlation (Figure 9B, P
< 0.05), indicating that the magnitude of the shift depends on both the stimulus orientation and the difference in orientation preference of the cell pair. These results are consistent with those reported by König et al. (König et al., 1995b
), where a similar dependence on stimulus orientation was found among neuronal pairs recorded in cat striate cortex.
|
In our sample of significant CCHs, we observed that nearly two-thirds of the long-range pairs and three-quarters of the short-range pairs exhibited oscillatory modulation (Table 1E). This suggested that synchronization might be related to the occurrence of oscillatory firing, as has been recently demonstrated in the cat (König et al., 1995a
) and the squirrel monkey (Livingstone, 1996
). To test this conjecture, we examined the ACHs of each pair of significantly correlated cells in our sample to determine if the occurrence of oscillatory firing was related to response synchronization. For each neuronal pair, one of three relationships could be true: both ACHs could be oscillatory, one cell of a pair could be oscillatory and the other nonoscillatory, or both could be nonoscillatory. If response synchronization occurred with equal probability among oscillatory and nonoscillatory cells, we would expect the relation between oscillatory firing and synchronization to simply reflect the incidence of rhythmic (56.3%) and nonrhythmic (43.7%) firing in our sample of single-unit and multiunit recordings (Friedman-Hill et al., 2000
). Specifically, among those pairs of cells that exhibit synchronization, ~32% (0.563 x 0.563) would be derived from cells that are both oscillatory, ~49% [(0.563 x 0.437) x 2] would contain one site that is oscillatory and ~19% [1 (0.32 + 0.49)] would come from cell pairs that are both nonoscillatory. On the other hand, if synchrony were more likely to occur when the cells exhibit oscillatory firing, we would expect the incidence of synchronization to be highest when the cells are firing rhythmically and lowest when they show no rhythmic modulation.
Our findings were consistent with the latter interpretation (Figure 10A). The incidence of significant synchronization (52%) was highest when both members of a pair were oscillatory, decreased to 38% when only one member of a pair was oscillatory and fell to 10% when neither site exhibited oscillatory firing (black-filled bars). Consistent with this finding, uncorrelated firing was least likely to occur when both members of a pair were oscillatory, and was most commonly observed when one or neither member of a pair was firing rhythmically (hatched bars). Using a Chi-square test, we found the correlation between oscillation and synchronization to be highly significant (P < 0.00001), and a similar result held when we restricted the analysis to pairs of single units only (P < 0.024). These data reveal that oscillatory firing, as detected in the ACH, was a strong predictor of the occurrence of significant peaks in the CCH. Significant crosscorrelation occurred more often than predicted by chance if both ACHs in a pair were oscillatory (i.e. 52% versus the predicted 32%).
|
Our finding that oscillatory firing was a strong predictor of synchronization raised the possibility that independent oscillators could, in principle, show significant correlation simply because they have similar frequencies. To address this issue, we modified the previous analysis by first selecting those pairs where both cells displayed oscillatory firing and then determined the incidence of significant correlation among these pairs. For the correlograms meeting this condition, 71% (64/90) of the short-range pairs and 30% (28/93) of the long-range pairs exhibited significant synchronous firing. This demonstrates that, although correlation probability is greatest when the participating cells fire rhythmically, oscillatory responses are by no means obligatorily coupled. They can occur simultaneously and remain temporally independent even when their frequencies are very similar. An example illustrating this result is shown in Figure 11. The data were taken from four recordings where the cells were simultaneously activated by two gratings drifting in the same direction (A). The cells recorded on channels 2, 3 and 4 were sampled from either the same or adjacent tetrodes and had overlapping receptive fields, while the activity recorded on channel 1 was sampled from a separate tetrode located 4 mm distant and its receptive field did not overlap with the others. All of the cells had preferred orientations that spanned a range of <30°. Strong oscillatory responses of similar mean frequency were evoked in all four spike trains during the presentation of the stimuli (B). However, while the activity on channels 2, 3 and 4 showed clear synchronization (D), the oscillatory firing between these cells and the activity on channel 1 was independent (C).
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Methodological Considerations
At least two factors may have affected our ability to assess the properties of synchronous activity. First, several forms of uncontrolled sampling bias might have affected our data, as discussed in the comanion paper (Friedman-Hill et al., 2000). Secondly, our measure of synchronization may also have introduced errors in our analysis. In designing this measure (the significance ratio), we attempted to minimize this influence in several ways. First, we made no assumptions regarding the temporal strucuture of the underlying activity. We simply searched for joint firing within the central 20 ms of the CCH that was greater and/or less than that estimated from the confidence limits set by the trialshuffled CCH. This provided a way to estimate the joint firing exceeding that introduced by the stimulus and by any common temporal structure of the signals. We then determined the significance of the SR values by comparing them to distributions of surrogate SR values computed from equivalent pseudorandom and interval-shuffled spike trains. Consequently, two strict criteria had to be met for an SR value to be considered significant. When we applied these methods, the quantitative assignment of significant correlations closely matched the assignment that each of us made by visually inspecting the cross-correlograms. We are therefore reasonably confident that our measure of statistical significance is rigorous and accurate. By limiting our estimate of correlated firing to the central 20 ms of the CCH, we might have easily missed correlation peaks whose time lags were outside this window, and we would have incorrectly estimated the magnitude of correlation peaks whose widths exceeded this window. However, careful visual inspection of all the CCHs in our data set showed no evidence of the former, and our post hoc analysis of the peak widths revealed no evidence of the latter within the time lags (±128 ms) of our CCH calculation.
Comparison to Other Studies of Synchronization in Monkeys
In several respects, our findings are both similar to and different from previous studies of neuronal synchronization in primate visual cortex. We have summarized these studies in Table 4, which reveals a surprisingly broad spectrum of findings for such a small number of reports. Our results most closely resemble those of Frien et al. and Livingstone (Frien et al., 1994
; Livingstone, 1996
), who reported abundant neuronal synchronization on a fast time scale that was both stimulus dependent and dominated by rhythmic firing in the gamma frequency band. Our findings differ, however, from several of the other studies of V1. For example, we found no evidence for synchronization at the intermediate time scale (20100 ms), while several studies have not reported evidence for synchronous gamma-band activity. Although it is difficult to fully account for these differences, there are several factors that are likely to contribute. The absence of intermediate time scale correlations in our data might simply reflect the abundance of gamma-band activity we observed (Friedman-Hill et al., 2000
). Synchronous oscillatory firing at 3060 Hz would naturally lead to correlation peak widths on the order of ~15 ms or less. If for some reason our data were biased in favor of gamma-band activity [see discussion in (Friedman-Hill et al., 2000
)], we would expect to see less of the slower time scale correlations.
|
Differences in analysis techniques may provide another source of variation. In some studies, cross-correlograms were computed from long epochs that combined periods of spontaneous and evoked activity across different stimuli (Aiple and Krüger, 1988; Krüger and Aiple, 1988
; Ts'o and Gilbert, 1988
). Several studies indicate that pooling data in this way can lead to erroneous estimates of the properties of synchronization (Engel et al., 1990
; Kreiter and Singer, 1996
; Livingstone, 1996
; Gray and Viana Di Prisco, 1997
; Nowak et al., 1999
). We attempted to reduce these factors by analyzing spontaneous and stimulus-evoked activity separately. Differences in the methods used to estimate the statistical significance of correlations are likely to provide a second source of analytic variation. Here the number of methods employed have been surprisingly diverse for such a small number of studies (Table 4
). Some studies simply reported the incidence of peaked cross-correlograms without providing a measure of their statistical significance. Other studies applied statistical confidence limits to assess correlated firing, but the details of these methods vary widely across studies. Needless to say, these factors make direct comparisons difficult, and given the wide variety of both qualitative and quantitative techniques, it is perhaps not surprising that such a wide range of results would occur.
Two other factors may contribute to the differences between studies. First, cross-correlation measurements from a number of cortical areas indicate that the temporal dynamics of synchronous activity vary across cortical areas (Gochin et al., 1991; Frien et al., 1994
; Freiwald et al., 1997
; Roe and Ts'o, 1997
; Nowak et al., 1999
). Thus, it is not entirely unexpected that the properties of correlated firing would differ between V1 and IT, for example. Second, variation in behavioral state is also likely to influence the properties of synchronous activity. One striking example of this is apparent between two studies of area MT. In alert fixating macaques, Kreiter and Singer (Kreiter and Singer, 1996
) reported that synchronization was stimulus dependent. They showed changes in synchronized firing that depended on whether the cells were stimulated by one or two drifiting bars. In contrast, de Oliveira et al. (de Oliveira et al., 1997
) showed that neuronal synchronization was suppressed by visual stimulation during a motion discrimination task. This suggests that task demands may have a profound impact on synchronization dynamics in area MT.
When each of these factors are taken into consideration, it is not possible to reconcile all the reported differences in results between the studies listed in Table 4. There are common findings in some of the studies, but there are also widespread differences. The broad divergence of findings suggests that agreement between studies will only come when efforts are made to apply common experimental and analytical methods to the investigation of specific cortical areas.
Comparison of Response Synchronization in Cat and Monkey
In spite of the complexity of correlation dynamics in the monkey, as well as in the cat, a number of our findings are comparable to those of previous studies. As has been reported numerous times, we find that the incidence and magnitude of correlated firing decreases with the spatial separation between the cells (Toyama et al., 1981a,b
; Michalski et al., 1983
; Ts'o et al., 1986
; Ts'o and Gilbert, 1988
; Gray et al., 1989
; Engel et al., 1990
; Schwarz and Bolz, 1991
; Hata et al., 1993
; Livingstone, 1996
), and that short-range synchronization occurs with roughly equal probability between cells of all orientation preferences (Engel et al., 1990
, 1991c
; Das and Gilbert, 1999
). These data are roughly consistent with the patterns of intracortical connections in striate cortex, where the density of horizontal connections falls off sharply with distance, and where local connections tend to lack the bias for connecting regions of similar orientation preference that is present in the long range connections (Gilbert and Wiesel, 1983
, 1989
; Rockland and Lund, 1983
; Kisvarday and Eysel, 1992
; Lund et al., 1993
; Malach et al., 1993
; Yoshioka et al., 1996
; Levitt et al., 1996
; Kisvarday et al., 1997
).
We also find that the properties of synchronous gamma-band activity in monkey striate cortex are very similar to those observed in the same area of the cat. In both species, rhythmically firing cells tend to fire in high-frequency bursts (Gray et al., 1990; Friedman-Hill et al., 1996
; Livingstone, 1996
). The probability of occurrence, the distributions of frequency and phase, and the stimulus dependence of synchronous oscillations are similar in the two species (Engel et al., 1990
; Fries et al., 1997
; Gray and Viana Di Prisco, 1997
). Synchronous gamma-band activity rarely occurs in the absence of visual stimulation, can occur over a wide range of stimulus orientations, directions and speeds (Engel et al., 1990
; Freiwald et al., 1995
; Gray and Viana Di Prisco, 1997
; Maldonado and Gray, 1997
), and exhibits specific changes in phase related to stimulus orientation and the orientation preferences of the cells (König et al., 1995b
). Together, these data suggest a close correspondence in the mechanisms generating synchronous gamma-band activity in the cat and monkey (Gray and McCormick, 1996
), and imply that this form of activity may be a general feature of mammalian striate cortex, including the human (Sem-Jacobsen et al., 1956
; Chatrian et al., 1960
; Tallon-Baudry et al., 1996
, 1997a
, Tallon-Baudry et al., b
).
Correlation between Synchrony and Oscillations
Finally, as has been reported in both the cat (König et al., 1995a) and the squirrel monkey (Livingstone, 1996
), we have found that the probability and strength of response synchronization is directly correlated with rhythmic firing in the gamma frequency band. These data indicate that, although oscillatory firing is neither necessary nor sufficient for synchronization, it is a good predictor of its occurrence. This raises the possibility that oscillatory firing may contribute to the generation of synchronous activity. How this might be mediated is largely unknown, but the recent discovery of superficial pyramidal neurons in striate cortex that are capable of intrinisically generating oscillatory activity suggests one possible mechanism (Gray and McCormick, 1996
). When depolarized above firing threshold, these neurons, referred to as chattering cells, fire brief bursts of 25 action potentials that are followed by a powerful afterhyperpolarization (Wang, 1999
). Both events may be particularly effective in coordinating the activity of other neurons. Bursts of spikes can lead to temporal summation postsynaptically (Miles and Wong, 1986
) and an increase in presynaptic release probability (Thomson et al., 1993
; Stevens and Wang, 1995
), each of which is likely to produce an enhanced postsynaptic response. The afterhyperpolarization following each burst may act synergistically with these effects by reducing the effectiveness of nonsynchronous inputs. The additional finding that chattering cells have horizontal axon collaterals that project within superficial cortical layers (Gray and McCormick, 1996
) places these neurons in a key position for generating synchronized activity. Evidence for chattering cells in the striate cortex of monkeys is limited and indirect (Friedman-Hill et al., 1996
, 2000
), thus their role in synchronization remains speculative until further evidence can be obtained.
![]() |
Notes |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Address correspondence to Charles M. Gray, Ph.D., Center for Computational Biology, Lewis Hall, Montana State University, Bozeman, MT 59717-3505, USA. Email: cmgray{at}nervana.montana.edu.
![]() |
Footnotes |
---|
2 Present address: Laboratory of Neuropsychology, NIMH, NIH, Bldg 49, Room 1B80, 49 Convent Drive, MSC 4415, Bethesda MD 208924415, USA
3 Present address: Center for Computational Biology, Lewis Hall, Montana State University, Bozeman, MT 59717-3505, USA
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Chatrian GE, Bickford RG, Uihlein A (1960) Depth electrographic study of a fast rhythm evoked from the calcarine region by steady illumination. Electroenceph Clin Neurophysiol 12:167176.[ISI]
Das A, Gilbert CD (1999) Topography of contextual modulations mediated by short-range interactions in primary visual cortex. Nature 399:655661.[ISI][Medline]
de Oliveira SC, ThieleA, Hoffmann KP (1997) Synchronization of neuronal activity during stimulus expectation in a direction discrimination task. J Neurosci 17:92489260.
Eckhorn R, Frien A, Bauer R, Woelbern T, Kehr H (1993) High frequency (6090 Hz) oscillations in primary visual cortex of awake monkey. NeuroReport 4:243246.[ISI][Medline]
Engel AK, König P, Gray CM, Singer W (1990) Stimulus-dependent neuronal oscillations in cat visual cortex: inter-columnar interaction as determined by cross-correlation analysis. Eur J Neurosci 2:588606.[ISI][Medline]
Engel AK, Kreiter AK, König P, Singer W (1991a) Synchronization of oscillatory neuronal responses between striate and extrastriate visual cortical areas of the cat. Proc Natl Acad Sci USA 88:60486052.[Abstract]
Engel AK, König P, Kreiter AK, Singer W (1991b) Interhemispheric synchronization of oscillatory neuronal responses in cat visual cortex. Science 252:11771179.[ISI][Medline]
Engel AK, König P, Singer W (1991c) Direct physiological evidence for scene segmentation by temporal coding. Proc Natl Acad Sci USA 88: 91369140.[Abstract]
Freiwald WA, Kreiter AK, Singer W (1995) Stimulus dependent intercolumnar synchormization of single unit responses in cat area 17. NeuroReport 6:23482352.[ISI][Medline]
Freiwald WA, Kreiter AK, Singer W (1997) Synchronous and oscillatory activity in the macaque inferotemporal cortex. In: Göttingen Neurobiology Report 1997, p. 550. Stuttgart: Thieme.
Friedman-Hill S, Maldonado PE, Gray CM (1995) Short and long-range response synchronization in macaque striate cortex. Soc Neurosci Abstr 21:592.3.
Friedman-Hill SR, Gray CM, Maldonado PE (1996) Incidence of oscillations and synchrony as a function of bursting properties of macaque striate neurons. Soc Neurosci Abstr 22:117.6.
Friedman-Hill S, Maldonado PE, Gray CM (2000) Temporal dynamics of neuronal activity in the striate cortex of alert macaque: I. Incidence and stimulus-dependence of oscillations. Cereb Cortex 10:11051116.
Frien A, Eckhorn R, Bauer R, Woelbern T, Kehr H (1994) Stimulus-specific fast oscillations at zero phase between visual areas V1 and V2 of awake monkey. NeuroReport 5:22732277.[ISI][Medline]
Fries P, Roelfsema PR, Engel AK, König P, Singer W (1997) Synchronization of oscillatory responses in visual cortex correlates with perception in interocular rivalry. Proc Natl Acad Sci USA 94: 1269912704.
Gawne TJ, Richmond BJ (1993) How independent are the messages carried by adjacent inferior temporal cortical neurons? J Neurosci 13:27582771.[Abstract]
Ghose GM, Maunsell J (1999) Specialized representations in visual cortex: a role for binding? Neuron 24:7985.[ISI][Medline]
Gilbert CD, Wiesel TN (1983) Clustered intrinsic connections in cat visual cortex. J Neurosci 3:11161133.[Abstract]
Gilbert CD, Wiesel TN (1989) Columnar specificity of intrinsic horizontal and corticocortical connections in cat visual cortex. J Neurosci 9: 24322442.[Abstract]
Gochin PM, Miller EK, Gross CG, Gerstein GL (1991) Functional interactions among neurons in inferior temporal cortex of the awake macaque. Exp Brain Res 84:505516.[ISI][Medline]
Gray CM (1999) The temporal correlation hypothesis of visual feature integration: still alive and well. Neuron 24:3147.[ISI][Medline]
Gray CM, McCormick DA (1996) Chattering cells: superficial pyramidal neurons contributing to the generation of synchronous oscillations in the visual cortex. Science 274:109113.
Gray CM, Viana Di Prisco G (1997) Stimulus-dependent neuronal oscillations and local synchronization in striate cortex of the alert cat. J Neurosci 17:32393253.
Gray CM, König P, Engel AK, Singer W (1989) Stimulus-specific neuronal oscillations in cat visual cortex exhibit inter-columnar synchronization which reflects global stimulus properties. Nature 338:334337.[ISI][Medline]
Gray CM, Engel AK, König P, Singer W (1990) Stimulus-dependent neuronal oscillations in cat visual cortex: receptive field properties and feature dependence. Eur J Neurosci 2:607619.[ISI][Medline]
Gray CM, Engel AK, König P, Singer W (1992) Synchronization of oscillatory neuronal responses in cat striate cortex: temporal properties. Vis Neurosci 8:337347.[ISI][Medline]
Gray CM, Maldonado PE, Wilson M, McNaughton B (1995) Tetrodes markedly improve the reliability and yield of multiple single unit isolation from multiunit recordings in cat striate cortex. J Neurosci Methods 63:4354.[ISI][Medline]
Hata Y, Tsumoto T, Sato H, Tamura T (1993) Horizontal interactions between visual cortical neurons studied by cross-correlation analysis in the cat. J Physiol 441:593614.[Abstract]
Herculano-Houzel S, Munk MHJ, Neuenschwander S, Singer W (1999) Precisely synchronized oscillatory firing patterns require electroencephalographic activation. J Neurosci 19:39924010.
Kisvarday AF, Eysel UT (1992) Cellular organization of reciprocal patchy networks in layer III of cat visual cortex (Area 17). Neuroscience 46:275286.[ISI][Medline]
Kisvarday ZF, Toth E, Rausch M, Eysel UT (1997) Orientation-specific relationship between populations of excitatory and inhibitory lateral connections in the visual cortex of the cat. Cereb Cortex 7:605618.[Abstract]
König P, Engel AK, Singer W (1995a) The relation between oscillatory activity and long-range synchronization in cat visual cortex. Proc Natl Acad Sci USA 92:290294.[Abstract]
König P, Engel AK, Roelfsema PR, Singer W (1995b) How precise is neuronal synchronization. Neural Computat 7:469485.[Abstract]
Kreiter AK, Singer W (1996) Stimulus-dependent synchronization of neuronal responses in the visual cortex of the awake macaque monkey. J Neurosci 16:23812396.[Abstract]
Krüger J, Aiple F (1988) Multimicroelectrode investigation of monkey striate cortex: spike train correlations in the infragranular layers. J Neurophysiol 60:798828.
Lamme VAF, Spekreijse H (1998) Neuronal synchrony does not represent texture segregation. Nature 396:362366.[ISI][Medline]
Levitt JB, Lund JS, Yoshioka T (1996) Anatomical substrates for early stages in cortical processing of visual information in the macaque monkey. Behav Brain Res 76:519.[ISI][Medline]
Livingstone MS (1996) Oscillatory firing and interneuronal correlations in squirrel monkey striate cortex. J Neurophyiol 75:24672485.
Lund JS, Yoshioka T, Levitt JB (1993) Comparison of intrinsic connectivity in different areas of macaque monkey cerebral cortex. Cereb Cortex 3:148162.[Abstract]
Malach R, Amir Y, Harel M, Grinvald A (1993) Relationship between intrinsic connections and functional architecture revealed by optical imaging and in vivo targeted biocytin injections in primate striate cortex. Proc Natl Acad Sci USA 90:1046910473.[Abstract]
Maldonado PE, Gray CM (1996) Heterogeinity in local distributions of orientation-selective neurons in the cat primary visual cortex. Vis Neurosci 13:509516.[ISI][Medline]
Maldonado PE, Gray CM (1997) Frequency of neuronal oscillations in monkey striate cortex is modulated by stimulus velocity. Soc Neurosci Abstr 23:13.13.
Michalski A, Gerstein GL, Czarkowska J, Tarnecki R (1983) Interactions between cat striate cortex neurons. Exp Brain Res 51:97107.[ISI][Medline]
Miles R, Wong RKS (1986) Excitatory synaptic interactions between CA3 neurones in the guinea-pig hippocampus. J Physiol 373:397418.[Abstract]
Milner P (1974) A model for visual shape recognition. Psychol Rev 81:521535.[ISI][Medline]
Munk MH, Roelfsema PR, König P, Engel AK, Singer W (1996) Role of reticular activation in the modulation of intracortical synchronization. Science 272:271274.[Abstract]
Nelson JI, Salin PA, Munk MHJ, Arzi M, Bullier J (1992) Spatial and temporal coherence in cortico-cortical connections: a crosscorrelation study in areas 17 and 18 of the cat. Vis Neurosci 9:2137.[ISI][Medline]
Nowak LG, Munk MHJ, Chounlamountri N, Bullier J (1994) Temporal aspects of information processing in areas V1 and V2 of the macaque monkey. In: Oscillatory event-related brain dynamics (Pantev C, ed.), pp. 8598. New York: Plenum Press.
Nowak LG, Munk MH, Nelson JI, James AC, Bullier J (1995) Structural basis of cortical synchronization. I. Three types of interhemispheric coupling. J Neurophysiol 74:23792400.
Nowak LG, Munk MHJ, James AC, Girard P, Bullier J (1999) Crosscorrelation study of the temporal interactions between areas V1 and V2 of the macaque monkey. J Neurophysiol 81:10571074.
Rockland KS, Lund JS (1983) Intrinsic laminar lattice connections in primate visual cortex. J Comp Neurol 216:303318.[ISI][Medline]
Roe AW, Ts'o DY (1997) The functional architecture of area V2 in the macaque monkey: physiology, topography and connectivity. In: Cerebral cortex: extrastriate cortex in primates (Kaas JH, Rockland KL, Peters AL, eds). New York: Plenum Press.
Roelfsema PR, Engel AK, König P, Singer W (1997) Visuomotor integration is associated with zero time-lag synchronization among cortical areas. Nature 385:157161.[ISI][Medline]
Schwarz C, Bolz J (1991) Functional specificity of the long-range horizontal connections in cat visual cortex: a cross-correlation study. J Neurosci 11:29953007.[Abstract]
Sem-Jacobsen CW, Petersen MC, Dodge HW, Lazarte JA, Holman CB (1956) Electroencephalographic rhythms from the depths of the parietal, occipital and temporal lobes in man. Electroenceph Clin Neurophysiol 8:263278.[ISI]
Shadlen MN, Movshon JA (1999) Synchrony unbound: a critical evaluation of the temporal binding hypothesis. Neuron 24:6777.[ISI][Medline]
Singer W (1999) Neuronal synchrony: a versatile code for the definition of relations. Neuron 24:4965.[ISI][Medline]
Singer W, Gray CM (1995) Visual feature integration and the temporal correlation hypothesis. Annu Rev Neurosci 18:555586.[ISI][Medline]
Stevens CF, Wang Y (1995) Facilitation and depression at single central synapses. Neuron 14:795802.[ISI][Medline]
Tallon-Baudry C, Bertrand O, Delpuech C, Pernier J (1996) Stimulus specificity of phase-locked and non-phase-locked 40 Hz visual responses in human. J Neurosci 16:42404249.
Tallon-Baudry C, Bertrand O, Delpuech C, Permier J (1997a) Oscillatory gamma-band (3070 Hz) activity induced by a visual search task in humans. J Neurosci 17:722734.
Tallon-Baudry C, Bertrand O, Wienbruch C, Ross B, Pantev C (1997b) Combined EEG and MEG recordings of visual 40 Hz responses to illusory triangles in human. NeuroReport 8:11031107.
Thomson AM, Deuchars J, West DC (1993) Single axon excitatory postsynaptic potentials in neocortical interneurons exhibit pronounced paired pulse facilitation. Neuroscience 54:347360.[ISI][Medline]
Toyama K, Kimura M, Tanaka K (1981a) Cross-correlation analysis of interneuronal connectivity in cat visual cortex. J Neurophysiol 46: 191201.
Toyama K, Kimura M, Tanaka K (1981b) Organization of cat visual cortex as investigated by cross-correlation techniques. J Neurophysiol 46: 202214.
Ts'o DY, Gilbert CG (1988) The organization of chromatic and spatial interactions in the primate striate cortex. J Neurosci 8:17121727.[Abstract]
Ts'o D, Gilbert C, Wiesel TN (1986) Relationship between horizontal interactions and functional architecture in cat striate cortex as revealed by cross-correlation analysis. J Neurosci 6:11601170.[Abstract]
von der Malsburg C (1981) The correlation theory of brain function. Internal report, Max-Planck-Institute for Biophysical Chemistry, Göttingen.
von der Malsburg C. (1985) Nervous structures with dynamical links. Ber Bunsenges Phys Chem 89:703710.[ISI]
Wang X-J (1999) Fast burst firing and short-term synaptic plasticity: a model of neocortical chattering neurons. Neuroscience 89:347362.[ISI][Medline]
Yoshioka T, Blasdel GG, Levitt JB, Lund JS (1996) Relation between patterns of intrinsic lateral connectivity, ocular dominance, and cytochrome oxidase-reactive regions in macaque monkey striate cortex. Cereb Cortex 6:297310.[Abstract]