No safety in the trees: Local and species-level adaptation of an arboreal squirrel to the venom of sympatric rattlesnakes

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.

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.

Widespread convergent evolution of alpha-neurotoxin resistance in African mammals

Convergent evolution is central to the study of adaptation and has been used to understand both the limits of evolution and the diverse patterns and processes which result in adaptive change. Resistance to snake venom alpha-neurotoxins (αNTXs) is a case of widespread convergence having evolved several times in snakes, lizards and mammals. Despite extreme toxicity of αNTXs, substitutions in its target, the nicotinic acetylcholine receptor (nAChR), prevent αNTX binding and render species resistant. Recently, the published meerkat (Herpestidae) genome revealed that meerkats have the same substitutions in nAChR as the venom-resistant Egyptian mongoose (Herpestidae), suggesting that venom-resistant nAChRs may be ancestral to Herpestids. Like the mongoose, many other species of feliform carnivores prey on venomous snakes, though their venom resistance has never been explored. To evaluate the prevalence and ancestry of αNTX resistance in mammals, we generate a dataset of mammalian nAChR using museum specimens and public datasets. We find five instances of convergent evolution within feliform carnivores, and an additional eight instances across all mammals sampled. Tests of selection show that these substitutions are evolving under positive selection. Repeated convergence suggests that this adaptation played an important role in the evolution of mammalian physiology and potentially venom evolution.

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

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.

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.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

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.

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.

Evidence for Snake Venom Plasticity in a Long-Term Study with Individual Captive Bothrops atrox

Variability in snake venom composition has been frequently reported and correlated to the adaptability of snakes to environmental conditions. Previous studies report plasticity for the venom phenotype. However, these observations are not conclusive, as the results were based on pooled venoms, which present high individual variability. Here we tested the hypothesis of plasticity by influence of confinement and single diet type in the venom composition of 13 adult specimens of Bothrops atrox snakes, maintained under captivity for more than three years. Individual variability in venom composition was observed in samples extracted just after the capture of the snakes. However, composition was conserved in venoms periodically extracted from nine specimens, which presented low variability restricted to the less abundant components. In a second group, composed of four snakes, drastic changes were observed in the venom samples extracted at different periods, mostly related to snake venom metalloproteinases (SVMPs), the core function toxins of B. atrox venom, which occurred approximately between 400 and 500 days in captivity. These data show plasticity in the venom phenotype during the lifetime of adult snakes maintained under captive conditions. Causes or functional consequences involved in the phenotype modification require further investigations.

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.

Snake venom potency and yield are associated with prey-evolution, predator metabolism and habitat structure

Snake venom is well known for its ability to incapacitate and kill prey. Yet, potency and the amount of venom available varies greatly across species, ranging from the seemingly harmless to those capable of killing vast numbers of potential prey. This variation is poorly understood, with comparative approaches confounded by the use of atypical prey species as models to measure venom potency. Here, we account for such confounding issues by incorporating the phylogenetic similarity between a snake's diet and the species used to measure its potency. In a comparative analysis of 102 species we show that snake venom potency is generally prey-specific. We also show that venom yields are lower in species occupying three dimensional environments and increases with body size corresponding to metabolic rate, but faster than predicted from increases in prey size. These results underline the importance of physiological and environmental factors in the evolution of predator traits.

Evidence for divergent patterns of local selection driving venom variation in Mojave Rattlesnakes (Crotalus scutulatus)

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

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.

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.