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

Combined venomics, venom gland transcriptomics, bioactivities, and antivenomics of two Bothrops jararaca populations from geographic isolated regions within the Brazilian Atlantic rainforest

Bothrops jararaca is a slender and semi-arboreal medically relevant pit viper species endemic to tropical and subtropical forests in southern Brazil, Paraguay, and northern Argentina (Misiones). Within its geographic range, it is often abundant and is an important cause of snakebite. Although no subspecies are currently recognized, geographic analyses have revealed the existence of two well-supported B. jararaca clades that diverged during the Pliocene ~3.8Mya and currently display a southeastern (SE) and a southern (S) Atlantic rainforest (Mata Atlântica) distribution. The spectrum, geographic variability, and ontogenetic changes of the venom proteomes of snakes from these two B. jararaca phylogroups were investigated applying a combined venom gland transcriptomic and venomic analysis. Comparisons of the venom proteomes and transcriptomes of B. jararaca from the SE and S geographic regions revealed notable interpopulational variability that may be due to the different levels of population-specific transcriptional regulation, including, in the case of the southern population, a marked ontogenetic venom compositional change involving the upregulation of the myotoxic PLA2 homolog, bothropstoxin-I. This population-specific marker can be used to estimate the proportion of venom from the southern population present in the B. jararaca venom pool used for the Brazilian soro antibotrópico (SAB) antivenom production. On the other hand, the southeastern population-specific D49-PLA2 molecules, BinTX-I and BinTX-II, lend support to the notion that the mainland ancestor of Bothrops insularis was originated within the same population that gave rise to the current SE B. jararaca phylogroup, and that this insular species endemic to Queimada Grande Island (Brazil) expresses a pedomorphic venom phenotype. Mirroring their compositional divergence, the two geographic B. jararaca venom pools showed distinct bioactivity profiles. However, the SAB antivenom manufactured in Vital Brazil Institute neutralized the lethal effect of both venoms to a similar extent. In addition, immobilized SAB antivenom immunocaptured most of the venom components of the venoms of both B. jararaca populations, but did not show immunoreactivity against vasoactive peptides. The Costa Rican bothropic-crotalic-lachesic (BCL) antivenom showed the same lack of reactivity against vasoactive peptides but, in addition, was less efficient immunocapturing PI- and PIII-SVMPs from the SE venom, and bothropstoxin-I, a CRISP molecule, and a D49-PLA2 from the venom of the southern B. jararaca phylogroup. The remarkable paraspecificity exhibited by the Brazilian and the Costa Rican antivenoms indicates large immunoreactive epitope conservation across the natural history of Bothrops, a genus that has its roots in the middle Miocene. This article is part of a Special Issue entitled: Omics Evolutionary Ecolog.

Statistical analysis of the individual variability of 1D protein profiles as a tool in ecology: an application to parasitoid venom

Understanding the forces that shape eco-evolutionary patterns often requires linking phenotypes to genotypes, allowing characterization of these patterns at the molecular level. DNA-based markers are less informative in this aim compared to markers associated with gene expression and, more specifically, with protein quantities. The characterization of eco-evolutionary patterns also usually requires the analysis of large sample sizes to accurately estimate interindividual variability. However, the methods used to characterize and compare protein samples are generally expensive and time-consuming, which constrains the size of the produced data sets to few individuals. We present here a method that estimates the interindividual variability of protein quantities based on a global, semi-automatic analysis of 1D electrophoretic profiles, opening the way to rapid analysis and comparison of hundreds of individuals. The main original features of the method are the in silico normalization of sample protein quantities using pictures of electrophoresis gels at different staining levels, as well as a new method of analysis of electrophoretic profiles based on a median profile. We demonstrate that this method can accurately discriminate between species and between geographically distant or close populations, based on interindividual variation in venom protein profiles from three endoparasitoid wasps of two different genera (Psyttalia concolor, Psyttalia lounsburyi and Leptopilina boulardi). Finally, we discuss the experimental designs that would benefit from the use of this method.

The red seaweed Plocamium brasiliense shows anti-snake venom toxic effects

BACKGROUND

Snakebite is considered a neglected tropical disease by the World Health Organization. In Brazil, about 70% of the envenomation cases are caused by Bothrops snakes. Its venom may provoke hemorrhage, pain, necrosis, hemolysis, renal or cardiac failure and even death in victims. Since commercial antivenom does not efficiently neutralize the local toxic effects of venoms, natural products have been tested in order to provide alternative or complementary treatment to serum therapy. Therefore, the present study aimed to evaluate the ability of the seaweed Plocamium brasiliense and its active derivatives to neutralize hemorrhagic, edematogenic, hemolytic, coagulant and proteolytic activities of B. jararaca venom.

METHODS

Specimens of P. brasiliense were collected in Rio de Janeiro state, Brazil, dried and submitted to oil extraction using four solvents of increasing polarities, n-hexane (HEX), dichloromethane (DCM), ethyl acetate (ETA) and hydroalcoholic solution (HYD). The solvents were evaporated, yielding HEX, DCM, ETA and HYD extracts. Further, all extracts were dissolved in dimethylsulfoxide. In addition, two monoterpenes (8-bromo-3,4,7-trichloro-3,7-dimethyl-1E, 5E-octadiene and 1,8-dibromo-3,4,7-trichloro-3,7-dimethyl-1E, 5E-octadiene) and a cholesterol fraction were isolated from the extract of P. brasiliense prepared in hexane. Algal samples were incubated for 30 minutes with B. jararaca venom, and then tested for lethality; hemorrhagic, edematogenic, hemolytic, coagulant and proteolytic effects.

RESULTS

Most of the algal extracts inhibited the toxic effects with different potencies. The DCM extract was the most effective, since it inhibited all types of toxic activity. On the other hand, the HYD extract failed to inhibit any effect. Moreover, the isolated products inhibited proteolysis and protected mice from hemorrhage in 30% of the cases, whereas 8-bromo-3,4,7-trichloro-3,7-dimethyl-1E, 5E-octadiene inhibited 100% and 20% of the hemorrhagic and proteolytic activities, respectively. None of the algal products were toxic to mice.

CONCLUSION

Seaweeds may be a promising source of inhibitors against toxic effects caused by B. jararaca envenomation, which may contribute to antivenom treatment.

Omics meets biology: application to the design and preclinical assessment of antivenoms

Snakebite envenoming represents a neglected tropical disease that has a heavy public health impact worldwide, mostly affecting poor people involved in agricultural activities in Africa, Asia, Latin America and Oceania. A key issue that complicates the treatment of snakebite envenomings is the poor availability of the only validated treatment for this disease, antivenoms. Antivenoms can be an efficacious treatment for snakebite envenoming, provided they are safe, effective, affordable, accessible and administered appropriately. The shortage of antivenoms in various regions, particularly in Sub-Saharan Africa and some parts of Asia, can be significantly alleviated by optimizing the use of current antivenoms and by the generation of novel polyspecific antivenoms having a wide spectrum of efficacy. Complementing preclinical testing of antivenom efficacy using in vivo and in vitro functional neutralization assays, developments in venomics and antivenomics are likely to revolutionize the design and preclinical assessment of antivenoms by being able to test new antivenom preparations and to predict their paraspecific neutralization to the level of species-specific toxins.

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

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

Extraction of venom and venom gland microdissections from spiders for proteomic and transcriptomic analyses

Venoms are chemically complex secretions typically comprising numerous proteins and peptides with varied physiological activities. Functional characterization of venom proteins has important biomedical applications, including the identification of drug leads or probes for cellular receptors. Spiders are the most species rich clade of venomous organisms, but the venoms of only a few species are well-understood, in part due to the difficulty associated with collecting minute quantities of venom from small animals. This paper presents a protocol for the collection of venom from spiders using electrical stimulation, demonstrating the procedure on the Western black widow (Latrodectus hesperus). The collected venom is useful for varied downstream analyses including direct protein identification via mass spectrometry, functional assays, and stimulation of venom gene expression for transcriptomic studies. This technique has the advantage over protocols that isolate venom from whole gland homogenates, which do not separate genuine venom components from cellular proteins that are not secreted as part of the venom. Representative results demonstrate the detection of known venom peptides from the collected sample using mass spectrometry. The venom collection procedure is followed by a protocol for dissecting spider venom glands, with results demonstrating that this leads to the characterization of venom-expressed proteins and peptides at the sequence level.

Immunological cross-reactivity and neutralisation of European viper venoms with the monospecific Vipera berus antivenom ViperaTAb

Medically important cases of snakebite in Europe are predominately caused by European vipers of the genus Vipera. The mainstay of snakebite therapy is polyclonal antibody therapy, referred to as antivenom. Here we investigate the capability of the monospecific V. berus antivenom, ViperaTAb^®, to cross-react with, and neutralise lethality induced by, a variety of European vipers. Using ELISA and immunoblotting, we find that ViperaTAb^® antibodies recognise and bind to the majority of toxic components found in the venoms of the Vipera species tested at comparably high levels to those observed with V. berus. Using in vivo pre-clinical efficacy studies, we demonstrate that ViperaTAb^® effectively neutralises lethality induced by V. berus, V. aspis, V. ammodytes and V. latastei venoms and at much higher levels than those outlined by regulatory pharmacopoeial guidelines. Notably, venom neutralisation was found to be superior to (V. berus, V. aspis and V. latastei), or as equally effective as (V. ammodytes), the monospecific V. ammodytes “Zagreb antivenom”, which has long been successfully used for treating European snake envenomings. This study suggests that ViperaTAb^® may be a valuable therapeutic product for treating snakebite by a variety of European vipers found throughout the continent.

A neutralizing recombinant single chain antibody, scFv, against BaP1, A P-I hemorrhagic metalloproteinase from Bothrops asper snake venom

BaP1 is a P-I class snake venom metalloproteinase (SVMP) relevant in the local tissue damage associated with envenomings by Bothrops asper, a medically important snake species in Central America and parts of South and North America. The main treatment for these accidents is the passive immunotherapy using antibodies raised in horses. In order to obtain more specific and batch-to-batch consistent antivenons, recombinant antibodies are considered a good option compared to animal immunization. We constructed a recombinant single chain variable fragment (scFv) from a monoclonal antibody against BaP1 (MABaP1) formerly secreted by a hybridoma clone. This recombinant antibody was cloned into pMST3 vector in fusion with SUMO protein and contains VH and VL domains linked by a flexible (G4S)3 polypeptide (scFvBaP1). The aim of this work was to produce scFvBaP1 and to evaluate its potential concerning the neutralization of biologically important activities of BaP1. The cytoplasmic expression of this construct was successfully achieved in C43 (DE3) bacteria. Our results showed that scFvBaP1-SUMO fusion protein presented an electrophoretic band of around 43 kDa from which SUMO alone corresponded to 13.6 kDa, and only the scFv was able to recognize BaP1 as well as the whole venom by ELISA. In contrast, neither an irrelevant scFv anti-LDL nor its MoAb partner recognized it. BaP1-induced fibrinolysis was significantly neutralized by scFvBaP1, but not by SUMO, in a concentration-dependent manner. In addition, scFvBaP1, as well as MaBaP1, completely neutralized in vivo hemorrhage, muscle necrosis, and inflammation induced by the toxin. Docking analyses revealed possible modes of interaction of the recombinant antibody with BaP1. Our data showed that scFv recognized BaP1 and whole B. asper venom, and neutralized biological effects of this SVMP. This scFv antibody can be used for understanding the molecular mechanisms of neutralization of SVMPs, and for exploring the potential of recombinant antibody fragments for improving the neutralization of local tissue damage in snakebite envenoming.

Two color morphs of the pelagic yellow-bellied sea snake, Pelamis platura, from different locations of Costa Rica: snake venomics, toxicity, and neutralization by antivenom

The yellow-bellied sea snake, Pelamis platura, is the most broadly distributed snake species. Despite being endowed with a highly lethal venom, a proteomic analysis of its toxin composition was unavailable. The venoms of specimens collected in Golfo de Papagayo and Golfo Dulce (Costa Rica), where two distinctive color morphs occur, were chromatographically compared. The latter inhabits a fjord-like gulf where the transit of oceanic sea snakes into and from the basin is restricted, thus possibly affecting gene flow. RP-HPLC evidenced a conserved venom protein profile in both populations, despite their divergent color phenotypes. Following a trend observed in other sea snakes, P. platura venom is relatively simple, being composed of proteins of the three-finger toxin (3FTx), phospholipase A2 (PLA2), cysteine-rich secretory protein (CRISP), 5'-nucleotidase, and metalloproteinase families. The first three groups represent 49.9%, 32.9%, and 9.1% of total venom protein, respectively. The most abundant component (~26%) is pelamitoxin (P62388), a short-chain 3FTx, followed by a major basic PLA2 (~20%) and a group of three isoforms of CRISPs (~9%). Whereas isolated pelamitoxin was highly lethal to mice, neither the PLA2 nor the CRISP fraction caused death. However, the PLA2 rapidly increased plasma creatine kinase activity after intramuscular injection, indicating its myotoxic action. Differing from myotoxic PLA2s of viperids, this PLA2 was not cytolytic to murine myogenic cells in vitro, suggesting possible differences in its mechanism of action. The median lethal dose (LD50) estimates for P. platura crude venom in mice and in three species of fishes did not differ significantly. The sea snake antivenom manufactured by CSL Ltd. (Australia), which uses Enhydrina schistosa as immunogen, cross-recognized the three major components of P. platura venom and, accordingly, neutralized the lethal activity of crude venom and pelamitoxin, therefore being of potential usefulness in the treatment of envenomations by this species.

BIOLOGICAL SIGNIFICANCE

Integrative analyses of animal venoms that combine the power of proteomics (venomics) with the characterization of their functional and immunological properties are significantly expanding knowledge on these remarkable bioweapons, both from a basic and a medical perspective. Costa Rica harbors a unique population of the yellow-bellied sea snake, Pelamis platura, that is restricted to a fjord-like gulf (Golfo Dulce). This population differs markedly from oceanic populations found elsewhere along the Pacific coast of this country, by presenting a patternless bright yellow coloration, instead of the typical bicolored or tricolored pattern of this species. It has been suggested that the dominance of this yellow-morph in Golfo Dulce might reflect gene flow restrictions, caused by the oceanographic conditions at this location. The present study demonstrates that the remarkable phenotypic variation between the two color morphs inhabiting Golfo Dulce and Golfo de Papagayo, respectively, is not associated with differences in the expression of venom components, as shown by their conserved RP-HPLC profiles. Proteomic analysis revealed the relatively simple toxin composition of P. platura venom, which contains three predominant types of proteins: three-finger toxins (protein abundance: 49.9%), phospholipases A2 (32.9%), and cysteine-rich secretory proteins (9.1%), together with few minor components. Further, the involvement of these most abundant proteins in the toxic effects of the venom, and their cross-recognition and neutralization by a sea snake antivenom produced against the venom of Enhydrina schistosa, were analyzed.

Venomous snakes of Costa Rica: biological and medical implications of their venom proteomic profiles analyzed through the strategy of snake venomics

In spite of its small territory of ~50,000km(2), Costa Rica harbors a remarkably rich biodiversity. Its herpetofauna includes 138 species of snakes, of which sixteen pit vipers (family Viperidae, subfamily Crotalinae), five coral snakes (family Elapidae, subfamily Elapinae), and one sea snake (Family Elapidae, subfamily Hydrophiinae) pose potential hazards to human and animal health. In recent years, knowledge on the composition of snake venoms has expanded dramatically thanks to the development of increasingly fast and sensitive analytical techniques in mass spectrometry and separation science applied to protein characterization. Among several analytical strategies to determine the overall protein/peptide composition of snake venoms, the methodology known as 'snake venomics' has proven particularly well suited and informative, by providing not only a catalog of protein types/families present in a venom, but also a semi-quantitative estimation of their relative abundances. Through a collaborative research initiative between Instituto de Biomedicina de Valencia (IBV) and Instituto Clodomiro Picado (ICP), this strategy has been applied to the study of venoms of Costa Rican snakes, aiming to obtain a deeper knowledge on their composition, geographic and ontogenic variations, relationships to taxonomy, correlation with toxic activities, and discovery of novel components. The proteomic profiles of venoms from sixteen out of the 22 species within the Viperidae and Elapidae families found in Costa Rica have been reported so far, and an integrative view of these studies is hereby presented. In line with other venomic projects by research groups focusing on a wide variety of snakes around the world, these studies contribute to a deeper understanding of the biochemical basis for the diverse toxic profiles evolved by venomous snakes. In addition, these studies provide opportunities to identify novel molecules of potential pharmacological interest. Furthermore, the establishment of venom proteomic profiles offers a fundamental platform to assess the detailed immunorecognition of individual proteins/peptides by therapeutic or experimental antivenoms, an evolving methodology for which the term 'antivenomics' was coined (as described in an accompanying paper in this special issue).

BIOLOGICAL SIGNIFICANCE

Venoms represent an adaptive trait and an example of both divergent and convergent evolution. A deep understanding of the composition of venoms and of the principles governing the evolution of venomous systems is of applied importance for exploring the enormous potential of venoms as sources of chemical and pharmacological novelty but also to fight the consequences of snakebite envenomings. Key to this is the identification of evolutionary and ecological trends at different taxonomical levels. However, the evolution of venomous species and their venoms do not always follow the same course, and the identification of structural and functional convergences and divergences among venoms is often unpredictable by a phylogenetic hypothesis. Snake venomics is a proteomic-centered strategy to deconstruct the complex molecular phenotypes the venom proteomes. The proteomic profiles of venoms from sixteen out of the 22 venomous species within the Viperidae and Elapidae families found in Costa Rica have been completed so far. An integrative view of their venom composition, including the identification of geographic and ontogenic variations, is hereby presented. Venom proteomic profiles offer a fundamental platform to assess the detailed immunorecognition of individual venom components by therapeutic or experimental antivenoms. This aspect is reviewed in the companion paper. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Proteomic and biochemical analyses of short-tailed pit viper (Gloydius brevicaudus) venom: age-related variation and composition-activity correlation

We conducted an in-depth analysis of the proteomic and biochemical profiles of the venom of neonate and adult short-tailed pit vipers (Gloydius brevicaudus). Identified proteins were assigned to a few main toxin families. Disintegrin, phospholipase A2 (PLA2), serine proteinase, cysteine-rich secretory protein, C-type lectin-like protein, l-amino acid oxidase and snake venom metalloproteinase (SVMP) were detected in both venoms, while 5'-nucleotidase was detected only in the adult venom. SVMP was the predominant protein family in both venoms (neonate: 65.7%; adult: 64.4%), followed by PLA2 (neonate: 13.4%; adult: 25.0%). Antivenomic analysis revealed that commercial G. brevicaudus antivenom almost neutralized the chromatographic peaks with medium and high molecular masses in both venoms, but did not completely recognize peaks with low molecular mass. Toxicological and enzymatic activities show remarkable age-related variation in G. brevicaudus venom, probably resulting from variation in venom composition. Our data demonstrate age-related variation across venomics, antivenomics and biochemical profiles of G. brevicaudus venom, and have implications for the management of G. brevicaudus bites, including improving antivenom preparation by combining both venoms.

BIOLOGICAL SIGNIFICANCE

This study investigates the composition and biochemical activity of neonate and adult Gloydius brevicaudus venoms. We found remarkable age-related variation in venom biological activity, likely the result of variation in venom composition. Antivenomics analysis was used to explore difference in neonate and adult G. brevicaudus venoms. Our findings have implications for the diagnosis and clinical management of G. brevicaudus bites, and the design of venom mixtures that will increase the efficacy of commercial antivenom. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Snake venomics of Bothrops punctatus, a semiarboreal pitviper species from Antioquia, Colombia

Bothrops punctatus is an endangered, semi-arboreal pitviper species distributed in Panamá, Colombia, and Ecuador, whose venom is poorly characterized. In the present work, the protein composition of this venom was profiled using the ‘snake venomics’ analytical strategy. Decomplexation of the crude venom by RP-HPLC and SDS-PAGE, followed by tandem mass spectrometry of tryptic digests, showed that it consists of proteins assigned to at least nine snake toxin families. Metalloproteinases are predominant in this secretion (41.5% of the total proteins), followed by C-type lectin/lectin-like proteins (16.7%), bradykinin-potentiating peptides (10.7%), phospholipases A₂ (93%), serine proteinases (5.4%), disintegrins (38%), L-amino acid oxidases (3.1%), vascular endothelial growth factors (17%), and cysteine-rich secretory proteins (1.2%). Altogether, 6.6% of the proteins were not identified. In vitro, the venom exhibited proteolytic, phospholipase A₂, and L-amino acid oxidase activities, as well as angiotensin-converting enzyme (ACE)-inhibitory activity, in agreement with the obtained proteomic profile. Cytotoxic activity on murine C2C12 myoblasts was negative, suggesting that the majority of venom phospholipases A₂ likely belong to the acidic type, which often lack major toxic effects. The protein composition of B. punctatus venom shows a good correlation with toxic activities here and previously reported, and adds further data in support of the wide diversity of strategies that have evolved in snake venoms to subdue prey, as increasingly being revealed by proteomic analyses.

Quantitative high-throughput profiling of snake venom gland transcriptomes and proteomes (Ovophis okinavensis and Protobothrops flavoviridis)

Background

Advances in DNA sequencing and proteomics have facilitated quantitative comparisons of snake venom composition. Most studies have employed one approach or the other. Here, both Illumina cDNA sequencing and LC/MS were used to compare the transcriptomes and proteomes of two pit vipers, Protobothrops flavoviridis and Ovophis okinavensis, which differ greatly in their biology.

Results

Sequencing of venom gland cDNA produced 104,830 transcripts. The Protobothrops transcriptome contained transcripts for 103 venom-related proteins, while the Ovophis transcriptome contained 95. In both, transcript abundances spanned six orders of magnitude. Mass spectrometry identified peptides from 100% of transcripts that occurred at higher than contaminant (e.g. human keratin) levels, including a number of proteins never before sequenced from snakes. These transcriptomes reveal fundamentally different envenomation strategies. Adult Protobothrops venom promotes hemorrhage, hypotension, incoagulable blood, and prey digestion, consistent with mammalian predation. Ovophis venom composition is less readily interpreted, owing to insufficient pharmacological data for venom serine and metalloproteases, which comprise more than 97.3% of Ovophis transcripts, but only 38.0% of Protobothrops transcripts. Ovophis venom apparently represents a hybrid strategy optimized for frogs and small mammals.

Conclusions

This study illustrates the power of cDNA sequencing combined with MS profiling. The former quantifies transcript composition, allowing detection of novel proteins, but cannot indicate which proteins are actually secreted, as does MS. We show, for the first time, that transcript and peptide abundances are correlated. This means that MS can be used for quantitative, non-invasive venom profiling, which will be beneficial for studies of endangered species.

Proteomic and functional profiling of the venom of Bothrops ayerbei from Cauca, Colombia, reveals striking interspecific variation with Bothrops asper venom

Bothrops ayerbei, a pitviper inhabiting the Patía River's basin (Valle Alto del Río Patía) in the Southwestern Department of Cauca, Colombia, was considered as a variant form of Bothrops asper prior to being proposed as a new species in 2010, on the basis of subtle morphological differences. This study reports the proteomic and functional profiling of B. ayerbei venom. Its most striking feature is an almost complete absence (0.7%) of phospholipases A2 (PLA2), which is in contrast to the high proportion of these enzymes (25.3%) in the venom of B. asper from Cauca, as well as in other species of Bothrops. The predominant proteins in B. ayerbei venom are metalloproteinases (53.7%), in agreement with its higher hemorrhagic and lethal activities compared to B. asper venom. Moreover, the negligible content of PLA2s in B. ayerbei venom correlates with its weaker myotoxic effect, in contrast to B. asper venom, here shown to contain abundant Asp49- and Lys49-type PLA2s responsible for its strong myotoxic activity. Other components identified in B. ayerbei venom include bradykinin-potentiating-like peptides and proteins belonging to the C-type lectin/lectin-like, serine proteinase, l-amino acid oxidase, disintegrin, cysteine-rich secretory protein, nerve growth factor, and phosphodiesterase families. The venom composition of B. ayerbei resembles that of neonate specimens of B. asper, which shows a predominance of metalloproteinases, with only low amounts of PLA2s. Therefore, the present findings suggest that the expression of venom proteins in B. ayerbei, in contrast to B. asper, might retain a marked 'paedomorphic' condition. Altogether, the proteomic and toxicological characterization of the venom of B. ayerbei here reported argues in favor of its taxonomical separation from B. asper in Cauca, Colombia.

BIOLOGICAL SIGNIFICANCE

B. ayerbei, a pitviper found in Cauca, Colombia, had been considered as a variant form of B. asper, but was recently described as a new species on the basis of subtle morphological differences. Our study provides the first detailed proteomic and functional analysis of the venom of B. ayerbei, revealing striking interspecific variation from B. asper, thus arguing in favor of their taxonomical separation. In addition, the observed composition of the venom of B. ayerbei correlates well with its functional and toxicological properties, helping to predict the main clinical manifestations in envenomings by this species, which inflicts a considerable number of snakebites in the Southwestern regions of Colombia.

Comparison of phylogeny, venom composition and neutralization by antivenom in diverse species of bothrops complex

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

Snakebite envenomation is a serious health issue in Latin America, particularly in the Amazon, where antivenom administration may be delayed due to logistic constraints. Bothrops snakes are involved in most of the snakebite-related accidents in Brazil. This work reports a comparative study of the toxin composition and antigenicity of the Bothrops venoms used to prepare the commercial antivenom and its effectiveness against the venom from Bothrops atrox, a prevalent Amazon species that is not included in the pool. Our data show a lack of connection between Bothrops taxonomic identity and venom composition. We also show that different toxins display distinct reactivity with the tested antivenom. However, the antivenom reacted similarly with each class of toxin present in the venoms of the different snakes studied. Important evidence was the neutralization of the major toxic effects of B. atrox venom, not included in the mixture of antigens used to produce the antivenom. Based on the observed antigenicity of the distinct protein classes of toxins, we suggest that it is possible to obtain pan-specific and efficient Bothrops antivenoms via immunization with venoms from a few species of snakes that are representative of the protein composition of a large number of targeted species.

The genesis of an exceptionally lethal venom in the timber rattlesnake (Crotalus horridus) revealed through comparative venom-gland transcriptomics

Background

Snake venoms generally show sequence and quantitative variation within and between species, but some rattlesnakes have undergone exceptionally rapid, dramatic shifts in the composition, lethality, and pharmacological effects of their venoms. Such shifts have occurred within species, most notably in Mojave (Crotalus scutulatus), South American (C. durissus), and timber (C. horridus) rattlesnakes, resulting in some populations with extremely potent, neurotoxic venoms without the hemorrhagic effects typical of rattlesnake bites.

Results

To better understand the evolutionary changes that resulted in the potent venom of a population of C. horridus from northern Florida, we sequenced the venom-gland transcriptome of an animal from this population for comparison with the previously described transcriptome of the eastern diamondback rattlesnake (C. adamanteus), a congener with a more typical rattlesnake venom. Relative to the toxin transcription of C. adamanteus, which consisted primarily of snake-venom metalloproteinases, C-type lectins, snake-venom serine proteinases, and myotoxin-A, the toxin transcription of C. horridus was far simpler in composition and consisted almost entirely of snake-venom serine proteinases, phospholipases A₂, and bradykinin-potentiating and C-type natriuretic peptides. Crotalus horridus lacked significant expression of the hemorrhagic snake-venom metalloproteinases and C-type lectins. Evolution of shared toxin families involved differential expansion and loss of toxin clades within each species and pronounced differences in the highly expressed toxin paralogs. Toxin genes showed significantly higher rates of nonsynonymous substitution than nontoxin genes. The expression patterns of nontoxin genes were conserved between species, despite the vast differences in toxin expression.

Conclusions

Our results represent the first complete, sequence-based comparison between the venoms of closely related snake species and reveal in unprecedented detail the rapid evolution of snake venoms. We found that the difference in venom properties resulted from major changes in expression levels of toxin gene families, differential gene-family expansion and loss, changes in which paralogs within gene families were expressed at high levels, and higher nonsynonymous substitution rates in the toxin genes relative to nontoxins. These massive alterations in the genetics of the venom phenotype emphasize the evolutionary lability and flexibility of this ecologically critical trait.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Effect of monospecific antibodies against baltergin in myotoxicity induced by Bothrops alternatus venom from northeast of Argentina. Role of metalloproteinases in muscle damage

Myotoxicity, one of the most relevant local manifestations in envenomation by Bothrops genus, may result from a direct action of myotoxins or be due to an indirect vascular degeneration and ischemia. Baltergin, a snake venom metalloproteinase (SVMP), isolated from Bothrops alternatus venom has been used to obtain monospecific IgG, in order to determine the relative role of toxin in myotoxicity induced by whole venom. Bothrops diporus venom, another medical relevant genus of the northeastern region of Argentina, was also studied. Anti-baltergin IgG was able to neutralize completely the hemorrhagic activity of B. alternatus venom at an antibodies:venom ratio of 30:1 (w:w). However, mice injected with B. diporus venom showed a small spot remaining even at the highest ratio of IgG:venom assayed (50:1; w:w). Specific antibodies were efficient to neutralize the myotoxicity of B. alternatus venom at ratio 30:1 (w:w) but did not neutralize the same effects in B. diporus venom. Anti-baltergin polyclonal antibodies were useful tools for revealing the central role of SVMPs in the development of myotoxicity of B. alternatus venom, as well as, helping to suggest indirectly presence of potent myotoxic phospholipases A2 (PLA2s) in B. diporus venom.

In vitro antiplasmodial activity of phospholipases A2 and a phospholipase homologue isolated from the venom of the snake Bothrops asper

The antimicrobial and antiparasite activity of phospholipase A₂ (PLA₂) from snakes and bees has been extensively explored. We studied the antiplasmodial effect of the whole venom of the snake Bothrops asper and of two fractions purified by ion-exchange chromatography: one containing catalytically-active phospholipases A₂ (PLA₂) (fraction V) and another containing a PLA₂ homologue devoid of enzymatic activity (fraction VI). The antiplasmodial effect was assessed on in vitro cultures of Plasmodium falciparum. The whole venom of B. asper, as well as its fractions V and VI, were active against the parasite at 0.13 ± 0.01 µg/mL, 1.42 ± 0.56 µg/mL and 22.89 ± 1.22 µg/mL, respectively. Differences in the cytotoxic activity on peripheral blood mononuclear cells between the whole venom and fractions V and VI were observed, fraction V showing higher toxicity than total venom and fraction VI. Regarding toxicity in mice, the whole venom showed the highest lethal effect in comparison to fractions V and VI. These results suggest that B. asper PLA₂ and its homologue have antiplasmodial potential.

Snake venom metalloproteinases

Recent proteomic analyses of snake venoms show that metalloproteinases represent major components in most of the Crotalid and Viperid venoms. In this chapter we discuss the multiple activities of the SVMPs. In addition to hemorrhagic activity, members of the SVMP family also have fibrin(ogen)olytic activity, act as prothrombin activators, activate blood coagulation factor X, possess apoptotic activity, inhibit platelet aggregation, are pro-inflammatory and inactivate blood serine proteinase inhibitors. Clearly the SVMPs have multiple functions in addition to their well-known hemorrhagic activity. The realization that there are structural variations in the SVMPs and the early studies that led to their classification represents an important event in our understanding of the structural forms of the SVMPs. The SVMPs were subdivided into the P-I, P-II and P-III protein classes. The noticeable characteristic that distinguished the different classes was their size (molecular weight) differences and domain structure: Class I (P-I), the small SVMPs, have molecular masses of 20-30 kDa, contain only a pro domain and the proteinase domain; Class II (P-II), the medium size SVMPs, molecular masses of 30-60 kDa, contain the pro domain, proteinase domain and disintegrin domain; Class III (P-III), the large SVMPs, have molecular masses of 60-100 kDa, contain pro, proteinase, disintegrin-like and cysteine-rich domain structure. Another significant advance in the SVMP field was the characterization of the crystal structure of the first P-I class SVMP. The structures of other P-I SVMPs soon followed and the structures of P-III SVMPs have also been determined. The active site of the metalloproteinase domain has a consensus HEXXHXXGXXHD sequence and a Met-turn. The "Met-turn" structure contains a conserved Met residue that forms a hydrophobic basement for the three zinc-binding histidines in the consensus sequence.

Combined snake venomics and venom gland transcriptomic analysis of Bothropoides pauloensis

Unraveling the repertoire of venom toxins of Bothropoides pauloensis was assessed by snake venomics and venom gland transcriptomic surveys. Both approaches yielded converging overall figures, pointing to metalloproteinases (~37%), PLA(2)s (26-32%), and vasoactive (bradykinin-potentiating) peptides (12-17%) as the major toxin classes. The high occurrence of SVMPs, PLA(2) molecules, vasoactive peptides, along with serine proteinases, explains the local and systemic effects observed in envenomations by B. pauloensis. Minor (<3%) C-type lectin, serine proteinase, L-amino acid oxidase, nerve growth factor, and CRISP molecules were also identified in the transcriptome and the proteome. Low abundance (0.3%) EST singletons coding for vascular endothelial growth factor (svVEGF), ohanin, hyaluronidase, and 5' nucleotidase were found only in the venom gland cDNA library. At the molecular level, the transcriptomic and proteomic datasets display low compositional concordance. In particular, although there is good agreement between transcriptome and proteome in the identity of BPPs, PLA(2) molecules and L-amino acid oxidase, both datasets strongly depart in their C-type lectin and SVMP complements. These data support the view that venom composition is influenced by transcriptional and translational mechanisms and emphasize the value of combining proteomic and transcriptomic approaches to acquire a more complete understanding of the toxinological profile and natural history of the snake venom.

A transcriptomic view of the proteome variability of newborn and adult Bothrops jararaca snake venoms

Background

Snake bite is a neglected public health problem in communities in rural areas of several countries. Bothrops jararaca causes many snake bites in Brazil and previous studies have demonstrated that the pharmacological activities displayed by its venom undergo a significant ontogenetic shift. Similarly, the venom proteome of B. jararaca exhibits a considerable variation upon neonate to adult transition, which is associated with changes in diet from ectothermic prey in early life to endothermic prey in adulthood. Moreover, it has been shown that the Brazilian commercial antibothropic antivenom, which is produced by immunization with adult venom, is less effective in neutralizing newborn venom effects. On the other hand, venom gland transcripts of newborn snakes are poorly known since all transcriptomic studies have been carried out using mRNA from adult specimens.

Methods/Principal Findings

Here we analyzed venom gland cDNA libraries of newborn and adult B. jararaca in order to evaluate whether the variability demonstrated for its venom proteome and pharmacological activities was correlated with differences in the structure of toxin transcripts. The analysis revealed that the variability in B. jararaca venom gland transcriptomes is quantitative, as illustrated by the very high content of metalloproteinases in the newborn venom glands. Moreover, the variability is also characterized by the structural diversity of SVMP precursors found in newborn and adult transcriptomes. In the adult transcriptome, however, the content of metalloproteinase precursors considerably diminishes and the number of transcripts of serine proteinases, C-type lectins and bradykinin-potentiating peptides increase. Moreover, the comparison of the content of ESTs encoding toxins in adult male and female venom glands showed some gender-related differences.

Conclusions/Significance

We demonstrate a substantial shift in toxin transcripts upon snake development and a marked decrease in the metalloproteinase P-III/P-I class ratio which are correlated with changes in the venom proteome complexity and pharmacological activities.

Bothrops jararaca is one of the most abundant venomous snake species in Brazil. It is primarily a nocturnal and generalist animal, however, it exhibits a notable ontogenetic shift in diet, feeding mainly on arthropods, lizards, and amphibians (ectothermic prey) through its juvenile phase and on small mammals (endothermic animals) during adult life. Due to its broad geographical distribution, this species is responsible for the majority of the accidents by Bothrops genus in Brazil. Studies on envenomation cases with newborn and adult B. jararaca snakes have shown distinct patterns, mainly related to blood coagulation disorders, which seems to be prominent in accidents with newborn specimens. Moreover, it has been demonstrated that the Brazilian commercial antibothropic antivenom, which is produced by immunization with adult venom, is less effective in neutralizing newborn venom effects. In this study we analyzed the venom gland transcriptome of newborn snake specimens and compared the content of toxin transcripts with that of adult specimens. We demonstrate that upon B. jararaca development, its repertoire of mRNAs encoding toxins changes both qualitatively and quantitatively and these alterations are associated with the venom proteome profiles and pharmacological activities displayed by newborn and adult specimens.

Snake-venom resistance as a mammalian trophic adaptation: lessons from didelphid marsupials

Mammals that prey on venomous snakes include several opossums (Didelphidae), at least two hedgehogs (Erinaceidae), several mongooses (Herpestidae), several mustelids, and some skunks (Mephitidae). As a group, these taxa do not share any distinctive morphological traits. Instead, mammalian adaptations for ophiophagy seem to consist only in the ability to resist the toxic effects of snake venom. Molecular mechanisms of venom resistance (as indicated by biochemical research on opossums, mongooses, and hedgehogs) include toxin-neutralizing serum factors and adaptive changes in venom-targeted molecules. Of these, toxin-neutralizing serum factors have received the most research attention to date. All of the toxin-neutralizing serum proteins discovered so far in both opossums and mongooses are human α1B-glycoprotein homologs that inhibit either snake-venom metalloproteinases or phospholipase A(2) myotoxins. By contrast, adaptive changes in venom-targeted molecules have received far less attention. The best-documented examples include amino-acid substitutions in mongoose nicotinic acetylcholine receptor that inhibit binding by α-neurotoxins, and amino-acid substitutions in opossum von Willebrand factor (vWF) that are hypothesized to weaken the bond between vWF and coagulopathic C-type lectins. Although multiple mechanisms of venom resistance are known from some species, the proteomic complexity of most snake venoms suggests that the evolved biochemical defences of ophiophagous mammals are likely to be far more numerous than currently recognized. Whereas most previous research in this field has been motivated by the potential for medical applications, venom resistance in ophiophagous mammals is a complex adaptation that merits attention from comparative biologists. Unfortunately, evolutionary inference is currently limited by ignorance about many relevant facts that can only be provided by future research.

Improving antivenom availability and accessibility: science, technology, and beyond

Snakebite envenomings constitute a serious and neglected public health problem. Despite the fact that effective treatment exists, i.e. administration of animal-derived antivenoms, the availability and accessibility of these life-saving immunobiologicals is deficitary in various parts of the world, particularly in sub-Saharan Africa and some regions of Asia. This article discusses some of the problems that need to be circumvented in order to improve the availability and accessibility of antivenoms. The conglomerate of antivenom manufacturers is highly heterogeneous in terms of technological base, qualification of staff, implementation of Good Manufacturing Practices (GMPs), and volume of production. Therefore, improvements in antivenom quality and availability should be based on strategies tailored to the situation of each region or country; in this context, three different scenarios are discussed. Accessibility of antivenoms demands concerted efforts at multiple levels, including raising the awareness of public health authorities on the relevance of the problem, implementing innovative antivenom purchasing schemes, strengthening national distribution channels on the basis of robust epidemiological information, improving the cold chain and the provision of health services in remote rural settings, supporting the correct use of antivenoms, and promoting the involvement of local community organizations in various aspects of prevention and management. These tasks should be envisaged in terms of synergistic, interprogrammatic and intersectorial interventions, with the participation of many players.

Exploring the proteomes of the venoms of the Peruvian pit vipers Bothrops atrox, B. barnetti and B. pictus

We report the comparative proteomic characterization of the venoms of Bothrops atrox, B. barnetti and B. pictus. The venoms were subjected to RP-HPLC and the resulting fractions analyzed by SDS-PAGE. The proteins were cut from the gels, digested with trypsin and identified via peptide mass fingerprint and manual sequencing of selected peptides by MALDI-TOF/TOF mass spectrometry. Around 20-25 proteins were identified belonging to only 6-7 protein families. Metalloproteinases of the classes P-I and P-III were the most abundant proteins in all venoms (58-74% based on peak area A214 nm), followed by phospholipases-A(2) (6.4-14%), disintegrins (3.2-9%) and serine proteinases (7-11%), and some of these proteins occurred in several isoforms. In contrast cysteine-rich secretory proteins and L-amino acid oxidases appeared only as single isoforms and were found only in B. atrox and B. barnetti. C-type lectins were also detected in all venoms but at low levels (~ 5%). Furthermore, the venoms contain variable numbers of peptides (<3 kDa) and non-protein compounds which were not considered in this work. The protein composition of the investigated Bothrops species is in agreement with their pharmacological and pathological effects.