Abstract
Many of the parvocellular pathway (PC) cells in primates show red–green spectral selectivity (cone opponency), but PC ganglion cells in the retina show no anatomical signs of cone selectivity. Here we asked whether responses of PC cells are compatible with “random wiring” of cone inputs. We measured long-wavelength-sensitive (L) and medium-wavelength-sensitive (M) cone inputs to PC receptive fields in the dorsal lateral geniculate of marmosets, using discrete stimuli (apertures and annuli) to achieve functional segregation of center and surround. Receptive fields between the fovea and 30° eccentricity were measured.
We show that, in opponent PC cells, the center is dominated by one (L or M) cone type, with normally <20% contribution from the other cone type (high “cone purity”), whereas non-opponent cells have mixed L and M cone inputs to the receptive field center. Furthermore, opponent response strength depends on the overall segregation of L and M cone inputs to center and surround rather than exclusive input from one cone type to either region. These data are consistent with random wiring. The majority of PC cells in both foveal (<8°) and peripheral retina nevertheless show opponent responses. This arises because cone purity in the receptive field surround is at least as high as in the center, and the surround in nearly all opponent PC cells is dominated by the opposite cone type to that which dominates the center. These functional biases increase the proportion of opponent PC cells, but their anatomical basis is unclear.
Introduction
In trichromatic (“normal color vision”) primates, most parvocellular pathway (PC) cells show functional antagonism of inputs from long-wavelength-sensitive (L) and medium-wavelength-sensitive (M) cones to the receptive field (DeMonasterio et al., 1975; Derrington et al., 1984; Lee et al., 1987, 1990; Yeh et al., 1995a; Benardete and Kaplan, 1999). Most foveal PC cells derive input from a single L or M cone via a single midget bipolar cell (Kolb and Dekorver, 1991; Calkins et al., 1994; Jusuf et al., 2006a) and thus should get a dominant functional input to the receptive field center from a single cone type. Under the “random wiring” hypothesis (Mollon et al., 1984; Shapley and Perry, 1986; Lennie et al., 1991), this one-to-one connectivity is the dominant source of spectral bias in PC fields and thus produces the afferent signals serving red–green color vision. Consistently, anatomical analyses of the foveal retina have so far revealed no other sign of selective wiring for L and M cones (Wässle et al., 1989; Calkins and Sterling, 1996, 1999; Jusuf et al., 2006a).
The random wiring hypothesis predicts that, in peripheral retina, there will be mixing of input from L and M cones to the receptive field center (Lennie et al., 1991; Mullen and Kingdom, 1996; Calkins and Sterling, 1999; Diller et al., 2004), and this is supported by anatomical and physiological observations (Diller et al., 2004; Jusuf et al., 2006b) (but see Solomon et al., 2005). It is therefore unclear why many PC cells in peripheral retina show cone-opponent signals (DeMonasterio and Gouras, 1975; Solomon et al., 2005). Furthermore, the question whether the surround of PC cells shows mixed input from L and M cones has not been resolved conclusively (Smith et al., 1992; Lankheet et al., 1998b; Lee et al., 1998; Reid and Shapley, 2002). Finally, the fundamental question of how chromatic sensitivity depends on functional segregation of L and M cone inputs has not been answered definitively.
Cone inputs to PC receptive fields have been measured using spatial- or temporal-frequency modulation (Kaplan et al., 1989; Smith et al., 1992; Yeh et al., 1995b; Lankheet et al., 1998b; Benardete and Kaplan, 1999) or spatially discrete stimuli such as pseudorandom checkerboards and counterphase edges (Lee et al., 1998; Reid and Shapley, 2002). In these methods, the receptive field center and surround are modulated simultaneously. This makes it difficult to measure cone mixing in the surround, because the overall response of the cell is dominated by the center mechanism. Here, we minimized this problem by using annuli to stimulate preferentially the receptive field surround and apertures to stimulate the receptive field center (Westheimer, 1967; Kilavik et al., 2003). We used cone-isolating (“silent substitution”) modulation to measure the contribution of L and M cones to center and surround. Our experimental goals were to relate the amplitude of cone-opponent responses to the functional weight of cone inputs and to compare the cone-opponent signal in foveal and peripheral retina.
Materials and Methods
Female marmosets (Callithrix jacchus) were obtained from the Australian National Health and Medical Research Council combined breeding facility. Three trichromatic adult females, carrying the genes for the 543 and 563 nm opsins [identified by PCR–restriction fragment length polymorphism (Blessing et al., 2004)] were selected for physiological recordings. Animals were anesthetized and prepared for extracellular single-unit recordings using experimental protocols described in detail previously (Blessing et al., 2004). Procedures were approved by institutional Animal Experimentation and Ethics Committee and conform to the Society for Neuroscience policy on the use of animals in neuroscience research. Action potential time series from individual PC cells in the dorsal lateral geniculate nucleus were digitized at 10 kHz. Responses were subjected to Fourier analysis. The first harmonic amplitude and phase were used as response measures.
Stimuli were presented on a computer monitor at 80 Hz frame refresh rate and centered on the receptive field via a gimbaled mirror. Stimulus–response curves for contrast, aperture, and annulus tuning were recorded using modulation along four color axes: L and M cone isolating, isoluminant chromatic (L − M), and luminance (L + M). Temporal frequency was 4 Hz. This value was chosen because it yielded robust responses in nearly all PC cells recorded but is low enough to avoid the masking of opponent responses that occurs at higher temporal frequencies (cf. Diller et al., 2004; Solomon et al., 2005). Mean luminance was close to 25 cd · m−2. Cone-isolating stimuli (calculated by convolution of the monitor phosphor spectra with photoreceptor spectral templates) were confirmed as accurate in separate experiments (Blessing et al., 2004; Forte et al., 2006). Contrast for L + M tuning functions was set to evoke close to 25 action potentials (impulses) per second on a non saturating region of the contrast–response curve. Contrast gain (impulses · s−1 · %−1 Michelson contrast) was estimated as the initial slope of the best-fit Naka–Rushton function (Croner and Kaplan, 1995; Yeh et al., 1995a). If the response at maximum stimulus contrast was <10 impulses/s, we estimated response gain as response amplitude at maximum stimulus contrast/100. The stimuli were generated using only the red and green monitor phosphors. Cells were also tested for signs of functional input from short-wavelength-sensitive (S) cones as described previously (Forte et al., 2006). None of the PC cells showed significant input from S cones.
The maximum Michelson cone contrast was close to 36% for L and M cone-isolating stimuli; the cone contrast for the chromatic (L − M) stimulus was adjusted to be half this value for each cone class. Accordingly, we categorized cells as non-opponent if the gain ratio (L − M)/(L + M) for large-field stimuli was ≤0.36 (that is, the sum of L and M cone contrasts) or if the maximum response to the L − M stimulus was below 10 impulses/s. Thus, an opponent cell response to L − M modulation is greater than the response to L + M modulation at the same cone contrast. However, the reader should note that these criteria do not imply the existence of discrete populations of PC cells, and the results show a continuum of response properties on all stimulus dimensions that we measured. Our purpose here was to relate our study to categorizations made by others and to the qualitative pattern of responses elicited by hand-held stimuli.
Results
Predictions of the random wiring hypothesis
The random wiring hypothesis is illustrated in Figure 1, which shows how a PC receptive field could draw functional inputs from a trichromatic cone mosaic in the perifoveal retina. Lennie et al. (1991) showed that red–green opponent responses of foveal PC cells could be explained by functional input from a single L or M receptor to the receptive field excitatory mechanism, opposed by a Gaussian-weighted mixture of L and M cones in the inhibitory mechanism. A simple elaboration of this model (Mullen and Kingdom, 1996) can be used to illustrate the effect of convergent input from L and M cones to both center and surround. Here (Fig. 1), the cone inputs are expressed as a sum of four Gaussian-weighted components (Rodieck and Stone, 1965; Derrington and Lennie, 1984). Each “pair” of Gaussians, (L and M excitation) and (L and M inhibition), is set to have the same width, but the L and M components of each pair can have different amplitudes. The center draws from a smaller number of cones than the surround, and thus the relative weight of L and M cone inputs [“cone purity” (Mullen and Kingdom, 1996)] shows greater imbalance in the center than in the surround. In the example shown, the excitatory (center) mechanism is dominated by M cones (Fig. 1C), whereas the inhibitory (surround) mechanism gets mixed L/M cone input close to the average L/M cone ratio (Fig. 1B). With increasing size of the receptive field, the degree of asymmetry in center and surround weights [“opponent purity” (Mullen and Kingdom, 1996)] is predicted to fall because the center draws from larger numbers of cones. These two predictions of the random wiring hypothesis (low cone purity in the surround at all eccentricities, and decreasing center cone purity with increasing eccentricity) can be tested using stimuli that differentially activate the center and surround but are insensitive to the exact position of each M or L cone in the mosaic. We show below that these requirements are met by temporal modulation of apertures and annuli and that both the center and surround of the majority of PC cells shows high cone purity.
Comparison of responses in foveal and peripheral visual field
Complete datasets were obtained for 32 PC cells. A complete dataset comprised triplicate measurement of responses to 12 aperture (“spot”) and 13 annulus (“donut”) sizes on L − M, L + M, L-isolating, and M-isolating cone axes, together with spatial frequency and contrast tuning curves, yielding a total of ∼360 measurement points for each cell. Data from five additional cells with incomplete aperture and annulus tuning data and two cells with incomplete spatial frequency tuning data were consistent with the complete datasets and are included when possible in the results reported below.
Figure 2 shows responses to chromatic (L − M) and luminance (L + M) modulation for opponent PC cells recorded with receptive fields near the fovea (Fig. 2A) and near 15° eccentricity (Fig. 2B) and one non-opponent PC cell (Fig. 2C) recorded near 15° visual field eccentricity. The pattern of responses is consistent with that reported previously for foveal and peripheral PC ganglion cells in macaque (Solomon et al., 2005) and foveal PC cells in marmosets (Yeh et al., 1995a; Blessing et al., 2004). All cells show mildly saturating contrast–response relationship for L + M modulation, and the opponent cells additionally display vigorous responses to L − M modulation. Figure 2D shows L − M and L + M gain for all PC cells recorded, together with the criterion (L − M)/(L + M) gain ratio (Fig. 2D, dashed line).
Figure 3 shows that the increase in receptive field center radius with eccentricity (Fig. 3A) is not associated with a large decrease in the proportion of opponent PC cells (Fig. 3B). Mean center radius for PC cells recorded above 8° eccentricity (0.127 ± 0.042°, mean ± SD; n = 21) was approximately three times greater than for PC cells recorded below 8° (0.039 ± 0.015°, mean ± SD; n = 13; p < 0.01, unpaired Wilcoxon's test), but the proportion of opponent cells recorded above 8° (13 of 23, 56.5%) is only slightly lower than the proportion recorded below 8° (10 of 14, 71%). Because we only classified cells as opponent if the L − M gain was >0.1 impulses · s−1 · %−1, there is some overlap of opponent and non-opponent cells near the equal cone contrast ratio of 0.36 (Fig. 3B). However, even under these conservative criteria, it is clear that chromatic opponent PC cells are not restricted to the central visual field. The log gain ratio for cells recorded below 8° eccentricity (0.005 ± 0.432, mean ± SD; n = 21) is slightly greater than for PC cells recorded above 8° (−0.220 ± 0.480, mean ± SD; n = 23; p = 0.11, unpaired Wilcoxon's test). On the basis of these variance measures, a sample of >40 cells in each group would be required to detect a decline of 50% in mean gain ratio with a power of 0.8 (Lenth, 2006). In summary, there is little evidence for a systematic decline in opponency with eccentricity in the sample of PC cells we measured.
The known anatomical properties of the marmoset retina allow the expected convergence of cone photoreceptors to the receptive field center of extrafoveal PC cells to be estimated. Figure 4A shows a drawing of a midget ganglion cell at 13° eccentricity (Jusuf et al., 2006b). The cell drawing is superimposed on a drawing of the cone pedicle mosaic close to this eccentricity (Fig. 4A, gray patches). The average excitatory center radius of PC cells (1 SD of the Gaussian weighting function) at this position is also indicated (Fig. 4A, red circle). This shows first that the physiological estimate of center radius of PC cells is consistent with the dendritic field radius of midget cells in peripheral marmoset retina, and second that the peripheral receptive fields would be expected to receive convergent input from >15 cones.
Goodchild et al. (1996) showed that, in extrafoveal primate retina, the dendritic field area of PC cells increases in an approximately inverse relationship to cone photoreceptor density. This allows a simple (albeit indirect) estimate of convergence to each PC cell in our dataset. Figure 4B shows data points replotted from Goodchild et al. (1996), their Figure 13. The average midget ganglion cell dendritic field area (square millimeters) divided by the average cone photoreceptor density (cones per square millimeters) is shown at a set of eccentricities above 5°. These data are well fit by linear regression Y = −22.1281 + 3.3488 × X, where Y is number of cones per cell, and X is eccentricity in degrees (r2 = 0.99) (Fig. 4B, solid line). The convergence measured for the cell illustrated in A (Fig. 4B, blue asterisk) is consistent with this prediction. The linear regression was applied to the eccentricity of extrafoveal (>5°) PC cells in our dataset, as shown in Figure 4C. For simplicity, we did not incorporate our measured variation in receptive field center size in this estimation. Convergence for cells below 5° eccentricity was set to 1, corresponding to the one-to-one connections of foveal midget ganglion cells (Kolb and Dekorver, 1991; Calkins et al., 1994; Jusuf et al., 2006a). It can be seen that the majority of PC cells is expected to receive convergent input from cone photoreceptors to the receptive field center. We conclude that the presence of chromatic opponent responses in marmoset PC cells is not dependent on one-to-one connectivity of the foveal midget pathway. We next asked how the cone inputs to PC pathway cells are organized to produce opponent responses, by a spatiochromatic analysis of the receptive field.
Spatiochromatic analysis of the receptive field
In the following, the reader should keep in mind a subtle distinction between the receptive field regions of a PC cell and the physiological mechanisms that give rise to these regions. The receptive field surround is a more-or-less annular-shaped inhibitory region, which is concentric with an excitatory center region (Kuffler, 1953; Wiesel and Hubel, 1966). However, as illustrated in Figures 1 and 5A, the process that gives rise to these regions is well described as a combination of a (weak but spatially broad) zone of inhibition, which overlaps with a (powerful but spatially restricted) zone of excitation (Rodieck and Stone, 1965; Dreher et al., 1976; Derrington and Lennie, 1984). Here we use the terms center and surround to refer to discrete regions of the receptive field and excitation and inhibition to refer to the corresponding mechanisms.
It is well known that the difference-of-Gaussians (DOG) receptive field model (Rodieck and Stone, 1965; Enroth-Cugell and Robson, 1966; Frishman et al., 1987) can describe the responses of PC cells to drifting sinusoid luminance gratings (Derrington and Lennie, 1984; Croner and Kaplan, 1995; Kremers and Weiss, 1997; White et al., 2001). Kilavik et al. (2003) and Kremers et al. (2004) showed that this model also can account for PC cell responses to counterphase luminance modulation, in apertures or annuli that are centered on the receptive field. We therefore explored the potential of such spatially discrete stimuli to stimulate selectively the center and surround and thus enable measurement of L and M cone contributions to each.
A prerequisite for such measurements is that the stimulus remains concentric with the receptive field, so it is important to ensure that drifts in eye position do not affect the recordings (Forte et al., 2002). Accurate centering of the stimulus on the receptive field was confirmed before and (in nearly all cases) after data collection, by monitoring the position of maximal response to a small (<0.05° radius) flashed spot. On rare occasions, drifts in eye position became evident during or after data collection, because such drifts caused high response variability for small stimuli and characteristic changes in the shape of annulus tuning curves. These data were discarded.
For a subset of cells (n = 7), we also quantified eye stability by mapping the spatial response profile for temporal square-wave flickered (4 Hz) horizontal and vertical bar stimuli (length, 8°; width, 0.05–0.2°). Figure 5B shows excitatory response profiles of four sequentially measured PC cells. Two sets of isoresponse contours are shown for each cell; these measurements were taken before and after data collection. The contours show a high degree of overlap, indicating that eye position has remained stable. The mean receptive field drift over this period (between 0.5 and 2 h) was 0.078 ± 0.052° (mean ± SD; n = 7), which is close to the diameter of foveal receptive field centers (compare Figs. 3A, 5B). In summary, we found that, with careful monitoring, the position of small stimuli can be reliably centered on the receptive field of PC cells for periods of at least 2 h. The filled circles in Figure 5B show the diameter of the aperture stimulus used to stimulate preferentially the receptive excitatory mechanism. This stimulus was normally set to be slightly larger than the excitatory field. The following analysis shows that this configuration gives a high degree of functional isolation yet minimizes the effects of eye drift, especially with foveal receptive fields (Forte et al., 2002).
Figure 5, C and D, shows “side-view” sketches of aperture and annulus stimuli in relation to the difference-of-Gaussians receptive field model. The responses of the model to these stimuli were calculated as follows: where R is the response modulation amplitude of the cell, and Rexc and Rinh are the response amplitudes of the excitatory (center) and inhibitory (surround) subunits, respectively. The response of each subunit is calculated as the integral of a two-dimensional Gaussian under the stimulus shape (aperture or annulus) following Kilavik et al. (2003). For an aperture stimulus of radius rAP, where ai is the subunit sensitivity (in impulses · s−1 · degrees−2), and σi is the standard deviation of the Gaussian (used as a measure of subunit radius) for subunits i ϵ (exc; inh).
For annular stimuli, we assumed that the outer diameter was sufficient to cover the entire receptive field. Using this simplification, the response to an annulus of inner radius rAN is modeled (using the notation of Eq. 2) as follows: The curves in Figure 5, E and F, shows the predicted responses of each subunit and the cell for aperture (Fig. 5E) and annulus (Fig. 5F) stimuli. Responses are normalized to the maximum cell response. These predictions are based on the average receptive field parameters we obtained (Fig. 3A) for our sample of receptive fields above 8° eccentricity (n = 17): center radius (σexc) = 0.12°; surround radius (σinh) = 1.00°; and sensitivity ratio (aexc/ainh) = 16.7. With increasing aperture diameter, the cell response (exc − inh) rises sharply as the center is recruited and then falls gradually as the surround contribution increases more rapidly than the center contribution. Consistently, as annulus inner diameter increases (Fig. 5F), the dominant drive from the receptive field center falls rapidly until it is balanced by the inhibitory subunit to yield a response minimum. The inhibitory mechanism dominates the cell response for annuli with larger inner diameters. Thus, although both the excitatory and inhibitory mechanisms contribute to the cell response for all stimuli, the relative contribution of these mechanisms varies dramatically with stimulus spatial configuration.
Figure 5G shows the proportional contribution of the center mechanism to the cell response as a function of aperture diameter. Over the range of small apertures in which cell response amplitude is rising, the surround contribution is negligible, because the surround sensitivity is much lower than that of the center. An aperture of close to 0.5° diameter thus produces highly selective activation of the excitatory mechanism. In the same way, annuli with inner diameter close to 0.75° (Fig. 5H) can be used to activate selectively the inhibitory mechanism (although at a lower absolute cell response amplitude; Eq. 3). We refer to these stimuli as “center-optimized aperture” (Fig. 5G, filled star) and “surround-optimized annulus” (Fig. 5H, open star) to indicate this functional segregation.
Figure 6 shows responses of a PC cell recorded at 17.7° eccentricity to illustrate how center-optimized and surround-optimized stimuli were determined. This cell, which has a relatively strong inhibitory surround (Table 1), illustrates the key features of our measurements. The spatial configuration of the stimulus is sketched in Figure 6A. As expected (Kilavik et al., 2003), response amplitude for luminance (L + M) modulation is captured well by the DOG model (Fig. 6B, solid lines). The response phase for aperture tuning is dominated by the center mechanism (Fig. 6C, left), whereas there is an abrupt phase shift of ∼180° at the amplitude null point (Fig. 6C, right) on the annulus tuning curve. This result conforms nicely with the predictions shown in Figure 5F: the dominant drive to the cell response shifts from excitatory (at small inner diameters) to inhibitory (at large inner diameters).
The tuning curves for chromatic modulation (Fig. 6D) support previous evidence that chromatic opponent responses in parvocellular cells are enhanced by synergy between excitatory and inhibitory mechanisms (Wiesel and Hubel, 1966; DeMonasterio et al., 1975; Dreher et al., 1976; Smith et al., 1992; Lankheet et al., 1998b; Lee et al., 2000; Reid and Shapley, 2002; Saito et al., 2004; Solomon et al., 2005). The peak response amplitude for L − M modulation is approximately half the amplitude of response to L + M modulation, but the cone contrast (18%) is much lower for the L − M stimulus than for the L + M stimulus (100%). Thus, L − M response gain is ∼2.8 times greater than L + M gain. Responses to large L − M apertures show little sign of surround attenuation, and response phase is primarily independent of stimulus size (Fig. 6D, filled gray histograms).
The aperture tuning curves obtained using L and M cone-isolating stimuli (Fig. 6E,F, left) show bandpass characteristic, suggesting that both cone types contribute to both (excitatory and inhibitory) components of the receptive field. There is, however, a functional bias: excitation is dominated by the M cone, whereas inhibition is dominated by the L cone. Accordingly, at all aperture diameters, the response phase for M cone modulation is dominated by the excitatory mechanism (Fig. 6E, left, histograms), and the annulus tuning curve for M cone modulation (Fig. 6E, right) shows a poorly defined null, with feeble responses to larger annuli. Conversely, the response phase for L cone modulation shows an ∼180° shift between small and large apertures, and for all annuli the effect of L cone modulation is inhibitory (Fig. 6F, histograms). Thus, there is substantial, but incomplete, segregation of M and L cone inputs to excitatory and inhibitory inputs to this PC cell. Responses to all stimulus configurations were highly reliable across trials, as indicated by the small error bars in Figure 6. This is additional evidence that eye position has remained stable.
In principle, the Gaussian integral analysis described above (Eqs. 2, 3) (Fig. 3) can be used to specify the aperture and annulus diameters that segregate optimally the excitatory and inhibitory inputs. Under the Gaussian model, this segregation can never be complete, because these two mechanisms are overlapping and have infinite spatial extent. Furthermore, PC cell receptive fields show deviation from circular symmetry (Smith et al., 1990; White et al., 2001; Passaglia et al., 2002; Forte et al., 2006), which together with unavoidable small drifts in receptive field position (Fig. 3) and the poorly predictable effects of chromatic aberrations at the stimulus edges (Forte et al., 2006), make the theoretical goal of ideal segregation impossible to achieve in practice. For our analysis, we therefore combined the predictions from the difference-of-Gaussians model with empirical analysis of response phase, to choose a center-optimized aperture and a surround-optimized annulus for each cell, as follows.
First, the excitatory response phase (“On” or “Off”) was measured at the peak of the L + M aperture tuning curve (On-center cells, 308.4 ± 7.4°, mean ± SD, n = 17; Off-center cells, 136.0 ± 17.3°, mean ± SD, n = 17). Second, the center-optimized aperture was taken as the aperture closest to the peak of the L + M tuning curve, with the additional constraint that the response phase for the cone-isolating (M or L) conditions was within 90° of the excitatory response phase (Fig. 6A, C–F, filled stars). Finally, the surround-optimized annulus was taken as the annulus closest to the peak of the L + M annulus tuning curve, with the constraint that the response phase for the cone-isolating conditions was greater than ±90° of the excitatory response phase. The upshot of this process for nearly all PC cells was to leave a “gap” between the outer border of the center-optimized aperture and the inner border of the surround-optimized annulus. Previous analyses (Derrington et al., 1984; Smith et al., 1992; Lankheet et al., 1998b) show that, in macaques, PC cells combine L and M cone inputs linearly at low cone contrasts, so observations made in noncontiguous regions of the receptive field can be reasonably extrapolated to explain cell behavior. Analysis of our data (described below) supports this assumption.
The cone inputs to center and surround were estimated from response amplitude to L and M cone-isolating conditions, for these optimized stimuli. In some cases (4 of 34 aperture and 11 of 34 annulus tuning curves), the phase constraints could not be met for both of the cone-isolating conditions (for example, if a cell receives pure M cone excitation, the phase of L cone response will never be close to the excitatory response phase). In these cases, the response amplitude for that cone was set to 0 for that condition.
As noted above, under the Gaussian receptive field model, the excitatory and inhibitory mechanisms are spatially overlapping and have infinite extent, so no spatial configuration can segregate these mechanisms completely. We therefore estimated the functional segregation achieved by our optimized stimuli by integrating the two Gaussian components using the receptive field parameters σexc, σinh, aexc, and ainh derived from the (M + L) spatial frequency tuning function for each cell (Eqs. 2, 3). These calculations showed maximum 4% contribution of the inhibitory Gaussian to center-optimized stimuli (0.024 ± 0.045, mean ± SD; n = 34). The excitatory Gaussian contribution to surround-optimized stimuli was also low (mean of 0.034) but more variable (SD of 0.076; range of 0–0.33). These values (together with the response phase criteria described above) give us confidence that our method gives a greater degree of segregation of excitatory and inhibitory input to PC cells than has been achieved previously. We therefore used this method to estimate the L and M cone contributions to PC cell center and surround.
Phase and amplitude of cone inputs
Cone weights and response phase for the 34 PC cells studied are shown as polar plots in Figure 7. Response phase for each cell is referenced to the response phase for center-optimized L + M modulation (vertical axis). Amplitudes are shown relative to the dominant cone for each (excitatory or inhibitory) mechanism. It can be seen that the phases of cone inputs to each (excitatory center and inhibitory surround) region are closely aligned and that center and surround inputs are approximately in opposite temporal phase at the contrast modulation frequency we used (4 Hz).
The radius of the gray circle at the center of each plot represents 20% of the dominant cone weight. Thus, it becomes apparent that two-thirds of the opponent PC cells (16 of 21) show <20% cone mixing in the receptive field center. Furthermore, the surround in nearly all (20 of 21) opponent cells is dominated by the opposite cone type to that which dominates the center. These results show that strong bias in cone weights is associated with opponent responses, but that complete segregation of cone types is not essential. In contrast to opponent cells, nearly all (12 of 13) of the non-opponent cells show >20% cone mixing to the center mechanism, and the surround of many non-opponent cells (8 of 13) is dominated by the same cone type that dominates the center. When considered as an entire population, PC cells show a full range of cone weights, and there are cases in which the distinction between opponent and non-opponent categories is not clearly reflected in differences in cone weight and temporal phase. We interpret this to mean that the high proportion of opponent PC cells is not the result of all-or-none selective wiring. Rather, it is the result of a more subtle process that biases cone weights in favor of opponent inputs to PC cells.
Both center and surround of opponent PC cells shows cone bias
We quantified the degree of cone segregation by the relative weight [L/(L + M)], or by cone purity, which describes the ratio of L − M difference to the total subunit weight (Mullen and Kingdom, 1996). A cone purity value of 0 means equal weight of L and M cones, whereas 1 means functional input from either L or M cones alone. Opponent purity describes the overall segregation of L and M cones between center and surround and is calculated as the absolute difference in L cone ratio between center and surround (Mullen and Kingdom, 1996): [L/(L + M)]ce − [L/(L + M)]su, where ce is center and su is surround. An opponent purity value of 0 means identical cone ratios, whereas 1 indicates a perfectly cone-opponent receptive field (pure L center and pure M surround or vice versa). Figure 8, A and B, shows the result.
For opponent cells, the cone purity of both the center and surround is greater (median center purity, 0.65; median surround purity, 0.93; n = 21) than for non-opponent cells (median center purity, 0.29; median surround purity, 0.35; n = 13; p < 0.01, unpaired Wilcoxon's test). Furthermore, for opponent cells, the surround cone purity is at least as high as the center cone purity (n = 21; p = 0.204, unpaired Wilcoxon's test), despite the fact that the surround draws input from a larger number of cones: in our sample, the surround diameter is, on average, 6.29 times larger than the center diameter. Finally, we found significant negative correlation (r = −0.79; p < 0.001) between center and surround cone weights, so cells with L cone-dominated centers tend to have M cone-dominated surrounds and vice versa (Fig. 8A). These results all are consistent with preferential, but not exclusive, selection of opponent cone types in the receptive fields of opponent PC cells.
Figure 8B shows that opponent PC cells are characterized by overall imbalance of cone types driving the centers and surround rather than exclusive cone inputs to either mechanism. In Figure 8C, the ratio of chromatic to luminance gain is plotted against opponent purity. The strength of opponent response is in approximately exponential relationship to opponent purity (correlation between log gain ratio [L − M]/[L + M] and opponent purity, r = 0.58; p < 0.01). Our criterion value for classifying cells as opponent therefore represents an opponent purity close to 0.5. Cells above this opponent purity value show average (L + M)/(L − M) gain ratio 5.7-fold greater than cells below this value. The opponent purity of green-On and green-Off cells (0.78 ± 0.14, mean ± SD; n = 10) is slightly greater than opponent purity of red-On and red-Off cells (0.67 ± 0.26, mean ± SD; n = 11), but the difference is not significant (p = 0.47, unpaired Wilcoxon's test). In summary, a cone weight difference of ∼2:1 between center and surround is sufficient to produce robust chromatic opponent responses in PC cells.
Opponent purity is greater than predicted by random wiring
A “hit-and-miss” binomial model of L and M cone inputs to PC receptive fields (Mullen and Kingdom, 1996) can be used to predict opponent purity of a randomly wired PC cell given knowledge of the number of cones converging to the receptive field center mechanism. We therefore used our estimates of cone convergence to PC cells in our dataset (see above, Fig. 4) to compare the predictions of the random wiring hypothesis with our measured values. We simplified the comparison by using the same L/M cone ratio (1:1) and ratio of center and surround cone numbers (1:6) as shown by Mullen and Kingdom (1996), their Figure 1B. This choice is justified on three grounds. First, the average L/M cone ratio in marmoset retina is expected to be close to 1:1 (Lyon, 1972; Travis et al., 1988; Tovée, 1993; Bowmaker et al., 2003) (see Discussion). Second, even twofold changes in the L/M ratio have only small effects on average opponent purity (Mullen and Kingdom, 1996). Third, on average, the radius of surround that we measured is ∼6 times larger than center radius, consistent with a surround to center cone number ratio of at least 6:1. Values from the published figure were fit to the form where c1 = 1.0813, c2 = 0.2723, λ1 = 1.5079, λ2 = 0.0647, and n is the number of cones in the receptive field center. This enables the opponent purity for a given value of cone convergence to be estimated.
The result of this procedure is shown in Figure 9A. As expected, the maximum opponent purity (0.5) is predicted for foveal PC cells, because they show one-to-one cone convergence to the receptive field center and, on average, equal weight of L and M cones in the surround. Opponent purity is predicted to fall rapidly with increasing convergence, so that cells outside the fovea [which receive convergent input from at least eight cones to the receptive field center (Fig. 4C)] are predicted to have opponent purity below 0.2.
Measured and predicted opponent purity are compared for each cell in Figure 9B. For nearly all PC cells (31 of 34, 91%) the measured opponent purity is greater than predicted by the random wiring (hit-and-miss) model. On average, measured opponent purity (0.55 ± 0.32, mean ± SD; n = 34) is more than double the value predicted by random wiring (0.21 ± 0.22, mean ± SD; p < 0.01, paired Wilcoxon's test). This result is consistent with the results shown in Figures 7 and 8. The functional bias for either L or M cones to dominate in the surround and the negative correlation between center and surround cone weight both serve to “push” opponent purity above the average predicted by random wiring. A more extensive modeling study could incorporate stochastic variation in cone weight to predict not only the average value but also the proportion of opponent PC cells at each eccentricity. We did not attempt this exercise.
Contrast saturation could enhance PC responses to low-contrast chromatic stimuli
The approximately exponential relationship between chromatic response amplitude and opponent purity (Fig. 8C) suggests that responses to chromatic contrast are selectively enhanced in PC cell receptive fields. A possible basis for this effect is shown in Figure 10, in which we compare responses to cone-isolating and L − M modulation across the PC population and across all aperture sizes. Figure 10A compares for one green-Off cell the aperture tuning curves for L cone (red circles) and M cone (green squares) modulation. The cone contrast for both conditions is close to 36%. Figure 10B shows the response of this cell to L − M modulation (gray diamonds), in which the M and L cone contrast is close to 18% but the cones are activated in opposite temporal phase. The solid line shows the vector (phase and amplitude) sum of the M and L cone stimuli from Figure 10A, scaled for cone contrast. The L − M response amplitude falls close to the vector sum prediction for apertures below 1° diameter but lies above the vector sum prediction for larger apertures. Thus, for these apertures, the response is enhanced by out-of-phase combination of two low-contrast (∼18%) components relative to the prediction given by higher-contrast (∼36%) stimuli.
Figure 10C shows that, overall, the magnitude of the enhancement effect is small; across all cells and stimulus conditions, there is tight correlation (r2 = 0.82) between predicted and actual response. This is consistent with other evidence that PC cells show close to linear combination of M and L cone inputs (Derrington et al., 1984; Smith et al., 1992; Reid and Shapley, 2002). The effect is greatest at apertures close to 1° (Fig. 10D). There are at least two simple explanations for the nonlinear relationship between opponent purity and opponent response amplitude. First, the combination of low-contrast components of the L − M stimulus produces negligible response saturation at a point before the opponent site (Fig. 10D, inset). Second, response rectification at high contrasts could lead to underestimation of response amplitude by the fundamental Fourier component. More comprehensive tests of these hypotheses require systematic variation of cone contrast and maintained response amplitude.
Discussion
Is the random wiring hypothesis still tenable?
When the entire PC population is considered, one key prediction of the random wiring hypothesis is supported by our data. That is, cone-opponent responses are degraded by mixed cone inputs to the receptive field center. However, our data also show that functional input from a single (L or M) cone to the receptive field center may be sufficient, but is not necessary, for cone opponency. Furthermore, our results do not support the prediction that the inhibitory input to PC cells is constituted of an average of the local L/M cone ratio. This suggests that, in the absence of selective wiring for the receptive field center mechanism, additional factors operate to bias the functional segregation of cone inputs to PC cells. In the following, we suggest that a “random wiring, functional bias” model may give a better account of PC cell responses than a purely random model.
Our foveal data are consistent with previous studies of marmoset foveal receptive fields, as well as results from four species of macaque (Macacca fascicularis, M. mulatta, M. iris, and M. nemestrina) and two species of New World monkeys (Cebus apella and Alouatta caraya). These results were obtained using several different techniques: temporal and/or spatial analysis of cone signatures (Derrington et al., 1984; Kaplan et al., 1989; Smith et al., 1992; Kremers et al., 1998; Lankheet et al., 1998a), pseudorandom (“M-sequence”) modulation of gratings or cone-isolating checkerboards (Gielen et al., 1982; Benardete and Kaplan, 1999; Reid and Shapley, 2002), and counterphase modulation of bipartite fields (Lee et al., 1998). All of these studies are consistent with the presence of cone-selective inputs to the surround mechanism of opponent PC cells. In common with previous recordings from macaque retina and geniculate (Wiesel and Hubel, 1966; DeMonasterio et al., 1975; Dreher et al., 1976), we found a small number of non-opponent PC cells with receptive fields below 5° eccentricity (4 of 13), although this is well within the eccentricity range (<10°) of one-to-one connections between midget bipolar and midget ganglion cells in marmoset (Wilder et al., 1996; Chan et al., 2001). We never encountered non-opponent PC cells under 2° eccentricity, suggesting that the PC cells in the centralmost fovea may be exclusively opponent, but our sample is too small to draw a firm conclusion.
Our data from peripheral retina are primarily consistent with results obtained from PC ganglion cells in peripheral macaque retina (Martin et al., 2001; Solomon et al., 2005). In contrast, Diller et al. (2004) measured light-evoked responses of macaque ganglion cells in vitro and found little sign of cone opponency in PC cells at eccentricities above 30°. It was later shown that the high (∼10 Hz) modulation frequency of stimulation used by these authors would not detect opponent PC cells in the peripheral retina, in which opponent responses are limited to low temporal frequencies (Solomon et al., 2005). Despite this limitation, the study by Diller et al. and by others in macaque retina (DeMonasterio and Gouras, 1975; DeMonasterio et al., 1975; Diller et al., 2004; Solomon et al., 2005) do agree that the proportion of opponent PC cells declines at eccentricities above 40°.
Determinants of chromatic gain
Early studies categorized PC cells as possessing either non-opponent or opponent receptive fields (DeValois et al., 1966; Wiesel and Hubel, 1966; DeMonasterio and Gouras, 1975), but this dichotomy was questioned by later studies that measured the relative cone weights of inputs to PC cells (Derrington et al., 1984; Lee et al., 1987; Smith et al., 1992). Our results suggest a simple explanation for this apparent discrepancy. Because chromatic response amplitude and opponent purity show an approximately exponential relationship, measurements of chromatic response amplitude (that is, the y-axis in Fig. 8C) would naturally yield a “cluster” of opponent cells, whereas measurement of relative cone weight (that is, the x-axis in Fig. 8C) would yield a broad distribution, consistent with a continuum of opponent response property. For most PC cells in our dataset, the cone bias of the center is normally opposite to that in the surround (Figs. 7, 8C). Although the reason for this differential bias is not clear, it should nevertheless work to preserve chromatic opponent responses, even as the size of PC receptive grows larger in the peripheral retina.
Source of functional bias in excitation
Direct visualization of the cone mosaic has shown for humans and macaque monkeys that the L and M cones are distributed randomly and that there is substantial intersubject variation in the average L/M cone ratio for human observers (Roorda et al., 2001; Hofer et al., 2005). The L and M mosaic in marmosets is governed by the process of X-inactivation or “Lyonization” (Lyon, 1972; Travis et al., 1988; Tovée, 1993) and comprises patches made up of clonal descendents of cells expressing either the maternal or paternal photopigment. The size of these patches in peripheral marmoset retina is not known, but limited microphotometric recordings from the fovea are consistent with a random arrangement (Bowmaker et al., 2003). Thus, different mechanisms determine the L/M mosaic in marmosets and macaques (Hunt et al., 1993; Tovée, 1993; Nathans, 1999; Wang et al., 1999), yet the processes both result in a matrix that is much closer to random arrangement than, say, the crystalline chromatic lattice of cones in the retina of fishes (Lythgoe and Partridge, 1989). Anatomical analysis has shown directly that marmoset PC cells are not selectivity wired to a subset of bipolar cells in the dendritic field (Jusuf et al., 2006a,b), and data from other primate retinas are consistent with random wiring at the level of ganglion cell dendrites (Dacey, 1993; Yamada et al., 1996; Calkins and Sterling, 1999). This all suggests that the high center cone purity in opponent PC cells arises from functional input from a single “patch” of like-type (L or M) cones in the periphery rather than selective wiring of ganglion cells for individual members of the L or M cone mosaic, as previously speculated (Dacey, 1993; Martin et al., 2001).
In many other respects studied so far, the anatomical and physiological properties of marmoset and macaque visual system are quantitatively comparable (Yeh et al., 1995a; Ghosh et al., 1996; White et al., 1998; Kilavik et al., 2003), and behavioral tests show that trichromatic New World monkeys display chromatic color discrimination capacity (Jacobs et al., 1987; Tovée et al., 1992). We cannot rule out the possibility that New and Old World primates have achieved post-receptoral cone selectivity by different routes, but it is reasonable to assume that our results would apply to the PC pathway in Old World primates. Indeed, our results can reconcile apparent contradictions in previous studies in macaque by showing that chromatic opponency depends on an overall imbalance of L and M inputs to center and surround rather than exclusive input from one cone type to either region.
Source of functional bias in inhibition
Although our results are broadly consistent with random wiring of inputs to the receptive field center mechanism, the high cone purity in PC receptive field surrounds (Fig. 7) and the negative correlation between center and surround cone weight (Fig. 8A) suggest that inhibitory mechanisms, and their interaction at the cone-opponent site, are at least as important as the properties of the center mechanism in generating opponent PC receptive fields.
We measured surround weights using an annular stimulus, so our measurements do not include the contribution of the inhibitory mechanism to the receptive field center (Fig. 5F,H). Although under the Gaussian receptive field model this contribution is very small (<5%), an in vitro analysis of the H1 class of horizontal cell receptive fields was consistent with a sensitivity profile, which shows a sharper peak than a Gaussian profile (Packer and Dacey, 2002). The fact that H1 cells contact both L and M cones (Dacey et al., 1996; Goodchild et al., 1996) implies that this inhibitory contribution to the centralmost cones in the receptive field would thus be mixed. Regardless of this possibility, our data show that, for the majority of PC cells, cone-selective inhibition can be elicited by annuli that are spatially remote from the receptive field center. This is existent proof for functional cone bias of inhibitory inputs to PC cells, but there may be some mixed L/M cone inhibition in which the inhibitory mechanism overlaps with the center. In summary, the local inhibitory influence of horizontal cells may contribute to both spatial and chromatic opponency yet may not be the critical locus for cone-selective inhibition. A network involving amacrine cells or amacrine–ganglion cell interactions becomes an equally likely candidate (Mills and Sperling, 1990; Taylor, 1999; Flores-Herr et al., 2001; Kenyon et al., 2004), but the details of this anatomical wiring are not established. Our results suggest that, whatever mechanism is at work, it is most likely common to central and mid-peripheral retina.
Footnotes
-
This work was supported by Australian National Health and Medical Research Council Grant 253621 and Australian Research Council Grant A00104053. We thank A. Lara and D. Matin for technical assistance, T. FitzGibbon and P. Jusuf for assistance with recordings, J. Forte, U. Grünert, S. Solomon, and T. R. Vidyasagar for helpful comments, and U. Grünert, P. Jusuf, and S. C. S. Lee for the drawing shown in Figure 4.
- Correspondence should be addressed to Paul R. Martin, National Vision Research Institute of Australia, Corner of Keppel and Cardigan Streets, Carlton, Victoria 3053, Australia. prmartin{at}unimelb.edu.au