Reproductive isolation related to mimetic divergence in the poison frog Ranitomeya imitator

Evidence for individual vocal recognition in a pair-bonding poison frog, Ranitomeya imitator

Individually distinctive vocalizations are widespread in nature, although the ability of receivers to discriminate these signals has only been explored through limited taxonomic and social lenses. Here, we asked whether anuran advertisement calls, typically studied for their role in territory defense and mate attraction, facilitate recognition and preferential association with partners in a pair-bonding poison frog (Ranitomeya imitator). Combining no- and two-stimulus choice playback experiments, we evaluated behavioral responses of females to male acoustic stimuli. Virgin females oriented to and approached speakers broadcasting male calls independent of caller identity, implying that females are generally attracted to male acoustic stimuli outside the context of a pair bond. When pair-bonded females were presented with calls of a mate and a stranger, they showed significant preference for calls of their mate. Moreover, behavioral responses varied with breeding status: females with eggs were faster to approach stimuli than females that were pair bonded but did not currently have eggs. Our study suggests a potential role for individual vocal recognition in the formation and maintenance of pair bonds in a poison frog and raises new questions about how acoustic signals are perceived in the context of monogamy and biparental care.

A histological analysis of coloration in the Peruvian mimic poison frog (Ranitomeya imitator)

Aposematism continues to be a phenomenon of central interest in evolutionary biology. The life history of the mimic poison frog, Ranitomeya imitator, relies heavily on aposematism. In order for aposematic signals to be effective, predators must be able to learn to avoid the associated phenotype. However, in R. imitator, aposematism is associated with four different color phenotypes that mimic a complex of congeneric species occurring across the mimic frog's geographic range. Investigations of the underlying mechanics of color production in these frogs can provide insights into how and why these different morphs evolved. We used histological samples to examine divergence in the color production mechanisms used by R. imitator to produce effective aposematic signals across its geographic range. We measured the coverage of melanophores and xanthophores (the area covered by chromatophores divided by total area of the skin section) in each color morph. We find that morphs that produce orange skin exhibit a higher coverage of xanthophores and lower coverage of melanophores than those that produce yellow skin. In turn, morphs that produce yellow skin exhibit a higher coverage of xanthophores and lower coverage of melanophores than those that produce green skin. Generally, across the morphs, a high ratio of xanthophores to melanophores is associated with colors of brighter spectral reflectance. Together, our results contribute to the understanding of color production in amphibians and document divergence in the histology of a species that is subject to divergent selection associated with aposematism.

Sex chromosomes as supergenes of speciation: why amphibians defy the rules?

As reflected by the two rules of speciation (Haldane's rule and the large X-/Z-effect), sex chromosomes are expected to behave like supergenes of speciation: they recombine only in one sex (XX females or ZZ males), supposedly recruit sexually antagonistic genes and evolve faster than autosomes, which can all contribute to pre-zygotic and post-zygotic isolation. While this has been mainly studied in organisms with conserved sex-determining systems and highly differentiated (heteromorphic) sex chromosomes like mammals, birds and some insects, these expectations are less clear in organismal groups where sex chromosomes repeatedly change and remain mostly homomorphic, like amphibians. In this article, we review the proposed roles of sex-linked genes in isolating nascent lineages throughout the speciation continuum and discuss their support in amphibians given current knowledge of sex chromosome evolution and speciation modes. Given their frequent recombination and lack of differentiation, we argue that amphibian sex chromosomes are not expected to become supergenes of speciation, which is reflected by the rarity of empirical studies consistent with a 'large sex chromosome effect' in frogs and toads. The diversity of sex chromosome systems in amphibians has a high potential to disentangle the evolutionary mechanisms responsible for the emergence of sex-linked speciation genes in other organisms. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

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.

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.

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.

Diversity in warning coloration: selective paradox or the norm?

Aposematic theory has historically predicted that predators should select for warning signals to converge on a single form, as a result of frequency-dependent learning. However, widespread variation in warning signals is observed across closely related species, populations and, most problematically for evolutionary biologists, among individuals in the same population. Recent research has yielded an increased awareness of this diversity, challenging the paradigm of signal monomorphy in aposematic animals. Here we provide a comprehensive synthesis of these disparate lines of investigation, identifying within them three broad classes of explanation for variation in aposematic warning signals: genetic mechanisms, differences among predators and predator behaviour, and alternative selection pressures upon the signal. The mechanisms producing warning coloration are also important. Detailed studies of the genetic basis of warning signals in some species, most notably Heliconius butterflies, are beginning to shed light on the genetic architecture facilitating or limiting key processes such as the evolution and maintenance of polymorphisms, hybridisation, and speciation. Work on predator behaviour is changing our perception of the predator community as a single homogenous selective agent, emphasising the dynamic nature of predator-prey interactions. Predator variability in a range of factors (e.g. perceptual abilities, tolerance to chemical defences, and individual motivation), suggests that the role of predators is more complicated than previously appreciated. With complex selection regimes at work, polytypisms and polymorphisms may even occur in Müllerian mimicry systems. Meanwhile, phenotypes are often multifunctional, and thus subject to additional biotic and abiotic selection pressures. Some of these selective pressures, primarily sexual selection and thermoregulation, have received considerable attention, while others, such as disease risk and parental effects, offer promising avenues to explore. As well as reviewing the existing evidence from both empirical studies and theoretical modelling, we highlight hypotheses that could benefit from further investigation in aposematic species. Finally by collating known instances of variation in warning signals, we provide a valuable resource for understanding the taxonomic spread of diversity in aposematic signalling and with which to direct future research. A greater appreciation of the extent of variation in aposematic species, and of the selective pressures and constraints which contribute to this once-paradoxical phenomenon, yields a new perspective for the field of aposematic signalling.

Testing the adaptive hypothesis of Batesian mimicry among hybridizing North American admiral butterflies

Batesian mimicry is characterized by phenotypic convergence between an unpalatable model and a palatable mimic. However, because convergent evolution may arise via alternative evolutionary mechanisms, putative examples of Batesian mimicry must be rigorously tested. Here, we used artificial butterfly facsimiles (N = 4000) to test the prediction that (1) palatable Limenitis lorquini butterflies should experience reduced predation when in sympatry with their putative model, Adelpha californica, (2) protection from predation on L. lorquini should erode outside of the geographical range of the model, and (3) mimetic color pattern traits are more variable in allopatry, consistent with relaxed selection for mimicry. We find support for these predictions, implying that this convergence is the result of selection for Batesian mimicry. Additionally, we conducted mark-recapture studies to examine the effect of mimicry and found that mimics survive significantly longer at sites where the model is abundant. Finally, in contrast to theoretical predictions, we found evidence that the Batesian model (A. californica) is protected from predation outside of its geographic range. We discuss these results considering the ongoing hybridization between L. lorquini and its sister species, L. weidemeyerii, and growing evidence that selection for mimicry predictably leads to a reduction in gene flow between nascent species.

Radiation of the polymorphic Little Devil poison frog (Oophaga sylvatica) in Ecuador

Some South American poison frogs (Dendrobatidae) are chemically defended and use bright aposematic colors to warn potential predators of their unpalatability. Aposematic signals are often frequency-dependent where individuals deviating from a local model are at a higher risk of predation. However, extreme diversity in the aposematic signal has been documented in poison frogs, especially in Oophaga. Here, we explore the phylogeographic pattern among color-divergent populations of the Little Devil poison frog Oophaga sylvatica by analyzing population structure and genetic differentiation to evaluate which processes could account for color diversity within and among populations. With a combination of PCR amplicons (three mitochondrial and three nuclear markers) and genome-wide markers from a double-digested RAD (ddRAD) approach, we characterized the phylogenetic and genetic structure of 199 individuals from 13 populations (12 monomorphic and 1 polymorphic) across the O. sylvatica distribution. Individuals segregated into two main lineages by their northern or southern latitudinal distribution. A high level of genetic and phenotypic polymorphism within the northern lineage suggests ongoing gene flow. In contrast, low levels of genetic differentiation were detected among the southern lineage populations and support recent range expansions from populations in the northern lineage. We propose that a combination of climatic gradients and structured landscapes might be promoting gene flow and phylogenetic diversification. Alternatively, we cannot rule out that the observed phenotypic and genomic variations are the result of genetic drift on near or neutral alleles in a small number of genes.

Polymorphism at a mimicry supergene maintained by opposing frequency-dependent selection pressures

Explaining the maintenance of adaptive diversity within populations is a long-standing goal in evolutionary biology, with important implications for conservation, medicine, and agriculture. Adaptation often leads to the fixation of beneficial alleles, and therefore it erodes local diversity so that understanding the coexistence of multiple adaptive phenotypes requires deciphering the ecological mechanisms that determine their respective benefits. Here, we show how antagonistic frequency-dependent selection (FDS), generated by natural and sexual selection acting on the same trait, maintains mimicry polymorphism in the toxic butterfly Heliconius numata Positive FDS imposed by predators on mimetic signals favors the fixation of the most abundant and best-protected wing-pattern morph, thereby limiting polymorphism. However, by using mate-choice experiments, we reveal disassortative mate preferences of the different wing-pattern morphs. The resulting negative FDS on wing-pattern alleles is consistent with the excess of heterozygote genotypes at the supergene locus controlling wing-pattern variation in natural populations of H. numata The combined effect of positive and negative FDS on visual signals is sufficient to maintain a diversity of morphs displaying accurate mimicry with other local prey, although some of the forms only provide moderate protection against predators. Our findings help understand how alternative adaptive phenotypes can be maintained within populations and emphasize the need to investigate interactions between selective pressures in other cases of puzzling adaptive polymorphism.

Poison frog color morphs express assortative mate preferences in allopatry but not sympatry

The concurrent divergence of mating traits and preferences is necessary for the evolution of reproductive isolation via sexual selection, and such coevolution has been demonstrated in diverse lineages. However, the extent to which assortative mate preferences are sufficient to drive reproductive isolation in nature is less clear. Natural contact zones between lineages divergent in traits and preferences provide exceptional opportunities for testing the predicted evolutionary consequences of such divergence. The strawberry poison frog (Oophaga pumilio) displays extreme color polymorphism in and around the young Bocas del Toro archipelago. In a transition zone between red and blue allopatric lineages, we asked whether female preferences diverged along with coloration, and whether any divergent preferences persist in a zone of sympatry. When choosing among red, blue and phenotypically intermediate males, females from monomorphic red and monomorphic blue populations both expressed assortative preferences. However, red, blue, and intermediate females from the contact zone all preferred red males, suggesting that divergent preferences may be insufficient to effect behavioral isolation. Our results highlight the complexity of behavioral isolation, and the need for studies that can reveal the circumstances under which divergent preferences do and do not contribute to speciation.

Mimetic Divergence and the Speciation Continuum in the Mimic Poison Frog Ranitomeya imitator

While divergent ecological adaptation can drive speciation, understanding the factors that facilitate or constrain this process remains a major goal in speciation research. Here, we study two mimetic transition zones in the poison frog Ranitomeya imitator, a species that has undergone a Müllerian mimetic radiation to establish four morphs in Peru. We find that mimetic morphs are strongly phenotypically differentiated, producing geographic clines with varying widths. However, distinct morphs show little neutral genetic divergence, and landscape genetic analyses implicate isolation by distance as the primary determinant of among-population genetic differentiation. Mate choice experiments suggest random mating at the transition zones, although certain allopatric populations show a preference for their own morph. We present evidence that this preference may be mediated by color pattern specifically. These results contrast with an earlier study of a third transition zone, in which a mimetic shift was associated with reproductive isolation. Overall, our results suggest that the three known mimetic transition zones in R. imitator reflect a speciation continuum, which we have characterized at the geographic, phenotypic, behavioral, and genetic levels. We discuss possible explanations for variable progress toward speciation, suggesting that multifarious selection on both mimetic color pattern and body size may be responsible for generating reproductive isolation.

Number of genes controlling a quantitative trait in a hybrid zone of the aposematic frog Ranitomeya imitator

The number of genes controlling mimetic traits has been a topic of much research and discussion. In this paper, we examine a mimetic, dendrobatid frog Ranitomeya imitator, which harbours extensive phenotypic variation with multiple mimetic morphs, not unlike the celebrated Heliconius system. However, the genetic basis for this polymorphism is unknown, and not easy to determine using standard experimental approaches, for this hard-to-breed species. To circumvent this problem, we first develop a new protocol for automatic quantification of complex colour pattern phenotypes from images. Using this method, which has the potential to be applied in many other systems, we define a phenotype associated with differences in colour pattern between different mimetic morphs. We then proceed to develop a maximum-likelihood method for estimating the number of genes affecting a quantitative trait segregating in a hybrid zone. This method takes advantage of estimates of admixture proportions obtained using genetic data, such as microsatellite markers, and is applicable to any other system where a phenotype has been quantified in an admixture/introgression zone. We evaluate the method using extensive simulations and apply it to the R. imitator system. We show that probably one or two, or at most three genes, control the mimetic phenotype segregating in a R. imitator hybrid zone identified using image analyses.