Next Article in Journal
Cells Deficient in the Fanconi Anemia Protein FANCD2 are Hypersensitive to the Cytotoxicity and DNA Damage Induced by Coffee and Caffeic Acid
Next Article in Special Issue
A Tricky Trait: Applying the Fruits of the “Function Debate” in the Philosophy of Biology to the “Venom Debate” in the Science of Toxinology
Previous Article in Journal
Distribution and Metabolism of Bt-Cry1Ac Toxin in Tissues and Organs of the Cotton Bollworm, Helicoverpa armigera
Previous Article in Special Issue
Tempo and Mode of the Evolution of Venom and Poison in Tetrapods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Venom Evolution Lab, School of Biological Sciences, University of Queensland, St Lucia, QLD 4072, Australia
2
Department of Earth and Biological Sciences, Loma Linda University, Loma Linda, CA 92350, USA
3
Department of Biochemistry & Molecular Biology, Biomedical Discovery Institute, Monash University, Clayton, VIC 3800, Australia
4
School of Chemistry and Molecular Biosciences, University of Queensland, St Lucia, QLD 4072, Australia
5
Baker IDI Heart and Diabetes Institute, 75 Commercial Road, Melbourne, Victoria 3004, Australia
6
Department of Physiology, Biomedical Discovery Institute, Monash University, Clayton, VIC 3800, Australia
7
Aquatic Animal Health, School of Biological Sciences, University of Queensland, St Lucia, QLD 4072 Australia
8
Institute for Molecular Bioscience, University of Queensland, St Lucia, QLD 4072, Australia
9
HEJ Research Institute of Chemistry, ICCBS, University of Karachi, Karachi-75270, Pakistan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Toxins 2016, 8(7), 210; https://doi.org/10.3390/toxins8070210
Submission received: 22 August 2015 / Revised: 28 May 2016 / Accepted: 16 June 2016 / Published: 8 July 2016
(This article belongs to the Collection Evolution of Venom Systems)

Abstract

:
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.

1. Introduction

Snake venoms are complex secretions of various enzymes, peptides and small organic molecules that have diverse modes of action on both prey and accidental human bite victims [1,2,3,4,5,6]. Though targeting a wide array of physiological processes and targets, the proteins that make up snake venom come from relatively few protein families. This is due to core protein scaffolds being maintained with accelerated rates of evolution at the secondary through quaternary levels facilitating the evolution of multiple derived functions within each protein family [3,7,8,9]. The use of mass spectrometry in venom proteomics [7,10,11,12,13,14,15,16,17,18,19] and venom gland transcriptome analysis [9,20,21,22,23,24,25,26,27,28] has improved knowledge of venom composition tremendously. Variation in venom profiles has been shown interspecifically [7,13,18,19,29,30,31,32,33,34] and intraspecifically, with intraspecific differences found among geographic locales [12,35,36,37,38,39,40,41,42,43], between sexes [36,37,44,45], and between age classes [10,15,46,47,48,49,50,51,52]. Some authors have argued that diversity in venom composition is the result of neutral evolutionary processes and is not subject to natural selection [53,54], whereas others have argued that adaptation has been driven, by strong natural selection, to specific prey species [19,29,36,37,55,56,57,58,59,60].
New World pitvipers are believed to be descendants of an ancestral Asian pitviper species, by way of a single New World colonization event, via the Bering land bridge [61,62,63,64]. In regards to the current study particular interest is given to the members of the genera Bothriechis and Bothrops that have independently evolved to occupy arboreal niches. That Bothriechis is a distinctive genus from the “bothropoid” genera [61] has been asserted by multiple authors [61,62,63,64,65,66,67]; however, phylogenetic relationships within the bothropoids, in particular the validity of Bothriopsis, have proven more difficult due to limited sampling of taxa. Numerous studies have presented evidence of the paraphyly of Bothrops with respect to Bothriopsis [61,62,63,64,65,66,67,68,69,70,71,72] leading some authors to advocate for the synonymizing of Bothriopsis and Bothrocophias with Bothrops [70,72] or Bothriechis [73], while others have maintained that Bothriopsis and Bothrocophias should remain distinct [61,62,65,74,75] (Figure 1). Regardless of taxonomical divisions preferred, it is clear that the arboreal specialization has occurred upon two separate occasions.
Recent studies utilizing a combined morphological/molecular approach and increased sampling of taxa have disagreed on generic assignment as well. Fenwick et al. (2009), sampling 90% of the known “bothropoid” taxa and utilizing both morphological characters and mtDNA sequence data found support for several monophyletic clades within Bothrops (sensu lato) and proposed taxonomic revisions, including recognition of Bothriopsis, based on the distinct species groups recovered [76]. Carrasco et al. (2012) [77] incorporated additional morphological characters not utilized by Fenwick et al. (2009) [76] in addition to previously published molecular sequences and ecological characters. Citing possible paraphyly of Bothrops with respect to Bothriopsis (with the inclusion of Bothrops sanctaecrucis), they proposed following the recommendation of Salomão et al. (1997) [70] and Wüster et al. (2002) [72] and synonymized Bothriopsis with Bothrops [77]. We have followed the proposals of Salomão et al. (1997) [70] and Wüster et al. (2002) [72] in recognizing the bothropoids as a single genus, Bothrops, as the weight of current evidence (cited above) favors this classification.
Despite the medical significance of pitviper envenomations in the New World [78,79,80,81,82,83], the majority of toxinological research published to date has been on North American members of the genus Crotalus [14,25,26,35,38,39,40,41,42,43,49,50,82,83,84,85,86,87,88,89,90], and the larger members of the genus Bothrops [10,91,92,93,94,95,96,97,98,99,100,101,102,103,104]. While investigations into Bothriechis, in particular B. schlegelii, venom [31,105,106,107] and some of the smaller bothropoid pitvipers [108,109,110,111] have occurred, knowledge gaps remain on the proteomic make up of the venom of the majority of New World pitviper species. This knowledge gap is unfortunate as snake venoms have proved valuable sources of investigational ligands [112,113,114,115,116] and potential therapeutics [117,118,119,120,121,122,123,124] due to the high degree of target specificity exhibited by their constituent toxins. Bothrops and Bothriechis are both medically relevant species across their range and together are responsible for a large number of severe envenomations [78,79,80,81,125,126]. American pitviper venoms are typically rich in two types of enzymes: Snake venom metalloprotease (SVMP) and phospholipase A2 (PLA2) [3,42,127]. Due to variance in retention of ancestral domains, SVMPs span a broad range of molecular weights from ~30–70 kDa and exhibit a variety of activities. Structural diversity within SVMPs is the basis for the development of a variety of symptoms in victims of envenomation including; haemorrhage, edema, hypotension, hypovolemia, inflammation, and necrosis [128,129]. The realities of inter/intraspecific venom variation necessitate a thorough understanding of species venom proteomes as the information allows for improved clinical treatment of envenomations by way of informing antivenom selection for heterologous venoms [127,130]. Determining a species proteome is often the first step of antivenomic studies undertaken to confirm antivenom neutralization efficacy of species/populations/age classes [12,84,131] and determination of toxins unable to be neutralized is valuable information not only for clinicians in charge of managing envenomations but also for future antivenom production.
In this study we investigate the venom proteomes of arboreal members of the Bothriechis and Bothrops genera that have independently invaded the tree canopies of South America [7,31,64,71,77,107]. In both genera variation in adult size has been documented, with Bothriechis aurifer, Bothriechis lateralis and Bothriechis marchi typically exceeding Bothriechis schlegelli in adult size, while Bothrops bilineata is a much smaller species than Bothrops taeniata. For comparison, we have included two terrestrial Bothrops species that also differ markedly in maximum adult size, with Bothrops asper reaching adult sizes in excess of twice that of Bothrops neuwiedi bolivianus. While B. lateralis and B. schlegelli have been proteomically previously profiled, as has Bothriechis nigroviridis [31,132], their venom variation has not been considered in relation to size and prey preference.

2. Results

Shotgun mass spectrometry recovered proteins of known pitviper toxin types (Table 1) in agreement with previous proteomic [10,11,13,14,30,32,83,84,92] and transcriptomic analyses [22,24,29,59,77,89,90,93,95,97,103].
1D gel analysis revealed greater complexity in all venoms than indicated by the shotgun results. The venoms of the larger Bothriechis species (B. aurifer, B. lateralis and B. marchi) contained more P-III, P-II and P-I SVMP than the venom of the smaller sized species B. schlegelii (Figure 2). In contrast, B. schlegelii venom contained more PLA2 (lower molecular weights of synovial PLA2 12–15 kDa). A similar correlation between body size and SVMP/PLA2 content was found in arboreal Bothrops. However this trend was reversed in the terrestrial Bothrops. These results were complemented by the 2D gel analyses (Figure 3).
The relative straightforward proteomics trends were not reflected by the most complex differences in venom composition revealed by the bioactivity testing, which is consistent with venoms being complex mixtures of bioactive substances. While B. aurifer was indeed much more potent than B. schlegelli in a fluorescence-based metalloprotease activity assay (Figure 4), and B. lateralis also more active, B. marchi was only weakly active, being equipotent to B. schlegelli. This result was despite B. marchi having very similar 1D and 2D gel patterns to the more potent B. lateralis, especially in the known metalloprotease regions. It must be noted that these trends are of course specific only to this particular substrate and that other assays (such as the gel based zymography) may present differential results. Indeed, further investigation of metalloprotease activity yielded results that contrasted those gathered with the aforementioned analyses (Figure 4). Zymography gels with a casein substrate revealed active P-III SVMP in all four Bothriechis venoms in the 75 kDa region (Figure 2). The region of digestion in the casein zymography gel (Figure 5) corresponds not with the 50 kDa heavy bands seen in the 1D gels (Figure 2), but a 75 kDa region above that, with the 50 kDa 1D bands clearly evident in the zymography gels below the zone of digestion. This suggests that this digestive activity is due to the presence of heavily glycosylated enzymes that are resistant to staining, while the well-stained 50 kDa serine protease bands do not have this type of activity. In contrast, zymography gels with gelatin as the substrate produced bands of digestion in the Bothriechis venoms at the 50 kDa region (Figure 5), except for B. aurifer despite it having seemingly homologous serine protease band in the same region.
The arboreal Bothrops venoms were both weakly active on the casein substrate zymography gel (Figure 5). However, these two venoms differed sharply in the gelatin substrate zymography gel, with B. taeniata having two discrete zones of digestion, one of which was similar to that of the Bothriechis species, while the other was higher in molecular weight (Figure 5). Neither arboreal Bothrops venom, however, displayed significant activity in the fluorescence based assay. Both terrestrial Bothrops species were only weakly active in either the casein or gelatin (Figure 5) substrate zymography gels. However both were active in the fluorescence based metalloprotease assay, with B. asper displaying high levels of activity.
PLA2 activity levels, however, were more congruent with the relative presence of venom components in the corresponding range (12–15 kDa), with B. aurifer and B. marchi having only weak activity while B. lateralis and B. schlegelli were very active (Figure 6). This is consistent with the latter two species having darker staining bands at ~12–15 kDa than the former two species. Thus envenomations by B. lateralis and B. schlegelli would be expected to produce more myotoxicity than envenomations by B. aurifer and B. marchi. Despite Bothrops bilineata and B. taeniata both containing abundant PLA2, they were only weakly active in this region, suggesting that myotoxicity would not be as significant a complication but that other non-enzymatic PLA2 activities (e.g., antiplatelet aggregation or neurotoxicity) may be more potent in these two species. Similarly, despite both being rich in PLA2, the two terrestrial Bothrops species displayed significantly different levels of PLA2 enzymatic activity. Bothrops asper displayed high levels of PLA2 enzymatic activity, consistent with myotoxic envenomation effects, while Bothrops neuwiedi bolivianus was much less potent in this regard.

3. Discussion

For the Bothrops species tested, SVMP levels are correlated with consumed prey class percentages, reported by Martins et al. (2002) [132], with species containing higher percentage of anurans in their diet (B. bilineata and B. taeniata) possessing lower levels of SVMP activity than species that contained non-volant mammals as a higher percentage of their diet (B. asper and B. neuwiendi). Pitviper venoms have been dichotomously classified into type 1 (proteolytic or “tenderizer”) and type II (more toxic) venoms [133], as anurans have a high surface area to volume ratio and skin lacking tough dermal scales, it is unlikely that strong selection pressures for venoms with high metalloprotease activity exist in species that contain anurans as a large percentage of their diet as the need for such cleavage enzymes is reduced. The increased levels of SVMP activity for non-volant mammal generalist B. asper (59.4% mammals) in regards to B. neuwiedi (93.1% mammals) may be due to size class of mammal taken, as pitvipers are gape-limited predators that consume their prey whole [134,135,136]. Pitviper SVL (snout vent length) in the North American pitviper, Crotalus atrox, was found as the best predictor of overall maximum gape [137]. Although B. asper is recorded as consuming less percentage of mammals than B. neuwiedi, B. asper obtains an average SVL length roughly twice that of B. neuwiedi [131], suggesting the requirement for increased amounts of SVMP, as larger prey requires increased levels of SVMP which impairs homeostasis during predation. In addition our venom samples were taken from adult individuals displaying the predicted adult requirement of toxin activity. Other studies have documented that pitvipers consume prey species with lower surface area to volume ratios as they grow [50,132,138,139] and the lytic activities of SVMP may aid in predigestion of acquired prey. A detailed dietary analysis including prey species/size would be necessary in order to determine if B. asper does indeed consume larger prey/ prey of lower surface area to volume ratio than B. neuwiedi.
Detailed dietary studies for members of Bothriechis are lacking. As a whole the genus appears to be a generalist with an ontogenetic shift toward mammalian prey, with B. lateralis and B. schlegelii even documented as occasionally predating on bats and small birds (Campbell and Lamar 2004 and references therein) [131]. We detected increased levels of SVMP activity for B. aurifer in relation to other members of Bothriechis. This increased level of SVMP activity in combination with Campbell and Lamar (2004) [131] noting that the snake is frequently encountered on the ground is suggestive that B. aurifer may contain a higher percentage of non-volant mammals in its diet than the other members of Bothriechis tested. In contrast to high activity levels of SVMP recorded in B. aurifer, B. schleglii and B. lateralis both present higher levels of sPLA2 activity in relation to other members of Bothriechis. The increased levels of sPLA2 activity present in the venom profiles of B. lateralis and B. schleglii is suggestive of their arboreal encounters with volant animals. As sPLA2 is a neurotoxin affecting the nervous system, recruiting such a toxin would aid in quickly immobilizing flying prey whilst ambushing from a tree branch high up in the forest canopy, therefore reducing the risk of escape and injury.
Detailed knowledge of venom components present in Bothrops and Bothriechis venoms has the ability to improve treatment for those receiving an envenomation [12,15,84,126,130,140,141]. This is an important task as the two genera are responsible for a considerable number of bites [78,79,80,81,82,83] and knowledge of inter and intraspecific venom variation is crucial for optimizing neutralization capabilities of antivenoms cf. [84,126].
Further research into the identification and activity of toxins present in the venoms of snakes within this clade will be conducted in order to elucidate the selection pressures that have shaped these fine-tuned weapon systems. This research will hopefully contribute towards filling existing knowledge gaps concerning the evolution of venom composition and activity in the Viperidae.

4. Experimental Section

4.1. Venom Collection

A series of 8 snake species were chosen due to their phylogenetic positioning, genetic relatedness, body size/length variation, and lack of existing knowledge of their venom composition for proteomic analyses. Pooled venom samples were collected from both from adult wild and captive specimens residing in private collections. Venom samples were collected by inducing the snake to bite a sterile container and collected venom was immediately stored on dry ice or in liquid nitrogen during transport. Once in the laboratory, samples were stored at −80 °C after lyophilization. Species investigated were Bothriechis aurifer (captive snakes of unknown locality), Bothriechis lateralis (captive snakes of unknown locality), Bothriechis marchi (captive snakes of unknown locality), Bothriechis schlegelii (Costa Rica), Bothrops bilineata (French Guyana), Bothrops taeniata (captive snakes of unknown locality), Bothrops asper (Costa Rica), and Bothrops neuwiedi bolivianus (Bolivia).

4.2. Proteomics

4.2.1. 1D and 2D Gels Using SDS PAGE Electrophoresis Methodology

1D SDS-PAGE

Reduced 1D Tris-glycine SDS-PAGE methodology followed that of the study conducted by Ali et al., 2013 (cf. [11]). Specific changes were as follows: samples (20 ug) loaded were 15 µL total volume, gels were run at 90 V for 20 mins and then 120 V for 40 min, and gels were de-stained in 100 mL Milli Q.

Mini 2D SDS-PAGE

First Dimension

A sample of 125 µL solubilization buffer was added to a pre-washed sample of 300 µg of protein from required venom stock, prepared in a 15 mL Falcon tube. A solution of 0.75 µL of Bio-Lyte Ampholyte was added to the buffer and gently mixed. A sterile IPG strip (Bio-Rad Ready Strip Oakland, California non-linear pH 3–10, 7 cm, Hercules, CA, USA) was then added to the Corning Falcon Conical tube (Tamaulipas, Mexico) and allowed to absorb the solution overnight.
IPG strips were removed from Falcon tubes and placed gel side up an Isoelectric Focusing (IEF) machine (PROTEAN i12 IEF CELL Bio-RadLab Oakland, Hercules, CA, USA). Strips were covered with mineral oil and run overnight. Running conditions followed that of Ali et al., 2013 (cf. [11]).

Second Dimension

Required mini gels (7 cm × 6 cm) were cast using the same protocol for 1D Tris-glycine SDS-PAGE resolving layers (without the stacking layer). A layer of gel-overlay buffer (1 mL) was placed on top of the resolving layer once dry. A solution of 38 mg dithiothreitol (DTT) and 2.5 mL Equilibration buffer was added to each strip and gently mixed for 10 min. DTT solution was removed and a solution of 50 mg iodoacetamide (IAA) and 2.5 mL Equilibration buffer, gently mixed for 20 min.
Gel-overlay buffer was removed and 1 mL liquid agarose was laid over the gel. IPG strips were then embedded into the agarose layer, positive end on the left. A molecular marker gel bit was then inserted on the negative end. Gels were then placed in the electrophoresis container with a combination of 10× electrode buffer and Milli Q (80:800 mL). Power pack was set at 20 mA for 60 min and run at room temperature.
Once gels were finished they were removed from their glass cast and placed in 10 mL Coomassie brilliant blue stain, and set on a rocker overnight.

Liquid Chromatography-tandem Mass Spectrometry (LC-MS/MS) Analyses

In order to identify the types of toxins present in each gel, bands/spots from 1D SDS-PAGE glycine gels were picked, digested and analysed by LC-MS/MS. Visible protein bands were cut and stored at −80 °C. De-stain solution (500 µL 50 mM Ammonium bicarbonate (ABC)/50% Acetonitrile (ACN)) was added to each picked band/spot and left overnight. The destain was removed and repeated.

Reduction and Alkylation

Once blue stain had disappeared from picked bands, destain was removed. For 1D bands, a solution of 40 µL of 10 mM DTT was added to each band. Bands were then incubated at 60 °C for 30 min reducing disulphide bonds, and DTT removed. A solution of 40 µL of 55 mM IAA was added to the bands and incubated for 30 min at room temperature in the dark, alkylating free cysteine residues. Once IAA was removed, 100 µL of 50 mM ABC wash was added. Tubes were then spun down using a vortex for 2 min, ABC wash was removed and process repeated, washing away any excess IAA (removing IAA from the samples was important, as the solution can create adducts on proteins of interest). Remaining IAA was removed, and 100 µL of 100% ACN was added and left for minimum 5 min. ACN dehydrates the gel pieces, turning them white. Once gel pieces had turned white, ACN was removed.

Enzymatic Digestion

Gel pieces were rehydrated with 8 µL trypsin, brought up in hydrochloric acid (10 ng/µL, further diluted in 50mM ABC). This was left for 10–20 min at 4 °C. Depending on gel piece size, an additional 6 µL to 16 µL 50 mM ABC buffer was added. Gel pieces were kept moist overnight during digestion, incubated at 37 °C.

Peptide Extraction

Digested peptides were resuspended in 10 µL of 2.5% ACN/0.1% FA. Samples were filtered via ziptips using the following protocol; tips were wetted with 100% ACN (3 × 10 µL), and equilibrated with 10 µL 5% ACN/0.1% TFA, also repeated 3 times. Samples were loaded and pipetted up and down at least twice the volume of sample. Tip was washed with 5% ACN/0.1% TFA (3 × 10 µL). Samples were eluted with 80% ACN/0.1% TFA (10 µL) in a new tube. Combined samples were transferred into glass vials and diluted with 0.1% FA (final loading volume 20 µL).
For LC-MS/MS analysis, parameters are as follows; samples were separated using RP-chromatography on a Dionex Ultimate 3000 RSLC nano-system (Lifetech, Carlsbad, CA, USA). Samples were desalted on a Thermo PepMap 100 C18 trap (Lifetech, Carlsbad, CA, USA) (0.3 × 5 mm, 5 µm) for 5 min with a flow rate of 30 µL/min. This was followed by separation on an Acclaim PepMap RSLC C18 (Lifetech, Carlsbad, CA, USA) (150 mm × 75 µm) column at a flow rate of 300 nL/min. A gradient of 10%–70% buffer B over 7 min where buffer A (1% ACN/0.1% FA) and buffer B (80% ACN/0.1% FA) was used to separate peptides. Eluted peptides were directly analysed on an Orbitap Elite mass spectrometer (Thermo Scientific, Carlsbad, CA, USA) using an NSI electrospray interface. Source parameters included a capillary temperature of 275 °C; S-Lens RF level at 60%; source voltage of 2 kV and maximum injections times of 200 ms for MS and 150 ms for MS2. Instrument parameters included an FTMS scan across m/z range 350–1800 at 60,000 resolution followed by information dependent acquisition of the top 10 peptides across m/z 40-1800. Dynamic ion exclusion was employed using a 15 s interval. Charge state screening was enabed with rejection of +1 charged ions and monoisotopic precursor selection enabled. Data was converted to mascot generic format (mgf) using the msConvert software (ProteoWizard v2.0) and searched using Protein Pilot™ v5.0 (Sciex).

4.2.2. Shotgun Analysis

Reduction, Alkylation and Trypsinization

Pre washed 10 µg protein venom samples, solubilised with 40 µL of Milli Q, was added to 5 µL of 1 M Ammonium Carbonate. A total of 55 µL per sample (97.5% acetonitrile, 2% iododethanol and 0.5% ttriethylphosphene) was added, with a total combined volume of 100 µL. Samples were then incubated at 37 °C for 2 h, and then freeze-dried. Samples were re-suspended in 25 µL of 40 mM ammonium bicarbonate. Trypsin (1 ng) was added to each re-suspended sample and incubated at 37 °C overnight. Samples were dried and re-suspended in 40 µL 5% acetonitrile and 1% formic acid in preparation for LC-MS-MS analysis according to methods described in detail elsewhere [11].

4.3. Bioactivity Testing

4.3.1. Fluorescent Determination of Metalloprotease Activity

Freeze dried venom was reconstituted in a buffer containing 150 mM NaCl and 50 mM Tri-HCl (pH 6.3). Metalloprotease activity in venom (10 ng/µL) was measured by adding quenched fluorescent substrate (10 µM final; Fluorogenic Peptide Substrate, R & D systems, Cat#ES001, Minneapolis, Minnesota). Fluorescence was monitored (excitation at 320 nm and emission at 405 nm) over 100 min, and rate of substrate cleavage calculated from a standard curve of known fluorophore (7-methoxycoumarin) concentrations. The specificity of metalloprotease activity was confirmed by incubating venom for 24 h at 37 °C with EDTA (50 mM).

4.3.2. Zymogram PAGE 10% and 12% Gels Bio-Rad

Dried samples of 30 µg were prepared for each species and reconstituted in 15 µL of Bio-Rad Sample Buffer (#161-0764 supplied). Bio-Rad Zymogram PAGE ready-made 10% gelatin and 12% casein gels (1,611,113 and 1,611,114 respectively) were run using standard PAGE protocol at 100V for 90 min in 1× Tris-Glycine SDS Running Buffer (working concentration 25 mM Tris, 192 mM glycine, 0.1% SDS). Gel was removed from case and incubated at room temperature for 30 min in 100 mL 1× Zymogram Renaturing Buffer (#161-0765, 2.5% Triton ×100, supplied). Renaturing Buffer was replaced with 100 mL 1x zymogram Development Buffer (#161-0766, 50 mM Tris, 200 mM NaCl, 5 mM CaCl2, 0.02% Brij-35 pH 7.5, supplied) and equilibrated for 30 min at room temperature. Development buffer was replaced with renaturing buffer and incubated overnight at 37 °C for maximum sensitivity.
Zymogram gels test for proteolytic activity when performing protein characterization. The varying gel compositions act as differing substrates for proteases separated in the gel. Gels were stained in R-250 Coomassie Brilliant Blue R-250 for 1 h and stored in 100 mL Milli Q water. Gelatinases, Matrix Metalloprotease (MMP) 2 and MMP 9 can be detected via 10% Zymography, while MMP 1, MMP 7, MMP 12, and MMP 13 are detected via 12% Zymography. By performing both Zymography tests with the identical venom concentrations under identical conditions, direct comparison between venoms can be made.

4.3.3. sPLA2 Assay Kit Cayman Activity

Triplicate venom concentrations of 1 µg were prepared, along with the of positive control bee venom supplied with the Cayman 765001 assay kit. Triplicates were used instead of duplicates to strengthen the statistical power. Supplied methods consisted of non-enzymatic controls, positive controls and sample wells. Non-enzymatic controls were prepared using 10 µL DTNB (10 mM DTNB in 0.4 M Tris-HCl, pH 8.0) (supplied) and 15 µL Assay Buffer (25 mM Tris-HCl, pH 7.5, 10 mM CaCl2, 100 mM KCl, 0.3 mM Triton X-100). Positive control wells were prepared using 10 uL DTNB, 10 uL provided bee venom (100 µg/mL) and 5 µL Assay buffer. Sample wells consisted of adding 10 µL DTNB, 10 µL sample and 5 µL Assay Buffer. Reactions were initiated by adding 200 µL Substrate Solution (supplied) to all wells. Samples were prepared in a 96 well plate and run at A-405 nm every minute for 10 min using Fluostar Optima absorbance (BMG Labtech, Ortenberg Germany). Data was analysed using MARS data analysis software (BMG Labtech, Ortenberg Germany) and results were published using PRISM software graphs.

Supplementary Materials

The following are available online at www.mdpi.com/2072-6651/8/7/210/s1.

Acknowledgments

BGF was funded by an Australian Research Council Future Fellowship and by the University of Queensland. SAA was the recipient of postdoctoral fellowship (PDRF Phase II Batch-V) from Higher Education Commission (HEC Islamabad) Pakistan. TNWJ was funded by an Australian Postgraduate Award. MJ and IK acknowledge funding from the University of Queensland (International Postgraduate Research Scholarship. CC was supported by the National Science Foundation Graduate Re- search Fellowship under Grant No. 2012134810 and therefore must include the statement “Any opinion, findings, and con- clusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation.”

Author Contributions

Jordan Debono—experimental work, data analysis, manuscript writing; Chip Cochran—data analysis, manuscript writing; Sanjaya Kuruppu—experimental work, data analysis; Amanda Nouwens—experimental work, data analysis; Niwanthi W. Rajapakse—experimental work, data analysis; Minami Kawaski—experimental work, data analysis; Kelly Wood—experimental work; James Dobson—experimental work; Kate Baumann—experimental work; Mahdokht Jouiaei—experimental work; Timothy N.W. Jackson—experimental work and manuscript preparation; Ivan Koludarov—experimental work; Dolyce Low—experimental work, data analysis, manuscript writing; Syed A. Ali—experimental work, data analysis, manuscript writing; Andrew Barnes—experimental work, data analysis; A. Ian Smith—experimental work, data analysis; Bryan G. Fry—study concept and design, experimental work, data analysis, manuscript writing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Anaya, M.; Rael, E.D.; Lieb, C.S.; Perez, J.C.; Salo, R.J. Antibody Detection of Venom Protein Variation within a Population of the Rattlesnake Crotalus v. viridis. J. Herpetol. 1992, 26, 473–482. [Google Scholar] [CrossRef]
  2. Bush, S.P.; Green, S.M.; Moynihan, J.A.; Hayes, W.K.; Cardwell, M.D. Crotalidae polyvalent immune Fab (ovine) antivenom is efficacious for envenomations by Southern Pacific rattlesnakes (Crotalus helleri). Ann. Emerg. Med. 2002, 40, 619–624. [Google Scholar] [CrossRef] [PubMed]
  3. Casewell, N.R.; Wüster, W.; Vonk, F.J.; Harrison, R.A.; Fry, B.G. Complex cocktails: The evolutionary novelty of venoms. Trends Ecol. Evol. 2012, 28, 219–229. [Google Scholar] [CrossRef] [PubMed]
  4. Chippaux, J.-P.; Williams, V.; White, J. Snake venom variability: Methods of study, results and interpretation. Toxicon 1991, 29, 1279–1303. [Google Scholar] [CrossRef]
  5. Fry, B.G.; Casewell, N.R.; Wuster, W.; Vidal, N.; Young, B.; Jackson, T.N. The structural and functional diversification of the Toxicofera reptile venom system. Toxicon 2012, 60, 434–448. [Google Scholar] [CrossRef] [PubMed]
  6. Mackessy, S.P.; Baxter, L.M. Bioweapons synthesis and storage: The venom gland of front-fanged snakes. Zool. Anz. J. Comp. Zool. 2006, 245, 147–159. [Google Scholar] [CrossRef]
  7. Angulo, Y.; Escolano, J.; Lomonte, B.; Gutierrez, J.M.; Sanz, L.; Calvete, J.J. Snake venomics of Central American pitvipers: Clues for rationalizing the distinct envenomation profiles of Atropoides nummifer and Atropoides picadoi. J. Proteome Res. 2008, 7, 708–719. [Google Scholar] [CrossRef] [PubMed]
  8. Fry, B.G. From genome to “venome”: Molecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins. Genome Res. 2005, 15, 403–420. [Google Scholar] [CrossRef] [PubMed]
  9. Fry, B.G.; Undheim, E.A.; Ali, S.A.; Jackson, T.N.; Debono, J.; Scheib, H.; Ruder, T.; Morgenstern, D.; Cadwallader, L.; Whitehead, D.; et al. Squeezers and Leaf-cutters: Differential Diversification and Degeneration of the Venom System in Toxicoferan Reptiles. Mol. Cell Proteom. 2013, 12, 1881–1899. [Google Scholar] [CrossRef] [PubMed]
  10. Alape-Girón, A.; Sanz, L.; Escolano, J.; Flores-Díaz, M.; Madrigal, M.; Sasa, M.; Calvete, J.J. Snake venomics of the lancehead pitviper Bothrops asper: Geographic, individual, and ontogenetic variations. J. Proteome Res. 2008, 7, 3556–3571. [Google Scholar] [CrossRef] [PubMed]
  11. Ali, S.A.; Yang, D.; Jackson, T.N.; Undheim, E.A.; Koludarov, I.; Wood, K.; Jones, A.; Hodgson, W.C.; McCarthy, S.; Ruder, T. Venom proteomic characterization and relative antivenom neutralization of two medically important Pakistani elapid snakes (Bungarus sindanus and Naja naja). J. Proteom. 2013, 89, 15–23. [Google Scholar] [CrossRef] [PubMed]
  12. Boldrini-Franca, J.; Correa-Netto, C.; Silva, M.M.; Rodrigues, R.S.; De La Torre, P.; Perez, A.; Soares, A.M.; Zingali, R.B.; Nogueira, R.A.; Rodrigues, V.M.; et al. Snake venomics and antivenomics of Crotalus durissus subspecies from Brazil: Assessment of geographic variation and its implication on snakebite management. J. Proteom. 2010, 73, 1758–1776. [Google Scholar] [CrossRef] [PubMed]
  13. Calvete, J.J.; Escolano, J.; Sanz, L. Snake venomics of Bitis species reveals large intragenus venom toxin composition variation: Application to taxonomy of congeneric taxa. J. Proteome Res. 2007, 6, 2732–2745. [Google Scholar] [CrossRef] [PubMed]
  14. Calvete, J.J.; Pérez, A.; Lomonte, B.; Sánchez, E.E.; Sanz, L. Snake venomics of Crotalus tigris: The minimalist toxin arsenal of the deadliest neartic rattlesnake venom. Evolutionary clues for generating a pan-specific antivenom against crotalid type II venoms. J. Proteome Res. 2012, 11, 1382–1390. [Google Scholar] [CrossRef] [PubMed]
  15. Calvete, J.J.; Sanz, L.; Cid, P.; de la Torre, P.; Flores-Díaz, M.; Dos Santos, M.C.; Borges, A.; Bremo, A.; Angulo, Y.; Lomonte, B. Snake venomics of the Central American rattlesnake Crotalus simus and the South American Crotalus durissus complex points to neurotoxicity as an adaptive paedomorphic trend along Crotalus dispersal in South America. J. Proteome Res. 2009, 9, 528–544. [Google Scholar] [CrossRef] [PubMed]
  16. Fernández, J.; Alape-Girón, A.; Angulo, Y.; Sanz, L.; Gutiérrez, J.M.; Calvete, J.J.; Lomonte, B. Venomic and antivenomic analyses of the Central American coral snake, Micrurus nigrocinctus (Elapidae). J. Proteome Res. 2011, 10, 1816–1827. [Google Scholar] [CrossRef] [PubMed]
  17. Fry, B.G.; Lumsden, N.G.; Wüster, W.; Wickramaratna, J.C.; Hodgson, W.C.; Kini, R.M. Isolation of a neurotoxin (α-colubritoxin) from a nonvenomous colubrid: Evidence for early origin of venom in snakes. J. Mol. Evol. 2003, 57, 446–452. [Google Scholar] [CrossRef] [PubMed]
  18. Fry, B.G.; Wickramaratna, J.C.; Hodgson, W.C.; Alewood, P.F.; Kini, R.; Ho, H.; Wüster, W. Electrospray liquid chromatography/mass spectrometry fingerprinting of Acanthophis (death adder) venoms: Taxonomic and toxinological implications. Rapid commun. Mass Spectrum. 2002, 16, 600–608. [Google Scholar] [CrossRef] [PubMed]
  19. Fry, B.G.; Wüster, W.; Ramjan, R.; Fadil, S.; Jackson, T.; Martelli, P.; Kini, R.M. Analysis of Colubroidea snake venoms by liquid chromatography with mass spectrometry: Evolutionary and toxinological implications. Rapid Commun. Mass Spectrom. 2003, 17, 2047–2062. [Google Scholar] [CrossRef] [PubMed]
  20. Aird, S.D.; Watanabe, Y.; Villar-Briones, A.; Roy, M.C.; Terada, K.; Mikheyev, A.S. Quantitative high-throughput profiling of snake venom gland transcriptomes and proteomes (Ovophis okinavensis and Protobothrops flavoviridis). BMC Genom. 2013, 14, 790. [Google Scholar] [CrossRef] [PubMed]
  21. Fry, B.G.; Roelants, K.; Winter, K.; Hodgson, W.C.; Griesman, L.; Kwok, H.F.; Scanlon, D.; Karas, J.; Shaw, C.; Wong, L. Novel venom proteins produced by differential domain-expression strategies in beaded lizards and gila monsters (genus Heloderma). Mol. Biol. Evol. 2010, 27, 395–407. [Google Scholar] [CrossRef] [PubMed]
  22. Fry, B.G.; Scheib, H.; de Azevedo, I.D.L.J.; Silva, D.A.; Casewell, N.R. Novel transcripts in the maxillary venom glands of advanced snakes. Toxicon 2012, 59, 696–708. [Google Scholar] [CrossRef] [PubMed]
  23. Fry, B.G.; Wroe, S.; Teeuwisse, W.; van Osch, M.J.; Moreno, K.; Ingle, J.; McHenry, C.; Ferrara, T.; Clausen, P.; Scheib, H. A central role for venom in predation by Varanus komodoensis (Komodo Dragon) and the extinct giant Varanus (Megalania) priscus. Proc. Natl. Acad. Sci. USA 2009, 106, 8969–8974. [Google Scholar] [CrossRef] [PubMed]
  24. Margres, M.J.; Aronow, K.; Loyacano, J.; Rokyta, D.R. The venom-gland transcriptome of the eastern coral snake (Micrurus fulvius) reveals high venom complexity in the intragenomic evolution of venoms. BMC Genom. 2013, 14, 531. [Google Scholar] [CrossRef] [PubMed]
  25. Rokyta, D.R.; Wray, K.P.; Lemmon, A.R.; Lemmon, E.M.; Caudle, S.B. A high-throughput venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) and evidence for pervasive positive selection across toxin classes. Toxicon 2011, 57, 657–671. [Google Scholar] [CrossRef] [PubMed]
  26. Rokyta, D.R.; Wray, K.P.; Margres, M.J. The genesis of an exceptionally lethal venom in the timber rattlesnake (Crotalus horridus) revealed through comparative venom-gland transcriptomics. BMC Genom. 2013, 14, 394. [Google Scholar] [CrossRef] [PubMed]
  27. Terrat, Y.; Sunagar, K.; Fry, B.G.; Jackson, T.N.; Scheib, H.; Fourmy, R.; Verdenaud, M.; Blanchet, G.; Antunes, A.; Ducancel, F. Atractaspis aterrima toxins: The first insight into the molecular evolution of venom in side-stabbers. Toxins 2013, 5, 1948–1964. [Google Scholar] [CrossRef] [PubMed]
  28. Wagstaff, S.C.; Sanz, L.; Juárez, P.; Harrison, R.A.; Calvete, J.J. Combined snake venomics and venom gland transcriptomic analysis of the ocellated carpet viper, Echis ocellatus. J. Proteom. 2009, 71, 609–623. [Google Scholar] [CrossRef] [PubMed]
  29. Fry, B.G.; Scheib, H.; van der Weerd, L.; Young, B.; McNaughtan, J.; Ramjan, S.R.; Vidal, N.; Poelmann, R.E.; Norman, J.A. Evolution of an arsenal structural and functional diversification of the venom system in the advanced snakes (Caenophidia). Mol. Cell. Proteom. 2008, 7, 215–246. [Google Scholar] [CrossRef] [PubMed]
  30. Gutiérrez, J.M.; Sanz, L.; Escolano, J.; Fernández, J.; Lomonte, B.; Angulo, Y.; Rucavado, A.; Warrell, D.A.; Calvete, J.J. Snake Venomics of the Lesser Antillean Pit Vipers Bothrops caribbaeus and Bothrops lanceolatus: Correlation with Toxicological Activities and Immunoreactivity of a Heterologous Antivenom†. J. Proteome Res. 2008, 7, 4396–4408. [Google Scholar] [CrossRef] [PubMed]
  31. Lomonte, B.; Escolano, J.; Fernandez, J.; Sanz, L.; Angulo, Y.; Gutierrez, J.M.; Calvete, J.J. Snake venomics and antivenomics of the arboreal neotropical pitvipers Bothriechis lateralis and Bothriechis schlegelii. J. Proteome Res. 2008, 7, 2445–2457. [Google Scholar] [CrossRef] [PubMed]
  32. Sanz, L.; Gibbs, H.L.; Mackessy, S.P.; Calvete, J.J. Venom proteomes of closely related Sistrurus rattlesnakes with divergent diets. J. Proteome Res. 2006, 5, 2098–2112. [Google Scholar] [CrossRef] [PubMed]
  33. Tashima, A.K.; Sanz, L.; Camargo, A.C.; Serrano, S.M.; Calvete, J.J. Snake venomics of the Brazilian pitvipers Bothrops cotiara and Bothrops fonsecai. Identification of taxonomy markers. J. Proteom. 2008, 71, 473–485. [Google Scholar] [CrossRef] [PubMed]
  34. van der Weyden, L.; Hains, P.G.; Broady, K.W. Characterisation of the biochemical and biological variations from the venom of the death adder species (Acanthophis antarcticus, A. praelongus and A. pyrrhus). Toxicon 2000, 38, 1703–1713. [Google Scholar] [CrossRef]
  35. Castro, E.N.; Lomonte, B.; del Carmen Gutiérrez, M.; Alagón, A.; Gutiérrez, J.M. Intraspecies variation in the venom of the rattlesnake Crotalus simus from Mexico: Different expression of crotoxin results in highly variable toxicity in the venoms of three subspecies. J. Proteom. 2013, 87, 103–121. [Google Scholar] [CrossRef] [PubMed]
  36. Daltry, J.C.; Ponnudurai, G.; Shin, C.K.; Tan, N.-H.; Thorpe, R.S.; Wolfgang, W. Electrophoretic profiles and biological activities: Intraspecific variation in the venom of the Malayan pit viper (Calloselasma rhodostoma). Toxicon 1996, 34, 67–79. [Google Scholar] [CrossRef]
  37. Daltry, J.C.; Wüster, W.; Thorpe, R.S. Diet and snake venom evolution. Nature 1996, 379, 537–540. [Google Scholar] [CrossRef] [PubMed]
  38. Forstner, M.; Hilsenbeck, R.; Scudday, J. Geographic variation in whole venom profiles from the mottled rock rattlesnake (Crotalus lepidus lepidus) in Texas. J. Herpetol. 1997, 31, 277–287. [Google Scholar] [CrossRef]
  39. French, W.J.; Hayes, W.K.; Bush, S.P.; Cardwell, M.D.; Bader, J.O.; Rael, E.D. Mojave toxin in venom of Crotalus helleri (Southern Pacific Rattlesnake): Molecular and geographic characterization. Toxicon 2004, 44, 781–791. [Google Scholar] [CrossRef] [PubMed]
  40. Mackessy, S.P. Evolutionary trends in venom composition in the western rattlesnakes (Crotalus viridis sensu lato): Toxicity vs. tenderizers. Toxicon 2010, 55, 1463–1474. [Google Scholar] [CrossRef] [PubMed]
  41. Salazar, A.M.; Guerrero, B.; Cantu, B.; Cantu, E.; Rodríguez-Acosta, A.; Pérez, J.C.; Galán, J.A.; Tao, A.; Sánchez, E.E. Venom variation in hemostasis of the southern Pacific rattlesnake (Crotalus oreganus helleri): Isolation of hellerase. Comp. Biochem. Physiol. C Toxicol. Pharmacol. 2009, 149, 307–316. [Google Scholar] [CrossRef] [PubMed]
  42. Sunagar, K.; Undheim, E.A.; Scheib, H.; Gren, E.C.; Cochran, C.; Person, C.E.; Koludarov, I.; Kelln, W.; Hayes, W.K.; King, G.F. Intraspecific venom variation in the medically significant Southern Pacific Rattlesnake (Crotalus oreganus helleri): Biodiscovery, clinical and evolutionary implications. J. Proteom. 2014, 99, 68–83. [Google Scholar] [CrossRef] [PubMed]
  43. Wilkinson, J.A.; Glenn, J.L.; Straight, R.C.; Sites, J.W., Jr. Distribution and Genetic Variation in Venom A and B Populations of the Mojave Rattlesnake (Crotalus scutulatus scutulatus in Arizona. Herpetologica 1991, 47, 54–68. [Google Scholar]
  44. Menezes, M.C.; Furtado, M.F.; Travaglia-Cardoso, S.R.; Camargo, A.C.; Serrano, S.M. Sex-based individual variation of snake venom proteome among eighteen Bothrops jararaca siblings. Toxicon 2006, 47, 304–312. [Google Scholar] [CrossRef] [PubMed]
  45. Pimenta, D.C.; Prezoto, B.C.; Konno, K.; Melo, R.L.; Furtado, M.F.; Camargo, A.; Serrano, S.M. Mass spectrometric analysis of the individual variability of Bothrops jararaca venom peptide fraction. Evidence for sex-based variation among the bradykinin-potentiating peptides. Rapid Commun. Mass Spectrum. 2007, 21, 1034–1042. [Google Scholar] [CrossRef] [PubMed]
  46. Durban, J.; Pérez, A.; Sanz, L.; Gómez, A.; Bonilla, F.; Chacón, D.; Sasa, M.; Angulo, Y.; Gutiérrez, J.M.; Calvete, J.J. Integrated “omics” profiling indicates that miRNAs are modulators of the ontogenetic venom composition shift in the Central American rattlesnake, Crotalus simus simus. BMC Genom. 2013, 14, 234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Guércio, R.A.; Shevchenko, A.; Shevchenko, A.; López-Lozano, J.L.; Paba, J.; Sousa, M.V.; Ricart, C.A. Ontogenetic variations in the venom proteome of the Amazonian snake Bothrops atrox. Proteome Sci. 2006, 4, 11. [Google Scholar] [CrossRef] [PubMed]
  48. López-Lozano, J.L.; de Sousa, M.V.; Ricart, C.A.O.; Chávez-Olortegui, C.; Sanchez, E.F.; Muniz, E.G.; Bührnheim, P.F.; Morhy, L. Ontogenetic variation of metalloproteinases and plasma coagulant activity in venoms of wild Bothrops atrox specimens from Amazonian rain forest. Toxicon 2002, 40, 997–1006. [Google Scholar] [CrossRef]
  49. Mackessy, S.P. Venom ontogeny in the Pacific rattlesnakes Crotalus viridis helleri and C. v. oreganus. Copeia 1988, 1988, 92–101. [Google Scholar] [CrossRef]
  50. Mackessy, S.P.; Williams, K.; Ashton, K.G. Ontogenetic variation in venom composition and diet of Crotalus oreganus concolor: A case of venom paedomorphosis? J. Inf. 2003, 2003, 769–782. [Google Scholar] [CrossRef]
  51. Saldarriaga, M.M.A.; Otero, R.; Núñez, V.; Toro, M.F.; Dı́az, A.; Gutiérrez, J.M.A. Ontogenetic variability of Bothrops atrox and Bothrops asper snake venoms from Colombia. Toxicon 2003, 42, 405–411. [Google Scholar] [CrossRef]
  52. Zelanis, A.; Tashima, A.K.; Rocha, M.M.; Furtado, M.F.; Camargo, A.C.; Ho, P.L.; Serrano, S.M. Analysis of the ontogenetic variation in the venom proteome/peptidome of Bothrops jararaca reveals different strategies to deal with prey. J. Proteome Res. 2010, 9, 2278–2291. [Google Scholar] [CrossRef] [PubMed]
  53. Mebs, D. Toxicity in animals. Trends in evolution? Toxicon 2001, 39, 87–96. [Google Scholar] [CrossRef]
  54. Sasa, M. Diet and snake venom evolution: Can local selection alone explain intraspecific venom variation? Toxicon 1999, 37, 249–252. [Google Scholar] [PubMed]
  55. Aird, S.D. Ophidian envenomation strategies and the role of purines. Toxicon 2002, 40, 335–393. [Google Scholar] [CrossRef]
  56. Barlow, A.; Pook, C.E.; Harrison, R.A.; Wüster, W. Coevolution of diet and prey-specific venom activity supports the role of selection in snake venom evolution. Proc. R. Soc. B Biol. Sci. 2009, 276, 2443–2449. [Google Scholar] [CrossRef] [PubMed]
  57. Brust, A.; Sunagar, K.; Undheim, E.A.; Vetter, I.; Yang, D.C.; Casewell, N.R.; Jackson, T.N.; Koludarov, I.; Alewood, P.F.; Hodgson, W.C. Differential evolution and neofunctionalization of snake venom metalloprotease domains. Mol. Cell. Proteom. 2013, 12, 651–663. [Google Scholar] [CrossRef] [PubMed]
  58. Gibbs, H.L.; Mackessy, S.P. Functional basis of a molecular adaptation: Prey-specific toxic effects of venom from Sistrurus rattlesnakes. Toxicon 2009, 53, 672–679. [Google Scholar] [CrossRef] [PubMed]
  59. Pawlak, J.; Mackessy, S.P.; Fry, B.G.; Bhatia, M.; Mourier, G.; Fruchart-Gaillard, C.; Servent, D.; Ménez, R.; Stura, E.; Ménez, A. Denmotoxin, a three-finger toxin from the colubrid snake Boiga dendrophila (Mangrove Catsnake) with bird-specific activity. J. Biol. Chem. 2006, 281, 29030–29041. [Google Scholar] [CrossRef] [PubMed]
  60. Sunagar, K.; Johnson, W.E.; O’Brien, S.J.; Vasconcelos, V.; Antunes, A. Evolution of CRISPs associated with toxicoferan-reptilian venom and mammalian reproduction. Mol. Biol. Evol. 2012, 29, 1807–1822. [Google Scholar] [CrossRef] [PubMed]
  61. Parkinson, C.L. Molecular systematics and biogeographical history of pitvipers as determined by mitochondrial ribosomal DNA sequences. Copeia 1999, 1999, 576–586. [Google Scholar] [CrossRef]
  62. Parkinson, C.L.; Campbell, J.A.; Chippindale, P.T.; Schuett, G. Multigene phylogenetic analysis of pitvipers, with comments on their biogeography. In Biology of the Vipers; Schuett, G.W., Höggren, M., Eds.; Eagle Mountain Publishing: Eagle Mountain, UT, USA, 2002; pp. 93–110. [Google Scholar]
  63. Kraus, F.; Mink, D.G.; Brown, W.M. Crotaline intergeneric relationships based on mitochondrial DNA sequence data. Copeia 1996, 1996, 763–773. [Google Scholar] [CrossRef]
  64. Wüster, W.; Peppin, L.; Pook, C.E.; Walker, D.E. A nesting of vipers: Phylogeny and historical biogeography of the Viperidae (Squamata: Serpentes). Mol. Phylogenetics Evol. 2008, 49, 445–459. [Google Scholar] [CrossRef] [PubMed]
  65. Castoe, T.A.; Parkinson, C.L. Bayesian mixed models and the phylogeny of pitvipers (Viperidae: Serpentes). Mol. Phylogenetics Evol. 2006, 39, 91–110. [Google Scholar] [CrossRef] [PubMed]
  66. Vidal, N.; Lecointre, G.; Vié, J.C.; Gasc, J.P. Molecular systematics of pitvipers: Paraphyly of the Bothrops complex. Comptes Rendus de l'Académie des Sciences-Series III-Sciences de la Vie 1997, 320, 95–101. [Google Scholar] [CrossRef]
  67. Castoe, T.A.; Daza, J.M.; Smith, E.N.; Sasa, M.M.; Kuch, U.; Campbell, J.A.; Chippindale, P.T.; Parkinson, C.L. Comparative phylogeography of pitvipers suggests a consensus of ancient Middle American highland biogeography. J. Biogeogr. 2009, 36, 88–103. [Google Scholar] [CrossRef]
  68. Cadle, J.E. Phylogenetic relationships among vipers: Immunological evidence. In Biology of the Pitvipers; Selva Natural History Book Publishers: Tyler, TX, USA, 1992; pp. 41–48. [Google Scholar]
  69. Salomão, M.; Wüster, W.; Thorpe, R. BBBSP MtDNA phylogeny of Neotropical pitvipers of the genus Bothrops (Squamata: Serpentes: Viperidae). Kaupia 1999, 8, 127–134. [Google Scholar]
  70. Salomão, M.; Wüster, W.; Thorpe, R.; Touzet, J. DNA evolution of South American pitvipers of the genus Bothrops (Reptilia: Serpentes: Viperidae). Venom. Snakes Ecol. Evol. Snakebite 1997, 70, 89–98. [Google Scholar]
  71. Werman, S.D. Phylogenetic Relationships of Central and South American Pitvipers of the Genus Bothrops (Sensu Lato): Cladistic Analyses of Biochemical and Anatomical Characters. In Biology Of the Pitvipers; Campbell, J.A., Brodie, E.D., Jr., Eds.; Selva: Tyler, TX, USA, 1992; pp. 21–40. [Google Scholar]
  72. Wüster, W.; Salomão, M.; Quijada-Mascareñas, J.A.; Thorpe, R.; Duckett, G.; Puorto, M.; Theakston, R.; Warrell, D. Origins and evolution of the South American pitviper fauna: Evidence from mitochondrial DNA sequence analysis. In Biology of the Vipers; Schuett, G.W., Höggren, M., Eds.; Eagle Mountain Publishing: Eagle Mountain, UT, USA, 2002; pp. 111–128. [Google Scholar]
  73. Schatti, B.; Kramer, E.; Touzet, J. Systematic remarks on a rare Crotalid snake from Ecuador, Bothriechis albocarinata (Shreve), with some comments on generic arrangement of Neotropical pitvipers. Ibid 1990, 97, 877–885. [Google Scholar] [CrossRef]
  74. Burger, W.L. Genera of Pitvipers (Serpentes: Crotalidae). PhD Thesis, University of Kansas, Lawrence, KS, USA, 1971. [Google Scholar]
  75. Campbell, J.A.; Lamar, W.W. Taxonomic status of miscellaneous Neotropical viperids, with the descriptions of a new genus. Occ. Pap. Mus. Tex. Tech. Univ. 1992, 153, 1–31. [Google Scholar]
  76. Fenwick, A.M.; Gutberlet, R.L., Jr.; Evans, J.A.; Parkinson, C.L. Morphological and molecular evidence for phylogeny and classification of South American pitvipers, genera Bothrops, Bothriopsis, and Bothrocophias (Serpentes: Viperidae). Zool. J. Linn. Soc. 2009, 156, 617–640. [Google Scholar] [CrossRef]
  77. Carrasco, P.A.; Mattoni, C.I.; Leynaud, G.C.; Scrocchi, G.J. Morphology, phylogeny and taxonomy of South American bothropoid pitvipers (Serpentes, Viperidae). Zool. Scr. 2012, 41, 109–124. [Google Scholar] [CrossRef]
  78. Otero, R.; Gutiérrez, J.; Mesa, M.A.B.; Duque, E.; Rodrı́guez, O.; Arango, J.L.; Gómez, F.; Toro, A.; Cano, F.; Rodrı́guez, L.M.A. Complications of Bothrops, Porthidium, and Bothriechis snakebites in Colombia. A clinical and epidemiological study of 39 cases attended in a university hospital. Toxicon 2002, 40, 1107–1114. [Google Scholar] [CrossRef]
  79. Rosenfeld, G. Symptomatology, pathology and treatment of snake bites in South America. Venom. Anim. Their Venoms 1971, 2, 345–384. [Google Scholar]
  80. Torrez, P.; Said, R.; Quiroga, M.; Duarte, M.; França, F. Forest pit viper (Bothriopsis bilineata bilineata) bite in the Brazilian Amazon with acute kidney injury and persistent thrombocytopenia. Toxicon 2014, 85, 27–30. [Google Scholar] [CrossRef] [PubMed]
  81. Warrell, D.A. Snake bite. Lancet 2010, 375, 77–88. [Google Scholar] [CrossRef]
  82. Calvete, J.J.; Sanz, L.; Angulo, Y.; Lomonte, B.; Gutierrez, J.M. Venoms, venomics, antivenomics. FEBS Lett. 2009, 583, 1736–1743. [Google Scholar] [CrossRef] [PubMed]
  83. Gibbs, H.L.; Sanz, L.; Chiucchi, J.E.; Farrell, T.M.; Calvete, J.J. Proteomic analysis of ontogenetic and diet-related changes in venom composition of juvenile and adult Dusky Pigmy rattlesnakes (Sistrurus miliarius barbouri). J. Proteom. 2011, 74, 2169–2179. [Google Scholar] [CrossRef] [PubMed]
  84. Juárez, P.; Sanz, L.; Calvete, J.J. Snake venomics: Characterization of protein families in Sistrurus barbouri venom by cysteine mapping, N-terminal sequencing, and tandem mass spectrometry analysis. Proteomics 2004, 4, 327–338. [Google Scholar] [CrossRef] [PubMed]
  85. Margres, M.J.; McGivern, J.J.; Seavy, M.; Wray, K.P.; Facente, J.; Rokyta, D.R. Comparative genotype-phenotype mapping reveals distinct modes of venom expression evolution in the sympatric eastern diamondback rattlesnake (Crotalus adamanteus) and eastern coral snake (Micrurus fulvius). Genetics 2014, 199, 165. [Google Scholar] [CrossRef] [PubMed]
  86. Margres, M.J.; McGivern, J.J.; Wray, K.P.; Seavy, M.; Calvin, K.; Rokyta, D.R. Linking the transcriptome and proteome to characterize the venom of the eastern diamondback rattlesnake (Crotalus adamanteus). J. Proteom. 2014, 96, 145–158. [Google Scholar] [CrossRef] [PubMed]
  87. Martínez-Romero, G.; Rucavado, A.; Lazcano, D.; Gutiérrez, J.M.; Borja, M.; Lomonte, B.; Garza-García, Y.; Zugasti-Cruz, A. Comparison of venom composition and biological activities of the subspecies Crotalus lepidus lepidus, Crotalus lepidus klauberi and Crotalus lepidus morulus from Mexico. Toxicon 2013, 71, 84–95. [Google Scholar] [CrossRef] [PubMed]
  88. Massey, D.J.; Calvete, J.J.; Sánchez, E.E.; Sanz, L.; Richards, K.; Curtis, R.; Boesen, K. Venom variability and envenoming severity outcomes of the Crotalus scutulatus scutulatus (Mojave rattlesnake) from Southern Arizona. J. Proteom. 2012, 75, 2576–2587. [Google Scholar] [CrossRef] [PubMed]
  89. Pahari, S.; Mackessy, S.P.; Kini, R.M. The venom gland transcriptome of the Desert Massasauga Rattlesnake (Sistrurus catenatus edwardsii): Towards an understanding of venom composition among advanced snakes (Superfamily Colubroidea). BMC Mol. Biol. 2007, 8, 115. [Google Scholar] [CrossRef] [PubMed]
  90. Rokyta, D.R.; Wray, K.P.; McGivern, J.J.; Margres, M.J. The transcriptomic and proteomic basis for the evolution of a novel venom phenotype within the Timber Rattlesnake (Crotalus horridus). Toxicon 2015, 98, 34–48. [Google Scholar] [CrossRef] [PubMed]
  91. Correa-Netto, C.; Teixeira-Araujo, R.; Aguiar, A.S.; Melgarejo, A.R.; De-Simone, S.G.; Soares, M.R.; Foguel, D.; Zingali, R.B. Immunome and venome of Bothrops jararacussu: A proteomic approach to study the molecular immunology of snake toxins. Toxicon 2010, 55, 1222–1235. [Google Scholar] [CrossRef] [PubMed]
  92. Dias, G.S.; Kitano, E.S.; Pagotto, A.H.; Sant’anna, S.V.S.; Rocha, M.M.; Zelanis, A.; Serrano, S.M. Individual variability in the venom proteome of juvenile Bothrops jararaca specimens. J. Proteome Res. 2013, 12, 4585–4598. [Google Scholar] [CrossRef] [PubMed]
  93. Durban, J.; Juárez, P.; Angulo, Y.; Lomonte, B.; Flores-Diaz, M.; Alape-Girón, A.; Sasa, M.; Sanz, L.; Gutiérrez, J.M.; Dopazo, J. Profiling the venom gland transcriptomes of Costa Rican snakes by 454 pyrosequencing. BMC Genom. 2011, 12, 259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Homsi-Brandeburgo, M.; Queiroz, L.; Santo-Neto, H.; Rodrigues-Simioni, L.; Giglio, J. Fractionation of Bothrops jararacussu snake venom: Partial chemical characterization and biological activity of bothropstoxin. Toxicon 1988, 26, 615–627. [Google Scholar] [CrossRef]
  95. Junqueira-de-Azevedo, I.L.; Bastos, C.M.V.; Ho, P.L.; Luna, M.S.; Yamanouye, N.; Casewell, N.R. Venom-related transcripts from Bothrops jararaca tissues provide novel molecular insights into the production and evolution of snake venom. Mol. Biol. Evol. 2015, 32, 754–766. [Google Scholar] [CrossRef] [PubMed]
  96. Kashima, S.; Roberto, P.G.; Soares, A.M.; Astolfi-Filho, S.; Pereira, J.O.; Giuliati, S.; Faria, M., Jr.; Xavier, M.A.; Fontes, M.R.; Giglio, J.R. Analysis of Bothrops jararacussu venomous gland transcriptome focusing on structural and functional aspects: I—gene expression profile of highly expressed phospholipases A 2. Biochimie 2004, 86, 211–219. [Google Scholar] [CrossRef] [PubMed]
  97. Mora-Obando, D.; Guerrero-Vargas, J.A.; Prieto-Sánchez, R.; Beltrán, J.; Rucavado, A.; Sasa, M.; Gutiérrez, J.M.; Ayerbe, S.; Lomonte, B. Proteomic and functional profiling of the venom of Bothrops ayerbei from Cauca, Colombia, reveals striking interspecific variation with Bothrops asper venom. J. Proteom. 2014, 96, 159–172. [Google Scholar] [CrossRef] [PubMed]
  98. Moura-da-Silva, A.M.; Furlan, M.S.; Caporrino, M.C.; Grego, K.F.; Portes-Junior, J.A.; Clissa, P.B.; Valente, R.H.; Magalhães, G.S. Diversity of metalloproteinases in Bothrops neuwiedi snake venom transcripts: Evidences for recombination between different classes of SVMPs. BMC Genet. 2011, 12, 94. [Google Scholar] [CrossRef] [PubMed]
  99. Serrano, S.M.; Oliveira, A.K.; Menezes, M.C.; Zelanis, A. The proteinase-rich proteome of Bothrops jararaca venom. Toxin Rev. 2014, 33, 169–184. [Google Scholar] [CrossRef]
  100. Sousa, L.F.; Nicolau, C.A.; Peixoto, P.S.; Bernardoni, J.L.; Oliveira, S.S.; Portes-Junior, J.A.; Mourão, R.H.V.; Lima-dos-Santos, I.; Sano-Martins, I.S.; Chalkidis, H.M. Comparison of phylogeny, venom composition and neutralization by antivenom in diverse species of Bothrops complex. PLoS Negl. Trop. Dis. 2013, 7, e2442. [Google Scholar] [CrossRef] [PubMed]
  101. Tashima, A.K.; Zelanis, A.; Kitano, E.S.; Ianzer, D.; Melo, R.L.; Rioli, V.; Sant’anna, S.S.; Schenberg, A.C.; Camargo, A.C.; Serrano, S.M. Peptidomics of three Bothrops snake venoms: Insights into the molecular diversification of proteomes and peptidomes. Mol. Cell. Proteom. 2012, 11, 1245–1262. [Google Scholar] [CrossRef] [PubMed]
  102. Valente, R.H.; Guimaraes, P.R.; Junqueira, M.; Neves-Ferreira, A.G.C.; Soares, M.R.; Chapeaurouge, A.; Trugilho, M.R.; Leon, I.R.; Rocha, S.L.; Oliveira-Carvalho, A.L. Bothrops insularis venomics: A proteomic analysis supported by transcriptomic-generated sequence data. J. Proteom. 2009, 72, 241–255. [Google Scholar] [CrossRef] [PubMed]
  103. Zelanis, A.; Andrade-Silva, D.; Rocha, M.M.; Furtado, M.F.; Serrano, S.M.; Junqueira-de-Azevedo, I.L.; Ho, P.L. A transcriptomic view of the proteome variability of newborn and adult Bothrops jararaca snake venoms. PLoS Negl. Trop. Dis. 2012, 6, e1554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Angulo, Y.; Chaves, E.; Alape, A.; Rucavado, A.; Marı́, J.; Lomonte, B. Isolation and Characterization of a Myotoxic Phospholipase A 2 from the Venom of the Arboreal Snake Bothriechis (Bothrops) schlegelii from Costa Rica. Arch. Biochem. Biophys. 1997, 339, 260–266. [Google Scholar] [CrossRef] [PubMed]
  105. Fernández, J.; Lomonte, B.; Sanz, L.; Angulo, Y.; Gutiérrez, J.M.; Calvete, J.J. Snake venomics of Bothriechis nigroviridis reveals extreme variability among palm pitviper venoms: Different evolutionary solutions for the same trophic purpose. J. Proteome Res. 2010, 9, 4234–4241. [Google Scholar] [CrossRef] [PubMed]
  106. Kuch, U.; Mebs, D.; Gutiérrez, J.M.; Freire, A. Biochemical and biological characterization of Ecuadorian pitviper venoms (genera Bothriechis, Bothriopsis, Bothrops and Lachesis). Toxicon 1996, 34, 714–717. [Google Scholar] [CrossRef]
  107. Corasolla Carregari, V.; Stuani Floriano, R.; Rodrigues-Simioni, L.; Winck, F.V.; Baldasso, P.A.; Ponce-Soto, L.A.; Marangoni, S. Biochemical, Pharmacological, and Structural Characterization of New Basic Bbil-TX from Bothriopsis bilineata Snake Venom. BioMed Res. Int. 2013, 2013. [Google Scholar] [CrossRef] [PubMed]
  108. Floriano, R.S.; Carregari, V.C.; de Abreu, V.A.; Kenzo-Kagawa, B.; Ponce-Soto, L.A.; da Cruz-Hoefling, M.A.; Hyslop, S.; Marangoni, S.; Rodrigues-Simioni, L. Pharmacological study of a new Asp49 phospholipase A 2 (Bbil-TX) isolated from Bothriopsis bilineata smargadina (forest viper) venom in vertebrate neuromuscular preparations. Toxicon 2013, 69, 191–199. [Google Scholar] [CrossRef] [PubMed]
  109. Porto, B.N.; Telli, C.A.; Dutra, T.P.; Alves, L.S.; Bozza, M.T.; Fin, C.A.; Thiesen, F.V.; Renner, M.F. Biochemical and biological characterization of the venoms of Bothriopsis bilineata and Bothriopsis taeniata (Serpentes: Viperidae). Toxicon 2007, 50, 270–277. [Google Scholar] [CrossRef] [PubMed]
  110. Rodrigues-Simioni, L.; Floriano, R.S.; Rostelato-Ferreira, S.; Sousa, N.C.; Marangoni, S.; Ponce-Soto, L.A.; Carregari, V.C.; Hyslop, S. Presynaptic action of Bothriopsis bilineata smargadina (forest viper) venom in vitro. Toxicon 2011, 58, 140–145. [Google Scholar] [CrossRef] [PubMed]
  111. Changeux, J.-P.; Kasai, M.; Lee, C.-Y. Use of a snake venom toxin to characterize the cholinergic receptor protein. Proc. Natl. Acad. Sci. USA 1970, 67, 1241–1247. [Google Scholar] [CrossRef] [PubMed]
  112. Kamiguti, A.; Zuzel, M.; Theakston, R. Snake venom metalloproteinases and disintegrins: Interactions with cells. Braz. J. Med. Biol. Res. 1998, 31, 7. [Google Scholar] [CrossRef]
  113. Serrano, S.M.; Kim, J.; Wang, D.; Dragulev, B.; Shannon, J.D.; Mann, H.H.; Veit, G.; Wagener, R.; Koch, M.; Fox, J.W. The Cysteine-rich Domain of Snake Venom Metalloproteinases Is a Ligand for von Willebrand Factor A Domains role in substrate targeting*. J. Biol. Chem. 2006, 281, 39746–39756. [Google Scholar] [CrossRef] [PubMed]
  114. Singh, G.; Kuc, R.E.; Maguire, J.J.; Fidock, M.; Davenport, A.P. Novel Snake Venom Ligand Dendroaspis Natriuretic Peptide Is Selective for Natriuretic Peptide Receptor-A in Human Heart Down regulation of Natriuretic Peptide Receptor-A in Heart Failure. Circ. Res. 2006, 99, 183–190. [Google Scholar] [CrossRef] [PubMed]
  115. Singh, G.; Maguire, J.; Kuc, R.; Skepper, J.; Fidock, M.; Davenport, A. Characterization of the snake venom ligand [125I]-DNP binding to natriuretic peptide receptor-A in human artery and potent DNP mediated vasodilatation. Br. J. Pharmacol. 2006, 149, 838–844. [Google Scholar] [CrossRef] [PubMed]
  116. Cushman, D.W.; Ondetti, M.A. History of the design of captopril and related inhibitors of angiotensin converting enzyme. Hypertension 1991, 17, 589–592. [Google Scholar] [CrossRef] [PubMed]
  117. Kemparaju, K.; Girish, K. Snake venom hyaluronidase: A therapeutic target. Cell Biochem. Funct. 2006, 24, 7–12. [Google Scholar] [CrossRef] [PubMed]
  118. Koh, D.; Armugam, A.; Jeyaseelan, K. Snake venom components and their applications in biomedicine. Cell. Mol. Life Sci. CMLS 2006, 63, 3030–3041. [Google Scholar] [CrossRef] [PubMed]
  119. Marcinkiewicz, C. Functional characteristic of snake venom disintegrins: Potential therapeutic implication. Curr. Pharm. Des. 2005, 11, 815–827. [Google Scholar] [CrossRef] [PubMed]
  120. Pal, S.K.; Gomes, A.; Dasgupta, S.; Gomes, A. Snake venom as therapeutic agents: From toxin to drug development. Indian J. Exp. Biol. 2002, 40, 1353–1358. [Google Scholar] [PubMed]
  121. Patlak, M. From viper’s venom to drug design: Treating hypertension. FASEB J. 2004, 18, 421e–421e. [Google Scholar] [CrossRef]
  122. Swenson, S.; Ramu, S.; Markland, F.S. Anti-angiogenesis and RGD-containing snake venom disintegrins. Curr. Pharm. Des. 2007, 13, 2860–2871. [Google Scholar] [CrossRef] [PubMed]
  123. Fox, J.W.; Serrano, S.M. Approaching the golden age of natural product pharmaceuticals from venom libraries: An overview of toxins and toxin-derivatives currently involved in therapeutic or diagnostic applications. Curr. Pharm. Des. 2007, 13, 2927–2934. [Google Scholar] [CrossRef] [PubMed]
  124. Calvete, J.J. Antivenomics and venom phenotyping: A marriage of convenience to address the performance and range of clinical use of antivenoms. Toxicon 2010, 56, 1284–1291. [Google Scholar] [CrossRef] [PubMed]
  125. Chippaux, J.P. Snake-bites: Appraisal of the global situation. Bull. World Health Organ. 1998, 76, 515. [Google Scholar] [PubMed]
  126. Gutiérrez, J.M.; Theakston, R.D.G.; Warrell, D.A. Confronting the neglected problem of snake bite envenoming: The need for a global partnership. PLoS Med. 2006, 3, e150. [Google Scholar] [CrossRef] [PubMed]
  127. Fox, J.W.; Serrano, S.M. Structural considerations of the snake venom metalloproteinases, key members of the M12 reprolysin family of metalloproteinases. Toxicon 2005, 45, 969–985. [Google Scholar] [CrossRef] [PubMed]
  128. Casewell, N.R.; Wagstaff, S.C.; Harrison, R.A.; Renjifo, C.; Wüster, W. Domain loss facilitates accelerated evolution and neofunctionalization of duplicate snake venom metalloproteinase toxin genes. Mol. Biol. Evol. 2011, 28, 2637–2649. [Google Scholar] [CrossRef] [PubMed]
  129. Gutierrez, J.M.; Lomonte, B.; Leon, G.; Alape-Giron, A.; Flores-Diaz, M.; Sanz, L.; Angulo, Y.; Calvete, J.J. Snake venomics and antivenomics: Proteomic tools in the design and control of antivenoms for the treatment of snakebite envenoming. J. Proteom. 2009, 72, 165–182. [Google Scholar] [CrossRef] [PubMed]
  130. Calvete, J.J. Proteomic tools against the neglected pathology of snake bite envenoming. Expert Rev. Proteom. 2011, 8, 739–758. [Google Scholar] [CrossRef] [PubMed]
  131. Campbell, J.A.; Lamar, W.W.; Brodie, E.D. The Venomous Reptiles of the Western Hemisphere; Comstock Publishing Associates: Ithaca, NY, USA, 2004. [Google Scholar]
  132. Martins, M.; Marques, O.A.; Sazima, I. Ecological and phylogenetic correlates of feeding habits in Neotropical pitvipers of the genus Bothrops. In Biology of the Vipers; Schuett, G.W., Höggren, M., Eds.; Eagle Mountain Publishing: Eagle Mountain, UT, USA, 2002; pp. 307–328. [Google Scholar]
  133. Mackessy, S. Venom composition in rattlesnakes: Trends and biological significance. In The Biology of Rattlesnakes; Hayes, W.K., Cardwell, M.D., Beaman, K.R., Bush, S.P., Eds.; Loma Linda University Press: Loma Linda, CA, USA, 2008; pp. 495–510. [Google Scholar]
  134. Gans, C. The feeding mechanism of snakes and its possible evolution. Am. Zool. 1961, 217–227. [Google Scholar] [CrossRef]
  135. Greene, H.W. Dietary correlates of the origin and radiation of snakes. Am. Zool. 1983, 23, 431–441. [Google Scholar] [CrossRef]
  136. Pough, F.H.; Groves, J.D. Specializations of the body form and food habits of snakes. Am. Zool. 1983, 23, 443–454. [Google Scholar] [CrossRef]
  137. Hampton, P.M.; Moon, B.R. Gape size, its morphological basis, and the validity of gape indices in western diamond-backed rattlesnakes (Crotalus atrox). J. Morphol. 2013, 274, 194–202. [Google Scholar] [CrossRef] [PubMed]
  138. Andrade, D.V.; Abe, A.S. Relationship of venom ontogeny and diet in Bothrops. Herpetologica 1999, 55, 200–204. [Google Scholar]
  139. MacArtney, J.M. Diet of the northern Pacific rattlesnake, Crotalus viridis oreganus, in British Columbia. Herpetologica 1989, 45, 299–304. [Google Scholar]
  140. Casasola, A.; Ramos-Cerrillo, B.; de Roodt, A.R.; Saucedo, A.C.; Chippaux, J.-P.; Alagón, A.; Stock, R.P. Paraspecific neutralization of the venom of African species of cobra by an equine antiserum against Naja melanoleuca: A comparative study. Toxicon 2009, 53, 602–608. [Google Scholar] [CrossRef] [PubMed]
  141. Ciscotto, P.H.; Rates, B.; Silva, D.A.; Richardson, M.; Silva, L.P.; Andrade, H.; Donato, M.F.; Cotta, G.A.; Maria, W.S.; Rodrigues, R.J. Venomic analysis and evaluation of antivenom cross-reactivity of South American Micrurus species. J. Proteom. 2011, 74, 1810–1825. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Taxonomical relationships of American pitvipers, based on Castoe 2006 [65].
Figure 1. Taxonomical relationships of American pitvipers, based on Castoe 2006 [65].
Toxins 08 00210 g001
Figure 2. 1D SDS page. MW = molecular weight marker; aur = Bothriechis aurifer; lat = Bothriechis lateralis; mar = Bothriechis marchi; sch = Bothriechis schlegelli; bil = Bothrops bilineata; tan = Bothrops taeniata; asp = Bothrops asper; neu = Bothrops neuwiedi bolivianus. Annotation indicates the dominant type in a region. However, other toxin types may also be present (see Supplementary File 1 for full annotation).
Figure 2. 1D SDS page. MW = molecular weight marker; aur = Bothriechis aurifer; lat = Bothriechis lateralis; mar = Bothriechis marchi; sch = Bothriechis schlegelli; bil = Bothrops bilineata; tan = Bothrops taeniata; asp = Bothrops asper; neu = Bothrops neuwiedi bolivianus. Annotation indicates the dominant type in a region. However, other toxin types may also be present (see Supplementary File 1 for full annotation).
Toxins 08 00210 g002
Figure 3. 2D SDS page analysis of (A) Bothriechis aurifer; (B) Bothriechis lateralis; (C) Bothriechis marchi; (D) Bothriechis schlegelli; (E) Bothrops bilineata; (F) Bothrops taeniata; (G) Bothrops asper and (H) Bothrops neuwiedi bolivianus. pI range is 3–10 (left to right) and molecular weight markers are as for Figure 2.
Figure 3. 2D SDS page analysis of (A) Bothriechis aurifer; (B) Bothriechis lateralis; (C) Bothriechis marchi; (D) Bothriechis schlegelli; (E) Bothrops bilineata; (F) Bothrops taeniata; (G) Bothrops asper and (H) Bothrops neuwiedi bolivianus. pI range is 3–10 (left to right) and molecular weight markers are as for Figure 2.
Toxins 08 00210 g003
Figure 4. Metalloprotease activity of venom—Metalloprotease activity of venom (10 ng/µL) was measured based on its ability to cleave a fluorogenic peptide substrate (Mca-PLGL-Dpa-AR-NH2, 10 µM final). aur = Bothriechis aurifer; lat = Bothriechis lateralis; mar = Bothriechis marchi; sch = Bothriechis schlegelii; bil = Bothrops bilineata; tan = Bothrops taeniata; asp = Bothrops asper; neu = Bothrops neuwiedi bolivianus.
Figure 4. Metalloprotease activity of venom—Metalloprotease activity of venom (10 ng/µL) was measured based on its ability to cleave a fluorogenic peptide substrate (Mca-PLGL-Dpa-AR-NH2, 10 µM final). aur = Bothriechis aurifer; lat = Bothriechis lateralis; mar = Bothriechis marchi; sch = Bothriechis schlegelii; bil = Bothrops bilineata; tan = Bothrops taeniata; asp = Bothrops asper; neu = Bothrops neuwiedi bolivianus.
Toxins 08 00210 g004
Figure 5. Crude venom zymography gel analysis using casein (A) and gelatin (B) as protein substrates. MW = molecular weight marker; aur = Bothriechis aurifer; lat = Bothriechis lateralis; mar = Bothriechis marchi; sch = Bothriechis schlegelii; bil = Bothrops bilineata; tan = Bothrops taeniata; asp = Bothrops asper; neu = Bothrops neuwiedi bolivianus.
Figure 5. Crude venom zymography gel analysis using casein (A) and gelatin (B) as protein substrates. MW = molecular weight marker; aur = Bothriechis aurifer; lat = Bothriechis lateralis; mar = Bothriechis marchi; sch = Bothriechis schlegelii; bil = Bothrops bilineata; tan = Bothrops taeniata; asp = Bothrops asper; neu = Bothrops neuwiedi bolivianus.
Toxins 08 00210 g005
Figure 6. Phospholipase A2 enzymatic activity profiling measured by means of absorbance over time at 1 µg for (A) Bothriechis species; (B) arboreal Bothrops species and (C) terrestrial Bothrops species.
Figure 6. Phospholipase A2 enzymatic activity profiling measured by means of absorbance over time at 1 µg for (A) Bothriechis species; (B) arboreal Bothrops species and (C) terrestrial Bothrops species.
Toxins 08 00210 g006
Table 1. Toxin types recovered by shotgun mass spectrometry, arboreality, and typical total body lengths of species under study. Body length information is from Campbell and Lamar (2004) [131].
Table 1. Toxin types recovered by shotgun mass spectrometry, arboreality, and typical total body lengths of species under study. Body length information is from Campbell and Lamar (2004) [131].
Toxin typesaurlatmarschbiltanaspneu
5’ nuc X XX
BPP/CNPX X X X
CRiSPX X XX
KallikreinXX XXXXX
Kunitz X
Lectin X X
LAAO XXXXXXX
PLA2 X X X
PLB X
PD XX
SVMP/DisX XX
Usual TL (cm)<70<80>80<60<70<100120–18060–70
Max TL (cm)10110096.897.9123175250100
ArborealYYYYYYNN
aur = Bothriechis aurifer; lat = Bothriechis lateralis; mar = Bothriechis marchi; sch = Bothriechis schlegelii; bil = Bothrops bilineata; tan = Bothrops taeniata; asp = Bothrops asper; neu = Bothrops neuwiedi bolivianus.

Share and Cite

MDPI and ACS Style

Debono, J.; Cochran, C.; Kuruppu, S.; Nouwens, A.; Rajapakse, N.W.; Kawasaki, M.; Wood, K.; Dobson, J.; Baumann, K.; Jouiaei, M.; et al. Canopy Venom: Proteomic Comparison among New World Arboreal Pit-Viper Venoms. Toxins 2016, 8, 210. https://doi.org/10.3390/toxins8070210

AMA Style

Debono J, Cochran C, Kuruppu S, Nouwens A, Rajapakse NW, Kawasaki M, Wood K, Dobson J, Baumann K, Jouiaei M, et al. Canopy Venom: Proteomic Comparison among New World Arboreal Pit-Viper Venoms. Toxins. 2016; 8(7):210. https://doi.org/10.3390/toxins8070210

Chicago/Turabian Style

Debono, Jordan, Chip Cochran, Sanjaya Kuruppu, Amanda Nouwens, Niwanthi W. Rajapakse, Minami Kawasaki, Kelly Wood, James Dobson, Kate Baumann, Mahdokht Jouiaei, and et al. 2016. "Canopy Venom: Proteomic Comparison among New World Arboreal Pit-Viper Venoms" Toxins 8, no. 7: 210. https://doi.org/10.3390/toxins8070210

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop