Is Hybridization a Source of Adaptive Venom Variation in Rattlesnakes? A Test, Using a Crotalus scutulatus × viridis Hybrid Zone in Southwestern New Mexico

A Comparative Analysis of the Cytotoxic and Vascular Activity Effects of Western Diamondback Rattlesnake (Crotalus atrox) and Eastern Diamondback Rattlesnake (Crotalus adamanteus) Venoms Using a Chick Embryo Model

Crotalus snakebites induce various toxicological effects, encompassing neurological, myotoxic, and cytotoxic symptoms, with potentially fatal outcomes. Investigating venom toxicity is essential for public health, and developing new tools allows for these effects to be studied more comprehensively. The research goals include the elucidation of the physiological consequences of venom exposure and the assessment of toxicity using animal models. Chicken embryos serve as valuable models for assessing venom toxicity through the chick embryotoxicity screening test (CHEST) and the chick chorioallantoic membrane (CAM) assay, particularly useful for evaluating vascular impacts. C. adamanteus venom application resulted in higher embryotoxicity and morphological abnormalities, such as Siamese twins. The CAM assay demonstrated the hemorrhagic effects of venom, varying with venom type and concentration. The irritant potential of both venom types was classified as slight or moderate depending on their concentration. Additionally, acetylcholinesterase (AChE) activity was performed to receive information about organ toxicity. The results show that both venoms induced changes in the whole embryo, heart, and liver weights, but the C. adamanteus venom was identified as more toxic. Specific venom concentrations affected AChE activity in embryonic tissues. These findings underscore the embryotoxic and vasoactive properties of Crotalus venoms, providing valuable insights into their mechanisms of toxicity and potential applications in biomedicine.

A Review of Rattlesnake Venoms

Venom components are invaluable in biomedical research owing to their specificity and potency. Many of these components exist in two genera of rattlesnakes, Crotalus and Sistrurus, with high toxicity and proteolytic activity variation. This review focuses on venom components within rattlesnakes, and offers a comparison and itemized list of factors dictating venom composition, as well as presenting their known characteristics, activities, and significant applications in biosciences. There are 64 families and subfamilies of proteins present in Crotalus and Sistrurus venom. Snake venom serine proteases (SVSP), snake venom metalloproteases (SVMP), and phospholipases A2 (PLA2) are the standard components in Crotalus and Sistrurus venom. Through this review, we highlight gaps in the knowledge of rattlesnake venom; there needs to be more information on the venom composition of three Crotalus species and one Sistrurus subspecies. We discuss the activity and importance of both major and minor components in biomedical research and drug development.

Hunting behavior and feeding ecology of Mojave rattlesnakes (Crotalus scutulatus), prairie rattlesnakes (Crotalus viridis), and their hybrids in southwestern New Mexico

Predators must contend with numerous challenges to successfully find and subjugate prey. Complex traits related to hunting are partially controlled by a large number of co-evolved genes, which may be disrupted in hybrids. Accordingly, research on the feeding ecology of animals in hybrid zones has shown that hybrids sometimes exhibit transgressive or novel behaviors, yet for many taxa, empirical studies of predation and diet across hybrid zones are lacking. We undertook the first such field study for a hybrid zone between two snake species, the Mojave rattlesnake (Crotalus scutulatus) and the prairie rattlesnake (Crotalus viridis). Specifically, we leveraged established field methods to quantify the hunting behaviors of animals, their prey communities, and the diet of individuals across the hybrid zone in southwestern New Mexico, USA. We found that, even though hybrids had significantly lower body condition indices than snakes from either parental group, hybrids were generally similar to non-hybrids in hunting behavior, prey encounter rates, and predatory attack and success. We also found that, compared to C. scutulatus, C. viridis was significantly more active while hunting at night and abandoned ambush sites earlier in the morning, and hybrids tended to be more viridis-like in this respect. Prey availability was similar across the study sites, including within the hybrid zone, with kangaroo rats (Dipodomys spp.) as the most common small mammal, both in habitat surveys and the frequency of encounters with hunting rattlesnakes. Analysis of prey remains in stomachs and feces also showed broad similarity in diets, with all snakes preying primarily on small mammals and secondarily on lizards. Taken together, our results suggest that the significantly lower body condition of hybrids does not appear to be driven by differences in their hunting behavior or diet and may instead relate to metabolic efficiency or other physiological traits we have not yet identified.

The best of both worlds? Rattlesnake hybrid zones generate complex combinations of divergent venom phenotypes that retain high toxicity

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

Evidence that genomic incompatibilities and other multilocus processes impact hybrid fitness in a rattlesnake hybrid zone

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

Mohave Rattlesnake (Crotalus scutulatus) Identification Revisited

Crotalus scutulatus (Mohave rattlesnake) is a clinically significant pit viper broadly distributed across much of the arid southwestern United States and mainland Mexico. Identification of C scutulatus is a concern among emergency medical service and emergency department personnel owing to its reputation for severe envenomations and difficulty in visually differentiating between C scutulatus and other species, primarily Crotalus atrox (western diamond-backed rattlesnake). We contrast distinctive characteristics of C scutulatus, C atrox, and 3 other sympatric species: Crotalus molossus, Crotalus ornatus, and Crotalus viridis (western and eastern black-tailed rattlesnakes and prairie rattlesnake, respectively). Greenish coloration eliminates C atrox but does not confirm C scutulatus. Obvious coarse and fine speckling of the dorsal pattern and a pale postocular stripe intersecting the mouth characterize C atrox. Dorsal speckling is insignificant or absent in the other species, whereas the pale postocular stripe passes above the mouth in C scutulatus and C viridis and is absent in C molossus and C ornatus. Tails boldly ringed with alternating black and white or contrasting shades of gray are shared by C atrox and C scutulatus, respectively, but a lack of boldly ringed tails characterizes the other species. The proximal rattle segment is yellow and black, or entirely yellow, in C scutulatus but black in the others. The most reliable visual identifications are based on evaluations of multiple traits, all of which are variable to some extent. Traits such as tail ring width and the size and number of crown scales have frequently been overemphasized in the past.

Hybridization between Crotalus aquilus and Crotalus polystictus Species: A Comparison of Their Venom Toxicity and Enzymatic Activities

Hybridization is defined as the interbreeding of individuals from two populations distinguishable by one or more heritable characteristics. Snake hybridization represents an interesting opportunity to analyze variability and how genetics affect the venom components between parents and hybrids. Snake venoms exhibit a high degree of variability related to biological and biogeographical factors. The aim of this work is to analyze the protein patterns and enzymatic activity of some of the main hemotoxic enzymes in snake venoms, such as serine proteases (trypsin-like, chymotrypsin-like, and elastase-like), metalloproteases, hyaluronidases, and phospholipase A₂. The lethal dose of 50 (LD₅₀) of venom from the Crotalus aquilus (Cabf) and Crotalus polystictus (Cpbm) parents and their hybrids in captivity was determined, and phenetic analysis is also conducted, which showed a high similarity between the hybrids and C. polystictus. The protein banding patterns and enzymatic activity analyze by zymography resulted in a combination of proteins from the parental venoms in the hybrids, with variability among them. In some cases, the enzymatic activity is higher in the hybrids with a lower LD₅₀ than in the parents, indicating higher toxicity. These data show the variability among snake venoms and suggest that hybridization is an important factor in changes in protein concentration, peptide variability, and enzymatic activity that affect toxicity and lethality.

Characterization of the venom and external morphology of a natural hybrid between Crotalus atrox and Crotalus mictlantecuhtli

Here we report, for the first time, a natural hybrid between Crotalus atrox and C. mictlantecuhtli based on intermediate characteristics of the external morphology and venom. Morphologically, the individual had characteristics of both parent species. The hybrid's venom exhibited an intermediate composition including the presence of crotoxin which has never been documented in C. atrox but is well documented in C. mictlantecuhtli. The hybrid's venom was highly toxic and showed an intermediate proteolytic activity between the parental species. The two Mexican antivenoms were able to neutralize the hybrid's venom's lethality.

Varying Intensities of Introgression Obscure Incipient Venom-Associated Speciation in the Timber Rattlesnake (Crotalus horridus)

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

Genome-wide data reveal extensive gene flow during the diversification of the western rattlesnakes (Viperidae: Crotalinae: Crotalus)

Hybridization and introgression are important, but often overlooked processes when inferring phylogenies. When these processes are not accounted for and a strictly diverging phylogenetic model is applied to groups with a history of hybridization, phylogenetic inference and parameter estimation can be inaccurate. Recent developments in phylogenetic network approaches coupled with the increasing availability of genomic data allow inferences of reticulate evolutionary histories across the tree of life. The western rattlesnake species group (C. viridis species complex, C. mitchellii species complex, C. scutulutas, and C. tigris) is an iconic snake lineage that is widespread across western North America. This group is composed of several species complexes with unclear species limits, likely the result of ongoing gene flow among nascent lineages. Here I generate reduced representation genomic data and test for a history of reticulation within this group. I demonstrate that all species have undergone hybridization with at least one other lineage, suggesting introgression is widespread in this group. Topologies differ between phylogenies estimated under the multispecies coalescent and multispecies network coalescent methods, indicating that gene flow has obscured phylogenetic relationships within this group. These past introgression events are predominantly restricted to species that co-occur geographically. However, within species that have a history of introgression, this signature is detected regardless of specimen sampling across geography. Overall, my results suggest the accumulation of reproductive isolating barriers occurs slowly in rattlesnakes which likely leads to the difficulty in delimiting species, furthermore, the results of this study have implications for trait evolution in this group.

Biological and Biochemical Characterization of Coronado Island Rattlesnake (Crotalus helleri caliginis) Venom and Antivenom Neutralization

The Baja California Peninsula has over 250 islands and islets with many endemic species. Among them, rattlesnakes are the most numerous but also one of the least studied groups. The study of island rattlesnake venom could guide us to a better understanding of evolutionary processes and the description of novel toxins. Crotalus helleri caliginis venom samples were analyzed to determine possible ontogenetic variation with SDS-PAGE in one and two dimensions and with RP-HPLC. Western Blot, ELISA, and amino-terminal sequencing were used to determine the main components of the venom. The biological and biochemical activities demonstrate the similarity of C. helleri caliginis venom to the continental species C. helleri helleri, with both having low proteolytic and phospholipase A₂ (PLA₂) activity but differing due to the absence of neurotoxin (crotoxin-like) in the insular species. The main components of the snake venom were metalloproteases, serine proteases, and crotamine, which was the most abundant toxin group (30–35% of full venom). The crotamine was isolated using size-exclusion chromatography where its functional effects were tested on mouse phrenic nerve–hemidiaphragm preparations in which a significant reduction in muscle twitch contractions were observed. The two Mexican antivenoms could neutralize the lethality of C. helleri caliginis venom but not the crotamine effects.

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.

A Meta-Analysis of the Protein Components in Rattlesnake Venom

The specificity and potency of venom components give them a unique advantage in developing various pharmaceutical drugs. Though venom is a cocktail of proteins, rarely are the synergy and association between various venom components studied. Understanding the relationship between various components of venom is critical in medical research. Using meta-analysis, we observed underlying patterns and associations in the appearance of the toxin families. For Crotalus, Dis has the most associations with the following toxins: PDE; BPP; CRL; CRiSP; LAAO; SVMP P-I and LAAO; SVMP P-III and LAAO. In Sistrurus venom, CTL and NGF have the most associations. These associations can predict the presence of proteins in novel venom and understand synergies between venom components for enhanced bioactivity. Using this approach, the need to revisit the classification of proteins as major components or minor components is highlighted. The revised classification of venom components is based on ubiquity, bioactivity, the number of associations, and synergies. The revised classification can be expected to trigger increased research on venom components, such as NGF, which have high biomedical significance. Using hierarchical clustering, we observed that the genera's venom compositions were similar, based on functional characteristics rather than phylogenetic relationships.

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.

When one phenotype is not enough: divergent evolutionary trajectories govern venom variation in a widespread rattlesnake species

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

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.

In vitro effects of Crotalus atrox snake venom on chick and mouse neuromuscular preparations

The neuromuscular effect of venoms is not a major clinical manifestation shared between rattlesnakes native to the Americas, which showed two different venom phenotypes. Taking into account this dichotomy, nerve muscle preparations from mice and chicks were used to investigate the ability of Crotalus atrox venom to induce in vitro neurotoxicity and myotoxicity. Unlike crotalic venoms of South America, low concentrations of C. atrox venom did not result in significant effects on mouse neuromuscular preparations. The venom was more active on avian nerve-muscle, showing reduction of twitch heights after 120 min of incubation with 10, 30 and 100 μg/mL of venom with diminished responses to agonists and KCl. Histological analysis highlighted that C. atrox was myotoxic in both species of experimental animals; as evidenced by degenerative events, including edematous cells, delta lesions, hypercontracted fibers and muscle necrosis, which can lead to neurotoxic action. These results provide key insights into the myotoxicity and low neurotoxicity of C. atrox in two animal models, corroborating with previous genomic and proteomic findings and would be useful for a deeper understanding of venom evolution in snakes belonging to the genus Crotalus.

Proteopeptidomic, Functional and Immunoreactivity Characterization of Bothrops moojeni Snake Venom: Influence of Snake Gender on Venom Composition

Venom composition varies across snakes from all taxonomic levels and is influenced by the snakes’ age, habitat, diet, and sexual dimorphism. The present study reports the first in-depth investigation of venom composition in male and female Bothrops moojeni (B. moojeni) snakes (BmooM and BmooF, respectively) through three proteomics approaches associated with functional, cytotoxic, and immunoreactivity characterization. Compared with BmooM venom, BmooF venom exhibited weaker hyaluronidase, metalloproteinase, and phospholipase activity; stronger recognition by anti-bothropic serum; 1.4-fold stronger cytotoxicity; and greater number of peptides. The increased L-amino acid oxidase expression probably accounted for the stronger immunoreactivity and cytotoxicity of BmooF venom. BmooF and BmooM venom shared only 19% peptides. Some venom components were gender-specific, such as phospholipases B, phospholipase inhibitor, and hyaluronidases in BmooM, and cysteine-rich secretory proteins in BmooF. In conclusion, we describe herein the first proteomics study of B. moojeni snake venom and an in-depth characterization of gender-specific differences in venom composition. Altogether, our findings not only stress the importance of considering the snake’s gender during antivenom production, but also help to identify new potential drugs and biotechnological tools.

Toxin-resolved antivenomics-guided assessment of the immunorecognition landscape of antivenoms

Snakebite envenoming represents a major issue in rural areas of tropical and subtropical regions across sub-Saharan Africa, South to Southeast Asia, Latin America and Oceania. Antivenoms constitute the only scientifically validated therapy for snakebite envenomings, provided they are safe, effective, affordable, accessible and administered appropriately. However, the lack of financial incentives in a technology that has remained relatively unchanged for more than a century, has contributed to some manufacturers leaving the market and others downscaling production or increasing the prices, leading to a decline in the availability and accessibility for these life-saving antidotes to millions of rural poor most at risk from snakebites in low income countries. The shortage of antivenoms can be significantly alleviated by optimizing the use of current antivenoms (through the assessment of their specific and paraspecific efficacy against the different medically relevant homologous and heterologous snake venoms) and by generating novel polyspecific antivenoms exhibiting broad clinical spectrum and wide geographic distribution range. Research on venoms has been continuously enhanced by advances in technology. Particularly, the last decade has witnessed the development of omics strategies for unravelling the toxin composition of venoms ("venomics") and to assess the immunorecognition profile of antivenoms ("antivenomics"). Here, we review recent developments and reflect on near future innovations that promise to revolutionize the mutually enlightening relationship between evolutionary and translational venomics.

Phenotypic Variation in Mojave Rattlesnake (Crotalus scutulatus) Venom Is Driven by Four Toxin Families

Phenotypic diversity generated through altered gene expression is a primary mechanism facilitating evolutionary response in natural systems. By linking the phenotype to genotype through transcriptomics, it is possible to determine what changes are occurring at the molecular level. High phenotypic diversity has been documented in rattlesnake venom, which is under strong selection due to its role in prey acquisition and defense. Rattlesnake venom can be characterized by the presence (Type A) or absence (Type B) of a type of neurotoxic phospholipase A2 (PLA2), such as Mojave toxin, that increases venom toxicity. Mojave rattlesnakes (Crotalus scutulatus), represent this diversity as both venom types are found within this species and within a single panmictic population in the Sonoran Desert. We used comparative venom gland transcriptomics of nine specimens of C. scutulatus from this region to test whether expression differences explain diversity within and between venom types. Type A individuals expressed significantly fewer toxins than Type B individuals owing to the diversity of C-type lectins (CTLs) and snake venom metalloproteinases (SVMPs) found in Type B animals. As expected, both subunits of Mojave toxin were exclusively found in Type A individuals but we found high diversity in four additional PLA2s that was not associated with a venom type. Myotoxin a expression and toxin number variation was not associated with venom type, and myotoxin a had the highest range of expression of any toxin class. Our study represents the most comprehensive transcriptomic profile of the venom type dichotomy in rattlesnakes and C. scutulatus. Even intra-specifically, Mojave rattlesnakes showcase the diversity of snake venoms and illustrate that variation within venom types blurs the distinction of the venom dichotomy.

Extremely Divergent Haplotypes in Two Toxin Gene Complexes Encode Alternative Venom Types within Rattlesnake Species

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

Biological and Proteolytic Variation in the Venom of Crotalus scutulatus scutulatus from Mexico

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

Venom On-a-Chip: A Fast and Efficient Method for Comparative Venomics

Venom research has attracted an increasing interest in disparate fields, from drug development and pharmacology, to evolutionary biology and ecology, and rational antivenom production. Advances in “-omics” technologies have allowed the characterization of an increasing number of animal venoms, but the methodology currently available is suboptimal for large-scale comparisons of venom profiles. Here, we describe a fast, reproducible and semi-automated protocol for investigating snake venom variability, especially at the intraspecific level, using the Agilent Bioanalyzer on-chip technology. Our protocol generated a phenotype matrix which can be used for robust statistical analysis and correlations of venom variation with ecological correlates, or other extrinsic factors. We also demonstrate the ease and utility of combining on-chip technology with previously fractionated venoms for detection of specific individual toxin proteins. Our study describes a novel strategy for rapid venom discrimination and analysis of compositional variation at multiple taxonomic levels, allowing researchers to tackle evolutionary questions and unveiling the drivers of the incredible biodiversity of 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.