Molecular evolution of HoxA13 and the multiple origins of limbless morphologies in amphibians and reptiles

Developmental processes and their results, morphological characters, are inherited through transmission of genes regulating development. While there is ample evidence that cis-regulatory elements tend to be modular, with sequence segments dedicated to different roles, the situation for proteins is less clear, being particularly complex for transcription factors with multiple functions. Some motifs mediating protein-protein interactions may be exclusive to particular developmental roles, but it is also possible that motifs are mostly shared among different processes. Here we focus on HoxA13, a protein essential for limb development. We asked whether the HoxA13 amino acid sequence evolved similarly in three limbless clades: Gymnophiona, Amphisbaenia and Serpentes. We explored variation in ω (dN/dS) using a maximum-likelihood framework and HoxA13sequences from 47 species. Comparisons of evolutionary models provided low ω global values and no evidence that HoxA13 experienced relaxed selection in limbless clades. Branch-site models failed to detect evidence for positive selection acting on any site along branches of Amphisbaena and Gymnophiona, while three sites were identified in Serpentes. Examination of alignments did not reveal consistent sequence differences between limbed and limbless species. We conclude that HoxA13 has no modules exclusive to limb development, which may be explained by its involvement in multiple developmental processes

Claims of Potential Expansion throughout the U.S. by Invasive Python Species Are Contradicted by Ecological Niche Models

BACKGROUND: Recent reports from the United States Geological Survey (USGS) suggested that invasive Burmese pythons in the Everglades may quickly spread into many parts of the U.S. due to putative climatic suitability. Additionally, projected trends of global warming were predicted to significantly increase suitable habitat and promote range expansion by these snakes. However, the ecological limitations of the Burmese python are not known and the possible effects of global warming on the potential expansion of the species are also unclear. METHODOLOGY/PRINCIPAL FINDINGS: Here we show that a predicted continental expansion is unlikely based on the ecology of the organism and the climate of the U.S. Our ecological niche models, which include variables representing climatic extremes as well as averages, indicate that the only suitable habitat in the U.S. for Burmese pythons presently occurs in southern Florida and in extreme southern Texas. Models based on the current distribution of the snake predict suitable habitat in essentially the only region in which the snakes are found in the U.S. Future climate models based on global warming forecasts actually indicate a significant contraction in suitable habitat for Burmese pythons in the U.S. as well as in their native range. CONCLUSIONS/SIGNIFICANCE: The Burmese python is strongly limited to the small area of suitable environmental conditions in the United States it currently inhabits due to the ecological niche preferences of the snake. The ability of the Burmese python to expand further into the U.S. is severely limited by ecological constraints. Global warming is predicted to significantly reduce the area of suitable habitat worldwide, underscoring the potential negative effects of climate change for many species

The Molecular Phylogenetic Signature of Clades in Decline

Molecular phylogenies have been used to study the diversification of many clades. However, current methods for inferring diversification dynamics from molecular phylogenies ignore the possibility that clades may be decreasing in diversity, despite the fact that the fossil record shows this to be the case for many groups. Here we investigate the molecular phylogenetic signature of decreasing diversity using the most widely used statistic for inferring diversity dynamics from molecular phylogenies, the γ statistic. We show that if a clade is in decline its molecular phylogeny may show evidence of the decrease in the diversification rate that occurred between its diversification and decline phases. The ability to detect the change in diversification rate depends largely on the ratio of the speciation rates of the diversification and decline phases, the higher the ratio the stronger the signal of the change in diversification rate. Consequently, molecular phylogenies of clades in relative rapid decline do not carry a signature of their decreasing diversification. Further, the signal of the change in diversification rate, if present, declines as the diversity drop. Unfortunately, the molecular signature of clades in decline is the same as the signature produced by diversity dependent diversification. Given this similarity, and the inability of current methods to detect declining diversity, it is likely that some of the extant clades that show a decrease in diversification rate, currently interpreted as evidence for diversity dependent diversification, are in fact in decline. Unless methods can be developed that can discriminate between the different modes of diversification, specifically diversity dependent diversification and declining diversity, we will need the fossil record, or data from some other source, to distinguish between these very different diversity trajectories

Shedding light on the ‘dark side’ of phylogenetic comparative methods

ORCID: 0000-0003-4919-8655© 2016 The Authors. Methods in Ecology and Evolution published by John Wiley & Sons Ltd on behalf of British Ecological Society This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. The attached file is the published version of the article

Macroevolutionary diversification of glands for chemical communication in squamate reptiles

Chemical communication plays a central role in social, sexual and ecological interactions among animals. However, the macroevolutionary diversification of traits responsible for chemical signaling remains fundamentally unknown. Most research investigating evolutionary diversification of glands responsible for the production of chemical signals has focused on arthropods, while its study among vertebrates remains neglected. Using a global-scale dataset covering > 80% (7,904 species) of the living diversity of lizards and snakes (squamates), we investigate rates, trajectories and phylogenetic patterns of diversification of their follicular glands for chemical communication. We observed these glands in 13.66% of species, that their expression has varying phylogenetic signal among lineages, and that the crown squamate ancestor lacked follicular glands, which therefore originated and diversified subsequently during their evolutionary history. Additionally, our findings challenge the longstanding view that within squamates the Iguania are visually oriented while Scleroglossa are chemically-oriented, given that Iguania doubles Scleroglossa in the frequency of glands. Our phylogenetic analyses identified stabilizing selection as the best model describing follicular gland diversification, and revealed high rates of disparity. We provide the first global-scale analysis investigating the diversification of one of the main forms of communication among reptiles, presenting a macroevolutionary angle to questions traditionally explored at microevolutionary scale

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the \u27venom-ome\u27 and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 \u27venom-ome-specific toxins\u27 (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery

Bayesian molecular clock dating of species divergences in the genomics era

It has been five decades since the proposal of the molecular clock hypothesis, which states that the rate of evolution at the molecular level is constant through time and among species. This hypothesis has become a powerful tool in evolutionary biology, making it possible to use molecular sequences to estimate the geological ages of species divergence events. With recent advances in Bayesian clock dating methodology and the explosive accumulation of genetic sequence data, molecular clock dating has found widespread applications, from tracking virus pandemics, to studying the macroevolutionary process of speciation and extinction, to estimating a timescale for Life on Earth

Island survivors: population genetic structure and demography of the critically endangered giant lizard of La Gomera, Gallotia bravoana

Background: The giant lizard of La Gomera (Gallotia bravoana), is an endemic lacertid of this Canary Island that lives confined to a very restricted area of occupancy in a steep cliff, and is catalogued as Critically Endangered by IUCN. We present the first population genetic analysis of the wild population as well as of captive-born individuals (for which paternity data are available) from a recovery center. Current genetic variability, and inferred past demographic changes were determined in order to discern the relative contribution of natural versus human-mediated effects on the observed decline in population size. Results: Genetic analyses indicate that the only known natural population of the species shows low genetic diversity and acts as a single evolutionary unit. Demographic analyses inferred a prolonged decline of the species for at least 230 generations. Depending on the assumed generation time, the onset of the decline was dated between 1200-13000 years ago. Pedigree analyses of captive individuals suggest that reproductive behavior of the giant lizard of La Gomera may include polyandry, multiple paternity and female long-term sperm retention. Conclusions: The current low genetic diversity of G. bravoana is the result of a long-term gradual decline. Because generation time is unknown in this lizard and estimates had large credibility intervals, it is not possible to determine the relative contribution of humans in the collapse of the population. Shorter generation times would favor a stronger influence of human pressure whereas longer generation times would favor a climate-induced origin of the decline. In any case, our analyses show that the wild population has survived for a long period of time with low levels of genetic diversity and a small effective population size. Reproductive behavior may have acted as an important inbreeding avoidance mechanism allowing the species to elude extinction. Overall, our results suggest that the species retains its adaptive potential and could restore its ancient genetic diversity under favorable conditions. Therefore, management of the giant lizard of La Gomera should concentrate efforts on enhancing population growth rates through captive breeding of the species as well as on restoring the carrying capacity of its natural habitat.Spanish Ministry of Education; European Life Project [LIFE 02 NAT-E-008614]; Ministerio de Ciencia e Innovacion [REN 2001- 1514/GLO, CGL 2010-18216]info:eu-repo/semantics/publishedVersio

Phylogeography of Sardinian Cave Salamanders (Genus Hydromantes) Is Mainly Determined by Geomorphology

Detecting the factors that determine the interruption of gene flow between populations is key to understanding how speciation occurs. In this context, caves are an excellent system for studying processes of colonization, differentiation and speciation, since they represent discrete geographical units often with known geological histories. Here, we asked whether discontinuous calcareous areas and cave systems represent major barriers to gene flow within and among the five species of Sardinian cave salamanders (genus Hydromantes) and whether intraspecific genetic structure parallels geographic distance within and among caves. We generated mitochondrial cytochrome b gene sequences from 184 individuals representing 48 populations, and used a Bayesian phylogeographic approach to infer possible areas of cladogenesis for these species and reconstruct historical and current dispersal routes among distinct populations. Our results show deep genetic divergence within and among all Sardinian cave salamander species, which can mostly be attributed to the effects of mountains and discontinuities in major calcareous areas and cave systems acting as barriers to gene flow. While these salamander species can also occur outside caves, our results indicate that there is a very poor dispersal of these species between separate cave systems

Chromosome evolution in Cophomantini (Amphibia, Anura, Hylinae)

The hylid tribe Cophomantini is a diverse clade of Neotropical treefrogs composed of the genera Aplastodiscus, Boana, Bokermannohyla, Hyloscirtus, and Myersiohyla. The phylogenetic relationships of Cophomantini have been comprehensively reviewed in the literature, providing a suitable framework for the study of chromosome evolution. Employing different banding techniques, we studied the chromosomes of 25 species of Boana and 3 of Hyloscirtus; thus providing, for the first time, data for Hyloscirtus and for 15 species of Boana. Most species showed karyotypes with 2n = 2x = 24 chromosomes; some species of the B. albopunctata group have 2n = 2x = 22, and H. alytolylax has 2n = 2x = 20. Karyotypes are all bi-armed in most species presented, with the exception of H. larinopygion (FN = 46) and H. alytolylax (FN = 38), with karyotypes that have a single pair of small telocentric chromosomes. In most species of Boana, NORs are observed in a single pair of chromosomes, mostly in the small chromosomes, although in some species of the B. albopunctata, B. pulchella, and B. semilineata groups, this marker occurs on the larger pairs 8, 1, and 7, respectively. In Hyloscirtus, NOR position differs in the three studied species: H. alytolylax (4p), H. palmeri (4q), and H. larinopygion (1p). Heterochromatin is a variable marker that could provide valuable evidence, but it would be necesserary to understand the molecular composition of the C-bands that are observed in different species in order to test its putative homology. In H. alytolylax, a centromeric DAPI+ band was observed on one homologue of chromosome pair 2. The band was present in males but absent in females, providing evidence for an XX/XY sex determining system in this species. We review and discuss the importance of the different chromosome markers (NOR position, C-bands, and DAPI/CMA3 patterns) for their impact on the taxonomy and karyotype evolution in Cophomantini