Bothrops jararaca venom proteome rearrangement upon neonate to adult transition

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.

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.

Acute Lonomia obliqua caterpillar envenomation-induced physiopathological alterations in rats: evidence of new toxic venom activities and the efficacy of serum therapy to counteract systemic tissue damage

The clinical manifestations of Lonomia obliqua caterpillar envenomation are systemic hemorrhage and acute kidney injury. In an effort to better understand the physiopathological mechanisms of envenomation, a rat model was established to study systemic tissue damage during L. obliqua envenomation. An array of acute venom effects was characterized, including biochemical, hematological, histopathological, myotoxic and genotoxic alterations. Rapid increases in serum alanine and aspartate transaminases, γ-glutamyl transferase, lactate dehydrogenase, hemoglobin, bilirubin, creatinine, urea and uric acid were observed, indicating that intravascular hemolysis and liver and kidney damage had occurred. Treatment with a specific antivenom (antilonomic serum) for up to 2 h post-venom injection neutralized the biochemical alterations. However, treatment after 6 h post-venom injection failed to normalize all biochemical parameters, despite its efficacy in reversing coagulation dysfunction. The hematological findings were consistent with hemolytic anemia and neutrophilic leukocytosis. The histopathological alterations were mainly related to hemorrhage and inflammation in the subcutaneous tissue, lung, heart and kidneys. Signs of congestion and hemosiderosis were evident in the spleen, and hemoglobin and/or myoglobin casts were also detected in the renal tubules. Increased levels of creatine kinase and creatine kinase-MB were correlated with the myocardial necrosis observed in vivo and confirmed the myotoxicity detected in vitro in isolated extensor digitorum longus muscles. Significant DNA damage was observed in the kidneys, heart, lung, liver and lymphocytes. The majority of the DNA lesions in the kidney were due to oxidative damage. The results presented here will aid in understanding the pathology underlying Lonomia's envenomation.

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.

Jararhagin, a hemorrhagic snake venom metalloproteinase from Bothrops jararaca

Jararhagin is a metalloproteinase isolated from Bothrops jararaca snake venom, which has been extensively studied. These studies showed its involvement on most of the systemic and local damaging effects of snakebite envenomings. In this review we comment on the major targets of jararhagin as the vascular endothelium, platelets and coagulation factors and also its action on other cell systems as inflammatory cells and their mediators, cancer and cell signaling. The mechanisms of jararhagin action are discussed together with structural features essential for the expression of its biological activities. The studies reviewed here denote jararhagin as a prototype for studies of snake venom metalloproteinases, bringing new insights into cellular-matrix interactions and adding for the improvement of snakebite treatment.

Peptidomics of three Bothrops snake venoms: insights into the molecular diversification of proteomes and peptidomes

Snake venom proteomes/peptidomes are highly complex and maintenance of their integrity within the gland lumen is crucial for the expression of toxin activities. There has been considerable progress in the field of venom proteomics, however, peptidomics does not progress as fast, because of the lack of comprehensive venom sequence databases for analysis of MS data. Therefore, in many cases venom peptides have to be sequenced manually by MS/MS analysis or Edman degradation. This is critical for rare snake species, as is the case of Bothrops cotiara (BC) and B. fonsecai (BF), which are regarded as near threatened with extinction. In this study we conducted a comprehensive analysis of the venom peptidomes of BC, BF, and B. jararaca (BJ) using a combination of solid-phase extraction and reversed-phase HPLC to fractionate the peptides, followed by nano-liquid chromatography-tandem MS (LC-MS/MS) or direct infusion electrospray ionization-(ESI)-MS/MS or MALDI-MS/MS analyses. We detected marked differences in the venom peptidomes and identified peptides ranging from 7 to 39 residues in length by de novo sequencing. Forty-four unique sequences were manually identified, out of which 30 are new peptides, including 17 bradykinin-potentiating peptides, three poly-histidine-poly-glycine peptides and interestingly, 10 L-amino acid oxidase fragments. Some of the new bradykinin-potentiating peptides display significant bradykinin potentiating activity. Automated database search revealed fragments from several toxins in the peptidomes, mainly from l-amino acid oxidase, and allowed the determination of the peptide bond specificity of proteinases and amino acid occurrences for the P4-P4' sites. We also demonstrate that the venom lyophilization/resolubilization process greatly increases the complexity of the peptidome because of the imbalance caused to the venom proteome and the consequent activity of proteinases on venom components. The use of proteinase inhibitors clearly showed different outcomes in the peptidome characterization and suggested that degradomic-peptidomic analysis of snake venoms is highly sensitive to the conditions of sampling procedures.

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.

Intraspecific variation of biological activities in venoms from wild and captive Bothrops jararaca

The venom of Bothrops jararaca is composed of complex mixture of molecules, mainly lectins, metalloproteinases, serinoproteinases, desintegrins, phospholipases, and peptides. This composition may vary according to the snake's age, gender, and region of origin. The aim of the was to determine individual variation in Bothrops jararaca venom in the Botucatu region, Sao Paulo State, Brazil, by means of enzymatic, biochemical, and pharmacological characterization, utilizing in vitro tests and biological assays. The activities were compared with those of Brazilian Reference Venom (BRV). Protein concentration varied between adult and juvenile groups. The electrophoretic profiles were similar, with molecular masses ranging between 25 and 50 kD, but with intraspecific variations. Reverse-phase high-performance liquid chromatography (RP-HPLC) revealed protein concentration differences. Coagulant activity did not differ significantly among adult groups, but there was a large variation between juvenile venom and BRV, which coagulated more extensively. Venoms from adults displayed greater hemorrhagic activity, especially in males recently obtained from the wild. In contrast, juveniles kept in captivity and adult males showed higher values. Edematogenic activity displayed an increase in edema in all groups. At the mean lethal dose (LD₅₀), toxicity varied significantly between groups, with venom from captive females being threefold more toxic than juvenile venom. Data illustrate the intra- and interspecific complexity that occurs in snake venoms, which may be attributed to ontogenetic, sexual, and environmental factors that affect variability in Bothrops jararaca venom. Further, it is proposed that Brazilian public health authorities document the constitution of pooled venom employed in the immunization of serum-producing animals due to this variability in venom properties. Given the large Brazilian territory, this variability requires regional monitoring and evaluation of the efficacy of bothropic antivenom in treatment of snakebite and consequent permanent sequelae observed.