Visual mate choice in poison frogs

Anti-predator defences are linked with high levels of genetic differentiation in frogs

Predator-prey interactions have been suggested as drivers of diversity in different lineages, and the presence of anti-predator defences in some clades is linked to higher rates of diversification. Warning signals are some of the most widespread defences in the animal world, and there is evidence of higher diversification rates in aposematic lineages. The mechanisms behind such species richness, however, are still unclear. Here, we test whether lineages that use aposematism as anti-predator defence exhibit higher levels of genetic differentiation between populations, leading to increased opportunities for divergence. We collated from the literature more than 3000 pairwise genetic differentiation values across more than 700 populations from over 60 amphibian species. We find evidence that over short geographical distances, populations of species of aposematic lineages exhibit greater genetic divergence relative to species that are not aposematic. Our results support a scenario where the use of warning signals could restrict gene flow, and suggest that anti-predator defences could impact divergence between populations and potentially have effects at a macro-evolutionary scale.

Selection on Visual Opsin Genes in Diurnal Neotropical Frogs and Loss of the SWS2 Opsin in Poison Frogs

Amphibians are ideal for studying visual system evolution because their biphasic (aquatic and terrestrial) life history and ecological diversity expose them to a broad range of visual conditions. Here, we evaluate signatures of selection on visual opsin genes across Neotropical anurans and focus on three diurnal clades that are well-known for the concurrence of conspicuous colors and chemical defense (i.e., aposematism): poison frogs (Dendrobatidae), Harlequin toads (Bufonidae: Atelopus), and pumpkin toadlets (Brachycephalidae: Brachycephalus). We found evidence of positive selection on 44 amino acid sites in LWS, SWS1, SWS2, and RH1 opsin genes, of which one in LWS and two in RH1 have been previously identified as spectral tuning sites in other vertebrates. Given that anurans have mostly nocturnal habits, the patterns of selection revealed new sites that might be important in spectral tuning for frogs, potentially for adaptation to diurnal habits and for color-based intraspecific communication. Furthermore, we provide evidence that SWS2, normally expressed in rod cells in frogs and some salamanders, has likely been lost in the ancestor of Dendrobatidae, suggesting that under low-light levels, dendrobatids have inferior wavelength discrimination compared to other frogs. This loss might follow the origin of diurnal activity in dendrobatids and could have implications for their behavior. Our analyses show that assessments of opsin diversification in across taxa could expand our understanding of the role of sensory system evolution in ecological adaptation.

Sexual signal evolution and patterns of assortative mating across an intraspecific contact zone

Contact zones provide important insights into the evolutionary processes that underlie lineage divergence and speciation. Here, we use a contact zone to ascertain speciation potential in the red-eyed treefrog (Agalychnis callidryas), a brightly coloured and polymorphic frog that exhibits unusually high levels of intraspecific variation. Populations of A. callidryas differ in a number of traits, several of which are known sexual signals that mediate premating reproductive isolation in allopatric populations. Along the Caribbean coast of Costa Rica, a ~100 km contact zone, situated between two phenotypically and genetically divergent parent populations, contains multiple colour pattern phenotypes and late-generation hybrids. This contact zone provides the opportunity to examine processes that are important in the earliest stages of lineage divergence. We performed analyses of colour pattern variation in five contact zone sites and six parental sites and found complex, continuous colour variation along the contact zone. We found discordance between the geographic distribution of colour pattern and previously described genomic population structure. We then used a parental site and contact zone site to measure assortative mating and directional selection from naturally-occurring amplectant mating pairs. We found assortative mating in a parental population, but no assortative mating in the contact zone. Furthermore, we uncovered evidence of directional preference towards the adjacent parental phenotype in the contact zone population, but no directional preference in the parent population. Combined, these data provide insights into potential dynamics at the contact zone borders and indicate that incipient speciation between parent populations will be slowed.

Outcomes of multifarious selection on the evolution of visual signals

Multifarious sources of selection shape visual signals and can produce phenotypic divergence. Theory predicts that variance in warning signals should be minimal due to purifying selection, yet polymorphism is abundant. While in some instances divergent signals can evolve into discrete morphs, continuously variable phenotypes are also encountered in natural populations. Notwithstanding, we currently have an incomplete understanding of how combinations of selection shape fitness landscapes, particularly those which produce polymorphism. We modelled how combinations of natural and sexual selection act on aposematic traits within a single population to gain insights into what combinations of selection favours the evolution and maintenance of phenotypic variation. With a rich foundation of studies on selection and phenotypic divergence, we reference the poison frog genus Oophaga to model signal evolution. Multifarious selection on aposematic traits created the topology of our model's fitness landscape by approximating different scenarios found in natural populations. Combined, the model produced all types of phenotypic variation found in frog populations, namely monomorphism, continuous variation and discrete polymorphism. Our results afford advances into how multifarious selection shapes phenotypic divergence, which, along with additional modelling enhancements, will allow us to further our understanding of visual signal evolution.

Mate choice in a promiscuous poison frog

In many animal species, members of one sex, most often females, exhibit a strong preference for mating partners with particular traits or resources. However, when females sequentially mate with multiple partners, strategies underlying female choice are not very well understood. Particularly, little is known if under such sequential polyandry females mate truly randomly, or if they actively try to spread mating events across multiple partners. In the present study, we used the highly promiscuous poison frog Allobates femoralis to investigate whether promiscuity could result from a preference for novel mates. Furthermore, we examined the importance of call characteristics for mate choice. We conducted mate choice experiments in a laboratory setup, by presenting females with recent mating partners or novel males. We recorded call characteristics of both males and the time females spent close to each male. In our trials, females preferred previous mating partners over novel males and also males with shorter advertisement calls. Results from previous studies on A. femoralis suggest that females in our trials recognized previous partners based on individual call characteristics. While mating decisions in the wild and in the laboratory might differ, our study provides first evidence for female mate choice in a poison frog with sequential polyandry.

Joint effects of female preference intensity and frequency-dependent predation on the polymorphism maintenance in aposematic sexual traits

Maintenance of variation in aposematic traits within and among populations is paradoxical because aposematic species are normally under positive frequency‐dependent predation (PFD), which is expected to erode variation. Aposematic traits can evolve in an ecological context where aposematic traits are simultaneously under mate choice. Here, we examine how the mate preference intensity affects the permissiveness of polymorphism in sexually selected aposematic traits under different PFD regimes. We use the haploid version of the classical sexual selection model and show that strong mate preferences can substantially increase the permissiveness of polymorphism in aposematic traits under different PFD regimes. The Fisher process can interact with PFD, and their interaction can promote the maintenance of polymorphism within populations when mate preferences are strong. We show that the same selective conditions that promote the maintenance of polymorphism within populations reduce the likelihood of divergence in aposematic traits among populations.

Reproductive behaviors promote ecological and phenotypic sexual differentiation in the critically endangered Lehmann’s poison frog

Territoriality and parental care are complex reproductive behaviors found in many taxa from insects to mammals. Parental care can be carried out by the female, the male, or both, depending on the species. Territoriality, in contrast, is predominantly displayed by males. Different selective pressures imposed on individuals from the sex performing territorial or parental care behaviors may also lead to sexual differentiation in other life-history traits. Due to their territorial behavior and their diversity of parental care behaviors, Neotropical poison frogs are an excellent study system to investigate whether behavioral traits can influence sexual differentiation in intrinsic or extrinsic traits of individuals. Here, we evaluate whether territorial and parental care behaviors mediate sexual differentiation in ecological (habitat use) and phenotypic (coloration, morphology) traits in the critically endangered Lehmann’s poison frog (Oophaga lehmanni), a species in which males defend territories while females provide parental care. We found sex differences in habitat use and morphological traits, but not in coloration. Males use trunks and green leaves as perches more frequently and are found on higher substrates, than females. We found no sex differences in body size, but females have longer arms than males, which is probably associated with their parental duties (climbing trees to feed the tadpoles). Altogether, our results provide evidence that selection pressures act differently on male and female traits, and that territoriality and parental care may promote the evolution of sexual differentiation in dendrobatids. Long-term wildlife observations are essential to identify important life-history traits and to evaluate hypotheses about the behavioral ecology and conservation of this and other vertebrate species.

Poison frog social behaviour under global change: potential impacts and future challenges

The current and cascading effects of global change challenges the interactions both between animal individuals (i.e. social and sexual behaviour) and the environment they inhabit. Amphibians are an ecologically diverse class with a wide range of social and sexual behaviours, making them a compelling model to understand the potential adaptations of animals faced with the effects of human-induced rapid environmental changes (HIREC). Poison frogs (Dendrobatoidea) are a particularly interesting system, as they display diverse social behaviours that are shaped by conspecific and environmental interactions, thus offering a tractable system to investigate how closely related species may respond to the impacts of HIREC. Here, we discuss the potential impacts of global change on poison frog behaviour, and the future challenges this group may face in response to such change. We pay special attention to parental care and territoriality, which are emblematic of this clade, and consider how different species may flexibly respond and adapt to increasingly frequent and diverse anthropogenic stress. More specifically, we hypothesise that some parents may increase care (i.e. clutch attendance and distance travelled for tadpole transport) in HIREC scenarios and that species with more generalist oviposition and tadpole deposition behaviours may fare more positively than their less flexible counterparts; we predict that the latter may either face increased competition for resources limited by HIREC or will be forced to adapt and expand their natural preferences. Likewise, we hypothesise that human-driven habitat alteration will disrupt the acoustic and visual communication systems due to increased noise pollution and/or changes in the surrounding light environment. We highlight the need for more empirical research combining behavioural ecology and conservation to better predict species’ vulnerability to global change and efficiently focus conservation efforts.

A natural experiment identifies an impending ecological trap for a neotropical amphibian in response to extreme weather events

Extreme weather events are predicted to increase as a result of climate change, yet amphibian responses to extreme disturbance events remain understudied, especially in the Neotropics. Recently, an unprecedented windstorm within a protected Costa Rican rainforest opened large light gaps in sites where we have studied behavioral responses of diurnal strawberry poison frogs (Oophaga pumilio) to ultraviolet radiation for nearly two decades. Previous studies demonstrate that O. pumilio selects and defends perches where ultraviolet radiation (UV‐B) is relatively low, likely because of the lethal and sublethal effects of UV‐B. In this natural experiment, we quantified disturbance to O. pumilio habitat, surveyed for the presence of O. pumilio in both high‐disturbance and low‐disturbance areas of the forest, and assessed UV‐B levels and perch selection behavior in both disturbance levels. Fewer frogs were detected in high‐disturbance habitat than in low‐disturbance habitat. In general, frogs were found vocalizing at perches in both disturbance levels, and in both cases, in significantly lower UV‐B levels relative to ambient adjacent surroundings. However, frogs at perches in high‐disturbance areas were exposed to UV‐B levels nearly 10 times greater than males at perches in low‐disturbance areas. Thus, behavioral avoidance of UV‐B may not reduce the risks associated with elevated exposure under these novel conditions, and similarly, if future climate and human‐driven land‐use change lead to sustained analogous environments.

Star finches Neochmia ruficauda have a visual preference for white dot patterns: a possible case of trypophilia

Many animals have polka dot patterns on their body surface, some of which are known to have signalling functions; however, their evolutionary origins remain unclear. Dot patterns can trigger a fear response (trypophobia) in humans and are known to function as aposematic signals in non-human animals, suggesting that dots may deserve attention for biological reasons. Interestingly in many birds, plumage dot patterns serve for social/sexual signalling. To understand their evolution, we have focused on the sensory bias hypothesis, which predicts the role of pre-existing sensory preference driven by natural selection in shaping signal design. Our previous phylogenetic comparative study supported the hypothesis and showed that diet-driven visual preference promoted the evolution of plumage patterns, as there was an evolutionary correlation between termite-eating (white roundish gregarious prey) and the presence of plumage dot patterns in species of the family Estrildidae. This suggests that these species possess an intrinsic preference for dots. To test this, we compared the responses of an Estrildid species with dot plumage pattern (star finch Neochmia ruficauda) towards simultaneously presented monochrome-printed white dot vs white stripe patterns under both food-deprived and -supplied conditions. Overall, star finches preferred dots to stripes. They showed foraging-like behaviours almost only toward dots when hungry and gazed at dots frequently even when food was available, suggesting both hunger-related and hunger-neutral dot preferences. These results are rather surprising, given how strongly the subjects were attracted to abstract dot patterns without organic structure, but provided good support for the sensory bias hypothesis.

Continuous Variation in an Aposematic Pattern Affects Background Contrast, but Is Not Associated With Differences in Microhabitat Use

Variation in aposematic signals was once predicted to be rare, yet in recent years it has become increasingly well documented. Despite increases in the frequency with which polytypism and polymorphism have been suggested to occur, population-wide variance is rarely quantified. We comprehensively sampled a subpopulation of the poison frog Oophaga sylvatica, a species which is polytypic across its distribution and also shows considerable within-population polymorphism. On one hand, color pattern polymorphism could be the result of multifarious selection acting to balance different signaling functions and leading to the evolution of discrete sub-morphs which occupy different fitness peaks. Alternatively, variance could simply be due to relaxed selection, where variation would be predicted to be continuous. We used visual modeling of conspecific and heterospecific observers to quantify the extent of within population phenotypic variation and assess whether this variation produced distinct signals. We found that, despite considerable color pattern variation, variance could not be partitioned into distinct groups, but rather all viewers would be likely to perceive variation as continuous. Similarly, we found no evidence that frog color pattern contrast was either enhanced or diminished in the frogs’ chosen microhabitats compared to alternative patches in which conspecifics were observed. Within population phenotypic variance therefore does not seem to be indicative of strong selection toward multiple signaling strategies, but rather pattern divergence has likely arisen due to weak purifying selection, or neutral processes, on a signal that is highly salient to both conspecifics and predators.

The impact of life stage and pigment source on the evolution of novel warning signal traits

Our understanding of how novel warning color traits evolve in natural populations is largely based on studies of reproductive stages and organisms with endogenously produced pigmentation. In these systems, genetic drift is often required for novel alleles to overcome strong purifying selection stemming from frequency‐dependent predation and positive assortative mating. Here, we integrate data from field surveys, predation experiments, population genomics, and phenotypic correlations to explain the origin and maintenance of geographic variation in a diet‐based larval pigmentation trait in the redheaded pine sawfly (Neodiprion lecontei), a pine‐feeding hymenopteran. Although our experiments confirm that N. lecontei larvae are indeed aposematic—and therefore likely to experience frequency‐dependent predation—our genomic data do not support a historical demographic scenario that would have facilitated the spread of an initially deleterious allele via drift. Additionally, significantly elevated differentiation at a known color locus suggests that geographic variation in larval color is currently maintained by selection. Together, these data suggest that the novel white morph likely spread via selection. However, white body color does not enhance aposematic displays, nor is it correlated with enhanced chemical defense or immune function. Instead, the derived white‐bodied morph is disproportionately abundant on a pine species with a reduced carotenoid content relative to other pine hosts, suggesting that bottom‐up selection via host plants may have driven divergence among populations. Overall, our results suggest that life stage and pigment source can have a substantial impact on the evolution of novel warning signals, highlighting the need to investigate diverse aposematic taxa to develop a comprehensive understanding of color variation in nature.

Contrasting environmental drivers of genetic and phenotypic divergence in an Andean poison frog (Epipedobates anthonyi)

Phenotypic and genetic divergence are shaped by the homogenizing effects of gene flow and the differentiating processes of genetic drift and local adaptation. Herein, we examined the mechanisms that underlie phenotypic (size and color) and genetic divergence in 35 populations (535 individuals) of the poison frog Epipedobates anthonyi along four elevational gradients (0-1800 m asl) in the Ecuadorian Andes. We found phenotypic divergence in size and color despite relatively low genetic divergence at neutral microsatellite loci. Genetic and phenotypic divergence were both explained by landscape resistance between sites (isolation-by-resistance, IBR), likely due to a cold and dry mountain ridge between the northern and southern elevational transects that limits dispersal and separates two color morphs. Moreover, environmental differences among sites also explained genetic and phenotypic divergence, suggesting isolation-by-environment (IBE). When northern and southern transects were analyzed separately, genetic divergence was predicted either by distance (isolation-by-distance, IBD; northern) or environmental resistance between sites (IBR; southern). In contrast, phenotypic divergence was primarily explained by environmental differences among sites, supporting the IBE hypothesis. These results indicate that although distance and geographic barriers are important drivers of population divergence, environmental variation has a two-fold effect on population divergence. On the one hand, landscape resistance between sites reduces gene flow (IBR), while on the other hand, environmental differences among sites exert divergent selective pressures on phenotypic traits (IBE). Our work highlights the importance of studying both genetic and phenotypic divergence to better understand the processes of population divergence and speciation along ecological gradients.

The geometry of decision-making in individuals and collectives

Almost all animals must make decisions on the move. Here, employing an approach that integrates theory and high-throughput experiments (using state-of-the-art virtual reality), we reveal that there exist fundamental geometrical principles that result from the inherent interplay between movement and organisms’ internal representation of space. Specifically, we find that animals spontaneously reduce the world into a series of sequential binary decisions, a response that facilitates effective decision-making and is robust both to the number of options available and to context, such as whether options are static (e.g., refuges) or mobile (e.g., other animals). We present evidence that these same principles, hitherto overlooked, apply across scales of biological organization, from individual to collective decision-making.

Choosing among spatially distributed options is a central challenge for animals, from deciding among alternative potential food sources or refuges to choosing with whom to associate. Using an integrated theoretical and experimental approach (employing immersive virtual reality), we consider the interplay between movement and vectorial integration during decision-making regarding two, or more, options in space. In computational models of this process, we reveal the occurrence of spontaneous and abrupt “critical” transitions (associated with specific geometrical relationships) whereby organisms spontaneously switch from averaging vectorial information among, to suddenly excluding one among, the remaining options. This bifurcation process repeats until only one option—the one ultimately selected—remains. Thus, we predict that the brain repeatedly breaks multichoice decisions into a series of binary decisions in space–time. Experiments with fruit flies, desert locusts, and larval zebrafish reveal that they exhibit these same bifurcations, demonstrating that across taxa and ecological contexts, there exist fundamental geometric principles that are essential to explain how, and why, animals move the way they do.

The evolution of age-specific choosiness and reproductive isolation in a model with overlapping generations

The strength of mate choice (choosiness) often varies with age, but theory to understand this variation is scarce. Additionally, theory has investigated the evolution of choosiness in speciation scenarios but has ignored that most organisms have overlapping generations. We investigate whether speciation can result in variation of choosiness with age, and whether such variation can in turn affect speciation. We develop a population-genetic model of the evolution of choosiness in organisms with overlapping generations in the context of secondary contact between two divergent populations. We assume that females choose males that match their phenotype, such that choosiness evolves by sexual selection. We demonstrate that speciation can result in the evolution of age-specific choosiness when the mating trait is under divergent ecological selection and age is not used as a mating cue. The cause of this result is that allele frequencies differ between choosy females and males. However, we find that the evolution of age-specific choosiness does not affect the overall level of reproductive isolation compared to a case without age-structure, supporting previous speciation theory. Overall, our results connect life history and speciation theory, and the mechanisms that we highlight have implications for the understanding of the role of sex-specific selection in the evolution of choosiness.

The influence of ultraviolet reflectance differs between conspicuous aposematic signals in neotropical butterflies and poison frogs

Warning signals are often characterized by highly contrasting, distinctive, and memorable colors. Greater chromatic (hue) and achromatic (brightness) contrast have both been found to contribute to greater signal efficacy, making longwave colored signals (e.g., red and yellow), that are perceived by both chromatic and achromatic visual pathways, particularly common. Conversely, shortwave colors (e.g., blue and ultraviolet) do not contribute to luminance perception yet are also commonly found in warning signals. Our understanding of the role of UV in aposematic signals is currently incomplete as UV perception is not universal, and evidence for its utility is at best mixed. We used visual modeling to quantify how UV affects signal contrast in aposematic heliconiian butterflies and poison frogs both of which reflect UV wavelengths, occupy similar habitats, and share similar classes of predators. Previous work on butterflies has found that UV reflectance does not affect predation risk but is involved in mate choice. As the butterflies, but not the frogs, have UV-sensitive vision, the function of UV reflectance in poison frogs is currently unknown. We found that despite showing up strongly in UV photographs, UV reflectance only appreciably affected visual contrast in the butterflies. As such, these results support the notion that although UV reflectance is associated with intraspecific communication in butterflies, it appears to be nonfunctional in frogs. Consequently, our data highlight that we should be careful when assigning a selection-based benefit to the presence of UV reflectance.

Aposematism facilitates the diversification of parental care strategies in poison frogs

Many organisms have evolved adaptations to increase the odds of survival of their offspring. Parental care has evolved several times in animals including ectotherms. In amphibians, ~ 10% of species exhibit parental care. Among these, poison frogs (Dendrobatidae) are well-known for their extensive care, which includes egg guarding, larval transport, and specialized tadpole provisioning with trophic eggs. At least one third of dendrobatids displaying aposematism by exhibiting warning coloration that informs potential predators about the presence of defensive skin toxins. Aposematism has a central role in poison frog diversification, including diet specialization, and visual and acoustic communication; and it is thought to have impacted their reproductive biology as well. We tested the latter association using multivariate phylogenetic methods at the family level. Our results show complex relationships between aposematism and certain aspects of the reproductive biology in dendrobatids. In particular, aposematic species tend to use more specialized tadpole-deposition sites, such as phytotelmata, and ferry fewer tadpoles than non-aposematic species. We propose that aposematism may have facilitated the diversification of microhabitat use in dendrobatids in the context of reproduction. Furthermore, the use of resource-limited tadpole-deposition environments may have evolved in tandem with an optimal reproductive strategy characterized by few offspring, biparental care, and female provisioning of food in the form of unfertilized eggs. We also found that in phytotelm-breeders, the rate of transition from cryptic to aposematic phenotype is 17 to 19 times higher than vice versa. Therefore, we infer that the aposematism in dendrobatids might serve as an umbrella trait for the evolution and maintenance of their complex offspring-caring activities.

Speciation rates are positively correlated with the rate of plumage color evolution in hummingbirds

A fascinating pattern in nature is the uneven distribution of biodiversity among clades, some with low species richness and phenotypic variation in contrast to others with remarkable species richness and phenotypic diversity. In animals, communication signals are crucial for intra- and interspecific interactions and are likely an important factor in speciation. However, evidence for the association between the evolution of such signals and speciation is mixed. In hummingbirds, plumage coloration is an important communication signal, particularly for mate selection. Here, using reflectance data for 237 hummingbird species (∼66% of total diversity), we demonstrate that color evolution rates are associated with speciation rates, and that differences among feather patches are consistent with an interplay between natural and sexual selection. We found that female color evolution rates of multiple plumage elements, including the gorget, were similar to those of males. Although male color evolution in this patch was associated with speciation, female gorget color evolution was not. In other patches, the relationship between speciation and color evolution rates was pervasive between sexes. We anticipate that future studies on animal communication will likely find that evolution of signaling traits of both sexes has played a vital role in generating signal and species diversity.

Habitat disturbance alters color contrast and the detectability of cryptic and aposematic frogs

Animals use color both to conceal and signal their presence, with patterns that match the background, disrupt shape recognition, or highlight features important for communication. The forms that these color patterns take are responses to the visual systems that observe them and the environments within which they are viewed. Increasingly, however, these environments are being affected by human activity. We studied how pattern characteristics and habitat change may affect the detectability of three frog color patterns from the Bocas del Toro archipelago in Panama: Beige-Striped Brown Allobates talamancae and two spotted morphs of Oophaga pumilio, Black-Spotted Green and Black-Spotted Red. To assess detectability, we used visual modeling of conspecifics and potential predators, along with a computer-based detection experiment with human participants. Although we found no evidence for disruptive camouflage, we did find clear evidence that A. talamancae stripes are inherently more cryptic than O. pumilio spots regardless of color. We found no evidence that color pattern polytypism in O. pumilio is related to differences in the forest floor between natural sites. We did, however, find strong evidence that human disturbance affects the visual environment and modifies absolute and rank order frog detectability. Human-induced environmental change reduces the effectiveness of camouflage in A. talamancae, reduces detectability of Black-Spotted Green O. pumilio, and increases chromatic contrast, but not detectability, in Black-Spotted Red O. pumilio. Insofar as predators may learn about prey defenses and make foraging decisions based on relative prey availability and suitability, such changes may have wider implications for predator–prey dynamics.

Light environment influences mating behaviours during the early stages of divergence in tropical butterflies

Speciation is facilitated when traits under divergent selection also act as mating cues. Fluctuations in sensory conditions can alter signal perception independently of adaptation to the broader sensory environment, but how this fine-scale variation may constrain or promote behavioural isolation has received little attention. The warning patterns of Heliconius butterflies are under selection for aposematism and act as mating cues. Using computer vision, we extracted behavioural data from 1481 h of video footage, for 387 individuals. We show that the putative hybrid species H. heurippa and its close relative H. timareta linaresi differ in their response to divergent warning patterns, but that these differences are strengthened with increased local illuminance. Trials with live individuals reveal low-level assortative mating that is sufficiently explained by differences in visual attraction. Finally, results from hybrid butterflies are consistent with linkage between a major warning pattern gene and the corresponding behaviour, though the differences in behaviour we observe are unlikely to cause rapid reproductive isolation as predicted under a model of hybrid trait speciation. Overall, our results reveal that the contribution of ecological mating cues to reproductive isolation may depend on the immediate sensory conditions during which they are displayed to conspecifics.

Behavioral Variation in the Pygmy Halfbeak Dermogenys collettei: Comparing Shoals With Contrasting Ecologies

Variation in biotic and abiotic factors among populations affects individual behaviors by transforming the social landscape and shaping mating systems. Consequently, describing behaviors in natural populations requires consideration of the biological and physical factors that different individuals face. Here, we examined variation in socio-sexual and locomotor behaviors in a small, livebearing, freshwater fish, the pygmy halfbeak Dermogenys collettei, across natural populations in Singapore. The pygmy halfbeak is a surface feeding fish that spends most of the time near the water surface, making it ideal for non-invasive behavioral observations. We compared behaviors between sexes among 26 shoals while simultaneously accounting for environmental variation. We demonstrated that sexual interactions and locomotor behaviors differed among shoals with varying levels of canopy cover and water flow. Specifically, in areas with greater canopy cover, sexual interactions decreased, whereas time spent in a stationary position increased. Sexual interactions were more numerous in still water, where fish spent less time swimming. Variation in the expression of socio-sexual and locomotor behaviors were not associated with differences in the amount of aquatic vegetation, water depth or halfbeak shoal size. Agonistic interactions were robust to environmental effects, showing little variation among environments. However, there were strong sex effects, with males performing more agonistic behaviors and spending less time in a stationary position compared to females, regardless of the environment. Moreover, sexual interactions, measured as actively performed by males and passively received by females, were on average more frequent in males than in females. Our findings help us explore the proximal causes of intraspecific behavioral variation and suggest that fundamental information on socio-sexual behaviors from wild populations can lead to a better understanding of how sexual selection operates when the strength of natural selection varies across environments.

Ultraviolet components offer minimal contrast enhancement to an aposematic signal

Aposematic and sexual signals are often characterized by bright, highly contrasting colors. Many species can see colors beyond the human visible spectrum, and ultraviolet (UV) reflection has been found to play an important role in communication and sexual selection. However, the role of UV in aposematic signals is poorly explored. Poison frogs frequently produce high-contrast signals that have been linked to both aposematism and intraspecific communication. Yet despite considerable efforts studying interspecific and intraspecific diversity in color, poison frogs are not known to perceive UV, and UV reflection of the integument has not been described. We report UV-reflective spots in a population of Oophaga sylvatica and quantify the effect of UV on visual contrast with models of avian vision. We found that the frogs are highly contrasting, but UV had a minimal effect on signal saliency. These data highlight the importance of considering UV reflectance within aposematic signals, but that UV should not necessarily be regarded as an independent signal.

Shifts in the developmental rate of spadefoot toad larvae cause decreased complexity of post-metamorphic pigmentation patterns

Amphibian larvae are plastic organisms that can adjust their growth and developmental rates to local environmental conditions. The consequences of such developmental alterations have been studied in detail, both at the phenotypic and physiological levels. While largely unknown, it is of great importance to assess how developmental alterations affect the pigmentation pattern of the resulting metamorphs, because pigmentation is relevant for communication, mate choice, and camouflage and hence influences the overall fitness of the toads. Here we quantify the variation in several aspects of the pigmentation pattern of juvenile spadefoot toads experimentally induced to accelerate their larval development in response to decreased water level. It is known that induced developmental acceleration comes at the cost of reduced size at metamorphosis, higher metabolic rate, and increased oxidative stress. In this study, we show that spadefoot toads undergoing developmental acceleration metamorphosed with a less complex, more homogeneous, darker dorsal pattern consisting of continuous blotches, compared to the more contrasted pattern with segregated blotches and higher fractal dimension in normally developing individuals, and at a smaller size. We also observed a marked effect of population of origin in the complexity of the pigmentation pattern. Complexity of the post-metamorphic dorsal pigmentation could therefore be linked to pre-metamorphic larval growth and development.

Male–male contest limits the expression of assortative mate preferences in a polymorphic poison frog

Codivergence of sexual traits and mate preferences can lead to assortative mating and subsequently reproductive isolation. However, mate choice rarely operates without intrasexual competition, and the effects of the latter on speciation are often overlooked. Maintaining trait polymorphisms despite gene flow and limiting assortative female preferences for less-competitive male phenotypes are two important roles that male–male competition may play in the speciation process. Both roles rely on the assumption that male–male competition limits the expression of divergent female preferences. We tested this assumption in the highly color-polymorphic strawberry poison frog (Oophaga pumilio). Females prefer males of the local color, suggesting that reproductive isolation may be evolving among color morphs. However, this inference does not account for male–male competition, which is also color-mediated. We housed females with two differently colored males, and compared reproductive patterns when the more attractive male was the territory holder versus when he was the nonterritorial male. Females mated primarily with the territory winner, regardless of coloration, suggesting that when a choice must be made between the two, male territoriality overrides female preferences for male coloration. Our results highlight the importance of considering the combined effects of mate choice and intrasexual competition in shaping phenotypic divergence and speciation.

Divergence, gene flow, and the origin of leapfrog geographic distributions: The history of colour pattern variation in Phyllobates poison-dart frogs

The geographic distribution of phenotypic variation among closely related populations is a valuable source of information about the evolutionary processes that generate and maintain biodiversity. Leapfrog distributions, in which phenotypically similar populations are disjunctly distributed and separated by one or more phenotypically distinct populations, represent geographic replicates for the existence of a phenotype, and are therefore especially informative. These geographic patterns have mostly been studied from phylogenetic perspectives to understand how common ancestry and divergent evolution drive their formation. Other processes, such as gene flow between populations, have not received as much attention. Here, we investigate the roles of divergence and gene flow between populations in the origin and maintenance of a leapfrog distribution in Phyllobates poison frogs. We found evidence for high levels of gene flow between neighbouring populations but not over long distances, indicating that gene flow between populations exhibiting the central phenotype may have a homogenizing effect that maintains their similarity, and that introgression between 'leapfroging' taxa has not played a prominent role as a driver of phenotypic diversity in Phyllobates. Although phylogenetic analyses suggest that the leapfrog distribution was formed through independent evolution of the peripheral (i.e. leapfrogging) populations, the elevated levels of gene flow between geographically close populations poise alternative scenarios, such as the history of phenotypic change becoming decoupled from genome-averaged patterns of divergence, which we cannot rule out. These results highlight the importance of incorporating gene flow between populations into the study of geographic variation in phenotypes, both as a driver of phenotypic diversity and as a confounding factor of phylogeographic inferences.

A ketocarotenoid-based colour polymorphism in the Sira poison frog Ranitomeya sirensis indicates novel gene interactions underlying aposematic signal variation

The accumulation of red ketocarotenoids is an important component of coloration in many organisms, but the underlying mechanisms are poorly understood. In some organisms, ketocarotenoids are sequestered from the diet and can accumulate when enzymes responsible for carotenoid breakdown are disrupted. In other organisms, ketocarotenoids are formed endogenously from dietary precursors via oxidation reactions carried out by carotenoid ketolase enzymes. Here, we study the genetic basis of carotenoid coloration in an amphibian. We demonstrate that a red/yellow polymorphism in the dendrobatid poison frog Ranitomeya sirensis is due to the presence/absence of ketocarotenoids. Using whole-transcriptome sequencing of skins and livers, we found that a transcript encoding a cytochrome P450 enzyme (CYP3A80) is expressed 3.4-fold higher in livers of red frogs versus yellow. As CYP3A enzymes are known carotenoid ketolases in other organisms, our results point to CYP3A80 as a strong candidate for a carotenoid ketolase in amphibians. Furthermore, in red frogs, the transcript encoding the carotenoid cleavage enzyme BCO2 is expressed at a low level or as a splice variant lacking key catalytic amino acids. This suggests that BCO2 function may be disrupted in red frogs, providing a mechanism whereby the accumulation of ketocarotenoids and their dietary precursors may be enhanced.

Plastic sexual ornaments: Assessing temperature effects on color metrics in a color-changing reptile

Conspicuous coloration is an important subject in social communication and animal behavior, and it can provide valuable insight into the role of visual signals in social selection. However, animal coloration can be plastic and affected by abiotic factors such as temperature, making its quantification problematic. In such cases, careful consideration is required so that metric choices are consistent across environments and least sensitive to abiotic factors. A detailed assessment of plastic trait in response to environmental conditions could help identify more robust methods for quantifying color. Temperature affects sexual ornamentation of eastern fence lizards, Sceloporus undulatus, with ventral coloration shifting from green to blue hues as temperatures rise, making the calculation of saturation (color purity) difficult under conditions where temperatures vary. We aimed to characterize how abiotic factors influence phenotypic expression and to identify a metric for quantifying animal color that is either independent from temperature (ideally) or best conserves individual's ranks. We compared the rates of change in saturation across two temperature treatments using seven metrics: three that are based on fixed spectral ranges (with two of them designed by us specifically for this system) and three that track the expressed hue (with one of them designed by us to circumvent spurious results in unornamented individuals). We also applied a lizard visual sensitivity model to understand how temperature-induced color changes may be perceived by conspecifics. We show that the rate of change in saturation between two temperatures is inconsistent across individuals, increasing at a higher rate in individuals with higher baseline saturation at lower temperatures. In addition, the relative color rank of individuals in a population varies with the temperature standardized by the investigator, but more so for some metrics than others. While we were unable to completely eliminate the effect of temperature, current tools for quantifying color allowed us to use spectral data to estimate saturation in a variety of ways and to largely preserve saturation ranks of individuals across temperatures while avoiding erroneous color scores. We describe our approaches and suggest best-practices for quantifying and interpreting color, particularly in cases where color changes in response to environmental factors.

Being red, blue and green: the genetic basis of coloration differences in the strawberry poison frog (Oophaga pumilio)

Background

Animal coloration is usually an adaptive attribute, under strong local selection pressures and often diversified among species or populations. The strawberry poison frog (Oophaga pumilio) shows an impressive array of color morphs across its distribution in Central America. Here we quantify gene expression and genetic variation to identify candidate genes involved in generating divergence in coloration between populations of red, green and blue O. pumilio from the Bocas del Toro archipelago in Panama.

Results

We generated a high quality non-redundant reference transcriptome by mapping the products of genome-guided and de novo transcriptome assemblies onto a re-scaffolded draft genome of O. pumilio. We then measured gene expression in individuals of the three color phenotypes and identified color-associated candidate genes by comparing differential expression results against a list of a priori gene sets for five different functional categories of coloration – pteridine synthesis, carotenoid synthesis, melanin synthesis, iridophore pathways (structural coloration), and chromatophore development. We found 68 candidate coloration loci with significant expression differences among the color phenotypes. Notable upregulated examples include pteridine synthesis genes spr, xdh and pts (in red and green frogs); carotenoid metabolism genes bco2 (in blue frogs), scarb1 (in red frogs), and guanine metabolism gene psat1 (in blue frogs). We detected significantly higher expression of the pteridine synthesis gene set in red and green frogs versus blue frogs. In addition to gene expression differences, we identified 370 outlier SNPs on 162 annotated genes showing signatures of diversifying selection, including eight pigmentation-associated genes.

Conclusions

Gene expression in the skin of the three populations of frogs with differing coloration is highly divergent. The strong signal of differential expression in pteridine genes is consistent with a major role of these genes in generating the coloration differences among the three morphs. However, the finding of differentially expressed genes across pathways and functional categories suggests that multiple mechanisms are responsible for the coloration differences, likely involving both pigmentary and structural coloration. In addition to regulatory differences, we found potential evidence of differential selection acting at the protein sequence level in several color-associated loci, which could contribute to the color polymorphism.

Mechanisms for Color Convergence in a Mimetic Radiation of Poison Frogs

In animals, bright colors often evolve to mimic other species when a resemblance is selectively favored. Understanding the proximate mechanisms underlying such color mimicry can give insights into how mimicry evolves-for example, whether color convergence evolves from a shared set of mechanisms or through the evolution of novel color production mechanisms. We studied color production mechanisms in poison frogs (Dendrobatidae), focusing on the mimicry complex of Ranitomeya imitator. Using reflectance spectrometry, skin pigment analysis, electron microscopy, and color modeling, we found that the bright colors of these frogs, both within and outside the mimicry complex, are largely structural and produced by iridophores but that color production depends crucially on interactions with pigments. Color variation and mimicry are regulated predominantly by iridophore platelet thickness and, to a lesser extent, concentration of the red pteridine pigment drosopterin. Compared with each of the four morphs of model species that it resembles, R. imitator displays greater variation in both structural and pigmentary mechanisms, which may have facilitated phenotypic divergence in this species. Analyses of nonmimetic dendrobatids in other genera demonstrate that these mechanisms are widespread within the family and that poison frogs share a complex physiological "color palette" that can produce diverse and highly reflective colors.

Complex patterns of differentiation and gene flow underly the divergence of aposematic phenotypes in Oophaga poison frogs

Hybridization and introgression can have complex consequences for both species evolution and conservation. Here, we investigated the origin and characteristics of a putative hybrid zone between two South American poison dart frog species, Oophaga anchicayensis and the critically endangered Oophaga lehmanni, which are heavily sought after on the illegal pet market. Using a combination of phenotypic (49 traits) and genomic (ddRADseq) data, we found that the putative hybrids are morphologically distinct from their parental species and confirmed genomic signatures of admixture in these populations. Several lines of evidence (hybrid indices, interspecific hybrid heterozygosity, genomic clines, comparisons with simulated hybrids and demographic modelling) support the conclusion that these populations are not comprised of early-generation hybrids and thus, they probably did not arise as a result of illegal translocations associated with wildlife trafficking. Instead, they probably represent an independent lineage which has persisted through isolation and has only relatively recently re-established gene flow with both parental species. Furthermore, we detected signals of differential introgression from parental species into these hybrid populations which suggest relaxed stabilizing selection on these aposematic colour morphs, potentially via context-dependent female choice. These populations thus provide a fascinating window into the role of hybridization, isolation and female choice in the diversification of South American poison dart frogs. In addition, our results underline the importance of landscape conservation measures to protect, not only known localities of nominal species, but also the phenotypic and genomic variation harbored by admixed lineages which represent crucial repositories for the impressive diversity in this system.

Reproductive State Modulates Retinal Sensitivity to Light in Female Túngara Frogs

Visual cues are often a vital part of animal communication and courtship. While a plethora of studies have focused on the role that hormones play in acoustic communication of anurans, relatively few have explored hormonal modulation of vision in these animals. Much of what we do know comes from behavioral studies, which show that a frog’s hormonal state can significantly affect both its visual behavior and mating decisions. However, to fully understand how frogs use visual cues to make these mating decisions, we must first understand how their visual system processes these cues, and how hormones affect these processes. To do this, we performed electroretinograms (ERGs) to measure retinal sensitivity of túngara frogs (Physalaemus pustulosus), a neotropical species whose mating behavior includes previously described visual cues. To determine the effect of hormonal state on visual sensitivity, ERGs were recorded under scotopic and photopic conditions in frogs that were either non-reproductive or hormone-treated with human chorionic gonadotropin (hCG) prior to testing. Additionally, measurements of optical anatomy determined how túngara frog eye and retina morphology related to physiological sensitivity. As expected, we found that both sexes display higher visual sensitivity under scotopic conditions compared to photopic conditions. However, hormone injections significantly increased retinal sensitivity of females under scotopic conditions. These results support the hypothesis that hormonal modulation of neural mechanisms, such as those mediating visually guided reproductive behavior in this species, include modulation of the receptor organ: the retina. Thus, our data serve as a starting point for elucidating the mechanism of hormonal modulation of visual sensitivity.

Links between prey assemblages and poison frog toxins: A landscape ecology approach to assess how biotic interactions affect species phenotypes

Ecological studies of species pairs showed that biotic interactions promote phenotypic change and eco-evolutionary feedbacks. However, it is unclear how phenotypes respond to synergistic interactions with multiple taxa. We investigate whether interactions with multiple prey species explain spatially structured variation in the skin toxins of the neotropical poison frog Oophaga pumilio. Specifically, we assess how dissimilarity (i.e., beta diversity) of alkaloid-bearing arthropod prey assemblages (68 ant species) and evolutionary divergence between frog populations (from a neutral genetic marker) contribute to frog poison dissimilarity (toxin profiles composed of 230 different lipophilic alkaloids sampled from 934 frogs at 46 sites). We find that models that incorporate spatial turnover in the composition of ant assemblages explain part of the frog alkaloid variation, and we infer unique alkaloid combinations across the range of O. pumilio. Moreover, we find that alkaloid variation increases weakly with the evolutionary divergence between frog populations. Our results pose two hypotheses: First, the distribution of only a few prey species may explain most of the geographic variation in poison frog alkaloids; second, different codistributed prey species may be redundant alkaloid sources. The analytical framework proposed here can be extended to other multitrophic systems, coevolutionary mosaics, microbial assemblages, and ecosystem services.

Preferences for and use of light microhabitats differ among and within populations of a polytypic poison frog

Anti-predator strategies can influence trade-offs governing other activities important to fitness. Crypsis, for example, might make conspicuous sexual display especially costly, whereas aposematism might reduce or remove such costs. We tested for correlates of anti-predator strategy in Oophaga pumilio, a polytypic poison frog with morphs spanning the crypsis–aposematism continuum. In the wild, males of visually conspicuous morphs display from conspicuous perches and behave as if they perceive predation risk to be low. We thus predicted that, given a choice of ambient light microhabitats, these males would use high ambient light conditions the most and be most likely to perch in high-light conditions. We found no evidence that differently colored male O. pumilio preferentially used bright microhabitats or that ambient light influenced perching in a morph-specific manner. Independent of light conditions, males from the most conspicuous population perched the least, but the most conspicuous individuals from a polymorphic population perched the most. These patterns suggest that preferences do not necessarily underlie among-morph differences observed in the wild. This could be explained, and remain consistent with theory, if risk aversion is shaped, in part, by experience.

Coevolution of male and female mate choice can destabilize reproductive isolation

Sexual interactions play an important role in the evolution of reproductive isolation, with important consequences for speciation. Theoretical studies have focused on the evolution of mate preferences in each sex separately. However, mounting empirical evidence suggests that premating isolation often involves mutual mate choice. Here, using a population genetic model, we investigate how female and male mate choice coevolve under a phenotype matching rule and how this affects reproductive isolation. We show that the evolution of female preferences increases the mating success of males with reciprocal preferences, favouring mutual mate choice. However, the evolution of male preferences weakens indirect selection on female preferences and, with weak genetic drift, the coevolution of female and male mate choice leads to periodic episodes of random mating with increased hybridization (deterministic 'preference cycling' triggered by stochasticity). Thus, counterintuitively, the process of establishing premating isolation proves rather fragile if both male and female mate choice contribute to assortative mating.

Phylogenetic relationships and systematics of the Amazonian poison frog genus Ameerega using ultraconserved genomic elements

The Amazonian poison frog genus Ameerega is one of the largest yet most understudied of the brightly colored genera in the anuran family Dendrobatidae, with 30 described species ranging throughout tropical South America. Phylogenetic analyses of Ameerega are highly discordant, lacking consistency due to variation in data types and methods, and often with limited coverage of species diversity in the genus. Here, we present a comprehensive phylogenomic reconstruction of Ameerega, utilizing state-of-the-art sequence capture techniques and phylogenetic methods. We sequenced thousands of ultraconserved elements from over 100 tissue samples, representing almost every described Ameerega species, as well as undescribed cryptic diversity. We generated topologies using maximum likelihood and coalescent methods and compared the use of maximum likelihood and Bayesian methods for estimating divergence times. Our phylogenetic inference diverged strongly from those of previous studies, and we recommend steps to bring Ameerega taxonomy in line with the new phylogeny. We place several species in a phylogeny for the first time, as well as provide evidence for six potential candidate species. We estimate that Ameerega experienced a rapid radiation approximately 7-11 million years ago and that the ancestor of all Ameerega was likely an aposematic, montane species. This study underscores the utility of phylogenomic data in improving our understanding of the phylogeny of understudied clades and making novel inferences about their evolution.

From habitat use to social behavior: natural history of a voiceless poison frog, Dendrobates tinctorius

Descriptive studies of natural history have always been a source of knowledge on which experimental work and scientific progress rely. Poison frogs are a well-studied group of small Neotropical frogs with diverse parental behaviors, distinct calls, and bright colors that warn predators about their toxicity; and a showcase of advances in fundamental biology through natural history observations. The dyeing poison frog, Dendrobates tinctorius, is emblematic of the Guianas region, widespread in the pet trade, and increasingly popular in research. This species shows several unusual behaviors, such as the lack of advertisement calls and the aggregation around tree-fall gaps, which remain poorly described and understood. Here, we summarize our observations from a natural population of D. tinctorius in French Guiana collected over various field trips between 2009 and 2017; our aim is to provide groundwork for future fundamental and applied research spanning parental care, animal dispersal, disease spread, habitat use in relation to color patterns, and intra-specific communication, to name a few. We report sex differences in habitat use and the striking invasion of tree-fall gaps; describe their courtship and aggressive behaviors; document egg development and tadpole transport; and discuss how the knowledge generated by this study could set the grounds for further research on the behavior, ecology, and conservation of this species.

Phylogenomic Reconstruction of the Neotropical Poison Frogs (Dendrobatidae) and Their Conservation

The evolutionary history of the Dendrobatidae, the charismatic Neotropical poison frog family, remains in flux, even after a half-century of intensive research. Understanding the evolutionary relationships between dendrobatid genera and the larger-order groups within Dendrobatidae is critical for making accurate assessments of all aspects of their biology and evolution. In this study, we provide the first phylogenomic reconstruction of Dendrobatidae with genome-wide nuclear markers known as ultraconserved elements. We performed sequence capture on 61 samples representing 33 species across 13 of the 16 dendrobatid genera, aiming for a broadly representative taxon sample. We compare topologies generated using maximum likelihood and coalescent methods and estimate divergence times using Bayesian methods. We find most of our dendrobatid tree to be consistent with previously published results based on mitochondrial and low-count nuclear data, with notable exceptions regarding the placement of Hyloxalinae and certain genera within Dendrobatinae. We also characterize how the evolutionary history and geographic distributions of the 285 poison frog species impact their conservation status. We hope that our phylogeny will serve as a backbone for future evolutionary studies and that our characterizations of conservation status inform conservation practices while highlighting taxa in need of further study.

Genomic Takeover by Transposable Elements in the Strawberry Poison Frog

We sequenced the genome of the strawberry poison frog, Oophaga pumilio, at a depth of 127.5× using variable insert size libraries. The total genome size is estimated to be 6.76 Gb, of which 4.76 Gb are from high copy number repetitive elements with low differentiation across copies. These repeats encompass DNA transposons, RNA transposons, and LTR retrotransposons, including at least 0.4 and 1.0 Gb of Mariner/Tc1 and Gypsy elements, respectively. Expression data indicate high levels of gypsy and Mariner/Tc1 expression in ova of O. pumilio compared with Xenopus laevis. We further observe phylogenetic evidence for horizontal transfer (HT) of Mariner elements, possibly between fish and frogs. The elements affected by HT are present in high copy number and are highly expressed, suggesting ongoing proliferation after HT. Our results suggest that the large amphibian genome sizes, at least partially, can be explained by a process of repeated invasion of new transposable elements that are not yet suppressed in the germline. We also find changes in the spliceosome that we hypothesize are related to permissiveness of O. pumilio to increases in intron length due to transposon proliferation. Finally, we identify the complement of ion channels in the first genomic sequenced poison frog and discuss its relation to the evolution of autoresistance to toxins sequestered in the skin.

Phenotypic and genetic diversity in aposematic Malagasy poison frogs (genus Mantella)

Intraspecific color variation has long fascinated evolutionary biologists. In species with bright warning coloration, phenotypic diversity is particularly compelling because many factors, including natural and sexual selection, contribute to intraspecific variation. To better understand the causes of dramatic phenotypic variation in Malagasy poison frogs, we quantified genetic structure and color and pattern variation across three closely related species, Mantella aurantiaca, Mantella crocea, and Mantella milotympanum. Although our restriction site-associated DNA (RAD) sequencing approach identified clear genetic clusters, they do not align with current species designations, which has important conservation implications for these imperiled frogs. Moreover, our results suggest that levels of intraspecific color variation within this group have been overestimated, while species diversity has been underestimated. Within major genetic clusters, we observed distinct patterns of variation including: populations that are phenotypically similar yet genetically distinct, populations where phenotypic and genetic breaks coincide, and populations that are genetically similar but have high levels of within-population phenotypic variation. We also detected admixture between two of the major genetic clusters. Our study suggests that several mechanisms-including hybridization, selection, and drift-are contributing to phenotypic diversity. Ultimately, our work underscores the need for a reevaluation of how polymorphic and polytypic populations and species are classified, especially in aposematic organisms.

Mate Choice versus Mate Preference: Inferences about Color-Assortative Mating Differ between Field and Lab Assays of Poison Frog Behavior

Codivergence of mating traits and mate preferences can lead to behavioral isolation among lineages in early stages of speciation. However, mate preferences limit gene flow only when expressed as mate choice, and numerous factors might be more important than preferences in nature. In the extremely color polytypic strawberry poison frog (Oophaga pumilio), female mate preferences have codiverged with color in most allopatric populations tested. Whether these lab-assayed preferences predict mating (gene flow) in the wild remains unclear. We observed courting pairs in a natural contact zone between red and blue lineages until oviposition or courtship termination. We found color-assortative mating in a disturbed habitat with high population density but not in a secondary forest with lower density. Our results suggest color-assortative O. pumilio mate choice in the wild but also mating patterns that do not match those predicted by lab-assayed preferences.

Quantifying the relationship between optical anatomy and retinal physiological sensitivity: A comparative approach

Light intensity varies 1 million-fold between night and day, driving the evolution of eye morphology and retinal physiology. Despite extensive research across taxa showing anatomical adaptations to light niches, surprisingly few empirical studies have quantified the relationship between such traits and the physiological sensitivity to light. In this study, we employ a comparative approach in frogs to determine the physiological sensitivity of eyes in two nocturnal (Rana pipiens, Hyla cinerea) and two diurnal species (Oophaga pumilio, Mantella viridis), examining whether differences in retinal thresholds can be explained by ocular and cellular anatomy. Scotopic electroretinogram (ERG) analysis of relative b-wave amplitude reveals 10- to 100-fold greater light sensitivity in nocturnal compared to diurnal frogs. Ocular and cellular optics (aperture, focal length, and rod outer segment dimensions) were assessed via the Land equation to quantify differences in optical sensitivity. Variance in retinal thresholds was overwhelmingly explained by Land equation solutions, which describe the optical sensitivity of single rods. Thus, at the b-wave, stimulus-response thresholds may be unaffected by photoreceptor convergence (which create larger, combined collecting areas). Follow-up experiments were conducted using photopic ERGs, which reflect cone vision. Under these conditions, the relative difference in thresholds was reversed, such that diurnal species were more sensitive than nocturnal species. Thus, photopic data suggest that rod-specific adaptations, not ocular anatomy (e.g., aperture and focal distance), drive scotopic thresholds differences. To the best of our knowledge, these data provide the first quantified relationship between optical and physiological sensitivity in vertebrates active in different light regimes.

Uniformity in premating reproductive isolation along an intraspecific cline

Premating reproductive isolation (RI) may reduce gene flow across populations that have differentiated in traits important for mate choice. Examining RI across genetic and phenotypic clines can inform the fundamental evolutionary processes that underlie population and lineage differentiation. We conducted female mate-choice studies across an intraspecific red-eyed treefrog cline in Costa Rica and Panama with 2 specific aims: (1) to characterize RI across the cline and examine the relationship between premating RI and genetic and phenotypic distance and (2) to evaluate our results within a broader evolutionary and taxonomic perspective through examination of other RI studies. We found that female red-eyed treefrogs prefer local males relative to non-local males, indicating that some premating RI has evolved in this system, but that preference strength is not associated with phenotypic or geographic distance. Our analysis of 65 other studies revealed no clear pattern between the strength of RI and geographic distribution (allopatry, parapatry, cline) or phenotypic distance, but revealed extreme variation and overlap in levels of intra- and interspecific levels of RI. This work contributes to a growing body of literature that examines intraspecific RI across a cline to understand the selective processes that shape evolutionary patterns at the earliest stages of divergence.