Mechanisms for Color Convergence in a Mimetic Radiation of Poison Frogs

Expression patterns of melanin-related genes are linked to crypsis and conspicuousness in a pumpkin toadlet

Colour signals play pivotal roles in different communication systems, and the evolution of these characters has been associated with behavioural ecology, integumentary production processes and perceptual mechanisms of the species involved. Here, we present the first insight into the molecular and histological basis of skin colour polymorphism within a miniaturized species of pumpkin toadlet, potentially representing the lowest size threshold for colour polytypism in tetrapods. Brachycephalus actaeus exhibits a coloration ranging from cryptic green to conspicuous orange skin, and our findings suggest that colour morphs differ in their capability to be detected by potential predators. We also found that the distribution and abundance of chromatophores are variable in the different colour morphs. The expression pattern of coloration related genes was predominantly associated with melanin synthesis (including dct, edn1, mlana, oca2, pmel, slc24a5, tyrp1 and wnt9a). Up-regulation of melanin genes in grey, green and brown skin was associated with higher melanophore abundance than in orange skin, where xanthophores predominate. Our findings provide a significant foundation for comparing and understanding the diverse pathways that contribute to the evolution of pigment production in the skin of amphibians.

What Makes a Mimic? Orange, Red, and Black Color Production in the Mimic Poison Frog (Ranitomeya imitator)

Aposematic organisms rely on their conspicuous appearance to signal that they are defended and unpalatable. Such phenotypes are strongly tied to survival and reproduction. Aposematic colors and patterns are highly variable; however, the genetic, biochemical, and physiological mechanisms producing this conspicuous coloration remain largely unidentified. Here, we identify genes potentially affecting color variation in two color morphs of Ranitomeya imitator: the orange-banded Sauce and the redheaded Varadero morphs. We examine gene expression in black and orange skin patches from the Sauce morph and black and red skin patches from the Varadero morph. We identified genes differentially expressed between skin patches, including those that are involved in melanin synthesis (e.g. mlana, pmel, tyrp1), iridophore development (e.g. paics, ppat, ak1), pteridine synthesis (e.g. gch1, pax3-a, xdh), and carotenoid metabolism (e.g. dgat2, rbp1, scarb2). In addition, using weighted correlation network analysis, we identified the top 50 genes with high connectivity from the most significant network associated with gene expression differences between color morphs. Of these 50 genes, 13 were known to be related to color production (gch1, gmps, gpr143, impdh1, mc1r, pax3-a, pax7, ppat, rab27a, rlbp1, tfec, trpm1, xdh).

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which species share a conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which mimicry is produced in poison frogs. We assembled a 6.02-Gbp genome with a contig N50 of 310 Kbp, a scaffold N50 of 390 Kbp and 85% of expected tetrapod genes. We leveraged this genome to conduct gene expression analyses throughout development of four colour morphs of Ranitomeya imitator and two colour morphs from both R. fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes differentially expressed throughout development, many of them related to melanophores/melanin, iridophore development and guanine synthesis. We also identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species, and identify several plausible candidates for colouration in vertebrates (e.g. cd36, ep-cadherin and perlwapin). Finally, we hypothesise that keratin genes (e.g. krt8) are important for producing different structural colours within these frogs.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomics resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomics resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, antipredator strategies, and resilience and adaptive responses. They also serve as critical models for understanding widespread genomic characteristics, including evolutionary genome expansions and contractions given they have the largest range in genome sizes of any animal taxon and multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The advent of long-read sequencing technologies, along with computational techniques that enhance scaffolding capabilities and streamline computational workload is now enabling the ability to overcome some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC) in early 2023. This burgeoning community already has more than 282 members from 41 countries (6 in Africa, 131 in the Americas, 27 in Asia, 29 in Australasia, and 89 in Europe). The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and outline how the AGC can enable amphibian genomics research to "leap" to the next level.

Pervasive mimicry in flight behavior among aposematic butterflies

Flight was a key innovation in the adaptive radiation of insects. However, it is a complex trait influenced by a large number of interacting biotic and abiotic factors, making it difficult to unravel the evolutionary drivers. We investigate flight patterns in neotropical heliconiine butterflies, well known for mimicry of their aposematic wing color patterns. We quantify the flight patterns (wing beat frequency and wing angles) of 351 individuals representing 29 heliconiine and 9 ithomiine species belonging to ten color pattern mimicry groupings. For wing beat frequency and up wing angles, we show that heliconiine species group by color pattern mimicry affiliation. Convergence of down wing angles to mimicry groupings is less pronounced, indicating that distinct components of flight are under different selection pressures and constraints. The flight characteristics of the Tiger mimicry group are particularly divergent due to convergence with distantly related ithomiine species. Predator-driven selection for mimicry also explained variation in flight among subspecies, indicating that this convergence can occur over relatively short evolutionary timescales. Our results suggest that the flight convergence is driven by aposematic signaling rather than shared habitat between comimics. We demonstrate that behavioral mimicry can occur between lineages that have separated over evolutionary timescales ranging from <0.5 to 70 My.

Multiple Routes to Color Convergence in a Radiation of Neotropical Poison Frogs

Convergent evolution is defined as the independent evolution of similar phenotypes in different lineages. Its existence underscores the importance of external selection pressures in evolutionary history, revealing how functionally similar adaptations can evolve in response to persistent ecological challenges through a diversity of evolutionary routes. However, many examples of convergence, particularly among closely related species, involve parallel changes in the same genes or developmental pathways, raising the possibility that homology at deeper mechanistic levels is an important facilitator of phenotypic convergence. Using the genus Ranitomeya, a young, color-diverse radiation of Neotropical poison frogs, we set out to 1) provide a phylogenetic framework for this group, 2) leverage this framework to determine if color phenotypes are convergent, and 3) to characterize the underlying coloration mechanisms to test whether color convergence occurred through the same or different physical mechanisms. We generated a phylogeny for Ranitomeya using ultraconserved elements and investigated the physical mechanisms underlying bright coloration, focusing on skin pigments. Using phylogenetic comparative methods, we identified several instances of color convergence, involving several gains and losses of carotenoid and pterin pigments. We also found a compelling example of nonparallel convergence, where, in one lineage, red coloration evolved through the red pterin pigment drosopterin, and in another lineage through red ketocarotenoids. Additionally, in another lineage, "reddish" coloration evolved predominantly through structural color mechanisms. Our study demonstrates that, even within a radiation of closely related species, convergent evolution can occur through both parallel and nonparallel mechanisms, challenging the assumption that similar phenotypes among close relatives evolve through the same mechanisms.

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.

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.

Phylogenomic analysis of evolutionary relationships in Ranitomeya poison frogs (Family Dendrobatidae) using ultraconserved elements

The use of genome-scale data in phylogenetics has enabled recent strides in determining the relationships between taxa that are taxonomically problematic because of extensive morphological variation. Here, we employ a phylogenomic approach to infer evolutionary relationships within Ranitomeya (Anura: Dendrobatidae), an Amazonian lineage of poison frogs consisting of 16 species with remarkable diversity in color pattern, range size, and parental care behavior. We infer phylogenies with all described species of Ranitomeya from ultraconserved nuclear genomic elements (UCEs) and also estimate divergence times. Our results differ from previous analyses regarding interspecific relationships. Notably, we find that R. toraro and R. defleri are not sister species but rather distantly related, contrary to previous analyses based on smaller genetic datasets. We recover R. uakarii as paraphyletic, designate certain populations formerly assigned to R. fantastica from Peru as R. summersi, and transfer the French Guianan and eastern Brazilian R. amazonica populations to R. variabilis. By clarifying both inter- and intraspecific relationships within Ranitomeya, our study paves the way for future tests of hypotheses on color pattern evolution and historical biogeography.

Environmental gradients predict the ratio of environmentally acquired carotenoids to self-synthesised pteridine pigments

Carotenoids are important pigments producing integument colouration; however, their dietary availability may be limited in some environments. Many species produce yellow to red hues using a combination of carotenoids and self-synthesised pteridine pigments. A compelling hypothesis is that pteridines replace carotenoids in environments where carotenoid availability is limited. To test this hypothesis, we quantified concentrations of five carotenoid and six pteridine pigments in multiple skin colours and individuals from 27 species of agamid lizards. We show that environmental gradients predict the ratio of carotenoids to pteridines; carotenoid concentrations are lower and pteridine concentrations higher in arid environments with low vegetation productivity. Both carotenoid and pteridine pigments were present in all species, but only pteridine concentrations explained colour variation among species and there were no correlations between carotenoid and pteridine pigments with a similar hue. These results suggest that in arid environments, where carotenoids are likely limited, species may compensate by synthesising more pteridines but do not necessarily replace carotenoids with pteridines of similar hue.

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which multiple species share the same conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which vivid colour and pattern are produced in a Müllerian mimicry complex of poison frogs. We did this by first assembling a high-quality de novo genome assembly for the mimic poison frog Ranitomeya imitator. This assembled genome is 6.8 Gbp in size, with a contig N50 of 300 Kbp R. imitator and two colour morphs from both Ranitomeya fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes that are differentially expressed throughout development, many of them related to melanocyte development, melanin synthesis, iridophore development and guanine synthesis. Polytypic differences within species may be the result of differences in expression and/or timing of expression, whereas convergence for colour pattern between species does not appear to be due to the same changes in gene expression. In addition, we identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species. Finally, we hypothesize that genes in the keratin family are important for producing different structural colours within these frogs.

Transcriptomic Analysis of Skin Color in Anole Lizards

Color and color pattern are critical for animal camouflage, reproduction, and defense. Few studies, however, have attempted to identify candidate genes for color and color pattern in squamate reptiles, a colorful group with over 10,000 species. We used comparative transcriptomic analyses between white, orange, and yellow skin in a color-polymorphic species of anole lizard to 1) identify candidate color and color-pattern genes in squamates and 2) assess if squamates share an underlying genetic basis for color and color pattern variation with other vertebrates. Squamates have three types of chromatophores that determine color pattern: guanine-filled iridophores, carotenoid- or pteridine-filled xanthophores/erythrophores, and melanin-filled melanophores. We identified 13 best candidate squamate color and color-pattern genes shared with other vertebrates: six genes linked to pigment synthesis pathways, and seven genes linked to chromatophore development and maintenance. In comparisons of expression profiles between pigment-rich and white skin, pigment-rich skin upregulated the pteridine pathway as well as xanthophore/erythrophore development and maintenance genes; in comparisons between orange and yellow skin, orange skin upregulated the pteridine and carotenoid pathways as well as melanophore maintenance genes. Our results corroborate the predictions that squamates can produce similar colors using distinct color-reflecting molecules, and that both color and color-pattern genes are likely conserved across vertebrates. Furthermore, this study provides a concise list of candidate genes for future functional verification, representing a first step in determining the genetic basis of color and color pattern in anoles.

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.