Transcriptome analysis of venom glands from a single fishing spider Dolomedes mizhoanus

Unveiling hidden toxin diversity: Discovery of novel venom components through manual curation of highly expressed sequences annotated as "no hits" in Phoneutria nigriventer spider venom gland transcriptome

Spider venoms have evolved over thousands of years, optimizing feeding and defense mechanisms. Venom components show pharmacological and biotechnological potential, rising interest in their study. However, the isolation of spider toxins for experimental evaluation poses significant challenges. To address this, transcriptomic analysis combined with computational tools has emerged as an appealing approach to characterizing spider venoms. However, many sequences remain unidentified after automatic annotation. In this study, we manually curated a subset of previously unannotated sequences from the Phoneutria nigriventer transcriptome and identified new putative venom components. Our manual analysis revealed 29 % of the analyzed sequences were potential venom components, 29 % hypothetical/uncharacterized proteins, and 17 % cellular function proteins. Only 25 % of the originally unannotated dataset remained without any identification. Most reclassified components were cysteine-rich peptides, including 23 novel putative toxins. We also found glycine-rich peptides (GRP), corroborating the previous description of GRPs in Phoneutria pertyi venom glands. Furthermore, to emphasize the recurrence of the lack of annotation in spider venom glands transcripts, we provide a survey of the percentage of unidentified sequences in several published spider venom transcriptomics studies. In conclusion, our study highlights the importance of manual curation in uncovering novel venom components and underscores the need for improved annotation strategies to fully exploit the medical and biotechnological potential of spider venoms.

Isolation of Peptide Inhibiting SGC-7901 Cell Proliferation from Aspongopus chinensis Dallas

Aspongopus chinensis Dallas is used as a traditional Chinese medicine as well as an edible insect. Although its anti-tumor effects have been observed, the anti-tumor active component(s) in the hemolymph of A. chinensis remains unknown. In this study, a combination usage of ultrafiltration, gel filtration chromatography, FPLC and RP-HPLC to separate and purify active peptides was performed based on the proliferation of the human gastric cancer SGC-7901 cell line treated with candidates. One peptide (MW = 2853.3 Da) was isolated from the hemolymph of A. chinensis. A total of 24 amino acid residues were continuously determined for the active peptide: N′-ECGYCAEKGIRCDDIHCCTGLKKK-C′. In conclusion, a peptide that can inhibit the proliferation of gastric cancer SGC-7901 cells in the hemolymph of A. chinensis was purified in this study, which is homologous to members of the spider toxin protein family. These results should facilitate further works for this peptide, such as the cloning of genes, expression in vitro by prokaryotic or eukaryotic systems, more specific tests of anti-tumor activity, and so on.

A Deep Learning Approach with Data Augmentation to Predict Novel Spider Neurotoxic Peptides

As major components of spider venoms, neurotoxic peptides exhibit structural diversity, target specificity, and have great pharmaceutical potential. Deep learning may be an alternative to the laborious and time-consuming methods for identifying these peptides. However, the major hurdle in developing a deep learning model is the limited data on neurotoxic peptides. Here, we present a peptide data augmentation method that improves the recognition of neurotoxic peptides via a convolutional neural network model. The neurotoxic peptides were augmented with the known neurotoxic peptides from UniProt database, and the models were trained using a training set with or without the generated sequences to verify the augmented data. The model trained with the augmented dataset outperformed the one with the unaugmented dataset, achieving accuracy of 0.9953, precision of 0.9922, recall of 0.9984, and F1 score of 0.9953 in simulation dataset. From the set of all RNA transcripts of Callobius koreanus spider, we discovered neurotoxic peptides via the model, resulting in 275 putative peptides of which 252 novel sequences and only 23 sequences showing homology with the known peptides by Basic Local Alignment Search Tool. Among these 275 peptides, four were selected and shown to have neuromodulatory effects on the human neuroblastoma cell line SH-SY5Y. The augmentation method presented here may be applied to the identification of other functional peptides from biological resources with insufficient data.

A Novel Insecticidal Spider Peptide that Affects the Mammalian Voltage-Gated Ion Channel hKv1.5

Spider venoms include various peptide toxins that modify the ion currents, mainly of excitable insect cells. Consequently, scientific research on spider venoms has revealed a broad range of peptide toxins with different pharmacological properties, even for mammal species. In this work, thirty animal venoms were screened against hKv1.5, a potential target for atrial fibrillation therapy. The whole venom of the spider Oculicosa supermirabilis, which is also insecticidal to house crickets, caused voltage-gated potassium ion channel modulation in hKv1.5. Therefore, a peptide from the spider O. supermirabilis venom, named Osu1, was identified through HPLC reverse-phase fractionation. Osu1 displayed similar biological properties as the whole venom; so, the primary sequence of Osu1 was elucidated by both of N-terminal degradation and endoproteolytic cleavage. Based on its primary structure, a gene that codifies for Osu1 was constructed de novo from protein to DNA by reverse translation. A recombinant Osu1 was expressed using a pQE30 vector inside the E. coli SHuffle expression system. recombinant Osu1 had voltage-gated potassium ion channel modulation of human hKv1.5, and it was also as insecticidal as the native toxin. Due to its novel primary structure, and hypothesized disulfide pairing motif, Osu1 may represent a new family of spider toxins.

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.

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.

An overview of Phoneutria nigriventer spider venom using combined transcriptomic and proteomic approaches

Phoneutria nigriventer is one of the largest existing true spiders and one of the few considered medically relevant. Its venom contains several neurotoxic peptides that act on different ion channels and chemical receptors of vertebrates and invertebrates. Some of these venom toxins have been shown as promising models for pharmaceutical or biotechnological use. However, the large diversity and the predominance of low molecular weight toxins in this venom have hampered the identification and deep investigation of the less abundant toxins and the proteins with high molecular weight. Here, we combined conventional and next-generation cDNA sequencing with Multidimensional Protein Identification Technology (MudPIT), to obtain an in-depth panorama of the composition of P. nigriventer spider venom. The results from these three approaches showed that cysteine-rich peptide toxins are the most abundant components in this venom and most of them contain the Inhibitor Cysteine Knot (ICK) structural motif. Ninety-eight sequences corresponding to cysteine-rich peptide toxins were identified by the three methodologies and many of them were considered as putative novel toxins, due to the low similarity to previously described toxins. Furthermore, using next-generation sequencing we identified families of several other classes of toxins, including CAPs (Cysteine Rich Secretory Protein-CRiSP, antigen 5 and Pathogenesis-Related 1-PR-1), serine proteinases, TCTPs (translationally controlled tumor proteins), proteinase inhibitors, metalloproteinases and hyaluronidases, which have been poorly described for this venom. This study provides an overview of the molecular diversity of P. nigriventer venom, revealing several novel components and providing a better basis to understand its toxicity and pharmacological activities.

De novo characterization of venom apparatus transcriptome of Pardosa pseudoannulata and analysis of its gene expression in response to Bt protein

Background

Pardosa pseudoannulata is a prevailing spider species, and has been regarded as an important bio-control agent of insect pests in farmland of China. However, the available genomic and transcriptomic databases of P. pseudoannulata and their venom are limited, which severely hampers functional genomic analysis of P. pseudoannulata. Recently high-throughput sequencing technology has been proved to be an efficient tool for profiling the transcriptome of relevant non-target organisms exposed to Bacillus thuringiensis (Bt) protein through food webs.

Results

In this study, the transcriptome of the venom apparatus was analyzed. A total of 113,358 non-redundant unigenes were yielded, among which 34,041 unigenes with complete or various length encoding regions were assigned biological function annotations and annotated with gene ontology and karyotic orthologous group terms. In addition, 3726 unigenes involved in response to stimulus and 720 unigenes associated with immune-response pathways were identified. Furthermore, we investigated transcriptomic changes in the venom apparatus using tag-based DGE technique. A total of 1724 differentially expressed genes (DEGs) were detected, while 75 and 372 DEGs were functionally annotated with KEGG pathways and GO terms, respectively. qPCR analyses were performed to verify the DEGs directly or indirectly related to immune and stress responses, including genes encoding heat shock protein, toll-like receptor, GST and NADH dehydrogenase.

Conclusion

This is the first study conducted to specifically investigate the venom apparatus of P. pseudoannulata in response to Bt protein exposure through tritrophic chain. A substantial fraction of transcript sequences was generated by high-throughput sequencing of the venom apparatus of P. pseudoannulata. Then a comparative transcriptome analysis showing a large number of candidate genes involved in immune response were identified by the tag-based DGE technology. This transcriptome dataset will provide a comprehensive sequence resource for furture molecular genetic research of the venom apparatus of P. pseudoannulata.

Electronic supplementary material

The online version of this article (10.1186/s12896-017-0392-z) contains supplementary material, which is available to authorized users.

Comparative transcriptomics of Entelegyne spiders (Araneae, Entelegynae), with emphasis on molecular evolution of orphan genes

Next-generation sequencing technology is rapidly transforming the landscape of evolutionary biology, and has become a cost-effective and efficient means of collecting exome information for non-model organisms. Due to their taxonomic diversity, production of interesting venom and silk proteins, and the relative scarcity of existing genomic resources, spiders in particular are excellent targets for next-generation sequencing (NGS) methods. In this study, the transcriptomes of six entelegyne spider species from three genera (Cicurina travisae, C. vibora, Habronattus signatus, H. ustulatus, Nesticus bishopi, and N. cooperi) were sequenced and de novo assembled. Each assembly was assessed for quality and completeness and functionally annotated using gene ontology information. Approximately 100 transcripts with evidence of homology to venom proteins were discovered. After identifying more than 3,000 putatively orthologous genes across all six taxa, we used comparative analyses to identify 24 instances of positively selected genes. In addition, between ~ 550 and 1,100 unique orphan genes were found in each genus. These unique, uncharacterized genes exhibited elevated rates of amino acid substitution, potentially consistent with lineage-specific adaptive evolution. The data generated for this study represent a valuable resource for future phylogenetic and molecular evolutionary research, and our results provide new insight into the forces driving genome evolution in taxa that span the root of entelegyne spider phylogeny.

Biotechnological Trends in Spider and Scorpion Antivenom Development

Spiders and scorpions are notorious for their fearful dispositions and their ability to inject venom into prey and predators, causing symptoms such as necrosis, paralysis, and excruciating pain. Information on venom composition and the toxins present in these species is growing due to an interest in using bioactive toxins from spiders and scorpions for drug discovery purposes and for solving crystal structures of membrane-embedded receptors. Additionally, the identification and isolation of a myriad of spider and scorpion toxins has allowed research within next generation antivenoms to progress at an increasingly faster pace. In this review, the current knowledge of spider and scorpion venoms is presented, followed by a discussion of all published biotechnological efforts within development of spider and scorpion antitoxins based on small molecules, antibodies and fragments thereof, and next generation immunization strategies. The increasing number of discovery and development efforts within this field may point towards an upcoming transition from serum-based antivenoms towards therapeutic solutions based on modern biotechnology.

Identification and characterization of toxins in the venom gland of the Chinese bird spider, Haplopelma hainanum, by transcriptomic analysis

Tarantula venoms provide a model system for studying toxin selectivity, structure-activity relationships and molecular evolution of peptide toxins. Previous studies have identified a large number of peptide toxins in the venom of the Chinese bird spider Haplopelma hainanum, generally regarded as a highly venomous spider. However, the lack of available RNA-seq transcriptomic and genomic data is an obstacle to understanding its venom at the molecular level. In this study, we investigated the venom gland transcriptome of H. hainanum by RNA-seq, in the absence of an available genomic sequence. We identified 201 potential toxins among 57 181 de novo assembled transcripts, including knottins, Kunitz-type toxins, enzymes and other proteins. We systematically identified most of the knottins and Kunitz-type toxins, some of which showed strongly biased expression in the venom gland, including members of the huwentoxin-1, huwentoxin-2 and magi-1 families. We also discovered several novel potential toxins. These data demonstrate the high molecular and structural diversity in the venom toxins of H. hainanum. This study offers a useful strategy for exploring the complex components of spider venoms.

ω-Tbo-IT1-New Inhibitor of Insect Calcium Channels Isolated from Spider Venom

Novel disulfide-containing polypeptide toxin was discovered in the venom of the Tibellus oblongus spider. We report on isolation, spatial structure determination and electrophysiological characterization of this 41-residue toxin, called ω-Tbo-IT1. It has an insect-toxic effect with LD50 19 μg/g in experiments on house fly Musca domestica larvae and with LD50 20 μg/g on juvenile Gromphadorhina portentosa cockroaches. Electrophysiological experiments revealed a reversible inhibition of evoked excitatory postsynaptic currents in blow fly Calliphora vicina neuromuscular junctions, while parameters of spontaneous ones were not affected. The inhibition was concentration dependent, with IC50 value 40 ± 10 nM and Hill coefficient 3.4 ± 0.3. The toxin did not affect frog neuromuscular junctions or glutamatergic and GABAergic transmission in rat brains. Ca(2+) currents in Calliphora vicina muscle were not inhibited, whereas in Periplaneta americana cockroach neurons at least one type of voltage gated Ca(2+) current was inhibited by ω-Tbo-IT1. Thus, the toxin apparently acts as an inhibitor of presynaptic insect Ca(2+) channels. Spatial structure analysis of the recombinant ω-Tbo-IT1 by NMR spectroscopy in aqueous solution revealed that the toxin comprises the conventional ICK fold containing an extended β-hairpin loop and short β-hairpin loop which are capable of making "scissors-like mutual motions".

A Comparative Analysis of the Venom Gland Transcriptomes of the Fishing Spiders Dolomedes mizhoanus and Dolomedes sulfurous

Dolomedes sulfurous and Dolomedes mizhoanus are predaceous arthropods catching and feeding on small fish. They live in the same area and have similar habits. Their venoms exhibit some similarities and differences in biochemical and electrophysiological properties. In the present work, we first performed a transcriptomic analysis by constructing the venom gland cDNA library of D. sulfurous and 127 novel putative toxin sequences were consequently identified, which were classified into eight families. This venom gland transcriptome was then compared with that of D. mizhoanus, which revealed that the putative toxins from both spider venoms might have originated from the same gene ancestors although novel toxins were evolved independently in the two spiders. The putative toxins from both spiders contain 6-12 cysteine residues forming seven cysteine patterns. As revealed by blast search, the two venoms are rich in neurotoxins targeting ion channels with pharmacological and therapeutic significance. This study provides insight into the venoms of two closely related species of spider, which will be of use for future investigations into the structure and function of their toxins.

Transcriptomics and venomics: implications for medicinal chemistry

Over the last three decades, transcriptomic studies of venom gland cells have continuously evolved, opening up new possibilities for exploring the molecular diversity of animal venoms, a prerequisite for the discovery of new drug candidates and molecular phylogenetics. The molecular complexity of animal venoms is much greater than initially thought. In this review, we describe the different technologies available for transcriptomic studies of venom, from the original individual cloning approaches to the more recent global Next Generation Sequencing strategies. Our understanding of animal venoms is evolving, with the discovery of complex and diverse bio-optimized cocktails of compounds, including mostly peptides and proteins, which are now beginning to be studied by academic and industrial researchers.

Analysis of transcriptomes of three orb-web spider species reveals gene profiles involved in silk and toxin

As an ancient arthropod with a history of 390 million years, spiders evolved numerous morphological forms resulting from adaptation to different environments. The venom and silk of spiders, which have promising commercial applications in agriculture, medicine and engineering fields, are of special interests to researchers. However, little is known about their genomic components, which hinders not only understanding spider biology but also utilizing their valuable genes. Here we report on deep sequenced and de novo assembled transcriptomes of three orb-web spider species, Gasteracantha arcuata, Nasoonaria sinensis and Gasteracantha hasselti which are distributed in tropical forests of south China. With Illumina paired-end RNA-seq technology, 54 871, 101 855 and 75 455 unigenes for the three spider species were obtained, respectively, among which 9 300, 10 001 and 10 494 unique genes are annotated, respectively. From these annotated unigenes, we comprehensively analyzed silk and toxin gene components and structures for the three spider species. Our study provides valuable transcriptome data for three spider species which previously lacked any genetic/genomic data. The results have laid the first fundamental genomic basis for exploiting gene resources from these spiders.

Peptidomics combined with cDNA library unravel the diversity of centipede venom

Centipedes are one of the oldest venomous arthropods using toxin as their weapon to capture prey. But little attention was focused on them and only few centipede toxins were demonstrated with activity on ion channels. Therefore, more deep works are needed to understand the diversity of centipede venom. In the present study, we use peptidomics combined with cDNA library to uncover the diversity of centipede Scolopendra subspinipes mutilans L. Koch. 192 peptides were identified by LC-MS/MS and 79 precursors were deduced by cDNA library. Surprisingly, the signal peptides of centipede toxins were more complicated than any other animal toxins and even exhibited large differences in homologues. Meanwhile, a large number of variants generated by alternative cleavage sites were detected by mass spectra. Odd number of cystein (3, 5, 7) found in the mature peptides were seldom seen in peptide toxins. In additional, two novel cysteine frameworks (C-C-C-CCC, C-C-C-C-CC-CC) were identified from 16 different cysteine frameworks from centipede peptides. Only 29 precursors have clear targets, while others may provide a potential diversity function for centipede. These findings highlight the extensive diversity of centipede toxins and provide powerful tools to understand the capture and defense weapon of centipede.

BIOLOGICAL SIGNIFICANCE

Peptide toxins from venomous animal have attracted increasing attentions due to their extraordinary chemical and pharmacological diversity. Centipedes are one of the most used Chinese traditional medicines, but little was known about the active components. The venom of Scolopendra subspinipes mutilans L. Koch is first deeply analyzed in this work and most of peptides were never discovered before. Interestingly, the number and arrangement of cysteine showed a larger different to known peptide toxins such spider or scorpion toxins. Moreover, only 29 peptides from this centipede venom were identified with known function. It suggested that our work not only important to understand the composition of centipede venom, but also provide many valuable peptides for potential biological functions.

Dragline silk: a fiber assembled with low-molecular-weight cysteine-rich proteins

Dragline silk has been proposed to contain two main protein constituents, MaSp1 and MaSp2. However, the mechanical properties of synthetic spider silks spun from recombinant MaSp1 and MaSp2 proteins have yet to approach natural fibers, implying the natural spinning dope is missing critical factors. Here we report the discovery of novel molecular constituents within the spinning dope that are extruded into dragline silk. Protein studies of the liquid spinning dope from the major ampullate gland, coupled with the analysis of dragline silk fibers using mass spectrometry, demonstrate the presence of a new family of low-molecular-weight cysteine-rich proteins (CRPs) that colocalize with the MA fibroins. Expression of the CRP family members is linked to dragline silk production, specifically MaSp1 and MaSp2 mRNA synthesis. Biochemical data support that CRP molecules are secreted into the spinning dope and assembled into macromolecular complexes via disulfide bond linkages. Sequence analysis supports that CRP molecules share similarities to members that belong to the cystine slipknot superfamily, suggesting that these factors may have evolved to increase fiber toughness by serving as molecular hubs that dissipate large amounts of energy under stress. Collectively, our findings provide molecular details about the components of dragline silk, providing new insight that will advance materials development of synthetic spider silk for industrial applications.

Comprehensive analysis of the venom gland transcriptome of the spider Dolomedes fimbriatus

A comprehensive transcriptome analysis of an expressed sequence tag (EST) database of the spider Dolomedes fimbriatus venom glands using single-residue distribution analysis (SRDA) identified 7,169 unique sequences. Mature chains of 163 different toxin-like polypeptides were predicted on the basis of well-established methodology. The number of protein precursors of these polypeptides was appreciably numerous than the number of mature polypeptides. A total of 451 different polypeptide precursors, translated from 795 unique nucleotide sequences, were deduced. A homology search divided the 163 mature polypeptide sequences into 16 superfamilies and 19 singletons. The number of mature toxins in a superfamily ranged from 2 to 49, whereas the diversity of the original nucleotide sequences was greater (2-261 variants). We observed a predominance of inhibitor cysteine knot toxin-like polypeptides among the cysteine-containing structures in the analyzed transcriptome bank. Uncommon spatial folds were also found.

Fish predation by semi-aquatic spiders: a global pattern

More than 80 incidences of fish predation by semi-aquatic spiders--observed at the fringes of shallow freshwater streams, rivers, lakes, ponds, swamps, and fens--are reviewed. We provide evidence that fish predation by semi-aquatic spiders is geographically widespread, occurring on all continents except Antarctica. Fish predation by spiders appears to be more common in warmer areas between 40° S and 40° N. The fish captured by spiders, usually ranging from 2-6 cm in length, are among the most common fish taxa occurring in their respective geographic area (e.g., mosquitofish [Gambusia spp.] in the southeastern USA, fish of the order Characiformes in the Neotropics, killifish [Aphyosemion spp.] in Central and West Africa, as well as Australian native fish of the genera Galaxias, Melanotaenia, and Pseudomugil). Naturally occurring fish predation has been witnessed in more than a dozen spider species from the superfamily Lycosoidea (families Pisauridae, Trechaleidae, and Lycosidae), in two species of the superfamily Ctenoidea (family Ctenidae), and in one species of the superfamily Corinnoidea (family Liocranidae). The majority of reports on fish predation by spiders referred to pisaurid spiders of the genera Dolomedes and Nilus (>75% of observed incidences). There is laboratory evidence that spiders from several more families (e.g., the water spider Argyroneta aquatica [Cybaeidae], the intertidal spider Desis marina [Desidae], and the 'swimming' huntsman spider Heteropoda natans [Sparassidae]) predate fish as well. Our finding of such a large diversity of spider families being engaged in fish predation is novel. Semi-aquatic spiders captured fish whose body length exceeded the spiders' body length (the captured fish being, on average, 2.2 times as long as the spiders). Evidence suggests that fish prey might be an occasional prey item of substantial nutritional importance.