Multilocus phylogeography of the sea snakeHydrophis curtusreveals historical vicariance and cryptic lineage diversity

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.

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.

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.

A comparative analysis of the proteomes and biological activities of the venoms from two sea snakes, Hydrophis curtus and Hydrophis cyanocinctus, from Hainan, China

We characterized and compared the venom protein profiles of Hydrophis curtus (synonyms: Lapemis hardwickii, Lapemis curtus and Hydrophis hardwickii) and Hydrophis cyanocinctus, the two representatives of medically important venomous sea snakes in Chinese waters using proteomic approaches. A total of 47 and 38 putative toxins were identified in H. curtus venom (HcuV) and H. cyanocinctus venom (HcyV), respectively, and these toxins could be grouped into 15 functional categories, mainly proteinases, phospholipases, three-finger toxins (3FTxs), lectins, protease inhibitors, ion channel inhibitors, cysteine-rich venom proteins (CRVPs) and snake venom metalloproteases (SVMPs). The constituent ratio of each toxin category varied between HcuV and HcyV with 3FTx (54% in HcuV/69% in HcyV) and PLA2 (38% in HcuV/22% in HcyV) unanimously ranked as the top two most abundant families. Both HcuV and HcyV exhibited relatively high lethality (LD₅₀ values in mice of 0.34 μg/g and 0.24 μg/g, respectively), specific PLA2 activity and hemolytic activity. On the basis of several previous reports of HcuV and HcyV collected from other areas, these findings greatly expand our understanding of geographical variation and interspecies diversity of the two sea snake venoms and can provide a scientific basis for the development of specific sea snake antivenom in the future.

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.

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.