JN Miami Valley Hospital
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 87: 2715-2725, 2002;
0022-3077/02 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (34)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Brosch, M.
Right arrow Articles by Scheich, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Brosch, M.
Right arrow Articles by Scheich, H.

The Journal of Neurophysiology Vol. 87 No. 6 June 2002, pp. 2715-2725
Copyright ©2002 by the American Physiological Society

Stimulus-Related Gamma Oscillations in Primate Auditory Cortex

Michael Brosch, Eike Budinger, and Henning Scheich

Leibniz-Institut für Neurobiologie, 39118 Magdeburg, Germany


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Brosch, Michael, Eike Budinger, and Henning Scheich. Stimulus-Related Gamma Oscillations in Primate Auditory Cortex. J. Neurophysiol. 87: 2715-2725, 2002. With a multielectrode system, we explored neuronal activity in the gamma  range (>40 Hz) in the primary and caudomedial auditory cortex of six anesthetized macaque monkeys. Stimuli were tone bursts of 100- to 500-ms duration that were presented at sound pressure levels of 40-60 dB and were varied over a wide range of frequencies. These stimuli induced gamma  oscillations, not phase-locked to the onset of stimulation, in 465 of 616 multiunit clusters and at 321 of 422 sites at which field potentials were recorded. Occurrence of gamma  activity was stimulus dependent. It was mostly seen when the stimulus was at the units' preferred frequency. The incidence of gamma  activity decreased with increasing difference between stimulus frequency and preferred frequency. gamma  activity emerged 100-900 ms after stimulus onset with highest incidence ~120 ms. Amplitudes of stimulus-induced gamma  oscillations in field potentials were, on average, almost twice the amplitude of spontaneously occurring gamma  oscillations. gamma  activity at different sites within the primary and the caudomedial auditory field could be synchronized at near-zero phase. Synchrony depended on the spatial distance and on the receptive fields similarity of pairs of units. It decreased with increasing distance between recording sites and increased with similarity of preferred frequencies of the pairs of units. The results indicate that stimulus-induced gamma  oscillations originate from sources in the auditory cortex. They further suggest that gamma  oscillations may provide a mechanism utilized in many parts of the sensory cortex, including the auditory cortex, to integrate neurons according to the similarity of their receptive fields.


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

High-frequency brain rhythms >30 Hz (gamma  oscillations) are thought to be involved in visual, auditory, tactile, and motor processing (for review, see Eckhorn 1999; Freeman 1998; Gray 1999; Singer 1999; Tallon-Baudry and Bertrand 1999). Two types of gamma  oscillations have been distinguished, evoked and induced gamma  oscillations (Galambos 1992). Evoked gamma  oscillations are precisely phase-locked to the onset of sensory stimulation and are present during the initial 100-ms period after stimulus onset. A particular type of evoked gamma  oscillation provides the steady-state response, which is elicited at greatest amplitudes by periodic stimulation at a frequency of ~40 Hz (Galambos et al. 1981). Evoked gamma  oscillations have been implicated with attentional processes (Sheer 1989; Tiitinen et al. 1993), vigilance, and consciousness (Llinas and Ribary 1993; May et al. 1994) and with the temporal binding of successive sensory events (Joliot et al. 1994). Induced gamma  oscillations are characterized by poor temporal locking to stimulus onset and occur ~200-400 ms past stimulus onset. Induced gamma  oscillations are task dependent (Bertrand et al. 1998; Marshall et al. 1996; Sheer 1989) and are related to associative learning (Miltner et al. 1999), sensory/motor integration (Murthy and Fetz 1992; Salenius et al. 1996; Sanes and Donoghue 1993), feature binding (Eckhorn 1999; Eckhorn et al. 1988; Gray 1999; Gray and Singer 1989; Singer 1999), and object representation (Pantev 1995; Tallon-Baudry and Bertrand 1999).

The neural basis of gamma  oscillations has been investigated most extensively in the visual system (Eckhorn 1999; Gray 1999; Singer 1999). In these reviews, it is stated that gamma  oscillations can be found in discharges and local slow wave field potentials recorded in the retina, thalamus, and various areas of the visual cortex in anesthetized and awake cats and monkeys. The occurrence, amplitude, and frequency of gamma  oscillations can be modified by the visual stimulus and by electrical stimulation of the mesencephalic reticular formation and are related to behavioral tasks (Amzica et al. 1997). gamma  oscillations at different parts of the visual system can be synchronous, and their synchrony is related to the similarity of stimulus elements.

Compared with the visual system, as well as to the olfactory system (Freeman 1998), little is known about gamma  oscillations in the auditory cortex. Although there is an increasing number of studies reporting acoustically evoked gamma  oscillations in human electro- and magnetoencephalographic (EEG and MEG) signals, there are relatively few reports that have actually observed gamma  oscillations in the neural activity in the auditory cortex. Except for a recent observation in the human electrocorticogram (Crone et al. 2001), most of our current knowledge is based on epidural recordings in rodents. In a series of experiments in lightly anesthetized rats, Barth and coworkers (Barth and MacDonald 1996; Franowicz and Barth 1995; MacDonald et al. 1996, 1998; Sukov and Barth 2001) have found gamma  oscillations with a frequency of ~40 Hz in the primary (AI) and secondary (AII) auditory cortex. The gamma  oscillations at different sites within the same area were highly synchronized with no phase lag. whereas there was a phase lag of ~2 ms between AI and AII. gamma  oscillation occurred in the ongoing activity but could also be evoked by electrical stimulation of the posterior intralaminar nucleus of thalamus. Moreover, ongoing gamma  oscillations were inhibited by click stimuli but reappeared after the evoked potential at enhanced amplitudes. Recently gamma  oscillations have also been found in slice preparations of rat and mouse auditory cortex, where they could be evoked by electrical stimulation of thalamic afferents (Metherate and Cruikshank 1999).

The goal of the present study is to elucidate properties of gamma  oscillations in the auditory cortex in more detail. This is important because of the increasing interest in the functional role of gamma  oscillations in EEG and MEG signals of humans. Because of their similarity to humans in terms of anatomy and physiology, experiments were carried out in macaque monkeys in which we addressed the following questions: are gamma  oscillations induced by specific acoustic stimuli? What is the relation between gamma  oscillations occurring in slow wave field potentials and in spike trains? Are gamma  oscillations at different cortical sites correlated, and if so, is correlation related to the frequency selectivity of units? Are gamma  oscillations precisely or loosely locked to acoustic stimuli?


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Animal preparation

Data are from six experiments performed on Macaca fascicularis (2.5-4.2 kg). The surgical procedure was initiated by injections of atropine (0.5 mg/kg) and a mixture of ketamine HCl (4 mg/kg) and xylazine (5 mg/kg im). Thereafter we performed a tracheotomy and an extensive craniotomy over the left auditory cortex, which included the lateral inferotemporal sulcus and parts of the central sulcus. In three animals, we aspirated parts of the parietal cortex overlaying the auditory core and belt. Animals received prophylactic injections of an antibiotic (Anicilin; 200 mg · kg-1 · day-1) and dexamethasone (0.2 mg · kg-1 · day-1). During the recordings, animals were continuously given, through an intraperitoneal catheter, ketamine (~30 mg · kg-1 · h-1) and xylazine (~ 25 mg · kg-1 · h-1) in a Ringer/glucose solution (~5 ml · kg-1 · h-1). In one subject (monkey E), xylazine was replaced by diazepam (~ 15 mg · kg-1 · h-1). The depth of anesthesia was checked by periodically monitoring heart beat, inspiration and expiration, body temperature, EEG, and state of eyelid and withdrawal reflexes. Recording sessions lasted 35-85 h. In three animals, 100 nl fluorescent latex beads (0.03-0.5 µm; Microprobe) were injected close to electrophysiologically characterized recording sites to aid localization of recording sites relative to the parvalbumin staining characteristic for AI. Experiments were approved by the authority for animal care and ethics of the federal state Saxony Anhalt (No. 43.2-42502/2-253 IfN) and conformed with the rules for animal experimentation of the European Communities Council Directive (86/609/EEC).

Neural recording

Experiments were conducted in an electrically shielded sound-attenuated double-walled room (IAC, model 1202-A). Simultaneous recordings of neural signals were made with a linear array of seven-fiber microelectrodes, which had an inter electrode separation of 330 µm and in which the electrodes could be moved independently from each other (Thomas Recording). Electrode impedance ranged between 2 and 4 MOmega (measured at 1 kHz). Signals on each electrode were referred to a contact mounted in the frontal bone, amplified, and passed through two parallel filter banks (Thomas Recording), set to yield action potentials (0.5-5 kHz) and intracortical slow wave field potentials filtered at 1 Hz (roll-off 12 dB/octave) and at 140 Hz (roll-off 30 dB/octave). All signals could be monitored on oscilloscopes and on an audio monitor. For quantitative analyses, the signals were linked to a 32-channel A/D data-acquisition system (Brain Wave; DataWave Technologies), which was programmed to detect the action potentials of individual neurons or a few neurons from each electrode. Action potentials were accepted when their amplitude was more than three times above the background signal and when their duration was between 50 and 500 µs. In the following, we did not distinguish between single- and multiunit data and referred to both as "unit." The data-acquisition system also stored the time of the occurrence of action potentials, created peristimulus dot rastergrams from them in real-time, and sampled field potentials at a rate of 1/651 per second.

In animals in which parts of the parietal cortex were aspirated, electrodes penetrated the supratemporal plane almost at a right angle. In the other animals electrodes first entered the parietal cortex at an angle of about 45° before they arrived in the auditory cortex. Most recordings were made 200-600 µm after the first observation of neural discharges in the supratemporal plane and thus were presumably made from upper cortical layers. The first electrode penetrations were from AI. Consecutive penetrations were done from adjacent locations in AI and the caudomedial auditory field (CM) by moving the electrode array along the mediolateral axis in steps of 330 µm or along the rostrocaudal direction in steps of 2.33 mm. Areal locations of recording sites was assessed by inspecting the gradient of best frequency (see following text) as well as the cytoarchitectonic features and the parvalbumin staining of brain sections relative to dye stained recording sites (Kosaki et al. 1997).

Acoustic stimuli

Acoustic search stimuli (tone pips, frequency sweeps, and noise bursts) were produced with a waveform generator (Tucker-Davis Technologies, model WG1). For quantitative analyses, acoustic signals were generated digitally in a computer (Pentium-PC interfaced with an array processor AP2-card, Tucker-Davis Technologies) at a sampling rate of 100 kHz and with a dynamic range of 88 dB and D/A converted to an analog signal (Tucker-Davis Technologies, model DA1). Aliasing was reduced with a low-pass filter set at a cutoff frequency of 35 kHz (Tucker-Davis Technologies, model FT5). Signal amplitude could be controlled over a 100-dB range by passive attenuators (Tucker-Davis Technologies, model PA4). Signals were then passed through an equalizer (SEA 4500; Conrad Electronics) to compensate for the frequency response of the sound system, amplified (Pioneer, model A202), and coupled to a free-field loudspeaker (Manger), which was located on the right frontal side of the animal at a distance of 1 m. The sound pressure level (SPL) was measured with a free-field 0.5-in microphone (GRAS, 40AC) located close to the monkey's ear and a spectrum analyzer (Rion, SA 77). The output of the sound delivering system varied 10 dB in the frequency range of 0.2-35 kHz. At SPL <90 dB, harmonic distortion was >= 36 dB below the signal level.

The frequency tuning of units was measured by presenting 400 single pure tones at 40 frequencies (each repeated 10 times) with a duration of 100 ms in pseudo-random order. Frequencies were equidistantly spaced on a logarithmic scale and covered a range of two to eight octaves, centered on the best frequencies of the neurons recorded simultaneously on the seven electrodes. Best frequency (BF) was defined as the tone that elicited the highest number of action potentials during the tone presentation; it ranged between 0.07 and 31.7 kHz in different units. For each recording, all tones were presented at the same SPL, which varied between 40 and 60 dB. Tone duration was 100 ms (including 5-ms rise and fall time) and intertone interval was 1,375 ms, except for monkey E, where it was 1,078 ms. Often measurements were repeated with tones of longer duration (300 or 500 ms).

Data analysis

Off-line data analyses were carried out with Brain Wave Common Processing version 5.0 (DataWave Technologies) and MATLAB 4.0 (Mathworks). Spectral and temporal properties of gamma  oscillations were analyzed from field potentials. The stimulus-specificity of gamma  oscillations was examined from the spike trains of units.

FIELD POTENTIALS. To assess the effect of pure tone stimulation on the spectral composition of field potentials, we compared different short-time spectra with the spectrum of ongoing activity. This analysis was performed on 70 trials. In 10 of these trials, the tone was at the frequency that evoked the greatest middle-latency potential, which was measured as the amplitude between the first positive and the first negative deflection in this potential. To increase the power of our statistical tests, we included another 60 trials for this analysis in which we presented tones with the six frequencies closest to the tone frequency evoking the largest middle-latency potential. The field potential record of each of the 70 trials was divided into 59 overlapping 64-point (98.3 ms) time windows with a shift interval of 16 points. In monkey E, the trial length of 1,178 ms accommodated 47 time windows. Each time window was tapered with a three-term Blackman-Harris window (Harris 1978) to reduce the high-frequency artifacts induced by the sectioning of the data. A discrete Fourier transform algorithm then yielded the complex spectrum, providing amplitude and phase information of different spectral components of the signal. The frequency resolution of the spectral estimates was 10.17 Hz. The lowest frequency bin ranged from 0 to 10.17 Hz with a center frequency of 5.08 Hz. The zero time bin started 49.15 ms before and ended 49.15 ms after stimulus onset. For an easier reading, values of time and frequency bins are rounded in the remainder of this article.

From the 70 complex spectra of each time window, we next calculated two different average amplitude spectra. In the first spectrum, phase differences between the 70 spectra were neglected. This was done by calculating the absolute value from each spectrum before they were averaged. This procedure, hence, revealed the mean amplitude of different frequency components of the field potential whether or not they had the same phase relation to stimulus onset in each trial. In the second spectrum, phase information of the different spectra was preserved by averaging the complex spectra without taking the absolute value prior to calculating the amplitude spectrum. This second procedure, hence, yielded only those frequency components that were phase-locked to stimulus onset in each trial.

For each 98.3-ms time window, we then assessed whether stimulation had modified the amplitude in individual frequency bins. This was done by employing t-tests separately for all frequency bins in which we compared the distribution of the amplitudes of the 70 trials of a given time window with the corresponding distribution of ongoing activity. The latter was estimated by analyzing the time window immediately before onset of stimulation, which corresponded to the time window 1,361-1,475 ms after onset of stimulation because of the cyclic stimulation. This time window provided a good estimate of ongoing or spontaneous activity because the field potential was no longer influenced by the preceding stimulus, in accordance with previous studies of the auditory-evoked potential (Franowicz and Barth 1995; Ohl et al. 2000). To account for the multiple comparisons performed on the 20 frequency bins between 5.1 and 198.4 Hz, a frequency bin was deemed to be significantly modified by the stimulus if the P value of the t-test was <0.0001.

The similarity of pairs of simultaneously recorded field potentials was assessed by compiling average coherence spectra for three time windows, namely 0-300 ms before and 0-100 and 100-400 ms after stimulus onset. This analysis was performed on the responses to the seven tones closest to the mean of the tones that evoked the greatest amplitude of the middle-latency auditory-evoked potential at the two recordings sites. The signal in each time window was sectioned into overlapping 64-point (98.3 ms) periods with a shift interval of 16 points, tapered with a three-term Blackman-Harris window, and Fourier transformed. Subsequently the cross spectrum was calculated by multiplying the complex spectra of the two field potentials. Last the squared coherence was calculated by dividing the square of the mean of all cross spectra through the average power spectra of the signal pair. With this method, identical signals are found to have a coherence value of 1, whereas statistically independent signals have a coherence of maximally 1 - (1 - alpha ).and [1/(N - 1)] (Rosenberg et al. 1989), where N is the number of time windows and alpha  was set to 0.99. The phase relation between frequency components of two field potential recordings was determined by the inverse tangent of the ratio of the imaginary and real part of the average cross spectrum. This value was divided by the frequency to give the average phase shift in milliseconds of a specific frequency component of two signals.

DISCHARGES. The temporal pattern of spike trains was examined for the period 100-400 ms after stimulus onset, separately for all 40 tone stimuli. The period was the same for all units because most field potential recordings exhibited increased power in the gamma  range during this period as shown later. We calculated the autocorrelation from spike trains that were binned at 2-ms resolution and summed the autocorrelations from the 10 presentations of the same frequency. Then the summed autocorrelation was smoothed by it convoluting with a five-point triangular kernel and Fourier transformed to yield the amplitude spectrum. A unit was considered to have a periodic firing pattern if the amplitude of at least one frequency bin >40 Hz was >4 SDs (P < 0.0001) above the amplitude of the respective frequency bin found in the spectrum of the autocorrelogram of a spike train whose interspike intervals were Poissonian distributed and whose mean rate was equal to that of the spike train under consideration.

To assess the similarity of the temporal patterns of pairs of spike trains, we calculated the average cross-correlation with a bin size of 2 ms from the period 100-400 ms after onset of all repetitions of all stimuli (e.g., Fig. 8). The cross-correlation was convoluted with a five-point triangular kernel and an amplitude spectrum was computed. Those frequency bins of the amplitude spectrum were considered to indicate common periodic firing patterns that significantly (P < 0.0001) exceeded the corresponding frequency bins of the amplitude spectrum of two uncorrelated spike trains. The latter amplitude spectrum was calculated by generating two independent Poissonian spike trains with the same mean rates as the spike trains under consideration.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

gamma oscillations were observed both in field potentials and in spike trains recorded from the auditory cortex of the macaque. Figure 1A shows individual field potential traces from the caudomedial field in response to 10 presentations of a 760-Hz tone burst. Inspection of these traces revealed a period with high-frequency oscillations that emerged after the middle-latency-evoked potential complex and lasted for several hundred milliseconds. Onset, frequency, and amplitude as well as duration of these high-frequency oscillations varied considerably from trial to trial, indicating that these oscillations were not phase-locked to stimulus onset. This was also reflected in the average of the 10 field potential traces (Fig. 1A, bottom), which showed little indication of high-frequency oscillations. The average was rather dominated by a middle-latency-evoked potential, indicating that this component of the field potential was precisely locked to the onset of stimulation.



View larger version (49K):
[in this window]
[in a new window]
 
Fig. 1. Stimulus-induced gamma  oscillations in the auditory cortex of the macaque monkey. A: field potentials in the caudomedial field recorded during and after 10 presentations of a 760-Hz tone burst (presence indicated by gray shading). Note the oscillatory patterns commencing ~100 ms after tone onset. Bottom: the average of the 10 traces shown above. B: dot display (top) and peristimulus time histogram (bottom) of the discharges recorded in parallel with the field potentials shown in A. Dashed line, the rate equal to 4 SDs above the mean spontaneous discharge rate. Note that, unlike many other units in the auditory cortex, this unit did not exhibit a transient response. C: average amplitude spectra of the field potentials shown in A observed in different time windows. The time windows were 0-300 ms before (gray line) and 0-100 ms (thin line) and 100-400 ms after tone onset (black line). Note the stimulus-induced enhancement of field potential frequencies between 40.7 and 71.2 Hz in the latter time window. D: autocorrelogram of the discharges shown in B calculated from the time window 100-400 ms after tone onset. Note the oscillatory structure of the autocorrelogram indicating periodic neuronal firing. The decrease of correlation with increasing delay is due to the rectangular analysis window. E: amplitude spectrum of the autocorrelogram shown in D. The 2 gray lines indicate the mean (bottom) and the mean plus 4 SDs (top) of the amplitude spectrum calculated from the autocorrelogram of a Poissonian spike train with the same mean rate as the 1 shown in B.

Because high-frequency oscillations varied considerably in frequency over time and were not phase-locked to stimulus onset, the spectral composition of field potentials was quantitatively assessed by adding the amplitude spectra of short data segments of the 10 individual trials (see METHODS). The resulting average amplitude spectra showed that, during the period 100-400 ms after stimulus onset, the field potential had considerable power in the range of 41-71 Hz, i.e., in the gamma  range (Fig. 1C). This was in contrast to the initial 100 ms after stimulus onset and to ongoing activity (estimated by analyzing the 300-ms period before stimulus onset), when the amplitude of the spectral components decreased monotonically with increasing frequency.

The period of enhanced gamma  oscillations in the field potential coincided with the late part of the spike response of a unit that was recorded in parallel with the field potential (Fig. 1B). The temporal structure of the spike train was analyzed with the autocorrelation technique for the period 100-400 ms after stimulus onset. The oscillatory structure of the autocorrelogram indicated that the unit had a periodic discharge pattern, which was characterized by a significantly increased number of interspike intervals in the range of 12.5 and 22.2 ms (corresponding to frequencies between 45 and 80 Hz; Fig. 1D). This was shown by calculating the amplitude spectrum of the autocorrelation (Fig. 1E) and comparing this spectrum to the spectrum of the autocorrelation of a spike train with Poisson-distributed interspike intervals and a mean equal to that of the spike train shown in Fig. 1B.

gamma oscillations in field potentials

Figure 2 depicts the average time course of frequency-specific increases of field potentials seen at 422 sites in the auditory fields AI and CM. The time course was determined by dividing the 1,475-ms-long field potential traces into 59 overlapping time windows of 98-ms duration and finding in each time window which spectral components of the field potential were significantly increased over corresponding components of ongoing activity (see METHODS). During ongoing activity, there was little gamma  activity; the spectral power decreased monotonically with increasing frequency (Fig. 2, inset). Our analysis was performed on the responses to the six tones closest to the largest middle-latency-evoked potential and the tone that elicited the largest middle-latency-evoked potential. The gray level of each pixel is proportional to the number of cases with a significant amplitude increase in the frequency bin specified on the ordinate and in the time bin specified on the abscissa.



View larger version (53K):
[in this window]
[in a new window]
 
Fig. 2. Compound time course of stimulus-related increases of field potential amplitude in different frequency bins. Increases are referred to amplitudes of ongoing field potentials. The gray level of each pixel is proportional to the number of cases with a significant amplitude increase in the frequency bin specified on the ordinate and in the time bin specified on the abscissa. Total number of recording sites was 422. Each time window had a duration of 98.3 ms and successive time windows overlapped by 3/4 of their length. The 0 time window was centered on the stimulus onset. The time window at 98.3 ms was centered on the stimulus offset. Inset: the average amplitude spectrum of ongoing field potentials, estimated from 98.3-ms time windows, immediately preceding the onset of stimulation (see METHODS).

This analysis revealed that acoustic stimulation increased field potential amplitudes in the time bins 0-909 ms after stimulus onset. During this period, the spectral distribution of increases of field potential amplitudes varied considerably. Initially, i.e., 49 ms after stimulus onset, most recording sites exhibited increased field potential amplitudes in the frequency bins between 20 and 41 Hz. These frequencies corresponded to the temporal structure of the middle-latency-evoked potential (cf. Fig. 1). At a number of sites, the amplitude increase in the 20- to 41-Hz band was accompanied by increases in other frequency bands, especially in the range >41 Hz. After this period, the spectral distribution of increased field potential amplitudes changed notably. There were markedly fewer sites with a low-frequency increase; rather most sites exhibited increases at frequencies between 41 and 102 Hz. This period with high frequencies started ~100 ms after stimulus onset and lasted maximally 800 ms. In the present report, we therefore considered a cortical site to exhibit stimulus-induced gamma  oscillations if the neural signal contained a frequency component >40 Hz.

Results similar to those shown in Fig. 2 were found for the dominant frequency of the field potential, i.e., the frequency bin whose amplitude was maximally increased over that of ongoing activity (Fig. 3). During the initial 100 ms after stimulus onset, the dominant frequency was mostly <= 40 Hz. Thereafter, dominant frequencies moved to higher values of 60-90 Hz. During the period 98-909 ms after stimulus onset, the amplitude of dominant frequencies in the band >40 Hz was increased by factors of 1.31-7.80 over the corresponding amplitude of ongoing field potentials with a median increase of 1.95. Figure 3 also shows that the time course of increased field potential amplitudes was largely similar to the time course of the response probability of the units that were recorded in parallel with the field potentials. The initial period of low-frequency increase of field potential amplitudes was associated with the transient part of the spike response of neurons, whereas gamma  oscillations coincided with the sustained part of their responses.



View larger version (27K):
[in this window]
[in a new window]
 
Fig. 3. Spectral and temporal characteristics of neural response parameters. A: average time course of response probability of units in comparison to the time course of the dominant frequency in field potentials. The dominant frequency was defined as the frequency bin with the greatest significant increase in power.  and , dominant frequencies more than and <40 Hz, respectively. A unit was considered to respond in a 1-ms time bin if its spike rate in response to pure tone stimulation was significantly larger than the spontaneous rate. B: distribution of the dominant frequency of field potentials in the time window 0-74 ms (, left axis) and 98-909 ms (, right axis) after stimulus onset. Note that the number of cases is larger in the latter time window because the analysis window was longer.

Acoustic stimulation with pure tones of 100-ms duration induced increased field potentials with frequency components >40 Hz in 321 of 422 (76.1%) recording sites in AI and CM in at least one time window >= 98 ms after stimulus onset (Table 1). Results were not much different when tones with a duration of 300 or 500 ms instead of 100 ms were tested; they were similar to those shown in Figs. 2 and 3 except that, at some recording sites, epochs with increased field potential amplitudes occurred for a slightly longer duration. This suggests that gamma  oscillations were induced by the onset or by the steady state part of acoustic stimulation.


                              
View this table:
[in this window]
[in a new window]
 
Table 1. Numbers of stimulus-induced and synchronized gamma  oscillations

Despite some differences in the prevalence, amplitudes, and frequencies of gamma  oscillations between subjects (Table 1), the stimulus dependence of gamma  oscillations and the relation between receptive field properties and gamma  oscillations were similar in all subjects and in AI and CM. In two monkeys in which recordings were made from both auditory fields, there was also little difference between AI and CM with regard to prevalence, amplitudes, and frequencies of gamma  oscillations. Therefore data from the two auditory fields were combined for all group results and figures (Figs. 2, 4-7, and 9) given in this report.

To examine the phase-locking to stimulus onset of different frequency components, we analyzed the spectrotemporal composition of the average field potential complex. Like in the previous analysis, we used the 70 trials in which we presented the six tones closest to the greatest middle-latency-evoked potential and the tone that elicited the greatest middle-latency-evoked potential (Fig. 4). This analysis revealed that there were phase-locked field potential components up to a period of maximally 246 ms after stimulus onset. Most of them were in the frequency range between 20 and 41 Hz, which corresponded to the temporal structure of the middle-latency-evoked potential (Fig. 1C shows the evoked potential at a single recording site). This indicated that most of the gamma  oscillations that were observed during the period 98-909 ms after stimulus onset were not phase-locked to stimulus onset. Thus they can be identified as induced gamma  oscillations (see INTRODUCTION).



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 4. Compound time course of field potential components in different frequency bins, which were phase-locked to stimulus onset and increased in amplitude by acoustic stimulation. For presentational purposes, the gray level of each pixel is proportional to the square root of the number of cases. Same conventions as in Fig. 2.

gamma oscillations in neuronal discharges

The temporal pattern of spike trains was assessed with the autocorrelation technique for 616 units. This analysis revealed that 465 units (75.5%) exhibited periodic discharge patterns with an excessive number of interspike intervals in the gamma  range (>40 Hz) during the time window 100-400 ms after stimulus onset for >= 1 of the 40 tones tested on each unit (Table 1). This time window was chosen because the prevalence of gamma  oscillations was highest in field potential recordings during this period. The oscillatory modulation of spike trains was weaker than in field potentials. The probability of observing gamma  activity was associated with the relation of the stimulus frequency to the spectral receptive field of a unit. Stimuli inside the spectral receptive field were more effective in evoking gamma  activity in a unit than were stimuli outside the spectral receptive field (Fig. 5). For stimuli inside the receptive field, the probability of inducing gamma  activity increased the closer a tone was to the BF of a unit and became maximal when the stimulus was at the BF.



View larger version (21K):
[in this window]
[in a new window]
 
Fig. 5. Stimulus-dependence of the occurrence of spike trains with gamma  oscillations (>40 Hz) in auditory cortex.  and , data for tones inside and outside the spectral receptive field, respectively. In this plot, the distribution of best frequency (BF) distance was divided into 20 bins of different widths but of almost equal number of cases. In each bin, the number of units with an oscillatory discharge patterns was divided by the total number of units within this bin to yield the oscillation probability. Analysis window was 100-400 ms after tone onset. Only units were included whose response to a BF tone was >= 2 times larger than the spontaneous rate. Total number of units was 494.

Correlation of gamma  oscillations at different cortical sites

The linear seven-electrode array with an interelectrode spacing of 330 µm allowed us to analyze the synchronization of neural signals that were simultaneously recorded at different sites within AI or CM. Because we were interested in the synchronization in specific frequency bands and time periods, we used short-term coherence spectra instead of the cross-correlation. This analysis was performed on the responses to the seven tones closest to the mean of the tones that evoked the greatest amplitude of the middle-latency-evoked potential at the two recordings sites (see METHODS).

Figure 6A depicts median coherence spectra for the time window 100-400 ms after stimulus onset for 215 recordings pairs in which amplitudes of gamma  oscillations in field potentials were increased by a factor >= 2. In this time window, most gamma  oscillations were observed, as shown in Fig. 3. The coherence spectra demonstrate that gamma  oscillations at different sites of AI or CM could be synchronized. Synchronization was strongest for cortical sites separated by 330 µm and decreased monotonically with increasing separation of sites (Fig. 6B). A linear regression analysis of the frequency bin centered at 45.8 Hz indicated that gamma  oscillations at different sites could be synchronous up to a distance of 2.9 mm, at which coherence attained nonsignificant values (not shown). The synchronization of field potentials, however, was not restricted to the gamma  range but also occurred at lower frequencies. Generally, coherence of field potentials was highest at the lowest frequency and decreased monotonically with increasing frequency to nonsignificant values between 76 and 158 Hz, depending on the spatial separation of recording pairs. When the gamma  oscillations at two sites were correlated, they had, on average, zero phase lag [mean = 0.3 ± 1.4 (SD) ms; Fig. 6C]. This was found by computing the phase relation of the 215 pairs at the frequency >40 Hz at which the maximal coherence was found.



View larger version (18K):
[in this window]
[in a new window]
 
Fig. 6. Synchrony of field potentials in auditory cortex. A: median coherence spectra for different spatial separations between recording pairs (330-2,000 µm, number of cases given in parenthesis), computed from the time window 100-400 ms after tone onset. - - -, 4 SDs above the mean chance coherence. Data were from 215 recording pairs in which the power in a frequency bin >40 Hz of both field potentials was increased >= 2 over the corresponding power of ongoing activity. B: spatial decrease of coherence of field potentials in different frequency bins. C: temporal relation between significantly correlated dominant frequencies in the gamma  range of pairs of simultaneously recorded field potentials. Mean delay was 0.3 ± 1.4 ms.

During the period 100-400 ms after stimulus onset, coherence of field potentials was selectively enhanced in the gamma  range over spontaneous coherence. Figure 7 compares the coherence spectrum 100-400 ms after stimulus onset with the spontaneous coherence spectrum (calculated from the period 0-300 ms before stimulus onset), pooled over all 215 recording pairs shown in Fig. 6A. During the time window of 100-400 ms after stimulus onset, significant (P < 0.0001) increases in coherence over values of ongoing activity were restricted to frequency bins between 56 and 87 Hz, as revealed by performing separate t-tests (P < 0.0001) for individual frequency bins between the distribution of coherence spectra of the two time windows. Enhanced coherence of field potentials was also found in the time window shortly after stimulus onset (0-100 ms). In contrast to the 100- to 400-ms period, however, coherence in the early time window was also increased over spontaneous coherence in frequency bins below the gamma  range.



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 7. Stimulus-induced changes of synchrony of field potentials in auditory cortex. A: average coherence spectra 100-400 ms after () and 0-300 ms before (open circle ) tone onset, pooled over all spatial separations of the recording sites shown in Fig. 6A. Bars give standard error means. B: ratio of the coherence spectra 100-400 ms after and 0-300 ms before tone onset shown in A. , significant differences of coherence values in the 2 time windows (P < 0.0001).

The periodic discharge of units simultaneously recorded at different sites of AI or CM could be synchronized. This is illustrated in Fig. 8, which shows a typical cross-correlation of the discharges of two units in CM during the time window 100-400 ms after stimulus onset. The cross-correlogram had an oscillatory structure with a central and several satellite peaks, which suggested that the periodic discharges of the two units were synchronized. This was verified by calculating the spectrum of the cross-correlation and by finding that the spectral peak in the gamma  range was significantly above the corresponding spectral peak of chance correlation.



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 8. Typical example of the correlation of 2 spike trains recorded simultaneously from the caudomedial auditory field. Separation of recording sites was 1,330 µm. A: cross-correlation calculated from a 300-ms long period commencing immediately after the cessation of a 600-Hz tone of 100-ms duration. B: amplitude spectrum of the cross-correlogram shown in A. Dashed line, the amplitude spectrum of the chance correlation of 2 spike trains, each with the same rate as the spike trains under investigation but with independent Poissonian interspike interval distributions. The thin line represents 4 SDs above the chance amplitude spectrum.

The synchronization of the spike trains of pairs of units could be examined for 1,803 recording pairs from the response to a tone equal to the mean BF of the two units under investigation. Of them 417 (23.1%) exhibited significantly synchronized gamma  activity (>40 Hz) during the time window 100-400 ms after stimulus onset. Generally, gamma  synchronization was related to the similarity of the receptive fields of two units. The probability of two units to discharge synchronously was highest for pairs with the same BF and decreased with increasing difference of their BFs (chi 2-test, P < 0.001; Fig. 9). This relation was found for units with small and large cortical separations. Results were similar when synchrony was analyzed for tones of different frequency.



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 9. Relation between the ratio of BFs of pairs of units and correlation of gamma  oscillations in auditory cortex. The distribution of BF ratio was divided into 8 classes such that all classes, except the 1st class, consisted of the same number of units (n = 201; 1st class: n = 246). The ordinate shows the percentage of units with significantly correlated spike trains in the gamma  range. Correlation between units was analyzed for the tone whose frequency was equal to the mean BF of 2 units during the time window 100-400 ms after stimulus onset.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Acoustic stimulation with pure tones could augment amplitudes of field potentials in auditory cortex for a period of <= 900 ms. In this period, increases of amplitudes of field potentials occurred in a frequency-specific manner. During the initial 100 ms, low-frequency components between 20 and 40 Hz were increased at most recordings sites. These components were mostly phase-locked to stimulus onset and corresponded to the temporal structure of the stimulus-evoked potential complex. At some sites, the low-frequency components were associated with increased frequency components >40 Hz that sometimes were phase-locked to the stimulus. The early increase of power in the gamma  band by acoustic stimulation could be interpreted as evidence for the existence of stimulus-evoked gamma  oscillations in auditory cortex. However, it could also be due merely to a broadband increase of stimulus-evoked field potential power. The early response period differed from the late response period in which most sites exhibited increases of high-frequency components between 40 and 100 Hz that rarely were associated with increases in lower-frequency bands. Therefore the high-frequency oscillations commencing ~100 ms after stimulus onset were probably not an artifact of a broadband increase of field potential power but rather because of their latency and their lack of phase-locking to stimulus onset, represent gamma  oscillations induced by sensory stimulation.

The stimulus-induced gamma  oscillations seen in the present study resemble in many respects gamma  oscillations observed previously in different types of neural signals recorded in auditory cortex of rodents (Barth and MacDonald 1996; Franowicz and Barth 1995; MacDonald et al. 1996,1998; Metherate and Cruikshank 1999; Sukov and Barth 2001). Nevertheless there appear to be some differences that might be accounted for, in part, by differences between preparation, acoustic stimulation, or species. The latency of gamma  oscillations in rat auditory cortex in response to clicks was 300-350 ms and, thus longer than the latency in monkey auditory cortex in response to pure tones. The frequency of induced gamma  oscillations was higher in the auditory cortex of monkeys than they were in rodents, where they ranged between 20 and 80 Hz. In rat auditory cortex, spontaneous oscillations, as well as oscillations induced by acoustic click stimuli or by electric stimulation of unspecific thalamic nuclei were observed over all regions of AI and AII with amplitudes being largest at the border of AI and AII. In monkey auditory cortex, gamma  oscillations could be induced by pure tone stimuli at many sites of AI and CM. In contrast to rat auditory cortex, the present study found no indication that amplitudes of gamma  oscillations were largest at any site of the fields AI and CM in monkeys.

The frequency and stimulus dependence of gamma  oscillations seen in this and other animal studies (Barth and MacDonald 1996; Franowicz and Barth 1995; MacDonald et al. 1996, 1998; Metherate and Cruikshank 1999) are comparable to gamma  oscillations observed in humans after acoustic stimulation (Bertrand et al. 1998; Crone et al. 2001; Joliot et al. 1994; Knief et al. 2000; Llinas and Ribary 1993; Marshall et al. 1996; May et al. 1994; Pantev 1995; Pantev et al. 1991; Tiitinen et al. 1993). The findings in animals that gamma  oscillations can be induced over wide regions of auditory cortex and that gamma  oscillations at different sites can be synchronized suggest that the signals of a large number of local oscillators can superimpose to yield a macroscopic signal of a size that can be measured outside the skull. This corroborates recent findings in the human electrocorticogram (Crone et al. 2001) that evoked and induced gamma  signals seen in human EEG and MEG recordings originate, at least in part, from neural activity emanating in the auditory cortex.

A prerequisite for the observation of oscillatory neural activity appears to be the existence of late sustained responses to sensory stimulation (Fig. 3). This is corroborated by experiments in a slice preparation of the auditory forebrain (Metherate and Cruikshank 1999) that demonstrated that gamma  oscillations never occurred in the absence of a slow potential, evoked by electric stimulation of the thalamic radiation. The absence of late responses might be a reason why acoustically induced firing patterns in the gamma  range have not been observed in most previous studies of auditory cortex. The emergence of late responses and, thus gamma  oscillations, appears to depend, at least partly, on the sensory stimulus. After stimulation with acoustic clicks, power of field potentials in the gamma  band was suppressed for a period of ~300-350 ms that was followed by a rebound of gamma  power (Franowicz and Barth 1995). Pure tones, by contrast, elicited gamma  activity at markedly shorter latencies (Fig. 3). The proper selection of the acoustic stimuli for the induction of gamma  oscillations is further stressed by the finding that, in visual cortex, gamma  oscillations occurred at highest amplitudes in response to smooth visual stimuli, whereas transient-rich visual stimulation evoked fewer gamma  oscillations (Kruse and Eckhorn 1996; but see Friedman-Hill et al. 2000).

Generally, there are striking similarities between induced gamma  oscillations observed in auditory cortex (Barth and MacDonald 1996; Franowicz and Barth 1995; MacDonald et al. 1996, 1998; Metherate and Cruikshank 1999; Sukov and Barth 2001; present study) and in visual cortex (Eckhorn 1999; Gray 1999; Singer 1999). This includes the frequency range, the temporal relation to the evoked potential, their lack of phase-locking to stimulus onset, the prevalence, the stimulus selectivity, and the finding that oscillatory activity recorded within the same or in different cortical fields could be synchronized with each other at near-zero phase lag. This suggests that different parts of the sensory cortex utilize similar discharge patterns to code stimulus information.

Most of the differences between visual and auditory cortex were seen with respect to amplitudes and synchrony of gamma  oscillations. In monkey auditory cortex, acoustic stimulation increased amplitudes of stimulus-induced gamma  oscillations in field potentials, on average, by a factor of 1.95 over respective amplitudes of ongoing activity. The increase was smaller than in the visual cortex of the cat (2.2) (Bauer et al. 1995) and of the monkey (2.8) (Eckhorn 1999). Figure 7 shows that in auditory cortex, synchrony of gamma  oscillations was restricted to a range <3 mm. This range is larger than one would expect from passive volume conduction for which the decay constant in the gamma  range was estimated to be maximally 200 µm in visual cortex (Engel et al. 1990). Nevertheless, the range is smaller than in the visual cortex of cats where coherence of gamma  oscillations extended over wider regions (Brosch et al. 1995; Engel et al. 1990). It was still ~0.5 at a separation of 2 mm and synchronous oscillations were found up to separations of >= 6 mm.

The differences with regard to amplitudes and synchrony of gamma  oscillations in field potentials might simply result from differences in the effect of sensory stimulation. In most of the studies in visual cortex, amplitudes and synchrony of gamma  oscillations were measured with stimuli extending over wide regions of the visual field and, hence, excited a large fraction of visual neurons (Eckhorn 1999; Gray 1999; Singer 1999). Reduction of stimulus size decreased amplitudes of induced gamma  oscillations considerably (Bauer et al. 1995). Therefore it is possible that gamma  oscillations with larger amplitudes and at more sites of auditory cortex can be observed with acoustic stimuli more complex than pure tones, which excite only a limited part of the auditory system (Crone et al. 2001). In addition it could be that stationary stimuli evoke more gamma  oscillations than transient-rich stimuli (Kruse and Eckhorn 1996; but see Friedman-Hill et al. 2000). The differences in oscillation amplitudes and synchrony in the auditory and visual cortex could also arise from dissimilarities in animal preparation or from physiological dissimilarities of the cortical tissue serving the visual and auditory modality. In visual cortex, oscillation amplitudes were considerably greater in awake monkeys (Eckhorn 1999) than in anesthetized cats (Bauer et al. 1995). In the current study, we used ketamine, whereas studies in visual cortex preferentially used gas anesthesia. Although ketamine is known to alter neurotransmission (e.g., Zurita et al. 1994), the effects on gamma  oscillations seem to be different in auditory and visual cortex. In cat visual cortex, ketamine has been reported to induce spontaneous gamma  oscillations at high amplitudes (R. Eckhorn, personal communication). In auditory cortex, neither the present study in monkeys nor previous studies in cats (e.g., Brosch and Schreiner 1999; Eggermont 1994) has found bursts of spontaneous gamma  oscillations with ketamine.

The occurrence of gamma  oscillations was related to the frequency specificity of neurons in auditory cortex. gamma  activity was most often found when neurons were stimulated with a BF tone. Comparable findings have recently been reported for the orientation tuning of neurons in visual cortex (Friedman-Hill et al. 2000; Frien et al. 2000). This indicates that gamma  activity is not only a simple arousal effect (Barth and MacDonald 1996) and occurs whenever neurons are activated but that it is related to the specific excitation of neurons. To our knowledge, the present finding provides the first evidence that, in addition to conventional rate and latency measures, high-frequency tones are represented by the temporal patterning of neural discharges in auditory cortex. It is interesting that the frequency of gamma  oscillations is above the highest frequency of modulated sounds to which most cortical neurons can phase-lock their response (Langner 1992). This offers additional coding capacity to cortical cells during the late part of the stimulus-evoked response when discharge rates are low.

It has been hypothesized that synchronized gamma  activity might be used for the definition of functional neuronal assemblies (Eckhorn 1999; Gray 1999; Singer 1999). Although our study was performed on anesthetized animals, present results are compatible with this hypothesis. Like in the visual cortex, the assembly formation in the auditory cortex by synchronizing the firing of neurons seems to be governed, in part, by the similarity of the receptive fields of neurons. Such assemblies might contribute to the cortical representation of auditory objects. This speculation does not exclude that assemblies might also be formed according to rules other than similarity of receptive fields. It neither excludes that assemblies might also be established by synchronization of stochastic activity or activity in other frequency ranges (von der Malsburg 1983; Reitböck 1983). In the present study, there was also an increase of synchrony of neural signals during the initial 100 ms after stimulus onset, which was most pronounced in the low-frequency range, consistent with recent findings of correlated activity in auditory cortex (Brosch and Schreiner 1999). These results suggest that different neuronal assemblies, dynamically formed by synchronization of neural discharges, contribute to the processing of acoustic signals in auditory cortex.


    ACKNOWLEDGMENTS

The authors thank A. Schulz for participating in the experiments, C. Bucks for help in data analysis, and Dr. Peter Heil for valuable suggestions on the manuscript.


    FOOTNOTES

Address for reprint requests: M. Brosch, Leibniz-Institut für Neurobiologie, Brenneckestrabeta e 6, 39118 Magdeburg, Germany (E-mail: brosch{at}ifn-magdeburg.de).

Received 16 July 2001; accepted in final form 23 January 2002.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES