- Research article
- Open Access
The venom gland transcriptome of the parasitoid wasp Nasonia vitripennis highlights the importance of novel genes in venom function
© The Author(s). 2016
- Received: 26 January 2016
- Accepted: 7 July 2016
- Published: 8 August 2016
Prior to egg laying the parasitoid wasp Nasonia vitripennis envenomates its pupal host with a complex mixture of venom peptides. This venom induces several dramatic changes in the host, including developmental arrest, immunosuppression, and altered metabolism. The diverse and potent bioactivity of N. vitripennis venom provides opportunities for the development of novel acting pharmaceuticals based on these molecules. However, currently very little is known about the specific functions of individual venom peptides or what mechanisms underlie the hosts response to envenomation. Many of the venom peptides also lack bioinformatically derived annotations because no homologs can be identified in the sequences databases. The RNA interference system of N. vitripennis provides a method for functional characterisation of venom protein encoding genes, however working with the current list of 79 candidates represents a daunting task. For this reason we were interested in determining the expression levels of venom encoding genes in the venom gland, as this information could be used to rank candidates for further study. To do this we carried out deep transcriptome sequencing of the venom gland and ovary tissue and used RNA-seq to rank the venom protein encoding genes by expression level. The generation of a specific venom gland transcriptome dataset also provides further opportunities to investigate novel features of this specialised organ.
RNA-seq revealed that the highest expressed venom encoding gene in the venom gland was ‘Venom protein Y’. The highest expressed annotated gene in this tissue was serine protease Nasvi2EG007167, which has previously been implicated in the apoptotic activity of N. vitripennis venom. As expected the RNA-seq confirmed that venom encoding genes are almost exclusively expressed in the venom gland relative to the neighbouring ovary tissue. Novel genes appear to perform key roles in N. vitripennis venom function, with over half of the 15 highest expressed venom encoding loci lacking bioinformatic annotations. The high throughput sequencing data also provided evidence for the existence of an additional 472 previously undescribed transcribed regions in the N. vitripennis genome. Finally, metatranscriptomic analysis of the venom gland transcriptome finds little evidence for the role of Wolbachia in the venom system.
The expression level information provided here for the N. vitripennis venom protein encoding genes represents a valuable dataset that can be used by the research community to rank candidates for further functional characterisation. These candidates represent bioactive peptides valuable in the development of new pharmaceuticals.
- Serine protease
- Drug development
- Venom gene
Nasonia vitripennis (Hymenoptera: Pteromalidae) (here after Nasonia) is a parasitoid wasp that utilises species of cyclorrhaphous diptera, such as house flies and blow flies, as hosts to support the growth and development of their offspring. Envenomation of the host increases the overall nutritional value of the food source available to the feeding wasp larvae. Some of these venom induced changes include developmental arrest, immune suppression and changes to the overall metabolite profile of the host [1–3]. The diverse and potent bioactivities of Nasonia venom peptides indicate that they may be useful in the development of novel pharmaceuticals [1, 4–7]. The venom peptides themselves are synthesised in the venom gland (VG) before being secreted into the venom reservoir ready for injection into the host via the ovipositor . A recent proteomic study on these venom reservoir extracts was able to identify 79 peptides [4, 9]. Although several of the peptides had similarity to sequences found in the venom of other animals, a majority had either not been associated with venom before or had no homology to sequences in the National Centre for Bioinformatic Information (NCBI) databases. Bioinformatic based annotation categories that could be applied to the venom peptides included immunity, metabolism, esterases, proteases, and recognition proteins . The most prevalent functional annotation was serine proteases, which is one of the largest gene families found in insect genomes [4, 10].
The response to Nasonia envenomation by the common laboratory host Sarcophaga bullata (Diptera: Sarcophagidae) has been extensively studied at both the molecular and physiological level. Experiments show that the respiration rate of S. bullata remains stable for up to 2 weeks following envenomation (although at lower levels than unstung hosts) indicating that the venom is not immediately lethal to this host . Other physiological changes to the host associated with envenomation include an increase in lipids and free amino acids, as well as the initiation of developmental arrest through an as yet an unknown mechanism. S. bullata immune processes are also affected, with the abundance and potency of circulating plasmatocytes declining dramatically (due to cell death) 60 min following envenomation by Nasonia . Importantly, Nasonia venom has been shown to be able to interfere with mammalian immunity pathways. For example, experiments using mouse cell lines demonstrated that Nasonia venom has both antimicrobial and anti-inflammatory properties . Reporter assays provided evidence that the observed immune responses were mediated through venom specific modifications to the steroid and Toll pathways . The activity of Nasonia venom in mammalian systems may indicate that venom peptides have evolved to target highly conserved core pathway components in order to limit the ability of the host to escape envenomation. In this situation evolution has already done much of the work in designing highly specialised peptides that target cellular pathways involved in diseases such as hypertension, diabetes, and cancer. Thus given the incredible species richness of the parasitoid wasps their venoms may represent a diverse untapped reservoir of molecules ready for incorporation into drug development pipelines.
In contrast to the extensive body of work exploring the host response to envenomation, relatively little is known about the individual venom peptides responsible for inducing these changes. Based on sequence homology a Nasonia venom metalloproteinase has been proposed to initiate developmental arrest by interfering with Notch signalling [1, 7, 13]. Supporting this proposal is the observation that several notch pathway genes are differentially expressed in the S. bullata hosts following envenomation. However, currently there is no direct evidence showing that this venom metalloprotease is specifically responsible for the observed expression changes in the Notch pathway . Several other venom peptides are thought to be involved in modifying host metabolic pathways based on homology to sequences from model systems. For example, a venom lipase could be responsible for the increase in phospholipid degradation observed in envenomated hosts [2, 14]. A venom trehalase has been proposed to convert the abundant host trehalose sugars into glucose . Additionally, several proteases found in Nasonia venom are homologous to peptides used by the ectoparasitoid Euplectrus separatae and the tick Haemaphysalis longicornis for blood meal digestion [15, 16].
The most direct functional characterisation of a venom encoding gene comes from a recent RNA interference (RNAi) study by Siebert et al. (2015). Knockdown of Nasonia venom calreticulin (Nasvi2EG037342) resulted in an increase in melanisation occurring at the oviposition wound on the host following envenomation . This observation suggested that venom calreticulin functions to interfere with normal wound healing in the host, perhaps by competing with the endogenous calreticulin in the host for binding of cofactors .
Several endoparasitoids from the Braconidae and Ichneumonidae wasp families inject viruses or virus like particles along with their venom in order to prevent the host immune system from attacking eggs deposited in the haemolymph. The virus like particles appear to provide protection to the developing parasitoid by induction apoptosis in host hemocytes [18–22]. As yet viral-like particles have not been implicated in Nasonia venom function and viral transcripts could not be identified in cDNA libraries synthesised from envenomated hosts . The intracellular symbiont Wolbachia is also present in high concentrations in the ovaries and host envenomation has been proposed as a possible mechanism for horizontal transfer of strains between parasitoid species .
Currently there is little direct experimental evidence supporting the proposed functions of the 79 Nasonia venom proteins and many of the bioinformatic annotations are based on homology to sequences from non-venomous model organisms. The RNAi knockdown system of Nasonia clearly provides an opportunity to begin to functionally characterise the venom protein encoding genes . However, given the number of genes involved, an important first step would be the identification of high value candidates. For this reason we carried out deep sequencing of the Nasonia VG and ovary in order to use the relative gene expression as a proxy for protein abundance in the venom reservoir. The RNA-seq information generated in this study also provides a valuable dataset for future studies into the reproductive system of this important developing model system.
RNA-seq analysis of the Nasonia VG and ovary transcriptomes
RNA-seq statistics of Nasonia vitripennis ovary and VG tissue
Multiple Alignments (%)b
No Feature (%)c
2.89E + 07
3.24E + 07
3.14E + 07
3.31E + 07
3.42E + 07
3.06E + 07
Overlap between the mapped reads and annotated features in the Nasonia official gene set version 2 (OGS2) is high, with an average of 81 % of mapped reads aligning to at least one gene body (Table 1). Of the 36,329 annotated genes in OGS2, 54 % (19,854) had detectable expression in either the VG or ovary. Based on normalised read counts, tissue specific expression (reads for a gene model detected in only one tissue) was observed for 3475 genes in the ovary and 1424 genes in the VG. Next we utilised the R-package DESeq to identify differentially expressed genes between the VG and ovary RNA-seq datasets . Significantly different expression was observed for 5794 (16 % of OGS2 loci) genes at an adjusted p-value threshold of 0.05 (Benjamini-Hochberg adjusted for multiple testing) (Figs. 1b and c; Additional file 1). Of the differentially expressed genes, 2524 were up and 3270 were down regulated in the VG relative to the ovary. Overall the fold changes of upregulated genes were significantly larger in the VG versus the ovary (Student’s t-test p-value < 0.001) (Fig. 1b).
The highly expressed venom protein encoding genes
Summary of VG and ovary RNA-seq data for the 79 Nasonia venom protein encoding genes
Ovary read countsb
VG read countsb
Log2 Fold Changec
Venom protein Y
Venom protein G
Venom protein Z
Venom protein Q
Venom protein X
Serine protease inhibitor 2
Venom protein L
Venom protein H
Venom protein P
Venom protein V
Venom protein I
Venom protein K
Venom protein M
Venom protein N
Serine protease homologue
Serine proteinase inhibitor
Venom protein D
Venom protein T
Chitin binding protein-like
Serine protease inhibitor 1
Venom protein J
Venom protein R
Venom protein F
Venom protein E
Venom protein S
γ-Glutamyl transpeptidase 2
Venom protein U
Venom protein O
C1q-like venom protein
1.00E + 00
Antigen 5-like protein
Venom protein W
1.00E + 00
Highest 15 expressed genes in the VG and ovary sorted in descending order by normalised read counts
Venom gland data
Ovary normalised countsb
VG normalised countsb
Log2 Fold Changec
Venom protein Y
Venom protein G
Venom protein Z
Venom protein Q
Venom protein X
Serine protease inhibitor 2
Venom protein L
Venom protein H
HEAT shock 70
CCHC-type zinc finger
ATP-dependent RNA helicase
40S ribosomal protein
Elongation factor 1-alpha
40S ribosomal protein
60S acidic ribosomal protein
60S ribosomal protein
60S acidic ribosomal protein
40S ribosomal protein
60S ribosomal protein
40S ribosomal protein
60S ribosomal protein
RNA-seq identifies 472 novel genes
Gene annotations of novel detected transcripts from Nasonia vitripennis VG and Ovary samples
Atp-dependent dna helicase q4
metal ion binding, helicase activity, heterocyclic compound binding, organic cyclic compound binding
Signal peptidase complex catalytic subunit sec11a
signal peptide processing, integral to membrane, proteolysis, serine-type peptidase activity
lyase activity, metabolic process
NADH dehydrogenase subunit 2
mitochondrial electron transport, NADH to ubiquinone, mitochondrial inner membrane, respiratory chain, respiratory chain, NADH dehydrogenase (ubiquinone) activity, integral to membrane
NADH dehydrogenase subunit partial
integral to membrane, mitochondrion, oxidation-reduction process, NADH dehydrogenase (ubiquinone) activity
Uncharacterized aarf domain-containing protein kinase 1
protein phosphorylation,oxidation-reduction process, protein kinase activity, ubiquinone biosynthetic process, ATP binding, flavin adenine dinucleotide binding oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen, NAD(P)H as one donor, and incorporation of one atom of oxygen
60s ribosomal protein
ribosome, structural constituent of ribosome, translation
28 s ribosomal protein mitochondrial
ribosome, structural constituent of ribosome, translation
signal peptidase complex subunit 3
peptidase activity, signal peptide processing, integral to membrane, signal peptidase complex
inosine-5 -monophosphate dehydrogenase 1b
oxidation-reduction process, adenyl nucleotide binding, metal ion binding, GMP biosynthetic process, cytoplasm, IMP dehydrogenase activity
Expression of venom protein encoding genes previously identified as high value candidates
Using our data we were then able to examine the expression pattern of genes previously identified as high value venom candidates based on bioinformatic annotations or results from studies in other venomous species [7, 9]. To differentiate the genes based broadly on their expression level we placed the 79 venom protein encoding genes into high, medium or low expression categories based on three equal divisions of normalised counts across our entire dataset (see methods). These categories were defined by normalised read count ranges of <181700, 181701-1218700,>1218701 for low, medium and high, respectively.
The highly expressed category was dominated by novel venom genes (Venom proteins G, Q, X, Y, and Z), with the single annotated loci being the trypsin-like serine protease Nasvi2EG007167 (Table 2). Serine proteases play broad functional roles across insect physiology, including apoptosis, immunity, and development [30, 31]. The medium expression category includes the gamma-glutamyltranspeptidase and endonuclease-like venom genes, which have been previously implicated in initiating venom induced apoptosis by interfering with normal metabolism of glutathione and by degrading nuclear DNA, respectively [7, 32, 33]. Similarly, the aminotransferase-like venom gene encodes an enzyme predicted to produce kynurenic acid from kynurenine and this may also be involved in apoptosis [34, 35]. Other medium expression category genes of note include IMPL-L2, enzymes involved in carbohydrate metabolism, and several serine proteases/proteinase inhibitor genes. Insulin binding protein (Nasvi2EG001168) is possibly involved in initiating the developmental arrest by blocking insulin signalling [7, 36].
The majority of venom protein encoding genes fit into the low expressed category based on our designations. Notable venom protein encoding genes in the low expression group include apyrase, beta-1-3-Glucan recognition protein, dipeptidyl peptidase, metalloprotease, trehalase, chitinase, and all five of the venom cystein-rich/Pacifastin protease inhibitors (Table 2). The first three of these aforementioned genes have been proposed to be involved in modifying the host immune response [7, 37, 38]. The metalloprotease has been suggested to be involved in venom induced developmental arrest by interfering with normal Notch signalling . The trehalase gene important for metabolising carbohydrates is also one of the few such genes expressed at higher levels in the ovary in our dataset. As noted by  the observation of proteins normally associated with breakdown of the chitin rich cuticle is somewhat puzzling given Nasonia larvae are able to mechanically access the host hemolymph using their mandibles [7, 39]. Antigen 5-like protein and the cystein-rich/pacifastin protease inhibitors have been identified in the venom of the Apis mellifera and Pimpla hypochondriaca [40, 41].
Metatranscriptomic analysis of Nasonia VG and ovary tissue
The expression levels of venom protein encoding genes obtained in this study provides additional information that can be used to select candidates for further functional analysis using either RNAi or recombinant techniques. The ovary RNA-seq also provides a proxy for the tissue specificity of the venom encoding gene expression that can be used to further refine candidate selection. It should be noted that high venom gene expression does not automatically imply that the encoded peptide is responsible for important envenomation phenotypes. Other biological features, such as translational regulators, post-translational modifications and mRNA/protein stability could also be critical to the function of gene products in the venom context. Indeed, it has been argued that as a result of extensive posttranslational modifications, final venom protein composition is largely independent of transcriptional control [45, 46]. The legitimacy of this view has been challenged by a more recent study showing strong correlations between the venom transcriptomes and proteomes of several snake species . Ultimately, additional high quality quantitative proteomic and transcriptomic data will be required to more fully appreciate the general importance of posttranslational mechanisms on venom composition. Finally, it is also worth considering that highly expressed venom genes may perform accessory functions, such as protecting the venom reservoir of the wasp from being attacked by its own venom, as has been observed in snakes .
The highest expressed annotated venom encoding gene was a trypsin-like serine protease. Unfortunately, serine proteases are a very common functional category in insect genomes making it difficult to propose a role for this highly expressed gene (especially in the venom context). Previous studies using Sf21 cells have shown that the apoptotic activity of Nasonia venom could be reduced by the addition of serine protease inhibitors . Interestingly, amongst the ten most highly expressed venom protein encoding genes with annotations, all three point to roles in apoptosis [7, 32, 33]. This may highlight the importance of apoptosis for the early envenomation effects, perhaps by assisting with the distribution of venom throughout the host or by suppressing the immune system by destroying circulating hemocytes . Nasonia venom has also been shown to induce apoptosis in S. bullata neural tissue and this has been suggested as a possible explanation for developmental arrest in envenomated this host .
Perhaps the most dramatic finding from this study was that over half of the top 15 expressed venom protein encoding genes have no functional annotation, including the highest expressed gene ‘Venom protein Y’. In retrospect maybe this should not be surprising given the complex evolutionary processes that drive specialisation of venom peptides . However, this observation does highlight the value of candidate selection methods (such as RNA-seq) that are not biased towards genes with existing annotation information. Even for those genes that are bioinformatically annotated, it should be recognised that most of the functional data comes from experiments in model systems that are non-venomous.
Another interesting finding from this study is that several venom peptides that have been previously implicated in venom function were shown to be only lowly expressed in our data [1, 7, 9]. Perhaps the best example is the metalloprotease (Nasvi2EG010516) predicted to be involved in initiating developmental arrest, was one of the lowest expressed of the venom protein encoding genes in our dataset [1, 52]. We also observed that all five cysteine rich protease inhibitors suggested to be involved in disrupting host immunity by inactivating the pro-phenol oxidase cascade were also expressed at low levels . The low expression of chitinase-type venom genes also suggests they may not be directly involved in venom function, especially given Nasonia larvae are able to mechanically disrupt the host integument during feeding.
Cufflink based transcript assembly of the VG and ovary RNA-seq data allowed us to identify 472 previously undescribed transcribed regions in the Nasonia genome. The low expression and lack of homology to other known genes might explain why they have not been previously described. Alternatively these observations may indicate they represent artefacts of the cufflinks methodology. Re-evaluation of the raw venom proteomic data may enable the identification of additional peptide fragments based on these new gene models. It is worth noting that the OGS2 models were developed based on RNA-seq data from the Wolbachia free Asymcx strain of Nasonia. As some of the new gene models identified here could be induced by Wolbachia infection they represent candidates for being involved in the cytoplasmic incompatibility phenotypes induced by this protobacteria.
Finally, the metatranscriptomic study revealed that viruses represented a large proportion of non-Nasonia reads in the VG and ovary transcripts. Particularly ‘Nasonia vitripennis virus’, which has a positive strand ssRNA genome and is part of the Picornavirales order and the Iflaviridae family but remains unassigned to lower taxonomic ranks . ‘Nasonia vitripennis virus’ has a ssRNA genome so reads may be derived from either messenger RNA or the ssRNA genome. The latter, opens up the possibility of performing a de novo assembly of the viral genome from our transcriptome data. Studies have shown that this virus causes no observed detrimental effect on Nasonia and it remains an open question as to whether these transcripts end up in the venom itself . As expected the ovary metatranscriptome was dominated by Wolbachia transcripts. In contrast, we identified relatively few Wolbachia sequences in the VG, suggesting this bacteria is unlikely to play an important role in venom function. However, it is important to note that the RNA-seq library preparation methods used in our study enriched for poly A+ reads and thus many bacterial sequences would not have been sequenced in this experiment. Follow up experiments using total RNA library preparation methods are required to conclusively rule out a possible role for Wolbachia in Nasonia venom function.
Animal venoms represent an important source of natural products with pharmaceutical value; as demonstrated by the development of important drugs, such as Exenatide from the Gila monster lizard, and Captopril isolated from the venom of the lancehead viper . More recently molecules isolated from reptile and arthropod venom have been shown to have therapeutic value in the treatment of cancer [54, 55]. With this in mind the venom of parasitoid hymenoptera holds particular promise, given the incredible species richness of this group of animals, a potentially vast and unique source of molecules awaits discovery. We also predict that venoms target highly conserved core cellular pathways, as a mechanism to limit the ability of the host to evolve detoxification strategies. Indeed, Nasonia venom has been shown to have bioactivity against mammalian immune pathways, further highlighting the pharmaceutical value of these molecules . The results reported here point to an important role for novel genes in venom function, and this information should assist with the selection of candidates in future functional studies.
In this study we used RNA-seq to generate a comprehensive dataset of expression information for VG and ovary in the developing model system Nasonia. Information on the expression level of the 79 venom protein encoding genes provides an unbiased approach for selecting RNAi candidates, especially as bioinformatic annotations are likely to be unreliable. As our knowledge of venom gene function increases our ability to understand the fundamental molecular mechanisms underlying envenomation will also be improved. This latter knowledge will be important in future efforts to use Nasonia venom in drug development pipelines.
Nasonia and tissue collection
The Nasonia (LABII strain) was a gift from Prof. Jack Werren, University of Rochester, USA. The Nasonia were hosted on Lucilia sericata pupae obtained from a commercial insectary (Biosupplies Ltd). Venom gland (VG) and ovary tissue from mated, host exposed, 1–3 day old females, was dissected under aseptic conditions and placed immediately into ice cold Trizol (Ambion). Each replicate contained VG or ovaries tissue pooled from between 40 and 50 individuals.
RNA-seq and functional annotation
RNA was extracted using the modified RNeasy (Qiagen) protocol. Briefly, tissue was initially homogenised in 1 ml of Trizol before the addition of 0.2 ml chloroform (Merck). The sample was centrifuged at 10,000 rpm for 10 min at room temperature and the supernatant was transferred to a Qiagen RNeasy column. The column was centrifuged at 10,000 rpm for 15 s before the bound RNA was washed twice with 0.5 ml Qiagen RPE buffer. The RNA was eluted in 30 ul of RNase free water as described in the RNeasy protocol. The quality of the purified RNA was verified using a Bioanalyzer (Agilent) with all samples having a RNA integrity score >7. Poly-T beads were used to enrich for mRNA and TruSeq (Illumina) stranded cDNA sequencing libraries were constructed by the Otago Genome Service (NZGL). The 100 bp paired end sequencing run was performed on the Illumina HiSeq2500 platform.
The raw RNA-seq data was processed with fastq-mcf to remove sequencing adapters and primers . The reads were quality trimmed to a Phred score of >20 using SolexaQA v3.1.3 . The paired-end reads were then mapped to the N. vitripennis genome  using Tophat2 (version 2.1) with the default settings, except for library type set as “fr-firststrand” based on the sequencing chemistry used during the library preparation . Read counts for gene feature in the N. vitripennis official gene set version 2 were generated using HTSeq-count (version 0.6.1p) with ‘union mode’ on exon features and incorporating strand information .
Differential expression between the ovary and VG read counts was determined using DESeq as described in the package vignette [25, 60]. A gene was considered differentially expressed if it had an adjusted p-value < 0.05. Principal component analysis was performed using R with the prcomp function. The heatmaps were generated using heatmap.2 (ggplots CRAN package) and clustered based on euclidean distance. A one-tailed enrichment analysis for enrichment of functional categories was carried out with Blast2GO using both the newly annotated venom and ovary expressed genes as the target and the entire gene set as the reference . Soley for the purposes of ranking genes based on expression level, we divided the DESeq normalised read counts by the transcript length (in kilobases) to compensate for the effects of gene length on counts. These length normalised counts are shown in Table 2 and Table 3 and also provided in Additional files 1, 2, 3 and 4. The expression level boundaries of high, medium, and lowly expressed venom genes were determined by binning all genes into windows of expression of 100. A bin was classified as low expressed if the value of the bin was multiplied by its number of containing genes and this value was less than a third of the sum of these values, high if equal or greater than two thirds, and medium elsewise.
The Cufflinks pipeline was used to identify putative alternative transcripts as well as unannotated genes based on comparisons to the current Nasonia gene models as described in the package documentation . Novel transcripts were annotated using Blast2GO v2.8 with the default parameters . Normalised read counts were generated for the new gene models by merging the cufflinks models with those of OGS version 2 and repeating the DESeq based analysis described above.
Reads that did not align with the Nasonia genome using Tophat2 were used in a Diamond BlastX search against the NCBI non-redundant (nr) database. Results presented do not include data from unmapped reads that were discarded by Tophat2 during the generation of non-redundant splice junctions. The high number of unmapped reads in the third VG sample necessitated that only a random sample of 10 % of the reads from this particular sample were used due to memory limitations. Megan5 was used for visualisation of the taxonomic distribution of reads in the data .
Complementary DNA was generated from 1 ug of RNA extracted from dissected ovary and VG tissue using the transcriptor first strand cDNA synthesis kit (Roche) as described by the manufacturer. The PCR was then performed using primers for the housekeeping gene RP49: 5’-CTTCCGCAAAGTCCTTGTTC-3’, 5’-TTTATTCATTCTCCTCAGAACG-3’. The primers specific for Nasvi2EG007167 were: 5’-TGGCTGTCATCAGATTGACG-3’, 5’-TATCCTGGAGCCAGTGTAG-3’. Reaction tubes were removed at 24 and 30 cycles, at a denaturing temperature of 94 °C for 30 s, annealing temperature of 55 °C for 30 s, and extension temperature of 72 °C for 45 s with one unit of taq polymerase (Roche).
We would like to thank Prof Jack Werren (University of Rochester) and the Werren lab for the gift of the Nasonia LABII strain. We also wish to thank Drs Patrick Biggs and Matthew Campbell for discussions and comments that improved the draft manuscript. We are grateful for the valuable feedback received from two anonymous reviewers. We also wish to thank Dr Helen Fitzsimons, Dr Natisha Magan, and James Connell for vital assistance with the import of Nasonia into New Zealand.
This work was supported by Massey Honours Scholarship to AS, Massey University start-up and MURF grants (RM17937) to DW.
Availability of data and materials
The raw RNA-seq reads are available at the SRA under accession: SRP067692. The RNA-seq data supporting the results of this article is available in the GEO repository under accession number GSE76257.
AS performed the bioinformatics and the molecular biology and helped in the preparation of this manuscript. DW conceived of the study, assisted with the bioinformatic analysis, and wrote the manuscript. Both authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This invertebrates research was exempt from animal ethics requirements by section 2 of the New Zealand Animal Welfare Act 1999.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Martinson EO, Wheeler D, Wright J, Mrinalini, Siebert AL, Werren JH. Nasonia vitripennis venom causes targeted gene expression changes in its fly host. Mol Ecol. 2014;23:5918–30.View ArticlePubMedGoogle Scholar
- Mrinalini, Siebert AL, Wright J, Martinson E, Wheeler D, Werren JH. Parasitoid venom induces metabolic cascades in fly hosts. Metabolomics. 2014;11:350–66.View ArticleGoogle Scholar
- Rivers DB, Denlinger DL. Venom-induced alterations in Fly lipid metabolism and its impact on larval development of the ectoparasitoid nasonia vitripennis (Walker) (Hymenoptera: Pteromalidae). J Invertebr Pathol. 1995;66:104–10.View ArticleGoogle Scholar
- Werren JH, Richards S, Desjardins CA, Niehuis O, Gadau J, Colbourne JK, et al. Functional and evolutionary insights from the genomes of three parasitoid Nasonia species. Science. 2010;327:343–8.View ArticlePubMedGoogle Scholar
- Aisha LS, Wheeler D, Werren JH. Combined parasitoid wasp RNA interference knockdown / fly host RNA sequencing to test and generate hypotheses about the function of venom calreticulin. Toxicon. 2015;100:304–16.Google Scholar
- Danneels EL, Gerlo S, Heyninck K, Van Craenenbroeck K, De Bosscher K, Haegeman G, et al. How the venom from the ectoparasitoid Wasp nasonia vitripennis exhibits anti-inflammatory properties on mammalian cell lines. PLoS One. 2014;9:e96825.View ArticlePubMedPubMed CentralGoogle Scholar
- Danneels EL, Rivers DB, De Graaf DC. Venom proteins of the parasitoid wasp Nasonia vitripennis: recent discovery of an untapped pharmacopee. Toxins. 2010;2:494–516.View ArticlePubMedPubMed CentralGoogle Scholar
- Ratcliffe NA, King PE. The “venom” system of Nasonia vitripennis (Walker) (Hymenoptera: Pteromalidae). Proceedings of the Royal Entomological Society of London Series A, General Entomology Blackwell Publishing Ltd. 1967;42:49–61.View ArticleGoogle Scholar
- De Graaf DC, Aerts M, Brunain M, Desjardins CA, Jacobs FJ, Werren JH, et al. Insights into the venom composition of the ectoparasitoid wasp Nasonia vitripennis from bioinformatic and proteomic studies. Insect Mol Biol. 2010;19:11–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Ross J, Jiang H, Kanost MR, Wang Y. Serine proteases and their homologs in the Drosophila melanogaster genome: an initial analysis of sequence conservation and phylogenetic relationships. Gene. 2003;304:117–31.View ArticlePubMedGoogle Scholar
- Rivers DB, Denlinger DL. Redirection of metabolism in the flesh fly, Sarcophaga bullata, following envenomation by the ectoparasitoid Nasonia vitripennis and correlation of metabolic effects with the diapause status of the host. J Insect Physiol. 1994;40:207–15.View ArticleGoogle Scholar
- Rivers DB, Ruggiero L, Hayes M. The ectoparasitic wasp Nasonia vitripennis (Walker)(Hymenoptera: Pteromalidae) differentially affects cells mediating the immune response of its flesh fly host, Sarcophaga bullata Parker (Diptera: Sarcophagidae). J Insect Physiol. 2002;48:1053–64.View ArticlePubMedGoogle Scholar
- Price DRG, Bell HA, Hinchliffe G, Fitches E, Weaver R, Gatehouse JA. A venom metalloproteinase from the parasitic wasp Eulophus pennicornis is toxic towards its host, tomato moth (Lacanobia oleracae). Insect Mol Biol. 2009;18:195–202.View ArticlePubMedGoogle Scholar
- Rivers DB, Rocco MM, Frayha AR. Venom from the ectoparasitic wasp Nasonia vitripennis increases Na+ influx and activates phospholipase C and phospholipase A 2 dependent signal transduction pathways in cultured insect cells. Toxicon. 2002;40:9-21.Google Scholar
- Miyoshi T, Tsuji N, Islam KM, Kamio T, Fujisaki K. Enzymatic characterization of a cubilin-related serine proteinase from the hard tick Haemaphysalis longicornis. J Vet Med Sci. 2004;66:1195–8.View ArticlePubMedGoogle Scholar
- Nakamatsu Y, Tanaka T. The function of a trypsin-like enzyme in the saliva of Euplectrus separatae larvae. J Insect Physiol. 2004;50:847–54.View ArticlePubMedGoogle Scholar
- Abt M, Rivers DB. Characterization of phenoloxidase activity in venom from the ectoparasitoid Nasonia vitripennis (Walker) (Hymenoptera: Pteromalidae). J Invertebr Pathol. 2007;94:108–18.View ArticlePubMedGoogle Scholar
- Hayakawa Y, Yazaki K. Envelope protein of parasitic wasp symbiont virus, polydnavirus, protects the wasp eggs from cellular immune reactions by the host insect. Eur J Biochem. 1997;246:820–6.View ArticlePubMedGoogle Scholar
- Huw Davies D, Strand MR, Vinson SB. Changes in differential haemocyte count and in vitro behaviour of plasmatocytes from host Heliothis virescens caused by Campoletis sonorensis polydnavirus. J Insect Physiol. 1987;33:143–53.View ArticleGoogle Scholar
- Labrosse C, Carton Y, Dubuffet A, Drezen JM, Poirie M. Active suppression of D. melanogaster immune response by long gland products of the parasitic wasp Leptopilina boulardi. J. Insect Physiol. 2003;49:513–22.View ArticleGoogle Scholar
- Asgari S, Hellers M, Schmidt O. Host haemocyte inactivation by an insect parasitoid: transient expression of a polydnavirus gene. J Gen Virol. 1996;77(Pt 10):2653–62.View ArticlePubMedGoogle Scholar
- Shelby KS, Webb BA. Polydnavirus-mediated suppression of insect immunity. J Insect Physiol. 1999;45:507–14.View ArticlePubMedGoogle Scholar
- Oliveira DCSG, Hunter WB, Ng J, Desjardins CA, Dang PM, Werren JH. Data mining cDNAs reveals three new single stranded RNA viruses in Nasonia (Hymenoptera: Pteromalidae). Insect Mol Biol. 2010;19 Suppl 1:99–107.View ArticlePubMedPubMed CentralGoogle Scholar
- Vavre F, Fleury F, Lepetit D, Fouillet P, Boulétreau M. Phylogenetic evidence for horizontal transmission of Wolbachia in host-parasitoid associations. Mol Biol Evol. 1999;16:1711–23.View ArticlePubMedGoogle Scholar
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.View ArticlePubMedPubMed CentralGoogle Scholar
- Mesaeli N, Nakamura K, Zvaritch E, Dickie P, Dziak E, Krause KH, et al. Calreticulin is essential for cardiac development. J Cell Biol. 1999;144:857–68.View ArticlePubMedPubMed CentralGoogle Scholar
- Pockley AG. Heat shock proteins as regulators of the immune response. Lancet. 2003;362:469–76.View ArticlePubMedGoogle Scholar
- Rivers DB, Brogan A. Venom glands from the ectoparasitoid Nasonia vitripennis (Walker)(Hymenoptera: Pteromalidae) produce a calreticulin-like protein that functions in developmental arrest and cell death in the flesh fly host, Sarcophaga bullata Parker (Diptera: Sarcophagidae). In: Maes RP, editor. Insect Physiology: New Research. New York: Nova; 2008. p. 259–78.Google Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Page MJ, Di Cera E. Evolution of peptidase diversity. J Biol Chem. 2008;283:30010–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Cônsoli FL, Tian H-S, Vinson SB, Coates CJ. Differential gene expression during wing morph differentiation of the ectoparasitoid Melittobia digitata (Hym., Eulophidae). Comp Biochem Physiol A Mol Integr Physiol. 2004;138:229–39.View ArticlePubMedGoogle Scholar
- Li LY, Luo X, Wang X. Endonuclease G is an apoptotic DNase when released from mitochondria. Nature. 2001;412:95–9.View ArticlePubMedGoogle Scholar
- Falabella PE, et al. A gamma-glutamyl transpeptidase of Aphidius ervi venom induces apoptosis in the ovaries of host aphids. Insect Biochem Mol Bio. 2007;37(5):453–65.View ArticleGoogle Scholar
- Cerstiaens A, Huybrechts J, Kotanen S, Lebeau I, Meylaers K, De Loof A, et al. Neurotoxic and neurobehavioral effects of kynurenines in adult insects. Biochem Biophys Res Commun. 2003;312:1171–7.View ArticlePubMedGoogle Scholar
- Fang J, Han Q, Li J. Isolation, characterization, and functional expression of kynurenine aminotransferase cDNA from the yellow fever mosquito, Aedes aegypti1. Insect Biochem Mol Biol. 2002;32:943–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Sloth Andersen A, Hertz Hansen P, Schaffer L, Kristensen C. A new secreted insect protein belonging to the immunoglobulin superfamily binds insulin and related peptides and inhibits their activities. J Biol Chem. 2000;275:16948–53.View ArticlePubMedGoogle Scholar
- Sarkis JJ, Guimarães JA, Ribeiro JM. Salivary apyrase of Rhodnius prolixus. Kinetics and purification Biochem J. 1986;233:885–91.PubMedGoogle Scholar
- Kim Y-S, Ryu J-H, Han S-J, Choi K-H, Nam K-B, Jang I-H, et al. Gram-negative bacteria-binding protein, a pattern recognition receptor for lipopolysaccharide and β-1,3-glucan that mediates the signaling for the induction of innate immune genes in drosophila melanogaster cells. J Biol Chem. 2000;275:32721–7.View ArticlePubMedGoogle Scholar
- Cutler JR. The Morphology of the head and the final instar lavra of Nasonia vitripennis Walker (Hymenoptera: Chalicdoidea). Proceedings of the Royal Entomological Society of London. Series A. General Entomology. 1955;30:73–81.Google Scholar
- Van Vaerenbergh M, Cardoen D, Formesyn EM, Brunain M, Van Driessche G, Blank S, et al. Extending the honey bee venome with the antimicrobial peptide apidaecin and a protein resembling wasp antigen 5. Insect Mol Biol. 2013;22:199–210.View ArticlePubMedGoogle Scholar
- Parkinson NM, Conyers C, Keen J, MacNicoll A, Smith I, Audsley N, et al. Towards a comprehensive view of the primary structure of venom proteins from the parasitoid wasp Pimpla hypochondriaca. Insect Biochem Mol Biol. 2004;34:565–71.View ArticlePubMedGoogle Scholar
- Huson DH, Auch AF, Qi J, Schuster SC. MEGAN analysis of metagenomic data. Genome Res. 2007;17:377–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.View ArticlePubMedGoogle Scholar
- Sironi M, Bandi C, Sacchi L, Di Sacco B, Damiani G, Genchi C. Molecular evidence for a close relative of the arthropod endosymbiont Wolbachia in a filarial worm. Mol Biochem Parasitol. 1995;74:223–7.View ArticlePubMedGoogle Scholar
- Casewell NR, Wagstaff SC, Wüster W, Cook DAN, Bolton FMS, King SI, Davinia P, Sanz L, Calvete JJ, Harrison RA. Medically important differences in snake venom composition are dictated by distinct postgenomic mechanisms. Proc Natl Acad Sci U S A. 2014;111:9205–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Durban J, Pérez A, Sanz L, Gómez A, Bonilla F, Rodríguez S, Chacón D, Sasa M, Angulo Y, Gutiérrez JM, Calvete JJ. Integrated “omics” profiling indicates that miRNAs are modulators of the ontogenetic venom composition shift in the Central American rattlesnake, Crotalus simus simus. BMC Genomics. 2013;14:234.View ArticlePubMedPubMed CentralGoogle Scholar
- Rokyta DR, Margres MJ, Calvin K. Post-transcriptional Mechanisms Contribute Little to Phenotypic Variation in Snake Venoms. G3. 2015;5:2375–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Fox JW. A brief review of the scientific history of several lesser-known snake venom proteins: l-amino acid oxidases, hyaluronidases and phosphodiesterases. Toxicon. 2013;62:75–82.View ArticlePubMedGoogle Scholar
- Formesyn EM, Heyninck K, de Graaf DC. The role of serine- and metalloproteases in Nasonia vitripennis venom in cell death related processes towards a Spodoptera frugiperda Sf21 cell line. J Insect Physiol. 2013;59:795–803.View ArticlePubMedGoogle Scholar
- Rivers DB, Keefer DA, Ergin E, Uçkan F. Morphology and Ultrastructure of Brain Tissue and Fat Body from the Flesh Fly, Sarcophaga bullata Parker (Diptera: Sarcophagidae), Envenomated by the Ectoparasitic Wasp Nasonia vitripennis (Walker) (Hymenoptera: Pteromalidae). Psyche 2011; Article ID 520875.Google Scholar
- Kordis D, Gubensek F. Adaptive evolution of animal toxin multigene families. Gene. 2000;261:43–52.View ArticlePubMedGoogle Scholar
- Jia LG, Shimokawa K, Bjarnason JB, Fox JW. Snake venom metalloproteinases: structure, function and relationship to the ADAMs family of proteins. Toxicon. 1996;34:1269–76.View ArticlePubMedGoogle Scholar
- Vetter I, Davis JL, Rash LD, Anangi R, Mobli M, Alewood PF, Lewis RL, King GF. Venomics: a new paradigm for natural products-based drug discovery. Amino Acids. 2011;40(1):15–28.View ArticlePubMedGoogle Scholar
- Heinen TE, da Veiga ABG. Arthropod venoms and cancer. Toxicon. 2011;57(4):497–511.View ArticlePubMedGoogle Scholar
- Zambelli VO, Pasqualoto KFM, Picolo G, Chudzinski-Tavassi AM, Cury Y. Harnessing the knowledge of animal toxins to generate drugs. Pharmacol Res. 2016. doi:10.1016/j.phrs.2016.01.009.PubMedGoogle Scholar
- Aronesty E. Comparison of Sequencing Utility Programs. Open Biotechnol J. 2013;7:1–8.View ArticleGoogle Scholar
- Cox MP, Peterson DA, Biggs PJ. SolexaQA: At-a-glance quality assessment of Illumina second-generation sequencing data. BMC Bioinformatics. 2010;11:485.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.View ArticlePubMedPubMed CentralGoogle Scholar
- Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.View ArticlePubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.View ArticlePubMedPubMed CentralGoogle Scholar
- Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36:3420–35.View ArticlePubMedPubMed CentralGoogle Scholar