California Ground Squirrel (Spermophilus beecheyi) Defenses Against Rattlesnake Venom Digestive and Hemostatic Toxins

A genome assembly for the southern Pacific rattlesnake, Crotalus oreganus helleri, in the western rattlesnake species complex

Rattlesnakes play important roles in their ecosystems by regulating prey populations, are involved in complex coevolutionary dynamics with their prey, and exhibit a variety of unusual adaptations, including maternal care, heat-sensing pit organs, hinged fangs, and medically-significant venoms. The western rattlesnake (Crotalus oreganus) is one of the widest ranging rattlesnake species, with a distribution from British Columbia, where it is listed as threatened, to Baja California and east across the Great Basin to western Wyoming, Colorado and New Mexico. Here, we report a new reference genome assembly for one of six currently recognized subspecies, C. oreganus helleri, as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomic sequencing strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 698 scaffolds spanning 1,564,812,557 base pairs, has a contig N50 of 64.7 Mb, a scaffold N50 of 110.8 Mb, and BUSCO complete score of 90.5%. This reference genome will be valuable for studies on the genomic basis of venom evolution and variation within Crotalus, in resolving the taxonomy of C. oreganus and its relatives, and for the conservation and management of rattlesnakes in general.

Stings on wings: Proteotranscriptomic and biochemical profiling of the lesser banded hornet (Vespa affinis) venom

Distinct animal lineages have convergently recruited venoms as weaponry for prey capture, anti-predator defence, conspecific competition, or a combination thereof. Most studies, however, have been primarily confined to a narrow taxonomic breadth. The venoms of cone snails, snakes, spiders and scorpions remain particularly well-investigated. Much less explored are the venoms of wasps (Order: Hymenoptera) that are infamous for causing excruciating and throbbing pain, justifying their apex position on Schmidt's pain index, including some that are rated four on four. For example, the lesser banded wasp (V. affinis) is clinically important yet has only been the subject of a few studies, despite being commonly found across tropical and subtropical Asia. Stings from these wasps, especially from multiple individuals of a nest, often lead to clinically severe manifestations, including mastocytosis, myasthenia gravis, optic neuropathy, and life-threatening pathologies such as myocardial infarction and organ failure. However, their venom composition and activity remain unexplored in the Indian subcontinent. Here, we report the proteomic composition, transcriptomic profile, and biochemical and pharmacological activities of V. affinis venom from southern India. Our findings suggest that wasp venoms are rich in diverse toxins that facilitate antipredator defence. Biochemical and pharmacological assessments reveal that these toxins can exhibit significantly higher activities than their homologues in medically important snakes. Their ability to exert potent effects on diverse molecular targets makes them a treasure trove for discovering life-saving therapeutics. Fascinatingly, wasp venoms, being evolutionarily ancient, exhibit a greater degree of compositional and sequence conservation across very distant populations/species, which contrasts with the patterns of venom evolution observed in evolutionarily younger lineages, such as advanced snakes and cone snails.

The roles of balancing selection and recombination in the evolution of rattlesnake venom

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.

Genetic characterization of potential venom resistance proteins in California ground squirrels ( Otospermophilus beecheyi ) using transcriptome analyses

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.

Phenotypic and functional variation in venom and venom resistance of two sympatric rattlesnakes and their prey

Predator-prey interactions often lead to the coevolution of adaptations associated with avoiding predation and, for predators, overcoming those defences. Antagonistic coevolutionary relationships are often not simple interactions between a single predator and prey but rather a complex web of interactions between multiple coexisting species. Coevolution between venomous rattlesnakes and small mammals has led to physiological venom resistance in several mammalian taxa. In general, viperid venoms contain large quantities of snake venom metalloproteinase toxins (SVMPs), which are inactivated by SVMP inhibitors expressed in resistant mammals. We explored variation in venom chemistry, SVMP expression, and SVMP resistance across four co-distributed species (California Ground Squirrels, Bryant's Woodrats, Southern Pacific Rattlesnakes, and Red Diamond Rattlesnakes) collected from four different populations in Southern California. Our aim was to understand phenotypic and functional variation in venom and venom resistance in order to compare coevolutionary dynamics of a system involving two sympatric predator-prey pairs to past studies that have focused on single pairs. Proteomic analysis of venoms indicated that these rattlesnakes express different phenotypes when in sympatry, with Red Diamonds expressing more typical viperid venom (with a diversity of SVMPs) and Southern Pacifics expressing a more atypical venom with a broader range of non-enzymatic toxins. We also found that although blood sera from both mammals were generally able to inhibit SVMPs from both rattlesnake species, inhibition depended strongly on the snake population, with snakes from one geographic site expressing SVMPs to which few mammals were resistant. Additionally, we found that Red Diamond venom, rather than woodrat resistance, was locally adapted. Our findings highlight the complexity of coevolutionary relationships between multiple predators and prey that exhibit similar offensive and defensive strategies in sympatry.

Bacterial expression of a snake venom metalloproteinase inhibitory protein from the North American opossum (D.virginiana)

A variety of opossum species are resistant to snake venoms due to the presence of antihemorrhagic and antimyotoxic acidic serum glycoproteins that inhibit several toxic venom components. Two virtually identical antihemorrhagic proteins isolated from either the North American opossum (D. virginiana) or the South American big-eared opossum (D. aurita), termed oprin or DM43 respectively, inhibit specific snake venom metalloproteinases (SVMPs). A better understanding of the structure of these proteins may provide useful insight to determine their mechanism of action and for the development of therapeutics against the global health concern of snake-bite envenomation. The aim of this work is to produce a recombinant snake venom metalloproteinase inhibitor (SVMPI) similar to the above opossum proteins in Escherichia coli and determine if this bacterially produced protein inhibits the proteolytic properties of Western Diamondback rattlesnake (C. atrox) venom. The resulting heterologous SVMPI was produced with either a 6-Histidine or maltose binding protein (MBP) affinity tag on either the C-terminus or N-terminus of the protein, respectively. The presence of the solubility enhancing MBP affinity tag resulted in significantly more soluble protein expression. The inhibitory activity was measured using two complementary assays and the MBP labeled SVMPI showed 7-fold less activity as compared to the 6-Histidine labeled SVMPI. Thus, the bacterially derived SVMPI with an unlabeled N-terminus showed high inhibitory activity (IC₅₀ = 4.5 μM). The use of a solubility enhancing MBP fusion protein construct appears to be a productive way to express sufficient quantities of this mammalian protein in E. coli for further study.

The molecular basis of venom resistance in a rattlesnake-squirrel predator-prey system

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.

Investigating the cytotoxic effects of the venom proteome of two species of the Viperidae family (Cerastes cerastes and Cryptelytrops purpureomaculatus) from various habitats

Animal secretions are of great interest in terms of drug development due to their complex protein and peptide composition. Especially, in the field of therapeutic medications such as anti-cancer drugs snake venoms receive attention. In this study, we address two Viperidae species from various habitats with a particular focus on the cytotoxic potential along with the decomplexation of the venom proteome: the horned desert viper (Cerastes cerastes), native to desert regions of North Africa and the mangrove pit viper (Cryptelytrops purpureomaculatus), found in coastal forests of Southeast Asia. Initial cytotoxic screenings of the crude venoms revealed diverse activity, with the highest effect against SHSY5Y human glioblastoma carcinoma cells compared to other cancerous and non-cancerous cell lines. In-depth cytotoxicity studies of SHSY5Y cells with purified venom fractions revealed heterodimeric disintegrins from C. cerastes venom, which exerted a high cytotoxic activity with IC₅₀ values from 0.11 to 0.58 μM and a disintegrin-like effect on SHSY5Y morphology was observed due to cell detachment. Furthermore, two polyproline BPP-related peptides, one PLA₂ and a peptide-rich fraction were determined for C. purpureomaculatus with moderate IC₅₀ values between 3 and 51 μM. Additionally, the decryption of the venom proteomes by snake venomic mass spectrometry and comparison of the same species from different habitats revealed slight differences in the composition.

Serum-based inhibition of pitviper venom by eastern indigo snakes (Drymarchon couperi)

When organisms possess chemical defenses, their predators may eventually evolve resistance to their toxins. Eastern indigo snakes (Drymarchon couperi; EIS) prey on pitvipers and are suspected to possess physiological resistance to their venom. In this study, we formally investigated this hypothesis using microassays that measured the ability of EIS blood sera to inhibit (A) hemolytic and (B) snake venom metalloproteinase (SVMP) activity of copperhead (Agkistrodon contortrix) venom. To serve as controls, we also tested the inhibitory ability of sera from house mice (Mus musculus) and checkered gartersnakes (Thamnophis marcianus), a snake that does not feed on pitvipers. Sera from both EIS and gartersnakes inhibited over 60% of SVMP activity, while only EIS sera also inhibited venom hemolytic activity (78%). Our results demonstrate that EIS serum is indeed capable of inhibiting two of the primary classes of toxins found in copperhead venom, providing the first empirical evidence suggesting that EIS possess physiological resistance to venom upon injection. Because we documented resistance to hemolytic components of pitviper venom within EIS but not gartersnakes, we speculate this resistance may be driven by selection from feeding on pitvipers while resistance to SVMP may be relatively widespread among snakes.

Summary: We confirm eastern indigo snakes possess a physiological resistance to pitviper venom and speculate that venom resistance in snakes may be a product of antagonistic interactions and phylogenetically conserved traits.

Venom Resistance as a Model for Understanding the Molecular Basis of Complex Coevolutionary Adaptations

SynopsisVenom and venom resistance are molecular phenotypes widely considered to have diversified through coevolution between predators and prey. However, while evolutionary and functional studies on venom have been extensive, little is known about the molecular basis, variation, and complexity of venom resistance. We review known mechanisms of venom resistance and relate these mechanisms to their predicted impact on coevolutionary dynamics with venomous enemies. We then describe two conceptual approaches which can be used to examine venom/resistance systems. At the intraspecific level, tests of local adaptation in venom and resistance phenotypes can identify the functional mechanisms governing the outcomes of coevolution. At deeper evolutionary timescales, the combination of phylogenetically informed analyses of protein evolution coupled with studies of protein function promise to elucidate the mode and tempo of evolutionary change on potentially coevolving genes. We highlight case studies that use each approach to extend our knowledge of these systems as well as address larger questions about coevolutionary dynamics. We argue that resistance and venom are phenotypic traits which hold exceptional promise for investigating the mechanisms, dynamics, and outcomes of coevolution at the molecular level. Furthermore, extending the understanding of single gene-for-gene interactions to the whole resistance and venom phenotypes may provide a model system for examining the molecular and evolutionary dynamics of complex multi-gene interactions.

Therapeutic potential of snake venom in cancer therapy: current perspectives

Many active secretions produced by animals have been employed in the development of new drugs to treat diseases such as hypertension and cancer. Snake venom toxins contributed significantly to the treatment of many medical conditions. There are many published studies describing and elucidating the anti-cancer potential of snake venom. Cancer therapy is one of the main areas for the use of protein peptides and enzymes originating from animals of different species. Some of these proteins or peptides and enzymes from snake venom when isolated and evaluated may bind specifically to cancer cell membranes, affecting the migration and proliferation of these cells. Some of substances found in the snake venom present a great potential as anti-tumor agent. In this review, we presented the main results of recent years of research involving the active compounds of snake venom that have anticancer activity.

Interacting impacts of invasive plants and invasive toads on native lizards

The ecological impacts of an invasive species may be reduced by prior invasions if selective pressures imposed by earlier events preadapt the native biota to deal with the newer arrival. In northwestern Australia, invasion of the cane toad (Rhinella marina) kills many native predators if they ingest the highly toxic toads. Remarkably, the toads' defensive toxins (bufadienolides) are chemically similar to those of another invasive species: an ornamental plant from Madagascar, Bryophyllum spp. (Crassulaceae, mother-of-millions). Omnivorous lizards (bluetongue skinks, Tiliqua scincoides) are imperiled by the invasion of toads in northwestern Australia, but conspecifics from other areas of the continent (those where exotic plants were introduced and including areas where toads have yet to invade) are less affected because they exhibit higher physiological tolerance of toad toxins (and also of plant toxins). The willingness of captive bluetongues to consume both toads and these plants and the high correlation in the lizards' sensitivity to toad toxins versus plant toxins suggest that exotic plants may have imposed strong selection on the lizards' physiological tolerance of bufadienolides. As a result, populations of lizards from areas previously exposed to these alien plants may be preadapted to deal with the toxins of the more recent anuran invader.

A limited role for gene duplications in the evolution of platypus venom

Gene duplication followed by adaptive selection is believed to be the primary driver of venom evolution. However, to date, no studies have evaluated the importance of gene duplications for venom evolution using a genomic approach. The availability of a sequenced genome and a venom gland transcriptome for the enigmatic platypus provides a unique opportunity to explore the role that gene duplication plays in venom evolution. Here, we identify gene duplication events and correlate them with expressed transcripts in an in-season venom gland. Gene duplicates (1,508) were identified. These duplicated pairs (421), including genes that have undergone multiple rounds of gene duplications, were expressed in the venom gland. The majority of these genes are involved in metabolism and protein synthesis not toxin functions. Twelve secretory genes including serine proteases, metalloproteinases, and protease inhibitors likely to produce symptoms of envenomation such as vasodilation and pain were detected. Only 16 of 107 platypus genes with high similarity to known toxins evolved through gene duplication. Platypus venom C-type natriuretic peptides and nerve growth factor do not possess lineage-specific gene duplicates. Extensive duplications, believed to increase the potency of toxic content and promote toxin diversification, were not found. This is the first study to take a genome-wide approach in order to examine the impact of gene duplication on venom evolution. Our findings support the idea that adaptive selection acts on gene duplicates to drive the independent evolution and functional diversification of similar venom genes in venomous species. However, gene duplications alone do not explain the "venome" of the platypus. Other mechanisms, such as alternative splicing and mutation, may be important in venom innovation.

A rapid and sensitive fluorometric method for the quantitative analysis of snake venom metalloproteases and their inhibitors

Metalloproteases are responsible for the hemorrhagic effects of many snake venoms and contribute to other pathways that lead to local tissue damage. Methods that quantify snake venom metalloproteases (SVMP) are therefore valuable tools in research on the clinical, physiological, and biochemical effects of envenomation. Comparative analysis of individual, population, and species differences requires screening of large numbers of samples and treatments, and therefore require a method of quantifying SVMP activity that is simple, rapid, and sensitive. This paper demonstrates the properties of a new fluorometric assay of SVMP activity that can provide a measure of metalloprotease activity in 1 h. The assay is reliable, with variation among replicates sufficiently small to reliably detect differences in between species (F(19,60) = 2924, p < 0.001), even for those venoms with low overall activity. It is also sensitive enough to detect differences among venoms using <2 ng of whole venom protein. We provide an example use of this assay to detect the presence of natural SVMP inhibitors in minute samples of blood plasma from rock squirrels (S. variegatus), a natural prey species for North American rattlesnakes. We propose this assay is a useful addition to the set of tools used to characterize venoms, as well as high-throughput screening of natural or synthetic inhibitors, or other novel therapeutic agents against SVMP effects.

Venom Proteomes of Closely RelatedSistrurusRattlesnakes with Divergent Diets

The protein composition of the venoms of the three subspecies of Sistrurus catenatus (S. c. catenatus, tergeminus, and edwardsii) and a basal species, Sistrurus miliarius barbouri, were analyzed by RP-HPLC, N-terminal sequencing, MALDI-TOF peptide mass fingerprinting, and CID-MS/MS. The venoms of the four Sistrurus taxa contain proteins from 11 families. The protein family profile and the relative abundance of each protein group in the different venoms are not conserved. Myotoxins and 2-chain PLA2s were detected only in S.c. catenatus and S.c. tergeminus, whereas C-type BPP and Kunitz-type inhibitors were exclusively found in S.c. edwardsii and Sistrurus miliarius barbouri. Among major protein families, taxa were most similar in their metalloproteases (protein similarity coefficient value: 34%) and most divergent in PLA2s (12%), with values for disintegrins and serine proteases lying between these extremes (25 and 20%, respectively). The patterns of venom diversity points to either a gain in complexity in S. catenatus taxa or a loss of venom diversity occurring early on in the evolution of the group involving the lineage connecting S. milarius to the other taxa. The high degree of differentiation in the venom proteome among recently evolved congeneric taxa emphasizes the uniqueness of the venom composition of even closely related species that have different diets. Comparative proteomic analysis of Sistrurus venoms provides a comprehensive catalog of secreted proteins, which may contribute to a deeper understanding of the biology and ecology of these North American snakes and may also serve as a starting point for studying structure-function correlations of individual toxins.