Phylogenomic Reconstruction of the Neotropical Poison Frogs (Dendrobatidae) and Their Conservation

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

Distribution and diversification of Adelphobates, emblematic poison frogs from Brazilian Amazonia

Adelphobates contains three species, and the inaccurate identification of A. quinquevittatus and the scarcity of records of A. castaneoticus complicate inference of their distributions; the latter species occurs in sympatry with A. galactonotus. Our objective was to revise the distributions of Adelphobates by compiling data and modeling habitat suitability, as range limits may be shaped by landscape features and biotic interactions. We initially analyzed the existence of operational taxonomic units within the nominal species and subsequently inferred the observed and potential distributions, taking into account the possible independent lineages for the three species, and we also generated a molecular timetree to understand the chronology of interspecific diversification events. Adelphobates quinquevittatus was found to have a more easterly distribution than previously described, and specimens with phenotypic variation were found to occur in areas inconsistent with the modeling, and A. castaneoticus was concentrated in the Tapajós-Xingu interfluve, surrounded by A. galactonotus. Models indicated that the right bank of the Xingu River is suitable for both species, indeed, both were found there. Despite Adelphobates species having their distributions delimited by major Amazonian rivers, estimated divergence times predate the formation of the modern river network, suggesting that other mechanisms were involved in their diversification.

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 novel probe set for the phylogenomics and evolution of RTA spiders

Spiders are important models for evolutionary studies of web building, sexual selection and adaptive radiation. The recent development of probes for UCE (ultra-conserved element)-based phylogenomic studies has shed light on the phylogeny and evolution of spiders. However, the two available UCE probe sets for spider phylogenomics (Spider and Arachnida probe sets) have relatively low capture efficiency within spiders, and are not optimized for the retrolateral tibial apophysis (RTA) clade, a hyperdiverse lineage that is key to understanding the evolution and diversification of spiders. In this study, we sequenced 15 genomes of species in the RTA clade, and using eight reference genomes, we developed a new UCE probe set (41 845 probes targeting 3802 loci, labelled as the RTA probe set). The performance of the RTA probes in resolving the phylogeny of the RTA clade was compared with the Spider and Arachnida probes through an in-silico test on 19 genomes. We also tested the new probe set empirically on 28 spider species of major spider lineages. The results showed that the RTA probes recovered twice and four times as many loci as the other two probe sets, and the phylogeny from the RTA UCEs provided higher support for certain relationships. This newly developed UCE probe set shows higher capture efficiency empirically and is particularly advantageous for phylogenomic and evolutionary studies of RTA clade and jumping spiders.

Identification of the NA+/K+-ATPase α-Isoforms in Six Species of Poison Dart Frogs and their Sensitivity to Cardiotonic Steroids

Cardiotonic steroids (CTS) are a group of compounds known to be toxic due to their ability to inhibit the Na⁺/K⁺-ATPase (NKA), which is essential to maintain the balance of ions in animal cells. An evolutionary strategy of molecular adaptation to avoid self-intoxication acquired by CTS defended organisms and their predators is the structural modification of their NKA where specific amino acid substitutions confer resistant phenotypes. Several lineages of poison dart frogs (Dendrobatidae) are well known to sequester a wide variety of lipophilic alkaloids from their arthropod diet, however there is no evidence of CTS-sequestration or dietary exposure. Interestingly this study identified the presence of α-NKA isoforms (α₁ and α₂) with amino acid substitutions indicative of CTS-resistant phenotypes in skeletal muscle transcriptomes obtained from six species of dendrobatids: Phyllobates aurotaenia, Oophaga anchicayensis, Epipedobates boulengeri, Andinobates bombetes, Andinobates minutus, and Leucostethus brachistriatus, collected in the Valle del Cauca (Colombia). P. aurotaenia, A. minutus, and E. boulengeri presented two variants for α₁-NKA, with one of them having these substitutions. In contrast, O. anchicayensis and A. bombetes have only one α₁-NKA isoform with an amino acid sequence indicative of CTS susceptibility and an α₂-NKA with one substitution that could confer a reduced affinity for CTS. The α₁ and α₂ isoforms of L. brachistriatus do not contain substitutions imparting CTS resistance. Our findings indicate that poison dart frogs express α-NKA isoforms with different affinities for CTS and the pattern of this expression might be influenced by factors related to evolutionary, physiological, ecological, and geographical burdens.

Size-dependent colouration balances conspicuous aposematism and camouflage

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

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.

Red List assessment of amphibian species of Ecuador: A multidimensional approach for their conservation

Ecuador is one of the most biodiverse countries in the world, but faces severe pressures and threats to its natural ecosystems. Numerous species have declined and require to be objectively evaluated and quantified, as a step towards the development of conservation strategies. Herein, we present an updated National Red List Assessment for amphibian species of Ecuador, with one of the most detailed and complete coverages for any Ecuadorian taxonomic group to date. Based on standardized methodologies that integrate taxonomic work, spatial analyses, and ecological niche modeling, we assessed the extinction risk and identified the main threats for all Ecuadorian native amphibians (635 species), using the IUCN Red List Categories and Criteria. Our evaluation reveals that 57% (363 species) are categorized as Threatened, 12% (78 species) as Near Threatened, 4% (26 species) as Data Deficient, and 27% (168 species) as Least Concern. Our assessment almost doubles the number of threatened species in comparison with previous evaluations. In addition to habitat loss, the expansion of the agricultural/cattle raising frontier and other anthropogenic threats (roads, human settlements, and mining/oil activities) amplify the incidence of other pressures as relevant predictors of ecological integrity. Potential synergic effects with climate change and emergent diseases (apparently responsible for the sudden declines), had particular importance amongst the threats sustained by Ecuadorian amphibians. Most threatened species are distributed in montane forests and paramo habitats of the Andes, with nearly 10% of them occurring outside the National System of Protected Areas of the Ecuadorian government. Based on our results, we recommend the following actions: (i) An increase of the National System of Protected Areas to include threatened species. (ii) Supporting the ex/in-situ conservation programs to protect species considered like Critically Endangered and Endangered. (iii) Focalizing research efforts towards the description of new species, as well as species currently categorized as Data Deficient (DD) that may turn out to be threatened. The implementation of the described actions is challenging, but urgent, given the current conservation crisis faced by amphibians.

Reproductive behavior drives female space use in a sedentary Neotropical frog

Longer-range movements of anuran amphibians such as mass migrations and habitat invasion have received a lot of attention, but fine-scale spatial behavior remains largely understudied. This gap is especially striking for species that show long-term site fidelity and display their whole behavioral repertoire in a small area. Studying fine-scale movement with conventional capture-mark-recapture techniques is difficult in inconspicuous amphibians: individuals are hard to find, repeated captures might affect their behavior and the number of data points is too low to allow a detailed interpretation of individual space use and time budgeting. In this study, we overcame these limitations by equipping females of the Brilliant-Thighed Poison Frog (Allobates femoralis) with a tag allowing frequent monitoring of their location and behavior. Neotropical poison frogs are well known for their complex behavior and diverse reproductive and parental care strategies. Although the ecology and behavior of the polygamous leaf-litter frog Allobates femoralis is well studied, little is known about the fine-scale space use of the non-territorial females who do not engage in acoustic and visual displays. We tracked 17 females for 6 to 17 days using a harmonic direction finder to provide the first precise analysis of female space use in this species. Females moved on average 1 m per hour and the fastest movement, over 20 m per hour, was related to a subsequent mating event. Traveled distances and activity patterns on days of courtship and mating differed considerably from days without reproduction. Frogs moved more on days with lower temperature and more precipitation, but mating seemed to be the main trigger for female movement. We observed 21 courtships of 12 tagged females. For seven females, we observed two consecutive mating events. Estimated home ranges after 14 days varied considerably between individuals and courtship and mating associated space use made up for ∼30% of the home range. Allobates femoralis females spent large parts of their time in one to three small centers of use. Females did not adjust their time or space use to the density of males in their surroundings and did not show wide-ranging exploratory behavior. Our study demonstrates how tracking combined with detailed behavioral observations can reveal the patterns and drivers of fine-scale spatial behavior in sedentary species.

Complex patterns of differentiation and gene flow underly the divergence of aposematic phenotypes in Oophaga poison frogs

Hybridization and introgression can have complex consequences for both species evolution and conservation. Here, we investigated the origin and characteristics of a putative hybrid zone between two South American poison dart frog species, Oophaga anchicayensis and the critically endangered Oophaga lehmanni, which are heavily sought after on the illegal pet market. Using a combination of phenotypic (49 traits) and genomic (ddRADseq) data, we found that the putative hybrids are morphologically distinct from their parental species and confirmed genomic signatures of admixture in these populations. Several lines of evidence (hybrid indices, interspecific hybrid heterozygosity, genomic clines, comparisons with simulated hybrids and demographic modelling) support the conclusion that these populations are not comprised of early-generation hybrids and thus, they probably did not arise as a result of illegal translocations associated with wildlife trafficking. Instead, they probably represent an independent lineage which has persisted through isolation and has only relatively recently re-established gene flow with both parental species. Furthermore, we detected signals of differential introgression from parental species into these hybrid populations which suggest relaxed stabilizing selection on these aposematic colour morphs, potentially via context-dependent female choice. These populations thus provide a fascinating window into the role of hybridization, isolation and female choice in the diversification of South American poison dart frogs. In addition, our results underline the importance of landscape conservation measures to protect, not only known localities of nominal species, but also the phenotypic and genomic variation harbored by admixed lineages which represent crucial repositories for the impressive diversity in this system.

Phylogenetic relationships and systematics of the Amazonian poison frog genus Ameerega using ultraconserved genomic elements

The Amazonian poison frog genus Ameerega is one of the largest yet most understudied of the brightly colored genera in the anuran family Dendrobatidae, with 30 described species ranging throughout tropical South America. Phylogenetic analyses of Ameerega are highly discordant, lacking consistency due to variation in data types and methods, and often with limited coverage of species diversity in the genus. Here, we present a comprehensive phylogenomic reconstruction of Ameerega, utilizing state-of-the-art sequence capture techniques and phylogenetic methods. We sequenced thousands of ultraconserved elements from over 100 tissue samples, representing almost every described Ameerega species, as well as undescribed cryptic diversity. We generated topologies using maximum likelihood and coalescent methods and compared the use of maximum likelihood and Bayesian methods for estimating divergence times. Our phylogenetic inference diverged strongly from those of previous studies, and we recommend steps to bring Ameerega taxonomy in line with the new phylogeny. We place several species in a phylogeny for the first time, as well as provide evidence for six potential candidate species. We estimate that Ameerega experienced a rapid radiation approximately 7-11 million years ago and that the ancestor of all Ameerega was likely an aposematic, montane species. This study underscores the utility of phylogenomic data in improving our understanding of the phylogeny of understudied clades and making novel inferences about their evolution.