Individual Variability in the Venom Proteome of Juvenile Bothrops jararaca Specimens

A snapshot of Bothrops jararaca snake venom gland subcellular proteome

Snake venom protein synthesis undergoes finely regulated processes in the specialized secretory epithelium within the venom gland. Such processes occur within a defined period in the cell and at specific cellular locations. Thus, the determination of subcellular proteomes allows the characterization of protein groups for which the site may be relevant to their biological roles, thereby allowing the deconvolution of complex biological circuits into functional information. In this regard, we performed subcellular fractionation of proteins from B. jararaca venom gland, focusing on nuclear proteins since this cellular compartment comprises key effectors that shape gene expression. Our results provided a snapshot of B. jararaca's subcellular venom gland proteome and pointed to a 'conserved' proteome core among different life stages (newborn and adult) and between sexes (adult male and female). Overall, the top 15 highly abundant proteins identified in B. jararaca venom glands mirrored the panel of highly expressed genes in human salivary glands. Therefore, the expression profile observed for such a protein set could be considered a conserved core signature of salivary gland secretory epithelium. Moreover, the newborn venom gland displayed a unique expression signature of transcription factors involved in regulating transcription and biosynthetic processes and may mirror biological constraints of the ontogenetic development of B. jararaca, contributing to venom proteome diversity.

Who are the most affected by Bothrops snakebite envenoming in Brazil? A Clinical-epidemiological profile study among the regions of the country

Snakebite envenoming represents an important Neglected Tropical Disease (NTD) that mainly affects tropical and subtropical developing countries according to the World Health Organization (WHO). As a priority issue in the tropics, it is estimated that accidental encounter between snakes and humans is the leading cause of morbidity and mortality among all NTDs in the world. In Brazil, an extremely diverse country with continental dimensions, snakebite envenoming is the second leading cause of reported human envenoming. Treating the disease has been an unprecedented challenge for Brazilian Health Systems for decades. Despite access to Antivenom therapy and distributing it free of charge across the country, Brazil faces numerous issues regarding the notification process and accurate treatment targeting for at-risk populations. Thus, this study aimed to identify the temporal epidemiological dynamics of accidents caused by Bothrops snakes in Brazil, the country's major group of venomous snakes, based on secondary information from the online database provided by The Brazilian Notifiable Diseases Information System (SINAN). For this purpose, reported Bothrops snakebites between 2012 and 2021 were counted, then the data were analyzed. We looked at the frequency, occurrence, mortality rates, case fatality rate (CFR), age and gender distribution, and the time lapse between the incident and the initiation of Antivenom therapy. The data were also organized considering regional variations of the country. Throughout the studied period, a total of 202,604 cases of envenoming caused by Bothrops spp. were notified, resulting in 766 fatalities. These accidents were found to occur in variable proportions across different regions in Brazil, with notable concentrations observed in the North, Northeast, and Southeast regions. The epidemiological profile of patients varied greatly between the regions, revealing that snake envenoming is much more a social, economic, and ecological problem than a medical one. In conclusion, our study provides an overview of the clinical and epidemiological profile of envenoming by Bothrops snakes in Brazil. Notably, this is the first study to present such information in a country as vast and diverse as Brazil, encompassing a comparative analysis of its regions using SINAN data, that proves to be a very useful national tool to improve the control and management of envenoming.

Mouse skin peptidomic analysis of the hemorrhage induced by a snake venom metalloprotease

Hemorrhage induced by snake venom metalloproteases (SVMPs) results from proteolysis, capillary disruption, and blood extravasation. HF3, a potent SVMP of Bothrops jararaca, induces hemorrhage at pmol doses in the mouse skin. To gain insight into the hemorrhagic process, the main goal of this study was to analyze changes in the skin peptidome generated by injection of HF3, using approaches of mass spectrometry-based untargeted peptidomics. The results revealed that the sets of peptides found in the control and HF3-treated skin samples were distinct and derived from the cleavage of different proteins. Peptide bond cleavage site identification in the HF3-treated skin showed compatibility with trypsin-like serine proteases and cathepsins, suggesting the activation of host proteinases. Acetylated peptides, which originated from the cleavage at positions in the N-terminal region of proteins in both samples, were identified for the first time in the mouse skin peptidome. The number of peptides acetylated at the residue after the first Met residue, mostly Ser and Ala, was higher than that of peptides acetylated at the initial Met. Proteins cleaved in the hemorrhagic skin participate in cholesterol metabolism, PPAR signaling, and in the complement and coagulation cascades, indicating the impairment of these biological processes. The peptidomic analysis also indicated the emergence of peptides with potential biological activities, including pheromone, cell penetrating, quorum sensing, defense, and cell-cell communication in the mouse skin. Interestingly, peptides generated in the hemorrhagic skin promoted the inhibition of collagen-induced platelet aggregation and could act synergistically in the local tissue damage induced by HF3.

Sialic acid-containing glycans play a role in the activity of snake venom proteases

Structural variability is a feature of snake venom proteins, and glycosylation is a post-translational modification that contributes to the diversification of venom proteomes. Studies by our group have shown that Bothrops venoms are distinctly defined by their glycoprotein content, and that most hybrid/complex N-glycans identified in these venoms contain sialic acid. Considering that metalloproteases and serine proteases are abundant components of Bothrops venoms and essential in the envenomation process, and that these enzymes contain several glycosylation sites, the role of sialic acid in venom proteolytic activity was evaluated. Here we show that removal of sialic acid by treatment of nine Bothrops venoms with neuraminidase (i) altered the pattern of gelatinolysis in zymography of most venoms and reduced the gelatinolytic activity of all venoms, (ii) decreased the proteolytic activity of some venoms on fibrinogen and the clotting activity of human plasma of all venoms, and (iii) altered the proteolysis profile of plasma proteins by B. jararaca venom, suggesting that sialic acid may play a role in the interaction of proteases with their protein substrates. In contrast, the profile of venom amidolytic activity on Bz-Arg-pNA did not change after removal of sialic acid, indicating that this monosaccharide is not essential in N-glycans of serine proteases acting on small substrates. In summary, these results expand the knowledge about the variability of the subproteomes of Bothrops venom proteases, and for the first time point to the importance of carbohydrate chains containing sialic acid in the enzymatic activities of venom proteases relevant in human envenomation.

Extracellular Vesicles from Bothrops jararaca Venom Are Diverse in Structure and Protein Composition and Interact with Mammalian Cells

Snake venoms are complex cocktails of non-toxic and toxic molecules that work synergistically for the envenoming outcome. Alongside the immediate consequences, chronic manifestations and long-term sequelae can occur. Recently, extracellular vesicles (EVs) were found in snake venom. EVs mediate cellular communication through long distances, delivering proteins and nucleic acids that modulate the recipient cell's function. However, the biological roles of snake venom EVs, including possible cross-organism communication, are still unknown. This knowledge may expand the understanding of envenoming mechanisms. In the present study, we isolated and characterized the EVs from Bothrops jararaca venom (Bj-EVs), giving insights into their biological roles. Fresh venom was submitted to differential centrifugation, resulting in two EV populations with typical morphology and size range. Several conserved EV markers and a subset of venom related EV markers, represented mainly by processing enzymes, were identified by proteomic analysis. The most abundant protein family observed in Bj-EVs was 5'-nucleotidase, known to be immunosuppressive and a low abundant and ubiquitous toxin in snake venoms. Additionally, we demonstrated that mammalian cells efficiently internalize Bj-EVs. The commercial antibothropic antivenom partially recognizes Bj-EVs and inhibits cellular EV uptake. Based on the proteomic results and the in vitro interaction assays using macrophages and muscle cells, we propose that Bj-EVs may be involved not only in venom production and processing but also in host immune modulation and long-term effects of envenoming.

From Snake Venoms to Therapeutics: A Focus on Natriuretic Peptides

Snake venom is a cocktail of multifunctional biomolecules that has evolved with the purpose of capturing prey and for defense. These biomolecules are classified into different classes based on their functions. They include three-finger toxins, natriuretic peptides, phospholipases and metalloproteinases. The focus for this review is on the natriuretic peptide (NP), which is an active component that can be isolated from the venoms of vipers and mambas. In these venoms, NPs contribute to the lowering of blood pressure, causing a rapid loss of consciousness in the prey such that its mobility is reduced, paralyzing the prey, and often death follows. Over the past 30 years since the discovery of the first NP in the venom of the green mamba, venom NPs have shown potential in the development of drug therapy for heart failure. Venom NPs have long half-lives, different pharmacological profiles, and may also possess different functions in comparison to the mammalian NPs. Understanding their mechanisms of action provides the strategies needed to develop new NPs for treatment of heart failure. This review summarizes the venom NPs that have been identified over the years and how they can be useful in drug development.

Profilings of subproteomes of lectin-binding proteins of nine Bothrops venoms reveal variability driven by different glycan types

Snake venom proteomes have long been investigated to explore a multitude of biologically active components that are used for prey capture and defense, and are involved in the pathological effects observed upon mammalian envenomation. Glycosylation is a major protein post-translational modification in venoms and contributes to the diversification of proteomes. We have shown that Bothrops venoms are markedly defined by their content of glycoproteins, and that most N-glycan structures of eight Bothrops venoms contain sialic acid, while bisected N-acetylglucosamine was identified in Bothrops cotiara venom. To further investigate the mechanisms involved in the generation of different venoms by related snakes, here the glycoproteomes of nine Bothrops venoms (Bothrops atrox, B. cotiara, Bothrops erythromelas, Bothrops fonsecai, B. insularis, Bothrops jararaca, Bothrops jararacussu, Bothrops moojeni and Bothrops neuwiedi) were comparatively analyzed by enrichment with three lectins of different specificities, recognizing bisecting N-acetylglucosamine- and sialic acid-containing glycoproteins, and mass spectrometry. The lectin capture strategy generated venom fractions enriched with several glycoproteins, including metalloprotease, serine protease, and L- amino acid oxidase, in addition to various types of low abundant enzymes. The different contents of lectin-enriched proteins underscore novel aspects of the variability of the glycoprotein subproteomes of Bothrops venoms and point to the role of distinct types of glycan chains in generating different venoms by closely related snake species.

Venom proteomic analysis of medically important Nigerian viper Echis ocellatus and Bitis arietans snake species

Snakebite envenoming remains a neglected tropical disease which poses severe health hazard, especially for the rural inhabitants in Africa. In Nigeria, vipers are responsible for the highest number of deaths. Hydrophilic interaction liquid chromatography coupled with LC-MS/MS was used to analyze the crude venoms of Echis ocellatus (Carpet viper) and Bitis arietans (Puff adder) in order to understand their venom proteomic identities. Results obtained revealed that gel-free proteomic analysis of the crude venoms led to the identification of 85 and 79 proteins, respectively. Seventy-eight (78) proteins were common between the two snake species with a 91.8% similarity score. The identified proteins belong to 18 protein families in E. ocellatus and 14 protein families in B. arietans. Serine proteases (22.31%) and metalloproteinases (21.06%) were the dominant proteins in the venom of B. arietans; while metalloproteinases (34.84%), phospholipase A₂s (21.19%) and serine proteases (15.50%) represent the major toxins in the E. ocellatus venom. Other protein families such as three-finger toxins and cysteine-rich venom proteins were detected in low proportions. This study provides an insight into the venom proteomic analysis of the two Nigerian viper species, which could be useful in identifying the toxin families to be neutralized in case of envenomation.

•

Venom proteomic of Nigeria's most medically important snakes is presented.

•

SVMP, SVSP and PLA₂ were the major toxin families in E. ocellatus and B. arietans.

•

The venom proteomes of these vipers displayed 91.8% similarity in composition.

From birth to adulthood: An analysis of the Brazilian lancehead (Bothrops moojeni) venom at different life stages

The Brazilian lancehead (Bothrops moojeni) has a wide distribution in Brazil and represents a serious public health hazard. Previous works reported that the symptoms of snakebites caused by B. moojeni juveniles’ bites were mainly related to coagulation, while those caused by adults’ bites had a more prominent local damage. In this work, we analyzed the venoms of B. moojeni at different life stages to better understand the ontogeny shift in this species. Snakes were grouped by age and sex, and venom pools were formed accordingly. Compositional analyses by one-dimensional electrophoresis (1-DE), chromatography, and mass spectrometry revealed that ontogenetic changes might be mostly related to phospholipase A₂ (PLA₂) and metalloproteases. Regarding the venoms functional aspect, proteolytic, L-amino acid oxidase, PLA₂, and coagulant in vitro activities were assayed, but only the first and the last ones showed age-related changes, with the venom of snakes up to 1 year-old displaying lower proteolytic and higher coagulant activities, while those from 2 years-old onward presented the opposite relation. The venoms of 3 years-old snakes were exceptions to the compositional and functional pattern of adults as both venoms presented profiles similar to neonates. Sex-related differences were observed in specific groups and were not age-related. In vivo experiments (median lethal dose and hemorrhagic activity) were statistically similar between neonates and adults, however we verified that the adult venom killed mice faster comparing to the neonates. All venoms were mostly recognized by the antibothropic serum and displayed similar profiles to 1-DE in western blotting. In conclusion, the Brazilian lancehead venom showed ontogenetic shift in its composition and activities. Furthermore, this change occurred in snakes from 1 to 2 years-old, and interestingly the venom pools from 3 years-old snakes had particular characteristics, which highlights the importance of comprehensive studies to better understand venom variability.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Proteomic analysis of three medically important Nigerian Naja (Naja haje, Naja katiensis and Naja nigricollis) snake venoms

Proteomics technologies enable a comprehensive study of complex proteins and their functions. The venom proteomes of three medically important Nigerian Elapidae snakes Naja haje, Naja katiensis and Naja nigricollis was studied using HILIC coupled with LC-MS/MS analysis. Results revealed a total of 57, 55, and 46 proteins in the venoms of N. haje, N. katiensis, and N. nigricollis, respectively, with molecular mass ranging between 5 and 185  kDa. These snakes have 38 common proteins in addition to 3 uncommon proteins: actiflagelin, cathelicidin, and cystatin identified in their venoms. The identified proteins belonged to 14 protein families in N. haje and N. katiensis, and 12 protein families in N. nigricollis. Of the total venom proteins, 3FTx was the most abundant protein family, constituting 52% in N. haje and N. katiensis, and 41% in N. nigricollis, followed by PLA₂, constituting 37% in N. nigricollis, 26% in N. haje, and 24% in N. katiensis. Other protein families, including LAAO, CRISPs, VEGF, PLB, CVF, SVMP, SVH, AMP, PI, Globin, Actin, and C-type lectins, were also detected, although, at very low abundances. Quantification of the relative abundance of each protein revealed that alpha and beta fibrinogenase and PLA_2, which constituted 18-26% of the total proteome, were the most abundant. The 3 uncommon proteins have no known function in snake venom. However, actiflagelin activates sperm motility; cystatin inhibits angiogenesis, while cathelicidin exerts antimicrobial effects. The three Nigerian Naja genus proteomes displayed 70% similarity in composition, which suggests the possibility of formulating antivenom that may cross-neutralise the venoms of cobra species found in Nigeria. These data provide insights into clinically relevant peptides/proteins present in the venoms of these snakes. Data are available via ProteomeXchange with identifier PXD024627.

The distinct N-terminomes of Bothrops jararaca newborn and adult venoms

Using approaches of transcriptomics and proteomics we have shown that the phenotype of Bothrops jararaca venom undergoes a significant rearrangement upon neonate to adult transition. Most regulatory processes in biology are intrinsically related to modifications of protein structure, function, and abundance. However, it is unclear to which extent intrinsic proteolysis affects toxins and snake venom phenotypes upon ontogenesis. Here we assessed the natural N-terminome of Bothrops jararaca newborn and adult venoms and explored the degree of N-terminal protein truncation in ontogenetic-based proteome variation. To this end we applied the Terminal Amine Isotopic Labeling of Substrates (TAILS) technology to characterize venom collected in the presence of proteinase inhibitors. We identified natural N-terminal sequences in the newborn (71) and adult (84) venoms, from which only 37 were common to both. However, truncated toxins were found in higher number in the newborn (212) than in the adult (140) venom. Moreover, sequences N-terminally blocked by pyroglutamic acid were identified in the newborn (55) and adult (49) venoms. Most toxin classes identified by their natural N-terminal sequences showed a similar number of unique peptides in the newborn and adult venoms, however, those of serine proteinases and C-type lectins were more abundant in the adult venom. Truncated sequences from at least ten toxin classes were detected, however the catalytic and cysteine-rich domains of metalloproteinases were the most prone to proteolysis, mainly in the newborn venom. Our results underscore the pervasiveness of truncations in most toxin classes and highlight variable post-translational events in newborn and adult venoms.

Venom Profiling of the Insular Species Bothrops alcatraz: Characterization of Proteome, Glycoproteome, and N-Terminome Using Terminal Amine Isotopic Labeling of Substrates

Bothrops alcatraz, a species endemic to Alcatrazes Islands, is regarded as critically endangered due to its small area of occurrence and the declining quality of its habitat. We recently reported the identification of N-glycans attached to toxins of Bothrops species, showing similar compositions in venoms of the B. jararaca complex (B. jararaca, B. insularis, and B. alcatraz). Here, we characterized B. alcatraz venom using electrophoretic, proteomic, and glycoproteomic approaches. Electrophoresis showed that B. alcatraz venom differs from B. jararaca and B. insularis; however, N-glycan removal revealed similarities between them, indicating that the occupation of N-glycosylation sites contributes to interspecies variability in the B. jararaca complex. Metalloproteinase was the major toxin class identified in the B. alcatraz venom proteome followed by serine proteinase and C-type lectin, and overall, the adult B. alcatraz venom resembles that of B. jararaca juvenile specimens. The comparative glycoproteomic analysis of B. alcatraz venom with B. jararaca and B. insularis indicated that there may be differences in the utilization of N-glycosylation motifs among their different toxin classes. Furthermore, we prospected for the first time the N-terminome of a snake venom using the terminal amine isotopic labeling of substrates (TAILS) approach and report the presence of ∼30% of N-termini corresponding to truncated toxin forms and ∼37% N-terminal sequences blocked by pyroglutamic acid in B. alcatraz venom. These findings underscore a low correlation between venom gland transcriptomes and proteomes and support the view that post-translational processes play a major role in shaping venom phenotypes.

Highlights of animal venom research on the geographical variations of toxin components, toxicities and envenomation therapy

Geographical variation of animal venom is common among venomous animals. This kind of intraspecific variation based on geographical location mainly concerned from envenomation cases and brought new problems in animal venom studies, including venom components regulatory mechanisms, differentiation of venom activities, and clinical treatment methods. At present, food is considered as the most related factor influencing venom development. Related research defined the variational venomous animal species by the comparison of venom components and activities in snakes, jellyfish, scorpions, cone snails, ants, parasitoid wasps, spiders and toads. In snake venom studies, researchers found that antivenom effectiveness was variated to different located venom samples. As described in some snake venom research, developing region-specific antivenom is the development trend. The difficulties of developing region-specific antivenom and theoretical solutions have been discussed. This review summarized biological studies of animal venom geographical variation by species, compared venom components and major biological activities of the vary venom from the same species, and listed the basic methods in comparing venom protein compositions and major toxicity differences to provide a comprehensive reference.

Toxicological comparison of Crotalus ruber lucasensis venom from different ecoregions of the Baja California Peninsula

The Baja California Peninsula possesses a mosaic of ecoregions that offers a wide variety of environments for the species that here inhabit. Here we report biological variations in. Crotalus ruber lucasensis venom from arid, semiarid and tropical eco-regions. Lethal (1.4-6.8 mg/kg), edematogenic (0.3-0.5 μg) and defibrinogenating (from non-detectable to 20 μg) activities were found to have significant differences among eco-regions.

Intra-Specific Venom Variation in the Australian Coastal Taipan Oxyuranus scutellatus

Intra-specific venom variation has the potential to provide important insights into the evolution of snake venom, but remains a relatively neglected aspect of snake venom studies. We investigated the venom from 13 individual coastal taipans Oxyuranus scutellatus from four localities on the north-east coast of Australia, spanning a distance of 2000 km. The intra-specific variation in taipan venom was considerably less than the inter-specific variation between it and the other Australian elapids to which it was compared. The electrophoretic venom profile of O. scutellatus was visually different to six other genera of Australian elapids, but not to its congener inland taipan O. microlepidotus. There was minimal geographical variation in taipan venom, as the intra-population variation exceeded the inter-population variation for enzymatic activity, procoagulant activity, and the abundance of neurotoxins. The pre-synaptic neurotoxin (taipoxin) was more abundant than the post-synaptic neurotoxins (3FTx), with a median of 11.0% (interquartile range (IQR): 9.7% to 18.3%; range: 6.7% to 23.6%) vs. a median of 3.4% (IQR: 0.4% to 6.7%; range: 0% to 8.1%). Three taipan individuals almost completely lacked post-synaptic neurotoxins, which was not associated with geography and occurred within two populations. We found no evidence of sexual dimorphism in taipan venom. Our study provides a basis for evaluating the significance of intra-specific venom variation within a phylogenetic context by comparing it to the inter-specific and inter-generic variation. The considerable intra-population variation we observed supports the use of several unpooled individuals from each population when making inter-specific comparisons.

Bite by a juvenile Bothrops venezuelensis (Venezuelan lancehead) resulting in severe envenomation: A case report

Bothrops venezuelensis is a venomous snake of the Viperidae family. It is associated with a high snakebite-related morbidity and mortality in Venezuela, although clinical case descriptions are scarce. Bites by other Bothrops sp. can result in coagulopathy and acute kidney injury. We describe a bite by a captive juvenile B. venezuelensis that caused local swelling, severe pain, endothelial damage, excessive fibrinolysis (INR >12, aPTT 136s, fibrinogen 0.3g/l) and incoagulable blood within 1.5 hours after the bite. The patient was treated with prothrombin complex factors concentrate, fibrinogen and antivenom (Antivipmyn®, Instituto Bioclon, Mexico) 4.5 h after the bite, which improved coagulation parameters progressively. Subsequently signs of compensated disseminated intravascular coagulation manifested and the patient received fresh frozen plasma and erythrocyte concentrate. The patient developed acute kidney injury with macroscopic hematuria. Fluid overload resulted in pulmonary edema requiring intermittent ventilation and diuretic treatment with furosemide. He was discharged with moderately elevated creatinine 16 days after hospitalization. Creatinine level normalized within another week. This case displays the life-threatening toxicity even after juvenile B. venezuelensis bites and the comparability to bites by other Bothrops sp.

Intact protein mass spectrometry reveals intraspecies variations in venom composition of a local population of Vipera kaznakovi in Northeastern Turkey

We report on the variable venom composition of a population of the Caucasus viper (Vipera kaznakovi) in Northeastern Turkey. We applied a combination of venom gland transcriptomics, de-complexing bottom-up and top-down venomics. In contrast to sole bottom-up venomics approaches and gel or chromatography based venom comparison, our combined approach enables a faster and more detailed comparison of venom proteomes from multiple individuals. In total, we identified peptides and proteins from 15 toxin families, including snake venom metalloproteinases (svMP; 37.8%), phospholipases A2 (PLA2; 19.0%), snake venom serine proteinases (svSP; 11.5%), C-type lectins (CTL; 6.9%) and cysteine-rich secretory proteins (CRISP; 5.0%), in addition to several low abundant toxin families. Furthermore, we identified intraspecies variations of the venom composition of V. kaznakovi, and find these were mainly driven by the age of the animals, with lower svSP abundance detected in juveniles. On the proteoform level, several small molecular weight toxins between 5 and 8 kDa in size, as well as PLA2s, drove the differences observed between juvenile and adult individuals. This study provides novel insights into the venom variability of V. kaznakovi and highlights the utility of intact mass profiling for fast and detailed comparison of snake venom. BIOLOGICAL SIGNIFICANCE: Population level and ontogenetic venom variation (e.g. diet, habitat, sex or age) can result in a loss of antivenom efficacy against snakebites from wide ranging snake populations. The current state of the art for the analysis of snake venoms are de-complexing bottom-up proteomics approaches. While useful, these have the significant drawback of being time-consuming and following costly protocols, and consequently are often applied to pooled venom samples. To overcome these shortcomings and to enable rapid and detailed profiling of large numbers of individual venom samples, we integrated an intact protein analysis workflow into a transcriptomics-guided bottom-up approach. The application of this workflow to snake individuals of a local population of V. kaznakovi revealed intraspecies variations in venom composition, which are primarily explained by the age of the animals, and highlighted svSP abundance to be one of the molecular drivers for the compositional differences observed.

Insights into individual variations in nematocyst venoms from the giant jellyfish Nemopilema nomurai in the Yellow Sea

The giant jellyfish, Nemopilema nomurai, is widely distributed from the Eastern China Sea to the northern part of the Yellow Sea and has resulted in numerous hospitalizations in coastal areas of China, especially in Northern China. Our previous studies have revealed sting-related proteins in the venom of the jellyfish N. nomurai by using experimental and omics-based approaches; however, the variable symptoms of patients who have been stung by N. nomurai are not fully understood. This limited knowledge led to an examination of whether intraspecific variations occur in the venom of different N. nomurai. In the present study, 13 specimens of N. nomurai were collected from the Yellow Sea, and their venom was characterized by profiling differences in biochemical properties and biological activities. SDS-PAGE analysis presented recognizable differences in the number, intensity and presence of some protein bands. Moreover, enzymatic assays revealed considerable quantitative variations in metalloproteinase activity and PLA₂-like activity. In particular, zymography assays of proteases demonstrated the general presence of abundant metalloproteinases in jellyfish nematocyst venom; however, the catalytic activities varied greatly among some specific metalloproteinases in the 28–46 kDa or 57–83 kDa range. Hemolytic assays using sheep erythrocytes suggested a predominant variance in the toxicities of different individual jellyfish venoms, with the difference between the most hemolytic and the least hemolytic venom as large as 77-fold. The current data suggested remarkable variations in the nematocyst venoms of individual N. nomurai jellyfish. These observations will provide a new understanding of the clinical manifestations induced by N. nomurai jellyfish stings and will therefore have important implications for preventing and treating jellyfish envenomations.

Bradykinin-potentiating PEPTIDE-10C, an argininosuccinate synthetase activator, protects against H2O2-induced oxidative stress in SH-SY5Y neuroblastoma cells

Bradykinin-potentiating peptides (BPPs - 5a, 7a, 9a, 10c, 11e, and 12b) of Bothrops jararaca (Bj) were described as argininosuccinate synthase (AsS) activators, improving l-arginine availability. Agmatine and polyamines, which are l-arginine metabolism products, have neuroprotective properties. Here, we investigated the neuroprotective effects of low molecular mass fraction from Bj venom (LMMF) and two synthetic BPPs (BPP-10c, <ENWPHPQIPP; BPP-12b, <EWGRPPGPPIPP) in the SH-SY5Y cell line against H₂O₂-induced oxidative stress. The neuroprotective effects against H₂O₂-induced were analyzed by reactive oxygen species (ROS - DCFH) production; lipid peroxidation (TBARS); intracellular GSH; AsS, iNOS, and NF-kB expressions; nitrite levels (Griess); mitochondrial membrane potential (TMRM); and antioxidant activity (DPPH). Analysis of variance followed by Tukey's post hoc test were calculated for statistical comparisons. Pre-treatment with both BPPs significantly reduced cell death induced by H₂O₂, but BPP-10c showed higher protective capacity than BPP-12b. LMMF pretreatment was unable to prevent the reduction of cell viability caused by H₂O₂. The neuroprotective mechanism of BPP-10c against oxidative stress was investigated. BPP-10c reduced ROS generation and lipid peroxidation in relation to cells treated only with H₂O₂. BBP-10c increased AsS expression and was not neuroprotective in the presence of MDLA, a specific inhibitor of AsS. BPP-10c reduced iNOS expression and nitrate levels but decreased NF-kB expression. Furthermore, BPP-10c protected the mitochondrial membrane against oxidation. Overall, we demonstrated for the first time neuroprotective mechanisms of BPPs against oxidative stress, opening new perspectives to the study and application of these peptides for the treatment of neurodegenerative diseases.

Proteopeptidomic, Functional and Immunoreactivity Characterization of Bothrops moojeni Snake Venom: Influence of Snake Gender on Venom Composition

Venom composition varies across snakes from all taxonomic levels and is influenced by the snakes’ age, habitat, diet, and sexual dimorphism. The present study reports the first in-depth investigation of venom composition in male and female Bothrops moojeni (B. moojeni) snakes (BmooM and BmooF, respectively) through three proteomics approaches associated with functional, cytotoxic, and immunoreactivity characterization. Compared with BmooM venom, BmooF venom exhibited weaker hyaluronidase, metalloproteinase, and phospholipase activity; stronger recognition by anti-bothropic serum; 1.4-fold stronger cytotoxicity; and greater number of peptides. The increased L-amino acid oxidase expression probably accounted for the stronger immunoreactivity and cytotoxicity of BmooF venom. BmooF and BmooM venom shared only 19% peptides. Some venom components were gender-specific, such as phospholipases B, phospholipase inhibitor, and hyaluronidases in BmooM, and cysteine-rich secretory proteins in BmooF. In conclusion, we describe herein the first proteomics study of B. moojeni snake venom and an in-depth characterization of gender-specific differences in venom composition. Altogether, our findings not only stress the importance of considering the snake’s gender during antivenom production, but also help to identify new potential drugs and biotechnological tools.

Venomics of Tropidolaemus wagleri, the sexually dimorphic temple pit viper: Unveiling a deeply conserved atypical toxin arsenal

Tropidolaemus wagleri (temple pit viper) is a medically important snake in Southeast Asia. It displays distinct sexual dimorphism and prey specificity, however its venomics and inter-sex venom variation have not been thoroughly investigated. Applying reverse-phase HPLC, we demonstrated that the venom profiles were not significantly affected by sex and geographical locality (Peninsular Malaya, insular Penang, insular Sumatra) of the snakes. Essentially, venoms of both sexes share comparable intravenous median lethal dose (LD₅₀) (0.56-0.63 μg/g) and cause neurotoxic envenomation in mice. LCMS/MS identified six waglerin forms as the predominant lethal principles, comprising 38.2% of total venom proteins. Fourteen other toxin-protein families identified include phospholipase A₂, serine proteinase, snaclec and metalloproteinase. In mice, HPLC fractions containing these proteins showed insignificant contribution to the overall venom lethality. Besides, the unique elution pattern of approximately 34.5% of non-lethal, low molecular mass proteins (3-5 kDa) on HPLC could be potential biomarker for this primitive crotalid species. Together, the study unveiled the venom proteome of T. wagleri that is atypical among many pit vipers as it comprises abundant neurotoxic peptides (waglerins) but little hemotoxic proteinases. The findings also revealed that the venom is relatively well conserved intraspecifically despite the drastic morphological differences between sexes.

Neutralization of toxicological activities of medically-relevant Bothrops snake venoms and relevant toxins by two polyvalent bothropic antivenoms produced in Peru and Brazil

Snakebite envenoming is a neglected public pathology, affecting especially rural communities or isolated areas of tropical and subtropical Latin American countries. The parenteral administration of antivenom is the mainstay and the only validated treatment of snake bite envenoming. Here, we assess the efficacy of polyspecific anti-Bothrops serum (α-BS) produced in the Instituto Nacional de Salud (INS, Peru) and at the Fundação Ezequiel Dias (FUNED, Brazil), to neutralize the main toxic activities induced by five medically-relevant venoms of: Bothrops atrox, B. barnetti, and B. pictus from Peru, and the Brazilian B. jararaca and B. leucurus, all of them inhabiting different geographical locations. Protein electrophoretic patterns of these venoms showed significant differences in composition, number and intensity of bands. Another goal was to evaluate the efficacy and safety of lyophilized α-BS developed at INS to neutralize the detrimental effects of these venoms using in vivo and in vitro assays. The availability of lyophilized α-BS has relevant significance in its distribution to distant rural communities where the access to antivenom in health facilities is more difficult. Despite the fact that different antigen mixtures were used for immunization during antivenom production, our data showed high toxin-neutralizing activity of α-BS raised against Bothrops venoms. Moreover, the antivenom cross-reacted even against venoms not included in the immunization mixture. Furthermore, we have evaluated the efficacy of both α-BS to neutralize key toxic compounds belonging to the predominant protein families of Bothrops snakes. Most significantly, both α-BS cross-specifically neutralized the main toxicological activities e.g. lethality and hemorrhage induced by these venoms. Thus, our data indicate that both α-BS are equally effective to treat snake bite victims inflicted by Bothrops snakes particularly B. atrox, responsible for the largest numbers of human envenomations in the Amazon regions of some South American countries including Peru and Brazil.

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.

Advances in venomics

The term “venomics” was coined to describe the global study of venom and venom glands, targeting comprehensive characterization of the whole toxin profile of a venomous animal by means of proteomics, transcriptomics, genomics and bioinformatics studies.

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 venoms are integrated systems, but abundant venom proteins evolve more rapidly

Background

While many studies have shown that extracellular proteins evolve rapidly, how selection acts on them remains poorly understood. We used snake venoms to understand the interaction between ecology, expression level, and evolutionary rate in secreted protein systems. Venomous snakes employ well-integrated systems of proteins and organic constituents to immobilize prey. Venoms are generally optimized to subdue preferred prey more effectively than non-prey, and many venom protein families manifest positive selection and rapid gene family diversification. Although previous studies have illuminated how individual venom protein families evolve, how selection acts on venoms as integrated systems, is unknown.

Results

Using next-generation transcriptome sequencing and mass spectrometry, we examined microevolution in two pitvipers, allopatrically separated for at least 1.6 million years, and their hybrids. Transcriptomes of parental species had generally similar compositions in regard to protein families, but for a given protein family, the homologs present and concentrations thereof sometimes differed dramatically. For instance, a phospholipase A₂ transcript comprising 73.4 % of the Protobothrops elegans transcriptome, was barely present in the P. flavoviridis transcriptome (<0.05 %). Hybrids produced most proteins found in both parental venoms. Protein evolutionary rates were positively correlated with transcriptomic and proteomic abundances, and the most abundant proteins showed positive selection. This pattern holds with the addition of four other published crotaline transcriptomes, from two more genera, and also for the recently published king cobra genome, suggesting that rapid evolution of abundant proteins may be generally true for snake venoms. Looking more broadly at Protobothrops, we show that rapid evolution of the most abundant components is due to positive selection, suggesting an interplay between abundance and adaptation.

Conclusions

Given log-scale differences in toxin abundance, which are likely correlated with biosynthetic costs, we hypothesize that as a result of natural selection, snakes optimize return on energetic investment by producing more of venom proteins that increase their fitness. Natural selection then acts on the additive genetic variance of these components, in proportion to their contributions to overall fitness. Adaptive evolution of venoms may occur most rapidly through changes in expression levels that alter fitness contributions, and thus the strength of selection acting on specific secretome components.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1832-6) contains supplementary material, which is available to authorized users.

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.

VTBuilder: a tool for the assembly of multi isoform transcriptomes

Background

Within many research areas, such as transcriptomics, the millions of short DNA fragments (reads) produced by current sequencing platforms need to be assembled into transcript sequences before they can be utilized. Despite recent advances in assembly software, creating such transcripts from read data harboring isoform variation remains challenging. This is because current approaches fail to identify all variants present or they create chimeric transcripts within which relationships between co-evolving sites and other evolutionary factors are disrupted. We present VTBuilder, a tool for constructing non-chimeric transcripts from read data that has been sequenced from sources containing isoform complexity.

Results

We validated VTBuilder using reads simulated from 54 Sanger sequenced transcripts (SSTs) expressed in the venom gland of the saw scaled viper, Echis ocellatus. The SSTs were selected to represent genes from major co-expressed toxin groups known to harbor isoform variants. From the simulated reads, VTBuilder constructed 55 transcripts, 50 of which had a greater than 99% sequence similarity to 48 of the SSTs. In contrast, using the popular assembler tool Trinity (r2013-02-25), only 14 transcripts were constructed with a similar level of sequence identity to just 11 SSTs. Furthermore VTBuilder produced transcripts with a similar length distribution to the SSTs while those produced by Trinity were considerably shorter. To demonstrate that our approach can be scaled to real world data we assembled the venom gland transcriptome of the African puff adder Bitis arietans using paired-end reads sequenced on Illumina’s MiSeq platform. VTBuilder constructed 1481 transcripts from 5 million reads and, following annotation, all major toxin genes were recovered demonstrating reconstruction of complex underlying sequence and isoform diversity.

Conclusion

Unlike other approaches, VTBuilder strives to maintain the relationships between co-evolving sites within the constructed transcripts, and thus increases transcript utility for a wide range of research areas ranging from transcriptomics to phylogenetics and including the monitoring of drug resistant parasite populations. Additionally, improving the quality of transcripts assembled from read data will have an impact on future studies that query these data. VTBuilder has been implemented in java and is available, under the GPL GPU V0.3 license, from http:// http://www.lstmed.ac.uk/vtbuilder.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0389-8) contains supplementary material, which is available to authorized users.