Phylogeny-Based Comparative Analysis of Venom Proteome Variation in a Clade of Rattlesnakes (Sistrurus sp.)

A long-standing question in evolutionary studies of snake venoms is the extent to which phylogenetic divergence and diet can account for between-species differences in venom composition. Here we apply phylogeny-based comparative methods to address this question. We use data on venom variation generated using proteomic techniques for all members of a small clade of rattlesnakes (Sistrurus sp.) and two outgroups for which phylogenetic and diet information is available. We first complete the characterization of venom variation for all members of this clade with a “venomic” analysis of pooled venoms from two members of this genus, S. milarius streckeri and S. m. milarius. These venoms exhibit the same general classes of proteins as those found in other Sistrurus species but differ in their relative abundances of specific protein families. We then test whether there is significant phylogenetic signal in the relative abundances of major venom proteins across species and if diet (measured as percent mammals and lizards among all prey consumed) covaries with venom composition after phylogenetic divergence is accounted for. We found no evidence for significant phylogenetic signal in venom variation: K values for seven snake venom proteins and two composite venom variables [PC 1 and 2]) were all nonsignificant and lower (mean = 0.11+0.06 sd) than mean K values (>0.35) previously reported for a wide range of morphological, life history, physiological and behavioral traits from other species. Finally, analyses based on Phylogenetic Generalized Least Squares (PGLS) methods reveal that variation in abundance of some venom proteins, most strongly CRISP is significantly related to snake diet. Our results demonstrate that venom variation in these snakes is evolutionarily a highly labile trait even among very closely-related taxa and that natural selection acting through diet variation may play a role in molding the relative abundance of specific venom proteins.     

Eggs-Only Diet: Its Implications for the Toxin Profile Changes and Ecology of the Marbled Sea Snake (Aipysurus eydouxii)

Studies so far have correlated the variation in the composition of snake venoms with the target prey population and snakes diet. Here we present the first example of an alternative evolutionary link between venom composition and dietary adaptation of snakes. We describe a dinucleotide deletion in the only three finger toxin gene expressed in the sea snake Aipysurus eydouxii (Marbled Sea Snake) venom and how it may have been the result of a significant change in dietary habits. The deletion leads to a frame shift and truncation with an accompanying loss of neurotoxicity. Due to the remarkable streamlining of sea snake venoms, a mutation of a single toxin can have dramatic effects on the whole venom, in this case likely explaining the 50- to 100-fold decrease in venom toxicity in comparison to that of other species in the same genus. This is a secondary result of the adaptation of A. eydouxii to a new dietary habit — feeding exclusively on fish eggs and, thus, the snake no longer using its venom for prey capture. This was parallel to greatly atrophied venom glands and loss of effective fangs. It is interesting to note that a potent venom was not maintained for use in defense, thus reinforcing that the primary use of snake venom is for prey capture.Nucleotide sequence data reported here have been deposited in the GenBank database under accession number AY559317.Reviewing Editor: Dr. Martin Kreitman     

The molecular basis of venom resistance in a rattlesnake‐squirrel predator‐prey system

Understanding how interspecific interactions mould the molecular basis of adaptations in coevolving species is a long‐sought goal of evolutionary biology. Venom in predators and venom resistance proteins in prey are coevolving molecular phenotypes, and while venoms are highly complex mixtures it is unclear if prey respond with equally complex resistance traits. Here, we use a novel molecular methodology based on protein affinity columns to capture and identify candidate blood serum resistance proteins (“venom interactive proteins” [VIPs]) in California Ground Squirrels (Otospermophilus beecheyi) that interact with venom proteins from their main predator, Northern Pacific Rattlesnakes (Crotalus o. oreganus). This assay showed that serum‐based resistance is both population‐ and species‐specific, with serum proteins from ground squirrels showing higher binding affinities for venom proteins of local snakes compared to allopatric individuals. Venom protein specificity assays identified numerous and diverse candidate prey resistance VIPs but also potential targets of venom in prey tissues. Many specific VIPs bind to multiple snake venom proteins and, conversely, single venom proteins bind multiple VIPs, demonstrating that a portion of the squirrel blood serum “resistome” involves broad‐based inhibition of nonself proteins and suggests that resistance involves a toxin scavenging mechanism. Analyses of rates of evolution of VIP protein homologues in related mammals show that most of these proteins evolve under purifying selection possibly due to molecular constraints that limit the evolutionary responses of prey to rapidly evolving snake venom proteins. Our method represents a general approach to identify specific proteins involved in co‐evolutionary interactions between species at the molecular level.     

Identification and characterization of a taxon-specific three-finger toxin from the venom of the Green Vinesnake (Oxybelis fulgidus; family Colubridae)

Snake venoms contain a variety of protein and peptide toxins, and the three-finger toxins (3FTxs) are among the best characterized family of venom proteins. The compact nature and highly conserved molecular fold of 3FTxs, together with their abundance in many venoms, has contributed to their utility in structure-function studies. Although many target the nicotinic acetylcholine receptor of vertebrate skeletal muscle, often binding with nanomolar K_ds, several non-conventional 3FTxs show pronounced taxon-specific neurotoxic effects. Here we describe the purification and characterization of fulgimotoxin, a monomeric 3FTx from the venom of Oxybelis fulgidus, a neotropical rear-fanged snake. Fulgimotoxin retains the canonical 5 disulfides of the non-conventional 3FTxs and is highly neurotoxic to lizards; however, mice are unaffected, demonstrating that this toxin is taxon-specific in its effects. Analysis of structural features of fulgimotoxin and other colubrid venom 3FTxs indicate the presence of a “colubrid toxin motif” (CYTLY) and a second conserved segment (WAVK) found in Boiga and Oxybelis taxon-specific 3FTxs, both in loop II. Because specific residues in loop II conventional α-neurotoxic 3FTxs are intimately associated with receptor binding, we hypothesize that this loop, with its highly conserved substitutions, confers taxon-specific neurotoxicity. These findings underscore the importance of rear-fanged snake venoms for understanding the evolution of toxin molecules and demonstrate that even among well-characterized toxin families, novel structural and functional motifs may be found.     

Snake venomics: comparative analysis of the venom proteomes of the Tunisian snakes Cerastes cerastes, Cerastes vipera and Macrovipera lebetina

The protein composition of the crude venoms of the three most important vipers of Tunisia was analyzed by RP-HPLC, N-terminal sequence analysis, MALDI-TOF mass determination, and in-gel tryptic digestion followed by PMF and CID-MS/MS of selected peptide ions in a quadrupole-linear IT instrument. Our results show that the venom proteomes of Cerastes cerastes, Cerastes vipera, and Macrovipera lebetina are composed of proteins belonging to a few protein families. However, each venom showed distinct degree of protein composition complexity. The three venoms shared a number of protein classes though the relative occurrence of these toxins was different in each snake species. On the other hand, the venoms of the Cerastes species and Macrovipera lebetina each contained unique components. The comparative proteomic analysis of Tunisian snake venoms provides a comprehensible catalogue of secreted proteins, which may contribute to a deeper understanding of the biological effects of the venoms, and may also serve as a starting point for studying structure-function correlations of individual toxins.     

Rapid Evolution by Positive Selection and Gene Gain and Loss: PLA<Subscript>2</Subscript> Venom Genes in Closely Related <Emphasis Type="Italic">Sistrurus</Emphasis> Rattlesnakes with Divergent Diets

Rapid evolution of snake venom genes by positive selection has been reported previously but key features of this process such as the targets of selection, rates of gene turnover, and functional diversity of toxins generated remain unclear. This is especially true for closely related species with divergent diets. We describe the evolution of PLA₂ gene sequences isolated from genomic DNA from four taxa of Sistrurus rattlesnakes which feed on different prey. We identified four to seven distinct PLA₂ sequences in each taxon and phylogenetic analyses suggest that these sequences represent a rapidly evolving gene family consisting of both paralogous and homologous loci with high rates of gene gain and loss. Strong positive selection was implicated as a driving force in the evolution of these protein coding sequences. Exons coding for amino acids that make up mature proteins have levels of variation two to three times greater than those of the surrounding noncoding intronic sequences. Maximum likelihood models of coding sequence evolution reveal that a high proportion (∼30%) of all codons in the mature protein fall into a class of codons with an estimated d _N /d _S (ω) ratio of at least 2.8. An analysis of selection on individual codons identified nine residues as being under strong (p < 0.01) positive selection, with a disproportionately high proportion of these residues found in two functional regions of the PLA₂ protein (surface residues and putative anticoagulant region). This is direct evidence that diversifying selection has led to high levels of functional diversity due to structural differences in proteins among these snakes. Overall, our results demonstrate that both gene gain and loss and protein sequence evolution via positive selection are important evolutionary forces driving adaptive divergence in venom proteins in closely related species of venomous snakes.     

Identification and phylogeny of Arabian snakes: Comparison of venom chromatographic profiles versus 16S rRNA gene sequences

Identification of snake species is important for various reasons including the emergency treatment of snake bite victims. We present a simple method for identification of six snake species using the gel filtration chromatographic profiles of their venoms. The venoms of Echis coloratus, Echis pyramidum, Cerastes gasperettii, Bitis arietans, Naja arabica, and Walterinnesia aegyptia were milked, lyophilized, diluted and centrifuged to separate the mucus from the venom. The clear supernatants were filtered and chromatographed on fast protein liquid chromatography (FPLC). We obtained the 16S rRNA gene sequences of the above species and performed phylogenetic analysis using the neighbor-joining method. The chromatograms of venoms from different snake species showed peculiar patterns based on the number and location of peaks. The dendrograms generated from similarity matrix based on the presence/absence of particular chromatographic peaks clearly differentiated Elapids from Viperids. Molecular cladistics using 16S rRNA gene sequences resulted in jumping clades while separating the members of these two families. These findings suggest that chromatographic profiles of snake venoms may provide a simple and reproducible chemical fingerprinting method for quick identification of snake species. However, the validation of this methodology requires further studies on large number of specimens from within and across species.     

Inhibition of Nicotinic Acetylcholine Receptors,a Novel Facet in the Pleiotropic Activities of Snake Venom Phospholipases A2

Phospholipases A₂ represent the most abundant family of snake venom proteins. They manifest an array of biological activities, which is constantly expanding. We have recently shown that a protein bitanarin, isolated from the venom of the puff adder Bitis arietans and possessing high phospholipolytic activity, interacts with different types of nicotinic acetylcholine receptors and with the acetylcholine-binding protein. To check if this property is characteristic to all venom phospholipases A₂, we have studied the capability of these enzymes from other snakes to block the responses of Lymnaea stagnalis neurons to acetylcholine or cytisine and to inhibit α-bungarotoxin binding to nicotinic acetylcholine receptors and acetylcholine-binding proteins. Here we present the evidence that phospholipases A₂ from venoms of vipers Vipera ursinii and V. nikolskii, cobra Naja kaouthia, and krait Bungarus fasciatus from different snake families suppress the acetylcholine- or cytisine-elicited currents in L. stagnalis neurons and compete with α-bungarotoxin for binding to muscle- and neuronal α7-types of nicotinic acetylcholine receptor, as well as to acetylcholine-binding proteins. As the phospholipase A₂ content in venoms is quite high, under some conditions the activity found may contribute to the deleterious venom effects. The results obtained suggest that the ability to interact with nicotinic acetylcholine receptors may be a general property of snake venom phospholipases A₂, which add a new target to the numerous activities of these enzymes.     

Extraction and protein component analysis of venom from the dissected venom glands of Latrodectus tredecimguttatus

Black widow spiders (genus Latrodectus) have attracted increasing attention due to frequently reported human injuries caused by them and the potential applications of biologically active components in their venoms. Although a number of studies have described the biological properties and structures of several venomous proteins such as latrotoxins, a comprehensive analysis of protein component of the venom from the spider is not available. We used combinative proteomic strategies to assess the protein components of the crude venom collected from Latrodectus tredecimguttatus by extracting the dissected venom glands. The experiments demonstrated that the crude venom of L. tredecimguttatus has a high abundance of acidic proteins with molecular masses greater than 15 kDa, and the content of proteins and peptides of below 15 kDa is low. 86 unique proteins were identified, part of which were contaminations of cellular components during the extraction, determined in comparison with venom obtained by electrostimulation. Except for members of latrotoxin family that were commonly considered as the primary toxic components of the venom, several other special enzymes and proteins were detected such as protease, phosphatase, lysozyme, inhibitory protein, and so on. These protein components, particularly the proteases, were speculated to play important roles in the action of L. tredecimguttatus venom.     

Cloning and characterisation of novel cystatins from elapid snake venom glands

Snake venoms contain a complex mixture of polypeptides that modulate prey homeostatic mechanisms through highly specific and targeted interactions. In this study we have identified and characterised cystatin-like cysteine-protease inhibitors from elapid snake venoms for the first time. Novel cystatin sequences were cloned from 12 of 13 elapid snake venom glands and the protein was detected, albeit at very low levels, in a total of 22 venoms. One highly conserved isoform, which displayed close sequence identity with family 2 cystatins, was detected in each elapid snake. Crude Austrelaps superbus (Australian lowland copperhead) snake venom inhibited papain, and a recombinant form of A. superbus cystatin inhibited cathepsin L ≅ papain > cathepsin B, with no inhibition observed for calpain or legumain. While snake venom cystatins have truncated N-termini, sequence alignment and structural modelling suggested that the evolutionarily conserved Gly-11 of family 2 cystatins, essential for cysteine protease inhibition, is conserved in snake venom cystatins as Gly-3. This was confirmed by mutagenesis at the Gly-3 site, which increased the dissociation constant for papain by 10⁴-fold. These data demonstrate that elapid snake venom cystatins are novel members of the type 2 family. The widespread, low level expression of type 2 cystatins in snake venom, as well as the presence of only one highly conserved isoform in each species, imply essential housekeeping or regulatory roles for these proteins.     

Intraspecific variation in the venom of the vermivorous cone snail Conus vexillum

A combination of proteomic and biochemical assays was used to examine variations in the venom of Conus vexillum taken from two locations (Hurgada and Sharm El-Shaikh) in the Red Sea, Egypt. Using MALDI/TOF-MS, a remarkable degree of intra-species variation between venom samples from both locations was identified. To evaluate variability in the cytotoxic effects of Conus venom, mice were injected with the same dose from each location. The oxidative stress biomarkers [malondialdehyde (MDA), protein carbonyl content (PCC)], antioxidants [glutathione (GSH), superoxide dismutase (SOD), catalase (CAT)], total antioxidant capacity (TAC) and nitric oxide (NO), were measured 3, 6, 9 and 12 h post venom injection. The venoms induced a significant increase in the levels of PCC, MDA, NO, GSH and CAT. The venoms significantly inhibited the activity of SOD and reduced the TAC. Toxicological data showed that the venom obtained from Hurgada was more potent than that obtained from Sharm El-Shaikh. It can be concluded that: (1) the venom of the same Conus species from different regions is highly diversified (2) the venoms from different locations reflect clear differences in venom potency and (3) the cytotoxic effects of C. vexillum venom can be attributed to its ability to induce oxidative stress.     

Phenotypic integration in the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus)

Selection can vary geographically across environments and temporally over the lifetime of an individual. Unlike geographic contexts, where different selective regimes can act on different alleles, age‐specific selection is constrained to act on the same genome by altering age‐specific expression. Snake venoms are exceptional traits for studying ontogeny because toxin expression variation directly changes the phenotype; relative amounts of venom components determine, in part, venom efficacy. Phenotypic integration is the dependent relationship between different traits that collectively produce a complex phenotype and, in venomous snakes, may include traits as diverse as venom, head shape and fang length. We examined the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus) across environments and over the lifetime of individuals and used a genotype–phenotype map approach, protein expression data and morphological data to demonstrate that: (i) ontogenetic effects explained more of the variation in toxin expression variation than geographic effects, (ii) both juveniles and adults varied geographically, (iii) toxin expression variation was a result of directional selection and (iv) different venom phenotypes covaried with morphological traits also associated with feeding in temporal (ontogenetic) and geographic (functional) contexts. These data are the first to demonstrate, to our knowledge, phenotypic integration between multiple morphological characters and a biochemical phenotype across populations and age classes. We identified copy number variation as the mechanism driving the difference in the venom phenotype associated with these morphological differences, and the parallel mitochondrial, venom and morphological divergence between northern and southern clades suggests that each clade may warrant classification as a separate evolutionarily significant unit.     

Snake venom components and their cross-reactivity: a review

Snake venoms are complex mixtures of organic and inorganic compounds, many of which display biological activity. It has been demonstrated that antisera raised against whole venom or a single purified venom protein from one species of snake will react with proteins in the venom of other species. This cross-reactivity between species may have applications in determining snake phylogeny, but recent studies on the variation of venom components within a species make these evolutionary conclusions questionable.     

Venom proteomic analysis of medically important Nigerian viper Echis ocellatus and Bitis arietans snake species

Snakebite envenoming remains a neglected tropical disease which poses severe health hazard, especially for the rural inhabitants in Africa. In Nigeria, vipers are responsible for the highest number of deaths. Hydrophilic interaction liquid chromatography coupled with LC-MS/MS was used to analyze the crude venoms of Echis ocellatus (Carpet viper) and Bitis arietans (Puff adder) in order to understand their venom proteomic identities. Results obtained revealed that gel-free proteomic analysis of the crude venoms led to the identification of 85 and 79 proteins, respectively. Seventy-eight (78) proteins were common between the two snake species with a 91.8% similarity score. The identified proteins belong to 18 protein families in E. ocellatus and 14 protein families in B. arietans. Serine proteases (22.31%) and metalloproteinases (21.06%) were the dominant proteins in the venom of B. arietans; while metalloproteinases (34.84%), phospholipase A₂s (21.19%) and serine proteases (15.50%) represent the major toxins in the E. ocellatus venom. Other protein families such as three-finger toxins and cysteine-rich venom proteins were detected in low proportions. This study provides an insight into the venom proteomic analysis of the two Nigerian viper species, which could be useful in identifying the toxin families to be neutralized in case of envenomation.     

Evaluation of cytotoxic activities of snake venoms toward breast (MCF-7) and skin cancer (A-375) cell lines

Snake venoms are mixtures of bioactive proteins and peptides that exhibit diverse biochemical activities. This wide array of pharmacologies associated with snake venoms has made them attractive sources for research into potentially novel therapeutics, and several venom-derived drugs are now in use. In the current study we performed a broad screen of a variety of venoms (61 taxa) from the major venomous snake families (Viperidae, Elapidae and “Colubridae”) in order to examine cytotoxic effects toward MCF-7 breast cancer cells and A-375 melanoma cells. MTT cell viability assays of cancer cells incubated with crude venoms revealed that most venoms showed significant cytotoxicity. We further investigated venom from the Red-bellied Blacksnake (Pseudechis porphyriacus); venom was fractionated by ion exchange fast protein liquid chromatography and several cytotoxic components were isolated. SDS-PAGE and MALDI-TOF mass spectrometry were used to identify the compounds in this venom responsible for the cytotoxic effects. In general, viper venoms were potently cytotoxic, with MCF-7 cells showing greater sensitivity, while elapid and colubrid venoms were much less toxic; notable exceptions included the elapid genera Micrurus, Naja and Pseudechis, which were quite cytotoxic to both cell lines. However, venoms with the most potent cytotoxicity were often not those with low mouse LD₅₀s, including some dangerously venomous viperids and Australian elapids. This study confirmed that many venoms contain cytotoxic compounds, including catalytic PLA₂s, and several venoms also showed significant differential toxicity toward the two cancer cell lines. Our results indicate that several previously uncharacterized venoms could contain promising lead compounds for drug development.     

Comparison of Phylogeny,Venom Composition and Neutralization by Antivenom in Diverse Species of Bothrops Complex

In Latin America, Bothrops snakes account for most snake bites in humans, and the recommended treatment is administration of multispecific Bothrops antivenom (SAB – soro antibotrópico). However, Bothrops snakes are very diverse with regard to their venom composition, which raises the issue of which venoms should be used as immunizing antigens for the production of pan-specific Bothrops antivenoms. In this study, we simultaneously compared the composition and reactivity with SAB of venoms collected from six species of snakes, distributed in pairs from three distinct phylogenetic clades: Bothrops, Bothropoides and Rhinocerophis. We also evaluated the neutralization of Bothrops atrox venom, which is the species responsible for most snake bites in the Amazon region, but not included in the immunization antigen mixture used to produce SAB. Using mass spectrometric and chromatographic approaches, we observed a lack of similarity in protein composition between the venoms from closely related snakes and a high similarity between the venoms of phylogenetically more distant snakes, suggesting little connection between taxonomic position and venom composition. P-III snake venom metalloproteinases (SVMPs) are the most antigenic toxins in the venoms of snakes from the Bothrops complex, whereas class P-I SVMPs, snake venom serine proteinases and phospholipases A₂ reacted with antibodies in lower levels. Low molecular size toxins, such as disintegrins and bradykinin-potentiating peptides, were poorly antigenic. Toxins from the same protein family showed antigenic cross-reactivity among venoms from different species; SAB was efficient in neutralizing the B. atrox venom major toxins. Thus, we suggest that it is possible to obtain pan-specific effective antivenoms for Bothrops envenomations through immunization with venoms from only a few species of snakes, if these venoms contain protein classes that are representative of all species to which the antivenom is targeted.