Search

Both the metabolic theory of ecology and dynamic energy budget theory predict that climate influences body size through its effects on first-order determinants of energetics: reactive temperatures, carbon resources and oxygen availability. Although oxygen is seldom limiting in terrestrial systems, temperature and resources vary spatially. We used redundancy analyses and variation partitioning to evaluate the influence of climatic temperature, precipitation and their seasonalities on multivariate body size across the distributions of four species of the western rattlesnake group in North America ( Crotalus pyrrhus , C. scutulatus , C. oreganus and C. viridis ). Most species showed a pattern of increased body size in cooler, mesic climates and decreased body size in warmer, xeric climates. Exceptions to the pattern provided additional context through climatic idiosyncrasies in the distributions of each species. For example, the general pattern of a negative influence of temperature on body size was not apparent for C. oreganus , which ranges across the mildest climates overall among the four species. In contrast to previous studies, we found that seasonality had negligible effects on body size. We suggest that precipitation gradients correlate positively with resource availability in driving intraspecific body size and that temperature compounds this gradient by increasing baseline metabolic demands and restricting activity in particularly warm or otherwise extreme climates., Competing Interests: We declare we have no competing interests., (© 2024 The Authors.)

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interests with the contents of this article., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Snakebite envenoming remains a devastating and neglected tropical disease, claiming over 100,000 lives annually and causing severe complications and long-lasting disabilities for many more 1,2 . Three-finger toxins (3FTx) are highly toxic components of elapid snake venoms that can cause diverse pathologies, including severe tissue damage 3 and inhibition of nicotinic acetylcholine receptors (nAChRs) resulting in life-threatening neurotoxicity 4 . Currently, the only available treatments for snakebite consist of polyclonal antibodies derived from the plasma of immunized animals, which have high cost and limited efficacy against 3FTxs 5,6,7 . Here, we use deep learning methods to de novo design proteins to bind short- and long-chain α-neurotoxins and cytotoxins from the 3FTx family. With limited experimental screening, we obtain protein designs with remarkable thermal stability, high binding affinity, and near-atomic level agreement with the computational models. The designed proteins effectively neutralize all three 3FTx sub-families in vitro and protect mice from a lethal neurotoxin challenge. Such potent, stable, and readily manufacturable toxin-neutralizing proteins could provide the basis for safer, cost-effective, and widely accessible next-generation antivenom therapeutics. Beyond snakebite, our computational design methodology should help democratize therapeutic discovery, particularly in resource-limited settings, by substantially reducing costs and resource requirements for development of therapies to neglected tropical diseases.

Background: Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris)., Results: Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis., Conclusions: Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype., (© 2024. The Author(s).)

Predation has the potential to impart strong selective pressures on organisms within their environments, resulting in adaptive changes in prey that minimize risk of predation. Pressures from venomous snakes present an exceptional challenge to prey, as venom represents a unique chemical arsenal evolutionarily tailored to incapacitate prey. In response, venom resistance has been detected in various snake prey species, and to varying degrees. This study analyzes venom resistance in an eastern Colorado grassland habitat, where the Prairie Rattlesnake ( Crotalus viridis ) and Desert Massasauga Rattlesnake ( Sistrurus tergeminus edwardsii ) co-occur with a suite of grassland rodents. We test for venom resistance across rodent and snake pairings using two geographically distant field sites to determine the role of 1) predation pressure and trophic ecology, and 2) sympatric and allopatric patterns of venom resistance. Resistance was measured using serum-based metalloproteinase inhibition assays to determine potential inhibition of proteolytic activity, augmented by median lethal dose (LD 50 ) assays on rodent species to assess toxicity of crude venoms. Resistance is present in several rodent species, with strong resistance present in populations of Eastern Woodrat ( Neotoma floridana ), Ord's Kangaroo Rat ( Dipodomys ordii ), and Northern Grasshopper Mouse ( Onychomys leucogaster ). Resistance is less developed in other species, including the House Mouse ( Mus musculus ) and Plains Pocket Mouse ( Perognathus flavescens ). An unexpected differential is present, where Lincoln County Kangaroo Rats are highly resistant to venom of co-occurring Prairie Rattlesnakes yet are sensitive to an allopatric population of Prairie Rattlesnakes in Weld County. Lincoln Co. Northern Grasshopper Mice also demonstrate extremely elevated resistance to Weld Co. Prairie Rattlesnake venoms, and they may possess resistance mechanisms for myotoxin a, an abundant component of Weld Co. C. v viridis venoms. This study illustrates the complexity of venom resistance in biological communities that can exist when incorporating multiple species interactions. Future studies aimed at characterizing resistance mechanisms at the molecular level will provide a more detailed physiological context for understanding mechanisms by which resistance to venoms occurs., Competing Interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Neil Balchan reports financial support was provided by 10.13039/100007666University of Northern Colorado., (© 2023 The Authors.)

Studying the consequences of hybridization between closely related species with divergent traits can reveal patterns of evolution that shape and maintain extreme trophic adaptations. Snake venoms are an excellent model system for examining the evolutionary and ecological patterns that underlie highly selected polymorphic traits. Here we investigate hybrid venom phenotypes that result from natural introgression between two rattlesnake species that express highly divergent venom phenotypes: Crotalus o. concolor and C. v. viridis. Though not yet documented, interbreeding between these species may lead to novel venom phenotypes with unique activities that break the typical trends of venom composition in rattlesnakes. The characteristics of these unusual phenotypes could unveil the roles of introgression in maintaining patterns of venom composition and variation, including the near ubiquitous dichotomy between neurotoxic or degradative venoms observed across rattlesnakes. We use RADseq data to infer patterns of gene flow and hybrid ancestry between these diverged lineages and link these genetic data with analyses of venom composition, biological activity, and whole animal model toxicity tests to understand the impacts of introgression on venom composition. We find that introgressed populations express admixed venom phenotypes that do not sacrifice biological activity (lethal toxicity) or overall abundance of dominant toxins compared to parental venoms. These hybridized venoms therefore do not represent a trade-off in functionality between the typical phenotypic extremes but instead represent a unique combination of characters whose expression appears limited to the hybrid zone., Competing Interests: Declaration of competing interest The authors have no competing interests to declare., (Copyright © 2023 Elsevier B.V. and Société Française de Biochimie et Biologie Moléculaire (SFBBM). All rights reserved.)

Background: Snake venoms are trophic adaptations that represent an ideal model to examine the evolutionary factors that shape polymorphic traits under strong natural selection. Venom compositional variation is substantial within and among venomous snake species. However, the forces shaping this phenotypic complexity, as well as the potential integrated roles of biotic and abiotic factors, have received little attention. Here, we investigate geographic variation in venom composition in a wide-ranging rattlesnake (Crotalus viridis viridis) and contextualize this variation by investigating dietary, phylogenetic, and environmental variables that covary with venom., Results: Using shotgun proteomics, venom biochemical profiling, and lethality assays, we identify 2 distinct divergent phenotypes that characterize major axes of venom variation in this species: a myotoxin-rich phenotype and a snake venom metalloprotease (SVMP)-rich phenotype. We find that dietary availability and temperature-related abiotic factors are correlated with geographic trends in venom composition., Conclusions: Our findings highlight the potential for snake venoms to vary extensively within species, for this variation to be driven by biotic and abiotic factors, and for the importance of integrating biotic and abiotic variation for understanding complex trait evolution. Links between venom variation and variation in biotic and abiotic factors indicate that venom variation likely results from substantial geographic variation in selection regimes that determine the efficacy of venom phenotypes across populations and snake species. Our results highlight the cascading influence of abiotic factors on biotic factors that ultimately shape venom phenotype, providing evidence for a central role of local selection as a key driver of venom variation., (© 2023. The Author(s).)

The ubiquitous cellular heterogeneity underlying many organism-level phenotypes raises questions about what factors drive this heterogeneity and how these complex heterogeneous systems evolve. Here, we use single-cell expression data from a Prairie rattlesnake (Crotalus viridis) venom gland to evaluate hypotheses for signaling networks underlying snake venom regulation and the degree to which different venom gene families have evolutionarily recruited distinct regulatory architectures. Our findings suggest that snake venom regulatory systems have evolutionarily co-opted trans-regulatory factors from extracellular signal-regulated kinase and unfolded protein response pathways that specifically coordinate expression of distinct venom toxins in a phased sequence across a single population of secretory cells. This pattern of co-option results in extensive cell-to-cell variation in venom gene expression, even between tandemly duplicated paralogs, suggesting this regulatory architecture has evolved to circumvent cellular constraints. While the exact nature of such constraints remains an open question, we propose that such regulatory heterogeneity may circumvent steric constraints on chromatin, cellular physiological constraints (e.g., endoplasmic reticulum stress or negative protein-protein interactions), or a combination of these. Regardless of the precise nature of these constraints, this example suggests that, in some cases, dynamic cellular constraints may impose previously unappreciated secondary constraints on the evolution of gene regulatory networks that favors heterogeneous expression., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Biofluid proteomics is a sensitive and high throughput technique that provides vast amounts of molecular data for biomarker discovery. More recently, dried blood spots (DBS) have gained traction as a stable, noninvasive, and relatively cheap source of proteomic data for biomarker identification in disease and injury. Snake envenomation is responsible for significant morbidity and mortality worldwide; however, much remains unknown about the systemic molecular response to envenomation and acquiring biological samples for analysis is a major hurdle. In this study, we utilized DBS acquired from a case of lethal rattlesnake envenomation to determine the feasibility of discovering biomarkers associated with human envenomation. We identified proteins that were either unique or upregulated in envenomated blood compared to non-envenomated blood and evaluated if physiological response pathways and protein markers that correspond to the observed syndromes triggered by envenomation could be detected. We demonstrate that DBS provide useful proteomic information on the systemic processes that resulted from envenomation in this case and find evidence for a massive and systemic inflammatory cascade, combined with coagulation dysregulation, complement system activation, hypoxia response activation, and apoptosis. We also detected potential markers indicative of lethal anaphylaxis, cardiac arrest, and brain death. Ultimately, DBS proteomics has the potential to provide stable and sensitive molecular data on envenomation syndromes and response pathways, which is particularly relevant in low-resource areas which may lack the materials for biofluid processing and storage., Competing Interests: The authors declare no conflict of interest., (© 2023 The Authors. Analytical Science Advances published by Wiley‐VCH GmbH.)

Infections caused by multi-drug-resistant (MDR) bacteria are a global threat to human health. As venoms are the source of biochemically diverse bioactive proteins and peptides, we investigated the antimicrobial activity and murine skin infection model-based wound healing efficacy of a 13 kDa protein. The active component PaTx-II was isolated from the venom of Pseudechis australis (Australian King Brown or Mulga Snake). PaTx-II inhibited the growth of Gram-positive bacteria in vitro, with moderate potency (MICs of 25 µM) observed against S. aureus , E. aerogenes, and P. vulgaris . The antibiotic activity of PaTx-II was associated with the disruption of membrane integrity, pore formation, and lysis of bacterial cells, as evidenced by scanning and transmission microscopy. However, these effects were not observed with mammalian cells, and PaTx-II exhibited minimal cytotoxicity (CC 50 > 1000 µM) toward skin/lung cells. Antimicrobial efficacy was then determined using a murine model of S. aureus skin infection. Topical application of PaTx-II (0.5 mg/kg) cleared S. aureus with concomitant increased vascularization and re-epithelialization, promoting wound healing. As small proteins and peptides can possess immunomodulatory effects to enhance microbial clearance, cytokines and collagen from the wound tissue samples were analyzed by immunoblots and immunoassays. The amounts of type I collagen in PaTx-II-treated sites were elevated compared to the vehicle controls, suggesting a potential role for collagen in facilitating the maturation of the dermal matrix during wound healing. Levels of the proinflammatory cytokines interleukin-1β (IL-1β), interleukin-6 (IL-6) and tumor necrosis factor-α (TNF-α), cyclooxygenase-2 (COX-2) and interleukin-10 (IL-10), factors known to promote neovascularization, were substantially reduced by PaTx-II treatment. Further studies that characterize the contributions towards efficacy imparted by in vitro antimicrobial and immunomodulatory activity with PaTx-II are warranted.

The Crotalus intermedius group is a clade of rattlesnakes consisting of several species adapted to a high elevation habitat, primarily in México. Crotalus tancitarensis was previously classified as C. intermedius , until individuals occurring on Cerro Tancítaro in Michoacán, México, were reevaluated and classified as a new species ( C. tancitarensis ) based on scale pattern and geographic location. This study aimed to characterize the venom of C. tancitarensis and compare the venom profile to those of other species within the Crotalus intermedius group using gel electrophoresis, biochemical assays, reverse-phase high performance liquid chromatography, mass spectrometry, and lethal toxicity (LD 50 ) assays. Results show that the venom profiles of species within the Crotalus intermedius group are similar, but with distinct differences in phospholipase A 2 (PLA 2 ), metalloproteinase PI (SVMP PI), and kallikrein-like serine proteinase (SVSP) activity and relative abundance. Proteomic analysis indicated that the highland forms produce venoms with 50-60 protein isoforms and a composition typical of type I rattlesnake venoms (abundant SVMPs, lack of presynaptic PLA 2 -based neurotoxins), as well as a diversity of typical Crotalus venom components such as serine proteinases, PLA 2 s, C-type lectins, and less abundant toxins (LAAOs, CRiSPs, etc.). The overall venom profile of C. tancitarensis appears most similar to C. transversus , which is consistent with a previous mitochondrial DNA analysis of the Crotalus intermedius group. These rattlesnakes of the Mexican highlands represent a radiation of high elevation specialists, and in spite of divergence of species in these Sky Island habitats, venom composition of species analyzed here has remained relatively conserved. The majority of protein family isoforms are conserved in all members of the clade, and as seen in other more broadly distributed rattlesnake species, differences in their venoms are largely due to relative concentrations of specific components.

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system., (© 2022 The Authors. Evolution © 2022 The Society for the Study of Evolution.)

Convergence is the phenomenon whereby similar phenotypes evolve independently in different lineages. One example is resistance to toxins in animals. Toxins have evolved many times throughout the tree of life. They disrupt molecular and physiological pathways in target species, thereby incapacitating prey or deterring a predator. In response, molecular resistance has evolved in many species exposed to toxins to counteract their harmful effects. Here, we review current knowledge on the convergence of toxin resistance using examples from a wide range of toxin families. We explore the evolutionary processes and molecular adaptations driving toxin resistance. However, resistance adaptations may carry a fitness cost if they disrupt the normal physiology of the resistant animal. Therefore, there is a trade-off between maintaining a functional molecular target and reducing toxin susceptibility. There are relatively few solutions that satisfy this trade-off. As a result, we see a small set of molecular adaptations appearing repeatedly in diverse animal lineages, a phenomenon that is consistent with models of deterministic evolution. Convergence may also explain what has been called 'autoresistance'. This is often thought to have evolved for self-protection, but we argue instead that it may be a consequence of poisonous animals feeding on toxic prey. Toxin resistance provides a unique and compelling model system for studying the interplay between trophic interactions, selection pressures and the molecular mechanisms underlying evolutionary novelties., (© 2022 The Authors. Biological Reviews published by John Wiley & Sons Ltd on behalf of Cambridge Philosophical Society.)

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a "refugium" for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Crotamine, myotoxin a and homologs are short peptides that often comprise major fractions of rattlesnake venoms and have been extensively studied for their bioactive properties. These toxins are thought to be important for rapidly immobilizing mammalian prey and are implicated in serious, and sometimes fatal, responses to envenomation in humans. While high quality reference genomes for multiple venomous snakes are available, the loci that encode myotoxins have not been successfully assembled in any existing genome assembly. Here, we integrate new and existing genomic and transcriptomic data from the Prairie Rattlesnake (Crotalus viridis viridis) to reconstruct, characterize, and infer the chromosomal locations of myotoxin-encoding loci. We integrate long-read transcriptomics (Pacific Bioscience's Iso-Seq) and short-read RNA-seq to infer gene sequence diversity and characterize patterns of myotoxin and paralogous β-defensin expression across multiple tissues. We also identify two long non-coding RNA sequences which both encode functional myotoxins, demonstrating a newly discovered source of venom coding sequence diversity. We also integrate long-range mate-pair chromatin contact data and linked-read sequencing to infer the structure and chromosomal locations of the three myotoxin-like loci. Further, we conclude that the venom-associated myotoxin is located on chromosome 1 and is adjacent to non-venom paralogs. Consistent with this locus contributing to venom composition, we find evidence that the promoter of this gene is selectively open in venom gland tissue and contains transcription factor binding sites implicated in broad trans-regulatory pathways that regulate snake venoms. This study provides the best genomic reconstruction of myotoxin loci to date and raises questions about the physiological roles and interplay between myotoxin and related genes, as well as the genomic origins of snake venom variation., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

Jatropha mollissima is endemic to Brazil and is used for traditional medicinal purposes, including the treatment of snakebite. In this study, latex obtained from this plant was fractioned using reversed-phase chromatography, and the fractions were then screened for peptides. A 755 g/mol peptide was obtained, and MS/MS analyses indicated it had a cyclic sequence (Pro-Leu-Gly-Val-Leu-Leu-Tyr). This peptide sequence was present in the Jatropha genome database, and an identity value of 90.71%, an E-value of 0.0, and a score of 883 with NO-associated protein 1/chloroplastic/mitochondria of Jatropha curcas were obtained from the NCBI nonredundant protein sequence (nr) database. Molecular docking analyses performed with the peptide against a metalloendopeptidase belonging to Crotalus adamanteus snake venom suggested the cyclic peptide establishes favorable interactions with the catalytic site of the enzyme. Therefore, it could inhibit enzyme catalysis. This belief was corroborated by the formation of 6 hydrogen bonds with the linear form of the peptide. Tighter complexation of the cyclic form (41 kcal/mol more energetic) revealed better spatial blocking. The linear form outperformed the cyclic form in complexing the required energy, recruiting more catalytic residues (6/2), and in establishing more hydrogen bonds (6/3). However, cyclic folding provided a more significant spatial block within the catalytic site. The set of results suggests that the cycle peptide, here called Jatromollistatin, which was previously described as jatrophidin and pohlianin A in two other species of Jatropha, is a promising candidate to inhibit venom proteases. This belief is corroborated by the topical use of the latex for initial treatment of snakebites., (© 2022 John Wiley & Sons Ltd.)

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis -regulatory sequences (enhancers and promoters), trans -regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis -regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade., (© 2022 Perry et al.; Published by Cold Spring Harbor Laboratory Press.)

The venom glands of reptiles, particularly those of front-fanged advanced snakes, must satisfy conflicting biological demands: rapid synthesis of potentially labile and highly toxic proteins, storage in the gland lumen for long periods, stabilization of the stored secretions, immediate activation of toxins upon deployment and protection of the animal from the toxic effects of its own venom. This dynamic system could serve as a model for the study of a variety of different phenomena involving exocrine gland activation, protein synthesis, stabilization of protein products and secretory mechanisms. However, these studies have been hampered by a lack of a long-term model that can be propagated in the lab (as opposed to whole-animal studies). Numerous attempts have been made to extend the lifetime of venom gland secretory cells, but only recently has an organoid model been shown to have the requisite qualities of recapitulation of the native system, self-propagation and long-term viability (>1 year). A tractable model is now available for myriad cell- and molecular-level studies of venom glands, protein synthesis and secretion. However, venom glands of reptiles are not identical, and many differ very extensively in overall architecture, microanatomy and protein products produced. This Review summarizes the similarities among and differences between venom glands of helodermatid lizards and of rear-fanged and front-fanged snakes, highlighting those areas that are well understood and identifying areas where future studies can fill in significant gaps in knowledge of these ancient, yet fascinating systems., Competing Interests: Competing interests The author declares no competing or financial interests., (© 2022. Published by The Company of Biologists Ltd.)

Introduction: Snake venoms contain many protein and peptide isoforms with high levels of sequence variation, even within a single species., Areas Covered: In this review, we highlight several examples, from both published and unpublished work in our lab, demonstrating how a combined venom gland transcriptome and proteome methodology allows for comprehensive characterization of venoms, including those from understudied rear-fanged snake species, and we provide recommendations for using these approaches., Expert Opinion: When characterizing venoms, peptide mass fingerprinting using databases built predominately from protein sequences originating from model organisms can be disadvantageous, especially when the intention is to document protein diversity. Therefore, the use of species-specific venom gland transcriptomes corrects for the absence of these unique peptide sequences in databases. The integration of transcriptomics and proteomics improves the accuracy of either approach alone for venom profiling.

Epidemiological data on snakebite in the Brazilian state of Ceará are scarce, as the only report on this subject was last published in 1997. However, according to the Brazilian system of recording disease incidents (Sistema de Informação de Agravos de Notificação [SINAN]), more than 13,000 snakebites have been registered since 2001 in the state of Ceará, making this disease an important public health issue. In the present study, we evaluate the influence of environmental changes, including drought and desertification, on the risk of snakebite envenoming in the Brazilian northeastern state of Ceará. We compare public data on snakebites from Brazilian Epidemiological Surveillance System (DATASUS), rainfall records, advanced desertification maps, pastures and socioeconomic information of the 184 municipals in Ceará between 2001 and 2017. During the period of investigation, 8,945 snakebites were recorded, the majority (93.8%) of which involved venomous snakes. Almost half of the municipals (48%) had 100 incidences or more per 100,000 inhabitants. Data collected also highlight month-to-month occurrences of snakebites, with trends to rise shortly after the onset of precipitation, peaking in July and then trending downward as rainfall decreases, reaching the lowest level in December. We deduce an inverse relationship between Human Development Index (HDI) and snakebites per area. Spearman correlation and principal component analysis support the hypothesis that water scarcity and desertification are linked to increased risk of snakebite envenoming. Our study indicates that besides poverty, dry and desertified areas represent risk factors associated with increased incidence of snakebite envenoming in the state of Ceará., (© 2021 John Wiley & Sons Ltd.)

In the Southern Asian countries, snakebite takes a substantial toll in terms of human life, inflicts acute morbidity and long term disability both physical and psychological, and therefore represents a neglected socio-economic problem and severe health issue that requires immediate medical attention. The 'Big Four' venomous snakes, viz. Daboia russelii, Naja naja, Bungarus caeruleus and Echis carinatus, are prominent, medically important species and are the most dangerous snakes of this region; therefore, the commercial polyvalent antivenom (PAV) contains antibodies against the venoms of these snakes. However, envenomations by species other than the 'Big Four' snakes are grossly neglected, and PAV is only partially effective in neutralizing the venom of these snakes. Many issues confounding effective treatment of snakebite are discussed in this review, and these hurdles preventing successful treatment of snakebite must be addressed. However, in South Asian countries, the pre-hospital treatment and appropriate first aid are equally important to mitigate the problem of snakebite and therefore, these issues are also highlighted here. Further, this review suggests a roadmap and guidelines for the prevention of snakebite and improvement of hospital management of snakebite in these Southern Asian countries., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Male-biased mutation rates occur in a diverse array of organisms. The ratio of male-to-female mutation rate may have major ramifications for evolution across the genome, and for sex-linked genes in particular. In ZW species, the Z chromosome is carried by males two-thirds of the time, leading to the prediction that male-biased mutation rates will have a disproportionate effect on the evolution of Z-linked genes relative to autosomes and the W chromosome. Colubroid snakes (including colubrids, elapids, and viperids) have ZW sex determination, yet male-biased mutation rates have not been well studied in this group. Here we analyze a population genomic dataset from rattlesnakes to quantify genetic variation within and genetic divergence between species. We use a new method for unbiased estimation of population genetic summary statistics to compare variation between the Z chromosome and autosomes and to calculate net nucleotide differentiation between species. We find evidence for a 2.03-fold greater mutation rate in male rattlesnakes relative to females, corresponding to an average μZ/μA ratio of 1.1. Our results from snakes are quantitatively similar to birds, suggesting that male-biased mutation rates may be a common feature across vertebrate lineages with ZW sex determination., (© The American Genetic Association. 2021. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Species of Oxybelis are extremely elongate arboreal snakes that are broadly distributed in the Americas, from extreme southeastern Arizona (USA) to central South America. Primarily feeding on lizards and birds, Oxybelis venoms are poorly known in general, but a prominent taxon-specific three-finger toxin (fulgimotoxin) was isolated from and is a prominent component of O. fulgidus venom; a homolog is also present in O. aeneus venom. As part of ongoing characterization of venoms from rear-fanged snakes, we describe here the composition of two broadly distributed species, O. aeneus and O. fulgidus. Venom proteomes were of very low complexity, and four protein families (LAAO, PIII SVMP, CRiSP and 3FTx) account for more than 90% of total protein composition. Venoms from both species are moderately toxic to mice and to Hemidactylus geckos, but they are nearly an order of magnitude more toxic to Anolis lizards (a native prey species). These results reflect a trend in colubrid venom composition that is becoming increasingly more common: the presence of taxon-specific toxins, specifically three-finger toxins, preferentially targeting lizards and/or birds., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

It is assumed that toxin expression is equivalent between left and right glands of a single snake. In the current study, we report venoms that differ in enzyme functionality and overall composition between the left and right gland of a single snake. The right gland produced venom of comparable composition to venom previously extracted from the same individual; however, the left gland produced venom with overall lower protein content and considerably less enzyme activity., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Despite the extensive body of research on snake venom, many facets of snake venom systems, such as the physiology and regulation of the venom gland itself, remain virtually unstudied. Here, we use time series gene expression analyses of the rattlesnake venom gland in comparison with several non-venom tissues to characterize physiological and cellular processes associated with venom production and to highlight key distinctions of venom gland cellular and physiological function. We find consistent evidence for activation of stress response pathways in the venom gland, suggesting that mitigation of cellular stress is a crucial component of venom production. Additionally, we demonstrate evidence for an unappreciated degree of cellular and secretory activity in the steady state venom gland relative to other secretory tissues and identify vacuolar ATPases as the likely mechanisms driving acidification of the venom gland lumen during venom production and storage.

A case of midget-faded rattlesnake (Crotalus oreganus concolor) envenomation of an adult male professional herpetologist occurred in a rural setting and resulted in an array of venom induced myoneurologic symptoms. The patient experienced blurry vision, total body paresthesia, dyspnea, chest tightness, and waves of spastic muscle movements of the hands and feet that resembled tetany. It was not apparent whether these symptoms were potentially venom induced or were related to stress-induced physiologic responses. Local envenomation effects were minimal, and coagulation parameters remained within normal limits. Antivenom was not administered per patient concerns related to a history of acute allergic reactions to antivenom. Venom was collected from the Crotalus oreganus concolor responsible for the bite, and analysis revealed the presence of high levels of myotoxins (SR calcium pump antagonists) and concolor toxin, a presynaptic neurotoxin that can have myotoxic effects and cause respiratory paralysis; several serine proteinases associated with coagulopathies were also present in the venom profile., (Copyright © 2020 Wilderness Medical Society. Published by Elsevier Inc. All rights reserved.)

The genera Ophiophagus and Naja comprise part of a clade of snakes referred to as cobras, dangerously venomous front-fanged snakes in the family Elapidae responsible for significant human mortality and morbidity throughout Asia and Africa. We evaluated venom enzyme variation for eleven cobra species and three N. kaouthia populations using SDS-PAGE venom fingerprinting and numerous enzyme assays. Acetylcholinesterase and PLA 2 activities were the most variable between species, and PLA 2 activity was significantly different between Malaysian and Thailand N. kaouthia populations. Venom metalloproteinase activity was low and significantly different among most species, but levels were identical for N. kaouthia populations; minor variation in venom L-amino acid oxidase and phosphodiesterase activities were seen between cobra species. Naja siamensis venom lacked the α-fibrinogenolytic activity common to other cobra venoms. In addition, venom from N. siamensis had no detectable metalloproteinase activity and exhibited an SDS-PAGE profile with reduced abundance of higher mass proteins. Venom profiles from spitting cobras (N. siamensis, N. pallida, and N. mossambica) exhibited similar reductions in higher mass proteins, suggesting the evolution of venoms of reduced complexity and decreased enzymatic activity among spitting cobras. Generally, the venom proteomes of cobras show highly abundant three-finger toxin diversity, followed by large quantities of PLA 2 s. However, PLA 2 bands and activity were very reduced for N. haje, N. annulifera and N. nivea. Venom compositionalenzy analysis provides insight into the evolution, diversification and distribution of different venom phenotypes that complements venomic data, and this information is critical for the development of effective antivenoms and snakebite treatment., Competing Interests: Declaration of competing interest The authors declare that they have no conflict of interest associated with this manuscript., (Copyright © 2020. Published by Elsevier Inc.)

Rear-fanged colubrid snakes include hundreds of species globally that possess a Duvernoy's venom gland and often one-several enlarged rear maxillary teeth. We investigated the venom proteome of the Central American Lyre Snake (Trimorphodon quadruplex), a moderate-sized rear-fanged colubrid snake and the southernmost Trimorphodon, using a bottom-up proteomic approach coupled with enzyme and inhibitor assays, cytotoxicity assays and lethal toxicity assays. Several enzymes uncommonly observed in colubrid venoms were purified and characterized further. Trimorphodon quadruplex has a rather low complexity venome, typical of many rear-fanged snakes, but its venom contains L-amino acid oxidase, phospholipase A 2 , and a dimeric 3FTx, and 3FTxs dominate the proteome. Its PLA 2 is catalytically quite active, but it lacks myotoxicity or acute toxicity; LAAO exhibits conserved structure and appears to be highly labile. Several P-III metalloproteinases are present and hydrolyze azocasein and the α-subunit of fibrinogen but lack hemorrhagic activity. Trimorphodon quadruplex produces venom and retains constriction, utilizing both chemically-mediated and mechanical feeding modes. SIGNIFICANCE: We demonstrate that T. quadruplex venom proteins are similar to those found in front-fanged snake species are present but show different biological activities. Our results underscore the importance of considering the biological roles of venoms from more than a mammal-centric perspective., Competing Interests: Declaration of Competing Interest The authors have declared no conflict of interest., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions-a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions-a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Studying animal toxin evolution requires sequences of these proteins and peptides, and transcript sequences allow for the construction of cladograms and evaluation of selection pressures from nonsynonymous and synonymous nucleotide mutation ratios. In addition, these translated sequences can be useful as custom databases for peptide identifications within venoms and for better proteomic quantification. Obtaining these transcripts is achieved by sequencing cDNA originating from venom gland tissue or venom. This chapter provides the methodology for (1) targeted sequencing of transcripts from a single venom protein family (RNA isolation and 3'RACE [rapid amplification of cDNA ends]), (2) generation of a venom gland transcriptome with next-generation sequencing (NGS) technology (de novo transcriptome assembly, toxin transcript identification, quantification, and positive selection analysis), and (3) combined high-throughput proteomics to identify secreted venom components. Transcriptomics has become fundamental for studying toxin evolution, but it creates many challenges for scientists who are unfamiliar with working with RNA, managing large NGS datasets and executing the required programs, particularly considering that there is an overabundance of available software in this field and not all perform optimally for venom gland transcriptome assembly. This chapter provides one pipeline for the integration of both low- and high-throughput transcriptomics with proteomics to characterize venoms.

The Twin-spotted Rattlesnake (Crotalus pricei) is a small lizard specialist restricted to higher montane habitat in the Sky Islands of Arizona and México. Though this species is restricted to high elevations and dispersal between mountaintops is impossible, few studies have investigated venom composition or the predator-prey relationship between C. p. pricei and its primary prey source, Yarrow's Spiny Lizard (Sceloporus jarrovii). Because of current isolation of populations in disjunct mountain ranges, it is possible that populations show local adaptation to prey and/or environmental conditions, as reflected by distinct venom phenotypes. Here we characterize venom composition of C. pricei from several different Sky Island mountain ranges and its relationship with S. jarrovii by using various analytical techniques and comparative toxicity tests. Results of venom analyses indicate that there is limited geographic variation in venom composition, occurring primarily in venoms of C. pricei from Durango, México, which had unusual PLA 2 isoforms that are lacking from venoms of U.S. populations. Toxicity assays reveal that S. jarrovii has not developed resistance specific to C. p. pricei venom but does display a general tolerance to venom of several snakes in the genus Crotalus. These results provide insight into the evolutionary relationship between a lizard specialist and its natural prey, in addition to novel information on the venom composition of a little-studied species with a narrow range in the United States., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

The Mohave Rattlesnake (Crotalus scutulatus) is a highly venomous pitviper inhabiting the arid interior deserts, grasslands, and savannas of western North America. Currently two subspecies are recognized: the Northern Mohave Rattlesnake (C. s. scutulatus) ranging from southern California to the southern Central Mexican Plateau, and the Huamantla Rattlesnake (C. s. salvini) from the region of Tlaxcala, Veracruz, and Puebla in south-central Mexico. Although recent studies have demonstrated extensive geographic variation in venom composition and cryptic genetic diversity in this species, no modern studies have focused on geographic variation in morphology. Here we analyzed a series of qualitative, meristic, and morphometric traits from 347 specimens of C. scutulatus and show that this species is phenotypically cohesive without discrete subgroups, and that morphology follows a continuous cline in primarily color pattern and meristic traits across the major axis of its expansive distribution. Interpreted in the context of previously published molecular evidence, our morphological analyses suggest that multiple episodes of isolation and secondary contact among metapopulations during the Pleistocene were sufficient to produce distinctive genetic populations, which have since experienced gene flow to produce clinal variation in phenotypes without discrete or diagnosable distinctions among these original populations. For taxonomic purposes, we recommend that C. scutulatus be retained as a single species, although it is possible that C. s. salvini, which is morphologically the most distinctive population, could represent a peripheral isolate in the initial stages of speciation.

Several snake species that feed infrequently in nature have evolved the ability to massively upregulate intestinal form and function with each meal. While fasting, these snakes downregulate intestinal form and function, and upon feeding restore intestinal structure and function through major increases in cell growth and proliferation, metabolism and upregulation of digestive function. Previous studies have identified changes in gene expression that underlie this regenerative growth of the python intestine, but the unique features that differentiate this extreme regenerative growth from non-regenerative post-feeding responses exhibited by snakes that feed more frequently remain unclear. Here, we leveraged variation in regenerative capacity across three snake species-two distantly related lineages ( Crotalus and Python) that experience regenerative growth, and one ( Nerodia) that does not-to infer molecular mechanisms underlying intestinal regeneration using transcriptomic and proteomic approaches. Using a comparative approach, we identify a suite of growth, stress response and DNA damage response signalling pathways with inferred activity specifically in regenerating species, and propose a hypothesis model of interactivity between these pathways that may drive regenerative intestinal growth in snakes.

Most colubrid snake venoms have been poorly studied, despite the fact that they represent a great resource for biological, ecological, toxinological and pharmacological research. Herein, we explore the venom delivery system of the Aesculapian False Coral Snake Erythrolamprus aesculapii as well as some biochemical and toxicological properties of its venom. Its Duvernoy's venom gland is composed of serous secretory cells arranged in densely packed secretory tubules, and the most striking feature of its fang is their double-curved shape, exhibiting a beveled bladelike appearance near the tips. Although E. aesculapii resembles elapid snakes of the genus Micrurus in color pattern, this species produces a venom reminiscent of viperid venoms, containing mainly tissue-damaging toxins such as proteinases. Prominent hemorrhage developed both locally and systemically in mice injected with the venom, and the minimum hemorrhagic dose was found to be 18.8 μg/mouse; the lethal dose, determined in mice, was 9.5 ± 3.7 μg/g body weight. This work has toxicological implications that bites to humans by E. aesculapii could result in moderately severe local (and perhaps systemic) hemorrhage and gives insight into future directions for research on the venom of this species., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Many species of snakes produce venom as a chemical means of procuring potentially fractious prey. Studies have increasingly focused on venom compositional variation between and within individual snakes of the same species/subspecies, with significant differences often being observed. This variation in composition has been attributed to differences in age, season, diet, and environment, suggesting that these factors could help explain the inter- and intra-specific variation found in some snake venoms, perhaps via some type of feedback mechanism(s). To address several of these possible sources of variation, this study utilized wild-caught Western Diamondback Rattlesnakes (Crotalus atrox) from Cochise Co., AZ. Sixteen adult C. atrox were maintained in the lab on a diet of NSA mice for eight months to determine whether venom composition changed in captivity under a static diet in a stable environment. Reducing 1-D SDS-PAGE, fibrinogen degradation assays, reversed-phase HPLC, and MALDI-TOF mass spectrometry revealed only minor differences over time within individuals. Venom L-amino acid oxidase (LAAO) and phosphodiesterase activities significantly increased over the course of captivity, with no changes occurring in azocasein metalloproteinase, kallikrein-like serine proteinase (KLSP), or thrombin-like serine proteinase (TLSP) activities. Snake total length was positively correlated with TLSP activity and negatively correlated with LAAO and KLSP activity. There was typically a much higher degree of variation between individuals than within individuals for all analyses performed and measurements collected. Because the overall "fingerprint" of each snake's venom remained more/less constant, it is concluded that biologically significant changes in venom composition did not occur within individual C. atrox as a function of captivity/diet. However, this study does indicate that differences in activity levels do occur in minor venom enzyme components, but the differences observed are likely to be of minimal significance to the production of antivenom or to subsequent treatment of human envenomations., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

The genus Trimeresurus comprises a group of venomous pitvipers endemic to Southeast Asia and the Pacific Islands. Of these, Trimeresurus insularis, the White-lipped Island Pitviper, is a nocturnal, arboreal species that occurs on nearly every major island of the Lesser Sunda archipelago. In the current study, venom phenotypic characteristics of T. insularis sampled from eight Lesser Sunda Islands (Flores, Lembata, Lombok, Pantar, Sumba, Sumbawa, Timor, and Wetar) were evaluated via SDS-PAGE, enzymatic activity assays, fibrinogenolytic assays, gelatin zymography, and RP-HPLC, and the Sumbawa sample was characterized by venomic analysis. For additional comparative analyses, venoms were also examined from several species in the Trimeresurus complex, including T. borneensis, T. gramineus, T. puniceus, T. purpureomaculatus, T. stejnegeri, and Protobothrops flavoviridis. Despite the geographical isolation, T. insularis venoms from all eight islands demonstrated remarkable similarities in gel electrophoretic profiles and RP-HPLC patterns, and all populations had protein bands in the mass ranges of phosphodiesterases (PDE), l-amino acid oxidases (LAAO), P-III snake venom metalloproteinases (SVMP), serine proteases, cysteine-rich secretory proteins (CRISP), phospholipases A 2 (PLA 2 ), and C-type lectins. An exception was observed in the Lombok sample, which lacked protein bands in the mass range of serine protease and CRISP. Venomic analysis of the Sumbawa venom also identified these protein families, in addition to several proteins of lesser abundance (<1%), including glutaminyl cyclase, aminopeptidase, PLA 2 inhibitor, phospholipase B, cobra venom factor, 5'-nucleotidase, vascular endothelial growth factor, and hyaluronidase. All T. insularis venoms exhibited similarities in thrombin-like and PDE activities, while significant differences were observed for LAAO, SVMP, and kallikrein-like activities, though these differences were only observed for a few islands. Slight but noticeable differences were also observed with fibrinogen and gelatin digestion activities. Trimeresurus insularis venoms exhibited overall similarity to the other Trimeresurus complex species examined, with the exception of P. flavoviridis venom, which showed the greatest overall differentiation. Western blot analysis revealed that all major T. insularis venom proteins were recognized by Green Pitviper ( T. albolabris) antivenom, and reactivity was also seen with most venom proteins of the other Trimeresurus species, but incomplete antivenom-venom recognition was observed against P. flavoviridis venom proteins. These results demonstrate significant conservation in the venom composition of T. insularis across the Lesser Sunda archipelago relative to the other Trimeresurus species examined.

Here we use a chromosome-level genome assembly of a prairie rattlesnake ( Crotalus viridis ), together with Hi-C, RNA-seq, and whole-genome resequencing data, to study key features of genome biology and evolution in reptiles. We identify the rattlesnake Z Chromosome, including the recombining pseudoautosomal region, and find evidence for partial dosage compensation driven by an evolutionary accumulation of a female-biased up-regulation mechanism. Comparative analyses with other amniotes provide new insight into the origins, structure, and function of reptile microchromosomes, which we demonstrate have markedly different structure and function compared to macrochromosomes. Snake microchromosomes are also enriched for venom genes, which we show have evolved through multiple tandem duplication events in multiple gene families. By overlaying chromatin structure information and gene expression data, we find evidence for venom gene-specific chromatin contact domains and identify how chromatin structure guides precise expression of multiple venom gene families. Further, we find evidence for venom gland-specific transcription factor activity and characterize a complement of mechanisms underlying venom production and regulation. Our findings reveal novel and fundamental features of reptile genome biology, provide insight into the regulation of snake venom, and broadly highlight the biological insight enabled by chromosome-level genome assemblies., (© 2019 Schield et al.; Published by Cold Spring Harbor Laboratory Press.)

Context: The vast majority of the 2.5 million annual worldwide venomous snakebites are attributed to Viperidae or Elapidae envenomations. Of the nearly 2000 Colubridae species described, only a handful are known to cause medically significant envenomations. Considered medically insignificant, Heterodon nasicus (Western Hognose Snake) is a North American rear-fanged colubrid common in the legal pet trading industry. Previously reported cases of envenomations describe local pain, swelling, edema, and blistering. However, there are no reported cases of systemic or hematologic toxicity., Case Details: A 20-year-old female sustained a bite while feeding a captive H. nasicus causing local symptoms and thrombocytopenia. On day three after envenomation, the patient was seen in the emergency department for persistent pain, swelling, and blistering. At that time, she was found to have a platelet count of 90 × 10 9 /L. Previous routine platelet counts ranged from 315 to 373 × 10 9 /L during the prior two years. Local symptoms peaked on day seven post envenomation. Her local symptoms and thrombocytopenia improved on evaluation four months after envenomation., Discussion: We report the first Heterodon nasicus envenomation causing both local toxicity and thrombocytopenia. Potential mechanisms based on H. nasicus venom composition are discussed in detail. Treatment is largely supportive. Bites by H. nascius should be evaluated by a toxicologist familiar with Colubridae species. This represents the first reported case of hematologic toxicity from envenomation by a North American colubrid snake., (Copyright © 2018. Published by Elsevier Ltd.)

Snake venoms represent an enriched system for investigating the evolutionary processes that lead to complex and dynamic trophic adaptations. It has long been hypothesized that natural selection may drive geographic variation in venom composition, yet previous studies have lacked the population genetic context to examine these patterns. We leverage range-wide sampling of Mojave Rattlesnakes (Crotalus scutulatus) and use a combination of venom, morphological, phylogenetic, population genetic, and environmental data to characterize the striking dichotomy of neurotoxic (Type A) and hemorrhagic (Type B) venoms throughout the range of this species. We find that three of the four previously identified major lineages within C. scutulatus possess a combination of Type A, Type B, and a 'mixed' Type A + B venom phenotypes, and that fixation of the two main venom phenotypes occurs on a more fine geographic scale than previously appreciated. We also find that Type A + B individuals occur in regions of inferred introgression, and that this mixed phenotype is comparatively rare. Our results support strong directional local selection leading to fixation of alternative venom phenotypes on a fine geographic scale, and are inconsistent with balancing selection to maintain both phenotypes within a single population. Our comparisons to biotic and abiotic factors further indicate that venom phenotype correlates with fang morphology and climatic variables. We hypothesize that links to fang morphology may be indicative of co-evolution of venom and other trophic adaptations, and that climatic variables may be linked to prey distributions and/or physiology, which in turn impose selection pressures on snake venoms.

DNA barcoding is a simple technique used to develop a large-scale system of classification that is broadly applicable across a wide variety of taxa. DNA-based analysis of snake venoms can provide a system of classification independent of currently accepted taxonomic relationships by generating DNA barcodes specific to each venom sample. DNA purification from dried snake venoms has previously required large amounts of starting material, has resulted in low yields and inconsistent amplification, and was possible with front-fanged snakes only. Here, we present a modified DNA extraction protocol applied to venoms of both front- and rear-fanged snakes that requires significantly less starting material (1 mg) and yields sufficient amounts of DNA for successful PCR amplification of regions commonly used for DNA barcoding. [Formula: see text].

The Mojave rattlesnake (Crotalus scutulatus) inhabits deserts and arid grasslands of the western United States and Mexico. Despite considerable interest in its highly toxic venom and the recognition of two subspecies, no molecular studies have characterized range-wide genetic diversity and population structure or tested species limits within C. scutulatus. We used mitochondrial DNA and thousands of nuclear loci from double-digest restriction site associated DNA sequencing to infer population genetic structure throughout the range of C. scutulatus, and to evaluate divergence times and gene flow between populations. We find strong support for several divergent mitochondrial and nuclear clades of C. scutulatus, including splits coincident with two major phylogeographic barriers: the Continental Divide and the elevational increase associated with the Central Mexican Plateau. We apply Bayesian clustering, phylogenetic inference, and coalescent-based species delimitation to our nuclear genetic data to test hypotheses of population structure. We also performed demographic analyses to test hypotheses relating to population divergence and gene flow. Collectively, our results support the existence of four distinct lineages within C. scutulatus, and genetically defined populations do not correspond with currently recognized subspecies ranges. Finally, we use approximate Bayesian computation to test hypotheses of divergence among multiple rattlesnake species groups distributed across the Continental Divide, and find evidence for co-divergence at this boundary during the mid-Pleistocene., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Introduction: The Russell's Viper (RV) (Daboia russelii), a category I medically important snake, is responsible for a significant level of morbidity and mortality in the Indian sub-continent. Areas covered: The current review highlights the variation in RV venom (RVV) composition from different geographical locales on the Indian sub-continent, as revealed by biochemical and proteomic analyses. A comparison of these RVV proteomes revealed significant differences in the number of toxin isoforms and relative toxin abundances, highlighting the impact of geographic location on RVV composition. Antivenom efficacy studies have shown differential neutralization of toxicity and enzymatic activity of different RVV samples from the Indian sub-continent by commercial polyvalent antivenom (PAV). The proteome analysis has provided deeper insights into the variation of RVV composition leading to differences in antivenom efficacy and severity of clinical manifestations post RV-envenomation across the Indian sub-continent. Expert commentary: Variation in RVV antigenicity due to geographical differences and poor recognition of low molecular mass (<20 kDa) RVV toxins by PAV are serious concerns for effective antivenom treatment against RV envenomation. Improvements in immunization protocols that take into account the poorly immunogenic components and geographic variation in RVV composition, can lead to better hospital management of RV bite patients.

High-throughput technologies were used to identify venom gland toxin expression and to characterize the venom proteomes of two rear-fanged snakes, Ahaetulla prasina (Asian Green Vine Snake) and Borikenophis portoricensis (Puerto Rican Racer). Sixty-nine complete toxin-coding transcripts from 12 venom protein superfamilies (A. prasina) and 50 complete coding transcripts from 11 venom protein superfamilies (B. portoricensis) were identified in the venom glands. However, only 18% (A. prasina) and 32% (B. portoricensis) of the translated protein isoforms were detected in the proteome of these venoms. Both venom gland transcriptomes and venom proteomes were dominated by P-III metalloproteinases. Three-finger toxins, cysteine-rich secretory proteins, and C-type lectins were present in moderate amounts, but other protein superfamilies showed very low abundances. Venoms contained metalloproteinase activity comparable to viperid snake venom levels, but other common venom enzymes were absent or present at negligible levels. Western blot analysis showed metalloproteinase and cysteine-rich secretory protein epitopes shared with the highly venomous Boomslang (Dispholidus typus). The abundance of metalloproteinases emphasizes the important trophic role of these toxins. Comprehensive, transcriptome-informed definition of proteomes and functional characterization of venom proteins in rear-fanged snake families help to elucidate toxin evolution and provide models for protein structure-function analyses., (Copyright © 2018 Elsevier B.V. All rights reserved.)

The use of -omics technologies allows for the characterization of snake venom composition at a fast rate and at high levels of detail. In the present study, we investigated the protein content of Red-headed Krait ( Bungarus flaviceps ) venom. This analysis revealed a high diversity of snake venom protein families, as evidenced by high-throughput mass spectrometric analysis. We found all six venom protein families previously reported in a transcriptome study of the venom gland of B. flaviceps , including phospholipases A₂ (PLA₂s), Kunitz-type serine proteinase inhibitors (KSPIs), three-finger toxins (3FTxs), cysteine-rich secretory proteins (CRISPs), snaclecs, and natriuretic peptides. A combined approach of automated database searches and de novo sequencing of tandem mass spectra, followed by sequence similarity searches, revealed the presence of 12 additional toxin families. De novo sequencing alone was able to identify 58 additional peptides, and this approach contributed significantly to the comprehensive description of the venom. Abundant protein families comprise 3FTxs (22.3%), KSPIs (19%), acetylcholinesterases (12.6%), PLA₂s (11.9%), venom endothelial growth factors (VEGFs, 8.4%), nucleotidases (4.3%), and C-type lectin-like proteins (snaclecs, 3.3%); an additional 11 toxin families are present at significantly lower concentrations, including complement depleting factors, a family not previously detected in Bungarus venoms. The utility of a multifaceted approach toward unraveling the proteome of snake venoms, employed here, allowed detection of even minor venom components. This more in-depth knowledge of the composition of B. flaviceps venom facilitates a better understanding of snake venom molecular evolution, in turn contributing to more effective treatment of krait bites.

Venom proteins evolve rapidly, and as a trophic adaptation are excellent models for predator-prey evolutionary studies. The key to a deeper understanding of venom evolution is an integrated approach, combining prey assays with analysis of venom gene expression and venom phenotype. Here, we use such an approach to study venom evolution in the Amazon puffing snake, Spilotes sulphureus , a generalist feeder. We identify two novel three-finger toxins: sulditoxin and sulmotoxin 1. These new toxins are not only two of the most abundant venom proteins, but are also functionally intriguing, displaying distinct prey-specific toxicities. Sulditoxin is highly toxic towards lizard prey, but is non-toxic towards mammalian prey, even at greater than 22-fold higher dosage. By contrast, sulmotoxin 1 exhibits the reverse trend. Furthermore, evolutionary analysis and structural modelling show highest sequence variability in the central loop of these proteins, probably driving taxon-specific toxicity. This is, to our knowledge, the first case in which a bimodal and contrasting pattern of toxicity has been shown for proteins in the venom of a single snake in relation to diet. Our study is an example of how toxin gene neofunctionalization can result in a venom system dominated by one protein superfamily and still exhibit flexibility in prey capture efficacy., (© 2018 The Author(s).)

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.