Molecular mechanisms underlying intraspecific variation in snake venom.

Elucidating the molecular mechanisms underlying snake venom variability provides important clues for understanding how the biological functions of this powerful toxic arsenal evolve. We analyzed in detail individual transcripts and venom protein isoforms produced by five specimens of a venomous snake (Bothrops atrox) from two nearby but genetically distinct populations from the Brazilian Amazon rainforest which show functional similarities in venom properties. Individual variation was observed among the venoms of these specimens, but the overall abundance of each general toxin family was conserved both in transcript and in venom protein levels. However, when expression of independent paralogues was analyzed, remarkable differences were observed within and among each toxin group, both between individuals and between populations. Transcripts for functionally essential venom proteins ("core function" proteins) were highly expressed in all specimens and showed similar transcription/translation rates. In contrast, other paralogues ("adaptive" proteins) showed lower expression levels and the toxins they coded for varied among different individuals. These results provide support for the inferences that (a) expression and translational differences play a greater role in defining adaptive variation in venom phenotypes than does sequence variation in protein coding genes and (b) convergent adaptive venom phenotypes can be generated through different molecular mechanisms. SIGNIFICANCE: Analysis of individual transcripts and venom protein isoforms produced by specimens of a venomous snake (Bothrops atrox), from the Brazilian Amazon rainforest, revealed that transcriptional and translational mechanisms contribute to venom phenotypic variation. Our finding of evidence for high expression of toxin proteins with conserved function supports the hypothesis that the venom phenotype consists of two kinds of proteins: conserved "core function" proteins that provide essential functional activities with broader relevance and less conserved "adaptive" proteins that vary in expression and may permit customization of protein function. These observations allowed us to suggest that genetic mechanisms controlling venom variability are not restricted to selection of gene copies or mutations in structural genes but also to selection of the mechanisms controlling gene expression, contributing to the plasticity of this important phenotype for venomous snakes.

Functional proteomic analyses of Bothrops atrox venom reveals phenotypes associated with habitat variation in the Amazon.

Venom variability is commonly reported for venomous snakes including Bothrops atrox. Here, we compared the composition of venoms from B. atrox snakes collected at Amazonian conserved habitats (terra-firme upland forest and várzea) and human modified areas (pasture and degraded areas). Venom samples were submitted to shotgun proteomic analysis as a whole or compared after fractionation by reversed-phase chromatography. Whole venom proteomes revealed a similar composition among the venoms with predominance of SVMPs, CTLs, and SVSPs and intermediate amounts of PLA2s and LAAOs. However, when distribution of particular isoforms was analyzed by either method, the venom from várzea snakes showed a decrease in hemorrhagic SVMPs and an increase in SVSPs, and procoagulant SVMPs and PLA2s. These differences were validated by experimental approaches including both enzymatic and in vivo assays, and indicated restrictions in respect to antivenom efficacy to variable components. Thus, proteomic analysis at the isoform level combined to in silico prediction of functional properties may indicate venom biological activity. These results also suggest that the prevalence of functionally distinct isoforms contributes to the variability of the venoms and could reflect the adaptation of B. atrox to distinct prey communities in different Amazon habitats. BIOLOGICAL SIGNIFICANCE: In this report, we compared isoforms present in venoms from snakes collected at different Amazonian habitats. By means of a species venom gland transcriptome and the in silico functional prediction of each isoform, we were able to predict the principal venom activities in vitro and in animal models. We also showed remarkable differences in the venom pools from snakes collected at the floodplain (várzea habitat) compared to other habitats. Not only was this venom less hemorrhagic and more procoagulant, when compared to the venom pools from the other three habitats studied, but also this enhanced procoagulant activity was not efficiently neutralized by Bothrops antivenom. Thus, using a functional proteomic approach, we highlighted intraspecific differences in B. atrox venom that could impact both in the ecology of snakes but also in the treatment of snake bite patients in the region.

An in-depth snake venom proteopeptidome characterization: Benchmarking Bothrops jararaca.

A large-scale proteomic approach was devised to advance the understanding of venom composition. Bothrops jararaca venom was fractionated by OFFGEL followed by chromatography, generating peptidic and proteic fractions. The latter was submitted to trypsin digestion. Both fractions were separately analyzed by reversed-phase nanochromatography coupled to high resolution mass spectrometry. This strategy allowed deeper and joint characterizations of the peptidome and proteome (proteopeptidome) of this venom. Our results lead to the identification of 46 protein classes (with several uniquely assigned proteins per class) comprising eight high-abundance bona fide venom components, and 38 additional classes in smaller quantities. This last category included previously described B. jararaca venom proteins, common Elapidae venom constituents (cobra venom factor and three-finger toxin), and proteins typically encountered in lysosomes, cellular membranes and blood plasma. Furthermore, this report is the most complete snake venom peptidome described so far, both in number of peptides and in variety of unique proteins that could have originated them. It is hypothesized that such diversity could enclose cryptides, whose bioactivities would contribute to envenomation in yet undetermined ways. Finally, we propose that the broad range screening of B. jararaca peptidome will facilitate the discovery of bioactive molecules, eventually leading to valuable therapeutical agents. BIOLOGICAL SIGNIFICANCE: Our proteopeptidomic strategy yielded unprecedented insights into the remarkable diversity of B. jararaca venom composition, both at the peptide and protein levels. These results bring a substantial contribution to the actual pursuit of large-scale protein-level assignment in snake venomics. The detection of typical elapidic venom components, in a Viperidae venom, reinforces our view that the use of this approach (hand-in-hand with transcriptomic and genomic data) for venom proteomic analysis, at the specimen-level, can greatly contribute for venom toxin evolution studies. Furthermore, data were generated in support of a previous hypothesis that venom gland secretory vesicles are specialized forms of lysosomes. Two testable hypotheses also emerge from the results of this work. The first is that a nucleobindin-2-derived protein could lead to prey disorientation during envenomation, aiding in its capture by the snake. The other being that the venom's peptidome might contain a population of cryptides, whose biological activities could lead to the development of new therapeutical agents.

Viperid Envenomation Wound Exudate Contributes to Increased Vascular Permeability via a DAMPs/TLR-4 Mediated Pathway.

Viperid snakebite envenomation is characterized by inflammatory events including increase in vascular permeability. A copious exudate is generated in tissue injected with venom, whose proteomics analysis has provided insights into the mechanisms of venom-induced tissue damage. Hereby it is reported that wound exudate itself has the ability to induce increase in vascular permeability in the skin of mice. Proteomics analysis of exudate revealed the presence of cytokines and chemokines, together with abundant damage associated molecular pattern molecules (DAMPs) resulting from both proteolysis of extracellular matrix and cellular lysis. Moreover, significant differences in the amounts of cytokines/chemokines and DAMPs were detected between exudates collected 1 h and 24 h after envenomation, thus highlighting a complex temporal dynamic in the composition of exudate. Pretreatment of mice with Eritoran, an antagonist of Toll-like receptor 4 (TLR4), significantly reduced the exudate-induced increase in vascular permeability, thus suggesting that DAMPs might be acting through this receptor. It is hypothesized that an "Envenomation-induced DAMPs cycle of tissue damage" may be operating in viperid snakebite envenomation through which venom-induced tissue damage generates a variety of DAMPs which may further expand tissue alterations.

Processing of Snake Venom Metalloproteinases: Generation of Toxin Diversity and Enzyme Inactivation.

Snake venom metalloproteinases (SVMPs) are abundant in the venoms of vipers and rattlesnakes, playing important roles for the snake adaptation to different environments, and are related to most of the pathological effects of these venoms in human victims. The effectiveness of SVMPs is greatly due to their functional diversity, targeting important physiological proteins or receptors in different tissues and in the coagulation system. Functional diversity is often related to the genetic diversification of the snake venom. In this review, we discuss some published evidence that posit that processing and post-translational modifications are great contributors for the generation of functional diversity and for maintaining latency or inactivation of enzymes belonging to this relevant family of venom toxins.

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