A global reptile assessment highlights shared conservation needs of tetrapods.

Comprehensive assessments of species' extinction risks have documented the extinction crisis1 and underpinned strategies for reducing those risks2. Global assessments reveal that, among tetrapods, 40.7% of amphibians, 25.4% of mammals and 13.6% of birds are threatened with extinction3. Because global assessments have been lacking, reptiles have been omitted from conservation-prioritization analyses that encompass other tetrapods4-7. Reptiles are unusually diverse in arid regions, suggesting that they may have different conservation needs6. Here we provide a comprehensive extinction-risk assessment of reptiles and show that at least 1,829 out of 10,196 species (21.1%) are threatened-confirming a previous extrapolation8 and representing 15.6 billion years of phylogenetic diversity. Reptiles are threatened by the same major factors that threaten other tetrapods-agriculture, logging, urban development and invasive species-although the threat posed by climate change remains uncertain. Reptiles inhabiting forests, where these threats are strongest, are more threatened than those in arid habitats, contrary to our prediction. Birds, mammals and amphibians are unexpectedly good surrogates for the conservation of reptiles, although threatened reptiles with the smallest ranges tend to be isolated from other threatened tetrapods. Although some reptiles-including most species of crocodiles and turtles-require urgent, targeted action to prevent extinctions, efforts to protect other tetrapods, such as habitat preservation and control of trade and invasive species, will probably also benefit many reptiles.

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.

First evidence of hemiclitores in snakes.

Female genitalia are conspicuously overlooked in comparison to their male counterparts, limiting our understanding of sexual reproduction across vertebrate lineages. This study is the first complete description of the clitoris (hemiclitores) in female snakes. We describe morphological variation in size and shape (n = 9 species, 4 families) that is potentially comparable to the male intromittent organs in squamate reptiles (hemipenes). Dissection, diffusible iodine contrast-enhanced micro-CT and histology revealed that, unlike lizard hemiclitores, the snake hemiclitores are non-eversible structures. The two individual hemiclitores are separated medially by connective tissue, forming a triangular structure that extends posteriorly. Histology of the hemiclitores in Australian death adders (Acanthophis antarcticus) showed erectile tissue and strands/bundles of nerves, but no spines (as is found in male hemipenes). These histological features suggest the snake hemiclitores have functional significance in mating and definitively show that the hemiclitores are not underdeveloped hemipenes or scent glands, which have been erroneously indicated in other studies. Our discovery supports that hemiclitores have been retained across squamates and provides preliminary evidence of differences in this structure among snake species, which can be used to further understand systematics, reproductive evolution and ecology across squamate reptiles.

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.

Horizontal transfer and subsequent explosive expansion of a DNA transposon in sea kraits (Laticauda).

Transposable elements (TEs) are self-replicating genetic sequences and are often described as important 'drivers of evolution'. This driving force is because TEs promote genomic novelty by enabling rearrangement, and through exaptation as coding and regulatory elements. However, most TE insertions potentially lead to neutral or harmful outcomes, therefore host genomes have evolved machinery to suppress TE expansion. Through horizontal transposon transfer (HTT) TEs can colonize new genomes, and since new hosts may not be able to regulate subsequent replication, these TEs may proliferate rapidly. Here, we describe HTT of the Harbinger-Snek DNA transposon into sea kraits (Laticauda), and its subsequent explosive expansion within Laticauda genomes. This HTT occurred following the divergence of Laticauda from terrestrial Australian elapids approximately 15-25 Mya. This has resulted in numerous insertions into introns and regulatory regions, with some insertions into exons which appear to have altered UTRs or added sequence to coding exons. Harbinger-Snek has rapidly expanded to make up 8-12% of Laticauda spp. genomes; this is the fastest known expansion of TEs in amniotes following HTT. Genomic changes caused by this rapid expansion may have contributed to adaptation to the amphibious-marine habitat.

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.

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.

Multiple rod-cone and cone-rod photoreceptor transmutations in snakes: evidence from visual opsin gene expression.

In 1934, Gordon Walls forwarded his radical theory of retinal photoreceptor 'transmutation'. This proposed that rods and cones used for scotopic and photopic vision, respectively, were not fixed but could evolve into each other via a series of morphologically distinguishable intermediates. Walls' prime evidence came from series of diurnal and nocturnal geckos and snakes that appeared to have pure-cone or pure-rod retinas (in forms that Walls believed evolved from ancestors with the reverse complement) or which possessed intermediate photoreceptor cells. Walls was limited in testing his theory because the precise identity of visual pigments present in photoreceptors was then unknown. Subsequent molecular research has hitherto neglected this topic but presents new opportunities. We identify three visual opsin genes, rh1, sws1 and lws, in retinal mRNA of an ecologically and taxonomically diverse sample of snakes central to Walls' theory. We conclude that photoreceptors with superficially rod- or cone-like morphology are not limited to containing scotopic or photopic opsins, respectively. Walls' theory is essentially correct, and more research is needed to identify the patterns, processes and functional implications of transmutation. Future research will help to clarify the fundamental properties and physiology of photoreceptors adapted to function in different light levels.

Dynamic Expansions and Retinal Expression of Spectrally Distinct Short-Wavelength Opsin Genes in Sea Snakes.

The photopigment-encoding visual opsin genes that mediate colour perception show great variation in copy number and adaptive function across vertebrates. An open question is how this variation has been shaped by the interaction of lineage-specific structural genomic architecture and ecological selection pressures. We contribute to this issue by investigating the expansion dynamics and expression of the duplicated Short-Wavelength-Sensitive-1 opsin (SWS1) in sea snakes (Elapidae). We generated one new genome, 45 resequencing datasets, 10 retinal transcriptomes, and 81 SWS1 exon sequences for sea snakes, and analysed these alongside 16 existing genomes for sea snakes and their terrestrial relatives. Our analyses revealed multiple independent transitions in SWS1 copy number in the marine Hydrophis clade, with at least three lineages having multiple intact SWS1 genes: the previously studied Hydrophis cyanocinctus and at least two close relatives of this species; H. atriceps-H. fasciatus; and an individual H. curtus. In each lineage, gene copy divergence at a key spectral tuning site resulted in distinct UV and Violet/Blue-sensitive SWS1 subtypes. Both spectral variants were simultaneously expressed in the retinae of H. cyanocinctus and H. atriceps, providing the first evidence that these SWS1 expansions confer novel phenotypes. Finally, chromosome annotation for nine species revealed shared structural features in proximity to SWS1 regardless of copy number. If these features are associated with SWS1 duplication, expanded opsin complements could be more common in snakes than is currently recognised. Alternatively, selection pressures specific to aquatic environments could favour improved chromatic distinction in just some lineages.

Anonymous nuclear loci in non-model organisms: making the most of high-throughput genome surveys.

MOTIVATION: When working with non-model organisms, few if any species-specific markers are available for phylogenetic, phylogeographic and population studies. Therefore, researchers often try to adapt markers developed in distantly related taxa, resulting in poor amplification and ascertainment bias in their target taxa. Markers can be developed de novo and anonymous nuclear loci (ANL) are proving to be a boon for researchers seeking large numbers of fast-evolving, independent loci. However, the development of ANL can be laboratory intensive and expensive. A workflow is described to identify suitable low-copy anonymous loci from high-throughput shotgun sequences, dramatically reducing the cost and time required to develop these markers and produce robust multilocus datasets. RESULTS: By successively removing repetitive and evolutionary conserved sequences from low coverage shotgun libraries, we were able to isolate thousands of potential ANL. Empirical testing of loci developed from two reptile taxa confirmed that our methodology yields markers with comparable amplification rates and nucleotide diversities to ANLs developed using other methodologies. Our approach capitalizes on next-generation sequencing technologies to enable the development of phylogenetic, phylogeographic and population markers for taxa lacking suitable genomic resources.

Recent rapid speciation and ecomorph divergence in Indo-Australian sea snakes.

The viviparous sea snakes (Hydrophiinae) are a young radiation of at least 62 species that display spectacular morphological diversity and high levels of local sympatry. To shed light on the mechanisms underlying sea snake diversification, we investigated recent speciation and eco-morphological differentiation in a clade of four nominal species with overlapping ranges in Southeast Asia and Australia. Analyses of morphology and stomach contents identified the presence of two distinct ecomorphs: a 'macrocephalic' ecomorph that reaches >2 m in length, has a large head and feeds on crevice-dwelling eels and gobies; and a 'microcephalic' ecomorph that rarely exceeds 1 m in length, has a small head and narrow fore-body and hunts snake eels in burrows. Mitochondrial sequences show a lack of reciprocal monophyly between ecomorphs and among putative species. However, individual assignment based on newly developed microsatellites separated co-distributed specimens into four significantly differentiated clusters corresponding to morphological species designations, indicating limited recent gene flow and progress towards speciation. A coalescent species tree (based on mitochondrial and nuclear sequences) and isolation-migration model (mitochondrial and microsatellite markers) suggest between one and three transitions between ecomorphs within the last approximately 1.2 million to approximately 840,000 years. In particular, the macrocephalic 'eastern' population of Hydrophis cyanocinctus and microcephalic H. melanocephalus appear to have diverged very recently and rapidly, resulting in major phenotypic differences and restriction of gene flow in sympatry. These results highlight the viviparous sea snakes as a promising system for speciation studies in the marine environment.

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

The viviparous sea snakes (Hydrophiinae: Hydrophiini) comprise a young but morphologically and ecologically diverse clade distributed throughout the Indo-Pacific. Despite presenting a very promising model for marine diversification studies, many relationships among the 62 species and 16 genera in Hydrophiini remain unresolved. Here, we extend previous taxonomic and genomic sampling for Hydrophiini using three mitochondrial fragments and five nuclear loci for multiple individuals of 39 species in 15 genera. Our results highlight many of the impediments to inferring phylogenies in recent rapid radiations, including low variation at all five nuclear markers, and conflicting relationships supported by mitochondrial and nuclear trees. However, concatenated Bayesian and likelihood analyses, and a multilocus coalescent tree, recovered concordant support for primary clades and several previously unresolved inter-specific groupings. The Aipysurus group is monophyletic, with egg-eating specialists forming separate, early-diverging lineages. All three monotypic semi-aquatic genera (Ephalophis, Parahydrophis and Hydrelaps) are robustly placed as early diverging lineages along the branch leading to the Hydrophis group, with Ephalophis recovered as sister to Parahydrophis. The molecular phylogeny implies extensive evolutionary convergence in feeding adaptations within the Hydrophis group, especially the repeated evolution of small-headed (microcephalic) forms. Microcephalophis (Hydrophis) gracilis is robustly recovered as a relatively distant sister lineage to all other sampled Hydrophis group species, here termed the 'core Hydrophis group'. Within the 'core Hydrophis group', Hydrophis is recovered as broadly paraphyletic, with several other genera nested within it (Pelamis, Enhydrina, Astrotia, Thalassophina, Acalyptophis, Kerilia, Lapemis, Disteira). Instead of erecting multiple new genera, we recommend dismantling the latter (mostly monotypic) genera and recognising a single genus, Hydrophis Latreille 1802, for the core Hydrophis group. Estimated divergence times suggest that all Hydrophiini last shared a common ancestor ∼6million years ago, but that the majority of extant lineages diversified over the last ∼3.5million years. The core Hydrophis group is a young and rapidly speciating clade, with 26 sampled species and 9 genera and dated at only ∼1.5-3million years old.

Molecular evidence that the deadliest sea snake Enhydrina schistosa (Elapidae: Hydrophiinae) consists of two convergent species.

We present a striking case of phenotypic convergence within the speciose and taxonomically unstable Hydrophis group of viviparous sea snakes. Enhydrina schistosa, the 'beaked sea snake', is abundant in coastal and inshore habitats throughout the Asian and Australian regions, where it is responsible for the large majority of recorded deaths and injuries from sea snake bites. Analyses of five independent mitochondrial and nuclear loci for populations spanning Australia, Indonesia and Sri Lanka indicate that this 'species' actually consists of two distinct lineages in Asia and Australia that are not closest relatives. As a result, Australian "E. schistosa" are elevated to species status and provisionally referred to Enhydrinazweifeli. Convergence in the characteristic 'beaked' morphology of these species is probably associated with the wide gape required to accommodate their spiny prey. Our findings have important implications for snake bite management in light of the medical importance of beaked sea snakes and the fact that the only sea snake anti-venom available is raised against Malaysian E. schistosa.

Functional Duplication of the Short-Wavelength-Sensitive Opsin in Sea Snakes: Evidence for Reexpanded Color Sensitivity Following Ancestral Regression.

Color vision is mediated by ancient and spectrally distinct cone opsins. Yet, while there have been multiple losses of opsin genes during the evolution of tetrapods, evidence for opsin gains via functional duplication is extremely scarce. Previous studies have shown that some secondarily marine elapid snakes have acquired expanded "UV-blue" sensitivity via changes at key spectral tuning amino acid sites of the Short-Wavelength Opsin 1 (SWS1) gene. Here, we use elapid reference genomes to show that the molecular origin of this adaptation involved repeated, proximal duplications of the SWS1 gene in the fully marine Hydrophis cyanocinctus. This species possesses four intact SWS1 genes; two of these genes have the ancestral UV sensitivity, and two have a derived sensitivity to the longer wavelengths that dominate marine habitats. We suggest that this remarkable expansion of the opsin repertoire of sea snakes functionally compensates for the ancestral losses of two middle-wavelength opsins in the earliest (dim-light adapted) snakes. This provides a striking contrast to the evolution of opsins during ecological transitions in mammals. Like snakes, early mammals lost two cone photopigments; however, lineages such as bats and cetaceans underwent further opsin losses during their adaptation to dim-light environments.

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.

Horizontal Transposon Transfer and Its Implications for the Ancestral Ecology of Hydrophiine Snakes.

Transposable elements (TEs), also known as jumping genes, are sequences able to move or copy themselves within a genome. As TEs move throughout genomes they often act as a source of genetic novelty, hence understanding TE evolution within lineages may help in understanding environmental adaptation. Studies into the TE content of lineages of mammals such as bats have uncovered horizontal transposon transfer (HTT) into these lineages, with squamates often also containing the same TEs. Despite the repeated finding of HTT into squamates, little comparative research has examined the evolution of TEs within squamates. Here we examine a diverse family of Australo-Melanesian snakes (Hydrophiinae) to examine if the previously identified, order-wide pattern of variable TE content and activity holds true on a smaller scale. Hydrophiinae diverged from Asian elapids ~30 Mya and have since rapidly diversified into six amphibious, ~60 marine and ~100 terrestrial species that fill a broad range of ecological niches. We find TE diversity and expansion differs between hydrophiines and their Asian relatives and identify multiple HTTs into Hydrophiinae, including three likely transferred into the ancestral hydrophiine from fish. These HTT events provide the first tangible evidence that Hydrophiinae reached Australia from Asia via a marine route.