Comparative Venom Multiomics Reveal the Molecular Mechanisms Driving Adaptation to Diverse Predator-Prey Ecosystems in Closely Related Sea Snakes.

Predator-prey arms races are ideal models for studying the natural selection and adaptive evolution that drive the formation of biological diversity. For venomous snakes, venom is a key bridge linking snakes with their prey, but whether and how venom evolves under the selection of diet remains unclear. Here, we focused on two closely related sea snakes, Hydrophis cyanocinctus and Hydrophis curtus, which show significant differences in prey preferences. Data-independent acquisition (DIA)-based proteomic analysis revealed different degrees of homogeneity in the venom composition of the two snakes, which was consistent with the differential phylogenetic diversity of their prey. By investigating the sequences and structures of three-finger toxins (3FTx), a predominant toxin family in elapid venom, we identified significant differences between the two sea snakes in the binding activity of 3FTx to receptors from different prey populations, which could explain the trophic specialization of H. cyanocinctus. Furthermore, we performed integrated multiomic profiling of the transcriptomes, microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and proteomes of the venom glands; constructed venom-related mRNA-miRNA-lncRNA networks; and identified a series of noncoding RNAs involved in the regulation of toxin gene expression in the two species. These findings are highly informative for elucidating the molecular basis and regulatory mechanisms that account for discrepant venom evolution in response to divergent diets in closely related snakes, providing valuable evidence for the study of coselection and coevolution in predator-prey ecosystems.

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.

Marine Envenomation in Okinawa: Overview and Treatment Concept.

Okinawa prefecture is a popular tourist destination due to its beaches and reefs. The reefs host a large variety of animals, including a number of venomous species. Because of the popularity of the reefs and marine activities, people are frequently in close contact with dangerous venomous species and, thus, are exposed to potential envenomation. Commonly encountered venomous animals throughout Okinawa include the invertebrate cone snail, sea urchin, crown-of-thorns starfish, blue-ringed octopus, box jellyfish, and fire coral. The vertebrates include the stonefish, lionfish, sea snake, and moray eel. Treatment for marine envenomation can involve first aid, hot water immersion, antivenom, supportive care, regional anesthesia, and pharmaceutical administration. Information on venomous animals, their toxins, and treatment should be well understood by prehospital care providers and physicians practicing in the prefecture.

A Bite by Shaw's Short Sea Snake (Hydrophis curtus): A Case of Mild Myotoxicity or a Dry Bite?

Although sea snakes (Elapidae) are commonly encountered by fishermen, accurately authenticated envenomings by them are uncommon in clinical literature. We report an authenticated case of Shaw's short, or spine-bellied, sea snake (Hydrophis curtus) bite in a young fisherman from northern Sri Lanka. The patient had clinical and biochemical evidence of mild transient myotoxicity but no evidence of neuromuscular paralysis or significant renal injury. Consideration of the clinical manifestations suggests either a mild envenoming or a dry bite. The patient completely recovered without any antivenom therapy and was discharged on the fourth day. Prolonged observation may be beneficial to exclude complications of sea snake envenoming.

A cathelicidin antimicrobial peptide from Hydrophis cyanocinctus inhibits Zika virus infection by downregulating expression of a viral entry factor.

Zika virus (ZIKV) is a re-emerging flavivirus that causes conditions such as microcephaly and testis damage. The spread of ZIKV has become a major public health concern. Recent studies indicated that antimicrobial peptides are an ideal source for screening antiviral candidates with broad-spectrum antiviral activities, including against ZIKV. We herein found that Hc-CATH, a cathelicidin antimicrobial peptide identified from the sea snake Hydrophis cyanocinctus in our previous work, conferred protection against ZIKV infection in host cells and showed preventative efficacy and therapeutic efficacy in C57BL/6J mice, Ifnar1-/- mice, and pregnant mice. Intriguingly, we revealed that Hc-CATH decreased the susceptibility of host cells to ZIKV by downregulating expression of AXL, a TAM (TYRO3, AXL and MERTK) family kinase receptor that mediates ZIKV infection, and subsequently reversed the negative regulation of AXL on host's type I interferon response. Furthermore, we showed that the cyclo-oxygenase-2/prostaglandin E2/adenylyl cyclase/protein kinase A pathway was involved in Hc-CATH-mediated AXL downregulation, and Hc-CATH in addition directly inactivated ZIKV particles by disrupting viral membrane. Finally, while we found Hc-CATH did not act on the late stage of ZIKV infection, structure-function relationship studies revealed that α-helix and phenylalanine residues are key structural requirements for its protective efficacy against initial ZIKV infection. In summary, we demonstrate that Hc-CATH provides prophylactic and therapeutic efficacy against ZIKV infection via downregulation of AXL, as well as inactivating the virion. Our findings reveal a novel mechanism of cathelicidin against viral infection and highlight the potential of Hc-CATH to prevent and treat ZIKV infection.

Identification and Characterization of a Novel Cathelicidin from Hydrophis cyanocinctus with Antimicrobial and Anti-Inflammatory Activity.

The abuse of antibiotics and lack of new antibacterial drugs has led to the emergence of superbugs that raise fears of untreatable infections. The Cathelicidin family of antimicrobial peptide (AMP) with varying antibacterial activities and safety is considered to be a promising alternative to conventional antibiotics. In this study, we investigated a novel Cathelicidin peptide named Hydrostatin-AMP2 from the sea snake Hydrophis cyanocinctus. The peptide was identified based on gene functional annotation of the H. cyanocinctus genome and bioinformatic prediction. Hydrostatin-AMP2 showed excellent antimicrobial activity against both Gram-positive and Gram-negative bacteria, including standard and clinical Ampicillin-resistant strains. The results of the bacterial killing kinetic assay demonstrated that Hydrostatin-AMP2 had faster antimicrobial action than Ampicillin. Meanwhile, Hydrostatin-AMP2 exhibited significant anti-biofilm activity including inhibition and eradication. It also showed a low propensity to induce resistance as well as low cytotoxicity and hemolytic activity. Notably, Hydrostatin-AMP2 apparently decreased the production of pro-inflammatory cytokines in the LPS-induced RAW264.7 cell model. To sum up, these findings indicate that Hydrostatin-AMP2 is a potential peptide candidate for the development of new-generation antimicrobial drugs fighting against antibiotic-resistant bacterial infections.

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.

The banded colour patterns of sea snakes discourage attack by predatory fishes, enabling Batesian mimicry by harmless species.

The evolution of bright 'warning' colours in nontoxic animals often is attributed to mimicry of toxic species, but empirical tests of that hypothesis must overcome the logistical challenge of quantifying differential rates of predation in nature. Populations of a harmless sea snake species (Emydocephalus annulatus) in New Caledonia exhibit colour polymorphism, with around 20% of individuals banded rather than melanic. Stability in that proportion over 20 years has been attributed to Batesian mimicry of deadly snake species by banded morphs of the harmless taxon. This hypothesis requires that banded colours reduce a snake's vulnerability to predation. We tested that idea by pulling flexible snake-shaped models through the water and recording responses by predatory fish. Black and banded lures attracted similar numbers of following fish, but attacks were directed almost exclusively to black lures. Our methods overcome several ambiguities associated with experimental studies on mimicry in terrestrial snakes and support the hypothesis that banded colour patterns reduce a non-venomous marine snake's vulnerability to predation.

Venom-gland transcriptomic, venomic, and antivenomic profiles of the spine-bellied sea snake (Hydrophis curtus) from the South China Sea.

BACKGROUND: A comprehensive evaluation of the -omic profiles of venom is important for understanding the potential function and evolution of snake venom. Here, we conducted an integrated multi-omics-analysis to unveil the venom-transcriptomic and venomic profiles in a same group of spine-bellied sea snakes (Hydrophis curtus) from the South China Sea, where the snake is a widespread species and might generate regionally-specific venom potentially harmful to human activities. The capacity of two heterologous antivenoms to immunocapture the H. curtus venom was determined for an in-depth evaluation of their rationality in treatment of H. curtus envenomation. In addition, a phylogenetic analysis by maximum likelihood was used to detect the adaptive molecular evolution of full-length toxin-coding unigenes. RESULTS: A total of 90,909,384 pairs of clean reads were generated via Illumina sequencing from a pooled cDNA library of six specimens, and yielding 148,121 unigenes through de novo assembly. Sequence similarity searching harvested 63,845 valid annotations, including 63,789 non-toxin-coding and 56 toxin-coding unigenes belonging to 22 protein families. Three protein families, three-finger toxins (3-FTx), phospholipase A2 (PLA2), and cysteine-rich secretory protein, were detected in the venom proteome. 3-FTx (27.15% in the transcriptome/41.94% in the proteome) and PLA2 (59.71%/49.36%) were identified as the most abundant families in the venom-gland transcriptome and venom proteome. In addition, 24 unigenes from 11 protein families were shown to have experienced positive selection in their evolutionary history, whereas four were relatively conserved throughout evolution. Commercial Naja atra antivenom exhibited a stronger capacity than Bungarus multicinctus antivenom to immunocapture H. curtus venom components, especially short neurotoxins, with the capacity of both antivenoms to immunocapture short neurotoxins being weaker than that for PLA2s. CONCLUSIONS: Our study clarified the venom-gland transcriptomic and venomic profiles along with the within-group divergence of a H. curtus population from the South China Sea. Adaptive evolution of most venom components driven by natural selection appeared to occur rapidly during evolutionary history. Notably, the utility of commercial N. atra and B. multicinctus antivenoms against H. curtus toxins was not comprehensive; thus, the development of species-specific antivenom is urgently needed.

The Genome of Shaw's Sea Snake (Hydrophis curtus) Reveals Secondary Adaptation to Its Marine Environment.

The transition of terrestrial snakes to marine life ∼10 Ma is ideal for exploring adaptive evolution. Sea snakes possess phenotype specializations including laterally compressed bodies, paddle-shaped tails, valvular nostrils, cutaneous respiration, elongated lungs, and salt glands, yet, knowledge on the genetic underpinnings of the transition remains limited. Herein, we report the first genome of Shaw's sea snake (Hydrophis curtus) and use it to investigate sea snake secondary marine adaptation. A hybrid assembly strategy obtains a high-quality genome. Gene family analyses date a pulsed coding-gene expansion to ∼20 Ma, and these genes associate strongly with adaptations to marine environments. Analyses of selection pressure and convergent evolution discover the rapid evolution of protein-coding genes, and some convergent features. Additionally, 108 conserved noncoding elements appear to have evolved quickly, and these may underpin the phenotypic changes. Transposon elements may contribute to adaptive specializations by inserting into genomic regions around functionally related coding genes. The integration of genomic and transcriptomic analyses indicates independent origins and different components in sea snake and terrestrial snake venom; the venom gland of the sea snake harbors the highest PLA2 (17.23%) expression in selected elapids and these genes may organize tandemly in the genome. These analyses provide insights into the genetic mechanisms that underlay the secondary adaptation to marine and venom production of this sea snake.

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 (Ne ) 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 Ne and less pronounced fluctuations in Ne , 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 Ne , while each terrestrial species shows a prominent upturn in Ne 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 Ne 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.

Fatal Sea Snake Envenomation.

ABSTRACT: A 23-year-old man working on a prawn trawler off the Northern Australian coast was bitten on the right hand by a black-banded sea snake (Laticauda colubrina), resulting in the rapid onset of ptosis, blurred vision, and respiratory difficulties with convulsions, loss of consciousness, and death. Resuscitation was unsuccessful. No antivenom was available because of the remote location. At autopsy, a small puncture wound was identified on the middle interphalangeal joint of the right little finger. No other significant abnormalities were detected. Elevations in creatine kinase (1558 µ/L) and myoglobin (24,800 µg/L) levels were found. Death was caused by respiratory failure after sea snake envenomation. This case demonstrates potential difficulties in establishing the diagnosis of lethal sea snake envenomation, including the inconspicuous nature of the bite marks with no local reaction and the underlying neurotoxic nature of the lethal mechanisms, which leaves very nonspecific pathological findings. Myonecrosis may occur but is variable. Unless reliably witnessed, sea snake envenomation may not be suspected. Failing vision or unexplained paresis in fishermen or individuals who have been in the sea may point to the diagnosis.

[Identification of Chinese medicine sea snake based on cytochrome B (Cytb)].

The identification of species primordium has been one of the hot issues in the identification of traditional Chinese medicine. Sea snake is one of the most valuable Chinese medicinal materials in China. In order to understand the origin and varieties of sea snake in the market, we studied the molecular identification of 46 sea snakes by cytochrome B(Cytb). After comparison and manual correction, the sequence length was 582 bp, and the content of A+T(58.9%) was higher than that of G+C(41.1%). There exist 197 variable sites and 179 parsimony-informative sites of the sequence. There are 44 kinds of sequence alignment with consistency equal to 100%, and 2 kinds equal to 96%. A total of 408 Cytb effective sequences were downloaded from GenBank database, with a total of 68 species. Phylogenetic tree of a total of 454 sea snake sequences with the samples in this study were constructed by neighbor-joining trees and Bayesian inference method, respectively, which can identify 42 samples of medicinal materials, while 4 samples can not be identified because of their low node support. The results showed that the species of the sea snake medicine were at least from 2 genera and 5 species, namely, Aipysurus eydouxii, Hydrophis curtus, H. caerulescen, H. curtus, H. ornatus and H. spiralis. This study suggested that the original species of commercial sea snake are very complex and can provide insight into the identification of sea snakes.

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.

Patterns of intracolumnar size variation inform the heterochronic mechanisms underlying extreme body shape divergence in microcephalic sea snakes.

Sea snakes (Hydrophiinae) that specialize on burrowing eel prey have repeatedly evolved tiny heads and reduced forebody relative to hindbody girths. Previous research has found that these "microcephalic" forms have higher counts of precaudal vertebrae, and postnatal ontogenetic changes cause their hindbodies to reach greater girths relative to their forebodies. We examine variation in vertebral size along the precaudal axis of neonates and adults of three species. In the nonmicrocephalic Hydrophis curtus, these intracolumnar patterns take the form of symmetrical curved profiles, with longer vertebrae in the midbody (50% of body length) relative to distal regions. In contrast, intracolumnar profiles in the microcephalic H. macdowelli and H. obscurus are strongly asymmetrical curves (negative skewness) due to the presence of numerous, smaller-sized vertebrate in the forebody (anterior to the heart). Neonate and adult H. macdowelli and H. obscurus specimens all exhibit this pattern, implying an onset of fore- versus hindbody decoupling in the embryo stage. Based on this, we suggest plausible developmental mechanisms involving the presence and positioning of Hox boundaries and heterochronic changes in segmentation. Tests of our hypotheses would give new insights into the drivers of rapid convergent shifts in evolution, but will ultimately require studies of gene expression in the embryos of relevant taxa.

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.

Loss of olfaction in sea snakes provides new perspectives on the aquatic adaptation of amniotes.

Marine amniotes, a polyphyletic group, provide an excellent opportunity for studying convergent evolution. Their sense of smell tends to degenerate, but this process has not been explored by comparing fully aquatic species with their amphibious relatives in an evolutionary context. Here, we sequenced the genomes of fully aquatic and amphibious sea snakes and identified repertoires of chemosensory receptor genes involved in olfaction. Snakes possess large numbers of the olfactory receptor (OR) genes and the type-2 vomeronasal receptor (V2R) genes, and expression profiling in the olfactory tissues suggests that snakes use the ORs in the main olfactory system (MOS) and the V2Rs in the vomeronasal system (VNS). The number of OR genes has decreased in sea snakes, and fully aquatic species lost MOS which is responsible for detecting airborne odours. By contrast, sea snakes including fully aquatic species retain a number of V2R genes and a well-developed VNS for smelling underwater. This study suggests that the sense of smell also degenerated in sea snakes, particularly in fully aquatic species, but their residual olfactory capability is distinct from that of other fully aquatic amniotes. Amphibious species show an intermediate status between terrestrial and fully aquatic snakes, implying their importance in understanding the process of aquatic adaptation.

Phototactic tails: Evolution and molecular basis of a novel sensory trait in sea snakes.

Dermal phototaxis has been reported in a few aquatic vertebrate lineages spanning fish, amphibians and reptiles. These taxa respond to light on the skin of their elongate hind-bodies and tails by withdrawing under cover to avoid detection by predators. Here, we investigated tail phototaxis in sea snakes (Hydrophiinae), the only reptiles reported to exhibit this sensory behaviour. We conducted behavioural tests in 17 wild-caught sea snakes of eight species by illuminating the dorsal surface of the tail and midbody skin using cold white, violet, blue, green and red light. Our results confirmed phototactic tail withdrawal in the previously studied Aipysurus laevis, revealed this trait for the first time in A. duboisii and A. tenuis, and suggested that tail photoreceptors have peak spectral sensitivities between blue and green light (457-514 nm). Based on these results, and an absence of photoresponses in five Aipysurus and Hydrophis species, we tentatively infer that tail phototaxis evolved in the ancestor of a clade of six Aipysurus species (comprising 10% of all sea snakes). Quantifying tail damage, we found that the probability of sustaining tail injuries was not influenced by tail phototactic ability in snakes. Gene profiling showed that transcriptomes of both tail skin and body skin lacked visual opsins but contained melanopsin (opn4x) in addition to key genes of the retinal regeneration and phototransduction cascades. This work suggests that a nonvisual photoreceptor (e.g., Gq rhabdomeric) signalling pathway underlies tail phototaxis, and provides candidate gene targets for future studies of this unusual sensory innovation in reptiles.

Underwater hearing in sea snakes (Hydrophiinae): first evidence of auditory evoked potential thresholds.

The viviparous sea snakes (Hydrophiinae) are a secondarily aquatic radiation of more than 60 species that possess many phenotypic adaptations to marine life. However, virtually nothing is known of the role and sensitivity of hearing in sea snakes. This study investigated the hearing sensitivity of the fully marine sea snake Hydrophis stokesii by measuring auditory evoked potential (AEP) audiograms for two individuals. AEPs were recorded from 40 Hz (the lowest frequency tested) up to 600 Hz, with a peak in sensitivity identified at 60 Hz (163.5 dB re. 1 µPa or 123 dB re. 1 µm s-2). Our data suggest that sea snakes are sensitive to low-frequency sounds but have relatively low sensitivity compared with bony fishes and marine turtles. Additional studies are required to understand the role of sound in sea snake life history and further assess these species' vulnerability to anthropogenic noise.

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.