Deep divergences among inconspicuously colored clades of Epipedobates poison frogs.

Poison frogs (Dendrobatidae) are famous for their aposematic species, having a combination of diverse color patterns and defensive skin toxins, yet most species in this family are inconspicuously colored and considered non-aposematic. Epipedobates is among the youngest genus-level clades of Dendrobatidae that includes both aposematic and inconspicuous species. Using Sanger-sequenced mitochondrial and nuclear markers, we demonstrate deep genetic divergences among inconspicuous species of Epipedobates but relatively shallow genetic divergences among conspicuous species. Our phylogenetic analysis includes broad geographic sampling of the inconspicuous lineages typically identified as E. boulengeri and E. espinosai, which reveals two putative new species, one in west-central Colombia (E. sp. 1) and the other in north-central Ecuador (E. aff. espinosai). We conclude that E. darwinwallacei is a junior subjective synonym of E. espinosai. We also clarify the geographic distributions of inconspicuous Epipedobates species including the widespread E. boulengeri. We provide a qualitative assessment of the phenotypic diversity in each nominal species, with a focus on the color and pattern of inconspicuous species. We conclude that Epipedobates contains eight known valid species, six of which are inconspicuous. A relaxed molecular clock analysis suggests that the most recent common ancestor of Epipedobates is âˆ¼11.1 million years old, which nearly doubles previous estimates. Last, genetic information points to a center of species diversity in the Chocó at the southwestern border of Colombia with Ecuador. A Spanish translation of this text is available in the supplementary materials.

Selection on Visual Opsin Genes in Diurnal Neotropical Frogs and Loss of the SWS2 Opsin in Poison Frogs.

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.

2b or not 2b? 2bRAD is an effective alternative to ddRAD for phylogenomics.

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.

ExplorATE: a new pipeline to explore active transposable elements from RNA-seq data.

MOTIVATION: Transposable elements (TEs) are ubiquitous in genomes and many remain active. TEs comprise an important fraction of the transcriptomes with potential effects on the host genome, either by generating deleterious mutations or promoting evolutionary novelties. However, their functional study is limited by the difficulty in their identification and quantification, particularly in non-model organisms. RESULTS: We developed a new pipeline [explore active transposable elements (ExplorATE)] implemented in R and bash that allows the quantification of active TEs in both model and non-model organisms. ExplorATE creates TE-specific indexes and uses the Selective Alignment (SA) to filter out co-transcribed transposons within genes based on alignment scores. Moreover, our software incorporates a Wicker-like criteria to refine a set of target TEs and avoid spurious mapping. Based on simulated and real data, we show that the SA strategy adopted by ExplorATE achieved better estimates of non-co-transcribed elements than other available alignment-based or mapping-based software. ExplorATE results showed high congruence with alignment-based tools with and without a reference genome, yet ExplorATE required less execution time. Likewise, ExplorATE expands and complements most previous TE analyses by incorporating the co-transcription and multi-mapping effects during quantification, and provides a seamless integration with other downstream tools within the R environment. AVAILABILITY AND IMPLEMENTATION: Source code is available at https://github.com/FemeniasM/ExplorATEproject and https://github.com/FemeniasM/ExplorATE_shell_script. Data available on request. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Computational Modeling of the Interaction of Silver Clusters with Carbohydrates.

Silver nanoparticles are recognized for their numerous physical, biological, and pharmaceutical applications. In the present study, the interaction of silver clusters with monosaccharide molecules is examined to identify which molecule works better as a reducing agent in the application of a green synthesis approach. Geometry optimization of clusters containing one, three, and five silver atoms is performed along with the optimization of α-d-glucose, α-d-ribose, d-erythrose, and glyceraldehyde using density functional theory. Optimized geometries allow identifying the interaction formed in the silver cluster and monosaccharide complexes. An electron localization function analysis is performed to further analyze the interaction found and explain the reduction process in the formation of silver nanoparticles. The overall results indicate that glyceraldehyde presents the best characteristics to serve as the most efficient reducing agent.

Pursuing High-Resolution Structures of Nicotinic Acetylcholine Receptors: Lessons Learned from Five Decades.

Since their discovery, nicotinic acetylcholine receptors (nAChRs) have been extensively studied to understand their function, as well as the consequence of alterations leading to disease states. Importantly, these receptors represent pharmacological targets to treat a number of neurological and neurodegenerative disorders. Nevertheless, their therapeutic value has been limited by the absence of high-resolution structures that allow for the design of more specific and effective drugs. This article offers a comprehensive review of five decades of research pursuing high-resolution structures of nAChRs. We provide a historical perspective, from initial structural studies to the most recent X-ray and cryogenic electron microscopy (Cryo-EM) nAChR structures. We also discuss the most relevant structural features that emerged from these studies, as well as perspectives in the field.

Aposematism facilitates the diversification of parental care strategies in poison frogs.

Many organisms have evolved adaptations to increase the odds of survival of their offspring. Parental care has evolved several times in animals including ectotherms. In amphibians, ~ 10% of species exhibit parental care. Among these, poison frogs (Dendrobatidae) are well-known for their extensive care, which includes egg guarding, larval transport, and specialized tadpole provisioning with trophic eggs. At least one third of dendrobatids displaying aposematism by exhibiting warning coloration that informs potential predators about the presence of defensive skin toxins. Aposematism has a central role in poison frog diversification, including diet specialization, and visual and acoustic communication; and it is thought to have impacted their reproductive biology as well. We tested the latter association using multivariate phylogenetic methods at the family level. Our results show complex relationships between aposematism and certain aspects of the reproductive biology in dendrobatids. In particular, aposematic species tend to use more specialized tadpole-deposition sites, such as phytotelmata, and ferry fewer tadpoles than non-aposematic species. We propose that aposematism may have facilitated the diversification of microhabitat use in dendrobatids in the context of reproduction. Furthermore, the use of resource-limited tadpole-deposition environments may have evolved in tandem with an optimal reproductive strategy characterized by few offspring, biparental care, and female provisioning of food in the form of unfertilized eggs. We also found that in phytotelm-breeders, the rate of transition from cryptic to aposematic phenotype is 17 to 19 times higher than vice versa. Therefore, we infer that the aposematism in dendrobatids might serve as an umbrella trait for the evolution and maintenance of their complex offspring-caring activities.

Noise constrains the evolution of call frequency contours in flowing water frogs: a comparative analysis in two clades.

BACKGROUND: The acoustic adaptation hypothesis (AAH) states that signals should evolve towards an optimal transmission of the intended information from senders to intended receivers given the environmental constraints of the medium that they traverse. To date, most AAH studies have focused on the effect of stratified vegetation on signal propagation. These studies, based on the AAH, predict that acoustic signals should experience less attenuation and degradation where habitats are less acoustically complex. Here, we explored this effect by including an environmental noise dimension to test some AAH predictions in two clades of widespread amphibians (Bufonidae and Ranidae) that actively use acoustic signals for communication. By using data from 106 species in these clades, we focused on the characterization of the differences in dominant frequency (DF) and frequency contour (i.e., frequency modulation [FM] and harmonic performances) of mating calls and compared them between species that inhabit flowing-water or still-water environments. RESULTS: After including temperature, body size, habitat type and phylogenetic relationships, we found that DF differences among species were explained mostly by body size and habitat structure. We also showed that species living in lentic habitats tend to have advertisement calls characterized by well-defined FM and harmonics. Likewise, our results suggest that flowing-water habitats can constrain the evolutionary trajectories of the frequency-contour traits of advertisement calls in these anurans. CONCLUSIONS: Our results may support AAH predictions in frogs that vocalize in noisy habitats because flowing-water environments often produce persistent ambient noise. For instance, these anurans tend to generate vocalizations with less well-defined FM and harmonic traits. These findings may help us understand how noise in the environment can influence natural selection as it shapes acoustic signals in affected species.

Pincho: A Modular Approach to High Quality De Novo Transcriptomics.

Transcriptomic reconstructions without reference (i.e., de novo) are common for data samples derived from non-model biological systems. These assemblies involve massive parallel short read sequence reconstructions from experiments, but they usually employ ad-hoc bioinformatic workflows that exhibit limited standardization and customization. The increasing number of transcriptome assembly software continues to provide little room for standardization which is exacerbated by the lack of studies on modularity that compare the effects of assembler synergy. We developed a customizable management workflow for de novo transcriptomics that includes modular units for short read cleaning, assembly, validation, annotation, and expression analysis by connecting twenty-five individual bioinformatic tools. With our software tool, we were able to compare the assessment scores based on 129 distinct single-, bi- and tri-assembler combinations with diverse k-mer size selections. Our results demonstrate a drastic increase in the quality of transcriptome assemblies with bi- and tri- assembler combinations. We aim for our software to improve de novo transcriptome reconstructions for the ever-growing landscape of RNA-seq data derived from non-model systems. We offer guidance to ensure the most complete transcriptomic reconstructions via the inclusion of modular multi-assembly software controlled from a single master console.

Nutritional niches reveal fundamental domestication trade-offs in fungus-farming ants.

During crop domestication, human farmers traded greater productivity for higher crop vulnerability outside specialized cultivation conditions. We found a similar domestication trade-off across the major co-evolutionary transitions in the farming systems of attine ants. First, the fundamental nutritional niches of cultivars narrowed over ~60 million years of naturally selected domestication, and laboratory experiments showed that ant farmers representing subsequent domestication stages strictly regulate protein harvest relative to cultivar fundamental nutritional niches. Second, ants with different farming systems differed in their abilities to harvest the resources that best matched the nutritional needs of their fungal cultivars. This was assessed by quantifying realized nutritional niches from analyses of items collected from the mandibles of laden ant foragers in the field. Third, extensive field collections suggest that among-colony genetic diversity of cultivars in small-scale farms may offer population-wide resilience benefits that species with large-scale farming colonies achieve by more elaborate and demanding practices to cultivate less diverse crops. Our results underscore that naturally selected farming systems have the potential to shed light on nutritional trade-offs that shaped the course of culturally evolved human farming.

A Cytochrome c-Chlorotoxin Hybrid Protein as a Possible Antiglioma Drug.

Malignant gliomas are the most lethal form of primary brain tumors. Despite advances in cancer therapy, the prognosis of glioma patients has remained poor. Cytochrome c (Cytc), an endogenous heme-based protein, holds tremendous potential to treat gliomas because of its innate capacity to trigger apoptosis. To this end, a hybrid cytochrome c-chlorotoxin (Cytc-CTX) protein was biosynthesized to enable cellular uptake of the cell impenetrable Cytc using CTX transporters. A nucleotide sequence containing 1 : 1 Cytc and CTX was constructed and separated by a hexahistidine-tag and an enterokinase cleavage site. The sequence was cloned into a pBTR1 plasmid, expressed in Escherichia coli, purified via 2-dimensional chromatography. The identity and size of the protein were determined by Western blot and mass spectrometry. Cytc in this soluble hybrid protein has similar structure and stability as human Cytc and the hybrid protein is endocytosed into a glioma cell line, while displaying potent cytotoxicity and a favorable therapeutic index. Its facile, low-cost, and high yield synthesis, biocompatibility, and robustness suggest that the hybrid protein is a promising candidate for antiglioma drug evaluation.

Phylogenetic Analyses of Chemokine Receptors from Sequence Retrieval to Phylogenetic Trees.

Phylogenetic trees are an essential requisite for comparative biology studies where hypotheses regarding the evolution of genes can be investigated. Trees provide visual and statistical guides to characterize the degree of relatedness among biological entities from genes to species. In a tree, ancestor-descendant relationships are represented by connections, and closely related entities share most of these links. In this chapter, I outlined a method to retrieve and label amino acid and nucleotide sequences of chemokine receptors, align them in sequence matrices, determine their best-model of molecular evolution, and estimate the corresponding phylogenetic trees with distance and maximum likelihood approaches. Most of these analyses are performed within the R environment, and all of these methods use open-source software.

The Characterization of Monoclonal Antibodies to Mouse TLT-1 Suggests That TLT-1 Plays a Role in Wound Healing.

Platelets play a vital role in hemostasis and inflammation. The membrane receptor TREM-like transcript-1 (TLT-1) is involved in platelet aggregation, bleeding, and inflammation, and it is localized in the α-granules of platelets. Upon platelet activation, TLT-1 is released from α-granules both in its transmembrane form and as a soluble fragment (sTLT-1). Higher levels of sTLT-1 have been detected in the plasma of patients with acute inflammation or sepsis, suggesting an important role for TLT-1 during inflammation. However, the roles of TLT-1 in hemostasis and inflammation are not well understood. We are developing the mouse model of TLT-1 to mechanistically test clinical associations of TLT-1 in health and disease. To facilitate our studies, monoclonal murine TLT-1 (mTLT-1) antibodies were produced by the immunization of a rabbit using the negatively charged region of the mTLT-1 extracellular domain 122PPVPGPREGEEAEDEK139. In the present study, we demonstrate that two selected clones, 4.6 and 4.8, are suitable for the detection of mTLT-1 by western blot, immunoprecipitation, immunofluorescent staining, flow cytometry and inhibit platelet aggregation in aggregometry assays. In addition, we found that the topical administration of clone 4.8 delayed the wound healing process in an experimental burn model. These results suggest that TLT-1 plays an important role in wound healing and because both clones specifically detect mTLT-1, they are suitable to further develop TLT-1 based models of inflammation and hemostasis in vivo.

Thermal physiology of Amazonian lizards (Reptilia: Squamata).

We summarize thermal-biology data of 69 species of Amazonian lizards, including mode of thermoregulation and field-active body temperatures (Tb). We also provide new data on preferred temperatures (Tpref), voluntary and thermal-tolerance ranges, and thermal-performance curves (TPC's) for 27 species from nine sites in the Brazilian Amazonia. We tested for phylogenetic signal and pairwise correlations among thermal traits. We found that species generally categorized as thermoregulators have the highest mean values for all thermal traits, and broader ranges for Tb, critical thermal maximum (CTmax) and optimal (Topt) temperatures. Species generally categorized as thermoconformers have large ranges for Tpref, critical thermal minimum (CTmin), and minimum voluntary (VTmin) temperatures for performance. Despite these differences, our results show that all thermal characteristics overlap between both groups and suggest that Amazonian lizards do not fit into discrete thermoregulatory categories. The traits are all correlated, with the exceptions of (1) Topt, which does not correlate with CTmax, and (2) CTmin, and correlates only with Topt. Weak phylogenetic signals for Tb, Tpref and VTmin indicate that these characters may be shaped by local environmental conditions and influenced by phylogeny. We found that open-habitat species perform well under present environmental conditions, without experiencing detectable thermal stress from high environmental temperatures induced in lab experiments. For forest-dwelling lizards, we expect warming trends in Amazonia to induce thermal stress, as temperatures surpass the thermal tolerances for these species.

Diversity within diversity: Parasite species richness in poison frogs assessed by transcriptomics.

Symbionts (e.g., endoparasites and commensals) play an integral role in their host's ecology, yet in many cases their diversity is likely underestimated. Although endoparasites are traditionally characterized using morphology, sequences of conserved genes, and shotgun metagenomics, host transcriptomes constitute an underused resource to identify these organisms' diversity. By isolating non-host transcripts from host transcriptomes, individual host tissues can now simultaneously reveal their endoparasite species richness (i.e., number of different taxa) and provide insights into parasite gene expression. These approaches can be used in host taxa whose endoparasites are mostly unknown, such as those of tropical amphibians. Here, we focus on the poison frogs (Dendrobatidae) as hosts, which are a Neotropical clade known for their bright coloration and defensive alkaloids. These toxins are an effective protection against vertebrate predators (e.g., snakes and birds), bacteria, and skin-biting ectoparasites (e.g., mosquitoes); however, little is known about their deterrence against eukaryotic endoparasites. With de novo transcriptomes of dendrobatids, we developed a bioinformatics pipeline for endoparasite identification that uses host annotated RNA-seq data and set of a priori parasite taxonomic terms, which are used to mine for specific endoparasites. We found a large community of helminths and protozoans that were mostly restricted to the digestive tract and a few systemic parasites (e.g., Trypanosoma). Contrary to our expectations, all dendrobatid frogs regardless of the presence of alkaloid defenses have endoparasites, with their highest species richness located in the frog digestive tract. Some of these organisms (e.g., roundworms) might prove to be generalists, as they were not found to be co-diversifying with their frog hosts. We propose that endoparasites may escape poison frogs' chemical defenses by colonizing tissues with fewer alkaloids than the frog's skin, where most toxins are stored.

Extinction risks forced by climatic change and intraspecific variation in the thermal physiology of a tropical lizard.

Temperature increases can impact biodiversity and predicting their effects is one of the main challenges facing global climate-change research. Ectotherms are sensitive to temperature change and, although predictions indicate that tropical species are highly vulnerable to global warming, they remain one of the least studied groups with respect to the extent of physiological variation and local extinction risks. We model the extinction risks for a tropical heliothermic teiid lizard (Kentropyx calcarata) integrating previously obtained information on intraspecific phylogeographic structure, eco-physiological traits and contemporary species distributions in the Amazon rainforest and its ecotone to the Cerrado savannah. We also investigated how thermal-biology traits vary throughout the species' geographic range and the consequences of such variation for lineage vulnerability. We show substantial variation in thermal tolerance of individuals among thermally distinct sites. Thermal critical limits were highly correlated with operative environmental temperatures. Our physiological/climatic model predicted relative extinction risks for local populations within clades of K. calcarata for 2050 ranging between 26.1% and 70.8%, while for 2070, extinction risks ranged from 52.8% to 92.8%. Our results support the hypothesis that tropical-lizard taxa are at high risk of local extinction caused by increasing temperatures. However, the thermo-physiological differences found across the species' distribution suggest that local adaptation may allow persistence of this tropical ectotherm in global warming scenarios. These results will serve as basis to further research to investigate the strength of local adaptation to climate change. Persistence of Kentropyx calcarata also depends on forest preservation, but the Amazon rainforest is currently under high deforestation rates. We argue that higher conservation priority is necessary so the Amazon rainforest can fulfill its capacity to absorb the impacts of temperature increase on tropical ectotherms during climate change.

Radiation of the polymorphic Little Devil poison frog (Oophaga sylvatica) in Ecuador.

Some South American poison frogs (Dendrobatidae) are chemically defended and use bright aposematic colors to warn potential predators of their unpalatability. Aposematic signals are often frequency-dependent where individuals deviating from a local model are at a higher risk of predation. However, extreme diversity in the aposematic signal has been documented in poison frogs, especially in Oophaga. Here, we explore the phylogeographic pattern among color-divergent populations of the Little Devil poison frog Oophaga sylvatica by analyzing population structure and genetic differentiation to evaluate which processes could account for color diversity within and among populations. With a combination of PCR amplicons (three mitochondrial and three nuclear markers) and genome-wide markers from a double-digested RAD (ddRAD) approach, we characterized the phylogenetic and genetic structure of 199 individuals from 13 populations (12 monomorphic and 1 polymorphic) across the O. sylvatica distribution. Individuals segregated into two main lineages by their northern or southern latitudinal distribution. A high level of genetic and phenotypic polymorphism within the northern lineage suggests ongoing gene flow. In contrast, low levels of genetic differentiation were detected among the southern lineage populations and support recent range expansions from populations in the northern lineage. We propose that a combination of climatic gradients and structured landscapes might be promoting gene flow and phylogenetic diversification. Alternatively, we cannot rule out that the observed phenotypic and genomic variations are the result of genetic drift on near or neutral alleles in a small number of genes.

Sexual dimorphism, phenotypic integration, and the evolution of head structure in casque-headed lizards.

Sexes can differ in features associated with differential reproduction, which can be used during courtship or aggressive encounters. Some traits tend to evolve independently between sexes and emerge as sexually dimorphic within the organismal phenotype. We characterize such a relationship by estimating the phenotypic integration of the head morphology and modularity of the crest in the casque-headed lizards (Corytophanidae). In this clade, some species show extreme sexual dimorphism (e.g., head crests in the genus Basiliscus) while in others, both sexes are monomorphic. To characterize these patterns, we define phenotypic integration at the interspecific level as a pattern or network of traits evidenced by phylogenetically adjusted correlations that persist among species. At this level, modularity is an increased connectedness (e.g., higher correlation) among sections of these networks that persist in a lineage during the evolution of complex phenotypes. To test both concepts, we used phylogenetic geomorphometrics to characterize the head structure of corytophanid lizards, based on a time-calibrated phylogeny that includes candidate fossil ancestors. We found evidence of an older diversification of corytophanids than previously reported (~67 vs. ~23.5 MYA) and show that this clade includes two morphological head architectures: (1) Sexually dimorphic crests present in males that are evolving independently from the rest of the head structure, and (2) full integration of the head morphology in monomorphic species. We propose that both architectures are optimal evolutionary trajectories of the parietal crest bones in the head of these lizards. In sexually dimorphic species, these bones are elongated and thinner, and gave rise to the extended crest used in male courtship displays. In monomorphic species, the parietal crest grew thicker in both sexes to allow for a better insertion of muscles associated with a stronger bite.

Interacting amino acid replacements allow poison frogs to evolve epibatidine resistance.

Animals that wield toxins face self-intoxication. Poison frogs have a diverse arsenal of defensive alkaloids that target the nervous system. Among them is epibatidine, a nicotinic acetylcholine receptor (nAChR) agonist that is lethal at microgram doses. Epibatidine shares a highly conserved binding site with acetylcholine, making it difficult to evolve resistance yet maintain nAChR function. Electrophysiological assays of human and frog nAChR revealed that one amino acid replacement, which evolved three times in poison frogs, decreased epibatidine sensitivity but at a cost of acetylcholine sensitivity. However, receptor functionality was rescued by additional amino acid replacements that differed among poison frog lineages. Our results demonstrate how resistance to agonist toxins can evolve and that such genetic changes propel organisms toward an adaptive peak of chemical defense.