QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

HOME  |   SEARCH  |   ARCHIVE  |   SUBSCRIBE  |   CONTACT  |   HELP

This Article Abstract Full Text (PDF) Submit an eLetter Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (26) Citing Articles via Google Scholar Google Scholar Articles by Katz, D. B. Articles by Nicolelis, M. A. L. Search for Related Content PubMed PubMed Citation Articles by Katz, D. B. Articles by Nicolelis, M. A. L.

 Previous Article  |  Next Article 

The Journal of Neuroscience, March 1, 2002, 22(5):1850-1857

Taste-Specific Neuronal Ensembles in the Gustatory Cortex of Awake Rats

Donald B. Katz¹, S. A. Simon^{1, 2}, and Miguel A. L. Nicolelis^{1, 2}

Departments of ¹ Neurobiology and ² Biomedical Engineering, Duke University, Durham, North Carolina 27710

	ABSTRACT

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

In gustatory cortex, single-neuron activity reflects the multimodal processing of taste stimuli. Little is known, however, about the interactions between gustatory cortical (GC) neurons during tastant processing. Here, these interactions were characterized. It was found that 36% (85 of 237) of neuron pairs, including many (61%) in which one or both single units were not taste specific, produced significant cross-correlations (CCs) to a subset of tastants across a hundreds of milliseconds timescale. Significant CCs arose from the coupling between the firing rates of neurons as those rates changed through time. Such coupling significantly increased the amount of tastant-specific information contained in ensembles. These data suggest that taste-specific GC assemblies may transiently form and coevolve on a behaviorally appropriate timescale, contributing to rats' ability to discriminate tastants.

Key words: insular; taste; cross-correlation; multiple electrode; coding; dynamics

	INTRODUCTION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

In awake rats, gustatory cortical (GC) neurons produce time-varying responses to gustatory stimuli (Katz et al., 2001a). Such time-varying responses may arise because these neurons are embedded in neural networks that process tastants interactively, as has been suggested for other sensory systems (Bazhenov et al., 2001; Ghazanfar et al., 2001). In recordings from pairs of neurons in anesthetized rats, the occurrence of tight (±20 msec) cross-correlations (CCs) has revealed functional interactions between GC (Yokota et al., 1996, 1997; Nakamura and Ogawa, 1997) or gustatory brainstem (Adachi et al., 1989; Di Lorenzo and Monroe, 1997) neurons. These previous studies, however, did not determine whether neural interactions across hundreds of milliseconds, the timescale of both single-neuron GC response dynamics (Katz et al., 2001a) and gustatory decision making (Halpern, 1985), might be related to system processing (Nowak and Bullier, 2000; Salinas and Sejnowski, 2001).

To determine whether time-varying GC responses are related to network processing, it is necessary to characterize GC neuronal interactions across this broader timescale by recording the simultaneous activity of small ensembles of GC neurons while delivering tastants to the tongues of awake rats. After application of a subset of tastants, pairs of neurons produced peaked CCs (mean half-width 300 msec), even when single-unit analyses suggested that one or both neurons were not taste specific. Taste-specific CCs were explicable as coupling in response rates of neurons as they changed through poststimulus time and suggest the emergence of taste-specific neural assemblies that interact on a behaviorally relevant timescale. The formation of such assemblies facilitates tastant processing in a manner consistent with population theories of taste perception (Erickson, 2001).

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Subjects. All procedures accorded with the National Institutes of Health guidelines for the treatment of animal subjects. Male (n = 3) and female (n = 9) Long-Evans rats (275-300 gm) served as subjects for this study. Because no sex-related differences were observed, the data from male and female rats were pooled.

Surgery. Anesthetized rats [5% halothane and then either pentobarbital (50 mg/kg, i.p.) or ketamine and xylazine (100 and 10 mg/kg, i.m., respectively)] were implanted unilaterally or bilaterally with electrode bundles in the GC [anteroposterior, 1.2-1.5 mm; mediolateral, 5.2 mm; dorsoventral, approximately 4.5 mm from dura] (Kosar et al., 1986). Intraoral cannulas (IOCs) were implanted bilaterally as well (Phillips and Norgren, 1970; Katz et al., 2001a).

Each electrode bundle included 16 Formvar-coated nichrome wires (diameter of 25-µm), cut flat with carbide-tipped scissors. The impedance of the wires was 200-500 K at 1 kHz. The wires were glued to a small microdrive such that they could be advanced through the brain in the weeks after surgery (Katz et al., 2001b).

Behavioral procedure. After recovery from surgery, water-restricted rats were adapted to mild restraint, after which 90-120 min sessions were run in which 40 µl of a randomly selected tastant was delivered every 90 sec, interspersed with 80 µl water rinses. The stimuli were delivered via IOCs or a nozzle situated directly in front of the mouth. The two delivery methods produced similar responses (Nishijo and Norgren, 1991; Katz et al., 2001a). The tastants were citric acid [CA (0.02 M)], NaCl [Na (0.1 M)], sucrose [Suc (0.1 M)], quinine HCl [Q (0.001 M)], nicotine [Nic (100 µM)], and water (W) (Katz et al., 2001a). One rat was not tested with nicotine.

The triggering of tastant delivery solenoids was transmitted to the data acquisition computer (Plexon Inc., Dallas, TX). Off-line, the data were adjusted for physical delays between this signal and the time at which fluid hit the tongue, in accordance with the following test. The stimulus delivery apparatus was placed an appropriate distance from an "artificial rat tongue" (two bare wire tips, separated by 1 mm of air, that were the ends of an open circuit including a battery and oscilloscope); the second input to the oscilloscope monitored the transistor-transistor logic pulse to the fluid-delivery solenoid. NaCl was delivered, connecting and completing the battery circuit, and the resultant delay between solenoid opening and stimulus hitting the tongue could be viewed on the oscilloscope. In the case of delivery via an intraoral cannula, this delay was reliably 45 msec (±4 msec range). In the case of nozzle delivery (which had to be placed further from the tongue), the delay was 90 msec (±8 msec range).

Electrophysiology. Single neurons of >3:1 signal-to-noise ratio were isolated using a waveform template algorithm (Nicolelis et al., 1997a) and corrected off-line using cluster cutting software (Plexon Inc.). All isolations lasted >2 hr, allowing for the delivery of 10-30 trials of each tastant. Figure 1 presents an example of a raw oscilloscope trace, as well as an example of extracted waveforms and associated interspike interval statistics of units isolated from individual microwires. The larger unit in Figure 1A is Neuron 1 from Figure 2A, and the two units in Figure 1, B and C, are Neurons 1 and 2 from Figure 2B.

View larger version (31K): [in this window] [in a new window]   Figure 1.   Examples of GC recordings. A, An oscilloscope trace showing the signal from a single microwire placed in GC. One larger and one smaller spike can be seen. The larger spike is Neuron 1 from Figure 2A. The trace has a time base of 5 msec per division. B, Extracted waveforms from two channels in a different animal. The signal-to-noise ratio of the left neuron was 6:1; that of the right neuron was 4:1. These are Neurons 1 and 2, respectively, used in Figure 2B. C, The interspike interval plots for the neurons in B. Note the presence of a refractory period in the records of both neurons. k = 1000.

View larger version (54K): [in this window] [in a new window]   Figure 2.   Taste-specific CCs between pairs of GC neurons. A, The CCs (left) and PSTHs (right) of a pair of GC neurons to the battery of tastants. The CCs of this pair to NaCl (mauve) is highly significant (significant portions are shown as bold lines); the half-height of the peak is ~300 msec. The abscissa is lag (in seconds); the ordinate is correlation. The PSTHs for Neuron 1 show an positive response to NaCl and subtle negative responses to quinine HCl and acid; Neuron 2 did not respond in a taste-specific manner. B, Similar CCs and PSTHs for another pair of neurons, showing significant negative interactions to nicotine and acid; only mild specificity of response can be seen in the PSTHs. C, CCs and PSTHs for two neurons recorded from a single wire. These neurons cohered in the presence of NaCl and quinine on a very short timescale, visible in the sharp peak just offset from 0 lag, as well as on the longer timescale. Neither neuron showed taste specificity of response in PSTHs. D, CCs and PSTHs for a between-hemisphere pair of neurons (same rat as in Fig. 1B), one of which was mildly taste specific according to PSTH, that produced significant positive CC peaks to nicotine, NaCl, and quinine and showed a biphasic negative-positive interaction to acid. Green, CA; mauve, Na; orange, Nic; dark blue, Q; black, Suc; light blue, W.

Paired-unit analyses. CCs were calculated using spikes produced during the 4 sec after tastant administration. Analyses were performed to ascertain which of various types of coupling were responsible for significant CC peaks and to determine the nature of GC interactions (Perkel et al., 1967; Aertsen et al., 1989; Brody, 1999a). These analyses involved sequentially removing various aspects of the observed responses (see below).

To remove the portions of CCs attributable to common input to both neurons, the "optimal shift predictors" [cross-correlations made after the trials of one of the neurons are shuffled; essentially cross-correlations between the two peristimulus histograms (PSTH) with trial-to-trial coherence ignored] were subtracted from CCs made from the original data (Perkel et al., 1967). The "shift predictor" contains only information that is time locked to stimulus presentation, that is, that pertains to stimulus-driven activation. This subtraction procedure therefore leaves only cross-correlation structure that is related to actual interdependence between the neurons. In fact, this procedure leads to an underestimation of the actual interdependence between neurons, because some true between-neuron interactions may in fact be stimulus locked.

To test whether direct spike-to-spike coupling (synchrony) was responsible for the CC peaks, raw data were compared with bootstrapped simulations (Seidemann et al., 1996; Baker and Lemon, 2000; Bair et al., 2001). First, single-trial spike trains were turned into firing probability vectors, determined at each time point by the interspike interval between previous and succeeding action potentials. These vectors were then smoothed through convolution with a  of 50 msec Gaussian distribution. Once the smoothed vectors (one per trial) were created, spike trains were simulated from each. Ten such "smoothed simulations" were produced, and the average CC for each simulated neuron pair was calculated. Neuronal activity was modeled as gamma processes, which is to say as Poisson processes with refractory periods and "fat tails." The orders of the gamma distributions were chosen to match the interspike interval plots (minimizing least-squared error) of the real data (Baker and Lemon, 2000).

Surrogate spike data were generated using the millisecond-by-millisecond firing probability estimation and the appropriate gamma distribution (Baker and Lemon, 2000). These simulations provided estimates of the expected variability in cross-correlogram height, taking into account the observed single-unit spiking distributions and broad single-trial response patterns but ignoring spike-to-spike interactions in the real data. CC peaks generated from spike trains simulated in this way solely reflected rate-related (i.e., not spike-to-spike) coupling.

If the firing rates of a neuron pair vary together between trials across a session, this coupling may cause the appearance of a peak in their cross-correlogram (Brody, 1999b). To test whether such coupled magnitude could account for CC peaks observed in this study, the above simulation procedure was repeated, but the computed probability vectors were normalized so that each trial had the same mean before spike-train simulation ("normalized smoothed simulation"). Such normalization removed any coupled magnitude effects by removing any between-trial structure in the responses and leaves only within-trial rate coupling.

To test whether coupled rate changes, the tendency for significant within-trial changes in the time-varying response of one neuron to be coupled to those of another neuron, played a role in producing CC peaks, the above data were compared with the results of yet another simulation: the "rate-change simulation." Trials were simulated as above, with one difference: no smoothing was done across times of significant within-trial firing rate changes (Pauluis and Baker, 2000). The algorithm used to determine times of firing rate change worked locally on sets of three successive action potentials. The interspike interval between the first two spikes served as the momentary average of the gamma distribution describing the firing of that neuron, the order having been estimated as described above. This provided the means to determine whether the third was unlikely (p < 0.01) to have been produced by the same gamma-distributed firing rate. Such "unlikely spikes" were determined to represent significant changes in firing rate (Pauluis and Baker, 2000).

CCs based on these vectors emphasized coupling of firing rate change times. Comparisons with smoothed simulations thus revealed the importance of coupled rate changes. These CCs were also compared with CCs calculated from "normalized rate-change simulations."

To test whether any coupled rate changes reflected purely coupled latency, the tendency for the initial response onsets of the neuron pairs to vary together between trials (Brody, 1999a), the initial response changes, determined as part of the analysis described in the above paragraphs, for interacting neuron pairs were examined. Any coupled latency will be reflected in a simple Pearson's product moment correlation between the initial latencies of the two neurons.

Because interpretations of CCs are complicated by peaks in the autocorrelograms of the neurons (Rosenberg et al., 1989; Nowak and Bullier, 2000), we bolstered our analyses with several supporting measures. First, the autocorrelograms of all units were examined, and only data for which the autocorrelograms did not differ significantly between tastants (according to ² tests) was further analyzed. Second, significance tests were constructed based on one final simulation, in which independent pairs of spike trains were constructed based on the observed single-unit firing rates and spiking distributions, but ignoring any single-trial interactions in the data. The PSTH of each neuron was used as the underlying firing probability vector in this analysis, and a number of trials equal to that actually collected were generated from this one probability vector. Ten such simulations were prepared for each pair of neurons, and binwise confidence intervals for chance interactions were constructed from the means and SDs of these samples. Actual and simulated CCs were smoothed with a 20 msec Gaussian to limit the effect of outliers.

To ensure conservatism of interpretation, the significance criterion was set at p < 0.001. Furthermore, only sections of CC that were significant for 100 msec, or for which the peak was greater than twice the amplitude of the significance criterion, were counted as significant. These four criteria, set after examining real and simulated CCs, made the occurrence of appreciable rates of spurious significance very unlikely.

The results differed little for significance criteria set between 0.01 and 0.001, a fact that stands testament to the robustness of the phenomena. Furthermore, the fact that the number of significantly interacting pairs did not scale with size of recorded ensemble [the correlation between these two variables did not reach significance (r² = 0.24; Fisher's z (10) = 1.77; NS)] adds credence, above and beyond the stringent significance criterion used, to our claim that the observed interactions did not occur randomly. Results using a variety of alternative techniques for building CC confidence intervals (Brillinger, 1992; Politis et al., 1992) also led to similar conclusions. The simulations were used because they facilitated additional analysis, as described above.

Between-simulation comparisons used the maximum height of the peak of the produced CCs, normalized to the peak height of the empirical CC. Widths at half-height were calculated as the distance between the sides of a peak at the height halfway between the peak and the average value of the correlation across all lags (or between the peak and zero, if the CC showed biphasic peaks). The similarities between cross-correlations was calculated using Pearson's product Moment correlation (r), computed using the values of one cross-correlation for lags between ±1 sec as one vector of numbers and the values of the second cross-correlation for the identical lags as the paired vector. This analysis permitted the estimation of similarities in the timing and direction of cross-correlation peaks. Two cross-correlations with similarly timed peaks in the same direction (positive or negative) will have an r > 0, two with similarly timed peaks in different directions (one positive and one negative) will have an r < 0, and two with unrelated peak timings will have an r = 0.

To further examine the information content of GC cross-correlations, the responses of whole ensembles were subjected to linear discriminant analysis (Gochin et al., 1994; Schoenbaum and Eichenbaum, 1995; Nicolelis et al., 1997b; Chapin, 1999). This multivariate technique looks for structure in the covariance matrix of the variables and uses this covariance structure to build optimal linear estimators for the responses to each stimulus. In this case, the covariance structure reflects information related to that which is exposed in the cross-correlations. The firing rate vectors of individual trials (described above) were the input to the algorithm (the input vectors were n × t numbers long, with n indicating number of neurons and t indicating number of time bins); the number of correctly and incorrectly classified trials was the output. Each trial was classified on the basis of the distance between that trial and the estimator of each stimulus, with the smallest distance determining the "predicted" tastant. An identical analysis was then run on the same data, with an additional random time jitter (of up to ±6 bins) introduced to each trial. This jitter destroyed the between-neuron structure, leaving only individual responses intact. Thus, the difference in error rates between the jittered and unjittered data provided a test of the usefulness of between-neuron patterns for coding tastants. The pattern of correct and incorrect classifications was transformed into "bits," a basic measure of the information content of a signal (Krippendorff, 1986).

Single-unit analysis. Details of the analysis of single-unit firing rates can be found by Katz et al. (2001).

Histology. After the end of recording, rats were killed with sodium pentobarbital (150 mg/kg) and perfused first with PBS and then 5% formalin in PBS. Electrolytic lesions made after perfusion (70 µA for 7 sec) in fixed, 80 µm coronal slices stained with cresyl violet revealed that all recordings were made in GC (dysgranular insular cortex).

	RESULTS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

GC neurons fire coherently on the timescale of hundreds of milliseconds

The sample consisted of 237 neuron pairs (mean yield per rat was 20 pairs; range of six to 45), 167 within a single hemisphere and 70 between hemispheres.

Stimulus-specific neural interactions occurred frequently at a relatively broad timescale (half-widths of 75-500 msec). Representative examples, revealing the relative likelihood (at a given lag) that a neuron 2 spike follows or precedes a neuron 1 spike, are shown in Figure 2A-D. The neuron pair in Figure 1A interacted significantly only in the presence of Na. The CCs for other tastants were not significant. The corresponding color-coded PSTHs are shown to their right (the color key is at the bottom). Neuron 1 produced an excitatory firing rate change to Na and a significant inhibitory response to CA. Neuron 2 did not respond in a statistically significant manner to any of the proffered tastants. This pattern was not uncommon: 67% (57 of 85) of the interacting pairs contained one (46%; 39 of 85) or two (21%; 18 of 85) neurons that were either unresponsive or that produced similarly shaped responses to all of the tested tastants (Table 1).

View this table: [in this window] [in a new window]   Table 1.   Summary of cross-correlations

The CCs in Figure 2B show a significant negative correlation peak in the presence of Nic and CA. The occurrence of a "negative" correlation means that neuron 2 tended not to fire in the vicinity of a neuron 1 spike. This pattern of interaction is not reflected by significant inhibition of the single-unit responses to these tastants.

Figure 2C shows one of the relatively few (17.6%; 15 of 85) sharp interactions (width of a few milliseconds) that was observed. In the presence of Na and Q, the CCs of this pair contained a sharp peak (centered on a broader peak) near 0 sec lag time. This pair, like 33% (5 of 15) of those showing such sharp CCs, was isolated from a single wire (9.8% of the total number of significant peaks came from single-wire pairs). None of the taste specificity of the CCs was paralleled by taste specificity in the PSTHs of the neurons.

Sharp CCs were not observed between any interhemispheric pairs, but many such pairs (22%;19 of 85) showed broad interactions. An example is shown in Figure 2D. This pair showed significant positive correlation peaks to Nic, NA, and Q and a biphasic negative-positive interaction to CA. The PSTHs of these neurons were not significantly different for any applied tastants. In all other regards, interhemispheric interactions were generally similar to interactions in intrahemispheric pairs.

Overall, 36% (85 of 237) of the GC pairs interacted in the presence of at least one tastant (mean half-width of 300 ± 15 msec). Of these, 59% (50 of 85) interacted significantly in the presence of more than one tastant (that is, the leftmost bar in Fig. 3A is 41% of the total number of pairs represented in the histogram). In total, 13% of the total number of CCs (177 of 1401) were significant. Most tastants induced similar numbers of significant CCs (Table 1); nicotine produced significantly fewer, perhaps because of relative concentration differences or adaptation (Dessirier et al., 2000). Palatability of the tastants did not determine the pattern of significant CCs (Table 1).

View larger version (53K): [in this window] [in a new window]   Figure 3.   Summary of cross-correlation patterns in GC neuron pairs. A, Frequency histogram showing the number of significant cross-covariances, of a total possible six (one for each tastant), for each neuron pair. Most pairs only covaried in response to one or two tastants, but a substantial number covaried in response to four or five. B, Frequency histograms showing that CCs produced to different tastants by the same neuron pair tend to be similar (white bars), in that the Pearson's r between such CCs tends to be greater than zero. For comparison, the black bars show the Pearson's r between cross-covariances that share only one neuron in common; note that this distribution is centered on zero.

When neuron pairs did interact to a subset of tastants (e.g., two to five), the subset of CCs with significant peaks tended to have similar "morphologies;" that is, they tended to show similarly timed peaks in the same direction (Fig. 2B-D). This observation, summarized in Figure 3B, shows the similarity between pairs of cross-correlations at the same set of lags. The white bars show the similarity (Pearson's product moment correlation r) between the shapes of significant CCs involving the same pair of neurons (for instance, the CCs to Na and Q for the pair shown in Fig. 2C); for this subgroup of pairs, the mean r = 0.36. Black bars show the similarity between any two significantly peaked CCs that have only one neuron in common; for these comparisons, the average r = 0.02. The two distributions are significantly different from each other (²₍₂₀₎ = 703.5; p < 0.000001), as are the mean correlations (Fisher's z (163) = 14.15; p < 0.000001).

Taste-specific assemblies of GC neurons change their firing rates together

The fact that many neuron pairs interacted similarly to a subset of tastants suggests that taste specificity can be found in the makeup of the groups of GC neurons interacting in response to each tastant. Figure 4 shows such taste-specific assemblies recorded for three of the 12 rats. Circles at the intersections in each half-matrix represent the pairwise interactions. The size of a circle represents the significance of the interaction (white for negative and dark gray for positive); small black circles represent nonsignificant interactions (for the neurons showing biphasic patterns of interaction, the size of the circle is keyed to the strongest interaction, and the circle is shown half-white and half-gray). Although positive interactions predominated, both were plentiful. For each tastant, a distinct but overlapping set of neurons interacted in response to presentation of each tastant. Many pairs (Fig. 4C, neurons 4 and 7) only interacted in the presence of a single tastant (in this case, Na), whereas others interacted in the presence of many to a range of degrees (Fig. 4B, neurons 1 and 2).

View larger version (21K): [in this window] [in a new window]   Figure 4.   Taste-specific assemblies of interacting neurons. A-C, The half-matrix of the neuron pairs in three separate ensembles, showing which pairs produced significant CC peaks in response to the different tastants. Small black circles denote lack of significant interaction, and the sizes of shaded circles denote relative significance level. White circles denote excitatory interactions, and dark gray circles denote inhibitory interactions (for the rare neuron pairs that showed biphasically significant CCs, the size of the circle refers to the most significant peak, and half of the circle is shown in each shade). In virtually all cases, the assemblies for different tastants are distinct (some pairs interact to only one tastant) but overlap (some pairs interact for several tastants).

We next determined whether such coherence structure within ensembles, measured at the timescale of CC widths, could be used to identify applied tastants. This possibility was tested using linear discriminant analysis (LDA). Single-trial firing rates of all neurons in ensembles (including those not involved in significant cross-correlations) were turned into rate vectors and binned at a range of timescales (100, 125, 160, 200, and 250 msec, chosen to produce whole numbers of bins). Across-ensemble vectors were used as input to the LDA algorithm, which produced a linear function that allowed maximal discrimination between the tastants, measured as the error rate of response classification.

Figure 5 presents the results of the LDA. Figure 5A shows, for the range of bin sizes tested, the significance of the difference in percentage of correct using normal and shuffled (±1-6 bins) data. Eliminating the covariance structure between neurons in the ensemble impairs performance for bin sizes 160 and 200 msec (both p values < 0.02), a range that is consistent with the timescale of the observed cross-correlations. Figure 5B shows the bits of information (that is, the amount of information useful for identifying tastants from the responses) for all ensembles at a bin width of 160 msec. For the vast majority of ensembles (n = 10), performance with shuffled data (white bars) was worse than performance with normal data (black bars); even with the small number of trials and neurons, the difference was significant for five individual ensembles. Figure 5C summarizes these latter data, showing the average amount of information (in bits) for normal and shuffled data. The shuffled data contained significantly lower amounts of information (t₍₁₁₎ = 2.73; p < 0.02).

View larger version (26K): [in this window] [in a new window]   Figure 5.   Use of the GC ensembles to discriminate between administered tastants. A, Quantification of the reduction in discriminative performance of LDA after time shifting (±6 bins) of single-unit spike trains within GC ensembles. The abscissa is bin size for the analysis, and the ordinate is significance level (p value of the t test comparing error rates for normal and shifted data). The dashed horizontal line shows p = 0.05. When the time bins used were between 160 and 200 msec, LDA performed significantly worse when time shifting eliminated the true between-neuron coherence in the data. B, Information contained within ensembles obtained from each rat, using a bin size of 160 msec. Black bars show the available information (in bits) using the real data, and white bars show the bits using time-shifted data. C, Summary of the data in B, using real and time-shifted data. *p < 0.02.

As noted, significant CCs can arise from a number of sources. We therefore undertook a series of analyses to determine what the observed interactions may "mean" in terms of the processing of gustatory information. The analytic process required obtaining the times of significant firing rate changes in individual trials. After quantification of these times, initial latencies for every trial involved in significant interactions were correlated. The average correlation (0.02) suggests that coupled latencies did not contribute to the CCs.

To test for other possible sources of CCs, simulations that excluded aspects of the data were produced from individual trials of the single-unit responses, and CCs calculated from these various simulations were compared. Figure 6A presents CCs constructed from the Na interaction of the neuron pair shown in Figure 2A. The empirical CC is shown in a thick solid line, the CC calculated from normalized rate-change simulations is shown in a thinner solid line, from normalized smoothed simulations in a dashed line, and from control simulations (all single-trial effects excluded) in a dot-dash line. The strongest interaction is seen in the real data, but the rate-change simulation produced a substantial CC, and the smoothed simulation produced an even lower CC. The corresponding PSTHs for neuron 1 can be seen to the right; each is remarkably similar to the others, clearly demonstrating that the CCs do not depend on the summed single-unit responses.

View larger version (18K): [in this window] [in a new window]   Figure 6.   CCs in GC are related to coupled within-trial rate changes. A, A set of CCs for the same neuron pair shown in Figure 1A. The thick solid line shows the observed CC (i.e., from real data) in response to NaCl. CCs calculated between simulations of the responses of that pair of neurons to NaCl are also shown: one in which all single-trial information has been discarded (dot-dash line), one in which only slow single-trial firing rate variations have been included (dashed line), and one in which times of significant firing rate change are sharp (thin solid line; see Results). The PSTHs for the neuron that produced a significant single-unit response to NaCl, for each condition, are shown to the right. The vast differences in CC are not reflected in appreciable differences between PSTHs for the real and variously simulated spike trains. B, Comparison between the maximum heights of the significant CCs produced from rate-change vector simulations and smoothed vector simulations. Proportion of the real data maximum peak height is on the ordinate; the error bars are SEM. The difference between the two simulations is highly significant. C, Similar comparison between CCs produced by scaled rate-change vector simulations and unscaled rate-change vector simulations. The difference is not significant.

The results of these analyses are summarized in Figure 6B. The difference between the rate change and smoothed simulations was highly significant (t₍₁₇₅₎ = 6.73; p < 0.000001). Coupled shifts in firing rate contribute significantly to the CCs. Figure 6C compares the rate-change and normalized rate-change simulations. Had magnitude coupling contributed to the interactions, the latter should have been significantly smaller than the former. Such was not the case: the size of the peaks are virtually identical for the two simulations (t₍₁₇₅₎ < 1).

In summary, the information reflected in the taste-specific coherent assemblies is preserved in the broad time course of individual trial responses, and significant amounts of that information can be traced to sudden changes in instantaneous firing rates, even when these changes do not show up in PSTHs.

	DISCUSSION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

GC neurons cohere during tasting

Information related to gustation evolves in GC responses (Katz et al., 2001a) on the same timescale as gustatory perception (Halpern, 1985). We investigated interactions between GC neurons during tastant responses and explored how such interactions might be related to the single-unit processing of gustatory stimuli. These data reveal that pairs of GC neurons interact across hundreds of milliseconds, that the interactions are taste specific, that they do not depend on the shapes of the PSTH or on stimulus palatability, and that they define distinct, but overlapping, neural assemblies that respond to the presence of each tastant by undergoing coupled changes in firing rate. Our analyses demonstrate that these couplings could, in principle, be used to discriminate between tastants.

Previous studies reported taste-specific CCs between GC neurons on a narrow timescale (±20 msec) (Nakamura and Ogawa, 1997; Yokota et al., 1997). We observed a small number of such CCs (Fig. 2C), the properties of which were consistent with these previous reports. The fact that sharp CCs were not found in interhemispheric pairs and the relatively high likelihood that they were found between neurons isolated from a single wire suggests that monosynaptic connections are formed preferentially between neurons in close spatial proximity to each other.

However, our findings differ from those reported previously in two ways. First, the effects of anesthesia, which may change system dynamics (Nowak and Bullier, 2000; Gaese and Ostwald, 2001), were avoided through the use of awake rats. Second, the broad (hundreds of milliseconds) interactions observed here (Fig. 2) are widespread and systemic, even appearing between neurons in opposite hemispheres (Fig. 2D). The systemic nature of this phenomenon, together with the well known millisecond-to-millisecond variability of neural responses (Shadlen and Newsome, 1998), suggests that broad interactions may better represent the processing of tastants than do tight (indicative of monosynaptic connections) interactions.

Gustatory processing via coherent rate changes in taste-specific assemblies

What do broad CCs have to do with chemosensory processing? Research in olfaction suggests that specific information may be carried in patterns of transient synchrony between neurons (Wehr and Laurent, 1996). In this study, we found that the shapes of CCs provide some information about tastant identity, in that most pairs of neurons only interact after application of only a subset of tastants. Figures 2B-D and 3 suggest, however, that significant CCs in GC may simply be a "flag" to the membership of a neuron pair in a processing group. Specifically, significant CCs can serve to flag the membership of a neuron pair in a stimulus-specific neural assembly (Fig. 4). The formation of such assemblies has been related to system function by a number of researchers studying other systems (Singer, 1990; Abeles et al., 1993; Vaadia et al., 1995; Stopfer et al., 1997; Nicolelis et al., 1998; Laurent, 1999; Christensen et al., 2000). Although it is difficult to ascertain whether interactions between neurons provide "more" coding-relevant information than single-neuron firing patterns (as they seem to in olfaction), elimination of this information via random time shifting of individual responses significantly degraded the ability of the ensembles to discriminate between applied tastants (Fig. 5). Although it should be remembered that such analyses represent proof only of an "observer's" ability to use neural interactions (and not that the rat does so), this supports the suggestion that these interactions between GC neurons can play a role in the processing of tastant stimuli.

The performed analyses permit the exclusion of common source, coupled latency, and coupled magnitude as sources of the interactions between taste-specific subset of coherent neurons (Fig. 6). Other possibilities can also be eliminated. Between-trial variations in the rats' handling of taste stimuli, for instance, could give rise to peaked CCs, but those CCs would be taste specific only if tastant delivery was similarly variable for a subset of stimuli and not at all variable for others. The fact that different but overlapping subsets of neurons interacted to different tastants (Fig. 2) further reduces the feasibility of this hypothesis. Alternatively, it might be argued that palatability-specific orofacial behaviors drove the between-neuron coupling. The fact that pairwise interactions seldom fell out according to palatability (Table 1), however, eliminates this explanation.

Therefore, we conclude that coupled changes in firing rate (a portion of which must "live" at a timescale smaller than the 50 msec filter used to produce our simulations) are the underlying source of GC interactions. Subsets of neurons in GC became coupled after the presentation of particular tastants, and the responses of neurons in that ensemble changed in concert with those of others. The action of gustatory assemblies is thus intimately related to rate changes on a timescale similar to that of taste behavior (Halpern, 1985). The failure of PSTHs to reflect the significant interactions (Figs. 2, 6) suggests that the firing rate coupling was not locked to stimulus onset. Rather, it is an internal, single-trial process.

The conclusions offered here are consistent with recent work from a variety of preparations. Researchers have suggested that networks of neurons form functional ensembles (Welsh et al., 1995; Laubach et al., 2000) and that such ensembles undergo coupled firing rate changes (Seidemann et al., 1996). Recent simulations suggest explicit mechanisms of between-neuron coupling to explain such a process in olfaction (Bazhenov et al., 2001). The current results are also consistent with research suggesting that the taste specificity of neural responses is partially determined by interconnections among neurons (Ogawa et al., 1998; Smith and Li, 1998).

In summary, the time-varying tastant responses produced by GC neurons in awake rats, even those that do not translate into taste-specific PSTHs, can be attributed to coherent rate changes in taste-specific assemblies. These data suggest that taste perception may involve the formation and action of such assemblies and to the way in which ensembles of neurons work together.

	FOOTNOTES

Received Oct. 15, 2001; revised Nov. 27, 2001; accepted Dec. 3, 2001.

This research was supported by National Institutes of Health Grants DC-00403 (D.B.K.), DC-01065 (S.A.S.), and DE-11121 (M.A.L.N.) and by a grant from the Philip Morris Research Center. We gratefully acknowledge the advice of Drs. Stuart Baker, Robert Erickson, Asif Ghazanfar, Rebecka Jornsten, Ben Rubin, Geoff Schoenbaum, Marshall Shuler, Alan Spector, and Johan Wessberg.

Correspondence should be addressed to Donald B. Katz, Box 2309, Duke University Medical Center, Durham, NC 27710. E-mail: dkatz{at}neuro.duke.edu.

	REFERENCES

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

• Abeles M, Bergman H, Margalit E, Vaadia E (1993) Spatiotemporal firing patterns in the frontal cortex of behaving monkeys. J Neurophysiol 70:1629-1638[Abstract/Free Full Text].
• Adachi M, Ohshima T, Yamada S, Satoh T (1989) Cross-correlation analysis of taste neuron pairs in rat solitary tract nucleus. J Neurophysiol 62:501-509[Abstract/Free Full Text].
• Aertsen AM, Gerstein GL, Habib MK, Palm G (1989) Dynamics of neuronal firing correlation: modulation of "effective connectivity." J Neurophysiol 61:900-917[Abstract/Free Full Text].
• Bair W, Zohary E, Newsome WT (2001) Correlated firing in macaque visual area MT: time scales and relationship to behavior. J Neurosci 21:1676-1697[Abstract/Free Full Text].
• Baker SN, Lemon RN (2000) Precise spatiotemporal repeating patterns in monkey primary and supplementary motor areas occur at chance levels. J Neurophysiol 84:1770-1780[Abstract/Free Full Text].
• Bazhenov M, Stopfer M, Rabinovich M, Abarbanel HD, Sejnowski TJ, Laurent G (2001) Model of cellular and network mechanisms for odor-evoked temporal patterning in the locust antennal lobe. Neuron 30:569-581[ISI][Medline].
• Brillinger DR (1992) Nerve cell spike train data analysis: a progression of technique. J Am Stat Assoc 87:260-271.
• Brody CD (1999a) Correlations without synchrony. Neural Comput 11:1537-1551[Abstract/Free Full Text].
• Brody CD (1999b) On artefactual spike train cross-correlations. Neurocomputing 26-27:957-962.
• Chapin JK (1999) Population-level analysis of multi-single neuron recording data: multivariate statistical methods. In: Methods for neural ensemble recordings (Nicolelis MAL, ed), pp 193-228. Boca Raton, FL: CRC.
• Christensen TA, Pawlowski VM, Lei H, Hildebrand JG (2000) Multi-unit recordings reveal context-dependent modulation of synchrony in odor-specific neural ensembles. Nat Neurosci 3:927-931[ISI][Medline].
• Dessirier J-M, Simons CT, Sudo M, Sudo S, Carstens E (2000) Sensitization, desensitization, stimulus-induced recovery (SIR) of rat trigeminal caudalis neurons to repeated oral application of capsaicin, nicotine. J Neurophysiol 84:1851-1862[Abstract/Free Full Text].
• Di Lorenzo PM, Monroe S (1997) Transfer of information about taste from the nucleus of the solitary tract to the parabrachial nucleus of the pons. Brain Res 763:167-181[ISI][Medline].
• Erickson RP (2001) The evolution and implications of population and modular neural coding ideas. Prog Brain Res 130:9-29[Medline].
• Gaese BH, Ostwald J (2001) Anesthesia changes frequency tuning of neurons in the rat primary auditory cortex. J Neurophysiol 86:1062-1066[Abstract/Free Full Text].
• Ghazanfar AA, Krupa DJ, Nicolelis MA (2001) Role of cortical feedback in the receptive field structure and nonlinear properties of somatosensory thalamic neurons. Exp Brain Res 141:88-100[ISI][Medline].
• Gochin PM, Colombo M, Dorfman GA, Gerstein GL, Gross CG (1994) Neural ensemble coding in inferior temporal cortex. J Neurophysiol 71:2325-2337[Abstract/Free Full Text].
• Halpern BP (1985) Time as a factor in gustation: temporal patterns of taste stimulation and response. In: Taste, olfaction, and the central nervous system (Pfaff DW, ed), pp 181-209. New York: Rockefeller UP.
• Katz DB, Simon SA, Nicolelis MA (2001a) Dynamic and multimodal responses of gustatory cortical neurons in awake rats. J Neurosci 21:4478-4489[Abstract/Free Full Text].
• Katz DB, Simon SA, Nicolelis MAL (2001b) Electrophysiological studies of gustation in awake rats. In: Methods and frontiers in the chemical senses (Simon SA, Nicolelis MAL, eds), pp 339-357. Boca Raton, FL: CRC.
• Kosar E, Grill HJ, Norgren R (1986) Gustatory cortex in the rat. I. Physiological properties and cytoarchitecture. Brain Res 379:329-341[ISI][Medline].
• Krippendorff K (1986) In: Information theory: structural models for qualitative data. Beverly Hills, CA: Sage UP.
• Laubach M, Wessberg J, Nicolelis MA (2000) Cortical ensemble activity increasingly predicts behaviour outcomes during learning of a motor task. Nature 405:567-571[Medline].
• Laurent G (1999) A systems perspective on early olfactory coding. Science 286:723-728[Abstract/Free Full Text].
• Nakamura T, Ogawa H (1997) Neural interaction between cortical taste neurons in rats: a cross-correlation analysis. Chem Senses 22:517-528[Abstract/Free Full Text].
• Nicolelis MA, Katz D, Krupa DJ (1998) Potential circuit mechanisms underlying concurrent thalamic and cortical plasticity. Rev Neurosci 9:213-224[ISI][Medline].
• Nicolelis MAL, Ghazanfar AA, Faggin BM, Votaw S, Oliveira LMO (1997a) Reconstructing the engram: simultaneous, multisite, many single neuron recordings. Neuron 18:529-537[ISI][Medline].
• Nicolelis MAL, Lin RCS, Chapin JK (1997b) Neonatal whisker removal reduces the discrimination of tactile stimuli by thalamic ensembles in adult rats. J Neurophysiol 78:1691-1706[Abstract/Free Full Text].
• Nishijo H, Norgren R (1991) Parabrachial gustatory neural activity during licking by rats. J Neurophysiol 66:974-985[Abstract/Free Full Text].
• Nowak LG, Bullier J (2000) Cross-correlograms for neuronal spike trains. Different types of temporal correlation in neocortex, their origin and significance. In: Time and the brain (Miller R, ed), pp 53-96. Sidney: Harwood.
• Ogawa H, Hasegawa K, Otawa S, Ikeda I (1998) GABAergic inhibition and modifications of taste responses in the cortical taste area in rats. Neurosci Res 32:85-95[ISI][Medline].
• Pauluis Q, Baker SN (2000) An accurate measure of the instantaneous discharge probability, with application to unitary joint-even analysis. Neural Comput 12:647-669[Abstract/Free Full Text].
• Perkel DH, Gerstein GL, Moore GP (1967) Neuronal spike trains and stochastic point processes. II. Simultaneous spike trains. Biophys J 7:419-440.
• Phillips MI, Norgren R (1970) A rapid method for permanent implantation of an intraoral fistula in rats. Behav Res Methods Instrum Comput 2:124.
• Politis DN, Romano JP, Lai TL (1992) Bootsrap confidence bands for spectra and cross-spectra. IEEE Trans Signal Proc 40:1206-1215.
• Rosenberg JR, Amjad AM, Breeze P, Brillinger DR, Halliday DM (1989) The Fourier approach to the identification of functional coupling between neuronal spike trains. Prog Biophys Mol Biol 53:1-31[ISI][Medline].
• Salinas E, Sejnowski TJ (2001) Correlated neuronal activity and the flow of neural information. Nat Rev Neurosci 2:539-550[ISI][Medline].
• Schoenbaum G, Eichenbaum H (1995) Information coding in the rodent prefrontal cortex. II. Ensemble activity in orbitofrontal cortex. J Neurophysiol 74:751-762[Abstract/Free Full Text].
• Seidemann E, Meilijson I, Abeles M, Bergman H, Vaadia E (1996) Simultaneously recorded single units in the frontal cortex go through sequences of discrete and stable states in monkeys performing a delayed localization task. J Neurosci 16:752-768[Abstract/Free Full Text].
• Shadlen MN, Newsome WT (1998) The variable discharge of cortical neurons: implications for connectivity, computation, and information coding. J Neurosci 18:3870-3896[Abstract/Free Full Text].
• Singer W (1990) The formation of cooperative cell assemblies in the visual cortex. J Exp Biol 153:177-197[Abstract/Free Full Text].
• Smith DV, Li CS (1998) Tonic GABAergic inhibition of taste-responsive neurons in the nucleus of the solitary tract. Chem Senses 23:159-169[Abstract].
• Stopfer M, Bhagavan S, Smith BH, Laurent G (1997) Impaired odour discrimination on desynchronization of odour-encoding neural assemblies. Nature 390:70-74[Medline].
• Vaadia E, Haalman I, Abeles M, Bergman H, Prut Y, Slovin H, Aertsen A (1995) Dynamics of neuronal interactions in monkey cortex in relation to behavioural events. Nature 373:515-518[Medline].
• Wehr M, Laurent G (1996) Odour encoding by temporal sequences of firing in oscillating neural assemblies. Nature 384:162-166[Medline].
• Welsh JP, Lang EJ, Sugihara I, Llinás R (1995) Dynamic organization of motor control within the olivocerebellar system. Nature 374:453-457[Medline].
• Yokota T, Eguchi K, Satoh T (1996) Correlated discharges of two neurons in rat gustatory cortex during gustatory stimulation. Neurosci Lett 209:204-206[Medline].
• Yokota T, Eguchi K, Satoh T (1997) Sensitivity of rat cortical neurons in distinguishing taste qualities by individual and correlative activities. Chem Senses 22:363-373[Abstract/Free Full Text].

---

Copyright © 2002 Society for Neuroscience  0270-6474/02/2251850-08$05.00/0

This article has been cited by other articles:

				S. Tsujimoto, A. Genovesio, and S. P. Wise Transient Neuronal Correlations Underlying Goal Selection and Maintenance in Prefrontal Cortex Cereb Cortex, March 20, 2008; (2008) bhn033v1. [Abstract] [Full Text] [PDF]

				S. E. Grossman, A. Fontanini, J. S. Wieskopf, and D. B. Katz Learning-Related Plasticity of Temporal Coding in Simultaneously Recorded Amygdala-Cortical Ensembles J. Neurosci., March 12, 2008; 28(11): 2864 - 2873. [Abstract] [Full Text] [PDF]

				J. V. Verhagen and D. B. Katz More Time to Taste. Focus on "Variability in Responses and Temporal Coding of Tastants of Similar Quality in the Nucleus of the Solitary Tract of the Rat" J Neurophysiol, February 1, 2008; 99(2): 413 - 414. [Full Text] [PDF]

				L. M. Jones, A. Fontanini, B. F. Sadacca, P. Miller, and D. B. Katz Natural stimuli evoke dynamic sequences of states in sensory cortical ensembles PNAS, November 20, 2007; 104(47): 18772 - 18777. [Abstract] [Full Text] [PDF]

				M. G. Veldhuizen, G. Bender, R. T. Constable, and D. M. Small Trying to Detect Taste in a Tasteless Solution: Modulation of Early Gustatory Cortex by Attention to Taste Chem Senses, July 1, 2007; 32(6): 569 - 581. [Abstract] [Full Text] [PDF]

				R. L. Rennaker, C.-F. F. Chen, A. M. Ruyle, A. M. Sloan, and D. A. Wilson Spatial and Temporal Distribution of Odorant-Evoked Activity in the Piriform Cortex J. Neurosci., February 14, 2007; 27(7): 1534 - 1542. [Abstract] [Full Text] [PDF]

				R. Accolla, B. Bathellier, C. C. H. Petersen, and A. Carleton Differential Spatial Representation of Taste Modalities in the Rat Gustatory Cortex J. Neurosci., February 7, 2007; 27(6): 1396 - 1404. [Abstract] [Full Text] [PDF]

				A. Fontanini and D. B. Katz State-Dependent Modulation of Time-Varying Gustatory Responses J Neurophysiol, December 1, 2006; 96(6): 3183 - 3193. [Abstract] [Full Text] [PDF]

				J. I. Glendinning, A. Davis, and M. Rai Temporal coding mediates discrimination of "bitter" taste stimuli by an insect. J. Neurosci., August 30, 2006; 26(35): 8900 - 8908. [Abstract] [Full Text] [PDF]

				C. H. Lemon and D. V. Smith Influence of response variability on the coding performance of central gustatory neurons. J. Neurosci., July 12, 2006; 26(28): 7433 - 7443. [Abstract] [Full Text] [PDF]

				J. R. Stapleton, M. L. Lavine, R. L. Wolpert, M. A. L. Nicolelis, and S. A. Simon Rapid Taste Responses in the Gustatory Cortex during Licking J. Neurosci., April 12, 2006; 26(15): 4126 - 4138. [Abstract] [Full Text] [PDF]