Multilocus phylogeny and recent rapid radiation of the viviparous sea snakes (Elapidae: Hydrophiinae)

New chromosome-scale genomes provide insights into marine adaptations of sea snakes (Hydrophis: Elapidae)

BACKGROUND

Sea snakes underwent a complete transition from land to sea within the last ~ 15 million years, yet they remain a conspicuous gap in molecular studies of marine adaptation in vertebrates.

RESULTS

Here, we generate four new annotated sea snake genomes, three of these at chromosome-scale (Hydrophis major, H. ornatus and H. curtus), and perform detailed comparative genomic analyses of sea snakes and their closest terrestrial relatives. Phylogenomic analyses highlight the possibility of near-simultaneous speciation at the root of Hydrophis, and synteny maps show intra-chromosomal variations that will be important targets for future adaptation and speciation genomic studies of this system. We then used a strict screen for positive selection in sea snakes (against a background of seven terrestrial snake genomes) to identify genes over-represented in hypoxia adaptation, sensory perception, immune response and morphological development.

CONCLUSIONS

We provide the best reference genomes currently available for the prolific and medically important elapid snake radiation. Our analyses highlight the phylogenetic complexity and conserved genome structure within Hydrophis. Positively selected marine-associated genes provide promising candidates for future, functional studies linking genetic signatures to the marine phenotypes of sea snakes and other vertebrates.

Polyphyly of the Dinurinae Looss, 1907 (Digenea: Hemiuridae) and resurrection of the Mecoderinae Skrjabin & Guschanskaja, 1954 based on novel collection of Tubulovesicula laticaudi Parukhin, 1969 from marine elapid snakes in Sri Lanka

With one exception, the only known hemiurid trematodes that do not use teleost fishes as definitive hosts instead occur in marine elapid snakes. These comprise six species across four genera and three subfamilies, and so presumably indicate at least three independent invasions of marine snakes from teleost fishes. Here, one of these taxa, Tubulovesicula laticaudi Parukhin, 1969 (= T. orientalis Chattopadhyaya, 1970 n. syn.) is reported from Sri Lanka, collected from Shaw's sea snake Hydrophis curtus (Shaw) (Elapidae: Hydrophiinae: Hydrophinii), the annulated sea snake H. cyanocinctus Daudin and the yellow sea snake H. spiralis (Shaw) off Nayaru in the Bay of Bengal, and from H. spiralis in Portugal Bay, Gulf of Mannar. Novel molecular data, for COI mtDNA and ITS2 and 28S rDNA, are the first for a species of Tubulovesicula Yamaguti, 1934. Nominally, Tubulovesicula belongs in the Dinurinae Looss, 1907, but in phylogenetic analyses based on 28S rDNA, our sequences for T. laticaudi resolved relatively distant from that for representatives of Dinurus Looss, 1907, the type-genus, rendering the subfamily polyphyletic. Tubulovesicula laticaudi resolved closest to data for the type-species of the Plerurinae Gibson & Bray, 1979, but that subfamily is also polyphyletic. These findings lead us to re-evaluate an alternative classification considered by Gibson & Bray (1979). We propose restricting the Dinurinae for forms with a permanent sinus-organ (Dinurus, Ectenurus Looss, 1907; Erilepturus Woolcock, 1935; Paradinurus Vigueras, 1958; Qadriana Bilqees, 1971) and resurrect the Mecoderinae Skrjabin & Guschanskaja, 1954 for forms with a temporary sinus-organ (Mecoderus Manter, 1940, Allostomachicola Yamaguti, 1958, Stomachicola Yamaguti, 1934 and Tubulovesicula).

Sexual dimorphism in aipysurine sea snakes (Elapidae, Hydrophiinae)

The transition from terrestrial to aquatic life by hydrophiine elapid snakes modified targets of natural selection and likely affected sexual selection also. Thus, the shift to marine life also might have affected sexual dimorphism. Our measurements of 419 preserved specimens of six species of aipysurine snakes (genera Emydocephalus and Aipysurus) revealed sexual dimorphism in mean adult snout-vent length (SVL), body width relative to SVL, lengths and widths of heads and tails relative to SVL, and eye diameter relative to head length. Females averaged larger than males in all taxa, and generally were wider-bodied with shorter and wider tails and smaller eyes. For other traits, sexual dimorphism varied among species: for example, relative head length ranged from male-biased to female-biased, and head shape (width relative to length) was highly dimorphic only in A. laevis. The transition to marine life may have eliminated male-male combat (reducing selection for large males) and favoured visual rather than pheromone-based mate-searching (favouring larger eyes in males). Variation in head-size dimorphism may reflect intersexual niche partitioning, with different taxa following different trajectories. Repeated evolutionary transitions from terrestrial to aquatic life in snakes provide a powerful opportunity to explore selective forces on sexually dimorphic traits.

Singapore's herpetofauna: updated and annotated checklist, history, conservation, and distribution

Given Singapore's location at the confluence of important maritime trading routes, and that it was established as a British East India Company trading post in 1819, it is unsurprising that Singapore has become one of the centres of natural history collecting and research in Southeast Asia. Despite its small size, Singapore is home to a diverse herpetofauna assemblage and boasts a rich herpetological history. The first systematic studies of Singapore's herpetofauna (within the Linnaean binomial framework) date back to Stamford Raffles and the naturalists hired by him who first came to the island in 1819. Specimens that were collected during and after this time were deposited in museums worldwide. Over time, 39 species from Singapore were described as new to science. Due to the entrepôt nature of Singapore with its associated purchasing and trading of specimens (both alive and dead), poor record-keeping, and human introductions, numerous extraneous species from outside of Singapore were reported to occur on the island. Such issues have left a complicated legacy of ambiguous records and taxonomic complications concerning the identity of Singapore's species-rich herpetofauna, many of which were only resolved in the past 30-40 years. By compiling a comprehensive collection of records and publications relating to the herpetofauna of Singapore, we construct an updated and more accurate listing of the herpetofauna of Singapore. Our investigation culminated in the evaluation of 309 species, in which we compiled a final species checklist recognising 166 species (149 native and 17 non-native established species). Among the 149 native species are two caecilians, 24 frogs, one crocodilian, 13 turtles (three visitors), 34 lizards, and 75 snakes. Of the 17 non-native species are five frogs, four turtles, six lizards, and two snakes. The remaining 143 species represent species to be excluded from Singapore's herpetofauna species checklist. For each of the 309 species examined, we provide species accounts and explanatory annotations. Furthermore, we discuss Singapore's herpetofauna from a historical and conservation perspective. Immediate deforestation and nationwide urbanisation following colonisation completely eliminated many species from throughout much of the country and restricted them to small, degraded forest patches. We hope this publication highlights the importance of publishing observations and serves as a valuable resource to future researchers, naturalists, biological consultants, and policy makers in initiating studies on species ecology, distribution, status, and promoting conservation efforts to safeguard Singapore's herpetofauna.

Medically important snakes and snakebite envenoming in Iran

Snakebite is a common health condition in Iran with a diverse snake fauna, especially in tropical southern and mountainous western areas of the country with plethora of snake species. The list of medically important snakes, circumstances and effects of their bite, and necessary medical care require critical appraisal and should be updated regularly. This study aims to review and map the distributions of medically important snake species of Iran, re-evaluate their taxonomy, review their venomics, describe the clinical effects of envenoming, and discuss medical management and treatment, including the use of antivenom. Nearly 350 published articles and 26 textbooks with information on venomous and mildly venomous snake species and snakebites of Iran, were reviewed, many in Persian (Farsi) language, making them relatively inaccessible to an international readership. This has resulted in a revised updated list of Iran's medically important snake species, with taxonomic revisions of some, compilation of their morphological features, remapping of their geographical distributions, and description of species-specific clinical effects of envenoming. Moreover, the antivenom manufactured in Iran is discussed, together with treatment protocols that have been developed for the hospital management of envenomed patients.

Genetic diversity among sea snakes of the genus Hydrophis (Elapidae, Reptilia) in the Persian Gulf and Gulf of Oman

Sea snakes of the genus Hydrophis are important components of animal diversity in Iranian waters of the Persian Gulf and Gulf of Oman. Ten species of Hydrophis have been identified from the these waters and, in this study, genetic structure of seven species was compared with other populations in the eastern Indian Ocean and the West Pacific. We found that six species (H.platurus, H.cyanocinctus, H.spiralis, H.schistosus, H.gracilis, and H.lapemiodes) show high genetic similarity with conspecific populations in the Indian Ocean and Australia. However, H.curtus from southern Iran shows a high level of genetic differentiation from conspecific populations in Sri Lanka and Indonesia (0.6% and 6% genetic distance from Sri Lankan samples for 16S and COI gene fragments, respectively). Variation between Iranian and Southeast Asian populations may reflect new genetic lineages and suggest the need of further morphological evaluations to re-evaluate their taxonomic position.

Macroevolution in axial morphospace: innovations accompanying the transition to marine environments in elapid snakes

Sea snakes in the Hydrophis-Microcephalophis clade (Elapidae) show exceptional body shape variation along a continuum from similar forebody and hindbody girths, to dramatically reduced girths of the forebody relative to hindbody. The latter is associated with specializations on burrowing prey. This variation underpins high sympatric diversity and species richness and is not shared by other marine (or terrestrial) snakes. Here, we examined a hypothesis that macroevolutionary changes in axial development contribute to the propensity, at clade level, for body shape change. We quantified variation in the number and size of vertebrae in two body regions (pre- and post-apex of the heart) for approximately 94 terrestrial and marine elapids. We found Hydrophis-Microcephalophis exhibit increased rates of vertebral evolution in the pre- versus post-apex regions compared to all other Australasian elapids. Unlike other marine and terrestrial elapids, axial elongation in Hydrophis-Microcephalophis occurs via the preferential addition of vertebrae pre-heart apex, which is the region that undergoes concomitant shifts in vertebral number and size during transitions along the relative fore- to hindbody girth axis. We suggest that this macroevolutionary developmental change has potentially acted as a key innovation in Hydrophis-Microcephalophis by facilitating novel (especially burrowing) prey specializations that are not shared with other marine snakes.

Multilocus phylogeny, natural history traits and classification of natricine snakes (Serpentes: Natricinae)

Natricine snakes are geographically widespread, species rich (with ~250 extant species) and both morphologically and ecologically diverse. We present a multilocus DNA sequence phylogeny for 249 natricine specimens representing 189 named species, including 69 specimens and 21 species not previously sampled. Our inferred Bayesian and maximum likelihood trees form the basis for evaluations of genus-level classification, historical biogeography, lineage diversification, and dietary, habit and reproductive-mode diversity and evolution, although several, mostly deeper, relationships remain poorly resolved. The optimal trees support natricine origins in Asia, with dispersals to Australo-Melanesia, sub-Saharan Africa (including Seychelles Archipelago, excluding Aldabra), Europe and North Africa and into North and Central America. Viviparity appears to have evolved independently three times in Natricinae but was not significantly associated with an aquatic habit. We found limited associations between habit and diet categories. We propose generic reallocations for four natricine species and highlight other points of uncertainty in natricine classification.

Morphology, Diet, and Reproduction of Coastal Hydrophis Sea Snakes (Elapidae: Hydrophiinae) at Their Northern Distribution Limit

The Ryukyu Archipelago represents the northern distribution limit for hydrophiine sea snakes, the largest group of marine reptiles. Ryukyuan sea snakes may have developed distinct local adaptations in morphology and ecology, but they have been poorly studied. We examined preserved specimens of 111 Hydrophismelanocephalusand 61 Hydrophis ornatusfrom the Ryukyu Archipelago to obtain data on morphology, diet, and reproduction. Sexual size dimorphism was detected in H. melanocephalus (mean ± standard deviation of adult snout-vent length: SVL, females 1062 ± 141 mm vs. males 959 ± 96 mm) but not in H. ornatus. Female H. melanocephalus had larger head widths and shorter tail lengths relative to SVL compared to males. Relative girth was low in neonates of both species (1.0-1.3), but increased in adults to about 1.7-2.6 in H. melanocephalus and 1.3-1.8 in female H. ornatus. Stomach contents of H. melanocephalus consisted of ophichthid and congrid eels, a sand diver, and gobies, whereas in H. ornatus, gobies and a goat fish were found. Litter size of three reproductive H. melanocephalus ranged from five to seven, and parturition seems to occur from August to October. Litter size of six H. ornatus ranged from two to seven, and was correlated with maternal SVL. Parturition in H. ornatus probably occurs around November. Different selective forces related to locomotion, feeding and predation risk, which influence the pregnant mother and neonates, may have resulted in having few, long but slender offspring that show positive allometric growth in hind-body girth.

Two reference-quality sea snake genomes reveal their divergent evolution of adaptive traits and venom systems

True sea snakes (Hydrophiini) are among the last and most successful clades of vertebrates that show secondary marine adaptation, exhibiting diverse phenotypic traits and lethal venom systems. To better understand their evolution, we generated the first chromosome-level genomes of two representative Hydrophiini snakes, Hydrophis cyanocinctus and H. curtus. Through comparative genomics we identified a great expansion of the underwater olfaction-related V2R gene family, consisting of more than 1,000 copies in both snakes. A series of chromosome rearrangements and genomic structural variations were recognized, including large inversions longer than 30 megabase (Mb) on sex chromosomes which potentially affect key functional genes associated with differentiated phenotypes between the two species. By integrating multiomics we found a significant loss of the major weapon for elapid predation, three-finger toxin genes, which displayed a dosage effect in H. curtus. These genetic changes may imply mechanisms that drove the divergent evolution of adaptive traits including prey preferences between the two closely related snakes. Our reference-quality sea snake genomes also enrich the repositories for addressing important issues on the evolution of marine tetrapods, and provide a resource for discovering marine-derived biological products.

Sexual dimorphism and skull size and shape in the highly specialized snake species, Aipysurus eydouxii (Elapidae: Hydrophiinae)

Background

Snakes exhibit sexual dimorphism in both head size and shape. Such differences are often attributed to different reproductive roles and feeding habits. We aim to investigate how sexual dimorphism is displayed in the highly specialised fish-egg-eating snake, Aipysurus eydouxii, by analysing two complementary features: body size and skull morphology.

Methods

We used data on body length, weight, and skull shape from 27 measurements of 116 males and females of A. eydouxii. We investigated both sexual dimorphism and allometric (multivariate and bi-variate) properties of skull growth in the analysed data set.

Results

We found that although there was female-biased sexual size dimorphism in body length, females were not heavier than males, contrary to what is commonly observed pattern among snakes. Moreover, females tend to possess relatively smaller heads than males. However, we only found very subtle differences in skull shape reflected in nasal width, mandibular fossa, quadrate crest and quadrate length.

Discussion

We suggest that the feeding specialisation in A. eydouxii does not allow for an increase in body thickness and the size of the head above a certain threshold. Our results may be interpreted as support for prey-size divergence as a factor driving skull dimorphism since such species in which the sexes do not differ in prey size also shows very subtle or no differences in skull morphology.

De Novo Venom-Gland Transcriptomics of Spine-Bellied Sea Snake (Hydrophis curtus) from Penang, Malaysia-Next-Generation Sequencing, Functional Annotation and Toxinological Correlation

Envenomation resulted from sea snake bite is a highly lethal health hazard in Southeast Asia. Although commonly caused by sea snakes of Hydrophiinae, each species is evolutionarily distinct and thus, unveiling the toxin gene diversity within individual species is important. Applying next-generation sequencing, this study investigated the venom-gland transcriptome of Hydrophis curtus (spine-bellied sea snake) from Penang, West Malaysia. The transcriptome was de novo assembled, followed by gene annotation and sequence analyses. Transcripts with toxin annotation were only 96 in number but highly expressed, constituting 48.18% of total FPKM in the overall transcriptome. Of the 21 toxin families, three-finger toxins (3FTX) were the most abundantly expressed and functionally diverse, followed by phospholipases A₂. Lh_FTX001 (short neurotoxin) and Lh_FTX013 (long neurotoxin) were the most dominant 3FTXs expressed, consistent with the pathophysiology of envenomation. Lh_FTX001 and Lh_FTX013 were variable in amino acid compositions and predicted epitopes, while Lh_FTX001 showed high sequence similarity with the short neurotoxin from Hydrophis schistosus, supporting cross-neutralization effect of Sea Snake Antivenom. Other toxins of low gene expression, for example, snake venom metalloproteinases and L-amino acid oxidases not commonly studied in sea snake venom were also identified, enriching the knowledgebase of sea snake toxins for future study.

Demographic analyses of marine and terrestrial snakes (Elapidae) using whole genome sequences

The question of whether spatial aspects of evolution differ in marine versus terrestrial realms has endured since Ernst Mayr's 1954 essay on marine speciation. Marine systems are often suggested to support larger and more highly connected populations, but quantitative comparisons with terrestrial systems have been lacking. Here, we compared the population histories of marine and terrestrial elapid snakes using the pairwise sequentially Markovian coalescent (PSMC) model to track historical fluctuations in species' effective population sizes (N_e ) from individual whole-genome sequences. To do this we generated a draft genome for the olive sea snake (Aiysurus laevis) and analysed this alongside six published elapid genomes and their sequence reads (marine species Hydrophis curtus, H. melanocephalus and Laticauda laticaudata; terrestrial species Pseudonaja textilis, Naja Naja and Notechis scutatus). Counter to the expectation that marine species should show higher overall N_e and less pronounced fluctuations in N_e , our analyses reveal demographic patterns that are highly variable among species and do not clearly correspond to major ecological divisions. At deeper time intervals, the four marine elapids appear to have experienced relatively stable N_e , while each terrestrial species shows a prominent upturn in N_e starting at ~4 million years ago (Ma) followed by an equally strong decline. However, over the last million years, all seven species show strong and divergent fluctuations. Estimates of N_e in the most recent intervals (~10 kya) are lowest in two of four marine species (H. melanocephalus and Laticauda), and do not correspond to contemporary range sizes in marine or terrestrial taxa.

Visual adaptation of opsin genes to the aquatic environment in sea snakes

BACKGROUND

Evolutionary transitions from terrestrial to aquatic life history cause drastic changes in sensory systems. Indeed, the drastic changes in vision have been reported in many aquatic amniotes, convergently. Recently, the opsin genes of the full-aquatic sea snakes have been reported. However, those of the amphibious sea snakes have not been examined in detail.

RESULTS

Here, we investigated opsin genes and visual pigments of sea snakes. We determined the sequences of SWS1, LWS, and RH1 genes from one terrestrial, three amphibious and four fully-aquatic elapids. Amino acid replacements at four and one spectra-tuning positions were found in LWS and RH1, respectively. We measured or predicted absorption of LWS and RH1 pigments with A1-derived retinal. During their evolution, blue shifts of LWS pigments have occurred stepwise in amphibious sea snakes and convergently in both amphibious and fully-aquatic species.

CONCLUSIONS

Blue shifted LWS pigments may have adapted to deep water or open water environments dominated by blue light. The evolution of opsins differs between marine mammals (cetaceans and pinnipeds) and sea snakes in two fundamental ways: (1) pseudogenization of opsins in marine mammals; and (2) large blue shifts of LWS pigments in sea snakes. It may be possible to explain these two differences at the level of photoreceptor cell composition given that cone and rod cells both exist in mammals whereas only cone cells exist in fully-aquatic sea snakes. We hypothesize that the differences in photoreceptor cell compositions may have differentially affected the evolution of opsins in divergent amniote lineages.

New Environment, New Invaders-Repeated Horizontal Transfer of LINEs to Sea Snakes

Although numerous studies have found horizontal transposon transfer (HTT) to be widespread across metazoans, few have focused on HTT in marine ecosystems. To investigate potential recent HTTs into marine species, we searched for novel repetitive elements in sea snakes, a group of elapids which transitioned to a marine habitat at most 18 Ma. Our analysis uncovered repeated HTTs into sea snakes following their marine transition. The seven subfamilies of horizontally transferred LINE retrotransposons we identified in the olive sea snake (Aipysurus laevis) are transcribed, and hence are likely still active and expanding across the genome. A search of 600 metazoan genomes found all seven were absent from other amniotes, including terrestrial elapids, with the most similar LINEs present in fish and marine invertebrates. The one exception was a similar LINE found in sea kraits, a lineage of amphibious elapids which independently transitioned to a marine environment 25 Ma. Our finding of repeated horizontal transfer events into marine snakes greatly expands past findings that the marine environment promotes the transfer of transposons. Transposons are drivers of evolution as sources of genomic sequence and hence genomic novelty. We identified 13 candidate genes for HTT-induced adaptive change based on internal or neighboring HTT LINE insertions. One of these, ADCY4, is of particular interest as a part of the KEGG adaptation pathway “Circadian Entrainment.” This provides evidence of the ecological interactions between species influencing evolution of metazoans not only through specific selection pressures, but also by contributing novel genomic material.

Exploring the functional meaning of head shape disparity in aquatic snakes

Phenotypic diversity, or disparity, can be explained by simple genetic drift or, if functional constraints are strong, by selection for ecologically relevant phenotypes. We here studied phenotypic disparity in head shape in aquatic snakes. We investigated whether conflicting selective pressures related to different functions have driven shape diversity and explore whether similar phenotypes may give rise to the same functional output (i.e., many-to-one mapping of form to function). We focused on the head shape of aquatically foraging snakes as they fulfill several fitness-relevant functions and show a large amount of morphological variability. We used 3D surface scanning and 3D geometric morphometrics to compare the head shape of 62 species in a phylogenetic context. We first tested whether diet specialization and size are drivers of head shape diversification. Next, we tested for many-to-one mapping by comparing the hydrodynamic efficiency of head shape characteristic of the main axes of variation in the dataset. We 3D printed these shapes and measured the forces at play during a frontal strike. Our results show that diet and size explain only a small amount of shape variation. Shapes did not fully functionally converge as more specialized aquatic species evolved a more efficient head shape than others. The shape disparity observed could thus reflect a process of niche specialization.

The life aquatic: an association between habitat type and skin thickness in snakes

An aquatic animal faces challenges not encountered by its terrestrial counterparts, promoting adaptive responses in multiple traits. For example, a thicker dermis might protect snakes when they are pushed against sharp objects by water currents, and might enable a snake to shed fouling organisms attached to its skin. We thus predicted that marine snakes should have thicker skin than terrestrial species, and that smaller sea snakes should have relatively thicker skin (because absolute, not relative, thickness determines vulnerability to fouling). Measurements of 192 snakes of 44 species supported those predictions. Many (but not all) sea snakes have skins 50% thicker than those of terrestrial and amphibious snake species, representing multiple independent evolutionary origins of thicker skin (in acrochordids, Laticauda sea kraits and both main clades of hydrophiine sea snakes). Marine snakes showed different allometries of skin thickness compared with their terrestrial counterparts; larger snakes had thicker skin within and among species of amphibious and terrestrial snakes, but larger aquatic snake species had thinner skin compared with smaller taxa. Interspecific variation in skin thickness was primarily due to increased collagen in the deep dermis, a physical barrier well suited to protecting against physical injury and to resisting penetration by epibionts.

Novel vascular plexus in the head of a sea snake (Elapidae, Hydrophiinae) revealed by high-resolution computed tomography and histology

Novel phenotypes are often linked to major ecological transitions during evolution. Here, we describe for the first time an unusual network of large blood vessels in the head of the sea snake Hydrophis cyanocinctus. MicroCT imaging and histology reveal an intricate modified cephalic vascular network (MCVN) that underlies a broad area of skin between the snout and the roof of the head. It is mostly composed of large veins and sinuses and converges posterodorsally into a large vein (sometimes paired) that penetrates the skull through the parietal bone. Endocranially, this blood vessel leads into the dorsal cerebral sinus, and from there, a pair of large veins depart ventrally to enter the brain. We compare the condition observed in H. cyanocinctus with that of other elapids and discuss the possible functions of this unusual vascular network. Sea snakes have low oxygen partial pressure in their arterial blood that facilitates cutaneous respiration, potentially limiting the availability of oxygen to the brain. We conclude that this novel vascular structure draining directly to the brain is a further elaboration of the sea snakes' cutaneous respiratory anatomy, the most likely function of which is to provide the brain with an additional supply of oxygen.

Heterochronic Shifts Mediate Ecomorphological Convergence in Skull Shape of Microcephalic Sea Snakes

Morphological variation among the viviparous sea snakes (Hydrophiinae), a clade of fully aquatic elapid snakes, includes an extreme “microcephalic” ecomorph that has a very small head atop a narrow forebody, while the hind body is much thicker (up to three times the forebody girth). Previous research has demonstrated that this morphology has evolved at least nine times as a consequence of dietary specialization on burrowing eels, and has also examined morphological changes to the vertebral column underlying this body shape. The question addressed in this study is what happens to the skull during this extreme evolutionary change? Here we use X-ray micro-computed tomography and geometric morphometric methods to characterize cranial shape variation in 30 species of sea snakes. We investigate ontogenetic and evolutionary patterns of cranial shape diversity to understand whether cranial shape is predicted by dietary specialization, and examine whether cranial shape of microcephalic species may be a result of heterochronic processes. We show that the diminutive cranial size of microcephalic species has a convergent shape that is correlated with trophic specialization to burrowing prey. Furthermore, their cranial shape is predictable for their size and very similar to that of juvenile individuals of closely related but non-microcephalic sea snakes. Our findings suggest that heterochronic changes (resulting in pedomorphosis) have driven cranial shape convergence in response to dietary specializations in sea snakes.

Structure and function of skin in the pelagic sea snake, Hydrophis platurus

We describe and interpret the functional morphology of skin of the Yellow-bellied sea snake, Hydrophis platurus. This is the only pelagic sea snake, and its integument differs from what is known for other species of snakes. In gross appearance, the scales of H. platurus consist of non-overlapping, polygonal knobs with flattened outer surfaces bearing presumptive filamentous sensillae. The deep recesses between scales ('hinge') entrap and wick water over the body surface, with mean retention of 5.1 g/cm of skin surface, similar to that determined previously for the roughened, spiny skin of marine file snakes, Acrochordus granulatus. This feature possibly serves to maintain the skin wet when the dorsal body protrudes above water while floating on calm oceanic slicks where they forage. In contrast with other snakes, including three species of amphibious, semi-marine sea kraits (Laticauda spp.), the outer corneous β-protein layer consists of a syncytium that is thinner than seen in most other species. The subjacent α-layer is also thin, and lipid droplets and lamellar bodies are seen among the immature, cornifying α-cells. A characteristic mesos layer, comprising the water permeability barrier, is either absent or very thin. These features are possibly related to (1) permeability requirements for cutaneous gas exchange, (2) reduced gradient for water efflux compared with terrestrial environments, (3) less need for physical protection in water compared with terrestrial ground environments, and (4) increased frequency of ecdysis thought to be an anti-fouling mechanism. The lipogenic features of the α-layer possibly compensate for the reduced or absent mesos layer, or produce layers of cells that comprise what functionally might be termed a mesos layer, but where the organization of barrier lipids nonetheless appears less robust than what is characteristically seen in squamates.

Venom Proteome of Spine-Bellied Sea Snake (Hydrophis curtus) from Penang, Malaysia: Toxicity Correlation, Immunoprofiling and Cross-Neutralization by Sea Snake Antivenom

The venom proteome of Hydrophis curtus (synonym: Lapemis hardwickii) from Penang, Malaysia was investigated with nano-electrospray ionization-liquid chromatography tandem mass spectrometry (ESI-LCMS/MS) of the reverse-phase high-performance liquid chromatography (HPLC) venom fractions. Thirty distinct protein forms were identified as toxins from ten families. The three major protein families were phospholipase A₂ (PLA₂, 62.0% of total venom proteins), three-finger toxin (3FTX, 26.33%) and cysteine-rich secretory protein (CRiSP, 9.00%). PLA₂ comprises diverse homologues (11 forms), predominantly the acidic subtypes (48.26%). 3FTX composed of one short alpha-neurotoxin (SNTX, 22.89%) and four long alpha-neurotoxins (LNTX, 3.44%). Both SNTX and LNTX were lethal in mice (intravenous LD₅₀ = 0.10 and 0.24 μg/g, respectively) but the PLA₂ were non-lethal (LD₅₀ >1 μg/g). The more abundant and toxic SNTX appeared to be the main driver of venom lethality (holovenom LD₅₀ = 0.20 μg/g). The heterologous Sea Snake Antivenom (SSAV, Australia) effectively cross-neutralized the venom (normalized potency = 9.35 mg venom neutralized per g antivenom) and the two neurotoxins in vivo, with the LNTX being neutralized more effectively (normalized potency = 3.5 mg toxin/g antivenom) than SNTX (normalized potency = 1.57 mg/g). SSAV immunorecognition was strong toward PLA₂ but moderate-to-weak toward the alpha-neurotoxins, indicating that neutralization of the alpha-neurotoxins should be further improved.

Rates of population differentiation and speciation are decoupled in sea snakes

Comparative phylogeography can inform many macroevolutionary questions, such as whether species diversification is limited by rates of geographical population differentiation. We examined the link between population genetic structure and species diversification in the fully aquatic sea snakes (Hydrophiinae) by comparing mitochondrial phylogeography across northern Australia in 16 species from two closely related clades that show contrasting diversification dynamics. Contrary to expectations from theory and several empirical studies, our results show that, at the geographical scale studied here, rates of population differentiation and speciation are not positively linked in sea snakes. The eight species sampled from the rapidly speciating Hydrophis clade have weak population differentiation that lacks geographical structure. By contrast, all eight sampled Aipysurus-Emydocephalus species show clear geographical patterns and many deep intraspecific splits, but have threefold slower speciation rates. Alternative factors, such as ecological specialization, species duration and geographical range size, may underlie rapid speciation in sea snakes.

Spine-bellied sea snake (Hydrophis curtus) venom shows greater skeletal myotoxicity compared with cardiac myotoxicity

For the first time the impedance-based xCELLigence real-time cell analysis system was used to measure the myotoxicity of sea snake venom. With a focus on the spine-bellied sea snake (Hydrophis curtus), the venom of four sea snake species and three terrestrial snake species were compared for myotoxicity against a human skeletal muscle cell line (HSkMC). Hydrophis curtus venom was also tested on a human cardiac muscle cell line (HCM). Surprisingly, all four sea snake venoms tested on HSkMC produced an initial 100-280% rise in xCELLigence cell index that peaked within the first two hours before falling. The cell index rise of H. curtus venom was correlated with the WST-1 cell proliferation assay, which demonstrated an increase in mitochondrial metabolism. The myotoxicity of H. curtus was 4.7-8.2 fold less potent than the other sea snakes tested, the Australian beaked sea snake (Hydrophis zweifeli), the elegant sea snake (Hydrophis elegans) and the olive sea snake (Aipysurus laevis). If our cell-based results translate to H. curtus envenomations, this implies that H. curtus would be less myotoxic than the other three. Yet the myotoxicity of H. curtus venom to cardiac muscle cells was nine times weaker than for skeletal muscle cells, providing evidence that the venom has a selective effect on skeletal muscle cells. This evidence, combined with the slow-acting nature of the venom, supports a digestive role for sea snake myotoxins.

The Venom of the Spine-Bellied Sea Snake (Hydrophis curtus): Proteome, Toxin Diversity and Intraspecific Variation

The spine-bellied sea snake (Hydrophis curtus) is known to cause human deaths, yet its venom composition has not yet been proteomically characterised. An in-depth proteomic analysis was performed on H. curtus venom from two different seasons, January and June, corresponding to adults and subadults, respectively. Venoms from adult and subadult H. curtus individuals were compared using reversed-phase high-performance liquid chromatography (RP-HPLC), matrix-assisted laser desorption ionisation-time of flight (MALDI-TOF) mass spectrometry and liquid chromatography electrospray ionisation mass spectrometry (LC-ESI-MS) to detect intraspecific variation, and the molecular weight data obtained with ESI-MS were used to assess toxin diversity. RP-HPLC and LC-ESI-MS/MS were used to characterise the venom proteome and estimate the relative abundances of protein families present. The most abundant protein family in January and June venoms is phospholipase A₂ (PLA₂: January 66.7%; June 54.5%), followed by three-finger toxins (3FTx: January 30.4%; June 40.4%) and a minor component of cysteine-rich secretory proteins (CRISP: January 2.5%; June 5%). Trace amounts of snake venom metalloproteinases (SVMP), C-type lectins and housekeeping and regulatory proteins were also found. Although the complexity of the venom is low by number of families present, each family contained a more diverse set of isoforms than previously reported, a finding that may have implications for the development of next-generation sea snake antivenoms. Intraspecific variability was shown to be minor with one obvious exception of a 14,157-Da protein that was present in some January (adult) venoms, but not at all in June (subadult) venoms. There is also a greater abundance of short-chain neurotoxins in June (subadult) venom compared with January (adult) venom. These differences potentially indicate the presence of seasonal, ontogenetic or sexual variation in H. curtus venom.

Oceanic circulation models help to predict global biogeography of pelagic yellow-bellied sea snake

It is well recognized that most marine vertebrates, and especially tetrapods, precisely orient and actively move in apparently homogeneous oceanic environments. Here, we investigate the presumptive role of oceanic currents in biogeographic patterns observed in a secondarily marine tetrapod, the yellow-bellied sea snake (Hydrophis [Pelamis] platurus). State-of-the-art world ocean circulation models show how H. platurus, the only pelagic species of sea snake, can potentially exploit oceanic currents to disperse and maintain population mixing between localities that spread over two-thirds of the Earth's circumference. The very close association of these snakes with surface currents seems to provide a highly efficient dispersal mechanism that allowed this species to range extensively and relatively quickly well beyond the central Indo-Pacific area, the centre of origin, abundance and diversity of sea snakes. Our results further suggest that the pan-oceanic population of this species must be extraordinarily large.

The evolution of scale sensilla in the transition from land to sea in elapid snakes

Scale sensilla are small tactile mechanosensory organs located on the head scales of many squamate reptiles (lizards and snakes). In sea snakes and sea kraits (Elapidae: Hydrophiinae), these scale organs are presumptive scale sensilla that purportedly function as both tactile mechanoreceptors and potentially as hydrodynamic receptors capable of sensing the displacement of water. We combined scanning electron microscopy, silicone casting of the skin and quadrate sampling with a phylogenetic analysis to assess morphological variation in sensilla on the postocular head scale(s) across four terrestrial, 13 fully aquatic and two semi-aquatic species of elapids. Substantial variation exists in the overall coverage of sensilla (0.8-6.5%) among the species sampled and is broadly overlapping in aquatic and terrestrial lineages. However, two observations suggest a divergent, possibly hydrodynamic sensory role of sensilla in sea snake and sea krait species. First, scale sensilla are more protruding (dome-shaped) in aquatic species than in their terrestrial counterparts. Second, exceptionally high overall coverage of sensilla is found only in the fully aquatic sea snakes, and this attribute appears to have evolved multiple times within this group. Our quantification of coverage as a proxy for relative 'sensitivity' represents the first analysis of the evolution of sensilla in the transition from terrestrial to marine habitats. However, evidence from physiological and behavioural studies is needed to confirm the functional role of scale sensilla in sea snakes and sea kraits.

A Review and Database of Snake Venom Proteomes

Advances in the last decade combining transcriptomics with established proteomics methods have made possible rapid identification and quantification of protein families in snake venoms. Although over 100 studies have been published, the value of this information is increased when it is collated, allowing rapid assimilation and evaluation of evolutionary trends, geographical variation, and possible medical implications. This review brings together all compositional studies of snake venom proteomes published in the last decade. Compositional studies were identified for 132 snake species: 42 from 360 (12%) Elapidae (elapids), 20 from 101 (20%) Viperinae (true vipers), 65 from 239 (27%) Crotalinae (pit vipers), and five species of non-front-fanged snakes. Approximately 90% of their total venom composition consisted of eight protein families for elapids, 11 protein families for viperines and ten protein families for crotalines. There were four dominant protein families: phospholipase A₂s (the most common across all front-fanged snakes), metalloproteases, serine proteases and three-finger toxins. There were six secondary protein families: cysteine-rich secretory proteins, l-amino acid oxidases, kunitz peptides, C-type lectins/snaclecs, disintegrins and natriuretic peptides. Elapid venoms contained mostly three-finger toxins and phospholipase A₂s and viper venoms metalloproteases, phospholipase A₂s and serine proteases. Although 63 protein families were identified, more than half were present in <5% of snake species studied and always in low abundance. The importance of these minor component proteins remains unknown.

A new subspecies of sea snake, Hydrophis platurus xanthos, from Golfo Dulce, Costa Rica

We describe a distinctive new subspecies of sea snake from the occasionally anoxic inner-basin waters of Golfo Dulce, Costa Rica, based on combined data garnered between 2010 and 2017 for 154 specimens, 123 free-ranging and 31 museum-held. The yellow sea snake, Hydrophisplaturusxanthos Bessesen & Galbreath, subsp. n., is diagnosed by a notably smaller body size and nearly uniform yellow coloration, which contrasts with the black and yellow striae and tail spots or bands typical of the species. Within the modest geographic range (circa 320 km²), nearly all specimens possess both diagnostic character states. Bathymetrics appear to restrict genetic flow between this allopatric population and conspecifics in the broader Eastern Pacific. In perspicuous contrast to typical H.platurus, H.p.xanthos shows no association with drift lines, and feeds at night in turbulent waters, assuming a sinusoidal ambush posture never previously reported for the species. This evolutionarily significant unit (ESU) warrants taxonomic recognition and active protection.

Postmortem examination of Australian sea snakes ( Hydrophiinae): Anatomy and common pathologic conditions

There is limited published information about disease in wild sea snakes and no standardized guideline for postmortem examination of sea snakes. Identifying causes of morbidity and mortality of marine vertebrate species has been pivotal to understanding disease factors implicated in stranding events and assisting with the formulation of conservation plans. Additionally, postmortem findings can provide valuable information on life history traits and the ecology of these reclusive species. Sick, moribund, or dead sea snakes are intermittently washed ashore along Australian and international beaches and provide an opportunity to examine a subset of the population and identify causes of population decline. We present an illustrated description of sea snake anatomy and describe a systematic approach to postmortem examination of sea snakes. We describe common pathologic conditions identified from clinical and postmortem examinations of stranded Australian sea snakes from southeast Queensland. Notable pathologic conditions include traumatic injury, inflammatory conditions, parasitic infections, and neoplasia.

Ontogenetic shifts of heart position in snakes

Heart position relative to total body length (TL) varies among snakes, with anterior hearts in arboreal species and more centrally located hearts in aquatic or ground-dwelling species. Anterior hearts decrease the cardiac work associated with cranial blood flow and minimize drops in cranial pressure and flow during head-up climbing. Here, we investigate whether heart position shifts intraspecifically during ontogenetic increases in TL. Insular Florida cottonmouth snakes, Agkistrodon conanti, are entirely ground-dwelling and have a mean heart position that is 33.32% TL from the head. In contrast, arboreal rat snakes, Pantherophis obsoleta, of similar lengths have a mean heart position that is 17.35% TL from the head. In both species, relative heart position shifts craniad during ontogeny, with negative slopes = -.035 and -.021% TL/cm TL in Agkistrodon and Pantherophis, respectively. Using a large morphometric data set available for Agkistrodon (N = 192 individuals, 23-140 cm TL), we demonstrate there is an anterior ontogenetic shift of the heart position within the trunk (= 4.56% trunk length from base of head to cloacal vent), independent of head and tail allometry which are both negative. However, in longer snakes > 100 cm, the heart position reverses and shifts caudally in longer Agkistrodon but continues toward the head in longer individuals of Pantherophis. Examination of data sets for two independent lineages of fully marine snakes (Acrochordus granulatus and Hydrophis platurus), which do not naturally experience postural gravity stress, demonstrate both ontogenetic patterns for heart position that are seen in the terrestrial snakes. The anterior migration of the heart is greater in the terrestrial species, even if TL is standardized to that of the longer P. obsoleta, and compensates for about 5 mmHg gravitational pressure head if they are fully upright.

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.

Sea snakes (Elapidae, Hydrophiinae) in their westernmost extent: an updated and illustrated checklist and key to the species in the Persian Gulf and Gulf of Oman

The Persian Gulf is known as the westernmost distribution limit for sea snakes, except for Hydrophis platurus (Linnaeus, 1766) that reaches southeastern Africa. Previous identification guides for sea snakes of the Persian Gulf and its adjacent waters in the Gulf of Oman were based on old data and confined mostly to written descriptions. Therefore, a series of field surveys were carried out in 2013 and 2014 through Iranian coastal waters of both gulfs to provide a comprehensive sampling of sea snakes in the area. This paper presents an illustrated and updated checklist and identification tool for sea snakes in the Persian Gulf and Gulf of Oman, which are based on new material and a review of the literature. This checklist includes ten species of marine hydrophiines, of which one, Microcephalophis cantoris (Günther, 1864), is a new record for the area. All specimens examined herein are deposited and available at the Zoological Museum of Shahid Bahonar University of Kerman, Kerman province, Iran.

A Species-Level Phylogeny of Extant Snakes with Description of a New Colubrid Subfamily and Genus

BACKGROUND

With over 3,500 species encompassing a diverse range of morphologies and ecologies, snakes make up 36% of squamate diversity. Despite several attempts at estimating higher-level snake relationships and numerous assessments of generic- or species-level phylogenies, a large-scale species-level phylogeny solely focusing on snakes has not been completed. Here, we provide the largest-yet estimate of the snake tree of life using maximum likelihood on a supermatrix of 1745 taxa (1652 snake species + 7 outgroup taxa) and 9,523 base pairs from 10 loci (5 nuclear, 5 mitochondrial), including previously unsequenced genera (2) and species (61).

RESULTS

Increased taxon sampling resulted in a phylogeny with a new higher-level topology and corroborate many lower-level relationships, strengthened by high nodal support values (> 85%) down to the species level (73.69% of nodes). Although the majority of families and subfamilies were strongly supported as monophyletic with > 88% support values, some families and numerous genera were paraphyletic, primarily due to limited taxon and loci sampling leading to a sparse supermatrix and minimal sequence overlap between some closely-related taxa. With all rogue taxa and incertae sedis species eliminated, higher-level relationships and support values remained relatively unchanged, except in five problematic clades.

CONCLUSION

Our analyses resulted in new topologies at higher- and lower-levels; resolved several previous topological issues; established novel paraphyletic affiliations; designated a new subfamily, Ahaetuliinae, for the genera Ahaetulla, Chrysopelea, Dendrelaphis, and Dryophiops; and appointed Hemerophis (Coluber) zebrinus to a new genus, Mopanveldophis. Although we provide insight into some distinguished problematic nodes, at the deeper phylogenetic scale, resolution of these nodes may require sampling of more slowly-evolving nuclear genes.

Venomic Analysis of the Poorly Studied Desert Coral Snake, Micrurus tschudii tschudii, Supports the 3FTx/PLA₂ Dichotomy across Micrurus Venoms

The venom proteome of the poorly studied desert coral snake Micrurus tschudii tschudii was unveiled using a venomic approach, which identified ≥38 proteins belonging to only four snake venom protein families. The three-finger toxins (3FTxs) constitute, both in number of isoforms (~30) and total abundance (93.6% of the venom proteome), the major protein family of the desert coral snake venom. Phospholipases A₂ (PLA₂s; seven isoforms, 4.1% of the venom proteome), 1–3 Kunitz-type proteins (1.6%), and 1–2 l-amino acid oxidases (LAO, 0.7%) complete the toxin arsenal of M. t. tschudii. Our results add to the growing evidence that the occurrence of two divergent venom phenotypes, i.e., 3FTx- and PLA₂-predominant venom proteomes, may constitute a general trend across the cladogenesis of Micrurus. The occurrence of a similar pattern of venom phenotypic variability among true sea snake (Hydrophiinae) venoms suggests that the 3FTx/PLA₂ dichotomy may be widely distributed among Elapidae venoms.

Diversification rates and phenotypic evolution in venomous snakes (Elapidae)

The relationship between rates of diversification and of body size change (a common proxy for phenotypic evolution) was investigated across Elapidae, the largest radiation of highly venomous snakes. Time-calibrated phylogenetic trees for 175 species of elapids (more than 50% of known taxa) were constructed using seven mitochondrial and nuclear genes. Analyses using these trees revealed no evidence for a link between speciation rates and changes in body size. Two clades (Hydrophis, Micrurus) show anomalously high rates of diversification within Elapidae, yet exhibit rates of body size evolution almost identical to the general elapid 'background' rate. Although correlations between speciation rates and rates of body size change exist in certain groups (e.g. ray-finned fishes, passerine birds), the two processes appear to be uncoupled in elapid snakes. There is also no detectable shift in diversification dynamics associated with the colonization of Australasia, which is surprising given that elapids appear to be the first clade of venomous snakes to reach the continent.