New findings from the first transcriptome of the Bothrops moojeni snake venom gland

Three snake venoms from Bothrops genus induced apoptosis and cell cycle arrest in K562 human leukemic cell line

Cancer is indisputably one of the leading causes of death worldwide. Snake venoms are a potential source of bioactive compounds, complex mixtures constituted mainly of proteins and peptides with several pharmacological possibilities, including the potential to inhibit tumoral cell growth. In the present study, it was evaluated the antitumor effect of crude venom of Bothrops erythromelas (BeV), Bothrops jararaca (from Southern and Southeastern- BjsV and BjsdV, respectively) and Bothrops alternatus (BaV) in in vitro Chronic myeloid leukemia (CML) cancer cell line model. After 24 h of cell exposure to 10 and 50 μg/mL, BjsV, BjsdV, and BaV exerted a decrease in cell viability in both concentrations. BeV was not cytotoxic and, therefore wasn't chosen for further mechanism of action investigation. Furthermore, morphological alterations show modification typical of apoptosis. Also, was observes a significant cell cycle arrest in the S phase by BjsdV and BaV treatment. Flow cytometry evidenced the involvement of changes in the cell membrane permeability and the mitochondrial function by BjsV and BjsdV, corroborating with the triggering of the apoptotic pathway by the venom administration. BjsV, BjsdV, and BaV also led to extensive DNA damage and were shown to modulate the gene expression of transcripts related to the cell cycle progression and suppress the expression of the BCR-ABL1 oncogene. Altogether, these findings suggest that the venoms trigger the apoptosis pathway due to mitochondrial damage and cell cycle arrest, with modulation of intracellular pathways important for CML progression. Thus, indicating the pharmacological potential of these venoms in the development of new antitumoral compounds.

Effects of Mauritia flexuosa L. f. buriti oil on symptoms induced by Bothrops moojeni snake envenomation

ETHNOPHARMACOLOGICAL RELEVANCE

In Brazil, there are species of snakes that become involved in accidents and cause serious health problems to the inhabitants, highlighting the genus Bothrops for being responsible for approximately 90% of accidents reported annually. In the northern region of the country, this genus is responsible for the largest number of accidents, especially among rural dwellers. These populations invest in alternative treatments for with the purpose of improving the symptoms caused by snakebites. The species Mauritia flexuosa L. f., known as buriti, is traditionally used for the treatment of envenomation by snakes.

AIM OF THE STUDY

This study aimed to evaluate the antiophidic potential of the oil of Mauritia flexuosa L. f. for Bothrops moojeni H. venom, confronting cultural and scientific knowledge.

MATERIALS AND METHODS

The physicochemical properties were determined, and the components present in the oil, extracted from fruit pulp, were analyzed by Gas Chromatography Coupled with Mass Spectrometry. The in vitro inhibitory capacity of the oil for phospholipase, metalloprotease and serine protease activities was investigated. In the in vivo studies, male Swiss mice were used to evaluate the effect of oil on lethality and toxicity, and hemorrhagic, myotoxic and edematogenic activities were assessed.

RESULTS

GC‒MS analysis identification of 90.95% of the constituents of the oil, with the main components being 9-eicosenoic acid, (Z)- (34.54%), n-hexadecanoic acid (25.55%) and (E)-9-octadecenoic acid ethyl ester (12.43%). For the substrates, the outcomes indicate that the oil inhibited the activity of the main classes of toxins present in Bothrops moojeni H. venom (VBm) at the highest dose tested (0.5 μL), with inhibition of 84% for the hydrolysis of the selective substrate for serine protease and inhibition of 60% for the hydrolysis of substrates for PLA₂ and metalloproteases. The antiophidic activity in vivo was evaluated with two concentrations of the oil: 1.5 mg, the dosage the population, diluted in mineral oil to a volume of 1 tablespoon and 15 mg, administered by gavage 30 min before poisoning and at time zero (concomitant to poisoning), and both concentrations administered by gavage in combination with topical use at time zero. The bleeding time in the group treated with oil at a concentration of 15 mg administered at time zero was significantly lower than that in the control group (p < 0.05). However, a greater inhibition of bleeding time was observed when local application was combined with the gavage treatment at both concentrations tested at time zero (p < 0.05). In the myotoxicity test, oil was efficient in reducing the myotoxic effects induced by the venom at the two concentrations tested, with gavage administration at time zero and gavage plus topical administration at time zero (p < 0.05).

CONCLUSIONS

The data obtained show that the oil is safe to use at the concentrations studied and contains fatty acids that may collaborate for cellular-level repair of the injuries caused by Bm poisoning. The in vitro and in vivo experiments showed that oil inhibits the main proteolytic enzymes present in the venom and that it has important activities to control the local effects caused by bothropic venom.

Isolation and characterization of the first phosphodiesterase (Bj-PDE) from the venom of Bothrops jararacussu snake

Phosphodiesterases are exonucleases that sequentially hydrolyse phosphodiester bonds of polynucleotides from the 3'-end and release 5-mononucleotides. After more than one decade without any advance in the study of Bothropic phosphodiesterases, we described here the isolation of the first phosphodiesterase from Bothrops jararacussu, which we named Bj-PDE. A five-step column chromatography procedure (size exclusion, hydrophobic interaction, cation exchange, lentil lectin affinity, and blue sepharose affinity) enabled isolation of Bj-PDE with preserved and stable enzymatic activity (bis(p-nitrophenyl) phosphate substrate), K_m = 6.9 mM (± 0.7 mM), k_cat/K_m = 1.7 × 10⁴ M^-1 s^-1 (± 0.2 × 10⁴ M^-1 s^-1), MW = 116 kDa (SDS-PAGE), optimum activity around 45 °C at pH 8.0, and stability for 81 days at different storage temperatures (8, -20, and - 80 °C). Ca²⁺ and Mg²⁺ ions positively influenced Bj-PDE activity, while EDTA had the opposite action. Zn²⁺ restored >50 % of enzyme activity after its inhibition by EDTA. The Bj-PDE partial sequence identified by mass spectrometry was very similar to the sequence of BATXPDE1 from Bothrops atrox, which was evolutionarily close to this new PDE. Therefore, our study represents an important progress on the isolation of this minor toxin and sheds new lights on the properties and bioprospection of bothropic phosphodiesterases.

Current situation of snakebites envenomation in the Neotropics: Biotechnology, a versatile tool in the production of antivenoms

Snakebite envenomation is a neglected tropical disease that affects millions of people around the world with a great impact on health and the economy. Unfortunately, public health programs do not include this kind of disease as a priority in their social programs. Cases of snakebite envenomations in the Neotropics are inaccurate due to inadequate disease management from medical records to the choice of treatments. Victims of snakebite envenomation are primarily found in impoverished agricultural areas where remote conditions limit the availability of antivenom. Antivenom serum is the only Food and Drug Administration-approved treatment used up to date. However, it has several disadvantages in terms of safety and effectiveness. This review provides a comprehensive insight dealing with the current epidemiological status of snakebites in the Neotropics and technologies employed in antivenom production. Also, modern biotechnological tools such as transcriptomic, proteomic, immunogenic, high-density peptide microarray and epitope mapping are highlighted for producing new-generation antivenom sera. These results allow us to propose strategic solutions in the Public Health Sector for managing this disease. Keywords: antivenom, biotechnology, neglected tropical disease, omics, recombinant antibody.

Differences in PLA2 Constitution Distinguish the Venom of Two Endemic Brazilian Mountain Lanceheads, Bothrops cotiara and Bothrops fonsecai

Interspecific differences in snake venom compositions can result from distinct regulatory mechanisms acting in each species. However, comparative analyses focusing on identifying regulatory elements and patterns that led to distinct venom composition are still scarce. Among venomous snakes, Bothrops cotiara and Bothrops fonsecai represent ideal models to complement our understanding of the regulatory mechanisms of venom production. These recently diverged species share a similar specialized diet, habitat, and natural history, but each presents a distinct venom phenotype. Here, we integrated data from the venom gland transcriptome and miRNome and the venom proteome of B. fonsecai and B. cotiara to better understand the regulatory mechanisms that may be acting to produce differing venom compositions. We detected not only the presence of similar toxin isoforms in both species but also distinct expression profiles of phospholipases A2 (PLA2) and some snake venom metalloproteinases (SVMPs) and snake venom serine proteinases (SVSPs) isoforms. We found evidence of modular expression regulation of several toxin isoforms implicated in venom divergence and observed correlated expression of several transcription factors. We did not find strong evidence for miRNAs shaping interspecific divergence of the venom phenotypes, but we identified a subset of toxin isoforms whose final expression may be fine-tuned by specific miRNAs. Sequence analysis on orthologous toxins showed a high rate of substitutions between PLA2s, which indicates that these toxins may be under strong positive selection or represent paralogous toxins in these species. Our results support other recent studies in suggesting that gene regulation is a principal mode of venom evolution across recent timescales, especially among species with conserved ecotypes.

Unraveling the structure and function of CdcPDE: A novel phosphodiesterase from Crotalus durissus collilineatus snake venom

This study reports the isolation, structural, biochemical, and functional characterization of a novel phosphodiesterase from Crotalus durissus collilineatus venom (CdcPDE). CdcPDE was successfully isolated from whole venom using three chromatographic steps and represented 0.7% of total protein content. CdcPDE was inhibited by EDTA and reducing agents, demonstrating that metal ions and disulfide bonds are necessary for its enzymatic activity. The highest enzymatic activity was observed at pH 8-8.5 and 37 °C. Kinetic parameters indicated a higher affinity for the substrate bis(p-nitrophenyl) phosphate compared to others snake venom PDEs. Its structural characterization was done by the determination of the protein primary sequence by Edman degradation and mass spectrometry, and completed by the building of molecular and docking-based models. Functional in vitro assays showed that CdcPDE is capable of inhibiting platelet aggregation induced by adenosine diphosphate in a dose-dependent manner and demonstrated that CdcPDE is cytotoxic to human keratinocytes. CdcPDE was recognized by the crotalid antivenom produced by the Instituto Butantan. These findings demonstrate that the study of snake venom toxins can reveal new molecules that may be relevant in cases of snakebite envenoming, and that can be used as molecular tools to study pathophysiological processes due to their specific biological activities.

Proteomic and Transcriptomic Techniques to Decipher the Molecular Evolution of Venoms

Nature's library of venoms is a vast and untapped resource that has the potential of becoming the source of a wide variety of new drugs and therapeutics. The discovery of these valuable molecules, hidden in diverse collections of different venoms, requires highly specific genetic and proteomic sequencing techniques. These have been used to sequence a variety of venom glands from species ranging from snakes to scorpions, and some marine species. In addition to identifying toxin sequences, these techniques have paved the way for identifying various novel evolutionary links between species that were previously thought to be unrelated. Furthermore, proteomics-based techniques have allowed researchers to discover how specific toxins have evolved within related species, and in the context of environmental pressures. These techniques allow groups to discover novel proteins, identify mutations of interest, and discover new ways to modify toxins for biomimetic purposes and for the development of new therapeutics.

A comparative study on the leishmanicidal activity of the L-amino acid oxidases BjussuLAAO-II and BmooLAAO-II isolated from Brazilian Bothrops snake venoms

This study aims to examine whether two L-amino acid oxidases isolated from Bothrops snake venom (SV-LAAOs) were cytotoxic to Leishmania (Leishmania) amazonensis and Leishmania (Viannia) braziliensis, two causative agents of leishmaniasis, which is an endemic disease in tropical and subtropical countries. The SV-LAAOs BjussuLAAO-II and BmooLAAO-II were isolated from Bothrops jararacussu and Bothrops moojeni venom, respectively, through a three-step chromatography process that used molecular exclusion, hydrophobic interaction, and affinity columns. BmooLAAO-II is a new SV-LAAO isoform that we isolated in this study. The purified BjussuLAAO-II and BmooLAAO-II had high L-amino acid oxidase-specific activity: 3481.17 and 4924.77 U/mg/min, respectively. Both SV-LAAOs were strongly cytotoxic to the two Leishmania species, even at low concentrations. At the same concentration, BjussuLAAO-II and BmooLAAO-II exerted different cytotoxic effects on the parasites. We reported for the first time that the SV-LAAOs suppressed cell proliferation and altered the mitochondrial membrane potential of the two Leishmania species. Surprisingly, BjussuLAAO-II increased the intracellular reactive oxygen species production only in L. (L.) amazonensis, while BmooLAAO-II increased the intracellular reactive oxygen species production only in L. (V.) braziliensis, indicating that these SV-LAAOs had a certain specificity of action.

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Background:

Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design.

Methods:

RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment.

Results:

We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin.

Conclusion:

It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.

The Sequence and Three-Dimensional Structure Characterization of Snake Venom Phospholipases B

Snake venom phospholipases B (SVPLBs) are the least studied enzymes. They constitute about 1% of Bothrops crude venoms, however, in other snake venoms, it is present in less than 1%. These enzymes are considered the most potent hemolytic agent in the venom. Currently, no structural information is available about these enzymes from snake venom. To better understand its three-dimensional structure and mechanisms of envenomation, the current work describes the first model-based structure report of this enzyme from Bothrops moojeni venom named as B. moojeni phospholipase B (PLB_Bm). The structure model of PLB_Bm was generated using model building software like I-TESSER, MODELLER 9v19, and Swiss-Model. The build PLB_Bm model was validated using validation tools (PROCHECK, ERRAT, and Verif3D). The analysis of the PLB_Bm modeled structure indicates that it contains 491 amino acid residues that form a well-defined four-layer αββα sandwich core and has a typical fold of the N-terminal nucleophile aminohydrolase (Ntn-hydrolase). The overall structure of PLB_Bm contains 18 β-strands and 17 α-helices with many connecting loops. The structure divides into two chains (A and B) after maturation. The A chain is smaller and contains 207 amino acid residues, whereas the B chain is larger and contains 266 amino acid residues. The sequence and structural comparison among homologous snake venom, bacterial, and mammals PLBs indicate that differences in the length and sequence composition may confer variable substrate specificity to these enzymes. Moreover, the surface charge distribution, average volume, and depth of the active site cavity also vary in these enzymes. The present work will provide more information about the structure-function relationship and mechanism of action of these enzymes in snakebite envenomation.

Photobiomodulation reduces cell death and cytokine production in C2C12 cells exposed to Bothrops venoms

Snakebites caused by the genus Bothrops are often associated with severe and complex local manifestations such as edema, pain, hemorrhage, and myonecrosis. Conventional treatment minimizes the systemic effects of venom; however, their local action is not neutralized. The purpose of this study was to evaluate the effect of photobiomodulation (PBM) on C2C12 muscle cells exposed to B. jararaca, B. jararacussu, and B. moojeni venoms on events involved in cell death and the release of inflammatory mediators. Cells were exposed to venoms and immediately irradiated with low-level laser (LLL) application in continuous wave at the wavelength of 660 nm, energy density of 4.4 J/cm², power of 10 mW, area of 0.045 cm², and time of 20 s. Cell integrity was analyzed by phase contrast microscope and cell death was performed by flow cytometry. In addition, interleukin IL1-β, IL-6, and IL-10 levels were measured in the supernatant. Our results showed that the application of PBM increases cell viability and decreases cell death by apoptosis and necrosis. Moreover, the release of pro-inflammatory interleukins was also reduced. The data reported here indicate that PBM resulted in cytoprotection on myoblast C2C12 cells after venom exposure. This protection involves the modulation of cell death mechanism and decreased pro-inflammatory cytokine release.

Proteome of fraction from Tityus serrulatus venom reveals new enzymes and toxins

Background

Tityus serrulatus venom (Ts venom) is a complex mixture of several compounds with biotechnological and therapeutical potentials, which highlights the importance of the identification and characterization of these components. Although a considerable number of studies have been dedicated to the characterization of this complex cocktail, there is still a limitation of knowledge concerning its venom composition. Most of Ts venom studies aim to isolate and characterize their neurotoxins, which are small, basic proteins and are eluted with high buffer concentrations on cation exchange chromatography. The first and largest fraction from carboxymethyl cellulose-52 (CMC-52) chromatography of Ts venom, named fraction I (Fr I), is a mixture of proteins of high and low molecular masses, which do not interact with the cation exchange resin, being therefore a probable source of components still unknown of this venom. Thus, the present study aimed to perform the proteome study of Fraction I from Ts venom, by high resolution mass spectrometry, and its biochemical characterization, by the determination of several enzymatic activities.

Methods

Fraction I was obtained by a cation exchange chromatography using 50 mg of crude venom. This fraction was subjected to a biochemical characterization, including determination of L-amino acid oxidase, phospholipase, hyaluronidase, proteases activities and inhibition of angiotensin converting enzyme (ACE) activity. Fraction I was submitted to reduction, alkylation and digestion processes, and the tryptic digested peptides obtained were analyzed in a Q-Exactive Orbitrap mass spectrometer. Data analysis was performed by PEAKS 8.5 software against NCBI database.

Results

Fraction I exhibits proteolytic activity and it was able to inhibit ACE activity. Its proteome analysis identified 8 different classes of venom components, among them: neurotoxins (48%), metalloproteinases (21%), hypotensive peptides (11%), cysteine-rich venom protein (9%), antimicrobial peptides (AMP), phospholipases and other enzymes (chymotrypsin and lysozymes) (3%) and phosphodiesterases (2%).

Conclusions

The combination of a proteomic and biochemical characterization strategies leads us to identify new components in the T. serrulatus scorpion venom. The proteome of venom´s fraction can provide valuable direction in the obtainment of components in their native forms in order to perform a preliminary characterization and, consequently, to promote advances in biological discoveries in toxinology.

Biochemical and molecular characterization of the hyaluronidase from Bothrops atrox Peruvian snake venom

Snake venoms are a rich source of enzymes such as metalloproteinases, serine proteinases phospholipases A2 and myotoxins, that have been well characterized structurally and functionally. However, hyaluronidases (E.C.3.2.1.35) have not been studied extensively. In this study, we describe the biochemical and molecular features of a hyaluronidase (Hyal-Ba) isolated from the venom of the Peruvian snake Bothrops atrox. Hyal-Ba was purified by a combination of ion-exchange and gel filtration chromatography. Purified Hyal-Ba is a 69-kDa (SDS-PAGE) monomeric glycoprotein with an N-terminal amino acid sequence sharing high identity with homologous snake venom hyaluronidases. Detected associated carbohydrates were hexoses (16.38%), hexosamines (2.7%) and sialic acid (0.69%). Hyal-Ba selectively hydrolyzed only hyaluronic acid (HA; specific activity = 437.5 U/mg) but it did not hydrolyze chondroitin sulfate or heparin. The optimal pH and temperature for maximum activity were 6.0 and 40 °C, respectively, and its Km was 0.31 μM. Its activity was inhibited by EDTA, iodoacetate, 2-mercaptoethanol, TLCK and dexamethasone. Na⁺ and K⁺ (0.2 M) positively affect hyaluronidase activity; while Mg²⁺, Br²⁺, Ba²⁺, Cu²⁺, Zn²⁺, and Cd²⁺ reduced catalytic activity. Hyal-Ba potentiates the hemorrhagic and hemolytic activity of whole venom, but decreased subplantar edema caused by an l-amino acid oxidase (LAAO). The Hyal-Ba cDNA sequence (2020 bp) encodes 449 amino acid residues, including the catalytic site residues (Glu135, Asp133, Tyr206, Tyr253 and Trp328) and three functional motifs for N-linked glycosylation, which are conserved with other snake hyaluronidases. Spatial modeling of Hyal-Ba displayed a TIM-Barrel (α/β) fold and an EGF-like domain in the C-terminal portion. The phylogenetic analysis of Hyal-Ba with other homologous Hyals showed the monophyly of viperids. Further, Hyal-Ba studies may extend our knowledge of B. atrox toxinology and provides insight to improve the neutralizing strategies of therapeutic antivenoms.

New Insights on Moojase, a Thrombin-Like Serine Protease from Bothrops moojeni Snake Venom

Snake venom serine proteases (SVSPs) are enzymes that are capable of interfering in various parts of the blood coagulation cascade, which makes them interesting candidates for the development of new therapeutic drugs. Herein, we isolated and characterized Moojase, a potent coagulant enzyme from Bothrops moojeni snake venom. The toxin was isolated from the crude venom using a two-step chromatographic procedure. Moojase is a glycoprotein with N-linked glycans, molecular mass of 30.3 kDa and acidic character (pI 5.80⁻6.88). Sequencing of Moojase indicated that it is an isoform of Batroxobin. Moojase was able to clot platelet-poor plasma and fibrinogen solutions in a dose-dependent manner, indicating thrombin-like properties. Moojase also rapidly induced the proteolysis of the Aα chains of human fibrinogen, followed by the degradation of the Bβ chains after extended periods of incubation, and these effects were inhibited by PMSF, SDS and DTT, but not by benzamidine or EDTA. RP-HPLC analysis of its fibrinogenolysis confirmed the main generation of fibrinopeptide A. Moojase also induced the fibrinolysis of fibrin clots formed in vitro, and the aggregation of washed platelets, as well as significant amidolytic activity on substrates for thrombin, plasma kallikrein, factor Xia, and factor XIIa. Furthermore, thermofluor analyses and the esterase activity of Moojase demonstrated its very high stability at different pH buffers and temperatures. Thus, studies such as this for Moojase should increase knowledge on SVSPs, allowing their bioprospection as valuable prototypes in the development of new drugs, or as biotechnological tools.

Evaluating the Performance of De Novo Assembly Methods for Venom-Gland Transcriptomics

Venom-gland transcriptomics is a key tool in the study of the evolution, ecology, function, and pharmacology of animal venoms. In particular, gene-expression variation and coding sequences gained through transcriptomics provide key information for explaining functional venom variation over both ecological and evolutionary timescales. The accuracy and usefulness of inferences made through transcriptomics, however, is limited by the accuracy of the transcriptome assembly, which is a bioinformatic problem with several possible solutions. Several methods have been employed to assemble venom-gland transcriptomes, with the Trinity assembler being the most commonly applied among them. Although previous evidence of variation in performance among assembly software exists, particularly regarding recovery of difficult-to-assemble multigene families such as snake venom metalloproteinases, much work to date still employs a single assembly method. We evaluated the performance of several commonly used de novo assembly methods for the recovery of both nontoxin transcripts and complete, high-quality venom-gene transcripts across eleven snake and four scorpion transcriptomes. We varied k-mer sizes used by some assemblers to evaluate the impact of k-mer length on transcript recovery. We showed that the recovery of nontoxin transcripts and toxin transcripts is best accomplished through different assembly software, with SDT at smaller k-mer lengths and Trinity being best for nontoxin recovery and a combination of SeqMan NGen and a seed-and-extend approach implemented in Extender as the best means of recovering a complete set of toxin transcripts. In particular, Extender was the only means tested capable of assembling multiple isoforms of the diverse snake venom metalloproteinase family, while traditional approaches such as Trinity recovered at most one metalloproteinase transcript. Our work demonstrated that traditional metrics of assembly performance are not predictive of performance in the recovery of complete and high quality toxin genes. Instead, effective venom-gland transcriptomic studies should combine and quality-filter the results of several assemblers with varying algorithmic strategies.

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.