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 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.

Shifts in microbiome community composition can have large effects on host health. It is therefore important to understand how perturbations, like those caused by the introduction of exogenous chemicals, modulate microbiome community composition. In poison frogs within the family Dendrobatidae, the skin microbiome is exposed to the alkaloids that the frogs sequester from their diet and use for defense. Given the demonstrated antimicrobial effects of these poison frog alkaloids, these compounds may be structuring the skin microbial community. To test this, we first characterized microbial communities from chemically defended and closely related non-defended frogs from Ecuador. Then we conducted a laboratory experiment to monitor the effect of the alkaloid decahydroquinoline (DHQ) on the microbiome of a single frog species. In both the field and lab experiments, we found that alkaloid-exposed microbiomes are more species rich and phylogenetically diverse, with an increase in rare taxa. To better understand the strain-specific behavior in response to alkaloids, we cultured microbial strains from poison frog skin and found the majority of strains exhibited either enhanced growth or were not impacted by the addition of DHQ. Additionally, stable isotope tracing coupled to nanoSIMS suggests that some of these strains are able to metabolize DHQ. Taken together, these data suggest that poison frog chemical defenses open new niches for skin-associated microbes with specific adaptations, including the likely metabolism of alkaloids, that enable their survival in this toxic environment. This work helps expand our understanding of how exposure to exogenous compounds like alkaloids can impact host microbiomes.

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.

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.

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.

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.

Colour is an important component of many different defensive strategies, but signal efficacy and detectability will also depend on the size of the coloured structures, and how pattern size interacts with the background. Consequently, size-dependent changes in colouration are common among many different species as juveniles and adults frequently use colour for different purposes in different environmental contexts. A widespread strategy in many species is switching from crypsis to conspicuous aposematic signalling as increasing body size can reduce the efficacy of camouflage, while other antipredator defences may strengthen. Curiously, despite being chemically defended, the gold-striped frog (Lithodytes lineatus, Leptodactylidae) appears to do the opposite, with bright yellow stripes found in smaller individuals, whereas larger frogs exhibit dull brown stripes. Here, we investigated whether size-dependent differences in colour support distinct defensive strategies. We first used visual modelling of potential predators to assess how colour contrast varied among frogs of different sizes. We found that contrast peaked in mid-sized individuals while the largest individuals had the least contrasting patterns. We then used two detection experiments with human participants to evaluate how colour and body size affected overall detectability. These experiments revealed that larger body sizes were easier to detect, but that the colours of smaller frogs were more detectable than those of larger frogs. Taken together our data support the hypothesis that the primary defensive strategy changes from conspicuous aposematism to camouflage with increasing size, implying size-dependent differences in the efficacy of defensive colouration. We discuss our data in relation to theories of size-dependent aposematism and evaluate the evidence for and against a possible size-dependent mimicry complex with sympatric poison frogs (Dendrobatidae).

Sequestration of chemical defenses from dietary sources is dependent on the availability of compounds in the environment and the mechanism of sequestration. Previous experiments have shown that sequestration efficiency varies among alkaloids in poison frogs, but little is known about the underlying mechanism. The aim of this study was to quantify the extent to which alkaloid sequestration and modification are dependent on alkaloid availability and/or sequestration mechanism. To do this, we administered different doses of histrionicotoxin (HTX) 235A and decahydroquinoline (DHQ) to captive-bred Adelphobates galactonotus and measured alkaloid quantity in muscle, kidney, liver, and feces. HTX 235A and DHQ were detected in all organs, whereas only DHQ was present in trace amounts in feces. For both liver and skin, the quantity of alkaloid accumulated increased at higher doses for both alkaloids. Accumulation efficiency in the skin increased at higher doses for HTX 235A but remained constant for DHQ. In contrast, the efficiency of HTX 235A accumulation in the liver was inversely related to dose and a similar, albeit statistically nonsignificant, pattern was observed for DHQ. We identified and quantified the N-methylation of DHQ in A. galactonotus, which represents a previously unknown example of alkaloid modification in poison frogs. Our study suggests that variation in alkaloid composition among individuals and species can result from differences in sequestration efficiency related to the type and amount of alkaloids available in the environment.

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.