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

HOME  |   SEARCH  |   ARCHIVE  |   SUBSCRIBE  |   CONTACT  |   HELP

The Journal of Neuroscience, October 18, 2006, 26(42):10826-10846; doi:10.1523/JNEUROSCI.2091-06.2006

This Article Abstract Full Text (PDF) Supplemental data Submit an eLetter View eLetters for this paper 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 (11) Citing Articles via Google Scholar Google Scholar Articles by Conway, B. R. Articles by Livingstone, M. S. Search for Related Content PubMed PubMed Citation Articles by Conway, B. R. Articles by Livingstone, M. S.

 Previous Article  |  Next Article 

Behavioral/Systems/Cognitive
Spatial and Temporal Properties of Cone Signals in Alert Macaque Primary Visual Cortex

Bevil R. Conway^1,2 and Margaret S. Livingstone¹

¹Department of Neurobiology, Harvard Medical School, Boston, Massachusetts 02115, and ²Wellesley College, Wellesley, Massachusetts 02481

	Abstract

Top Abstract Introduction Materials and Methods Results Discussion References

 
Neurons in the lateral geniculate nucleus cannot perform the spatial color calculations necessary for color contrast and color constancy. Under neutral-adapting conditions, we mapped the cone inputs (L, M, and S) to 83 cone-opponent cells representing the central visual field of the next stage of visual processing, primary visual cortex (V1), to determine how the color signals are spatially transformed. Cone-opponent cells, constituting 10% of V1 cells, formed two populations, red–green (L vs M; 66 of 83) and blue–yellow (S vs L+M; 17 of 83). Many cone-opponent cells (48 of 83) were double-opponent, with circular receptive-field centers and crescent-shaped surrounds (0.63° offset) that had opposite chromatic tuning to the centers and a time-to-peak 11 ms later than the centers. The remaining cone-opponent cells were either spatially opponent in only one cone system (20 of 83) or lacked spatial opponency (15 of 83). Cells lacking spatial opponency had smaller receptive fields (0.5–0.7°) than spatial-opponent cell centers (1°). We found that red–green cells received S-cone input, which aligned with M input, and, unlike blue–yellow cells, red–green cells gave push–pull responses: receptive-field centers of red-ON cells were excited by both L increments (bright red) and M decrements (dark red) and were suppressed by both L decrements (dark green) and M increments (bright green). Excitatory responses to decrements were slightly larger than to increments, which may account for the lower detection and discrimination thresholds of decrements shown psychophysically. By virtue of their specialized receptive fields, the neurons described here spatially transform the cone signals and represent the first stage in the visual system at which spatially opponent color calculations are made.

Key words: receptive-field; primate; V1; color; double-opponent; color constancy

	Introduction

Top Abstract Introduction Materials and Methods Results Discussion References

 
Color perception is achieved by comparing the responses of the three cone types (L, M, and S; loosely red, green, and blue). This calculation is typified by retinal blue-ON neurons, which are excited by S-cone activity and suppressed by (L+M)-cone activity (Dacey and Lee, 1994). Blue-ON neurons likely connect to blue–yellow type II cells in the lateral geniculate nucleus (LGN) (Martin et al., 1997). Type I neurons in the LGN also show cone opponency but between L and M cones, so they are often called red–green (De Valois et al., 1958). It remains unclear how the cortex uses these signals to make spatial color comparisons, a step required for color contrast and color constancy (Land, 1977).

Spatially opponent color calculations are not achieved by the LGN. Type II cells do not have center-surround receptive fields (Wiesel and Hubel, 1966; Chichilnisky and Baylor, 1999), as would be required, and the opponent systems of type I cells are spatially opponent in luminance, not color (Wiesel and Hubel, 1966) (see Fig. 1A). The opponent systems of a given type 1 cell are actually centered at the same location (Reid and Shapley, 2002), a configuration captured by a difference of Gaussians (Rodieck, 1965): the center is narrower but stronger than the surround (see Fig. 1A).

View larger version (46K): [in this window] [in a new window]   Figure 1. Receptive field of a parvocellular type I cell in the lateral geniculate nucleus (A), two theoretical receptive fields capable of spatial color calculations (B), and predicted response maps. A, Receptive field and profile view as a difference of Gaussians (left diagrams) and responses to cone-isolating stimuli (right). Cone-isolating stimuli can either increase (+) or decrease (–) the activity of a single cone class. This cell shows excitation to L+ and a broader region of excitation to M– (white regions in the response maps). Note that the optimal stimulus is not a red spot on a green background but rather a bright red spot on a dark red background. The colored squares behind L+, L–, M+, and M– provide a rough idea of what cone-isolating stimuli look like. B, Model 1 shows a set of Gaussians that produce spatial and chromatic opponency. The peaks of all predicted response maps are centered at the same location. Model 2, which also produces chromatic opponency, is built up from many red-ON type I cells: cells feeding the center are wired with excitatory connections and cells feeding the surround with inhibitory connections (adapted from Billock, 1991). The peak response to the L– stimulus will be offset from the peak response to the L+ stimulus.

 
Spatially opponent color calculations could be performed by receptive fields that, like type I cells, can also be represented as a difference of Gaussians (see Fig. 1B, model 1) (Lennie and D’Zmura, 1988). However, the inputs required by this model would need to have receptive fields not seen in type 1 cells or any LGN cell: an L+ center would need to be accompanied by an L– surround rather than an M– surround (Hubel and Livingstone, 1990a).

Billock (1991) has proposed an alternative color-contrast model that combines inputs from many type I cells, low-pass filtering the luminance signal of each input while retaining its chromatic opponency (see Fig. 1B, model 2). The model samples one pool of type 1 cells for the center and one for the surround, a "difference of offset Gaussians" (Young, 1987; Billock, 1995), and predicts that resulting receptive fields will (1) be large compared with those of type I inputs, (2) have spatially offset center and surround excitatory regions, and (3) show residual cross-talk responses to luminance.

Here we test these predictions by studying the receptive fields of cone-opponent neurons in the next stage of visual processing, primary visual cortex (V1). We used neutral-adapting stimuli to overcome criticisms (Shapley and Hawken, 2002) of previous cone-isolating methods (Conway, 2001). The spatiotemporal receptive-field maps generated here can be directly compared with LGN cell maps of Reid and Shapley (2002) and support models like that of Billock.

Finally, it is known that light increments are less visible than light decrements (Whittle, 1986; Bowen et al., 1989). By measuring the responses to stimuli that increase the activity of a given cone (+ stimuli) separately from the responses to stimuli that decrease the activity of a given cone (– stimuli), we tested for asymmetries in V1 responses that might account for this phenomenon.

	Materials and Methods

Top Abstract Introduction Materials and Methods Results Discussion References

 
General design
Experiments were conducted in alert adult male macaque monkeys using chronic recording procedures described previously (Livingstone et al., 1996; Conway and Livingstone, 2003). All experiments were performed according to National Institutes of Health guidelines for the use of animals and with the approval of the Harvard Medical School Standing Committee on the use of animals. Macaques are a useful model for human color vision because psychophysical data from them match those of humans (De Valois et al., 1974; Sandell et al., 1979); the psychophysical results on human color matching are well predicted from the spectral sensitivities of the macaque cones (Baylor et al., 1987). Moreover, using alert animals overcomes the effects of anesthesia, which alter motion tuning (Pack et al., 2001) and may interfere with color processing (Solomon et al., 2004).

Well isolated single units were recorded from V1 using tungsten microelectrodes (Hubel, 1957) (Frederick Haer Company, Bowdoinham, ME). Units were isolated using a dual-window discriminator (BAK Electronics, Germantown, MD) after they were amplified and bandpass filtered (1–10 kHz). Ambiguous single units were discarded if >2% total spikes for a given run were in the first millisecond bin of the interspike interval histogram.

Spikes were collected at 1 ms resolution. Stimuli were presented (in a dark recording chamber) on a 21 inch monitor (Barco Display Systems, Kortrijk, Belgium) with a 60 Hz refresh rate (non-interlaced). The monitor screen was 57 cm in front of the animal for some experiments and 39 cm in others. The animal was given a juice reward for keeping its gaze within 1° of a fixation spot, and spikes were rejected from analysis if they were collected while the animal’s gaze was outside 1° of the fixation spot. Eye position was monitored using either a search coil in a magnetic field (Judge et al., 1980) or the ISCAN (Burlington, MA) system; the field-coil system was from CNC Engineering (Enfield, CT). Data were not corrected for eye position. Eye monitoring was only done to ensure that the monkeys maintained fixation. Microsaccadic eye movements do not have a significant effect on measured receptive fields in well trained animals, as were used here (Tsao et al., 2003).

From a continuous history of spike timing and stimulus position (at 1 ms precision), the activity was reverse correlated to produce a spatiotemporal record of the neural response (Jones and Palmer, 1987; Conway, 2001; Conway and Livingstone, 2003; Livingstone and Conway, 2003). The response maps provide a complete description of the first-order receptive fields of the neurons.

This approach is similar to that used by Reid and Shapley (2002) to map LGN neurons with three minor differences. First, the stimuli used here were sparse noise as opposed to checkerboards that modulate each element in the entire stimulation window. Second, the present stimuli were not constrained to be displayed within a grid, allowing higher-resolution maps (Martinez et al., 2005). Third, the responses to stimuli that increase the activity of a given cone class ("bright" or + stimuli) were analyzed separately from the responses to stimuli that decrease the activity of a given cone class ("dark" or – stimuli).

Screening for cone-opponent cells
Quantitative receptive-field analysis was performed only on neurons that gave some indication of cone opponency to an initial screening test. A description of the cone-isolating stimuli (L+, L–, M+, M–, S+, S–) (Fig. 1A) is given below (see below, Visual stimuli). Neurons were first hand mapped using a computer mouse to adjust the stimulus location on the stimulus monitor while listening to the responses on an audio monitor; hand mapping included small spots of cone-isolating light that were flashed (1 Hz) in and around the receptive fields. For screening, a brief run of the quantitative mapping stimulus (see below, Visual stimuli) was used to test for opposite patterns (excitation vs suppression) in the difference responses (see Eq. 3) to the three cone classes. Data was analyzed on-line; cells showing excitation to one cone type and suppression to another, in any region of their receptive fields, were subsequently mapped with longer stimulus runs. Neurons that were excited or suppressed by both L and M but were nonresponsive to S were classified as luminance selective. The screening method did not require that the neurons show excitatory rebound responses to the cessation of a suppressive stimulus. This is important because color-coding neurons often lack excitatory discharges after cessation of suppressive stimuli (Friedman et al., 2003). Note that the screening methods did not presuppose discrete categories of cone-opponent cells.

Visual stimuli
Cone-isolating stimuli were used to evaluate the cone inputs. These stimuli were presented under neutral-adapting conditions produced by a neutral gray background. The gray-adapting field maintained a constant level of activation of all three cone classes, allowing a direct comparison of the responses to different cone-isolating stimuli. In addition, because the cone-isolating stimuli had spatial structure (they were small spots), they allowed a direct evaluation of the spatial organization of the cone inputs. Small patches of the background were replaced with a color that altered the activity of a single cone class at a time, a method called silent substitution. A + stimulus increased the activity of the respective cone; a – stimulus decreased it. The six cone-isolating stimuli used were as follows (Commission Internationale de l’Eclairage coordinates: x, y, luminance in cd/m²): L+ (0.413, 0.262, 30.7), L– (0.209, 0.277, 21.3), M+ (0.223, 0.336, 28.1), M– (0.378, 0.224, 22.5), S+ (0.254, 0.155, 23.1), and S– (0.376, 0.382, 25.9). Two additional achromatic luminance stimuli were also used in some experiments. These were generated by maximal (white: 70 cd/m²) and minimal (black: <1 cd/m²) gun activities. The adapting background for all experiments was (0.318, 0.271, 25); the mean luminance (25 cd/m²) ensured the stimuli were well within the photopic range. The cone fundamentals of Stockman and Sharpe (2000) were used to estimate cone activity. L and M stimuli were calculated to have equal Michelson’s contrast, as follows:

where L+ is the L cone activity produced by the L+ stimulus, and L_bk is the activity of the L cone produced by the neutral background. The M+ and L+ stimuli had +17% contrast; the M– and L– stimuli had –17% contrast. The S+ and S– stimuli each had 41% contrast. This resulted in the following Weber’s contrasts: [(L+) – (L_bk)/(L_bk)]: L+ was 0.40; L– was –0.29, M+ was 0.40, M– was –0.29; S+ was 1.41, and S– was –0.60. According to colorimetric mixing principles, this produces a + and – sum that is slightly brighter than the background, although the stimuli used were constrained so that the overlap was the background gray. Detailed methods for generating cone-isolating stimuli were given previously (Conway, 2001; Reid and Shapley, 2002).

Note that both L+ and M– look red, but L+ is brighter than the background and M– is darker (Fig. 1A). Therefore, + and – also constitute a luminance stimulus; this is a feature of all cone-isolating stimuli. Some studies have attempted to avoid this luminance component by modulating two cone classes at once, for example by combining an L+ and an M– (Derrington et al., 1984; Lennie et al., 1990). We chose not to use such stimuli so that we could directly assay the inputs of the different cone types and of the + and – states of each cone type, separately.

Stimuli were sparse noise. During any given frame, only two small patches of the neutral background were replaced with a cone-isolating color: one of the spots increased the activity of the given cone type (a + stimulus) and the other decreased it (a – stimulus). The patches were between 0.25 and 0.75° square; the patches had no fixed spatial relationship to each other from one frame to the next. The size of the patches was chosen to be the smallest that still reliably elicited responses. Both patches were the same size and were the same size for all cone stimuli, for a given cell.

In any given frame, the two patches could be displayed anywhere relative to each other, at 0.06° pixel resolution within a stimulation window between 2.25 and 4.5° square centered on the receptive field. If the patches overlapped, the overlap region appeared as the adapting gray. Enough stimuli were presented so that all locations within the stimulation window were sampled evenly; on average, five presentations were delivered to each pixel location; stimulus runs were on an average 17.5 min per cone map (range of 2–55 min). In most cases, the stimulation window extended well beyond both the center and the surround of the receptive field. The same stimulus run (and spike train) was used to generate the response maps for both spots, which was possible because the presentation position of the two spots was independent (the two spots had no consistent spatial relationship from frame to frame). While generating the response map for one of the spots, the activity produced by the other spot was averaged uniformly throughout the map. This average activity, which is common to both + and – maps, is removed during analysis by considering only the difference in activity between the two maps (see below, Data analysis).

The stimulus durations were either 34, 51, or 68 ms (two, three, or four frames) and are stated in the figure legends. The same stimulus duration was used to generate all maps for a given cell. The spatial maps were unaffected by the choice of stimulus duration. In a few cases, a 17 ms stimulus duration (one frame) was used to maximize data collection, although such brief flashes did not always reliably elicit spikes.

Chromatic aberration
As in all studies of color vision, the spatial distribution of short-wavelength light stimuli, including S-cone stimuli, is confounded by chromatic aberration. The effects of chromatic aberration and the preretinal absorption characteristics of the eye are especially significant for neurons recorded at eccentricities >10° (Cottaris, 2003). All neurons in the present study were between 2 and 8°, making it unlikely that chromatic aberration had a significant impact.

Data analysis
Response maps show neural activity, in spikes per second, according to a color scale bar. Reverse correlation typically is reported as a probability of a stimulus preceding a spike, but this can be converted into firing rate according to Bayes’ Rule:

where P(spike/stim) is the probability of a spike given a stimulus (the firing rate), P(stim/spike) is the probability of a stimulus preceding a spike, P_spike is the probability of a spike and P_stim is the probability of a stimulus (the ratio is a constant).

The response maps, which were smoothed with a 0.1° wide Gaussian filter in both spatial dimensions, are depicted with a 0.75° overlying grid to assist in comparing the patterns in different maps for a given cell. The locations of the stimuli were not constrained by this grid. Stimuli were positioned with 1 pixel resolution (0.06 or 0.08° depending on the distance between monkey and monitor).

The first-order reverse correlation for each stimulus run produced two maps, one for each stimulus: a + map, showing the response to an increase in activity of the given cone class, and a – map, showing the response to a decrease in the activity of that cone class. In these "raw" response maps, black indicates 0 firing, and red indicates maximal firing (see Fig. 2, top, + and – maps).

View larger version (41K): [in this window] [in a new window]   Figure 2. Response maps (top) and the spike-triggered average histograms (bottom) to cone-isolating stimuli for a single cone-opponent neuron in primary visual cortex of an alert macaque. The response maps show the spatial structure of cone inputs to the receptive field at the optimal reverse-correlation delay of the receptive-field center (93 ms). The difference maps are the – response subtracted from the + response and show excitation in red and suppression in blue. The neuron showed both spatial and chromatic opponency, i.e., double-opponency. For this neuron, the peak center is aligned with the center of the response maps. The spike-triggered average histograms show the temporal evolution of the response of the peak center and peak surround locations (method 1 was used to index the center and surround; see Materials and Methods). Two SDs above and below the mean background firing rate are indicated by the horizontal dotted lines. Solid arrowheads show the time-to-peak of the center; open arrowheads show the time-to-peak of the surround.

 
These response maps can be analyzed separately or they can be subtracted:

This subtraction discounts the baseline firing rate, which is common to both + and – responses for a given cone map. The baseline firing rate is determined by the average activity of the neuron throughout the stimulus run so that, for stimuli that cause considerably more activity, the baseline is higher (e.g., compare the baseline response to S-cone stimuli, to which the neuron was less sensitive, with that to the M- or L-cone stimuli in Fig. 10). Subtraction was used to quantify the data shown in Figures 5, 7, 11, 13, and 16 and supplemental Figure 2 (available at www.jneurosci.org as supplemental material). In the subtracted maps (see Fig. 2, top, difference map), white indicates no difference between the + and – maps, black indicates maximal suppression, which could result from any combination of real suppression in the + map and excitation in the – map, and red indicates maximal excitation, which could result from any combination of excitation in the + map and suppression in the – map. SDs were taken as the sum of the SD of the + map and that of the – map, divided by 2.

To generate fraction-of-background responses (see Figs. 3, 17), we used raw maps, comparing responses of each stimulus (+ or –) with the baseline firing rate, defined as the average firing rate at a reverse-correlation delay between 0 and 10 ms, before the visual latency of the cell was reached:

To compare ON and OFF responses (see Fig. 4), we subtracted the background response from the responses of the stimuli (+ or –) that produced center excitation:

Responses were the average throughout the significantly active region (see below, Method 3, area response).

View larger version (13K): [in this window] [in a new window]   Figure 3. Push–pull organization of cone inputs to cone-opponent cells in V1. Most L-ON-center neurons (red dots) were both excited by L+ stimuli and suppressed by L– stimuli, indicated by the location of the responses within quadrant IV (left plot). Similarly, M-ON-center neurons (green spots) were excited by L– and suppressed by L+, indicated by the location of the responses within quadrant II (left plot). Push–pull functional organization was also seen for L-ON and M-ON cells to the M stimulus (middle plot). Blue–yellow neurons showed push–pull to the L and M stimuli less often (most blue and yellow dots in quadrants I and III, L and M plots). Peak responses at the center location were determined using Equation 4 (see Materials and Methods). Error bars show the SD of the firing rate of each neuron. For 2 values describing the distribution of the cells within the plots, see supplemental Table 1A–C.

 

View larger version (24K): [in this window] [in a new window]   Figure 4. ON and OFF responses in cone-opponent neurons reveal asymmetries in responses to light increments and decrements. The response was determined as the average significant response throughout the center receptive-field subregion (method 3; see Materials and Methods). OFF responses were stronger than ON responses (paired t test, p < 0.001). For L-ON cells, the ON response was the response to L+ and the OFF response was that to M–; for M-ON cells, the ON response was the response to M+ and the OFF response was that to L–. For S-ON cells, the ON response was that to S+ and the OFF response was the average response to M– and L–; for (L+M)-ON cells, the ON response was the response to L or M (whichever was stronger), and the OFF response was that to S–. The displacement above the x = y diagonal is significant for the red–green cells and for the total population but not for the blue–yellow cells alone.

 

View larger version (15K): [in this window] [in a new window]   Figure 5. The balance of cone inputs to the receptive-field centers of cone-opponent neurons. Units are spikes per second. Error bars are SDs of the baseline firing rate of each neuron. A, M versus L input to red–green cells. The regression line combining L-ON and M-ON cells showed balance opponency of L and M (M = –0.96 x L; r2 = 0.9; measured using the peak difference method, method 2; see Materials and Methods). Area responses (method 4; Materials and Methods) revealed a subtle difference in cone ratios of the two populations, reflected as a kink about the origin when the average L versus M input to M-ON cells is plotted separately from that to L-ON cells (inset; L-ON cells, M = –1.13 x L; M-ON cells, M = –0.71 x L; regression lines are carried into the adjacent quadrant as dotted lines). This difference reflects the ON/OFF asymmetries shown in Figure 4. B, S versus M input to red–green cells. The majority of red–green neurons were driven by the S-cone stimulus, in a manner that pitted M+S against L. C, D, M versus L input (C) and S versus M input (D) to blue–yellow cells and achromatic cells.

 
Defining receptive-field center and surround
The analysis required identification of the center and the surround subregions of the response maps. Several methods were used to define these regions, to facilitate comparison with other studies. Note that the peak response of the "surround" of LGN cells is centered at the same visual field location as the center. Surround signals of LGN cells are distinguished because they cover a wider area and are generally weaker than the centers (Fig. 1A), but it is only by subtracting the surround response from the center response that the stereotypical center-surround LGN receptive field is generated. The receptive-field structure of cone-opponent cortical neurons is fundamentally different: as we show, the peak response of the surround is centered at a different receptive-field location. To directly assay this spatial offset, we analyzed the unsubtracted maps. To compare the responses with the LGN responses, which are based on difference maps, we also show an analysis of subtracted maps. There is no consensus on whether analysis should be based on peak responses or average responses throughout the receptive field. We therefore analyzed the data in both ways, yielding four separate analyses determined by the following: (1) the peak (i.e., point) unsubtracted response (Eq. 5); (2) the peak difference response (Eq. 3); (3) the average significant (i.e., area) unsubtracted response; and (4) the average significant difference response. In all major conclusions, the methods produced the same results. A comparison of the results of the different methods is shown in supplemental Figure 2 (available at www.jneurosci.org as supplemental material). The specific method used for each analysis is given in each figure legend.

Quantification and analysis was done using programs written in Matlab (Mathworks, Natick, MA).

Defining center and surround maps
The first step of all four methods introduced above was to define, for a given cell, one of the maps (+ or –) of each cone stimulus as the center map (the other is then the surround map). Thus, for a neuron in which all three cone response maps were determined, there were three center maps and three surround maps.

The + or – map containing the earliest excitatory response was defined as the center map. For the cell shown in Figure 2, the L+ map and the M– map would be center maps and the L– and M+ maps would be surround maps; these responses were categorized as ON-center (excitation to the + stimulus) and OFF-center (excitation to the – stimulus) accordingly. It was necessary to use the earliest excitatory response because some neurons gave temporally biphasic responses: the first phase, suppression, corresponding to the onset of the stimulus; the second phase, excitation, corresponding to the cessation of stimulation (see Figs. 11, 16). Center and surround maps were defined by stimulus-onset responses, not rebound responses. The spatial offset of the peak center location within the center maps for a given cell was <0.2° and usually <0.1°.

Determining receptive-field subregions within center and surround maps
Method 1, peak response. The center location was defined as the absolute peak response within the center map. The center location has spatial coordinates, corresponding to degrees of visual angle, and temporal coordinates, corresponding to the optimal reverse-correlation delay. The peak responses at the center location were determined in the separate + and – maps.

The peak surround also has spatial and temporal coordinates. The spatial coordinates were defined by the peak response in the surround map, within ±10 ms of the peak center response. The temporal coordinates were defined by the reverse-correlation delay of the peak response at this location, within the surround map, at any reverse-correlation delay. The temporal coordinates of the peak surround response are therefore not restricted to a fixed temporal window. The response at the spatial and temporal coordinates of the surround were then determined in the separate + and – maps.

Note that the temporal window used to define the spatial location of the surround was not long enough in duration to include excitatory rebound responses to the cessation of the stimulus because the stimuli durations were 34 ms or longer and the temporal window was 10 ms.

When analyzed using this method, a neuron theoretically could produce a response having an identical center and surround location (both spatial and temporal coordinates). This was crucial to be able to evaluate the spatial (see Fig. 9) and temporal (see Fig. 12) offset of center and surround. Figure 7A and supplemental Figure 2 (peak) (available at www.jneurosci.org as supplemental material) show data analyzed using this method.

Method 2, peak difference response. The location of the peak excitatory difference response defined the center location for ON-center maps; the location of the peak suppressive difference response defined the center location for OFF-center maps. The responses were determined at this location in the separate + and – maps.

The surround region was determined in two steps. First, the response at a temporal delay corresponding to the peak center difference response ±10 ms was identified. The spatial coordinates of the peak excitatory difference response within this temporal window defined the location of the surround for OFF-center maps, and the spatial coordinates of the peak suppressive difference response within this temporal window defined the location of the surround for ON-center maps. The surround was then defined as the maximum difference response at this spatial location at any reverse-correlation delay, for the OFF-center maps, and the minimum difference response at this spatial location at any reverse-correlation delay, for the ON-center maps. Responses at the surround location were then measured in the separate + and – maps.

This method constrains the surround response to be in a spatial location offset from the center location. To be significant, the surround response had to be different from the average difference response by >2.5 SDs of the mean difference response. This method is similar to those that use peak responses to cone-isolating gratings in which one phase of the grating is bright and the other dark. See supplemental Figure 2 (peak^d) (available at www.jneurosci.org as supplemental material).

Method 3, area response. The spatial extent of the receptive-field center was defined as that occupied by a response >2.5 SDs above the background for the center maps, at the optimal reverse-correlation delay. To be considered significant, the region had to cover an area of at least 0.2 x 0.2°. The average response of this region was determined in the separate + and – maps.

The temporal coordinates of the surround were the same as for method 1. The spatial coordinates were determined at this delay, as the response >2.5 SDs of the background, at least 0.2 x 0.2°, in the surround maps. The average responses of this region were determined in the separate + and – maps.

This analysis method (and method 4; see below) is most similar to that of Reid and Shapley (2002), who used area responses. See supplemental Figure 2 (area response) (available at www.jneurosci.org as supplemental material).

Method 4, area difference response. The spatial extent of the receptive-field center was defined as that occupied by a difference response >2.5 SDs above the background difference response for the ON-center maps and <2.5 SDs below the background difference response for the OFF-center maps, at the optimal reverse-correlation delay in the difference map. To be considered significant, the region had to cover an area of at least 0.2 x 0.2°. The average response of this region was determined in the separate + and – maps.

The temporal coordinates of the surround were the same as for method 2. The spatial coordinates were determined at this delay, as the difference response <2.5 SDs of the background difference response, at least 0.2 x 0.2° for ON-center maps and as the difference response >2.5 SDs of the background difference response, at least 0.2 x 0.2°, for OFF-center maps. The average responses of this region were determined in the separate + and – maps. See Figure 7B and supplemental Figure 2 (area response^d) (available at www.jneurosci.org as supplemental material).

Determining the spatial extent of the receptive field and the spatial offset of center and surround (see Figs. 8, 9, 14B)
The total area of the receptive-field center and surround were defined as all parts of the response map that produced responses >2.5 SDs of the background, at the optimal reverse-correlation delay of the center and the surround. The extent of the center and surround was then defined as the square root of the total area, to give a sense of the width of the receptive-field subregion. The square root was used for both the center and surround to allow a direct comparison, although the spatial configurations of the center and surround were usually quite different.

The offset (O) of center and surround (for cone maps with significant surrounds), in terms of percentage extent of the receptive-field center, was determined as follows:

where S is the separation (in degrees) between the peak in the center map and the peak in the surround map (see Results, Defining center and surround maps), and RF_c is the extent of the receptive-field center (in degrees). Offsets were only measured in neurons that had significant responses in both center and surround maps, which was the case for the majority of neurons for the L and M maps but less so for the S map. An offset of 50% represents half the width of the receptive-field center.

Rebound responses measurements (see Figs. 16B, 17)
Rebound responses to stimulus cessation were quantified for the receptive-field center; results shown in Figure 16B were determined as the peak response after the zero crossing in the spike-triggered average (STA) histograms, using center locations defined by method 2 (see above). These peak difference responses are given as a fraction of the baseline response. In neurons that showed no significant rebound response, the "rebound" was the response of the center at a reverse-correlation delay set by the average time-to-rebound-peak of the population of neurons whose responses were measured using stimuli of the same duration. Rebound responses shown in Figure 17 were determined using Equation 4 above and show the rebound to + stimuli separately from the rebound to – stimuli. Quadrants II and IV contain neurons that showed opposite response patterns (excitation vs suppression) for + stimuli versus – stimuli; quadrant III contains neurons that were suppressed by both + and – stimuli at long delays.

	Results

Top Abstract Introduction Materials and Methods Results Discussion References

 
The spatial and temporal organization of the cone inputs to cone-opponent cells in V1 was determined using small patches of cone-isolating light presented under neutral-adapting conditions. The cone contrast of the L stimulus was the same as the contrast of the M stimulus, so the relative strength of L and M inputs can be directly determined. The S contrast was higher.

All neurons had receptive-field eccentricities between 2 and 8° and were recorded in alert macaques trained to fixate a spot presented in the center of a computer monitor. Each neuron encountered was first hand mapped with patches and bars of cone-isolating stimuli, flashed at 1 Hz, to determine the location of the receptive-field, assay orientation selectivity and screen for cone opponency (see Materials and Methods). This screening procedure did not presuppose discrete categories of cone-opponent neurons: to pass the screen, neurons simply needed to show excitation to one or two cone types (and not all three). More than 1152 neurons were screened; 165 cells were studied further. From this, we make a conservative estimate that the number of cone-opponent cells in V1 is 10%. Complete maps of the L- and M-cone inputs were obtained for 83 neurons; L, M, and S maps were obtained for 75 neurons. Sixty-six neurons showed dominant opponency between L and M cones (with variable S-cone input), and 17 showed dominant S-cone input, pitted against some combination of L+M. This population of 83 neurons is the basis for the present study (Tables 1–3). Throughout the paper, we will refer to the L versus M cells (with their variable S-cone input) as red–green and the S versus L+M cells as blue–yellow. This is not intended to indicate that these neurons encode specific hues but hopefully makes the text easier to read while following traditional naming convention.

View this table: [in this window] [in a new window]   Table 1. Types of L/M opponent cells

 

View this table: [in this window] [in a new window]   Table 2. Types of S-(L + M) opponent cells

 

View this table: [in this window] [in a new window]   Table 3. Summary of cone-opponent cells

 
Push–pull responses of cone-opponent neurons
Figure 2 (top) shows the spatial distribution of the responses of a red-ON-center neuron in primary visual cortex to cone-isolating stimuli. The temporal pattern, also called the STA histogram or temporal weighting function, is shown at the bottom of Figure 2. The spatial maps show the response at the optimal reverse-correlation delay of the center, 93 ms (the complete spatiotemporal movie of this neuron is given in supplemental Fig. 1, available at www.jneurosci.org as supplemental material). The maps show non-uniform response distributions, indicating that the receptive field was not homogeneous. To quantify the response, we identified different regions of the receptive field as center and surround. For this neuron, this was easy: the neuron gave a stronger response to one of each +/– pair of stimuli. These maps, L+, M–, and S–, were therefore defined as center maps, and the location of the peak of the responses within these maps defined the receptive-field center. The location of the center was identical in all center maps. The remaining maps, L–, M+, and S+, were defined as surround maps. The peak response within these maps defined the receptive-field surround. The time course of the response at the center location and surround location was then extracted from the spike train to produce the spike-triggered averages shown in Figure 2 (bottom).

The neuron whose responses are shown in Figure 2 was chromatically opponent. At the location of the receptive-field center, the cell was excited by L+, a stimulus that looks bright red, and suppressed by M+, a stimulus that looks bright green. (Figure 1A gives an approximate idea of what the stimuli colors looked like.) This chromatic opponency was not defined by luminance: the center was also excited by M–, a stimulus that looks dark red, and suppressed by L–, a stimulus that looks dark green. In both cases, the suppression was evident as a reduction in baseline firing, which can be seen in the response maps as a black region at the location of the center in the M+ and L– maps. This suppression is also clear in the STA histograms: the black traces in the M+ and L– STAs dip below baseline.

Previous studies have not analyzed separately the responses to the two states of a given cone-isolating stimulus, called + and –, or bright and dark. If the cortex were entirely linear in the way it combined inputs from the two states, it would be superfluous to explicitly distinguish them. However, there is evidence that V1 neurons, especially blue–yellow cells (Horwitz et al., 2005), do not sum signals linearly; moreover, psychophysical research indicates that sensitivity is not equal for light increments and light decrements (see Results, Responses to light increments and decrements).

Distinguishing the responses to the two states allowed us to investigate two issues: first, are the responses of cone-opponent neurons push–pull? Second, do cone-opponent neurons show asymmetries between the ON and OFF responses? Push–pull is a term used to describe the functional organization of simple cells in cat primary visual cortex: the light-excitatory, or ON, subregion is both excited by white and inhibited by dark; the dark-excitatory, or OFF, subregion is both excited by dark and inhibited by light (Ferster, 1988; Hirsch et al., 1998). Each component of a single subregion may be attributed to a distinct retinal process, relayed to the cortex by the LGN. For example, the ON subregion of a cortical simple cell may receive excitation from ON-center retinal ganglion cells and inhibition from OFF-center retinal ganglion cells. We sought to evaluate whether cone-opponent neurons in primary visual cortex are analogously push–pull: is the receptive-field center of a red-ON cell both excited by L+, an ON stimulus, and also suppressed by L–, an OFF stimulus? Although extracellular recording does not definitively characterize the mechanism of suppression, intracellular recordings of cat simple cells have shown that extracellular suppression is produced by inhibition (Hirsch et al., 1998).

The neuron shown in Figure 2 was indeed push–pull: L+ caused excitation at the location of the receptive-field center; L– caused suppression, evident as a reduction of baseline firing rate. This push–pull was also present in the M and S maps: at the receptive-field center, the neuron was not only suppressed by M+ and S+ but also excited by M– and S–.

In Figure 3, we analyze the push–pull organization of all of the color-opponent cells recorded by plotting the degree of excitation versus suppression to each state of each cone class. Suppression is not always apparent with extracellular recording when the baseline firing rate is low, but we found evidence for suppression in the majority of cone-opponent neurons we recorded. Push-pull responses, which were calculated using Equation 4 (see Materials and Methods), are indicated below the x-axis or to the left of the y-axis (quadrants II and IV). In Figure 3 and throughout this paper, we distinguish different types of neurons, using a naming convention that depends on the identity of the cone whose + stimulus elicited the strongest center excitatory response. The neuron shown in Figure 2 is called an L-ON-center cell and is depicted by a red dot in Figure 3.

The red–green cells in Figure 3 were distributed differently from the blue–yellow cells in the L and M plots but not in the S plot: red–green cells were more frequently located in quadrants II or IV in the L and M plots, indicating push–pull opponency to + and – stimuli; blue–yellow cells were more frequently located in quadrants I or III, indicating a lack of push–pull (see Fig. 15). The distribution was evaluated using a ² test. In the L plot, 17 red–green cells are located in quadrants I or III and 49 in quadrants II or IV; 12 blue–yellow cells are located in quadrants I and III and 4 in quadrants II or IV. This yields a ² value of 13.66, and a probability of arising by chance p < 0.001 (supplemental Table 1A–C, available at www.jneurosci.org as supplemental material). The ² value of the distribution of red–green and blue–yellow cells for the M plot is 9.4 (p < 0.01). These differences in the distribution of red–green versus blue–yellow cells were not found for responses to the S stimulus, in which both populations of cells showed the same distribution (S plot, ² = 0.23; p < 1.0). Neurons that do not show inverted responses (opponency) to + and – stimuli are not push–pull and nonlinear. Therefore, this result shows that on average (1) blue–yellow neurons gave nonlinear responses to L and M stimuli and (2) this nonlinearity quantitatively distinguished red–green from blue–yellow cells.

The suppression shown by many neurons to some stimuli resulted in a near complete reduction in background firing rate (Figs. 2, 3), leading one to suspect that the stimulus was capable of more suppression than could by measured extracellularly. To overcome this rectification nonlinearity in simple cells, the full extent of the response is often inferred by subtracting the responses to opposite-contrast stimuli (Ferster, 1994). The assumption that the excitatory response to a – stimulus is the same magnitude, but opposite sign, as the suppression to a + stimulus has been shown to be valid for LGN cells in the cat (Martinez et al., 2005). We made this assumption, too, to represent the responses succinctly by difference maps (Fig. 2) (see Materials and Methods, Eq. 3), which allow a direct quantitative calculation of the cone weights. Note that a positive difference can result from any combination of real excitation in the + map and suppression in the – map; similarly, a negative difference can result from any combination of real suppression in the + map and excitation in the – map.

Responses to light increments and decrements
Light decrements (OFF signals) are both more detectable and more readily discriminated than light increments (Walraven, 1977; Krauskopf, 1980; Whittle, 1986; Bowen et al., 1989; Kremers et al., 1993; Chichilnisky and Wandell, 1996), although the subjective contrast of both, as a function of the background luminance, is essentially symmetric (Burkhardt et al., 1984; Whittle, 1986). The higher salience of decrements is thought to reflect an adjustment by the visual system for the fact that, in the natural world, "the range of increment intensities exceeds the range of decrement intensities (increments may be arbitrarily intense while decrements are bounded by zero)" (Chichilnisky and Wandell, 1996). Asymmetries of ON and OFF responses have been described in parasol retinal ganglion cells (Chichilnisky and Kalmar, 2002) and cone-opponent midget retinal ganglion cells (E. J. Chichilnisky, unpublished observations), but asymmetries in cone-opponent neurons in primary visual cortex have not been explored.

We were able to address this issue because the cone-isolating stimuli generate a luminance signal: the M+ and L+ stimuli are increments (ON) and the M– and L– stimuli are decrements (OFF). Luminance values for each of the stimuli are given in Materials and Methods. L-ON and M-ON cells were each excited by one state of the L stimulus and the opposite state of the M stimulus: each neuron received both ON excitation and OFF excitation (Fig. 3). For example, the L-ON cell in Figure 2 was excited by both L+ and M–. We investigated whether the ON and OFF inputs to these neurons were balanced for spatial extent and response magnitude. In general, the ON receptive fields for L versus M cells were slightly larger than the OFF receptive fields (OFF_{rf_size} = 0.93° x ON_{rf_size}; paired t test, p < 0.003), a finding that might not be surprising given that the increments (+ stimuli) had a slightly higher Weber contrast than the decrements (– stimuli; see Materials and Methods). Next we tested for differences in the area response magnitude to increments versus decrements (method 4; see Materials and Methods). OFF area response was stronger than ON area response, indicated by the displacement of the data in Figure 4 (red and green dots) above the x = y diagonal (OFF_{area_response} = 1.2 x ON_{area_response}; r² = 0.54; paired t test, p < 0.001). This bias may be even stronger than that shown here because the Weber contrast was lower for OFF stimuli than for ON stimuli.

Blue–yellow neurons were also excited by ON and OFF stimuli, although they did not show consistent push–pull: S-ON cells were excited by S+ and also by L– or M– (Fig. 3). The blue–yellow cells did not, however, show a significant bias in area response magnitude (Fig. 4, blue and yellow dots), which may be attributed to the higher variability of the responses of this population of neurons; additionally, the S-cone stimulus does not have the same luminance relationship as the M and L stimuli (S– is a luminance increment and S+ is a decrement).

How does the magnitude of the asymmetry shown in Figure 4 compare with psychophysical findings? Because the psychophysical effects vary from person to person and depend on the size of the increment/decrement, the stimulus duration, the brightness of the background, and whether the stimuli are colored (Boynton et al., 1964; Chichilnisky and Wandell, 1996), a direct comparison is difficult. However, consider the finding of Bowen et al. (1989) using sawtooth stimuli that there is "a sensitivity advantage of 0.1 log unit for decremental stimuli at medium to low temporal frequencies (13 Hz)." By interpolating the results of Whittle (1986), who used step stimuli, we infer that a luminance decrement of –5 cd/m² on a background of 25 cd/m² (corresponding to the – stimuli of the L and M cones) would need to be an additional 0.6 cd/m² dimmer to be discriminated. Conversely, a +5 cd/m² increment would need to be an additional 4 cd/m² to be discriminated. The difference (4 – 0.6) is 14% of the background, or 68% of the step size. By comparison, decrements of the red–green cells elicited an average response 20% greater than the response to the increments.

Absolute and relative strength of cone inputs to cone-opponent cells
Figure 5, A and C, quantifies the center response to L and M stimuli for the population of cone-opponent neurons. As a population, the peak responses (Fig. 5A) and area responses (Fig. 5A, inset) of red–green neurons showed balanced L and M input. The area responses suggest that M-ON cells and L-ON cells have slightly different ratios of L versus M input, reflected as a kink about the origin in the regression lines (Fig. 5A, inset); this difference is attributed to the asymmetry in ON and OFF responses (Fig. 4). By comparison, the whole population of type I cells in the LGN receive balanced, equal and opposite, L and M input (Derrington et al., 1984)^a; however, when considered separately, the two kinds of red LGN cells (L-ON and M-OFF) show slightly different ratios of L and M input, as do the two kinds of green LGN cells (M-ON, L-OFF) (Reid and Shapley, 2002).

The majority of red–green neurons we studied also responded to the S-cone stimulus; these responses aligned with responses to the M stimulus in all but two neurons (Figs. 2, 5B). Because the M and S cones are often loosely referred to as green and blue cones and because the combination of green and blue light appears cyan, we suggested calling these cells red–cyan (Conway, 2001). The coupling of M and S cones in opposition to L-cone inputs in V1 neurons is supported by other studies (Wachtler et al., 2003) (G. D. Horwitz, E. J. Chichilnisky, and T. D. Albright, unpublished observations).

Unlike the red–green neurons, the blue–yellow neurons showed considerable variation in the magnitude of the input from L and M cones (Fig. 5C,D), with some neurons showing virtually no M input and others little or no L input. The balance of cone inputs to blue–yellow type II cells in the LGN has not been thoroughly investigated, although the few cases studied show grossly unequal ratios of L and M input (Reid and Shapley, 2002) (but see Chichilnisky and Baylor, 1999).

Spatially opponent receptive-field structure in cone-opponent neurons
In addition to chromatic opponency, each cone input usually showed spatial opponency, resulting in opponent responses (excitation vs suppression) to stimulation in adjacent regions of visual space. A total of 82% of cone-opponent neurons showed spatial opponency in at least one of the chromatically opponent systems (either the red or the green of red–green cells; either the blue or the yellow of blue–yellow cells); 58% showed spatial structure in both opponent systems (Table 3). Spatial opponency was clear in both the raw maps and the difference maps, which have the advantage of subtracting out the common baseline response, leaving a white, zero-difference background (Fig. 2). Against this white background, the L map of the cell in Figure 2 shows excitation in the center and suppression in the surround, depicted as a red region surrounded by a blue crescent. The cell also showed spatial opponency in the other cone system: the M-cone (and S-cone) map shows suppression in the center surrounded by excitation. Like the receptive-field centers, the surrounds of all the cone maps for this cell coincide spatially.

Spatial structure in the difference maps comes about because of the complementary spatial distributions of the responses to the two states of a given cone stimulus. The spatial distribution of one state predicts the spatial distribution of the other: the two responses interdigitate. The interdigitation of + and – responses of a single cone map is analogous to ON and OFF subregions of simple cells, to light and dark stimuli (Hubel and Wiesel, 1962; Conway and Livingstone, 2003; Tsao et al., 2003), although the structure of the responses in cone-opponent neurons was not luminance dependent. The M+ and L+ stimuli are both brighter than the background, yet elicit inverse patterns of responses in each receptive-field subregion.

The combined spatial and chromatic opponency of the neuron shown in Figure 2 earns it the distinction double-opponent. Figure 6 shows the difference maps for 11 more neurons, which are grouped into L versus M+S cells (red–green cells) and S versus L+M cells (blue–yellow cells). Also shown are the cone maps for a luminance simple cell (Fig. 6, Achrom.). Most of the cone maps show regions of both excitation (red) and suppression (blue), indicating spatial opponency; the spatial opponency between the maps of most of the cone-opponent cells are complementary, which was not the case for the achromatic luminance simple cell. For example, the peak response of cell 2 was excited by L and suppressed by M and S (shown in the response maps as a red blob at the center of the L map but a blue blob at this location in the M and S maps); the surround (to the right of center) was excited by M and S and also suppressed by L. The achromatic cell showed the same spatial pattern of opponency in each map (it received negligible S input, as shown by the noisy S-cone map).

View larger version (52K): [in this window] [in a new window]   Figure 6. Spatially opponent receptive fields of cone-opponent neurons in V1. Each row represents the cone difference maps for a single cell. Red indicates excitation, and blue indicates suppression. With the exception of the achromatic cell (bottom), all cells showed chromatic opponency. The L versus M+S cells (red–green cells) showed opposite patterns of response (excitation or suppression) between the L and M maps (and often S maps) within a given spatial location. The S versus L+M cells (blue–yellow cells) showed opposite patterns of response within a given spatial location between S and L or M maps. In addition, each single map shows spatial opponency, indicated by areas of excitation and areas of suppression. Spatial scale of the overlying grid is 0.75°. The reverse-correlation delays (in milliseconds) of the response maps are given in the bottom left corner of the L map for each cell. Scale bar is in spikes per second. Complete spatiotemporal maps are shown for cells 2, 3, 4, 8, 10, and 11 in supplemental Figure 7 (available at www.jneurosci.org as supplemental material), supplemental Figure 3 (available at www.jneurosci.org as supplemental material), Figure 10, supplemental Figure 6 (available at www.jneurosci.org as supplemental material), supplemental Figure 5 (available at www.jneurosci.org as supplemental material), and Figure 15, respectively.

 
These cone maps reveal several features about the cone-opponent neurons. Few of the double-opponent neurons had completely annular surround responses; rather than being doughnut shaped, as hypothesized previously (Hubel and Wiesel, 1968; Gouras, 1974; Michael, 1978; Livingstone and Hubel, 1984; Billock, 1991), the surrounds were often crescents, or adjacent parallel subfields (Fig. 6, cell 6), much the same shape as the ON and OFF subfields of luminance simple cells. These shapes resemble color basis functions, or filters, generated from Independent Component Analysis of Munsell Color Chips and colored natural scenes (Tailor et al., 2000; Buchsbaum and Bloch, 2002; Doi et al., 2003; Caywood et al., 2004). The asymmetry of the surround also predicts that the neurons will show some orientation preference, and many did as assayed by initial hand mapping (see also Conway, 2001).

For the red–green cells, the surround in one cone map aligned spatially with the surround of the remaining cone maps. For example, the L-OFF surround of cell 2, which sits to the right of the receptive-field center, is in the same spatial location as the M-ON surround. This spatial coincidence was not evident in all blue–yellow cells. Cell 10 (Fig. 6) (complete maps in supplemental Fig. 5, available at www.jneurosci.org as supplemental material) shows a crescent-shaped M-OFF surround that sits below the receptive-field center, but the L-OFF and S-ON surrounds straddle the receptive-field center left and right.

Many cells showed reproducible hot spots within the surround (Fig. 6, M-cone map for cell 4). Moreover, the surrounds showed variable spatial extent and strength, relative to the strength of the centers, with some cells showing surrounds of equal size or larger than the center response (Fig. 6, L-cone map, cell 6) and others showing barely significant or insignificant surround responses (for a summary, see Table 3). There was no correlation between the magnitude of the response (using any of the methods to define center and surround) and the spatial extent of the centers or surrounds (data not shown). Twenty-four percent of cone-opponent neurons had significant surround responses in only one map (Table 3), a pattern of response that has been defined as -opponent (Livingstone and Hubel, 1984). Three neurons showed a variation of -opponency in which the surround, but not the center, of one of the cone maps was significant. This is seen in cell 5 (Fig. 6), which showed no significant S– response but a strong S+ response located exclusively in the surround region defined by the L and M maps. These patterns suggest that double-opponent cells are wired by sampling separately the ON and OFF signals from each of the cone types.

In Figure 7, we plot the magnitude of the peak center response against the peak surround response, separately for each cone type. The peak responses of the surrounds were one-fifth as strong as the centers, for each of the cone maps (Fig. 7A), but this analysis gives an underestimate of the strength of the surrounds of spatial-opponent cells because cone-opponent neurons, having insignificant surrounds and which form a discrete population of neurons (see Results, Cone-opponent neurons without spatial opponency), were included. Excluding cells with insignificant surrounds gives surround estimates ranging from approximately one-third as strong as the centers (if method 2, peak difference responses, is used to define center and surround) to approximately two-thirds as strong as the centers (if method 4, area difference responses, is used to define center and surround). A comparison of the results of the different methods is shown in supplemental Figure 2 (available at www.jneurosci.org as supplemental material), in which cells with insignificant surrounds are plotted along the x-axes for methods 2 and 4.

View larger version (23K): [in this window] [in a new window]   Figure 7. Quantification of center/surround opponency and cone weights. A, Spatial opponency of the L-cone maps (left), M-cone maps (center), and S-cone maps (right). The peak response (method 1; see Materials and Methods) was used to measure the strength of the center versus the strength of the surround. Quadrants II and IV contain neurons that had spatial opponency, indicated by the schematic "Mexican hat" receptive-field profile. Units are spikes per second. Lsurround = –0.21 x Lcenter, r2 = 0.2; Msurround = –0.22 x Mcenter, r2 = 0.3; Ssurround = –0.32 x 32 Scenter, r2= –0.1. Insets show LGN responses to similar stimuli; axes limits are 40 spikes/(second x C), where C is the contrast of the stimulus (data from Reid and Shapley, 2002). Spatial opponency is evident using any method of indexing the center and surround (supplemental Fig. 2, available at www.jneurosci.org as supplemental material). B, Cone weights. The sum of the absolute values of the responses (method 4; see Materials and Methods) to each cone stimulus were normalized to a total of 1. The normalized values of the L- and M-cone inputs are plotted, preserving the sign of the response. The S response is indicated by the distance from the diagonal lines toward the origin. Neurons in which S-cone responses were not measured are depicted outside the diagonal lines; neurons with negligible S-cone input lie on the diagonal lines. A red-ON cell having balanced, opponent L and M input, with no S input, would be plotted as L = 0.5, M = –0.5, at the center of the diagonal in quadrant IV.

 
Regardless of the method of analysis, the center versus surround responses of all cone-opponent cells are distributed predominantly in quadrants II and IV for all cone maps, indicating that cone-opponent neurons on average show spatial opponency (Fig. 7A) (supplemental Fig. 2, available at www.jneurosci.org as supplemental material). All methods of analysis also show that the centers and surrounds have opposite chromatic tuning: the L-ON cells are almost exclusively distributed in quadrant IV in the L cone map and almost exclusively in quadrant II in the M cone map; the M-ON cells show the reverse pattern. Spatial chromatic opponency is underscored in cone-weight plots (Fig. 7B): in the center plot, red-ON cells fall in quadrant II and green-ON cells in quadrant IV, but in the surround plot, red-ON cells fall in quadrant IV and green-ON cells in quadrant II. Blue-ON and yellow-ON cells also swap quadrants. Cone-weight plots also support the conclusion that the red–green neurons receive approximately balanced L and M input: both center and surround responses cluster around the midpoints on the diagonals (Fig. 7B) (see also Fig. 5A). This clustering is similar to that of type I cells plotted in cone-weight space (Derrington et al., 1984). Blue–yellow neurons, conversely, did not on average receive equal input from L and M cones or show a consistent ratio of L and M input, as indicated by the scatter of blue and yellow dots in quadrants I and III (Fig. 7B) (see also Fig. 5C,D). That the cone weights for red–green and blue–yellow cells are nonoverlapping is not an artifact of the selection criteria because neurons were not selected for the type of opponency or for the cone weights they receive.

The spatial organization of the responses of V1 cone-opponent neurons can be directly compared with the responses of LGN type I cells. Type I cells fall almost exclusively in quadrants I and III (Fig. 7A, insets) (data from Reid and Shapley, 2002); in contrast, cortical cone-opponent cells fall almost exclusively in quadrants II and IV. This shows that the striate cortex makes a spatial transformation of the cone signals, which produces neurons capable of spatial color comparisons.

Although weaker in magnitude than the centers, the surrounds occupied a relatively large region of visual space (L and M inputs, average of 0.7° wide), over half the extent of the receptive-field center (L and M average of 1° wide) (Fig. 8A). Center and surround sizes were determined by the entire significantly activated region of the response maps. It may be surprising that the surround is actually smaller than the center (compare the sizes of the red and blue blobs in each panel of Fig. 6), unlike LGN neurons in which the surround is usually larger than the center. This is partly attributed to a difference in definition of surround required by cortical neurons (see Materials and Methods, Defining receptive-field center and surround). Moreover, some cortical cone-opponent neurons lack a surround (see Results, Cone-opponent neurons without spatial opponency) yet were included in the analysis of Figure 8A; even excluding these cells, however, the surrounds were still smaller than the centers (Fig. 8B). There was a weak correlation between the spatial extent of the receptive-field center and the spatial extent of the surround for the L and M plots but not the S plot (L plot, r² = 0.35, p = 0.005; M plot, r² = 0.26, p = 0.04; S plot, r² = 0.11, p = 0.6); there was no correlation between the magnitude of the response and the spatial extent of the receptive field for either center or surround (analysis not shown). Despite the smaller surround-to-center size ratio of cortical cells compared with LGN cells, both the centers and surrounds of cortical cells were still larger than centers and surrounds of LGN cells, whose average receptive-field width at comparable eccentricity is 0.24–0.5° (Reid and Shapley, 2002