Venom-gland transcriptomics and venom proteomics of the Hentz striped scorpion (Centruroides hentzi; Buthidae) reveal high toxin diversity in a harmless member of a lethal family

Revisiting toxins with transcriptomics-informed proteomics of venom glands and crude venom from Centruroides bicolor from Panama

The sting of the scorpion Centruroides bicolor causes a large morbidity in Panama. To characterize its venom, transcriptomic and proteomic analyses of the venom glands and the crude venom were performed. These two approaches utilized high-throughput sequencing to enhance the likelihood of detecting a wide range of venom proteins correlated with the venom proteome. After RNA venom gland extraction, a cDNA library was constructed and sequenced by RNA-seq. Also, the crude venom was digested using trypsin and chymotrypsin, and the resulting peptides were analyzed using a nano-LC-MS/MS. Notably, transcriptomic and proteomic venom approaches identified a hyaluronidase, alpha- and beta-neurotoxins that affect Na+ channels, CRISP proteins, metalloproteinases, transferrin, monooxygenase alpha-peptidyl-glycine, serine proteases, alpha pancreatic amylase, lysozyme, neurotoxins targeting K+ channels, neprilysin, scorpine, angiotensin-converting enzyme, insulin-like growth factor-binding domain proteins, nucleobindin-like proteins, and uncharacterized proteins. Interestingly, some of the venom proteins such as nucleobindin and angiotensin-converting enzymes have been not reported in the proteome, their predicted presence has only been previously derived from the genomic sequence of Centruroides sculpturatus and C. vittatus. These newly identified components enhance the understanding of the venomous nature of C. bicolor. SIGNIFICANCE: The proteins and peptides found in Centruroides bicolor venom by transcriptomic and proteomic analyses were assessed according to the protein and toxin databases available on public domains. Notably, some of the venom proteins such as nucleobindin and angiotensin-converting enzymes have been not reported in the proteome, their predicted presence has only been previously derived from the genomic sequence of Centruroides sculpturatus and C. vittatus. Moreover, enzymatic assays, including hyaluronidase, phospholipase A2, and proteolytic activity were conducted to confirm the presence or absence of those enzymes. Interestingly, neurotoxins from C. limbatus, a related species in the region, were found in the proteome but no mRNAs were identified in the transcriptome. These newly identified components enhance the understanding of the venomous nature of Centruroides bicolor.

Optimizing Scorpion Toxin Processing through Artificial Intelligence

Scorpion toxins are relatively short cyclic peptides (<150 amino acids) that can disrupt the opening/closing mechanisms in cell ion channels. These peptides are widely studied for several reasons including their use in drug discovery. Although improvements in RNAseq have greatly expedited the discovery of new scorpion toxins, their annotation remains challenging, mainly due to their small size. Here, we present a new pipeline to annotate toxins from scorpion transcriptomes using a neural network approach. This pipeline implements basic neural networks to sort amino acid sequences to find those that are likely toxins and thereafter predict the type of toxin represented by the sequence. We anticipate that this pipeline will accelerate the classification of scorpion toxins in forthcoming scorpion genome sequencing projects and potentially serve a useful role in identifying targets for drug development.

Taking the sting out of scorpions: Electrophysiological investigation of the relative efficacy of three antivenoms against medically significant Centruroides species

In this study, we report the innovative application of whole-cell patch-clamp electrophysiology in assessing broad-spectrum neutralisation by three different antivenoms, of venoms from the medically significant scorpion genus Centruroides. Envenomations by as many as 21 species from the Centruroides genus result in up to 300,000 envenomations per year in Mexico, which poses significant and potentially life-threatening pathophysiology. We first evaluated the in vitro manifestation of envenomation against two human voltage-gated sodium (hNaV) channel subtypes: hNaV1.4 and hNaV1.5, which are primarily expressed in skeletal muscles and cardiomyocytes, respectively. The neutralisation of venom activity was then characterised for three different antivenoms using a direct competition model against the more potent target, hNaV1.4. While broad-spectrum neutralisation was identified, variation in neutralisation arose for Centruroides elegans, C. limpidus, C. noxius and C. suffusus venoms, despite the presence of a number of these venoms within the immunising mixture. This raises questions regarding the truly "broad" neutralisation capacity of the antivenoms. This study not only extends previous validation of the in vitro investigation of antivenom efficacy utilising the whole-cell patch-clamp technique but also underscores the potential of this animal-free model in exploring cross-reactivity, experimental scalability, and most importantly, informing clinical management practices regarding the administration of antivenom in Mexico.

Varying Modes of Selection Among Toxin Families in the Venoms of the Giant Desert Hairy Scorpions (Hadrurus)

Venoms are primarily believed to evolve under strong diversifying selection resulting from persistent coevolution between predator and prey. Recent research has challenged this hypothesis, proposing that venoms from younger venomous lineages (e.g., snakes and cone snails) are governed predominantly by diversifying selection, while venoms from older venomous lineages (e.g., centipedes, scorpions, and spiders) are under stronger purifying selection. However, most research in older lineages has tested selection at more diverse phylogenetic scales. Although these tests are important for evaluating broad macroevolutionary trends underlying venom evolution, they are less equipped to detect species-level evolutionary trends, which likely have large impacts on venom variation seen at more diverse phylogenetic scales. To test for selection among closely related species from an older venomous lineage, we generated high-throughput venom-gland transcriptomes and venom proteomes for four populations of Giant Desert Hairy Scorpions (Hadrurus), including three Hadrurus arizonensis populations and one Hadrurus spadix population. We detected significant episodic and pervasive diversifying selection across a highly abundant toxin family that likely has a major role in venom function ([Formula: see text]KTxs), providing a contrast to the stronger purifying selection identified from other studies on scorpion venoms. Conversely, we detected weak episodic diversifying and/or stronger purifying selection in four toxin families (non-disulfide bridged peptides, phospholipase A2s, scorpine-like antimicrobial peptides, and serine proteases), most of which were less abundant and likely have ancillary functional roles. Finally, although we detected several major toxin families at disproportionate transcriptomic and/or proteomic abundances, we did not identify significant sex-based variation in Hadrurus venoms.

The arylsulfatase- and phospholipase-rich venom of the plutoniumid centipede Theatops posticus

Research on centipede venoms has led to the discovery of a diverse array of novel proteins and peptides, including those with homology to previously discovered toxin families (e.g., phospholipase A2s and pM12a metalloproteases) and novel toxin families not previously detected in venoms (e.g., β-pore forming toxins and scoloptoxins). Most of this research has focused on centipedes in the order Scolopendromorpha, particularly those in the families Scolopendridae, Cryptopidae, and Scolopocryptopidae. To generate the first high-throughput venom characterization for a centipede in the scolopendromorph family Plutoniumidae, we performed venom-gland transcriptomics and venom proteomics on two Theatops posticus. We identified a total of 64 venom toxins, 60 of which were detected in both the venom-gland transcriptome and venom proteome and four of which were only detected transcriptomically. We detected a single highly abundant arylsulfatase B (ARSB) toxin, the first ARSB toxin identified from centipede venoms. As ARSBs have been detected in other venomous species (e.g., scorpions), ARSBs in T. posticus highlights a new case of convergent evolution across venoms. Theatops posticus venom also contained a much higher abundance and diversity of phospholipase A2 toxins compared to other characterized centipede venoms. Conversely, we detected other common centipedes toxins, such as CAPs and scoloptoxins, at relatively low abundances and diversities. Our observation of a diverse set of toxins from T. posticus venom, including those from novel toxin families, emphasizes the importance of studying unexplored centipede taxonomic groups and the continued potential of centipede venoms for novel toxin discovery and unraveling the molecular mechanisms underlying trait evolution.

Scorpion Venom as a Source of Antimicrobial Peptides: Overview of Biomolecule Separation, Analysis and Characterization Methods

Scorpion venoms have long captivated scientific researchers, primarily due to the potency and specificity of the mechanism of action of their derived components. Among other molecules, these venoms contain highly active compounds, including antimicrobial peptides (AMPs) and ion channel-specific components that selectively target biological receptors with remarkable affinity. Some of these receptors have emerged as prime therapeutic targets for addressing various human pathologies, including cancer and infectious diseases, and have served as models for designing novel drugs. Consequently, extensive biochemical and proteomic investigations have focused on characterizing scorpion venoms. This review provides a comprehensive overview of the key methodologies used in the extraction, purification, analysis, and characterization of AMPs and other bioactive molecules present in scorpion venoms. Noteworthy techniques such as gel electrophoresis, reverse-phase high-performance liquid chromatography, size exclusion chromatography, and "omics" approaches are explored, along with various combinations of methods that enable bioassay-guided venom fractionation. Furthermore, this review presents four adapted proteomic workflows that lead to the comprehensive dissection of the scorpion venom proteome, with an emphasis on AMPs. These workflows differ based on whether the venom is pre-fractionated using separation techniques or is proteolytically digested directly before further proteomic analyses. Since the composition and functionality of scorpion venoms are species-specific, the selection and sequence of the techniques for venom analyses, including these workflows, should be tailored to the specific parameters of the study.

Unveiling the Protein Components of the Secretory-Venom Gland and Venom of the Scorpion Centruroides possanii (Buthidae) through Omic Technologies

Centruroides possanii is a recently discovered species of "striped scorpion" found in Mexico. Certain species of Centruroides are known to be toxic to mammals, leading to numerous cases of human intoxications in the country. Venom components are thought to possess therapeutic potential and/or biotechnological applications. Hence, obtaining and analyzing the secretory gland transcriptome and venom proteome of C. possanii is relevant, and that is what is described in this communication. Since this is a newly described species, first, its LD50 to mice was determined and estimated to be 659 ng/g mouse weight. Using RNA extracted from this species and preparing their corresponding cDNA fragments, a transcriptome analysis was obtained on a Genome Analyzer (Illumina) using the 76-base pair-end sequencing protocol. Via high-throughput sequencing, 19,158,736 reads were obtained and ensembled in 835,204 sequences. Of them, 28,399 transcripts were annotated with Pfam. A total of 244 complete transcripts were identified in the transcriptome of C. possanii. Of these, 109 sequences showed identity to toxins that act on ion channels, 47 enzymes, 17 protease inhibitors (PINs), 11 defense peptides (HDPs), and 60 in other components. In addition, a sample of the soluble venom obtained from this scorpion was analyzed using an Orbitrap Velos apparatus, which allowed for identification by liquid chromatography followed by mass spectrometry (LC-MS/MS) of 70 peptides and proteins: 23 toxins, 27 enzymes, 6 PINs, 3 HDPs, and 11 other components. Until now, this work has the highest number of scorpion venom components identified through omics technologies. The main novel findings described here were analyzed in comparison with the known data from the literature, and this process permitted some new insights in this field.

Venom phenotype conservation suggests integrated specialization in a lizard-eating snake

Biological specialization reduces the size of niche space while increasing efficiency in the use of available resources. Specialization often leads to phenotypic changes via natural selection aligning with niche space constraints. Commonly observed changes are in size, shape, behavior, and traits associated with feeding. One often selected trait for dietary specialization is venom, which, in snakes, often shows variation dependent on diet across and within species. The Neotropical Blunt-headed Treesnake (Imantodes cenchoa) is a highly specialized, rear-fanged, arboreal, lizard hunter that displays a long thin body, enlarged eyes, and a large Duvernoy's gland. However, toxin characterization of I. cenchoa has never been completed. Here, we use RNA-seq and mass spectrometry to assemble, annotate, and analyze the venom gland transcriptomes of four I. cenchoa from across their range. We find a lack of significant venom variation at the sequence and expression levels, suggesting venom conservation across the species. We propose this conservation provides evidence of a specialized venom repertoire, adapted to maximize efficiency of capturing and processing lizards. Importantly, this study provides the most complete venom gland transcriptomes of I. cenchoa and evidence of venom specialization in a rear-fanged snake, giving insight into selective pressures of venom across all snake species.

The phylogenetic relationship among two species of genus Nebo (Scorpiones: Diplocentridae) from Saudi Arabia and Middle East

BACKGROUND

The genus Nebo has been identified as a medically important scorpion species distributed across Arabia and the Middle East. However, its taxonomic status remains unclear.

AIM

The molecular phylogeny of two Nebo species from Saudi Arabia and comparative sequences from Palestine is presented based on the mitochondrial cytochrome oxidase subunit I (COI) gene.

METHODOLOGY

Scorpion specimens were collected from two different localities, mainly the Southern part of Saudi Arabia. Then, DNA was extracted, amplified using invertebrate universal primers, and sequenced to identify the COI gene. The obtained sequences were analyzed, and phylogenetic trees based on maximum parsimony, neighbor-joining, and Bayesian inference were constructed.

RESULTS

The inferred phylogeny indicates the monophyletic status of the family Diplocentridae and its subfamily Nebinae and Diplocentrinae. Also, the phylogenetic analyses support the existence of interspecific and intraspecific variations among/ within Nebo hierichonticus and Nebo yemenensis which may indicate distinct species.

CONCLUSION

Further morphological studies with additional specimens from the Arabian Peninsula may reveal possible undiscovered and cryptic species.

The remarkably enzyme-rich venom of the Big Bend Scorpion (Diplocentrus whitei)

Scorpion venoms have long been studied for their peptide discovery potential, with modern high-throughput venom-characterization techniques paving the way for the discovery of thousands of novel putative toxins. Research into these toxins has provided insight into the pathology and treatment of human diseases, even resulting in the development of one compound with Food and Drug Administration (FDA) approval. Although most of this research has focused on the toxins of scorpion species considered medically significant to humans, the venom of harmless scorpion species possess toxins that are homologous to those from medically significant species, indicating that harmless scorpion venoms may also serve as valuable sources of novel peptide variants. Furthermore, as harmless scorpions represent a vast majority of scorpion species diversity, and therefore venom toxin diversity, venoms from these species likely contain entirely new toxin classes. We sequenced the venom-gland transcriptome and venom proteome of two male Big Bend scorpions (Diplocentrus whitei), providing the first high-throughput venom characterization for a member of this genus. We identified a total of 82 toxins in the venom of D. whitei, 25 of which were identified in both the transcriptome and proteome, and 57 of which were only identified in the transcriptome. Furthermore, we identified a unique, enzyme-rich venom dominated by serine proteases and the first arylsulfatase B toxins identified in scorpions.

An overview of Tityus cisandinus scorpion venom: Transcriptome and mass fingerprinting reveal conserved toxin homologs across the Amazon region and novel lipolytic components

Tityus cisandinus, a neglected medically important scorpion in Ecuadorian and Peruvian Amazonia, belongs to a complex of species related to the eastern Amazon endemic Tityus obscurus, spanning a distribution of ca. 4000 km. Despite high morbidity and mortality rates, no effective scorpion antivenom is currently available in the Amazon region. Knowledge of the structural/functional relationships between T. cisandinus venom components and those from related Amazonian species is crucial for designing region-specific therapeutic antivenoms. In this work, we carried out the first venom gland transcriptomic study of an Amazonian scorpion outside Brazil, T. cisandinus. We also fingerprinted its total venom through MALDI-TOF MS, which supported our transcriptomic findings. We identified and calculated the expression level of 94 components: 60 toxins, 25 metalloproteases, five disulfide isomerases, three amidating enzymes, one hyaluronidase, and also uncovered transcripts encoding novel lipolytic beta subunits produced by New World buthid scorpions. This study demonstrates the high similarity between T. cisandinus and T. obscurus venoms, reinforcing the existence of a neglected complex of genetically and toxinologically related Amazonian scorpions of medical importance. Finally, we demonstrated the low recognition of currently available therapeutic sera against T. cisandinus and T. obscurus venoms, and concluded that these should be improved to protect against envenomation by Amazonian Tityus spp.

French Scorpionism (Mainland and Oversea Territories): Narrative Review of Scorpion Species, Scorpion Venom, and Envenoming Management

Sixty-seven scorpion species have been described in France and its territories, where they have been found to be heterogeneously distributed. Indeed, only one species can be found on Réunion Island, while 38 species exist in French Guiana. The number of stings is also heterogenous, with up to 90 stings per 100,000 inhabitants occurring annually. Scorpion species can frequently be determined through simple visual factors, including species of medical importance (i.e., Buthus, Centruroides and Tityus). Scorpion venom is composed of local enzymes and peptides with a cysteine-stabilized α/β motif (NaTxs, Ktxs, Calcines), which allow for venom diffusion and the prey's incapacitation, respectively. Harmful scorpion species are limited to Centruroides pococki in the French West Indies, which can induce severe envenoming, and the Tityus obscurus and Tityus silvestris in French Guiana, which can cause fatalities in children and can induce severe envenoming, respectively. Envenomation by one of these scorpions requires hospital monitoring as long as systemic symptoms persist. Typical management includes the use of a lidocaine patch, pain killers, and local antiseptic. In the case of heart failure, the use of dobutamine can improve survival, and pregnant women must consult an obstetrician because of the elevated risk of preterm birth or stillbirth. France does not have scorpion antivenom, as scorpion stings are generally not fatal.

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

BACKGROUND

The Androctonus crassicauda, belonging to the genus Androctonus of the family Buthidae, is the most venomous scorpion in Middle East countries. However, the venom gland transcriptome profile of A. crassicauda scorpion has not yet been studied. In this study, we elucidated and compared the venom gland gene expression profiles of adult and juvenile male scorpion A. crassicauda using high-throughput transcriptome sequencing. This is the first report of transcriptional analysis of the venom glands of scorpions in different growth stages, with insights into the identification of the key genes during venom gland development.

RESULTS

A total of 209,951 mRNA transcripts were identified from total RNA-seq data, of which 963 transcripts were differentially expressed (DE) in adult and juvenile scorpions (p < 0.01). Overall, we identified 558 up-regulated and 405 down-regulated transcripts in the adult compared to the juvenile scorpions, of which 397 and 269 unique unigenes were annotated, respectively. GO and KEGG enrichment analyses indicated that the metabolic, thermogenesis, cytoskeleton, estrogen signaling, GnRH signaling, growth hormone signaling, and melanogenesis pathways were affected by two different growth conditions and the results suggested that the DE genes related to those pathways are important genes associated with scorpion venom gland development, in which they may be important in future studies, including Chs, Elovl, MYH, RDX, ACTN, VCL, PIP5K, PP1C, FGFR, GNAS, EGFR, CREB, CoA, PLCB, CALM, CACNA, PKA and CAMK genes.

CONCLUSIONS

These findings broadened our knowledge of the differences between adult and juvenile scorpion venom and opened new perspectives on the application of comparative transcriptome analysis to identify the special key genes.

Evolution, expression patterns and distribution of novel ribbon worm predatory and defensive toxins

Ribbon worms are active predators that use an eversible proboscis to inject venom into their prey and defend themselves with toxic epidermal secretions. Previous work on nemertean venom has largely focused on just a few species and has not investigated the different predatory and defensive secretions in detail. Consequently, our understanding of the composition and evolution of ribbon worm venoms is still very limited. Here, we present a comparative study of nemertean venom combining RNA-seq differential gene expression analyses of venom-producing tissues, tandem mass spectrometry-based proteomics of toxic secretions, and mass spectrometry imaging of proboscis sections, to shed light onto the composition and evolution of predatory and defensive toxic secretions in Antarctonemertes valida. Our analyses reveal a wide diversity of putative defensive and predatory toxins with tissue-specific gene expression patterns and restricted distributions to the mucus and proboscis proteomes respectively, suggesting that ribbon worms produce distinct toxin cocktails for predation and defense. Our results also highlight the presence of numerous lineage-specific toxins, indicating that venom evolution is highly divergent across nemerteans, producing toxin cocktails that might be finely tuned to subdue different prey. Our data also suggest that the hoplonemertean proboscis is a highly specialized predatory organ that seems to be involved in a variety of biological functions besides predation, including secretion and sensory perception. Overall, our results advance our knowledge into the diversity and evolution of nemertean venoms and highlight the importance of combining different types of data to characterize toxin composition in understudied venomous organisms.

The Enzymatic Core of Scorpion Venoms

Enzymes are an integral part of animal venoms. Unlike snakes, in which enzymes play a primary role in envenomation, in scorpions, their function appears to be ancillary in most species. Due to this, studies on the diversity of scorpion venom components have focused primarily on the peptides responsible for envenomation (toxins) and a few others (e.g., antimicrobials), while enzymes have been overlooked. In this work, a comprehensive study on enzyme diversity in scorpion venoms was performed by transcriptomic and proteomic techniques. Enzymes of 63 different EC types were found, belonging to 330 orthogroups. Of them, 24 ECs conform the scorpion venom enzymatic core, since they were determined to be present in all the studied scorpion species. Transferases and lyases are reported for the first time. Novel enzymes, which can play different roles in the venom, including direct toxicity, as venom spreading factors, activators of venom components, venom preservatives, or in prey pre-digestion, were described and annotated. The expression profile for transcripts coding for venom enzymes was analyzed, and shown to be similar among the studied species, while being significantly different from their expression pattern outside the telson.

Snake Venomics: Fundamentals, Recent Updates, and a Look to the Next Decade

Venomic research, powered by techniques adapted from proteomics, transcriptomics, and genomics, seeks to unravel the diversity and complexity of venom through which knowledge can be applied in the treatment of envenoming, biodiscovery, and conservation. Snake venom proteomics is most extensively studied, but the methods varied widely, creating a massive amount of information which complicates data comparison and interpretation. Advancement in mass spectrometry technology, accompanied by growing databases and sophisticated bioinformatic tools, has overcome earlier limitations of protein identification. The progress, however, remains challenged by limited accessibility to samples, non-standardized quantitative methods, and biased interpretation of -omic data. Next-generation sequencing (NGS) technologies enable high-throughput venom-gland transcriptomics and genomics, complementing venom proteomics by providing deeper insights into the structural diversity, differential expression, regulation and functional interaction of the toxin genes. Venomic tissue sampling is, however, difficult due to strict regulations on wildlife use and transfer of biological materials in some countries. Limited resources for techniques and funding are among other pertinent issues that impede the progress of venomics, particularly in less developed regions and for neglected species. Genuine collaboration between international researchers, due recognition of regional experts by global organizations (e.g., WHO), and improved distribution of research support, should be embraced.

Contrasting patterns of venom regeneration in a centipede (Scolopendra viridis) and a scorpion (Centruroides hentzi)

As biochemical traits with clear fitness consequences, venoms serve a critical ecological role for the animals that produce them. Understanding how venoms are maintained and regenerated after use will, therefore, provide valuable insight into the ecology of venomous animals. Furthermore, most studies on venomous organisms often require removing animals from the wild and waiting extended periods of time between venom extractions. Uncovering the patterns of venom regeneration across different species will likely lead to the development of more efficient venom extraction protocols, reducing both experimental time and the number of animals required. Using reversed-phase high-performance liquid chromatography, we identified asynchronous regeneration of venom protein component abundances in the centipede Scolopendra viridis but found no evidence for asynchronous venom regeneration in the scorpion Centruroides hentzi. We also observed high levels of intraspecific venom variation in C. hentzi, emphasizing the importance of testing for intraspecific venom variation in studies evaluating the synchronicity of venom regeneration. Although the regeneration of relative venom protein component abundances is an asynchronous process in S. viridis, we provide evidence that the presence-absence of major venom components is not an asynchronous process and suggest that studies relying on just the presence/absence of individual proteins (e.g. bioprospecting, drug discovery) could use catch-and-release methods of venom extraction to reduce the number of animals removed from the wild.

Toxinology in the proteomics era: a review on arachnid venom proteomics

The word venomics was coined to acknowledge the studies that use omics to investigate venom proteins and peptides. Venomics has evolved considerably over the last 20 years. The first works on scorpion or spider venomics were published in the early 2000's. Such studies relied on peptide mass fingerprinting (PMF) to characterize venom complexity. After the introduction of new mass spectrometers with higher resolution, sensitivity and mass accuracy, and the next-generation nucleotide sequencing, the complexity of data reported in research on scorpion and spider venomics increased exponentially, which allowed more comprehensive studies. In the present review article, we covered key publications on scorpion venomics and spider venomics, presenting historical grounds and implemented technologies over the last years. The literature presented in this review was selected after searching the PubMed database using the terms "(scorpion venom) AND (proteome)" for scorpion venomics, and "(spider venom) AND (proteome)" for publications on spider venomics. We presented the key aspects related to proteomics in the covered papers including, but not restricted to, the employed proteomic strategy (i.e., PMF, two-dimensional gel electrophoresis, shotgun/bottom-up and/or top-down/peptidome), and the type of mass spectrometer used. Some conclusions can be drawn from the present study. For example, the scorpion genus Tityus is the most studied concerning venomics, followed by Centruroides; whereas for spiders the studied genera were found more equally distributed. Another interesting conclusion is the lack of high throughput studies on post-translational modifications (PTMs) of scorpion and spider proteins. In our opinion, PTMs should be more studied as they can modulate the activity of scorpion and spider toxins.

Transcriptomic and proteomic analyses of venom glands from scorpions Liocheles australasiae, Mesobuthus martensii, and Scorpio maurus palmatus

Scorpion venom contains a cocktail of differing peptides and proteins. Previous studies focused on the identification of species-specific components in scorpion venoms, and whether there could be peptides and/or proteins conserved in the venom gland of a scorpion ancestor has been rarely investigated. Here, using a combination of transcriptomic and proteomic approaches, putative conserved toxins from the venom glands of scorpions Liocheles australasiae, Mesobuthus martensii, and Scorpio maurus palmatus were identified and compared. Similar to other studies, more than half of the conserved toxins are predominantly proteins including proteases. On the other hand, unique venom peptides, including ion channel toxins were revealed specifically in the M. martensii. The sodium channel toxin peptides revealed in M. martensii consolidated that scorpions in the Buthidae are able to envenomate their prey wih highly neurotoxic venom. This study suggested that these conserved proteins had already formed part of the arsenal in the venom gland of the common ancestor of scorpions, and likely perform important functional roles in envenomation during scorpion evolution.

Terrestrial venomous animals, the envenomings they cause, and treatment perspectives in the Middle East and North Africa

The Middle East and Northern Africa, collectively known as the MENA region, are inhabited by a plethora of venomous animals that cause up to 420,000 bites and stings each year. To understand the resultant health burden and the key variables affecting it, this review describes the epidemiology of snake, scorpion, and spider envenomings primarily based on heterogenous hospital data in the MENA region and the pathologies associated with their venoms. In addition, we discuss the venom composition and the key medically relevant toxins of these venomous animals, and, finally, the antivenoms that are currently in use to counteract them. Unlike Asia and sub-Saharan Africa, scorpion stings are significantly more common (approximately 350,000 cases/year) than snakebites (approximately 70,000 cases/year) and present the most significant contributor to the overall health burden of envenomings, with spider bites being negligible. However, this review also indicates that there is a substantial lack of high-quality envenoming data available for the MENA region, rendering many of these estimates speculative. Our understanding of the venoms and the toxins they contain is also incomplete, but already presents clear trends. For instance, the majority of snake venoms contain snake venom metalloproteinases, while sodium channel-binding toxins and potassium channel-binding toxins are the scorpion toxins that cause most health-related challenges. There also currently exist a plethora of antivenoms, yet only few are clinically validated, and their high cost and limited availability present a substantial health challenge. Yet, some of the insights presented in this review might help direct future research and policy efforts toward the appropriate prioritization of efforts and aid the development of future therapeutic solutions, such as next-generation antivenoms.

Reduced Toxicity of Centruroides vittatus (Say, 1821) May Result from Lowered Sodium β Toxin Gene Expression and Toxin Protein Production

Body tissue and venom glands from an eastern population of the scorpion Centruroides vittatus (Say, 1821) were homogenized and molecular constituents removed to characterize putative sodium β toxin gene diversity, RT-qPCR, transcriptomic, and proteomic variation. We cloned sodium β toxins from genomic DNA, conducted RT-qPCR experiments with seven sodium β toxin variants, performed venom gland tissue RNA-seq, and isolated venom proteins for mass spectrophotometry. We identified >70 putative novel sodium β toxin genes, 111 toxin gene transcripts, 24 different toxin proteins, and quantified sodium β toxin gene expression variation among individuals and between sexes. Our analyses contribute to the growing evidence that venom toxicity among scorpion taxa and their populations may be associated with toxin gene diversity, specific toxin transcripts variation, and subsequent protein production. Here, slight transcript variation among toxin gene variants may contribute to the major toxin protein variation in individual scorpion venom composition.

A non-lethal method for studying scorpion venom gland transcriptomes, with a review of potentially suitable taxa to which it can be applied

Scorpion venoms are mixtures of proteins, peptides and small molecular compounds with high specificity for ion channels and are therefore considered to be promising candidates in the venoms-to-drugs pipeline. Transcriptomes are important tools for studying the composition and expression of scorpion venom. Unfortunately, studying the venom gland transcriptome traditionally requires sacrificing the animal and therefore is always a single snapshot in time. This paper describes a new way of generating a scorpion venom gland transcriptome without sacrificing the animal, thereby allowing the study of the transcriptome at various time points within a single individual. By comparing these venom-derived transcriptomes to the traditional whole-telson transcriptomes we show that the relative expression levels of the major toxin classes are similar. We further performed a multi-day extraction using our proposed method to show the possibility of doing a multiple time point transcriptome analysis. This allows for the study of patterns of toxin gene activation over time a single individual, and allows assessment of the effects of diet, season and other factors that are known or likely to influence intraindividual venom composition. We discuss the gland characteristics that may allow this method to be successful in scorpions and provide a review of other venomous taxa to which this method may potentially be successfully applied.

Identification and Characterization of Novel Proteins from Arizona Bark Scorpion Venom That Inhibit Nav1.8, a Voltage-Gated Sodium Channel Regulator of Pain Signaling

The voltage-gated sodium channel Nav1.8 is linked to neuropathic and inflammatory pain, highlighting the potential to serve as a drug target. However, the biophysical mechanisms that regulate Nav1.8 activation and inactivation gating are not completely understood. Progress has been hindered by a lack of biochemical tools for examining Nav1.8 gating mechanisms. Arizona bark scorpion (Centruroides sculpturatus) venom proteins inhibit Nav1.8 and block pain in grasshopper mice (Onychomys torridus). These proteins provide tools for examining Nav1.8 structure–activity relationships. To identify proteins that inhibit Nav1.8 activity, venom samples were fractioned using liquid chromatography (reversed-phase and ion exchange). A recombinant Nav1.8 clone expressed in ND7/23 cells was used to identify subfractions that inhibited Nav1.8 Na⁺ current. Mass-spectrometry-based bottom-up proteomic analyses identified unique peptides from inhibitory subfractions. A search of the peptides against the AZ bark scorpion venom gland transcriptome revealed four novel proteins between 40 and 60% conserved with venom proteins from scorpions in four genera (Centruroides, Parabuthus, Androctonus, and Tityus). Ranging from 63 to 82 amino acids, each primary structure includes eight cysteines and a “CXCE” motif, where X = an aromatic residue (tryptophan, tyrosine, or phenylalanine). Electrophysiology data demonstrated that the inhibitory effects of bioactive subfractions can be removed by hyperpolarizing the channels, suggesting that proteins may function as gating modifiers as opposed to pore blockers.

An integrated transcriptomic and proteomic approach to identify the main Torymus sinensis venom components

During oviposition, ectoparasitoid wasps not only inject their eggs but also a complex mixture of proteins and peptides (venom) in order to regulate the host physiology to benefit their progeny. Although several endoparasitoid venom proteins have been identified, little is known about the components of ectoparasitoid venom. To characterize the protein composition of Torymus sinensis Kamijo (Hymenoptera: Torymidae) venom, we used an integrated transcriptomic and proteomic approach and identified 143 venom proteins. Moreover, focusing on venom gland transcriptome, we selected additional 52 transcripts encoding putative venom proteins. As in other parasitoid venoms, hydrolases, including proteases, phosphatases, esterases, and nucleases, constitute the most abundant families in T. sinensis venom, followed by protease inhibitors. These proteins are potentially involved in the complex parasitic syndrome, with different effects on the immune system, physiological processes and development of the host, and contribute to provide nutrients to the parasitoid progeny. Although additional in vivo studies are needed, initial findings offer important information about venom factors and their putative host effects, which are essential to ensure the success of parasitism.

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.

Exploring the biological activities and proteome of Brazilian scorpion Rhopalurus agamemnon venom

Scorpion venoms are formed by toxins harmful to various organisms, including humans. Several techniques have been developed to understand the role of proteins in animal venoms, including proteomics approach. Rhopalurus agamemnon (Koch, 1839) is the largest scorpion in the Buthidae family in the Brazilian Cerrado, measuring up to 110 mm in total length. The accident with R. agamemnon is painful and causes some systemic reactions, but the specie's venom remains uninvestigated. We explore the venom protein composition using a proteomic and a biological-directed approach identifying 230 protein compounds including enzymes like Hyaluronidase, metalloproteinase, L-amino acid oxidase and amylase, the last two are first reported for scorpion venoms. Some of those new reports are important to demonstrate how distant we are from a total comprehension of the diversity about venoms in general, due to their diversity in composition and function. BIOLOGICAL SIGNIFICANCE: In this study, we explored the composition of venom proteins from the scorpion Rhopalurus agamemnon. We identified 230 proteins from the venom including new enzyme reports. These data highlight the unique diversity of the venom proteins from the scorpion R. agamemnon, provide insights into new mechanisms of envenomation and enlarge the protein database of scorpion venoms. The discovery of new proteins provides a new scenario for the development of new drugs and suggests molecular targets to venom components.

De novo transcriptomic and proteomic analysis and potential toxin screening of Mesobuthus martensii samples from four different provinces

ETHNOPHARMACOLOGICAL RELEVANCE

As well-known medicinal materials in traditional Chinese medicine, scorpions, commonly called as Quanxie () in Chinese, have been widely used to treat several diseases such as rheumatoid arthritis, apoplexy, epilepsy and chronic pain for more than a thousand years. Not only in the ancient times, the scorpions have also been recorded nowadays in the Pharmacopoeia of the People's Republic of China since 1963.

AIM OF STUDY

This study aims to explore the differences in composition of the venom of scorpions from different regions by using the method of transcriptomics and proteomics.

MATERIALS AND METHODS

Whole de novo transcriptomes, proteomics and their bioinformatic analyses were performed on samples of the scorpion Mesobuthus martensii and their venoms from four different provinces with clear geographical boundaries, including Hebei, Henan, Shandong and Shanxi.

RESULTS

The four captured samples had the same morphology, and the conserved CO-1 sequence matched that of M. martensii. A total of 141,003 of 174,653 transcripts were identified as unigenes, of which we successfully annotated 51,627 (36.61%), 21,970 (15.58%), 7,168 (5.08%), and 45,263 (32.10%) unigenes with the NR, GO, KEGG and SWISSPROT databases, respectively, while a total of 427 proteins were collected from the protein extracted from venoms. Both GO and KEGG annotations exhibited only slight differences among the four samples while the expression level of gene and protein was quite different. A total of 249 toxin-related unigenes were successfully screened, including 41 serine proteases and serine protease inhibitors, 39 potassium channel toxins, 38 phospholipases, 16 host defense peptides, 9 metalloproteases, and 50 other toxins. Although the toxin species were similar among the four samples, the gene expression of each toxin varied considerably, for example, the scorpion from HB province has the most abundant expression quality in sequences c48391_g1, c55239_g1 and c47749_g1 while the lowest expressions of c51178_g1, c62033_g3 and c63754_g2.

CONCLUSION

The regional differences in the transcriptomes and proteomes of M. martensii are mainly from expression levels e.g. toxins rather than expression species, of which the method can be further extended to evaluate the qualities of traditional Chinese medicines obtained from different regions.

Identification of an antiviral component from the venom of the scorpion Liocheles australasiae using transcriptomic and mass spectrometric analyses

Scorpion venom contains a variety of biologically active peptides. Among them, neurotoxins are major components in the venom, but it also contains peptides that show antimicrobial activity. Previously, we identified three insecticidal peptides from the venom of the Liocheles australasiae scorpion, but activities and structures of other venom components remained unknown. In this study, we performed a transcriptome analysis of the venom gland of the scorpion L. australasiae to gain a comprehensive understanding of its venom components. The result shows that potassium channel toxin-like peptides were the most diverse, whereas only a limited number of sodium channel toxin-like peptides were observed. In addition to these neurotoxin-like peptides, many non-disulfide-bridged peptides were identified, suggesting that these components have some critical roles in the L. australasiae venom. In this study, we also isolated a component with antiviral activity against hepatitis C virus using a bioassay-guided fractionation approach. By integrating mass spectrometric and transcriptomic data, we successfully identified LaPLA₂-1 as an anti-HCV component. LaPLA₂-1 is a phospholipase A₂ having a heterodimeric structure that is N-glycosylated at the N-terminal region. Since the antiviral activity of LaPLA₂-1 was inhibited by a PLA₂ inhibitor, the enzymatic activity of LaPLA₂-1 is likely to be involved in its antiviral activity.

An Economic Dilemma Between Molecular Weapon Systems May Explain an Arachno-atypical Venom in Wasp Spiders (Argiope bruennichi)

Spiders use venom to subdue their prey, but little is known about the diversity of venoms in different spider families. Given the limited data available for orb-weaver spiders (Araneidae), we selected the wasp spider Argiope bruennichi for detailed analysis. Our strategy combined a transcriptomics pipeline based on multiple assemblies with a dual proteomics workflow involving parallel mass spectrometry techniques and electrophoretic profiling. We found that the remarkably simple venom of A. bruennichi has an atypical composition compared to other spider venoms, prominently featuring members of the cysteine-rich secretory protein, antigen 5 and pathogenesis-related protein 1 (CAP) superfamily and other, mostly high-molecular-weight proteins. We also detected a subset of potentially novel toxins similar to neuropeptides. We discuss the potential function of these proteins in the context of the unique hunting behavior of wasp spiders, which rely mostly on silk to trap their prey. We propose that the simplicity of the venom evolved to solve an economic dilemma between two competing yet metabolically expensive weapon systems. This study emphasizes the importance of cutting-edge methods to encompass the lineages of smaller venomous species that have yet to be characterized in detail, allowing us to understand the biology of their venom systems and to mine this prolific resource for translational research.

Transcriptomic and proteomic analyses of the venom and venom glands of Centruroides hirsutipalpus, a dangerous scorpion from Mexico

Centruroides hirsutipalpus (Scorpiones: Buthidae) is related to the "striped scorpion" group inhabiting the western Pacific region of Mexico. Human accidents caused by this species are medically important due to the great number of people stung and the severity of the resulting intoxication. This communication reports an extensive venom characterization using high-throughput proteomic and Illumina transcriptomic sequencing performed with RNA purified from its venom glands. 2,553,529 reads were assembled into 44,579 transcripts. From these transcripts, 23,880 were successfully annoted using Trinotate. Using specialized databases and by performing bioinformatic searches, it was possible to identify 147 putative venom protein transcripts. These include α- and β-type sodium channel toxins (NaScTx), potassium channel toxins (KScTx) (α-, β-, δ-, γ- and λ-types), enzymes (metalloproteases, hyaluronidases, phospholipases, serine proteases, and monooxygenases), protease inhibitors, host defense peptides (HDPs) such as defensins, non-disulfide bridge peptides (NDBPs), anionic peptides, superfamily CAP proteins, insulin growth factor-binding proteins (IGFBPs), orphan peptides, and other venom components (La1 peptides). De novo tandem mass spectrometric sequencing of digested venom identificatied 50 peptides. The venom of C. hirsutipalpus contains the highest reported number (77) of transcripts encoding NaScTxs, which are the components responsible for human fatalities.

Venom Peptides with Dual Modulatory Activity on the Voltage-Gated Sodium Channel NaV1.1 Provide Novel Leads for Development of Antiepileptic Drugs

Voltage-gated sodium (Na_V) channels play a fundamental role in normal neurological function, especially via the initiation and propagation of action potentials. The Na_V1.1 subtype is found in inhibitory interneurons of the brain and it is essential for maintaining a balance between excitation and inhibition in neuronal networks. Heterozygous loss-of-function mutations of SCN1A, the gene encoding Na_V1.1, underlie Dravet syndrome (DS), a severe pediatric epilepsy. We recently demonstrated that selective inhibition of Na_V1.1 inactivation prevents seizures and premature death in a mouse model of DS. Thus, selective modulators of Na_V1.1 might be useful therapeutics for treatment of DS as they target the underlying molecular deficit. Numerous scorpion-venom peptides have been shown to modulate the activity of Na_V channels, but little is known about their activity at Na_V1.1. Here we report the isolation, sequence, three-dimensional structure, recombinant production, and functional characterization of two peptidic modulators of Na_V1.1 from venom of the buthid scorpion Hottentotta jayakari. These peptides, Hj1a and Hj2a, are potent agonists of Na_V1.1 (EC₅₀ of 17 and 32 nM, respectively), and they present dual α/β activity by modifying both the activation and inactivation properties of the channel. NMR studies of rHj1a indicate that it adopts a cystine-stabilized αβ fold similar to known scorpion toxins. Although Hj1a and Hj2a have only limited selectivity for Na_V1.1, their unusual dual mode of action provides an alternative approach to the development of selective Na_V1.1 modulators for the treatment of DS.

Venom characterization of the bark scorpion Centruroides edwardsii (Gervais 1843): Composition, biochemical activities and in vivo toxicity for potential prey

In this study, we characterize the venom of Centruroides edwardsii, one of the most abundant scorpions in urban and rural areas of Costa Rica, in terms of its biochemical constituents and their biological activities. C. edwardsii venom is rich in peptides but also contains some higher molecular weight protein components. No phospholipase A₂, hemolytic or fibrinogenolytic activities were found, but the presence of proteolytic and hyaluronidase enzymes was evidenced by zymography. Venom proteomic analysis indicates the presence of a hyaluronidase, several cysteine-rich secretory proteins, metalloproteinases and a peptidylglycine α-hydroxylating monooxygenase like-enzyme. It also includes peptides similar to the K⁺-channel blocker margatoxin, a dominant toxin in the venom of the related scorpion C. margaritatus. MS and N-terminal sequencing analysis also reveals the presence of Na⁺-channel-modulating peptides with sequence similarity to orthologs present in other scorpion species of the genera Centruroides and Tityus. We purified the hyaluronidase (which co-eluted with an allergen 5-like CRiSP) and sequenced ~60% of this enzyme. We also sequenced some venom gland transcripts that include other cysteine-containing peptides and a Non-Disulfide Bridged Peptide (NDBP). Our in vivo experiments characterizing the effects on potential predators and prey show that C. edwardsii venom induces paralysis in several species of arthropods and geckos; crickets being the most sensitive and cockroaches and scorpions the most resistant organisms tested. Envenomation signs were also observed in mice, but no lethality was reached by intraperitoneal administration of this venom up to 120 μg/g body weight.

Transcriptome annotation and characterization of novel toxins in six scorpion species

Background

Venom has evolved in parallel in multiple animals for the purpose of self-defense, prey capture or both. These venoms typically consist of highly complex mixtures of toxins: diverse bioactive peptides and/or proteins each with a specific pharmacological activity. Because of their specificity, they can be used as experimental tools to study cell mechanisms and develop novel medicines and drugs. It is therefore potentially valuable to explore the venoms of various animals to characterize their toxins and identify novel toxin-families. This study focuses on the annotation and exploration of the transcriptomes of six scorpion species from three different families. The transcriptomes were annotated with a custom-built automated pipeline, primarily consisting of Basic Local Alignment Search Tool searches against UniProt databases and filter steps based on transcript coverage.

Results

We annotated the transcriptomes of four scorpions from the family Buthidae, one from Iuridae and one from Diplocentridae using our annotation pipeline. We found that the four buthid scorpions primarily produce disulfide-bridged ion-channel targeting toxins, while the non-buthid scorpions have a higher abundance of non-disulfide-bridged toxins. Furthermore, analysis of the “unidentified” transcripts resulted in the discovery of six novel putative toxin families containing a total of 37 novel putative toxins. Additionally, 33 novel toxins in existing toxin-families were found. Lastly, 19 novel putative secreted proteins without toxin-like disulfide bonds were found.

Conclusions

We were able to assign most transcripts to a toxin family and classify the venom composition for all six scorpions. In addition to advancing our fundamental knowledge of scorpion venomics, this study may serve as a starting point for future research by facilitating the identification of the venom composition of scorpions and identifying novel putative toxin families.

Electronic supplementary material

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

The Dual α-Amidation System in Scorpion Venom Glands

Many peptides in scorpion venoms are amidated at their C-termini. This post-translational modification is paramount for the correct biological function of ion channel toxins and antimicrobial peptides, among others. The discovery of canonical amidation sequences in transcriptome-derived scorpion proproteins suggests that a conserved enzymatic α-amidation system must be responsible for this modification of scorpion peptides. A transcriptomic approach was employed to identify sequences putatively encoding enzymes of the α-amidation pathway. A dual enzymatic α-amidation system was found, consisting of the membrane-anchored, bifunctional, peptidylglycine α-amidating monooxygenase (PAM) and its paralogs, soluble monofunctional peptidylglycine α-hydroxylating monooxygenase (PHMm) and peptidyl-α-hydroxyglycine α-amidating lyase (PALm). Independent genes encode these three enzymes. Amino acid residues responsible for ion coordination and enzymatic activity are conserved in these sequences, suggesting that the enzymes are functional. Potential endoproteolytic recognition sites for proprotein convertases in the PAM sequence indicate that PAM-derived soluble isoforms may also be expressed. Sequences potentially encoding proprotein convertases (PC1 and PC2), carboxypeptidase E (CPE), and other enzymes of the α-amidation pathway, were also found, confirming the presence of this pathway in scorpions.

Convergent recruitment of adamalysin-like metalloproteases in the venom of the red bark centipede (Scolopocryptops sexspinosus)

Many venom proteins have presumably been convergently recruited by taxa from diverse venomous lineages. These toxic proteins have characteristics that allow them to remain stable in solution and have a high propensity for toxic effects on prey and/or potential predators. Despite this well-established convergent toxin recruitment, some toxins seem to be lineage specific. To further investigate the toxic proteins found throughout venomous lineages, venom proteomics and venom-gland transcriptomics were performed on two individual red bark centipedes (Scolopocryptops sexspinosus). Combining the protein phenotype with the transcript genotype resulted in the first in-depth venom characterization of S. sexspinosus, including 72 venom components that were identified in both the transcriptome and proteome and 1468 nontoxin transcripts identified in the transcriptome. Ten different toxin families were represented in the venom and venom gland with the majority of the toxins belonging to metalloproteases, CAPS (cysteine-rich secretory protein, antigen 5, and pathogenesis-related 1 proteins), and β-pore-forming toxins. Nine of these toxin families shared a similar proteomic structure to venom proteins previously identified from other centipedes. However, the most highly expressed toxin family, the adamalysin-like metalloproteases, has until now only been observed in the venom of snakes. We confirmed adamalysin-like metalloprotease activity by means of in vivo functional assays. The recruitment of an adamalysin-like metalloprotease into centipede venom represents a striking case of convergent evolution.

Dissecting Toxicity: The Venom Gland Transcriptome and the Venom Proteome of the Highly Venomous Scorpion Centruroides limpidus (Karsch, 1879)

Venom glands and soluble venom from the Mexican scorpion Centruroides limpidus (Karsch, 1879) were used for transcriptomic and proteomic analyses, respectively. An RNA-seq was performed by high-throughput sequencing with the Illumina platform. Approximately 80 million reads were obtained and assembled into 198,662 putative transcripts, of which 11,058 were annotated by similarity to sequences from available databases. A total of 192 venom-related sequences were identified, including Na⁺ and K⁺ channel-acting toxins, enzymes, host defense peptides, and other venom components. The most diverse transcripts were those potentially coding for ion channel-acting toxins, mainly those active on Na⁺ channels (NaScTx). Sequences corresponding to β- scorpion toxins active of K⁺ channels (KScTx) and λ-KScTx are here reported for the first time for a scorpion of the genus Centruroides. Mass fingerprint corroborated that NaScTx are the most abundant components in this venom. Liquid chromatography coupled to mass spectometry (LC-MS/MS) allowed the identification of 46 peptides matching sequences encoded in the transcriptome, confirming their expression in the venom. This study corroborates that, in the venom of toxic buthid scorpions, the more abundant and diverse components are ion channel-acting toxins, mainly NaScTx, while they lack the HDP diversity previously demonstrated for the non-buthid scorpions. The highly abundant and diverse antareases explain the pancreatitis observed after envenomation by this species.

The Dual Prey-Inactivation Strategy of Spiders-In-Depth Venomic Analysis of Cupiennius salei

Most knowledge of spider venom concerns neurotoxins acting on ion channels, whereas proteins and their significance for the envenomation process are neglected. The here presented comprehensive analysis of the venom gland transcriptome and proteome of Cupiennius salei focusses on proteins and cysteine-containing peptides and offers new insight into the structure and function of spider venom, here described as the dual prey-inactivation strategy. After venom injection, many enzymes and proteins, dominated by α-amylase, angiotensin-converting enzyme, and cysteine-rich secretory proteins, interact with main metabolic pathways, leading to a major disturbance of the cellular homeostasis. Hyaluronidase and cytolytic peptides destroy tissue and membranes, thus supporting the spread of other venom compounds. We detected 81 transcripts of neurotoxins from 13 peptide families, whereof two families comprise 93.7% of all cysteine-containing peptides. This raises the question of the importance of the other low-expressed peptide families. The identification of a venom gland-specific defensin-like peptide and an aga-toxin-like peptide in the hemocytes offers an important clue on the recruitment and neofunctionalization of body proteins and peptides as the origin of toxins.

Biochemical characterization of the venom of Central American scorpion Didymocentrus krausi Francke, 1978 (Diplocentridae) and its toxic effects in vivo and in vitro

Venoms of medically important scorpions from Buthidae family have been intensively studied, in contrast to non-buthid venoms, for which knowledge is scarce. In this work, we characterized the venom of a Diplocentridae species, Didymocentrus krausi, a small fossorial scorpion that inhabits the Tropical Dry Forest of Central America. D. krausi venom soluble fraction contains proteases with enzymatic activity on gelatin and casein. Mass spectrometry and venomic analysis confirmed the presence of elastase-like, cathepsin-O-like proteases and a neprilysin-like metalloproteinase. We did not detect phospholipase A₂, C or D, nor hyaluronidase activity in the venom. By homology-based venom gland transcriptomic analysis, NDBPs, a β-KTx-like peptide, and other putative toxin transcripts were found, which, together with a p-benzoquinone compound present in the venom, could potentially explain its direct hemolytic and cytotoxic effects in several mammalian cell lines. Cytotoxicity of D. krausi venom was higher than the effect of venoms from two buthid scorpion species distributed in Costa Rica, Centruroides edwardsii and Tityus pachyurus. Even though D. krausi venom was not lethal to mice or crickets, when injected in mouse gastrocnemius muscle at high doses it induced pathological effects at 24 h, which include myonecrosis, weak hemorrhage, and inflammatory infiltration. We observed an apparent thrombotic effect in the skin blood vessels, but no in vitro fibrinogenolytic activity was detected. In crickets, D. krausi venom induced toxicity and paralysis in short periods of time.

Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (Crotalus cerastes) lineages reveal little differential expression despite individual variation

Changes in gene expression can rapidly influence adaptive traits in the early stages of lineage diversification. Venom is an adaptive trait comprised of numerous toxins used for prey capture and defense. Snake venoms can vary widely between conspecific populations, but the influence of lineage diversification on such compositional differences are unknown. To explore venom differentiation in the early stages of lineage diversification, we used RNA-seq and mass spectrometry to characterize Sidewinder Rattlesnake (Crotalus cerastes) venom. We generated the first venom-gland transcriptomes and complementary venom proteomes for eight individuals collected across the United States and tested for expression differences across life history traits and between subspecific, mitochondrial, and phylotranscriptomic hypotheses. Sidewinder venom was comprised primarily of hemorrhagic toxins, with few cases of differential expression attributable to life history or lineage hypotheses. However, phylotranscriptomic lineage comparisons more than doubled instances of significant expression differences compared to all other factors. Nevertheless, only 6.4% of toxins were differentially expressed overall, suggesting that shallow divergence has not led to major changes in Sidewinder venom composition. Our results demonstrate the need for consensus venom-gland transcriptomes based on multiple individuals and highlight the potential for discrepancies in differential expression between different phylogenetic hypotheses.

Transcriptomic and Proteomic Analyses Reveal the Diversity of Venom Components from the Vaejovid Scorpion Serradigitus gertschi

To understand the diversity of scorpion venom, RNA from venomous glands from a sawfinger scorpion, Serradigitus gertschi, of the family Vaejovidae, was extracted and used for transcriptomic analysis. A total of 84,835 transcripts were assembled after Illumina sequencing. From those, 119 transcripts were annotated and found to putatively code for peptides or proteins that share sequence similarities with the previously reported venom components of other species. In accordance with sequence similarity, the transcripts were classified as potentially coding for 37 ion channel toxins; 17 host defense peptides; 28 enzymes, including phospholipases, hyaluronidases, metalloproteases, and serine proteases; nine protease inhibitor-like peptides; 10 peptides of the cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 protein superfamily; seven La1-like peptides; and 11 sequences classified as "other venom components". A mass fingerprint performed by mass spectrometry identified 204 components with molecular masses varying from 444.26 Da to 12,432.80 Da, plus several higher molecular weight proteins whose precise masses were not determined. The LC-MS/MS analysis of a tryptic digestion of the soluble venom resulted in the de novo determination of 16,840 peptide sequences, 24 of which matched sequences predicted from the translated transcriptome. The database presented here increases our general knowledge of the biodiversity of venom components from neglected non-buthid scorpions.

Female-biased population divergence in the venom of the Hentz striped scorpion (Centruroides hentzi)

Sex-biased genes are expressed at higher levels in one sex and contribute to phenotypic differences between males and females, as well as overall phenotypic variation within and among populations. Venom has evolved primarily for predation and defense, making venom expression a highly variable phenotype as a result of local adaptation. Several scorpion species have shown both intraspecific and intersexual venom variation, and males have been observed using venom in courtship and mating, suggesting the existence of venom-specific, sex-biased genes that may contribute to population divergence. We used reversed-phase high-performance liquid chromatography (RP-HPLC), Agilent protein bioanalyzer chips, nano-liquid chromatography mass spectrometry (nLC/MS/MS), and median lethal dose (LD₅₀) assays in fruit flies (Drosophila melanogaster) and banded crickets (Gryllodes sigillatus) to investigate proteomic and functional venom variation within and among three Florida populations of the Hentz striped scorpion (Centruroides hentzi). We found significant venom variation among populations, with females, not males, being responsible for this divergence. We also found significant variation in venom expression within populations, with males contributing more to within population variation than females. Our results provide evidence that male and female scorpions experience different natural and sexual selective pressures that have led to the expression of sex-biased venom genes and that these genes may be consequential in population divergence.

Venom-gland transcriptomics and venom proteomics of the giant Florida blue centipede, Scolopendra viridis

The limited number of centipede venom characterizations have revealed a rich diversity of toxins, and recent work has suggested centipede toxins may be more rapidly diversifying than previously considered. Additionally, many identified challenges in venomics research, including assembly and annotation methods, toxin quantification, and the ability to provide biological or technical replicates, have yet to be addressed in centipede venom characterizations. We performed high-throughput, quantifiable transcriptomic and proteomic methods on two individual Scolopendra viridis centipedes from North Florida. We identified 39 toxins that were proteomically confirmed, and 481 nontoxins that were expressed in the venom gland of S. viridis. The most abundant toxins expressed in the venom of S. viridis belonged to calcium and potassium ion-channel toxins, venom allergens, metalloproteases, and β-pore forming toxins. We compared our results to the previously characterized S. viridis from Morelos, Mexico, and found only five proteomically confirmed toxins in common to both localities, suggesting either extreme toxin divergence within S. viridis, or that these populations may represent entirely different species. By using multiple assembly and annotation methods, we generated a comprehensive and quantitative reference transcriptome and proteome of a Scolopendromorpha centipede species, while overcoming some of the challenges present in venomics research.

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.

Phenotypic Variation in Mojave Rattlesnake (Crotalus scutulatus) Venom Is Driven by Four Toxin Families

Phenotypic diversity generated through altered gene expression is a primary mechanism facilitating evolutionary response in natural systems. By linking the phenotype to genotype through transcriptomics, it is possible to determine what changes are occurring at the molecular level. High phenotypic diversity has been documented in rattlesnake venom, which is under strong selection due to its role in prey acquisition and defense. Rattlesnake venom can be characterized by the presence (Type A) or absence (Type B) of a type of neurotoxic phospholipase A2 (PLA2), such as Mojave toxin, that increases venom toxicity. Mojave rattlesnakes (Crotalus scutulatus), represent this diversity as both venom types are found within this species and within a single panmictic population in the Sonoran Desert. We used comparative venom gland transcriptomics of nine specimens of C. scutulatus from this region to test whether expression differences explain diversity within and between venom types. Type A individuals expressed significantly fewer toxins than Type B individuals owing to the diversity of C-type lectins (CTLs) and snake venom metalloproteinases (SVMPs) found in Type B animals. As expected, both subunits of Mojave toxin were exclusively found in Type A individuals but we found high diversity in four additional PLA2s that was not associated with a venom type. Myotoxin a expression and toxin number variation was not associated with venom type, and myotoxin a had the highest range of expression of any toxin class. Our study represents the most comprehensive transcriptomic profile of the venom type dichotomy in rattlesnakes and C. scutulatus. Even intra-specifically, Mojave rattlesnakes showcase the diversity of snake venoms and illustrate that variation within venom types blurs the distinction of the venom dichotomy.