Mojave rattlesnakes (Crotalus scutulatus scutulatus) lacking the acidic subunit DNA sequence lack Mojave toxin in their venom

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.

Neurologic and hematologic sequalae following a timber rattlesnake (Crotalus horridus) envenomation in a dachshund

A 2-year-old female Dachshund had a witnessed timber rattlesnake envenomation. Although rattlesnake envenomations are a common, potentially life-threatening event in companion animals, timber rattlesnake envenomations in the dog are rarely reported. This dog described in this case report had significant hematologic and neurologic clinical derangements consistent with Types A and B rattlesnake venom and a suspected hypersensitivity reaction to the venom. This patient was treated aggressively with antivenom and fully recovered without any persistent neurologic signs at follow-up.

Severe, persistent thrombocytopenia in Crotalus horridus envenomation despite antivenom: A retrospective review

BACKGROUND

In the late nineties, Bond and Burkhardt described a severe thrombocytopenic phenomenon from envenomation by Crotalus horridus. This thrombocytopenia persisted despite administration of platelets and antivenom. Questions remain regarding the clinical significance and time to resolution of this thrombocytopenia. In addition, as new antivenoms are available in North America, the response to current treatment is not well reported. The purpose of this study is to provide further insight into the approach to treatment of Crotalus horridus envenomation.

METHODS

This is a retrospective chart review of 21 cases of presumed envenomation by C. horridus. Data collected included age, sex, antivenom administration, laboratory data, length of hospital stay, blood products administered, and general clinical course. We also evaluated platelet response to antivenom, bleeding outcomes, and complications from envenomation.

RESULTS

Patients' ages ranged from 19 to 71 years. All patients were men. Most patients presented with thrombocytopenia and all had limb swelling. Patients responded initially to antivenom treatment, however subsequently developed a profound thrombocytopenia, including fourteen with platelet counts less than 20 × 10⁹/L. Abnormalities in thromboelastography (TEG) were noted in conjunction with thrombocytopenia. Patients displayed persistent thrombocytopenia despite administration of Crotalidae polyvalent immune Fab or Crotalidae immune F(ab')2. Median time to rebound platelet count greater than 20 × 10⁹/L was ten days (range 6-12 days) from envenomation. Complications included a partial finger amputation in one patient, bleeding gums in four patients, bloody stools in two patients, bloody nasogastric output in one patient. No patients required red blood cell transfusion and no deaths occurred.

CONCLUSION

Practitioners treating C. horridus should recognize the possibility of severe thrombocytopenia and its persistence despite antivenom. They should counsel patients on appropriate abstention from activities that could lead to trauma, as well as the importance of follow up for repeat laboratory studies to ensure the resolution of thrombocytopenia.

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.

A Single Mutation Unlocks Cascading Exaptations in the Origin of a Potent Pitviper Neurotoxin

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

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.

Ontogenetic Change in the Venom of Mexican Black-Tailed Rattlesnakes (Crotalus molossus nigrescens)

Ontogenetic changes in venom composition have important ecological implications due the relevance of venom in prey acquisition and defense. Additionally, intraspecific venom variation has direct medical consequences for the treatment of snakebite. However, ontogenetic changes are not well documented in most species. The Mexican Black-tailed Rattlesnake (Crotalus molossus nigrescens) is large-bodied and broadly distributed in Mexico. To document venom variation and test for ontogenetic changes in venom composition, we obtained venom samples from twenty-seven C. m. nigrescens with different total body lengths (TBL) from eight states in Mexico. The primary components in the venom were detected by reverse-phase HPLC, western blot, and mass spectrometry. In addition, we evaluated the biochemical (proteolytic, coagulant and fibrinogenolytic activities) and biological (LD₅₀ and hemorrhagic activity) activities of the venoms. Finally, we tested for recognition and neutralization of Mexican antivenoms against venoms of juvenile and adult snakes. We detected clear ontogenetic venom variation in C. m. nigrescens. Venoms from younger snakes contained more crotamine-like myotoxins and snake venom serine proteinases than venoms from older snakes; however, an increase of snake venom metalloproteinases was detected in venoms of larger snakes. Venoms from juvenile snakes were, in general, more toxic and procoagulant than venoms from adults; however, adult venoms were more proteolytic. Most of the venoms analyzed were hemorrhagic. Importantly, Mexican antivenoms had difficulties recognizing low molecular mass proteins (<12 kDa) of venoms from both juvenile and adult snakes. The antivenoms did not neutralize the crotamine effect caused by the venom of juveniles. Thus, we suggest that Mexican antivenoms would have difficulty neutralizing some human envenomations and, therefore, it may be necessary improve the immunization mixture in Mexican antivenoms to account for low molecular mass proteins, like myotoxins.

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.

Rattling the border wall: Pathophysiological implications of functional and proteomic venom variation between Mexican and US subspecies of the desert rattlesnake Crotalus scutulatus

While some US populations of the Mohave rattlesnake (Crotalus scutulatus scutulatus) are infamous for being potently neurotoxic, the Mexican subspecies C. s. salvini (Huamantlan rattlesnake) has been largely unstudied beyond crude lethality testing upon mice. In this study we show that at least some populations of this snake are as potently neurotoxic as its northern cousin. Testing of the Mexican antivenom Antivipmyn showed a complete lack of neutralisation for the neurotoxic effects of C. s. salvini venom, while the neurotoxic effects of the US subspecies C. s. scutulatus were time-delayed but ultimately not eliminated. These results document unrecognised potent neurological effects of a Mexican snake and highlight the medical importance of this subspecies, a finding augmented by the ineffectiveness of the Antivipmyn antivenom. These results also influence our understanding of the venom evolution of Crotalus scutulatus, suggesting that neurotoxicity is the ancestral feature of this species, with the US populations which lack neurotoxicity being derived states.

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.

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

Venomous snakes often display extensive variation in venom composition both between and within species. However, the mechanisms underlying the distribution of different toxins and venom types among populations and taxa remain insufficiently known. Rattlesnakes (Crotalus, Sistrurus) display extreme inter- and intraspecific variation in venom composition, centered particularly on the presence or absence of presynaptically neurotoxic phospholipases A₂ such as Mojave toxin (MTX). Interspecific hybridization has been invoked as a mechanism to explain the distribution of these toxins across rattlesnakes, with the implicit assumption that they are adaptively advantageous. Here, we test the potential of adaptive hybridization as a mechanism for venom evolution by assessing the distribution of genes encoding the acidic and basic subunits of Mojave toxin across a hybrid zone between MTX-positive Crotalus scutulatus and MTX-negative C. viridis in southwestern New Mexico, USA. Analyses of morphology, mitochondrial and single copy-nuclear genes document extensive admixture within a narrow hybrid zone. The genes encoding the two MTX subunits are strictly linked, and found in most hybrids and backcrossed individuals, but not in C. viridis away from the hybrid zone. Presence of the genes is invariably associated with presence of the corresponding toxin in the venom. We conclude that introgression of highly lethal neurotoxins through hybridization is not necessarily favored by natural selection in rattlesnakes, and that even extensive hybridization may not lead to introgression of these genes into another species.

Full-Length Venom Protein cDNA Sequences from Venom-Derived mRNA: Exploring Compositional Variation and Adaptive Multigene Evolution

Envenomation of humans by snakes is a complex and continuously evolving medical emergency, and treatment is made that much more difficult by the diverse biochemical composition of many venoms. Venomous snakes and their venoms also provide models for the study of molecular evolutionary processes leading to adaptation and genotype-phenotype relationships. To compare venom complexity and protein sequences, venom gland transcriptomes are assembled, which usually requires the sacrifice of snakes for tissue. However, toxin transcripts are also present in venoms, offering the possibility of obtaining cDNA sequences directly from venom. This study provides evidence that unknown full-length venom protein transcripts can be obtained from the venoms of multiple species from all major venomous snake families. These unknown venom protein cDNAs are obtained by the use of primers designed from conserved signal peptide sequences within each venom protein superfamily. This technique was used to assemble a partial venom gland transcriptome for the Middle American Rattlesnake (Crotalus simus tzabcan) by amplifying sequences for phospholipases A₂, serine proteases, C-lectins, and metalloproteinases from within venom. Phospholipase A₂ sequences were also recovered from the venoms of several rattlesnakes and an elapid snake (Pseudechis porphyriacus), and three-finger toxin sequences were recovered from multiple rear-fanged snake species, demonstrating that the three major clades of advanced snakes (Elapidae, Viperidae, Colubridae) have stable mRNA present in their venoms. These cDNA sequences from venom were then used to explore potential activities derived from protein sequence similarities and evolutionary histories within these large multigene superfamilies. Venom-derived sequences can also be used to aid in characterizing venoms that lack proteomic profiles and identify sequence characteristics indicating specific envenomation profiles. This approach, requiring only venom, provides access to cDNA sequences in the absence of living specimens, even from commercial venom sources, to evaluate important regional differences in venom composition and to study snake venom protein evolution.

This work demonstrates that full-length venom protein messenger RNAs are present in secreted venoms and can be used to acquire full-length protein sequences of toxins from both front-fanged (Elapidae, Viperidae) and rear-fanged (Colubridae) snake venoms, eliminating the need to use venom glands. Full-length transcripts were obtained from venom samples that were fresh, newly lyophilized, old, field desiccated or commercially prepared, representing a significant advance over previous attempts which produced only partial sequence transcripts. Transcripts for all major venom protein families (metalloproteinases, serine proteases, C-type lectins, phospholipases A₂ and three-finger toxins) responsible for clinically significant snakebite symptoms were obtained from venoms. These sequences aid in the identification and characterization of venom proteome profiles, allowing for the identification of peptide sequences, specific isoforms, and novel venom proteins. The application of this technique will help to provide venom protein sequences for many snake species, including understudied rear-fanged snakes. Venom protein transcripts offer important insights into potential snakebite envenomation profiles and the molecular evolution of venom protein multigene families. By requiring only venom to obtain venom protein cDNAs, the approach detailed here will provide access to cDNA-based protein sequences from commercial and other venom sources, facilitating study of snake venom protein composition and evolution.

Combined venomics, venom gland transcriptomics, bioactivities, and antivenomics of two Bothrops jararaca populations from geographic isolated regions within the Brazilian Atlantic rainforest

Bothrops jararaca is a slender and semi-arboreal medically relevant pit viper species endemic to tropical and subtropical forests in southern Brazil, Paraguay, and northern Argentina (Misiones). Within its geographic range, it is often abundant and is an important cause of snakebite. Although no subspecies are currently recognized, geographic analyses have revealed the existence of two well-supported B. jararaca clades that diverged during the Pliocene ~3.8Mya and currently display a southeastern (SE) and a southern (S) Atlantic rainforest (Mata Atlântica) distribution. The spectrum, geographic variability, and ontogenetic changes of the venom proteomes of snakes from these two B. jararaca phylogroups were investigated applying a combined venom gland transcriptomic and venomic analysis. Comparisons of the venom proteomes and transcriptomes of B. jararaca from the SE and S geographic regions revealed notable interpopulational variability that may be due to the different levels of population-specific transcriptional regulation, including, in the case of the southern population, a marked ontogenetic venom compositional change involving the upregulation of the myotoxic PLA2 homolog, bothropstoxin-I. This population-specific marker can be used to estimate the proportion of venom from the southern population present in the B. jararaca venom pool used for the Brazilian soro antibotrópico (SAB) antivenom production. On the other hand, the southeastern population-specific D49-PLA2 molecules, BinTX-I and BinTX-II, lend support to the notion that the mainland ancestor of Bothrops insularis was originated within the same population that gave rise to the current SE B. jararaca phylogroup, and that this insular species endemic to Queimada Grande Island (Brazil) expresses a pedomorphic venom phenotype. Mirroring their compositional divergence, the two geographic B. jararaca venom pools showed distinct bioactivity profiles. However, the SAB antivenom manufactured in Vital Brazil Institute neutralized the lethal effect of both venoms to a similar extent. In addition, immobilized SAB antivenom immunocaptured most of the venom components of the venoms of both B. jararaca populations, but did not show immunoreactivity against vasoactive peptides. The Costa Rican bothropic-crotalic-lachesic (BCL) antivenom showed the same lack of reactivity against vasoactive peptides but, in addition, was less efficient immunocapturing PI- and PIII-SVMPs from the SE venom, and bothropstoxin-I, a CRISP molecule, and a D49-PLA2 from the venom of the southern B. jararaca phylogroup. The remarkable paraspecificity exhibited by the Brazilian and the Costa Rican antivenoms indicates large immunoreactive epitope conservation across the natural history of Bothrops, a genus that has its roots in the middle Miocene. This article is part of a Special Issue entitled: Omics Evolutionary Ecolog.

A case of neurotoxicity following envenomation by the Sidewinder rattlesnake, Crotalus cerastes

INTRODUCTION

North American rattlesnake envenomations typically result in local tissue injury and hematologic derangements. Neurotoxicity is uncommon but when present often manifests as fasciculations and paresthesias. Neurotoxicity following Sidewinder (Crotalus cerastes) envenomation has not been previously reported.

CASE REPORT

A 56-year-old man bitten on the right foot developed painful paresthesias, weakness and fasciculations of the right lower extremity, and involuntary muscle contractions of the anterior thigh. Local tissue effects and hemotoxicity never developed. The patient was discharged 5 days after the bite with resolution of fasciculations but continued to have right-sided weakness. The snake was identified as a Sidewinder, C. cerastes, by the patient and two independent herpetologists.

CONCLUSION

This is the first reported case of a Sidewinder rattlesnake envenomation resulting in neurotoxicity.

Extreme venom variation in Middle Eastern vipers: a proteomics comparison of Eristicophis macmahonii, Pseudocerastes fieldi and Pseudocerastes persicus

Venoms of the viperid sister genera Eristicophis and Pseudocerastes are poorly studied despite their anecdotal reputation for producing severe or even lethal envenomations. This is due in part to the remote and politically unstable regions that they occupy. All species contained are sit and wait ambush feeders. Thus, this study examined their venoms through proteomics techniques in order to establish if this feeding ecology, and putatively low levels of gene flow, have resulted in significant variations in venom profile. The techniques indeed revealed extreme venom variation. This has immediate implications as only one antivenom is made (using the venom of Pseudocerastes persicus) yet the proteomic variation suggests that it would be of only limited use for the other species, even the sister species Pseudocerastes fieldi. The high degree of variation however also points toward these species being rich resources for novel compounds which may have use as lead molecules in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results show extreme venom variation between these closely related snakes. These results have direct implications for the treatment of the envenomed patient.

Medically important differences in snake venom composition are dictated by distinct postgenomic mechanisms

Variation in venom composition is a ubiquitous phenomenon in snakes and occurs both interspecifically and intraspecifically. Venom variation can have severe outcomes for snakebite victims by rendering the specific antibodies found in antivenoms ineffective against heterologous toxins found in different venoms. The rapid evolutionary expansion of different toxin-encoding gene families in different snake lineages is widely perceived as the main cause of venom variation. However, this view is simplistic and disregards the understudied influence that processes acting on gene transcription and translation may have on the production of the venom proteome. Here, we assess the venom composition of six related viperid snakes and compare interspecific changes in the number of toxin genes, their transcription in the venom gland, and their translation into proteins secreted in venom. Our results reveal that multiple levels of regulation are responsible for generating variation in venom composition between related snake species. We demonstrate that differential levels of toxin transcription, translation, and their posttranslational modification have a substantial impact upon the resulting venom protein mixture. Notably, these processes act to varying extents on different toxin paralogs found in different snakes and are therefore likely to be as important as ancestral gene duplication events for generating compositionally distinct venom proteomes. Our results suggest that these processes may also contribute to altering the toxicity of snake venoms, and we demonstrate how this variability can undermine the treatment of a neglected tropical disease, snakebite.

Molecular models of the Mojave rattlesnake (Crotalus scutulatus scutulatus) venom metalloproteinases reveal a structural basis for differences in hemorrhagic activities

Rattlesnake venom can differ in composition and in metalloproteinase-associated activities. The molecular basis for this intra-species variation in Crotalus scutulatus scutulatus (Mojave rattlesnake) remains an enigma. To understand the molecular basis for intra-species variation of metalloproteinase-associated activities, we modeled the three-dimensional structures of four metalloproteinases based on the amino acid sequence of four variations of the proteinase domain of the C. s. scutulatus metalloproteinase gene (GP1, GP2, GP3, and GP4). For comparative purposes, we modeled the atrolysin metalloproteinases of C. atrox as well. All molecular models shared the same topology. While the atrolysin metalloproteinase molecular models contained highly conserved substrate binding sites, the Mojave rattlesnake metalloproteinases showed higher structural divergence when superimposed onto each other. The highest structural divergence among the four C. s. scutulatus molecular models was located at the northern cleft wall and the S'1-pocket of the substrate binding site, molecular regions that modulate substrate selectivity. Molecular dynamics and field potential maps for each C. s. scutulatus metalloproteinase model demonstrated that the non-hemorrhagic metalloproteinases (GP2 and GP3) contain highly basic molecular and field potential surfaces while the hemorrhagic metalloproteinases GP1 and atrolysin C showed extensive acidic field potential maps and shallow but less dynamic active site pockets. Hence, differences in the spatial arrangement of the northern cleft wall, the S'1-pocket, and the physico-chemical environment surrounding the catalytic site contribute to differences in metalloproteinase activities in the Mojave rattlesnake. Our results provide a structural basis for variation of metalloproteinase-associated activities in the rattlesnake venom of the Mojave rattlesnake.

Genetic Basis for Variation of Metalloproteinase-Associated Biochemical Activity in Venom of the Mojave Rattlesnake (Crotalus scutulatus scutulatus)

The metalloproteinase composition and biochemical profiles of rattlesnake venom can be highly variable among rattlesnakes of the same species. We have previously shown that the neurotoxic properties of the Mojave rattlesnake (Crotalus scutulatus scutulatus) are associated with the presence of the Mojave toxin A subunit suggesting the existence of a genetic basis for rattlesnake venom composition. In this report, we hypothesized the existence of a genetic basis for intraspecies variation in metalloproteinase-associated biochemical properties of rattlesnake venom of the Mojave rattlesnake. To address this question, we PCR-amplified and compared the genomic DNA nucleotide sequences that code for the mature metalloproteinase domain of fourteen Mojave rattlesnakes captured from different geographical locations across the southwest region of the United States. In addition, the venoms from the same rattlesnakes were tested for their ability to hydrolyze fibrinogen, fibrin, casein, and hide powder azure and for induction of hemorrhage in mice. Overall, based on genomic sequencing and biochemical data, we classified Mojave rattlesnake venom into four distinct groups of metalloproteinases. These findings indicate that differences in nucleotide sequences encoding the mature proteinase domain and noncoding regions contribute to differences in venom metalloproteinase activities among rattlesnakes of the same species.

The genesis of an exceptionally lethal venom in the timber rattlesnake (Crotalus horridus) revealed through comparative venom-gland transcriptomics

Background

Snake venoms generally show sequence and quantitative variation within and between species, but some rattlesnakes have undergone exceptionally rapid, dramatic shifts in the composition, lethality, and pharmacological effects of their venoms. Such shifts have occurred within species, most notably in Mojave (Crotalus scutulatus), South American (C. durissus), and timber (C. horridus) rattlesnakes, resulting in some populations with extremely potent, neurotoxic venoms without the hemorrhagic effects typical of rattlesnake bites.

Results

To better understand the evolutionary changes that resulted in the potent venom of a population of C. horridus from northern Florida, we sequenced the venom-gland transcriptome of an animal from this population for comparison with the previously described transcriptome of the eastern diamondback rattlesnake (C. adamanteus), a congener with a more typical rattlesnake venom. Relative to the toxin transcription of C. adamanteus, which consisted primarily of snake-venom metalloproteinases, C-type lectins, snake-venom serine proteinases, and myotoxin-A, the toxin transcription of C. horridus was far simpler in composition and consisted almost entirely of snake-venom serine proteinases, phospholipases A₂, and bradykinin-potentiating and C-type natriuretic peptides. Crotalus horridus lacked significant expression of the hemorrhagic snake-venom metalloproteinases and C-type lectins. Evolution of shared toxin families involved differential expansion and loss of toxin clades within each species and pronounced differences in the highly expressed toxin paralogs. Toxin genes showed significantly higher rates of nonsynonymous substitution than nontoxin genes. The expression patterns of nontoxin genes were conserved between species, despite the vast differences in toxin expression.

Conclusions

Our results represent the first complete, sequence-based comparison between the venoms of closely related snake species and reveal in unprecedented detail the rapid evolution of snake venoms. We found that the difference in venom properties resulted from major changes in expression levels of toxin gene families, differential gene-family expansion and loss, changes in which paralogs within gene families were expressed at high levels, and higher nonsynonymous substitution rates in the toxin genes relative to nontoxins. These massive alterations in the genetics of the venom phenotype emphasize the evolutionary lability and flexibility of this ecologically critical trait.

Venom variability and envenoming severity outcomes of the Crotalus scutulatus scutulatus (Mojave rattlesnake) from Southern Arizona

Twenty-one Mojave rattlesnakes, Crotalus scutulatus scutulatus (C. s. scutulatus), were collected from Arizona and New Mexico U.S.A. Venom proteome of each specimen was analyzed using reverse-phase HPLC and SDS-PAGE. The toxicity of venoms was analyzed using lethal dose 50 (LD(50)). Health severity outcomes between two Arizona counties U.S.A., Pima and Cochise, were determined by retrospective chart review of the Arizona Poison and Drug Information Center (APDIC) database between the years of 2002 and 2009. Six phenotypes (A-F) were identified based on three venom protein families; Mojave toxin, snake venom metalloproteinases PI and PIII (SVMP), and myotoxin-A. Venom changed geographically from SVMP-rich to Mojave toxin-rich phenotypes as you move from south central to southeastern Arizona. Phenotypes containing myotoxin-A were only found in the transitional zone between the SVMP and Mojave toxin phenotypes. Venom samples containing the largest amounts of SVMP or Mojave toxin had the highest and lowest LD(50s), respectively. There was a significant difference when comparing the presence of neurotoxic effects between Pima and Cochise counties (p=0.001). No significant difference was found when comparing severity (p=0.32), number of antivenom vials administered (p=0.17), days spent in a health care facility (p=0.23) or envenomation per 100,000 population (p=0.06). Although not part of the original data to be collected, death and intubations, were also noted. There is a 10× increased risk of death and a 50× increased risk of intubations if envenomated in Cochise County.

Snake venomics of Crotalus tigris: the minimalist toxin arsenal of the deadliest Nearctic rattlesnake venom. Evolutionary Clues for generating a pan-specific antivenom against crotalid type II venoms [corrected]

We report the proteomic and antivenomic characterization of Crotalus tigris venom. This venom exhibits the highest lethality for mice among rattlesnakes and the simplest toxin proteome reported to date. The venom proteome of C. tigris comprises 7-8 gene products from 6 toxin families; the presynaptic β-neurotoxic heterodimeric PLA(2), Mojave toxin, and two serine proteinases comprise, respectively, 66 and 27% of the C. tigris toxin arsenal, whereas a VEGF-like protein, a CRISP molecule, a medium-sized disintegrin, and 1-2 PIII-SVMPs each represent 0.1-5% of the total venom proteome. This toxin profile really explains the systemic neuro- and myotoxic effects observed in envenomated animals. In addition, we found that venom lethality of C. tigris and other North American rattlesnake type II venoms correlates with the concentration of Mojave toxin A-subunit, supporting the view that the neurotoxic venom phenotype of crotalid type II venoms may be described as a single-allele adaptation. Our data suggest that the evolutionary trend toward neurotoxicity, which has been also reported for the South American rattlesnakes, may have resulted by pedomorphism. The ability of an experimental antivenom to effectively immunodeplete proteins from the type II venoms of C. tigris, Crotalus horridus , Crotalus oreganus helleri, Crotalus scutulatus scutulatus, and Sistrurus catenatus catenatus indicated the feasibility of generating a pan-American anti-Crotalus type II antivenom, suggested by the identification of shared evolutionary trends among South and North American Crotalus species.

Protein complexes in snake venom

Snake venom contains mixture of bioactive proteins and polypeptides. Most of these proteins and polypeptides exist as monomers, but some of them form complexes in the venom. These complexes exhibit much higher levels of pharmacological activity compared to individual components and play an important role in pathophysiological effects during envenomation. They are formed through covalent and/or non-covalent interactions. The subunits of the complexes are either identical (homodimers) or dissimilar (heterodimers; in some cases subunits belong to different families of proteins). The formation of complexes, at times, eliminates the non-specific binding and enhances the binding to the target molecule. On several occasions, it also leads to recognition of new targets as protein-protein interaction in complexes exposes the critical amino acid residues buried in the monomers. Here, we describe the structure and function of various protein complexes of snake venoms and their role in snake venom toxicity.

Complement inactivating proteins and intraspecies venom variation in Crotalus oreganus helleri

Complement inactivating properties were detected in venom from the southern California distribution of Crotalus oreganus helleri (Southern Pacific Rattlesnake). This activity showed strong geographic bias to the San Bernardino Mountain range, and venom from this area reacted strongly with Fraction 5 antiserum (AF5). However, venoms from the San Jacinto Mountain range, which have been previously shown to contain Mojave toxin, did not inhibit complement and did not react with AF5. AF5 has been previously shown to recognize a protease in C. scutulatus venom that inactivates complement, but the identity of this protein has remained unknown. Using a functional venomic approach, utilizing two-dimensional gel electrophoresis coupled with liquid chromatography and tandem mass spectrometry (LC/MS/MS), we have identified catrocollastatin and hemorrhagic toxin II (HT-2) as the primary proteins recognized by AF5. The information we present within this manuscript further illustrates the now well-known reality of intraspecies venom variation and the challenges faced in providing comprehensive polyvalent antivenoms.

Disintegrin, hemorrhagic, and proteolytic activities of Mohave rattlesnake, Crotalus scutulatus scutulatus venoms lacking Mojave toxin

Venom from the Mohave rattlesnake, Crotalus scutulatus scutulatus, has been reported to be either: (1) neurotoxic; (2) hemorrhagic, or both (3) neurotoxic and hemorrhagic. In this study, 14 Mohave rattlesnakes from Arizona and Texas (USA) were analyzed for the presence of disintegrins and Mojave toxin. All venom samples were analyzed for the presence of hemorrhagic, proteolytic and disintegrin activities. The venoms were each chromatographed by reverse phase and their fractions tested for disintegrin activity. All specimens containing Mojave toxin were the most toxic and lacked proteolytic, hemorrhagic and disintegrin activities. In contrast, the venoms containing these activities lacked Mojave toxin. Two disintegrin genes, scutustatin and mojavestatin, were identified by PCR of genomic sequences. Scutustatin is a highly conserved disintegrin, while mojavestatin shows low conservation to other known disintegrins. Venoms with the highest LD50 measurements lacked both disintegrin genes, while the specimens with intermediate and low LD50 contained both genes. The intermediate LD50 group contained Mojave toxin and both disintegrin genes, but lacked hemorrhagic and disintegrin activity. Our results raise the possibility that scutustatin and mojavestatin are not expressed in the intermediate LD50 group, or that they may not be the same disintegrins responsible for the disintegrin activity found in the venom. Therefore, it is possible that Mohave rattlesnakes may produce more than two disintegrins.

Mojave toxin in venom of Crotalus helleri (Southern Pacific Rattlesnake): molecular and geographic characterization

Mojave toxin (MT) was detected in five of 25 Crotalus helleri (Southern Pacific rattlesnake) sampled using anti-MT antibodies and nucleotide sequence analysis. All of the venoms that were positive for MT were collected from Mt San Jacinto in Riverside Co., California. Since this population is geographically isolated from C. scutulatus scutulatus (Mojave rattlesnake), it is unlikely that this finding is due to recent hybridization. MT concentration differences between C. helleri and C. s. scutulatus reflected the presence of 'isoforms' of the toxin in the venom. Whereas C. s. scutulatus generally has several isoforms of the toxin (detected by Western blotting), only one 'isoform' that focused at pI 5.1 was detected in C. helleri. Both acidic and basic subunits of MT sequences were obtained from C. helleri DNA with primers specific for MT, but only from snakes that had MT in their venom. The sequence identity of the C. helleri acidic subunit to the C. s. scutulatus subunit was 84.9%, whereas the sequence identity of the C. helleri basic subunit was 97% to the C. s. scutulatus basic subunit. Using casein, fibrin, and hide powder azure as substrates, assays for proteolytic activity suggested that C. helleri possesses several different types of metalloproteinases in their venom. However, proteolytic activity was not detected, or present in reduced amounts, in specimens having MT. Clinical neurotoxicity following envenomation by certain populations of C. helleri may be due to MT.

Molecular evolution and structure-function relationships of crotoxin-like and asparagine-6-containing phospholipases A2 in pit viper venoms

Some myotoxic or neurotoxic PLA2s (phospholipases A2) from pit viper venoms contain characteristic N6 substitutions. Our survey of the venoms of more than ten pit viper genera revealed that N6-PLA2s exist only in limited Asian pit vipers of two genera, Protobothrops and Gloydius, and exist as either monomers or the basic subunits of heterodimers in some New World pit vipers. For the newly identified N6-PLA2s, the neuromuscular blocking activities were assayed with the chick biventer cervicis neuromuscular tissue, whereas the increased serum creatine kinase level assessed their myotoxicities. The purified N6-PLA2s from Protobothrops mangshanensis and Gloydius intermedius saxatilis were found to be presynaptic neurotoxins. In contrast, all N6-PLA2s from the venoms of Sistrurus miliarius strackeri, S. m. barbouri, Crotalus viridis viridis, C. lepidus lepidus, Cerrophidion godmani and Bothreichis schlegelii were myotoxins without neurotoxicity even in the presence of crotoxin A. Crotoxin-like complexes were for the first time purified from the venoms of Sitrurus catenatus tergeminus, C. mitchelli mitchelli, C. horridus atricaudatus, C. basiliscus and C. durissus cumanensis. The cDNAs encoding six novel N6-PLA2s and subunits of the crotoxin-like complex from S. c. tergeminus were cloned and fully sequenced. Phylogeny analysis showed that two structural subtypes of N6-PLA2s with either F24 or S24 substitution have been evolved in parallel, possibly descended respectively from species related to present-day Protobothrops and Gloydius. Calmodulin binds all the N6-PLA2s but crotoxin A may inhibit its binding to crotoxin B and to other neurotoxic N6-PLA2s. Structure-activity relationships at various regions of the PLA2 molecules were extensively discussed.