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

Next-generation mitogenomics: A comparison of approaches applied to caecilian amphibian phylogeny

Mitochondrial genome (mitogenome) sequences are being generated with increasing speed due to the advances of next-generation sequencing (NGS) technology and associated analytical tools. However, detailed comparisons to explore the utility of alternative NGS approaches applied to the same taxa have not been undertaken. We compared a 'traditional' Sanger sequencing method with two NGS approaches (shotgun sequencing and non-indexed, multiplex amplicon sequencing) on four different sequencing platforms (Illumina's HiSeq and MiSeq, Roche's 454 GS FLX, and Life Technologies' Ion Torrent) to produce seven (near-) complete mitogenomes from six species that form a small radiation of caecilian amphibians from the Seychelles. The fastest, most accurate method of obtaining mitogenome sequences that we tested was direct sequencing of genomic DNA (shotgun sequencing) using the MiSeq platform. Bayesian inference and maximum likelihood analyses using seven different partitioning strategies were unable to resolve compellingly all phylogenetic relationships among the Seychelles caecilian species, indicating the need for additional data in this case

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

Cranial biomechanics in basal urodeles: the Siberian salamander (Salamandrella keyserlingii) and its evolutionary and developmental implications

Developmental changes in salamander skulls, before and after metamorphosis, afect the feeding capabilities of these animals. How changes in cranial morphology and tissue properties afect the function of the skull are key to decipher the early evolutionary history of the crown-group of salamanders. Here, 3D cranial biomechanics of the adult Salamandrella keyserlingii were analyzed under diferent tissue properties and ossifcation sequences of the cranial skeleton. This helped unravel that: (a) Mechanical properties of tissues (as bone, cartilage or connective tissue) imply a consensus between the stifness required to perform a function versus the fxation (and displacement) required with the surrounding skeletal elements. (b) Changes on the ossifcation pattern, producing fontanelles as a result of bone loss or failure to ossify, represent a trend toward simplifcation potentially helping to distribute stress through the skull, but may also imply a major destabilization of the skull. (c) Bone loss may be originated due to biomechanical optimization and potential reduction of developmental costs. (d) Hynobiids are excellent models for biomechanical reconstruction of extinct early urodeles

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 Development of Three Long Universal Nuclear Protein-Coding Locus Markers and Their Application to Osteichthyan Phylogenetics with Nested PCR

BACKGROUND: Universal nuclear protein-coding locus (NPCL) markers that are applicable across diverse taxa and show good phylogenetic discrimination have broad applications in molecular phylogenetic studies. For example, RAG1, a representative NPCL marker, has been successfully used to make phylogenetic inferences within all major osteichthyan groups. However, such markers with broad working range and high phylogenetic performance are still scarce. It is necessary to develop more universal NPCL markers comparable to RAG1 for osteichthyan phylogenetics. METHODOLOGY/PRINCIPAL FINDINGS: We developed three long universal NPCL markers (>1.6 kb each) based on single-copy nuclear genes (KIAA1239, SACS and TTN) that possess large exons and exhibit the appropriate evolutionary rates. We then compared their phylogenetic utilities with that of the reference marker RAG1 in 47 jawed vertebrate species. In comparison with RAG1, each of the three long universal markers yielded similar topologies and branch supports, all in congruence with the currently accepted osteichthyan phylogeny. To compare their phylogenetic performance visually, we also estimated the phylogenetic informativeness (PI) profile for each of the four long universal NPCL markers. The PI curves indicated that SACS performed best over the whole timescale, while RAG1, KIAA1239 and TTN exhibited similar phylogenetic performances. In addition, we compared the success of nested PCR and standard PCR when amplifying NPCL marker fragments. The amplification success rate and efficiency of the nested PCR were overwhelmingly higher than those of standard PCR. CONCLUSIONS/SIGNIFICANCE: Our work clearly demonstrates the superiority of nested PCR over the conventional PCR in phylogenetic studies and develops three long universal NPCL markers (KIAA1239, SACS and TTN) with the nested PCR strategy. The three markers exhibit high phylogenetic utilities in osteichthyan phylogenetics and can be widely used as pilot genes for phylogenetic questions of osteichthyans at different taxonomic levels

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

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