Abstract
Most mammals are dichromats, having short-wavelength-sensitive (S) and middle-wavelength-sensitive (M) cones. Smaller terrestrial species commonly express a dual gradient in opsins, with M opsin concentrated superiorly and declining inferiorly, and vice-versa for S opsin. Some ganglion cells in these retinas combine S- and M-cone inputs antagonistically, but no direct evidence links this physiological opponency with morphology; nor is it known whether opponency varies with the opsin gradients. By recording from >3000 ganglion cells in guinea pig, we identified small numbers of color-opponent cells. Chromatic properties were characterized by responses to monochromatic spots and/or spots produced by mixtures of two primary lights. Superior retina contained cells with strong S+/M− and M+/S− opponency, whereas inferior retina contained cells with weak opponency. In superior retina, the opponent cells had well-balanced M and S weights, while in inferior retina the weights were unbalanced, with the M weights being much weaker. The M and S components of opponent cell receptive fields had approximately the same diameter. Opponent cells injected with Lucifer yellow restricted their dendrites to the ON stratum of the inner plexiform layer and provided sufficient membrane area (∼2.1 × 104 μm2) to collect ∼3.9 × 103 bipolar synapses. Two bistratified cells studied were nonopponent. The apparent decline in S/M opponency from superior to inferior retina is consistent with the dual gradient and a model where photoreceptor signals in both superior and inferior retina are processed by the same postreceptoral circuitry.
Introduction
Most mammals have short-wavelength-sensitive (S) cones and one or two types of middle/long-wavelength-sensitive (M/L) cones (with some deviations, Jacobs, 1993; Peichl et al., 2001) (for review, see Peichl, 2005; also Arrese et al., 2006; Cowing et al., 2008). Inputs from S and M/L cones combine antagonistically to form the basis for blue-yellow color vision. In trichromatic and dichromatic primate retinas, S-ON:M/L-OFF opponency is conveyed by both small and large sparse bistratified ganglion cells, while M/L-ON:S-OFF opponency is conveyed by large sparse monostratified and S-midget ganglion cells (Dacey, 1993; Dacey and Lee, 1994; Ghosh et al., 1997; Calkins et al., 1998; Ghosh and Grünert, 1999; Silveira et al., 1999; Dacey and Packer, 2003; Klug et al., 2003; Schein et al., 2004; Dacey et al., 2005; Martin et al., 2007). In other dichromatic mammals, physiological measurements reveal S/M opponent ganglion cells (Cleland and Levick, 1974; Caldwell and Daw, 1978; Ekesten et al., 2000; Grivich et al., 2005). These measurements, however, do not link the physiologically opponent cells to morphologically defined cell types. Based on anatomy it has been suggested that for rabbit, bistratified ganglion cells such as the G3 might correspond to the bistratified opponent ganglion cells in primate and thus carry S-ON:M-OFF signals (Rockhill et al., 2002). To clarify the nature of color vision in a nonprimate dichromatic mammal, we studied the morphology and physiology of opponent ganglion cells in guinea pig.
In guinea pig superior retina, each cone expresses either S or M opsin exclusively, and S cones comprise ∼6% of all cones (Röhlich et al., 1994; Yin et al., 2006). This expression pattern is typical of many dichromatic mammals (Jacobs, 1993; Szél et al., 1996; Ahnelt and Kolb, 2000; Szél et al., 2000; Peichl, 2005; Schiviz et al., 2008). Therefore, we asked whether the S+/M− ganglion cells bistratify (as in primate), sending dendrites to both ON and OFF layers of the IPL. The inferior retina of guinea pig expresses mostly S opsin, and the small fraction of M opsin is found in cones that coexpress both opsins (Röhlich et al., 1994; Yin et al., 2006). Therefore, we asked whether this region contains color-opponent ganglion cells at all, and if so, whether the strength of opponency reflects the gradient of opsin expression.
We recorded from ganglion cells via loose-patch electrode and tested for color opponency. Opponency proved rare, occurring in <1% of all cells tested. We characterized opponent cell spectral properties by measuring the response to flashed monochromatic spots and/or spots produced by mixtures of primary lights. Then we replaced the recording electrode with a sharp electrode and filled the cell with Lucifer yellow to reveal the dendritic arbor. We found both S+/M− and M+/S− cells; both types were monostratified with flat arbors restricted to the IPL's ON stratum. Opponent cells were found in both superior and inferior retina. Across the retina, the strength of opponency changed with the change in proportion of S and M opsins.
Materials and Methods
Tissue preparation.
We used albino guinea pigs weighing 300–700 g. An animal was overdosed with anesthetic (ketamine 100 mg/kg; xylazine 20 mg/kg and pentobarbital 50 mg/kg) and both eyes were removed and hemisected. The eyecup was cut radially and flattened onto a membrane filter with ganglion cells up. This preparation was placed into a chamber on a microscope stage (BX50WI, Olympus America) and superfused with oxygenated (95% O2 and 5% CO2) Ames medium (Sigma) containing sodium bicarbonate (1.9 g/L) and glucose (0.8 g/L) at ∼5–7 ml/min. The chamber temperature was maintained between 34 and 37°C. The retina was dark adapted for 0.5 h before measurements were made.
To help find particular retinal regions during the experiment, holes were cut in the membrane filter (Fig. 1A). When the tissue was flattened onto the filter, the small central hole (0.8–1 mm diameter) was positioned over the optic disk and the filter was oriented so that the other two elongated holes were located above or below the optic disk. The middle of each region was aligned along the superior-inferior axis of the tissue through the optic disk (3 mm superior or inferior to the optic disk). Recordings were made mainly from within these regions. In the first experimental stage of the study (defined below), recordings were also made from regions in superior retina closer to the optic disk. Retinal locations of the recording sites were measured relative to the optic disk.
S/M opponent cells are present in superior and inferior retina. A. Right, Regions searched. Left, Locations of all opponent cells relative to optic disk. Cells in superior and inferior retina are plotted as solid and open circles respectively. The S+/M− and M+/S− opponent cells are plotted in blue and green. Three cells identified in superior retina lie closer to optic disk (<2 mm) than the shaded circles shown in the right panel. One was encountered during an earlier experiment and the other two were encountered in the first experimental stage, when we periodically broadened the search region (see Materials and Methods). The arrow indicates where two S+/M− cells from superior retina plot to the same location. One M+/S− cell from superior retina which was found in the first stage is not plotted, because data on its location were not available. B, Right, An S+/M− opponent cell from superior retina in loose-patch configuration (DIC optics). Left, Soma sizes (estimated along the minor axis) versus maintained discharge rates in steady light. The arrow indicates where two M+/S− cells from superior retina plot to the same location. Same cells as in A. C, Proportions of opponent cells encountered in the third and fourth experimental stages (see Materials and Methods). Of all the ganglion cells patched that met the anatomical screening criteria, the percentage of cells that also met the physiological screening criterion are shown as the combined height of light and dark bars. The light bars represent cells that were nonopponent, or responded weakly or not at all to contrast. The dark bars represent opponent cells. The numbers of cells for the light bars and the dark bars from left to right were as follows: 5 and 1, 17 and 3, and 6 and 3. Total numbers of cells that met the anatomical criterion for each population are indicated above the bars. Not indicated in the figure is the fact that we also patched 15 cells in inferior retina in stage four, none of which were opponent. One of these 15 cells met our physiological screening criteria.
Electrical recording.
Ganglion cells were identified for recording by observing the tissue through the infrared differential interference contrast (DIC) optics of the microscope (Fig. 1B, right panel) with a 60× (0.9 NA) objective. Cells were recorded extracellularly with loose patch electrodes (tip impedance = 2–4 MΩ). Responses were amplified (NeuroData IR-283, Cygnus Technology), digitized at a 5 kHz sampling rate with 12-bit precision (DigiData 1200, Molecular Devices), and stored for later analysis (AxoScope, Molecular Devices). Spikes were identified using a threshold detection method. Spike rates were calculated for ∼15 ms bins.
Searching for S/M opponent cells.
To search for S/M opponent cells, we went through four experimental stages for a total of 66 experiments, adopting different search strategies in each stage. In the first stage (n = 17 experiments), we began by patching ganglion cells in superior retina and testing every cell for opponency (see below for a description of the opponency test). In targeting cells for patching we excluded brisk-transient ganglion cells, which could be identified because of their large soma size, but otherwise patched all cells encountered. Later analysis revealed that the patched cells generally had soma sizes between 10 and 15 μm, with a few as large as 18 μm (Fig. 1B, left) (soma size was typically estimated along the minor axis). We patched >145 cells in the first stage and encountered three opponent cells. Before the first stage, in experiments conducted as part of a separate study, one additional opponent cell was encountered near the superior region of the transition zone, where some M cones coexpress S opsin (Röhlich et al., 1994). We include that cell here (Fig. 1A, blue solid circle closest to middle line in left panel; morphology in the third panel from top in Fig. 7A).
For the opponent cells encountered in the first stage, we observed certain anatomical and physiological properties in common. For example, they tended to have a medium soma size and high maintained discharge rate in light. We took those properties as predictive of whether a ganglion cell would be opponent in stage two of the experiments.
In stage two (n = 26 experiments), we searched for ganglion cells in superior retina that matched the properties described above and only took the time to test for opponency when a cell matched these qualitative selection criteria. As we encountered more opponent cells (n = 9 more cells in stage two, out of >1085 cells patched), we gradually refined the qualitative selection criteria. These observations led us to develop quantitative screening criteria for identifying S/M opponent cells in superior retina: (1) medium soma size, ranging between 10 and 14 μm and (2) high maintained discharge rates to a uniform background (main background, see below; >15 spikes/s), and also moderate to high maintained rates in the dark (> 5 spikes/s). All opponent cells encountered in this study satisfied these criteria (Fig. 1B, left panel; maintained rates in dark are not shown, as they were not measured for some of the opponent cells in the first stage). Criterion (1) was based on anatomy and was met by a high percentage of ganglion cells (visual observation; Do-Nascimento et al., 1991). Criterion (2) was based on physiology and was met only by a very small percentage of ganglion cells (see Fig. 1C). In our experience, most ganglion cells had a maintained rate to the main background of ≤ 10 spikes/s. We also studied >41 cells from inferior retina in stage 2, without encountering an opponent cell.
In stages one and two we were focused on finding opponent cells and developing a sense of their general properties. We did not in these exploratory stages keep explicit records of every nonopponent cell patched. The reported numbers (145 for superior retina in stage one, 1085 for superior retina in stage 2, 41 for inferior retina for stage 2) are the number of cells explicitly documented, and represent a lower bound on the number of cells actually patched.
In the third stage (n = 18 experiments), we applied the screening criteria developed in stage two for experiments conducted in both superior and inferior retina. We only patched ganglion cells with appropriate soma size. Of all the cells patched in superior retina (327 cells, 1 opponent), we tried to study in detail those that had appropriate maintained rates and were opponent (Fig. 1C, left and middle bars). Of all the cells patched in inferior retina (1392 cells, 3 opponent), we tried to study every cell with appropriate maintained rates, except those that responded weakly or not at all to contrast (Fig. 1C, middle bars). We occasionally studied cells whose maintained rates fell below the cutoffs we established and none of those cells were opponent.
The final experimental stage (stage four) was designed to check whether the physiological screening criteria might have excluded a major population of opponent cells. In this stage (n = 6 experiments), we tested every cell patched in superior retina for opponency, and studied all of the opponent cells encountered (n = 3, out of 122 cells patched). All of the opponent cells identified in this stage satisfied the criteria established in stage two (Fig. 1C, right bars). We also patched 15 cells in inferior retina during this stage, none of which were opponent.
Characterizing chromatic response properties of ganglion cells.
Chromatic response properties of ganglion cells were characterized by recording responses to spots of different spectral composition. “Flicker experiments” measured responses to spots on a spatially uniform background. These spots were produced by modulating mixtures of two primary lights. This method provided an efficient way to assess whether a cell was opponent, and to estimate the spatial extent of its receptive field. We used this method, as described above, to identify S/M opponent cells. “Flash experiments” were designed to provide a quantitative characterization of a cell's opponency. These experiments measured responses to flashed monochromatic spots on a spatially uniform background. This method was slow but allowed direct characterization of the spectral properties of a ganglion cell. We used this method not only to verify that a cell was opponent, but also to measure the relative strengths of M-cone, S-cone, and rod inputs (M and S and rod “photopigment weights”).
We used the same two optical systems as in our previous study (Yin et al., 2006). Briefly, (1) “LCD system,” light from an LCD projector (PowerLite 730c, Epson America) produced images that subtended ∼2.4 × 3.2 mm on the retina, with each pixel corresponding to ∼3.1 × 3.1 μm. (2) “Lamp system,” light from a xenon lamp (HLX 64642, Osram), collimated and fed through a tunable narrowband filter (VariSpec, Cambridge Research & Instrumentation) and an adjustable aperture produced a uniform spot of monochromatic light (∼10 nm bandwidth) on the tissue. The temporal profile of the spot was controlled by a mechanical shutter (VS25S2T1 shutter and 122-BP controller, Uniblitz). The outputs of the two optical systems were combined through a beamsplitter and delivered to the tissue through the camera port on the microscope and the objective (4×). Full details of the configurations of the two systems, and on spectral characterization and calibration, are reported in Yin et al. (2006). As the original lamp system did not produce sufficient light intensity at 430 nm for this study we expanded the lamp system with another channel, which was fed through a narrowband interference filter of 430 nm.
The flicker experiments used the LCD system alone to modulate a spot (∼100–1600 μm diameter) around the background at 2–4 Hz (sinusoidal temporal waveform). The main background produced nearly equal isomerization rates in M and S cones (4.06–4.17 and 4.01–4.12 log10 Rh*/photoreceptor/s) and a somewhat higher rate in rods (4.65–4.76 log10 Rh*/photoreceptor/s). The exact background intensities are noted in the text where necessary. In the earlier experiment (described above, one opponent cell), a CRT instead of the LCD projector was used in the optical system. The background used in that experiment produced isomerization rates of 3.99 (M), 3.53 (S), 4.44 (Rod) log10 Rh*/photoreceptor/s.
Reported isomerization rates were computed as described by Yin et al. (2006), based on parameters estimated for intact guinea pig retina. In our in vitro preparation, the basic orientation of the receptors with respect to the incident light was the same as for intact retina, but detailed differences in orientation could produce deviations between computed and actual isomerization rates. We did not correct for any such differences. Note that our analyses of cell properties are not sensitive to a change in scaling of the computed isomerization rates, since these were based on contrast.
Ganglion cells were adapted continuously throughout the recording. To screen a cell for opponency, a sequence of flickering spots modulated along M-cone and S-cone isolating directions as well as in M+S and M−S color directions was repeated (typical stimulus parameters: ∼800 μm in diameter, temporal frequencies of 2 or 4 Hz, duration of 8 s for each direction). The spike traces and voltage pulses marking the start and end of each modulation cycle were visualized using AxoScope (Molecular Devices). From the relative response amplitudes and phases to modulations along the four color directions, we could determine whether a cell was opponent. The characteristics of S/M opponent cells were that (1) they responded to modulations along both M- and S-cone isolating directions, but with opposite response polarities and (2) they responded strongly to modulation along the M−S opponent direction but weakly to the M+S nonopponent direction (Fig. 2A). Nonopponent ganglion cells responded with the same polarity to modulation along both M-cone and S-cone isolating directions, and more strongly to modulation along the M+S nonopponent direction than the M−S opponent direction. In some experiments, to test cells for opponency as efficiently as possible, the flickering stimulus was presented directly through the 60× objective instead of 4× objective. This reduced the size of the stimulus. We did not perform a separate spectral calibration for the 60× objective.
S/M opponent cells characterized by flickering and flashed monochromatic spots. A, Opponency test using flickering spots. All panels show the responses (PSTH, dark lines) of an S+/M− opponent cell in superior retina to flickering spots (diameter 800 μm; temporal frequency ∼4 Hz) modulated along M-cone and S-cone isolating directions, as well as in M+S and M−S color directions. Temporal waveforms of the four modulations (green and blue lines) are shown in the insets of each plot. The M and S cone contrasts, [CM, CS], of the four modulations were (from left to right): [35%, 0%], [0%, 35%], [25%, 25%], and [24%, −24%]. B, Spatial structure of the receptive field. Same cell as in A. The plot shows the signed response amplitude (F1 ± SEM) to modulations in the same four color directions used in A, measured in another block of trials, where the spot diameter was varied from 200 to 1600 μm. C, Response to flashed monochromatic spots. All panels show the responses (PSTHs, in blue) of an M+/S− opponent cell in superior retina to flashed monochromatic spots (diameter ∼850 μm; flash duration 300 ms). This cell responded to flashed monochromatic spots of short and middle/long wavelengths with opposite response polarities; the morphology of this cell shown in the third panel from the top in Figure 7B. Stimulus onset was at 100 ms. PSTHs were fitted (red lines) with a template (Appendix A). PSTHs and fits were shifted vertically to align the fitted baseline rate with zero spikes/s.
We recorded the spike responses to flickering spots of ganglion cells identified as opponent by the screening procedure (typical stimulus parameters: 800 μm in diameter, temporal frequency of 4 Hz, duration of 32 s for each direction). The number of contrasts presented for each color direction (M-isolating, S-isolating, M+S, M−S) varied from cell to cell. For some cells, we also recorded area summation data in the form of responses to flickering spots of diameters between 100 and 1600 μm (Figs. 2B, 9; supplemental Fig. 1, available at www.jneurosci.org as supplemental material). For quantitative analysis, the response was taken as the signed fundamental response amplitude (F1) of the peristimulus time histogram (PSTH), with the sign of the response indicating whether the fundamental was best described as in-phase or out-of-phase with the stimulus modulation (after accounting for response delay). Responses during the first and last stimulus cycles were excluded from the analysis.
The flash experiments used monochromatic spots from the lamp system on a background produced by the LCD system. The main background was the same that used for the flicker experiments, and cells were adapted continuously throughout the recording. A sequence of flashed monochromatic spots of wavelengths between 430 and 620 nm was presented. We refer to the presentation of a series of flashed monochromatic spots at one wavelength as a “trial,” and to the presentation of a sequence of trials for different wavelengths as a “block”; stimulus parameters: 800–900 μm diameter, flash duration of 300 ms at 1 s intervals (500 ms for two cells), 12–32 repetitions per trial, and 19–30 trials per block.
The spike traces and voltage traces marking the flash onset and offset were visualized using AxoScope (Molecular Devices) and fed to an audio-speaker for real time monitoring. From the response polarities to monochromatic flashes at short, middle, and long wavelengths we could decide whether a cell was opponent. The characteristics of S/M opponent cells were that (1) they had spectral neutral points, that is a wavelength range where they responded only weakly to monochromatic flashes and (2) they responded to monochromatic flashes of wavelengths longer or shorter than spectral neutral points with opposite response polarities.
As the PSTHs were noisy, we derived fits to each individual PSTH using a template method and extracted the response amplitudes from the fits (Appendix A). Responses were taken as the peak of spike rate increment from baseline at either the flash onset or offset (Fig. 2C). We report photopigment weights derived from spectral measurements for seven opponent ganglion cells and one possibly opponent ganglion cell (see Results) (Fig. 6). We report spectral neutral points for five additional cells where full spectral characterization was not possible (see Results) (Fig. 3). Stimulus parameters for these cells were the same, except that one was studied with a 2000-μm-diameter spot.
Spectral neutral points of S/M opponent cells vary across retina. Spectral neutral points are indicated by horizontal lines. At wavelengths within each horizontal line, the cell responded with indeterminate response polarity; at the wavelengths indicated by the lines' ends the cell responded with clear and opposite polarities. For cells in superior retina (solid lines), data were combined across multiple blocks. For cells in inferior retina (dashed lines), measurements from the individual blocks are shown separately. Measurements for S+/M− and M+/S− opponent cells are presented below and above the dashed black line, and are also distinguished by being shown with black and gray bars respectively. For S+/M− opponent cell 1, the second block of trials did not reveal a spectral neutral point within the wavelength range we measured (triangle). Photopigment weights for a subset of cells (S+/M−, index 1, 2, and 3; M+/S−, index 1, 2, 5, and 6) are shown in Figures 4C, 4F, 5C, and 6.
Estimating photopigment weights.
From the responses to monochromatic spots, we could estimate the relative strength of M, S, and rod inputs to a ganglion cell. Appendix B describes the linear–nonlinear (LNL) model we used to estimate contrast-response nonlinearities, action spectra, and photopigment weights. The responses of all nine ganglion cells where we had sufficient flash data to perform this analysis were well-described by the LNL model. For these cells, the maximal cone or rod contrasts presented (across wavelengths) were <120% for three cells, 150% for three cells, 200% for one cell, 250% for one cell and 350% for one cell.
Dye staining, immunostaining, and quantitative morphology.
After recording, some ganglion cells were penetrated with a sharp electrode and stained with Lucifer yellow (3%, in 0.1 m LiCl solution; Invitrogen) to visualize their morphology. To verify that the cell stained was the one recorded, DIC views of the cell soma were captured before the loose patch electrode was released and after the sharp electrode was positioned near the soma. After staining, the tissue was fixed in 4% paraformaldehyde in 0.1 m phosphate buffer (PB) for 1 h in dark and store in PB in 4 C° before immunostaining.
After a brief wash in PB, the retina was peeled from the pigment epithelium and incubated in blocking buffer (10% normal goat serum, 5% Triton X-100 in PB with 5% sucrose) for 1 h to reduce nonspecific binding. We reacted the retina with primary antibodies to Lucifer yellow to enhance the staining of the dendrites of the ganglion cell, and to S opsin to reveal the overlaying S-cone mosaic. For Lucifer yellow we used rabbit anti-Lucifer yellow (Invitrogen) at a dilution of 1:500 and for S opsin we used rabbit anti-S opsin (Millipore Bioscience Research Reagents) at a dilution of 1:600. The reaction was carried out initially at room temperature for 2 h, then at 4 C° for 35 h. After several washes in PB with 5% sucrose (SPB), we reacted the retina with secondary antibodies to visualize the primary antibody staining (4 h, at room temperature). Both the dendrites of the ganglion cell and S opsin were tagged with FITC or CY3 (Invitrogen). After several washes in SPB, we reacted the tissue with SYTO 61 or SYTO 13 nucleic acid stains (Invitrogen) to visualize the amacrine cell layer (ACL) and the ganglion cell layer (GCL), which outlined the boundaries of the IPL. Finally after several washes in SPB, we mounted the retina on a glass slide with the ganglion cell side up and covered it with mounting medium (Vectashield, Vector Laboratories) and a coverslip.
The tissue was examined using a confocal microscope (Leica TCSNT, Leica) with 10× (0.3 NA), 40× (1.25 NA), and 100× (1.4 NA) objectives. Image stacks were obtained at two z-depths, with one covering IPL, including the stained ganglion cell and another one covering the photoreceptor layer, including the S-cone outer segments. To quantify the dendritic morphology of a ganglion cell: (1) we estimated the diameter of the dendritic field from the tangential view at 20× (10× objective with 2× magnification; 500 × 500 μm). We drew a polygonal contour connecting the most distal dendritic tips and then fit the polygon with an ellipse using the NIH Image software (http://rsb.info.nih.gov/nih-image/). The estimated diameter was the average length of the major and minor axis of the fitted ellipse. (2) We measured the length of the ON or OFF dendrites from the tangential image stack at 40× (250 × 250 μm) by tracing all dendrites in 3D using the Volocity software (Improvision). The density of the ON dendritic processes was calculated by dividing the length of ON dendrites by the area of the polygon containing all the ON dendrites in the same view. This density is one of the parameters that proved effective for morphological classification of ganglion cells in a cluster analysis reported by Kong et al. (2005). (3) We measured the relative depths of distal dendrites within IPL at several locations from the tangential image stack at 40× and sometimes 100× (250 × 250 μm), also using the Volocity software. The depth of dendritic stratification was expressed as the percentage of the total depth of the IPL, starting from 0% at the boundary with the amacrine cell layer and ending with 100% at the boundary with the ganglion cell layer. To generate radial views, we resliced a section of ∼250 × 60 μm within the tangential image stack using the NIH Image software (see Figs. 7A–C, 8A).
We quantified the dendritic membrane area of opponent ganglion cells using a procedure similar to that described by Xu et al. (2008). First, we divided the tangential image stack obtained at 40× magnification into subregions: a circular region of 30 μm diameter centered around the cell soma and a series of outward expanding concentric rings with equal widths of 25 μm. All of the image stacks analyzed had x- and y-resolutions of 0.24 μm and a z-resolution of 0.49 μm. One image stack (that shown in Fig. 7A, the third panel from top) was excluded because its z-resolution was larger than the diameters of typical distal dendrites. Second, after smoothing the images with 3 × 3 median filter, we used the Volocity software to recognize and reconstruct the dendrites contained in each of the ring subregions. Finally, we measured the dendritic membrane area from the reconstructed dendrites and calculated the dendritic membrane density by dividing the dendritic membrane area by the retinal area of corresponding ring.
Results
Responses to modulations
Figure 2A shows responses from an S+/M− cell from superior retina to 4 Hz modulations of an 800 μm spot. The responses clearly reveal the opponent nature of this cell. Responses to modulation along M- or S-cone isolating directions had opposite response polarities. The response to modulation along the M−S opponent direction was strong, but that to the M+S nonopponent modulation was weak. M+/S− opponent cells showed a similar response pattern, but with opposite response polarities.
Figure 2B shows the dependence of the responses on spot size for the same cell. For M-cone, S-cone, and M−S modulations, response amplitude increased with spot diameter and then remained fairly constant, indicating a lack of spatial antagonism for these modulation directions. For the M+S modulation response first increased with diameter and then decreased, indicating a spatially antagonistic receptive field structure for this modulation direction. These measurements are qualitatively consistent with an S-ON center opposed by a slightly larger M-OFF surround (see Fig. 9A, second panel from top on the left). Figure 9 and supplemental Figure 1 (available at www.jneurosci.org as supplemental material) provide the full set of area summation measurements we made for opponent cells, along with quantitative model fits that provide estimates of their receptive field structure (see detailed discussion below; Appendix C). We checked for, and did not observe, changes in response phase with stimulus size (data not shown).
Spectral neutral points, action spectra, derived photopigment weights, and receptive field structure for cells identified as opponent by their responses to modulations are presented below. These include S+/M− and M+/S− cells in superior retina and S+/M− cells in inferior retina. We also encountered one M−/S+ cell in inferior retina, as classified by its responses to flickering M-isolating, S-isolating, M+S, and M−S modulations. We were not able to hold this cell long enough to characterize further its physiology, nor were we able to fill it and obtain its morphology.
Spectral neutral points of S/M opponent cells
The response polarities of S/M opponent cells to flashed monochromatic spots varied with wavelength, depending on the relative strengths of the excitatory input from one cone type and the inhibitory input from another cone type (Figs. 2C, 4A,D, 5A). At the spectral neutral points, where the excitatory and inhibitory inputs balanced each other, response polarities were indeterminate (Fig. 2C, PSTH for 480 nm). For S+/M− opponent cells, at wavelengths shorter than the spectral neutral points, the response polarities were ON, reflecting the dominant S-ON (S+) input. At wavelengths longer than the spectral neutral points, the response polarities were OFF, reflecting the dominant M-OFF (M−) input. The M+/S− opponent cells showed the reverse pattern.
Opponent cells in superior retina showed relatively balanced M/S antagonism. A, Response waveforms (PSTHs, in blue) of an S+/M− opponent cell from superior retina in response to flashed monochromatic spots (diameter ∼850 μm; flash duration 300 ms; flash onset at 100 ms in the plots). Same format as Figure 2C. The flash intensities were 3.99 × 104, 5.68 × 104, 3.23 × 104, and 2.40 × 105 quanta/μm2/s (440, 460, 530, and 610 nm). The response amplitudes were measured from the fitted waveforms at either flash onset or offset, depending on the polarity of the response (see Appendix A). The background produced 4.11 (M), 4.06 (S), and 4.70 (Rod) log10 Rh*/photoreceptor/s. Data shown are from the first block of trials for this cell. B and C also show data from this block of trials for the same cell and background intensity. The morphology of this cell is shown in the second panel from top in Figure 7A. B, Flash contrast-response function. Response amplitudes of either ON (solid gray circles) or OFF (solid blue circles) responses to monochromatic flashes of wavelengths from 420 to 620 nm were fitted with two static nonlinearities (solid red line for ON responses; dashed red line for OFF responses) as described in Appendix B. The x-axis of the plot is the weighted input contrast, Cinput, derived from the model fit (see Eq. B1 in Appendix B). C, Action spectrum of an S+/M− opponent cell in superior retina. The same data as in B are plotted as a spectral sensitivity. The solid and dashed red lines represent the magnitude of the weighted sum of M, S, and rod contrast sensitivities as determined by the model fit. The solid line represents the wavelength region where there was an ON response, while the dashed line represents the wavelength region where there was an OFF response. The relative photopigment weights used to produce the action spectrum were M: −0.41, S: 0.58, and rod: −0.02 and are shown by the horizontal bars in the inset (magnitudes of the M, S, and rod weights are represented by green, blue, and gray bars). Sensitivity dropped significantly between 460 and 480 nm, indicating the spectral neutral point (see also Fig. 3). The overall shape of the action spectrum here and in F resembles the guinea pig increment-threshold spectral sensitivity measured behaviorally by Jacobs and Deegan (1994). The vertical location of the data points on the spectral plots has been corrected for variation in stimulus intensity across wavelength and for the shape of the cells' ON and OFF static nonlinearities (see Appendix B) (Yin et al., 2006). D, Response waveforms (PSTHs, in blue) of an M+/S− opponent cell from superior retina to flashed monochromatic spots (diameter ∼850 μm; flash duration 300 ms). Same format and background intensity as A. The flash intensities were 3.42 × 104, 5.08 × 104, 2.96 × 104, and 2.12 × 105 quanta/μm2/s (440, 460, 530, and 610 nm). Data shown are from the first block of trials for this cell. E and F also show data from this block of trials. The morphology of this cell is shown in the second panel from top in Figure 7B. E, Flash contrast-response function. Same format as B. F, Action spectrum of an S+/M− opponent cell from superior retina. Same format as C. The relative weights were M: 0.42, S: −0.51, and rod: 0.07 and are shown in the inset. Because the response amplitude to the flash at 470 nm was zero, no data point is plotted for that wavelength.
Opponent cells in inferior retina were dominated by S input. A, Response waveforms (PSTHs, in blue) of an S+/M− opponent cell from inferior retina to flashed monochromatic spots (diameter of ∼800 μm; flash duration of 300 ms). Same format and background intensity as Figure 4A. The flash intensities were 5.69 × 104, 1.57 × 105, 8.74 × 104, 5.14 × 104, 4.62 × 104, 1.66 × 105, and 3.74 × 105 quanta/μm2/s (440, 460, 480, 500, 530, 590, and 610 nm). Data shown are from the first block of trials for this cell. B and C also show data from this block of trials. The morphology of this cell is shown in the top panel of Figure 7C. B, Flash contrast-response function. Same format as Figure 4B. For one data point (one of four measurements made at 520 nm) at low weight input contrast, the response polarity of the data differed from that of the model prediction. This point is shown in the plot as having a negative response value. Note that this sign reversal is of small magnitude and occurs near the spectral neutral point of the cell. C, Action spectrum of an S+/M− opponent cell from inferior retina. Same format as Figure 4C. The relative weights were M: −0.01, S: 0.99, and rod: 0 and are shown in the inset. That data point indicated in B as having a negative response is not shown on this logarithmic plot.
All S/M opponent cells in superior retina had similar spectral neutral points, for both S+/M− (n = 4) and M+/S− (n = 6) types (Fig. 3). The neutral points were similar to those determined behaviorally for guinea pig (Jacobs and Deegan, 1994).
S+/M− opponent cells in inferior retina (n = 2) (Fig. 3) had spectral neutral points at longer wavelengths than those in superior retina. For one cell, the opponency which had been apparent in the flicker screening measurements and in the first block of flash trials was absent in the second block of trials (Fig. 3, triangle). That the spectral neutral points occur at longer wavelengths for cells from inferior retina is consistent with a decline in the relative strength of M input for S/M opponent cells from superior to inferior retina. This decline was confirmed by the quantitative characterization of photopigment weights presented in the next section.
Photopigment weights of S/M opponent cells
We measured photopigment contributions to S/M opponent cells from their responses to monochromatic flashes (see Appendices A and B). Action spectra of opponent cells across retina were well fit with weighted sums of M, S, and rod spectral sensitivities (Figs. 4C,F, 5C). The nonlinearities in the contrast responses of each cell were also well described (Figs. 4B,E, 5B). The overall shape of the action spectra for S/M opponent cells in superior retina resembled the increment-threshold spectral sensitivity functions measured behaviorally against an achromatic background [Jacobs and Deegan (1994), their Figs. 6, 7]. This suggests that S/M opponent cells in superior retina subserve guinea pig color vision.
Photopigment weights of S/M opponent cells vary across retina and differ from those for horizontal cells and other types of ganglion cells. Photopigment weights of all ganglion cells measured in this study are shown. The x- and y-axes indicate the relative M and S weights. The position of each point along a third axis, extending from the origin along the diagonals at −45°, 45°, and 135°, indicates the absolute value of the rod weight. Cells from superior and inferior retina are plotted as solid and open symbols respectively. Opponent cells are plotted as solid and open circles in the second and fourth quadrants, with symbol colors of blue and green for S+/M− and M+/S− types. One possible S+/M− opponent cell is plotted as an open diamond in the second quadrant (see Results). For cells having opposite signs for S and rod weights, the edge of the corresponding symbol is shown in red, otherwise in black. All measurements were made at similar background intensities (see Materials and Methods). For cells studied for multiple blocks of trials, photopigment weights were estimated from each block; and the corresponding data points are connected with lines. For cases where there was a substantial shift in weights across blocks, the arrow on the connecting line points to the second block. Opponent cells in superior retina showed relatively balanced antagonism between M and S inputs, but in inferior retina S input dominated. For comparison, the mean photopigment weights (± SD) of horizontal cells in superior and inferior retina at a similar background intensity are plotted as solid and open triangles [Yin et al. (2006), their Fig. 3]. The photopigment weights of brisk-transient ganglion cells (see Fig. 11) are similar to those of horizontal cells (Yin et al., 2006).
Dendrites of S/M opponent cells were located in the IPL's ON stratum. A–C, Top panels, Tangential and radial views of all well stained opponent cells. In the radial view, the IPL is the region between the GCL and ACL (bands of red staining on the left and right). In the figure, the horizontal dimension of each radial view is stretched by 300% relative to the vertical dimension. Based on the staining pattern (upper and middle panels) or the distance between the dendritic layer of a cell (bottom panel) to the right edge of its soma (∼3 μm), all cells stratified close to the GCL, that is in the ON stratum. If available, the depths of dendritic stratification are labeled on the tangential views (bottom-right corners) and also visualized in the bottom. For one M+/S− opponent cell (the third panel from top in B), a small proportion of dendrites looped into the OFF stratum (at a depth of 15%; tangential view, arrow), then terminated back into the ON stratum (10% of the total length of the ON dendrites; data not shown). S+/M− opponent cells in inferior retina had many very short lateral branches in their dendrites (one example shown as arrowhead in the third panel from top in C, which was enlarged from the panel above). We tested seven of the eight cells for opponency using both flickering spots and flashed monochromatic spots (see also Figs. 2C, 3, and 4A,D and 5A), and further estimated the photopigment weights of six cells from their flash data (see Figs. 4C,F and 5C and 6). One S+/M− opponent cell (the third panel from top in A) was encountered during a different experiment and was only tested for opponency using flickering spots (see Materials and Methods). For another S+/M− opponent cell (top panel in A) we did not have sufficient data from the flash experiments to allow quantitative estimation of photopigment weights. A–C, Bottom panels, Summary of dendritic stratification for opponent cells. The means of the depth of dendritic stratification were similar across opponent types and retinal locations (A, 72%; B, 77%; and C, 77%). Error bars in red show the SD.
The action spectra can also be interpreted in terms of the S and M photopigment weights. S/M opponent cells in superior retina had balanced inputs from S and M cones, for both S+/M− and M+/S− types. Correspondingly, in Figure 6, the points representing photopigment weights of S/M opponent cells fell in the vicinity of the −45° or 135° diagonals. Because S cones are much sparser in superior retina than M cones, the balanced weights indicate that signals from S cones are selectively amplified.
S/M opponent cells in inferior retina received unbalanced inputs from M and S photopigments, with the relative strength of the M input was much weaker than that of the S input (Fig. 5). Correspondingly, in Figure 6, the data points representing the photopigments weights for S/M opponent cells in inferior retina were close to the y-axis, where the M weight equals zero. The cell that failed to show a spectral neutral point in the second block of trials (Fig. 3, triangle) still showed slight opponency for this block when the data were analyzed in terms of photopigment weights (Fig. 6, arrow).
One additional cell studied from inferior retina is worth comment. This cell met the anatomical screening criteria, but its maintained discharge rate fell below the screening cutoff we established. Still, we recorded from it, to learn whether it represented a novel cell type. This cell was monostratified in the ON stratum of the IPL. It was not opponent according to its responses to flickering stimuli, nor did it have a spectral neutral point in response to monochromatic flashed stimuli. For these reasons, we have not grouped this cell with the two S+/M− cells we identified in inferior retina. The quantitative analysis, however, did reveal weak M and S weights of opposite sign (Fig. 6, diamond) and it is possible that this cell represents a third exemplar of an S+/M− opponent cell.
As shown in Figure 6, the S weights of four S/M opponent cells decreased substantially during the course of experiments (connecting lines with arrows). It is possible that cone opsins were bleached under continuous light adaptation, although our retinal preparation was capable of regenerating bleached photopigments (see Materials and Methods). The finding in mouse that S cones are more susceptible to bleaching than M cones could then potentially explain our observation (Nikonov et al., 2006).
Photopigment weights of S/M opponent cells across retina (Fig. 6, solid and open circles) differed from those of horizontal cells measured at a similar background intensity [Yin et al. (2006), their Fig. 3], as well as from those of brisk-transient ganglion cells (see Fig. 11) [Yin et al. (2006), their Fig. 6], in two aspects. First, opponent cells in superior retina received much larger S input than horizontal cells or brisk-transient ganglion cells in superior retina, and for the opponent cells the S weight (absolute value: 0.37 ± 0.17 SD, n = 5 for both S+/M− and M+/S− types) (Fig. 11) was much larger than the percentage of S cones in superior retina (∼6%). This suggests that S input to the opponent cells is processed by a retinal pathway with high relative gain. Presumably this pathway involves S-ON bipolar cells (Mariani, 1984; Kouyama and Marshak, 1992; Calkins et al., 1998; Haverkamp et al., 2005; Li and DeVries, 2006; MacNeil and Gaul, 2008). Second, opponent cells generally received smaller amounts of rod input than horizontal cells from the same retinal regions. This would be explained if rod signals reach S/M opponent cells through both of the S and M antagonistic pathways that drive the ganglion cell and thus tend to cancel. For all opponent cells where we could quantify the inputs, the rod weight was either negligible (absolute value ≤7%) or else had the same sign as the S- cone input (Fig. 6).
Morphology of S/M opponent cells
Based on the morphologies of well-stained opponent cells (Fig. 7), S/M opponent cells across retina were monostratified, either completely or primarily, in the IPL's ON stratum. This held for both S+/M− (superior retina: n = 3; inferior retina: n = 2) and for M+/S− types (superior retina: n = 3) (Fig. 7). Opponent cells from both superior and inferior retina were similar in terms of the parameters we used to quantify cell morphology: the diameter of the dendritic field, the depth of dendritic stratification and the density of the ON dendritic processes. Values for these parameters are provided in Table 1. There may, however, be small differences in aspects of the morphology that we did not quantify. For example the presence of very short lateral dendritic branches appeared more prominent in the cells from inferior retina (Fig. 7C, the third panel from top). We also obtained incomplete morphologies (data not shown) of five more S/M opponent cells in superior retina (S+/M−: n = 3; M+/S−: n = 2). Four of these cells appeared to stratify in the ON stratum. A fifth appeared to stratify in the OFF stratum, based on staining of the primary and secondary dendrites. But the distal dendrites, which would provide more reliable information about the locus of dendritic stratification, were not stained and prevent us from attaching much significance to this observation.
Morphological and anatomical parameters for opponent ganglion cells
Table 1 also provides the number of S cones labeled within the dendritic field of opponent cells in superior retina, the density of S cones within the dendritic field, and estimates of the S-cone percentage within the dendritic field. The latter were obtained by assuming an overall cone density of 20,000 cones/mm2 (Yin et al., 2006). The mean value for the S-cone percentage (n = 5) was 4.5%, slightly smaller than previous estimates of ∼6% S cones in guinea pig superior retina (Röhlich et al., 1994; Yin et al., 2006). We regard the previous estimates as more reliable because they were obtained in a preparation where both S and M cones were labeled in the same retinal region, and thus not subject to error introduced by local variation in overall cone density.
In tangential view the dendritic pattern of opponent ganglion cells was similar to that of rabbit G9 cells, while their dendritic stratification resembled that of rabbit G10 cells (Rockhill et al., 2002). None of the known rabbit ganglion cells, however, have a morphology that corresponds completely to that of the opponent ganglion cells we studied.
The amount of information that a ganglion cell collects, and the area over which it does so, are key aspects of its function. Both are constrained by the number and distribution of its bipolar synapses. For all ganglion cell types studied to date synaptic density is constant at ∼19 synapses per 100 μm2 of dendritic membrane, independent of cell type (Xu et al., 2008). Therefore, measuring dendritic membrane area for the color-opponent cells provides an indirect estimate the number of bipolar synapses for comparison to other known types. For our cells, the total area of the dendritic membrane was 2.1 ± 0.3 (SD) × 104 μm2 (data available from seven of the eight cells in Fig. 7; see also supplemental Fig. 2, available at www.jneurosci.org as supplemental material). The average peak density of dendritic membrane area was 51 μm2 per 100 μm2 retinal area. Thus assuming the same density and distribution of synapses as for other ganglion types, the color-opponent cell would collect about ∼3900 synapses. This number is similar to those obtained for ON-OFF direction-selective (DS) and local-edge detecting (LE) ganglion cells, but smaller than that for brisk-sustained (BS) ganglion cells and larger than that for brisk-transient (BT) ganglion cells (supplemental Fig. 2, available at www.jneurosci.org as supplemental material) (Xu et al., 2008).
Some bistratified ganglion cells were nonopponent
We stained a few cells that met the anatomical and physiological screening criteria but that were not opponent as revealed by the flickering spot protocol. Two of these were bistratified ganglion cells. Their morphologies are shown in Figure 8A, their physiological responses in Figure 8B,C. The morphology, including the depth of dendritic stratification, differs from that of other well-studied bistratified ganglion cells in guinea pig, such as the ON-OFF direction-selective and ON-OFF ganglion cells (data not shown). A ganglion cell type with similar dendritic stratification, the G3 bistratified ganglion cell, has been reported in rabbit [Rockhill et al. (2002), their Figs. 2, 4]. Indeed, Rockhill et al. (2002) speculated that this bistratified type might be color-opponent. If the bistratified cells we studied are the homolog of the G3, however, this is not the case. Moreover, the fact that the bistratified cells we identified are not opponent also makes it unlikely that they are the homolog of the small bistratified ganglion cell that has been well-studied in primate (Dacey, 1993; Dacey and Lee, 1994; Ghosh et al., 1997; Ghosh and Grünert, 1999; Silveira et al., 1999; Dacey and Packer, 2003; Dacey et al., 2005).
Bistratified ganglion cells did not show explicit opponency. A, Bottom panels, Tangential views of bistratified nonopponent ganglion cells that met both the anatomical and physiological screening criteria. For each cell, tangential views for ON and OFF dendrites as well as a radial are shown from left to right along each row. Top panels show a cell from superior retina and bottom panels show a cell from inferior retina. In the radial views, the IPL is the region between the GCL and ACL (bands of red staining on the left and right). The total length of the OFF dendrites was 23% of the total length of the ON dendrites for the cell from superior retina. The corresponding number for the cell from inferior retina was 78%. A, Bottom, Summary of dendritic stratification for the bistratified ganglion cells. Error bars in red show the SD for each stratum. B, C, Responses to flickering spots (PSTHs, dark lines) of the two bistratified ganglion cells shown in A and B. Neither cell showed opponency. In B, the stimulus diameter was 600 μm and the stimulus was flickered at a temporal frequency of 4 Hz. Modulation directions shown in the plots from left to right were M-cone isolating, S-cone isolating, M+S and M−S. The M and S cone contrasts, [CM, CS], were (from left to right): [20%, 0%], [0%, 20%], [20%, 20%] and [20%, −20%]. In C, the stimulus diameter was 800 μm and the stimulus was flickered at a temporal frequency of 4 Hz. The M and S cone contrasts, [CM, CS], were (from left to right): [30%, 0%], [0%, 30%], [21%, 21%] and [21%, −21%]. In B, the response to M-cone isolating modulation was stronger than that to S-cone isolating modulation, vice versa for C. In both B and C, responses to modulation along the M+S nonopponent direction were stronger than those along the M−S opponent direction.
Receptive field structure
We collected area summation data of the sort shown in Figure 2B for six S+/M− and six M+/S− opponent cells from superior retina and one S+/M− cell from inferior retina (see Fig. 9 and supplemental Fig. 1, available at www.jneurosci.org as supplemental material). These data can be used to infer receptive field properties of the cells, by fitting them with a Difference of Gaussians (DOG) receptive field model, as described in Appendix C. Briefly, we fit each cell with a Gaussian center and surround. Both the center and surround took a weighted sum of S and M inputs; the weights for both center and surround were parameters of the model. Because the area summation data were collected using the LCD system, we did not have direct measurements of the rod contribution at each stimulus size. In fitting the data, we assumed a constant rod-cone coupling strength and fixed the value of this strength at the average value derived from measurements of spectral sensitivity (see Fig. 11 and Table 2). Derivation of this value is discussed below (Functional model; Appendix D). In addition, we used data from each cell collected at one spot size to constrain the static nonlinearity of that cell, as described in Appendix C.
S- and M-cone receptive field components derived from DOG model. A–C, Left panels, One dimensional view of S and M receptive field components (RF profiles; blue for S and green for M), normalized by the peak strength of the dominant input type (S or M). Actual RF profiles are circularly symmetric two-dimensional functions. The RF profiles were derived from the area summation data shown in the right panels (see Appendix C). In A, B, and C, RF profiles for S+/M− opponent cells from superior retina, M+/S− opponent cells from superior retina and an S+/M− opponent cell from inferior retina are shown respectively. (Data and fits for additional S/M opponent cells from superior retina are shown in supplemental Fig. 1, available at www.jneurosci.org as supplemental material.) The spatial extents of the S and M receptive field components are indicated by the horizontal blue and green bars. These extents were calculated as the radius of the circular region that contained 80% of the unsigned volume of each component, where we defined the unsigned volume as the integral of the absolute value of the two-dimensional RF profile. The dendritic field diameter for each cell is shown as a gray horizontal bar. From the top to bottom, the cells in each panel of this figure correspond to those shown in Figure 7 in the following order: the second and the first panels from top in A, the second, the third and the first panels from top in B, and the second panel from top in C. A–C, Right panels, Area summation data with DOG receptive field model fits. Plots in the same format as Figure 2B. The data from Figure 2B are replotted in the second panel from the top in A. The signed response amplitudes (F1) to modulations in four color directions are shown as green and blue crosses (for M- or S-cone isolating directions respectively) and red and green dots (for M+S and M−S color directions, respectively). The DOG model fits are plotted in the same color as the corresponding data points.
Parameters used to predict photopigment weights of retinal neurons
Figure 9 shows area summation data along with the receptive field structure derived from the fits for six cells: two S+/M− opponent cells from superior retina (Fig. 9A), three M+/S− opponent cells from superior retina (Fig. 9B), and one S+/M− opponent cell from inferior retina (Fig. 9C). For these cells, we also had morphological data that allowed us to estimate dendritic field diameter. Supplemental Figure 1 (available at www.jneurosci.org as supplemental material) shows area summation data and receptive field structure for an additional seven cells (four S+/M− and three M+/S− opponent cells from superior retina). For the most part, the area summation curves were well fit by the model. The receptive field structure is plotted as the spatial profile of the S and M components; within the model each of these is a Difference of Gaussians. For 11 of the 13 cells, the receptive field structure had opposed S- and M-cone receptive fields, each with a dominant center Gaussian component. This is the classic type II receptive field structure (Wiesel and Hubel, 1966) exhibited by primate S-cone opponent cells (Dacey and Lee, 1994; Dacey, 1996; Dacey and Packer, 2003; Field et al., 2007).
There are two exceptions to the general characterization above. (1) For one S+/M− cell from superior retina (supplemental Fig. 1A, second panel from the top, available at www.jneurosci.org as supplemental material), the M-cone receptive field showed more marked spatial antagonism than for the other cells. Correspondingly, the area summation data showed that this cell's responses had a nonopponent character for small spot sizes. The S+/M− opponency emerges in the area summation data only for large spot sizes. (2) The S+/M− cell from inferior retina does not appear opponent in fits or in the area summation data (Fig. 9C). We attribute this to drift in the rod input to the cell over time. The area summation data were collected after the monochromatic flash data presented in Figure 3 (S+/M− opponent cell, index 1). As shown in Figure 3, the second block for this cell did not reveal a spectral neutral point, and the cell had a large estimated rod weight for this block (Fig. 6, arrow). If we refit the area summation data for this cell assuming a rod-cone coupling constant consistent with this larger rod weight, the derived M receptive field component is very weakly antagonistic to the derived S component (fit not shown).
We quantified the spatial extent of the S- and M-cone components of the receptive fields for the 11 cells that showed type II receptive fields (that is, for all cells except the two noted in the previous paragraph) by finding the diameter that accounted for 80% of the unsigned volume of each components' receptive field. These diameters are indicated for each cell by horizontal blue and green bars (Fig. 9; supplemental Fig. 1, left panels, available at www.jneurosci.org as supplemental material). Figure 10 shows a scatterplot of the M-cone component diameters versus S-cone component diameters. The spatial extents of the M-cone component of the receptive fields were slightly larger than those of the S-cone components. A similar observation has been reported for primate S-cone opponent cells (Dacey, 2000; Field et al., 2007). For cells where we had morphological data, we also plot the dendritic field diameter (Fig. 9, gray bars in left panels). Dendritic field diameters are smaller than the physiologically characterized diameters, although clearly the exact relation depends on the criterion used to define the physiological diameters. In dichromatic mammals and primates, S-ON cone bipolar cells have long and meandering dendrites and axons (Mariani, 1984; Kouyama and Marshak, 1992; Haverkamp et al., 2005; MacNeil and Gaul, 2008), which could provide a mechanism for opponent ganglion cells to collect input from S cones located in regions adjacent to their dendritic fields. The relationship in the relative size between physiologically measured receptive field and anatomically measured dendritic field that we observed is roughly consistent with that found for other ganglion cell types, such as brisk-transient ganglion cells (Borghuis et al., 2008).
Comparison of M and S receptive field components. Scatter-plot of the spatial extents of M versus S receptive field components. S+/M− and M+/S− opponent cells are plotted in blue and green respectively for 11 opponent cells, all from superior retina. Two of the cells shown in Figure 9 and supplemental Figure 1 (available at www.jneurosci.org as supplemental material) are excluded from this plot, as noted in the text. The values of spatial extents plotted are the lengths of the blue and green horizontal bars for the corresponding cells in the left panels of Figure 9 and supplemental Figure 1, available at www.jneurosci.org as supplemental material. The average ratio of the spatial extents of M to S receptive field component for all 11 cells was 1.22, significantly >1 (two-tailed t test, p < 0.01). In superior retina, this ratio was slightly larger for S+/M− (average ratio 1.35) than for M+/S− cells (average ratio 1.12). This difference was marginally significant (two-tailed t test, p < 0.10).
Functional model
In the context of the evolution of color vision, it is of interest to know whether a change in cone opsin expression can enhance visual performance in the absence of a concomitant change in postreceptoral circuitry (Jacobs et al., 2007; Mancuso et al., 2007). We thus asked, can the differences in cone weights across the superior-inferior opsin gradient, both for the opponent cells reported here and for achromatic cells (Yin et al., 2006), be explained by assuming common postreceptoral circuitry in superior and inferior retina? To address this, we constructed a functional model that predicts the strength of photopigment inputs to each class of cells studied. In outline, the model is characterized by the following parameters: (1) the proportion of S cones in superior retina, which is also the proportion of genuine S cones (Haverkamp et al., 2005) in inferior retina, (2) the relative gain of S-cone signal conveyed to the S/M opponent ganglion cells, (3) the strength of rod coupling to cones, which is assumed to be indiscriminate across all cones (Ahnelt et al., 1990), and (4) the proportion of M pigment coexpressed in nongenuine S cones in inferior retina. Appendix D describes the relation between these parameters and predicted photopigment weights. We used a numerical search procedure to find the parameter values that best predicted the measured weights. The model also incorporates two additional assumptions, which we think are reasonable. These are: (1) rod signals reach ganglion cells only through rod-cone gap junctions at the light intensities of our study (Sharpe and Stockman, 1999; Bloomfield and Dacheux, 2001; Sterling, 2004), and (2) M- and S-opsins drive the phototransduction cascade with equal gain, independent of the cone type in which they are expressed (Nikonov et al., 2006).
The modeling results are provided in Table 1, and the corresponding predictions are plotted with the data in Figure 11. The predicted weights capture all of the broad features of the measurements and are in reasonable quantitative agreement. Given the simplicity of the model, we take this agreement as indicating that the model provides a plausible account of the physiology. The key conclusion is that a single change in opsin expression pattern across the retina, without any other change in circuitry, is consistent with all of the observations reported here and in our previous study.
Functional model predictions. Modeling results in the same format as Figure 6. Solid and open circles show the average photopigment weights for opponent cells from superior and inferior retina (for superior retina, weights were averaged across S+/M− and M+/S− opponent cells by inverting the sign of the weights from the M+/S− opponent cells). Triangles (± SD) show averaged data from horizontal cells in superior (solid) and inferior retina (open), with two triangles for each retinal location representing results from two measurement protocols (Yin et al., 2006). Squares (± SD) show averaged data from brisk-transient ganglion cells (superior retina, solid; inferior retina, open). The red plusses show the photopigment weights derived from the model fit (see Table 1).
Discussion
Properties of S/M opponent cells in guinea pig
Our study is the first to jointly characterize the physiology and morphology of S/M opponent cells in a nonprimate mammalian species. In guinea pig superior retina, we find that S/M opponency is conveyed by both S+/M− and M+/S− opponent cells, consistent with earlier physiological reports for ground squirrel, cat, rabbit, mouse, and guinea pig (Michael, 1966; Cleland and Levick, 1974; Caldwell and Daw, 1978; Jacobs et al., 1981; Ekesten et al., 2000; Hemmi et al., 2002; Grivich et al., 2005). Both opponent cell types share similar morphology, with a flat dendritic arbor that stratifies in the IPL's ON stratum. Both types also have similar chromatic properties, in that S and M photopigment weights are approximately balanced. S+/M− and M+/S− cells in superior retina have type II receptive fields: the S and M components of the receptive field are of approximately equal size, with the M component being slightly larger for both S+/M− and M+/S− types.
In inferior retina, we found two S+/M− cells and one M+/S− cell. The opponency of the S+/M− cells was weak, with the S input dominating the M input. We were not able to obtain sufficient physiological data on the M+/S− cell to characterize the strength of its opponency.
Proportion of S/M opponent cells
In all, we tested >1679 cells in superior retina for opponency in stages 1–4. Of these, 16 were opponent. Taking the number of cells measured at the lower bound of 1679, an upper bound on the observed proportion of opponent cells is 1%. The 99% Bayesian confidence interval (Lee, 1989) on this proportion, computed using an uninformative uniform prior, is 0.4% to 1.7%. This proportion and confidence interval are exclusive of brisk transient cells and, for many of the measurements, of cells that did not meet the screening criteria.
Borghuis et al. (2008) show that in guinea pig, as in many other mammals (review by Wassle, 2004), the coverage factor for ganglion cells, defined as the product of dendritic field area and cell density, is roughly constant (between 3 and 4) for cell types in which dendritic tips extend to the neighboring soma. For these cell types, the coverage factor is independent of cell size and eccentricity. The mean dendritic field diameter for the opponent cells was 344 ± 39 (SD) μm (n = 8). If we assume a coverage factor of 3.5, we obtain a density estimate for each of the S+/M− and M+/S− types of ∼40 cells/mm2. Given the total density for all ganglion cells of ∼1500 cells/mm2 (similar to Do-Nascimento et al., 1991; Kao and Sterling, 2006), the S+M− and M+/S− types would each represent ∼2–3% of the ganglion cell population. This is higher than the ∼1% proportion we observed in superior retina for the two types combined. Either our screening procedure missed some opponent cells or the coverage factor for the opponent cells is lower than for other cell types.
In inferior retina, 3 of 1407 cells that met our screening criteria showed opponency, for a proportion of 0.2% (99% confidence interval 0% to 0.7%). The lower hit rate in inferior retina could reflect that fact that S/M opponency there is weak. In particular, cells with stronger than average rod input might not show overt opponency when screened. None of the opponent cells in inferior retina had balanced S and M photopigment weights. If there are strongly opponent cells in inferior retina that meet the screening criteria, they are very rare.
Ekesten and Gouras (2005) reported that 2% of ganglion cells they studied (physiology only) in mouse retina were opponent, somewhat higher than in our study. Given differences in species and sampling methods, however, it seems plausible that our data are tapping the same population of opponent cells they sampled. An earlier study from the same group (Ekesten et al., 2000) reported a higher percentage of S/M opponent cells in both superior and inferior retina in mouse.
Relation to melanopsin-expressing ganglion cells
In primate, “giant” monostratified melanopsin-expressing ganglion cells are color opponent (Dacey et al., 2005). Although we did not explicitly investigate whether the S/M opponent cells we studied express melanopsin, we think it is unlikely that they do. First, the spectral sensitivity of melanopsin peaks at around 480 nm (Berson et al., 2002; Dacey et al., 2005). However, as shown in Figure 3, the spectral neutral point of the S/M opponent cells in superior retina is at around that wavelength. On this point, we acknowledge that our procedures would not maximize the contribution of the sluggish temporal response characteristic of melanopsin (Berson et al., 2002; Dacey et al., 2005). Second, melanopsin-expressing ganglion cells are known to encode absolute irradiance (Dacey et al., 2005). However, for 11 S/M opponent cells, we measured maintained rates in both light and dark. These were similar (average of 29 and 24 spikes/s, respectively). Third, Dacey et al. (2005) report that the “giant” melanopsin-expressing ganglion cells were (L+M)+/S−. Here we report both M+/S− and S+/M− types, with similar dendritic stratifications.
Implications for retinal circuitry
Our data constrain the circuit configurations for S+/M− and M+/S− ganglion opponent cells. First, the excitatory signals should flow through ON bipolar cells. Second, in superior retina at least, the circuit should allow selective amplification of signals from S cones to produce the observed balanced opponency. Finally, the circuitry should account for the fact that the S- and M-cone components of the receptive field have approximately the same diameter.
The S+/M− cells probably receive their input from S-ON bipolar cells, whose existence in mammals is well-established (Mariani, 1984; Kouyama and Marshak, 1992; Calkins et al., 1998; Herr et al., 2003; Schein et al., 2004; Haverkamp et al., 2005; Li and DeVries, 2006; MacNeil and Gaul, 2008). The M-antagonism might arise at the outer retina, as has been suggested for primate (Dacey et al., 1996; Dacey, 2000; Field et al., 2007; Packer et al., 2007; Davenport et al., 2008), and/or from amacrine inhibition driven by an M-ON bipolar cell (Li and DeVries, 2006). In either case, the S-ON bipolar cells could selectively amplify the signals from S cones.
For the M+/S− cells, the excitatory M input could arise directly from an M-ON bipolar cell (Li and DeVries, 2006). In this case, the most likely source of the amplified inhibitory S input would be through an amacrine cell that connected to the S-ON bipolar cell. Alternately, if the S-ON bipolar cell itself carries antagonistic M signals, the M+/S− opponent cell could simply receive the majority of its input from the S-ON bipolar cell through a sign-inverting amacrine cell. This latter possibility would make the S opponent circuit similar to that of the mammalian rod bipolar pathway, where a single type of rod ON bipolar cell drives both ON and OFF ganglion cells (Sterling, 1983; Müller et al., 1988; Bloomfield and Dacheux, 2001). Similar possibilities have been suggested for the inhibitory S input to M/L-ON:S-OFF ganglion cells in primate (Dacey and Packer, 2003).
Weak color vision in guinea pig inferior retina?
Our results indicate that S/M opponency is considerably weaker in inferior retina than in superior retina. As color vision is generally conceived as depending on antagonistic processing of signals from different cone types, the weak opponency of identified S/M opponent cells in inferior retina suggests that color vision may in turn be weak or nonexistent in inferior retina. Note that behavioral measurements of guinea pig color vision (Jacobs and Deegan, 1994) did not control the retinal location of the stimuli, and that the spectral sensitivities revealed by those measurements are similar to the spectral properties of S/M opponent cells in superior retina.
Functional significance of the dual gradient in cone opsin expression
We speculate that the dual gradient is an adaptation that sacrificed color vision in inferior retina, presumably in favor of improved performance for some other aspect of vision. The improvement might be in achromatic contrast detection against the blue sky, the background typically seen by inferior retina (Szél et al., 2000; Osorio and Vorobyev, 2005; Peichl, 2005; Yin et al., 2006). Such an adaptation could be quite important for a species such as guinea pig, which forages in low grass and is under strong pressure from avian predators (Cassini, 1991; Cassini and Galante, 1992; Asher et al., 2004). Previously we showed that the dual gradient in cone opsin expression is preserved in the physiology of the brisk-transient ganglion cells, which are thought to mediate achromatic contrast detection. That result plus our current data strengthens the speculation that the dual gradient is an adaptation to maximize contrast sensitivity in inferior retina.
Appendix A: measuring response amplitudes with template fitting
This appendix describes the template fitting method used to characterize the response amplitudes from the PSTHs of a ganglion cell.
Rationale for the template fitting method
We assumed that the individual PSTHs of a ganglion cell could be divided into one of two groups according to their ON or OFF response polarities, and that all the PSTHs within a group shared a common temporal response waveform that differed only in overall response amplitude and baseline rate. Under these assumptions, we could form a temporal response waveform template for each group by combining data from the PSTHs within that group. We fit each individual PSTH with a scaled and shifted version of the response template. This allowed us to obtain a more reliable estimate of response amplitude to each stimulus than if we had analyzed the PSTH for each stimulus independently.
Forming the response template
Each ganglion cell had its own templates. We first determined whether the PSTHs obtained from of a ganglion cell (all blocks of trials) showed both ON and OFF responses, or whether there was only a single group of responses. This was done by eye. If we could identify one or more PSTHs of clear ON polarity and one or more PSTHs of clear OFF polarity, we assumed there were two groups and formed two response templates for the cell, one ON and one OFF. Otherwise, we formed only one response template. For the S/M opponent cells we studied, we found that the OFF template was not a sign-inverted version of the ON template.
To form the ON or OFF response template for a ganglion cell, we started by manually selecting PSTHs that (1) had clearly ON or OFF response polarities and (2) rectified mildly or not at all. That is, we tried to avoid using PSTHs where the spike rate dropped near zero for any substantial period. We then formed the response template as the weighted sum of the selected PSTHs. Each PSTH was weighted by the corresponding number of stimulus cycles that it contained. For convenience in later calculations, we normalized the template response by its maximum value. For many ganglion cells, it was not possible to avoid including some PSTHs that showed rectification.
Fitting a response template and estimating response
For opponent cells, we fit each individual PSTH with the ON template, the OFF template, and with a baseline template. For nonopponent cells, only a single ON or OFF template plus the baseline template was used. Template fitting for ON and OFF templates was accomplished by scaling and shifting, followed by clipping of any negative predicted rates to zero. The scale factor was constrained so that the peak rate was at least 1 spike/s. Baseline template fitting was accomplished by choosing a single positive baseline rate. To obtain the best fit, we minimized root mean squared error between predicted and measured PSTH time bins, excluding times when the template shape was not reliably determined because of rectification in the PSTHs used to form it.
After fitting the ON, OFF, and baseline templates we then choose the best fitting template based on the overall fit to the PSTH. For this comparison, no time bins were excluded in computing the error. A PSTH that was best fit with the ON template was considered an ON response, and its magnitude was taken as the positive peak of the fitted ON template's onset response (Figs. 2C, 4A,D, 5A). A PSTH that was best fit with the OFF template was considered an OFF response, and its response magnitude was taken as the positive peak of the fitted OFF template's offset response (Fig. 2C, 4A,D, 5A). A PSTH that was best fit with the baseline rate template was considered to have zero response. We denote the response estimated through this process at each wavelength as RPSTH(λ). Response amplitudes were either positive or zero, and classified as ON or OFF.
Iterating the fitting procedure
The procedure described above could be sensitive to the PSTHs initially selected to form the response templates. To minimize this sensitivity, we iterated the procedure as follows. For each template (ON or OFF), we fit all the cell's PSTHs as described above and classified each PSTH as ON or OFF based on the fit. We then used the fits to the PSTHs to reselect those used to form each template, based on examination of the quality of the fit of each template to each PSTH and the polarity of the best fit template. In addition, we excluded PSTHs that were rectified in time bins inconsistent with the prediction from the best fitting template. This selection-template formation-check loop was iterated until we were satisfied that the resulting templates had stabilized.
Appendix B: estimating action spectra from flash data for ganglion cells
This appendix describes how we estimate action spectra and in turn the photopigment weights of a ganglion cell. The procedures are similar to those we used in Yin et al. (2006) to analyze spectral data obtained from horizontal cells. The modifications allow us to handle responses of opponent cells.
Linear–nonlinear model
As ganglion cells had response nonlinearities, we could not convert responses to action spectra directly by dividing the measured response amplitude at each wavelength λ by the flash intensity I(λ) at that wavelength. We used an LNL model to correct for response nonlinearities (Chichilnisky, 2001; Chatterjee and Callaway, 2002; Chichilnisky and Kalmar, 2002; Zaghloul et al., 2003).
The LNL model describes the responses of a ganglion cell as the result of passing an input contrast Cinput(λ) through a nonlinear contrast-response function. The input contrast was taken as the weighted sum of the contrasts CM(λ), CS(λ), and CRod(λ) seen by M cones, S cones, and rods in response to the flash at wavelength λ:
Here, wMc, wSc, and wRodc are the M, S, and rod weights and the sum of their absolute values is constrained to unity, |wMc| + |wSc| + |wRodc| = 1. The values of CM(λ), CS(λ), and CRod(λ) are calculated from the spectral sensitivities of the photoreceptors SMi(λ), SSi(λ), and SRodi(λ), the flash intensity I(λ), and the isomerization rates seen by each receptor class to the background RhM, bg*, RhS, bg*, and RhRod, bg* (Yin et al., 2006). Spectral sensitivities were those derived in the study by Yin et al. (2006) (see also Jacobs and Deegan, 1994; Govardovskii et al., 2000) and allowed computation of isomerizations/s per receptor (Rh*) from stimulus irradiance. Contrasts were computed as CM = I(λ)SMi(λ) / RhM, bg*, CS = I(λ)SSi(λ) / RhS, bg*, and CRod = I(λ)SRodi(λ) / RhRod, bg*.
When Cinput(λ) ≥ 0, the LNL model predicted an ON or zero response. When Cinput(λ) < 0, the LNL model predicted an OFF response. We used a two-branched nonlinearity that allowed the amplitude of the response to vary with contrast differently for ON and OFF responses:
The response predicted from Equation B2 was always positive or zero, and had polarity ON or OFF according to the sign of Cinput(λ). When computing the error between predicted and measured response, we inverted the sign of the predicted response for cases where the predicted and measured polarities differed.
To fit the LNL model to a cell's responses, we used numerical search to find the parameters wMc, wSc, and wRodc in Equation B1, and the parameters Rmax_ON, Rmax_OFF, n, and C0 in Equation B2 to minimize the root mean squared error between the predictions and measurements.
Computing action spectra from linear responses
The cell's action spectrum was taken as the weighted sum of the normalized sensitivities of the M cones, S cones, and rods: SMn(λ) = SMi(λ) / RhM, bg*, SSn(λ) = SSi(λ) / RhS, bg*, and SRodn(λ) = SRodi(λ) / RhRod, bg*. The normalization accounts for the fact that we take the cell's response to be a weighted sum of receptor contrasts rather than a weighted sum of receptor isomerization rates. To plot the measured responses against the action spectrum, we corrected each datum for the cell's nonlinearity:
In cases where the predicted polarity (based on the sign of Cinput(λ)) differed from the measured polarity, we inverted the sign of Rlinear(λ).
Appendix C: DOG receptive field model
We used a DOG receptive field model in two dimensions (Enroth-Cugell and Robson, 1966; Derrington and Lennie, 1984; Solomon et al., 2005) to describe the spatial structure of the receptive fields of S/M opponent cells. For each cell, we assumed the receptive field (RF) was the weighted difference between a concentric Gaussian center and surround. The normalized spatial profiles of the center and surround were as follows:
Here, σCenter and σSurround determine the spatial extent of the center and surround and r is distance from the center of the receptive field. Because we did not observe any shift in response phase with stimulus size, we did not include a parameter to describe a relative delay in the response of the surround with respect to the center (Enroth-Cugell et al., 1983).
We further assumed that both the receptive field center and surround received a weighted combination of M and S inputs, the relative strengths of which were described by four non-negative weights: wMCenter, wSCenter, wMSurround, and wSSurround. The constraint of non-negativity on wMCenter and wSCenter means that the M and S inputs contributed to center with the same sign as each other; the constraint of non-negativity on wMSurround and wSSurround means that the two inputs also had the same sign as each other for the surround. Although we could have used a more general form for the spatio-chromatic properties of the cell, the form we used captures the notion of center-surround antagonism and provided a good fit to the area summation data. The overall opponency of the cell was incorporated in the model by taking the difference between center and surround.
To account for the rod input to the cell, we assumed that rod signal flowed to cones through rod-cone gap junctions, with a coupling strength of KRod-Cone. Therefore the relative strengths of rod inputs to the center and surround were described as follows:
In fitting the data, we fixed the value of KRod-Cone to the value given in Table 1, derived using the functional model described in Appendix D. For convenience we added the constraints that wMCenter + wSCenter + wRodCenter = 1 and wMSurround + wSSurround + wRodSurround = 1, and then used a parameter KSurround/Center to describe the ratio of center-surround strength.
As with our basic analysis of cell responses, we included a static nonlinearity (see Appendix B). Because the area summation data by itself did not contain enough contrast variation to constrain the shape of this nonlinearity, we used additional data obtained with a fixed spot size for this purpose. For cells where we had data from monochromatic flashes (such as those shown in Fig. 4), we used these. For cells where we did not have monochromatic flash data, we instead used measurements to flickering spots with varying combinations of CM and CS, which we had obtained for these cells.
We used numerical search to find the estimates of the receptive field parameters that best predicted the cell's responses. The error measure minimized was the root mean squared error between the predicted and the measured response amplitudes.
Appendix D: functional model for photopigment weights
The model is characterized by 4 parameters. These are (1) fS, the proportion of genuine S cones (Haverkamp et al., 2005), taken to be constant across the retina; (2) gS, the gain of signals from the genuine S cones through the S cone bipolar to S/M opponent ganglion cells, expressed relative to the gain for these same signals as they reach nonopponent cells; (3) cR-C, the strength of rod coupling to cones (Lee et al., 2003; Hornstein et al., 2005; Verweij et al., 2008), taken to be independent of cone type (Ahnelt et al., 1990) and retinal location; and (4) β, the proportion of M-pigment coexpressed in the cones that are nongenuine S cones in inferior retina. (This proportion is fixed at β = 1 in superior retina.)
Model description for nonopponent cells
We start with the formula for the S weight for nonopponent cells (horizontal and brisk-transient ganglion cells) in superior retina. Here the strength of S input relative to M input (excluding rods) is given directly by fS/(1 − fS). This is because we assume that the gain from cones to nonopponent cells is the same for all cone types and because the S and M pigments are assumed to be completely segregated within separate cones. There is also a rod weight, which we take directly as the parameter cR-C. This convention defines the units of cR-C in terms of the strength of rod input per cone. This leads to the expression for the relative S weight as follows:
for nonopponent cells in superior retina.
We can generalize Equation D1 for the case where S and M pigment are both expressed in the nongenuine S cones, as is the case for inferior retina.
Here the quantity (1 − β) is the strength of S opsin input from nongenuine S cones. Note that Equation D2 reduces to Equation D1 when β = 1.
In analogy to Equation D2, we also have expressions for the relative M and rod weights to nonopponent cells in both superior and inferior retina.
and
Note that the weights obtained through Equations D2, D3, and D4 were normalized so that their sum was one.
Model description for opponent cells
The formulae for photopigment weights for opponent cells are based on the same ideas, but take into account the amplification of signals from genuine S cones. Note that this amplification is applied both to the S opsin signal for these cones and also to the rod signal that passes through these cones. We obtain for the magnitude of the cone weights the following:
and for the rod weight
The sign of wSopponent was positive for S+/M− cells and negative for M+/S− cells, and vice-versa for the sign of wMopponent. The sign of the rod weight given by Equation D7 is relative to the sign of wSopponent: the signed value computed through Equation D7 was multiplied by the sign of wSopponent. The sum of the magnitudes of the weights obtained through Equations D5–D7 is one. In fitting the model, we averaged the weights of S+/M− and M+/S− by inverting the sign of the weights measured for M+/S− cells.
Fitting the model parameters
Equations D2–D7 predict the photopigment weights of opponent and nonopponent cells in both superior and inferior retina as a function of the four model parameters. We fixed fS at 5.5% based on our immunostaining data from superior retina [this number is given as rounded to 6% in the work of Yin et al. (2006) and elsewhere in this study]. We used numerical search to find the values of the remaining three parameters that provided the best fit (in a least squares sense) to the mean values for each cell type, as measured here and in our previous study (Yin et al., 2006). In determining the mean values we combined S+/M− and M+/S− opponent cells, inverting the sign of weights for the M+/S− cells.
Footnotes
-
This work was supported by National Eye Institute Grants EY10016, EY08124, EY016607, and P30 EY001583. The data are also presented as part of the first author's Ph.D. thesis (L. Yin, 2008, unpublished Ph.D. thesis, Neuroscience Graduate Group, University of Pennsylvania). We thank Noga Vardi, Ying Xu, and Marie Fina for providing advice on immunostaining, confocal imaging, and quantification of dendritic morphology; Michael Freed for sharing data on the classification of guinea pig ganglion cells; Edward Pugh and Peter H. Li for useful discussions on a proposed quantitative model; John Andrews-Labenski, Michael Suplick, Fred Letterio, and Jian Li for technical assistance; and Sharron Fina for manuscript preparation.
- Correspondence should be addressed to Dr. David H. Brainard, Department of Psychology, University of Pennsylvania, Suite 302C, 3401 Walnut Street, Philadelphia, PA 19104. brainard{at}psych.upenn.edu