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 color 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, recql4, xdh), and carotenoid metabolism (e.g., dgat2, rbp1, scarb2). In addition, using weighted gene 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, 14 were known to be related to color production (gch1, gmps, gpr143, impdh1, mc1r, pax3-a, pax7, ppat, rab27a, rlbp1, tfec, trpm1, xdh).

Genome assembly of the foot-flagging frog, Staurois parvus: a resource for understanding mechanisms of behavior.

Elaborate and skilled movements of the body have been selected in a variety of species as courtship and rivalry signals. One roadblock in studying these behaviors has been a lack of resources for understanding how they evolved at the genetic level. The Bornean rock frog (Staurois parvus) is an ideal species in which to address this issue. Males wave their hindlimbs in a "foot-flagging" display when competing for mates. The evolution of foot flagging in S. parvus and other species is accompanied by increases in the expression of the androgen receptor gene within its neuromuscular system, but it remains unclear what genetic or transcriptional changes are associated with this behavioral phenotype. We have now assembled the genome of S. parvus, resulting in 3.98 Gbp of 22,402 contigs with an N50 of 611,229 bp. The genome will be a resource for finding genes related to the physiology underlying foot flagging and to adaptations of the neuromuscular system. As a first application of the genome, we also began work in comparative genomics and differential gene expression analysis. We show that the androgen receptor is diverged from other anuran species, and we identify unique expression patterns of genes in the spinal cord and leg muscle that are important for axial patterning, cell specification and morphology, or muscle contraction. This genome will continue to be an important tool for future -omics studies to understand the evolution of elaborate signaling behaviors in this and potentially related species.

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.

Investigating signal modalities of aposematism in a poison frog.

Aposematic species combine a conspicuous signal with a secondary defence, the majority of which are studied in the context of a visual signal. While multimodality of the aposematic signal appears to be common in invertebrate species, we know very little about the presence or absence of multimodality in vertebrates. Here, we examine the possibility of multimodality of aposematism in the green-and-black poison frog, Dendrobates auratus. Using a non-visual predator (the cat-eyed snake, Leptodeira annulata) and extractions of chemicals in frog skins, we test whether there is sufficient non-visual information for predators to avoid this aposematic species without using visual cues. We found that experienced predators avoid chemicals in this poison frog's skin by olfactory cues alone in trials with live frogs and extracts from captive poison frogs, whereas extracts from wild poison frogs did not lead to avoidance behaviours in predators. Further, in our limited sampling, naïve predators demonstrate no avoidance. This not only indicates that predators can make informed decisions from the frog's odour, but also indicates that avoidance based on olfactory cue is a learned response.

Diversity of cave Phlebotomines (Diptera: Psychodidae) from a Colombian cave.

Sandflies are vector species of Leishmania, among many other pathogens, with a global distribution and a variety of ecological niches. Previous samplings have found that karstic formations (i.e., caves, grottos, and folds formed by the erosion of limestone) serve as a natural habitat to sandfly species. The majority of samplings of cave sandfly diversity have occurred in Brazil and to date none have studied the species composition in a cave in the Northern Andes. We collected sandflies in the Cave "Los Guácharos", in the state of Antioquia, Colombia. The sampling was carried out during two consecutive nights in September 2019. CDC-type light traps were installed inside the cavern and in other surrounding karst systems (caves, rock-breaks, and folds). In total, we identified 17 species of sandfly from the cave and surrounding karst systems, including a new record for Colombia (Bichromomyia olmeca), and provide the first karstic reports for four other species (Lutzomyia gomezi, Lutzomyia hartmanni, Pintomyia ovallesi, and Psychodopygus panamensis). We then used the results of our survey and published literature to test two hypotheses. First, that sandfly diversity in Neotropical caves is richest nearer to the equator, and second that there is a phylogenetic signal of karstic habitat use in sandflies. Counter to our predictions, we found no evidence that diversity follows a latitudinal gradient. Further, we find no evidence of a phylogenetic signal of karstic habitat use, instead finding that the use of caves likely evolved multiple times across several genera. Our results highlight the importance of a wide sampling to understand the natural habitat of sandflies and other disease vectors.

Evolution: Environmental conditions determine how Wolbachia interacts with its host.

Wolbachia is a genus of globally distributed endosymbiotic bacteria found in thousands of species of arthropods and nematodes. One strain, wMel, has spread globally but varies dramatically in frequency. In a new paper, Hague et al. discovered that the interaction of host genotype, symbiont genotype, and environment explains natural wMel variation.

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.

Molecular mechanisms of local adaptation for salt-tolerance in a treefrog.

Salinization is a global phenomenon affecting ecosystems and forcing freshwater organisms to deal with increasing levels of ionic stress. However, our understanding of mechanisms that permit salt tolerance in amphibians is limited. This study investigates mechanisms of salt tolerance in locally adapted, coastal populations of a treefrog, Hyla cinerea. Using a common garden experiment, we (i) determine the extent that environment (i.e., embryonic and larval saltwater exposure) or genotype (i.e., coastal vs. inland) affects developmental benchmarks and transcriptome expression, and (ii) identify genes that may underpin differences in saltwater tolerance. Differences in gene expression, survival, and plasma osmolality were most strongly associated with genotype. Population genetic analyses on expressed genes also delineated coastal and inland groups based on genetic similarity. Coastal populations differentially expressed osmoregulatory genes including ion transporters (atp1b1, atp6V1g2, slc26a), cellular adhesion components (cdh26, cldn1, gjb3, ocln), and cytoskeletal components (odc1-a, tgm3). Several of these genes are the same genes expressed by euryhaline fish after exposure to freshwater, which is a novel finding for North American amphibians and suggests that these genes may be associated with local salinity adaptation. Coastal populations also highly expressed glycerol-3-phosphate dehydrogenase 1 (gpd1), which indicates they use glycerol as a compatible osmolyte to reduce water loss - another mechanism of saltwater tolerance previously unknown in frogs. These data signify that Hyla cinerea inhabiting coastal, brackish wetlands have evolved a salt-tolerant ecotype, and highlights novel candidate pathways that can lead to salt tolerance in freshwater organisms facing habitat salinization.

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

BACKGROUND: Color and pattern phenotypes have clear implications for survival and reproduction in many species. However, the mechanisms that produce this coloration are still poorly characterized, especially at the genomic level. Here we have taken a transcriptomics-based approach to elucidate the underlying genetic mechanisms affecting color and pattern in a highly polytypic poison frog. We sequenced RNA from the skin from four different color morphs during the final stage of metamorphosis and assembled a de novo transcriptome. We then investigated differential gene expression, with an emphasis on examining candidate color genes from other taxa. RESULTS: Overall, we found differential expression of a suite of genes that control melanogenesis, melanocyte differentiation, and melanocyte proliferation (e.g., tyrp1, lef1, leo1, and mitf) as well as several differentially expressed genes involved in purine synthesis and iridophore development (e.g., arfgap1, arfgap2, airc, and gart). CONCLUSIONS: Our results provide evidence that several gene networks known to affect color and pattern in vertebrates play a role in color and pattern variation in this species of poison frog.

Evolution and assessment of colour patterns in stream-resident and anadromous male threespine stickleback Gasterosteus aculeatus from three regions.

We compared the colour patterns of free swimming, reproductively active male threespine stickleback Gasterosteus aculeatus of the anadromous and stream ecotypes from three geographically distinct regions. Consistent with the hypothesis of environmentally mediated selection, our results indicate ecologically replicated differences in G. aculeatus coloration between anadromous and stream-resident populations, and that G. aculeatus probably have the visual acuity to discriminate colour pattern differences between anadromous and stream-resident fish.

Diversity in warning coloration: selective paradox or the norm?

Aposematic theory has historically predicted that predators should select for warning signals to converge on a single form, as a result of frequency-dependent learning. However, widespread variation in warning signals is observed across closely related species, populations and, most problematically for evolutionary biologists, among individuals in the same population. Recent research has yielded an increased awareness of this diversity, challenging the paradigm of signal monomorphy in aposematic animals. Here we provide a comprehensive synthesis of these disparate lines of investigation, identifying within them three broad classes of explanation for variation in aposematic warning signals: genetic mechanisms, differences among predators and predator behaviour, and alternative selection pressures upon the signal. The mechanisms producing warning coloration are also important. Detailed studies of the genetic basis of warning signals in some species, most notably Heliconius butterflies, are beginning to shed light on the genetic architecture facilitating or limiting key processes such as the evolution and maintenance of polymorphisms, hybridisation, and speciation. Work on predator behaviour is changing our perception of the predator community as a single homogenous selective agent, emphasising the dynamic nature of predator-prey interactions. Predator variability in a range of factors (e.g. perceptual abilities, tolerance to chemical defences, and individual motivation), suggests that the role of predators is more complicated than previously appreciated. With complex selection regimes at work, polytypisms and polymorphisms may even occur in Müllerian mimicry systems. Meanwhile, phenotypes are often multifunctional, and thus subject to additional biotic and abiotic selection pressures. Some of these selective pressures, primarily sexual selection and thermoregulation, have received considerable attention, while others, such as disease risk and parental effects, offer promising avenues to explore. As well as reviewing the existing evidence from both empirical studies and theoretical modelling, we highlight hypotheses that could benefit from further investigation in aposematic species. Finally by collating known instances of variation in warning signals, we provide a valuable resource for understanding the taxonomic spread of diversity in aposematic signalling and with which to direct future research. A greater appreciation of the extent of variation in aposematic species, and of the selective pressures and constraints which contribute to this once-paradoxical phenomenon, yields a new perspective for the field of aposematic signalling.

Alkaloid defenses of co-mimics in a putative Müllerian mimetic radiation.

BACKGROUND: Polytypism in aposematic species is unlikely according to theory, but commonly seen in nature. Ranitomeya imitator is a poison frog species exhibiting polytypic mimicry of three congeneric model species (R. fantastica, R. summersi, and two morphs of R. variabilis) across four allopatric populations (a "mimetic radiation"). In order to investigate chemical defenses in this system, a key prediction of Müllerian mimicry, we analyzed the alkaloids of both models and mimics from four allopatric populations. RESULTS: In this study we demonstrate distinct differences in alkaloid profiles between co-mimetic species within allopatric populations. We further demonstrate that R. imitator has a greater number of distinct alkaloid types than the model species and more total alkaloids in all but one population. CONCLUSIONS: Given that R. imitator is the more abundant species in these populations, R. imitator is likely driving the majority of predator-learned avoidance in these complexes. The success of Ranitomeya imitator as a putative advergent mimic may be a direct result of differences in alkaloid sequestration. Furthermore, we propose that automimicry within co-mimetic species is an important avenue of research.