The UniProtKB/Swiss-Prot Tox-Prot program: A central hub of integrated venom protein data

Peptidomic and transcriptomic profiling of four distinct spider venoms

Venom based research is exploited to find novel candidates for the development of innovative pharmacological tools, drug candidates and new ingredients for cosmetic and agrochemical industries. Moreover, venomics, as a well-established approach in systems biology, helps to elucidate the genetic mechanisms of the production of such a great molecular biodiversity. Today the advances made in the proteomics, transcriptomics and bioinformatics fields, favor venomics, allowing the in depth study of complex matrices and the elucidation even of minor compounds present in minute biological samples. The present study illustrates a rapid and efficient method developed for the elucidation of venom composition based on NextGen mRNA sequencing of venom glands and LC-MS/MS venom proteome profiling. The analysis of the comprehensive data obtained was focused on cysteine rich peptide toxins from four spider species originating from phylogenetically distant families for comparison purposes. The studied species were Heteropoda davidbowie (Sparassidae), Poecilotheria formosa (Theraphosidae), Viridasius fasciatus (Viridasiidae) and Latrodectus mactans (Theridiidae). This led to a high resolution profiling of 284 characterized cysteine rich peptides, 111 of which belong to the Inhibitor Cysteine Knot (ICK) structural motif. The analysis of H. davidbowie venom revealed a high richness in term of venom diversity: 95 peptide sequences were identified; out of these, 32 peptides presented the ICK structural motif and could be classified in six distinct families. The profiling of P. formosa venom highlighted the presence of 126 peptide sequences, with 52 ICK toxins belonging to three structural distinct families. V. fasciatus venom was shown to contain 49 peptide sequences, out of which 22 presented the ICK structural motif and were attributed to five families. The venom of L. mactans, until now studied for its large neurotoxins (Latrotoxins), revealed the presence of 14 cysteine rich peptides, out of which five were ICK toxins belonging to the CSTX superfamily. This in depth profiling of distinct ICK peptide families identified across the four spider species highlighted the high conservation of these neurotoxins among spider families.

Transcriptome Analysis to Understand the Toxicity of Latrodectus tredecimguttatus Eggs

Latrodectus tredecimguttatus is a kind of highly venomous black widow spider, with toxicity coming from not only venomous glands but also other parts of its body as well as newborn spiderlings and eggs. Up to date, although L. tredecimguttatus eggs have been demonstrated to be rich in proteinaceous toxins, there is no systematic investigation on such active components at transcriptome level. In this study, we performed a high-throughput transcriptome sequencing of L. tredecimguttatus eggs with Illumina sequencing technology. As a result, 53,284 protein-coding unigenes were identified, of which 14,185 unigenes produced significant hits in the available databases, including 280 unigenes encoding proteins or peptides homologous to known proteinaceous toxins. GO term and KEGG pathway enrichment analyses of the 280 unigenes showed that 375 GO terms and 18 KEGG pathways were significantly enriched. Functional analysis indicated that these unigene-coded toxins have the bioactivities to degrade tissue proteins, inhibit ion channels, block neuromuscular transmission, provoke anaphylaxis, induce apoptosis and hyperalgesia, etc. No known typical proteinaceous toxins in L. tredecimguttatus venomous glands, such as latrotoxins, were identified, suggesting that the eggs have a different toxicity mechanism from that of the venom. Our present transcriptome analysis not only helps to reveal the gene expression profile and toxicity mechanism of the L. tredecimguttatus eggs, but also provides references for the further related researches.

The antifungal plant defensin AtPDF2.3 from Arabidopsis thaliana blocks potassium channels

Scorpion toxins that block potassium channels and antimicrobial plant defensins share a common structural CSαβ-motif. These toxins contain a toxin signature (K-C₄-X-N) in their amino acid sequence, and based on in silico analysis of 18 plant defensin sequences, we noted the presence of a toxin signature (K-C₅-R-G) in the amino acid sequence of the Arabidopsis thaliana defensin AtPDF2.3. We found that recombinant (r)AtPDF2.3 blocks K_v1.2 and K_v1.6 potassium channels, akin to the interaction between scorpion toxins and potassium channels. Moreover, rAtPDF2.3[G36N], a variant with a KCXN toxin signature (K-C₅-R-N), is more potent in blocking K_v1.2 and K_v1.6 channels than rAtPDF2.3, whereas rAtPDF2.3[K33A], devoid of the toxin signature, is characterized by reduced K_v channel blocking activity. These findings highlight the importance of the KCXN scorpion toxin signature in the plant defensin sequence for blocking potassium channels. In addition, we found that rAtPDF2.3 inhibits the growth of Saccharomyces cerevisiae and that pathways regulating potassium transport and/or homeostasis confer tolerance of this yeast to rAtPDF2.3, indicating a role for potassium homeostasis in the fungal defence response towards rAtPDF2.3. Nevertheless, no differences in antifungal potency were observed between the rAtPDF2.3 variants, suggesting that antifungal activity and K_v channel inhibitory function are not linked.

Intragenome Diversity of Gene Families Encoding Toxin-like Proteins in Venomous Animals

The evolution of venoms is the story of how toxins arise and of the processes that generate and maintain their diversity. For animal venoms these processes include recruitment for expression in the venom gland, neofunctionalization, paralogous expansions, and functional divergence. The systematic study of these processes requires the reliable identification of the venom components involved in antagonistic interactions. High-throughput sequencing has the potential of uncovering the entire set of toxins in a given organism, yet the existence of non-venom toxin paralogs and the misleading effects of partial census of the molecular diversity of toxins make necessary to collect complementary evidence to distinguish true toxins from their non-venom paralogs. Here, we analyzed the whole genomes of two scorpions, one spider and one snake, aiming at the identification of the full repertoires of genes encoding toxin-like proteins. We classified the entire set of protein-coding genes into paralogous groups and monotypic genes, identified genes encoding toxin-like proteins based on known toxin families, and quantified their expression in both venom-glands and pooled tissues. Our results confirm that genes encoding toxin-like proteins are part of multigene families, and that these families arise by recruitment events from non-toxin genes followed by limited expansions of the toxin-like protein coding genes. We also show that failing to account for sequence similarity with non-toxin proteins has a considerable misleading effect that can be greatly reduced by comparative transcriptomics. Our study overall contributes to the understanding of the evolutionary dynamics of proteins involved in antagonistic interactions.

A combined proteomic and transcriptomic analysis of slime secreted by the southern bottletail squid, Sepiadarium austrinum (Cephalopoda)

Sepiadarium austrinum, the southern bottletail squid, is a small squid that inhabits soft sediments along Australia's south-east coast. When provoked, it rapidly secretes large volumes of slime, presumably as a form of chemical defense. We analyzed the proteomic composition of this slime using tandem mass spectrometry and transcriptomics and found that it was remarkably complex with 1735 identified protein groups (FDR:0.01). To investigate the chemical defense hypothesis we performed an Artemia toxicity assay and used sequence analysis to search for toxin-like molecules. Although the slime did not appear to be toxic to Artemia we found 13 proteins in slime with the hallmarks of toxins, namely cysteine richness, short length, a signal peptide and/or homology to known toxins. These included three short (80-130AA) cysteine rich secreted proteins with no homology to proteins on the NCBI or UniProt databases. Other protein families found included, CAP, phospholipase-B, ShKT-like peptides, peptidase S10, Kunitz BPTI and DNase II. Quantitative analysis using intensity based absolute quantification (iBAQ via MaxQuant) revealed 20 highly abundant proteins, accounting for 67% of iBAQ signal, and three of these were toxin-like. No mucin homologues were found suggesting that the structure of the slime gel may be formed by an unknown mechanism.

BIOLOGICAL SIGNIFICANCE

This study is the first known instance of a slime secretion from a cephalopod to be analyzed by proteomics methods and is the first investigation of a member of the family Sepiadariidae using proteomic methods. 1735 proteins were identified with 13 of these fitting criteria established for the identification of putative toxins. The slime is dominated by 20 highly abundant proteins with secreted, cysteine rich proteins. The study highlights the importance of 'omics approaches in understanding novel organisms.

Using a Novel Ontology to Inform the Discovery of Therapeutic Peptides from Animal Venoms

Venoms and venom-derived compounds constitute a rich and largely unexplored source of potentially therapeutic compounds. To facilitate biomedical research, it is necessary to design a robust informatics infrastructure that will allow semantic computation of venom concepts in a standardized, consistent manner. We have designed an ontology of venom-related concepts - named Venom Ontology - that reuses an existing public data source: UniProt's Tox-Prot database. In addition to describing the ontology and its construction, we have performed three separate case studies demonstrating its utility: (1) An exploration of venom peptide similarity networks within specific genera; (2) A broad overview of the distribution of available data among common taxonomic groups spanning the known tree of life; and (3) An analysis of the distribution of venom complexity across those same taxonomic groups. Venom Ontology is publicly available on BioPortal at http://bioportal.bioontology.org/ontologies/CU-VO.

Characterising the enzymatic profile of crude tentacle extracts from the South Atlantic jellyfish Olindias sambaquiensis (Cnidaria: Hydrozoa)

Jellyfish venoms are of medical and biotechnological importance, with toxins displaying antimicrobial, analgesic and anti-tumor activities. Although proteolytic enzymes have also been described, detailed characterisation of these proteins is scant in Olindias spp. High throughput mass spectrometry profiling of cnidarian venoms has become increasingly popular since the first description of the proteomic profile of putative toxins isolated from nematocysts of the hydrozoan jellyfish Olindias sambaquiensis describing the presence of orthologous enzymes as presented in venoms of advanced species as snakes. Rigorous bioinformatics analyses can aid functional annotation, but biochemical assays are prerequisite to unambiguously assign toxic function to a peptide or protein. Here we present results that experimentally confirm previously predicted proteomic analysis that crude venom extracts from tentacles of O. sambaquiensis are composed of polypeptides with metalloproteinase, serine proteinase and phospholipases A2 activities. Surprisingly, levels of serine proteinase and phospholipase A2 activities were comparable to those observed in venoms of Bothrops snakes which were used as positive controls in this study. Hence, these data offer new opportunities to explore serine proteinase and phospholipase A2 activities in the clinical sequelae following O. sambaquiensis envenomation, with future possible biopharmaceutical applications.

SpinachDB: A Well-Characterized Genomic Database for Gene Family Classification and SNP Information of Spinach

Spinach (Spinacia oleracea L.), which originated in central and western Asia, belongs to the family Amaranthaceae. Spinach is one of most important leafy vegetables with a high nutritional value as well as being a perfect research material for plant sex chromosome models. As the completion of genome assembly and gene prediction of spinach, we developed SpinachDB (http://222.73.98.124/spinachdb) to store, annotate, mine and analyze genomics and genetics datasets efficiently. In this study, all of 21702 spinach genes were annotated. A total of 15741 spinach genes were catalogued into 4351 families, including identification of a substantial number of transcription factors. To construct a high-density genetic map, a total of 131592 SSRs and 1125743 potential SNPs located in 548801 loci of spinach genome were identified in 11 cultivated and wild spinach cultivars. The expression profiles were also performed with RNA-seq data using the FPKM method, which could be used to compare the genes. Paralogs in spinach and the orthologous genes in Arabidopsis, grape, sugar beet and rice were identified for comparative genome analysis. Finally, the SpinachDB website contains seven main sections, including the homepage; the GBrowse map that integrates genome, genes, SSR and SNP marker information; the Blast alignment service; the gene family classification search tool; the orthologous and paralogous gene pairs search tool; and the download and useful contact information. SpinachDB will be continually expanded to include newly generated robust genomics and genetics data sets along with the associated data mining and analysis tools.

From Mollusks to Medicine: A Venomics Approach for the Discovery and Characterization of Therapeutics from Terebridae Peptide Toxins

Animal venoms comprise a diversity of peptide toxins that manipulate molecular targets such as ion channels and receptors, making venom peptides attractive candidates for the development of therapeutics to benefit human health. However, identifying bioactive venom peptides remains a significant challenge. In this review we describe our particular venomics strategy for the discovery, characterization, and optimization of Terebridae venom peptides, teretoxins. Our strategy reflects the scientific path from mollusks to medicine in an integrative sequential approach with the following steps: (1) delimitation of venomous Terebridae lineages through taxonomic and phylogenetic analyses; (2) identification and classification of putative teretoxins through omics methodologies, including genomics, transcriptomics, and proteomics; (3) chemical and recombinant synthesis of promising peptide toxins; (4) structural characterization through experimental and computational methods; (5) determination of teretoxin bioactivity and molecular function through biological assays and computational modeling; (6) optimization of peptide toxin affinity and selectivity to molecular target; and (7) development of strategies for effective delivery of venom peptide therapeutics. While our research focuses on terebrids, the venomics approach outlined here can be applied to the discovery and characterization of peptide toxins from any venomous taxa.

Prediction of Toxin Genes from Chinese Yellow Catfish Based on Transcriptomic and Proteomic Sequencing

Fish venom remains a virtually untapped resource. There are so few fish toxin sequences for reference, which increases the difficulty to study toxins from venomous fish and to develop efficient and fast methods to dig out toxin genes or proteins. Here, we utilized Chinese yellow catfish (Pelteobagrus fulvidraco) as our research object, since it is a representative species in Siluriformes with its venom glands embedded in the pectoral and dorsal fins. In this study, we set up an in-house toxin database and a novel toxin-discovering protocol to dig out precise toxin genes by combination of transcriptomic and proteomic sequencing. Finally, we obtained 15 putative toxin proteins distributed in five groups, namely Veficolin, Ink toxin, Adamalysin, Za2G and CRISP toxin. It seems that we have developed a novel bioinformatics method, through which we could identify toxin proteins with high confidence. Meanwhile, these toxins can also be useful for comparative studies in other fish and development of potential drugs.

Tentacle Transcriptome and Venom Proteome of the Pacific Sea Nettle, Chrysaora fuscescens (Cnidaria: Scyphozoa)

Jellyfish venoms are rich sources of toxins designed to capture prey or deter predators, but they can also elicit harmful effects in humans. In this study, an integrated transcriptomic and proteomic approach was used to identify putative toxins and their potential role in the venom of the scyphozoan jellyfish Chrysaora fuscescens. A de novo tentacle transcriptome, containing more than 23,000 contigs, was constructed and used in proteomic analysis of C. fuscescens venom to identify potential toxins. From a total of 163 proteins identified in the venom proteome, 27 were classified as putative toxins and grouped into six protein families: proteinases, venom allergens, C-type lectins, pore-forming toxins, glycoside hydrolases and enzyme inhibitors. Other putative toxins identified in the transcriptome, but not the proteome, included additional proteinases as well as lipases and deoxyribonucleases. Sequence analysis also revealed the presence of ShKT domains in two putative venom proteins from the proteome and an additional 15 from the transcriptome, suggesting potential ion channel blockade or modulatory activities. Comparison of these potential toxins to those from other cnidarians provided insight into their possible roles in C. fuscescens venom and an overview of the diversity of potential toxin families in cnidarian venoms.

Venom gland transcriptome analyses of two freshwater stingrays (Myliobatiformes: Potamotrygonidae) from Brazil

Stingrays commonly cause human envenoming related accidents in populations of the sea, near rivers and lakes. Transcriptomic profiles have been used to elucidate components of animal venom, since they are capable of providing molecular information on the biology of the animal and could have biomedical applications. In this study, we elucidated the transcriptomic profile of the venom glands from two different freshwater stingray species that are endemic to the Paraná-Paraguay basin in Brazil, Potamotrygon amandae and Potamotrygon falkneri. Using RNA-Seq, we identified species-specific transcripts and overlapping proteins in the venom gland of both species. Among the transcripts related with envenoming, high abundance of hyaluronidases was observed in both species. In addition, we built three-dimensional homology models based on several venom transcripts identified. Our study represents a significant improvement in the information about the venoms employed by these two species and their molecular characteristics. Moreover, the information generated by our group helps in a better understanding of the biology of freshwater cartilaginous fishes and offers clues for the development of clinical treatments for stingray envenoming in Brazil and around the world. Finally, our results might have biomedical implications in developing treatments for complex diseases.

Gene Expressing and sRNA Sequencing Show That Gene Differentiation Associates with a Yellow Acer palmatum Mutant Leaf in Different Light Conditions

Acer palmatum Thunb., like other maples, is a widely ornamental-use small woody tree for leaf shapes and colors. Interestingly, we found a yellow-leaves mutant “Jingling Huangfeng” turned to green when grown in shade or low-density light condition. In order to study the potential mechanism, we performed high-throughput sequencing and obtained 1,082 DEGs in leaves grown in different light conditions that result in A. palmatum significant morphological and physiological changes. A total of 989 DEGs were annotated and clustered, of which many DEGs were found associating with the photosynthesis activity and pigment synthesis. The expression of CHS and FDR gene was higher while the expression of FLS gene was lower in full-sunlight condition; this may cause more colorful substance like chalcone and anthocyanin that were produced in full-light condition, thus turning the foliage to yellow. Moreover, this is the first available miRNA collection which contains 67 miRNAs of A. palmatum, including 46 conserved miRNAs and 21 novel miRNAs. To get better understanding of which pathways these miRNAs involved, 102 Unigenes were found to be potential targets of them. These results will provide valuable genetic resources for further study on the molecular mechanisms of Acer palmatum leaf coloration.

Maize pan-transcriptome provides novel insights into genome complexity and quantitative trait variation

Gene expression variation largely contributes to phenotypic diversity and constructing pan-transcriptome is considered necessary for species with complex genomes. However, the regulation mechanisms and functional consequences of pan-transcriptome is unexplored systematically. By analyzing RNA-seq data from 368 maize diverse inbred lines, we identified almost one-third nuclear genes under expression presence and absence variation, which tend to play regulatory roles and are likely regulated by distant eQTLs. The ePAV was directly used as “genotype” to perform GWAS for 15 agronomic phenotypes and 526 metabolic traits to efficiently explore the associations between transcriptomic and phenomic variations. Through a modified assembly strategy, 2,355 high-confidence novel sequences with total 1.9 Mb lengths were found absent within reference genome. Ten randomly selected novel sequences were fully validated with genomic PCR, including another two NBS_LRR candidates potentially affect flavonoids and disease-resistance. A simulation analysis suggested that the pan-transcriptome of the maize whole kernel is approaching a maximum value of 63,000 genes, and through developing two test-cross populations and surveying several most important yield traits, the dispensable genes were shown to contribute to heterosis. Novel perspectives and resources to discover maize quantitative trait variations were provided to better understand the kernel regulation networks and to enhance maize breeding.

Advances in venomics

The term “venomics” was coined to describe the global study of venom and venom glands, targeting comprehensive characterization of the whole toxin profile of a venomous animal by means of proteomics, transcriptomics, genomics and bioinformatics studies.

Screening Method for the Discovery of Potential Bioactive Cysteine-Containing Peptides Using 3D Mass Mapping

Animal venoms and toxins are a valuable source of bioactive peptides with pharmacologic relevance as potential drug leads. A large subset of biologically active peptides discovered up till now contain disulfide bridges that enhance stability and activity. To discover new members of this class of peptides, we developed a workflow screening specifically for those peptides that contain inter- and intra-molecular disulfide bonds by means of three-dimensional (3D) mass mapping. Two intrinsic properties of the sulfur atom, (1) its relatively large negative mass defect, and (2) its isotopic composition, allow for differentiation between cysteine-containing peptides and peptides lacking sulfur. High sulfur content in a peptide decreases the normalized nominal mass defect (NMD) and increases the normalized isotopic shift (NIS). Hence in a 3D plot of mass, NIS, and NMD, peptides with sulfur appear in this plot with a distinct spatial localization compared with peptides that lack sulfur. In this study we investigated the skin secretion of two frog species; Odorrana schmackeri and Bombina variegata. Peptides from the crude skin secretions were separated by nanoflow LC, and of all eluting peptides high resolution zoom scans were acquired in order to accurately determine both monoisotopic mass and average mass. Both the NMD and the NIS were calculated from the experimental data using an in-house developed MATLAB script. Candidate peptides exhibiting a low NMD and high NIS values were selected for targeted de novo sequencing, and this resulted in the identification of several novel inter- and intra-molecular disulfide bond containing peptides. Graphical Abstract ᅟ.

Designer and natural peptide toxin blockers of the KcsA potassium channel identified by phage display

Peptide neurotoxins are powerful tools for research, diagnosis, and treatment of disease. Limiting broader use, most receptors lack an identified toxin that binds with high affinity and specificity. This paper describes isolation of toxins for one such orphan target, KcsA, a potassium channel that has been fundamental to delineating the structural basis for ion channel function. A phage-display strategy is presented whereby ∼1.5 million novel and natural peptides are fabricated on the scaffold present in ShK, a sea anemone type I (SAK1) toxin stabilized by three disulfide bonds. We describe two toxins selected by sorting on purified KcsA, one novel (Hui1, 34 residues) and one natural (HmK, 35 residues). Hui1 is potent, blocking single KcsA channels in planar lipid bilayers half-maximally (Ki) at 1 nM. Hui1 is also specific, inhibiting KcsA-Shaker channels in Xenopus oocytes with a Ki of 0.5 nM whereas Shaker, Kv1.2, and Kv1.3 channels are blocked over 200-fold less well. HmK is potent but promiscuous, blocking KcsA-Shaker, Shaker, Kv1.2, and Kv1.3 channels with Ki of 1-4 nM. As anticipated, one Hui1 blocks the KcsA pore and two conserved toxin residues, Lys21 and Tyr22, are essential for high-affinity binding. Unexpectedly, potassium ions traversing the channel from the inside confer voltage sensitivity to the Hui1 off-rate via Arg23, indicating that Lys21 is not in the pore. The 3D structure of Hui1 reveals a SAK1 fold, rationalizes KcsA inhibition, and validates the scaffold-based approach for isolation of high-affinity toxins for orphan receptors.

In silico hybridization enables transcriptomic illumination of the nature and evolution of Myxozoa

Background

The Myxozoa, a group of oligocellular, obligate endoparasites, has long been poorly understood in an evolutionary context. Recent genome-level sequencing techniques such as RNA-seq have generated large amounts of myxozoan sequence data, providing valuable insight into their evolutionary history. However, sequences from host tissue contamination are present in next-generation sequencing reactions of myxozoan tissue, and differentiating between the two has been inadequately addressed. In order to shed light on the genetic underpinnings of myxozoan biology, assembled contigs generated from these studies that derived from the myxozoan must be decoupled from transcripts derived from host tissue and other contamination. This study describes a pipeline for categorization of transcripts asmyxozoan based on similarity searching with known host and parasite sequences, explores the extent to which host contamination is present in previously existing myxozoan datasets, and implements this pipeline on a newly sequenced transcriptome of Myxobolus pendula, a parasite of the common creek chub gill arch.

Methods

The insilico hybridization pipeline uses iterative BLAST searching and database-driven e-value comparison to categorize transcripts as deriving from host, parasite, or other contamination. Functional genetic analysis of M. pendula was conducted using further BLAST searching, Hidden Markov Modeling, and sequence alignment and phylogenetic reconstruction.

Results

Three RNA libraries of encysted M. pendula plasmodia were sequenced and subjected to the method. Nearly half of the final set of contiguous assembly sequences (47.3 %) was identified as putative myxozoan transcripts. Putative contamination was also identified in at least 1/3^rd of previously published myxozoan transcripts. The set of M. pendula transcripts was mined for a range of biologically insightful genes, including taxonomically restricted nematocyst structural proteins and nematocyst proteins identified through mass tandem spectrometry of other cnidarians. Several novel findings emerged, including a fourth myxozoan minicollagen gene, putative myxozoan toxin proteins,and extracellular matrix glycoproteins.

Conclusions

This study serves as a model for the handling of next-generation myxozoan sequence. The need for careful categorization was demonstrated in both previous and new sets of myxozoan sequences. The final set of confidently assigned myxozoan transcripts can be mined for any biologically relevant gene or gene family without spurious misidentification of host contamination as a myxozoan homolog. As exemplified by M. pendula, the repertoire of myxozoan polar capsules may be more complex than previously thought, with an additional minicollagen homolog and putative expression of toxin proteins.

Electronic supplementary material

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

Gene duplications are extensive and contribute significantly to the toxic proteome of nematocysts isolated from Acropora digitifera (Cnidaria: Anthozoa: Scleractinia)

Background

Gene duplication followed by adaptive selection is a well-accepted process leading to toxin diversification in venoms. However, emergent genomic, transcriptomic and proteomic evidence now challenges this role to be at best equivocal to other processess . Cnidaria are arguably the most ancient phylum of the extant metazoa that are venomous and such provide a definitive ancestral anchor to examine the evolution of this trait.

Methods

Here we compare predicted toxins from the translated genome of the coral Acropora digitifera to putative toxins revealed by proteomic analysis of soluble proteins discharged from nematocysts, to determine the extent to which gene duplications contribute to venom innovation in this reef-building coral species. A new bioinformatics tool called HHCompare was developed to detect potential gene duplications in the genomic data, which is made freely available (https://github.com/rgacesa/HHCompare).

Results

A total of 55 potential toxin encoding genes could be predicted from the A. digitifera genome, of which 36 (65 %) had likely arisen by gene duplication as evinced using the HHCompare tool and verified using two standard phylogeny methods. Surprisingly, only 22 % (12/55) of the potential toxin repertoire could be detected following rigorous proteomic analysis, for which only half (6/12) of the toxin proteome could be accounted for as peptides encoded by the gene duplicates. Biological activities of these toxins are dominatedby putative phospholipases and toxic peptidases.

Conclusions

Gene expansions in A. digitifera venom are the most extensive yet described in any venomous animal, and gene duplication plays a significant role leading to toxin diversification in this coral species. Since such low numbers of toxins were detected in the proteome, it is unlikely that the venom is evolving rapidly by prey-driven positive natural selection. Rather we contend that the venom has a defensive role deterring predation or harm from interspecific competition and overgrowth by fouling organisms. Factors influencing translation of toxin encoding genes perhaps warrants more profound experimental consideration.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1976-4) contains supplementary material, which is available to authorized users.

Ecological venomics: How genomics, transcriptomics and proteomics can shed new light on the ecology and evolution of venom

Animal venom is a complex cocktail of bioactive chemicals that traditionally drew interest mostly from biochemists and pharmacologists. However, in recent years the evolutionary and ecological importance of venom is realized as this trait has direct and strong influence on interactions between species. Moreover, venom content can be modulated by environmental factors. Like many other fields of biology, venom research has been revolutionized in recent years by the introduction of systems biology approaches, i.e., genomics, transcriptomics and proteomics. The employment of these methods in venom research is known as 'venomics'. In this review we describe the history and recent advancements of venomics and discuss how they are employed in studying venom in general and in particular in the context of evolutionary ecology. We also discuss the pitfalls and challenges of venomics and what the future may hold for this emerging scientific field.

Whole-genome sequencing in outbreak analysis

In addition to the ever-present concern of medical professionals about epidemics of infectious diseases, the relative ease of access and low cost of obtaining, producing, and disseminating pathogenic organisms or biological toxins mean that bioterrorism activity should also be considered when facing a disease outbreak. Utilization of whole-genome sequencing (WGS) in outbreak analysis facilitates the rapid and accurate identification of virulence factors of the pathogen and can be used to identify the path of disease transmission within a population and provide information on the probable source. Molecular tools such as WGS are being refined and advanced at a rapid pace to provide robust and higher-resolution methods for identifying, comparing, and classifying pathogenic organisms. If these methods of pathogen characterization are properly applied, they will enable an improved public health response whether a disease outbreak was initiated by natural events or by accidental or deliberate human activity. The current application of next-generation sequencing (NGS) technology to microbial WGS and microbial forensics is reviewed.

Transcriptome and venom proteome of the box jellyfish Chironex fleckeri

BACKGROUND

The box jellyfish, Chironex fleckeri, is the largest and most dangerous cubozoan jellyfish to humans. It produces potent and rapid-acting venom and its sting causes severe localized and systemic effects that are potentially life-threatening. In this study, a combined transcriptomic and proteomic approach was used to identify C. fleckeri proteins that elicit toxic effects in envenoming.

RESULTS

More than 40,000,000 Illumina reads were used to de novo assemble ∼ 34,000 contiguous cDNA sequences and ∼ 20,000 proteins were predicted based on homology searches, protein motifs, gene ontology and biological pathway mapping. More than 170 potential toxin proteins were identified from the transcriptome on the basis of homology to known toxins in publicly available sequence databases. MS/MS analysis of C. fleckeri venom identified over 250 proteins, including a subset of the toxins predicted from analysis of the transcriptome. Potential toxins identified using MS/MS included metalloproteinases, an alpha-macroglobulin domain containing protein, two CRISP proteins and a turripeptide-like protease inhibitor. Nine novel examples of a taxonomically restricted family of potent cnidarian pore-forming toxins were also identified. Members of this toxin family are potently haemolytic and cause pain, inflammation, dermonecrosis, cardiovascular collapse and death in experimental animals, suggesting that these toxins are responsible for many of the symptoms of C. fleckeri envenomation.

CONCLUSIONS

This study provides the first overview of a box jellyfish transcriptome which, coupled with venom proteomics data, enhances our current understanding of box jellyfish venom composition and the molecular structure and function of cnidarian toxins. The generated data represent a useful resource to guide future comparative studies, novel protein/peptide discovery and the development of more effective treatments for jellyfish stings in humans. (Length: 300).

High-throughput production of two disulphide-bridge toxins

A quick and efficient production method compatible with high-throughput screening was developed using 36 toxins belonging to four different families of two disulphide-bridge toxins. Final toxins were characterized using HPLC co-elution, CD and pharmacological studies.

Chemoselective tarantula toxins report voltage activation of wild-type ion channels in live cells

Electrically excitable cells, such as neurons, exhibit tremendous diversity in their firing patterns, a consequence of the complex collection of ion channels present in any specific cell. Although numerous methods are capable of measuring cellular electrical signals, understanding which types of ion channels give rise to these signals remains a significant challenge. Here, we describe exogenous probes which use a novel mechanism to report activity of voltage-gated channels. We have synthesized chemoselective derivatives of the tarantula toxin guangxitoxin-1E (GxTX), an inhibitory cystine knot peptide that binds selectively to Kv2-type voltage gated potassium channels. We find that voltage activation of Kv2.1 channels triggers GxTX dissociation, and thus GxTX binding dynamically marks Kv2 activation. We identify GxTX residues that can be replaced by thiol- or alkyne-bearing amino acids, without disrupting toxin folding or activity, and chemoselectively ligate fluorophores or affinity probes to these sites. We find that GxTX-fluorophore conjugates colocalize with Kv2.1 clusters in live cells and are released from channels activated by voltage stimuli. Kv2.1 activation can be detected with concentrations of probe that have a trivial impact on cellular currents. Chemoselective GxTX mutants conjugated to dendrimeric beads likewise bind live cells expressing Kv2.1, and the beads are released by channel activation. These optical sensors of conformational change are prototype probes that can indicate when ion channels contribute to electrical signaling.

A Polychaete's powerful punch: venom gland transcriptomics of Glycera reveals a complex cocktail of toxin homologs

Glycerids are marine annelids commonly known as bloodworms. Bloodworms have an eversible proboscis adorned with jaws connected to venom glands. Bloodworms prey on invertebrates, and it is known that the venom glands produce compounds that can induce toxic effects in animals. Yet, none of these putative toxins has been characterized on a molecular basis. Here we present the transcriptomic profiles of the venom glands of three species of bloodworm, Glycera dibranchiata, Glycera fallax and Glycera tridactyla, as well as the body tissue of G. tridactyla. The venom glands express a complex mixture of transcripts coding for putative toxin precursors. These transcripts represent 20 known toxin classes that have been convergently recruited into animal venoms, as well as transcripts potentially coding for Glycera-specific toxins. The toxins represent five functional categories: Pore-forming and membrane-disrupting toxins, neurotoxins, protease inhibitors, other enzymes, and CAP domain toxins. Many of the transcripts coding for putative Glycera toxins belong to classes that have been widely recruited into venoms, but some are homologs of toxins previously only known from the venoms of scorpaeniform fish and monotremes (stonustoxin-like toxin), turrid gastropods (turripeptide-like peptides), and sea anemones (gigantoxin I-like neurotoxin). This complex mixture of toxin homologs suggests that bloodworms employ venom while predating on macroscopic prey, casting doubt on the previously widespread opinion that G. dibranchiata is a detritivore. Our results further show that researchers should be aware that different assembly methods, as well as different methods of homology prediction, can influence the transcriptomic profiling of venom glands.

Echidna venom gland transcriptome provides insights into the evolution of monotreme venom

Monotremes (echidna and platypus) are egg-laying mammals. One of their most unique characteristic is that males have venom/crural glands that are seasonally active. Male platypuses produce venom during the breeding season, delivered via spurs, to aid in competition against other males. Echidnas are not able to erect their spurs, but a milky secretion is produced by the gland during the breeding season. The function and molecular composition of echidna venom is as yet unknown. Hence, we compared the deeply sequenced transcriptome of an in-season echidna crural gland to that of a platypus and searched for putative venom genes to provide clues into the function of echidna venom and the evolutionary history of monotreme venom. We found that the echidna venom gland transcriptome was markedly different from the platypus with no correlation between the top 50 most highly expressed genes. Four peptides found in the venom of the platypus were detected in the echidna transcriptome. However, these genes were not highly expressed in echidna, suggesting that they are the remnants of the evolutionary history of the ancestral venom gland. Gene ontology terms associated with the top 100 most highly expressed genes in echidna, showed functional terms associated with steroidal and fatty acid production, suggesting that echidna “venom” may play a role in scent communication during the breeding season. The loss of the ability to erect the spur and other unknown evolutionary forces acting in the echidna lineage resulted in the gradual decay of venom components and the evolution of a new role for the crural gland.

Venom: the sharp end of pain therapeutics

Adequate pain control is still a significant challenge and largely unmet medical need in the 21st century. With many small molecules failing to reach required levels of potency and selectivity, drug discovery is once again turning to nature to replenish pain therapeutic pipelines. Venomous animals are frequently stereotyped as inflictors of pain and distress and have historically been vilified by mankind. Yet, ironically, the very venoms that cause pain when directly injected by the host animal may actually turn out to contain the next generation of analgesics when injected by the clinician. The last 12 months have seen dramatic discoveries of analgesic tools within venoms. Spiders, snakes and even centipedes are yielding peptides with immense therapeutic potential. Significant advances are also taking place in delivery methods that can improve bioavailability and pharmacokinetics of these exciting natural resources. Turning proteinaceous venom into pharmaceutical liquid gold is the goal of venomics and the focus of this article.

Atractaspis aterrima toxins: the first insight into the molecular evolution of venom in side-stabbers

Although snake venoms have been the subject of intense research, primarily because of their tremendous potential as a bioresource for design and development of therapeutic compounds, some specific groups of snakes, such as the genus Atractaspis, have been completely neglected. To date only limited number of toxins, such as sarafotoxins have been well characterized from this lineage. In order to investigate the molecular diversity of venom from Atractaspis aterrima—the slender burrowing asp, we utilized a high-throughput transcriptomic approach completed with an original bioinformatics analysis pipeline. Surprisingly, we found that Sarafotoxins do not constitute the major ingredient of the transcriptomic cocktail; rather a large number of previously well-characterized snake venom-components were identified. Notably, we recovered a large diversity of three-finger toxins (3FTxs), which were found to have evolved under the significant influence of positive selection. From the normalized and non-normalized transcriptome libraries, we were able to evaluate the relative abundance of the different toxin groups, uncover rare transcripts, and gain new insight into the transcriptomic machinery. In addition to previously characterized toxin families, we were able to detect numerous highly-transcribed compounds that possess all the key features of venom-components and may constitute new classes of toxins.

The first venomous crustacean revealed by transcriptomics and functional morphology: remipede venom glands express a unique toxin cocktail dominated by enzymes and a neurotoxin

Animal venoms have evolved many times. Venomous species are especially common in three of the four main groups of arthropods (Chelicerata, Myriapoda, and Hexapoda), which together represent tens of thousands of species of venomous spiders, scorpions, centipedes, and hymenopterans. Surprisingly, despite their great diversity of body plans, there is no unambiguous evidence that any crustacean is venomous. We provide the first conclusive evidence that the aquatic, blind, and cave-dwelling remipede crustaceans are venomous and that venoms evolved in all four major arthropod groups. We produced a three-dimensional reconstruction of the venom delivery apparatus of the remipede Speleonectes tulumensis, showing that remipedes can inject venom in a controlled manner. A transcriptomic profile of its venom glands shows that they express a unique cocktail of transcripts coding for known venom toxins, including a diversity of enzymes and a probable paralytic neurotoxin very similar to one described from spider venom. We screened a transcriptomic library obtained from whole animals and identified a nontoxin paralog of the remipede neurotoxin that is not expressed in the venom glands. This allowed us to reconstruct its probable evolutionary origin and underlines the importance of incorporating data derived from nonvenom gland tissue to elucidate the evolution of candidate venom proteins. This first glimpse into the venom of a crustacean and primitively aquatic arthropod reveals conspicuous differences from the venoms of other predatory arthropods such as centipedes, scorpions, and spiders and contributes valuable information for ultimately disentangling the many factors shaping the biology and evolution of venoms and venomous species.

Isolation of an orally active insecticidal toxin from the venom of an Australian tarantula

Many insect pests have developed resistance to existing chemical insecticides and consequently there is much interest in the development of new insecticidal compounds with novel modes of action. Although spiders have deployed insecticidal toxins in their venoms for over 250 million years, there is no evolutionary selection pressure on these toxins to possess oral activity since they are injected into prey and predators via a hypodermic needle-like fang. Thus, it has been assumed that spider-venom peptides are not orally active and are therefore unlikely to be useful insecticides. Contrary to this dogma, we show that it is possible to isolate spider-venom peptides with high levels of oral insecticidal activity by directly screening for per os toxicity. Using this approach, we isolated a 34-residue orally active insecticidal peptide (OAIP-1) from venom of the Australian tarantula Selenotypus plumipes. The oral LD50 for OAIP-1 in the agronomically important cotton bollworm Helicoverpa armigera was 104.2±0.6 pmol/g, which is the highest per os activity reported to date for an insecticidal venom peptide. OAIP-1 is equipotent with synthetic pyrethroids and it acts synergistically with neonicotinoid insecticides. The three-dimensional structure of OAIP-1 determined using NMR spectroscopy revealed that the three disulfide bonds form an inhibitor cystine knot motif; this structural motif provides the peptide with a high level of biological stability that probably contributes to its oral activity. OAIP-1 is likely to be synergized by the gut-lytic activity of the Bacillus thuringiensis Cry toxin (Bt) expressed in insect-resistant transgenic crops, and consequently it might be a good candidate for trait stacking with Bt.

KEGG orthology-based annotation of the predicted proteome of Acropora digitifera: ZoophyteBase - an open access and searchable database of a coral genome

BACKGROUND

Contemporary coral reef research has firmly established that a genomic approach is urgently needed to better understand the effects of anthropogenic environmental stress and global climate change on coral holobiont interactions. Here we present KEGG orthology-based annotation of the complete genome sequence of the scleractinian coral Acropora digitifera and provide the first comprehensive view of the genome of a reef-building coral by applying advanced bioinformatics.

DESCRIPTION

Sequences from the KEGG database of protein function were used to construct hidden Markov models. These models were used to search the predicted proteome of A. digitifera to establish complete genomic annotation. The annotated dataset is published in ZoophyteBase, an open access format with different options for searching the data. A particularly useful feature is the ability to use a Google-like search engine that links query words to protein attributes. We present features of the annotation that underpin the molecular structure of key processes of coral physiology that include (1) regulatory proteins of symbiosis, (2) planula and early developmental proteins, (3) neural messengers, receptors and sensory proteins, (4) calcification and Ca2+-signalling proteins, (5) plant-derived proteins, (6) proteins of nitrogen metabolism, (7) DNA repair proteins, (8) stress response proteins, (9) antioxidant and redox-protective proteins, (10) proteins of cellular apoptosis, (11) microbial symbioses and pathogenicity proteins, (12) proteins of viral pathogenicity, (13) toxins and venom, (14) proteins of the chemical defensome and (15) coral epigenetics.

CONCLUSIONS

We advocate that providing annotation in an open-access searchable database available to the public domain will give an unprecedented foundation to interrogate the fundamental molecular structure and interactions of coral symbiosis and allow critical questions to be addressed at the genomic level based on combined aspects of evolutionary, developmental, metabolic, and environmental perspectives.

Predicting function from sequence in a large multifunctional toxin family

Venoms contain active substances with highly specific physiological effects and are increasingly being used as sources of novel diagnostic, research and treatment tools for human disease. Experimental characterisation of individual toxin activities is a severe rate-limiting step in the discovery process, and in-silico tools which allow function to be predicted from sequence information are essential. Toxins are typically members of large multifunctional families of structurally similar proteins that can have different biological activities, and minor sequence divergence can have significant consequences. Thus, existing predictive tools tend to have low accuracy. We investigated a classification model based on physico-chemical attributes that can easily be calculated from amino-acid sequences, using over 250 (mostly novel) viperid phospholipase A₂ toxins. We also clustered proteins by sequence profiles, and carried out in-vitro tests for four major activities on a selection of isolated novel toxins, or crude venoms known to contain them. The majority of detected activities were consistent with predictions, in contrast to poor performance of a number of tested existing predictive methods. Our results provide a framework for comparison of active sites among different functional sub-groups of toxins that will allow a more targeted approach for identification of potential drug leads in the future.

Diversification of animal venom peptides-were jellyfish amongst the first combinatorial chemists?

An ocean of data: Huge numbers of protein toxins are found in animal venoms. This diversity is widely believed to have arisen by gene duplication events. However, recent data now challenges this tradition view. Here we highlight how jellyfish could hold the key to unravelling toxin diversification, with a view towards future combinatorial biosynthesis of toxin libraries.

Extensive inter- and intraspecific venom variation in closely related parasites targeting the same host: the case of Leptopilina parasitoids of Drosophila

The arms race between immune suppressive parasites that produce virulence factors and hosts that evolve resistance to these factors is suggested to be a key driver for the diversification of both partners. However, little is known regarding the diversity of virulence factors in closely related parasites or the mechanisms underlying the variation of virulence. One of the best-described model to address this issue is the interaction between Leptopilina parasitic wasps and their Drosophila hosts, in which variation of virulence is well documented. Thanks to a combined transcriptomic and proteomic approach, we have identified the main secreted proteins in the venom of Leptopilina heterotoma (Gotheron strain, 66 proteins) and of two well-characterized strains of Leptopilina boulardi, ISm and ISy (65 and 49 proteins, respectively). Results revealed significant quantitative differences in venom components between the L. boulardi strains, in agreement with their different virulence properties. Strikingly, the two related Leptopilina species did not share any abundant venom protein. The main identified proteins in L. boulardi were RhoGAPs and serpins while an aspartylglucosaminidase (AGA) was found abundant in L. heterotoma. The extensive quantitative variation observed between these species may be related with their use of different virulence strategies and/or to differences in their host range (specialist versus generalist). Altogether, our data suggests that parasitoid venom can quickly evolve, mainly through rapid changes in regulation of gene expression. It also evidences venom evolutionary processes largely described in other venomous animals i.e. the convergent recruitment of venom proteins between phylogenetically unrelated organisms, and the role of duplications in the emergence of multigenic families of virulence factors.

Proteomic characterisation of toxins isolated from nematocysts of the South Atlantic jellyfish Olindias sambaquiensis

Surprisingly little is known of the toxic arsenal of cnidarian nematocysts compared to other venomous animals. Here we investigate the toxins of nematocysts isolated from the jellyfish Olindias sambaquiensis. A total of 29 unique ms/ms events were annotated as potential toxins homologous to the toxic proteins from diverse animal phyla, including cone-snails, snakes, spiders, scorpions, wasp, bee, parasitic worm and other Cnidaria. Biological activities of these potential toxins include cytolysins, neurotoxins, phospholipases and toxic peptidases. The presence of several toxic enzymes is intriguing, such as sphingomyelin phosphodiesterase B (SMase B) that has only been described in certain spider venoms, and a prepro-haystatin P-IIId snake venom metalloproteinase (SVMP) that activates coagulation factor X, which is very rare even in snake venoms. Our annotation reveals sequence orthologs to many representatives of the most important superfamilies of peptide venoms suggesting that their origins in higher organisms arise from deep eumetazoan innovations. Accordingly, cnidarian venoms may possess unique biological properties that might generate new leads in the discovery of novel pharmacologically active drugs.

Mass landscapes of seven scorpion species: The first analyses of Australian species with 1,5-DAN matrix

Scorpion venoms have been studied for over fifty years; however, the majority of research has focussed primarily on medically important Buthidae species. Additionally, venoms of the estimated 200 species of scorpion native to Australia have received very little attention. The first venom mass profiles of six non-buthid and one buthid scorpion species are presented herein, four of which are endemic to Australia. While masses under 5 kDa dominated the venoms of all species, the buthid venom contained considerably more masses between 7 and 8 kDa than those of the non-buthids, corroborating the emergent trend that buthids are richer in long-chain neurotoxins than non-buthids. The Australian scorpion venom fractions were also analysed with the relatively new MALDI-ToF matrix 1,5-DAN. Over forty partial sequences were obtained, the majority of which are homologous to scorpion antimicrobials such as opistoporin and IsCT2. Overall, this study is the single most comprehensive mass spectrometric analysis of scorpion venom landscapes to date and provides an insight into untapped Australian species.