The Center for Neuroscience, University of California, 1544 Newton Ct, Davis, CA 95616, USA
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
These studies have led to vigorous debates concerning the origin, properties and functional significance of synchronous activity (Gray, 1999; Shadlen and Movshon, 1999
). 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 (2070 Hz). These findings have led to the proposal that gamma-band activity may contribute to the generation of neuronal synchrony (Gray et al., 1989
; König et al., 1995
; Gray and McCormick, 1996
; Livingstone, 1996
). 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, 1992
; Tovee and Rolls, 1992
; Young et al., 1992
; Bair et al., 1994
). 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., 2000
) 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, 1987, 1989
; Eckhorn et al., 1988
; Gray et al., 1990
; Fries et al., 1997
; Gray and Viana Di Prisco, 1997
; Herculano-Houzel et al., 1999
) [but see (Ghose and Freeman, 1992
)], 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, 1992
; Young et al., 1992
; Bair et al., 1994
; de Oliveira et al., 1997
; Nowak et al., 1999
). Other groups have found oscillatory activity, but of a higher frequency (6090 Hz) than in cat (3060 Hz) (Eckhorn et al., 1993
; Frien et al., 1994
). 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, 1992
, 1996
; Livingstone, 1996
).
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., 1992; Livingstone, 1996
) 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., 1993
; Frien et al., 1994
) 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 (2070 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., 2000) have been reported previously in abstract form (Friedman-Hill et al., 1995
; Gray et al., 1995a
).
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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., 1980) and a stainless steel post for head restraint (Gray and Viana Di Prisco, 1997
). 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 5080 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 (26), tungsten microelectrodes (12 M resistance) or custom-fabricated nichrome tetrodes (Gray et al., 1995b
), separated by either 250 µm or 34 mm. The signals were amplified (10K), bandpass filtered (0.6 kHz6kHz), 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, 1977) or by using the stereo properties of tetrodes to cluster spikes by peak-to-peak amplitude, peak time, spike width, etc. (Gray et al., 1995b
). 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.53.0 cycles/°, 812 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 530 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.51 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.250.75 s for a duration of 2.02.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., 1967). 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, 1977; Gattass et al., 1981
). 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, 1997). 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, 1996
). 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., 1967) 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.02.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 2080 Hz (Gray and Viana Di Prisco, 1997). 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, 1997
).
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 1A 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 1A
. The ISIH in Figure 1C
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 1
as chattering cells. We have documented the existence of these cells in extracellular (Gray et al., 1990
; Gray and Viana Di Prisco, 1997
) and intracellular (Gray and McCormick, 1996
) recordings from cats and in extracellular recordings from awake monkeys (Friedman-Hill et al., 1996
). The autocorrelation histogram of this cell exhibits robust oscillations at a frequency of 52 Hz (Figure 1D,E
).
|
|
|
|
|
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 4AD) 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, 1992, 1997
). We addressed this question in two ways. First, we computed ACHs for 200500 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: 2366 Hz) and a median peak/cutoff of 1.36 (range: 1.024.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 2). 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 5
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 5
with Figure 4
), but significant rhythmic activity was still present for single trials. Moreover, the distribution of oscillation frequencies for trial and session ACHs were also similar.
|
|
|
|
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, 1987, 1989
; Eckhorn et al., 1988
; Gray et al., 1990
), 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., 1993
). 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 8 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.
|
|
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 9C 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., 1988; Engel et al., 1990
; Gray et al., 1990
), 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,B
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 10
. 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 10
. Even though both correlograms are significant, they have spectral peaks at very different frequencies (i.e. 35 versus 47 Hz).
|
|
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., 1996; Mechler et al., 1996
). 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., 1988
, 1993
; Gray and Singer, 1989
) or by a tungsten light source (Livingstone, 1996
). We have frequently observed rhythmic firing when using gratings or bars employing a DC light source to roughly assess receptive field position. Figure 12
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 12
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 12
at several other sites.
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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., 1988; Gray et al., 1990
; Fries et al., 1997
; Gray and Viana Di Prisco, 1997
). Thus, contrary to the view espoused by Young et al. (Young et al., 1992
) [see also the discussion of Shastri and Ajjanagadde (Shastri and Ajjanagadde, 1993
)], 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, 1988; Young et al., 1992
; Nowak et al., 1999
) or have found oscillatory activity of a different frequency range (Eckhorn et al., 1993
; Frien et al., 1994
; Livingstone, 1996
). 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, 1997
). 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., 1996
; Herculano-Houzel et al., 1999
). 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., 2000
)].
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, 1996) 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, 1996
), 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, 1996; Gray and Viana Di Prisco, 1997
). Livingstone (Livingstone, 1996
) 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 2
). 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, 1988). Others have employed curve-fitting procedures to identify rhythmic components in correlation histograms (Engel et al., 1990
; Young et al., 1992
; Livingstone, 1996
), and still others have applied spectral analysis to auto- and cross-correlation histograms (Ghose and Freeman, 1992
; Gray and Viana Di Prisco, 1997
), to unit activity directly (Bair et al., 1994
) or to local field potential recordings (Eckhorn et al., 1988
, 1993
; Gray and Singer, 1989
; Engel et al., 1990
; Frien et al., 1994
). 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, 1997
), 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., 1988, 1993
; Gray and Singer, 1989
; Frien et al., 1994
). 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, 1994
; Gray and Viana Di Prisco, 1997
), 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., 1993; Frien et al., 1994
; Livingstone, 1996
). 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., 1988
; Gray et al., 1989
; Gray and Viana Di Prisco, 1997
). The results reported here were obtained with drift velocities in the range of 26°/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, 1997
). Although Livingstone (Livingstone, 1996
) did not report the stimulus velocities used in her study, Eckhorn et al. (Eckhorn et al., 1993
) used drift velocities in the range of 12°/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, 1992, 1997
) 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, 1967
; Ghose and Freeman, 1992
; Ito et al., 1994
; Neuenschwander and Singer, 1996
; Castelo-Branco et al., 1998
). 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 30100 Hz (Yen and Gray, 1999
). 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, 1995
; Gray, 1999
; 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., 1988; Gray et al., 1990
). 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., 1988
; Gray et al., 1990
) and awake cats (Gray and Viana Di Prisco, 1997
), and in awake monkeys (Eckhorn et al., 1993
; Maldonado and Gray, 1997
). 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., 1988; Gray et al., 1990
, 1992
) and provide strong evidence for a neural induction of oscillatory responses rather than a strictly stimulus-locked process.
![]() |
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: Instituto de Ciencias Biomedicas, Facultad de Medicina, Universidad de Chile, Santiago, Chile
3 Present address: Center for Computational Biology, Lewis Hall, Montana State University, Bozeman, MT 59717-3505, USA
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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:28702892.[Abstract]
Barlow HB, Blakemore C, Pettigrew JD (1967) The neural mechanism of binocular depth discrimination. J Physiol 193:327342.[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:63956410.
de Oliveira SC, Thiele A, Hoffmann KP (1997) Synchronization of neuronal activity during stimulus expectation in a direction discrimination task. J Neurosci 17:92489260.
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:121130.[ISI][Medline]
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: intercolumnar 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]
Freiwald WA, Kreiter AK, Singer W (1995) Stimulus dependent inter-columnar synchronization of single unit responses in cat area 17. NeuroReport 6:23482352.[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: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.
Gattass R, Gross CG, Sandell JH (1981) Visual topography of V2 in the macaque. J Comp Neurol 201:519539.[ISI][Medline]
Ghose GM, Freeman RD (1992) Oscillatory discharge in the visual system: does it have a functional role? J Neurophysiol 68:15581574.
Ghose GM, Freeman RD (1997) Intracortical connections are not required for oscillatory activity in the visual cortex. Vis Neurosci 14:963979.[ISI][Medline]
Gray CM (1994) Synchronous oscillations in neuronal systems: mechanisms and functions. J Comp Neurosci 1:1138.
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, 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:16981702.[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:32393253.
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: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, 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:4354.[ISI][Medline]
Herculano-Houzel S, Munk MHJ, Neuenschwander S, Singer W (1999) Precisely synchronized oscillatory firing patterns require electro-encephalographic activation. J Neurosci 19:39924010.
Hubel DH, Wiesel TN (1977) Uniformity of monkey striate cortex: a parallel relationship between field size, scatter, and magnification factor. J Comp Neurol 158:295306.
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:535538.[ISI][Medline]
Kiper DC, Gegenfurtner KR, Movshon JA (1996) Cortical oscillatory responses do not affect visual segmentation. Vision Res 36:539544.[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:290294.[Abstract]
Kreiter AK, Singer W (1992) Oscillatory neuronal responses in the visual cortex of the awake macaque monkey. Eur J Neurosci 4:369375.[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: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.
Laufer M, Verzeano M (1967) Periodic activity in the visual system of the cat. Vision Res 7:215229.[ISI][Medline]
Livingstone MS (1996) Oscillatory firing and interneuronal correlations in squirrel monkey striate cortex. J Neurophysiol 75:24672485.
Maldonado PE, Gray CM (1996) Heterogeneity 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.
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:11171131.
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:521535.
Munk MHJ, Roelfsema PR, König P, Engel AK, Singer W (1996) Role of reticular activation in the modulation of intracortical synchronization. Science 272:271274.[Abstract]
Neuenschwander,S. and Singer,W. (1996) Long-range synchronization of oscillatory light responses in the cat retina and lateral geniculate nucleus. Nature 379:728733.[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:10571074.
Perkel DH, Gerstein GL, Moore GP (1967) Neuronal spike trains and stochastic point processes. I. The single spike train. Biophys J 7: 391418.[ISI][Medline]
Shadlen MN, Movshon JA (1999) Synchrony unbound: a critical evaluation of the temporal binding hypothesis. Neuron 24:6777.[ISI][Medline]
Shastri L, Ajjanagadde V (1993) From simple associations to systematic reasoning. Behav Brain Sci 16:417494.[ISI]
Singer W (1993) Synchronization of cortical activity and its putative role in information processing and learning. Annu Rev Physiol 55: 349374.[ISI][Medline]
Singer W, Gray CM (1995) Visual feature integration and the temporal correlation hypothesis. Annu Rev Neurosci 18:555586.[ISI][Medline]
Tovee MJ, Rolls ET (1992) Oscillatory activity is not evident in the primate temporal visual cortex with static stimuli. NeuroReport 3:369372.[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:703710.[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:14641474.