Fossil-informed biogeographic analysis suggests Eurasian regionalization in crown Squamata during the early Jurassic

Background

Squamata (lizards, snakes, and amphisbaenians) is a Triassic lineage with an extensive and complex biogeographic history, yet no large-scale study has reconstructed the ancestral range of early squamate lineages. The fossil record indicates a broadly Pangaean distribution by the end- Cretaceous, though many lineages (e.g., Paramacellodidae, Mosasauria, Polyglyphanodontia) subsequently went extinct. Thus, the origin and occupancy of extant radiations is unclear and may have been localized within Pangaea to specific plates, with potential regionalization to distinct Laurasian and Gondwanan landmasses during the Mesozoic in some groups.

Methods

We used recent tectonic models to code extant and fossil squamate distributions occurring on nine discrete plates for 9,755 species, with Jurassic and Cretaceous fossil constraints from three extinct lineages. We modeled ancestral ranges for crown Squamata from an extant-only molecular phylogeny using a suite of biogeographic models accommodating different evolutionary processes and fossil-based node constraints from known Jurassic and Cretaceous localities. We hypothesized that the best-fit models would not support a full Pangaean distribution (i.e., including all areas) for the origin of crown Squamata, but would instead show regionalization to specific areas within the fragmenting supercontinent, likely in the Northern Hemisphere where most early squamate fossils have been found.

Results

Incorporating fossil data reconstructs a localized origin within Pangaea, with early regionalization of extant lineages to Eurasia and Laurasia, while Gondwanan regionalization did not occur until the middle Cretaceous for Alethinophidia, Scolecophidia, and some crown Gekkotan lineages. While the Mesozoic history of extant squamate biogeography can be summarized as a Eurasian origin with dispersal out of Laurasia into Gondwana, their Cenozoic history is complex with multiple events (including secondary and tertiary recolonizations) in several directions. As noted by previous authors, squamates have likely utilized over-land range expansion, land-bridge colonization, and trans-oceanic dispersal. Tropical Gondwana and Eurasia hold more ancient lineages than the Holarctic (Rhineuridae being a major exception), and some asymmetries in colonization (e.g., to North America from Eurasia during the Cenozoic through Beringia) deserve additional study. Future studies that incorporate fossil branches, rather than as node constraints, into the reconstruction can be used to explore this history further.

Phylogenomics of Psammodynastes and Buhoma (Elapoidea: Serpentes), with the description of a new Asian snake family

Asian mock vipers of the genus Psammodynastes and African forest snakes of the genus Buhoma are two genera belonging to the snake superfamily Elapoidea. The phylogenetic placements of Psammodynastes and Buhoma within Elapoidea has been extremely unstable which has resulted in their uncertain and debated taxonomy. We used ultraconserved elements and traditional nuclear and mitochondrial markers to infer the phylogenetic relationships of these two genera with other elapoids. Psammodynastes, for which a reference genome has been sequenced, were found, with strong branch support, to be a relatively early diverging split within Elapoidea that is sister to a clade consisting of Elapidae, Micrelapidae and Lamprophiidae. Hence, we allocate Psammodynastes to its own family, Psammodynastidae new family. However, the phylogenetic position of Buhoma could not be resolved with a high degree of confidence. Attempts to identify the possible sources of conflict in the rapid radiation of elapoid snakes suggest that both hybridisation/introgression during the rapid diversification, including possible ghost introgression, as well as incomplete lineage sorting likely have had a confounding role. The usual practice of combining mitochondrial loci with nuclear genomic data appears to mislead phylogeny reconstructions in rapid radiation scenarios, especially in the absence of genome scale data.

Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing

BACKGROUND

Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates.

METHODS

In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species.

RESULTS

A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

It's a trap?! Escape from an ancient, ancestral sex chromosome system and implication of Foxl2 as the putative primary sex-determining gene in a lizard (Anguimorpha; Shinisauridae)

Although sex determination is ubiquitous in vertebrates, mechanisms of sex determination vary from environmentally to genetically influenced. In vertebrates, genetic sex determination is typically accomplished with sex chromosomes. Groups like mammals maintain conserved sex chromosome systems, while sex chromosomes in most vertebrate clades are not conserved across similar evolutionary timescales. One group inferred to have an evolutionarily stable mode of sex determination is Anguimorpha, a clade of charismatic taxa including monitor lizards, Gila monsters, and crocodile lizards. The common ancestor of extant anguimorphs possessed a ZW system that has been retained across the clade. However, the sex chromosome system in the endangered, monotypic family of crocodile lizards (Shinisauridae) has remained elusive. Here, we analyze genomic data to demonstrate that Shinisaurus has replaced the ancestral anguimorph ZW system on LG7 with a novel ZW system on LG3. The linkage group, LG3, corresponds to chromosome 9 in chicken, and this is the first documented use of this syntenic block as a sex chromosome in amniotes. Additionally, this ~1 Mb region harbors approximately 10 genes, including a duplication of the sex-determining transcription factor, Foxl2, critical for the determination and maintenance of sexual differentiation in vertebrates, and thus a putative primary sex-determining gene for Shinisaurus.

Regionalization of the vertebral column and its correlation with heart position in snakes: Implications for evolutionary pathways and morphological diversification

Spinal regionalization has important implications for the evolution of vertebrate body plans. We determined the variation in the number and morphology of vertebrae across the vertebral column (i.e., vertebral formula) for 63 snake species representing 13 families using intracolumnar variation in vertebral shape. Vertebral counts were used to determine the position of the heart, pylorus, and left kidney for each species. Across all species we observed a conspicuous midthoracic transition in vertebral shape, indicating four developmental domains of the precloacal vertebral column (cervical, anterior thoracic, posterior thoracic, and lumbar). Using phylogenetic analyses, the boundary between the anterior and posterior thoracic vertebrae was correlated with heart position. No associations were found between shifts in morphology of the vertebral column and either the pylorus or left kidney. We observed that among taxa, the number of preapex and postapex vertebrae could change independently from one another and from changes in the total number of precloacal vertebrae. Ancestral state reconstruction of the preapex and postapex vertebrae illustrated several evolutionary pathways by which diversity in the vertebral column and heart position have been attained. In addition, no conspicuous pattern was observed among the heart, pylorus, or kidney indicating that their relative positions to each other evolve independently. We conclude that snakes exhibit four morphologically distinct regions of the vertebral column. We discuss the implications of the forebody and hindbody vertebral formula on the morphological diversification of snakes.

Molecular data reveals a new genus of blindsnakes within Asiatyphlopinae from India

The genus Indotyphlops has a widespread distribution in the Indian landmass and Southeast Asia, with 20 reported species. The current classification within the genus is based on morphology. In this study, we sampled all the reported Indotyphlops species from subcontinental India, to resolve relationships within this genus and to understand biogeographic patterns that resulted in the widespread distribution. We generated sequences for five nuclear markers which were used in the global typhlopoid phylogeny and built phylogenetic trees of the superfamily Typhlopoidea. We also carried out divergence time analysis and biogeographic analysis to understand the time and modes of dispersal and diversification of these species. The results show Indotyphlops sensu lato to be polyphyletic, with the clade consisting of I. porrectus and I. exiguus sister to a clade consisting of the southeast Asian typhlopid genera Ramphotyphlops, Anilios, Malayotyphlops, Acutotyphlops, Sundatyphlops, and Indotyphlops sensu stricto. The other clade consists of I. pammeces and I. braminus from the Indian subcontinent and I. albiceps from Southeast Asia. Biogeographical analysis suggests two dispersals from Asia to the Indian landmass-an earlier dispersal from Eurasia into India led to the lineage consisting of I. porrectus and I. exiguus, followed by a later dispersal that evolved into I. pammeces and I. braminus. These results necessitate a taxonomic revision. We propose the genus Pseudoindotyphlops gen. nov. for the clade currently consisting of the most recent common ancestor (MRCA) of I. porrectus and I. exiguus, and all descendants thereof.

Ecological factors and parity mode correlate with genome size variation in squamate reptiles

BACKGROUND

Evidence of correlation between genome size, the nuclear haploid DNA content of a cell, environmental factors and life-history traits have been reported in many animal species. Genome size, however, spans over three orders of magnitude across taxa and such a correlation does not seem to follow a universal pattern. In squamate reptiles, the second most species-rich order of vertebrates, there are currently no studies investigating drivers of genome size variability. We run a series of phylogenetic generalized least-squares models on 227 species of squamates to test for possible relationships between genome size and ecological factors including latitudinal distribution, bioclimatic variables and microhabitat use. We also tested whether genome size variation can be associated with parity mode, a highly variable life history trait in this order of reptiles.

RESULTS

The best-fitting model showed that the interaction between microhabitat use and parity mode mainly accounted for genome size variation. Larger genome sizes were found in live-bearing species that live in rock/sand ecosystems and in egg-laying arboreal taxa. On the other hand, smaller genomes were found in fossorial live-bearing species.

CONCLUSIONS

Environmental factors and species parity mode appear to be among the main parameters explaining genome size variation in squamates. Our results suggest that genome size may favour adaptation of some species to certain environments or could otherwise result from the interaction between environmental factors and parity mode. Integration of genome size and genome sequencing data could help understand the role of differential genome content in the evolutionary process of genome size variation in squamates.

Genomic evidence sheds new light on phylogeny of Rhabdophis nuchalis (sensu lato) complex (Serpentes: Natricidae)

Stable taxonomy and robust phylogeny are essential for the evolution and conservation of organisms. The Rhabdophis nuchalis (sensu lato) complex presently contains three species (R. nuchalis, R. chiwen, R. pentasupralabialis). Although several studies have explored the diversity and phylogeography of this group, certain issues related to systematics and taxonomy remain unresolved. Here, based on genome-wide data, including single nucleotide polymorphisms (SNPs) generated from ddRAD-seq and mitochondrial DNA (mtDNA), we re-evaluated the phylogenetic relationships and cryptic diversity of this species group. Our results are generally consistent with previous studies but provide some new insights. Phylogenetic relationship reconstruction based on SNPs and mtDNA revealed that three species in the R. nuchalis (sensu lato) complex did not form a monophyly but each species is well supported as monophyletic lineage in SNP-based analyses. Population structure analyses showed genetic admixture between several species pairs. Additionally, the population in eastern Yunnan, China, was identified as a potential cryptic species and thus described as a new species based on morphological data. From our results and previous studies, we redefined the distribution boundary for each species in the R. nuchalis (sensu lato) species complex.

Hiding on jagged karst pinnacles: A new microendemic genus and species of a limestone-dwelling agamid lizard (Squamata: Agamidae: Draconinae) from Khammouan Province, Laos

We describe a unique new species and genus of agamid lizard from the karstic massifs of Khammouan Province, central Laos. Laodracon carsticola Gen. et sp. nov. is an elusive medium-sized lizard (maximum snout-vent length 101 mm) specifically adapted to life on limestone rocks and pinnacles. To assess the phylogenetic position of the new genus amongst other agamids, we generated DNA sequences from two mitochondrial gene fragments (16S rRNA and ND2) and three nuclear loci ( BDNF, RAG1 and c-mos), with a final alignment comprising 7 418 base pairs for 64 agamid species. Phylogenetic analyses unambiguously place the new genus in the mainland Asia subfamily Draconinae, where it forms a clade sister to the genus Diploderma from East Asia and the northern part of Southeast Asia. Morphologically, the new genus is distinguished from all other genera in Draconinae by possessing a notably swollen tail base with enlarged scales on its dorsal and ventral surfaces. Our work provides further evidence that limestone regions of Indochina represent unique "arks of biodiversity" and harbor numerous relict lineages. To date, Laodracon carsticola Gen. et sp. nov. is known from only two adult male specimens and its distribution seems to be restricted to a narrow limestone massif on the border of Khammouan and Bolikhamxai provinces of Laos. Additional studies are required to understand its life history, distribution, and conservation status.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

Large-scale snake genome analyses provide insights into vertebrate development

Snakes are a remarkable squamate lineage with unique morphological adaptations, especially those related to the evolution of vertebrate skeletons, organs, and sensory systems. To clarify the genetic underpinnings of snake phenotypes, we assembled and analyzed 14 de novo genomes from 12 snake families. We also investigated the genetic basis of the morphological characteristics of snakes using functional experiments. We identified genes, regulatory elements, and structural variations that have potentially contributed to the evolution of limb loss, an elongated body plan, asymmetrical lungs, sensory systems, and digestive adaptations in snakes. We identified some of the genes and regulatory elements that might have shaped the evolution of vision, the skeletal system and diet in blind snakes, and thermoreception in infrared-sensitive snakes. Our study provides insights into the evolution and development of snakes and vertebrates.

Genetic mapping and molecular mechanism behind color variation in the Asian vine snake

BACKGROUND

Reptiles exhibit a wide variety of skin colors, which serve essential roles in survival and reproduction. However, the molecular basis of these conspicuous colors remains unresolved.

RESULTS

We investigate color morph-enriched Asian vine snakes (Ahaetulla prasina), to explore the mechanism underpinning color variations. Transmission electron microscopy imaging and metabolomics analysis indicates that chromatophore morphology (mainly iridophores) is the main basis for differences in skin color. Additionally, we assemble a 1.77-Gb high-quality chromosome-anchored genome of the snake. Genome-wide association study and RNA sequencing reveal a conservative amino acid substitution (p.P20S) in SMARCE1, which may be involved in the regulation of chromatophore development initiated from neural crest cells. SMARCE1 knockdown in zebrafish and immunofluorescence verify the interactions among SMARCE1, iridophores, and tfec, which may determine color variations in the Asian vine snake.

CONCLUSIONS

This study reveals the genetic associations of color variation in Asian vine snakes, providing insights and important resources for a deeper understanding of the molecular and genetic mechanisms related to reptilian coloration.

A systematic review of snake translocations to identify potential tactics for reducing postrelease effects

Advancements in the field of reintroduction biology are needed, but understanding of how to effectively conduct translocations, particularly with snakes, is lacking. We conducted a systematic review of snake translocation studies to identify potential tactics for reducing postrelease effects. We included studies on intentional, human-mediated, wild-wild, or captive-wild translocations to any location, regardless of motive or number of snakes translocated. Only studies that presented results for at least 1 of 4 outcomes (movement behavior, site fidelity, survival, or population establishment) were included. We systematically searched 4 databases for published studies and used 5 methods to search the gray literature. Our search and screening criteria yielded 121 data sources, representing 130 translocation cases. We quantified the association between 15 translocation tactics and short-term translocation outcomes by calculating odds ratios and used forest plots to display results. Snake translocations involved 47 species (from mainly 2 families), and most were motivated by research, were monitored for at least 6 months, occurred in North America, and took place from the 1990s onward. The odds of a positive snake translocation outcome were highest with release of captive reared or juvenile snakes, release of social groups together, delayed release, provision of environmental enrichment or social housing before release, or minimization of distance translocated. The odds of a positive outcome were lowest when snakes were released early in their active season. Our results do not demonstrate causation, but outcomes of snake translocation were associated with 8 tactics (4 of which were strongly correlated). In addition to targeted comparative studies, we recommend practitioners consider the possible influence of these tactics when planning snake translocations.

What snakes and caecilians have in common? Molecular interaction units and the independent origins of similar morphotypes in Tetrapoda

Developmental pathways encompass transcription factors and cis-regulatory elements that interact as transcription factor-regulatory element (TF-RE) units. Independent origins of similar phenotypes likely involve changes in different parts of these units, a hypothesis promisingly tested addressing the evolution of the rib-associated lumbar (RAL) morphotype that characterizes emblematic animals such as snakes and elephants. Previous investigation in these lineages identified a polymorphism in the Homology region 1 [H1] enhancer of the Myogenic factor-5 [Myf5], which interacts with HOX10 proteins to modulate rib development. Here we address the evolution of TF-RE units focusing on independent origins of RAL morphotypes. We compiled an extensive database for H1-Myf5 and HOX10 sequences with two goals: (i) evaluate if the enhancer polymorphism is present in amphibians exhibiting the RAL morphotype and (ii) test a hypothesis of enhanced evolutionary flexibility mediated by TF-RE units, according to which independent origins of the RAL morphotype might involve changes in either component of the interaction unit. We identified the H1-Myf5 polymorphism in lineages that diverged around 340 Ma, including Lissamphibia. Independent origins of the RAL morphotype in Tetrapoda involved sequence variation in either component of the TF-RE unit, confirming that different changes may similarly affect the phenotypic outcome of a given developmental pathway.

No link between population isolation and speciation rate in squamate reptiles

Rates of species formation vary widely across the tree of life and contribute to massive disparities in species richness among clades. This variation can emerge from differences in metapopulation-level processes that affect the rates at which lineages diverge, persist, and evolve reproductive barriers and ecological differentiation. For example, populations that evolve reproductive barriers quickly should form new species at faster rates than populations that acquire reproductive barriers more slowly. This expectation implicitly links microevolutionary processes (the evolution of populations) and macroevolutionary patterns (the profound disparity in speciation rate across taxa). Here, leveraging extensive field sampling from the Neotropical Cerrado biome in a biogeographically controlled natural experiment, we test the role of an important microevolutionary process-the propensity for population isolation-as a control on speciation rate in lizards and snakes. By quantifying population genomic structure across a set of codistributed taxa with extensive and phylogenetically independent variation in speciation rate, we show that broad-scale patterns of species formation are decoupled from demographic and genetic processes that promote the formation of population isolates. Population isolation is likely a critical stage of speciation for many taxa, but our results suggest that interspecific variability in the propensity for isolation has little influence on speciation rates. These results suggest that other stages of speciation-including the rate at which reproductive barriers evolve and the extent to which newly formed populations persist-are likely to play a larger role than population isolation in controlling speciation rate variation in squamates.

Phylogenomic analyses resolve relationships among garter snakes (Thamnophis: Natricinae: Colubridae) and elucidate biogeographic history and morphological evolution

Garter snakes (Thamnophis) are a successful group of natricines endemic to North America. They have become important natural models for ecological and evolutionary research, yet prior efforts to resolve phylogenetic relationships have resulted in conflicting topologies and weak support for certain relationships. Here, we use genomic data generated with a reduced representation double-digest RADseq approach to reassess evolutionary relationships across Thamnophis. We then use the resulting phylogeny to better understand how biogeography and feeding ecology have influenced lineage diversification and morphological evolution. We recovered highly congruent and strongly supported topologies from maximum likelihood and Bayesian analyses, but some discordance with a multispecies coalescent approach. All phylogenomic estimates split Thamnophis into two clades largely defined by northern and southern North American species. Divergence time estimates and biogeographic analyses indicate a mid-Miocene origin of Thamnophis in Mexico. In addition, historic vicariant events thought to explain biogeographic patterns in other lineages (e.g., Isthmus of Tehuantepec, Rocky Mountain Range, and Trans-Mexican Volcanic Belt) appear to have influenced patterns of diversification in Thamnophis as well. Analyses of morphological traits associated with feeding ecology showed moderate to strong phylogenetic signal. Nevertheless, phylogenetic ANOVA suggested significant differences in certain cranial morphologies between aquatic specialists and garter snakes that are terrestrial-aquatic generalists, independent of evolutionary history. Our new estimate of Thamnophis phylogeny yields an improved understanding of the biogeographic history and morphological evolution of garter snakes, and provides a robust framework for future research on these snakes.

Molecular phylogeny of Neotropical Parrot Snakes (Serpentes: Colubrinae: Leptophis) supports underestimated species richness

Tetrapod taxa with broad geographic distributions across the Neotropics are often composed of multiple evolutionary lineages. In this paper, we present the most complete phylogeny of Leptophis to date and assess morphology-based species limits within the broadly distributed green parrot snake Leptophis ahaetulla sensu lato, which occurs from Mexico to Argentina. Although L. ahaetulla sensu stricto, L. nigromarginatus and L. occidentalis were recovered as paraphyletic, tree topology tests failed to reject their monophyly. Monophyly of L. bocourti, L. coeruleodorsus, L. cupreus, L. depressirostris, L. marginatus, L. riveti and L. sp. nov. was strongly supported. Our phylogenetic trees support recognition of multiple species within Leptophis ahaetulla sensu lato and suggest that color evolution and the uplift of the Andes played an important role in the diversification of parrot snakes.

How mitonuclear discordance and geographic variation have confounded species boundaries in a widely studied snake

As DNA sequencing technologies and methods for delimiting species with genomic data become more accessible and numerous, researchers have more tools than ever to investigate questions in systematics and phylogeography. However, easy access to sophisticated computational tools is not without its drawbacks. Choosing the right approach for one's question can be challenging when presented with multitudinous options, some of which fail to distinguish between species and intraspecific population structure. Here, we employ a methodology that emphasizes intensive geographic sampling, particularly at contact zones between populations, with a focus on differentiating intraspecific genetic clusters from species in the Pantherophis guttatus complex, a group of North American ratsnakes. Using a mitochondrial marker as well as ddRADseq data, we find evidence of mitonuclear discordance which has contributed to historical confusion about the relationships within this group. Additionally, we identify geographically and genetically structured populations within the species Pantherophis emoryi that are congruent with previously described morphological variation. Importantly, we find that these structured populations within P. emoryi are highly admixed throughout the range of the species and show no evidence of any reproductive isolation. Our data support a revision of the taxonomy of this group, and we recognize two species within the complex and three subspecies within P. emoryi. This study illustrates the importance of thorough sampling of contact zones and consideration of gene flow when delimiting species in widespread complexes containing parapatric lineages.

ToxCodAn: a new toxin annotator and guide to venom gland transcriptomics

MOTIVATION

Next-generation sequencing has become exceedingly common and has transformed our ability to explore nonmodel systems. In particular, transcriptomics has facilitated the study of venom and evolution of toxins in venomous lineages; however, many challenges remain. Primarily, annotation of toxins in the transcriptome is a laborious and time-consuming task. Current annotation software often fails to predict the correct coding sequence and overestimates the number of toxins present in the transcriptome. Here, we present ToxCodAn, a python script designed to perform precise annotation of snake venom gland transcriptomes. We test ToxCodAn with a set of previously curated transcriptomes and compare the results to other annotators. In addition, we provide a guide for venom gland transcriptomics to facilitate future research and use Bothrops alternatus as a case study for ToxCodAn and our guide.

RESULTS

Our analysis reveals that ToxCodAn provides precise annotation of toxins present in the transcriptome of venom glands of snakes. Comparison with other annotators demonstrates that ToxCodAn has better performance with regard to run time ($>20x$ faster), coding sequence prediction ($>3x$ more accurate) and the number of toxins predicted (generating $>4x$ less false positives). In this sense, ToxCodAn is a valuable resource for toxin annotation. The ToxCodAn framework can be expanded in the future to work with other venomous lineages and detect novel toxins.

Replacement and Parallel Simplification of Nonhomologous Proteinases Maintain Venom Phenotypes in Rear-Fanged Snakes

Novel phenotypes are commonly associated with gene duplications and neofunctionalization, less documented are the cases of phenotypic maintenance through the recruitment of novel genes. Proteolysis is the primary toxic character of many snake venoms, and ADAM metalloproteinases, named snake venom metalloproteinases (SVMPs), are largely recognized as the major effectors of this phenotype. However, by investigating original transcriptomes from 58 species of advanced snakes (Caenophidia) across their phylogeny, we discovered that a different enzyme, matrix metalloproteinase (MMP), is actually the dominant venom component in three tribes (Tachymenini, Xenodontini, and Conophiini) of rear-fanged snakes (Dipsadidae). Proteomic and functional analyses of these venoms further indicate that MMPs are likely playing an "SVMP-like" function in the proteolytic phenotype. A detailed look into the venom-specific sequences revealed a new highly expressed MMP subtype, named snake venom MMP (svMMP), which originated independently on at least three occasions from an endogenous MMP-9. We further show that by losing ancillary noncatalytic domains present in its ancestors, svMMPs followed an evolutionary path toward a simplified structure during their expansion in the genomes, thus paralleling what has been proposed for the evolution of their Viperidae counterparts, the SVMPs. Moreover, we inferred an inverse relationship between the expression of svMMPs and SVMPs along the evolutionary history of Xenodontinae, pointing out that one type of enzyme may be substituting for the other, whereas the general (metallo)proteolytic phenotype is maintained. These results provide rare evidence on how relevant phenotypic traits can be optimized via natural selection on nonhomologous genes, yielding alternate biochemical components.

Snake W Sex Chromosome: The Shadow of Ancestral Amniote Super-Sex Chromosome

: Heteromorphic sex chromosomes, particularly the ZZ/ZW sex chromosome system of birds and some reptiles, undergo evolutionary dynamics distinct from those of autosomes. The W sex chromosome is a unique karyological member of this heteromorphic pair, which has been extensively studied in snakes to explore the origin, evolution, and genetic diversity of amniote sex chromosomes. The snake W sex chromosome offers a fascinating model system to elucidate ancestral trajectories that have resulted in genetic divergence of amniote sex chromosomes. Although the principal mechanism driving evolution of the amniote sex chromosome remains obscure, an emerging hypothesis, supported by studies of W sex chromosomes of squamate reptiles and snakes, suggests that sex chromosomes share varied genomic blocks across several amniote lineages. This implies the possible split of an ancestral super-sex chromosome via chromosomal rearrangements. We review the major findings pertaining to sex chromosomal profiles in amniotes and discuss the evolution of an ancestral super-sex chromosome by collating recent evidence sourced mainly from the snake W sex chromosome analysis. We highlight the role of repeat-mediated sex chromosome conformation and present a genomic landscape of snake Z and W chromosomes, which reveals the relative abundance of major repeats, and identifies the expansion of certain transposable elements. The latest revolution in chromosomics, i.e., complete telomere-to-telomere assembly, offers mechanistic insights into the evolutionary origin of sex chromosomes.

Genome mapping of a LYST mutation in corn snakes indicates that vertebrate chromatophore vesicles are lysosome-related organelles

Reptiles exhibit a spectacular diversity of skin colors and patterns brought about by the interactions among three chromatophore types: black melanophores with melanin-packed melanosomes, red and yellow xanthophores with pteridine- and/or carotenoid-containing vesicles, and iridophores filled with light-reflecting platelets generating structural colors. Whereas the melanosome, the only color-producing endosome in mammals and birds, has been documented as a lysosome-related organelle, the maturation paths of xanthosomes and iridosomes are unknown. Here, we first use 10x Genomics linked-reads and optical mapping to assemble and annotate a nearly chromosome-quality genome of the corn snake Pantherophis guttatus The assembly is 1.71 Gb long, with an N50 of 16.8 Mb and L50 of 24. Second, we perform mapping-by-sequencing analyses and identify a 3.9-Mb genomic interval where the lavender variant resides. The lavender color morph in corn snakes is characterized by gray, rather than red, blotches on a pink, instead of orange, background. Third, our sequencing analyses reveal a single nucleotide polymorphism introducing a premature stop codon in the lysosomal trafficking regulator gene (LYST) that shortens the corresponding protein by 603 amino acids and removes evolutionary-conserved domains. Fourth, we use light and transmission electron microscopy comparative analyses of wild type versus lavender corn snakes and show that the color-producing endosomes of all chromatophores are substantially affected in the LYST mutant. Our work provides evidence characterizing xanthosomes in xanthophores and iridosomes in iridophores as lysosome-related organelles.

Landscape Genomics of a Widely Distributed Snake, Dolichophis caspius (Gmelin, 1789) across Eastern Europe and Western Asia

Across the distribution of the Caspian whipsnake (Dolichophis caspius), populations have become increasingly disconnected due to habitat alteration. To understand population dynamics and this widespread but locally endangered snake’s adaptive potential, we investigated population structure, admixture, and effective migration patterns. We took a landscape-genomic approach to identify selected genotypes associated with environmental variables relevant to D. caspius. With double-digest restriction-site associated DNA (ddRAD) sequencing of 53 samples resulting in 17,518 single nucleotide polymorphisms (SNPs), we identified 8 clusters within D. caspius reflecting complex evolutionary patterns of the species. Estimated Effective Migration Surfaces (EEMS) revealed higher-than-average gene flow in most of the Balkan Peninsula and lower-than-average gene flow along the middle section of the Danube River. Landscape genomic analysis identified 751 selected genotypes correlated with 7 climatic variables. Isothermality correlated with the highest number of selected genotypes (478) located in 41 genes, followed by annual range (127) and annual mean temperature (87). We conclude that environmental variables, especially the day-to-night temperature oscillation in comparison to the summer-to-winter oscillation, may have an important role in the distribution and adaptation of D. caspius.

Revisiting the Karyotype Evolution of Neotropical Boid Snakes: A Puzzle Mediated by Chromosomal Fissions

The Boidae family is an ancient group of snakes widely distributed across the Neotropical region, where several biogeographic events contributed towards shaping their evolution and diversification. Most species of this family have a diploid number composed of 2n = 36; however, among Booidea families, the Boidae stands out by presenting the greatest chromosomal diversity, with 2n ranging between 36 and 44 chromosomes and an undifferentiated XY sex chromosome system. Here, we applied a comparative chromosome analysis using cross-species chromosome paintings in five species representing four Boidae genera, to decipher the evolutionary dynamics of some chromosomes in these Neotropical snakes. Our study included all diploid numbers (2n = 36, 40, and 44) known for this family and our comparative chromosomal mappings point to a strong evolutionary relationship among the genera Boa, Corallus, Eunectes, and Epicrates. The results also allowed us to propose the cytogenomic diversification that had occurred in this family: a process mediated by centric fissions, including fission events of the putative and undifferentiated XY sex chromosome system in the 2n = 44 karyotype, which is critical in solving the puzzle of the karyotype evolution of boid snakes.

The tuatara genome reveals ancient features of amniote evolution

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Interstitial Telomeric Repeats Are Rare in Turtles

Telomeres are nucleoprotein complexes protecting chromosome ends in most eukaryotic organisms. In addition to chromosome ends, telomeric-like motifs can be accumulated in centromeric, pericentromeric and intermediate (i.e., between centromeres and telomeres) positions as so-called interstitial telomeric repeats (ITRs). We mapped the distribution of (TTAGGG)_n repeats in the karyotypes of 30 species from nine families of turtles using fluorescence in situ hybridization. All examined species showed the expected terminal topology of telomeric motifs at the edges of chromosomes. We detected ITRs in only five species from three families. Combining our and literature data, we inferred seven independent origins of ITRs among turtles. ITRs occurred in turtles in centromeric positions, often in several chromosomal pairs, in a given species. Their distribution does not correspond directly to interchromosomal rearrangements. Our findings support that centromeres and non-recombining parts of sex chromosomes are very dynamic genomic regions, even in turtles, a group generally thought to be slowly evolving. However, in contrast to squamate reptiles (lizards and snakes), where ITRs were found in more than half of the examined species, and birds, the presence of ITRs is generally rare in turtles, which agrees with the expected low rates of chromosomal rearrangements and rather slow karyotype evolution in this group.

Combining genetic and demographic monitoring better informs conservation of an endangered urban snake

Conversion and fragmentation of wildlife habitat often leads to smaller and isolated populations and can reduce a species' ability to disperse across the landscape. As a consequence, genetic drift can quickly lower genetic variation and increase vulnerability to extirpation. For species of conservation concern, quantification of population size and connectivity can clarify the influence of genetic drift in local populations and provides important information for conservation management and recovery strategies. Here, we used genome-wide single nucleotide polymorphism (SNP) data and capture-mark-recapture methods to evaluate the genetic diversity and demography within seven focal sites of the endangered San Francisco gartersnake (Thamnophis sirtalis tetrataenia), a species affected by alteration and isolation of wetland habitats throughout its distribution. The primary goals were to determine the population structure and degree of genetic isolation among T. s. tetrataenia populations and estimate effective size and population abundance within sites to better understand the present and future importance of genetic drift. We also used temporally sampled datasets to examine the magnitude of genetic change over time. We found moderate population genetic structure throughout the San Francisco Peninsula that partitions sites into northern and southern regional clusters. Point estimates of both effective size and population abundance were generally small (≤ 100) for a majority of the sites, and estimates were particularly low in the northern populations. Genetic analyses of temporal datasets indicated an increase in genetic differentiation, especially for the most geographically isolated sites, and decreased genetic diversity over time in at least one site (Pacifica). Our results suggest that drift-mediated processes as a function of small population size and reduced connectivity from neighboring populations may decrease diversity and increase differentiation. Improving genetic diversity and connectivity among T. s. tetrataenia populations could promote persistence of this endangered snake.

Prevalence of inclusion body disease and associated comorbidity in captive collections of boid and pythonid snakes in Belgium

Inclusion body disease (IBD) is caused by reptarenaviruses and constitutes one of the most notorious viral diseases in snakes. Although central nervous system disease and various other clinical signs have been attributed to IBD in boid and pythonid snakes, studies that unambiguously reveal the clinical course of natural IBD and reptarenavirus infection are scarce. In the present study, the prevalence of IBD and reptarenaviruses in captive snake collections and the correlation of IBD and reptarenavirus infection with the clinical status of the sampled snakes were investigated. In three IBD positive collections, long-term follow-up during a three- to seven-year period was performed. A total of 292 snakes (178 boas and 114 pythons) from 40 collections in Belgium were sampled. In each snake, blood and buffy coat smears were evaluated for the presence of IBD inclusion bodies (IB) and whole blood was tested for reptarenavirus RNA by RT-PCR. Of all tested snakes, 16.5% (48/292) were positive for IBD of which all were boa constrictors (34.0%; 48/141) and 17.1% (50/292) were reptarenavirus RT-PCR positive. The presence of IB could not be demonstrated in any of the tested pythons, while 5.3% (6/114) were reptarenavirus positive. In contrast to pythons, the presence of IB in peripheral blood cells in boa constrictors is strongly correlated with reptarenavirus detection by RT-PCR (P<0.0001). Although boa constrictors often show persistent subclinical infection, long-term follow-up indicated that a considerable number (22.2%; 6/27) of IBD/reptarenavirus positive boas eventually develop IBD associated comorbidities.

Discovery of a deeply divergent new lineage of vine snake (Colubridae: Ahaetuliinae: Proahaetulla gen. nov.) from the southern Western Ghats of Peninsular India with a revised key for Ahaetuliinae

The Western Ghats are well known as a biodiversity hotspot, but the full extent of its snake diversity is yet to be uncovered. Here, we describe a new genus and species of vine snake Proahaetulla antiqua gen. et sp. nov., from the Agasthyamalai hills in the southern Western Ghats. It was found to be a member of the Ahaetuliinae clade, which currently comprises the arboreal snake genera Ahaetulla, Dryophiops, Dendrelaphis and Chrysopelea, distributed in South and Southeast Asia. Proahaetulla shows a sister relationship with all currently known taxa belonging to the genus Ahaetulla, and shares ancestry with Dryophiops. In addition to its phylogenetic position and significant genetic divergence, this new taxon is also different in morphology from members of Ahaetuliinae in a combination of characters, having 12-13 partially serrated keels on the dorsal scale rows, 20 maxillary teeth and 3 postocular scales. Divergence dating reveals that the new genus is ancient, dating back to the Mid-Oligocene, and is one of the oldest persisting monotypic lineages of snakes in the Western Ghats. This discovery adds to the growing list of ancient lineages endemic to the Agasthyamalai hills and underscores the biogeographic significance of this isolated massif in the southern Western Ghats.

Taxonomic and conservation implications of population genetic admixture, mito-nuclear discordance, and male-biased dispersal of a large endangered snake, Drymarchon couperi

Accurate species delimitation and description are necessary to guide effective conservation of imperiled species, and this synergy is maximized when multiple data sources are used to delimit species. We illustrate this point by examining Drymarchon couperi (Eastern Indigo Snake), a large, federally-protected species in North America that was recently divided into two species based on gene sequence data from three loci and heuristic morphological assessment. Here, we re-evaluate the two-species hypothesis for D. couperi by evaluating both population genetic and gene sequence data. Our analyses of 14 microsatellite markers revealed 6–8 genetic population clusters with significant admixture, particularly across the contact zone between the two hypothesized species. Phylogenetic analyses of gene sequence data with maximum-likelihood methods suggested discordance between mitochondrial and nuclear markers and provided phylogenetic support for one species rather than two. For these reasons, we place Drymarchon kolpobasileus into synonymy with D. couperi. We suggest inconsistent patterns between mitochondrial and nuclear DNA are driven by high dispersal of males relative to females. We advocate for species delimitation exercises that evaluate admixture and gene flow in addition to phylogenetic analyses, particularly when the latter reveal monophyletic lineages. This is particularly important for taxa, such as squamates, that exhibit strong sex-biased dispersal. Problems associated with over-delimitation of species richness can become particularly acute for threatened and endangered species, because of high costs to conservation when taxonomy demands protection of more individual species than are supported by accumulating data.

Karyotype Analysis of Four Blind Snake Species (Reptilia: Squamata: Scolecophidia) and Karyotypic Changes in Serpentes

The suborder Serpentes is divided into 2 infraorders, Scolecophidia and Alethinophidia, which diverged at an early stage of snake diversification. In this study, we examined karyotypes of 4 scolecophidian species (Letheobia simonii, Xerotyphlops vermicularis, Indotyphlops braminus, and Myriopholis macrorhyncha) and performed FISH with 18S-28S rDNA as well as microchromosomal and Z chromosome-linked genes of Elaphe quadrivirgata (Alethinophidia) to investigate the karyotype evolution in the scolecophidian lineage. Diploid chromosome numbers of X. vermicularis and L. simonii were 30 (16 macrochromosomes and 14 microchromosomes) and 32 (16 macrochromosomes and 16 microchromosomes), respectively. The karyotype of a female M. macrorhyncha consisted of 15 macrochromosomes and 19 microchromosomes, including a heterochromatic microchromosome, indicating the presence of a heteromorphic chromosome pair. E. quadrivirgata Z-linked genes mapped to chromosome 4 of M. macrorhyncha, not to the heteromorphic pair. Therefore, M. macrorhyncha may have differentiated ZW sex chromosomes which are not homologous to those of E. quadrivirgata. One of the E. quadrivirgata microchromosomal genes mapped to the terminal region of chromosome 4q in X. vermicularis, suggesting that fusions between microchromosomes and macrochromosomes occurred in this species. rDNA was localized in different macrochromosomal pairs in the 2 diploid scolecophidian snakes examined here, whereas the gene location in a microchromosomal pair was conserved in 5 alethinophidian species examined. These results might imply the occurrence of chromosome fusions in the scolecophidian lineages. In I. braminus, a unique parthenogenetic snake with a triploid karyotype (21 macrochromosomes and 21 microchromosomes), morphological heteromorphisms were identified in chromosomes 1 and 7. Such heteromorphisms in 2 chromosomes were also observed in individuals from distant locations in the broad distribution range of this species, suggesting that the heteromorphisms were fixed in the genome at an early stage of its speciation.

Phylogenetic origins for severe acetaminophen toxicity in snake species compared to other vertebrate taxa

While it has been known for a while that some snake species are extremely sensitive to acetaminophen, the underlying mechanism for this toxicity has not been reported. To investigate if essential detoxification enzymes are missing in snake species that are responsible for biotransformation of acetaminophen in other vertebrate species, livers were collected from a variety of snake species, together with samples from alligator, snapping turtle, cat, rat, and cattle. Subcellular fractions were analyzed for enzymatic activities of phenol-type sulfotransferase and UDP‑glucuronosyltransferase, total glutathione S‑transferase, and N‑acetyltransferase. The results showed that none of the snake species, together with the cat samples, had any phenol-type glucuronidation activity, and that this activity was much lower in alligator and turtle samples than in the mammalian species. Combined with the lack of N‑acetyltransferase activity in snakes and cats, this would explain the accumulation of the aminophenol metabolite, which induces methemoglobinemia and subsequent suffocation of snakes and cats after acetaminophen exposure. While previous investigations have concluded that in cats the gene for the phenol-type glucuronosyltransferase isoform has turned into a pseudogene because of several point mutations, evaluation of genomic information for snake species revealed that they have only 2 genes that may code for glucuronosyltransferase isoforms. Similarity of these genes with mammalian genes is <50%, and suggests that the expressed enzymes may act on other types of substrates than aromatic amines. This indicates that the extreme sensitivity for acetaminophen in snakes is based on a different phylogenetic origin than the sensitivity observed in cats.

Noise and biases in genomic data may underlie radically different hypotheses for the position of Iguania within Squamata

Squamate reptiles are a major component of vertebrate biodiversity whose crown-clade traces its origin to a narrow window of time in the Mesozoic during which the main subclades diverged in rapid succession. Deciphering phylogenetic relationships among these lineages has proven challenging given the conflicting signals provided by genomic and phenomic data. Most notably, the placement of Iguania has routinely differed between data sources, with morphological evidence supporting a sister relationship to the remaining squamates (Scleroglossa hypothesis) and molecular data favoring a highly nested position alongside snakes and anguimorphs (Toxicofera hypothesis). We provide novel insights by generating an expanded morphological dataset and exploring the presence of phylogenetic signal, noise, and biases in molecular data. Our analyses confirm the presence of strong conflicting signals for the position of Iguania between morphological and molecular datasets. However, we also find that molecular data behave highly erratically when inferring the deepest branches of the squamate tree, a consequence of limited phylogenetic signal to resolve this ancient radiation with confidence. This, in turn, seems to result from a rate of evolution that is too high for historical signals to survive to the present. Finally, we detect significant systematic biases, with iguanians and snakes sharing faster rates of molecular evolution and a similarly biased nucleotide composition. A combination of scant phylogenetic signal, high levels of noise, and the presence of systematic biases could result in the misplacement of Iguania. We regard this explanation to be at least as plausible as the complex scenario of convergence and reversals required for morphological data to be misleading. We further evaluate and discuss the utility of morphological data to resolve ancient radiations, as well as its impact in combined-evidence phylogenomic analyses, with results relevant for the assessment of evidence and conflict across the Tree of Life.

Phylogeography of the reticulated python (Malayopython reticulatus ssp.): Conservation implications for the worlds' most traded snake species

As an important economic natural resource in Southeast Asia, reticulated pythons (Malayopython reticulatus ssp.) are primarily harvested from the wild for their skins-which are prized in the luxury leather goods industry. Trade dynamics of this CITES Appendix II listed species are complex and management approaches on the country or regional level appear obscure. Little is known about the actual geographic point-of-harvest of snakes, how genetic diversity is partitioned across the species range, how current harvest levels may affect the genetic viability of populations, and whether genetic structure could (or should) be accounted for when managing harvest quotas. As an initial survey, we use mitochondrial sequence data to define the broad-scale geographic structure of genetic diversity across a significant portion of the reticulated python's native range. Preliminary results reveal: (1) prominent phylogenetic structure across populations east and west of Huxley's modification of Wallace's line. Thirty-four haplotypes were apportioned across two geographically distinct groups, estimated to be moderately (5.2%); (2) Philippine, Bornean and Sulawesian populations appear to cluster distinctly; (3) individuals from Ambon Island suggest recent human introduction. Malayopython reticulatus is currently managed as a single taxonomic unit across Southeast Asia yet these initial results may justify special management considerations of the Philippine populations as a phylogenetically distinct unit, that warrants further examination. In Indonesia, genetic structure does not conform tightly to political boundaries and therefore we advocate the precautionary designation and use of Evolutionary Significant Units within Malayopython reticulatus, to inform and guide regional adaptive management plans.

A Step-by-Step Guide to Assemble a Reptilian Genome

Multiple technologies and software are now available facilitating the de novo sequencing and assembly of any vertebrate genome. Yet the quality of most available sequenced genomes is substantially poorer than that of the golden standard in the field: the human genome. Here, we present a step-by-step protocol for the successful sequencing and assembly of a high-quality snake genome that can be applied to any other reptilian or avian species. We combine the great sequencing depth and accuracy of short reads with the use of different insert size libraries for extended scaffolding followed by optical mapping. We show that this procedure improved the corn snake scaffold N50 from 3.7 kbp to 1.4 Mbp, currently making it one of the snake genomes with the longest scaffolds. Short guidelines are also given on the extraction of long DNA molecules from reptilian blood and the necessary modifications in DNA extraction protocols. This chapter is accompanied by a website ( www.reptilomics.org/stepbystep.html ), where we provide links to the suggested software, examples of input and output files, and running parameters.

A Tricky Trait: Applying the Fruits of the "Function Debate" in the Philosophy of Biology to the "Venom Debate" in the Science of Toxinology

The "function debate" in the philosophy of biology and the "venom debate" in the science of toxinology are conceptually related. Venom systems are complex multifunctional traits that have evolved independently numerous times throughout the animal kingdom. No single concept of function, amongst those popularly defended, appears adequate to describe these systems in all their evolutionary contexts and extant variations. As such, a pluralistic view of function, previously defended by some philosophers of biology, is most appropriate. Venom systems, like many other functional traits, exist in nature as points on a continuum and the boundaries between "venomous" and "non-venomous" species may not always be clearly defined. This paper includes a brief overview of the concept of function, followed by in-depth discussion of its application to venom systems. A sound understanding of function may aid in moving the venom debate forward. Similarly, consideration of a complex functional trait such as venom may be of interest to philosophers of biology.

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.

Coral snakes predict the evolution of mimicry across New World snakes

Batesian mimicry, in which harmless species (mimics) deter predators by deceitfully imitating the warning signals of noxious species (models), generates striking cases of phenotypic convergence that are classic examples of evolution by natural selection. However, mimicry of venomous coral snakes has remained controversial because of unresolved conflict between the predictions of mimicry theory and empirical patterns in the distribution and abundance of snakes. Here we integrate distributional, phenotypic and phylogenetic data across all New World snake species to demonstrate that shifts to mimetic coloration in nonvenomous snakes are highly correlated with coral snakes in both space and time, providing overwhelming support for Batesian mimicry. We also find that bidirectional transitions between mimetic and cryptic coloration are unexpectedly frequent over both long- and short-time scales, challenging traditional views of mimicry as a stable evolutionary 'end point' and suggesting that insect and snake mimicry may have different evolutionary dynamics.

Pliocene-Pleistocene lineage diversifications in the Eastern Indigo Snake (Drymarchon couperi) in the Southeastern United States

Indigo Snakes (Drymarchon; with five currently recognized species) occur from northern Argentina, northward to the United States in southern Texas and eastward in disjunct populations in Florida and Georgia. Based on this known allopatry and a difference in supralabial morphology the two United States taxa previously considered as subspecies within D. corais (Boie 1827), the Western Indigo Snake, D. melanurus erebennus (Cope 1860), and Eastern Indigo Snake, D. couperi (Holbrook 1842), are currently recognized as separate species. Drymarchon couperi is a Federally-designated Threatened species by the United States Fish and Wildlife Service under the Endangered Species Act, and currently being incorporated into a translocation program. This, combined with its disjunct distribution makes it a prime candidate for studying speciation and genetic divergence. In this study, we (1) test the hypothesis that D. m. erebennus and D. couperi are distinct lineages by analyzing 2411 base pairs (bp) of two mitochondrial (mtDNA) loci and one single copy nuclear (scnDNA) locus; (2) estimate the timing of speciation using a relaxed phylogenetics method to determine if Milankovitch cycles during the Pleistocene might have had an influence on lineage diversifications; (3) examine historical population demography to determine if identified lineages have undergone population declines, expansions, or remained stable during the most recent Milankovitch cycles; and (4) use this information to assist in an effective and scientifically sound translocation program. Our molecular data support the initial hypothesis that D. melanurus and D. couperi should be recognized as distinct species, but further illustrate that D. couperi is split into two distinct genetic lineages that correspond to historical biogeography and sea level changes in peninsular Florida. These two well-supported genetic lineages (herein termed Atlantic and Gulf lineages) illustrate a common biogeographic distributional break previously identified for other plants and animals, suggesting that these organisms might have shared a common evolutionary history related to historic sea level changes caused by Milankovitch cycles. Our estimated divergence times suggest that the most recent common ancestor (MRCA) between D. melanurus and southeastern United States Drymarchon occurred ca. 5.9Ma (95% HPD=2.5-9.8Ma; during the late Blancan of the Pleistocene through the Hemphillian of the Miocene), whereas the MRCA between the Atlantic and Gulf lineages in the southeastern United States occurred ca. 2.0Ma (95% HPD=0.7-3.7Ma; during the Irvingtonian of the Pleistocene through the Blancan of the Pliocene). During one or more glacial intervals within these times, these two lineages must have become separated and evolved independently. Despite numerous Milankovitch cycles along with associated forming of physical barriers (i.e., sea level fluctuations, high elevation sand ridges, clayey soils, and/or insufficient habitats) since their initial lineage diversification, these two lineages have likely come in and out of contact with each other many times, yet today they still illustrate near discrete geographic distributions. Although the Atlantic and Gulf lineages appear to be cryptic, a thorough study examining morphological characters should be conducted. We believe that our molecular data is crucial and should be incorporated in making conscious decisions in the management of a translocation program. We suggest that source populations for translocations include maintaining the integrity of the known genetic lineages found herein, as well as those coming from the closest areas that currently support sizable Drymarchon populations.

Comparing species tree estimation with large anchored phylogenomic and small Sanger-sequenced molecular datasets: an empirical study on Malagasy pseudoxyrhophiine snakes

BACKGROUND

Using molecular data generated by high throughput next generation sequencing (NGS) platforms to infer phylogeny is becoming common as costs go down and the ability to capture loci from across the genome goes up. While there is a general consensus that greater numbers of independent loci should result in more robust phylogenetic estimates, few studies have compared phylogenies resulting from smaller datasets for commonly used genetic markers with the large datasets captured using NGS. Here, we determine how a 5-locus Sanger dataset compares with a 377-locus anchored genomics dataset for understanding the evolutionary history of the pseudoxyrhophiine snake radiation centered in Madagascar. The Pseudoxyrhophiinae comprise ~86 % of Madagascar's serpent diversity, yet they are poorly known with respect to ecology, behavior, and systematics. Using the 377-locus NGS dataset and the summary statistics species-tree methods STAR and MP-EST, we estimated a well-supported species tree that provides new insights concerning intergeneric relationships for the pseudoxyrhophiines. We also compared how these and other methods performed with respect to estimating tree topology using datasets with varying numbers of loci.

METHODS

Using Sanger sequencing and an anchored phylogenomics approach, we sequenced datasets comprised of 5 and 377 loci, respectively, for 23 pseudoxyrhophiine taxa. For each dataset, we estimated phylogenies using both gene-tree (concatenation) and species-tree (STAR, MP-EST) approaches. We determined the similarity of resulting tree topologies from the different datasets using Robinson-Foulds distances. In addition, we examined how subsets of these data performed compared to the complete Sanger and anchored datasets for phylogenetic accuracy using the same tree inference methodologies, as well as the program *BEAST to determine if a full coalescent model for species tree estimation could generate robust results with fewer loci compared to the summary statistics species tree approaches. We also examined the individual gene trees in comparison to the 377-locus species tree using the program MetaTree.

RESULTS

Using the full anchored dataset under a variety of methods gave us the same, well-supported phylogeny for pseudoxyrhophiines. The African pseudoxyrhophiine Duberria is the sister taxon to the Malagasy pseudoxyrhophiines genera, providing evidence for a monophyletic radiation in Madagascar. In addition, within Madagascar, the two major clades inferred correspond largely to the aglyphous and opisthoglyphous genera, suggesting that feeding specializations associated with tooth venom delivery may have played a major role in the early diversification of this radiation. The comparison of tree topologies from the concatenated and species-tree methods using different datasets indicated the 5-locus dataset cannot beused to infer a correct phylogeny for the pseudoxyrhophiines under any method tested here and that summary statistics methods require 50 or more loci to consistently recover the species-tree inferred using the complete anchored dataset. However, as few as 15 loci may infer the correct topology when using the full coalescent species tree method *BEAST. MetaTree analyses of each gene tree from the Sanger and anchored datasets found that none of the individual gene trees matched the 377-locus species tree, and that no gene trees were identical with respect to topology.

CONCLUSIONS

Our results suggest that ≥50 loci may be necessary to confidently infer phylogenies when using summaryspecies-tree methods, but that the coalescent-based method *BEAST consistently recovers the same topology using only 15 loci. These results reinforce that datasets with small numbers of markers may result in misleading topologies, and further, that the method of inference used to generate a phylogeny also has a major influence on the number of loci necessary to infer robust species trees.

Island colonisation and the evolutionary rates of body size in insular neonate snakes

Island colonisation by animal populations is often associated with dramatic shifts in body size. However, little is known about the rates at which these evolutionary shifts occur, under what precise selective pressures and the putative role played by adaptive plasticity on driving such changes. Isolation time played a significant role in the evolution of body size in island Tiger snake populations, where adaptive phenotypic plasticity followed by genetic assimilation fine-tuned neonate body and head size (hence swallowing performance) to prey size. Here I show that in long isolated islands (>6000 years old) and mainland populations, neonate body mass and snout-vent length are tightly correlated with the average prey body mass available at each site. Regression line equations were used to calculate body size values to match prey size in four recently isolated populations of Tiger snakes. Rates of evolution in body mass and snout-vent length, calculated for seven island snake populations, were significantly correlated with isolation time. Finally, rates of evolution in body mass per generation were significantly correlated with levels of plasticity in head growth rates. This study shows that body size evolution occurs at a faster pace in recently isolated populations and suggests that the level of adaptive plasticity for swallowing abilities may correlate with rates of body mass evolution. I hypothesise that, in the early stages of colonisation, adaptive plasticity and directional selection may combine and generate accelerated evolution towards an 'optimal' phenotype.

The Interplay between Environmental Filtering and Spatial Processes in Structuring Communities: The Case of Neotropical Snake Communities

Both habitat filters and spatial processes can influence community structure. Space alone affects species immigration from the regional species pool, whereas habitat filters affect species distribution and inter-specific interactions. This study aimed to understand how the interplay between environmental and geographical processes influenced the structure of Neotropical snake communities in different habitat types. We selected six studies that sampled snakes in forests, four conducted in savannas and two in grasslands (the latter two are grouped in a non-forest category). We used the net relatedness and nearest taxon indices to assess phylogenetic structure within forest and non-forest areas. We also used the phylogenetic fuzzy-weighting algorithm to characterize phylogenetic structure across communities and the relation of phylogenetic composition patterns to habitat type, structure, and latitude. Finally, we tested for morphological trait convergence and phylogenetic niche conservatism using four forest and four non-forest areas for which morphological data were available. Community phylogenetic composition changed across forest and non-forest areas suggesting that environmental filtering influences community structure. Species traits were affected by habitat type, indicating convergence at the metacommunity level. Tail length, robustness, and number of ventral scales maximized community convergence among forest and non-forest areas. The observed patterns suggested environmental filtering, indicating that less vertically structured habitats represent a strong filter. Despite the fact that phylogenetic structure was not detected individually for each community, we observed a trend towards communities composed by more closely related species in higher latitudes and more overdispersed compositions in lower latitudes. Such pattern suggests that the limited distribution of major snake lineages constrained species distributions. Structure indices for each community were also related to habitat type, showing that communities from non-forest areas tend to be more clustered. Our study showed that both environmental filtering and spatial gradients play important roles in shaping the composition of Neotropical snake communities.

Temperate snake community in South America: is diet determined by phylogeny or ecology?

Communities are complex and dynamic systems that change with time. The first attempts to explain how they were structured involve contemporary phenomena like ecological interactions between species (e.g., competition and predation) and led to the competition-predation hypothesis. Recently, the deep history hypothesis has emerged, which suggests that profound differences in the evolutionary history of organisms resulted in a number of ecological features that remain largely on species that are part of existing communities. Nevertheless, both phylogenetic structure and ecological interactions can act together to determine the structure of a community. Because diet is one of the main niche axes, in this study we evaluated, for the first time, the impact of ecological and phylogenetic factors on the diet of Neotropical snakes from the subtropical-temperate region of South America. Additionally, we studied their relationship with morphological and environmental aspects to understand the natural history and ecology of this community. A canonical phylogenetical ordination analysis showed that phylogeny explained most of the variation in diet, whereas ecological characters explained very little of this variation. Furthermore, some snakes that shared the habitat showed some degree of diet convergence, in accordance with the competition-predation hypothesis, although phylogeny remained the major determinant in structuring this community. The clade with the greatest variability was the subfamily Dipsadinae, whose members had a very different type of diet, based on soft-bodied invertebrates. Our results are consistent with the deep history hypothesis, and we suggest that the community under study has a deep phylogenetic effect that explains most of the variation in the diet.

Integrated analyses resolve conflicts over squamate reptile phylogeny and reveal unexpected placements for fossil taxa

Squamate reptiles (lizards and snakes) are a pivotal group whose relationships have become increasingly controversial. Squamates include >9000 species, making them the second largest group of terrestrial vertebrates. They are important medicinally and as model systems for ecological and evolutionary research. However, studies of squamate biology are hindered by uncertainty over their relationships, and some consider squamate phylogeny unresolved, given recent conflicts between molecular and morphological results. To resolve these conflicts, we expand existing morphological and molecular datasets for squamates (691 morphological characters and 46 genes, for 161 living and 49 fossil taxa, including a new set of 81 morphological characters and adding two genes from published studies) and perform integrated analyses. Our results resolve higher-level relationships as indicated by molecular analyses, and reveal hidden morphological support for the molecular hypothesis (but not vice-versa). Furthermore, we find that integrating molecular, morphological, and paleontological data leads to surprising placements for two major fossil clades (Mosasauria and Polyglyphanodontia). These results further demonstrate the importance of combining fossil and molecular information, and the potential problems of estimating the placement of fossil taxa from morphological data alone. Thus, our results caution against estimating fossil relationships without considering relevant molecular data, and against placing fossils into molecular trees (e.g. for dating analyses) without considering the possible impact of molecular data on their placement.

[DNA barcoding and its utility in commonly-used medicinal snakes]

Identification accuracy of traditional Chinese medicine is crucial for the traditional Chinese medicine research, production and application. DNA barcoding based on the mitochondrial gene coding for cytochrome c oxidase subunit I (COI), are more and more used for identification of traditional Chinese medicine. Using universal barcoding primers to sequence, we discussed the feasibility of DNA barcoding method for identification commonly-used medicinal snakes (a total of 109 samples belonging to 19 species 15 genera 6 families). The phylogenetic trees using Neighbor-joining were constructed. The results indicated that the mean content of G + C(46.5%) was lower than that of A + T (53.5%). As calculated by Kimera-2-parameter model, the mean intraspecies genetic distance of Trimeresurus albolabris, Ptyas dhumnades and Lycodon rufozonatus was greater than 2%. Further phylogenetic relationship results suggested that identification of one sample of T. albolabris was erroneous. The identification of some samples of P. dhumnades was also not correct, namely originally P. korros was identified as P. dhumnades. Factors influence on intraspecific genetic distance difference of L. rufozonatus need to be studied further. Therefore, DNA barcoding for identification of medicinal snakes is feasible, and greatly complements the morphological classification method. It is necessary to further study in identification of traditional Chinese medicine.

Heterochrony and early left-right asymmetry in the development of the cardiorespiratory system of snakes

Snake lungs show a remarkable diversity of organ asymmetries. The right lung is always fully developed, while the left lung is either absent, vestigial, or well-developed (but smaller than the right). A 'tracheal lung' is present in some taxa. These asymmetries are reflected in the pulmonary arteries. Lung asymmetry is known to appear at early stages of development in Thamnophis radix and Natrix natrix. Unfortunately, there is no developmental data on snakes with a well-developed or absent left lung. We examine the adult and developmental morphology of the lung and pulmonary arteries in the snakes Python curtus breitensteini, Pantherophis guttata guttata, Elaphe obsoleta spiloides, Calloselasma rhodostoma and Causus rhombeatus using gross dissection, MicroCT scanning and 3D reconstruction. We find that the right and tracheal lung develop similarly in these species. By contrast, the left lung either: (1) fails to develop; (2) elongates more slowly and aborts early without (2a) or with (2b) subsequent development of faveoli; (3) or develops normally. A right pulmonary artery always develops, but the left develops only if the left lung develops. No pulmonary artery develops in relation to the tracheal lung. We conclude that heterochrony in lung bud development contributes to lung asymmetry in several snake taxa. Secondly, the development of the pulmonary arteries is asymmetric at early stages, possibly because the splanchnic plexus fails to develop when the left lung is reduced. Finally, some changes in the topography of the pulmonary arteries are consequent on ontogenetic displacement of the heart down the body. Our findings show that the left-right asymmetry in the cardiorespiratory system of snakes is expressed early in development and may become phenotypically expressed through heterochronic shifts in growth, and changes in axial relations of organs and vessels. We propose a step-wise model for reduction of the left lung during snake evolution.

Early mesozoic coexistence of amniotes and hepadnaviridae

Hepadnaviridae are double-stranded DNA viruses that infect some species of birds and mammals. This includes humans, where hepatitis B viruses (HBVs) are prevalent pathogens in considerable parts of the global population. Recently, endogenized sequences of HBVs (eHBVs) have been discovered in bird genomes where they constitute direct evidence for the coexistence of these viruses and their hosts from the late Mesozoic until present. Nevertheless, virtually nothing is known about the ancient host range of this virus family in other animals. Here we report the first eHBVs from crocodilian, snake, and turtle genomes, including a turtle eHBV that endogenized >207 million years ago. This genomic “fossil” is >125 million years older than the oldest avian eHBV and provides the first direct evidence that Hepadnaviridae already existed during the Early Mesozoic. This implies that the Mesozoic fossil record of HBV infection spans three of the five major groups of land vertebrates, namely birds, crocodilians, and turtles. We show that the deep phylogenetic relationships of HBVs are largely congruent with the deep phylogeny of their amniote hosts, which suggests an ancient amniote–HBV coexistence and codivergence, at least since the Early Mesozoic. Notably, the organization of overlapping genes as well as the structure of elements involved in viral replication has remained highly conserved among HBVs along that time span, except for the presence of the X gene. We provide multiple lines of evidence that the tumor-promoting X protein of mammalian HBVs lacks a homolog in all other hepadnaviruses and propose a novel scenario for the emergence of X via segmental duplication and overprinting of pre-existing reading frames in the ancestor of mammalian HBVs. Our study reveals an unforeseen host range of prehistoric HBVs and provides novel insights into the genome evolution of hepadnaviruses throughout their long-lasting association with amniote hosts.

Viruses are not known to leave physical fossil traces, which makes our understanding of their evolutionary prehistory crucially dependent on the detection of endogenous viruses. Ancient endogenous viruses, also known as paleoviruses, are relics of viral genomes or fragments thereof that once infiltrated their host's germline and then remained as molecular “fossils” within the host genome. The massive genome sequencing of recent years has unearthed vast numbers of paleoviruses from various animal genomes, including the first endogenous hepatitis B viruses (eHBVs) in bird genomes. We screened genomes of land vertebrates (amniotes) for the presence of paleoviruses and identified ancient eHBVs in the recently sequenced genomes of crocodilians, snakes, and turtles. We report an eHBV that is >207 million years old, making it the oldest endogenous virus currently known. Furthermore, our results provide direct evidence that the Hepadnaviridae virus family infected birds, crocodilians and turtles during the Mesozoic Era, and suggest a long-lasting coexistence of these viruses and their amniote hosts at least since the Early Mesozoic. We challenge previous views on the origin of the oncogenic X gene and provide an evolutionary explanation as to why only mammalian hepatitis B infection leads to hepatocellular carcinoma.

[Application of rapid PCR to authenticate medicinal snakes]

To obtained an accurate, rapid and efficient method for authenticate medicinal snakes listed in Chinese Pharmacopoeia (Zaocysd humnades, Bungarus multicinctus, Agkistrodon acutus), a rapid PCR method for authenticate snakes and its adulterants was established based on the classic molecular authentication methods. DNA was extracted by alkaline lysis and the specific primers were amplified by two-steps PCR amplification method. The denatured and annealing temperature and cycle numbers were optimized. When 100 x SYBR Green I was added in the PCR product, strong green fluorescence was visualized under 365 nm UV whereas adulterants without. The whole process can complete in 30-45 minutes. The established method provides the technical support for authentication of the snakes on field.

Rediscovery of an endemic vertebrate from the remote Islas Revillagigedo in the eastern Pacific Ocean: the Clarión nightsnake lost and found

Vertebrates are currently going extinct at an alarming rate, largely because of habitat loss, global warming, infectious diseases, and human introductions. Island ecosystems are particularly vulnerable to invasive species and other ecological disturbances. Properly documenting historic and current species distributions is critical for quantifying extinction events. Museum specimens, field notes, and other archived materials from historical expeditions are essential for documenting recent changes in biodiversity. The Islas Revillagigedo are a remote group of four islands, 700–1100 km off the western coast of mainland México. The islands are home to many endemic plants and animals recognized at the specific- and subspecific-levels, several of which are currently threatened or have already gone extinct. Here, we recount the initial discovery of an endemic snake Hypsiglena ochrorhyncha unaocularus Tanner on Isla Clarión, the later dismissal of its existence, its absence from decades of field surveys, our recent rediscovery, and recognition of it as a distinct species. We collected two novel complete mitochondrial (mt) DNA genomes and up to 2800 base-pairs of mtDNA from several other individuals, aligned these with previously published mt-genome data from samples throughout the range of Hypsiglena, and conducted phylogenetic analyses to infer the biogeographic origin and taxonomic status of this population. We found the Isla Clarión population to be most closely related to populations in the Sonora–Sinaloa state border area of mainland México and Isla Santa Catalina, in the Gulf of California. Based on genetics, morphology, and geographic distributions, we also recognize these two other lineages as distinct species. Our study shows the importance of museum specimens, field notes, and careful surveys to accurately document biodiversity and brings these island endemics (Clarión and Santa Catalina nightsnakes) and mainland population near the Sonora–Sinaloa state border to the attention of conservation biologists currently monitoring biodiversity in these fragile subtropical ecosystems.