Snake Venomics of the Lancehead Pitviper Bothrops asper: Geographic, Individual, and Ontogenetic Variations

The habu genome reveals accelerated evolution of venom protein genes

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

Proteomic profile, biological activities and antigenic analysis of the venom from Bothriopsis bilineata smaragdina ("loro machaco"), a pitviper snake from Peru

In order to determine Bothriopsis bilineata smaragdina venom (BbsV) composition, proteomic approaches were performed. Venom components were analyzed by RP-HPLC, SDS- PAGE and nano LC on line with LTQ Orbitrap XL. Results showed a total of 189 identified proteins, grouped into 11 different subgroups, which include snake venom metalloproteinases (SVMPs, 54.67%), snake C-type lectins (Snaclecs, 15.78%), snake venom serine proteinases (SVSPs, 14.69%), cystein-rich secretory proteins (CRISP, 2.61%), phospholipases A₂ (PLA₂, 1.14%), phosphodiesterase (PDE, 1.17%), venom endothelial growth factor (VEGF, 1.06%) 5'nucleotidases (0.33%), L-amino acid oxidases (LAAOs, 0.28%) and other proteins. In vitro enzymatic activities (SVMP, SVSP, LAAO, Hyal and PLA₂) of BbsV were also analyzed. BbsV showed high SVSP activity but low PLA₂ activity, when compared to other Bothrops venoms. In vivo, BbsV induced hemorrhage and edema in mice and showed intraperitoneal median lethal dose (LD₅₀) of 92.74 (± 0.15) μg/20 g of mice. Furthermore, BbsV reduced cell viability when incubated with VERO cells. Peruvian and Brazilian bothropic antivenoms recognize BbsV proteins, as detected by ELISA and Western Blotting. Both antivenoms were able to neutralize in vivo edema and hemorrhage.

SIGNIFICANCE

In Peru, snakebite is a public health problem, especially in the rain forest, as a result of progressive colonization of this geographical area. This country is the second in Latin America, after Brazil, to exhibit the largest variety of venomous snakes. B. atrox and B. b. smaragdina snakes are sympatric species in Peruvian Amazon region and are responsible for approximately 95% of the envenomings reported in this region. B. b. smaragdina may cause a smaller share (3 to 38%) of those accidents, due to its arboreal habits, that make human encounters with these snakes less likely to happen. Despite B. b. smaragdina recognized medical importance, its venom composition and biological activities have been poorly studied. Furthermore, BbsV is not a component of the antigenic pool used to produce the corresponding Peruvian bothropic antivenom (P-BAV). Our results not only provide new insights on BbsV composition and biological activity, but also demonstrate that both P-BAV and B-BAV polyvalent antivenoms have a considerable recognition of proteins from BbsV and, more importantly, neutralized hemorrhage and edema, the main local effects of bothropic envenomation.

Venom Ontogeny in the Mexican Lance-Headed Rattlesnake (Crotalus polystictus)

As trophic adaptations, rattlesnake venoms can vary in composition depending on several intrinsic and extrinsic factors. Ontogenetic changes in venom composition have been documented for numerous species, but little is known of the potential age-related changes in many rattlesnake species found in México. In the current study, venom samples collected from adult and neonate Crotalus polystictus from Estado de México were subjected to enzymatic and electrophoretic analyses, toxicity assays (LD₅₀), and MALDI-TOF mass spectrometry, and a pooled sample of adult venom was analyzed by shotgun proteomics. Electrophoretic profiles of adult males and females were quite similar, and only minor sex-based variation was noted. However, distinct differences were observed between venoms from adult females and their neonate offspring. Several prominent bands, including P-I and P-III snake venom metalloproteinases (SVMPs) and disintegrins (confirmed by MS/MS) were present in adult venoms and absent/greatly reduced in neonate venoms. Age-dependent differences in SVMP, kallikrein-like, phospholipase A₂ (PLA₂), and L-amino acid oxidase (LAAO) activity levels were confirmed by enzymatic activity assays, and like many other rattlesnake species, venoms from adult snakes have higher SVMP activity than neonate venoms. Conversely, PLA₂ activity was approximately 2.5 × greater in venoms from neonates, likely contributing to the increased toxicity (neonate venom LD₅₀ = 4.5 μg/g) towards non-Swiss albino mice when compared to adult venoms (LD₅₀ = 5.5 μg/g). Thrombin-like (TLE) and phosphodiesterase activities did not vary significantly with age. A significant effect of sex (between adult male and adult female venoms) was also observed for SVMP, TLE, and LAAO activities. Analysis of pooled adult venom by LC-MS/MS identified 14 toxin protein families, dominated by bradykinin-inhibitory peptides, SVMPs (P-I, P-II and P-III), disintegrins, PLA₂s, C-type-lectins, CRiSPs, serine proteinases, and LAAOs (96% of total venom proteins). Neonate and adult C. polystictus in this population consume almost exclusively mammals, suggesting that age-based differences in composition are related to physical differences in prey (e.g., surface-to-volume ratio differences) rather than taxonomic differences between prey. Venoms from adult C. polystictus fit a Type I pattern (high SVMP activity, lower toxicity), which is characteristic of many larger-bodied rattlesnakes of North America.

Intraspecies geographical variability in the South American tigra mariposa (Bothrops venezuelensis Sandner 1952) snake venom activities

Bothrops venezuelensis snake venoms, from five localities in the North-Central Venezuelan regions, showed biochemical and haemostatic differences. In this study, bioactivities of B. venezuelensis venoms from different regions (Aragua state; Waraira Repano (Capital District); Baruta, La Boyera and Lagunetica (Miranda state)) were compared using both natural and synthetic substrates. The protein contents of these venoms were Lagunetica 89%, La Boyera 79%, Baruta 71%, Waraira Repano 68% and Aragua 64%. Toxic activities effects were: Intraperitoneal LD50s: Aragua-14 mg/kg; Waraira Repano-6.4 mg/kg; Baruta: 8.3 mg/kg; La Boyera-4.4 mg/kg; Lagunetica-16.2 mg/kg. The MHD results: Aragua-21.4 μg/mouse; Waraira Repano-2.5 μg/mouse; Baruta-1.2 μg/mouse; La Boyera-1.4 μg/mouse and Lagunetica-12 μg/mouse. The hide powder azure results: Aragua-1.24 U/mg; La Boyera-2.26 U/mg; Baruta-2.83 U/mg; Lagunetica-3.28 U/mg and Waraira Repano-5.77 U/mg. Esterase specific activity on BAEE results: Waraira Repano-666.66 U/mg; La Boyera-805.5 U/mg; Baruta-900.00 U/mg; Lagunetica-922.19 U/mg and Aragua-1960.67 U/mg. Casein zymography showed digestion bands in the molecular weight above 100 and at 66.2 and 21.5 kDa. Analysis of casein degradation by SDS-PAGE showed two different degradation patterns. Fibrinolytic activity (mm2/μg) on fibrin plates results: Aragua-6.07; Lagunetica-27.6; Waraira Repano-35.7; La Boyera-44.27 and Baruta-45.63. In the fibrinogenolytic assay, the five venoms completely degraded the α chain after 1 min of incubation. None of the venoms completely degraded the β and γ chains after 24 h incubation. The research indicated that venoms of B. venezuelensis of different geographic areas in Venezuela exhibit variances in composition and component concentrations; except the Aragua venom, all of them had high proteolytic activities.

Characterization of the Venom of C. d. cumanesis of Colombia: Proteomic Analysis and Antivenomic Study

The Colombian rattlesnake Crotalus durissus cumanensis is distributed in three geographic zones of the country: the Atlantic Coast, the upper valley of the Magdalena River, and the eastern plains of the Colombian Orinoquía. Its venom induces neurological symptoms, such as eyelid ptosis, myasthenic facies, and paralysis of the respiratory muscles, which can lead to death. Identification and analysis of C. d. cumanensis showed nine groups of proteins responsible for the neurotoxic effect, of which the crotoxin complex was the most abundant (64.71%). Immunorecognition tests of C. d. cumanensis showed that the use of a commercial antivenom manufactured in Mexico resulted in immunoreactivity.

Comparative Venomics of the Vipera ammodytes transcaucasiana and Vipera ammodytes montandoni from Turkey Provides Insights into Kinship

The Nose-horned Viper (Vipera ammodytes) is one of the most widespread and venomous snakes in Europe, which causes high frequent snakebite accidents. The first comprehensive venom characterization of the regional endemic Transcaucasian Nose-horned Viper (Vipera ammodytes transcaucasiana) and the Transdanubian Sand Viper (Vipera ammodytes montandoni) is reported employing a combination of intact mass profiling and bottom-up proteomics. The bottom-up analysis of both subspecies identified the major snake protein families of viper venoms. Furthermore, intact mass profiling revealed the presence of two tripeptidic metalloprotease inhibitors and their precursors. While previous reports applied multivariate analysis techniques to clarify the taxonomic status of the subspecies, an accurate classification of Vipera ammodytestranscaucasiana is still part of the ongoing research. The comparative analysis of the viper venoms on the proteome level reveals a close relationship between the Vipera ammodytes subspecies, which could be considered to clarify the classification of the Transcaucasian Nose-horned Viper. However, the slightly different ratio of some venom components could be indicating interspecific variations of the two studied subspecies or intraspecies alternations based on small sample size. Additionally, we performed a bioactivity screening with the crude venoms against several human cancerous and non-cancerous cell lines, which showed interesting results against a human breast adenocarcinoma epithelial cell line. Several fractions of Vipera a. transcaucasiana demonstrated a strong cytotoxic effect on triple negative MDA MB 231 breast cancer cells.

The Venom of the Spine-Bellied Sea Snake (Hydrophis curtus): Proteome, Toxin Diversity and Intraspecific Variation

The spine-bellied sea snake (Hydrophis curtus) is known to cause human deaths, yet its venom composition has not yet been proteomically characterised. An in-depth proteomic analysis was performed on H. curtus venom from two different seasons, January and June, corresponding to adults and subadults, respectively. Venoms from adult and subadult H. curtus individuals were compared using reversed-phase high-performance liquid chromatography (RP-HPLC), matrix-assisted laser desorption ionisation-time of flight (MALDI-TOF) mass spectrometry and liquid chromatography electrospray ionisation mass spectrometry (LC-ESI-MS) to detect intraspecific variation, and the molecular weight data obtained with ESI-MS were used to assess toxin diversity. RP-HPLC and LC-ESI-MS/MS were used to characterise the venom proteome and estimate the relative abundances of protein families present. The most abundant protein family in January and June venoms is phospholipase A₂ (PLA₂: January 66.7%; June 54.5%), followed by three-finger toxins (3FTx: January 30.4%; June 40.4%) and a minor component of cysteine-rich secretory proteins (CRISP: January 2.5%; June 5%). Trace amounts of snake venom metalloproteinases (SVMP), C-type lectins and housekeeping and regulatory proteins were also found. Although the complexity of the venom is low by number of families present, each family contained a more diverse set of isoforms than previously reported, a finding that may have implications for the development of next-generation sea snake antivenoms. Intraspecific variability was shown to be minor with one obvious exception of a 14,157-Da protein that was present in some January (adult) venoms, but not at all in June (subadult) venoms. There is also a greater abundance of short-chain neurotoxins in June (subadult) venom compared with January (adult) venom. These differences potentially indicate the presence of seasonal, ontogenetic or sexual variation in H. curtus venom.

Allopurinol attenuates acute kidney injury following Bothrops jararaca envenomation

Snakebites have been recognized as a neglected public health problem in several tropical and subtropical countries. Bothrops snakebites frequently complicate with acute kidney injury (AKI) with relevant morbidity and mortality. To date, the only treatment available for Bothrops envenomation is the intravenous administration of antivenom despite its several limitations. Therefore, the study of novel therapies in Bothrops envenomation is compelling. The aim of this study was to evaluate the protective effect of Allopurinol (Allo) in an experimental model of Bothrops jararaca venom (BJ)-associated AKI. Five groups of Wistar rats were studied: Sham, Allo, BJ, BJ+Allo, BJ+ipAllo. BJ (0.25 mg/kg) was intravenously injected during 40'. Saline at same dose and infusion rate was administered to Sham and Allo groups. Allo and BJ+Allo groups received Allo (300 mg/L) in the drinking water 7 days prior to Saline or BJ infusion respectively. BJ+ipAllo rats received intraperitoneal Allo (25 mg/Kg) 40' after BJ infusion. BJ rats showed markedly reduced glomerular filtration rate (GFR, inulin clearance) associated with intense renal vasoconstriction, hemolysis, hemoglobinuria, reduced glutathione and increased systemic and renal markers of nitro-oxidative stress (Nitrotyrosine). Allo ameliorated GFR, renal blood flow (RBF), renal vascular resistance and arterial lactate levels. In addition, Allo was associated with increased serum glutathione as well as reduced levels of plasma and renal Nitrotyrosine. Our data show that Allo attenuated BJ-associated AKI, reduced oxidative stress, improved renal hemodynamics and organ perfusion. It might represent a novel adjuvant approach for Bothrops envenomation, a new use for an old and widely available drug.

Cross-reactivity and cross-immunomodulation between venoms of the snakes Bothrops asper, Crotalus simus and Lachesis stenophrys, and its effect in the production of polyspecific antivenom for Central America

A mixture of the venoms of Bothrops asper, Crotalus simus and Lachesis stenophrys is used as immunogen to produce the polyspecific Central American antivenom (PoliVal-ICP). In this work, we studied the ability of each of these venoms to modulate the antibody response induced by the other two venoms included in the immunization mixture. For that, equine monospecific, bispecific and polyspecific antivenoms were prepared and compared regarding their ability to neutralize the phospholipase A₂, coagulant and lethal activities of each venom, and their anti-venom antibodies concentration. Results indicate that there is low cross-reactivity and cross-neutralization between venoms of B. asper, C. simus and L. stenophrys, hence justifying the use of all of them as immunogens for the production of the Central American antivenom. It was also found that the venom of B. asper reduces the anti-crotalic response while the venom of C. simus does not affect the anti-bothropic response. On the other hand, the venoms of B. asper and C. simus increase the anti-lachesic response, and L. stenoprhys venom reduced both the anti-bothropic and anti-crotalic responses. On the basis of these results, the immunization strategy can be adjusted by preventing or taking advantage of cross-immunomodulation between venoms, in order to maximize the antibody response towards all venoms. Immune responses can be improved by injecting horses with several immunogen mixtures, composed by one or two of the three venoms, and administering them at different times during the immunization, eventually generating a high titer against the three venoms. Our results suggest that addressing the issue of immunomodulation by venoms might improve antivenom manufacture worldwide.

Cross-recognition of a pit viper (Crotalinae) polyspecific antivenom explored through high-density peptide microarray epitope mapping

Snakebite antivenom is a 120 years old invention based on polyclonal mixtures of antibodies purified from the blood of hyper-immunized animals. Knowledge on antibody recognition sites (epitopes) on snake venom proteins is limited, but may be used to provide molecular level explanations for antivenom cross-reactivity. In turn, this may help guide antivenom development by elucidating immunological biases in existing antivenoms. In this study, we have identified and characterized linear elements of B-cell epitopes from 870 pit viper venom protein sequences by employing a high-throughput methodology based on custom designed high-density peptide microarrays. By combining data on antibody-peptide interactions with multiple sequence alignments of homologous toxin sequences and protein modelling, we have determined linear elements of antibody binding sites for snake venom metalloproteases (SVMPs), phospholipases A₂s (PLA₂s), and snake venom serine proteases (SVSPs). The studied antivenom antibodies were found to recognize linear elements in each of the three enzymatic toxin families. In contrast to a similar study of elapid (non-enzymatic) neurotoxins, these enzymatic toxins were generally not recognized at the catalytic active site responsible for toxicity, but instead at other sites, of which some are known for allosteric inhibition or for interaction with the tissue target. Antibody recognition was found to be preserved for several minor variations in the protein sequences, although the antibody-toxin interactions could often be eliminated completely by substitution of a single residue. This finding is likely to have large implications for the cross-reactivity of the antivenom and indicate that multiple different antibodies are likely to be needed for targeting an entire group of toxins in these recognized sites.

Although snakebite antivenom is a 120-year-old invention, saving lives and limbs of thousands of snakebite victims every year, little is known about the mechanisms and molecular interactions of how antivenoms neutralize snake toxins. Antivenoms are produced by immunizing large animals with cocktails of snake venoms resulting in antibodies recognizing toxic as well as non-toxic venom proteins to variable degrees. As a result, high doses of antivenom are needed for treating a snakebite victim, causing more severe adverse reactions due to a high burden of heterologous antivenom proteins. For the first time, we have characterized the antibody recognition sites on hundreds of pit viper toxins using high-throughput peptide microarray technology and an antivenom specific for three pit vipers inflicting a high number of bites in Central America. Most pit viper toxins are enzymes known to have a catalytic site important for toxicity. However, our results suggest that the employed antivenom generally does not target such sites, but instead inhibits toxicity by binding to alternative sites, possibly causing conformational shifts in the toxin structures or interference with toxin-target recognition. The identification of these toxin-specific recognition sites may explain why the antivenom is effective against certain snakebites from pit vipers whose venoms are not part of the immunization mixture.

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

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

Preclinical Evaluation of the Efficacy of Antivenoms for Snakebite Envenoming: State-of-the-Art and Challenges Ahead

Animal-derived antivenoms constitute the mainstay in the therapy of snakebite envenoming. The efficacy of antivenoms to neutralize toxicity of medically-relevant snake venoms has to be demonstrated through meticulous preclinical testing before their introduction into the clinical setting. The gold standard in the preclinical assessment and quality control of antivenoms is the neutralization of venom-induced lethality. In addition, depending on the pathophysiological profile of snake venoms, the neutralization of other toxic activities has to be evaluated, such as hemorrhagic, myotoxic, edema-forming, dermonecrotic, in vitro coagulant, and defibrinogenating effects. There is a need to develop laboratory assays to evaluate neutralization of other relevant venom activities. The concept of the 3Rs (Replacement, Reduction, and Refinement) in Toxinology is of utmost importance, and some advances have been performed in their implementation. A significant leap forward in the study of the immunological reactivity of antivenoms against venoms has been the development of “antivenomics”, which brings the analytical power of mass spectrometry to the evaluation of antivenoms. International partnerships are required to assess the preclinical efficacy of antivenoms against snake venoms in different regions of the world in order to have a detailed knowledge on the neutralizing profile of these immunotherapeutics.

Good vibrations: Assessing the stability of snake venom composition after researcher-induced disturbance in the laboratory

Phenotypic plasticity contributes to intraspecific variation in traits of many animal species. Venom is an integral trait to the success and survival of many snake species, and potential plasticity in venom composition is important to account for in the context of basic research as well as in human medicine for treating the various symptoms of snakebite and producing effective anti-venoms. Researchers may unknowingly induce changes in venom variation by subjecting snakes to novel disturbances and potential stressors. We explored phenotypic plasticity in snake venom composition over time in captive Pacific rattlesnakes (Crotalus oreganus) exposed to vibration treatment, compared to an undisturbed control group. Venom composition did not change significantly in response to vibration, nor was there a detectable effect of overall time in captivity, even though snakes re-synthesized venom stores while subjected to novel disturbance in the laboratory. This result indicates that venom composition is a highly repeatable phenotype over short time spans and that the composition of venom within adult individuals may be resistant to or unaffected by researcher-induced disturbance. On the other hand, the change in venom composition, measured as movement along the first principle component of venom phenotype space, was associated with baseline corticosterone (CORT) levels in the snakes. While differential forms of researcher-induced disturbance may not affect venom composition, significant changes in baseline CORT, or chronic stress, may affect the venom phenotype, and further investigations will be necessary to assess the nature of the relationship between CORT and venom protein expression.

Differential transcript profile of inhibitors with potential anti-venom role in the liver of juvenile and adult Bothrops jararaca snake

Background

Snakes belonging to the Bothrops genus are vastly distributed in Central and South America and are responsible for most cases of reported snake bites in Latin America. The clinical manifestations of the envenomation caused by this genus are due to three major activities—proteolytic, hemorrhagic and coagulant—mediated by metalloproteinases, serine proteinases, phospholipases A₂ and other toxic compounds present in snake venom. Interestingly, it was observed that snakes are resistant to the toxic effects of its own and other snake’s venoms. This natural immunity may occur due the absence of toxin target or the presence of molecules in the snake plasma able to neutralize such toxins.

Methods

In order to identify anti-venom molecules, we construct a cDNA library from the liver of B. jararaca snakes. Moreover, we analyzed the expression profile of four molecules—the already known anti-hemorrhagic factor Bj46a, one gamma-phospholipase A₂ inhibitor, one inter-alpha inhibitor and one C1 plasma protease inhibitor—in the liver of juvenile and adult snakes by qPCR.

Results

The results revealed a 30-fold increase of gamma-phospholipase A₂ inhibitor and a minor increase of the inter-alpha inhibitor (5-fold) and of the C1 inhibitor (3-fold) in adults. However, the Bj46a factor seems to be equally transcribed in adults and juveniles.

Discussion

The results suggest the up-regulation of different inhibitors observed in the adult snakes might be a physiological adaptation to the recurrent contact with their own and even other snake’s venoms throughout its lifespan. This is the first comparative analysis of ontogenetic variation of expression profiles of plasmatic proteins with potential anti-venom activities of the venomous snake B. jararaca. Furthermore, the present data contributes to the understanding of the natural resistance described in these snakes.

Discovery of small molecule inhibitors for the snake venom metalloprotease BaP1 using in silico and in vitro tests

Snakebites represent an important public health problem, with a great number of victims with permanent sequelae or fatal outcomes, particularly in rural, agriculturally active areas. The snake venom metalloproteases (SVMPs) are the principal proteins responsible for some clinically-relevant effects, such as local and systemic hemorrhage, dermonecrosis, and myonecrosis. Because of the difficulties in neutralizing them rapidly and locally by antivenoms, the search and design of small molecules as inhibitors of SVMPs are proposed. The Bothrops asper metalloprotease P1 (BaP1) is hereby used as a target protein and by High Throughput Virtual Screening (HTVS) approach, the free access virtual libraries: ZINC, PubChem and ChEMBL, were searched for potent small molecule inhibitors. Results from the aforementioned approaches provided strong evidences on the structural requirements for the efficient BaP1 inhibition such as the presence of the pyrimidine-2,4,6-trione moiety. The two proposed compounds have also shown excellent results in performed in vitro interaction studies against BaP1.

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.

Unraveling the Proteome Composition and Immuno-profiling of Western India Russell's Viper Venom for In-Depth Understanding of Its Pharmacological Properties, Clinical Manifestations, and Effective Antivenom Treatment

The proteome composition of western India (WI) Russell's viper venom (RVV) was correlated with pharmacological properties and pathological manifestations of RV envenomation. Proteins in the 5-19 and 100-110 kDa mass ranges were the most predominate (∼35.1%) and least abundant (∼3.4%) components, respectively, of WI RVV. Non-reduced SDS-PAGE indicated the occurrence of multiple subunits, non-covalent oligomers, self-aggregation, and/or interactions among the RVV proteins. A total of 55 proteins belonging to 13 distinct snake venom families were unambiguously identified by ESI-LC-MS/MS analysis. Phospholipase A2 (32.5%) and Kunitz-type serine protease inhibitors (12.5%) represented the most abundant enzymatic and non-enzymatic proteins, respectively. However, ATPase, ADPase, and hyaluronidase, detected by enzyme assays, were not identified by proteomic analysis owing to limitations in protein database deposition. Several biochemical and pharmacological properties of WI RVV were also investigated. Neurological symptoms exhibited by some RV-bite patients in WI may be correlated to the presence of neurotoxic phospholipase A2 enzymes and Kunitz-type serine protease inhibitor complex in this venom. Monovalent antivenom was found to be better than polyvalent antivenom in immuno-recognition and neutralization of the tested pharmacological properties and enzyme activities of WI RVV; nevertheless, both antivenoms demonstrated poor cross-reactivity and neutralization of pharmacological activities shown by low-molecular-mass proteins (<18 kDa) of this venom.

What killed Karl Patterson Schmidt? Combined venom gland transcriptomic, venomic and antivenomic analysis of the South African green tree snake (the boomslang), Dispholidus typus

Background

Non-front-fanged colubroid snakes comprise about two-thirds of extant ophidian species. The medical significance of the majority of these snakes is unknown, but at least five species have caused life-threatening or fatal human envenomings. However, the venoms of only a small number of species have been explored.

Methods

A combined venomic and venom gland transcriptomic approach was employed to characterise of venom of Dispholidus typus (boomslang), the snake that caused the tragic death of Professor Karl Patterson Schmidt. The ability of CroFab™ antivenom to immunocapture boomslang venom proteins was investigated using antivenomics.

Results

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families: three-finger toxin (3FTx); phospholipase A₂ (PLA₂); cysteine-rich secretory proteins (CRISP); snake venom (SV) serine proteinase (SP); C-type lectin-like (CTL); SV metalloproteinases (SVMPs); and disintegrin-like/cysteine-rich (DC) proteolytic fragments. CroFab™ antivenom efficiently immunodepleted some boomslang SVMPs.

Conclusions

The present work is the first to address the overall proteomic profile of D. typus venom. This study allowed us to correlate the toxin composition with the toxic activities of the venom. The antivenomic analysis suggested that the antivenom available at the time of the unfortunate accident could have exhibited at least some immunoreactivity against the boomslang SVMPs responsible for the disseminated intravascular coagulation syndrome that caused K.P. Schmidt's fatal outcome.

General significance

This study may stimulate further research on other non-front-fanged colubroid snake venoms capable of causing life-threatening envenomings to humans, which in turn should contribute to prevent fatal human accidents, such as that unfortunately suffered by K.P. Schmidt.

•

The venom proteome of Dispholidus typus (boomslang) is reported.

•

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families.

•

Boomslang venom proteome is dominated (75%) by snake venom PIII-metalloproteinases (PIII-SVMPs).

•

CroFab™ antivenom efficiently immunodepleted some boomslang PIII-SVMPs.

The binding effectiveness of anti-r-disintegrin polyclonal antibodies against disintegrins and PII and PIII metalloproteases: An immunological survey of type A, B and A+B venoms from Mohave rattlesnakes

Snake venoms are known to have different venom compositions and toxicity, but differences can also be found within populations of the same species contributing to the complexity of treatment of envenomated victims. One of the first well-documented intraspecies venom variations comes from the Mohave rattlesnake (Crotalus scutulatus scutulatus). Initially, three types of venoms were described; type A venom is the most toxic as a result of ~45% Mojave toxin in the venom composition, type B lacks the Mojave toxin but contains over 50% of snake venom metalloproteases (SVMPs). Also, type A+B venom contains a combination of Mojave toxin and SVMP. The use of an anti-disintegrin antibody in a simple Enzyme-Linked Immunosorbent Assay (ELISA) can be used to identify the difference between the venoms of the type A, B, and A+B Mohave rattlesnakes. This study implements the use of an anti-recombinant disintegrin polyclonal antibody (ARDPA) for the detection of disintegrins and ADAMs (a disintegrin and metalloproteases) in individual crude snake venoms of Mohave rattlesnakes (Crotalus scutulatus scutulatus) of varying geographical locations. After correlation with Western blots, coagulation activity and LD₅₀ data, it was determined that the antibody allows for a quick and cost-efficient identification of venom types.

Proteomic analysis of venom variability and ontogeny across the arboreal palm-pitvipers (genus Bothriechis)

Bothriechis is a genus of eleven currently recognized slender and arboreal venomous snakes, commonly called palm-pitvipers that range from southern Mexico to northern South America. Despite dietary studies suggesting that palm-pitvipers are generalists with an ontogenetic shift toward endothermic prey, venom proteomic analyses have revealed remarkable divergence between the venoms of the Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. To achieve a more complete picture of the venomic landscape across Bothriechis, the venom proteomes of biodiversity of the northern Middle American highland palm-pitvipers, B. thalassinus, B. aurifer, and B. bicolor from Guatemala, B. marchi from Honduras, and neonate Costa Rican B. lateralis and B. schlegelii, were investigated. B. thalassinus and B. aurifer venoms are comprised by similar toxin arsenals dominated by SVMPs (33-39% of the venom proteome), CTLs (11-16%), BPP-like molecules (10-13%), and CRISPs (5-10%), and are characterized by the absence of PLA₂ proteins. Conversely, the predominant (35%) components of B. bicolor are D49-PLA₂ molecules. The venom proteome of B. marchi is similar to B. aurifer and B. thalassinus in that it is rich in SVMPs and BPPs, but also contains appreciable amounts (14.3%) of PLA₂s. The major toxin family found in the venoms of both neonate B. lateralis and B. schlegelii, is serine proteinase (SVSP), comprising about 20% of their toxin arsenals. The venom of neonate B. schlegelii is the only palm-pitviper venom where relative high amounts of Kunitz-type (6.3%) and γPLA₂ (5.2%) inhibitors have been identified. Despite notable differences between their proteomes, neonate venoms are more similar to each other than to adults of their respective species. However, the ontogenetic changes taking place in the venom of B. lateralis strongly differ from those that occur in the venom of B. schlegelii. Thus, the ontogenetic change in B. lateralis produces a SVMP-rich venom, whereas in B. schlegelii the age-dependent compositional shift generates a PLA₂-rich venom. Overall, genus-wide venomics illustrate the high evolvability of palm-pitviper venoms. The integration of the pattern of venom variation across Bothriechis into a phylogenetic and biogeographic framework may lay the foundation for assessing, in future studies, the evolutionary path that led to the present-day variability of the venoms of palm-pitvipers.

SIGNIFICANCE

Bothriechis represents a monophyletic basal genus of eleven arboreal palm-pitvipers that range from southern Mexico to northern South America. Despite palm-pitvipers' putative status as diet generalists, previous proteomic analyses have revealed remarkable divergence between the venoms of Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. Our current proteomic study of Guatemalan species, B. thalassinus, B. aurifer, and B. bicolor, Honduran B. marchi, and neonate B. lateralis and B. schlegelii from Costa Rica was undertaken to deepen our understanding of the evolutionary pattern of venom proteome diversity across Bothriechis. Ancestral characters are often, but not always, preserved in an organism's development. Venoms of neonate B. lateralis and B. schlegelii are more similar to each other than to adults of their respective species, suggesting that the high evolvability of palm-pitviper venoms may represent an inherent feature of Bothriechis common ancestor. Our genus-wide data identified four nodes of venom phenotype differentiation across the phylogeny of Bothriechis. Integrated into a phylogenetic and biogeographic framework, the pattern of venom variation across Bothriechis may lay the groundwork to establish whether divergence was driven by selection for efficient resource exploitation in arboreal 'islands', thereby contributing to the ecological speciation of the genus.

Understanding Biological Roles of Venoms Among the Caenophidia: The Importance of Rear-Fanged Snakes

Snake venoms represent an adaptive trophic response to the challenges confronting a limbless predator for overcoming combative prey, and this chemical means of subduing prey shows several dominant phenotypes. Many front-fanged snakes, particularly vipers, feed on various vertebrate and invertebrate prey species, and some of their venom components (e.g., metalloproteinases, cobratoxin) appear to have been selected for "broad-brush" incapacitation of different prey taxa. Using proteomic and genomic techniques, the compositional diversity of front-fanged snakes is becoming well characterized; however, this is not the case for most rear-fanged colubroid snakes. Because these species consume a high diversity of prey, and because venoms are primarily a trophic adaptation, important clues for understanding specific selective pressures favoring venom component composition will be found among rear-fanged snake venoms. Rear-fanged snakes typically (but not always) produce venoms with lower complexity than front-fanged snakes, and there are even fewer dominant (and, arguably, biologically most relevant) venom protein families. We have demonstrated taxon-specific toxic effects, where lizards and birds show high susceptibility while mammals are largely unaffected, for both Old World and New World rear-fanged snakes, strongly indicating a causal link between toxin evolution and prey preference. New data are presented on myotoxin a, showing that the extremely rapid paralysis induced by this rattlesnake toxin is specific for rodents, and that myotoxin a is ineffectual against lizards. Relatively few rear-fanged snake venoms have been characterized, and basic natural history data are largely lacking, but directed sampling of specialized species indicates that novel compounds are likely among these specialists, particularly among those species feeding on invertebrate prey such as scorpions and centipedes. Because many of the more than 2200 species of colubroid snakes are rear-fanged, and many possess a Duvernoy's venom gland, understanding the nature of their venoms is foundational to understanding venom evolution in advanced snakes.

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

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

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

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

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

BACKGROUND

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

METHODOLOGY/PRINCIPAL FINDINGS

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

CONCLUSIONS/SIGNIFICANCE

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

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

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

Muscle Tissue Damage Induced by the Venom of Bothrops asper: Identification of Early and Late Pathological Events through Proteomic Analysis

The time-course of the pathological effects induced by the venom of the snake Bothrops asper in muscle tissue was investigated by a combination of histology, proteomic analysis of exudates collected in the vicinity of damaged muscle, and immunodetection of extracellular matrix proteins in exudates. Proteomic assay of exudates has become an excellent new methodological tool to detect key biomarkers of tissue alterations for a more integrative perspective of snake venom-induced pathology. The time-course analysis of the intracellular proteins showed an early presence of cytosolic and mitochondrial proteins in exudates, while cytoskeletal proteins increased later on. This underscores the rapid cytotoxic effect of venom, especially in muscle fibers, due to the action of myotoxic phospholipases A2, followed by the action of proteinases in the cytoskeleton of damaged muscle fibers. Similarly, the early presence of basement membrane (BM) and other extracellular matrix (ECM) proteins in exudates reflects the rapid microvascular damage and hemorrhage induced by snake venom metalloproteinases. The presence of fragments of type IV collagen and perlecan one hour after envenoming suggests that hydrolysis of these mechanically/structurally-relevant BM components plays a key role in the genesis of hemorrhage. On the other hand, the increment of some ECM proteins in the exudate at later time intervals is likely a consequence of the action of endogenous matrix metalloproteinases (MMPs) or of de novo synthesis of ECM proteins during tissue remodeling as part of the inflammatory reaction. Our results offer relevant insights for a more integrative and systematic understanding of the time-course dynamics of muscle tissue damage induced by B. asper venom and possibly other viperid venoms.

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

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

Ontogenetic Variation in Biological Activities of Venoms from Hybrids between Bothrops erythromelas and Bothrops neuwiedi Snakes

Lance-headed snakes are found in Central and South America, and they account for most snakebites in Brazil. The phylogeny of South American pitvipers has been reviewed, and the presence of natural and non-natural hybrids between different species of Bothrops snakes demonstrates that reproductive isolation of several species is still incomplete. The present study aimed to analyze the biological features, particularly the thrombin-like activity, of venoms from hybrids born in captivity, from the mating of a female Bothrops erythromelas and a male Bothrops neuwiedi, two species whose venoms are known to display ontogenetic variation. Proteolytic activity on azocoll and amidolytic activity on N-benzoyl-DL-arginine-p-nitroanilide hydrochloride (BAPNA) were lowest when hybrids were 3 months old, and increased over body growth, reaching values similar to those of the father when hybrids were 12 months old. The clotting activity on plasma diminished as hybrids grew; venoms from 3- and 6-months old hybrids showed low clotting activity on fibrinogen (i.e., thrombin-like activity), like the mother venom, and such activity was detected only when hybrids were older than 1 year of age. Altogether, these results point out that venom features in hybrid snakes are genetically controlled during the ontogenetic development. Despite the presence of the thrombin-like enzyme gene(s) in hybrid snakes, they are silenced during the first six months of life.

Snake Venomics and Antivenomics of Bothrops diporus, a Medically Important Pitviper in Northeastern Argentina

Snake species within genus Bothrops are responsible for more than 80% of the snakebites occurring in South America. The species that cause most envenomings in Argentina, B. diporus, is widely distributed throughout the country, but principally found in the Northeast, the region with the highest rates of snakebites. The venom proteome of this medically relevant snake was unveiled using a venomic approach. It comprises toxins belonging to fourteen protein families, being dominated by PI- and PIII-SVMPs, PLA₂ molecules, BPP-like peptides, L-amino acid oxidase and serine proteinases. This toxin profile largely explains the characteristic pathophysiological effects of bothropic snakebites observed in patients envenomed by B. diporus. Antivenomic analysis of the SAB antivenom (Instituto Vital Brazil) against the venom of B. diporus showed that this pentabothropic antivenom efficiently recognized all the venom proteins and exhibited poor affinity towards the small peptide (BPPs and tripeptide inhibitors of PIII-SVMPs) components of the venom.

A novel synthetic quinolinone inhibitor presents proteolytic and hemorrhagic inhibitory activities against snake venom metalloproteases

Metalloproteases play a fundamental role in snake venom envenomation inducing hemorrhagic, fibrigen(ogen)olytic and myotoxic effects in their victims. Several snake venoms, such as those from the Bothrops genus, present important local effects which are not efficiently neutralized by conventional serum therapy. Consequently, these accidents may result in permanent sequelae and disability, creating economic and social problems, especially in developing countries, leading the attention of the World Health Organization that considered ophidic envenomations a neglected tropical disease. Aiming to produce an efficient inhibitor against bothropic venoms, we synthesized different molecules classified as quinolinones - a group of low-toxic chemical compounds widely used as antibacterial and antimycobacterial drugs - and tested their inhibitory properties against hemorrhage caused by bothropic venoms. The results from this initial screening indicated the molecule 2-hydroxymethyl-6-methoxy-1,4-dihydro-4-quinolinone (Q8) was the most effective antihemorrhagic compound among all of the assayed synthetic quinolinones. Other in vitro and in vivo experiments showed this novel compound was able to inhibit significantly the hemorrhagic and/or proteolytic activities of bothropic crude venoms and isolated snake venom metalloproteases (SVMPs) even at lower concentrations. Docking and molecular dynamic simulations were also performed to get insights into the structural basis of Q8 inhibitory mechanism against proteolytic and hemorrhagic SVMPs. These structural studies demonstrated that Q8 may form a stable complex with SVMPs, impairing the access of substrates to the active sites of these toxins. Therefore, both experimental and structural data indicate that Q8 compound is an interesting candidate for antiophidic therapy, particularly for the treatment of the hemorrhagic and necrotic effects induced by bothropic venoms.

Effects of geographical heterogeneity in species interactions on the evolution of venom genes

Geographical heterogeneity in the composition of biotic interactions can create a mosaic of selection regimes that may drive the differentiation of phenotypes that operate at the interface of these interactions. Nonetheless, little is known about effects of these geographical mosaics on the evolution of genes encoding traits associated with species interactions. Predatory marine snails of the family Conidae use venom, a cocktail of conotoxins, to capture prey. We characterized patterns of geographical variation at five conotoxin genes of a vermivorous species, Conus ebraeus, at Hawaii, Guam and American Samoa, and evaluated how these patterns of variation are associated with geographical heterogeneity in prey utilization. All populations show distinct patterns of prey utilization. Three 'highly polymorphic' conotoxin genes showed significant geographical differences in allelic frequency, and appear to be affected by different modes of selection among populations. Two genes exhibited low levels of diversity and a general lack of differentiation among populations. Levels of diversity of 'highly polymorphic' genes exhibit a positive relationship with dietary breadth. The different patterns of evolution exhibited by conotoxin genes suggest that these genes play different roles in prey capture, and that some genes are more greatly affected by differences in predator-prey interactions than others. Moreover, differences in dietary breadth appear to have a greater influence on the differentiation of venoms than differences in the species of prey.

Label-Free (XIC) Quantification of Venom Procoagulant and Neurotoxin Expression in Related Australian Elapid Snakes Gives Insight into Venom Toxicity Evolution

This study demonstrates a direct role of venom protein expression alteration in the evolution of snake venom toxicity. Avian skeletal muscle contractile response to exogenously administered acetylcholine is completely inhibited upon exposure to South Australian and largely preserved following exposure to Queensland eastern brown snake Pseudonaja textilis venom, indicating potent postsynaptic neurotoxicity of the former and lack thereof of the latter venom. Label-free quantitative proteomics reveals extremely large differences in the expression of postsynaptic three-finger α-neurotoxins in these venoms, explaining the difference in the muscle contractile response and suggesting that the type of toxicity induced by venom can be modified by altered expression of venom proteins. Furthermore, the onset of neuromuscular paralysis in the rat phrenic nerve-diaphragm preparation occurs sooner upon exposure to the venom (10 μg/mL) with high expression of α-neurotoxins than the venoms containing predominately presynaptic β-neurotoxins. The study also finds that the onset of rat plasma coagulation is faster following exposure to the venoms with higher expression of venom prothrombin activator subunits. This is the first quantitative proteomic study that uses extracted ion chromatogram peak areas (MS1 XIC) of distinct homologous tryptic peptides to directly show the differences in the expression of venom proteins.

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

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

Post-transcriptional Mechanisms Contribute Little to Phenotypic Variation in Snake Venoms

Protein expression is a major link in the genotype–phenotype relationship, and processes affecting protein abundances, such as rates of transcription and translation, could contribute to phenotypic evolution if they generate heritable variation. Recent work has suggested that mRNA abundances do not accurately predict final protein abundances, which would imply that post-transcriptional regulatory processes contribute significantly to phenotypes. Post-transcriptional processes also appear to buffer changes in transcriptional patterns as species diverge, suggesting that the transcriptional changes have little or no effect on the phenotypes undergoing study. We tested for concordance between mRNA and protein expression levels in snake venoms by means of mRNA-seq and quantitative mass spectrometry for 11 snakes representing 10 species, six genera, and three families. In contrast to most previous work, we found high correlations between venom gland transcriptomes and venom proteomes for 10 of our 11 comparisons. We tested for protein-level buffering of transcriptional changes during species divergence by comparing the difference between transcript abundance and protein abundance for three pairs of species and one intraspecific pair. We found no evidence for buffering during divergence of our three species pairs but did find evidence for protein-level buffering for our single intraspecific comparison, suggesting that buffering, if present, was a transient phenomenon in venom divergence. Our results demonstrated that post-transcriptional mechanisms did not contribute significantly to phenotypic evolution in venoms and suggest a more prominent and direct role for cis-regulatory evolution in phenotypic variation, particularly for snake venoms.

Isolation and characterization of four medium-size disintegrins from the venoms of Central American viperid snakes of the genera Atropoides, Bothrops, Cerrophidion and Crotalus

Four disintegrins were isolated from the venoms of the Central American viperid snakes Atropoides mexicanus (atropoimin), Bothrops asper (bothrasperin), Cerrophidion sasai (sasaimin), and Crotalus simus (simusmin). Purifications were performed by reverse-phase HPLC. The four disintegrins have biochemical characteristics, i.e. molecular mass and location of Cys, which allow their classification within the group of medium-size disintegrins. All of them present the canonical RGD sequence, which determines their interaction with integrins in cell membranes. The disintegrins inhibited ADP and collagen-induced human platelet aggregation, with similar IC50s in the nM range. In addition, disintegrins inhibited the adhesion of an endothelial cell line and a melanoma cell line to the extracellular matrix proteins type I collagen, laminin, fibronectin, and vitronectin, albeit showing variable ability to exert this activity. This study expands the inventory of this family of viperid venom proteins, and reports, for the first time, disintegrins from the venoms of species of the genera Atropoides and Cerrophidion.