Phylogenetic relationships of elapid snakes based on cytochrome b mtDNA sequences

Published molecular phylogenetic studies of elapid snakes agree that the marine and Australo-Melanesian forms are collectively monophyletic. Recent studies, however, disagree on the relationships of the African, American, and Asian forms. To resolve the relationships of the African, American, and Asian species to each other and to the marine/Australo-Melanesian clade, we sequenced the entire cytochrome b gene for 28 elapids; 2 additional elapid sequences from GenBank were also included. This sample includes all African, American, and Asian genera (except for the rare African Pseudohaje), as well as a representative sample of marine/Australo-Melanesian genera. The data were analyzed by the methods of maximum-parsimony and maximum-likelihood. Both types of analyses yielded similar trees, from which the following conclusions can be drawn: (1) Homoroselaps falls outside a clade formed by the remaining elapids; (2) the remaining elapids are divisible into two broad sister clades, the marine/Australo-Melanesian species vs the African, American, and Asian species; (3) American coral snakes cluster with Asian coral snakes; and (4) the "true" cobra genus Naja is probably not monophyletic as the result of excluding such genera as Boulengerina and Paranaja.     

An immunological assessment of the phylogenetic position of New World coral snakes

The New World coral snakes (micrurines), genera Micrurus and Micruroides have recently been seen as derived from a lineage of South American colubrids, rather than from a common lineage with Old World elapids and sea snakes as traditionally accepted. We compared serum albumins of representative coral snakes, Old World elapids, sea snakes, and neotropical colubrids immunologically. Phylogenetic analysis of the biochemical data unambiguously allies the micrurines with the family Elapidae as it is currently understood. Using the albumin molecular clock calibration derived from other terrestrial vertebrates. we suggest a late Oligocene-early Miocene separation between the New and Old World elapid lineages. This requires a movement of elapid stocks from Asia into North America, and supporting evidence for this model is derived from several paleontological sources. We suggest that a number of extant micrurine lineages have had long independent histories.     

Fingerprint correspondence of hemoglobins and the relationships of sea snakes.

1. Peptide fingerprints of tryptic digests of the globins of sea snake species of Hydrophis, Pelamis, Aipysurus, Laticauda and the terrestrial elapid Naja were compared. 2. Globin divergence, as estimated from peptide fingerprints, paralleled closely transferrin divergence, as measured immunologically. 3. Taxonomic affinities, suggested by the fingerprint data, are concordant with McDowell's taxonomic system for sea snakes with the following exceptions: (a) Laticauda shows a closer affinity to the true sea snakes than to the terrestrial elapid Naja. (b) Sea snakes appear to be more widely divergent from terrestrial elapids than his scheme suggests.     

Venom glands and some associated muscles in sea snakes

The venom glands and related muscles of sea snakes conform in their general structure to those of the terrestrial elapids. The venom gland, however, is smaller in size and the accessory gland is considerably reduced. A similar pattern is found in the Australian elapid Notechis. The musculus compressor glandulae is well developed in the sea snakes and in some species its posterior-medial portion runs uninterruptedly from the origin to the insertion of the muscle. This might be considered as a primitive condition suggesting an early divergence of the sea snakes from an ancestral elapid stock. Three species of sea snakes, Aipysurus eydouxi, Emydocephalus annulatus, and E. ijimae, feed on fish eggs and have very small, but still functioning, venom glands. The reduced accessory gland of the sea snakes is apparently connected with their aquatic environment, as a similar condition is found also in the elapine Boulengerina annulata which lives in large lakes of Central Africa. The similarity in structure of the venom gland between sea snakes and Notechis scutatus may point to a possible phylogenetic relationship between this group of Australian elapids and hydrophiine snakes.     

Phylogeny of Australasian venomous snakes (Colubroidea, Elapidae, Hydrophiinae) based on phenotypic and molecular evidence

Scanlon, John D. & Lee, Michael S. Y. (2004). Phylogeny of Australasian venomous snakes (Colubroidea, Elapidae, Hydrophiinae) based on phenotypic and molecular evidence. — Zoologica Scripta , 33 , 335–366.
Phylogenetic relationships among Hydrophiinae (Australasian and marine elapid snakes) are inferred using 87 characters from external, skeletal, hemipenial and internal anatomy, ecology, and chromosomes as well as available sequences of two mitochondrial genes (cytochrome b and 16S rRNA). Parsimony analysis of the combined data retrieves many widely accepted clades; while observed bootstrap or branch (Bremer) support for these is often weak, most have never been corroborated previously by a rigorous numerical analysis. Sea kraits ( Laticauda ) and Solomon Islands elapids are basal to the remaining hydrophiines (Australian terrestrial forms and hydrophiin sea snakes). The latter clade includes three main lineages: a large-bodied oviparous lineage, a small-bodied oviparous lineage, and a viviparous lineage (which also includes the hydrophiin sea snakes, strongly reaffirmed as monophyletic). While the Solomons retain a relictual fauna, New Guinea has less endemism and has been invaded multiple times by Australian lineages, so there is no clear 'stepping stone' pattern supporting a northern (Asian, rather than Gondwanan) biogeographical origin.     

Molecular phylogeny and divergence dates for Australasian elapids and sea snakes (hydrophiinae): evidence from seven genes for rapid evolutionary radiations

One of the most prolific radiations of venomous snakes, the Australo-Melanesian Hydrophiinae includes approximately 100 species of Australasian terrestrial elapids plus all approximately 60 species of viviparous sea snakes. Here, we estimate hydrophiine relationships based on a large data set comprising 5800 bp drawn from seven genes (mitochondrial: ND4, cytb, 12S, 16S; nuclear: rag1, cmos, myh). These data were analysed using parsimony, likelihood and Bayesian methods to better resolve hydrophiine phylogeny and provide a timescale for the terrestrial and marine radiations. Among oviparous forms, Cacophis, Furina and Demansia are basal to other Australian elapids (core oxyuranines). The Melanesian Toxicocalamus and Aspidomorphus group with Demansia, indicating multiple dispersal events between New Guinea and Australia. Oxyuranus and Pseudonaja form a robust clade. The small burrowing taxa form two separate clades, one consisting of Vermicella and Neelaps calanotus, and the other including Simoselaps, Brachyurophis and Neelaps bimaculatus. The viviparous terrestrial elapids form three separate groups: Acanthophis, the Rhinoplocephalus group and the Notechis-Hemiaspis group. True sea snakes (Hydrophiini) are robustly united with the Notechis-Hemiaspis group. Many of the retrieved groupings are consistent with previous molecular and morphological analyses, but the polyphyly of the viviparous and burrowing groups, and of Neelaps, are novel results. Bayesian relaxed clock analyses indicate very recent divergences: the approximately 160 species of the core Australian radiation (including sea snakes) arose within the last 10 Myr, with most inter-generic splits dating to between 10 and 6 Ma. The Hydrophis sea snake lineage is an exceptionally rapid radiation, with > 40 species evolving within the last 5 Myr.     

Ecological and historical legacies on global diversity gradients in marine elapid snakes

Global diversity gradients have been extensively investigated for several biological groups. However, little is known whether the diversity drivers of clades that underwent major environmental transition (e.g., from land to sea) are equivalent across these different environmental settings. Here, we ask if the pattern of diversity of marine elapid snakes is determined by factors analogous to those previously found for terrestrial lineages. Through a model selection framework, we compare the effect of factors that represent five ecological and historical hypotheses. We found that both ecological and historical factors play significant roles, but habitat structure, which can be linked to historical climatic changes, was better supported. This result agrees with that previously found for terrestrial elapids, despite the different environmental pressures in terrestrial and marine contexts. Our findings suggest an equivalence of the underlying process of diversity gradient within the same lineage from land and sea, irrespective of the different physiological, spatial and historical constraints.     

Inferring Species Trees from Gene Trees: A Phylogenetic Analysis of the Elapidae (Serpentes) Based on the Amino Acid Sequences of Venom Proteins

Toward the goal of recovering the phylogenetic relationships among elapid snakes, we separately found the shortest trees from the amino acid sequences for the venom proteins phospholipase A₂and the short neurotoxin, collectively representing 32 species in 16 genera. We then applied a method we term gene tree parsimony for inferring species trees from gene trees that works by finding the species tree which minimizes the number of deep coalescences or gene duplications plus unsampled sequences necessary to fit each gene tree to the species tree. This procedure, which is both logical and generally applicable, avoids many of the problems of previous approaches for inferring species trees from gene trees. The results support a division of the elapids examined into sister groups of the Australian and marine (laticaudines and hydrophiines) species, and the African and Asian species. Within the former clade, the sea snakes are shown to be diphyletic, with the laticaudines and hydrophiines having separate origins. This finding is corroborated by previous studies, which provide support for the usefulness of gene tree parsimony.     

Moving in two worlds: aquatic and terrestrial locomotion in sea snakes (Laticauda colubrina,Laticaudidae)

Yellow‐lipped sea kraits (Laticauda colubrina) are amphibious in their habits. We measured their locomotor speeds in water and on land to investigate two topics: (1) to what degree have adaptations to increase swimming speed (paddle‐like tail etc.) reduced terrestrial locomotor ability in sea kraits?; and (2) do a sea krait’s sex and body size influence its locomotor ability in these two habitats, as might be expected from the fact that different age and sex classes of sea kraits use the marine and terrestrial environments in different ways? To estimate ancestral states for locomotor performance, we measured speeds of three species of Australian terrestrial elapids that spend part of their time foraging in water. The evolutionary modifications of Laticauda for marine life have enhanced their swimming speeds by about 60%, but decreased their terrestrial locomotor speed by about 80%. Larger snakes moved faster than smaller individuals in absolute terms but were slower in terms of body lengths travelled per second, especially on land. Male sea kraits were faster than females (independent of the body‐size effect), especially on land. Prey items in the gut reduced locomotor speeds both on land and in water. Proteroglyphous snakes may offer exceptional opportunities to study phylogenetic shifts in locomotor ability, because (1) they display multiple independent evolutionary shifts from terrestrial to aquatic habits, and (2) one proteroglyph lineage (the laticaudids) displays considerable intraspecific and interspecific diversity in terms of the degree to which they use terrestrial vs. aquatic habitats.     

A Simple and Novel Strategy for the Production of a Pan-specific Antiserum against Elapid Snakes of Asia

Snakebite envenomation is a serious medical problem in many tropical developing countries and was considered by WHO as a neglected tropical disease. Antivenom (AV), the rational and most effective treatment modality, is either unaffordable and/or unavailable in many affected countries. Moreover, each AV is specific to only one (monospecific) or a few (polyspecific) snake venoms. This demands that each country to prepare AV against its local snake venoms, which is often not feasible. Preparation of a ‘pan-specific’ AV against many snakes over a wide geographical area in some countries/regions has not been possible. If a ‘pan-specific’ AV effective against a variety of snakes from many countries could be prepared, it could be produced economically in large volume for use in many countries and save many lives. The aim of this study was to produce a pan-specific antiserum effective against major medically important elapids in Asia. The strategy was to use toxin fractions (TFs) of the venoms in place of crude venoms in order to reduce the number of antigens the horses were exposed to. This enabled inclusion of a greater variety of elapid venoms in the immunogen mix, thus exposing the horse immune system to a diverse repertoire of toxin epitopes, and gave rise to antiserum with wide paraspecificity against elapid venoms. Twelve venom samples from six medically important elapid snakes (4 Naja spp. and 2 Bungarus spp.) were collected from 12 regions/countries in Asia. Nine of these 12 venoms were ultra-filtered to remove high molecular weight, non-toxic and highly immunogenic proteins. The remaining 3 venoms were not ultra-filtered due to limited amounts available. The 9 toxin fractions (TFs) together with the 3 crude venoms were emulsified in complete Freund’s adjuvant and used to immunize 3 horses using a low dose, low volume, multisite immunization protocol. The horse antisera were assayed by ELISA and by in vivo lethality neutralization in mice. The findings were: a) The 9 TFs were shown to contain all of the venom toxins but were devoid of high MW proteins. When these TFs, together with the 3 crude venoms, were used as the immunogen, satisfactory ELISA antibody titers against homologous/heterologous venoms were obtained. b) The horse antiserum immunologically reacted with and neutralized the lethal effects of both the homologous and the 16 heterologous Asian/African elapid venoms tested. Thus, the use of TFs in place of crude venoms and the inclusion of a variety of elapid venoms in the immunogen mix resulted in antiserum with wide paraspecificity against elapid venoms from distant geographic areas. The antivenom prepared from this antiserum would be expected to be pan-specific and effective in treating envenomations by most elapids in many Asian countries. Due to economies of scale, the antivenom could be produced inexpensively and save many lives. This simple strategy and procedure could be readily adapted for the production of pan-specific antisera against elapids of other continents.     

Higher-level snake phylogeny inferred from mitochondrial DNA sequences of 12S rRNA and 16S rRNA genes

Portions of two mitochondrial genes (12S and 16S ribosomal RNA) were sequenced to determine the phylogenetic relationships among the major clades of snakes. Thirty-six species, representing nearly all extant families, were examined and compared with sequences of a tuatara and three families of lizards. Snakes were found to constitute a monophyletic group (confidence probability [CP] = 96%), with the scolecophidians (blind snakes) as the most basal lineages (CP = 99%). This finding supports the hypothesis that snakes underwent a subterranean period early in their evolution. Caenophidians (advanced snakes), excluding Acrochordus, were found to be monophyletic (CP = 99%). Among the caenophidians, viperids were monophyletic (CP = 98%) and formed the sister group to the elapids plus colubrids (CP = 94%). Within the viperids, two monophyletic groups were identified: true vipers (CP = 98%) and pit vipers plus Azemiops (CP = 99%). The elapids plus Atractaspis formed a monophyletic clade (CP = 99%). Within the paraphyletic Colubridae, the largely Holarctic Colubrinae was found to be a monophyletic assemblage (CP = 98%), and the Xenodontinae was found to be polyphyletic (CP = 91%). Monophyly of the henophidians (primitive snakes) was neither supported nor rejected because of the weak resolution of relationships among those taxa, except for the clustering of Calabaria with a uropeltid, Rhinophis (CP = 94%).      

Phylogenetic relationships within laticaudine sea snakes (Elapidae)

The sea snake subfamily Laticaudinae consists of a single genus with eight named species, based on morphological characters. We used microsatellite and mitochondrial DNA (mtDNA) data to clarify the adaptive radiation of these oviparous sea snakes in the South Pacific, with special reference to New Caledonia and Vanuatu. A mitochondrial DNA data set (ND4 gene 793 bp) was obtained from 345 individuals of the five species of Laticauda sp. sea snakes endemic to the region. Maximum likelihood and Bayesian approaches yielded the same optimal tree topology, identifying two major clades (yellow-banded and blue-banded sea snakes). Although all laticaudine sea snakes rely on small islands as oviposition sites, the two lineages differ in their use of marine vs. terrestrial habitats. A highly aquatic species (Laticauda laticaudata) shows a strong pattern of genetic isolation by distance, implying that the patchy distribution of terrestrial habitats has had little impact on gene flow. The more terrestrial clade (Laticauda colubrina, Laticauda frontalis, Laticauda guineai, Laticauda saintgironsi) shows stronger geographic differentiation in allelic frequencies, associated with island groups rather than with geographic distance. Microsatellites and mtDNA suggest that L. frontalis (restricted to Vanuatu) represents a recent founder-induced speciation event, from allopatric migrants of the New Caledonian taxon L. saintgironsi. A major divergence in speciation patterns between the two major clades of laticaudine snakes thus correlates with (and perhaps, is driven by) differences in the importance of terrestrial habitats in the species' ecology.     

Ecophysiological steps of marine adaptation in extant and extinct non-avian tetrapods

Marine reptiles and mammals are phylogenetically so distant from each other that their marine adaptations are rarely compared directly. We reviewed ecophysiological features in extant non-avian marine tetrapods representing 31 marine colonizations to test whether there is a common pattern across higher taxonomic groups, such as mammals and reptiles. Marine adaptations in tetrapods can be roughly divided into aquatic and haline adaptations, each of which seems to follow a sequence of three steps. In combination, these six categories exhibit five steps of marine adaptation that apply across all clades except snakes: Step M1, incipient use of marine resources; Step M2, direct feeding in the saline sea; Step M3, water balance maintenance without terrestrial fresh water; Step M4, minimized terrestrial travel and loss of terrestrial feeding; and Step M5, loss of terrestrial thermoregulation and fur/plumage. Acquisition of viviparity is not included because there is no known case where viviparity evolved after a tetrapod lineage colonized the sea. A similar sequence is found in snakes but with the haline adaptation step (Step M3) lagging behind aquatic adaptation (haline adaptation is Step S5 in snakes), most likely because their unique method of water balance maintenance requires a supply of fresh water. The same constraint may limit the maximum body size of fully marine snakes. Steps M4 and M5 in all taxa except snakes are associated with skeletal adaptations that are mechanistically linked to relevant ecophysiological features, allowing assessment of marine adaptation steps in some fossil marine tetrapods. We identified four fossil clades containing members that reached Step M5 outside of stem whales, pinnipeds, sea cows and sea turtles, namely Eosauropterygia, Ichthyosauromorpha, Mosasauroidea, and Thalattosuchia, while five other clades reached Step M4: Saurosphargidae, Placodontia, Dinocephalosaurus, Desmostylia, and Odontochelys. Clades reaching Steps M4 and M5, both extant and extinct, appear to have higher species diversity than those only reaching Steps M1 to M3, while the total number of clades is higher for the earlier steps. This suggests that marine colonizers only diversified greatly after they minimized their use of terrestrial resources, with many lineages not reaching these advanced steps. Historical patterns suggest that a clade does not advance to Steps M4 and M5 unless these steps are reached early in the evolution of the clade. Intermediate forms before a clade reached Steps M4 and M5 tend to become extinct without leaving extant descendants or fossil evidence. This makes it difficult to reconstruct the evolutionary history of marine adaptation in many clades. Clades that reached Steps M4 and M5 tend to last longer than other marine tetrapod clades, sometimes for more than 100 million years.     

Physical properties of some snake plasma albumins

The albumin-like proteins were purified from the plasma of three terrestrial elapids and two sea snakes. The albumin-like fraction averaged 25% (range: 21-30%) in concentration of the total plasma proteins as determined by densitometer. The physical properties of the albumin-like proteins purified from these snakes were compared. These properties, e.g. electrophoretic mobility, isoelectric point, extinction coefficient, and molecular weight, were shown to be strikingly similar to those of human plasma albumin. The physical properties of the plasma albumins of the snakes studied are similar to one another. This similarity does not explain our previous observation that Naja albumin is considerably remote immunologically from those of other elapids (Mao et al., 1983).     

Mid-Tertiary elapid snakes (Squamata, Colubroidea) from Riversleigh, northern Australia: early steps in a continent-wide adaptive radiation

Vertebral and cranial remains of elapid snakes have been collected from fossil assemblages at Riversleigh, north-west Queensland, Australia; most are Miocene but one may be late Oligocene and another as young as Pliocene. The oldest specimen (probably the oldest elapid yet known anywhere) is a vertebra that can be referred provisionally to the extant taxon Laticauda (Hydrophiinae, sensu Slowinski and Keogh, 2000), implying that the basal divergences among Australasian hydrophiine lineages had occurred by the early Miocene, in contrast to most previous estimates for the age of this geographically isolated adaptive radiation. Associated vertebrae and jaw elements from a Late Miocene deposit are described as Incongruelaps iteratus nov. gen. et sp., which has a unique combination of unusual derived characters otherwise found separately in several extant hydrophiine taxa that are only distantly related. Associated vertebrae from other sites, and two parietals from a possibly Pliocene deposit, suggest the presence of several other taxa distinct from extant forms, but the amount of material (and knowledge of variation in extant taxa) is currently insufficient to diagnose these forms. The Tertiary elapids of Riversleigh thus appear to be relatively diverse taxonomically, but low in abundance and, with one exception, not referable to extant taxa below the level of Hydrophiinae. This implies that the present diversity of hydrophiine elapids (31 recognized terrestrial genera, and approximately 16 marine) represents the result of substantial extinction as well as the “cone of increasing diversity” that could be inferred from phylogenetic studies on extant forms.     

Uncoupling ecological innovation and speciation in sea snakes (Elapidae, Hydrophiinae, Hydrophiini)

The viviparous sea snakes (Hydrophiini) are by far the most successful living marine reptiles, with ～ 60 species that comprise a prominent component of shallow-water marine ecosystems throughout the Indo-West Pacific. Phylogenetically nested within the ～ 100 species of terrestrial Australo-Melanesian elapids (Hydrophiinae), molecular timescales suggest that the Hydrophiini are also very young, perhaps only ～ 8-13 Myr old. Here, we use likelihood-based analyses of combined phylogenetic and taxonomic data for Hydrophiinae to show that the initial invasion of marine habitats was not accompanied by elevated diversification rates. Rather, a dramatic three to six-fold increase in diversification rates occurred at least 3-5 Myr after this transition, in a single nested clade: the Hydrophis group accounts for ～ 80% of species richness in Hydrophiini and ～ 35% of species richness in (terrestrial and marine) Hydrophiinae. Furthermore, other co-distributed lineages of viviparous sea snakes (and marine Laticauda, Acrochordus and homalopsid snakes) are not especially species rich. Invasion of the oceans has not (by itself) accelerated diversification in Hydrophiini; novelties characterizing the Hydrophis group alone must have contributed to its evolutionary and ecological success.     

Comparative analysis of prothrombin activators from the venom of Australian elapids

A key component of the venom of many Australian snakes belonging to the elapid family is a toxin that is structurally and functionally similar to that of the mammalian prothrombinase complex. In mammals, this complex is responsible for the cleavage of prothrombin to thrombin and is composed of factor Xa in association with its cofactors calcium, phospholipids, and factor Va. The snake prothrombin activators have been classified on the basis of their requirement for cofactors for activity. The two major subgroups described in Australian elapid snakes, groups C and D, are differentiated by their requirement for mammalian coagulation factor Va. In this study, we describe the cloning, characterization, and comparative analysis of the factor X- and factor V-like components of the prothrombin activators from the venom glands of snakes possessing either group C or D prothrombin activators. The overall domain arrangement in these proteins was highly conserved between all elapids and with the corresponding mammalian clotting factors. The deduced protein sequence for the factor X-like protease precursor, identified in elapids containing either group C or D prothrombin activators, demonstrated a remarkable degree of relatedness to each other (80%-97%). The factor V-like component of the prothrombin activator, present only in snakes containing group C complexes, also showed a very high degree of homology (96%-98%). Expression of both the factor X- and factor V-like proteins determined by immunoblotting provided an additional means of separating these two groups at the molecular level. The molecular phylogenetic analysis described here represents a new approach for distinguishing group C and D snake prothrombin activators and correlates well with previous classifications.     

Isolation and amino acid sequence of a short-chain neurotoxin from an Australian elapid snake, Pseudechis australis.

A short-chain neurotoxin Pseudechis australis a (toxin Pa a) was isolated from the venom of an Australian elapid snake Pseudechis australis (king brown snake) by sequential chromatography on CM-cellulose, Sephadex G-50 and CM-cellulose columns. Toxin Pa a has an LD50 (intravenous) value of 76 micrograms/kg body wt. in mice and consists of 62 amino acid residues. The amino acid sequence of Pa a shows considerable homology with those of short-chain neurotoxins of elapid snakes, especially of true sea snakes.