True Lies: Using Proteomics to Assess the Accuracy of Transcriptome-Based Venomics in Centipedes Uncovers False Positives and Reveals Startling Intraspecific Variation in Scolopendra Subspinipes

Enlightening the toxinological dark matter of spider venom enzymes

Spiders produce highly adapted venoms featuring a complex mixture of biomolecules used mainly for hunting and defense. The most prominent components are peptidic neurotoxins, a major focus of research and drug development, whereas venom enzymes have been largely neglected. Nevertheless, investigation of venom enzymes not only reveals insights into their biological functions, but also provides templates for future industrial applications. Here we compared spider venom enzymes validated at protein level contained in the VenomZone database and from all publicly available proteo-transcriptomic spider venom datasets. We assigned reported enzymes to cellular processes and known venom functions, including toxicity, prey pre-digestion, venom preservation, venom component activation, and spreading factors. Our study unveiled extensive discrepancy between public databases and publications with regard to enzyme coverage, which impedes the development of novel spider venom enzyme-based applications. Uncovering the previously unrecognized abundance and diversity of venom enzymes will open new avenues for spider venom biodiscovery.

Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

Transcriptome-based analyses reveal venom diversity in two araneomorph spiders, Psechrus triangulus and Hippasa lycosina

The spiders Psechrus triangulus and Hippasa lycosina are widely distributed in Yunnan Province, China, and are important natural enemies of agricultural pests, yet studies regarding the composition of their venom are lacking. In this study, cDNA libraries were constructed from venom gland tissue of P. triangulus and H. lycosina and used for transcriptomic analysis. From the analysis, 39 and 31 toxin-like sequences were predicted for P. triangulus and H. lycosina, respectively. The predicted neurotoxin sequences were categorized according to cysteine sequence motifs, and the predicted neurotoxin sequences of P. triangulus and H. lycosina could be classified into 9 and 6 toxin families, respectively. In addition, potential acetylcholinesterase, hyaluronidase, and astaxanthin-like metalloproteinases were identified through annotation. In summary, transcriptomic techniques were invaluable in mining the gene expression information from these two spider species to explore the toxin composition of their venom and determine how they differ. Studies of this type provide essential baseline data for studying the evolution and physiological activities of spider toxins and for the potential development of medicinal compounds.

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.

Re-evaluating and dating myriapod diversification with phylotranscriptomics under a regime of dense taxon sampling

Recent transcriptomic studies of myriapod phylogeny have been based on relatively small datasets with <40 myriapod terminals and variably supported or contradicted the traditional morphological groupings of Progoneata and Dignatha. Here we amassed a large dataset of 104 myriapod terminals, including multiple species for each of the four myriapod classes. Across the tree, most nodes are stable and well supported. Most analyses across a range of gene occupancy levels provide moderate to strong support for a deep split of Myriapoda into Symphyla + Pauropoda (=Edafopoda) and an uncontradicted grouping of Chilopoda + Diplopoda (=Pectinopoda nov.), as in other recent transcriptome-based analyses; no analysis recovers Progoneata or Dignatha as clades. As in all recent multi-locus and phylogenomic studies, chilopod interrelationships resolve with Craterostigmus excluded from Amalpighiata rather than uniting with other centipedes with maternal brood care in Phylactometria. Diplopod ordinal interrelationships are largely congruent with morphology-based classifications. Chilognathan clades that are not invariably advocated by morphologists include Glomerida + Glomeridesmida, such that the volvation-related characters of pill millipedes may be convergent, and Stemmiulida + Polydesmida more closely allied to Juliformia than to Callipodida + Chordeumatida. The latter relationship implies homoplasy in spinnerets and contradicts Nematophora. A time-tree with nodes calibrated by 25 myriapod and six outgroup fossil terminals recovers Cambrian-Ordovician divergences for the deepest splits in Myriapoda, Edafopoda and Pectinopoda, predating the terrestrial fossil record of myriapods as in other published chronograms, whereas age estimates within Chilopoda and Diplopoda overlap with or do not appreciably predate the calibration fossils. The grouping of Chilopoda and Diplopoda is recovered in all our analyses and is formalized as Pectinopoda nov., named for the shared presence of mandibular comb lamellae. New taxonomic proposals for Chilopoda based on uncontradicted clades are Tykhepoda nov. for the three blind families of Scolopendromorpha that share a "sieve-type" gizzard, and Taktikospina nov. for Scolopendromorpha to the exclusion of Mimopidae.

Function Prediction of Peptide Toxins with Sequence-Based Multi-Tasking PU Learning Method

Peptide toxins generally have extreme pharmacological activities and provide a rich source for the discovery of drug leads. However, determining the optimal activity of a new peptide can be a long and expensive process. In this study, peptide toxins were retrieved from Uniprot; three positive-unlabeled (PU) learning schemes, adaptive basis classifier, two-step method, and PU bagging were adopted to develop models for predicting the biological function of new peptide toxins. All three schemes were embedded with 14 machine learning classifiers. The prediction results of the adaptive base classifier and the two-step method were highly consistent. The models with top comprehensive performances were further optimized by feature selection and hyperparameter tuning, and the models were validated by making predictions for 61 three-finger toxins or the external HemoPI dataset. Biological functions that can be identified by these models include cardiotoxicity, vasoactivity, lipid binding, hemolysis, neurotoxicity, postsynaptic neurotoxicity, hypotension, and cytolysis, with relatively weak predictions for hemostasis and presynaptic neurotoxicity. These models are discovery-prediction tools for active peptide toxins and are expected to accelerate the development of peptide toxins as drugs.

Venom resistance mechanisms in centipede show tissue specificity

Venomous animals utilize venom glands to secrete and store powerful toxins for intraspecific and/or interspecific antagonistic interactions, implying that tissue-specific resistance is essential for venom glands to anatomically separate toxins from other tissues. Here, we show the mechanism of tissue-specific resistance in centipedes (Scolopendra subspinipes mutilans), where the splice variant of the receptor repels its own toxin. Unlike the well-known resistance mechanism by mutation in a given exon, we found that the KCNQ1 channel is highly expressed in the venom gland as a unique splice variant in which the pore domain and transmembrane domain six, partially encoded by exon 6 (rather than 7 as found in other tissues), contain eleven mutated residues. Such a splice variant is sufficient to gain resistance to SsTx (a lethal toxin for giant prey capture) in the venom gland due to a partially buried binding site. Therefore, the tissue-specific KCNQ1 modification confers resistance to the toxins, establishing a safe zone in the venom-storing/secreting environment.

Toxin diversity revealed by de novo transcriptome assembly for venom gland in two species of spiders (Trichonephila clavata and Sinopoda pengi)

During long-term predator-prey coevolution, spiders have generated a vast diversity of toxins. Trichonephila clavata is a web-spinning spider whose large, well-constructed webs and venomous arsenal facilitate prey capture. In contrast, Sinopoda pengi is an ambush predator with agile locomotion and strong chelicerae for hunting. In this study, transcriptomic analysis was performed to describe the predicted toxins of S. pengi and T. clavata. A total of 43 and 47 of these unigenes from S. pengi and T. clavata, respectively, were predicted to have toxin activity. Putative neurotoxins were classified to the family level according to cysteine arrangement; 4 and 6 toxin families were produced by S. pengi and T. clavata, respectively. In addition, potential metalloproteases, acetylcholinesterases, serine proteases, hyaluronidases and phospholipases were found by annotation in databases. In summary, molecular templates with potential application value for medical and biological fields were obtained by classifying and characterizing presumed venom components, which established a foundation for further study of venom.

Venomics of the Central European Myrmicine Ants Myrmica rubra and Myrmica ruginodis

Animal venoms are a rich source of novel biomolecules with potential applications in medicine and agriculture. Ants are one of the most species-rich lineages of venomous animals. However, only a fraction of their biodiversity has been studied so far. Here, we investigated the venom components of two myrmicine (subfamily Myrmicinae) ants: Myrmica rubra and Myrmica ruginodis. We applied a venomics workflow based on proteotranscriptomics and found that the venoms of both species are composed of several protein classes, including venom serine proteases, cysteine-rich secretory protein, antigen 5 and pathogenesis-related 1 (CAP) superfamily proteins, Kunitz-type serine protease inhibitors and venom acid phosphatases. Several of these protein classes are known venom allergens, and for the first time we detected phospholipase A1 in the venom of M. ruginodis. We also identified two novel epidermal growth factor (EGF) family toxins in the M. ruginodis venom proteome and an array of additional EGF-like toxins in the venom gland transcriptomes of both species. These are similar to known toxins from the related myrmicine ant, Manica rubida, and the myrmecine (subfamily Myrmeciinae) Australian red bulldog ant Myrmecia gullosa, and are possibly deployed as weapons in defensive scenarios or to subdue prey. Our work suggests that M.rubra and M. ruginodis venoms contain many enzymes and other high-molecular-weight proteins that cause cell damage. Nevertheless, the presence of EGF-like toxins suggests that myrmicine ants have also recruited smaller peptide components into their venom arsenal. Although little is known about the bioactivity and function of EGF-like toxins, their presence in myrmicine and myrmecine ants suggests they play a key role in the venom systems of the superfamily Formicoidea. Our work adds to the emerging picture of ant venoms as a source of novel bioactive molecules and highlights the need to incorporate such taxa in future venom bioprospecting programs.

Modern venomics-Current insights, novel methods, and future perspectives in biological and applied animal venom research

Venoms have evolved >100 times in all major animal groups, and their components, known as toxins, have been fine-tuned over millions of years into highly effective biochemical weapons. There are many outstanding questions on the evolution of toxin arsenals, such as how venom genes originate, how venom contributes to the fitness of venomous species, and which modifications at the genomic, transcriptomic, and protein level drive their evolution. These questions have received particularly little attention outside of snakes, cone snails, spiders, and scorpions. Venom compounds have further become a source of inspiration for translational research using their diverse bioactivities for various applications. We highlight here recent advances and new strategies in modern venomics and discuss how recent technological innovations and multi-omic methods dramatically improve research on venomous animals. The study of genomes and their modifications through CRISPR and knockdown technologies will increase our understanding of how toxins evolve and which functions they have in the different ontogenetic stages during the development of venomous animals. Mass spectrometry imaging combined with spatial transcriptomics, in situ hybridization techniques, and modern computer tomography gives us further insights into the spatial distribution of toxins in the venom system and the function of the venom apparatus. All these evolutionary and biological insights contribute to more efficiently identify venom compounds, which can then be synthesized or produced in adapted expression systems to test their bioactivity. Finally, we critically discuss recent agrochemical, pharmaceutical, therapeutic, and diagnostic (so-called translational) aspects of venoms from which humans benefit.

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.

Antimicrobial, Insecticidal and Cytotoxic Activity of Linear Venom Peptides from the Pseudoscorpion Chelifer cancroides

Linear cationic venom peptides are antimicrobial peptides (AMPs) that exert their effects by damaging cell membranes. These peptides can be highly specific, and for some, a significant therapeutic value was proposed, in particular for treatment of bacterial infections. A prolific source of novel AMPs are arthropod venoms, especially those of hitherto neglected groups such as pseudoscorpions. In this study, we describe for the first time pharmacological effects of AMPs discovered in pseudoscorpion venom. We examined the antimicrobial, cytotoxic, and insecticidal activity of full-length Checacin1, a major component of the Chelifer cancroides venom, and three truncated forms of this peptide. The antimicrobial tests revealed a potent inhibitory activity of Checacin1 against several bacteria and fungi, including methicillin resistant Staphylococcus aureus (MRSA) and even Gram-negative pathogens. All peptides reduced survival rates of aphids, with Checacin1 and the C-terminally truncated Checacin11-21 exhibiting effects comparable to Spinosad, a commercially used pesticide. Cytotoxic effects on mammalian cells were observed mainly for the full-length Checacin1. All tested peptides might be potential candidates for developing lead structures for aphid pest treatment. However, as these peptides were not yet tested on other insects, aphid specificity has not been proven. The N- and C-terminal fragments of Checacin1 are less potent against aphids but exhibit no cytotoxicity on mammalian cells at the tested concentration of 100 µM.

Recombinant Production and Characterization of a New Toxin from Cryptops iheringi Centipede Venom Revealed by Proteome and Transcriptome Analysis

Among the Chilopoda class of centipede, the Cryptops genus is one of the most associated with envenomation in humans in the metropolitan region of the state of São Paulo. To date, there is no study in the literature about the toxins present in its venom. Thus, in this work, a transcriptomic characterization of the Cryptops iheringi venom gland, as well as a proteomic analysis of its venom, were performed to obtain a toxin profile of this species. These methods indicated that 57.9% of the sequences showed to be putative toxins unknown in public databases; among them, we pointed out a novel putative toxin named Cryptoxin-1. The recombinant form of this new toxin was able to promote edema in mice footpads with massive neutrophils infiltration, linking this toxin to envenomation symptoms observed in accidents with humans. Our findings may elucidate the role of this toxin in the venom, as well as the possibility to explore other proteins found in this work.

A Pseudoscorpion's Promising Pinch: The venom of Chelifer cancroides contains a rich source of novel compounds

With pedipalps modified for venom injection, some pseudoscorpions possess a unique venom delivery system, which evolved independently from those of other arachnids like scorpions and spiders. Up to now, only a few studies have been focused on pseudoscorpion venom, which either identified a small fraction of venom compounds, or were based on solely transcriptomic approaches. Only one study addressed the bioactivity of pseudoscorpion venom. Here, we expand existing knowledge about pseudoscorpion venom by providing a comprehensive proteomic and transcriptomic analysis of the venom of Chelifer cancroides. We identified the first putative genuine toxins in the venom of C. cancroides and we showed that a large fraction of the venom comprises novel compounds. In addition, we tested the activity of the venom at specific ion channels for the first time. These tests demonstrate that the venom of C. cancroides causes inhibition of a voltage-gated insect potassium channel (Shaker IR) and modulates the inactivation process of voltage-gated sodium channels from Varroa destructor. For one of the smallest venomous animals ever studied, today's toolkits enabled a comprehensive venom analysis. This is demonstrated by allocating our identified venom compounds to more than half of the prominent ion signals in MALDI-TOF mass spectra of venom samples. The present study is a starting point for understanding the complex composition and activity of pseudoscorpion venom and provides a potential rich source of bioactive compounds useable for basic research and industrial application.

A Combined Transcriptomics and Proteomics Approach Reveals the Differences in the Predatory and Defensive Venoms of the Molluscivorous Cone Snail Cylinder ammiralis (Caenogastropoda: Conidae)

Venoms are complex mixtures of proteins that have evolved repeatedly in the animal kingdom. Cone snail venoms represent one of the best studied venom systems. In nature, this venom can be dynamically adjusted depending on its final purpose, whether to deter predators or hunt prey. Here, the transcriptome of the venom gland and the proteomes of the predation-evoked and defensive venoms of the molluscivorous cone snail Cylinder ammiralis were catalogued. A total of 242 venom-related transcripts were annotated. The conotoxin superfamilies presenting more different peptides were O1, O2, T, and M, which also showed high expression levels (except T). The three precursors of the J superfamily were also highly expressed. The predation-evoked and defensive venoms showed a markedly distinct profile. A total of 217 different peptides were identified, with half of them being unique to one venom. A total of 59 peptides ascribed to 23 different protein families were found to be exclusive to the predatory venom, including the cono-insulin, which was, for the first time, identified in an injected venom. A total of 43 peptides from 20 protein families were exclusive to the defensive venom. Finally, comparisons of the relative abundance (in terms of number of peptides) of the different conotoxin precursor superfamilies showed that most of them present similar abundance regardless of the diet.

Recruitment of toxin-like proteins with ancestral venom function supports endoparasitic lifestyles of Myxozoa

Cnidarians are the oldest lineage of venomous animals and use nematocysts to discharge toxins. Whether venom toxins have been recruited to support parasitic lifestyles in the Endocnidozoa (Myxozoa + Polypodium) is, however, unknown. To examine this issue we variously employed transcriptomic, proteomic, associated molecular phylogenies, and localisation studies on representative primitive and derived myxozoans (Malacosporea and Myxosporea, respectively), Polypodium hydriforme, and the free-living staurozoan Calvadosia cruxmelitensis. Our transcriptomics and proteomics analyses provide evidence for expression and translation of venom toxin homologs in myxozoans. Phylogenetic placement of Kunitz type serine protease inhibitors and phospholipase A2 enzymes reveals modification of toxins inherited from ancestral free-living cnidarian toxins, and that venom diversity is reduced in myxozoans concordant with their reduced genome sizes. Various phylogenetic analyses of the Kunitz-type toxin family in Endocnidozoa suggested lineage-specific gene duplications, which offers a possible mechanism for enhancing toxin diversification. Toxin localisation in the malacosporean Buddenbrockia plumatellae substantiates toxin translation and thus illustrates a repurposing of toxin function for endoparasite development and interactions with hosts, rather than for prey capture or defence. Whether myxozoan venom candidates are expressed in transmission stages (e.g. in nematocysts or secretory vesicles) requires further investigation.

ToxCodAn: a new toxin annotator and guide to venom gland transcriptomics

MOTIVATION

Next-generation sequencing has become exceedingly common and has transformed our ability to explore nonmodel systems. In particular, transcriptomics has facilitated the study of venom and evolution of toxins in venomous lineages; however, many challenges remain. Primarily, annotation of toxins in the transcriptome is a laborious and time-consuming task. Current annotation software often fails to predict the correct coding sequence and overestimates the number of toxins present in the transcriptome. Here, we present ToxCodAn, a python script designed to perform precise annotation of snake venom gland transcriptomes. We test ToxCodAn with a set of previously curated transcriptomes and compare the results to other annotators. In addition, we provide a guide for venom gland transcriptomics to facilitate future research and use Bothrops alternatus as a case study for ToxCodAn and our guide.

RESULTS

Our analysis reveals that ToxCodAn provides precise annotation of toxins present in the transcriptome of venom glands of snakes. Comparison with other annotators demonstrates that ToxCodAn has better performance with regard to run time ($>20x$ faster), coding sequence prediction ($>3x$ more accurate) and the number of toxins predicted (generating $>4x$ less false positives). In this sense, ToxCodAn is a valuable resource for toxin annotation. The ToxCodAn framework can be expanded in the future to work with other venomous lineages and detect novel toxins.

Phylogenetic analyses suggest centipede venom arsenals were repeatedly stocked by horizontal gene transfer

Venoms have evolved over a hundred times in animals. Venom toxins are thought to evolve mostly by recruitment of endogenous proteins with physiological functions. Here we report phylogenetic analyses of venom proteome-annotated venom gland transcriptome data, assisted by genomic analyses, to show that centipede venoms have recruited at least five gene families from bacterial and fungal donors, involving at least eight horizontal gene transfer events. These results establish centipedes as currently the only known animals with venoms used in predation and defence that contain multiple gene families derived from horizontal gene transfer. The results also provide the first evidence for the implication of horizontal gene transfer in the evolutionary origin of venom in an animal lineage. Three of the bacterial gene families encode virulence factors, suggesting that horizontal gene transfer can provide a fast track channel for the evolution of novelty by the exaptation of bacterial weapons into animal venoms.

Deadly Proteomes: A Practical Guide to Proteotranscriptomics of Animal Venoms

Animal venoms are renowned for their toxicity, biochemical complexity, and as a source of compounds with potential applications in medicine, agriculture, and industry. Polypeptides underlie much of the pharmacology of animal venoms, and elucidating these arsenals of polypeptide toxins-known as the venom proteome or venome-is an important step in venom research. Proteomics is used for the identification of venom toxins, determination of their primary structure including post-translational modifications, as well as investigations into the physiology underlying their production and delivery. Advances in proteomics and adjacent technologies has led to a recent upsurge in publications reporting venom proteomes. Improved mass spectrometers, better proteomic workflows, and the integration of next-generation sequencing of venom-gland transcriptomes and venomous animal genomes allow quicker and more accurate profiling of venom proteomes with greatly reduced starting material. Technologies such as imaging mass spectrometry are revealing additional insights into the mechanism, location, and kinetics of venom toxin production. However, these numerous new developments may be overwhelming for researchers designing venom proteome studies. Here, the field of venom proteomics is reviewed and some practical solutions for simplifying mass spectrometry workflows to study animal venoms are offered.

Proteo-Transcriptomic Analysis Identifies Potential Novel Toxins Secreted by the Predatory, Prey-Piercing Ribbon Worm Amphiporus lactifloreus

Nemerteans (ribbon worms) employ toxins to subdue their prey, but research thus far has focused on the small-molecule components of mucus secretions and few protein toxins have been characterized. We carried out a preliminary proteotranscriptomic analysis of putative toxins produced by the hoplonemertean Amphiporus lactifloreus (Hoplonemertea, Amphiporidae). No variants were found of known nemertean-specific toxin proteins (neurotoxins, cytotoxins, parbolysins or nemertides) but several toxin-like transcripts were discovered, expressed strongly in the proboscis, including putative metalloproteinases and sequences resembling sea anemone actitoxins, crown-of-thorn sea star plancitoxins, and multiple classes of inhibitor cystine knot/knottin family proteins. Some of these products were also directly identified in the mucus proteome, supporting their preliminary identification as secreted toxin components. Two new nemertean-typical toxin candidates could be described and were named U-nemertotoxin-1 and U-nemertotoxin-2. Our findings provide insight into the largely overlooked venom system of nemerteans and support a hypothesis in which the nemertean proboscis evolved in several steps from a flesh-melting organ in scavenging nemerteans to a flesh-melting and toxin-secreting venom apparatus in hunting hoplonemerteans.

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.

An Integrated Proteomic and Transcriptomic Analysis Reveals the Venom Complexity of the Bullet Ant Paraponera clavata

A critical hurdle in ant venom proteomic investigations is the lack of databases to comprehensively and specifically identify the sequence and function of venom proteins and peptides. To resolve this, we used venom gland transcriptomics to generate a sequence database that was used to assign the tandem mass spectrometry (MS) fragmentation spectra of venom peptides and proteins to specific transcripts. This was performed alongside a shotgun liquid chromatography-mass spectrometry (LC-MS/MS) analysis of the venom to confirm that these assigned transcripts were expressed as proteins. Through the combined transcriptomic and proteomic investigation of Paraponera clavata venom, we identified four times the number of proteins previously identified using 2D-PAGE alone. In addition to this, by mining the transcriptomic data, we identified several novel peptide sequences for future pharmacological investigations, some of which conform with inhibitor cysteine knot motifs. These types of peptides have the potential to be developed into pharmaceutical or bioinsecticide peptides.

Structural venomics reveals evolution of a complex venom by duplication and diversification of an ancient peptide-encoding gene

Spiders are one of the most successful venomous animals, with more than 48,000 described species. Most spider venoms are dominated by cysteine-rich peptides with a diverse range of pharmacological activities. Some spider venoms contain thousands of unique peptides, but little is known about the mechanisms used to generate such complex chemical arsenals. We used an integrated transcriptomic, proteomic, and structural biology approach to demonstrate that the lethal Australian funnel-web spider produces 33 superfamilies of venom peptides and proteins. Twenty-six of the 33 superfamilies are disulfide-rich peptides, and we show that 15 of these are knottins that contribute >90% of the venom proteome. NMR analyses revealed that most of these disulfide-rich peptides are structurally related and range in complexity from simple to highly elaborated knottin domains, as well as double-knot toxins, that likely evolved from a single ancestral toxin gene.

Solenodon genome reveals convergent evolution of venom in eulipotyphlan mammals

Venom systems are key adaptations that have evolved throughout the tree of life and typically facilitate predation or defense. Despite venoms being model systems for studying a variety of evolutionary and physiological processes, many taxonomic groups remain understudied, including venomous mammals. Within the order Eulipotyphla, multiple shrew species and solenodons have oral venom systems. Despite morphological variation of their delivery systems, it remains unclear whether venom represents the ancestral state in this group or is the result of multiple independent origins. We investigated the origin and evolution of venom in eulipotyphlans by characterizing the venom system of the endangered Hispaniolan solenodon (Solenodon paradoxus). We constructed a genome to underpin proteomic identifications of solenodon venom toxins, before undertaking evolutionary analyses of those constituents, and functional assessments of the secreted venom. Our findings show that solenodon venom consists of multiple paralogous kallikrein 1 (KLK1) serine proteases, which cause hypotensive effects in vivo, and seem likely to have evolved to facilitate vertebrate prey capture. Comparative analyses provide convincing evidence that the oral venom systems of solenodons and shrews have evolved convergently, with the 4 independent origins of venom in eulipotyphlans outnumbering all other venom origins in mammals. We find that KLK1s have been independently coopted into the venom of shrews and solenodons following their divergence during the late Cretaceous, suggesting that evolutionary constraints may be acting on these genes. Consequently, our findings represent a striking example of convergent molecular evolution and demonstrate that distinct structural backgrounds can yield equivalent functions.

Proteo-Transcriptomic Characterization of the Venom from the Endoparasitoid Wasp Pimpla turionellae with Aspects on Its Biology and Evolution

Within mega-diverse Hymenoptera, non-aculeate parasitic wasps represent 75% of all hymenopteran species. Their ovipositor dual-functionally injects venom and employs eggs into (endoparasitoids) or onto (ectoparasitoids) diverse host species. Few endoparasitoid wasps such as Pimpla turionellae paralyze the host and suppress its immune responses, such as encapsulation and melanization, to guarantee their offspring's survival. Here, the venom and its possible biology and function of P. turionellae are characterized in comparison to the few existing proteo-transcriptomic analyses on parasitoid wasp venoms. Multiple transcriptome assembly and custom-tailored search and annotation strategies were applied to identify parasitoid venom proteins. To avoid false-positive hits, only transcripts were finally discussed that survived strict filter settings, including the presence in the proteome and higher expression in the venom gland. P. turionella features a venom that is mostly composed of known, typical parasitoid enzymes, cysteine-rich peptides, and other proteins and peptides. Several venom proteins were identified and named, such as pimplin2, 3, and 4. However, the specification of many novel candidates remains difficult, and annotations ambiguous. Interestingly, we do not find pimplin, a paralytic factor in Pimpla hypochondriaca, but instead a new cysteine inhibitor knot (ICK) family (pimplin2), which is highly similar to known, neurotoxic asilid1 sequences from robber flies.

Parallel Evolution of Complex Centipede Venoms Revealed by Comparative Proteotranscriptomic Analyses

Centipedes are among the most ancient groups of venomous predatory arthropods. Extant species belong to five orders, but our understanding of the composition and evolution of centipede venoms is based almost exclusively on one order, Scolopendromorpha. To gain a broader and less biased understanding we performed a comparative proteotranscriptomic analysis of centipede venoms from all five orders, including the first venom profiles for the orders Lithobiomorpha, Craterostigmomorpha, and Geophilomorpha. Our results reveal an astonishing structural diversity of venom components, with 93 phylogenetically distinct protein and peptide families. Proteomically-annotated gene trees of these putative toxin families show that centipede venom composition is highly dynamic across macroevolutionary timescales, with numerous gene duplications as well as functional recruitments and losses of toxin gene families. Strikingly, not a single family is found in the venoms of representatives of all five orders, with 67 families being unique for single orders. Ancestral state reconstructions reveal that centipede venom originated as a simple cocktail comprising just four toxin families, with very little compositional evolution happening during the approximately 50 My before the living orders had diverged. Venom complexity then increased in parallel within the orders, with scolopendromorphs evolving particularly complex venoms. Our results show that even venoms composed of toxins evolving under the strong constraint of negative selection can have striking evolutionary plasticity on the compositional level. We show that the functional recruitments and losses of toxin families that shape centipede venom arsenals are not concentrated early in their evolutionary history, but happen frequently throughout.

The Diversity of Venom: The Importance of Behavior and Venom System Morphology in Understanding Its Ecology and Evolution

Venoms are one of the most convergent of animal traits known, and encompass a much greater taxonomic and functional diversity than is commonly appreciated. This knowledge gap limits the potential of venom as a model trait in evolutionary biology. Here, we summarize the taxonomic and functional diversity of animal venoms and relate this to what is known about venom system morphology, venom modulation, and venom pharmacology, with the aim of drawing attention to the importance of these largely neglected aspects of venom research. We find that animals have evolved venoms at least 101 independent times and that venoms play at least 11 distinct ecological roles in addition to predation, defense, and feeding. Comparisons of different venom systems suggest that morphology strongly influences how venoms achieve these functions, and hence is an important consideration for understanding the molecular evolution of venoms and their toxins. Our findings also highlight the need for more holistic studies of venom systems and the toxins they contain. Greater knowledge of behavior, morphology, and ecologically relevant toxin pharmacology will improve our understanding of the evolution of venoms and their toxins, and likely facilitate exploration of their potential as sources of molecular tools and therapeutic and agrochemical lead compounds.

Spider Venom: Components, Modes of Action, and Novel Strategies in Transcriptomic and Proteomic Analyses

This review gives an overview on the development of research on spider venoms with a focus on structure and function of venom components and techniques of analysis. Major venom component groups are small molecular mass compounds, antimicrobial (also called cytolytic, or cationic) peptides (only in some spider families), cysteine-rich (neurotoxic) peptides, and enzymes and proteins. Cysteine-rich peptides are reviewed with respect to various structural motifs, their targets (ion channels, membrane receptors), nomenclature, and molecular binding. We further describe the latest findings concerning the maturation of antimicrobial, and cysteine-rich peptides that are in most known cases expressed as propeptide-containing precursors. Today, venom research, increasingly employs transcriptomic and mass spectrometric techniques. Pros and cons of venom gland transcriptome analysis with Sanger, 454, and Illumina sequencing are discussed and an overview on so far published transcriptome studies is given. In this respect, we also discuss the only recently described cross contamination arising from multiplexing in Illumina sequencing and its possible impacts on venom studies. High throughput mass spectrometric analysis of venom proteomes (bottom-up, top-down) are reviewed.

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.

Venom collection and analysis in the pseudoscorpion Chelifer cancroides (Pseudoscorpiones: Cheliferidae)

Pseudoscorpions are very small arthropods with almost worldwide distribution. They possess a unique venom delivery system in the chelal hands of their pedipalps that has evolved independently from that of scorpions and spiders. Studies on the venom composition of pseudoscorpions are very rare. Recently, the potential venom composition of the pseudoscorpion Synsphyronus apimelus Harvey, 1987 (Pseudoscorpiones: Garypidae) has been studied by transcriptome analysis. However, a proteome analysis of venom to identify the genuine venom compounds of pseudoscorpions has not yet been performed. In our study, we have developed a non-invasive approach for extracting minute amounts of venom, which for the first time allowed collecting pure venom samples of pseudoscorpions with minimal contaminations and high reproducibility. These experiments first required a morphological investigation of the venom delivery system with a focus on the role of the lamina defensor in the release of venom. Likely, the venom delivery system of pseudoscorpions has a mechanism that prevents the release of venom if the prey is not successfully penetrated by a venom tooth. Electrical stimulation of a gland-containing chelal hand in combination with a mechanical stimulation of the lamina defensor at the base of the venom tooth resulted in an average of 5 nl of collected venom. The utility of the method was then validated by repeated venom extractions and subsequent analysis of the venom composition using MALDI-TOF mass fingerprinting. Subsequent proteomics analysis in combination with transcriptome analyses of chelal hand tissue has identified the first genuine venom compounds of pseudoscorpions with putative antimicrobial peptides. For our experiments, we used the house pseudoscorpion Chelifer cancroides (Linnaeus, 1758) (Pseudoscorpiones: Cheliferidae).

Can we resolve the taxonomic bias in spider venom research?

The rate of discovery of new spider species greatly exceeds the rate of spider venom characterisation, leading to an increasing number of species with unstudied venoms. However, recent advances in proteomics and genomics that enable the study of venoms from smaller species has expanded the accessible taxonomic range. Thus, although the number of unstudied spider venoms is likely to further increase, future research should focus on the characterisation of venoms and toxins from previously unstudied spider families.

•

The rate of discovery of new spider species by far exceeds the rate of spider venom characterisation.

•

Advances in proteomics and genomics techniques now allow the study of venoms from a much wider taxonomic range.

•

Closing the taxonomic gap at the family level should be the main focus for toxinologists over the next decade.

Shotgun Proteomics Analysis of Saliva and Salivary Gland Tissue from the Common Octopus Octopus vulgaris

The salivary apparatus of the common octopus ( Octopus vulgaris) has been the subject of biochemical study for over a century. A combination of bioassays, behavioral studies and molecular analysis on O. vulgaris and related species suggests that its proteome should contain a mixture of highly potent neurotoxins and degradative proteins. However, a lack of genomic and transcriptomic data has meant that the amino acid sequences of these proteins remain almost entirely unknown. To address this, we assembled the posterior salivary gland transcriptome of O. vulgaris and combined it with high resolution mass spectrometry data from the posterior and anterior salivary glands of two adults, the posterior salivary glands of six paralarvae and the saliva from a single adult. We identified a total of 2810 protein groups from across this range of salivary tissues and age classes, including 84 with homology to known venom protein families. Additionally, we found 21 short secreted cysteine rich protein groups of which 12 were specific to cephalopods. By combining protein expression data with phylogenetic analysis we demonstrate that serine proteases expanded dramatically within the cephalopod lineage and that cephalopod specific proteins are strongly associated with the salivary apparatus.

A comprehensive portrait of the venom of the giant red bull ant, Myrmecia gulosa, reveals a hyperdiverse hymenopteran toxin gene family

Ants (Hymenoptera: Formicidae) are diverse and ubiquitous, and their ability to sting is familiar to many of us. However, their venoms remain largely unstudied. We provide the first comprehensive characterization of a polypeptidic ant venom, that of the giant red bull ant, Myrmecia gulosa. We reveal a suite of novel peptides with a range of posttranslational modifications, including disulfide bond formation, dimerization, and glycosylation. One venom peptide has sequence features consistent with an epidermal growth factor fold, while the remaining peptides have features suggestive of a capacity to form amphipathic helices. We show that these peptides are derived from what appears to be a single, pharmacologically diverse, gene superfamily (aculeatoxins) that includes most venom peptides previously reported from the aculeate Hymenoptera. Two aculeatoxins purified from the venom were found to be capable of activating mammalian sensory neurons, consistent with the capacity to produce pain but via distinct mechanisms of action. Further investigation of the major venom peptide MIITX1-Mg1a revealed that it can also incapacitate arthropods, indicative of dual utility in both defense and predation. MIITX1-Mg1a accomplishes these functions by generating a leak in membrane ion conductance, which alters membrane potential and triggers neuronal depolarization. Our results provide the first insights into the evolution of the major toxin gene superfamily of the aculeate Hymenoptera and provide a new paradigm in the functional evolution of toxins from animal venoms.

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.

Studying Smaller and Neglected Organisms in Modern Evolutionary Venomics Implementing RNASeq (Transcriptomics)-A Critical Guide

Venoms are evolutionary key adaptations that species employ for defense, predation or competition. However, the processes and forces that drive the evolution of venoms and their toxin components remain in many aspects understudied. In particular, the venoms of many smaller, neglected (mostly invertebrate) organisms are not characterized in detail, especially with modern methods. For the majority of these taxa, even their biology is only vaguely known. Modern evolutionary venomics addresses the question of how venoms evolve by applying a plethora of -omics methods. These recently became so sensitive and enhanced that smaller, neglected organisms are now more easily accessible to comparatively study their venoms. More knowledge about these taxa is essential to better understand venom evolution in general. The methodological core pillars of integrative evolutionary venomics are genomics, transcriptomics and proteomics, which are complemented by functional morphology and the field of protein synthesis and activity tests. This manuscript focuses on transcriptomics (or RNASeq) as one toolbox to describe venom evolution in smaller, neglected taxa. It provides a hands-on guide that discusses a generalized RNASeq workflow, which can be adapted, accordingly, to respective projects. For neglected and small taxa, generalized recommendations are difficult to give and conclusions need to be made individually from case to case. In the context of evolutionary venomics, this overview highlights critical points, but also promises of RNASeq analyses. Methodologically, these concern the impact of read processing, possible improvements by perfoming multiple and merged assemblies, and adequate quantification of expressed transcripts. Readers are guided to reappraise their hypotheses on venom evolution in smaller organisms and how robustly these are testable with the current transcriptomics toolbox. The complementary approach that combines particular proteomics but also genomics with transcriptomics is discussed as well. As recently shown, comparative proteomics is, for example, most important in preventing false positive identifications of possible toxin transcripts. Finally, future directions in transcriptomics, such as applying 3rd generation sequencing strategies to overcome difficulties by short read assemblies, are briefly addressed.