True S-cones are concentrated in the ventral mouse retina and wired for color detection in the upper visual field
Abstract
Color, an important visual cue for survival, is encoded by comparing signals from photoreceptors with different spectral sensitivities. The mouse retina expresses a short wavelength-sensitive and a middle/long wavelength-sensitive opsin (S- and M-opsin), forming opposing, overlapping gradients along the dorsal-ventral axis. Here, we analyzed the distribution of all cone types across the entire retina for two commonly used mouse strains. We found, unexpectedly, that ‘true S-cones’ (S-opsin only) are highly concentrated (up to 30% of cones) in ventral retina. Moreover, S-cone bipolar cells (SCBCs) are also skewed towards ventral retina, with wiring patterns matching the distribution of true S-cones. In addition, true S-cones in the ventral retina form clusters, which may augment synaptic input to SCBCs. Such a unique true S-cone and SCBC connecting pattern forms a basis for mouse color vision, likely reflecting evolutionary adaptation to enhance color coding for the upper visual field suitable for mice’s habitat and behavior.
eLife digest
Many primates, including humans, can see color better than most other mammals. This difference is due to the variety of light-detecting proteins – called opsins – that are produced in the eye by cells known as cones. While humans have three, mice only have two different opsins, known as S and M, which detect blue/UV and green light, respectively. Mouse cones produce either S-opsins, M-opsins or both. Fewer than 10 percent of cone cells in mice produce just the S-opsin, and these cells are essential for color vision.
Mice are commonly used in scientific research, and so their vision has been well studied. However, previous research has produced conflicting results. Some studies report that cone cells that contain only S-opsin are evenly spread out across the retina. Other evidence suggests that color vision in mice exists only for the upper field of their vision, in other words, that mice can only distinguish colors that appeared above them.
Nadal-Nicolás et al. set out to understand how to reconcile these contrasting findings. Molecular tools were used to detect S- and M-opsin in the retina of mice and revealed large differences between the lower part, known as the ventral retina, and the upper part, known as the dorsal retina. The ventral retina detects light coming from above the animal, and about a third of cone cells in this region produced exclusively S-opsin, compared to only 1 percent of cones in the dorsal retina. These S-opsin cone cells in the ventral retina group into clusters, where they connect with a special type of nerve cells that transmit this signal.
To better understand these findings, Nadal-Nicolás et al. also studied albino mice. Although albino mice have a different distribution of S-opsin protein in the retina, the cone cells producing only S-opsin are similarly clustered in the ventral retina. This suggests that the concentration of S-opsin cone cells in the ventral retina is an important feature in mouse sight. This new finding corrects the misconception that S-opsin-only cone cells are evenly spread throughout the retina and supports the previous evidence that mouse color vision is greatest in the upper part of their field of vision. Nadal-Nicolás et al. suggest this arrangement could help the mice to detect predators that may attack them from above during the daytime. Together, these new findings could help to improve the design of future studies involving vision in mice and potentially other similar species.
Introduction
Topographic representation of the visual world in the brain originates from the light-sensitive photoreceptors in the retina (Rhim et al., 2017). Although the neuronal architecture of the retina is similar among different vertebrates, the numbers and distributions of photoreceptors vary considerably (Hunt and Peichl, 2014). Such patterns have been evolutionarily selected, adapting to the animal’s unique behavior (diurnal or nocturnal) and lifestyle (prey or predator) for better use of the visual information in the natural environment (Dominy and Lucas, 2001; Gerl and Morris, 2008; Peichl, 2005). Color, an important visual cue for survival, is encoded by comparing signals carried by photoreceptors with different spectral preferences (Baden and Osorio, 2019). While amongst mammals, trichromatic color vision is privileged for some primates (Jacobs et al., 1996; Nathans et al., 1986; Yokoyama and Yokoyama, 1989), most terrestrial mammals are dichromatic (Marshak and Mills, 2014; Puller and Haverkamp, 2011; Jacobs, 1993). The mouse retina expresses two types of cone opsins, S- and M-opsin, with peak sensitivities at 360 nm and 508 nm, respectively (Jacobs et al., 1991; Nikonov et al., 2006). The expression patterns of these two opsins form opposing and overlapping gradients along the dorsal-ventral axis, resulting in a majority of cones expressing both opsins (herein either ‘mixed cones’ or M+S+) (Applebury et al., 2000; Ng et al., 2001; Wang et al., 2011). Thus, S-opsin enrichment in the ventral retina better detects short-wavelength light from the sky, and M-opsin in the dorsal retina perceives the ground (e.g., a grassy field) (Baden et al., 2013; Gouras and Ekesten, 2004; Osorio and Vorobyev, 2005; Szél et al., 1992), while co-expression of both opsins (herein either mixed cones or M+S+) (Röhlich et al., 1994) broadens the spectral range of individual cones and improves perception under varying conditions of ambient light (Chang et al., 2013).
This unusual opsin expression pattern poses a challenge for color-coding, particularly so for mixed cones. However, it has been discovered that a small population of cones only expresses S-opsin (‘true S-cones’, or S+M-). These true S-cones are thought to be evenly distributed across the retina (Franke et al., 2019; Haverkamp et al., 2005; Szatko et al., 2019; Wang et al., 2011) and to be critical for encoding color, especially in the dorsal retina where they are quasi-evenly distributed in a sea of cones expressing only M-opsin (M+S-), a pattern akin to mammalian retinas in general (Haverkamp et al., 2005; Wang et al., 2011). Nonetheless, subsequent physiological studies revealed that color-opponent retinal ganglion cells (RGCs) are more abundant in the dorsal-ventral transition zone (Chang et al., 2013) and the ventral retina (Joesch and Meister, 2016). Recent large scale two-photon imaging results further demonstrated that color opponent cells were mostly located in the ventral retina (Szatko et al., 2019). Intriguingly, a behavior-based mouse study demonstrated that their ability to distinguish color is also restricted to the ventral retina (Denman et al., 2018). These results prompt us to study, at the single-cell level and across the whole retina, the spatial distributions of cone types with different opsin expression configurations and, more importantly, with regard to S-cone bipolar cell connections in order to better understand the anatomical base for the unique color-coding scheme of the mouse retina.
Results and discussion
True S-cones are highly concentrated in the ventral retina of pigmented mouse
In mouse retina, the gradients of S- and M-opsin expression along the dorsal-ventral axis have been well documented (Figure 1A–B; Applebury et al., 2000; Calderone and Jacobs, 1995; Chang et al., 2013; Haverkamp et al., 2005; Jelcick et al., 2011; Lyubarsky et al., 1999; Ortín-Martínez et al., 2014; Szél et al., 1992; Wang et al., 2011), but the distribution of individual cone types with different combinations of opsin expression across the whole retina has not been characterized (but see Baden et al., 2013; Eldred et al., 2020), which we discuss below). We developed a highly reliable algorithm to automatically quantify the different opsins (S and M) and cone types (M+S-, true S, and mixed cones, Figure 2, Figure 2—figure supplement 1) based on high-resolution images of entire flat-mount retinas immunolabeled with S- and M-opsin antibodies (Figure 2—figure supplement 1). As demonstrated in examples of opsin labeling from dorsal, medial, and ventral retinal areas of the pigmented mouse (Figure 1B, left), while M opsin-expressing cones (M+: M+S+ + M+S-) were relatively evenly distributed across three regions, S opsin-expressing cones (S+: M+S+ + S+M-) showed considerable anisotropy, with a high density in the ventral retina and a precipitous drop in the dorsal retina, confirming previous observations (Haverkamp et al., 2005; Jelcick et al., 2011; Ortín-Martínez et al., 2014). Surprisingly, instead of finding an even distribution of true S-cones as previously presumed (Baden et al., 2013; Haverkamp et al., 2005; Wang et al., 2011), we found the ventral region had much more numerous true S-cones (~30% of the local cone population; Figure 1C left, Supplementary file 1A) than did the dorsal region (~1%). This result is evident from density maps of cone types from three examples of pigmented mice, showing highly concentrated true S-cones in the ventral retina (Figure 2A, left column, bottom row). In addition, M+S--cones were concentrated in the dorsal retina, whereas mixed cones dominated the medial and ventral retina (Figure 1C left and Figure 2A, left column, 4th and 5th rows).
Despite the vast difference in S-opsin expression pattern, the distribution of true S-cones is strikingly similar between the pigmented and albino mouse
Such a highly skewed distribution of true S-cones conflicts with the general notion that true S-cones only account for ~5% of cones and are evenly distributed across the mouse retina (Baden et al., 2013; Franke et al., 2019; Haverkamp et al., 2005; Szatko et al., 2019; Wang et al., 2011); however, it is not unprecedented considering the diverse S-cone patterns seen in mammals (Ahnelt et al., 2000; Ahnelt and Kolb, 2000; Calderone et al., 2003; Hendrickson et al., 2000; Hendrickson and Hicks, 2002; Kryger et al., 1998; Müller and Peichl, 1989; Nadal-Nicolás et al., 2018; Ortín-Martínez et al., 2014; Ortín-Martínez et al., 2010; Peichl, 2005; Schiviz et al., 2008; Szél et al., 2000). Therefore, we also examined an albino mouse line to determine whether this observation persists across different mouse strains. Overall, albino retinas had slightly smaller cone populations (Figure 2B, Supplementary file 1B; Ortín-Martínez et al., 2014). Interestingly, while M-opsin expressing cones had similar distributions in both strains, S-opsin expression extended well into the dorsal retina of the albino mouse, exhibiting a greatly reduced gradient of S-opsin expression toward the dorsal retina compared to that seen in pigmented mice (Figure 1B–C, Figure 2A second row; Applebury et al., 2000; Ortín-Martínez et al., 2014). Consequently, most cones in the dorsal retina were mixed cones, and M+S- cones were very sparse (7%, compared to 97% in pigmented mouse, Figure 1C right, Supplementary file 1A, Figure 2A right). However, despite these differences, the percentage and distribution of true S-cones were remarkably conserved between strains. In both strains, true S-cones were extremely sparse in the dorsal retina (1%) but highly concentrated in the ventral retina (33% vs 29%, Figure 1C and Supplementary file 1A). Notably, the density maps of true S-cones are nearly identical in both strains (Figure 2A, bottom row). Evaluating the distribution of three main cone populations (mixed, M+S-, and true S-cone) in four retinal quadrants centered upon the optic nerve head reveals different profiles between pigmented and albino strain for mixed and M+S- cones (Figure 2C). For example, in the dorsotemporal (DT) quadrant, we observed an increase of M+S- cones from the center to the periphery (green line) in pigmented mice, compared to a majority of mixed cones (gray line) in albino mice. However, true S-cone profiles (magenta lines) were similar between the two strains in all quadrants, except for a slightly increased density along the edge of the ventronasal (VN) quadrant in pigmented mice. A recent study successfully modeled cone opsin expression and type determination according to graded thyroid hormone signaling in a pigmented mouse strain (C57BL/6) (Eldred et al., 2020). It would be interesting to see whether a different pattern of thyroid hormone and/or receptor distribution could recapitulate a similar true S-cone distribution with a very different form of S-opsin expression.
S-cone bipolar cells exhibit a dorsal-ventral gradient with a higher density in the ventral retina
One major concern regarding cone classification based on opsin immunolabeling is that some S+M- cones may instead be mixed cones with low M-opsin expression (Applebury et al., 2000; Baden et al., 2013; Nikonov et al., 2006; Röhlich et al., 1994). Even though a similar cone-type distributions have been observed in mouse retina, it has been assumed that only a fraction of the S+M- cones are ‘true’ S-cones (Baden et al., 2013; Eldred et al., 2020). Out of caution, S+M- cones were only referred to as ‘anatomical’ S-cones due to a lack of confirmation regarding their bipolar connections (Baden et al., 2013). Thus, both true S-cones and S-cone bipolar cells have been generally acknowledged to be evenly distributed across the retina (Haverkamp et al., 2005; Wang et al., 2011; Baden et al., 2013; Szatko et al., 2019; Franke et al., 2019; Eldred et al., 2020). In order to confirm the distribution of true S-cones, it is critical to uncover the distribution and dendritic contacts of S-cone bipolar cells (type 9, or SCBCs). Previously, SCBCs have only been identified among other bipolar, amacrine and ganglion cells in a Thy1-Clomeleon mouse line, rendering the quantification of their distribution across the entire retina impractical (Haverkamp et al., 2005). We generated a Copine9-Venus mouse line, in which SCBCs are specifically marked (Figure 2—figure supplement 1C), owing to the fact that Cpne9 is an SCBC-enriched gene (Shekhar et al., 2016). In retinal sections, these Venus+ bipolar cells have axon terminals narrowly ramified in sub-lamina 5 of IPL (Figure 3A), closely resembling type 9 BCs as identified in EM reconstructions (Behrens et al., 2016; Stabio et al., 2018a). In flat-mount view, these bipolar cells are often seen to extend long dendrites to reach true S-cones, bypassing other cone types (Figure 3B–C). The majority of dendritic endings formed enlarged terminals beneath true S-cones pedicles (Figure 3C–c’), but occasional slender ‘blind’ endings were present (arrow in Figure 3C–c”), which have been documented for S-cone bipolar cells in many species (Haverkamp et al., 2005; Herr et al., 2003; Kouyama and Marshak, 1992). Unexpectedly, we found that the distribution of SCBCs was also skewed toward VN retina, albeit with a shallower gradient (Figure 3D–E). To examine the connections between true S-cones and SCBCs, we immunolabeled S- and M-opsins in Copine9-Venus mouse retinas. Because M-opsin antibody signals did not label cone structures other than their outer segments, we first identified true S-cones at the outer segment level and then traced S-opsin labeling to their pedicles in the outer plexiform layer (OPL), where they connect with SCBCs (Figure 3C, for more details see material and methods). Although convergent as well divergent connections were found between true S-cones and SCBCs in both dorsal and ventral retina (see the source data), we noted different connectivity patterns. While in the dorsal retina, a single true S-cone connected to approximately 4 SCBCs (3.8 ± 0.2, see material and methods), in the ventral retina, a single SCBC contacted approximately 5 true S-cones (4.6 ± 0.4; Figure 3C, Supplementary file 2). These results agree well with the true S-cone to SCBC ratios calculated from cell densities in the DT and VN retina. Specifically, in the dorsal retina, the true S-cone to SCBC ratio was approximately 1:3.6, compared to 5.3:1 in the ventral retina (Supplementary file 3). Accordingly, both data sets support the presence of a prevalent divergence of true S-cone to SCBC connections in the dorsal retina, in comparison to a prominent convergence of contacts from true S-cones to SCBCs in VN retina. Critically, the specificity of wiring from true S-cones to SCBCs also confirms the identity of true S-cones as revealed by opsin labeling and further supports the finding that true S-cones are highly concentrated in VN mouse retina.
True S-cones in the ventral retina are not evenly distributed but form clusters
As demonstrated above, in the mouse retina, despite a large population of mixed cones, SCBCs precisely connect with true S-cones, preserving this fundamental mammalian color circuitry motif (Behrens et al., 2016; Breuninger et al., 2011; Haverkamp et al., 2005; Mills et al., 2014). However, the increased density of SCBCs in the ventral retina does not match that of true S-cones (compare Figure 3D and Figure 2a, last row). Thus, individual SCBCs in the ventral retina may be required to develop more dendrites to maximize the number of contacts made with different S-cone terminals (Supplementary file 2, graphs in Figure 3C). Intriguingly, we discovered in both strains that true S-cones in the ventral retina appeared to cluster together rather than forming an even distribution, as revealed by K-nearest neighbor analysis (Figure 4A–B, Supplementary file 2). Ideally, such true S-cone clustering may increase the availability of targets for individual SCBCs in a reduced space.
To quantify the spatial patterning of true S-cone populations (or their lack thereof), we compared the observed true S-cone distributions within 1 mm diameter VN and DT retinal samples to artificially generated alternative populations (Figure 4C). To this end, we considered two extreme patterning rules: First, one in which the space between true S-cone locations was maximized within the set of actual locations for all cones, creating a relatively uniform (evenly ‘distributed’) mosaic of true S-cones. At the other extreme, cone identities were permuted randomly (‘shuffled’) among observed cone locations (Figure 4C). Repetition of these algorithms generated distributions of patterning metrics for true S-cones (see below) that remain constrained by the observed cone locations and proportions of cone types for each 1 mm sample.
To quantitatively compare the patterning of real true S-cone populations to their artificial counterparts, we first computed two measures of regularity for true S-cones: nearest neighbor and Voronoi diagram regularity indices (NNRI and VDRI, respectively; Reese and Keeley, 2015; Figure 4C–D); larger values of these metrics indicate smaller variability in the spacing between cones and thus more regular patterns. Interestingly, far from being regularly distributed, true S-cone placement was quite irregular and nearly indistinguishable from shuffled populations (including a slight trend toward regularity measures lower than random, which may indicate a tendency toward clustering, Figure 4D; see Reese, 2008). To further probe the possibility of true S-cone clustering, we computed the ratios of true S-cone neighbors for each cone (denoted here as the S-cone neighbor ratio [SCNR]; see Methods for the calculation of the SCNR search radius for each retinal sample). Intriguingly, SCNRs were significantly larger for true S-cones than for other cone types, which were equal to expected ratios due to random chance—especially so in ventral retinas, further indicating a clustering of true S-cones in those areas (Figure 4E). Notably, a more extreme form of clustering of S-cones has been observed in the ‘wild’ mouse (Warwick et al., 2018) and with much lower densities in some felids (Ahnelt et al., 2000). Here, such clustering may reflect the mode of true S-cone development in the ventral retina, for example, by ‘clonal expansion’ to achieve unusually high densities (Bruhn and Cepko, 1996; Reese et al., 1999). It is tempting to speculate that it may also facilitate the wiring of true S-cones with sparsely distributed SCBCs, which were not observed to cluster in the ventral retina (Figure 3E). Indeed, we observed examples of groups of true S-cones forming clusters whose pedicles in the OPL were tightly congregated in a patch and contacted by a nearby SCBC (Figure 4F).
Enriched true S-cones in the ventral retina may provide an anatomical base for mouse color vision
Despite being nocturnal and having a rod-dominated retina (Carter-Dawson and LaVail, 1979; Jeon et al., 1998), mice can detect color (Denman et al., 2018; Jacobs et al., 2004). Although it remains uncertain whether the source of long-wavelength sensitive signals for color opponency arises in rods or M-cones (Baden and Osorio, 2019; Ekesten et al., 2000; Ekesten and Gouras, 2005; Joesch and Meister, 2016; Reitner et al., 1991), it is clear that true S-cones provide short-wavelength signals for color discrimination. Given the previously-held notion that true S-cones are evenly distributed across the retina (Baden et al., 2013; Franke et al., 2019; Haverkamp et al., 2005; Szatko et al., 2019; Wang et al., 2011), whereas M+S- cones are concentrated in the dorsal retina of pigmented mouse, it is intuitive to speculate that color coding is prevalent in the dorsal retina. However, previous physiological and behavioral studies indicate that, although luminance detection can occur across the mouse retina, color discrimination is restricted to the ventral retina (Breuninger et al., 2011; Denman et al., 2018; Szatko et al., 2019). Thus, our discovery of high enrichment of true S-cones in the ventral retina provides a previously missed anatomical feature for mouse color vision that could help to re-interpret these results. From projections mapping true S-cone densities into visual space (Figure 4—figure supplement 1; Sterratt et al., 2013), it is conceivable that high ventral true S-cone density will provide a much higher sensitivity of short-wavelength signals, thus facilitating color detection for the upper visual field. Although the true S-cone signals carried by SCBCs in the dorsal retina might not be significant for color detection, they could certainly participate in other functions, such as non-image forming vision, that are known to involve short-wavelength signals (Altimus et al., 2008; Doyle et al., 2008; Patterson et al., 2020). Interestingly, the overall true S-cone percentage in the mouse retina remains approximately 10% (Figure 2B), and the average true S-cone to SCBC ratio across the whole retina is about 1.7:1 (Supplementary file 1B-C), similar to what has been reported in other mammals (Ahnelt et al., 2006; Ahnelt and Kolb, 2000; Bumsted et al., 1997; Bumsted and Hendrickson, 1999; Curcio et al., 1991; Hendrickson and Hicks, 2002; Hunt and Peichl, 2014; Kryger et al., 1998; Lukáts et al., 2005; Müller and Peichl, 1989; Ortín-Martínez et al., 2010; Peichl et al., 2000; Schiviz et al., 2008; Shinozaki et al., 2010; Szél et al., 1988).
Such a spatial rearrangement of true S-cones and SCBCs likely reflects evolutionary adaptation to enhance short-wavelength signaling and color coding for the upper visual field as best suited for the habitat and behavior of mice (Baden et al., 2020). For example, it may facilitate aerial predator detection during daytime (Yilmaz and Meister, 2013). Similarly, skewed S-cone arrangement has been reported for other terrestrial prey mammals (Famiglietti and Sharpe, 1995; Juliusson et al., 1994; Röhlich et al., 1994), while zebrafish possess a UV-enriched ventral retina that enhances their predation (Zimmermann et al., 2018). In addition, we observed that the clustering of true S-cones in the ventral retina may allow several neighboring cones of the same type to converge onto the same SCBC (Figure 4F), which could potentially enhance signal-to-noise ratios for more accurate detection, as described recently in human fovea (Schmidt et al., 2019). It is also remarkable that despite the very different S-opsin expression patterns in both mouse strains, the true S-cone population and distribution are strikingly similar between pigmented and albino mice, suggesting a common functional significance.
Materials and methods
Animal generation, handling and ethic statement
Request a detailed protocolThree months old male pigmented (C57BL/6J, n = 5), albino (CD1, n = 5) mice were obtained from the National Eye Institute breeding colony. The Venus-Cpne9 mouse line (n = 5; based on previous single cell sequencing data [Shekhar et al., 2016]) carries a reporter (Venus) allele under the control of the mouse Cpne9 locus. The reporter allele was created directly in B6.SJL(F1) zygotes using CRISPR-mediated homologous recombination (HR) (Yang et al., 2013). Briefly, a HR targeting template was assembled with PCR fragments of 5’ and 3’ homology arms of 910 bp and 969 bp respectively, flanking exon one, and a Venus expression cassette carrying the bovine growth hormone polyadenylation (bGH-PolyA) signal sequence as the terminator. Homology arms were designed such that integration of the reporter cassette would be at the position right after the first codon of the Cpne9 gene in exon one. A pair of guide RNAs (gRNA), with outward orientation (38 bp apart), were synthesized by in vitro transcription as described (Yang et al., 2013) and tested for their efficiency and potential toxicity in a zygote differentiation assay where mouse fertilized eggs were electroporated with SpCas9 protein and gRNA ribonuclear particles. Eggs were cultured in vitro for 4 days in KSOM (Origio Inc, CT) until differentiated to blastocysts. Viability and indel formation were counted respectively. gRNA sequences are (1) Copine9_gRNA_L(73/25), 5’GAGACATGACTGGTCCAA3’; (2) Copine9_gRNA_R(62/4.40), 5’GCCTCGGAGCGTAGCGTCC3’. A mixture of the targeting plasmid (super coiled, 25 ng/µl) with two tested gRNAs (25 ng/µl each) and the SpCas9 protein (Life Science technology, 30 ng/µl) were microinjected into mouse fertilized eggs and transferred to pseudopregnant female recipients as described elsewhere (Yang et al., 2013). With a total of 15 F0 live births from 6 pseudopregnant females, 11 were found to carry the knockin allele by homologous recombination, a HR rate of 73%. F0 founders in B6.SJL F1 (50% C57BL6 genome) were crossed consecutively for 3 generations with C57BL6/J mice to reach near congenic state to C57BL6/J.
Mice were housed a 12:12 hr light/dark cycle. All experiments and animal care are conducted in accordance with protocols approved by the Animal Care and Use Committee of the National Institutes of Health and following the Association for Research in Vision and Ophthalmology guidelines for the use of animals in research.
Tissue collection
Request a detailed protocolAll animals were sacrificed with an overdose of CO2 and perfused transcardially with saline followed by 4% paraformaldehyde. To preserve retinal orientation, eight retinas per mouse strain/line were dissected as flat whole-mounts by making four radial cuts (the deepest one in the dorsal pole previously marked with a burn signal as described [Nadal-Nicolás et al., 2018; Stabio et al., 2018b]). The two remaining retinas were cut in dorso-ventral orientation (14 μm) after cryoprotection in increasing gradients of sucrose (Sigma-Aldrich SL) and embedding in optimal cutting temperature (OCT; Sakura Finetek).
Immunohistochemical labeling
Request a detailed protocolImmunodetection of flat-mounted retinas or retinal sections was carried out as previously described (Nadal-Nicolás et al., 2018). Importantly, the retinal pigmented epithelium was removed before the immunodetection. First, whole-retinas were permeated (4 × 10’) in PBS 0.5% Triton X-100 (Tx) and incubated by shaking overnight at room temperature with S-opsin (1:1200) and M-opsin (1:1000) or cone arrestin (1:300) primary antibodies diluted in blocking buffer (2% normal donkey serum). Cpne9-Venus retinas were additionally incubated with an anti-GFP antibody (1:100) to enhance the original Venus signal. Retinas were washed in PBS 0.5% Tx before incubating the appropriate secondary antibodies overnight (1:500). Finally, retinas were thoroughly washed prior to mounting with photoreceptor side up on slides and covered with anti-fading solution. Retinal sections were counterstained with DAPI.
Image acquisition
Request a detailed protocolRetinal whole-mounts were imaged with a 20x objective using a LSM 780 Zeiss confocal microscope equipped with computer-driven motorized stage controlled by Zen Lite software (Black edition, Zeiss). M- and S-opsins were imaged together to allow the identification and quantification of different cone types. Magnifications from flat mounts and retinal cross-sections (Figure 1) were taken from dorsal, medial and ventral areas using a 63x objective for opsin co-expression analysis. Images from retinal cross-sections were acquired ~1.5 mm dorsally or ventrally from the optic disc.
Sampling and opsin co-expression measurement
Request a detailed protocolIn four retinas per strain, we acquired images from three 135 × 135 μm samples (63x) per each area of interest (dorsal, medial and ventral). These areas were selected according to the S-opsin gradient in wholemount retinas (see scheme in Figure 1C). Cone outer segments were manually classified as M+S-, true S- (S+M-) or mixed (M+S+) cones depending on their opsin expression. Data representation was performed using GraphPad Prism 8.3 software.
Image processing: manual and automated whole quantification
Request a detailed protocolTo characterize the distribution of the different cone photoreceptor types in the mouse retina, we developed and validated an automatic routine (ImageJ, NIH) to identify, quantify the total number of outer segments and finally extract the location of each individual cone (Figure 2—figure supplement 1A). Briefly, maximum-projection images were background-subtracted and thresholded (background-noise mean value, 9.6 ± 1.2% and 15.2 ± 3.2% for S- and M-opsin respectively, the threshold was applied at 15.7%) to create a binary mask that was then processed using watershed and despeckle filters to isolate individual cones and reduce noise. The ‘3D Objects Counter’ plugin was applied to such images to count cones within fixed parameters (shape and size) and extract their xy coordinates for further analysis. This automation was validated by statistical comparison with manual counting performed by an experienced investigator (Pearson correlation coefficient R2 = 96–99% for M- or S-opsin respectively, Figure 2—figure supplement 1B). To count cone subtypes, images were pre-processed with image processing software (Adobe Photoshop CC) to isolate the desired subtype and then manually marked using Photoshop, or automatically counted using ImageJ as described above. Total cone populations were determined by combining M- and S-opsin channels, while mixed M+S+ cones were obtained by masking the M-opsin signal with the S-opsin channel. M+S- cones in pigmented mice were obtained by subtracting the S-opsin signal from the M-opsin photomontage. Finally, M+S- cones (in albino samples), true S-cones (both strains) (Figure 2—figure supplement 1C) and Venus+ SCBCs (Cpne9-Venus mouse line) were manually marked on the retinal photomontage (Adobe Photoshop CC).
Topographical distributions
Request a detailed protocolTopographical distributions of cone population densities were calculated from cone locations identified in whole-mount retinas using image processing (see above). From these populations, isodensity maps were created using Sigmaplot 13.0 (Systat Software). These maps are filled contour plots generated by assigning to each area of interest (83.3 × 83.3 μm) a color code according to its cone density, ranging from 0 (purple) to 17,300 cones/mm2 for all cone types except for true S-cones and M+S--cone in the albino strain (5000 cones/mm2), as represented in the last image of each row of Figure 2A, or 1,400 SCBCs/mm2 (Figure 3D) within a 10-step color-scale. These calculations allow as well, the illustration of the number of cones at a given position from the ON center. To analyze the relative opsin expression along the retinal surface, we have considered three cone populations (mixed, M+S-- and true S-cones) dividing the retina in four quadrants: dorsotemporal, dorsonasal, ventrotemporal and ventronasal (DT, DN, VT and VN respectively, scheme in Figure 2C). The relative percentage of cone-types is represented in line graphs from four retinas/strain (SigmaPlot 13.0).
SCBC sampling and ‘true S-cone’ connectivity
Request a detailed protocolTo characterize the connectivity of Venus+ S-cone bipolar cells (Venus+ SCBCs) with true S-cone terminals, we acquired images from the same area (260 × 260 μm, 63x) at two focal planes: First, we focused upon the INL+OPL, then the corresponding photoreceptor outer segment (OS) layer, respectively, for two areas of interest (DT and VN). To verify connectivity between Venus+ SCBC dendrites and true S-cone pedicles in the OPL, in addition to S-opsin immunodetection, we also labeled retinas using cone arrestin antibodies to discriminate mixed cone pedicles from true S-cone pedicles, because true S-cone pedicles contain either low or no cone arrestin (Figure 3B, Haverkamp et al., 2005). In other retinas, SCBC contacts were verified by tracking each cell body from cone pedicles to their respective OS to confirm S+M- opsin labeling (Figure 3C). In five retinas (with S- and M-opsin double immunodetection), we analyzed the connectivity between 186 Venus+ SCBCs (133 and 53 for DT and VN respectively) and 263 true S-cone pedicles (74 and 189, DT and VN respectively). The number of synaptic contacts was assessed by tracking manually each SCBC-branch from the cell body using the Zen lite black visualization package (Z-stack with 1 μm interval). Multiple branch contacts in one true S-cone pedicle from a single SCBC were considered a single contact and counted only once (Figure 3B), while secondary bifurcations were considered as multiple contacts (Figure 3c’). SCBC-blind endings were not counted (Figure 3c”). The average number of contacts per retina was used to calculate the DT and VN means (Supplementary file 2 and graphs in 3C).
Clustering analysis
K-neighbor maps and variance analysis of Voronoi dispersion
Request a detailed protocolTo assess the true S-cones and S-cone bipolar cell (SCBC) clustering, we performed two comparable sets of analyses. First, we extracted two circular areas (1 mm diameter) in the DT-VN axis at 1 mm from the optic disc center (scheme in 4B). A K-nearest neighbor algorithm (Nadal-Nicolás et al., 2014) was used to map the number of neighboring true S-cones within a 18 μm radius of each true S-cone to a color-code in its retinal position (Figure 4B). Regularity indices were computed for each retinal sample using Voronoi diagrams for cone positions as well as nearest neighbor distances (VDRI and NNRI, respectively [Reese and Keeley, 2015; Figure 4C–E]). NNRIs were computed as the ratio of the mean to the standard deviation for the distance from true S-cones to their nearest true S-cone neighbor. true S-cone neighbor ratios (SCNR) were calculated for each retinal sample as the average proportion of true s-cones within a given radius for each cone. This search radius was calculated separately for each sample to correct for sample-to-sample variations in total density: this radius (r) was calculated as r = 3√ (A / (√2 πN)), where A is the circular area of the 1 mm diameter retinal sample and N is the total number of cones in that sample. For a highly regular cell mosaic containing N cells filling an area A, this calculation estimates the location of the first minimum in the density recovery profile (Rodieck, 1991), providing the average radius of a circle centered upon a cone that will encompass its first tier of cone neighbors (but exclude the second tier) in an evenly distributed mosaic. To minimize edge effects from computations of NNRI, VDRI, SCNR, those values for cones closer to the outer edge of the sample than the SCNR search radius were discarded. To produce simulated cone mosaics for comparison with observed values, cone distributions with evenly ‘distributed’ true S-cones were generated by first using a simple mutual repulsion simulation to maximize the distances between true S-cones, followed by assigning the nearest positions among all cone locations as being ‘true S’. ‘Shuffled’ populations of true S-cones were generated by permuting cone identities randomly among all cone locations, holding the proportion of true S-cones constant. Voronoi diagrams, neighbor calculations, and mosaic generation and other computations were performed using MATLAB R2016b.
True S-cone cluster and SCBC synaptic contacts evaluation
Request a detailed protocolTo characterize the true S-cone cluster connectivity in the VN retina, retinal whole-mounts were imaged with a 63x objective, from the photoreceptor outer segments to the OPL, in a Z-stack image with 0.5 μm interval. To visualize the true S-cone clustering and Venus+ SCBC connectivity, we identified numerically, and color coded each true S-outer segment form a cluster. The corresponding true S-pedicles were identified by tracking the cell body from their S+M-OSs. Focusing on the outer plexiform layer (OPL), each individual true S-cone pedicle -that form a cluster- was manually outlined and color coded accordingly. Lastly, the SCBC synaptic terminals that belong to a single SCBC, were identified by their specific contacts to the respective true S-cone pedicle (Figure 4F).
Retinal reconstruction and visuotopic projection
Request a detailed protocolRetinal images were reconstructed and projected into visual space using R software v.3.5.2 for 64-bit Microsoft Windows using Retistruct v.0.6.2 as in Sterratt et al., 2013. Reconstruction parameters from that study were used: namely, a rim angle of 112° (phi0 = 22°), and eye orientation angles of 22° (elevation) and 64° (azimuthal). For Figure 4—figure supplement 1, true S-cone density contour lines and heatmaps were computed in MATLAB and overlaid onto flat-mount retina opsin labeling images using ImageJ prior to processing by Retistruct.
Statistical analysis
Request a detailed protocolStatistical comparisons for the percentage of cones/retinal location, the total cone quantifications (Supplementary file 1) and the DT or VN true S-cones and Venus+SCBCs (Supplementary file 2) were carried out using GraphPad Prism v8.3 for Microsoft Windows. Data are presented as mean ± standard deviation. All data sets passed the D'Agostino-Pearson test for normality, and the comparisons between strains were performed with Student’s t-test.
For each 1 mm retinal sample, VDRI, NNRI, and SCNR values were normalized and compared to the distributions of ‘shuffled’ cone populations. Such comparisons were not performed against ‘distributed’ populations, because in those populations, VDRI and NNRI values were consistently much higher—and SCNR much lower—than in real samples (see Figure 4D–E). The ‘shuffled’ populations for each retinal region produced measurements that were well described by normal distributions (Kolmogorov-Smirnov test, MATLAB). Thus, to allow comparisons across samples, we converted each measurement into a Z-score using the mean and standard deviation of those measures from shuffled populations. One-tailed Student’s t-tests were performed to compare the normalized measures to the distribution of ‘randomly shuffled’ cone population measures, and significance was determined at the p<0.05 level.
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files.
References
-
Irregular S-cone mosaics in felid retinas spatial interaction with axonless horizontal cells, revealed by cross correlationJournal of the Optical Society of America A 17:580–588.https://doi.org/10.1364/JOSAA.17.000580
-
The mammalian photoreceptor mosaic-adaptive designProgress in Retinal and Eye Research 19:711–777.https://doi.org/10.1016/S1350-9462(00)00012-4
-
Understanding the retinal basis of vision across speciesNature Reviews Neuroscience 21:5–20.https://doi.org/10.1038/s41583-019-0242-1
-
The retinal basis of vertebrate color visionAnnual Review of Vision Science 5:177–200.https://doi.org/10.1146/annurev-vision-091718-014926
-
Chromatic bipolar cell pathways in the mouse retinaJournal of Neuroscience 31:6504–6517.https://doi.org/10.1523/JNEUROSCI.0616-11.2011
-
Development of the pattern of photoreceptors in the chick retinaThe Journal of Neuroscience 16:1430–1439.https://doi.org/10.1523/JNEUROSCI.16-04-01430.1996
-
Spatial and temporal expression of cone opsins during monkey retinal developmentThe Journal of Comparative Neurology 378:117–134.https://doi.org/10.1002/(SICI)1096-9861(19970203)378:1<117::AID-CNE7>3.0.CO;2-7
-
Distribution and development of short-wavelength cones differ between macaca monkey and human foveaThe Journal of Comparative Neurology 403:502–516.https://doi.org/10.1002/(SICI)1096-9861(19990125)403:4<502::AID-CNE6>3.0.CO;2-N
-
Topography of photoreceptors and retinal ganglion cells in the spotted hyena (Crocuta crocuta)Brain, Behavior and Evolution 62:182–192.https://doi.org/10.1159/000073270
-
Regional variations in the relative sensitivity to UV light in the mouse retinaVisual Neuroscience 12:463–468.https://doi.org/10.1017/S0952523800008361
-
Rods and cones in the mouse retina. I. structural analysis using light and electron microscopyThe Journal of Comparative Neurology 188:245–262.https://doi.org/10.1002/cne.901880204
-
Distribution and morphology of human cone photoreceptors stained with anti-blue opsinThe Journal of Comparative Neurology 312:610–624.https://doi.org/10.1002/cne.903120411
-
Retinal pathways influence temporal nichePNAS 105:13133–13138.https://doi.org/10.1073/pnas.0801728105
-
Cone inputs to murine retinal ganglion cellsVision Research 40:2573–2577.https://doi.org/10.1016/S0042-6989(00)00122-X
-
Modeling binary and graded cone cell fate patterning in the mouse retinaPLOS Computational Biology 16:e1007691.https://doi.org/10.1371/journal.pcbi.1007691
-
The causes and consequences of color visionEvolution: Education and Outreach 1:476–486.https://doi.org/10.1007/s12052-008-0088-x
-
Why do mice have ultra-violet vision?Experimental Eye Research 79:887–892.https://doi.org/10.1016/j.exer.2004.06.031
-
The primordial, blue-cone color system of the mouse retinaJournal of Neuroscience 25:5438–5445.https://doi.org/10.1523/JNEUROSCI.1117-05.2005
-
Nocturnal tarsier retina has both short and long/medium-wavelength cones in an unusual topographyThe Journal of Comparative Neurology 424:718–730.https://doi.org/10.1002/1096-9861(20000904)424:4<718::AID-CNE12>3.0.CO;2-Z
-
Distribution and density of medium- and short-wavelength selective cones in the domestic pig retinaExperimental Eye Research 74:435–444.https://doi.org/10.1006/exer.2002.1181
-
Inner S-cone bipolar cells provide all of the central elements for S cones in macaque retinaThe Journal of Comparative Neurology 457:185–201.https://doi.org/10.1002/cne.10553
-
S cones: evolution, retinal distribution, development, and spectral sensitivityVisual Neuroscience 31:115–138.https://doi.org/10.1017/S0952523813000242
-
The distribution and nature of colour vision among the mammalsBiological Reviews 68:413–471.https://doi.org/10.1111/j.1469-185X.1993.tb00738.x
-
The major cell populations of the mouse retinaThe Journal of Neuroscience 18:8936–8946.https://doi.org/10.1523/JNEUROSCI.18-21-08936.1998
-
Complementary cone fields of the rabbit retinaInvestigative Ophthalmology & Visual Science 35:811–818.
-
Bipolar cells specific for blue cones in the macaque retinaThe Journal of Neuroscience 12:1233–1252.https://doi.org/10.1523/JNEUROSCI.12-04-01233.1992
-
Photopigment coexpression in mammals: comparative and developmental aspectsHistology and Histopathology 20:551–574.https://doi.org/10.14670/HH-20.551
-
Short-wavelength cone-opponent retinal ganglion cells in mammalsVisual Neuroscience 31:165–175.https://doi.org/10.1017/S095252381300031X
-
Three distinct blue-green color pathways in a mammalian retinaJournal of Neuroscience 34:1760–1768.https://doi.org/10.1523/JNEUROSCI.3901-13.2014
-
Topography of cones and rods in the tree shrew retinaJournal of Comparative Neurology 282:581–594.https://doi.org/10.1002/cne.902820409
-
Displaced retinal ganglion cells in albino and pigmented ratsFrontiers in Neuroanatomy 8:99.https://doi.org/10.3389/fnana.2014.00099
-
The aging rat retina: from function to anatomyNeurobiology of Aging 61:146–168.https://doi.org/10.1016/j.neurobiolaging.2017.09.021
-
Physiological features of the S- and M-cone photoreceptors of wild-type mice from single-cell recordingsJournal of General Physiology 127:359–374.https://doi.org/10.1085/jgp.200609490
-
Automated quantification and topographical distribution of the whole population of S- and L-Cones in adult albino and pigmented ratsInvestigative Opthalmology & Visual Science 51:3171–3183.https://doi.org/10.1167/iovs.09-4861
-
Photoreceptor sectral sensitivities in terrestrial animals: adaptations for luminance and colour visionProceedings of the Royal Society B: Biological Sciences 272:1745–1752.https://doi.org/10.1098/rspb.2005.3156
-
Photoreceptor types and distributions in the retinae of insectivoresVisual Neuroscience 17:937–948.https://doi.org/10.1017/S0952523800176138
-
Diversity of mammalian photoreceptor properties: adaptations to habitat and lifestyle?The Anatomical Record Part A: Discoveries in Molecular, Cellular, and Evolutionary Biology 287A:1001–1012.https://doi.org/10.1002/ar.a.20262
-
Bipolar cell pathways for color vision in non-primate dichromatsVisual Neuroscience 28:51–60.https://doi.org/10.1017/S0952523810000271
-
Clonal expansion and cell dispersion in the developing mouse retinaEuropean Journal of Neuroscience 11:2965–2978.https://doi.org/10.1046/j.1460-9568.1999.00712.x
-
Book1.22 - Mosaics, Tiling, and Coverage by Retinal NeuronsIn: Masland R. H, Albright T. D, Albright T. D, Masland R. H, Dallos P, Oertel D, Firestein S, Beauchamp G. K, Catherine Bushnell M, Basbaum A. I, Kaas J. H, Gardner E. P, editors. The Senses: A Comprehensive Reference. New York: Academic Press. pp. 439–456.https://doi.org/10.1016/B978-012370880-9.00278-4
-
Design principles and developmental mechanisms underlying retinal mosaicsBiological Reviews 90:854–876.https://doi.org/10.1111/brv.12139
-
Maps of cone opsin input to mouse V1 and higher visual areasJournal of Neurophysiology 117:1674–1682.https://doi.org/10.1152/jn.00849.2016
-
Retinal cone topography of artiodactyl mammals: influence of body height and habitatThe Journal of Comparative Neurology 507:1336–1350.https://doi.org/10.1002/cne.21626
-
Topography of ganglion cells and photoreceptors in the sheep retinaThe Journal of Comparative Neurology 518:2305–2315.https://doi.org/10.1002/cne.22333
-
A novel map of the mouse eye for orienting retinal topography in anatomical spaceJournal of Comparative Neurology 526:1749–1759.https://doi.org/10.1002/cne.24446
-
Identification of the blue-sensitive cones in the mammalian retina by anti-visual pigment antibodyThe Journal of Comparative Neurology 273:593–602.https://doi.org/10.1002/cne.902730413
-
Unique topographic separation of two spectral classes of cones in the mouse retinaThe Journal of Comparative Neurology 325:327–342.https://doi.org/10.1002/cne.903250302
-
Photoreceptor distribution in the retinas of subprimate mammalsJournal of the Optical Society of America A 17:568–579.https://doi.org/10.1364/josaa.17.000568
-
Spectral and temporal sensitivity of cone-mediated responses in mouse retinal ganglion cellsJournal of Neuroscience 31:7670–7681.https://doi.org/10.1523/JNEUROSCI.0629-11.2011
-
Inhomogeneous encoding of the visual field in the mouse retinaCurrent Biology 28:655–665.https://doi.org/10.1016/j.cub.2018.01.016
-
Rapid innate defensive responses of mice to looming visual stimuliCurrent Biology 23:2011–2015.https://doi.org/10.1016/j.cub.2013.08.015
-
Molecular evolution of human visual pigment genesMolecular Biology and Evolution 6:186–197.https://doi.org/10.1093/oxfordjournals.molbev.a040537
Article and author information
Author details
Funding
National Eye Institute (Intramural Research Program)
- Wei Li
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors would like to thank the NEI Animal Care team, especially Megan Kopera and Ashley Yedlicka.
Ethics
Animal experimentation: All experiments and animal care are conducted in accordance with protocols approved by the Animal Care and Use Committee of the National Institutes of Health and following the Association for Research in Vision and Ophthalmology guidelines for the use of animals in research. The protocol was approved by the Animal Care and Use Committee of the National Institutes of Health (ASP#606).
Copyright
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Metrics
-
- 4,585
- views
-
- 689
- downloads
-
- 97
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
- Genetics and Genomics
Evolutionary arms races can arise at the contact surfaces between host and viral proteins, producing dynamic spaces in which genetic variants are continually pursued. However, the sampling of genetic variation must be balanced with the need to maintain protein function. A striking case is given by protein kinase R (PKR), a member of the mammalian innate immune system. PKR detects viral replication within the host cell and halts protein synthesis to prevent viral replication by phosphorylating eIF2α, a component of the translation initiation machinery. PKR is targeted by many viral antagonists, including poxvirus pseudosubstrate antagonists that mimic the natural substrate, eIF2α, and inhibit PKR activity. Remarkably, PKR has several rapidly evolving residues at this interface, suggesting it is engaging in an evolutionary arms race, despite the surface’s critical role in phosphorylating eIF2α. To systematically explore the evolutionary opportunities available at this dynamic interface, we generated and characterized a library of 426 SNP-accessible nonsynonymous variants of human PKR for their ability to escape inhibition by the model pseudosubstrate inhibitor K3, encoded by the vaccinia virus gene K3L. We identified key sites in the PKR kinase domain that harbor K3-resistant variants, as well as critical sites where variation leads to loss of function. We find K3-resistant variants are readily available throughout the interface and are enriched at sites under positive selection. Moreover, variants beneficial against K3 were also beneficial against an enhanced variant of K3, indicating resilience to viral adaptation. Overall, we find that the eIF2α-binding surface of PKR is highly malleable, potentiating its evolutionary ability to combat viral inhibition.
-
- Evolutionary Biology
- Genetics and Genomics
It is well established that several Homo sapiens populations experienced admixture with extinct human species during their evolutionary history. Sometimes, such a gene flow could have played a role in modulating their capability to cope with a variety of selective pressures, thus resulting in archaic adaptive introgression events. A paradigmatic example of this evolutionary mechanism is offered by the EPAS1 gene, whose most frequent haplotype in Himalayan highlanders was proved to reduce their susceptibility to chronic mountain sickness and to be introduced in the gene pool of their ancestors by admixture with Denisovans. In this study, we aimed at further expanding the investigation of the impact of archaic introgression on more complex adaptive responses to hypobaric hypoxia evolved by populations of Tibetan/Sherpa ancestry, which have been plausibly mediated by soft selective sweeps and/or polygenic adaptations rather than by hard selective sweeps. For this purpose, we used a combination of composite-likelihood and gene network-based methods to detect adaptive loci in introgressed chromosomal segments from Tibetan WGS data and to shortlist those presenting Denisovan-like derived alleles that participate to the same functional pathways and are absent in populations of African ancestry, which are supposed to do not have experienced Denisovan admixture. According to this approach, we identified multiple genes putatively involved in archaic introgression events and that, especially as regards TBC1D1, RASGRF2, PRKAG2, and KRAS, have plausibly contributed to shape the adaptive modulation of angiogenesis and of certain cardiovascular traits in high-altitude Himalayan peoples. These findings provided unprecedented evidence about the complexity of the adaptive phenotype evolved by these human groups to cope with challenges imposed by hypobaric hypoxia, offering new insights into the tangled interplay of genetic determinants that mediates the physiological adjustments crucial for human adaptation to the high-altitude environment.