Electrophoretic profiles and biological activities: intraspecific variation in the venom of the Malayan pit viper (Calloselasma rhodostoma)

Indian green pit vipers: A lesser-known snake group of north-east India

Green pit vipers are one of the most widely distributed group of venomous snakes in south-east Asia. In Indian, green pit vipers are found in the Northern and North-eastern states spreading across eastern and central India and one of the lesser studied venoms. High morphological similarity among them has been a long-established challenge for species identification, however, a total of six species of Indian green pit viper belonging to genus Trimeresurus, Popeia and Viridovipera has been reported from North-east India. Biochemical and biological studies have revealed that venom exhibits substantial variation in protein expression level along with functional variability. The symptoms of envenomation are painful swelling at bite site, bleeding, necrosis along with systemic toxicity such as prolonged coagulopathy. Clinical data of green pit viper envenomated patients from Demow community health centre, Assam advocated against the use of Indian polyvalent antivenom pressing the need for a suitable antivenom for the treatment of green pit viper envenomation. To design effective and specific antivenom for green pit vipers, unveiling the proteome profile of these snakes is needed. In this study, a comparative venomic of green pit vipers of Northern and North-eastern India, their clinical manifestation as well as treatment protocol has been reviewed.

Snakes on a plain: biotic and abiotic factors determine venom compositional variation in a wide-ranging generalist rattlesnake

BACKGROUND

Snake venoms are trophic adaptations that represent an ideal model to examine the evolutionary factors that shape polymorphic traits under strong natural selection. Venom compositional variation is substantial within and among venomous snake species. However, the forces shaping this phenotypic complexity, as well as the potential integrated roles of biotic and abiotic factors, have received little attention. Here, we investigate geographic variation in venom composition in a wide-ranging rattlesnake (Crotalus viridis viridis) and contextualize this variation by investigating dietary, phylogenetic, and environmental variables that covary with venom.

RESULTS

Using shotgun proteomics, venom biochemical profiling, and lethality assays, we identify 2 distinct divergent phenotypes that characterize major axes of venom variation in this species: a myotoxin-rich phenotype and a snake venom metalloprotease (SVMP)-rich phenotype. We find that dietary availability and temperature-related abiotic factors are correlated with geographic trends in venom composition.

CONCLUSIONS

Our findings highlight the potential for snake venoms to vary extensively within species, for this variation to be driven by biotic and abiotic factors, and for the importance of integrating biotic and abiotic variation for understanding complex trait evolution. Links between venom variation and variation in biotic and abiotic factors indicate that venom variation likely results from substantial geographic variation in selection regimes that determine the efficacy of venom phenotypes across populations and snake species. Our results highlight the cascading influence of abiotic factors on biotic factors that ultimately shape venom phenotype, providing evidence for a central role of local selection as a key driver of venom variation.

Multilevel Comparison of Indian Naja Venoms and Their Cross-Reactivity with Indian Polyvalent Antivenoms

Snake envenoming is caused by many biological species, rather than a single infectious agent, each with a multiplicity of toxins in their venom. Hence, developing effective treatments is challenging, especially in biodiverse and biogeographically complex countries such as India. The present study represents the first genus-wide proteomics analysis of venom composition across Naja species (N. naja, N. oxiana, and N. kaouthia) found in mainland India. Venom proteomes were consistent between individuals from the same localities in terms of the toxin families present, but not in the relative abundance of those in the venom. There appears to be more compositional variation among N. naja from different locations than among N. kaouthia. Immunoblotting and in vitro neutralization assays indicated cross-reactivity with Indian polyvalent antivenom, in which antibodies raised against N. naja are present. However, we observed ineffective neutralization of PLA2 activities of N. naja venoms from locations distant from the source of immunizing venoms. Antivenom immunoprofiling by antivenomics revealed differential antigenicity of venoms from N. kaouthia and N. oxiana, and poor reactivity towards 3FTxs and PLA2s. Moreover, there was considerable variation between antivenoms from different manufacturers. These data indicate that improvements to antivenom manufacturing in India are highly desirable.

A comparison of two different analytical workflows to determine the venom proteome composition of Naja kaouthia from North-East India and immunological profiling of venom against commercial antivenoms

The Indian monocled cobra (Naja kaouthia) is one of the most prevalent venomous snakes in northeast India (NEI) and is the cause of many fatalities. The composition of NEI N. kaouthia venom (NkV) was deciphered using two different proteomic approaches: (i) 1D SDS-PAGE coupled to label-free quantification of protein bands using stringent identification criteria and (ii) reversed-phase high-performance liquid chromatography (RP-HPLC) followed by quantification based on area under the RP-HPLC peaks. The proteomic data from both strategies were compared. Proteomic analyses from both workflows identified 32 proteins (toxins) distributed over 10-14 snake venom protein families in NEI NkV. The relative abundances of the venom proteins determined from the analytical workflows coincided with the densitometry band intensities of the NEI NkV. Phospholipase A₂ (13.1-16.0%) and three-finger toxins (58.5-64.2%) represented the most abundant enzymatic and non-enzymatic proteins in NEI NkV, respectively. Immuno-cross-reactivity studies by enzyme-linked immunoassay and immunoblot analyses pointed to the poor efficacy of commercial PAVs in recognizing the low molecular mass (<15 kDa) toxins of NEI NkV. Spectrofluorometric titration determined the presence of NEI NkV-specific antibodies in commercial PAV, at a level that was higher than that previously reported for eastern India NkV-specific antibodies in commercial antivenom.

Venom of several Indian green pit vipers: Comparison of biochemical activities and cross-reactivity with antivenoms

Green pit vipers, a name that can refer to several unrelated species, comprise a large group of venomous snakes found across the humid areas of tropical and sub-tropical Asia, and are responsible for most of the bite cases across this region. In India, green pit vipers belonging to several genera are prevalent in the northern and north-eastern hilly region, unrelated to species present in the peninsular region. In the present study, crude venom of representative species of green pit vipers present in the north and north-eastern hilly region of India (Trimeresurus erythrurus, T. septentrionalis, Viridovipera medoensis, and Popiea popieorum) were characterized to elucidate venom composition and venom variation. Profiling of crude venoms using SDS-PAGE and RP-HPLC methods revealed quantitative differences among the species. Further, in vitro biochemical assays reveal variable levels of phospholipase activity, coagulation activity, thrombin-like activity, fibrinogenolytic and haemolytic activity. This correlates with the pseudo-procoagulant effects on the haemostatic system of victims, which causes consumptive coagulopathy, frequently observed in patients bitten by green pit vipers. The immunoreactivity of Indian polyvalent antivenom and Thai green pit viper antivenom towards crude venoms were also evaluated by western blotting and inhibition of biochemical activities. The results exhibited poor efficacy of Indian polyvalent antivenom in neutralizing the venom toxins of crude venoms; however, Thai green pit viper antivenin (raised against the venom of Trimeresurus allbolabris, not present in India) showed higher immunoreactivity towards congeneric venoms tested. Analysis of green pit viper bite patients records from a community health centre in Assam, India, further revealed the inability of Indian polyvalent antivenom to reverse the extended coagulopathy featured.

Proteomic comparison of adult and juvenile Santa Catalina rattlesnake (Crotalus catalinensis) venom

Rattlesnake's venom constitutes an important ecological trait that dynamically changes over time. Venoms of adult and juvenile rattleless rattlesnakes, Crotalus catalinensis, an endemic insular species from the Gulf of California, were compared by electrophoretic profile, fibrinogenolytic activity, and proteomic composition to assess ontogenetic variability. The SDS-PAGE profiles show important differences at 12, 22, and 45 kDa, which were prominent in adult samples and absent in juvenile samples, while bands around 20, 25, and 70 kDa are almost absent in adults. Both venoms hydrolyze Aa and Bb chains of fibrinogen generating different patterns of degradation products. This activity was partially inhibited by EDTA and PMSF and completely abolished only in the presence of both inhibitors. More than 260 proteins were identified and quantified in both venoms by proteomic analysis. Metalloproteinases (more than 60%), serine proteinases (14.5% in adult venom and 17.7% in juvenile venom), and C-type lectins (7.1 and 5.9%) represent the three most abundant toxin-related protein families. Bradykinin inhibitor peptides and L-amino acid oxidases were not detected in juvenile venom. A protein-specific comparison shows that adult and juvenile venom share about 30.5% of total toxin-related proteins, while 32% and 35% are exclusively present in adult and juvenile venoms, respectively. This work represents one of the first efforts to understand phenotypic diversity in the venom composition of insular rattlesnake species from Mexico.

Proteomic Investigations of Two Pakistani Naja Snake Venoms Species Unravel the Venom Complexity, Posttranslational Modifications, and Presence of Extracellular Vesicles

Latest advancement of omics technologies allows in-depth characterization of venom compositions. In the present work we present a proteomic study of two snake venoms of the genus Naja i.e., Naja naja (black cobra) and Naja oxiana (brown cobra) of Pakistani origin. The present study has shown that these snake venoms consist of a highly diversified proteome. Furthermore, the data also revealed variation among closely related species. High throughput mass spectrometric analysis of the venom proteome allowed to identify for the N. naja venom 34 protein families and for the N. oxiana 24 protein families. The comparative evaluation of the two venoms showed that N. naja consists of a more complex venom proteome than N. oxiana venom. Analysis also showed N-terminal acetylation (N-ace) of a few proteins in both venoms. To the best of our knowledge, this is the first study revealing this posttranslational modification in snake venom. N-ace can shed light on the mechanism of regulation of venom proteins inside the venom gland. Furthermore, our data showed the presence of other body proteins, e.g., ankyrin repeats, leucine repeats, zinc finger, cobra serum albumin, transferrin, insulin, deoxyribonuclease-2-alpha, and other regulatory proteins in these venoms. Interestingly, our data identified Ras-GTpase type of proteins, which indicate the presence of extracellular vesicles in the venom. The data can support the production of distinct and specific anti-venoms and also allow a better understanding of the envenomation and mechanism of distribution of toxins. Data are available via ProteomeXchange with identifier PXD018726.

Toxicological comparison of Crotalus ruber lucasensis venom from different ecoregions of the Baja California Peninsula

The Baja California Peninsula possesses a mosaic of ecoregions that offers a wide variety of environments for the species that here inhabit. Here we report biological variations in. Crotalus ruber lucasensis venom from arid, semiarid and tropical eco-regions. Lethal (1.4-6.8 mg/kg), edematogenic (0.3-0.5 μg) and defibrinogenating (from non-detectable to 20 μg) activities were found to have significant differences among eco-regions.

Intra-Specific Venom Variation in the Australian Coastal Taipan Oxyuranus scutellatus

Intra-specific venom variation has the potential to provide important insights into the evolution of snake venom, but remains a relatively neglected aspect of snake venom studies. We investigated the venom from 13 individual coastal taipans Oxyuranus scutellatus from four localities on the north-east coast of Australia, spanning a distance of 2000 km. The intra-specific variation in taipan venom was considerably less than the inter-specific variation between it and the other Australian elapids to which it was compared. The electrophoretic venom profile of O. scutellatus was visually different to six other genera of Australian elapids, but not to its congener inland taipan O. microlepidotus. There was minimal geographical variation in taipan venom, as the intra-population variation exceeded the inter-population variation for enzymatic activity, procoagulant activity, and the abundance of neurotoxins. The pre-synaptic neurotoxin (taipoxin) was more abundant than the post-synaptic neurotoxins (3FTx), with a median of 11.0% (interquartile range (IQR): 9.7% to 18.3%; range: 6.7% to 23.6%) vs. a median of 3.4% (IQR: 0.4% to 6.7%; range: 0% to 8.1%). Three taipan individuals almost completely lacked post-synaptic neurotoxins, which was not associated with geography and occurred within two populations. We found no evidence of sexual dimorphism in taipan venom. Our study provides a basis for evaluating the significance of intra-specific venom variation within a phylogenetic context by comparing it to the inter-specific and inter-generic variation. The considerable intra-population variation we observed supports the use of several unpooled individuals from each population when making inter-specific comparisons.

Research into the Causes of Venom-Induced Mortality and Morbidity Identifies New Therapeutic Opportunities

Snakebite primarily affects rural subsistent farming populations in underdeveloped and developing nations. The annual number of deaths (100,000) and physical disabilities (400,000) of snakebite victims is a societal tragedy that poses a significant added socioeconomic burden to the society. Antivenom therapy is the treatment of choice for snakebite but, as testified by the continuing high rates of mortality and morbidity, too many rural tropical snakebite victims fail to access effective treatment. Here, we advocate for more basic research to better understand the pathogenesis of systemic and local envenoming and describe how research outcomes can identify novel snakebite therapeutic strategies with the potential to be more accessible and affordable to victims than current treatment.

Activity of two key toxin groups in Australian elapid venoms show a strong correlation to phylogeny but not to diet

Background

The relative influence of diet and phylogeny on snake venom activity is a poorly understood aspect of snake venom evolution. We measured the activity of two enzyme toxin groups – phospholipase A₂ (PLA₂), and L-amino acid oxidase (LAAO) – in the venom of 39 species of Australian elapids (40% of terrestrial species diversity) and used linear parsimony and BayesTraits to investigate any correlation between enzyme activity and phylogeny or diet.

Results

PLA₂ activity ranged from 0 to 481 nmol/min/mg of venom, and LAAO activity ranged from 0 to 351 nmol/min/mg. Phylogenetic comparative methods, implemented in BayesTraits showed that enzyme activity was strongly correlated with phylogeny, more so for LAAO activity. For example, LAAO activity was absent in both the Vermicella and Pseudonaja/Oxyuranus clade, supporting previously proposed relationships among these disparate taxa. There was no association between broad dietary categories and either enzyme activity. There was strong evidence for faster initial rates of change over evolutionary time for LAAO (delta parameter mean 0.2), but no such pattern in PLA₂ (delta parameter mean 0.64). There were some exceptions to the phylogenetic patterns of enzyme activity: different PLA₂ activity in the ecologically similar sister-species Denisonia devisi and D. maculata; large inter-specific differences in PLA₂ activity in Hoplocephalus and Austrelaps.

Conclusions

We have shown that phylogeny is a stronger influence on venom enzyme activity than diet for two of the four major enzyme families present in snake venoms. PLA₂ and LAAO activities had contrasting evolutionary dynamics with the higher delta value for PLA₂ Some species/individuals lacked activity in one protein family suggesting that the loss of single protein family may not incur a significant fitness cost.

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

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

Venom composition of adult Western Diamondback Rattlesnakes (Crotalus atrox) maintained under controlled diet and environmental conditions shows only minor changes

Many species of snakes produce venom as a chemical means of procuring potentially fractious prey. Studies have increasingly focused on venom compositional variation between and within individual snakes of the same species/subspecies, with significant differences often being observed. This variation in composition has been attributed to differences in age, season, diet, and environment, suggesting that these factors could help explain the inter- and intra-specific variation found in some snake venoms, perhaps via some type of feedback mechanism(s). To address several of these possible sources of variation, this study utilized wild-caught Western Diamondback Rattlesnakes (Crotalus atrox) from Cochise Co., AZ. Sixteen adult C. atrox were maintained in the lab on a diet of NSA mice for eight months to determine whether venom composition changed in captivity under a static diet in a stable environment. Reducing 1-D SDS-PAGE, fibrinogen degradation assays, reversed-phase HPLC, and MALDI-TOF mass spectrometry revealed only minor differences over time within individuals. Venom L-amino acid oxidase (LAAO) and phosphodiesterase activities significantly increased over the course of captivity, with no changes occurring in azocasein metalloproteinase, kallikrein-like serine proteinase (KLSP), or thrombin-like serine proteinase (TLSP) activities. Snake total length was positively correlated with TLSP activity and negatively correlated with LAAO and KLSP activity. There was typically a much higher degree of variation between individuals than within individuals for all analyses performed and measurements collected. Because the overall "fingerprint" of each snake's venom remained more/less constant, it is concluded that biologically significant changes in venom composition did not occur within individual C. atrox as a function of captivity/diet. However, this study does indicate that differences in activity levels do occur in minor venom enzyme components, but the differences observed are likely to be of minimal significance to the production of antivenom or to subsequent treatment of human envenomations.

Basal but divergent: Clinical implications of differential coagulotoxicity in a clade of Asian vipers

Envenomations by Asian pitvipers can induce multiple clinical complications resulting from coagulopathic and neuropathic effects. While intense research has been undertaken for some species, functional coagulopathic effects have been neglected. As these species' venoms affect the blood coagulation cascade we investigated their effects upon the human clotting cascade using venoms of species from the Azemiops, Calloselasma, Deinagkistrodon and Hypnale genera. Calloselasma rhodostoma, Deinagkistrodon acutus, and Hypnale hypnale produced net anticoagulant effects through pseudo-procoagulant clotting of fibrinogen, resulting in weak, unstable, transient fibrin clots. Tropidolaemus wagleri was only weakly pseudo-procoagulant, clotting fibrinogen with only a negligible net anticoagulant effect. Azemiops feae and Tropidolaemus subannulatus did not affect clotting. This is the first study to examine in a phylogenetic context the coagulotoxic effects of related genera of basal Asiatic pit-vipers. The results reveal substantial variation between sister genera, providing crucial information about clinical effects and implications for antivenom cross-reactivity.

The habu genome reveals accelerated evolution of venom protein genes

Evolution of novel traits is a challenging subject in biological research. Several snake lineages developed elaborate venom systems to deliver complex protein mixtures for prey capture. To understand mechanisms involved in snake venom evolution, we decoded here the ~1.4-Gb genome of a habu, Protobothrops flavoviridis. We identified 60 snake venom protein genes (SV) and 224 non-venom paralogs (NV), belonging to 18 gene families. Molecular phylogeny reveals early divergence of SV and NV genes, suggesting that one of the four copies generated through two rounds of whole-genome duplication was modified for use as a toxin. Among them, both SV and NV genes in four major components were extensively duplicated after their diversification, but accelerated evolution is evident exclusively in the SV genes. Both venom-related SV and NV genes are significantly enriched in microchromosomes. The present study thus provides a genetic background for evolution of snake venom composition.

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.

Correlation between ontogenetic dietary shifts and venom variation in Australian brown snakes (Pseudonaja)

Venom is a key evolutionary trait, as evidenced by its widespread convergent evolution across the animal kingdom. In an escalating prey-predator arms race, venoms evolve rapidly to guarantee predatory or defensive success. Variation in venom composition is ubiquitous among snakes. Here, we tested variation in venom activity on substrates relevant to blood coagulation among Pseudonaja (brown snake) species, Australian elapids responsible for the majority of medically important human envenomations in Australia. A functional approach was employed to elucidate interspecific variation in venom activity in all nine currently recognised species of Pseudonaja. Fluorometric enzymatic activity assays were performed to test variation in whole venom procoagulant activity among species. Analyses confirmed the previously documented ontogenetic shift from non-coagulopathic venom in juveniles to coagulopathic venom as adults, except for the case of P. modesta, which retains non-coagulopathic venom as an adult. These shifts in venom activity correlate with documented ontogenetic shifts in diet among brown snakes from specialisation on reptilian prey as juveniles (and throughout the life cycle of P. modesta), to a more generalised diet in adults that includes mammals. The results of this study bring to light findings relevant to both clinical and evolutionary toxinology.

Comparison of proteomic profiles of the venoms of two of the 'Big Four' snakes of India, the Indian cobra (Naja naja) and the common krait (Bungarus caeruleus), and analyses of their toxins

Snake venoms are mixtures of biologically-active proteins and peptides, and several studies have described the characteristics of some of these toxins. However, complete proteomic profiling of the venoms of many snake species has not yet been done. The Indian cobra (Naja naja) and common krait (Bungarus caeruleus) are elapid snake species that are among the 'Big Four' responsible for the majority of human snake envenomation cases in India. As understanding the composition and complexity of venoms is necessary for successful treatment of envenomation in humans, we utilized three different proteomic profiling approaches to characterize these venoms: i) one-dimensional SDS-PAGE coupled with in-gel tryptic digestion and electrospray tandem mass spectrometry (ESI-LC-MS/MS) of individual protein bands; ii) in-solution tryptic digestion of crude venoms coupled with ESI-LC-MS/MS; and iii) separation by gel-filtration chromatography coupled with tryptic digestion and ESI-LC-MS/MS of separated fractions. From the generated data, 81 and 46 different proteins were identified from N. naja and B. caeruleus venoms, respectively, belonging to fifteen different protein families. Venoms from both species were found to contain a variety of phospholipases A₂ and three-finger toxins, whereas relatively higher numbers of snake venom metalloproteinases were found in N. naja compared to B. caeruleus venom. The analyses also identified less represented venom proteins including L-amino acid oxidases, cysteine-rich secretory proteins, 5'-nucleotidases and venom nerve growth factors. Further, Kunitz-type serine protease inhibitors, cobra venom factors, phosphodiesterases, vespryns and aminopeptidases were identified in the N. naja venom, while acetylcholinesterases and hyaluronidases were found in the B. caeruleus venom. We further analyzed protein coverage (Lys/Arg rich and poor regions as well as potential glycosylation sites) using in-house software. These studies expand our understanding of the proteomes of the venoms of these two medically-important species.

Venomics of Tropidolaemus wagleri, the sexually dimorphic temple pit viper: Unveiling a deeply conserved atypical toxin arsenal

Tropidolaemus wagleri (temple pit viper) is a medically important snake in Southeast Asia. It displays distinct sexual dimorphism and prey specificity, however its venomics and inter-sex venom variation have not been thoroughly investigated. Applying reverse-phase HPLC, we demonstrated that the venom profiles were not significantly affected by sex and geographical locality (Peninsular Malaya, insular Penang, insular Sumatra) of the snakes. Essentially, venoms of both sexes share comparable intravenous median lethal dose (LD₅₀) (0.56-0.63 μg/g) and cause neurotoxic envenomation in mice. LCMS/MS identified six waglerin forms as the predominant lethal principles, comprising 38.2% of total venom proteins. Fourteen other toxin-protein families identified include phospholipase A₂, serine proteinase, snaclec and metalloproteinase. In mice, HPLC fractions containing these proteins showed insignificant contribution to the overall venom lethality. Besides, the unique elution pattern of approximately 34.5% of non-lethal, low molecular mass proteins (3-5 kDa) on HPLC could be potential biomarker for this primitive crotalid species. Together, the study unveiled the venom proteome of T. wagleri that is atypical among many pit vipers as it comprises abundant neurotoxic peptides (waglerins) but little hemotoxic proteinases. The findings also revealed that the venom is relatively well conserved intraspecifically despite the drastic morphological differences between sexes.

Sex Differences in Defensive Behavior and Venom of The Striped Bark Scorpion Centruroides vittatus (Scorpiones: Buthidae)

Studies of venom variability have advanced from describing the mechanisms of action and relative potency of medically important toxins to understanding the ecological and evolutionary causes of the variability itself. While most studies have focused on differences in venoms among taxa, populations, or age-classes, there may be intersexual effects as well. Striped bark scorpions (Centruroides vittatus) provide a good model for examining sex differences in venom composition and efficacy, as this species exhibits dramatic sexual dimorphism in both size and defensive behavior; when threatened by an enemy, larger, slower females stand and fight while smaller, fleeter males prefer to run. We here add evidence suggesting that male and female C. vittatus indeed have different defensive propensities; when threatened via an electrical stimulus, females were more likely to sting than were males. We reasoned that intersexual differences in defensive phenotypes would select for venoms with different functions in the two sexes; female venoms should be effective at predator deterrence, whereas male venoms, less utilized defensively, might be better suited to capturing prey or courting females. This rationale led to our predictions that females would inject more venom and/or possess more painful venom than males. We were wrong. While females do inject more venom than males in a defensive sting, females are also larger; when adjusted for body size, male and female C. vittatus commit equal masses of venom in a sting to a potential enemy. Additionally, house mice (Mus musculus) find an injection of male venom more irritating than an equal amount of female venom, likely because male venom contains more of the toxins that induce pain. Taken together, our results suggest that identifying the ultimate causes of venom variability will, as we move beyond adaptive storytelling, be hard-won.

Venom Resistance as a Model for Understanding the Molecular Basis of Complex Coevolutionary Adaptations

SynopsisVenom and venom resistance are molecular phenotypes widely considered to have diversified through coevolution between predators and prey. However, while evolutionary and functional studies on venom have been extensive, little is known about the molecular basis, variation, and complexity of venom resistance. We review known mechanisms of venom resistance and relate these mechanisms to their predicted impact on coevolutionary dynamics with venomous enemies. We then describe two conceptual approaches which can be used to examine venom/resistance systems. At the intraspecific level, tests of local adaptation in venom and resistance phenotypes can identify the functional mechanisms governing the outcomes of coevolution. At deeper evolutionary timescales, the combination of phylogenetically informed analyses of protein evolution coupled with studies of protein function promise to elucidate the mode and tempo of evolutionary change on potentially coevolving genes. We highlight case studies that use each approach to extend our knowledge of these systems as well as address larger questions about coevolutionary dynamics. We argue that resistance and venom are phenotypic traits which hold exceptional promise for investigating the mechanisms, dynamics, and outcomes of coevolution at the molecular level. Furthermore, extending the understanding of single gene-for-gene interactions to the whole resistance and venom phenotypes may provide a model system for examining the molecular and evolutionary dynamics of complex multi-gene interactions.

Canopy Venom: Proteomic Comparison among New World Arboreal Pit-Viper Venoms

Central and South American pitvipers, belonging to the genera Bothrops and Bothriechis, have independently evolved arboreal tendencies. Little is known regarding the composition and activity of their venoms. In order to close this knowledge gap, venom proteomics and toxin activity of species of Bothriechis, and Bothrops (including Bothriopsis) were investigated through established analytical methods. A combination of proteomics and bioactivity techniques was used to demonstrate a similar diversification of venom composition between large and small species within Bothriechis and Bothriopsis. Increasing our understanding of the evolution of complex venom cocktails may facilitate future biodiscoveries.

Stabilising the Integrity of Snake Venom mRNA Stored under Tropical Field Conditions Expands Research Horizons

BACKGROUND

Snake venoms contain many proteinaceous toxins that can cause severe pathology and mortality in snakebite victims. Interestingly, mRNA encoding such toxins can be recovered directly from venom, although yields are low and quality is unknown. It also remains unclear whether such RNA contains information about toxin isoforms and whether it is representative of mRNA recovered from conventional sources, such as the venom gland. Answering these questions will address the feasibility of using venom-derived RNA for future research relevant to biomedical and antivenom applications.

METHODOLOGY/PRINCIPAL FINDINGS

Venom was extracted from several species of snake, including both members of the Viperidae and Elapidae, and either lyophilized or immediately added to TRIzol reagent. TRIzol-treated venom was incubated at a range of temperatures (4-37°C) for a range of durations (0-48 hours), followed by subsequent RNA isolation and assessments of RNA quantity and quality. Subsequently, full-length toxin transcripts were targeted for PCR amplification and Sanger sequencing. TRIzol-treated venom yielded total RNA of greater quantity and quality than lyophilized venom, and with quality comparable to venom gland-derived RNA. Full-length sequences from multiple Viperidae and Elapidae toxin families were successfully PCR amplified from TRIzol-treated venom RNA. We demonstrated that venom can be stored in TRIzol for 48 hours at 4-19°C, and 8 hours at 37°C, at minimal cost to RNA quality, and found that venom RNA encoded multiple toxin isoforms that seemed homologous (98-99% identity) to those found in the venom gland.

CONCLUSIONS/SIGNIFICANCE

The non-invasive experimental modifications we propose will facilitate the future investigation of venom composition by using venom as an alternative source to venom gland tissue for RNA-based studies, thus obviating the undesirable need to sacrifice snakes for such research purposes. In addition, they expand research horizons to rare, endangered or protected snake species and provide more flexibility to performing fieldwork on venomous snakes in tropical conditions.

An analysis of venom ontogeny and prey-specific toxicity in the Monocled Cobra (Naja kaouthia)

Venoms of snakes of the family Elapidae (cobras, kraits, mambas, and relatives) are predominantly composed of numerous phospholipases A2 (PLA2s) and three-finger toxins (3FTxs), some of which are lethal while others are not significantly toxic. Currently, the only identified prey-specific toxins are several nonconventional 3FTxs, and given the large diversity of 3FTxs within Monocled Cobra (Naja kaouthia) venom, it was hypothesized that several 3FTxs, previously found to be non-toxic or weakly toxic 3FTxs in murine models, could potentially be toxic towards non-murine prey. Additionally, it was hypothesized that ontogenetic dietary shifts will be correlated with observable changes in specific 3FTx isoform abundance. Adult and juvenile N. kaouthia venom composition was investigated using ion-exchange FPLC, 1D and 2D SDS-PAGE, mass spectrometry, and various enzymatic and LD50 assays. Alpha-cobratoxin (α-elapitoxin) was the only significantly toxic (LD50 < 1 μg/g) 3FTx found in N. kaouthia venom and was equally toxic toward both lizard and mouse models. The abundance and diversity of 3FTxs and most enzyme activities did not vary between adult and juvenile cobra venoms; however, total venom PLA2 activity and specific PLA2 isoforms did vary, with juveniles lacking several of the least acidic PLA2s, and these differences could have both biological (related to predation) and clinical (antivenom efficacy) implications. Nevertheless, the ubiquitous presence of α-cobratoxin in both adult and juvenile cobra venoms, with high toxicity toward both reptiles and mammals, represents a venom compositional strategy wherein a single potent toxin effectively immobilizes a variety of prey types encountered across life history stages.

Venom and Purified Toxins of the Spectacled Cobra (Naja naja) from Pakistan: Insights into Toxicity and Antivenom Neutralization

Geographical variations of snake venoms can result in suboptimal effectiveness of Indian antivenoms that are currently used in most South Asian countries. This study investigated the toxicity and neutralization profile of the venom and toxins from Pakistani spectacled cobra, Naja naja, using VINS polyvalent antivenom (VPAV, India), Naja kaouthia monovalent antivenom (NKMAV, Thailand), and neuro bivalent antivenom (NBAV, Taiwan). Cation-exchange and reverse-phase high-performance liquid chromatography fractionations followed by toxin identification through liquid chromatography-mass spectrometry (MS)/MS indicated that the venom comprised mainly of postsynaptic neurotoxins (NTXs) (long neurotoxins [LNTXs], 28.3%; short neurotoxins [SNTXs], 8%), cytotoxins (CTXs) (31.2%), and acidic phospholipases A2 (12.3%). NKMAV is the most effective in neutralizing the lethal effect of the venom (potency = 1.1 mg venom/mL) and its LNTX (potency = 0.5 mg toxin/mL), consistent with the high content of LNTX in N. kaouthia venom. VPAV was effective in neutralizing the CTX (potency = 0.4 mg toxin/mL), in agreement with the higher CTX abundance in Indian cobra venom. All the three antivenoms were weak in neutralizing the SNTX (potency = 0.03-0.04 mg toxin/mL), including NBAV that was raised from the SNTX-rich Taiwanese cobra venom. In a challenge-rescue experiment, envenomed mice were prevented from death by a maximal dose of VPAV (intravenous 200 μL) but the recovery from paralysis was slow, indicating the need for higher or repeated doses of VPAV. Our results suggest that optimal neutralization for Pakistani N. naja venom may be achieved by improving the formulation of antivenom production to enhance antivenom immunoreactivity against long and SNTXs.

Ecological venomics: How genomics, transcriptomics and proteomics can shed new light on the ecology and evolution of venom

Animal venom is a complex cocktail of bioactive chemicals that traditionally drew interest mostly from biochemists and pharmacologists. However, in recent years the evolutionary and ecological importance of venom is realized as this trait has direct and strong influence on interactions between species. Moreover, venom content can be modulated by environmental factors. Like many other fields of biology, venom research has been revolutionized in recent years by the introduction of systems biology approaches, i.e., genomics, transcriptomics and proteomics. The employment of these methods in venom research is known as 'venomics'. In this review we describe the history and recent advancements of venomics and discuss how they are employed in studying venom in general and in particular in the context of evolutionary ecology. We also discuss the pitfalls and challenges of venomics and what the future may hold for this emerging scientific field.

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.

Cutaneous loxoscelism caused by Loxosceles similis venom and neutralization capacity of its specific antivenom

Members of the spider genus Loxosceles pose a marked health risk to humans because of the seriousness of the necrotic and systemic effects of their bite, known as loxoscelism. The recent confirmation of Loxosceles similis in residences of Belo Horizonte in Minas Gerais Province, Brazil increases the local potential risk of loxoscelism at higher levels. The first characterization of the venom from this species showed that its main biological effects had a similar intensity as other species (e.g. Loxosceles intermedia, Loxosceles laeta, and Loxosceles gaucho). Therefore, we wished to further analyse the biological activity of the L. similis venom as well as the capacity of anti-L. similis-venom serum to reduce dermonecrotic effects to rabbit skin. Histological analysis of rabbit skin 2, 4 and 8h after intradermal injection of L. similis venom demonstrated a dense inflammatory infiltrate, edema, degeneration and necrosis of the skin muscle, dissociation of collagen fibers, and disruption of reticular fibers. Importantly, pre-incubation of the venom with anti-L. similis-venom serum significantly decreased all of these effects. Anti-L. similis antivenom generated antibodies that were strongly reactive to L. similis venom and capable of neutralizing the dermonecrotic effects in rabbits caused by this venom. Moreover, the antivenom significantly reduced the sphingomyelinase activity of L. similis crude venom. Venoms produced by male and female spiders were equally reactive towards anti-L. similis and anti-L. intermedia antivenoms, but female venom induced larger lesions on rabbits. In contrast, female venom acted as an immunization enhancer and protected animals from L. similis envenomation to a greater degree than male venom. In conclusion, the results shown in this study for L. similis antivenom merits a more in depth study of its properties, which may become a valuable tool against loxoscelism.

Intraspecies variability in Vipera ammodytes ammodytes venom related to its toxicity and immunogenic potential

Vipera ammodytes is the most venomous European snake, whose venom has been used as antigen for immunization of antivenom-producing animals. Same as venom of any other snake, it is a complex mixture of proteins, peptides and other compounds which biochemical and pharmacological variability has been demonstrated at interspecies and intraspecies level. In this work we demonstrated intraspecific variability between 8 venom production batches using both the conventional and the new methodology. Moreover, in contrast to the literature on different venoms' variability, for the first time we were able to select those biochemical differences that are related to and give information on the venom's toxicity and immunogenicity. We have shown that methods quantifying ammodytoxin (the most toxic compound identified so far in the Vipera ammodytes ammodytes venom) content of the venom clearly distinguish between high and low immunogenic venoms.

Intraspecific variation in the Egyptian scorpion Scorpio maurus palmatus venom collected from different biotopes

The present study was conducted to explore the following hypotheses: (i) do scorpions (Scorpio maurus palmatus) from different biotopes exhibit intraspecific diversity in their venom? (ii) if so, is this variation associated with ecological or genetic factors, geographical distance, and/or multiple interrelated parameters? To address these questions, scorpions were collected from four geographically isolated localities in Egypt. Three of these locations are from mutually isolated pockets in the arid biotope of Southern Sinai (Wadi Sahab, El-Agramia and Rahaba plains). The fourth population was sampled from the semiarid biotope of Western Mediterranean Costal Desert (WMCD). Using reducing gel electrophoresis (SDS-PAGE), we have shown biotope-specific variation in the expression of peptides from scorpions collected from these distinct areas. WMCD sourced venom samples contain higher molecular weight protein components (219, 200, 170, 139, 116 kDa) than Southern Sinai scorpion venom samples. The Southern Sinai venom is characterized by the presence of 11 protein bands (93-0.58 kDa) that are not mirrored in the individual venom samples of WMCD. Bands of 33 and 3.4 kDa were characteristics of all individual venom samples of the scorpion populations. Even within Southern Sinai area, Sahab venom contains five fractions that are not detected in both El-Agramia and Rahaba venom samples. Moreover, male and female venom analysis revealed some sex-related proteomic similarities and differences between scorpion populations. Female venom appears to be more complicated than the male venom. Female venom samples showed bands of 219, 200, 77.5, 55.5, 45, 39, 37, 24 and 16 kDa which were absent in the male venom. The random amplified polymorphic DNA (RAPD) technique was used to estimate the genetic distance between the four scorpion populations. The RAPD data confirmed the genetic diversity at molecular level among the sampled populations. More than 77 RAPD bands (ranging in size from 125 to 15,000 bp) were defined from the four scorpion populations. Of the 77 bands, 57 (76.2%) were polymorphic and 20 were monomorphic among the populations. The similarity coefficient data of venom and DNA were used to construct separate dendrograms, which grouped together the Southern Sinai populations and these were some distance away from the WMCD population. Taken together, we suspect that a combination of local environmental conditions, geographical separation and genetic separation may play a major role in the intraspecific variation of venom of S. m. palmatus.

Neutralization of lethality and proteolytic activities of Malayan pit viper (Calloselasma rhodostoma) venom with North American Virginia opossum (Didelphis virginiana) serum

Malayan pit viper (Calloselasma rhodostoma) envenomation is a major health problem in South East Asia. During envenomation, venom components mainly affect the hemostatic system. The sera from the North American Virginia opossums (Didelphis virginiana) were able to neutralize the venom of the Malayan pit viper. These natural inhibitors could be explored as potential therapeutics against envenomations of a variety of venomous snake species in different geographical habitats.

Reappraisal of Vipera aspis venom neurotoxicity

Background

The variation of venom composition with geography is an important aspect of intraspecific variability in the Vipera genus, although causes of this variability remain unclear. The diversity of snake venom is important both for our understanding of venomous snake evolution and for the preparation of relevant antivenoms to treat envenomations. A geographic intraspecific variation in snake venom composition was recently reported for Vipera aspis aspis venom in France. Since 1992, cases of human envenomation after Vipera aspis aspis bites in south-east France involving unexpected neurological signs were regularly reported. The presence of genes encoding PLA₂ neurotoxins in the Vaa snake genome led us to investigate any neurological symptom associated with snake bites in other regions of France and in neighboring countries. In parallel, we used several approaches to characterize the venom PLA₂ composition of the snakes captured in the same areas.

Methodology/Principal Findings

We conducted an epidemiological survey of snake bites in various regions of France. In parallel, we carried out the analysis of the genes and the transcripts encoding venom PLA₂s. We used SELDI technology to study the diversity of PLA₂ in various venom samples. Neurological signs (mainly cranial nerve disturbances) were reported after snake bites in three regions of France: Languedoc-Roussillon, Midi-Pyrénées and Provence-Alpes-Côte d'Azur. Genomes of Vipera aspis snakes from south-east France were shown to contain ammodytoxin isoforms never described in the genome of Vipera aspis from other French regions. Surprisingly, transcripts encoding venom neurotoxic PLA₂s were found in snakes of Massif Central region. Accordingly, SELDI analysis of PLA₂ venom composition confirmed the existence of population of neurotoxic Vipera aspis snakes in the west part of the Massif Central mountains.

Conclusions/Significance

The association of epidemiological studies to genetic, biochemical and immunochemical analyses of snake venoms allowed a good evaluation of the potential neurotoxicity of snake bites. A correlation was found between the expression of neurological symptoms in humans and the intensity of the cross-reaction of venoms with anti-ammodytoxin antibodies, which is correlated with the level of neurotoxin (vaspin and/or ammodytoxin) expression in the venom. The origin of the two recently identified neurotoxic snake populations is discussed according to venom PLA₂ genome and transcriptome data.

Neuromuscular and phospholipase activities of venoms from three subspecies of Bothrops neuwiedi (B. n. goyazensis, B. n. paranaensis and B. n. diporus)

The Bothrops neuwiedi (Neuwied's lancehead) species complex consists of a variety of subspecies with a wide distribution in South America. In this work, we compared the neuromuscular blockade caused by venoms from three subspecies (B. n. goyazensis, B. n. paranaensis and B. n. diporus) of this complex using chick biventer cervicis (BC) and mouse phrenic nerve-diaphragm (PND) preparations and investigated their phospholipase A2 (PLA2) activities and electrophoretic profiles. The order of potency of PLA2 activity was B. n. diporus>B. n. paranaensis>B. n. goyazensis. In BC preparations, B. n. goyazensis venom (50 microg/mL) was significantly (p<0.05) more active than B. n. paranaensis and B. n. diporus venoms, which did not produce a significant blockade at this time interval; after 120 min, B. n. goyazensis, B. n. paranaensis and B. n. diporus venoms (100 microg/mL) produced blockades of 57.4+/-5%, 30+/-3% and 17.4+/-7% (n=3-6 each), respectively. The three venoms inhibited contractures in response to ACh, indicating interference with postsynaptic neurotransmission. Only B. n. goyazensis and B. n. paranaensis venoms caused a long-lasting, concentration-dependent muscle contracture prior to blockade. In PND preparations, all of the venoms blocked the twitch-tension responses within 45-100 min, indicating that these preparations were more sensitive than avian preparations. There was a correlation between PLA2 activity and the time for 50% blockade in PND but not in BC preparations. SDS-PAGE showed quantitative rather than qualitative differences among the venoms. These results indicate that the venoms of the three subspecies had similar profiles of neuromuscular activity, although the relationship with PLA2 activity varied with the preparation used.

Venom ontogeny, diet and morphology in Carukia barnesi, a species of Australian box jellyfish that causes Irukandji syndrome

Venom profiles of two age groups of the medically important Australian box jellyfish Carukia barnesi [Southcott, R.V., 1967. Revision of some Carybdeidae (Scyphozoa, Cubomedusae), including description of jellyfish responsible for the 'Irukandji' syndrome. Aust. J. Zool. 15, 651-657] were compared. Sodium dodecyl sulphate-polyacrylamide gel electrophoresis revealed differences in protein banding of tentacular venom between immature and mature animals. This correlates to a change in diet from invertebrate prey in immature C. barnesi medusae to vertebrate prey in mature medusae. Unlike other cubozoan studies, a change in venom did not equate to a change in nematocyst types or their relative frequencies. Additionally, comparison of tentacle structure and bell wart number showed developmental differences between the two age classes. Observations of prey capture in mature individuals and differences in bell warts between immature and mature medusae suggest different methods of prey capture are employed at different life stages of C. barnesi.

Mass spectrometric analysis of the individual variability of Bothrops jararaca venom peptide fraction. Evidence for sex-based variation among the bradykinin-potentiating peptides

Variation in the snake venom proteome is well documented and it is a ubiquitous phenomenon at all taxonomical levels. However, variation in the snake venom peptidome is so far not described. In this work we used mass spectrometry [liquid chromatography/mass spectrometry (LC/MS) and matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOFMS)] to explore sex-based differences among the venom peptides of eighteen sibling specimens of Bothrops jararaca of a single litter born and raised in the laboratory. MALDI-TOFMS analyses showed individual variability among the bradykinin-potentiating peptides (BPPs), and, interestingly, four new peptides were detected only in female venoms and identified by de novo sequencing as cleaved BPPs lacking the C-terminal Q-I-P-P sequence. Similar results were obtained with venom from wild-caught adult non-sibling specimens of B. jararaca and in this case we were able to identify the gender of the specimen by analyzing the MALDI-TOF profile of the peptide fraction and finding the cleaved peptides only in female venoms. Synthetic replicates of the cleaved BPPs were less potent than the full-length BPP-10c in potentiating the bradykinin hypotensive effect, suggesting that the C-terminus is critical for the interaction of the BPPs with their mammalian molecular targets. This work represents a comprehensive mass spectrometric analysis of the peptide fraction of B. jararaca venom and shows for the first time sex-based differences in the snake venom peptidome of sibling and non-sibling snakes and suggests that the BPPs may follow distinct processing pathways in female and male individuals.

Venom Proteomes of Closely RelatedSistrurusRattlesnakes with Divergent Diets

The protein composition of the venoms of the three subspecies of Sistrurus catenatus (S. c. catenatus, tergeminus, and edwardsii) and a basal species, Sistrurus miliarius barbouri, were analyzed by RP-HPLC, N-terminal sequencing, MALDI-TOF peptide mass fingerprinting, and CID-MS/MS. The venoms of the four Sistrurus taxa contain proteins from 11 families. The protein family profile and the relative abundance of each protein group in the different venoms are not conserved. Myotoxins and 2-chain PLA2s were detected only in S.c. catenatus and S.c. tergeminus, whereas C-type BPP and Kunitz-type inhibitors were exclusively found in S.c. edwardsii and Sistrurus miliarius barbouri. Among major protein families, taxa were most similar in their metalloproteases (protein similarity coefficient value: 34%) and most divergent in PLA2s (12%), with values for disintegrins and serine proteases lying between these extremes (25 and 20%, respectively). The patterns of venom diversity points to either a gain in complexity in S. catenatus taxa or a loss of venom diversity occurring early on in the evolution of the group involving the lineage connecting S. milarius to the other taxa. The high degree of differentiation in the venom proteome among recently evolved congeneric taxa emphasizes the uniqueness of the venom composition of even closely related species that have different diets. Comparative proteomic analysis of Sistrurus venoms provides a comprehensive catalog of secreted proteins, which may contribute to a deeper understanding of the biology and ecology of these North American snakes and may also serve as a starting point for studying structure-function correlations of individual toxins.