Morphological Plasticity of the Retina of Viperidae Snakes Is Associated With Ontogenetic Changes in Ecology and Behavior

Snakes of the Viperidae family have retinas adapted to low light conditions, with high packaging of rod-photoreceptors containing the rhodopsin photopigment (RH1), and three types of cone-photoreceptors, large single and double cones with long-wavelength sensitive opsins (LWS), and small single cones with short-wavelength sensitive opsins (SWS1). In this study, we compared the density and distribution of photoreceptors and ganglion cell layer (GCL) cells in whole-mounted retinas of two viperid snakes, the lancehead Bothrops jararaca and the rattlesnake Crotalus durissus, and we estimated the upper limits of spatial resolving power based on anatomical data. The ground-dwelling C. durissus inhabits savannah-like habitats and actively searches for places to hide before using the sit-and-wait hunting strategy to ambush rodents. B. jararaca inhabits forested areas and has ontogenetic changes in ecology and behavior. Adults are terrestrial and use similar hunting strategies to those used by rattlesnakes to prey on rodents. Juveniles are semi-arboreal and use the sit-and-wait strategy and caudal luring to attract ectothermic prey. Our analyses showed that neuronal densities were similar for the two species, but their patterns of distribution were different between and within species. In adults and juveniles of C. durissus, cones were distributed in poorly defined visual streaks and rods were concentrated in the dorsal retina, indicating higher sensitivity in the lower visual field. In adults of B. jararaca, both cones and rods were distributed in poorly defined visual streaks, while in juveniles, rods were concentrated in the dorsal retina and cones in the ventral retina, enhancing sensitivity in the lower visual field and visual acuity in the upper field. The GCL cells had peak densities in the temporal retina of C. durissus and adults of B. jararaca, indicating higher acuity in the frontal field. In juveniles of B. jararaca, the peak density of GCL cells in the ventral retina indicates better acuity in the upper field. The estimated visual acuity varied from 2.3 to 2.8 cycles per degree. Our results showed interspecific differences and suggest ontogenetic plasticity of the retinal architecture associated with changes in the niche occupied by viperid snakes, and highlight the importance of the retinal topography for visual ecology and behavior of snakes.


INTRODUCTION
The highly diverse group of Snakes, with more than 3,800 species (Uetz et al., 2020) has a fascinating diversity of retinal morphology, especially regarding the photoreceptor types (Walls, 1942;Underwood, 1967aUnderwood, , 1970Caprette, 2005;Hauzman et al., 2017;Hauzman, 2020). This group represents a valuable model to test hypotheses of correlation between the types of retinal specialization and species ecology and behavior. The vertebrate retina is formed by layers of cells and nerve plexuses organized in a highly conserved fashion that allows the vertical flow of luminous information from the photoreceptors in the outermost retina toward the ganglion cells in the innermost retina (Ramón y Cajal, 1983). The photoreceptors contain the visual pigments that absorb photons and trigger an enzymatic cascade within the cell. The light energy is converted into neural signals that are transmitted to bipolar cells and from those to the ganglion cells (GCs) that conduct the information to the brain (Ramón y Cajal, 1983).
The density and distribution of cells in the retinas are highly variable among species. Specific regions of higher cell density, the retinal specializations, reflect areas of the visual field that have greater importance for photon uptake, spatial resolution, or other visual functions depending on the cell type (Baden et al., 2020). Two main types of retinal specializations, visual streak and area centralis, were described in many vertebrates and are usually associated with the habitat occupied by the species (Hughes, 1977;Moore et al., 2017). Visual streaks are elongated regions of higher cell density that allow a wide screening of the surroundings without the constant need for head and eye movements and are usually associated with the use of open environments. Areae centrales are concentric regions of higher cell density usually found in species that occupy closed environments such as forests, where the horizon is obstructed by vegetation (Hughes, 1977;Collin, 2008;Moore et al., 2017).
In Snakes, despite the ecological diversity of the group, a very limited number of studies investigated the organization of neurons in the retinas (Wong, 1989;Hart et al., 2012;Hauzman et al., 2014Hauzman et al., , 2018. Different types of specializations were described even among sympatric and closely related species. In marine Elapidae snakes the GCs are arranged in horizontal streaks that might enable a better view of the open ocean environment , and in two out of three marine species analyzed by Hart et al. (2012), an additional area centralis in the ventral retina was associated with specific foraging strategies. In the arboreal Dipsadidae snake Philodryas olfersii, the photoreceptors and GC are arranged in horizontal streaks, while in the close-related ground-dwelling Philodryas patagoniensis, these neurons are concentrated in a ventral area centralis, indicating better spatial resolution of the upper visual field . In a comparative study of the distribution of GCs of diurnal and nocturnal Dipsadidae snakes, it was suggested that the type of specialization may also be associated with daily activity patterns and foraging strategies (Hauzman et al., 2018), in which diurnal species that actively forage during the day display visual streaks, while nocturnal species or those that feed on slow-moving prey have area centralis in different regions of the retinas (Hauzman et al., 2018). These studies revealed the variability of adaptations of the visual structures of snakes, and indicate that different selective forces may shape the retinal architecture irrespective of phylogenetic blueprints.
The Viperidae family represents a valuable model for investigating adaptations of the visual structures due to the diversity of species and habitats occupied, predation strategies with accurate strike performances (Reiserer, 2002;Chen et al., 2017;Schraft and Clark, 2019), and an elaborate thermosensitive sensory system in pitviper species (subfamily Crotalinae) integrated with inputs from visual neurons in the tectum (Hartline et al., 1978). Additionally, some viperid species have ontogenetic changes in the niche occupied and thus, represent a unique opportunity to explore how morphological adaptations of the retina might be associated with their visual ecology. Viperids are primarily nocturnal or crepuscular, and their retinas have a predominance of rods, highly sensitive photoreceptors adapted to low light (scotopic) conditions, and three types of cones, photoreceptors responsible for daylight (photopic) vision. The rods contain the typical rhodopsin (RH1) photopigment, and cones are classified as single cones and double cones sensitive to medium/long wavelengths, with the LWS photopigment, and single cones sensitive to short wavelengths, with the SWS1 photopigment (Bittencourt et al., 2019;Gower et al., 2019).
In this study, we compared the density and distribution of photoreceptors and ganglion cell layer (GCL) cells in the retinas of two pitvipers. The rattlesnake, Crotalus durissus (Figure 1), is a terrestrial snake that inhabits open fields of the Cerrado, a Brazilian savannah-like habitat, and actively search for places to hide before using the sit-and-wait hunting strategy to ambush rodents (Salomão et al., 1995;Sawaya et al., 2008;Tozetti andMartins, 2008, 2013). The lancehead, Bothrops jararaca (Figure 1), inhabits predominantly forested areas of the Atlantic Rain Forest and has ontogenetic changes in niche occupied and in behavior. Adults are terrestrial and use similar hunting strategies as rattlesnakes to prey on mammals. Juveniles are semi-arboreal and use the sitand-wait strategy and caudal luring to attract ectothermic vertebrates, mainly anurans (Sazima, 1991(Sazima, , 1992(Sazima, , 2006Marques and Sazima, 2004). We hypothesized that the differences in behavior and niche occupied by juveniles and adults of B. jararaca might be associated with rearrangements of the retinal architecture according to specific visual needs. With a stereological approach to quantify neurons in whole-mounted retinas, we observed differences in the density and distribution of cells between species, especially regarding the proportion and distribution of rods and cones, and differences in the retinal topography of juveniles and adults of B. jararaca that might reflect ontogenetic changes in the visual ecology. This study demonstrates that the habitat occupied by snakes and their foraging strategies are associated with different patterns of neuron distribution in the retina, which highlights the importance of their retinal specializations for visually guided behaviors.

Snakes
Snakes (n = 33) were provided by the Butantan Institute, São Paulo, Brazil, and were euthanized with a lethal injection of sodium thiopental (100 mg/kg). Animal procedures were in accordance with ethical principles of animal management and experimentation established by the Brazilian Animal Experiment College (COBEA). The project was approved by the Ethics Committee of Animal Research of the Psychology Institute, University of São Paulo, Brazil. Individuals were classified as adults or juveniles based on the snout-vent length (SVL) and body mass (Supplementary Table 1; Sazima and Abe, 1991;Almeida-Santos, 2005;de Moraes, 2008;Barros, 2011;Fiorillo et al., 2020a).

Tissue Processing
Following euthanasia, the eyes were enucleated, their axial lengths were measured, and a small radial incision was made in the dorsal region for later orientation. The corneas were removed and the eyecups were fixed in 4% paraformaldehyde (PFA) diluted in phosphate buffer (PB) 0.1 M, for 3 h. The retinas were carefully dissected and maintained in PB 0.1 M at 4 • C. When needed, the retinas were bleached in 10% hydrogen peroxide diluted in PB 0.1 M, for 24-48 h, at room temperature, prior to immunohistochemistry or Nissl procedures.

Antibody Characterization and Specificity
Immunohistochemistry procedures were performed with polyclonal antibodies raised in rabbits against the last 42 amino acids of the C-terminal of human blue opsin (Sigma-Aldrich; AB5407), against the last 38 amino acids of the C-terminal of human red/green opsins (Sigma-Aldrich; AB5405), or raised in goats against a synthetic peptide with 20 amino acids of human blue opsin (Santa Cruz Biotechnology; sc-14363) ( Table 1). The specificity of the antibodies for snakes was described previously (Hauzman et al., , 2017Bittencourt et al., 2019;Gower et al., 2019). Double labeling with the antibodies against SWS1 and against LWS opsins showed differential labeling of distinct photoreceptor populations (Figure 2), further indicating the specificity of both antibodies for particular types of cones. We also assessed the specificity of the two anti-SWS1 antibodies, by incubating 12 µm retinal sections of B. jararaca and C. durissus obtained at −25 • C with a cryostat (Leica, CM1100; Nussloch, Germany), with a mixture of both antibodies, rabbit anti-SWS1 (AB5407; 1:200) and goat anti-SWS1 (sc14363; 1:200). Immunofluorescence visualization showed a small number of small single cones labeled by both antibodies (data not shown), indicating the specificity of both antibodies against SWS1 cones.

Nissl Staining
We used Nissl staining technique in whole-mounted retinas to analyze the population of GCL cells. Small radial incisions were made in the dissected retinas to allow them to be flat-mounted onto gelatinized glass slides, with the GCL facing up. The retinas were incubated with paraformaldehyde vapors overnight, at room temperature, for enhancing the adherence to the slide and to increase the differentiation of ganglion cells during staining (Coimbra et al., 2006). The tissues were rehydrated in decreasing ethanol series (95, 70, 50%), rinsed in distilled water acidified with glacial acetic acid, stained in aqueous solution of 2% cresyl violet for approximately 1 min, dehydrated in increasing concentrations of ethanol, cleared in xylene, and coverslipped with DPX (Sigma-Aldrich. St. Louis, MO, EUA). To analyze the density and distribution of GCL cells, we used the cytological criteria proposed by Wong (1989) to distinguish ganglion cells from amacrine and glial cells. The ganglion cells  were differentiated by their larger polygonal soma, abundant Nissl substance in the cytoplasm, and a prominent nucleolus (Figure 3; Hart et al., 2012). Glial cells were differentiated by smaller soma size and a round and darkly stained profile, and amacrine cells were identified by a smaller and circular profile with a more darkly stained nucleus compared to ganglion cells . The glial and amacrine cells were not counted.

Stereological Assessment of the Density and Distribution of Retinal Neurons
The density and distribution of retinal neurons were accessed using a stereological approach based on the optical fractionator method (West et al., 1991), modified for retinal whole-mounts (Coimbra et al., 2009, using a motorized microscope (DM5500B, Leica Microsystems, Germany), connected to a computer running the Stereo Investigator software (MicroBrightField, Colchester, VT). The coordinates of the outer edges of the retinas were obtained using a 5x/NA 0.15 objective. Approximately 200 counting frames were positioned in a random fashion covering the entire area of the retina. Cells were counted when laying entirely within the counting frame or when intersected the acceptance lines, without touching the rejection lines (Gundersen, 1997). The area of the counting frame and the sampling grids varied according to the cell types and the retinal area, and were defined in pilot experiments, in order to obtain an acceptable Scheaffer coefficient of error (CE) (<0.10) (Scheaffer et al., 1996). The stereological parameters used to estimate the number of photoreceptors and GCL cells of each retina are described in Supplementary Tables 2, 3. To estimate the total population of neurons (N total ), we considered the area of sampling fraction (asf) according to the algorithm: N total = Q × 1/asf , where Q is the sum of the total number of neurons counted and the area of sampling fraction is the ratio between the counting frame and the sampling grid (Coimbra et al., 2009).

Anatomical Estimates of the Visual Acuity
We estimated the upper limits of spatial resolving power based on the peak density of GCL cells and the presumed focal distance of the eyes. The focal lengths were estimated by freezing and sectioning eyes from one adult and one juvenile individual of each species. The fresh eyes were enucleated and rapidly frozen, embedded in Tissue-Tek OCT compound (Sakura Finetechnical Co., Tokyo, Japan), and the blocks were sectioned at −25 • C on a cryostat (Leica CM1100; Nussloch, Germany). Photographs of the blocks were taken every 12 µm, with a camera (Axio CamMR, Carl ZeissVision, Germany) coupled to a stereomicroscope (SMZ775-T, NIKON, Japan), and a computed running the Axio Vision 4.1 software (Carl Zeiss, Germany). The areas of the lens and eyes were measured using the ImageJ software (NIH, Bethesda, United States) and the photographs that had the larger lens diameters were identified and used for optical measurements. Measurements of the axial length of the eye, lens axial diameter, and posterior nodal distance (PND) that represents the focal length and corresponds to the distance from the center of the lens to the retina-choroid border were taken along the optical axis, which was located by connecting the geometric centers of the optical components (Lisney and Collin, 2008). In B. jararaca, the focal length corresponded to 50% of the adult eye axial length of the adult ( Figure 4A) and 67% of the juvenile. In C. durissus the focal distance corresponded to 60% of eyes axial length of both adult ( Figure 4B) and juvenile. Those values were used to estimate the focal distance of the eyes used for GCL cell counts.
We estimated the upper limits of spatial resolution considering two possible arrangements of the ganglion cells, a hexagonal and a square array. The distance d subtended as 1 degree on the retina was determined from the PND and the equation d = (2πPND)/360 (Pettigrew et al., 1998). Considering a hexagonal array, we estimated the average spacing between cells (S) using the formula S 2 = 2/(D √ 3), where D is the peak density of GCL cells in cells/mm 2 . The maximum spatial   frequency (v) (Nyquist) of a sinusoidal grating (Snyder and Miller, 1977) that has resolution with this cellular arrangement was calculated as v = 1/S √ 3. This value was multiplied by the distance d, to obtain the spatial resolution in cycles per degree (cpd). In a second approach, considering that the ganglion cells might be organized in a square distribution, we estimated the linear density of GCL cells (cells/mm 2 ) from the square root of the peak density (D), and divided the linear cell density by 2 (because at least 2 cells are required to detect 1 cycle of a given spatial frequency). The resulting value was multiplied by the distance d to obtain the visual acuity in cpd (Pettigrew et al., 1998;Coimbra et al., 2013).
To infer the behavioral significance of the estimated visual acuity in an ecological context, we predicted the minimum size of objects that the viperid snakes can spatially resolve (Coimbra et al., 2017). To do so, we estimated the angular distance in the retina that corresponds to one cycle by calculating the inverse of the spatial resolving power (cycles/degree). This value was divided by 2 to obtain the minimum angle of resolution (MAR), which represents the angular distance of the smallest resolvable detail on the retina. Subsequently, using the trigonometric relationship between the MAR and a presumed distance (D), relevant, for instance, for foraging a prey or for predator detection, we estimated the minimum object size (obj), according to the equation: D = obj/tanMar. According to the Nyquist sampling theorem, because an object needs to be twice the threshold to be spatially resolved, we multiplied the minimum object size by 2 (Marshall, 2000).

Statistical Analysis
Statistical analyzes were performed with R (version 4.0.2) 1 and the RStudio software (1.3.959), to compare the total population and mean densities of retinal neurons among the four sampled groups, adults and juveniles of B. jararaca and of C. durissus, as well as the retinal area, eyes axial length and visual acuity. The normality of the distribution of values in each group was checked using the Shapiro-Wilk test, and homogeneities among groups were analyzed using the Levene test. The non-parametric test of Mann-Whitney for independent samples was applied for comparisons, even when values had normal distribution, due to the low sampling size in each group. Differences were considered significant when p < 0.05.

RESULTS
We analyzed the density and distribution of photoreceptors and GCL cells in 18 retinas of B. jararaca (adults: n = 8; juveniles: n = 10) and 17 retinas of C. durissus (adults: n = 9; juveniles: n = 8) (Tables 2-4). The different populations of cones were identified by immunohistochemistry labeling with antibodies against SWS1 and LWS opsins (Figure 2A). The total photoreceptor population was viewed under bright light and by adjusting the focus of the microscope into the photoreceptor's inner segments level ( Figure 2B). Cones and rods were differentiated by the larger diameters of the inner segments of cones compared to the small and highly packed inner segments of rods (Gower et al., 2019; Figure 2B). For cell counting in each sampling field, the photoreceptors labeled by the antibodies were counted first, then all photoreceptors (cones and rods) were counted under bright light to analyze the proportion of each photoreceptor type (rods, LWS, and SWS1 cones). Single and double LWS cones were not always easily differentiated from each other and were quantified together (Figure 2A).

Photoreceptors Topography
The distribution of photoreceptors showed differences between and within species. In adults of B. jararaca, rods and cones were distributed in poorly defined horizontal streaks (Figure 7 and Supplementary Figures 1, 2). In juveniles, rods were concentrated in an anisotropic area centralis in the dorsal retina, while cones were concentrated in the ventral retina (Figure 7 and Supplementary Figures 1, 2). Mean density values were estimated from retinal sectors (dorsal, ventral, temporal, and nasal). In adults and juveniles of B. jararaca, higher densities of rods were located in the temporal retina. Higher density of cones were located in the ventral retina of juveniles and in the temporal region in adults (Supplementary Figure 3 and Supplementary Table 4). In C. durissus, the isodensity maps of rods and cones were similar between adults and juveniles. Rods were concentrated in the dorsal retina in an anisotropic area centralis, and cones were organized in poorly defined visual streaks (Figure 7 and Supplementary Figures 1, 2). Density estimates of retinal sectors indicated higher density of rods in the temporal retina of juveniles and in the dorsal retina of adults.  Table 4). In both species, the SWS1 cones showed diffuse distributions, with peak densities located in the central or ventral retina (Figure 8 and Supplementary Figure 4). The distribution of LWS cones in the retinas of adults and juveniles of both species was similar to that observed for the distribution of total cones, with poorly defined visual streaks in adults of both species and in juveniles of C. durissus, and a higher concentration of LWS cones in the central and ventral retina of juveniles of B. jararaca (Figure 8 and Supplementary Figure 5).

Ganglion Cell Layer Cells and Estimates of the Spatial Resolving Power
The population of GCL cells was estimated from 11 retinas of B. jararaca (adults: n = 4, juveniles: n = 7) and 9 retinas of C. durissus (adults: n = 4, juveniles: n = 5) ( Table 4). The total population of GCL cells was similar between adults and juveniles of both species, ranging from 169,569 ± 13,730 (median: 163,733) cells in juveniles of B. jararaca and 215,141 ± 22,291 (median: 221,798) cells in adults of C. durissus (Figure 9 and Table 4). The average density of GCL cells was significantly higher in retinas of juveniles (9,871 ± 900 cells/mm 2 ; median: 9,496 cells/mm 2 in B. jararaca and 9,774 ± 1,078 cells/mm 2 ; median: 9,951 cells/mm 2 in C. durissus) compared with adults (4,880 ± 1,092 cells/mm 2 ; median: 4,657 cells/mm 2 in B. jararaca and 4,429 ± 448 cells/mm 2 ; median: 4,517 cells/mm 2 in C. durissus) (Figure 9 and Table 4). In the four groups analyzed, the isodensity maps of GCL cells did not show defined distribution patterns (Figure 10 and Supplementary  Figure 6). In retinas of adults and juveniles of C. durissus and adults of B. jararaca, we observed diffuse distributions and peak density of cells in the temporal retina (Figure 10 and  Figure 6). In juveniles of B. jararaca the GCL cells were concentrated in the ventral retina with a decreasing ventral-dorsal gradient and peak density located in the ventral area (Figure 10 and Supplementary Figure 6). Density estimates in retinal sectors showed higher values in the temporal retina of adults and juveniles of C. durissus and adults of B. jararaca (Supplementary Figure 3 and Supplementary Table 5). In juveniles of B. jararaca, higher densities were located in the ventral retinas (Supplementary Figure 3 and Supplementary  Table 5), in agreement with the observed isodensity maps (Figure 10 and Supplementary Figure 6).
The theoretical upper limits of spatial resolving power were estimated based on the ganglion cell peak densities using two approaches, one considering that the ganglion cells are organized in a hexagonal array and the other considering a square array.
From the estimated visual acuity values we calculated the minimum angle of resolution. The estimated distance relevant for predatory behavior was calculated considering the tangent of the minimum angle of resolution and minimum target size. For species with the lowest (2.2 cpd) and highest (2.8 cpd) spatial resolving power, we estimated a minimum angle of resolution of 0.45 • and 0.36 • . In an ecological context, snakes with lower estimated spatial resolution can observe an object with a minimum size of 10 cm, such as a small mammal, at a distance of approximately 13 m. Snakes with the higher presumed spatial resolution might be able to observe an object of the same size at a distance of about 17 m. These estimates indicate that at these predicted distances objects larger than the minimum target size can be spatially detected.

DISCUSSION
In this study, we analyzed the density and distribution of neurons in whole-mounted retinas of Viperidae snakes, considering an ontogenetic approach. Our analyzes revealed a predominance of rods in the outer retinas of Bothrops jararaca and Crotalus durissus, and three distinct populations of cones: single and double cones containing the LWS photopigment and single cones containing the SWS1 photopigment, as previously described for viperids (Walls, 1942;Underwood, 1967b;Bittencourt et al., 2019;Gower et al., 2019). Our results showed that the two species have similar density values of photoreceptors and GCL cells. However, the distribution of these neurons differed between species and between juveniles and adults of B. jararaca, pointing to a reorganization of the retinal architecture that might be associated with the ontogenetic changes in the niche occupied and hunting strategies, as summarized in Figure 11.

Photoreceptor Population
The retinas of B. jararaca and C. durissus had a high proportion of rods (about 90% of the photoreceptors) ( Figure 5 and Table 2), which indicates high sensitivity to light, in agreement with their nocturnal or crepuscular activity pattern (Martins et al., 2001;Fiorillo et al., 2020b). Previous studies described the predominance of rods in retinas of viperid snakes based on analysis of retinal sections (Walls, 1934(Walls, , 1942Underwood, 1967b;Bittencourt et al., 2019;Gower et al., 2019) and fragments of flat-mounted retinas (Gower et al., 2019). Snakes stand out among vertebrates by their highly variable patterns of photoreceptor morphology (Walls, 1942). In the Caenophidia group ("advanced" snakes), nocturnal species from different families, including viperids, have four types of photoreceptors, with a unique type of double cone that differs from double cones of other vertebrates, with a large principal member and an extremely slender and attached accessory member (Walls, 1942). On the other hand, diurnal caenophidian snakes have "purecone" retinas, with the absence of typical rods, the presence of a transmuted cone-like rod, and lower photoreceptor density compared to nocturnal species (Walls, 1942;Schott et al., 2016;Hauzman et al., 2017;Hauzman, 2020).
The comparison of the populations of cones and rods in B. jararaca and C. durissus showed a higher proportion of cones in retinas of C. durissus (10% in juveniles and adults) compared to B. jararaca (7.5% in juveniles and 6.5% in adults) ( Figure 5 and Table 2). This difference, although subtle, might reflect functional differences in temporal and spatial resolution between both species, which agree with the habitat occupied. C. durissus inhabit predominantly open areas of the Cerrado, with higher incidence of light (Miranda et al., 1997) compared to the closed forested areas predominantly occupied by B. jararaca (McWilliam et al., 1993;de Paula and Lemos Filho, 2001). In both species, the population of cones is dominated by large single cones and double cones that contain the LWS photopigment (about 80-90%) (Figure 6 and Table 3), with spectral sensitivity peak (λ max ) predicted at 555 nm (Bittencourt et al., 2019). Small single cones with the SWS1 photopigment comprise about 10% of the cone population (Figure 6 and Table 3), with λ max predicted at the UV range (360-370 nm) in both species (Bittencourt et al., 2019). Gower et al. (2019) identified double cones with the SWS1 photopigment in retinas of two viperid snakes, Echis coloratus and C. durissus, a unique type of photoreceptor described for the first time in vertebrates. However, in our analysis of wholemounted retinas, this type of cone was not identified, and all double cones observed were found to contain only the LWS photopigment (Figure 2).
The distribution of photoreceptors differed between C. durissus and B. jararaca and between juveniles and adults of B. jararaca (Figure 7). In C. durissus, the visual streak formed by cones might reflect a better panoramic view of the environment under photopic conditions (Figure 7). This specialization might benefit scanning of the environment while searching for shelters for body temperature control during the day, a frequent behavior of snakes that occupy open environments where they are subject to overheating (Tozetti and Martins, 2008). On the other hand, the higher density of rods in the dorsal retina (Figure 7) might improve light sensitivity in the lower visual field, possibly favoring foraging behavior and searching for rodents during twilight and at night. These patterns of distribution of cones and rods were found in both, juveniles and adults of C. durissus (Figures 7, 11).
In B. jararaca, ontogenetic changes in the niche occupied seem to be associated with plasticity of the retinal architecture. In juveniles, a higher density of cones was observed in the ventral retina (Figures 7, 11), a specialization that might provide higher acuity in the upper visual field. Oppositely, a higher density of rods in the dorsal retina might benefit the view of the lower visual field under scotopic conditions (Figure 7). Juveniles of B. jararaca occupy the arboreal stratum and use sit-and-wait and caudal luring as hunting strategies. Therefore, it is plausible to speculate that the difference in the distribution of cones and rods might be associated with the direction of predation pressure. Cones located in a ventral area might favor the view of aerial predators, such as diurnal birds approaching from above during the day (Figure 11) (Sazima, 1992;Costa et al., 2014). Higher concentration of rods in the dorsal retina might benefit the view of terrestrial predators, such as marsupials, approaching from below, during the night (Figure 11) (Emmons, 1990;Sazima, 1992;Jared et al., 1998;Oliveira and Santori, 1999). Compared to terrestrial snakes, arboreal species display a higher number of defensive tactics, which is likely associated with greater exposure to predators approaching from a variety of directions .
Adults of B. jararaca occupy terrestrial environments, lose the caudal luring, and actively forage for endothermic prey (Sazima, 1992(Sazima, , 2006. Their retinas have cones and rods distributed in poorly defined visual streaks, a specialization that benefits the panoramic view of the terrestrial forested stratum under photopic and scotopic conditions (Figures 7, 11). A similar distribution of cones was described in terrestrial, arboreal and semiaquatic colubrids (Hauzman, 2014;Hauzman et al., 2014), indicating that this specialization is widely observed in snakes and might contribute to active foraging behavior in different environments. As far as we are aware, this is the first description of a horizontal streak formed by rods in the retinas of snakes.
The distribution of LWS cones (single and double) was similar to the distribution of total cones in both species, as expected based on their high proportion (Figure 8). The distribution of SWS1 cones did not show a defined pattern of specialization. Higher densities of SWS1 cones were found in the ventral and central retina in adults and juveniles of both species (Figure 8), as described for diurnal colubrid snakes . This distribution might favor the view of potential aerial predators approaching from the upper visual field, as many bird species are important predators of snakes (Martins and Oliveira, 1998;Tozetti, 2006;Specht et al., 2008;Costa et al., 2014). Higher densities of UV cones were also described in the ventral retina of mammals (Szél et al., 1992(Szél et al., , 1994Famiglietti and Sharpe, 1995;Peichl et al., 2005;Huber et al., 2010). In mice, UV cones contribute to chromatic discrimination of the upper visual field (Szatko et al., 2020), and might be relevant for visualizing the silhouette of aerial predators against the blue sky background (Calderone and Jacobs, 1995;Szél et al., 1992).

Density and Topography of Ganglion Cell Layer Cells and Estimates of the Visual Acuity
The mean density of GCL cells was similar between B. jararaca and C. durissus. However, juveniles had higher average density values compared to adults (Figure 9 and Table 4). This difference can be attributed to the increase in the area of the retina. Juveniles have smaller retinas and thus, higher cell packaging. As the animals grow, the eye increases, such as the retinal area, but the cell population remains constant, resulting in lower density values. A higher GCL cell density in juveniles compared to adults was also described in fish (Hagedorn and Fernald, 1992;Bailes et al., 2006), birds (Straznicky and Chehade, 1987), amphibians (Nguyen and Straznicky, 1989), and mammals (Robinson et al., 1989). Considering there is no neuron loss or additional generation of cells in the GCL, we suggest the increase of retinal area followed by the decrease of cell density is associated with interstitial growth, as described for chicken retinas (Straznicky and Chehade, 1987). It is notable that in Bothrops jararaca these changes in cell density are associated with rearrangements of the cell distribution.
The distribution of GCL cells was similar between adults and juveniles of C. durissus, with no defined type of specialization, and with peak density of cells in the temporal retina (Figures 10, 11). This specialization might benefit the view of the frontal field and favor strike performance of endothermic prey. The same topographic pattern of GCL cells was found in retinas of adults of B. jararaca (Figures 10, 11). The visual information arriving in the midbrain, from the projections of this temporal specialization might be combined with infrared information from the loreal pit, a thermosensory organ of Crotalinae snakes located between the eyes and the nostrils (Noble and Schmidt, 1973). The nerve FIGURE 11 | Schematic diagram illustrating the niche occupied by adults and juveniles of Bothrops jararaca and Crotalus durissus, and the distribution of neurons in the retinas of each species. (Upper panel) B. jararaca inhabits forested areas and has ontogenetic changes in diet, behavior, and microhabitat occupied. Juveniles are arboreal, feed on ectothermic prey (attracted by caudal luring), and are susceptible to predation pressure coming from different directions. A higher density of cones and GCL cells in the ventral retina might benefit the view of aerial predators under photopic conditions. A higher density of rods in the ventral retina might favor the view of terrestrial predators in the lower visual field under scotopic conditions. Adults are terrestrial and actively hunt for endothermic prey. Cones and rods are distributed in poorly defined visual streaks, and GCL cells have higher densities in the temporal retina. These specializations might benefit the panoramic view of the environment under photopic and scotopic conditions and of strike performances. (Lower panel) Adults and juveniles of C. durissus occupy open environments (Cerrado) and feed on endothermic prey. In both, cones form poorly defined visual streaks that might benefit scanning the environment to search for shelters for body temperature control during the day to avoid overheating, a usual behavior in snakes that occupy open areas. The higher density of rods in the dorsal retina might favor foraging behavior at scotopic conditions. The peak density of GCL cells in the temporal retina might benefit the view of the frontal field and favor strike performances. N, nasal; V, ventral.
terminals of the loreal pit cells project to the optic tectum in the midbrain (Newman and Hartline, 1982), where the visual and thermosensory inputs are combined, allowing an integrated perception from ultra-violet to infrared wavelengths (Newman and Hartline, 1981;Moiseenkova et al., 2003;Goris, 2011), important in both scotopic and photopic conditions. During the night, thermal detection enables the perception of temperature changes. During the day, discrimination of infrared information may be compromised due to higher temperatures (Chen et al., 2012), yielding predominance of perceptual information to visual inputs.
In juveniles of B. jararaca the higher density of GCL cells in the ventral retina is in agreement with the distribution of cones (Figures 8, 10, 11), which indicates a higher convergence from cones to GCs in this region (from approximately 4:1 in the periphery, to ∼1:1 in the area centralis), enabling higher spatial resolution in the upper visual field. A similar distribution pattern of GCL cells, with higher density in the ventral retina was observed in the dipsadid snakes Dipsas albifrons and Sibynomorphus neuwiedi (Hauzman et al., 2018). Both species are nocturnal, semi-arboreal, and feed on ectothermic prey (gooeaters) (Maia et al., 2011;Sazima and Muscat, 2016), such as juveniles of B. jararaca. Thus, we suggest that this type of specialization might be associated with the niche occupied by snakes, the direction of predatory threats, and particularities of hunting strategies, which do not involve active foraging of fast moving prey (Figure 11). Reorganization of the distribution of GCL cells according to ontogenetic changes in ecology and behavior was also described in frogs (Dunlop and Beazley, 1981) and in fish (Shand, 1997;Bozzano and Catalan, 2002;Miyazaki et al., 2011). In the fish Acanthopagrus butcheri, the retinas of juveniles have a temporal area centralis, which favors object detection in the frontal visual field and might be important for feeding on plankton in pelagic waters. On the other hand, adults occupy the benthic environment. This ecological switch is accompanied by a displacement of the area toward the dorsal or dorso-temporal retina, which might increase acuity in the lowerfrontal visual field and benefit their hunting strategies for the capture of small fish, polychaete worms and detritus obtained from the benthos (Shand et al., 2000).
The estimated upper limits of spatial resolving power based on the peak cell density of GCL cells of the viperid snakes varied from 2 to 3 cpd ( Table 4). These values were similar to those estimated for diurnal dipsadid snakes of the genus Philodryas spp.  and for the marine elapid Aipysurus laevis , with spatial resolution between 2.3 and 2.6 cpd, and were higher than the values estimated for nocturnal dipsadids, with about 1.3 cpd (Hauzman et al., 2018). Our comparisons showed slightly higher estimated visual acuity in adults of C. durissus compared to the other groups analyzed (Figure 9). We predicted that the estimated visual acuity of ∼2.3 cpd of adults and juveniles of B. jararaca and juveniles of C. durissus would allow the view of potential prey such as rodents, at a minimum size of 10 cm, at a distance of about 13 m. On the other hand, adults of C. durissus with estimated visual acuity at ∼2.8 cpd, might be able to detect the same minimum object size at a higher distance, about 17 m. These estimates also indicate that these snakes can have a high capability to observe potential predators such as birds of prey (Martins and Oliveira, 1998;Tozetti, 2006;Specht et al., 2008;Costa et al., 2014) at longer distances. However, behavioral analyzes are necessary to confirm if the estimated values of spatial resolution based on anatomical data are reached, and to ascertain whether the differences imply functional significance, relevant to the species ecology and visual behavior.

CONCLUSION
Interspecific differences in the density and distribution of retinal neurons were identified in the viperid snakes B. jararaca and C. durissus. Rod-dominated retinas represent a functional advantage associated with the nocturnal activity pattern of both species. The higher proportion of cones in C. durissus might represent an adaptation to the open and brighter Cerrado environment. This study showed for the first time in snakes that ontogenetic changes in ecology and behavior are associated with morphological plasticity of the retinas. We suggest that the differences in the niche occupied, hunting strategies, prey type, and direction of approach of predatory threats are reflected in reorganizations of the distribution of retinal neurons between juveniles and adults of B. jararaca. In comparison, the absence of variation in neuron distribution patterns in retinas of C. durissus throughout life agrees with the maintenance of the same ecological and behavioral traits in juveniles and adults. These results highlight the importance of retinal specializations for the performance of visually guided behaviors, and how habitat use and hunting strategies might represent relevant ecological forces that shape the retinal architecture in snakes. Future behavioral studies should be applied to verify functional implications of the retinal organizations as well as the estimated upper limits of spatial resolving power, between 2 and 3 cpd. In addition, the patterns of neuronal connections in the retina rely on a highly intricate network for visual processing. The results of this study open up a new avenue for future analysis on the connectivity patterns of the inner retina of Viperidae snakes.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
The animal study was reviewed and approved by the Ethics Committee of Animal Research of the Psychology Institute, University of São Paulo, Brazil.
Hamasaki for providing reagents. We are grateful to Carola A. M. Yovanovich and Marcio R. C. Martins for valuable discussions, and Fernando A. F. Rocha and Marcelo F. Costa for methodological contributions. We also acknowledge the two reviewers for their valuable suggestions.