Dynamics of Striate Cortical Activity in the Alert Macaque: I. Incidence and Stimulus-dependence of Gamma-band Neuronal Oscillations

Stacia Friedman-Hill1, Pedro E. Maldonado2 and Charles M. Gray3

The Center for Neuroscience, University of California, 1544 Newton Ct, Davis, CA 95616, USA


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Notes
 References
 
Using single and multiunit recordings in the striate cortex of alert macaque monkeys, we find that gamma-band (20–70 Hz) oscillations in neuronal firing are a prominent feature of V1 neuronal activity. The properties of this rhythmic activity are very similar to those previously observed in the cat. Gamma-band activity is strongly dependent on visual stimulation, largely absent during spontaneous activity and, under the conditions of our experiment, not time-locked to the vertical refresh of the computer monitor (80 Hz) used to present the stimuli. In our sample, 61% of multiunit activity (MUA) and 46% of single-unit activity (SUA) was significantly oscillatory, with mean frequencies of 48 ± 9 and 42 ± 13 Hz, respectively. Gamma-band activity was most likely to occur when cells were activated by their optimal stimuli, but still occurred, although less often and with lower amplitude, in response to nonoptimal stimuli. The frequency of gamma-band activity also reflected stimulus properties, with drifting gratings evoking higher-frequency oscillations than stationary gratings. As in the cat, the spike trains of single cells showing gamma-band oscillations often displayed a pattern of repetitive burst firing, with intraburst firing rates of 300–800 Hz. The overall similarity of rhythmic neuronal activity in the primary visual cortex of cats and monkeys suggests that the phenomenon is not species-specific. The stimulus-dependence of the rhythmic activity is consistent with a functional role in visual perception.


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Notes
 References
 
Pattern recognition by the mammalian visual system is thought to be preceded by a rapid, preattentive process that groups related features in an image. Because most visual images evoke activity that is distributed over multiple cortical locations, it has been recognized that some mechanism is needed to identify members of a representation as belonging together. One hypothesized mechanism for this grouping process posits that cortical neurons responding to common features exhibit temporally correlated firing on a millisecond time scale (Milner, 1974Go; von der Malsburg, 1981Go, 1985Go; Singer, 1993Go; Singer and Gray, 1995Go). Support for this hypothesis has been obtained from experiments demonstrating that synchronous firing in the visual cortex is stimulus-dependent (Gray and Singer, 1987Go, 1989Go; Eckhorn et al., 1988Go; Gray et al., 1989Go; Engel et al., 1990Go, 1991aGo, Engel et al., bGo; Freiwald et al., 1995Go; Kreiter and Singer, 1996Go; Livingstone, 1996Go; Fries et al., 1997Go; Gray and Viana Di Prisco, 1997Go).

These studies have led to vigorous debates concerning the origin, properties and functional significance of synchronous activity (Gray, 1999Go; Shadlen and Movshon, 1999Go). However, much of the controversy surrounding this debate stems from a confusion of two separate, but related findings. First, the studies cited above have revealed that neurons, located within the same or different areas of the visual cortex, are capable of synchronizing their responses with millisecond precision when activated by appropriate visual stimuli. Second, the same studies show that synchronous activity is often associated with pronounced oscillations of firing probability in the gamma frequency range (20–70 Hz). These findings have led to the proposal that gamma-band activity may contribute to the generation of neuronal synchrony (Gray et al., 1989Go; König et al., 1995Go; Gray and McCormick, 1996Go; Livingstone, 1996Go). But they have also led to the erroneous notion that the properties of synchronous activity can be inferred from the measurement of gamma-band oscillations at single recording sites (Ghose and Freeman, 1992Go; Tovee and Rolls, 1992Go; Young et al., 1992Go; Bair et al., 1994Go). The ensuing debate has led to the emergence of two further questions: is gamma-band activity a general feature of neuronal firing in mammalian visual cortex, and if so, what are its properties? And what is the relationship between neuronal synchronization on a millisecond time scale and gamma-band activity? The current paper addresses the first of these questions, while the companion paper (Maldonado et al., 2000Go) focuses on the relationship between gamma-band oscillations and synchronous activity.

While a number of studies have reported common features of gamma-band activity in cat visual cortex (Gray and Singer, 1987Go, 1989Go; Eckhorn et al., 1988Go; Gray et al., 1990Go; Fries et al., 1997Go; Gray and Viana Di Prisco, 1997Go; Herculano-Houzel et al., 1999Go) [but see (Ghose and Freeman, 1992Go)], there is no general agreement on the nature of this activity in the visual cortex of nonhuman primates. Some researchers have explicitly reported a failure to find gamma-band activity in monkey visual cortex (Tovee and Rolls, 1992Go; Young et al., 1992Go; Bair et al., 1994Go; de Oliveira et al., 1997Go; Nowak et al., 1999Go). Other groups have found oscillatory activity, but of a higher frequency (60–90 Hz) than in cat (30–60 Hz) (Eckhorn et al., 1993Go; Frien et al., 1994Go). Still others have reported the presence of oscillatory firing in monkey visual cortical areas, with frequency ranges comparable to those documented in cat (Kreiter and Singer, 1992Go, 1996Go; Livingstone, 1996Go).

It is difficult to reconcile these divergent findings because the studies vary in (i) the species of monkey used; (ii) the experimental paradigm employed; (iii) the methods of data analysis; and (iv) the cortical area investigated. For example, while both Young et al. and Livingstone (Young et al., 1992Go; Livingstone, 1996Go) recorded from the striate cortex of anesthetized monkeys, the first study used macaque monkeys and analyzed cumulative autocorrelograms whereas the second study utilized squirrel monkeys and analyzed single-trial autocorrelograms. Eckhorn et al. and Frien et al. (Eckhorn et al., 1993Go; Frien et al., 1994Go) have been the only reports of gamma-band activity in the striate cortex of the alert monkey; however, these investigations employed only one animal, and the neural activity was sampled from large groups of multiple units or from the local field potential. In an effort to resolve these differences, we have analyzed single and multiunit activities recorded from the striate cortex of two alert macaque monkeys. Our aim was to determine the incidence and properties of gamma-band activity, and thereby clarify whether this form of activity is species-specific or is a more general phenomenon. We focused our investigations on the primary visual cortex because the original reports were based on data collected from areas 17 and 18 of the cat. And we utilized alert macaque monkeys because these animals are most commonly used in behavioral paradigms. In this study, we demonstrate that oscillatory firing in the gamma frequency band (20–70 Hz) is a robust property of visually evoked neuronal activity in the striate cortex of monkeys performing a visual fixation task.

The results of this study and the companion paper (Maldonado et al., 2000Go) have been reported previously in abstract form (Friedman-Hill et al., 1995Go; Gray et al., 1995aGo).


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Notes
 References
 
Subjects and Surgical Procedures

Two adult, female rhesus monkeys (Macaca mulatta) served as subjects for this study. The monkeys underwent two sterile surgical procedures to prepare them for training and recording. In the first procedure, we implanted a pair of scleral search coils for monitoring eye position (Judge, et al., 1980Go) and a stainless steel post for head restraint (Gray and Viana Di Prisco, 1997Go). Following the behavioral training described below, the monkeys underwent a second surgical procedure, in which a hard plastic recording chamber was mounted over the opercular surface of striate cortex and secured to the skull with orthopedic screws and dental acrylic. A craniotomy (5 x 10 mm) was made in the bone overlying one hemisphere. The animals were given ~20 days to recover from each surgical procedure before behavioral training or recording was initiated. All surgical procedures were performed at the California Regional Primate Research Center (CRPRC) and were supervised by the attending veterinary staff. Surgical and experimental techniques were in accordance with institutional and NIH guidelines.

Behavioral Training

The monkeys were trained to maintain their gaze within 1.0° of a fixation spot, in the presence of moving or stationary visual stimuli, for a period of up to 3 s. Successful trials were rewarded with a drop of dilute apple juice. During behavioral training and recording, the monkeys' access to water was restricted. On days the monkeys were working, they were given a minimum of 30 ml/kg/day of water or juice in the form of fixation rewards or supplemental water. In practice, the monkeys earned an average of 50–80 ml/kg/day. On Fridays and Saturdays, the monkeys received 80 ml/kg/day of water. Throughout the entire experimental period the animals were monitored by the veterinary staff of the CRPRC. The monkeys were given supplemental water if any of the routinely monitored physiological parameters (i.e. plasma osmolality, urine specific gravity, body weight) became unacceptable.

Recording Techniques

Neuronal signals were recorded with multiple (2–6), tungsten microelectrodes (1–2 M{Omega} resistance) or custom-fabricated nichrome tetrodes (Gray et al., 1995bGo), separated by either 250 µm or 3–4 mm. The signals were amplified (10K), bandpass filtered (0.6 kHz–6kHz), digitized (30 kHz/channel) and stored for off-line spike sorting and analysis. We recorded both multiunit activity (MUA) and isolated, single-unit activity (SUA) from both tetrodes and electrodes. On individual penetrations, units were recorded whenever the signal-to-noise ratio exceeded 3:1. After a recording was completed at a given site, the electrodes or tetrodes were always moved at least 200 µm before another recording was obtained. On successive penetrations occurring in the same guide tube, activity was always sampled at least 200 µm deeper than the previous recording session. This sampling procedure was continued until activity could no longer be measured at a given guide tube location, at which time the guide tubes were repositioned. Sometimes the sampling procedure across sessions involved recording from a second layer of V1 in the head of the calcarine sulcus. This was confirmed by changes in receptive field position.

After recording, individual units were resolved from the MUA through principle components analysis of waveforms (Abeles and Goldstein, 1977Go) or by using the stereo properties of tetrodes to cluster spikes by peak-to-peak amplitude, peak time, spike width, etc. (Gray et al., 1995bGo). The extracted spike trains were stored at 1.0 ms resolution. We carried out a separate analysis of the incidence of oscillatory activity for MUA and SUA. Because some of the single units were recorded simultaneously from a single electrode or tetrode and separated off-line, we also pooled these units and included the pooled data in our MUA analysis. In our discussion below, a ‘site’ refers to an electrode or tetrode position. Cells simultaneously recorded from different electrodes or tetrodes are considered to be different sites; cells simultaneously recorded from the same electrode or tetrode are considered to be the same site.

Visual Stimuli

Visual stimuli, consisting of single or multiple, drifting square-wave or sine-wave gratings (0.5–3.0 cycles/°, 8–12 cd/m2 luminance) presented on a dark background, were generated by a personal computer and displayed on a 19-inch color monitor (80 Hz noninterlaced, 1024 x 768 resolution). The vertical refresh of the video monitor was synchronized with the data acquisition at a resolution of 300 µs. The stimuli were presented 5–30 times during each recording session while the animals maintained ocular fixation on a central target. At the start of each trial, the monkey was given 0.5–1 s to look at the fixation point. Data collection began once the target was acquired and continued throughout the period of fixation. If the monkey failed to move its eyes to the fixation point within the required time or broke fixation early, the trial was aborted and replaced with a new trial following a delay. On each successful trial, the stimulus was presented at a latency of 0.25–0.75 s for a duration of 2.0–2.5 s. The interval between the end of one trial and the beginning of the next varied from 1.0 to 4.0 s in duration. The monkey was rewarded with apple juice for each complete trial, but was not penalized for incomplete trials.

Prior to the automated presentation of stimuli, receptive fields were hand-mapped with computer-generated light or dark bars or gratings, while the monkey maintained fixation. The minimum response field was determined following the methods of Barlow et al. (Barlow et al., 1967Go). We were usually able to make rough estimates of orientation, direction, velocity and spatial frequency preferences during the mapping routine; this information was used to choose a range of stimuli that was subsequently presented in automated tuning curve sequences.

Anatomical Localization

Cells in V1 and V2 were distinguished on the basis of electrode location and depth, as well as published reports of receptive field size versus eccentricity (Hubel and Wiesel, 1977Go; Gattass et al., 1981Go). Because of the angle and location of our penetrations, all cells initially encountered within a penetration were from the opercular surface of V1. As the electrodes were advanced, the amount of change in receptive field location and size was used to determine whether a second folding of V1 had been penetrated, or whether the electrode had moved into V2. All cells in our database were from the lower visual field and had receptive fields within 13° of the center of gaze.

Data Analysis

For each recording and stimulus, we computed the peri-stimulus time histogram (PSTH, bin width 50 ms). Using this measure, two epochs were defined to calculate spontaneous (win1) and stimulus-evoked (win2) activity (Gray and Viana Di Prisco, 1997Go). Mean firing rates were computed within these epochs for each stimulus and were used to construct tuning curves. Direction tuning curves were fit by a Gaussian function and deemed to be ‘tuned’ or ‘untuned’ using the methods employed by Maldonado and Gray (Maldonado and Gray, 1996Go). Interspike interval histograms (ISIH) were calculated with a bin width of 1 ms and a maximum time lag of 200 ms.

To study rhythmic neuronal firing, we computed autocorrelation histograms (ACH) for individual trials and for cumulative data (summed across trials) (Perkel et al., 1967Go) with time lags of ±128 ms. These calculations were performed separately for both spontaneous (win1) and stimulus-evoked (win2) activity. Correlated firing time-locked to the stimulus was estimated by computing the shift-predictor control correlogram. Each of these histograms (single trial, session and shift-predictor) was then corrected to enable a comparison of the data across cells and stimuli: (i) each bin in the ACH was normalized by dividing by the central bin; (ii) the central bin in the ACH was then replaced by the mean of the two adjacent bins, in order to eliminate a large, sharp peak at zero time lag; and (iii) a correction was applied to each ACH to compensate for the edge effects that arise when a correlation histogram is computed on windowed data. Each bin of the ACH was multiplied by the factor [window length/(window length – absolute value of the time lag)]. In general, this correction had very little effect on the ACHs computed from the stimulus-evoked data, because these data were collected with long window durations (i.e. 1.0–2.0 s) compared with the time-lag of the ACH. All histograms depicted in the figures accompanying this paper and the companion paper are uncorrected.

To quantify neuronal rhythmicity, we calculated the power spectrum of each ACH (±128 ms) and extracted the frequency and amplitude of the peak value in the frequency range of 20–80 Hz (Gray and Viana Di Prisco, 1997Go). The statistical significance of the spectral peaks was estimated using a Monte Carlo simulation. For this procedure, we generated an equivalent pseudorandom spike train for each window of data collected on each trial using a random sample from a uniform distribution and a 1 ms refractory period. The ACH and corresponding power spectrum were computed for each simulated data set and the peak value in the spectrum between 20 and 80 Hz was extracted. The simulation was repeated 500 times and any experimental value was considered significant if it was greater than the largest simulated value (the ‘cutoff’). The ratio of the experimental peak to the simulated cutoff served as a measure of oscillation strength or amplitude. We utilized this ratio (peak/cutoff) because it provided a scale that was easy to understand: any ratio >1.00 was statistically significant. A ratio of 3.00 exceeds the cutoff for significance by 300%; a ratio of 1.10 corresponds to a correlogram with very weak oscillations, just visible upon inspection. Spectral amplitudes are expressed as this ratio throughout the study and all statistics regarding oscillation strength were performed with these ratios. This enabled us to compare oscillation strength across different cells or stimuli independent of the total spike count in each correlogram. The increased stringency of the method also allowed us to include ACHs with very low spike counts and avoided the need for assigning an arbitrary spike count cutoff that excludes valuable data from the analysis (Gray and Viana Di Prisco, 1997Go).


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Notes
 References
 
Incidence of Gamma-band Activity

The data set included in this study consisted of 357 MUA and 169 SUA recordings taken from striate cortex, and a further 33 MUA and 12 SUA recordings sampled from V2. Rhythmic firing of single units or groups of cells was often obvious in the raw data. Figure 1AGo illustrates the oscillatory activity of a single neuron when stimulated by a drifting square wave grating of optimal orientation and spatial frequency. This cell fired repetitive bursts of spikes during visual stimulation, as can be seen in the expanded trace on the bottom of Figure 1AGo. The ISIH in Figure 1CGo has a peak at 2 ms, representing the intraburst firing rate of the cell, and a second peak at 14 ms, corresponding to intervals between the last spike in one burst and the first spike in a subsequent burst. Because the number of spikes within a burst and the length of the burst may vary, the second peak of the ISIH is only a rough estimate of the interburst firing rate. We refer to cells with bursting properties similar to the cell shown in Figure 1Go as ‘chattering cells’. We have documented the existence of these cells in extracellular (Gray et al., 1990Go; Gray and Viana Di Prisco, 1997Go) and intracellular (Gray and McCormick, 1996Go) recordings from cats and in extracellular recordings from awake monkeys (Friedman-Hill et al., 1996Go). The autocorrelation histogram of this cell exhibits robust oscillations at a frequency of 52 Hz (Figure 1D,EGo).



View larger version (27K):
[in this window]
[in a new window]
 
Figure 1. Oscillatory activity of a single, isolated V1 cell. (A) The activity recorded during a single stimulus presentation is shown at low (above) and high (below) time resolution. The brackets indicate the interval shown at the expanded time scale. (B) Peri-stimulus–time histogram of the average response. The binwidth for this and all other PSTHs in this paper is 50 ms. (C) Cumulative interspike interval histogram for all trials. (D) Cumulative ACH calculated from the spike trains during the response to the visual stimulus. The horizontal lines indicate the mean ± 2 SD of the control correlogram for this same data set. The number above the calibration bar on the left refers to spike count. The binwidth is 1 ms. All correlograms depicted in this paper use this convention. (E) The power spectrum of the ACH. The horizontal lines indicate the mean and the mean + 3 SD of the power spectrum of a correlogram computed from an equivalent random spike train. The number above the calibration bar in this and subsequent figures reflects the magnitude of power at each frequency as a percentage of the DC component. The numbers to the right side of the graph indicate the amplitude (A, expressed in units of the significance cutoff) and the frequency (F, in Hz) of the spectral peak. All other power spectra in this paper utilize the same convention.

 
In most recordings displaying oscillations, we found that the rhythmicity was easily detected and quantified from the data sampled on single trials. Figure 2Go shows 10 trials of activity recorded from the same cell as shown in Figure 1Go. The cumulative ACH and its shift-predictor control are shown in A and B, respectively. The absence of any peaks in the shift-predictor demonstrates that the oscillations are not time-locked to the stimulus. This excludes the possibility that the rhythmic firing is driven by the video refresh (80 Hz) of the computer monitor. Panels C and D show the ACHs and their associated power spectra computed from 10 repetitions of the same stimulus. Strong oscillatory activity near 50 Hz is present on each trial with some fluctuation of frequency and amplitude from trial to trial. Compared with the episodic and nonstationary neuronal oscillations previously seen in alert cats (Gray and Viana Di Prisco, 1997Go), we found much more robust and stable oscillatory activity in the alert monkey. The greater incidence and stability of oscillatory activity in the alert macaque may reflect genuine differences between cats and monkeys or may have resulted from differences in stimuli [drifting bars (Gray and Viana Di Prisco, 1997Go) versus drifting gratings (present study)].



View larger version (38K):
[in this window]
[in a new window]
 
Figure 2. Oscillatory activity is clearly present in individual trial data and is not due to screen refresh. (A) The cumulative ACH obtained from the same cell as shown in Figure 1DGo. (B) The trial-shuffled ACH from this cell is flat, indicating that oscillations are not time-locked to the stimulus, as would be expected if they reflected monitor flicker. (C) ACHs from 10 individual trials. (D) Power spectra of each of the single-trial ACHs shown in C. Spectral amplitudes >1.0 are considered statistically significant.

 
We visually inspected the ACH and power spectra of all the MUA and SUA sites in order to evaluate the descriptive power of the significance measure (i.e. peak/cutoff) derived from the Monte Carlo analysis. The ACHs and power spectra from five recordings are plotted in Figure 3Go, along with the oscillation amplitude, measured by the peak/cutoff, and the frequency. The correlograms range in strength from very weak oscillations (Figure 3AGo) to very robust rhythmic activity (Figure 3EGo). The corresponding peak/cutoff provides an accurate description of the rhythmic structure of the ACH. The weakly oscillatory ACH (Figure 3AGo) just barely passes the Monte Carlo cutoff, with a significance ratio of 1.09, while the cell with very strong oscillations (Figure 3EGo) exceeds the cutoff by almost 800%. Examples of moderate-to-strong oscillations (Fig 3GoB–D) also exhibit a good correspondence between the ACH and the significance measure. For comparison, we have included the peak/DC measure used in our earlier study in the alert cat (Gray and Viana Di Prisco, 1997Go). This parameter expresses the peak value in the power spectrum in the 20–70 Hz frequency range relative to the magnitude of the DC component.



View larger version (30K):
[in this window]
[in a new window]
 
Figure 3. ACHs and associated power spectra for 5 single or multiunit recordings illustrating gradations in the strength of oscillatory activity. The amplitude of oscillation (A = peak/cutoff) ranges from just greater than detection threshold (A) to very strong (E).

 
Of the total sample of 357 MUA records in V1, 218 (61%) showed significant oscillations in their cumulative ACHs. The median peak/cutoff was 2.88 and the mean oscillation frequency was 48 ± 9 Hz (Figure 4A,BGo and Table 1Go). Oscillations were less likely to occur in SUA autocorrelograms (46%; 78 of 169 cells) than in the MUA sample, with a median peak/cutoff of 1.68 and a mean frequency of 42 ± 13 Hz (Figure 4C,DGo). Thus, the data shown in Figure 3CGo are representative of the overall sample. It is also evident that the range of oscillation amplitudes and frequencies for MUA and SUA are similar (Figure 4GoA–D), although there are more MUA sites with oscillation amplitudes >5.0 and a greater percentage of SUA sites with oscillation frequencies <40 Hz.



View larger version (18K):
[in this window]
[in a new window]
 
Figure 4. Histograms of the amplitudes and peak frequencies for all significant cumulative V1 ACHs. MUA data are shown in A and B, and SUA data are shown in C and D. The numbers in the upper right of each graph indicate the median of the peak/cutoff (for graphs on the left side) or the mean and SD of the peak frequency (for graphs on the right side). The number of significant and total correlograms included in the analysis are also indicated in A, C.

 

View this table:
[in this window]
[in a new window]
 
Table 1 Summary statistics of the incidence of significant oscillatory firing among MUA and SUA recordings in areas V1 and V2
 
We found comparable results in our small sample of V2 recordings. Ten out of 33 (30%) MUA ACHs had significant oscillatory peaks. The mean frequency and median peak/cutoff were 43 ± 9 Hz and 2.10, respectively. Four of the 12 (33%) V2 units were significantly oscillatory, with a median peak/cutoff of 1.34 ± 0.16 and a mean frequency of 40 ± 5 Hz.

There were 571 shift-predictor correlograms, but only 27 (5%) of these had significant spectral peaks, and they came primarily from sessions with small numbers of trials and strongly oscillating cells. For some sessions with many stimuli (e.g. orientation tuning curves), the monkey completed only five trials per stimulus or a cell was lost abruptly because of movement of the animal. With small numbers of trials, the data are not adequately shuffled. This is particularly problematic for strongly oscillatory neurons. The significant shift-predictors had spectral peaks ranging in frequency from 39 to 54 Hz, well below the 80 Hz refresh rate of the monitor. These data strongly suggest that the oscillations seen in the experimental correlograms (Figure 4GoA–D) were not driven, in a time-locked fashion, by the visual stimulus.

It has been suggested that neuronal oscillations reflect spontaneous rhythmicity that is unrelated to visual stimulation (Ghose and Freeman, 1992Go, 1997Go). We addressed this question in two ways. First, we computed ACHs for 200–500 ms windows prior to the onset of the visual stimulus. Of 571 spontaneous ACHs, only 29 (5%) were significantly oscillatory. These spontaneous oscillations had a mean frequency of 36 ± 10 Hz (range: 23–66 Hz) and a median peak/cutoff of 1.36 (range: 1.02–4.85). Second, in a subset of our recordings, we included fixation trials in which no stimulus was presented. This enabled us to compare spontaneous and visually evoked activity recorded under comparable conditions. From this sample, 14 of the sites displayed vigorous gamma-band activity in response to an optimal drifting grating. None of these recordings, however, showed significant oscillations during an equivalent period of spontaneous firing. The mean peak/cutoff values for spontaneous and stimulus-evoked activity were 0.61 ± 0.41 and 4.57 ± 2.28, respectively (P < 0.0001, paired t-test). Together these data indicate that striate cortical neurons rarely exhibit gamma-band activity in the absence of a stimulus, and when they do it is usually of low amplitude.

In addition to the cumulative autocorrelograms, we also computed the ACHs for individual trials (Table 2Go). From the total sample obtained in V1, 25% of the MUA trials and 9% of the SUA trials displayed significant oscillations. For V2, 9% of the MUA trial correlograms and 1% of SUA trial correlograms were significant. Figure 5Go plots the amplitudes and frequencies for individual trials for striate MUA (A,B) and SUA (C,D). The strength of oscillations for single trial ACHs was somewhat weaker than the cumulative ACHs (compare Figure 5Go with Figure 4Go), but significant rhythmic activity was still present for single trials. Moreover, the distribution of oscillation frequencies for trial and session ACHs were also similar.


View this table:
[in this window]
[in a new window]
 
Table 2 Summary statistics of the incidence of significant oscillatory firing among MUA and SUA recordings in areas V1 and V2 measured on single trials
 


View larger version (19K):
[in this window]
[in a new window]
 
Figure 5. Histograms of the amplitudes and peak frequencies for all significant single-trial V1 ACHs. MUA data are displayed in A and B, and SUA data are shown in C and D. Numbers in the upper right of each graph follow the same convention as in Figure 4Go.

 
Throughout the recording phase of these experiments, we developed the subjective impression that oscillatory activity evoked by drifting stimuli was most prominent in the 40–50 Hz frequency range. This impression was confirmed by our finding that the mean oscillation frequency was near 48 Hz. Lowfrequency (20–30 Hz) oscillations were not usually significant for single trial correlograms, but when accumulated across many trials, weak oscillations sometimes reached statistical significance (Figure 4DGo). However, as illustrated in Figure 6Go, the strongest oscillatory activity in our sample was in the 40–50 Hz frequency range. These results reveal that the most robust oscillations, those with spectral amplitudes >5.0, correspond to gamma activity in the 30–60 Hz range. Thus, the rhythmic activity occurring at the limits of the frequency distribution in our sample (i.e. 20–30 and 60–80 Hz) tended to be lower in amplitude.



View larger version (18K):
[in this window]
[in a new window]
 
Figure 6. Oscillation amplitude is plotted as a function of oscillation frequency for all significant correlograms. The strongest oscillatory activity occurs within the 30–60 Hz range.

 
On the basis of earlier studies (Engel et al., 1990Go; Gray et al., 1990Go), and evidence presented below, it could be argued that oscillatory firing reflects activity in the cortical network when the cells are driven to fire at high rates. Using this line of reasoning, we should expect to see a greater incidence of significant oscillatory activity in those cells that fire at high rates. We examined this question by comparing the firing rates of the oscillatory and nonoscillatory single units in our sample of 169 V1 cells. The results are shown in Figure 7Go. Both oscillatory (Figure 7AGo) and nonoscillatory (Figure 7BGo) cells displayed a wide range of firing rates. Although oscillatory activity tended to occur in neurons firing at higher mean rates, significant oscillations were not limited solely to cells firing at high rates; there were many units that fired rhythmically at rates of <20 spikes/s. The distributions of firing rates for oscillatory and nonoscillatory units differed significantly (P < 0.0001, Mann–Whitney U-test), perhaps reflecting the fact that many of the oscillatory units are likely to be chattering cells with high intraburst firing rates (Friedman-Hill et al., 1996Go). Figure 7CGo plots oscillation strength (peak/cutoff) as a function of firing rate for significantly oscillatory single units. From this figure, it is clear that absolute firing rate and oscillation strength were largely independent: There were units with high response rates but very weak oscillatory activity, and there were units which fired very rhythmically but at low rates.



View larger version (13K):
[in this window]
[in a new window]
 
Figure 7. Histograms of firing rate distributions for SUA that exhibited significant oscillations (A) and cells with no significant oscillatory activity (B). Oscillation amplitude (peak/cutoff) was independent of firing rate for the population of units with significant ACHs (C).

 
Because other researchers have observed a difference in oscillation frequency between subjects (M.S. Livingstone, personal communication), we compared the data from the two monkeys that participated in this study. We found small, but significant differences in the incidence, strength and frequency of rhythmic activity in the two animals. In one animal, 66% of the recordings were significantly oscillatory, while in the other animal 45% of the recordings were rhythmic. The amplitude and frequency values were also different. In the monkey with the higher percentage of oscillatory cells, the ACHs had a mean peak/cutoff of 3.99 and a mean frequency of 45.4 Hz. In the second monkey, the mean peak/cutoff was 2.87 and the mean frequency was 46.5 Hz. Both of these differences were statistically significant (peak/cutoff P < 0.002; frequency P < 0.046, Mann–Whitney U-test). These data indicate that the properties of gamma-band activity can vary from animal to animal, but it is difficult to estimate the extent to which sampling biases may have contributed to the differences.

Oscillations Reflect Orientation and Direction Preference

Aside from documenting the existence of rhythmic neuronal activity in the visual cortex, another major aim of this study was to determine how the presence or absence of oscillations was related to orientation and direction preference. Although this issue has been previously investigated in the anesthetized cat (Gray and Singer, 1987Go, 1989Go; Eckhorn et al., 1988Go; Gray et al., 1990Go), it has received much less attention in the awake monkey. Only one study has examined the influence of stimulus direction or orientation on oscillatory activity in the alert monkey and this report included only 11 MUA sites recorded in a single monkey (Eckhorn et al., 1993Go). Our sample consisted of 188 directiontuning curves. For the majority of sessions, 16 directions of a drifting square-wave grating were presented in randomly interleaved trials. For one session, only eight directions of motion were tested.

The effect of changing stimulus orientation on oscillatory activity is illustrated in Figure 8Go for a direction-selective V1 single unit. The ACHs of the six stimuli that encompassed the preferred direction are shown on the left side of the figure (B); the corresponding power spectra appear on the right (C). Stimuli that caused a vigorous response (e.g. 3 and 4) evoked strong oscillations, while nonoptimal stimuli elicited fewer spikes (e.g. 1 and 6) and yielded weakly oscillatory or nonoscillatory correlograms. Interestingly, the strongest oscillations (as quantified by the peak/cutoff) occurred in response to the stimulus which was 22° away (stimulus 4) from the ‘optimal’ stimulus (i.e. the peak of the tuning curve). Thus, the stimulus that elicited the most spikes was not always the stimulus that evoked the strongest oscillations.



View larger version (32K):
[in this window]
[in a new window]
 
Figure 8. Direction-tuning curve for a direction-selective V1 single-unit recording (A). The numbers 1–6 above the horizontal axis indicate the stimuli to which the ACH (B) and power spectra (C) correspond.

 
In order to understand the effect of orientation preference on the probability and strength of oscillatory firing for our entire sample, we compared correlograms for preferred and nonpreferred stimuli, for sessions displaying significant oscillations in response to the optimal stimulus. If rhythmic firing is completely independent of the stimulus, then one might expect to find significant correlograms for all stimuli capable of evoking a visual response. On the other hand, if rhythmic neuronal activity reflects stimulus preferences, then one might expect that optimal stimuli would be associated with a higher incidence of oscillatory activity than nonoptimal stimuli. We identified sites that were orientation tuned and which were oscillatory for at least the orientation that elicited the maximal firing rate. Of the 180 tuning curves, 123 were tuned for orientation; of these, 76 (62%) demonstrated oscillations for the optimal orientation. Figure 9AGo illustrates that the probability of significant oscillatory firing in this latter population decreased with deviation from the preferred orientation. Stimuli at the preferred orientation were nearly seven times more likely to evoke oscillatory firing than at the orthogonal orientation.



View larger version (20K):
[in this window]
[in a new window]
 
Figure 9. Percentage and strength of oscillatory activity, as a function of orientation. (A) Incidence of oscillatory activity reflects orientation tuning preferences, with non-optimal stimuli yielding fewer significant correlograms than optimal or near-optimal stimuli. The analysis was limited to sites that were: (1) orientation tuned and (2) significantly oscillatory for the stimulus that elicited the maximal firing rate (n = 76). (B) Mean oscillation amplitude (peak/cutoff) as a function of normalized firing rate for stimuli yielding significant oscillations. Error bars indicate 2 SD. (C) Histogram comparing the orientation evoking the strongest oscillations with the orientation yielding the maximal firing rate.

 
We also examined whether the amplitude and/or frequency of oscillation was modulated by stimulus orientation. Figure 9BGo plots the strength of oscillation (mean peak/cutoff) as a function of normalized firing rate. It is evident that nonoptimal stimuli elicit weaker oscillations than do optimal or near-optimal stimuli. However, it is also worth noting that even when cells are firing at less than half their maximal firing rates, they may still exhibit significant oscillations. For our entire sample of direction tuning curves, we examined the oscillation frequency for all stimuli that yielded significant oscillations. Oscillation frequency did not vary in any systematic manner as a function of stimulus direction. This finding is consistent with similar results from the anesthetized cat (Gray et al., 1990Go).

Finally, we were interested in how well the tuning based on firing rate matches the tuning based on oscillation strength. As noted above, the amplitude of oscillation for some cells was greater for near-optimal stimuli than for stimuli at the peak of the tuning curve. Figure 9CGo shows a histogram of the orientation difference between the best stimulus in terms of firing rate versus optimal stimulus for oscillation magnitude. The majority of cells show a match between the two measures. However, there are quite a few cells that oscillate more strongly for a stimulus that is 22 or 180° away from optimal. Thus, while oscillatory activity tends to reflect orientation preference in general, the relationship between firing rate and the incidence and strength of oscillation is not simple.

Differences in Oscillation Frequency between Stationary and Moving Stimuli

In our preliminary recordings, we discovered an interesting difference in oscillation frequency between responses to stationary and moving gratings. Because we and others have previously reported that oscillations in area 17 of the anesthetized cat rarely occur in response to stationary stimuli (Eckhorn et al., 1988Go; Engel et al., 1990Go; Gray et al., 1990Go), we chose to pursue this initial finding in more depth. We recorded 84 sessions where presentations of a stationary bar or grating were interleaved with presentations of a bar or grating drifting at the preferred velocity. Figure 10A,BGo shows an example of MUA which exhibited significant oscillations in response to a stationary bar (A), but not a moving bar (B). It is worth noting that the firing rate in response to the moving stimulus was actually greater than that evoked by the stationary stimulus. An example of a single unit that fired rhythmically in response to a moving grating (C) but not to a stationary grating (D) is shown in the second row of Figure 10Go. Although this cell did not fire at a high rate, oscillations are clearly visible in the ACH for the drifting stimulus. Finally, some neurons exhibited oscillatory firing for both stationary (E) and drifting stimuli (F), as illustrated in the third row of Figure 10Go. Even though both correlograms are significant, they have spectral peaks at very different frequencies (i.e. 35 versus 47 Hz).



View larger version (28K):
[in this window]
[in a new window]
 
Figure 10. Autocorrelograms (columns 1 and 3) and corresponding power spectra (columns 2 and 4) computed from responses to stationary and drifting stimuli. (A,B) MUA displaying significant oscillations for static, but not moving stimuli. (C,D) A single V1 unit showing rhythmic firing in response to a drifting grating, but very few spikes in response to a stationary grating. (E,F) Oscillatory MUA recorded in response to both stationary and moving stimuli.

 
Of the entire sample, there were eight sessions where stationary stimuli elicited significant oscillations, but moving stimuli did not; six sessions where moving stimuli led to significant oscillations, but static stimuli did not; and 39 sessions where both types of stimuli evoked significant rhythmicity. For the latter 39 sessions, we found that oscillations evoked by stationary stimuli were lower in frequency (mean = 36 ± 6 Hz) than those evoked by stimuli moving at the preferred velocity (mean = 47 ± 8 Hz) (Wilcoxon Signed Rank, P << 0.001). This trend is illustrated in the upper plot of Figure 11Go, where nearly all the data points fall below the unity slope line. Although stimulus motion consistently led to higher oscillation frequencies, the two types of stimuli evoked equally strong rhythmic responses (P = 0.70). The mean peak/cutoff value for stationary stimuli was 4.33 ± 3.72; for moving stimuli the mean was 4.64 ± 4.10. In the bottom plot of Figure 11Go the data points fall equally above and below the unity slope line. These findings differ from those reported in anesthetized cat (Gray et al., 1990Go), where stationary stimuli often failed to elicit rhythmic activity or evoked oscillations of a smaller amplitude than moving stimuli. Whether this is a genuine difference between cat and monkey striate cells, whether it is due to the effects of anesthesia, or whether it reflects a difference in sampling or the effects of residual eye movements in the alert state remains to be explored.



View larger version (17K):
[in this window]
[in a new window]
 
Figure 11. Scatter plots of oscillation frequency (A) and amplitude (B) for stationary versus moving stimuli. Multiple data points with the same value are occluded. The solid diagonal lines indicate where a 1:1 relationship would fall. Oscillation frequencies (A) clearly fall below this line, indicating that dynamic stimuli elicited higher-frequency oscillations than did stationary stimuli. Oscillation strength (B) was not systematically related to the presence or absence of stimulus motion.

 
Possible Confounding Effects of Flicker Due to Monitor Refresh Rate

It has been suggested that oscillatory neuronal responses may be confounded by the rhythmic activity induced by the flicker of the video monitor employed for visual stimulation (Kiper et al., 1996Go; Mechler et al., 1996Go). Thus, it became important to assess the effect of video monitor refresh on the properties of oscillatory responses. Other studies have documented the existence of oscillatory activity using stimuli generated by a projector with a DC powered light source (Eckhorn et al., 1988Go, 1993Go; Gray and Singer, 1989Go) or by a tungsten light source (Livingstone, 1996Go). We have frequently observed rhythmic firing when using gratings or bars employing a DC light source to roughly assess receptive field position. Figure 12Go compares the response of a single unit stimulated by 20 presentations of a stationary square-wave grating generated by the computer (left) or by 40 presentations of a hand-held DC-generated grating (right). The hand-held stimulus consisted of a black and white slide, taped to a flashlight and placed in the receptive field by the experimenter. A fixation point was illuminated on the computer monitor, but aside from this small dot, the rest of the monitor screen was covered with opaque black paper. Shortly after the fixation point was turned on, the stimulus was presented until the fixation point was turned off. It is evident in Figure 12Go that significant oscillations can result when cells are stimulated by DC-generated stimuli. The frequency of oscillation (i.e. 35 Hz) is identical to the frequency observed when the cell was stimulated by the computer-generated stimulus. We have found results similar to those in Figure 12Go at several other sites.



View larger version (32K):
[in this window]
[in a new window]
 
Figure 12. Comparison of responses to a stationary grating generated by computer versus a hand-held DC light source. PSTHs (A,B), ACHs (C,D) and power spectra (E,F) for a single V1 unit stimulated with a computer-generated square-wave grating (left) and with a grating generated by a DC light source (right).

 

    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Notes
 References
 
Incidence and Frequency of Oscillatory Activity

Our results demonstrate that the properties of gamma-band activity in the striate cortex of the alert macaque monkey are very similar to those observed in the cat (Eckhorn et al., 1988Go; Gray et al., 1990Go; Fries et al., 1997Go; Gray and Viana Di Prisco, 1997Go). Thus, contrary to the view espoused by Young et al. (Young et al., 1992Go) [see also the discussion of Shastri and Ajjanagadde (Shastri and Ajjanagadde, 1993Go)], these data show that gamma-band oscillations are not a physiological curiosity of the cat, but rather are a robust and general feature of neuronal activity in the striate cortex of monkeys.

It is not clear, however, why some research groups have found little or no evidence for gamma-band activity (Krüger and Aiple, 1988Go; Young et al., 1992Go; Nowak et al., 1999Go) or have found oscillatory activity of a different frequency range (Eckhorn et al., 1993Go; Frien et al., 1994Go; Livingstone, 1996Go). One possible confounding factor is the use of anesthesia. Gray and Viana Di Prisco reported that gamma-band activity in cat striate cortex occurs with greater probability and larger amplitude when the animals are alert, suggesting that the state of arousal could affect the strength of oscillatory firing (Gray and Viana Di Prisco, 1997Go). Support for this conjecture comes from the finding that the probability and magnitude of stimulus-induced oscillations are closely correlated with the state of electroencephalogram (EEG) arousal and activation of the mesencephalic reticular formation (Munk et al., 1996Go; Herculano-Houzel et al., 1999Go). Thus, in experiments where the reported incidence of oscillations is low, it is possible that the bulk of the data were collected during states of EEG synchronization [see the discussion of Maldonado et al. (Maldonado et al., 2000Go)].

Variations in neuronal sampling may provide another explanation for the reported differences among studies. Sampling biases can arise in a variety of ways that depend on electrode properties, cellular responsiveness, laminar position and the bias of the experimenter. We attempted to reduce sampling bias by collecting data whenever the signal-to-noise ratio of our recordings was at least 3:1, regardless of the degree of visual responsiveness of the neurons. We also minimized the resampling of previously recorded neurons by restricting our recordings to sites that were at least 200 µm apart. Throughout our study, we found that the occurrence of oscillatory activity frequently varied with the recording depth, and we often encountered regions of strong oscillatory activity extending several hundred microns along a penetration. This apparent laminar specificity is supported by Livingstone's (Livingstone, 1996Go) finding that rhythmic firing in the monkey is most prevalent in superficial layers 2/3 and 4B, and by the study of Gray and McCormick (Gray and McCormick, 1996Go), who reported that rhythmically bursting neurons are localized in layer 2/3 of the cat. These findings suggest that oscillatory activity is more prevalent in the superficial layers of cortex. If this is the case, variations in the incidence of rhythmic activity might stem from disproportionate sampling from the superficial layers. Another form of bias may also have crept into our measurements, however. Because our aim was to assess the properties as well as the prevalence of gamma-band activity, we often spent additional time collecting data when this activity was encountered. Consequently, less recording time was available for sampling neuronal activity at other cortical depths and locations. We cannot exclude the possibility that each of these types of bias may have contributed to our estimates of the incidence and properties of gamma-band activity.

A further reason for variations in the reported incidence of rhythmic activity may stem from fluctuations in the frequency and probability of oscillations across trials. Intertrial frequency variability can dampen the cumulative ACH, even when strong oscillations are present on single trials (Livingstone, 1996Go; Gray and Viana Di Prisco, 1997Go). Livingstone (Livingstone, 1996Go) discussed this problem and chose to classify sessions as oscillatory if 25% of the single trial correlograms were significantly oscillatory. Intertrial variability of oscillatory firing appears to be less of a factor in our data, however. In many instances, we observed strong oscillations that were sustained throughout each trial and varied little in their frequency content across trials (see Figure 2Go). A possible explanation for this apparent difference between our findings and those of Livingstone may again relate to anesthesia or to differences in the stimuli or species of monkey used.

Finally, another potential source of variation between studies may stem from the differences in analysis methods used by different investigators. For example, some studies have relied on simple visual inspection of raw auto- and cross-correlation histograms (Krüger and Aiple, 1988Go). Others have employed curve-fitting procedures to identify rhythmic components in correlation histograms (Engel et al., 1990Go; Young et al., 1992Go; Livingstone, 1996Go), and still others have applied spectral analysis to auto- and cross-correlation histograms (Ghose and Freeman, 1992Go; Gray and Viana Di Prisco, 1997Go), to unit activity directly (Bair et al., 1994Go) or to local field potential recordings (Eckhorn et al., 1988Go, 1993Go; Gray and Singer, 1989Go; Engel et al., 1990Go; Frien et al., 1994Go). These differences in analysis methods are complicated by the fact that the measures used to ascribe statstical significance also vary widely across studies. Given that each method is associated with its own set of pitfalls, it is reasonable to surmise that these factors contribute to the disagreement between studies. In this study, we attempted to maintain consistency with our previous method (Gray and Viana Di Prisco, 1997Go), but at the same time introduced a new measure of statistical significance based on a Monte Carlo simulation of pseudorandom spike trains. While this method does not provide an adequate model to deal with such factors as nonstationarity and burst firing, it does provide a rigorous test of the null hypothesis that neuronal spike trains are Poisson-like in their behavior. Agreement among investigators, however, will probably not be obtained until a uniform set of experimental and analysis methods can be applied to all studies.

Another question raised by our findings is why MUA is more often rhythmic than SUA. One possible explanation is that MUA recordings may contain a mix of oscillatory and nonoscillatory units. If one or more simultaneously recorded units is firing rhythmically and at a high rate, the ACH of the entire group is likely to show significant oscillations, even though some of the units in the group are not oscillatory. This masking of nonoscillatory firing can also explain why rhythmic activity is even more robust for local field potentials than for MUA (Eckhorn et al., 1988Go, 1993Go; Gray and Singer, 1989Go; Frien et al., 1994Go). Another possible explanation for the greater incidence of oscillatory MUA could be that rhythmic activity is a population phenomenon. In this scenario, individual cells may display little or no evidence of oscillatory firing in their ACHs, but are correlated with a population of cells that collectively display rhythmic activity. This type of behavior, which is well documented in the olfactory system, hippocampus, motor cortex and cat visual cortex (Gray, 1994Go; Gray and Viana Di Prisco, 1997Go), can occur when individual cells fire at low rates and the population behavior is variable in frequency. In such instances, it is easy to underestimate the incidence of rhythmic activity when recording solely from single units.

In addition to variations in the incidence of oscillations, some studies have reported higher-frequency oscillations in the monkey compared with the cat or with the data reported here (Eckhorn et al., 1993Go; Frien et al., 1994Go; Livingstone, 1996Go). At present, we are unable to account for these differences, although it is possible that the relation between stimulus velocity and oscillation frequency could provide a partial explanation (Eckhorn et al., 1988Go; Gray et al., 1989Go; Gray and Viana Di Prisco, 1997Go). The results reported here were obtained with drift velocities in the range of 2–6°/s. In a separate study, we have used drift velocities as high as 9°/s and in no instance have we observed significant oscillatory responses that exceeded a frequency of 70 Hz (Maldonado and Gray, 1997Go). Although Livingstone (Livingstone, 1996Go) did not report the stimulus velocities used in her study, Eckhorn et al. (Eckhorn et al., 1993Go) used drift velocities in the range of 1–2°/s. Thus, it appears unlikely that significantly higher stimulus velocities can account for the higher oscillation frequencies observed by these authors compared with the results presented here.

Stimulus Dependence of Oscillatory Activity

The stimulus-dependence of oscillatory activity in V1 has been a further point of contention in the literature. Ghose and Freeman (Ghose and Freeman, 1992Go, 1997Go) have argued that cortical oscillations are not stimulus-dependent and may be generated by oscillatory input occurring spontaneously in the lateral geniculate nucleus (LGN). Several lines of evidence are inconsistent with this hypothesis. First, oscillatory firing in the LGN of the anesthetized cat occurs in a frequency range that only partially overlaps with that observed in the cortex (Laufer and Verzeano, 1967Go; Ghose and Freeman, 1992Go; Ito et al., 1994Go; Neuenschwander and Singer, 1996Go; Castelo-Branco et al., 1998Go). Second, under conditions identical to those employed here, we have recently recorded unit activity in the LGN of one of the monkeys used in this study. While we found numerous instances of rhythmic firing time-locked to the vertical refresh of our video monitor, only 9% (8/86) of our recordings showed evidence of weak stimulus-induced oscillations in the frequency range of 30–100 Hz (Yen and Gray, 1999Go). Third, this study, as well as numerous others, have consistently found that gamma-band activity is not only induced by visual stimulation, but its properties, such as frequency and probability of occurrence, are influenced by the properties of the stimuli [for reviews see (Singer and Gray, 1995Go; Gray, 1999Go; Singer, 1999)]. These findings clearly demonstrate that the occurrence and properties of gamma-band activity in V1 are dependent on visual stimulation.

In this context, we found on average that oscillatory firing occurred most often and with greatest magnitude with those stimuli that evoked the highest firing rates in each cell. However, we found no correlation between the occurrence of gamma-band oscillations and the absolute firing rate across the population of single units. We found many oscillatory cells that fired at rates below 30 Hz in response to their preferred stimuli. In contrast, neither the orientation nor the direction of motion of the stimulus bore a consistent relationship to the frequency of oscillation. Stimulus motion was the only parameter we tested that showed a systematic influence on oscillation frequency. Drifting stimuli consistently evoked higher-frequency oscillations than stationary stimuli, even though oscillation amplitude was the same for drifting and stationary gratings. Together these findings demonstrate that the magnitude and probability of oscillatory activity is directly correlated with the tuning preferences of the cells, and hence their relative firing rates. Oscillation frequency, on the other hand, appears to be dissociated from firing rate.

The finding that stationary stimuli evoke lower-frequency oscillations is interesting, on the one hand, because previous reports have suggested that static stimuli elicit nonrhythmic firing in the cat visual cortex (Eckhorn et al., 1988Go; Gray et al., 1990Go). However, the data in both these reports were taken from a small sample and could represent the subgroup of cells in our study that fire rhythmically only in response to moving gratings. On the other hand, velocity-sensitive oscillations have been documented in anesthetized (Eckhorn et al., 1988Go; Gray et al., 1990Go) and awake cats (Gray and Viana Di Prisco, 1997Go), and in awake monkeys (Eckhorn et al., 1993Go; Maldonado and Gray, 1997Go). In these studies, oscillation frequency was directly proportional to stimulus velocity. Thus, the low-frequency oscillations we observed in response to stationary stimuli may simply represent the lower end-point of the speed spectrum.

Finally, although oscillatory activity is stimulus-dependent, several lines of evidence indicate that, under the conditions of our experiment (i.e. low luminance intensity), it is not time-locked to the stimulus or the refresh rate (80 Hz) of our monitor (but see Mechler et al., 1996). Vigorous oscillatory responses could be evoked by a hand-held DC-generated stimulus, and the distribution of oscillation frequencies in response to computergenerated stimuli was well below the 80 Hz refresh rate of the monitor. Only a small percentage (~5%) of the trial-shuffled correlograms displayed significant spectral peaks. None of these had frequencies at or close to 80 Hz and their magnitudes were smaller than those measured from the unshuffled data. We also found that oscillatory responses displayed variations in frequency and amplitude across repeated trials to the same stimulus. Each of these findings closely parallel similar observations in the cat (Eckhorn et al., 1988Go; Gray et al., 1990Go, 1992Go) and provide strong evidence for a neural induction of oscillatory responses rather than a strictly stimulus-locked process.


    Notes
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Notes
 References
 
We thank Raul Aguilar for software support and Lee Rognlie-Howes for her excellent care of the animals and technical assistance. We also thank Mike Wehr, Rony Azouz, Jean-Philippe Lachaux and Shih-Cheng Yen for their helpful comments on an earlier version of the manuscript. This work was supported by a grant from the National Eye Institute (C.M.G.) and Fellowships from the Klingenstein Foundation (C.M.G.), NIMH (S.F.H.) and McDonnell-Pew Foundation (S.F.H. and P.E.M.).

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
 
1 Present address: Laboratory of Neuropsychology, NIMH, NIH, Bldg 49, Room 1B80, 49 Convent Drive, MSC 4415, Bethesda MD 20892–4415, USA Back

2 Present address: Instituto de Ciencias Biomedicas, Facultad de Medicina, Universidad de Chile, Santiago, Chile Back

3 Present address: Center for Computational Biology, Lewis Hall, Montana State University, Bozeman, MT 59717-3505, USA Back


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Notes
 References
 
Abeles M, Goldstein MH (1977) Multispike train analysis. Proc. IEEE 65:762–773.[ISI]

Bair W, Koch C, Newsome W, Britten K (1994) Power spectrum analysis of bursting cells in area MT in the behaving monkey. J Neurosci 14:2870–2892.[Abstract]

Barlow HB, Blakemore C, Pettigrew JD (1967) The neural mechanism of binocular depth discrimination. J Physiol 193:327–342.[ISI][Medline]

Castelo-Branco M, Neuenschwander S, Singer W (1998) Synchronization of visual responses between cortex, lateral geniculate nucleus, and retina in the anesthetized cat. J Neurosci 18:6395–6410.[Abstract/Free Full Text]

de Oliveira SC, Thiele A, Hoffmann KP (1997) Synchronization of neuronal activity during stimulus expectation in a direction discrimination task. J Neurosci 17:9248–9260.[Abstract/Free Full Text]

Eckhorn R, Bauer R, Jordan W, Brosch M, Kruse W, Munk M, Reitboeck HJ (1988) Coherent oscillations: a mechanism of feature linking in the visual cortex? Biol Cybern 60:121–130.[ISI][Medline]

Eckhorn R, Frien A, Bauer R, Woelbern T, Kehr H (1993) High frequency (60–90 Hz) oscillations in primary visual cortex of awake monkey. NeuroReport 4:243–246.[ISI][Medline]

Engel AK, König P, Gray CM, Singer W (1990) Stimulus-dependent neuronal oscillations in cat visual cortex: intercolumnar interaction as determined by cross-correlation analysis. Eur J Neurosci 2:588–606.[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:6048–6052.[Abstract]

Engel AK, König P, Kreiter AK, Singer W (1991b) Interhemispheric synchronization of oscillatory neuronal responses in cat visual cortex. Science 252:1177–1179.[ISI][Medline]

Freiwald WA, Kreiter AK, Singer W (1995) Stimulus dependent inter-columnar synchronization of single unit responses in cat area 17. NeuroReport 6:2348–2352.[ISI][Medline]

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 S, 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.

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:2273–2277.[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: 12699–12704.[Abstract/Free Full Text]

Gattass R, Gross CG, Sandell JH (1981) Visual topography of V2 in the macaque. J Comp Neurol 201:519–539.[ISI][Medline]

Ghose GM, Freeman RD (1992) Oscillatory discharge in the visual system: does it have a functional role? J Neurophysiol 68:1558–1574.[Abstract/Free Full Text]

Ghose GM, Freeman RD (1997) Intracortical connections are not required for oscillatory activity in the visual cortex. Vis Neurosci 14:963–979.[ISI][Medline]

Gray CM (1994) Synchronous oscillations in neuronal systems: mechanisms and functions. J Comp Neurosci 1:11–38.

Gray CM (1999) The temporal correlation hypothesis of visual feature integration: still alive and well. Neuron 24:31–47.[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:109–113.[Abstract/Free Full Text]

Gray CM, Singer W (1987) Stimulus-specific neuronal oscillations in the cat visual cortex: a cortical functional unit. Soc Neurosci Abstr 13:404.3.

Gray CM, Singer W (1989) Stimulus-specific neuronal oscillations in orientation columns of cat visual cortex. Proc Natl Acad Sci USA 86:1698–1702.[Abstract]

Gray CM, Viana Di Prisco G (1997) Stimulus-dependent neuronal oscillations and local synchronization in striate cortex of the alert cat. J Neurosci 17:3239–3253.[Abstract/Free Full Text]

Gray CM, König P, Engel AK, Singer W (1989) Oscillatory responses in cat visual cortex exhibit intercolumnar synchronization which reflects global stimulus properties. Nature 338:334–337.[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:607–619.[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:337–347.[ISI][Medline]

Gray CM, Friedman-Hill SR, Maldonado PE (1995a) Stimulus-dependent neuronal oscillations in areas V1 and V2 of the alert macaque monkey. Soc Neurosci Abstr 21:592.2.

Gray CM, Maldonado PE, Wilson M, McNaughton B (1995b) Tetrodes markedly improve the reliability and yield of multiple single-unit isolation from multiunit recordings in cat striate cortex. J Neurosci Methods 63:43–54.[ISI][Medline]

Herculano-Houzel S, Munk MHJ, Neuenschwander S, Singer W (1999) Precisely synchronized oscillatory firing patterns require electro-encephalographic activation. J Neurosci 19:3992–4010.[Abstract/Free Full Text]

Hubel DH, Wiesel TN (1977) Uniformity of monkey striate cortex: a parallel relationship between field size, scatter, and magnification factor. J Comp Neurol 158:295–306.

Ito H, Gray CM, Viana Di Prisco G (1994) Can oscillatory activity in the LGN account for the occurrence of synchronous oscillations in the visual cortex? Soc Neurosci Abstr 20:61.7.

Judge SJ, Richmond BJ, Chu FC (1980) Implantation of magnetic search coils for measurement of eye position: an improved method. Vision Res 20:535–538.[ISI][Medline]

Kiper DC, Gegenfurtner KR, Movshon JA (1996) Cortical oscillatory responses do not affect visual segmentation. Vision Res 36:539–544.[ISI][Medline]

König P, Engel AK, Singer W (1995) Relation between oscillatory activity and long-range synchronization in cat visual cortex. Proc Natl Acad Sci USA 92:290–294.[Abstract]

Kreiter AK, Singer W (1992) Oscillatory neuronal responses in the visual cortex of the awake macaque monkey. Eur J Neurosci 4:369–375.[ISI][Medline]

Kreiter AK, Singer W (1996) Stimulus-dependent synchronization of neuronal responses in the visual cortex of the awake macaque monkey. J Neurosci 16:2381–2396.[Abstract]

Krüger J, Aiple F (1988) Multimicroelectrode investigation of monkey striate cortex: spike train correlations in the infragranular layers. J Neurophysiol 60:798–828.[Abstract/Free Full Text]

Laufer M, Verzeano M (1967) Periodic activity in the visual system of the cat. Vision Res 7:215–229.[ISI][Medline]

Livingstone MS (1996) Oscillatory firing and interneuronal correlations in squirrel monkey striate cortex. J Neurophysiol 75:2467–2485.[Abstract/Free Full Text]

Maldonado PE, Gray CM (1996) Heterogeneity in local distributions of orientation selective neurons in the cat primary visual cortex. Vis Neurosci 13:509–516.[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.

Maldonado PE, Friedman-Hill SR, Gray CM (2000) Temporal dynamics of neuronal activity in the striate cortex of alert macaque: II. Short and long-range synchronization. Cereb Cortex 10:1117–1131.[Abstract/Free Full Text]

Mechler F, Shapley R, Hawken MJ, Ringach DL (1996) Video refresh entrains neurons in monkey V1 cortex. Soc Neurosci Abstr 22:672.10.

Milner P (1974) A model for visual shape recognition. Psychol Rev 816:521–535.

Munk MHJ, Roelfsema PR, König P, Engel AK, Singer W (1996) Role of reticular activation in the modulation of intracortical synchronization. Science 272:271–274.[Abstract]

Neuenschwander,S. and Singer,W. (1996) Long-range synchronization of oscillatory light responses in the cat retina and lateral geniculate nucleus. Nature 379:728–733.[ISI][Medline]

Nowak LG, Munk MH, James AC, Girard P, Bullier J (1999) Cross-correlation study of the temporal interactions between areas V1 and V2 of the macaque monkey. J Neurophysiol 81:1057–1074.[Abstract/Free Full Text]

Perkel DH, Gerstein GL, Moore GP (1967) Neuronal spike trains and stochastic point processes. I. The single spike train. Biophys J 7: 391–418.[ISI][Medline]

Shadlen MN, Movshon JA (1999) Synchrony unbound: a critical evaluation of the temporal binding hypothesis. Neuron 24:67–77.[ISI][Medline]

Shastri L, Ajjanagadde V (1993) From simple associations to systematic reasoning. Behav Brain Sci 16:417–494.[ISI]

Singer W (1993) Synchronization of cortical activity and its putative role in information processing and learning. Annu Rev Physiol 55: 349–374.[ISI][Medline]

Singer W, Gray CM (1995) Visual feature integration and the temporal correlation hypothesis. Annu Rev Neurosci 18:555–586.[ISI][Medline]

Tovee MJ, Rolls ET (1992) Oscillatory activity is not evident in the primate temporal visual cortex with static stimuli. NeuroReport 3:369–372.[ISI][Medline]

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:703–710.[ISI]

Yen S-C, Gray CM (1999) Gamma-band neuronal oscillations are largely absent in the lateral geniculate nucleus of the awake-behaving macaque monkey. Soc Neurosci Abstr 25:573.18.

Young MP, Tanaka K, Yamane S (1992) On oscillating neuronal responses in the visual cortex of the monkey. J Neurophysiol 67:1464–1474.[Abstract/Free Full Text]