Phylogenetic relationships of colubroid snakes based on mitochondrial DNA sequences

A New Species of Pareas (Squamata, Pareidae) from Guangxi Province, China

We described a new species of genus Pareas from Baise City, Guangxi Zhuang Autonomous Region, China, based on morphological and molecular evidence. Pareas baiseensis sp. nov. is distinguished from its congeners by the combination of (1) Yellowish-brown body colouration; (2) Frontal subhexagonal to diamond-shaped with its lateral sides converging posteriorly; (3) The anterior pair of chin shields is longer than it is broad; (4) Loreal not in contact with the eye, prefrontal in contact with the eye, two or three suboculars; (5) Rows of 15-15-15 dorsal scales, five rows of mid-dorsal scales keeled at the middle of the body, one vertebral scale row enlarged; (6) 187-191 ventrals, 89-97 subcaudals, all divided, cloacal plate single; (7) Two postocular stripes, the nuchal area forming a dark black four-pointed fork collar with the middle tines shorter than the outside tines. The genetic divergence (uncorrected p-distance) between the new species and other representatives of Pareas ranged from 13.9% to 24.4% for Cytochrome b (Cyt b) and 12.1% to 25.5% for NADH dehydrogenase subunit 4 (ND4). Phylogenetic analyses of mitochondrial DNA gene data recovered the new species from being the sister taxon to (P. boulengeri + P. chinensis) from China.

Ultraconserved elements-based phylogenomic systematics of the snake superfamily Elapoidea, with the description of a new Afro-Asian family

The highly diverse snake superfamily Elapoidea is considered to be a classic example of ancient, rapid radiation. Such radiations are challenging to fully resolve phylogenetically, with the highly diverse Elapoidea a case in point. Previous attempts at inferring a phylogeny of elapoids produced highly incongruent estimates of their evolutionary relationships, often with very low statistical support. We sought to resolve this situation by sequencing over 4,500 ultraconserved element loci from multiple representatives of every elapoid family/subfamily level taxon and inferring their phylogenetic relationships with multiple methods. Concatenation and multispecies coalescent based species trees yielded largely congruent and well-supported topologies. Hypotheses of a hard polytomy were not retained for any deep branches. Our phylogenies recovered Cyclocoridae and Elapidae as diverging early within Elapoidea. The Afro-Malagasy radiation of elapoid snakes, classified as multiple subfamilies of an inclusive Lamprophiidae by some earlier authors, was found to be monophyletic in all analyses. The genus Micrelaps was consistently recovered as sister to Lamprophiidae. We establish a new family, Micrelapidae fam. nov., for Micrelaps and assign Brachyophis to this family based on cranial osteological synapomorphy. We estimate that Elapoidea originated in the early Eocene and rapidly diversified into all the major lineages during this epoch. Ecological opportunities presented by the post-Cretaceous-Paleogene mass extinction event may have promoted the explosive radiation of elapoid snakes.

Ibrahim A, El-Abd E. Molecular phylogeny of Lytorhynchus diadema (Reptilia, Colubridae) populations in Saudi Arabia

This study presents the molecular phylogenetic relationships among Lytorhynchus diadema (Duméril, Bibron & Duméril, 1854) populations in Saudi Arabia relative to populations from Africa and Asia. This phylogenetic analysis was based on mitochondrial 16S and 12S rRNA partial gene fragments using Neighbor-joining, Maximum Parsimony, and Bayesian methods. The results strongly support the monophyly of Lytorhynchus based on two concatenated genes and the 12S rRNA gene separately. Also, a significant separation is observed between the Arabian samples from Saudi Arabia, Yemen, and Oman, and the African populations from Egypt, Tunisia, and Morocco.

A revised taxonomy of Asian snail-eating snakes Pareas (Squamata, Pareidae): evidence from morphological comparison and molecular phylogeny

The Asian snail-eating snakes Pareas is the largest genus of the family Pareidae (formerly Pareatidae), and widely distributed in Southeast Asia. However, potential diversity remains poorly explored due to their highly conserved morphology and incomplete samples. Here, on basis of more extensive sampling, interspecific phylogenetic relationships of the genus Pareas were reconstructed using two mitochondrial fragments (cyt b and ND4) and two nuclear genes (c-mos and Rag1), and multivariate morphometrics conducted for external morphological data. Both Bayesian Inference and Maximum Likelihood analyses consistently showed that the genus Pareas was comprised of two distinct, monophyletic lineages with moderate to low support values. Based on evidences from molecular phylogeny and morphological data, cryptic diversity of this genus was uncovered and two new species were described. In additional, the validity of P.macularius is confirmed.

Large-scale molecular phylogeny, morphology, divergence-time estimation, and the fossil record of advanced caenophidian snakes (Squamata: Serpentes)

Caenophidian snakes include the file snake genus Acrochordus and advanced colubroidean snakes that radiated mainly during the Neogene. Although caenophidian snakes are a well-supported clade, their inferred affinities, based either on molecular or morphological data, remain poorly known or controversial. Here, we provide an expanded molecular phylogenetic analysis of Caenophidia and use three non-parametric measures of support-Shimodaira-Hasegawa-Like test (SHL), Felsentein (FBP) and transfer (TBE) bootstrap measures-to evaluate the robustness of each clade in the molecular tree. That very different alternative support values are common suggests that results based on only one support value should be viewed with caution. Using a scheme to combine support values, we find 20.9% of the 1265 clades comprising the inferred caenophidian tree are unambiguously supported by both SHL and FBP values, while almost 37% are unsupported or ambiguously supported, revealing the substantial extent of phylogenetic problems within Caenophidia. Combined FBP/TBE support values show similar results, while SHL/TBE result in slightly higher combined values. We consider key morphological attributes of colubroidean cranial, vertebral and hemipenial anatomy and provide additional morphological evidence supporting the clades Colubroides, Colubriformes, and Endoglyptodonta. We review and revise the relevant caenophidian fossil record and provide a time-calibrated tree derived from our molecular data to discuss the main cladogenetic events that resulted in present-day patterns of caenophidian diversification. Our results suggest that all extant families of Colubroidea and Elapoidea composing the present-day endoglyptodont fauna originated rapidly within the early Oligocene-between approximately 33 and 28 Mya-following the major terrestrial faunal turnover known as the "Grande Coupure" and associated with the overall climate shift at the Eocene-Oligocene boundary. Our results further suggest that the caenophidian radiation originated within the Caenozoic, with the divergence between Colubroides and Acrochordidae occurring in the early Eocene, at ~ 56 Mya.

More evolution underground: Accelerated mitochondrial substitution rate in Australian burrowing freshwater crayfishes (Decapoda: Parastacidae)

To further understand the evolutionary history and mitogenomic features of Australia's highly distinctive freshwater crayfish fauna, we utilized a recently described rapid mitogenome sequencing pipeline to generate 24 new crayfish mitogenomes including a diversity of burrowing crayfish species and the first for Astacopsis gouldi, the world's largest freshwater invertebrate. Whole mitogenome-based phylogeny estimates using both Bayesian and Maximum Likelihood methods substantially strengthen existing hypotheses for systematic relationships among Australian freshwater crayfish with evidence of pervasive diversifying selection and accelerated mitochondrial substitution rate among the members of the clade representing strongly burrowing crayfish that may reflect selection pressures for increased energy requirement for adaptation to terrestrial environment and a burrowing lifestyle. Further, gene rearrangements are prevalent in the burrowing crayfish mitogenomes involving both tRNA and protein coding genes. In addition, duplicated control regions were observed in two closely related Engaeus species, together with evidence for concerted evolution. This study significantly adds to the understanding of Australian freshwater crayfish evolutionary relationships and suggests a link between mitogenome evolution and adaptation to terrestrial environments and a burrowing lifestyle in freshwater crayfish.

Colubrid Venom Composition: An -Omics Perspective

Snake venoms have been subjected to increasingly sensitive analyses for well over 100 years, but most research has been restricted to front-fanged snakes, which actually represent a relatively small proportion of extant species of advanced snakes. Because rear-fanged snakes are a diverse and distinct radiation of the advanced snakes, understanding venom composition among "colubrids" is critical to understanding the evolution of venom among snakes. Here we review the state of knowledge concerning rear-fanged snake venom composition, emphasizing those toxins for which protein or transcript sequences are available. We have also added new transcriptome-based data on venoms of three species of rear-fanged snakes. Based on this compilation, it is apparent that several components, including cysteine-rich secretory proteins (CRiSPs), C-type lectins (CTLs), CTLs-like proteins and snake venom metalloproteinases (SVMPs), are broadly distributed among "colubrid" venoms, while others, notably three-finger toxins (3FTxs), appear nearly restricted to the Colubridae (sensu stricto). Some putative new toxins, such as snake venom matrix metalloproteinases, are in fact present in several colubrid venoms, while others are only transcribed, at lower levels. This work provides insights into the evolution of these toxin classes, but because only a small number of species have been explored, generalizations are still rather limited. It is likely that new venom protein families await discovery, particularly among those species with highly specialized diets.

A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes

BACKGROUND

The extant squamates (>9400 known species of lizards and snakes) are one of the most diverse and conspicuous radiations of terrestrial vertebrates, but no studies have attempted to reconstruct a phylogeny for the group with large-scale taxon sampling. Such an estimate is invaluable for comparative evolutionary studies, and to address their classification. Here, we present the first large-scale phylogenetic estimate for Squamata.

RESULTS

The estimated phylogeny contains 4161 species, representing all currently recognized families and subfamilies. The analysis is based on up to 12896 base pairs of sequence data per species (average = 2497 bp) from 12 genes, including seven nuclear loci (BDNF, c-mos, NT3, PDC, R35, RAG-1, and RAG-2), and five mitochondrial genes (12S, 16S, cytochrome b, ND2, and ND4). The tree provides important confirmation for recent estimates of higher-level squamate phylogeny based on molecular data (but with more limited taxon sampling), estimates that are very different from previous morphology-based hypotheses. The tree also includes many relationships that differ from previous molecular estimates and many that differ from traditional taxonomy.

CONCLUSIONS

We present a new large-scale phylogeny of squamate reptiles that should be a valuable resource for future comparative studies. We also present a revised classification of squamates at the family and subfamily level to bring the taxonomy more in line with the new phylogenetic hypothesis. This classification includes new, resurrected, and modified subfamilies within gymnophthalmid and scincid lizards, and boid, colubrid, and lamprophiid snakes.

What are the closest relatives of the hot-spring snakes (Colubridae, Thermophis), the relict species endemic to the Tibetan Plateau?

We conducted phylogenetic analyses to identify the closest related living relatives of the Xizang and Sichuan hot-spring snakes (T. baileyi and T. zhaoermii) endemic to the Tibetan Plateau, using mitochondrial DNA sequences (cyt b, ND4) from eight specimens, together with sequences from 95 additional caenophidian and five henophidian genera that were downloaded from GenBank. Phylogenetic trees were obtained using Bayesian Inference and Maximum likelihood methods. Results suggest that hot-spring snakes, which are adapted to high and cold environments, were clustered in the monophyletic Xenodontinae. Xenodontinae is one of the largest subfamilies of colubrid snakes, with about 90 genera and more than 500 species known, and are primarily tropical snakes previously thought to be restricted to the New World. Our data failed to provide any evidence that the New World xenodontines diverged from Thermophis and dispersed into the New World, also failed to suggest a colonization of Asia by New World xenodontines by dispersal from the New World. An alternative plausible scenario may be that Thermophis and the New World xenodontines evolved independently in Asia and America, respectively, after the divergence of their common ancestor. The divergence of the two species in Thermophis was caused by the barrier of the Hengduan Mountains, and the speciation had almost occurred when Tibetan Plateau attained present elevation.

The evolution of venom-conducting fangs: Insights from developmental biology

The present study of the origin of the various types of fang represented among colubroid snakes (i.e., tubular, grooved, and ungrooved) attempts to reconcile the morphology of adult fangs with current phylogenetic hypotheses. Observations of growth series of developing tubular fangs were hypothesised to shed light on the evolutionary origin of fangs in snakes. While molecular phylogenies and evolutionary studies of venom proteins and of other anatomical components of the venom-delivery system reconstruct a consistent evolutionary scenario, the character of a tubular venom-conducting fang does not fit in this scenario. The present review offers a series of possible scenarios to resolve this anomaly. Of these, a new idea argues that a heterochronic mechanism (alteration of the timing of developmental events) may provide the answer that the ungrooved and grooved teeth of colubrid snakes evolved from an ancestral tubular fang by means of attachment of replacement tubular fangs to the maxilla at an earlier developmental stage than usual (precocial ankylosis).

Isolation and characterization of rufoxin, a novel protein exhibiting neurotoxicity from venom of the psammophiine, Rhamphiophis oxyrhynchus (Rufous beaked snake)

Colubrid snake venoms potentially represent a vast source of novel biological actives and structural motifs owing to their diverse phylogeny. The present study describes the identification of rufoxin, a neurotoxin from the venom of Rhamphiophis oxyrhynchus (Rufous beaked snake) which is a member of the African colubrid lineage, the psammophiines. Rufoxin (1 microM) displayed reversible post-synaptic neurotoxic activity as evidenced by significant inhibition of indirect twitches and responses to exogenous nicotinic agonists in the chick biventer cervicis nerve-muscle preparation. Rufoxin (0.1-1.0 microM) also caused a rightward parallel shift of cumulative concentration-response curves to carbachol (CCh; 0.6-80 microM) without a significant depression of the maximum response, suggestive of classical competitive antagonism at the skeletal muscle nicotinic receptor. Rufoxin lacks NH(2)-terminal sequence homology to previously identified snake venom toxins. This work indicates a wider distribution of neurotoxins across the advanced snake superfamily than previously described.

The phylogeny and classification of caenophidian snakes inferred from seven nuclear protein-coding genes

More than 80% of the approximately 3000 living species of snakes are placed in the taxon Caenophidia (advanced snakes), a group that includes the families Acrochordidae, Viperidae, Elapidae, Atractaspididae, and the paraphyletic 'Colubridae'. Previous studies using DNA sequences have involved few nuclear genes (one or two). Several nodes have therefore proven difficult to resolve with statistical significance. Here, we investigated the higher-level relationships of caenophidian snakes with seven nuclear protein-coding genes and obtained a well-supported topology. Accordingly, some adjustments to the current classification of Caenophidia are made to better reflect the relationships of the groups. The phylogeny also indicates that, ancestrally, caenophidian snakes are Asian and nocturnal in origin, although living species occur on nearly all continents and are ecologically diverse.

Relict snakes of North America and their relationships within Caenophidia, using likelihood-based Bayesian methods on mitochondrial sequences

This paper focuses on the phylogenetic relationships of eight North American caenophidian snake species (Carphophis amoena, Contia tenuis, Diadophis punctatus, Farancia abacura, Farancia erytrogramma, Heterodon nasicus, Heterodon platyrhinos, and Heterodon simus) whose phylogenetic relationships remain controversial. Past studies have referred to these "relict" North American snakes either as colubrid, or as Neotropical dipsadids and/or xenodontids. Based on mitochondrial DNA ribosomal gene sequences and a likelihood-based Bayesian analysis, our study suggests that these North American snakes are not monophyletic and are nested within a group (Dipsadoidea) that contains the Dipsadidae, Xenodontidae, and Natricidae. In addition, we use the relationships proposed here to highlight putative examples of parallel evolution of hemipenial morphology among snake clades.

Multiple colonization of Madagascar and Socotra by colubrid snakes: evidence from nuclear and mitochondrial gene phylogenies

Colubrid snakes form a speciose group of unclarified phylogeny. Their almost cosmopolitan distribution could be interpreted as a product of plate-tectonic vicariance. We used sequences of the nuclear c-mos, the mitochondrial cytochrome b and the 16S rRNA genes in 41 taxa to elucidate the relationships between the endemic colubrid genera found in Madagascar and in the Socotra archipelago. The well-resolved trees indicate multiple origins of both the Malagasy and the Socotran taxa. The Malagasy genus Mimophis was nested within the Psammophiinae, and the Socotran Hemerophis was closely related to Old World representatives of the former genus Coluber. The remaining 14 genera of Malagasy colubrids formed a monophyletic sister group of the Socotran Ditypophis (together forming the Pseudoxyrhophiinae). Molecular-clock estimates place the divergence of Malagasy and Socotran colubrids from their non-insular sister groups into a time-frame between the Eocene and Miocene. Over-seas rafting is the most likely hypothesis for the origin of at least the Malagasy taxa. The discovery of a large monophyletic clade of colubrids endemic to Madagascar indicates a need for taxonomic changes. The relationship of this radiation to the Socotran Ditypophis highlights the potential of the Indian Ocean islands to act as an evolutionary reservoir for lineages that have become extinct in Africa and Asia.

Analysis of Colubroidea snake venoms by liquid chromatography with mass spectrometry: evolutionary and toxinological implications

The evolution of the venomous function of snakes and the diversification of the toxins has been of tremendous research interest and considerable debate. It has become recently evident that the evolution of the toxins in the advanced snakes (Colubroidea) predated the evolution of the advanced, front-fanged delivery mechanisms. Historically, the venoms of snakes lacking front-fanged venom-delivery systems (conventionally grouped into the paraphyletic family Colubridae) have been largely neglected. In this study we used liquid chromatography with mass spectrometry (LC/MS) to analyze a large number of venoms from a wide array of species representing the major advanced snake clades Atractaspididae, Colubrinae, Elapidae, Homalopsinae, Natricinae, Psammophiinae, Pseudoxyrhophiinae, Xenodontinae, and Viperidae. We also present the first sequences of toxins from Azemiops feae as well as additional toxin sequences from the Colubrinae. The large body of data on molecular masses and retention times thus assembled demonstrates a hitherto unsuspected diversity of toxins in all lineages, having implications ranging from clinical management of envenomings to venom evolution to the use of isolated toxins as leads for drug design and development. Although definitive assignment of a toxin to a protein family can only be done through demonstrated structural studies such as N-terminal sequencing, the molecular mass data complemented by LC retention information, presented here, do permit formulation of reasonable hypotheses concerning snake venom evolution and potential clinical effects to a degree not possible till now, and some hypotheses of this kind are proposed here. The data will also be useful in biodiscovery.

Higher-level relationships of caenophidian snakes inferred from four nuclear and mitochondrial genes

Higher-level caenophidian snake relationships are inferred from sequence analyses of one nuclear gene (C-mos) and three mitochondrial genes (12S rRNA, 16S rRNA and ND4). Caenophidians, which are haenophidian closest relatives, have an Asiatic origin. An African clade comprising atractaspidids, psammophiines, 'lamprophiines' and 'pseudoxyrhophiines' is identified. We discern no evolutionary trend such as an improvement of the venom apparatus with a linear progression from the absence of a venom system to the presence of a front-fanged one. The venom apparatus is contemporary with the origin of colubroids and its absence in a few lineages results from secondary losses. The front-fanged venom system appeared three times independently. The active diurnal foraging mode (associated with a high metabolic rate) appears in a derived position among colubroids.

Molecular clocks in reptiles: life history influences rate of molecular evolution

Life history has been implicated as a determinant of variation in rate of molecular evolution amongst vertebrate species because of a negative correlation between body size and substitution rate for many molecular data sets. Both the generality and the cause of the negative body size trend have been debated, and the validity of key studies has been questioned (particularly concerning the failure to account for phylogenetic bias). In this study, a comparative method has been used to test for an association between a range of life-history variables-such as body size, age at maturity, and clutch size-and DNA substitution rate for three genes (NADH4, cytochrome b, and c-mos). A negative relationship between body size and rate of molecular evolution was found for phylogenetically independent pairs of reptile species spanning turtles, lizards, snakes, crocodile, and tuatara. Although this study was limited by the number of comparisons for which both sequence and life-history data were available, the results suggest that a negative body size trend in rate of molecular evolution may be a general feature of reptile molecular evolution, consistent with similar studies of mammals and birds. This observation has important implications for uncovering the mechanisms of molecular evolution and warns against assuming that related lineages will share the same substitution rate (a local molecular clock) in order to date evolutionary divergences from DNA sequences.

A molecular phylogeny for the frog genus Limnodynastes (Anura: myobatrachidae)

Members of the genus Limnodynastes are a prominent and widespread feature of the Australian frog fauna. Yet despite their potential to be informative about biogeographic history and mechanisms of speciation, the relationships among these taxa are not well known. We investigated phylogenetic relationships within the genus Limnodynastes via sequencing of mitochondrial (mt)DNA from current members of the genus Limnodynastes and the monotypic genus Megistolotis. a 450-bp fragment of the 16S rRNA gene and a 370-bp fragment of the protein-coding gene ND4 were used to infer a molecular phylogeny. We revise traditional species groupings and now recognize four species groups within Limnodynastes: the L. ornatus group (L. ornatus and L. spenceri), the L. peronii group (L. peronii, L. tasmaniensis, L. fletcheri, the L. depressus), the L. salmini group (L. salmini, L. convexiusculus, and L. lignarius), and the L. dorsalis group (L. dorsalis, L. terraereginae, L. dumerilii and L. interioris). The L. ornatus species group forms a highly distinctive clade that is a sister group to the other Limnodynastes groups. Pending broader phylogenetic studies it could be removed from the genus Limnodynastes. Our results concur with previous suggestions that Megistolotis lignarius is nested within Limnodynastes, and we therefore reclassify this species as Limnodynastes lignarius. Furthermore, specimens identified as L. depressus form a mtDNA lineage distinct from other species in the genus, confirming the validity of the species. Specimens of species from the L. dorsalis group (L. dorsalis, L. dumerilii, L. interioris, and L. terraereginae) are closely related such that L. dumerilii is paraphyletic with two other species. Finally, our study provides broad support for previous phylogenies based on microcomplement fixation.