Search

Understanding the proximate and ultimate causes of phenotypic variation is fundamental in evolutionary research, as such variation provides the substrate for selection to act upon. Although trait variation can arise due to selection, the importance of neutral processes is sometimes understudied. We presented the first reference-quality genome of the Red Diamond Rattlesnake (Crotalus ruber) and used range-wide 'omic data to estimate the degree to which neutral and adaptive evolutionary processes shaped venom evolution. We characterized population structure and found substantial genetic differentiation across two populations, each with distinct demographic histories. We identified significant differentiation in venom expression across age classes with substantially reduced but discernible differentiation across populations. We then used conditional redundancy analysis to test whether venom expression variation was best predicted by neutral divergence patterns or geographically variable (a)biotic factors. Snake size was the most significant predictor of venom variation, with environment, prey availability, and neutral sequence variation also identified as significant factors, though to a lesser degree. By directly including neutrality in the model, our results confidently highlight the predominant, yet not singular, role of life history in shaping venom evolution., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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.)

Disintegrins, a family of snake venom protein, which are capable of modulating the activity of integrins that play a fundamental role in the regulation of many physiological and pathological processes. The main purpose of this study is to obtain the recombinant disintegrin (r-DI) and evaluate its biological activity. In this study, we explored a high-level expression prokaryotic system and purification strategy for r-DI. Then, r-DI was treated to assay effects on cell growth, migration, and invasion. The affinity for the interactions of r-DI with integrin was determined using Surface plasmon resonance (SPR) analyses. The r-DI can be expressed in Escherichia coli and purified by one-step chromatography. The r-DI can inhibit B16F10 cells proliferation, migration, and invasion. Also, we found that r-DI could interact with the integrin αIIbβ3 (GPIIb/IIIa). The r-DI can be expressed, purified, characterized through functional assays, and can also maintain strong biological activities. Thus, this study showed potential therapeutic effects of r-DI for further functional and structural studies., (© 2024. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake ( Crotalus adamanteus ), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus . Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly., Competing Interests: Competing interests statement:The authors declare no competing interest.

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.)

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis., (© 2023. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

RNA purification and cDNA synthesis represents the starting point for molecular analyses of snake venom proteins-enzymes. Usually, the sacrifice of snakes is necessary for venom gland extraction to identify protein-coding transcripts; however, the venom can be used as a source of transcripts. Although there are methods for obtaining RNA from venom, no comparative analysis has been conducted in the Bothrops genus. In the present study, we compared four commercial methods for RNA purification and cDNA synthesis from venom (liquid, lyophilized, or long-term storage) of four clinically relevant species of Peruvian Bothrops . Our results show that the TRIzol method presents the highest yield of RNA purified from venom (59 ± 11 ng/100 µL or 10 mg). The SuperScript First-Strand Synthesis System kit produced high amounts of cDNA (3.2 ± 1.2 ng cDNA/ng RNA), and the highest value was from combination with the Dynabeads mRNA DIRECT kit (4.8 ± 2.0 ng cDNA/ng RNA). The utility of cDNA was demonstrated with the amplification of six relevant toxins: thrombin-like enzymes, P-I and P-III metalloproteinases, acid and basic phospholipases A2, and disintegrins. To our knowledge, this is the first comparative study of RNA purification and cDNA synthesis methodologies from Bothrops genus venom.

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).)

Understanding the molecular basis of adaptations in coevolving species requires identifying the genes that underlie reciprocally selected phenotypes, such as those involved in venom in snakes and resistance to the venom in their prey. In this regard, California ground squirrels (CGS; Otospermophilus beecheyi) are eaten by northern Pacific rattlesnakes (Crotalus oreganus oreganus), but individual squirrels may still show substantial resistance to venom and survive bites. A recent study using proteomics identified venom interactive proteins (VIPs) in the blood serum of CGS. These VIPs represent possible resistance proteins, but the sequences of genes encoding them are unknown despite the value of such data to molecular studies of coevolution. To address this issue, we analyzed a de novo assembled transcriptome from CGS liver tissue-where many plasma proteins are synthesized-and other tissues from this species. We then examined VIP sequences in terms of three characteristics that identify them as possible resistance proteins: evidence for positive selection, high liver expression, and nonsynonymous variation across CGS populations. Based on these characteristics, we identified five VIPs (i.e., α-2-macroglobulin, α-1-antitrypsin-like protein GS55-LT, apolipoprotein A-II, hibernation-associated plasma protein HP-20, and hibernation-associated plasma protein HP-27) as the most likely candidates for resistance proteins among VIPs identified to date. Four of these proteins have been previously implicated in conferring resistance to the venom in mammals, validating our approach. When combined with the detailed information available for rattlesnake venom proteins, these results set the stage for future work focused on understanding coevolutionary interactions at the molecular level between these species., (© 2022 The Authors. Journal of Experimental Zoology Part B: Molecular and Developmental Evolution published by Wiley Periodicals LLC.)

The complexity of snake venom composition reflects adaptation to the diversity of prey and may be driven at times by a coevolutionary arms race between snakes and venom-resistant prey. However, many snakes are also resistant to their own venom due to serum-borne inhibitors of venom toxins, which raises the question of how snake autoinhibitors maintain their efficacy as venom proteins evolve. To investigate this potential three-way arms race among venom, prey, and autoinhibitors, we have identified and traced the evolutionary origin of serum inhibitors of snake venom metalloproteinases (SVMPs) in the Western Diamondback rattlesnake Crotalus atrox which possesses the largest known battery of SVMP genes among crotalids examined. We found that C. atrox expresses five members of a Fetuin A-related metalloproteinase inhibitor family but that one family member, FETUA-3, is the major SVMP inhibitor that binds to approximately 20 different C. atrox SVMPs and inhibits activities of all three SVMP classes. We show that the fetua-3 gene arose deep within crotalid evolution before the origin of New World species but, surprisingly, fetua-3 belongs to a different paralog group than previously identified SVMP inhibitors in Asian and South American crotalids. Conversely, the C. atrox FETUA-2 ortholog of previously characterized crotalid SVMP inhibitors shows limited activity against C. atrox SVMPs. These results reveal that there has been a functional evolutionary shift in the major SVMP inhibitor in the C. atrox lineage as the SVMP family expanded and diversified in the Crotalus lineage. This broad-spectrum inhibitor may be of potential therapeutic interest.

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.)

Opossums in the tribe Didelphini are resistant to pit viper venoms and are hypothesized to be coevolving with venomous snakes. Specifically, a protein involved in blood clotting (von Willebrand factor [vWF] which is targeted by snake venom C-type lectins [CTLs]) has been found to undergo rapid adaptive evolution in Didelphini. Several unique amino acid changes in vWF could explain their resistance; however, experimental evidence that these changes disrupt binding to venom CTLs was lacking. Furthermore, without explicit testing of ancestral phenotypes to reveal the mode of evolution, the assertion that this system represents an example of coevolution rather than noncoevolutionary adaptation remains unsupported. Using expressed vWF proteins and purified venom CTLs, we quantified binding affinity for vWF proteins from all resistant taxa, their venom-sensitive relatives, and their ancestors. We show that CTL-resistant vWF is present in opossums outside clade Didelphini and likely across a wider swath of opossums (family Didelphidae) than previously thought. Ancestral reconstruction and in vitro testing of vWF phenotypes in a clade of rapidly evolving opossums reveal a pattern consistent with trench warfare coevolution between opossums and their venomous snake prey., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Proteins isolated from natural sources can be composed of a mixture of isoforms with similar physicochemical properties that coexist in the final steps of purification. Yet, even where unverified, the assumed sequence is enforced throughout the structural studies. Herein, we propose a novel perspective to address the usually neglected sequence heterogeneity of natural products by integrating biophysical, genetic and structural data in our program SEQUENCE SLIDER. The aim is to assess the evidence supporting chemical composition in structure determination. Locally, we interrogate the experimental map to establish which side chains are supported by the structural data, and the genetic information relating sequence conservation is integrated into this statistic. Hence, we build a constrained peptide database, containing most probable sequences to interpret mass spectrometry data (MS). In parallel, we perform MS de novo sequencing with genomic-based algorithms to detect point mutations. We calibrated SLIDER with Gallus gallus lysozyme, whose sequence is unequivocally established and numerous natural isoforms are reported. We used SLIDER to characterize a metalloproteinase and a phospholipase A2-like protein from the venom of Bothrops moojeni and a crotoxin from Crotalus durissus collilineatus. This integrated approach offers a more realistic structural descriptor to characterize macromolecules isolated from natural sources., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Interspecific differences in snake venom compositions can result from distinct regulatory mechanisms acting in each species. However, comparative analyses focusing on identifying regulatory elements and patterns that led to distinct venom composition are still scarce. Among venomous snakes, Bothrops cotiara and Bothrops fonsecai represent ideal models to complement our understanding of the regulatory mechanisms of venom production. These recently diverged species share a similar specialized diet, habitat, and natural history, but each presents a distinct venom phenotype. Here, we integrated data from the venom gland transcriptome and miRNome and the venom proteome of B. fonsecai and B. cotiara to better understand the regulatory mechanisms that may be acting to produce differing venom compositions. We detected not only the presence of similar toxin isoforms in both species but also distinct expression profiles of phospholipases A2 (PLA2) and some snake venom metalloproteinases (SVMPs) and snake venom serine proteinases (SVSPs) isoforms. We found evidence of modular expression regulation of several toxin isoforms implicated in venom divergence and observed correlated expression of several transcription factors. We did not find strong evidence for miRNAs shaping interspecific divergence of the venom phenotypes, but we identified a subset of toxin isoforms whose final expression may be fine-tuned by specific miRNAs. Sequence analysis on orthologous toxins showed a high rate of substitutions between PLA2s, which indicates that these toxins may be under strong positive selection or represent paralogous toxins in these species. Our results support other recent studies in suggesting that gene regulation is a principal mode of venom evolution across recent timescales, especially among species with conserved ecotypes.

Ecologically divergent selection can lead to the evolution of reproductive isolation through the process of ecological speciation, but the balance of responsible evolutionary forces is often obscured by an inadequate assessment of demographic history and the genetics of traits under selection. Snake venoms have emerged as a system for studying the genetic basis of adaptation because of their genetic tractability and contributions to fitness, and speciation in venomous snakes can be associated with ecological diversification such as dietary shifts and corresponding venom changes. Here, we explored the neurotoxic (type A)-hemotoxic (type B) venom dichotomy and the potential for ecological speciation among Timber Rattlesnake ( Crotalus horridus ) populations. Previous work identified the genetic basis of this phenotypic difference, enabling us to characterize the roles geography, history, ecology, selection, and chance play in determining when and why new species emerge or are absorbed. We identified significant genetic, proteomic, morphological, and ecological/environmental differences at smaller spatial scales, suggestive of incipient ecological speciation between type A and type B C. horridus . Range-wide analyses, however, rejected the reciprocal monophyly of venom type, indicative of varying intensities of introgression and a lack of reproductive isolation across the range. Given that we have now established the phenotypic distributions and ecological niche models of type A and B populations, genome-wide data are needed and capable of determining whether type A and type B C. horridus represent distinct, reproductively isolated lineages due to incipient ecological speciation or differentiated populations within a single species.

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca , a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A 2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option., Competing Interests: The authors declare no competing interest.

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network., Competing Interests: The authors declare no competing interest.

Ontogenetic changes in venom composition have been described in Bothrops snakes, but only a few studies have attempted to identify the targeted paralogues or the molecular mechanisms involved in modifications of gene expression during ontogeny. In this study, we decoded B. jararacussu venom gland transcripts from six specimens of varying sizes and analyzed the variability in the composition of independent venom proteomes from 19 individuals. We identified 125 distinct putative toxin transcripts, and of these, 73 were detected in venom proteomes and only 10 were involved in the ontogenetic changes. Ontogenetic variability was linearly related to snake size and did not correspond to the maturation of the reproductive stage. Changes in the transcriptome were highly predictive of changes in the venom proteome. The basic myotoxic phospholipases A 2 (PLA 2 s) were the most abundant components in larger snakes, while in venoms from smaller snakes, PIII-class SVMPs were the major components. The snake venom metalloproteinases (SVMPs) identified corresponded to novel sequences and conferred higher pro-coagulant and hemorrhagic functions to the venom of small snakes. The mechanisms modulating venom variability are predominantly related to transcriptional events and may consist of an advantage of higher hematotoxicity and more efficient predatory function in the venom from small snakes.

Ontogenetic shifts in venom occur in many snakes but establishing their nature as gradual or discrete processes required additional study. We profiled shifts in venom expression from the neonate to adult sizes of two rattlesnake species, the eastern diamondback and the timber rattlesnake. We used serial sampling and venom chromatographic profiling to test if ontogenetic change occurs gradually or discretely. We found evidence for gradual shifts in overall venom composition in six of eight snakes, which sometimes spanned more than two years. Most chromatographic peaks shift gradually, but one quarter shift in a discrete fashion. Analysis of published diet data showed gradual shifts in overall diet composition across the range of body sizes attained by our eight study animals, while the shifts in abundance of different prey classes varied in form from gradual to discrete. Testosterone concentrations were correlated with the change in venom protein composition, but the relationship is not strong enough to suggest causation. Venom research employing simple juvenile versus adult size thresholds may be failing to account for continuous variation in venom composition lifespan. Our results imply that venom shifts represent adaptive matches to dietary shifts and highlight venom for studies of alternative gene regulatory mechanisms.

Understanding how interspecific interactions mould the molecular basis of adaptations in coevolving species is a long-sought goal of evolutionary biology. Venom in predators and venom resistance proteins in prey are coevolving molecular phenotypes, and while venoms are highly complex mixtures it is unclear if prey respond with equally complex resistance traits. Here, we use a novel molecular methodology based on protein affinity columns to capture and identify candidate blood serum resistance proteins ("venom interactive proteins" [VIPs]) in California Ground Squirrels (Otospermophilus beecheyi) that interact with venom proteins from their main predator, Northern Pacific Rattlesnakes (Crotalus o. oreganus). This assay showed that serum-based resistance is both population- and species-specific, with serum proteins from ground squirrels showing higher binding affinities for venom proteins of local snakes compared to allopatric individuals. Venom protein specificity assays identified numerous and diverse candidate prey resistance VIPs but also potential targets of venom in prey tissues. Many specific VIPs bind to multiple snake venom proteins and, conversely, single venom proteins bind multiple VIPs, demonstrating that a portion of the squirrel blood serum "resistome" involves broad-based inhibition of nonself proteins and suggests that resistance involves a toxin scavenging mechanism. Analyses of rates of evolution of VIP protein homologues in related mammals show that most of these proteins evolve under purifying selection possibly due to molecular constraints that limit the evolutionary responses of prey to rapidly evolving snake venom proteins. Our method represents a general approach to identify specific proteins involved in co-evolutionary interactions between species at the molecular level., (© 2020 John Wiley & Sons Ltd.)

The genetic origins of novelty are a central interest of evolutionary biology. Most new proteins evolve from preexisting proteins but the evolutionary path from ancestral gene to novel protein is challenging to trace, and therefore the requirements for and order of coding sequence changes, expression changes, or gene duplication are not clear. Snake venoms are important novel traits that are comprised of toxins derived from several distinct protein families, but the genomic and evolutionary origins of most venom components are not understood. Here, we have traced the origin and diversification of one prominent family, the snake venom metalloproteinases (SVMPs) that play key roles in subduing prey in many vipers. Genomic analyses of several rattlesnake ( Crotalus ) species revealed the SVMP family massively expanded from a single, deeply conserved adam28 disintegrin and metalloproteinase gene, to as many as 31 tandem genes in the Western Diamondback rattlesnake ( Crotalus atrox ) through a number of single gene and multigene duplication events. Furthermore, we identified a series of stepwise intragenic deletions that occurred at different times in the course of gene family expansion and gave rise to the three major classes of secreted SVMP toxins by sequential removal of a membrane-tethering domain, the cysteine-rich domain, and a disintegrin domain, respectively. Finally, we show that gene deletion has further shaped the SVMP complex within rattlesnakes, creating both fusion genes and substantially reduced gene complexes. These results indicate that gene duplication and intragenic deletion played essential roles in the origin and diversification of these novel biochemical weapons., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

Small serum proteins (SSPs) are low-molecular-weight proteins in snake serum with affinities for various venom proteins. Five SSPs, Pf SSP-1 through Pf SSP-5, have been reported in Protobothrops flavoviridis ("habu", Pf ) serum so far. Recently, we reported that the five genes encoding these Pf SSPs are arranged in tandem on a single chromosome. However, the physiological functions and evolutionary origins of the five SSPs remain poorly understood. In a detailed analysis of the habu draft genome, we found a gene encoding a novel SSP, SSP-6. Structural analysis of the genes encoding SSPs and their genomic arrangement revealed the following: (1) SSP-6 forms a third SSP subgroup; (2) SSP-5 and SSP-6 were present in all snake genomes before the divergence of non-venomous and venomous snakes, while SSP-4 was acquired only by venomous snakes; (3) the composition of paralogous SSP genes in snake genomes seems to reflect snake habitat differences; and (4) the evolutionary emergence of SSP genes is probably related to the physiological functions of SSPs, with an initial snake repertoire of SSP-6 and SSP-5. SSP-4 and its derivative, SSP-3, as well as SSP-1 and SSP-2, appear to be venom-related and were acquired later.

Background: Modularity is the tendency for systems to organize into semi-independent units and can be a key to the evolution and diversification of complex biological systems. Snake venoms are highly variable modular systems that exhibit extreme diversification even across very short time scales. One well-studied venom phenotype dichotomy is a trade-off between neurotoxicity versus hemotoxicity that occurs through the high expression of a heterodimeric neurotoxic phospholipase A 2 (PLA 2 ) or snake venom metalloproteinases (SVMPs). We tested whether the variation in these venom phenotypes could occur via variation in regulatory sub-modules through comparative venom gland transcriptomics of representative Black-Speckled Palm-Pitvipers (Bothriechis nigroviridis) and Talamancan Palm-Pitvipers (B. nubestris)., Results: We assembled 1517 coding sequences, including 43 toxins for B. nigroviridis and 1787 coding sequences including 42 toxins for B. nubestris. The venom gland transcriptomes were extremely divergent between these two species with one B. nigroviridis exhibiting a primarily neurotoxic pattern of expression, both B. nubestris expressing primarily hemorrhagic toxins, and a second B. nigroviridis exhibiting a mixed expression phenotype. Weighted gene coexpression analyses identified six submodules of transcript expression variation, one of which was highly associated with SVMPs and a second which contained both subunits of the neurotoxic PLA 2 complex. The sub-module association of these toxins suggest common regulatory pathways underlie the variation in their expression and is consistent with known patterns of inheritance of similar haplotypes in other species. We also find evidence that module associated toxin families show fewer gene duplications and transcript losses between species, but module association did not appear to affect sequence diversification., Conclusion: Sub-modular regulation of expression likely contributes to the diversification of venom phenotypes within and among species and underscores the role of modularity in facilitating rapid evolution of complex traits.

(1) Background. Snake venom phosphodiesterases (SVPDEs) are among the least studied venom enzymes. In envenomation, they display various pathological effects, including induction of hypotension, inhibition of platelet aggregation, edema, and paralysis. Until now, there have been no 3D structural studies of these enzymes, thereby preventing structure-function analysis. To enable such investigations, the present work describes the model-based structural and functional characterization of a phosphodiesterase from Crotalus adamanteus venom, named PDE_ Ca . (2) Methods. The PDE_ Ca structure model was produced and validated using various software (model building: I-TESSER, MODELLER 9v19, Swiss-Model, and validation tools: PROCHECK, ERRAT, Molecular Dynamic Simulation, and Verif3D). (3) Results. The proposed model of the enzyme indicates that the 3D structure of PDE_ Ca comprises four domains, a somatomedin B domain, a somatomedin B-like domain, an ectonucleotide pyrophosphatase domain, and a DNA/RNA non-specific domain. Sequence and structural analyses suggest that differences in length and composition among homologous snake venom sequences may account for their differences in substrate specificity. Other properties that may influence substrate specificity are the average volume and depth of the active site cavity. (4) Conclusion. Sequence comparisons indicate that SVPDEs exhibit high sequence identity but comparatively low identity with mammalian and bacterial PDEs., Competing Interests: The authors declare no conflict of interest.

Snake venoms are complex mixtures of toxic proteins encoded by various gene families that function synergistically to incapacitate prey. A huge repertoire of snake venom genes and proteins have been reported, and alternative splicing is suggested to be involved in the production of divergent gene transcripts. However, a genome-wide survey of the transcript repertoire and the extent of alternative splicing still remains to be determined. In this study, the comprehensive analysis of transcriptomes in the venom gland was achieved by using PacBio sequencing. Extensive alternative splicing was observed in three venom protein gene families, metalloproteinase (MP), serine protease (SP), and vascular endothelial growth factors (VEGF). Eleven MP and SP genes and a VEGF gene are expressed as a total of 81, 61, and 8 transcript variants, respectively. In the MP gene family, individual genes are transcribed into different classes of MPs by alternative splicing. We also observed trans-splicing among the clustered SP genes. No other venom genes as well as non-venom counterpart genes exhibited alternative splicing. Our results thus indicate a potential contribution of mRNA alternative and trans-splicing in the production of highly variable transcripts of venom genes in the habu snake., Competing Interests: The authors declare no conflict of interest.

Traits can evolve rapidly through changes in gene expression or protein-coding sequences. However, these forms of genetic variation can be correlated and changes to one can influence the other. As a result, we might expect traits lacking differential expression to preferentially evolve through changes in protein sequences or morphological adaptation. Given the lack of differential expression across the distribution of sidewinder rattlesnakes ( Crotalus cerastes), we tested this hypothesis by comparing the coding regions of genes expressed in the venom gland transcriptomes and fang morphology. We calculated Tajima's D and F ST across four populations comparing toxin and nontoxin loci. Overall, we found little evidence of directional selection or differentiation between populations, suggesting that changes to protein sequences do not underlie the evolution of sidewinder venom or that toxins are under extremely variant selection pressures. Although low-expression toxins do not have higher sequence divergence between populations, they do have more standing variation on which selection can act. Additionally, we found significant differences in fang length among populations. The lack of differential expression and sequence divergence suggests sidewinders-given their generalist diet, moderate gene flow and environmental variation-are under stabilizing selection which functions to maintain a generalist phenotype. Overall, we demonstrate the importance of examining the relationship between gene expression and protein-coding changes to understand the evolution of complex traits.

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.)

Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.

The migration-selection interaction is the strongest determinant of whether a beneficial allele increases in frequency within a population. Results of empirical studies examining the role of gene flow in an adaptive context, however, have largely been equivocal, with examples of opposing outcomes being repeatedly documented (e.g., local adaptation with high levels of gene flow vs. gene swamping). We compared neutral genomic and venom expression divergence for three sympatric pit vipers with differing ecologies to determine if and how migration-selection disequilibria result in local adaptation. We specifically tested whether neutral differentiation predicted phenotypic differentiation within an isolation-by-distance framework. The decoupling of neutral and phenotypic differentiation would indicate selection led to adaptive divergence irrespective of migration, whereas a significant relationship between neutral and venom expression differentiation would provide evidence in favor of the constraining force of gene flow. Neutral differentiation and geographic distance predicted phenotypic differentiation only in the generalist species, indicating that selection was the predominant mechanism in the migration-selection balance underlying adaptive venom evolution in both specialists. Dispersal is thought to be a stronger influence on genetic differentiation than specialization, but our results suggest the opposite. Greater specialization may lead to greater diversification rates, and extreme spatial and temporal variation in selective pressures can favor generalist phenotypes evolving under strong stabilizing selection. Our results are consistent with these expectations, suggesting that the equivocal findings of studies examining the role of gene flow in an adaptive context may be explained by ecological specialization theory.

A mRNA transcript that codes for a phospholipase (PLA 2 ) was isolated from a single venom gland of the Bothrops ammodytoides viper. The PLA 2 transcript was cloned onto a pCR ® 2.1-TOPO vector and subsequently expressed heterologously in the E. coli strain M15, using the pQE30 vector. The recombinant phospholipase was named rBamPLA2_1, and is composed of an N-terminal fusion protein of 16 residues, along with 122 residues from the mature protein that includes 14 cysteines that form 7 disulfide bonds. Following bacterial expression, rBamPLA2_1 was obtained from inclusion bodies and extracted using a chaotropic agent. rBamPLA2_1 had an experimental molecular mass of 15,692.5 Da that concurred with its theoretical molecular mass. rBamPLA2_1 was refolded in in vitro conditions and after refolding, three main protein fractions with similar molecular masses, were identified. Although, the three fractions were considered to represent different oxidized cystine isoforms, their secondary structures were comparable. All three recombinant isoforms were active on egg-yolk phospholipid and recognized similar cell membrane phospholipids to be native PLA 2 s, isolated from B. ammodytoides venom. A mixture of the three rBamPLA2_1 cystine isoforms was used to immunize a horse in order to produce serum antibodies (anti-rBamPLA2_1), which partially inhibited the indirect hemolytic activity of B. ammodytoides venom. Although, anti-rBamPLA2_1 antibodies were not able to recognize crotoxin, a PLA 2 from the venom of a related but different viper genus, Crotalus durissus terrificus, they recognized PLA 2 s in other venoms from regional species of Bothrops., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Herein we evaluated the genotoxic effects of BnSP-6, a Lys-49 phospholipase A 2 (PLA 2 ) from Bothrops pauloensis, on breast cancer cells. BnSP-6 was able to induce a higher cytotoxic and genotoxic activity in MDA-MB-231 cells, when compared to MCF10A (a non-tumorigenic breast cell line), suggesting that this protein presented a possible preference for cancer cells. BnSP-6 inhibited MDA-MB-231 proliferation at 24, 48 and 72 h. In addition, BnSP-6 induced significant increase in the percentage of TUNEL-positive cells, a marker of DNA damage. To obtain novel insight into the direct DNA damage interference in MDA-MB-231 survival and proliferation, we evaluated cell cycle progression. BnSP-6 produced a significant decrease in 2N (G1) and an increase in the G2/M phase and this capacity is likely related to the modulation of expression of progression cell cycle-associated genes (CCND1, CCNE1, CDC25A, CHEK2, E2F1, CDH-1 and NF-kB). Taken together, these results indicate that BnSP-6 induces DNA damage in breast cancer cells and is an attractive model for developing innovative therapeutic agents against breast cancer., (Copyright © 2018 Elsevier B.V. All rights reserved.)

The first toxin isolated from the venomous pit viper Porthidium ophryomegas is a basic pentameric phospholipase A 2 (PophPLA 2 ). Elucidation of its amino acid sequence showed that it belongs to the group IIA of secreted PLA 2 s, with the presence of all 14 conserved cysteine positions. The toxin displayed catalytic activity, in agreement with the presence of Asp49 in its sequence of 121 residues. SDS-PAGE analysis revealed that this toxin is pentameric in non-reducing conditions, a structural organization that has not been described for any viperid PLA 2 . PophPLA 2 displayed moderate myotoxic (in vivo) and cytotoxic (in vitro) activities, as well as anticoagulant activity on human plasma (in vitro). PophPLA 2 was not lethal, and did not induce signs of toxicity or distress in mice, when administered intravenously at a dose of up to 100 μg (5.9 μg/g). The toxin showed highest sequence identity with other PLA 2 s from the venoms of ancestral Asian pit viper species., (Copyright © 2018. Published by Elsevier B.V.)

In order to study how acidic pro-peptides inhibit the antimicrobial activity of antimicrobial peptides, we introduce a simple model system, consisting of a 19 amino-acid long antimicrobial peptide, and an N-terminally attached, 10 amino-acid long acidic model pro-peptide. The antimicrobial peptide is a fragment of the crotalicidin peptide, a member of the cathelidin family, from rattlesnake venom. The model pro-peptide is a deca (glutamic acid). Attachment of the model pro-peptide only leads to a moderately large reduction in the binding to- and induced leakage of model liposomes, while the antimicrobial activity of the crotalicidin fragment is completely inhibited by attaching the model pro-peptide. Attaching the pro-peptide induces a conformational change to a more helical conformation, while there are no signs of intra- or intermolecular peptide complexation. We conclude that inhibition of antimicrobial activity by the model pro-peptide might be related to a conformational change induced by the pro-peptide domain, and that additional effects beyond induced changes in membrane activity must also be involved.

Venom-gland transcriptomics is a key tool in the study of the evolution, ecology, function, and pharmacology of animal venoms. In particular, gene-expression variation and coding sequences gained through transcriptomics provide key information for explaining functional venom variation over both ecological and evolutionary timescales. The accuracy and usefulness of inferences made through transcriptomics, however, is limited by the accuracy of the transcriptome assembly, which is a bioinformatic problem with several possible solutions. Several methods have been employed to assemble venom-gland transcriptomes, with the Trinity assembler being the most commonly applied among them. Although previous evidence of variation in performance among assembly software exists, particularly regarding recovery of difficult-to-assemble multigene families such as snake venom metalloproteinases, much work to date still employs a single assembly method. We evaluated the performance of several commonly used de novo assembly methods for the recovery of both nontoxin transcripts and complete, high-quality venom-gene transcripts across eleven snake and four scorpion transcriptomes. We varied k -mer sizes used by some assemblers to evaluate the impact of k -mer length on transcript recovery. We showed that the recovery of nontoxin transcripts and toxin transcripts is best accomplished through different assembly software, with SDT at smaller k -mer lengths and Trinity being best for nontoxin recovery and a combination of SeqMan NGen and a seed-and-extend approach implemented in Extender as the best means of recovering a complete set of toxin transcripts. In particular, Extender was the only means tested capable of assembling multiple isoforms of the diverse snake venom metalloproteinase family, while traditional approaches such as Trinity recovered at most one metalloproteinase transcript. Our work demonstrated that traditional metrics of assembly performance are not predictive of performance in the recovery of complete and high quality toxin genes. Instead, effective venom-gland transcriptomic studies should combine and quality-filter the results of several assemblers with varying algorithmic strategies.

Snake venoms are a mixture of different molecules that can be used in the design of drugs for various diseases. The study of these venoms has relied on strategies that use complete venom extracted from animals in captivity or from venom glands that require the sacrifice of the animals. Colombia, a country with political and geographical conflicts has difficult access to certain regions. A strategy that can prevent the sacrifice of animals and could allow the study of samples collected in the field is necessary. We report the use of lyophilized venom from Crotalus durissus cumanensis as a model to test, for the first time, a protocol for the amplification of complete toxins from Colombian venom samples collected in the field. In this protocol, primers were designed from conserved region from Crotalus sp. mRNA and EST regions to maximize the likelihood of coding sequence amplification. We obtained the sequences of Metalloproteinases II, Disintegrins, Disintegrin-Like, Phospholipases A 2, C-type Lectins and Serine proteinases from Crotalus durissus cumanensis and compared them to different Crotalus sp sequences available on databases obtaining concordance between the toxins amplified and those reported. Our strategy allows the use of lyophilized venom to obtain complete toxin sequences from samples collected in the field and the study of poorly characterized venoms in challenging environments., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Natural selection is generally expected to favor one form of a given trait within a population. The presence of multiple functional variants of traits involved in activities such as feeding, reproduction, or the defense against predators is relatively uncommon within animal species. The genetic architecture and evolutionary mechanisms underlying the origin and maintenance of such polymorphisms are of special interest. Among rattlesnakes, several instances of the production of biochemically distinct neurotoxic or hemorrhagic venom types within the same species are known. Here, we investigated the genetic basis of this phenomenon in three species and found that neurotoxic and hemorrhagic individuals of the same species possess markedly different haplotypes at two toxin gene complexes. For example, neurotoxic and hemorrhagic Crotalus scutulatus individuals differ by 5 genes at the phospholipase A2 (PLA2) toxin gene complex and by 11 genes at the metalloproteinase (MP) gene complex. A similar set of extremely divergent haplotypes also underlies alternate venom types within C. helleri and C. horridus. We further show that the MP and PLA2 haplotypes of neurotoxic C. helleri appear to have been acquired through hybridization with C. scutulatus-a rare example of the horizontal transfer of a potentially highly adaptive suite of genes. These large structural variants appear analogous to immunity gene complexes in host-pathogen arms races and may reflect the impact of balancing selection at the PLA2 and MP complexes for predation on different prey., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Evolutionary innovations and complex phenotypes seemingly require an improbable amount of genetic change to evolve. Rattlesnakes display two dramatically different venom phenotypes. Type I venoms are hemorrhagic with low systemic toxicity and high expression of tissue-destroying snake venom metalloproteinases. Type II venoms are highly neurotoxic and lack snake venom metalloproteinase expression and associated hemorrhagic activity. This dichotomy hinges on Mojave toxin (MTx), a phospholipase A2 (PLA2) based β-neurotoxin expressed in Type II venoms. MTx is comprised of a nontoxic acidic subunit that undergoes extensive proteolytic processing and allosterically regulates activity of a neurotoxic basic subunit. Evolution of the acidic subunit presents an evolutionary challenge because the need for high expression of a nontoxic venom component and the proteolytic machinery required for processing suggests genetic changes of seemingly little immediate benefit to fitness. We showed that MTx evolved through a cascading series of exaptations unlocked by a single nucleotide change. The evolution of one new cleavage site in the acidic subunit unmasked buried cleavage sites already present in ancestral PLA2s, enabling proteolytic processing. Snake venom serine proteases, already present in the venom to disrupt prey hemostasis, possess the requisite specificities for MTx acidic subunit proteolysis. The dimerization interface between MTx subunits evolved by exploiting a latent, but masked, hydrophobic interaction between ancestral PLA2s. The evolution of MTx through exaptation of existing functional and structural features suggests complex phenotypes that depend on evolutionary innovations can arise from minimal genetic change enabled by prior evolution.

Rattlesnake venoms may be classified according to the presence/absence and relative abundance of the neurotoxic phospholipases A 2 s (PLA 2 s), such as Mojave toxin, and snake venom metalloproteinases (SVMPs). In Mexico, studies to determine venom variation in Mojave Rattlesnakes ( Crotalus scutulatus scutulatus ) are limited and little is known about the biological and proteolytic activities in this species. Tissue (34) and venom (29) samples were obtained from C. s. scutulatus from different locations within their distribution in Mexico. Mojave toxin detection was carried out at the genomic (by PCR) and protein (by ELISA) levels for all tissue and venom samples. Biological activity was tested on representative venoms by measuring LD 50 and hemorrhagic activity. To determine the approximate amount of SVMPs, 15 venoms were separated by RP-HPLC and variation in protein profile and proteolytic activity was evaluated by SDS-PAGE ( n = 28) and Hide Powder Azure proteolytic analysis ( n = 27). Three types of venom were identified in Mexico which is comparable to the intraspecific venom diversity observed in the Sonoran Desert of Arizona, USA: Venom Type A (∼Type II), with Mojave toxin, highly toxic, lacking hemorrhagic activity, and with scarce proteolytic activity; Type B (∼Type I), without Mojave toxin, less toxic than Type A, highly hemorrhagic and proteolytic; and Type A + B, containing Mojave toxin, as toxic as venom Type A, variable in hemorrhagic activity and with intermediate proteolytic activity. We also detected a positive correlation between SVMP abundance and hemorrhagic and proteolytic activities. Although more sampling is necessary, our results suggest that venoms containing Mojave toxin and venom lacking this toxin are distributed in the northwest and southeast portions of the distribution in Mexico, respectively, while an intergradation in the middle of both zones is present., Competing Interests: The authors declare no conflict of interest. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Two acidic Asp 49 -PLA 2 s with Glu 6 substitution and a neutral Lys 49 -PLA (designated Gst-K49) were cloned from G. strauchii venom glands, their full amino acid sequences were deduced. The predominant acidic PLA 2 (designated Gst-E6a) contains 124 residues and the M 18 W 30 substitutions, while the minor acidic PLA 2 (designated Gst-E6b) contains 122 residues and the V 18 A 30 substitutions. Their sequences are most similar to those of the respective orthologous PLA 2 s of G. intermedius venom. Gst-E6a and Gst-E6b appear to be paralogs and possibly have different predatory targets or functions. The LC-MS/MS results indicate the presence of only three PLA 2 gene products in the crude venom, the relative expression levels were in the order of Gst-E6a ≫ Gst-E6b > Gst-K49, as confirmed by qPCR results. In contrast to other Gloydius, G. strauchii venom does not contain neurotoxic or basic anticoagulant Asp 49 -PLA 2 s, but Gst-K49 is the first Lys 49 -PLA 2 identified in Gloydius venoms. However, its venom content is relatively low and its pI value 7.3 is much lower than those of other Lys 49 -PLA 2 s and. The Lys 49 -PLA 2 genes appear to regress in the venom of most of Gloydius and related rattlesnake, and this evolutionary regression occurred before the dispersal of Asian pitvipers to the New World., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Snakebites are a serious health problem in tropical countries. In Brazil, the genus Bothrops (Viperidae family) causes most of the ophidic accidents, characterized by proteolysis and haemorrhage. Snake venoms are rich sources of toxins with great therapeutic and biotechnological potential and omics approaches is a valuable tool for identification of new bioactive components in the venom. In this study, we described the first transcriptome of the venom gland of Bothrops moojeni snake, using the next-generation sequencing with the Illumina platform. We identified: (i) 20 venom components classes, among which metalloproteases were the most expressed ones, followed by serine proteases and phospholipases; and (ii) the 33 full-length amino acid sequences of toxins that have never been reported before in B. moojeni venom, such as one cysteine-rich secretory protein (Moojin), two hyaluronidases (BmooHyal-1 and BmooHyal-2), and one three-finger toxin (Bmoo-3FTx). Altogether, the transcripts identified herein represent a starting point for the analysis of structure-function relationships of toxins, which shall help develop novel biological tools and therapeutic drugs., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Adult rattlesnakes within genus Crotalus express one of two distinct venom phenotypes, type I (hemorrhagic) and type II (neurotoxic). In Costa Rican Central American rattlesnake, ontogenetic changes in the concentration of miRNAs modulate venom type II to type I transition. Venomics and venom gland transcriptome analyses showed that adult C. simus and C. tzabcan expressed intermediate patterns between type II and type I venoms, whereas C. culminatus had a canonical type I venom. Neonate/juvenile and adult Mexican rattlesnakes showed notable inter- and intraspecific variability in the number, type, abundance and ontogenetic shifts of the transcriptional and translational venom gland activities. These results support a role for miRNAs in the ontogenetic venom compositional changes in the three congeneric Mexican rattlesnakes. It is worth noting the finding of dual-action miRNAs, which silence the translation of neurotoxic heterodimeric PLA 2 crotoxin and acidic PLA 2 mRNAs while simultaneously up-regulating SVMP-targeting mRNAs. Dual transcriptional regulation potentially explains the existence of mutually exclusive crotoxin-rich (type-II) and SVMP-rich (type-I) venom phenotypic dichotomy among rattlesnakes. Our results support the hypothesis that alterations of the distribution of miRNAs, modulating the translational activity of venom gland toxin-encoding mRNAs in response to an external cue, may contribute to the mechanism generating adaptive venom variability.

In our previous work, a thrombin-like enzyme (TLE), agkihpin, was successfully isolated, purified, cloned and named from the venom of Gloydius halys Pallas, having fibrinolytic, fibrinogenolytic and thrombosis-reduced activities, attenuating migration of liver cancer cell, and without bleeding risk. To explore the possibility of agkihpin as a thrombolytic and/or anti-metastasis agent in the future, in this study recombinant agkihpin was expressed and purified in Escherichia coli, and its biological activities investigated. Thus, r-agkihpin-2 was successfully expressed and purified and confirmed by Western blot and peptide mass fingerprinting. After purification and renaturation, 46 mg (399 U) of active r-agkihpin-2 was obtained from 1 L bacterial culture. The results of the arginine esterase activity assay, fibrin plate test fibrinogenolytic activity assay, thrombin-induced venous thrombosis assay, Scratch-Wound assay and bleeding assay showed that active r-agkihpin-2 had slightly lower TAME hydrolytic, fibrinolytic, fibrinogenolytic, thrombus-reduced and migration-attenuated activities than those of native agkihpin, and had no bleeding risk. These findings confirmed that, active r-agkihpin-2 could be further investigated for thrombolytic and/or anti-metastasis drug discovery in the future., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Nine distinct venom serine proteases (vSPs) of Gloydius intermedius were studied by transcriptomic, sub-proteomic and phylogenetic analyses. Their complete amino acid sequences were deduced after Expression Sequence Tag (EST) analyses followed by cDNA cloning and sequencing. These vSPs appear to be paralogs and contain the catalytic triads and 1-4 potential N-glycosylation sites. Their relative expression levels evaluated by qPCR were grossly consistent with their EST hit-numbers. The major vSPs were purified by HPLC and their N-terminal sequences matched well to the deduced sequences, while fragments of the minor vSPs were detected by LC-MS/MS identification. Specific amidolytic activities of the fractions from HPLC and anion exchange separation were assayed using four chromogenic substrates, respectively. Molecular phylogenetic tree based on the sequences of these vSPs and their orthologs revealed six major clusters, one of them covered four lineages of plasminogen activator like vSPs. N-glycosylation patterns and variations for the vSPs are discussed. The high sequence similarities between G. intermedius vSPs and their respective orthologs from American pitvipers suggest that most of the isoforms evolved before Asian pitvipers migrated to the New World. Our results also indicate that the neurotoxic venoms contain more kallikrein-like vSPs and hypotensive components than the hemorrhagic venoms., Significance: Full sequences and expression levels of nine paralogous serine proteases (designated as GiSPs) of Gloydius intermedius venom have been studied. A kallikrein-like enzyme is most abundant and four isoforms homologous to venom plasminogen-activators are also expressed in this venom. Taken together, the present and previous data demonstrate that the neurotoxic G. intermedius venoms contain more hypotensive vSPs relative to other hemorrhagic pitviper venoms and the pitviper vSPs are highly versatile and diverse. Their structure-function relationships remain to be explored and compared. A novel, simplified phylogenetic tree based on the sequences of GiSPs and their closely related orthologs from other pitvipers reveals six major subtypes and offers a better understanding of vSP duplication and evolution in pitvipers of both the Old and New Worlds. It is well known that specific vSPs are potential therapeutic or diagnostic agents that target the plasma proteins or coagulation factors. Our results not only render deeper insights into the variation and evolution of vSPs, but may help to choose right venoms for the development of better therapeutic leads., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Venom variability is commonly reported for venomous snakes including Bothrops atrox. Here, we compared the composition of venoms from B. atrox snakes collected at Amazonian conserved habitats (terra-firme upland forest and várzea) and human modified areas (pasture and degraded areas). Venom samples were submitted to shotgun proteomic analysis as a whole or compared after fractionation by reversed-phase chromatography. Whole venom proteomes revealed a similar composition among the venoms with predominance of SVMPs, CTLs, and SVSPs and intermediate amounts of PLA 2 s and LAAOs. However, when distribution of particular isoforms was analyzed by either method, the venom from várzea snakes showed a decrease in hemorrhagic SVMPs and an increase in SVSPs, and procoagulant SVMPs and PLA 2 s. These differences were validated by experimental approaches including both enzymatic and in vivo assays, and indicated restrictions in respect to antivenom efficacy to variable components. Thus, proteomic analysis at the isoform level combined to in silico prediction of functional properties may indicate venom biological activity. These results also suggest that the prevalence of functionally distinct isoforms contributes to the variability of the venoms and could reflect the adaptation of B. atrox to distinct prey communities in different Amazon habitats., Biological Significance: In this report, we compared isoforms present in venoms from snakes collected at different Amazonian habitats. By means of a species venom gland transcriptome and the in silico functional prediction of each isoform, we were able to predict the principal venom activities in vitro and in animal models. We also showed remarkable differences in the venom pools from snakes collected at the floodplain (várzea habitat) compared to other habitats. Not only was this venom less hemorrhagic and more procoagulant, when compared to the venom pools from the other three habitats studied, but also this enhanced procoagulant activity was not efficiently neutralized by Bothrops antivenom. Thus, using a functional proteomic approach, we highlighted intraspecific differences in B. atrox venom that could impact both in the ecology of snakes but also in the treatment of snake bite patients in the region., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Many animal toxins may target the same molecules that need to be controlled in certain pathologies; therefore, some toxins have led to the formulation of drugs that are presently used, and many other drugs are still under development. Nevertheless, collecting sufficient toxins from the original source might be a limiting factor in studying their biological activities. Thus, molecular biology techniques have been applied in order to obtain large amounts of recombinant toxins into Escherichia coli . However, most animal toxins are difficult to express in this system, which results in insoluble, misfolded, or unstable proteins. To solve these issues, toxins have been fused with tags that may improve protein expression, solubility, and stability. Among these tags, the SUMO (small ubiquitin-related modifier) has been shown to be very efficient and can be removed by the Ulp1 protease. However, removing SUMO is a labor- and time-consuming process. To enhance this system, here we show the construction of a bicistronic vector that allows the expression of any protein fused to both the SUMO and Ulp1 protease. In this way, after expression, Ulp1 is able to cleave SUMO and leave the protein interest-free and ready for purification. This strategy was validated through the expression of a new phospholipase D from the spider Loxosceles gaucho and a disintegrin from the Bothrops insularis snake. Both recombinant toxins showed good yield and preserved biological activities, indicating that the bicistronic vector may be a viable method to produce proteins that are difficult to express.

Tropidolaemus wagleri (temple pit viper) is a medically important snake in Southeast Asia. It displays distinct sexual dimorphism and prey specificity, however its venomics and inter-sex venom variation have not been thoroughly investigated. Applying reverse-phase HPLC, we demonstrated that the venom profiles were not significantly affected by sex and geographical locality (Peninsular Malaya, insular Penang, insular Sumatra) of the snakes. Essentially, venoms of both sexes share comparable intravenous median lethal dose (LD 50 ) (0.56-0.63 μg/g) and cause neurotoxic envenomation in mice. LCMS/MS identified six waglerin forms as the predominant lethal principles, comprising 38.2% of total venom proteins. Fourteen other toxin-protein families identified include phospholipase A 2 , serine proteinase, snaclec and metalloproteinase. In mice, HPLC fractions containing these proteins showed insignificant contribution to the overall venom lethality. Besides, the unique elution pattern of approximately 34.5% of non-lethal, low molecular mass proteins (3-5 kDa) on HPLC could be potential biomarker for this primitive crotalid species. Together, the study unveiled the venom proteome of T. wagleri that is atypical among many pit vipers as it comprises abundant neurotoxic peptides (waglerins) but little hemotoxic proteinases. The findings also revealed that the venom is relatively well conserved intraspecifically despite the drastic morphological differences between sexes.

Objective: To design a specific polyclonal antibody against Deinagkistrodon acutus venom (DA-pAb) by immunizating New Zealand white rabbits., Results: The IgG fraction was purified by affinity chromatography, and specific antibodies were purified by immunoaffinity chromatography. Polyclonal antibodies were subjected to ELISA and western blotting to evaluate their immune reactivity. We identified the mimotopes by screening a phage display 12-mer peptide library against D. acutus venom. After three rounds of biopanning with DA-pAb, 30 positive clones were identified. Eighteen phage clones were sequenced, and their corresponding amino acid sequences were deduced. Additional immunoassays with the peptides and DA-pAb identified five sequences as possible epitopes. Recombinant antigens synthesized with the five epitopes were used for the immunization of BALB/c mice., Conclusion: The antibodies induced by these peptides recognized the recombinant antigen and D. acutus venom and protected mice against the hemorrhagic effects of the venom.

The genetic origin of novel traits is a central but challenging puzzle in evolutionary biology. Among snakes, phospholipase A2 (PLA2)-related toxins have evolved in different lineages to function as potent neurotoxins, myotoxins, or hemotoxins. Here, we traced the genomic origin and evolution of PLA2 toxins by examining PLA2 gene number, organization, and expression in both neurotoxic and non-neurotoxic rattlesnakes. We found that even though most North American rattlesnakes do not produce neurotoxins, the genes of a specialized heterodimeric neurotoxin predate the origin of rattlesnakes and were present in their last common ancestor (∼22 mya). The neurotoxin genes were then deleted independently in the lineages leading to the Western Diamondback (Crotalus atrox) and Eastern Diamondback (C. adamanteus) rattlesnakes (∼6 mya), while a PLA2 myotoxin gene retained in C. atrox was deleted from the neurotoxic Mojave rattlesnake (C. scutulatus; ∼4 mya). The rapid evolution of PLA2 gene number appears to be due to transposon invasion that provided a template for non-allelic homologous recombination., (Copyright © 2016 Elsevier Ltd. All rights reserved.)