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

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.

Eco-Metabolomics Applied to the Chemical Ecology of Poison Frogs (Dendrobatoidea)

Amphibians are one of the most remarkable sources of unique natural products. Biogenic amines, peptides, bufodienolides, alkaloids, and volatile organic compounds have been characterized in different species. The superfamily Dendrobatoidea represents one of the most enigmatic cases of study in chemical ecology because their skin secretome is composed by a complex mixture (i.e. cocktail) of highly lethal and noxious unique alkaloid structures. While chemical defences from dendrobatoids (families Dendrobatidae and Aromobatidae) have been investigated employing ecological, behavioral, phylogenetic and evolutionary perspectives, studies about the analytical techniques needed to perform the chemical characterization have been neglected for many years. Therefore, our aim is to summarize the current methods applied for the characterization of chemical profiles in dendrobatoids and to illustrate innovative Eco-metabolomics strategies that could be translated to this study model. This approach could be extended to natural products other than alkaloids and implemented for the chemical analysis of different species of dendrobatoids employing both low- and high-resolution mass spectrometers. Here, we overview important biological features to be considered, procedures that could be applied to perform the chemical characterization, steps and tools to perform an Eco-metabolomic analysis, and a final discussion about future perspectives.

2b or not 2b? 2bRAD is an effective alternative to ddRAD for phylogenomics

Restriction-site-associated DNA sequencing (RADseq) has become an accessible way to obtain genome-wide data in the form of single-nucleotide polymorphisms (SNPs) for phylogenetic inference. Nonetheless, how differences in RADseq methods influence phylogenetic estimation is poorly understood because most comparisons have largely relied on conceptual predictions rather than empirical tests. We examine how differences in ddRAD and 2bRAD data influence phylogenetic estimation in two non-model frog groups. We compare the impact of method choice on phylogenetic information, missing data, and allelic dropout, considering different sequencing depths. Given that researchers must balance input (funding, time) with output (amount and quality of data), we also provide comparisons of laboratory effort, computational time, monetary costs, and the repeatability of library preparation and sequencing. Both 2bRAD and ddRAD methods estimated well-supported trees, even at low sequencing depths, and had comparable amounts of missing data, patterns of allelic dropout, and phylogenetic signal. Compared to ddRAD, 2bRAD produced more repeatable datasets, had simpler laboratory protocols, and had an overall faster bioinformatics assembly. However, many fewer parsimony-informative sites per SNP were obtained from 2bRAD data when using native pipelines, highlighting a need for further investigation into the effects of each pipeline on resulting datasets. Our study underscores the importance of comparing RADseq methods, such as expected results and theoretical performance using empirical datasets, before undertaking costly experiments.

Factors determining the dorsal coloration pattern of aposematic salamanders

Aposematic bright colors have a key role for animal defense and can be expressed through metabolic production or by acquiring pigments from diet. Aposematic coloration can be related to both local adaptations and availability of trophic resources. The European fire salamander (Salamandra salamandra) shows significant color variability and occurs across a broad range of habitats. Here we combined field observations with common rearing experiments to disentangle the role of environmental conditions and local adaptations in determining aposematic coloration of salamander populations. We assessed color variation and measured habitat features and food availability in adults from 25 populations. Furthermore, we reared newborn larvae from 10 populations under different food availability and analyzed color of metamorphs. To assess color pattern, we measured the percentage of yellow covering the body, and the Hue, Saturation and Value of yellow coloration. Adult showed strong variation of color pattern; variation was strongly related to the individual's size, to habitat productivity and to food availability. Under common garden conditions, differences between populations were not anymore evident, and coloration was only affected by resource availability during larval development. Our results suggest that environmental conditions and food availability are more important than local adaptations in determining differences in aposematic color pattern.

Genetic isolation by distance underlies colour pattern divergence in red-eyed treefrogs (Agalychnis callidryas)

Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of colour-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg colour patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple colour-pattern morphs. To explore possible causes of this variation, we generated a single nucleotide polymorphism data set to analyse population genetic structure, measure genetic diversity and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond to individual measures of colour pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that colour-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance indicates that colour-pattern variation was either retained through other processes such as ancestral colour polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on colour pattern at a small geographic scale.

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.

Structure and evolution of the sexually dimorphic integumentary swelling on the hands of dendrobatid poison frogs and their relatives (Amphibia: Anura: Dendrobatoidea)

The hands of adult males of many dendrobatid poison frogs and their relatives possess swelling formed by glandular tissue hypothesized to secrete courtship pheromones delivered to the female during cephalic amplexus. Variation in the occurrence and external morphology of the swelling of finger IV has provided important evidence for dendrobatoid systematics for decades, but its underlying structure has not been investigated. We undertook a detailed comparative analysis of the integument of the hand, including both external morphology and histology, of 36 species representing the diversity of dendrobatoid frogs and several close relatives. The swelling is caused by four densely packed, hypertrophic, morpho-histochemical types of specialized mucous glands (SMGs). We observed type I SMGs on fingers II-V and the wrist, including areas that are not swollen, types II and IV exclusively on finger IV, and type III on finger IV and the wrist. Type I SMGs occur either in isolation or together with types II, III, or IV; types II, III, and IV never occur together or without type I. We delimited 15 characters to account for the variation in external morphology and the occurrence of SMGs. Our data suggest that type I SMGs are a new synapomorphy for Dendrobatoidea and that type II SMGs originated in either the most recent common ancestor (MRCA) of Dendrobatidae or independently in the MRCAs of Aromobatidae and Colostethinae, respectively, while types III and IV are restricted to Anomaloglossus. The discovery of these SMGs adds a new dimension to studies of poison frog reproductive biology, which have investigated acoustic, visual, and tactile cues in courtship, mating, and parental care across the diversity of Dendrobatoidea for decades but have almost entirely overlooked the possible role of chemical cues.

First characterization of toxic alkaloids and volatile organic compounds (VOCs) in the cryptic dendrobatid Silverstoneia punctiventris

Background

Poison frogs are known for the outstanding diversity of alkaloid-based chemical defences with promising therapeutic applications. However, current knowledge about chemical defences in Dendrobatoidea superfamily has two sources of bias. First, cryptic, brown-colored species have been neglected in comparison to those conspicuously colored, and second, there has been little interest in characterizing metabolites other than alkaloids mediating defensive functions. In an effort to contribute to fill the gap of knowledge about cryptic species and broadening the spectrum of compounds analyzed we have applied head-space solid phase microextraction coupled to gas chromatography and mass spectrometry (HS-SPME/GC-MS) for extracting amphibian alkaloids and volatile organic compounds (VOCs) from Silverstoneia punctiventris.

Results

Using the skin from 8 specimens in 4 biological replicates we have found 33 different compounds. Twenty of them were classified as VOCs into 15 chemical classes including alkanes, alcohols, carbonyl compounds, methylpyridines, benzothiazoles, N-alkylpyrrolidines, pyrazines, and sesquiterpenoids, some of which were previously reported as repellents, defence compounds or defence pheromones in other organisms, and as sex pheromones in a treefrog. Interestingly, six of the remaining compounds were identified as alkaloids previously reported in other toxic/unpalatable dendrobatid frogs.

Conclusions

This is the first report of alkaloids and VOCs found in the Silverstoneia genus, which has been assumed for decades as non-chemically defended. This study establishes HS-SPME/GC-MS as a new application for a simultaneous approach to amphibian alkaloids and VOCs in poison frogs while opens up new research questions to assess the co-occurrence of both type of compounds and to investigate the evolutionary significance of a defence gradient that includes olfactory avoidance, unpalatability, and toxicity in dendrobatids. In addition, our results show that amphibian alkaloids could have a dual function (olfactory at distance, taste by contact) never explored before neither in Silverstonaeia nor in any other dendrobatid species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00420-1.

Is mimicry a diversification-driver in ants? Biogeography, ecology, ethology, genetics and morphology define a second West-Palaearctic Colobopsis species (Hymenoptera: Formicidae)

The West-Palaearctic Colobopsis ant populations have long been considered a single species (Colobopsis truncata). We studied the diversity of this species by employing a multidisciplinary approach and combining data from our surveys, museum and private collections, and citizen science platforms. As a result, we have revealed the existence of a second species, which we describe as Colobopsis imitans sp. nov., distributed allopatrically from Co. truncata and living in the Maghreb, Sicily and southern Iberia. While the pigmentation of Co. truncata is reminiscent of Dolichoderus quadripunctatus, that of Co. imitans is similar to Crematogaster scutellaris, with which Co. imitans lives in close spatial association, and whose foraging trails it habitually follows, similar to Camponotus lateralis and other ant-mimicking ants. The isolation between Co. imitans and Co. truncata seems to have occurred relatively recently because of significant, yet not extreme, morphometric differentiation, and to mtDNA polyphyly. Both Co. imitans and Co. truncata appear to employ mimicry of an unpalatable or aggressive ant species as an important defensive strategy; this ‘choice’ of a different model species is motivated by biogeographic reasons and appears to act as a critical evolutionary driver of their diversification.

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.

Integrating alpha, beta, and phylogenetic diversity to understand anuran fauna along environmental gradients of tropical forests in western Ecuador

The study of current distribution patterns of amphibian species in South America is of particular interest in areas such as evolutionary ecology and conservation biology. These patterns could be playing an important role in biological interactions, population size, and connectivity, and potential extinction risk in amphibians. Here, we tested the effects of spatial and environmental factors on the variation, turnover, and phylogenetic diversity of anuran amphibian species in tropical forests of western Ecuador. Data for presence/absence of 101 species of 34 genera and 10 families registered in 12 sites (nested in four biogeographic units) were obtained through fieldwork, museum collections, and literature records. We examined the influence of geographical, altitudinal, temperature, and precipitation distances on differences in anuran composition between sites. We found significant positive correlations among all of these variables with anuran distribution. The greatest alpha diversity (species richness) was found in the Equatorial Chocó biogeographic unit. Equatorial Pacific biogeographic unit could act as a transition zone between the Equatorial Chocó and Equatorial Tumbes. The western Andes (Western Cordillera biogeographic unit) was the most dissimilar and exhibited a higher species turnover rate than the other biogeographic units. Our results suggest that precipitation and elevation play a key role in maintaining the diversity of amphibian species in western Ecuador.

Transcriptomic Signatures of Experimental Alkaloid Consumption in a Poison Frog

In the anuran family Dendrobatidae, aposematic species obtain their toxic or unpalatable alkaloids from dietary sources, a process known as sequestering. To understand how toxicity evolved in this family, it is paramount to elucidate the pathways of alkaloid processing (absorption, metabolism, and sequestering). Here, we used an exploratory skin gene expression experiment in which captive-bred dendrobatids were fed alkaloids. Most of these experiments were performed with Dendrobates tinctorius, but some trials were performed with D. auratus, D. leucomelas and Allobates femoralis to explore whether other dendrobatids would show similar patterns of gene expression. We found a consistent pattern of up-regulation of genes related to muscle and mitochondrial processes, probably due to the lack of mutations related to alkaloid resistance in these species. Considering conserved pathways of drug metabolism in vertebrates, we hypothesize alkaloid degradation is a physiological mechanism of resistance, which was evidenced by a strong upregulation of the immune system in D. tinctorius, and of complement C2 across the four species sampled. Probably related to this strong immune response, we found several skin keratins downregulated, which might be linked to a reduction of the cornified layer of the epidermis. Although not conclusive, our results offer candidate genes and testable hypotheses to elucidate alkaloid processing in poison frogs.

Isolation and characterization of 12 microsatellite loci in Epipedobates anthonyi (Amphibia: Anura: Dendrobatidae) for population genetic analysis

Tropical anurans are among the most diverse and vulnerable organisms on Earth, yet the evolutionary mechanisms behind their diversity remain relatively unexplored. Epipedobates anthonyi is a poison frog that inhabits southern Ecuador and northern Peru along a broad elevational range (0-1800 m). Throughout its range, this species exhibits variation in phenotypic traits, such as color, advertisement calls, and alkaloid composition. The aim of this study is to isolate and characterize microsatellite loci to investigate patterns of genetic variation within the species. Using a next-generation sequencing approach to screen an enriched genomic library, we report twelve polymorphic microsatellite loci. The number of alleles per locus ranged from 7 to 15 per population. For the two populations tested, mean observed heterozygosity was 0.69 and 0.79, and mean expected heterozigosity was 0.84 and 0.85 respectively. Only locus EAN002 showed significant departure of HWE in both populations. None of the loci showed consistent null alleles in both populations. Also, no evidence of linkage disequilibrium was found across loci. In this paper, we report for the first time 12 microsatellite loci for E. anthonyi. These markers will be used to further elucidate evolutionary mechanisms underlying genetic and phenotypic variation across the species' range.

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.

Uncovering hidden specific diversity of Andean glassfrogs of the Centrolene buckleyi species complex (Anura: Centrolenidae)

The glassfrog Centrolene buckleyi has been recognized as a species complex. Herein, using coalescence-based species delimitation methods, we evaluate the specific diversity within this taxon. Four coalescence approaches (generalized mixed Yule coalescents, Bayesian general mixed Yule-coalescent, Poisson tree processes, and Bayesian Poisson tree processes) were consistent with the delimitation results, identifying four lineages within what is currently recognized as C. buckleyi. We propose three new candidate species that should be tested with nuclear markers, morphological, and behavioral data. In the meantime, for conservation purposes, candidate species should be considered evolutionary significant units, in light of observed population crashes in the C. buckleyi species complex. Finally, our results support the validity of C. venezuelense, formerly considered as a subspecies of C. buckleyi.

Biased generalization of salient traits drives the evolution of warning signals

The importance of receiver biases in shaping the evolution of many signalling systems is widely acknowledged. Here, we show that receiver bias can explain which traits evolve to become warning signals. For warning coloration, a generalization bias for a signalling trait can result from predators learning to discriminate unprofitable from profitable prey. However, because the colour patterns of prey are complex traits with multiple components, it is crucial to understand which of the many aspects of prey appearance evolve into signals. We provide experimental evidence that the more salient differences in prey traits give rise to greater generalization bias, corresponding to stronger selection towards trait exaggeration. Our results are based on experiments with domestic chickens as predators in a Skinner-box-like setting, and imply that the difference in appearance between profitable and unprofitable prey that is most rapidly learnt produces the greatest generalization bias. As a consequence, certain salient traits of unprofitable prey are selected towards exaggeration to even higher salience, driving the evolution of warning coloration. This general idea may also help to explain the evolution of many other striking signalling traits found in nature.

Diversity within diversity: Parasite species richness in poison frogs assessed by transcriptomics

Symbionts (e.g., endoparasites and commensals) play an integral role in their host's ecology, yet in many cases their diversity is likely underestimated. Although endoparasites are traditionally characterized using morphology, sequences of conserved genes, and shotgun metagenomics, host transcriptomes constitute an underused resource to identify these organisms' diversity. By isolating non-host transcripts from host transcriptomes, individual host tissues can now simultaneously reveal their endoparasite species richness (i.e., number of different taxa) and provide insights into parasite gene expression. These approaches can be used in host taxa whose endoparasites are mostly unknown, such as those of tropical amphibians. Here, we focus on the poison frogs (Dendrobatidae) as hosts, which are a Neotropical clade known for their bright coloration and defensive alkaloids. These toxins are an effective protection against vertebrate predators (e.g., snakes and birds), bacteria, and skin-biting ectoparasites (e.g., mosquitoes); however, little is known about their deterrence against eukaryotic endoparasites. With de novo transcriptomes of dendrobatids, we developed a bioinformatics pipeline for endoparasite identification that uses host annotated RNA-seq data and set of a priori parasite taxonomic terms, which are used to mine for specific endoparasites. We found a large community of helminths and protozoans that were mostly restricted to the digestive tract and a few systemic parasites (e.g., Trypanosoma). Contrary to our expectations, all dendrobatid frogs regardless of the presence of alkaloid defenses have endoparasites, with their highest species richness located in the frog digestive tract. Some of these organisms (e.g., roundworms) might prove to be generalists, as they were not found to be co-diversifying with their frog hosts. We propose that endoparasites may escape poison frogs' chemical defenses by colonizing tissues with fewer alkaloids than the frog's skin, where most toxins are stored.

Interacting amino acid replacements allow poison frogs to evolve epibatidine resistance

Animals that wield toxins face self-intoxication. Poison frogs have a diverse arsenal of defensive alkaloids that target the nervous system. Among them is epibatidine, a nicotinic acetylcholine receptor (nAChR) agonist that is lethal at microgram doses. Epibatidine shares a highly conserved binding site with acetylcholine, making it difficult to evolve resistance yet maintain nAChR function. Electrophysiological assays of human and frog nAChR revealed that one amino acid replacement, which evolved three times in poison frogs, decreased epibatidine sensitivity but at a cost of acetylcholine sensitivity. However, receptor functionality was rescued by additional amino acid replacements that differed among poison frog lineages. Our results demonstrate how resistance to agonist toxins can evolve and that such genetic changes propel organisms toward an adaptive peak of chemical defense.