Variation in pigmentation gene expression is associated with distinct aposematic color morphs in the poison frog Dendrobates auratus

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 genomic 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 genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, anti-predator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess 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 emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address 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, https://mvs.unimelb.edu.au/amphibian-genomics-consortium ) in early 2023. This burgeoning community already has more than 282 members from 41 countries. 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 call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

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 genomic 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 genomic 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 essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess 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 emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address 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, https://mvs.unimelb.edu.au/amphibian-genomics-consortium) in early 2023. This burgeoning community already has more than 282 members from 41 countries. 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 call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

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.

Correlated evolution of conspicuous colouration and burrowing in crayfish

Conspicuous colours have fascinated biologists for centuries, leading to research on the evolution and functional significance of colour traits. In many cases, research suggests that conspicuous colours are adaptive and serve a function in sexual or aposematic signalling. In other cases, a lack of evidence for the adaptive value of conspicuous colours garners interest from biologists, such as when organisms that live underground and are rarely exposed to the surface are nevertheless colourful. Here, we use phylogenetic comparative methods to investigate colour evolution throughout freshwater crayfishes that vary in burrowing ability. Within the taxa we analysed, conspicuous colours have evolved independently over 50 times, and these colours are more common in semi-terrestrial crayfishes that construct extensive burrows. The intuitive but not evolutionarily justified assumption when presented with these results is to assume that these colours are adaptive. But contrary to this intuition, we discuss the hypothesis that colouration in crayfish is neutral. Supporting these ideas, the small population sizes and reduced gene flow within semi-terrestrial burrowing crayfishes may lead to the fixation of colour-phenotype mutations. Overall, our work brings into question the traditional view of animal colouration as a perfectly adapted phenotype.

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.

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.

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.

Ionizing radiation and melanism in Chornobyl tree frogs

Human actions are altering ecosystems worldwide. Among human-released pollutants, ionizing radiation arises as a rare but potentially devastating threat to natural systems. The Chornobyl accident (1986) represents the largest release of radioactive material to the environment. Our aim was to examine how exposure to radiation from the Chornobyl accident influences dorsal skin coloration of Eastern tree frog (Hyla orientalis) males sampled across a wide gradient of radioactive contamination in northern Ukraine. We assessed the relationship between skin frog coloration (which can act as a protective mechanism against ionizing radiation), radiation conditions and oxidative stress levels. Skin coloration was darker in localities closest to areas with high radiation levels at the time of the accident, whereas current radiation levels seemed not to influence skin coloration in Chornobyl tree frogs. Tree frogs living within the Chornobyl Exclusion Zone had a remarkably darker dorsal skin coloration than frogs from outside the Zone. The maintenance of dark skin coloration was not linked to physiological costs in terms of frog body condition or oxidative status, and we did not detect short-term changes in frog coloration. Dark coloration is known to protect against different sources of radiation by neutralizing free radicals and reducing DNA damage, and, particularly melanin pigmentation has been proposed as a buffering mechanism against ionizing radiation. Our results suggest that exposure to high levels of ionizing radiation, likely at the time of the accident, may have been selected for darker coloration in Chornobyl tree frogs. Further studies are needed to determine the underlying mechanisms and evolutionary consequences of the patterns found here.

Population genetic differentiation and genomic signatures of adaptation to climate in an abundant lizard

Species distributed across climatic gradients will typically experience spatial variation in selection, but gene flow can prevent such selection from causing population genetic differentiation and local adaptation. Here, we studied genomic variation of 415 individuals across 34 populations of the common wall lizard (Podarcis muralis) in central Italy. This species is highly abundant throughout this region and populations belong to a single genetic lineage, yet there is extensive phenotypic variation across climatic regimes. We used redundancy analysis to, first, quantify the effect of climate and geography on population genomic variation in this region and, second, to test if climate consistently sorts specific alleles across the landscape. Climate explained 5% of the population genomic variation across the landscape, about half of which was collinear with geography. Linear models and redundancy analyses identified loci that were significantly differentiated across climatic regimes. These loci were distributed across the genome and physically associated with genes putatively involved in thermal tolerance, regulation of temperature-dependent metabolism and reproductive activity, and body colouration. Together, these findings suggest that climate can exercise sufficient selection in lizards to promote genetic differentiation across the landscape in spite of high gene flow.

Molecular physiology of pumiliotoxin sequestration in a poison frog

Poison frogs bioaccumulate alkaloids for chemical defense from their arthropod diet. Although many alkaloids are accumulated without modification, some poison frog species can metabolize pumiliotoxin (PTX 251D) into the more potent allopumiliotoxin (aPTX 267A). Despite extensive research characterizing the chemical arsenal of poison frogs, the physiological mechanisms involved in the sequestration and metabolism of individual alkaloids remain unclear. We first performed a feeding experiment with the Dyeing poison frog (Dendrobates tinctorius) to ask if this species can metabolize PTX 251D into aPTX 267A and what gene expression changes are associated with PTX 251D exposure in the intestines, liver, and skin. We found that D. tinctorius can metabolize PTX 251D into aPTX 267A, and that PTX 251D exposure changed the expression level of genes involved in immune system function and small molecule metabolism and transport. To better understand the functional significance of these changes in gene expression, we then conducted a series of high-throughput screens to determine the molecular targets of PTX 251D and identify potential proteins responsible for metabolism of PTX 251D into aPTX 267A. Although screens of PTX 251D binding human voltage-gated ion channels and G-protein coupled receptors were inconclusive, we identified human CYP2D6 as a rapid metabolizer of PTX 251D in a cytochrome P450 screen. Furthermore, a CYP2D6-like gene had increased expression in the intestines of animals fed PTX, suggesting this protein may be involved in PTX metabolism. These results show that individual alkaloids can modify gene expression across tissues, including genes involved in alkaloid metabolism. More broadly, this work suggests that specific alkaloid classes in wild diets may induce physiological changes for targeted accumulation and metabolism.

Pterin-based pigmentation in animals

Pterins are one of the major sources of bright coloration in animals. They are produced endogenously, participate in vital physiological processes and serve a variety of signalling functions. Despite their ubiquity in nature, pterin-based pigmentation has received little attention when compared to other major pigment classes. Here, we summarize major aspects relating to pterin pigmentation in animals, from its long history of research to recent genomic studies on the molecular mechanisms underlying its evolution. We argue that pterins have intermediate characteristics (endogenously produced, typically bright) between two well-studied pigment types, melanins (endogenously produced, typically cryptic) and carotenoids (dietary uptake, typically bright), providing unique opportunities to address general questions about the biology of coloration, from the mechanisms that determine how different types of pigmentation evolve to discussions on honest signalling hypotheses. Crucial gaps persist in our knowledge on the molecular basis underlying the production and deposition of pterins. We thus highlight the need for functional studies on systems amenable for laboratory manipulation, but also on systems that exhibit natural variation in pterin pigmentation. The wealth of potential model species, coupled with recent technological and analytical advances, make this a promising time to advance research on pterin-based pigmentation in animals.

Untangling the structural and molecular mechanisms underlying colour and rapid colour change in a lizard, Agama atra

With functions as diverse as communication, protection and thermoregulation, coloration is one of the most important traits in lizards. The ability to change colour as a function of varying social and environmental conditions is thus an important innovation. While colour change is present in animals ranging from squids, to fish and reptiles, not much is known about the mechanisms behind it. Traditionally, colour change was attributed to migration of pigments, in particular melanin. More recent work has shown that the changes in nanostructural configuration inside iridophores are able to produce a wide palette of colours. However, the genetic mechanisms underlying colour, and colour change in particular, remain unstudied. Here we use a combination of transcriptomic and microscopic data to show that melanin, iridophores and pteridines are the main colour-producing mechanisms in Agama atra, and provide molecular and structural data suggesting that rapid colour change is achieved via melanin dispersal in combination with iridophore organization. This work demonstrates the power of combining genotypic (gene expression) and phenotypic (microscopy) information for addressing physiological questions, providing a basis for future studies of colour change.

Ultraviolet components offer minimal contrast enhancement to an aposematic signal

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

Molecular Plasticity in Animal Pigmentation: Emerging Processes Underlying Color Changes

Animal coloration has been rigorously studied and has provided morphological implications for fitness with influences over social behavior, predator–prey interactions, and sexual selection. In vertebrates, its study has developed our understanding across diverse fields ranging from behavior to molecular biology. In the search for underlying molecular mechanisms, many have taken advantage of pedigree-based and genome-wide association screens to reveal the genetic architecture responsible for pattern variation that occurs in early development. However, genetic differences do not provide a full picture of the dynamic changes in coloration that are most prevalent across vertebrates at the molecular level. Changes in coloration that occur in adulthood via phenotypic plasticity rely on various social, visual, and dietary cues independent of genetic variation. Here, I will review the contributions of pigment cell biology to animal color changes and recent studies describing their molecular underpinnings and function. In this regard, conserved epigenetic processes such as DNA methylation play a role in lending plasticity to gene regulation as it relates to chromatophore function. Lastly, I will present African cichlids as emerging models for the study of pigmentation and molecular plasticity for animal color changes. I posit that these processes, in a dialog with environmental stimuli, are important regulators of variation and the selective advantages that accompany a change in coloration for vertebrate animals.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Mechanisms for Color Convergence in a Mimetic Radiation of Poison Frogs

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

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.