Search for: All records

Abstract Snake venoms are complex mixtures of toxic proteins that hold significant medical, pharmacological and evolutionary interest. To better understand the genetic diversity underlying snake venoms, we developed VenomCap, a novel exon‐capture probe set targeting toxin‐coding genes from a wide range of elapid snakes, with a particular focus on the ecologically diverse and medically important subfamily Hydrophiinae. We tested the capture success of VenomCap across 24 species, representing all major elapid lineages. We included snake phylogenomic probes in the VenomCap capture set, allowing us to compare capture performance between venom and phylogenomic loci and to infer elapid phylogenetic relationships. We demonstrated VenomCap's ability to recover exons from ~1500 target markers, representing a total of 24 known venom gene families, which includes the dominant gene families found in elapid venoms. We find that VenomCap's capture results are robust across all elapids sampled, and especially among hydrophiines, with respect to measures of target capture success (target loci matched, sensitivity, specificity and missing data). As a cost‐effective and efficient alternative to full genome sequencing, VenomCap can dramatically accelerate the sequencing and analysis of venom gene families. Overall, our tool offers a model for genomic studies on snake venom gene diversity and evolution that can be expanded for comprehensive comparisons across the other families of venomous snakes. 

Characterization of species distributions is a fundamental challenge in biodiversity science, with particular significance for downstream evolutionary studies, conservation efforts, field-based faunal studies and estimates of species diversity. Checklists and phylogenetic studies often focus on poorly known, rare taxa with limited ranges. However, studies of widely distributed, ecologically important species that are abundant in their preferred microhabitats are also important for systematics and local conservation efforts, but less often studied. We collected novel natural history data during fieldwork (2019–2023) for Philippine populations of bockadams (Homalopsidae:Cerberus), one of the most abundant vertebrates in Southeast Asian aquatic systems. Considered a coastal snake, many studies reportCerberusinland. We report the frequency of encounters ofCerberus schneiderii, and the IUCN data-deficient, Philippine-endemicCerberus microlepisduring six expeditions (62 days; 1041 person-hours). We report new occurrence data for 69C. schneideriiand6 C. microlepisfor coastal and inland populations, water measurements and dietary observations. Regression analyses and ecological niche models show the importance of coastal and mangrove habitats forCerberus. Our study is the most comprehensive assessment of PhilippineCerberuspopulations to date and provides critical baseline natural history data for downstream research on widespread and range-restricted species of Southeast Asian snakes. 

Abstract We present genome assemblies for 18 snake species representing 18 families (Serpentes: Caenophidia): Acrochordus granulatus, Aparallactus werneri, Boaedon fuliginosus, Calamaria suluensis, Cerberus rynchops, Grayia smithii, Imantodes cenchoa, Mimophis mahfalensis, Oxyrhabdium leporinum, Pareas carinatus, Psammodynastes pulverulentus, Pseudoxenodon macrops, Pseudoxyrhopus heterurus, Sibynophis collaris, Stegonotus admiraltiensis, Toxicocalamus goodenoughensis, Trimeresurus albolabris, and Tropidonophis doriae. From these new genome assemblies, we extracted thousands of loci commonly used in systematic and phylogenomic studies on snakes, including target-capture datasets composed of ultraconserved elements (UCEs) and anchored hybrid enriched loci (AHEs), as well as traditional Sanger loci. Phylogenies inferred from the two target-capture loci datasets were identical with each other and strongly congruent with previously published snake phylogenies. To show the additional utility of these non-model genomes for investigative evolutionary research, we mined the genome assemblies of two New Guinea island endemics in our dataset (S. admiraltiensis and T. doriae) for the ATP1a3 gene, a thoroughly researched indicator of resistance to toad toxin ingestion by squamates. We find that both these snakes possess the genotype for toad toxin resistance despite their endemism to New Guinea, a region absent of any toads until the human-mediated introduction of Cane Toads in the 1930s. These species possess identical substitutions that suggest the same bufotoxin resistance as their Australian congenerics (Stegonotus australis and Tropidonophis mairii) which forage on invasive Cane Toads. Herein, we show the utility of short-read high-coverage genomes, as well as improving the deficit of available squamate genomes with associated voucher specimens. 

Abstract Sunda-Papuan keelback snakes (Serpentes: Natricidae: Tropidonophis Jan 1863) include 20 species distributed from the Philippines south-east through the Moluccas to New Guinea and Australia. Diversity of this insular snake lineage peaks on the island of New Guinea. Previous phylogenetic studies incorporating Tropidonophis have been limited to multi-locus Sanger-sequenced datasets with broad squamate or family-level focus. We used a targeted-sequence capture approach to sequence thousands of nuclear ultraconserved elements (UCEs) to construct the most comprehensive sequence-based phylogenetic hypothesis for this genus and estimate ancestral biogeography. Phylogenies indicate the genus is monophyletic given recent taxonomic reassignment of Rhabdophis spilogaster to Tropidonophis. All UCE phylogenies recovered a monophyletic Tropidonophis with reciprocally monophyletic Philippine and New Guinean clades. Divergence dating and ancestral range estimation suggest dispersal to New Guinea from the Philippines to have occurred during the Mid-Miocene via the Oceanic Arc Terranes. From Late Miocene into the Pliocene the genus experienced rapid diversification from orogeny of the New Guinean Central Cordillera from Oceanic Arc Terrane accretion on the northern boundary of the Sahul Shelf. Future collecting of missing taxa from the Moluccas and Indonesian Papua will better the understanding of non-volant faunal biogeography and diversification in this tectonically complex Pacific arena. 

Abstract Speciation rates vary substantially across the tree of life. These rates should be linked to the rate at which population structure forms if a continuum between micro and macroevolutionary patterns exists. Previous studies examining the link between speciation rates and the degree of population formation in clades have been shown to be either correlated or uncorrelated depending on the group, but no study has yet examined the relationship between speciation rates and population structure in a young group that is constrained spatially to a single‐island system. We examine this correlation in 109 gemsnakes (Pseudoxyrhophiidae) endemic to Madagascar and originating in the early Miocene, which helps control for extinction variation across time and space. We find no relationship between rates of speciation and the formation rates of population structure over space in 33 species of gemsnakes. Rates of speciation show low variation, yet population structure varies widely across species, indicating that speciation rates and population structure are disconnected. We suspect this is largely due to the persistence of some lineages not susceptible to extinction. Importantly, we discuss how delimiting populations versus species may contribute to problems understanding the continuum between shallow and deep evolutionary processes. 

Over the past decade, museum genomics studies have focused on obtaining DNA of sufficient quality and quantity for sequencing from fluid-preserved natural history specimens, primarily to be used in systematic studies. While these studies have opened windows to evolutionary and biodiversity knowledge of many species worldwide, published works often focus on the success of these DNA sequencing efforts, which is undoubtedly less common than obtaining minimal or sometimes no DNA or unusable sequence data from specimens in natural history collections. Here, we attempt to obtain and sequence DNA extracts from 115 fresh and 41 degraded samples of homalopsid snakes, as well as from two degraded samples of a poorly known snake, Hydrablabes periops . Hydrablabes has been suggested to belong to at least two different families (Natricidae and Homalopsidae) and with no fresh tissues known to be available, intractable museum specimens currently provide the only opportunity to determine this snake’s taxonomic affinity. Although our aim was to generate a target-capture dataset for these samples, to be included in a broader phylogenetic study, results were less than ideal due to large amounts of missing data, especially using the same downstream methods as with standard, high-quality samples. However, rather than discount results entirely, we used mapping methods with references and pseudoreferences, along with phylogenetic analyses, to maximize any usable molecular data from our sequencing efforts, identify the taxonomic affinity of H. periops , and compare sequencing success between fresh and degraded tissue samples. This resulted in largely complete mitochondrial genomes for five specimens and hundreds to thousands of nuclear loci (ultra-conserved loci, anchored-hybrid enrichment loci, and a variety of loci frequently used in squamate phylogenetic studies) from fluid-preserved snakes, including a specimen of H. periops from the Field Museum of Natural History collection. We combined our H. periops data with previously published genomic and Sanger-sequenced datasets to confirm the familial designation of this taxon, reject previous taxonomic hypotheses, and make biogeographic inferences for Hydrablabes . A second H. periops specimen, despite being seemingly similar for initial raw sequencing results and after being put through the same protocols, resulted in little usable molecular data. We discuss the successes and failures of using different pipelines and methods to maximize the products from these data and provide expectations for others who are looking to use DNA sequencing efforts on specimens that likely have degraded DNA. Life Science Identifier ( Hydrablabes periops ) urn:lsid:zoobank.org :pub:F2AA44 E2-D2EF-4747-972A-652C34C2C09D. 

Our knowledge of the biodiversity of Asia and Australasia continues to expand with more focused studies on systematics of various groups and their biogeography. Historically, fluctuating sea levels and cyclic connection and separation of now-disjunct landmasses have been invoked to explain the accumulation of biodiversity via species pump mechanisms. However, recent research has shown that geological shifts of the mainland and species dispersal events may be better explanations of the biodiversity in these regions. We investigate these processes using the poorly studied and geographically widespread Mud Snakes (Serpentes: Homalopsidae) using a target capture approach of ~4,800 nuclear loci from fresh tissues and supplemental mitochondrial data from formalin tissues from museum specimens. We use these datasets to reconstruct the first resolved phylogeny of the group, identify their biogeographic origins, and test hypotheses regarding the roles of sea-level change and habitat selection on their diversification. Divergence dating and ancestral range estimation yielded support for an Oligocene origin and diversification from mainland Southeast Asia and Sundaland in the rear-fanged group ~20 million years ago, followed by eastward and westward dispersal. GeoHiSSE models indicate that niche expansion of ancestral, rear-fanged lineages into aquatic environments did not impact their diversification rates. Our results highlight that Pleistocene sea-level changes and habitat specificity did not primarily lead to the extant species richness of Homalopsidae and that, alternatively, geological shifts in mainland Southeast Asia may have been a major driver of diversity in this group. We also emphasize the importance of using fresh and degraded tissues, and both nuclear and mitochondrial DNA, for filling knowledge gaps in poorly known but highly diverse and conceptually important groups. Here, Homalopsidae represents a non-traditional but effective model study system for understanding transitions between terrestrial, marine, and freshwater environments. 

Abstract Divergence dating analyses in systematics provide a framework to develop and test biogeographic hypotheses regarding speciation. However, as molecular datasets grow from multilocus to genomic, sample sizes decrease due to computational burdens, and the testing of fine-scale biogeographic hypotheses becomes difficult. In this study, we use coalescent demographic models to investigate the diversification of poorly known rice paddy snakes from Southeast Asia (Homalopsidae:Hypsiscopus), which have conflicting dates of origin based on previous studies. We use coalescent modeling to test the hypothesis thatHypsiscopusdiversified 2.5 mya during the Khorat Plateau uplift in Thailand. Additionally, we use ecological niche analyses to identify potential differences in the niche space of the two most widely distributed species in the past and present. Our results suggestHypsiscopusdiversified ~ 2.4 mya, supporting that the Khorat Plateau may have initiated the diversification of rice paddy snakes. We also find significant niche differentiation and shifts between species ofHypsiscopus, indicating that environmental differences may have sustained differentiation of this genus after the Khorat Plateau uplift. Our study expands on the diversification history of snakes in Southeast Asia, and highlights how results from smaller multilocus datasets can be useful in developing and testing biogeographic hypotheses alongside genomic datasets. 

ABSTRACT AimThe Indo‐Australian Archipelago (IAA) is one of the most geologically complex and species‐rich regions on Earth. However, our knowledge of the geological processes and dispersal mechanisms that generate archipelago‐wide distributions across the IAA is limited to a few vertebrate groups and often solely attributed to Plio‐Pleistocene sea‐level fluctuations. In this study, we use phylogenomics‐based analyses to investigate two speciose and closely related genera of snakes,LycodonandStegonotus(Serpentes: Colubridae), which are widely distributed across the IAA to identify which biogeographic and environmental processes have shaped snake diversity in this region. LocationSouth Asia and the Indo‐Australian Archipelago (Indochina, Sundaland, Philippines, Wallacea, Australasia). TaxonSnakes (Colubridae:Lycodon, Stegonotus). MethodsWe inferred a phylogeny using a genomic dataset consisting of ultraconserved elements, anchored hybrid enrichment loci, and protein‐coding genes (~5400 nuclear loci) from 38 species (154 samples) ofLycodonandStegonotus. We used ancestral range estimation analysis to identify dispersal patterns across the IAA. Additionally, we implemented ensemble species distribution models to identify potential hotspots ofLycodonandStegonotusspecies richness and determined the environmental influence on geographic distributions and species diversity. ResultsWe find that these snakes comprise six deeply divergent lineages (genera) that initially originated in Mainland Southeast Asia during the Oligocene. Diversification of these lineages is influenced not by Plio‐Pleistocene sea‐level fluctuations but by multiple historical processes, including in situ diversification, island hopping, long‐distance rafting, possible microcontinental block drifting, Cenozoic land bridge migrations and founder events. The species distribution models do not consistently estimate lower or higher species richness in any particular region within the IAA, but precipitation overall is considered an important factor in estimated species richness. Main ConclusionAlthough Plio‐Pleistocene sea‐level fluctuations are notorious species pump diversification paradigms in the IAA, multiple and more ancient geological processes and environmental factors contributed to current diversity levels and distributions. Our approach expands future investigations of alternative hypotheses of biodiversity sources in the IAA and greatly expands the diversity of causal mechanisms for discussions of terrestrial Southeast Asian biodiversity beyond dispersal versus vicariance hypotheses.