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

Unique physiological and regulatory activity drives divergent toxin and non-toxin gene expression in rattlesnake accessory venom glands

Understanding the mechanisms by which organs and tissues evolve new physiological functions is central to understanding the evolution of novelty. This is particularly interesting in the context of related tissues that evolve specialized, yet complementary, functions. Snake venom glands are an attractive system to test hypotheses related to the evolution and specialization of novel physiological function, as these modified salivary glands have evolved over ∼60 MY to synthesize and store venom. Front-fanged venomous snakes (elapids and viperids) possess two types of venom glands: the main and accessory glands. The larger main gland produces greater quantities of venom toxins and has been studied extensively, while the smaller accessory gland has received less attention. Here, we explore gene expression differences between main and accessory venom glands across three rattlesnake species (Crotalus cerberus, C. oreganus concolor and C. viridis). Our findings indicate that accessory glands express most venom genes at significantly lower levels than the main gland, with a few exceptions that may represent biologically relevant contributions of accessory glands to venom. The two glands also exhibit distinct trans-regulatory environments that we link to key differences in their underlying physiology and secretory roles. Our results further suggest that two signaling pathways that regulate venom, the unfolded protein response (UPR) and extracellular signal-regulated kinase (ERK), show significantly lower activation in the accessory gland. These findings provide insight into the physiological and functional diversification of snake venom systems, highlighting how distinct glandular systems have evolved contrasting and complementary roles driven by distinct physiological and regulatory mechanisms.

Snake Venomics of the Arboreal Talamancan Palm-Pitviper, Bothriechis nubestris, Provides Clues on the Origin of a Phenotypic Dichotomy between Type-I and Type-II Venoms

We report the biochemical and proteomic characterization of the venom of the arboreal Talamancan palm-pitviper, Bothriechis nubestris, a species endemic to north and central parts of the Costa Rican Cordillera de Talamanca, at 2100 to 3000 m above sea level. The Talamancan palm-pitviper venom arsenal comprised the translated products of 26 unique transcripts into 10 toxin classes, where metalloproteinases PIIa and PIII represent the dominant components. In vitro, the venom proteolyzed azocasein and gelatin but showed no phospholipase A2 or human plasma coagulant activities. In vivo, B. nubestris venom exhibited an intravenous median lethal dose (LD50) of 21.5 (95% confidence interval: 15.6-29.5) μg/mouse, and a minimum hemorrhagic dose (MHD) of 1.85 μg. PoliVal-ICP antivenom neutralized the venom's lethal activity with a potency of 9.7 mg of venom/g of antivenom and significantly reduced the hemorrhagic effect. Comparison of venom gland transcriptomes and venom proteomes of B. nubestris and its closest congeneric relative, B. nigroviridis, revealed that highly conserved venom gland transcriptomes are differentially processed by each species to produce divergent (type-I vs type-II, respectively) venoms. This phenomenon contributes to the remarkable venom phenotypic variability found across the palm-pitviper phylogeny. A possible mechanism for the occurrence of type-I/type-II venom phenotypic dichotomy is discussed.

Ambient temperature and toxic diets constrain snake venom resistance in a desert rodent

Variation in ambient temperature and food availability is commonplace in nature and likely exerts several types of eco-evolutionary pressures that mediate species' interactions. In predator-prey interactions between snakes and rodents, the molecular interface is formed by snake venom and mammalian venom resistance proteins. However, the abiotic factors modulating resistance are poorly understood. Here, we measured serum-based venom resistance of desert woodrats (Neotoma lepida) maintained at either cool or warm ambient temperatures and fed creosote bush resin (native diet) or a control diet. Woodrat serum was collected and tested for its ability to inhibit the activity of rattlesnake venom metalloproteinases. Woodrats raised at cooler temperatures, as well as those consuming diets with creosote resin, were significantly less able to inhibit snake venom, suggesting that they would be more susceptible to snakebite. These results suggest that temperature and dietary variation across the distribution over which these rattlesnakes and woodrats interact could structure the outcomes of these predator-prey interactions. Additionally, these results may help explain why ambient temperatures, rather than dietary differences, predict the presence of neurotoxic versus proteolytic venom phenotypes in some rattlesnake species.

Raptors without resistance: No evidence for endogenous inhibition of rattlesnake venom metalloproteinases in a Great Plains raptor assemblage

Snake venoms are complex biochemical secretions under strong selection for prey subjugation, and venoms are tightly linked to the biotic communities that snakes inhabit. Physiological adaptations for venom resistance have been identified in various snake prey species, but fewer snake predators, with research in this area largely biased towards mammalian species. Fewer investigations have assayed for the presence of resistance mechanisms in avian systems. Birds of prey (hereafter "raptors"; orders Accipitriformes, Falconiformes, and Strigiformes) represent major sources of predation for snakes. Raptor dietary habits range from snake specialists to non-snake feeders, and this continuum of snake predation frequency among species creates the ideal system in which to explore the presence and strength of venom resistance. We assayed sera from a suite of Great Plains raptors against snake venom metalloproteinases (SVMPs) of the Prairie rattlesnake (Crotalus v. viridis) to test the general hypotheses that 1) raptor sera will display elevated SVMP inhibition compared to a naïve avian model (domestic chicken; Gallus gallus) and 2) raptor species with high levels of rattlesnake predation will more effectively inhibit SVMP activity than those that are not known to feed on rattlesnakes. We found that raptors do possess elevated SVMP inhibition in comparison to a naïve avian model, but this level of inhibition remains low and is unlikely to be biologically significant in detoxifying venoms. We found no evidence suggesting that inhibitory potential of different raptor sera corresponds to the level of rattlesnake predation associated with each species. The widespread lack of SVMP inhibition in diverse raptors underscores the complexity of venom resistance dynamics in natural systems and further suggests that physiological venom resistance mechanisms may be poorly developed in birds more broadly.

Tandem duplication of serpin genes yields functional variation and snake venom inhibitors

Tandem duplication of genes can play a critical role in the evolution of functional novelty, but our understanding is limited concerning gene duplication's role in coevolution between species. Much is known about the evolution and function of tandemly duplicated snake venom genes, however the potential of gene duplication to fuel venom resistance within prey species is poorly understood. In this study, we characterize patterns of gene duplication of the SERPINA subfamily of genes across in vertebrates and experimentally characterize functional variation in the SERPINA3-like paralogs of a wild rodent. We find the hallmarks of rapid birth-death evolution of SERPINA1-like and SERPINA3-like genes within and between rodent lineages. Next, we recombinantly expressed the 2 paralogous duplicates of SERPINA1 and 12 paralogous duplicates of SERPINA3 found in the genome of the big-eared woodrat (Neotoma macrotis), a species known to be resistant to protease-rich rattlesnake venoms. We found that two SERPINA3 paralogs inhibit snake venom serine proteases, indicating that these proteins have potential as resistance factors in SERPIN-mediated venom resistance. In addition, functional variation is apparent among paralogs, including neofunctionalization to inhibit both chymotrypsin-like and and trypsin-like proteases simultaneously for one venom-inhibiting paralog. Our results provide further evidence that the rapid evolution of SERPINA1 and SERPINA3 gene copy number across rodents has adaptive potential by producing functionally-diverse inhibitors.

Prey specificity of predatory venoms

Venom represents a key adaptation of many venomous predators, allowing them to immobilise prey quickly through chemical rather than physical warfare. Evolutionary arms races between prey and a predator are believed to be the main factor influencing the potency and composition of predatory venoms. Predators with narrowly restricted diets are expected to evolve specifically potent venom towards their focal prey, with lower efficacy on alternative prey. Here, we evaluate hypotheses on the evolution of prey-specific venom, focusing on the effect of restricted diet, prey defences, and prey resistance. Prey specificity as a potential evolutionary dead end is also discussed. We then provide an overview of the current knowledge on venom prey specificity, with emphasis on snakes, cone snails, and spiders. As the current evidence for venom prey specificity is still quite limited, we also overview the best approaches and methods for its investigation and provide a brief summary of potential model groups. Finally, possible applications of prey-specific toxins are discussed.

High-Voltage Toxin'Roll: Electrostatic Charge Repulsion as a Dynamic Venom Resistance Trait in Pythonid Snakes

The evolutionary interplay between predator and prey has significantly shaped the development of snake venom, a critical adaptation for subduing prey. This arms race has spurred the diversification of the components of venom and the corresponding emergence of resistance mechanisms in the prey and predators of venomous snakes. Our study investigates the molecular basis of venom resistance in pythons, focusing on electrostatic charge repulsion as a defense against α-neurotoxins binding to the alpha-1 subunit of the postsynaptic nicotinic acetylcholine receptor. Through phylogenetic and bioactivity analyses of orthosteric site sequences from various python species, we explore the prevalence and evolution of amino acid substitutions that confer resistance by electrostatic repulsion, which initially evolved in response to predatory pressure by Naja (cobra) species (which occurs across Africa and Asia). The small African species Python regius retains the two resistance-conferring lysines (positions 189 and 191) of the ancestral Python genus, conferring resistance to sympatric Naja venoms. This differed from the giant African species Python sebae, which has secondarily lost one of these lysines, potentially due to its rapid growth out of the prey size range of sympatric Naja species. In contrast, the two Asian species Python brongersmai (small) and Python bivittatus (giant) share an identical orthosteric site, which exhibits the highest degree of resistance, attributed to three lysine residues in the orthosteric sites. One of these lysines (at orthosteric position 195) evolved in the last common ancestor of these two species, which may reflect an adaptive response to increased predation pressures from the sympatric α-neurotoxic snake-eating genus Ophiophagus (King Cobras) in Asia. All these terrestrial Python species, however, were less neurotoxin-susceptible than pythons in other genera which have evolved under different predatory pressure as: the Asian species Malayopython reticulatus which is arboreal as neonates and juveniles before rapidly reaching sizes as terrestrial adults too large for sympatric Ophiophagus species to consider as prey; and the terrestrial Australian species Aspidites melanocephalus which occupies a niche, devoid of selection pressure from α-neurotoxic predatory snakes. Our findings underline the importance of positive selection in the evolution of venom resistance and suggest a complex evolutionary history involving both conserved traits and secondary evolution. This study enhances our understanding of the molecular adaptations that enable pythons to survive in environments laden with venomous threats and offers insights into the ongoing co-evolution between venomous snakes and their prey.

A plethora of rodents: Rattlesnake predators generate unanticipated patterns of venom resistance in a grassland ecosystem

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

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium

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.

Resistance Is Not Futile: Widespread Convergent Evolution of Resistance to Alpha-Neurotoxic Snake Venoms in Caecilians (Amphibia: Gymnophiona)

Predatory innovations impose reciprocal selection pressures upon prey. The evolution of snake venom alpha-neurotoxins has triggered the corresponding evolution of resistance in the post-synaptic nicotinic acetylcholine receptors of prey in a complex chemical arms race. All other things being equal, animals like caecilians (an Order of legless amphibians) are quite vulnerable to predation by fossorial elapid snakes and their powerful alpha-neurotoxic venoms; thus, they are under strong selective pressure. Here, we sequenced the nicotinic acetylcholine receptor alpha-1 subunit of 37 caecilian species, representing all currently known families of caecilians from across the Americas, Africa, and Asia, including species endemic to the Seychelles. Three types of resistance were identified: (1) steric hindrance from N-glycosylated asparagines; (2) secondary structural changes due to the replacement of proline by another amino acid; and (3) electrostatic charge repulsion of the positively charged neurotoxins, through the introduction of a positively charged amino acid into the toxin-binding site. We demonstrated that resistance to alpha-neurotoxins convergently evolved at least fifteen times across the caecilian tree (three times in Africa, seven times in the Americas, and five times in Asia). Additionally, as several species were shown to possess multiple resistance modifications acting synergistically, caecilians must have undergone at least 20 separate events involving the origin of toxin resistance. On the other hand, resistance in non-caecilian amphibians was found to be limited to five origins. Together, the mutations underlying resistance in caecilians constitute a robust signature of positive selection which strongly correlates with elapid presence through both space (sympatry with caecilian-eating elapids) and time (Cenozoic radiation of elapids). Our study demonstrates the extent of convergent evolution that can be expected when a single widespread predatory adaptation triggers parallel evolutionary arms races at a global scale.

One Size Fits All-Venomics of the Iberian Adder (Vipera seoanei, Lataste 1878) Reveals Low Levels of Venom Variation across Its Distributional Range

European vipers (genus Vipera) are medically important snakes displaying considerable venom variation, occurring at different levels in this group. The presence of intraspecific venom variation, however, remains understudied in several Vipera species. Vipera seoanei is a venomous snake endemic to the northern Iberian Peninsula and south-western France, presenting notable phenotypic variation and inhabiting several diverse habitats across its range. We analysed the venoms of 49 adult specimens of V. seoanei from 20 localities across the species' Iberian distribution. We used a pool of all individual venoms to generate a V. seoanei venom reference proteome, produced SDS-PAGE profiles of all venom samples, and visualised patterns of variation using NMDS. By applying linear regression, we then assessed presence and nature of venom variation between localities, and investigated the effect of 14 predictors (biological, eco-geographic, genetic) on its occurrence. The venom comprised at least 12 different toxin families, of which five (i.e., PLA2, svSP, DI, snaclec, svMP) accounted for about 75% of the whole proteome. The comparative analyses of the SDS-PAGE venom profiles showed them to be remarkably similar across the sampled localities, suggesting low geographic variability. The regression analyses suggested significant effects of biological and habitat predictors on the little variation we detected across the analysed V. seoanei venoms. Other factors were also significantly associated with the presence/absence of individual bands in the SDS-PAGE profiles. The low levels of venom variability we detected within V. seoanei might be the result of a recent population expansion, or of processes other than directional positive selection.

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.

Ancestrally Reconstructed von Willebrand Factor Reveals Evidence for Trench Warfare Coevolution between Opossums and Pit Vipers

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.

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.

Convergent evolution of toxin resistance in animals

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

A comparative study of endogenous phospholipase A2 inhibitors in the serum of Brazilian pit vipers

This work compared the presence of phospholipase A₂ inhibitors (PLIs) in the serum of 19 snake species maintained at Instituto Butantan to better understand the mechanisms of venom resistance in snakes and improve the treatment of snakebite. PLI was isolated from blood of 19 snake species by one-step chromatography and identified in all samples, besides its identity was confirmed through the interaction with both phospholipase A₂ and anti-γPLI. These findings highlight the diversity of snake serum PLIs and emphasize the importance of structure-function studies.

Venom gene sequence diversity and expression jointly shape diet adaptation in pitvipers

Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA₂s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms—such as higher absolute levels of expression—are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey.

A Strategy for Efficient Preparation of Genus-Specific Diagnostic Antibodies for Snakebites

As said by former United Nations Secretary-General Kofi Annan, "Snakebite is the most important tropical disease you've never heard of." Listed as a priority neglected tropical disease by the World Health Organization, snakebite envenoming (SBE) kills in excess of 125,000 people per year. However, due to the complexity and overlap of snake venom compositions, few reliable venom diagnostic methods for genus-/species-specific identification, which is crucial for successful SBE therapy, are available. Here, we develop a strategy to select and prepare genus-specific snake venom antibodies, which allows rapid and efficient clinical diagnosis of snakebite. Multi-omics approaches are used to choose candidate antigens from snake venoms and identify genus-specific antigenic epitope peptide fragments (GSAEPs) with ideal immunogenicity, specificity, and spatial accessibility. Double-antibody sandwich ELISA kit was established by matching a polyclonal antibody against a natural antigen and a monoclonal antibody that was prepared by natural protein as antigen and can specifically target the GSAEPs. The kit shows the ability to accurately identify venoms from similar genera of Trimeresurus and Protobothrops with a detection limit of 6.25 ng/ml on the snake venoms and a little cross-reaction, thus proving high feasibility and applicability.

Evolutionary, proteomic, and experimental investigations suggest the extracellular matrix of cumulus cells mediates fertilization outcomes†

Studies of fertilization biology often focus on sperm and egg interactions. However, before gametes interact, mammalian sperm must pass through the cumulus layer; in mice, this consists of several thousand cells tightly glued together with hyaluronic acid and other proteins. To better understand the role of cumulus cells and their extracellular matrix, we perform proteomic experiments on cumulus oophorus complexes (COCs) in house mice (Mus musculus), producing over 24,000 mass spectra to identify 711 proteins. Seven proteins known to stabilize hyaluronic acid and the extracellular matrix were especially abundant (using spectral counts as an indirect proxy for abundance). Through comparative evolutionary analyses, we show that three of these evolve rapidly, a classic signature of genes that influence fertilization rate. Some of the selected sites overlap regions of the protein known to impact function. In a follow-up experiment, we compared COCs from females raised in two different social environments. Female mice raised in the presence of multiple males produced COCs that were smaller and more resistant to dissociation by hyaluronidase compared to females raised in the presence of a single male, consistent with a previous study that demonstrated such females produced COCs that were more resistant to fertilization. Although cumulus cells are often thought of as enhancers of fertilization, our evolutionary, proteomic, and experimental investigations implicate their extracellular matrix as a potential mediator of fertilization outcomes.

Evidence for Resistance to Coagulotoxic Effects of Australian Elapid Snake Venoms by Sympatric Prey (Blue Tongue Skinks) but Not by Predators (Monitor Lizards)

Some Australian elapids possess potently procoagulant coagulotoxic venoms which activate the zymogen prothrombin into the functional enzyme thrombin. Although the activity of Australian elapid prothrombin-activators has been heavily investigated with respect to the mammalian, and in particular, human clotting cascades, very few studies have investigated the activity of their venom upon reptile plasmas. This is despite lizards representing both the primary diet of most Australian elapids and also representing natural predators. This study investigated the procoagulant actions of a diverse range of Australian elapid species upon plasma from known prey species within the genera Tiliqua (blue tongue skinks) as well as known predator species within the genera Varanus (monitor lizards). In addition to identifying significant variation in the natural responses of the coagulation cascade between species from the genera Tiliqua and Varanus relative to each other, as well as other vertebrate lineages, notable differences in venom activity were also observed. Within the genus Tiliqua, both T. rugosa and T. scincoides plasma displayed significant resistance to the procoagulant activity of Pseudechis porphyriacus venom, despite being susceptible to all other procoagulant elapid venoms. These results indicate that T. rugosa and T. scincoides have evolved resistance within their plasma to the coagulotoxic venom activity of the sympatric species P. porphyriacus. Other venoms were able to activate Tiliqua prothrombin, which suggests that the lessened activity of P. porphyriacus venom is not due to modifications of the prothrombin and may instead be due to a serum factor that specifically binds to P. porphyriacus toxins, as has been previously seen for squirrels resistant to rattlesnake venom. In contrast, none of the predatory lizards studied (Varanus giganteus, V. mertensi and V. varius) demonstrated resistance to the venom. This suggests that the mechanical protection afforded by thick osteodermic scales, and prey handling behaviour, removes a selection pressure for the evolution of resistance in these large predatory lizards. These results therefore reveal differential interactions between venoms of snakes with sympatric lizards that are on opposite sides of the predator-prey arms race.

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.

A symphony of destruction: Dynamic differential fibrinogenolytic toxicity by rattlesnake (Crotalus and Sistrurus) venoms

What factors influence the evolution of a heavily selected functional trait in a diverse clade? This study adopts rattlesnakes as a model group to investigate the evolutionary history of venom coagulotoxicity in the wider context of phylogenetics, natural history, and biology. Venom-induced clotting of human plasma and fibrinogen was determined and mapped onto the rattlesnake phylogenetic tree to reconstruct the evolution of coagulotoxicity across the group. Our results indicate that venom phenotype is often independent of phylogenetic relationships in rattlesnakes, suggesting the importance of diet and/or other environmental variables in driving venom evolution. Moreover, the striking inter- and intraspecific variability in venom activity on human blood highlights the considerable variability faced by physicians treating envenomation. This study is the most comprehensive effort to date to describe and characterize the evolutionary and biological aspects of coagulotoxins in rattlesnake venom. Further research at finer taxonomic levels is recommended to elucidate patterns of variation within species and lineages.

Mutual enlightenment: A toolbox of concepts and methods for integrating evolutionary and clinical toxinology via snake venomics and the contextual stance

Snakebite envenoming is a neglected tropical disease that may claim over 100,000 human lives annually worldwide. Snakebite occurs as the result of an interaction between a human and a snake that elicits either a defensive response from the snake or, more rarely, a feeding response as the result of mistaken identity. Snakebite envenoming is therefore a biological and, more specifically, an ecological problem. Snake venom itself is often described as a "cocktail", as it is a heterogenous mixture of molecules including the toxins (which are typically proteinaceous) responsible for the pathophysiological consequences of envenoming. The primary function of venom in snake ecology is pre-subjugation, with defensive deployment of the secretion typically considered a secondary function. The particular composition of any given venom cocktail is shaped by evolutionary forces that include phylogenetic constraints associated with the snake's lineage and adaptive responses to the snake's ecological context, including the taxa it preys upon and by which it is predated upon. In the present article, we describe how conceptual frameworks from ecology and evolutionary biology can enter into a mutually enlightening relationship with clinical toxinology by enabling the consideration of snakebite envenoming from an "ecological stance". We detail the insights that may emerge from such a perspective and highlight the ways in which the high-fidelity descriptive knowledge emerging from applications of -omics era technologies - "venomics" and "antivenomics" - can combine with evolutionary explanations to deliver a detailed understanding of this multifactorial health crisis.

BoaγPLI from Boa constrictor Blood is a Broad-Spectrum Inhibitor of Venom PLA2 Pathophysiological Actions

The use of venom in predation exerts a corresponding selection pressure for the evolution of venom resistance. One of the mechanisms related to venom resistance in animals (predators or prey of snakes) is the presence of molecules in the blood that can bind venom toxins, and inhibit their pharmacological effects. One such toxin type are venom phospholipase A₂s (PLA₂s), which have diverse effects including anticoagulant, myotoxic, and neurotoxic activities. BoaγPLI isolated from the blood of Boa constrictor has been previously shown to inhibit venom PLA₂s that induced myotoxic and edematogenic activities. Recently, in addition to its previously described and very potent neurotoxic effect, the venoms of American coral snakes (Micrurus species) have been shown to have anticoagulant activity via PLA₂ toxins. As coral snakes eat other snakes as a major part of their diet, neonate Boas could be susceptible to predation by this sympatric species. Thus, this work aimed to ascertain if BoaγPLI provided a protective effect against the anticoagulant toxicity of venom from the model species Micrurus laticollaris in addition to its ability shown previously against other toxin types. Using a STA R Max coagulation analyser robot to measure the effect upon clotting time, and TEG5000 thromboelastographers to measure the effect upon clot strength, we evaluated the ability of BoaγPLI to inhibit M. laticollaris venom. Our results indicate that BoaγPLI is efficient at inhibiting the M. laticollaris anticoagulant effect, reducing the time of coagulation (restoring them closer to non-venom control values) and increasing the clot strength (restoring them closer to non-venom control values). These findings demonstrate that endogenous PLA₂ inhibitors in the blood of non-venomous snakes are multi-functional and provide broad resistance against a myriad of venom PLA₂-driven toxic effects including coagulotoxicity, myotoxicity, and neurotoxicity. This novel form of resistance could be evidence of selective pressures caused by predation from venomous snakes and stresses the need for field-based research aimed to expand our understanding of the evolutionary dynamics of such chemical arms race.

Rapid toxin sequestration modifies poison frog physiology

Poison frogs sequester chemical defenses from their diet of leaf litter arthropods for defense against predation. Little is known about the physiological adaptations that confer this unusual bioaccumulation ability. We conducted an alkaloid-feeding experiment with the Diablito poison frog (Oophaga sylvatica) to determine how quickly alkaloids are accumulated and how toxins modify frog physiology using quantitative proteomics. Diablito frogs rapidly accumulated the alkaloid decahydroquinoline within 4 days, and dietary alkaloid exposure altered protein abundance in the intestines, liver and skin. Many proteins that increased in abundance with decahydroquinoline accumulation are plasma glycoproteins, including the complement system and the toxin-binding protein saxiphilin. Other protein classes that change in abundance with decahydroquinoline accumulation are membrane proteins involved in small molecule transport and metabolism. Overall, this work shows that poison frogs can rapidly accumulate alkaloids, which alter carrier protein abundance, initiate an immune response, and alter small molecule transport and metabolism dynamics across tissues.

The scales of coevolution: comparative phylogeography and genetic demography of a locally adapted venomous predator and its prey

Coevolutionary theory predicts that differences in the genetic demography of interacting species can influence patterns of local adaptation by affecting the potential of local populations to respond to selection. We conducted a comparative phylogeographical study of venomous rattlesnakes and their venom-resistant ground squirrel prey across California, and assessed how effective population size (Ne) estimates correspond with a previously documented pattern of rattlesnake local adaptation. Using RAD sequencing markers, we detected lineage relationships among both the rattlesnakes (Crotalus oreganus ssp.) and ground squirrels (Otospermophilus sp.) that are incongruent with previous phylogenetic hypotheses. Both rattlesnakes and squirrels share a deep divergence at the Sacramento–San Joaquin River Delta. At this broad phylogeographical scale, we found that the locally adapted rattlesnakes had higher Ne than squirrels. At the population scale, snakes also had larger Ne accompanied by larger values of several metrics of population genetic diversity. However, the specific magnitude of local adaptation of venom activity to ground squirrel venom resistance was not significantly correlated with local differences in Ne or other diversity statistics between predator and prey populations, suggesting that other factors in the geographic mosaic of coevolution contribute to the specific local-scale outcomes of this interaction. These results suggest an evolutionary mechanism that may explain some (but clearly not all) of rattlesnake local adaptation in this coevolutionary interaction – larger population sizes raise the adaptive potential of rattlesnakes compared to ground squirrels.

Physiological Stress Integrates Resistance to Rattlesnake Venom and the Onset of Risky Foraging in California Ground Squirrels

Using venom for predation often leads to the evolution of resistance in prey. Understanding individual variation in venom resistance is key to unlocking basic mechanisms by which antagonistic coevolution can sustain variation in traits under selection. For prey, the opposing challenges of predator avoidance and resource acquisition often lead to correlated levels of risk and reward, which in turn can favor suites of integrated morphological, physiological and behavioral traits. We investigate the relationship between risk-sensitive behaviors, physiological resistance to rattlesnake venom, and stress in a population of California ground squirrels. For the same individuals, we quantified foraging decisions in the presence of snake predators, fecal corticosterone metabolites (a measure of “stress”), and blood serum inhibition of venom enzymatic activity (a measure of venom resistance). Individual responses to snakes were repeatable for three measures of risk-sensitive behavior, indicating that some individuals were consistently risk-averse whereas others were risk tolerant. Venom resistance was lower in squirrels with higher glucocorticoid levels and poorer body condition. Whereas resistance failed to predict proximity to and interactions with snake predators, individuals with higher glucocorticoid levels and in lower body condition waited the longest to feed when near a snake. We compared alternative structural equation models to evaluate alternative hypotheses for the relationships among stress, venom resistance, and behavior. We found support for stress as a shared physiological correlate that independently lowers venom resistance and leads to squirrels that wait longer to feed in the presence of a snake, whereas we did not find evidence that resistance directly facilitates latency to forage. Our findings suggest that stress may help less-resistant squirrels avoid a deadly snakebite, but also reduces feeding opportunities. The combined lethal and non-lethal effects of stressors in predator–prey interactions simultaneously impact multiple key traits in this system, making environmental stress a potential contributor to geographic variation in trait expression of toxic predators and resistant prey.