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.

Structures of N-Glycans of Bothrops Venoms Revealed as Molecular Signatures that Contribute to Venom Phenotype in Viperid Snakes.

The complexity of snake venoms has long been investigated to explore a myriad of biologically active proteins and peptides that are used for immobilizing or killing prey, and are responsible for the pathological effects observed on envenomation. Glycosylation is the main post-translational modification (PTM) of viperid venoms but currently there is little understanding of how protein glycosylation impacts the variation of venom proteomes. We have previously reported that Bothrops venom glycoproteomes contain a core of components that markedly define their composition and parallel their phylogenetic classification. Here we extend those observations to eight Bothrops species evaluating the N-glycomes by LC-MS as assigned cartoon structures and detailing those structures separately as methylated analogs using ion-trap mass spectrometry (MSn). Following ion disassembly through multiple steps provided sequence and linkage isomeric details that characterized 52 unique compositions in Bothrops venoms. These occurred as 60 structures, of which 26 were identified in the venoms of the Jararaca Complex (B. alcatraz, B. insularis, and B. jararaca), 20 in B. erythromelas, B. jararacussu, B. moojeni and B. neuwiedi venoms, and 22 in B. cotiara venom. Further, quantitative analysis of these N-glycans showed variable relative abundances in the venoms. For the first time a comprehensive set of N-glycan structures present in snake venoms are defined. Despite the fact that glycosylation is not template-defined, the N-glycomes of these venoms mirror the phylogeny cladograms of South American bothropoid snakes reported in studies on morphological, molecular data and feeding habits, exhibiting distinct molecular signatures for each venom. Considering the complexity of N-glycan moieties generally found in glycoproteins, characterized by different degrees of branching, isomer structures, and variable abundances, our findings point to these factors as another level of complexity in Bothrops venoms, features that could dramatically contribute to their distinct biological activities.

Dynamic Rearrangement in Snake Venom Gland Proteome: Insights into Bothrops jararaca Intraspecific Venom Variation.

We carried out an analysis of the venom gland proteome of Bothrops jararaca taking into account two distinct phases of its ontogenetic development (i.e., newborn and adult) and the marked sexual dimorphism recently reported on its venom proteome. Proteomic data analysis showed a dynamic rearrangement in the proteome landscape of B. jararaca venom gland upon development and gender-related changes. Differentially expressed proteins covered a number of biological pathways related to protein synthesis, including proteins associated with transcription and translation, which were found to be significantly higher expressed in the newborn venom gland. Our results suggest that the variation in the expression levels of cellular proteins might give rise to an even higher variation in the levels of the expressed toxins. Upon aging, the venom gland proteome repertoire related to the protein synthesis together with ecological traits would have an impact on the toxin repertoire, which, in the case of B. jararaca species, would enable the species to deal with different prey types during its lifespan. Proteomic data are available via ProteomeXchange with identifier PXD004186.

Proteomic and Glycoproteomic Profilings Reveal That Post-translational Modifications of Toxins Contribute to Venom Phenotype in Snakes.

Snake venoms are biological weapon systems composed of secreted proteins and peptides that are used for immobilizing or killing prey. Although post-translational modifications are widely investigated because of their importance in many biological phenomena, we currently still have little understanding of how protein glycosylation impacts the variation and stability of venom proteomes. To address these issues, here we characterized the venom proteomes of seven Bothrops snakes using a shotgun proteomics strategy. Moreover, we compared the electrophoretic profiles of native and deglycosylated venoms and, in order to assess their subproteomes of glycoproteins, we identified the proteins with affinity for three lectins with different saccharide specificities and their putative glycosylation sites. As proteinases are abundant glycosylated toxins, we examined the effect of N-deglycosylation on their catalytic activities and show that the proteinases of the seven venoms were similarly affected by removal of N-glycans. Moreover, we prospected putative glycosylation sites of transcripts of a B. jararaca venom gland data set and detected toxin family related patterns of glycosylation. Based on our global analysis, we report that Bothrops venom proteomes and glycoproteomes contain a core of components that markedly define their composition, which is conserved upon evolution in parallel to other molecular markers that determine their phylogenetic classification.

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.

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.

Systemic toxicity of snake venom metalloproteinases: Multi-omics analyses of kidney and blood plasma disturbances in a mouse model.

Snakebite envenomation is classified as a Neglected Tropical Disease. Bothrops jararaca venom induces kidney injury and coagulopathy. HF3, a hemorrhagic metalloproteinase of B. jararaca venom, participates in the envenomation pathogenesis. We evaluated the effects of HF3 in mouse kidney and blood plasma after injection in the thigh muscle, mimicking a snakebite. Transcriptomic analysis showed differential expression of 31 and 137 genes related to kidney pathology after 2 h and 6 h, respectively. However, only subtle changes were observed in kidney proteome, with differential abundance of 15 proteins after 6 h, including kidney injury markers. N-terminomic analysis of kidney proteins showed 420 proteinase-generated peptides compatible with meprin specificity, indicating activation of host proteinases. Plasma analysis revealed differential abundance of 90 and 219 proteins, respectively, after 2 h and 6 h, including coagulation-cascade and complement-system components, and creatine-kinase, whereas a semi-specific search of N-terminal peptides indicated activation of endogenous proteinases. HF3 promoted host reactions, altering the gene expression and the proteolytic profile of kidney tissue, and inducing plasma proteome imbalance driven by changes in abundance and proteolysis. The overall response of the mouse underscores the systemic action of a hemorrhagic toxin that transcends local tissue damage and is related to known venom-induced systemic effects.

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.

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.

Community-based network analyses reveal emerging connectivity patterns of protein-protein interactions in murine melanoma secretome.

Protein-protein interaction networks (PPINs) are static representations of protein connections in which topological features such as subgraphs (communities) may contain proteins functionally related, revealing an additional layer of interactome complexity. We created two PPINs from the secretomes of a paired set of murine melanocytes (a normal melanocyte and its transformed phenotype). Community structures, identified by a graph clustering algorithm, resulted in the identification of subgraphs in both networks. Interestingly, the underlying structure of such communities revealed shared and exclusive proteins (core and exclusive nodes, respectively), in addition to proteins that changed their location within each community (rewired nodes). Functional enrichment analysis of core nodes revealed conserved biological functions in both networks whereas exclusive and rewired nodes in the tumoral phenotype network were enriched in cancer-related processes, including TGFß signaling. We found a remarkable shift in the tumoral interactome, resulting in an emerging pattern which was driven by the presence of exclusive nodes and may represent functional network motifs. Our findings suggest that the rearrangement in the tumoral interactome may be correlated with the malignant transformation of melanocytes associated with substrate adhesion impediment. The interactions found in core and new/rewired nodes might potentially be targeted for therapeutic intervention in melanoma treatment. SIGNIFICANCE: Malignant transformation is a result of synergistic action of multiple molecular factors in which genetic alterations as well as protein expression play paramount roles. During oncogenesis, cellular crosstalk through the secretion of soluble mediators modulates the phenotype of transformed cells which ultimately enables them to successfully disrupt important signaling pathways, including those related to cell growth and proliferation. Therefore, in this work we profiled the secretomes of a paired set of normal and transformed phenotypes of a murine melanocyte. After assembling the two interactomes, clusters of functionally related proteins (network communities) were observed as well as emerging patterns of network rewiring which may represent an interactome signature of transformed cells. In summary, the significance of this study relies on the understanding of the repertoire of 'normal' and 'tumoral' secretomes and, more importantly, the set of interacting proteins (the interactome) in both of these conditions, which may reveal key components that might be potentially targeted for therapeutic intervention.

Cell Cycle, Telomeres, and Telomerase in Leishmania spp.: What Do We Know So Far?

Leishmaniases belong to the inglorious group of neglected tropical diseases, presenting different degrees of manifestations severity. It is caused by the transmission of more than 20 species of parasites of the Leishmania genus. Nevertheless, the disease remains on the priority list for developing new treatments, since it affects millions in a vast geographical area, especially low-income people. Molecular biology studies are pioneers in parasitic research with the aim of discovering potential targets for drug development. Among them are the telomeres, DNA-protein structures that play an important role in the long term in cell cycle/survival. Telomeres are the physical ends of eukaryotic chromosomes. Due to their multiple interactions with different proteins that confer a likewise complex dynamic, they have emerged as objects of interest in many medical studies, including studies on leishmaniases. This review aims to gather information and elucidate what we know about the phenomena behind Leishmania spp. telomere maintenance and how it impacts the parasite's cell cycle.

Co-Exposure of Cardiomyocytes to IFN-γ and TNF-α Induces Mitochondrial Dysfunction and Nitro-Oxidative Stress: Implications for the Pathogenesis of Chronic Chagas Disease Cardiomyopathy.

Infection by the protozoan Trypanosoma cruzi causes Chagas disease cardiomyopathy (CCC) and can lead to arrhythmia, heart failure and death. Chagas disease affects 8 million people worldwide, and chronic production of the cytokines IFN-γ and TNF-α by T cells together with mitochondrial dysfunction are important players for the poor prognosis of the disease. Mitochondria occupy 40% of the cardiomyocytes volume and produce 95% of cellular ATP that sustain the life-long cycles of heart contraction. As IFN-γ and TNF-α have been described to affect mitochondrial function, we hypothesized that IFN-γ and TNF-α are involved in the myocardial mitochondrial dysfunction observed in CCC patients. In this study, we quantified markers of mitochondrial dysfunction and nitro-oxidative stress in CCC heart tissue and in IFN-γ/TNF-α-stimulated AC-16 human cardiomyocytes. We found that CCC myocardium displayed increased levels of nitro-oxidative stress and reduced mitochondrial DNA as compared with myocardial tissue from patients with dilated cardiomyopathy (DCM). IFN-γ/TNF-α treatment of AC-16 cardiomyocytes induced increased nitro-oxidative stress and decreased the mitochondrial membrane potential (ΔΨm). We found that the STAT1/NF-κB/NOS2 axis is involved in the IFN-γ/TNF-α-induced decrease of ΔΨm in AC-16 cardiomyocytes. Furthermore, treatment with mitochondria-sparing agonists of AMPK, NRF2 and SIRT1 rescues ΔΨm in IFN-γ/TNF-α-stimulated cells. Proteomic and gene expression analyses revealed that IFN-γ/TNF-α-treated cells corroborate mitochondrial dysfunction, transmembrane potential of mitochondria, altered fatty acid metabolism and cardiac necrosis/cell death. Functional assays conducted on Seahorse respirometer showed that cytokine-stimulated cells display decreased glycolytic and mitochondrial ATP production, dependency of fatty acid oxidation as well as increased proton leak and non-mitochondrial oxygen consumption. Together, our results suggest that IFN-γ and TNF-α cause direct damage to cardiomyocytes' mitochondria by promoting oxidative and nitrosative stress and impairing energy production pathways. We hypothesize that treatment with agonists of AMPK, NRF2 and SIRT1 might be an approach to ameliorate the progression of Chagas disease cardiomyopathy.

Heterotypic signaling between dermal fibroblasts and melanoma cells induces phenotypic plasticity and proteome rearrangement in malignant cells.

The signaling events triggered by soluble mediators released from both transformed and stromal cells shape the phenotype of tumoral cells and have significant implications in cancer development and progression. In this study we performed an in vitro heterotypic signaling assays by evaluating the proteome diversity of human dermal fibroblasts after stimulation with the conditioned media obtained from malignant melanoma cells. In addition, we also evaluated the changes in the proteome of melanoma cells after stimulation with their own conditioned media as well as with the conditioned medium from melanoma-stimulated fibroblasts. Our results revealed a clear rearrangement in the proteome of stromal and malignant cells upon crosstalk of soluble mediators. The main proteome signature of fibroblasts stimulated with melanoma conditioned medium was related to protein synthesis, which indicates that this process might be an early response of stromal cells. In addition, the conditioned medium derived from 'primed' stromal cells (melanoma-stimulated fibroblasts) was more effective in altering the functional phenotype (cell migration) of malignant cells than the conditioned medium from non-stimulated fibroblasts. Collectively, self- and cross-stimulation may play a key role in shaping the tumor microenvironment and enable tumoral cells to succeed in the process of melanoma progression and metastasis. Although the proteome landscape of cells participating in such a heterotypic signaling represents a snapshot of a highly dynamic state, understanding the diversity of proteins and enriched biological pathways resulting from stimulated cell states may allow for targeting specific cell regulatory motifs involved in melanoma progression and metastasis.

Comparative analysis of the high molecular mass subproteomes of eight Bothrops snake venoms.

Snake venoms are extremely active biological secretions composed primarily of various classes of enzymes. The genus Bothrops comprises various pit viper species that represent the most medically significant taxa in Central and South America, accounting for more human envenomations and fatalities than any other snakes in the region. Venom proteomes of many Bothrops species have been well-characterized but investigations have focused almost exclusively on proteins smaller than 100 kDa despite expression of larger components being documented in several Bothrops venoms. This study sought to achieve detailed identification of major components in the high molecular mass subproteome of venoms from eight Bothrops species (B. brazili, B. cotiara, B. insularis, B. jararaca, B. jararacussu, B. leucurus, B. moojeni and B. neuwiedi). Enzymes such as metalloproteinases and L-amino acid oxidases were the most prominent components identified in the first size-exclusion chromatography fractions of these venoms. Minor components also identified in the first peaks included 5'-nucleotidase, aminopeptidase, phosphodiesterase, and phospholipases A2 and B. Most of these components disappeared in electrophoretic profiles under reducing conditions, suggesting that they may be composed of more than one polypeptide chain. A significant shift in the molecular masses of these protein bands was observed following enzymatic N-deglycosylation, indicating that they may contain N-glycans. Furthermore, none of the identified high molecular mass proteins were shared by all eight species, revealing a high level of interspecific variability among these venom components.

Biochemical characterization of venom from Pseudoboa neuwiedii (Neuwied's false boa; Xenodontinae; Pseudoboini).

In this work, we examined the proteolytic and phospholipase A2 (PLA2) activities of venom from the opisthoglyphous colubrid Pseudoboa neuwiedii. Proteolytic activity (3 and 10 µg of venom) was comparable to that of Bothrops neuwiedii venom but less than Bothrops atrox. This activity was inhibited by EDTA and 1,10-phenanthroline but only slightly affected (≤30% inhibition) by PMSF and AEBSF, indicating it was mediated by snake venom metalloproteinases (SVMPs). The pH and temperature optima for proteolytic activity were 8.0 and 37 °C, respectively. The venom had no esterase activity, whereas PLA2 activity was similar to B. atrox, greater than B. neuwiedii but less than B. jararacussu. SDS-PAGE revealed venom proteins >100 kDa, 45-70 kDa, 21-24 kDa and ~15 kDa, and mass spectrometry of protein bands revealed SVMPs, cysteine-rich secretory proteins (CRISPs) and PLA2, but no serine proteinases. In gelatin zymography, the most active bands occurred at 65-68 kDa (seen with 0.05-0.25 µg of venom). Caseinolytic activity occurred at 50-66 kDa and was generally weaker than gelatinolytic activity. RP-HPLC of venom yielded 15 peaks, five of which showed gelatinolytic activity; peak 7 was the most active and apparently contained a P-III class SVMP. The venom showed α-fibrinogenase activity, without affecting the ß and γ chains; this activity was inhibited by EDTA and 1,10-phenanthroline. The venom did not clot rat citrated plasma but reduced the rate and extent of coagulation after plasma recalcification. In conclusion, P. neuwiedii venom is highly proteolytic and could potentially affect coagulation in vivo by degrading fibrinogen via SVMPs.

Functional and proteomic comparison of Bothrops jararaca venom from captive specimens and the Brazilian Bothropic Reference Venom.

Snake venom is a variable phenotypic trait, whose plasticity and evolution are critical for effective antivenom production. A significant reduction of the number of snake donations to Butantan Institute (São Paulo, Brazil) occurred in recent years, and this fact may impair the production of the Brazilian Bothropic Reference Venom (BBRV). Nevertheless, in the last decades a high number of Bothrops jararaca specimens have been raised in captivity in the Laboratory of Herpetology of Butantan Institute. Considering these facts, we compared the biochemical and biological profiles of B. jararaca venom from captive specimens and BBRV in order to understand the potential effects of snake captivity upon the venom composition. Electrophoretic analysis and proteomic profiling revealed few differences in venom protein bands and some differentially abundant toxins. Comparison of enzymatic activities showed minor differences between the two venoms. Similar cross-reactivity recognition pattern of both venoms by the antibothropic antivenom produced by Butantan Institute was observed. Lethality and neutralization of lethality for B. jararaca venom from captive specimens and BBRV showed similar values. Considering these results we suggest that the inclusion of B. jararaca venom from captive specimens in the composition of BBRV would not interfere with the quality of this reference venom. BIOLOGICAL SIGNIFICANCE: Snakebite envenomation is a neglected tropical pathology whose treatment is based on the use of specific antivenoms. Bothrops jararaca is responsible for the majority of snakebites in South and Southeastern Brazil. Its venom shows individual, sexual, and ontogenetic variability, however, the effect of animal captivity upon venom composition is unknown. Considering the reduced number of wild-caught snakes donated to Butantan Institute in the last decades, and the increased life expectancy of the snakes raised in captivity in the Laboratory of Herpetology, this work focused on the comparative profiling of B. jararaca venom from captive snakes and the Brazilian Bothropic Reference Venom (BBRV). BBRV is composed of venom obtained upon the first milking of wild-caught B. jararaca specimens, and used to assess the potency of all bothropic antivenoms produced by Brazilian suppliers. The use of proteomic strategies, added to biochemical and neutralization tests, allowed to conclude that, despite some subtle differences detected between these two venoms, venom from captive specimens could be used in the BBRV composition without affecting its quality in antivenom potency assays.

Trends in the Evolution of Snake Toxins Underscored by an Integrative Omics Approach to Profile the Venom of the Colubrid Phalotris mertensi.

Only few studies on snake venoms were dedicated to deeply characterize the toxin secretion of animals from the Colubridae family, despite the fact that they represent the majority of snake diversity. As a consequence, some evolutionary trends observed in venom proteins that underpinned the evolutionary histories of snake toxins were based on data from a minor parcel of the clade. Here, we investigated the proteins of the totally unknown venom from Phalotris mertensi (Dipsadinae subfamily), in order to obtain a detailed profile of its toxins and to appreciate evolutionary tendencies occurring in colubrid venoms. By means of integrated omics and functional approaches, including RNAseq, Sanger sequencing, high-resolution proteomics, recombinant protein production, and enzymatic tests, we verified an active toxic secretion containing up to 21 types of proteins. A high content of Kunitz-type proteins and C-type lectins were observed, although several enzymatic components such as metalloproteinases and an L-amino acid oxidase were also present in the venom. Interestingly, an arguable venom component of other species was demonstrated as a true venom protein and named svLIPA (snake venom acid lipase). This finding indicates the importance of checking the actual protein occurrence across species before rejecting genes suggested to code for toxins, which are relevant for the discussion about the early evolution of reptile venoms. Moreover, trends in the evolution of some toxin classes, such as simplification of metalloproteinases and rearrangements of Kunitz and Wap domains, parallel similar phenomena observed in other venomous snake families and provide a broader picture of toxin evolution.

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.