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.     

Independent innovation in the evolution of paddle-shaped tails in viviparous sea snakes (Elapidae: Hydrophiinae)

The viviparous sea snakes (Hydrophiinae) comprise ~90% of living marine reptiles and display many physical and behavioral adaptations for breathing, diving, and achieving osmotic balance in marine habitats. Among the most important innovations found in marine snakes are their paddle-shaped (dorsoventrally expanded) tails, which provide propulsive thrust in the dense aquatic medium. Here, we reconstruct the evolution of caudal paddles in viviparous sea snakes using a dated molecular phylogeny for all major lineages and computed tomography of internal osteological structures. Bayesian ancestral state reconstructions show that extremely large caudal paddles supported by elongated vertebral processes are unlikely to have been present in the most recent common ancestor of extant sea snakes. Instead, these characters appear to have been acquired independently in two highly marine lineages of relatively recent origin. Both the Aipysurus and Hydrophis lineages have elongated neural spines that support the dorsal edge of their large paddles. However, whereas in the Aipysurus lineage the ventral edge of the paddle is supported by elongated haemapophyses, this support is provided by elongated and ventrally directed pleurapophyses in the Hydrophis lineage. Three semi-marine lineages (Hydrelaps, Ephalophis, and Parahydrophis) form the sister group to the Hydrophis clade and have small paddles with poorly developed dorsal and ventral supports, consistent with their amphibious lifestyle. Overall, our results suggest that not only are the viviparous hydrophiines the only lineage of marine snakes to have acquired extremely large, skeletally supported caudal paddles but also that this innovation has occurred twice in the group in the past ~2-6 million years.     

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.     

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.     

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.     

Marine invasions by non-sea snakes, with thoughts on terrestrial-aquatic-marine transitions

Few species of snakes show extensive adaptations to aquatic environments and even fewer exploit the oceans. A survey of morphology, lifestyles, and habitats of 2552 alethenophidian snakes revealed 362 (14%) that use aquatic environments, are semi-aquatic, or aquatic; about 70 (2.7%) of these are sea snakes (Hydrophiinae and Laticaudinae). The ancient and aquatic family Acrochordidae contains three extant species, all of which have populations inhabiting brackish or marine environments, as well as freshwater. The Homalopsidae have the most ecologically diverse representatives in coastal habitats. Other families containing species exploiting saline waters with populations in freshwater environments include: the Dipsadidae of the western hemisphere, the cosmopolitan Natricidae, the African Grayinae, and probably a few Colubridae. Species with aquatic and semi-aquatic lifestyles are compared with more terrestrial (fossorial, cryptozoic, and arboreal) species for morphological traits and life histories that are convergent with those found in sea snakes; this may provide clues to the evolution of marine snakes and increase our understanding of snake diversity.     

Vertebral evolution and ontogenetic allometry: The developmental basis of extreme body shape divergence in microcephalic sea snakes

Snakes exhibit a diverse array of body shapes despite their characteristically simplified morphology. The most extreme shape changes along the precloacal axis are seen in fully aquatic sea snakes (Hydrophiinae): “microcephalic” sea snakes have tiny heads and dramatically reduced forebody girths that can be less than a third of the hindbody girth. This morphology has evolved repeatedly in sea snakes that specialize in hunting eels in burrows, but its developmental basis has not previously been examined. Here, we infer the developmental mechanisms underlying body shape changes in sea snakes by examining evolutionary patterns of changes in vertebral number and postnatal ontogenetic growth. Our results show that microcephalic species develop their characteristic shape via changes in both the embryonic and postnatal stages. Ontogenetic changes cause the hindbodies of microcephalic species to reach greater sizes relative to their forebodies in adulthood, suggesting heterochronic shifts that may be linked to homeotic effects (axial regionalization). However, microcephalic species also have greater numbers of vertebrae, especially in their forebodies, indicating that somitogenetic effects also contribute to evolutionary changes in body shape. Our findings highlight sea snakes as an excellent system for studying the development of segment number and regional identity in the snake precloacal axial skeleton.     

Sea snake diversity at the Entrecasteaux atolls,Coral Sea,as revealed by video observations at unbaited stations

Coral Reefs - The first survey of sea snakes (Elapidae, Hydrophiinae) at the remote Entrecasteaux atolls, Coral Sea, was conducted using remote unbaited 360° video cameras (RUV360), in 2021....     

Feeding Strategies in Marine Snakes: An Analysis of Evolutionary, Morphological, Behavioral and Ecological Relationships

Analysis of 1,063 stomach contents from 39 species of sea snakesindicates that about one-third of the shallow, warm, marine,Indo-Australian fish families are preyed upon by sea snakes.Families of eels and gobies are taken by the greatest numbersof snake species. Most species of sea snakes feed on fish familieswhose members are relatively sedentary, dwelling along the bottom,within burrows or reef crevices. With one exception, a fishegg-eating specialization found uniquely in the Aipysurus-Emydocephaluslineage, the dietary habits of sea snakes cannot be categorizedaccording to the snakes' three phylogenetic lineages. Eels,mullet-like, rabbitfish-like and goby-like fish forms are takenby all three lineages. Two or three snake species are generalists,and numerous ones specialize on eels, goby-like fish or catfish.There are differences among sea snake species in the relationshipbetween snake neck girth and the maximum diameter of the prey;in the relationships of both snake gape measurements and fanglength, to the type of prey taken; and in the relationship ofsnake shape and body proportions to the prey selected. Severalmodes of feeding have been observed among sea snakes: feedingin nooks and crannies in the bottom or in reefs, cruising nearthe bottom, and feeding in drift lines. Analysis of percentdigestion of stomach contents and projections backward to thetimes of prey capture provides evidence for feeding periodicity.The greatest amount of diet overlap is for two species of seasnakes which do not both occur at the same locality. Where speciesdo co-occur, diet overlap index values are lower. The numbersof species present as well as their relative abundances varyamong localities as does the relative importance of generalists,eel-eaters, egg-eaters and other specialized feeders.     

The invasion biology of the invasive earwig, <Emphasis Type="Italic">Forficula auricularia</Emphasis> in Australasian ecosystems

The European earwig, Forficula auricularia, is a cosmopolitan insect endemic to Europe, West Asia and North Africa, which has invaded many temperate regions of the world including Australia and New Zealand. F. auricularia has been shown to be a complex of morphologically identical, reproductively isolated lineages that possess two distinct clades of mitochondrial DNA. Entomological collection data, historical literature and further field collections were used to develop a greater understanding of Australian F. auricularia’s invasion biology and its current distribution. Genetic analysis of F. auricularia collected from Australia and New Zealand using two mitochondrial genes (COI and a fragment overlapping parts of the COI -COII genes) was also undertaken. To identify the possible source populations of the Australasian invasion these sequences were compared to those from 16 locations within Britain and continental Europe. All Australasian populations were shown to be of the clade B lineage. Tasmanian and New Zealand populations consist of a single subclade comprised of only 4 and 1 haplotypes respectively. The Australian mainland populations also contained a second subclade consisting of up to 11 haplotypes indicating that multiple introductions possibly occurred on the Australian mainland. Comparison of mitochondrial genomes from Australasian and European populations showed the Australian mainland subclade to be most closely related to Portuguese haplotypes, and the Tasmanian and New Zealand clade closely related to those in Brittany, France. No European haplotypes perfectly matched the Australasian sequences. Therefore, the original source populations are still to be identified with harbours on the Iberian Peninsula’s western coast and those on the English Channel likely candidate areas.     

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.     

First records of sea snakes (Elapidae: Hydrophiinae) diving to the mesopelagic zone (>200 m)

Viviparous sea snakes (Elapidae: Hydrophiinae) are fully marine reptiles distributed in the tropical and subtropical waters of the Indian and Pacific Oceans. Their known maximum diving depth ranges between 50 and 100 m and this is thought to limit their ecological ranges to shallow habitats. We report two observations, from industry‐owned remotely operated vehicles, of hydrophiine sea snakes swimming and foraging at depths of approximately 250 m in the Browse Basin on Australia's North West Shelf, in 2014 and 2017. These observations show that sea snakes are capable of diving to the dim‐lit, cold‐water mesopelagic zone, also known as the ‘twilight’ zone. These record‐setting dives raise new questions about the thermal tolerances, diving behaviour and ecological requirements of sea snakes. In addition to significantly extending previous diving records for sea snakes, these observations highlight the importance of university‐industry collaboration in surveying understudied deep‐sea habitats.     

Degeneration patterns of the olfactory receptor genes in sea snakes

The sense of smell relies on the diversity of olfactory receptor (OR) repertoires in vertebrates. It has been hypothesized that different types of ORs are required in terrestrial and marine environments. Here we show that viviparous sea snakes, which do not rely on a terrestrial environment, have significantly lost ORs compared with their terrestrial relatives, supporting the hypothesis. On the other hand, oviparous sea snakes, which rely on a terrestrial environment for laying eggs, still maintain their ORs, reflecting the importance of the terrestrial environment for them. Furthermore, we found one Colubroidea snake (including sea snakes and their terrestrial relatives)‐specific OR subfamily which had diverged widely during snake evolution after the blind snake–Colubroidea snake split. Interestingly, no pseudogenes are included in this subfamily in sea snakes, and this subfamily seems to have been expanding rapidly even in an underwater environment. These findings suggest that the Colubroidea‐specific ORs detect nonvolatile odorants.     

A phylogenetic analysis of variation in reproductive mode within an Australian lizard (Saiphos equalis, Scincidae)

Saiphos equalis , a semi-fossorial scincid lizard from south-eastern Australia, is one of only three reptile species world-wide that are known to display geographic variation in reproductive mode. Uniquely, Saiphos equalis includes populations with three reproductive modes: oviparous with long (15-day) incubation periods; oviparous with short (5-day) incubation periods; and viviparous (0-day incubation periods). No Saiphos populations show 'normal' scincid oviparity (> 30-day incubation period). We used mitochondrial nucleotide sequences ( ND2 and cytochrome b ) to reconstruct relationships among populations from throughout the species' distribution in New South Wales, Australia. Under the phylogenetic species concept, phylogenetic analyses are consistent with the oviparous and viviparous populations of S. equalis being conspecific. Phylogenetic analyses suggest that the long incubation period oviparous lineage is the sister group to all other populations; and that the viviparous populations belong to a cluster of weakly supported clades basal to the short-incubation-period oviparous clade. These clades correspond to variation in reproductive mode and geographic location.     

Molecular phylogeny of elapid snakes and a consideration of their biogeographic history

Evolutionary relationships among the major elapid clades, particularly the taxonomic position of the partially aquatic sea kraits (Latkauda) and the fully aquatic true sea snakes have been the subject of much debate. To discriminate among existing phylogenetic and biogeographic hypotheses, portions of both the 16S rRNA and cytochrome b mitochondrial DNA genes were sequenced from 16 genera and 17 species representing all major elapid snake clades from throughout the world and two non-elapid outgroups. This sequence data yielded 181 informative sites under parsimony. Parsimony analyses of the separate data sets produced trees of broad agreement although less well supported than the single most parsimonious tree resulting from the combined analyses. These results support the following hypotheses: (1) the Afro-Asian cobra radiation forms one or more sister groups to other elapids, (2) American and Asian coral snakes form a clade, corroborating morphological studies, (3) Bungarus forms a sister group to the hydrophiines comprised of Latkauda, terrestrial Australo-Papuan elapids and true sea snakes, (4) Latkauda and true sea snakes do not form a monophyletic group but instead each group shares an independent history with terrestrial Australo-Papuan elapids, corroborating previous studies, (5) a lineage of Melanesian elapids forms the sister group to Latkauda, terrestrial Australian species and true sea snakes. In agreement with previous morphologically based studies, the sequence data suggests that Bungarus and Latkauda represent transitional clades between the elapine 'palatine erectors' and hydrophiine 'palatine draggers'. Both intra and inter-clade genetic distances are considerable, implying that each of the major radiations have had long independent histories. I suggest an African, Asian, or Afro-Asian origin for elapids as a group, with independent Asian origins for American coral snakes and the hydrophiines.     

Systematic revision,molecular phylogeny and biogeography of the antlion tribe Acanthoplectrini (Neuroptera: Myrmeleontidae: Dendroleontinae), with emphasis on the Oriental lineage

The tribe Acanthoplectrini (Myrmeleontidae: Dendroleontinae) includes a group of antlion genera widely distributed across the Australasian and Oriental regions. The intergeneric and interspecific relationships between or within the Australian and Oriental lineages of this tribe as well as their historical biogeography remain largely unexplored. Here, we present a molecular phylogenetic and biogeographic analyses of Acanthoplectrini to infer the diversification history of this tribe, with emphasis on the Oriental lineage. Both the Oriental and Australian lineages are monophyletic and recovered as sister groups. Ancestral area reconstruction suggests that the ancestor of Acanthoplectrini might have been once widely distributed from Indochina to Australia and then split into the Oriental and Australian lineages during the early-Miocene. Our analyses recovered northeastern Indochina and south China as the ancestral range of the Oriental Acanthoplectrini. During the mid-Miocene to the mid-Pliocene, orographic events such as the rising of mountain ranges (including the Himalayas) and the formation of major islands in southeastern Asia triggered several dispersal and vicariance events in the Oriental Acanthoplectrini, driving their speciation. We revise the classification of the Oriental Acanthoplectrini, establishing the new genus Paralayahima gen. n. , which is recovered sister to Layahima Navás. Moreover, we describe four new species of Layahima, Layahima aspoeckorum sp. n. , Layahima monba sp. n. , Layahima lhoba sp. n. and Layahima xinliae sp. n. , and we reinstate two previously synonymized species, Layahima melanocoris (Yang) stat. rev. and comb. n. and Layahima nebulosa Navás stat. rev.     

Reproductive mode may determine geographic distributions in Australian venomous snakes (Pseudechis,Elapidae)

Summary Why are viviparous squamate reptiles more common in cold climates, and oviparous ones in warmer areas? The usual explanation is that (1) oviparous squamates cannot reproduce successfully in cold areas because soil temperatures are too low for embryonic development; and (2) viviparous squamates experience lower survivorship or reproductive success than oviparous taxa in warmer areas. These hypotheses suggest that the boundaries of geographic distributions of congeneric oviparous and viviparous squamates should be predictable from data on thermal tolerances of embryos, and estimated temperatures of soils and gravid female reptiles throughout the potential geographic range of the taxon. In large venomous Australian snakes of the genus Pseudechis, distributional boundaries of oviparous and viviparous taxa can be accurately predicted from such data. This predictive ability, if substantiated by studies of other reproductively biomodal squamate taxa, would support the putative role of reproductive mode as a direct determinant of reptilian geographic distributions.     

The amino acid sequence and position of the free thiol group of a short-chain neurotoxin from common-death-adder (Acanthophis antarcticus) venom.

The amino acid sequence of a short-chain neurotoxin Acanthophis antarcticus c (toxin Aa c) from the venom of an Australian elapid snake, the common death adder (Acanthophis antarcticus, subfamily Acanthophiinae) was elucidated. Toxin Aa c is composed of 62 amino acid residues, including eight half-cystine residues and a cysteine residue. The amino acid sequence of toxin Aa c is homologous with those of other short-chain neurotoxins found in snakes of the family Elapidae, especially with those from snakes of the subfamily Hydrophiinae. The single cysteine residue was located in position 4. Toxin Aa c has a lethal dose (LD50) of 0.08 micrograms/g body weight of mouse on intramuscular injection.     

Digging their own macroevolutionary grave: fossoriality as an evolutionary dead end in snakes

The tree of life is highly asymmetrical in its clade wise species richness, and this has often been attributed to variation in diversification rates either across time or lineages. Variations across lineages are usually associated with traits that increase lineage diversification. Certain traits can also hinder diversification by increasing extinction, and such traits are called evolutionary dead ends. Ecological specialization has usually been considered as an evolutionary dead end. However, recent analyses of specializations along single axes have provided mixed support for this model. Here, we test if fossoriality, a trait that forces specialization at multiple axes, acts as an evolutionary dead end in squamates (lizards and snakes) using recently developed phylogenetic comparative methods. We show that fossoriality is an evolutionary dead end in snakes but not in lizards. Fossorial snakes exhibit reduced speciation and increased extinction compared to nonfossorial snakes. Our analysis also indicates that transition rates from fossoriality to nonfossoriality in snakes are significantly lower than transition rates from nonfossoriality to fossoriality. Overall our results suggest that broad‐scale ecological interactions that lead to specialization at multiple axes limit diversification.