Hybridization promotes color polymorphism in the aposematic harlequin poison frog, Oophaga histrionica

Deep divergences among inconspicuously colored clades of Epipedobates poison frogs

Poison frogs (Dendrobatidae) are famous for their aposematic species, having a combination of diverse color patterns and defensive skin toxins, yet most species in this family are inconspicuously colored and considered non-aposematic. Epipedobates is among the youngest genus-level clades of Dendrobatidae that includes both aposematic and inconspicuous species. Using Sanger-sequenced mitochondrial and nuclear markers, we demonstrate deep genetic divergences among inconspicuous species of Epipedobates but relatively shallow genetic divergences among conspicuous species. Our phylogenetic analysis includes broad geographic sampling of the inconspicuous lineages typically identified as E. boulengeri and E. espinosai, which reveals two putative new species, one in west-central Colombia (E. sp. 1) and the other in north-central Ecuador (E. aff. espinosai). We conclude that E. darwinwallacei is a junior subjective synonym of E. espinosai. We also clarify the geographic distributions of inconspicuous Epipedobates species including the widespread E. boulengeri. We provide a qualitative assessment of the phenotypic diversity in each nominal species, with a focus on the color and pattern of inconspicuous species. We conclude that Epipedobates contains eight known valid species, six of which are inconspicuous. A relaxed molecular clock analysis suggests that the most recent common ancestor of Epipedobates is ∼11.1 million years old, which nearly doubles previous estimates. Last, genetic information points to a center of species diversity in the Chocó at the southwestern border of Colombia with Ecuador. A Spanish translation of this text is available in the supplementary materials.

Does testosterone underly the interplay between male traits and territorial behavior in neotropical poison frogs?

The ability of individual animals to defend a territory as well as various phenotypic and behavioral traits may be targets of sexual selection used by males to evaluate their competitors or by females to choose males. A frequent question in animal behavior is whether male traits and characteristics of their territory are correlated and what are the mechanisms that may mediate such associations when they exist. Because hormones link phenotype to behavior, by studying the role of testosterone in territoriality one may come closer to understanding the mechanisms mediating correlations or lack thereof between characteristics of territories and of males. We evaluated whether variation in characteristics of territories (size and quality) are correlated with variation in morphology, coloration, testosterone, heterozygosity, and calls in two species of poison frogs. The Amazonian frog Allobates aff. trilineatus exhibits male care and defends territories only during the breeding season, while the endangered frog Oophaga lehmanni displays maternal care and defends territories throughout the year. We found that morphological traits (body length, weight, thigh size), call activity, and testosterone levels correlated with size and various indicators of quality of the territory. However, the direction of these correlations (whether positive or negative) and which specific morphological, acoustic traits or testosterone level variables covaried depended on the species. Our findings highlight an endocrine pathway as part of the physiological machinery that may underlie the interplay between male traits and territorial behavior. We were able to identify some male traits related to territory attributes, but whether females choose males based on these traits requires further research.

recolorize: An R package for flexible colour segmentation of biological images

Colour pattern variation provides biological information in fields ranging from disease ecology to speciation dynamics. Comparing colour pattern geometries across images requires colour segmentation, where pixels in an image are assigned to one of a set of colour classes shared by all images. Manual methods for colour segmentation are slow and subjective, while automated methods can struggle with high technical variation in aggregate image sets. We present recolorize, an R package toolbox for human-subjective colour segmentation with functions for batch-processing low-variation image sets and additional tools for handling images from diverse (high-variation) sources. The package also includes export options for a variety of formats and colour analysis packages. This paper illustrates recolorize for three example datasets, including high variation, batch processing and combining with reflectance spectra, and demonstrates the downstream use of methods that rely on this output.

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.

Unexpected colour pattern variation in mimetic frogs: implication for the diversification of warning signals in the genus Ranitomeya

Predation is expected to promote uniformity in the warning coloration of defended prey, but also mimicry convergence between aposematic species. Despite selection constraining both colour-pattern and population divergence, many aposematic animals display numerous geographically structured populations with distinct warning signal. Here, we explore the extent of phenotypic variation of sympatric species of Ranitomeya poison frogs and test for theoretical expectations on variation and convergence in mimetic signals. We demonstrate that both warning signal and mimetic convergence are highly variable and are negatively correlated: some localities display high variability and no mimicry while in others the phenotype is fixed and mimicry is perfect. Moreover, variation in warning signals is always present within localities, and in many cases this variation overlaps between populations, such that variation is continuous. Finally, we show that coloration is consistently the least variable element and is likely of greater importance for predator avoidance compared to patterning. We discuss the implications of our results in the context of warning signal diversification and suggest that, like many other locally adapted traits, a combination of standing genetic variation and founding effect might be sufficient to enable divergence in colour pattern.

Skin transcriptional profiles in Oophaga poison frogs

Aposematic organisms advertise their defensive toxins to predators using a variety of warning signals, including bright coloration. While most Neotropical poison frogs (Dendrobatidae) rely on crypsis to avoid predators, Oophaga poison frogs from South America advertise their chemical defenses, a complex mix of diet-derived alkaloids, by using conspicuous hues. The present study aimed to characterize the skin transcriptomic profiles of South American Oophaga poison frogs. Our analyses showed very similar transcriptomic profiles for these closely related species in terms of functional annotation and relative abundance of gene ontology terms expressed. Analyses of expression profiles of Oophaga and available skin transcriptomes of cryptic anurans allowed us to propose initial hypotheses for the active sequestration of alkaloid-based chemical defenses and to highlight some genes that may be potentially involved in resistance mechanisms to avoid self-intoxication and skin coloration. In doing so, we provide an important molecular resource for the study of warning signals that will facilitate the assembly and annotation of future poison frog genomes.

Genetic and phenotypic evidence of a contact zone between divergent colour morphs of the iconic red-eyed treefrog

Hybrid zones act as natural laboratories where divergent genomes interact, providing powerful systems for examining the evolutionary processes underlying biological diversity. In this study, we characterized patterns of genomic and phenotypic variation resulting from hybridization between divergent intraspecific lineages of the Neotropical red-eyed treefrog (Agalychnis callidryas). We found genetic evidence of a newly discovered contact zone and phenotypic novelty in leg colour-a trait suspected to play a role in mediating assortative mating in this species. Analysis of hybrid ancestry revealed an abundance of later-generation F_n individuals, suggesting persistence of hybrids in the contact zone. Hybrids are predominantly of southern ancestry but are phenotypically more similar to northern populations. Genome-wide association mapping revealed QTL with measurable effects on leg-colour variation, but further work is required to dissect the architecture of this trait and establish causal links. Further, genomic cline analyses indicated substantial variation in patterns of introgression across the genome. Directional introgression of loci associated with different aspects of leg colour are inherited from each parental lineage, creating a distinct hybrid colour pattern. We show that hybridization can generate new phenotypes, revealing the evolutionary processes that potentially underlie patterns of phenotypic diversity in this iconic polytypic frog. Our study is consistent with a role of hybridization and sexual selection in lineage diversification, evolutionary processes that have been implicated in accelerating divergence in the most phenotypically diverse species.

The Influence of Environmental Variation on the Genetic Structure of a Poison Frog Distributed Across Continuous Amazonian Rainforest

Biogeographic barriers such as rivers have been shown to shape spatial patterns of biodiversity in the Amazon basin, yet relatively little is known about the distribution of genetic variation across continuous rainforest. Here, we characterize the genetic structure of the brilliant-thighed poison frog (Allobates femoralis) across an 880-km-long transect along the Purus-Madeira interfluve south of the Amazon river, based on 64 individuals genotyped at 7609 single-nucleotide polymorphism (SNP) loci. A population tree and clustering analyses revealed 4 distinct genetic groups, one of which was strongly divergent. These genetic groups were concomitant with femoral spot coloration differences, which was intermediate within a zone of admixture between two of the groups. The location of these genetic groups did not consistently correspond to current ecological transitions between major forest types. A multimodel approach to quantify the relative influence of isolation-by-geographic distance (IBD) and isolation-by-environmental resistance (IBR) nevertheless revealed that, in addition to a strong signal of IBD, spatial genetic differentiation was explained by IBR primarily linked to dry season intensity (r2 = 8.4%) and canopy cover (r2 = 6.4%). We show significant phylogenetic divergence in the absence of obvious biogeographical barriers and that finer-scaled measures of genetic structure are associated with environmental variables also known to predict the density of A. femoralis.

The evolution of polymorphism in the warning coloration of the Amazonian poison frog Adelphobates galactonotus

While intraspecific variation in aposematic signals can be selected for by different predatory responses, their evolution is also contingent on other processes shaping genetic variation. We evaluate the relative contributions of selection, geographic isolation, and random genetic drift to the evolution of aposematic color polymorphism in the poison frog Adelphobates galactonotus, distributed throughout eastern Brazilian Amazonia. Dorsal coloration was measured for 111 individuals and genetic data were obtained from 220 individuals at two mitochondrial genes (mtDNA) and 7963 Single Nucleotide Polymorphisms (SNPs). Four color categories were described (brown, blue, yellow, orange) and our models of frog and bird visual systems indicated that each color was distinguishable for these taxa. Using outlier and correlative analyses we found no compelling genetic evidence for color being under divergent selection. A time-calibrated mtDNA tree suggests that the present distribution of dorsal coloration resulted from processes occurring during the Pleistocene. Separate phylogenies based on SNPs and mtDNA resolved the same well supported clades, each containing different colored populations. Ancestral character state analysis provided some evidence for evolutionary transitions in color type. Genetic structure was more strongly associated with geographic features, than color category, suggesting that the distribution of color is explained by localized processes. Evidence for geographic isolation together with estimates of low effective population size implicates drift as playing a key role in color diversification. Our results highlight the relevance of considering the neutral processes involved with the evolution of traits with important fitness consequences.

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.

Phenotypic and genetic divergence in reed frogs across a mosaic hybrid zone on São Tomé Island

Although naturally heterogeneous environments can lead to mosaic hybrid zones, human-induced habitat fragmentation can also lead to environmental heterogeneity and hybridization. Here we quantify phenotypic and molecular divergence across a reed frog mosaic hybrid zone on São Tomé Island as a first step towards understanding the consequences of hybridization across this heterogeneous landscape. The São Tomé giant reed frog (Hyperolius thomensis) is strongly tied to cool, wet, forest habitats whereas the distribution of Moller’s reed frog (H. molleri) spans cool, wet, forests to warm, dry, disturbed habitats. Correspondingly, hybridization is concentrated in the more forested, cool, wet sites relative to non-forested, warmer, drier habitats. Four of six sites with hybrid frogs are artificial water bodies near the forest edge, indicating that both breeding habitat and broader scale environmental variation are probably important for understanding interspecific interactions and the extent of hybridization in this system. Phenotypic variation (body size and ventral coloration) largely tracks genetic and environmental variation across the hybrid zone with larger and more pigmented frogs occurring in forested, cool, wet habitats. Understanding whether human-induced changes in habitat break down reproductive barriers will be essential for conservation management of the less abundant, forest-associated H. thomensis in the face of rampant hybridization.

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.

Origin and cross-century dynamics of an avian hybrid zone

Background

Characterizations of the dynamics of hybrid zones in space and time can give insights about traits and processes important in population divergence and speciation. We characterized a hybrid zone between tanagers in the genus Ramphocelus (Aves, Thraupidae) located in southwestern Colombia. We evaluated whether this hybrid zone originated as a result of secondary contact or of primary differentiation, and described its dynamics across time using spatial analyses of molecular, morphological, and coloration data in combination with paleodistribution modeling.

Results

Models of potential historical distributions based on climatic data and genetic signatures of demographic expansion suggested that the hybrid zone likely originated following secondary contact between populations that expanded their ranges out of isolated areas in the Quaternary. Concordant patterns of variation in phenotypic characters across the hybrid zone and its narrow extent are suggestive of a tension zone, maintained by a balance between dispersal and selection against hybrids. Estimates of phenotypic cline parameters obtained using specimens collected over nearly a century revealed that, in recent decades, the zone appears to have moved to the east and to higher elevations, and may have become narrower. Genetic variation was not clearly structured along the hybrid zone, but comparisons between historical and contemporary specimens suggested that temporal changes in its genetic makeup may also have occurred.

Conclusions

Our data suggest that the hybrid zone likey resulted from secondary contact between populations. The observed changes in the hybrid zone may be a result of sexual selection, asymmetric gene flow, or environmental change.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12862-017-1096-7) contains supplementary material, which is available to authorized users.

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.

Shifting barriers and phenotypic diversification by hybridisation

The establishment of hybrid taxa relies on reproductive isolation from the parental forms, typically achieved by ecological differentiation. Here, we present an alternative mechanism, in which shifts in the strength and location of dispersal barriers facilitate diversification by hybridisation. Our case study concerns the highly diverse, stenotopic rock-dwelling cichlids of the African Great Lakes, many of which display geographic colour pattern variation. The littoral habitat of these fish has repeatedly been restructured in the course of ancient lake level fluctuations. Genetic data and an experimental cross support the hybrid origin of a distinct yellow-coloured variant of Tropheus moorii from ancient admixture between two allopatric, red and bluish variants. Deficient assortative mating preferences imply that reproductive isolation continues to be contingent on geographic separation. Linking paleolimnological data with the establishment of the hybrid variant, we sketch a selectively neutral diversification process governed solely by rearrangements of dispersal barriers.

Diversity in warning coloration is easily recognized by avian predators

Warning coloration is a widespread strategy to alert predators about prey unprofitability. The success of this strategy partly depends on predators being able to learn and recognize certain signals as indicators of toxicity, and theory predicts that this is easier if signals converge on similar colours. However, the diversity in warning signal form is astonishing, contradicting predictions. Here, we quantified ladybird signal diversity with respect to avian vision, measuring how unique and discernible each signal is from one another. In addition, we measured signal conspicuousness against a series of backgrounds, namely an average green, average brown, and where we collected each species, to determine whether signals are more contrasting against the ladybirds' local substrates than compared to average ones. This allowed us to establish whether there are local adaptations in conspicuousness that promote signal diversity. We found that while ladybird signals are unique and recognizable, specialist species are more contrasting against the background they are most commonly found on. However, overall our study suggests that warning signals have evolved to be effective against a wide range of natural backgrounds, partly explaining the success of this strategy in nature.

The birth of aposematism: High phenotypic divergence and low genetic diversity in a young clade of poison frogs

Rapid radiation coupled with low genetic divergence often hinders species delimitation and phylogeny estimation even if putative species are phenotypically distinct. Some aposematic species, such as poison frogs (Dendrobatidae), have high levels of intraspecific color polymorphism, which can lead to overestimation of species when phenotypic divergence primarily guides species delimitation. We explored this possibility in the youngest origin of aposematism (3-7 MYA) in poison frogs, Epipedobates, by comparing genetic divergence with color and acoustic divergence. We found low genetic divergence (2.6% in the 16S gene) despite substantial differences in color and acoustic signals. While chemical defense is inferred to have evolved in the ancestor of Epipedobates, aposematic coloration evolved at least twice or was lost once in Epipedobates, suggesting that it is evolutionarily labile. We inferred at least one event of introgression between two cryptically colored species with adjacent ranges (E. boulengeri and E. machalilla). We also find evidence for peripheral isolation resulting in phenotypic divergence and potential speciation of the aposematic E. tricolor from the non-aposematic E. machalilla. However, we were unable to estimate a well-supported species tree or delimit species using multispecies coalescent models. We attribute this failure to factors associated with recent speciation including mitochondrial introgression, incomplete lineage sorting, and too few informative molecular characters. We suggest that species delimitation within young aposematic lineages such as Epipedobates will require genome-level molecular studies. We caution against relying solely on DNA barcoding for species delimitation or identification and highlight the value of phenotypic divergence and natural history in delimiting species.

Recovering the Genetic Identity of an Extinct-in-the-Wild Species: The Puzzling Case of the Alagoas Curassow

The conservation of many endangered taxa relies on hybrid identification, and when hybrids become morphologically indistinguishable from the parental species, the use of molecular markers can assign individual admixture levels. Here, we present the puzzling case of the extinct in the wild Alagoas Curassow (Pauxi mitu), whose captive population descends from only three individuals. Hybridization with the Razor-billed Curassow (P. tuberosa) began more than eight generations ago, and admixture uncertainty affects the whole population. We applied an analysis framework that combined morphological diagnostic traits, Bayesian clustering analyses using 14 microsatellite loci, and mtDNA haplotypes to assess the ancestry of all individuals that were alive from 2008 to 2012. Simulated data revealed that our microsatellites could accurately assign an individual a hybrid origin until the second backcross generation, which permitted us to identify a pure group among the older, but still reproductive animals. No wild species has ever survived such a severe bottleneck, followed by hybridization, and studying the recovery capability of the selected pure Alagoas Curassow group might provide valuable insights into biological conservation theory.

Stream noise, hybridization, and uncoupled evolution of call traits in two lineages of poison frogs: Oophaga histrionica and Oophaga lehmanni

According to the acoustic adaptation hypothesis, communication signals are evolutionary shaped in a way that minimizes its degradation and maximizes its contrast against the background noise. To compare the importance for call divergence of acoustic adaptation and hybridization, an evolutionary force allegedly promoting phenotypic variation, we compared the mate recognition signal of two species of poison frogs (Oophaga histrionica and O. lehmanni) at five localities: two (one per species) alongside noisy streams, two away from streams, and one interspecific hybrid. We recorded the calls of 47 males and characterized the microgeographic variation in their spectral and temporal features, measuring ambient noise level, body size, and body temperature as covariates. As predicted, frogs living in noisy habitats uttered high frequency calls and, in one species, were much smaller in size. These results support a previously unconsidered role of noise on streams as a selective force promoting an increase in call frequency and pleiotropic effects in body size. Regarding hybrid frogs, their calls overlapped in the signal space with the calls of one of the parental lineages. Our data support acoustic adaptation following two evolutionary routes but do not support the presumed role of hybridization in promoting phenotypic diversity.