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

Functional variability of snake venom metalloproteinases: adaptive advantages in targeting different prey and implications for human envenomation

Snake venom metalloproteinases (SVMPs) are major components in most viperid venoms that induce disturbances in the hemostatic system and tissues of animals envenomated by snakes. These disturbances are involved in human pathology of snake bites and appear to be essential for the capture and digestion of snake's prey and avoidance of predators. SVMPs are a versatile family of venom toxins acting on different hemostatic targets which are present in venoms in distinct structural forms. However, the reason why a large number of different SVMPs are expressed in some venoms is still unclear. In this study, we evaluated the interference of five isolated SVMPs in blood coagulation of humans, birds and small rodents. P-III class SVMPs (fractions Ic, IIb and IIc) possess gelatinolytic and hemorrhagic activities, and, of these, two also show fibrinolytic activity. P-I class SVMPs (fractions IVa and IVb) are only fibrinolytic. P-III class SVMPs reduced clotting time of human plasma. Fraction IIc was characterized as prothrombin activator and fraction Ic as factor X activator. In the absence of Ca²⁺, a firm clot was observed in chicken blood samples with fractions Ic, IIb and partially with fraction IIc. In contrast, without Ca²⁺, only fraction IIc was able to induce a firm clot in rat blood. In conclusion, functionally distinct forms of SVMPs were found in B. neuwiedi venom that affect distinct mechanisms in the coagulation system of humans, birds and small rodents. Distinct SVMPs appear to be more specialized to rat or chicken blood, strengthening the current hypothesis that toxin diversity enhances the possibilities of the snakes for hunting different prey or evading different predators. This functional diversity also impacts the complexity of human envenoming since different hemostatic mechanisms will be targeted by SVMPs accounting for the complexity of the response of humans to venoms.

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.

Influence of post-starvation extraction time and prey-specific diet in Tityus serrulatus scorpion venom composition and hyaluronidase activity

The role of diet in venom composition has been a topic of intense research interest. This work presents evidence that the variation in the venom composition from the scorpion Tityus serrulatus (Ts) is closely associated with post-starvation extraction time and prey-specific diet. The scorpions were fed with cockroach, cricket, peanut beetle or giant Tenebrio. The venoms demonstrated a pronounced difference in the total protein and toxins composition, which was evaluated by electrophoresis, reversed-phase chromatography, densitometry, hyaluronidase activity and N-terminal sequencing. Indeed, many toxins and peptides, such as Ts1, Ts2, Ts4, Ts5, Ts6, Ts15, Ts19 frag. II, hypotensins 1 and 3, PAPE peptide and peptide 9797 (first described in Ts venom), were all identified in different proportions in the analyzed Ts venoms. This study is pioneer on assessing the influence of the starvation time and the prey diet on hyaluronidase activity as well as to describe a modification of Tricine-gel-electrophoresis to evaluate this enzyme activity. Altogether, this study reveal a large contribution of the extraction time and diet on Ts venom variability as well as present a background to recommend the cockroach diet to obtain higher protein content and the cricket diet to obtain higher hyaluronidase specific activity.

Immunological profile of antivenoms: preclinical analysis of the efficacy of a polyspecific antivenom through antivenomics and neutralization assays

Parenteral administration of animal-derived antivenoms constitutes the mainstay in the treatment of snakebite envenomings. Despite the fact that this therapy has been available for over a century, the detailed understanding of the neutralizing and immunoreactivity profiles of the majority of antivenoms is pending. Currently, a combination of preclinical neutralization tests and 'antivenomics', i.e. a proteomic-based assessment of antivenom immunoreactivity, provides a powerful analytical platform to investigate the preclinical efficacy of antivenoms. In this review, the studies performed on the polyvalent antivenom manufactured by Instituto Clodomiro Picado, Costa Rica, are summarized. This antivenom is prepared by immunizing horses with a mixture of the venoms of Bothrops asper, Crotalus simus and Lachesis stenophrys, and is used in Central America for the treatment of envenomings by viperid species. Overall, the antivenom shows a widespread pattern of immunological reactivity against homologous and heterologous venoms, which correlates with its ability to neutralize lethal, hemorrhagic, myotoxic, coagulant, defibrinogenating, phospholipase A2 and proteinase activities of viperid venoms. At the same time, antivenomics detected several venom components against which the antivenom shows only partial or negligible immunorecognition, such as low molecular mass vasoactive peptides, disintegrins, and some phospholipases A2, P-I metalloproteinases and serine proteinases. Such information can be used to design strategies for enhancing the antibody response of horses against poorly immunogenic, toxicologically-relevant venom components in order to further improve the efficacy of this antivenom.

BIOLOGICAL SIGNIFICANCE

The timely parenteral administration of an appropriate antivenom remains, more than a century after the development of the first serum antivenimeux by Calmette and Phisalix and Bertrand, the only currently effective treatment for snakebite envenomings. A key technical issue in the generation of novel antivenoms is the design of optimized immunization venom mixtures that ensure that the resulting antidotes will be effective against the highest number of venoms from snakes of medical concern across the geographical range where they will be used. Antivenomics is a proteomics-based protocol developed to complement in vitro and in vivo standard preclinical tests in the qualitative and quantitative characterization of the immunological profile and the extent of cross-reactivity of antivenoms against homologous and heterologous venoms. Antivenomics is translational venomics. The combination of antivenomics and neutralization assays represents a powerful analytical platform to investigate the efficacy of antivenoms at the molecular and preclinical levels. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

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.

Intraspecific venom variation in the medically significant Southern Pacific Rattlesnake (Crotalus oreganus helleri): biodiscovery, clinical and evolutionary implications

Due to the extreme variation of venom, which consequently results in drastically variable degrees of neutralization by CroFab antivenom, the management and treatment of envenoming by Crotalus oreganus helleri (the Southern Pacific Rattlesnake), one of the most medically significant snake species in all of North America, has been a clinician's nightmare. This snake has also been the subject of sensational news stories regarding supposed rapid (within the last few decades) evolution of its venom. This research demonstrates for the first time that variable evolutionary selection pressures sculpt the intraspecific molecular diversity of venom components in C. o. helleri. We show that myotoxic β-defensin peptides (aka: crotamines/small basic myotoxic peptides) are secreted in large amounts by all populations. However, the mature toxin-encoding nucleotide regions evolve under the constraints of negative selection, likely as a result of their non-specific mode of action which doesn't enforce them to follow the regime of the classic predator-prey chemical arms race. The hemorrhagic and tissue destroying snake venom metalloproteinases (SVMPs) were secreted in larger amounts by the Catalina Island and Phelan rattlesnake populations, in moderate amounts in the Loma Linda population and in only trace levels by the Idyllwild population. Only the Idyllwild population in the San Jacinto Mountains contained potent presynaptic neurotoxic phospholipase A2 complex characteristic of Mohave Rattlesnake (Crotalus scutulatus) and Neotropical Rattlesnake (Crotalus durissus terrificus). The derived heterodimeric lectin toxins characteristic of viper venoms, which exhibit a diversity of biological activities, including anticoagulation, agonism/antagonism of platelet activation, or procoagulation, appear to have evolved under extremely variable selection pressures. While most lectin α- and β-chains evolved rapidly under the influence of positive Darwinian selection, the β-chain lectin of the Catalina Island population appears to have evolved under the constraint of negative selection. Both lectin chains were conspicuously absent in both the proteomics and transcriptomics of the Idyllwild population. Thus, we not only highlight the tremendous biochemical diversity in C. o. helleri's venom-arsenal, but we also show that they experience remarkably variable strengths of evolutionary selection pressures, within each toxin class among populations and among toxin classes within each population. The mapping of geographical venom variation not only provides additional information regarding venom evolution, but also has direct medical implications by allowing prediction of the clinical effects of rattlesnake bites from different regions. Such information, however, also points to these highly variable venoms as being a rich source of novel toxins which may ultimately prove to be useful in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results have direct implications for the treatment of envenomed patients. The variable venom profile of Crotalus oreganus helleri underscores the biodiscovery potential of novel snake venoms.

Molecular evolution of vertebrate neurotrophins: co-option of the highly conserved nerve growth factor gene into the advanced snake venom arsenalf

Neurotrophins are a diverse class of structurally related proteins, essential for neuronal development, survival, plasticity and regeneration. They are characterized by major family members, such as the nerve growth factors (NGF), brain-derived neurotrophic factors (BDNF) and neurotrophin-3 (NT-3), which have been demonstrated here to lack coding sequence variations and follow the regime of negative selection, highlighting their extremely important conserved role in vertebrate homeostasis. However, in stark contrast, venom NGF secreted as part of the chemical arsenal of the venomous advanced snake family Elapidae (and to a lesser extent Viperidae) have characteristics consistent with the typical accelerated molecular evolution of venom components. This includes a rapid rate of diversification under the significant influence of positive-selection, with the majority of positively-selected sites found in the secreted β-polypeptide chain (74%) and on the molecular surface of the protein (92%), while the core structural and functional residues remain highly constrained. Such focal mutagenesis generates active residues on the toxin molecular surface, which are capable of interacting with novel biological targets in prey to induce a myriad of pharmacological effects. We propose that caenophidian NGFs could participate in prey-envenoming by causing a massive release of chemical mediators from mast cells to mount inflammatory reactions and increase vascular permeability, thereby aiding the spread of other toxins and/or by acting as proapoptotic factors. Despite their presence in reptilian venom having been known for over 60 years, this is the first evidence that venom-secreted NGF follows the molecular evolutionary pattern of other venom components, and thus likely participates in prey-envenomation.

Venomics of New World pit vipers: genus-wide comparisons of venom proteomes across Agkistrodon

We report a genus-wide comparison of venom proteome variation across New World pit vipers in the genus Agkistrodon. Despite the wide variety of habitats occupied by this genus and that all its taxa feed on diverse species of vertebrates and invertebrate prey, the venom proteomes of copperheads, cottonmouths, and cantils are remarkably similar, both in the type and relative abundance of their different toxin families. The venoms from all the eleven species and subspecies sampled showed relatively similar proteolytic and PLA2 activities. In contrast, quantitative differences were observed in hemorrhagic and myotoxic activities in mice. The highest myotoxic activity was observed with the venoms of A. b. bilineatus, followed by A. p. piscivorus, whereas the venoms of A. c. contortrix and A. p. leucostoma induced the lowest myotoxic activity. The venoms of Agkistrodon bilineatus subspecies showed the highest hemorrhagic activity and A. c. contortrix the lowest. Compositional and toxicological analyses agree with clinical observations of envenomations by Agkistrodon in the USA and Central America. A comparative analysis of Agkistrodon shows that venom divergence tracks phylogeny of this genus to a greater extent than in Sistrurus rattlesnakes, suggesting that the distinct natural histories of Agkistrodon and Sistrurus clades may have played a key role in molding the patterns of evolution of their venom protein genes.

BIOLOGICAL SIGNIFICANCE

A deep understanding of the structural and functional profiles of venoms and of the principles governing the evolution of venomous systems is a goal of venomics. Isolated proteomics analyses have been conducted on venoms from many species of vipers and pit vipers. However, making sense of these large inventories of data requires the integration of this information across multiple species to identify evolutionary and ecological trends. Our genus-wide venomics study provides a comprehensive overview of the toxic arsenal across Agkistrodon and a ground for understanding the natural histories of, and clinical observations of envenomations by, species of this genus.

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.

Integrated "omics" profiling indicates that miRNAs are modulators of the ontogenetic venom composition shift in the Central American rattlesnake, Crotalus simus simus

BACKGROUND

Understanding the processes that drive the evolution of snake venom is a topic of great research interest in molecular and evolutionary toxinology. Recent studies suggest that ontogenetic changes in venom composition are genetically controlled rather than environmentally induced. However, the molecular mechanisms underlying these changes remain elusive. Here we have explored the basis and level of regulation of the ontogenetic shift in the venom composition of the Central American rattlesnake, Crotalus s. simus using a combined proteomics and transcriptomics approach.

RESULTS

Proteomic analysis showed that the ontogenetic shift in the venom composition of C. s. simus is essentially characterized by a gradual reduction in the expression of serine proteinases and PLA2 molecules, particularly crotoxin, a β-neurotoxic heterodimeric PLA2, concominantly with an increment of PI and PIII metalloproteinases at age 9-18 months. Comparison of the transcriptional activity of the venom glands of neonate and adult C. s. simus specimens indicated that their transcriptomes exhibit indistinguisable toxin family profiles, suggesting that the elusive mechanism by which shared transcriptomes generate divergent venom phenotypes may operate post-transcriptionally. Specifically, miRNAs with frequency count of 1000 or greater exhibited an uneven distribution between the newborn and adult datasets. Of note, 590 copies of a miRNA targeting crotoxin B-subunit was exclusively found in the transcriptome of the adult snake, whereas 1185 copies of a miRNA complementary to a PIII-SVMP mRNA was uniquely present in the newborn dataset. These results support the view that age-dependent changes in the concentration of miRNA modulating the transition from a crotoxin-rich to a SVMP-rich venom from birth through adulthood can potentially explain what is observed in the proteomic analysis of the ontogenetic changes in the venom composition of C. s. simus.

CONCLUSIONS

Existing snake venom toxins are the result of early recruitment events in the Toxicofera clade of reptiles by which ordinary genes were duplicated, and the new genes selectively expressed in the venom gland and amplified to multigene families with extensive neofunctionalization throughout the approximately 112-125 million years of ophidian evolution. Our findings support the view that understanding the phenotypic diversity of snake venoms requires a deep knowledge of the mechanisms regulating the transcriptional and translational activity of the venom gland. Our results suggest a functional role for miRNAs. The impact of specific miRNAs in the modulation of venom composition, and the integration of the mechanisms responsible for the generation of these miRNAs in the evolutionary landscape of the snake's venom gland, are further challenges for future research.