Differential venom gland gene expression analysis of juvenile and adult scorpions Androctonus crassicauda

ScorpDb: A Novel Open-Access Database for Integrative Scorpion Toxinology

Scorpion stings are a significant public health concern globally, particularly in tropical and subtropical regions. Scorpion venoms contain a diverse array of bioactive peptides, and different scorpion species around the world typically exhibit varying venom profiles, resulting in a wide range of envenomation symptoms. Despite their harmful effects, scorpion venom peptides hold immense potential for drug development due to their unique characteristics. Therefore, the establishment of a comprehensive database that catalogs scorpions along with their known venom peptides and proteins is imperative in furthering research efforts in this research area. We hereby present ScorpDb, a novel database that offers convenient access to data related to different scorpion species, the peptides and proteins found in their venoms, and the symptoms they can cause. To this end, the ScorpDb database has been primarily advanced to accommodate data on the Iranian scorpion fauna. From there, we propose future community efforts to include a larger diversity of scorpions and scorpion venom components. ScorpDb holds the promise to become a valuable resource for different professionals from a variety of research fields, like toxinologists, arachnologists, and pharmacologists. The database is available at https://www.scorpdb.com/.

First transcriptome analysis of the venom glands of the scorpion Hottentotta zagrosensis (Scorpions: Buthidae) with focus on venom lipolysis activating peptides

Introduction

Scorpion venom is a rich source of biological active peptides and proteins. Transcriptome analysis of the venom gland provides detailed insights about peptide and protein venom components. Following the transcriptome analysis of different species in our previous studies, our research team has focused on the Hottentotta zagrosensis as one of the endemic scorpions of Iran to obtain information about its venom proteins, in order to develop biological research focusing on medicinal applications of scorpion venom components and antivenom production. To gain insights into the protein composition of this scorpion venom, we performed transcriptomic analysis.

Methods

Transcriptomic analysis of the venom gland of H. zagrosensis, prepared from the Khuzestan province, was performed through Illumina paired-end sequencing (RNA-Seq), Trinity de novo assembly, CD-Hit-EST clustering, and annotation of identified primary structures using bioinformatics approaches.

Results

Transcriptome analysis showed the presence of 96.4% of complete arthropod BUSCOs, indicating a high-quality assembly. From total of 45,795,108 paired-end 150 bp trimmed reads, the clustering step resulted in the generation of 101,180 de novo assembled transcripts with N50 size of 1,149 bp. 96,071 Unigenes and 131,235 transcripts had a significant similarity (E-value 1e-3) with known proteins from UniProt, Swissprot, Animal toxin annotation project, and the Pfam database. The results were validated using InterProScan. These mainly correspond to ion channel inhibitors, metalloproteinases, neurotoxins, protease inhibitors, protease activators, Cysteine-rich secretory proteins, phospholipase A enzymes, antimicrobial peptides, growth factors, lipolysis-activating peptides, hyaluronidase, and, phospholipase D. Our venom gland transcriptomic approach identified several biologically active peptides including five LVP1-alpha and LVP1-beta isoforms, which we named HzLVP1_alpha1, HzLVP1_alpha2, HzLVP1_alpha3, HzLVP1_beta1, and HzLVP1_beta and have extremely characterized here.

Discussion

Except for HzLVP1_beta1, all other identified LVP1s are predicted to be stable proteins (instability index <40). Moreover, all isoform of LVP1s alpha and beta subunits are thermostable, with the most stability for HzLVP1_alpha2 (aliphatic index = 71.38). HzLVP1_alpha2 has also the highest half-life. Three-dimensional structure of all identified proteins compacts with three disulfide bridges. The extra cysteine residue may allow the proteins to form a hetero- or homodimer. LVP1 subunits of H. zagrosensis potentially interact with adipose triglyceride lipase (ATGL) and hormone-sensitive lipase (HSL), two key enzymes in regulation of lipolysis in adipocytes, suggesting pharmacological properties of these identified proteins.

Dataset of PLA2 family identified from transcriptomic high-throughput sequencing of Androctonus crassicauda (Scorpionida: Buthidae) venom gland

Recently, RNA sequencing has been widely applied to deeply understand the molecular diversity of the venom compounds of various venomous animal species, including scorpions. Among the venomous scorpion species of the Buthidae family, there are many documents of stinging and severe envenoming of victims by the scorpion of Androctonus crassicauda. We present here a high-throughput RNA sequencing dataset of the venom glands from five A. crassicauda individuals, including male and female scorpions. Furthermore, the assembled data corresponding to annotated PLA2 transcripts are also presented. The dataset in this report is related to our research article entitled: "Whole transcriptome sequencing reveals the activity of the PLA2 family members in Androctonus crassicauda (Scorpionida: Buthidae) venom gland" [1]. Here, the venom gland transcriptome analysis of the A. crassicauda was performed. The analysis of concatenated clustered transcriptome assembly using TrinityStats.pl showed that de novo assembly of 517,799,704 clean read pairs generated 744,804 trinity transcripts representing 563,526 trinity genes. BUSCO score for the concatenated clustered transcriptome assembly against orthologs from Arachnida showed 96.7 % complete, 1.6 % fragmented, 1.7 % missing genes, and 2934 genes. Subsequently, the sequences represented PLA2 annotation were extracted from the transcriptome dataset using BLAST searches against the local PLA2 database. We found several cDNA sequences representing PLA2 annotations, which based on sequence similarity to previously found PLA2s, we named platelet-activating factor acetylhydrolases, calcium-dependent PLA2s, calcium-independent PLA2s, and secreted PLA2s. The PLA2 data significantly enrich KEGG pathways related to lipid metabolism. This manuscript complements the primary research article by providing additional data on the abundant estimation of PLA2s.

Venomics of Scorpion Ananteris platnicki (Lourenço, 1993), a New World Buthid That Inhabits Costa Rica and Panama

Ananteris is a scorpion genus that inhabits dry and seasonal areas of South and Central America. It is located in a distinctive morpho-group of Buthids, the 'Ananteris group', which also includes species distributed in the Old World. Because of the lack of information on venom composition, the study of Ananteris species could have biological and medical relevance. We conducted a venomics analysis of Ananteris platnicki, a tiny scorpion that inhabits Panama and Costa Rica, which shows the presence of putative toxins targeting ion channels, as well as proteins with similarity to hyaluronidases, proteinases, phospholipases A2, members of the CAP-domain family, and hemocyanins, among others. Venom proteolytic and hyaluronidase activities were corroborated. The determination of the primary sequences carried out by mass spectrometry evidences that several peptides are similar to the toxins present in venoms from Old World scorpion genera such as Mesobuthus, Lychas, and Isometrus, but others present in Tityus and Centruroides toxins. Even when this venom displays the characteristic protein families found in all Buthids, with a predominance of putative Na+-channel toxins and proteinases, some identified partial sequences are not common in venoms of the New World species, suggesting its differentiation into a distinctive group separated from other Buthids.

Whole transcriptome sequencing reveals the activity of the PLA2 family members in Androctonus crassicauda (Scorpionida: Buthidae) venom gland

Phospholipase A2 is the most abundant venom gland enzyme, whose activity leads to the activation of the inflammatory response by accumulating lipid mediators. This study aimed to identify, classify, and investigate the properties of venom PLA2 isoforms. Then, the present findings were confirmed by chemically measuring the activity of PLA2. The sequences representing PLA2 annotation were extracted from the Androctonus crassicauda transcriptome dataset using BLAS searches against the local PLA2 database. We found several cDNA sequences of PLA2 classified and named by conducting multiple searches as platelet-activating factor acetylhydrolases, calcium-dependent PLA2s, calcium-independent PLA2s, and secreted PLA2s. The largest and smallest isoforms of these proteins range between approximately 70.34 kDa (iPLA2) and 17.75 kDa (cPLA2). Among sPLA2 isoforms, sPLA2GXIIA and sPLA2G3 with ORF encoding 169 and 299 amino acids are the smallest and largest secreted PLA2, respectively. These results collectively suggested that A. crassicauda venom has PLA2 activity, and the members of this protein family may have important biological roles in lipid metabolism. This study also revealed the interaction between members of PLA2s in the PPI network. The results of this study would greatly help with the classification, evolutionary relationships, and interactions between PLA2 family proteins in the gene network.

Identification, classification, and characterization of alpha and beta subunits of LVP1 protein from the venom gland of four Iranian scorpion species

Scorpion venoms contain bioactive peptides and proteins. Some, can be used for pharmaceutical purposes. So, identification of venom proteins matters because, in addition to determining the function of the toxins can also be an excellent guide to developing new drugs. Here, we got transcriptome of venom glands from four Iranian scorpion species, including Hemsicorpius lepturus, Mesobuthus eupeus, Andructunus crassicuada, and Hottentotta saulcyi using cDNA library synthesis and high-throughput transcriptomic analysis of the venom glands. In a comparative way, we identified the cDNA encoding isoforms of subunits (alpha and beta) of BotLVP1/BmLVP1-like protein in the venom gland of three species except for H. lepturus. Characterization and structure determination of the LVP1_like proteins combined with gene map analysis provided evidence of the existence of some isoforms of LVP1_like proteins, encoded by genes with two exons and one intron, which can be classified in CSαβ superfamily in the venom gland of three Iranian scorpion species. According to the high similarity with BotLVP1 and BmLVP1, these proteins could also be potent to mediate cholesterol homeostasis. However, further research is needed to prove it, and this study just may lay the foundation lead to light up this way.

Genome-wide identification, structural homology analysis, and evolutionary diversification of the phospholipase D gene family in the venom gland of three scorpion species

BACKGROUND

Venom phospholipase D (PLDs), dermonecrotic toxins like, are the major molecules in the crude venom of scorpions, which are mainly responsible for lethality and dermonecrotic lesions during scorpion envenoming. The purpose of this study was fivefold: First, to identify transcripts coding for venom PLDs by transcriptomic analysis of the venom glands from Androctonus crassicauda, Hottentotta saulcyi, and Hemiscorpius lepturus; second, to classify them by sequence similarity to known PLDs and motif extraction method; third, to characterize scorpion PLDs; fourth to structural homology analysis with known dermonecrotic toxins; and fifth to investigate phylogenetic relationships of the PLD proteins.

RESULTS

We found that the venom gland of scorpions encodes two PLD isoforms: PLD1 ScoTox-beta and PLD2 ScoTox-alpha I. Two highly conserved regions shared by all PLD1s beta are GAN and HPCDC (HX2PCDC), and the most important conserved regions shared by all PLD2s alpha are two copies of the HKDG (HxKx4Dx6G) motif. We found that PLD1 beta is a 31-43 kDa acidic protein containing signal sequences, and PLD2 alpha is a 128 kDa basic protein without known signal sequences. The gene structures of PLD1 beta and PLD2 alpha contain 6 and 21 exons, respectively. Significant structural homology and similarities were found between the modeled PLD1 ScoTox-beta and the crystal structure of dermonecrotic toxins from Loxosceles intermedia.

CONCLUSIONS

This is the first report on identifying PLDs from A. crassicauda and H. saulcyi venom glands. Our work provides valuable insights into the diversity of scorpion PLD genes and could be helpful in future studies on recombinant antivenoms production.

Proteomic analyses of venom from a Spider Hawk, Pepsis decorata

Background

The composition of the venom from solitary wasps is poorly known, although these animals are considered sources of bioactive substances. Until the present moment, there is only one proteomic characterization of the venom of wasps of the family Pompilidae and this is the first proteomic characterization for the genus Pepsis.

Methods

To elucidate the components of Pepsis decorata venom, the present work sought to identify proteins using four different experimental conditions, namely: (A) crude venom; (B) reduced and alkylated venom; (C) trypsin-digested reduced and alkylated venom, and; (D) chymotrypsin-digested reduced and alkylated venom. Furthermore, three different mass spectrometers were used (Ion Trap-Time of Flight, Quadrupole-Time of Flight, and Linear Triple Quadruple).

Results

Proteomics analysis revealed the existence of different enzymes related to the insect's physiology in the venom composition. Besides toxins, angiotensin-converting enzyme (ACE), hyaluronidase, and Kunitz-type inhibitors were also identified.

Conclusion

The data showed that the venom of Pepsis decorata is mostly composed of proteins involved in the metabolism of arthropods, as occurs in parasitic wasps, although some classical toxins were recorded, and among them, for the first time, ACE was found in the venom of solitary wasps. This integrative approach expanded the range of compounds identified in protein analyses, proving to be efficient in the proteomic characterization of little-known species. It is our understanding that the current work will provide a solid base for future studies dealing with other Hymenoptera venoms.