- Research article
- Open Access
Venomics of the ectoparasitoid wasp Bracon nigricans
BMC Genomics volume 21, Article number: 34 (2020)
Venom is one of the most important sources of regulation factors used by parasitic Hymenoptera to redirect host physiology in favour of the developing offspring. This has stimulated a number of studies, both at functional and “omics” level, which, however, are still quite limited for ectophagous parasitoids that permanently paralyze and suppress their victims (i.e., idiobiont parasitoids).
Here we present a combined transcriptomic and proteomic study of the venom of the generalist idiobiont wasp Bracon nigricans, an ectophagous larval parasitoid of different lepidopteran species, for which we recently described the host regulation strategy and the functional role of the venom in the induction of physiological changes in parasitized hosts. The experimental approach used led to the identification of the main components of B. nigricans venom involved in host regulation. Enzymes degrading lipids, proteins and carbohydrates are likely involved in the mobilization of storage nutrients from the fat body and may concurrently be responsible for the release of neurotoxic fatty acids inducing paralysis, and for the modulation of host immune responses.
The present work contributes to fill the gap of knowledge on venom composition in ectoparasitoid wasps, and, along with our previous physiological study on this species, provides the foundation on which to develop a functional model of host regulation, based both on physiological and molecular data. This paves the way towards a better understanding of parasitism evolution in the basal lineages of Hymenoptera and to the possible exploitation of venom as source of bioinsecticidal molecules.
Parasitic Hymenoptera are part of one of the most speciose insect orders, which includes the largest number of insect natural enemies . This group of insects exhibit an impressive diversity of adaptive strategies to regulate host physiology . More basal lineages include ectoparasitoids, whose larvae feed externally on the host’s body . In most cases, these wasps have a fairly broad host range and are idiobionts, since they induce a rapid paralysis of the host which is quickly exploited as a static source of nutrients [4, 5]. In contrast, koinobionts, which in most cases develop as endophagous parasitoids (i.e., inside the body of the host), show more complex and subtle host regulation strategies, allowing a prolonged interaction of its juveniles with a living host [2, 6]. However, in both cases the growth and development of parasitoids is dependent on the physiological regulation of the host, mediated by a wide range of parasitic factors, some of which are present both in ectoparasitic and endoparasitic wasps (i.e., venom and larval secretions), while others only occur in endoparasitoids (e.g., polydnaviruses and teratocytes) [2, 7,8,9,10].
The most common and widely studied source of host regulation factors is the venom: a complex blend of proteinaceous and non-proteinaceous compounds injected at the oviposition in the host by the parasitoid female . The role and composition of venom reflect the different lifestyles of ecto- and endoparasitic wasps. In general, the venom from ectoparasitoids is involved in the rapid host immobilization, to facilitate food uptake by their larvae, while venom from endoparasitoids triggers a very diverse set of alterations, interfering with the host immune system and development or synergizing the effects of other maternal factors introduced into the host (e.g., polydnaviruses) .
The large number of studies on venom of endoparasitoid wasps has led to the identification and functional characterization of several molecules involved in the host regulation [12,13,14,15,16,17,18,19]. More recently, the advent of high-throughput technologies greatly contributed to this research area through a “multi-omic” approach often denoted as venomics, which is the integration of genomics, transcriptomics and proteomics [20,21,22,23,24,25,26,27]. Such an integrated approach provides a remarkable amount of molecular information and paves the way for the identification and exploitation of new biomolecules potentially useful for therapeutic and agricultural applications .
In contrast, only few venom components have been identified from a limited number of ectoparasitoids, despite the abundance of known species . A significant part of the research efforts are focused on host paralysis, which is the most evident and dramatic symptom induced by the venom of ectoparasitic wasps. The venom from Ampulex compressa contains components with both pre- and post-synaptic effects on GABA-gated chloride channels, determining host paralysis . The venom of Philanthus triangulum contains neurotoxic compounds (philantotoxins) which affect both the central and the peripheral nervous system, blocking the neuromuscular transmission . Envenomation by Liris niger causes host paralysis, due to blockage of synaptic transmission, but the venom components are still uncharacterized . The venoms of Bracon (=Habrobracon; =Microbracon) hebetor and Eupelmus orientalis cause a permanent paralysis, likely triggered by the neurotoxic activity of phospholipases [32,33,34].
Studies on other host regulation properties of ectoparasitoid venom are limited, with the exception of the pupal ectoparasitoid Nasonia vitripennis, for which a more comprehensive analysis of host alterations has been carried out, taking into consideration developmental arrest, metabolic changes and immunosuppression [35,36,37]. The main components of this venom  include, among others, serine- and metalloproteases, triggering apoptotic processes in Spodoptera frugiperda Sf21 cell line , two serpins and another protease inhibitor, interfering with prophenoloxidase activation in the host Musca domestica [40, 41], and a chitinase, inducing an upregulation of host genes involved in the immune response against fungi . Additional functional studies will be likely fostered in this research area since N. vitripennis is a powerful model system, for which the genome sequence and molecular tools are available [43, 44].
High-throughput technologies contributed to explore, to a limited extent, the venom composition of other ectoparasitic wasps, such as Anisopteromalus calandrae , Scleroderma guani  and B. hebetor . This latter generalist species that parasitizes a number of moth larvae was one of the first studied for its venom composition, which includes neurotoxic proteinaceus components, only partially characterized [33, 48,49,50,51].
Here we present a venomic study on a congeneric species, the generalist idiobiont wasp Bracon nigricans, which, like B. hebetor, is an ectophagous larval parasitoid of several lepidopteran species [52, 53]. This offers the opportunity to compare the molecular toolkit used by two related generalist wasp species, attacking lepidopteran hosts. The integration of transcriptomics and proteomics aims at identifying the venom components of B. nigricans responsible for the observed effects on the laboratory host Spodoptera littoralis, which include a reduced immune competence and fat body degradation associated with an enhanced host nutritional suitability for the feeding wasp larvae . The present study sheds light on the identity of venom molecules involved in the host regulation strategy adopted by B. nigricans, providing (i) the information required to develop the molecular details of the physiological model of host-parasitoid interaction recently proposed , (ii) evolutionary insights on venom blend diversification in related species with a similar host range, (iii) the opportunity to exploit venom components as potential bioinsecticides.
Assembly and analysis of venom glands transcriptome
In order to obtain the transcriptome profile of B. nigricans venom glands (Fig. 1), a cDNA library was constructed and sequenced adopting the Illumina paired-end reads sequencing. The sequencing output consisted of 25,252,591 read pairs that were reduced to 24,437,756 pairs and 796,318 single reads after adapter removal, trimming and quality check. De novo assembly of processed reads by Trinity software resulted in a total of 42,334 transcripts, with their length ranging from 201 to 29,885 bp, and a mean assembled length of 1206.93 bp (N50 = 2636 bp). The Trinity assembly output specifically consisted of 25,782 unigenes, each one representing a set of transcripts from the same locus. Main results and features of the assembly are presented in Table 1.
BLASTx similarity search results revealed that about 52.5% of the total assembled transcripts (22,218 sequences) have at least one match with the UniProtKB database. Among these sequences, 520 transcripts show similarities with venom and toxin related proteins deposited in the Tox Prot database. In particular, 455 and 65 transcripts resulted similar to putative venom proteins and putative toxin proteins, respectively.
The remaining 20,116 (47.5%) transcripts did not match with any other sequence in the UniProtKB database, and 2% of them (387 transcripts) showed a putative signal peptide. The functional annotation of the transcripts, performed using BLAST2GO, shows that 22,222 sequences (52.48% of the total assembly) share significant similarity to proteins with assigned Gene Ontology (GO) terms. These sequences are classified, at level two, into 37 functional subcategories belonging to the three main ontological categories: biological process, cellular component and molecular function. “Metabolic process” and “cellular process” are the dominant GO terms in the biological process category (Fig. 2). Highly represented terms within cellular component are “cell part” and “organelle”, whereas in the molecular function category “catalytic activity” and “binding” are the most represented terms (Fig. 2).
The SignalP analysis led to the identification of 3089 transcripts encoding proteins with a possible signal peptide at the N-terminus, suggesting that the encoded products are destined to the secretory pathway. The overall distribution of transcript abundance, expressed in RPKM, is shown in Fig. S1 (Additional file 1) and the 49 most abundant (RPKM> 100) transcripts corresponding to this subset, are listed in Additional file 2: Table S1. However, this information (i.e., the presence of signal peptide and the transcript abundance) has been subsequently integrated with the proteomic analysis, as described below.
Identification of venom proteins by proteomic approach
Upon separation of B. nigricans venom proteins by SDS-PAGE, bands with apparent molecular masses ranging from 15 to 300 kDa were observed. Figure 3 shows the 27 selected protein bands and the relative range of molecular weight (apparent molecular mass). The analyzed bands were resolved into 109 proteins by LC/MS-MS (Additional file 3: Table S2), identified by matching the resulting peptides to the full-length proteins predicted from the venom gland transcriptome assembly. Matches of at least two peptides were considered valid.
Integration of transcriptomic and proteomic analyses
To identify a robust set of abundant venom proteins, all data were analyzed under the assumption that venom proteins are secreted and highly expressed in the venom glands. This was performed by filtering the proteome for the presence of the signal peptide identified by SignalP tool (approximately 22% of the 109 total proteins) and, among these, we focused on those in the upper third of RPKM values distribution, exceeding the threshold value of 70. This allowed to generate a list of 18 venom proteins highly abundant and, as such, with important functions (Table 2 and Additional file 4: Table S3). Only 10 of these candidates showed similarity with protein sequences deposited in the UniProt/SwissProt database.
Specificity of expression in the venom gland of selected genes
To experimentally corroborate the validity of the approach used in the identification of the major functional players in the venom blend of B. nigricans, a qRT-PCR experiment was performed on total RNA extracted from venom glands, whole adult males and females deprived of the venom apparatus, focusing on a sample set of 8 genes encoding proteins selected among those showing similarity with sequences deposited in the UniProt/SwissProt database and/or highly represented in the venom gland transcriptome. Female body devoid of venom glands was used as calibrator sample. A significantly higher transcription level in venom glands was observed for all genes considered (Additional file 5: Table S4): phospholipase A2 (P < 0.001), lipase (P < 0.001), odorant binding protein (P < 0.001), trypsin-like serine protease (P < 0.001), carboxylesterase (P < 0.001), aminopeptidase (P < 0.001), protein disulfide-isomerase (P < 0.001) and lysosomal alpha-mannosidase (P < 0.005) (Fig. 4). Moreover, a much lower but significantly higher transcription level was found also in males, compared to females devoid of venom glands, for phospholipase A2 (P < 0.005), protein disulfide-isomerase (P < 0.005) and lysosomal alpha-mannosidase (P < 0.05) (Fig. 4).
Female parasitoid wasps inject complex cocktails of proteins into their hosts to induce physiological alterations which facilitate the successful development of their offspring . B. nigricans is an ectoparasitic wasp which determines a rapid host developmental arrest by envenomation. Venom injection immobilizes the host, inhibits immune response and triggers metabolism redirection, in order to enhance host nutritional suitability for the developing parasitoid larvae . To identify the major venom components involved in the host regulation, we performed an in-depth venomics approach, which combines high-throughput transcriptomics and proteomics .
Concurrent analysis of venom gland transcriptome and venom proteome revealed the presence of a large proportion of sequences (47.5% of the total transcripts) showing no similarity with those available in the UniProtKB database. Since parasitism factors can be very specific, displaying poor similarity among related evolutionary lineages [57,58,59], these unannotated proteins very likely include bioactive venom components that, however, are not easy to characterize from a functional point of view, even though they represent an untapped source of biological novelties worth of further studies.
Both transcriptome and proteome reported the presence of membrane, cytoskeletal or mitochondrial proteins, lacking the predicted signal peptide (e.g., myosin, calcium-transporting ATPase, tubulin). While these genes are expected to be found in the transcriptome of venom glands, their presence in the venom proteome very likely derives from a contamination due to cell breakage during sample collection .
Venom includes secreted proteins, which are expressed at different levels, and their abundance is often associated with important biological functions [15, 60]. In order to identify a robust set of proteins highly expressed in venom glands, we combined transcriptomic and proteomic data. In particular, the criteria for transcript selection were: (1) the presence of peptides identified by proteomic analysis; (2) the abundance at the transcript level; (3) the homology with known effectors in parasitoid-host interaction studies. This approach allowed the identification of a set of genes, very likely having an important role in the host regulation, which were subjected to qRT-PCR analysis, to assess the degree of specificity of their expression in the venom glands of B. nigricans. All these genes were highly expressed in venom glands, to a various extent, confirming their key-relevance in the venom blend.
Multiple alignments revealed that most venom proteins of B. nigricans are putatively functional enzymes, sharing the same conserved active sites of their homologs expressed in species lacking venom glands (Additional file 6: Figures S2-S9). Indeed, as already suggested for parasitic Hymenoptera, the mechanism of acquisition of venom components is likely related to co-option (i.e., new use of existing genes through expression alterations) rather than to duplication and neofunctionalization . However, most of venom proteins identified in this study clearly diverge from homologs present both in non-venomous species, and in species using venom for defense (e.g. ants, bees) (Additional file 7: Figures S10-S17); this suggests the occurrence of a divergent evolution driven by colonization of separate host niches and specialization to different parasitic life habits.
The main components of B. nigricans venom are hydrolases acting on ester bonds (EC 3.1), which include a large group of different proteins, such as phospholipases A2 (PLA2) and lipases, frequently reported as venom components of several parasitoid species [55, 62, 63].
The PLA2 (BnPLA2; comp22364_c1_seq2) identified by our analysis resulted highly expressed in venom glands, in accordance with the wide occurrence of this enzyme, which is one of the main venom components of Hymenoptera, including honeybees  and many parasitic wasps, such as E. orientalis , Toxoneuron nigriceps  and Psyttalia concolor . BnPLA2 shows 48% of identity, along 99% of the protein length, with two putative PLA2s of B. hebetor found by BLASTp in NCBI patented protein sequences database (GenBank: CAB42203.1; CAA03259.1). In B. hebetor the two sequenced PLA2s, named BrhTX-1(b) and BrhTX-1(c), are two of the four subunits of the paralyzing BrhTX-1 toxin . PLA2s (EC 126.96.36.199) are a large super-family of lipolytic enzymes that cleave the glycerol backbone of phospholipids, usually in a metal-dependent reaction, to release free fatty acids and lysophospholipids . These products of lipid hydrolysis (i.e., arachidonic and linoleic acid) are cytotoxic , as demonstrated by many studies on catalytically active venom PLA2s exerting neurotoxic [68, 69], myotoxic , anticoagulant  and antibacterial activities . A pharmacological effect independent of catalytic activity is also predicted for divergent putative homologs of BnPLA2 in other Hymenoptera, such as A. compressa, which lacks the catalytic site (H34) (Additional file 6: Figure S2 and Additional file 7: Figure S10) [73, 74].
Among transcripts highly expressed in venom glands, a lipase (BnLIP; comp22420_c0_seq1) exhibits 32% of identity, along 90% of the protein length, with a homolog found in N. vitripennis venom (GenBank: NP_001154991.1) . Lipases (EC 188.8.131.52) are also found in the venom of Pteromalus puparum , Pimpla hypochondriaca , Chelonus inanitus , Microctonus aethiopoides  and Ooencyrtus telenomicida . These enzymes perform essential roles in the digestion, transport and processing of dietary lipids in most living organisms and might participate in the breakdown of the energy stores contained in the fat body of the host, in order to increase the nutritional suitability of its body fluids ingested by parasitoid larvae [79, 80]. Moreover, these enzymes may contribute by lipid hydrolysis to the generation of the toxic molecules mentioned above.
A similar function is hypothesized for another esterase, a carboxylesterase (comp24797_c0_seq1), which results highly expressed in the venom glands of B. nigricans. Carboxylesterases (EC 184.108.40.206) are reported as components of diverse venoms of parasitoids, including the venoms of N. vitripennis , Hyposoter didymator  and A. calandrae . Interestingly, a putative carboxylesterase was isolated from the teratocytes of Dinocampus coccinellae, which are cells of embryonic origin that support the development of parasitoid larvae . Considering that de novo lipid synthesis is energetically expensive , the mobilization and consumption of host lipids, through the action of the abovementioned lipolytic venom hydrolases, could provide a selective advantage for B. nigricans larvae by enhancing the host nutritional suitability.
A trypsin-like serine protease (BnTRY; comp16818_c0_seq2; EC:3.4.21) is highly expressed in venom glands, compared to other tissues. BnTRY shows 52% of identity, along 100% of the protein length, with a putative trypsin-like serine protease of B. hebetor venom found by BLASTp in NCBI patented protein sequences database (GenBank: CAB42201.1) which corresponds to “Sequence 7 from Patent WO9744355” (GenBank: A67382.1) . Moreover, this protein exhibits a 23% of identity, along 90% of the sequence length, with a N. vitripennis venom serine protease (NP_001155042.1), which is the annotated gene in the venom glands of this pupal ectoparasitoid showing the highest expression level . Serine proteases are a very common functional category in insect genomes that can have very diverse roles in parasitoid physiology. A classical function exerted by serine proteases is digestion, as occurs in the larvae of the ectoparasitoid Euplectrus separatae, which releases a salivary secretion containing a trypsin-like enzyme to digest the host tissues . When highly expressed in venom glands, serine proteases may also play important roles in interfering with the immune response of the host by altering the proteolytic cascades activated by the detection of non-self intruders [85, 86]. Indeed, serine proteases were found in the venom of various parasitoid species, such as, for example, Aphidius ervi , C. inanitus , N. vitripennis , P. hypochondriaca , and P. puparum . However, the disruption of the proteolytic activating cascade involved in the melanization of host hemolymph can be also induced by a mutated serine protease with inhibitory effects, such as reported for the endoparasitoid Cotesia rubecula . Interestingly, B. nigricans envenomation leads to a reduced encapsulation and melanization response in the host S. littoralis , supporting the hypothesis that BnTRY can be involved in the immunosuppressive syndrome.
A leucyl-cystinil aminopeptidase (BnLCA; comp13820_c0_seq1; EC 220.127.116.11) shows 24% identity, along 98% of the sequence length, with an aminopeptidase N-like protein (GenBank: EFN65598.1) from Camponotus floridanus (Hymenoptera, Formicidae). BnLCA contains a peptidase_M1 domain (pfam01433), which characterizes a family of zinc-metalloenzymes involved in the cleavage of amino acids from the amino terminus (N-terminus) of proteins or peptides . Aminopeptidase activity was detected in the venom of Cotesia chilonis , P. hypochondriaca  and Psyttalia lounsboury , although no further functional data are available in the literature. On the contrary, the aminopeptidase N of Plasmodium falciparum, A-M1 (PfA-M1) is well characterized. PfA-M1 has various localizations and, in particular, in the food vacuole (a unique proteolytic organelle), where is involved in parasite metabolism and participates in the last steps of hemoglobin degradation . Thus, it is reasonable to assume that BnLCA may participate to the general degradation of host tissues, perhaps increasing their permeability to other venom components and facilitating nutritional exploitation of the host, as suggested for other venom zinc-metallopeptidases [92, 93].
A lysosomal alpha-mannosidase (BnLAM; comp22165_c0_seq2; EC 18.104.22.168), which takes part in the sequential degradation of complex, hybrid and high-mannose N-linked oligosaccharides , is highly expressed in venom glands of B. nigricans. A homolog of BnLAM has been recently identified by proteomic approach in the venom of an encyrtid endoparasitoid wasp, Diversinervus elegans . However, the role of lysosomal alpha-mannosidase (LAM) and, more in general, the role of carbohydrases in entomophagous insects has been poorly characterized. LAMs, as well as the abovementioned trypsin-like serine protease and leucyl-cystinil aminopeptidase, are present in the digestive fluids of entomophagous arachnids (scorpions and spiders), which feed on their insect preys through extra-oral digestion (EOD) [95,96,97]. Enzymatic digestion through the action of saliva and digestive fluids represents a common strategy adopted by predatory arthropods and parasitoid wasps to exploit nutritional resources [98, 99]. In parasitic Hymenoptera, the EOD has been described in different species, as, for example, in the case of the larvae of Trichogramma australicum and E. separatae [84, 100], as well as for the teratocytes of A. ervi, which determine the cytolytic degradation of the formed embryos of the aphid Acyrthosiphon pisum . However, as said above, these digestive enzymes identified in the venom of B. nigricans can be likely involved also in a more subtle regulation of the host immunity. Indeed, it is reasonable to propose that this enzyme could induce the structural disruption of glycoconjugants which mediate the recognition of altered self domains associated with feeding wounds . Then, by preventing the recruitment of hemocytes for wound healing, the feeding hole remains patent and the food uptake by wasp larvae is not impaired, as needed for all ectophagous parasitic arthropods .
A pheromone/general odorant binding protein (BnOBP; comp24796_c0_seq1; IPR006170) results one of the most abundant components in the venom of B. nigricans, as suggested by the strong intensity of the corresponding band (number 27, Fig. 3) in SDS-PAGE, high RPKM (the highest in the transcriptome) and qRT-PCR data. BnOBP, which has no matches in SwissProt database, shows 27% identity, along the 62% of the sequence length, with the translated transcript GECT01010095 (GenBank) from the endoparasitoid wasp P. puparum. OBPs participate in solubilization and transport of small hydrophobic odorant molecules and pheromones, and are characterized by six conserved cysteine residues, forming three disulfide bonds that stabilize the folded structure of the protein . BnOBP contains only four cysteine residues in the mature protein sequence (Additional file 6: Figure S8) and thus can be classified as Minus-C OBP (i.e., less than six cysteine residues), tough the relationship between function and structural changes for this protein family is still largely unknown . It has been suggested that with less cysteine residues, Minus-C OBPs might have more structural flexibility than the Classic OBPs , and thus higher binding affinity for a broader number of ligands . Phylogenetic reconstruction (Additional file 7: Figure S16) showed the presence of several paralogs in the venom gland transcriptome of B. nigricans, while few putative orthologs are found in the nr NCBI db, suggesting that OBPs are evolving through a rapid neofunctionalization . The presence of members of this protein family in the venoms of parasitoid wasps has been reported for N. vitripennis , C. inanitus , Leptopilina heterotoma , P. puparum  and A. calandrae . As hypothesized by these studies, the detection of OBPs in the venom may suggest their involvement in host selection [27, 110], when, after initial paralysis, B. nigricans female intensively probes the host with its ovipositor before egg laying . Alternatively, beyond their canonical chemosensory role, it has been proposed that OBPs can act as carriers, mediating the solubilization of hydrophobic molecules, such as free fatty acids released by lipases [111, 112]. Then, BnOBP may have a nutritional relevance for the developing parasitoid larvae acting in tandem with BnLIP and, thus, contributing to lipid mobilization in the parasitized host.
A protein disulfide-isomerase (BnPDI; comp21345_c1_seq3; EC 22.214.171.124) shows 73.8% identity, along the 100% of the sequence length, with a recently annotated PDI from the genome of the parasitoid wasp Diachasma alloeum (NCBI Reference Sequence: XP_015122546.1). This enzyme family has also been identified in the venom glands of A. ervi  and Psytallia species , as well as in the crude venom extract of P. puparum , D. elegans  and C. chilonis . PDIs are involved in the folding of proteins by catalysing the oxidation, isomerization and reduction of disulfide bonds, that covalently link specific cysteine residues and confer stability to proteins [113, 114]. In venomous cone snails, PDIs are strictly localized in the venom glands, where they guide the folding of cysteine-rich peptide toxins (conotoxins) into their native state, while they are absent in the secreted venom [115, 116]. Indeed, PDIs are mainly localized in the endoplasmic reticulum, and their presence in the extracellular space is considered rather rare . Intriguingly, BnPDI sequence lacks the C-terminal motif Lys-Asp-Glu-Leu (KDEL) (Additional file 6: Figure S9), which prevents the secretion of luminal proteins from endoplasmic reticulum . Therefore, the PDI presence in the venom of B. nigricans, if not due to the breakage of venom glands during venom collection, could represent a new interesting case of secretion, worth of further studies.
The integration of transcriptomics and proteomics used in the present work provides the first description of the main components of B. nigricans venom and contributes to the expansion of the limited information available for the venom of ectophagous parasitoids, that permanently paralyze and rapidly suppress their victims (i.e., idiobiont parasitoids). The most abundant venom components of B. nigricans, mainly esterases and proteases likely involved in paralysis (BnPLA2) and digestion of the host’s tissues, are in line with a typical ectoparasitic idiobiont life-style. Enzymes with lipolytic activity are likely involved in the mobilization of storage nutrients from fat bodies and/or in the release of neurotoxic fatty acids inducing paralysis, while enzymes related to carbohydrate catabolism could be responsible for the alteration of glycoconjugants which mediate the recognition of altered self domains mediating wound healing, and may induce the observed changes of the carbohydrate titer in the host hemolymph . The host regulation strategy emerging from this latter study is well corroborated by the discovery of venom component that can have immunosuppressive activity, such as a trypsin-like serine protease (BnTRY) and a lysosomal alpha-mannosidase (BnLAM).
The venom gland transcriptome of B. nigricans shows the occurrence only of a few transcripts in common with B. hebetor . The most significant similarity is the presence in both species of neurotoxic molecules (i.e., phospholipases), which were previously identified in B. hebetor and patented ; only a few other transcripts are shared (arginine kinase and venom acid phosphatase). Even though the data available for B. hebetor are somewhat limited, it is quite evident that these two closely related species, attacking the same group of insect hosts, rely upon a molecular toolkit for host regulation and exploitation which is only partly conserved, suggesting the occurrence of an ancestral broad molecular biodiversity of the venom blend, which has been one of the major pre-requisites allowing the evolutionary diversification of parasitic Hymenoptera.
The identification of these venom components can have important implications from an applied perspective, given their bioinsecticide activity . Indeed, the functional characterization of the venom proteins identified in the present study will offer new tools to develop bioinspired strategies for pest control, based on the use of natural antagonists beyond the organism level, as a source of insecticide molecules.
Parasitoid and host rearing
B. nigricans was reared on larvae of the noctuid moth S. littoralis, which were used as laboratory host, maintained on artificial diet as previously described . Adults of B. nigricans were maintained at 27 ± 1 °C, 70 ± 5% RH and 16:8 h (L:D) photoperiod, as previously described . The starting material of B. nigricans colony, used for transcriptomic and proteomic studies, was kindly donated by Lucia Zappalà (University of Catania, Catania, Italy) and derived from a laboratory population established in 2010 with specimens collected in Italy, Sicily region (Catania area), to which field material from the same area was added every year to refresh the colony. S. littoralis colony was established in 2010 and refreshed with insects provided by ISAGRO (Novara, Italy).
Venom gland collection and RNA isolation
Female wasps were anaesthetized on ice and dissected in sterile PBS (Phosphate Buffered Saline: 137 mM NaCl, 2.7 mM KCl, 10 mM phosphate buffer, pH 7.4) by grasping the ovipositor tip with fine forceps, under a stereoscope (Discovery v8, Zeiss). Reservoir and venom glands (Fig. 1) were separated using dissection needles and forceps. Venom glands were collected into a microcentrifuge tube containing 300 μl of TRIzol® Reagent (Thermo Fisher Scientific, Waltham, MA, USA) kept on ice and, then, stored at − 80 °C until RNA extraction. Total RNA extraction was performed according to the manufacturer’s directions, using TRIzol® Reagent, while the PureLink™ RNA Micro Kit was used to purify and concentrate RNA (Thermo Fisher Scientific). In order to eliminate genomic DNA, an on-column DNase treatment (PureLink®DNase, Thermo Fisher Scientific) was carried out. An amount of 2.4 μg of total RNA was obtained from venom glands of 80 female wasps. The quantity and the quality of total RNA were assessed using Varioskan Flash spectrophotometer (Thermo Fisher Scientific) and, prior to cDNA library construction, with Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).
cDNA library construction, sequencing and computational analysis
TruSeq (Illumina) stranded cDNA sequencing libraries were constructed. The 100 bp paired-end sequencing run was performed on the Illumina HiSeq1500 platform. Raw sequences were trimmed from adapters using Cutadapt . The remaining sequences were checked for quality (Q20 and error rate 0.01) using Trim Galore v. 0.4.0 [122, 123]. The de novo transcriptome assembly of the cleaned reads was performed with Trinity, setting the k-mer length at the default value (25) [124, 125]. To remove redundancy, the assembled transcripts were further collapsed using CD-HIT-EST, setting identity (−c option) and coverage (−s option) equal to 1 . The sequences obtained were annotated by BLASTx comparisons versus the UniProtKB/SwissProt protein database , using default parameters and setting the e-value threshold at 1e-5. GO terms assignments, encoded enzymes and functional domain identifications were supplemented to the functional annotation by the web based pipeline FastAnnotator . Signal peptides for translated transcripts were predicted using SignalP 4.0 software .
In order to detect the presence of putative toxins in the resulting data, the protein hits associated with each transcript by the BLASTx search versus the UniProtKB database were crossed with the sequences deposited in the UniProtKB/Swiss-Prot Tox-Prot db (http://www.uniprot.org/program/Toxins), which are categorized in toxin and venom proteins . The assembled and annotated venom gland transcriptome was used to generate a custom-made protein database, by translating the six reading frames of the nucleotide sequences in their corresponding amino acid sequences by SEQtools (http://www.seqtools.dk/). The relative expression level of each transcript was estimated calculating the RPKM (Reads per Kilobase per Million) based on the number of reads re-mapped on the assembled sequences using Bowtie2 alignment tool .
Extraction and SDS-PAGE fractionation of venom proteins
Adult females of B. nigricans were dissected in ice-cold Pringle’s solution  with 1 mM phenylmethylsulfonyl fluoride, as described above. After resection from venom glands, the reservoir was gently ruptured, using dissecting needles, and the exuding crude venom content transferred to a microcentrifuge tube. The crude extract was centrifuged at 5000×g, at 4 °C, for 5 min to remove cell debris, and the supernatant was thus filtered through a 0.22 μm filter (Millipore, Merck KGaA, Germany) and stored at − 80 °C until use. Protein concentration of the venom extracted was assessed by Bradford method .
In total, 12 μg of venom proteins were fractioned on 10% SDS-PAGE minigel and stained with colloidal Coomassie (0.1% (w/v) Coomassie Brilliant Blue G250 (Thermo Fisher Scientific), 2% (w/v) ortophosphoric acid, 10% (w/v) ammonium sulfate). Selected protein bands were thus excised from the gel and washed in a solution of 50 mM ammonium bicarbonate (pH 8.0) in 50% acetonitrile until complete destaining. Each gel slice was resuspended in 50 mM ammonium bicarbonate and incubated with 100 ng of trypsin for 2 h at 4 °C and overnight at 37 °C. The supernatant containing the resulting mixture of peptides was removed and the gel pieces were re-extracted with acetonitrile. The two fractions were then pooled and freeze-dried.
LC/MS-MS and protein identification
The peptide mixture from each gel slice was analyzed by LC/MS-MS, using the LC/MSD Trap XCT Ultra (Agilent Technologies, Palo Alto, CA), equipped with a 1100 HPLC system and a chip cube (Agilent Technologies). After loading, the peptide mixture (8 μl in 0.2% (v/v) formic acid in water) was first concentrated and washed at 4 μl/min in a 40 nl enrichment column (Agilent Technologies chip), with 0.1% formic acid as eluent (eluent A). The sample was then fractionated on a C18 reverse-phase capillary column (75 μm × 43 mm in the Agilent Technologies chip) at a flow rate of 300 nl/min, with a linear gradient of eluent B (0.1% (v/v) formic acid in acetonitrile) in A, from 5 to 60% in 50 min. Elution was monitored with the mass spectrometers without any splitting device. Peptide analysis was performed using data-dependent acquisition of one MS scan (m/z range from 400 to 1600 Da/e), followed by MS/MS scans of the three most abundant ions in each MS scan. Dynamic exclusion was used to acquire a more complete survey of the peptides by automatic recognition and temporary exclusion (2 min) of ions from which definitive mass spectral data had previously been acquired. Moreover, a permanent exclusion list of the most frequent peptide contaminants (keratins and trypsin peptides) was included in the acquisition method in order to focus the analyses on significant data.
The resulting MS/MS spectra were searched against the translated B. nigriceps transcriptome with Mascot 2.4 (Matrix Science, Boston, MA, USA), using the Ion Search option. Carbamidomethyl of cysteine was set as fixed modification, while Pyro-Carbamidomethyl of cysteine (N-term C), Pyro-Glu (N-term Q) and oxidation of methionine were set as variable modifications. The other parameters were set at 0.6 Da for mass tolerance, ± 660 ppm for precursor ion tolerance and 1 for maximum missed cleavages allowed.
Quantitative RT-PCR (qRT-PCR) analysis
To confirm in vivo the result of proteomic and transcriptomic analyses, the expression profiles of selected genes were assessed using qRT-PCR. Total RNA was separately extracted from (1) venom glands, (2) females deprived of the venom glands and (3) males, using TRIzol® reagent (Thermo Fisher Scientific) according to the manufacturer’s protocol. Each sample was formed by a pool of RNA obtained from 5 individuals and was analyzed in triplicate on a Step One Real Time PCR System (Applied Biosystems, Carlsbad, CA, USA). Differential relative expression of the target genes was measured by one-step qRT-PCR, using the SYBR Green PCR Kit (Applied Biosystems), according to the manufacturer’s instructions. B. nigricans 40S ribosomal protein S3 (RPS3) gene (GenBank accession number: MK631956) was used as endogenous control for RNA loading. All primers were designed using Primer Express, version 1.0 software (Applied Biosystems) and are reported in (Additional file 8: Table S5). Relative gene expression data were analyzed using the ∆∆Ct method [134, 135]. For validation of the ∆∆Ct method, the difference between the Ct value of the target and the Ct value of RPS3 transcripts [∆Ct = Ct(target gene)-Ct (RPS3)] was plotted versus the log of ten-fold serial dilutions (100, 10, 1, 0.1 and 0.01 ng) of the purified RNA samples. The plot of log total RNA input versus ∆Ct displayed a slope less than 0.1, indicating that the efficiencies of the two amplicons were approximately equal. The relative expression of the target genes in female body deprived of the venom glands was used as calibrator (relative expression = 1). The results are presented as mean fold changes of three independent biological replicates. ∆Cts were compared using one-way analysis of variance (ANOVA) and Tukey’s test with statistical significance set at P < 0.05. All statistical analyses were performed using the Statistical Analysis Systems software (Sigma Stat Statistical Software, SPSS Science, Chicago, IL, USA).
Sequence analysis and phylogeny reconstruction
Putative homologous sequences of the most representative venom proteins were identified by sequence similarity searches through a BlastP analysis versus the non-redundant NCBI database (nr NCBI, release October 2019) and Swiss-Prot (release 2019_10 of 13-Nov-19). Additional rounds of BLASTp analyses were performed selecting Arachnida (taxid:6854) and Hymenoptera (taxid:7399) groups in the nr NCBI collection. However, using the BnOBP protein sequence as a query resulted in only 3 hits below E-value threshold of 10− 1. Therefore a tBLASTn search against the transcriptome of B. nigricans venom glands was also performed. Most representative hits selected below the E-value threshold of 10− 5 were aligned using Muscle 3.8 , with default settings.
The alignments of venomous and non-venomous orthologs were plotted using Jalview 2 , and sequences were analyzed with ScanProsite  and HMMscan (HmmerWeb version 2.40.0) tools, in order to identify active sites and conserved patterns.
To reconstruct phylogeny, alignments were manually trimmed to avoid comparisons of non-conserved regions present only in a subset of the taxa. Best-fit model of amino acid substitution and phylogenetic reconstruction was performed using RAxML 8.2.12 . The maximum-likelihood tree was run for 1000 bootstrap replicates and the tree figure was plotted using FigTree v1.4.3.
Analysis of Variance
Basic Local Alignment Search Tool
Bracon nigricans lysosomal alpha-mannosidase
Bracon nigricans leucyl-cystinil aminopeptidase
Bracon nigricans lipase
Bracon nigricans odorant binding protein
Bracon nigricans protein disulfide-isomerase
Bracon nigricans phospholipase A2
Bracon nigricans trypsin-like serine protease
High Performance Liquid Chromatography
Liquid Chromatography/Mass Selective Detector
Liquid Chromatography coupled with tandem Mass Spectrometry
Odorant binding protein
Plasmodium falciparum aminopeptidase N
Quantitative Reverse Transcriptase - Polymerase Chain Reaction
Reads Per Kilobase per Million
Ribosomal Protein S3
Sodium Dodecyl Sulphate - PolyAcrylamide Gel Electrophoresis
Trypsin-like serine protease
Forbes AA, Bagley RK, Beer MA, Hippee AC, Widmayer HA. Quantifying the unquantifiable: why Hymenoptera, not Coleoptera, is the most speciose animal order. BMC Ecol. 2018;18(1):21.
Pennacchio F, Strand MR. Evolution of developmental strategies in parasitic hymenoptera. Annu Rev Entomol. 2006;51:233–58.
Whitfield JB. Phylogenetic insights into the evolution of parasitism in Hymenoptera. Advances in Parasitology, Volume 54; London: Academic Press. 2003;p. 69–100.
Gauld ID. Evolutionary patterns of host utilization by ichneumonoid parasitoids (Hymenoptera: Ichneumonidae and Braconidae). Biol J Linn Soc. 1988;35:351–77.
Hawkins BA, Askew RR, Shaw MR. Influences of host feeding-niche and foodplant type on generalist and specialist parasitoids. Ecol Entomol. 1990;15:275–80.
Quicke DLJ. Parasitic wasps. London: Chapman & Hall; 1997.
Beckage NE. Modulation of immune responses to parasitoids by polydnaviruses. Parasitology. 1998;116(Suppl):S57–64.
Dahlman DL. Teratocytes and host/parasitoid interactions. Biol Control. 1991;1:118–26.
Fritz RS. Selection for host modification by insect parasitoids. Evolution. 1982;36:283–8.
Vinson SB. Biochemical coevolution between parasitoids and their hosts. In: Price PW, editor. Evolutionary strategies of parasitic insects and mites. Boston, MA: Springer US; 1975. p. 14–48. doi:https://0-doi-org.brum.beds.ac.uk/10.1007/978-1-4615-8732-3_2.
Moreau SJM, Asgari S. Venom proteins from parasitoid wasps and their biological functions. Toxins. 2015;7:2385–412.
Asgari S, Reineke A, Beck M, Schmidt O. Isolation and characterization of a neprilysin-like protein from Venturia canescens virus-like particles. Insect Mol Biol. 2002;11:477–85.
Colinet D, Dubuffet A, Cazes D, Moreau S, Drezen J-M, Poirié M. A serpin from the parasitoid wasp Leptopilina boulardi targets the Drosophila phenoloxidase cascade. Dev Comp Immunol. 2009;33:681–9.
Dani MP, Richards EH. Cloning and expression of the gene for an insect haemocyte anti-aggregation protein (VPr3), from the venom of the endoparasitic wasp, Pimpla hypochondriaca. Arch Insect Biochem Physiol. 2009;71:191–204.
Falabella P, Riviello L, Caccialupi P, Rossodivita T, Teresa Valente M, Luisa De Stradis M, et al. A gamma-glutamyl transpeptidase of Aphidius ervi venom induces apoptosis in the ovaries of host aphids. Insect Biochem Mol Biol. 2007;37:453–65.
Labrosse C, Eslin P, Doury G, Drezen JM, Poirié M. Haemocyte changes in D. melanogaster in response to long gland components of the parasitoid wasp Leptopilina boulardi: a rho-GAP protein as an important factor. J Insect Physiol. 2005;51:161–70.
Parkinson N, Smith I, Audsley N, Edwards JP. Purification of pimplin, a paralytic heterodimeric polypeptide from venom of the parasitoid wasp Pimpla hypochondriaca, and cloning of the cDNA encoding one of the subunits. Insect Biochem Mol Biol. 2002;32:1769–73.
Shen X, Ye G, Cheng X, Yu C, Yao H, Hu C. Novel antimicrobial peptides identified from an endoparasitic wasp cDNA library. J Pept Sci. 2010;16:58–64.
Zhang G, Schmidt O, Asgari S. A novel venom peptide from an endoparasitoid wasp is required for expression of polydnavirus genes in host hemocytes. J Biol Chem. 2004;279:41580–5.
Goecks J, Mortimer NT, Mobley JA, Bowersock GJ, Taylor J, Schlenke TA. Integrative approach reveals composition of endoparasitoid wasp venoms. PLoS One. 2013;8:e64125.
Laurino S, Grossi G, Pucci P, Flagiello A, Bufo SA, Bianco G, et al. Identification of major Toxoneuron nigriceps venom proteins using an integrated transcriptomic/proteomic approach. Insect Biochem Mol Biol. 2016;76:49–61.
Lin Z, Wang R-J, Cheng Y, Du J, Volovych O, Han L-B, et al. Insights into the venom protein components of Microplitis mediator, an endoparasitoid wasp. Insect Biochem Mol Biol. 2019;105:33–42.
Liu N-Y, Wang J-Q, Zhang Z-B, Huang J-M, Zhu J-Y. Unraveling the venom components of an encyrtid endoparasitoid wasp Diversinervus elegans. Toxicon. 2017;136:15–26.
Liu N-Y, Xu Z-W, Yan W, Ren X-M, Zhang Z-Q, Zhu J-Y. Venomics reveals novel ion transport peptide-likes (ITPLs) from the parasitoid wasp Tetrastichus brontispae. Toxicon Off J Int Soc Toxinology. 2018;141:88–93.
Mathé-Hubert H, Colinet D, Deleury E, Belghazi M, Ravallec M, Poulain J, et al. Comparative venomics of Psyttalia lounsburyi and P. concolor, two olive fruit fly parasitoids: a hypothetical role for a GH1 β-glucosidase. Sci Rep. 2016;6:35873.
Teng Z-W, Xiong S-J, Xu G, Gan S-Y, Chen X, Stanley D, et al. Protein discovery: combined transcriptomic and proteomic analyses of venom from the endoparasitoid Cotesia chilonis (Hymenoptera: Braconidae). Toxins. 2017;9:135.
Yan Z, Fang Q, Wang L, Liu J, Zhu Y, Wang F, et al. Insights into the venom composition and evolution of an endoparasitoid wasp by combining proteomic and transcriptomic analyses. Sci Rep. 2016;6:19604.
Wilson D, Daly NL. Venomics: a mini-review. High-Throughput. 2018;7(3).
Moore EL, Haspel G, Libersat F, Adams ME. Parasitoid wasp sting: a cocktail of GABA, taurine, and β-alanine opens chloride channels for central synaptic block and transient paralysis of a cockroach host. J Neurobiol. 2006;66(8):811–20.
Eldefrawi AT, Eldefrawi ME, Konno K, Mansour NA, Nakanishi K, Oltz E, et al. Structure and synthesis of a potent glutamate receptor antagonist in wasp venom. Proc Natl Acad Sci. 1988;85(13):4910–3.
Ferber M, Hörner M, Cepok S, Gnatzy W. Digger wasp versus cricket: mechanisms underlying the total paralysis caused by the predator’s venom. J Neurobiol. 2001;47(3):207–22.
Doury G, Bigot Y, Periquet G. Physiological and biochemical analysis of factors in the female venom gland and larval salivary secretions of the ectoparasitoid wasp Eupelmus orientalis. J Insect Physiol. 1997;43(1):69–81.
Quistad GB, Nguyen Q, Bernasconi P, Leisy DJ. Purification and characterization of insecticidal toxins from venom glands of the parasitic wasp, Bracon hebetor. Insect Biochem Mol Biol. 1994;24(10):955–61.
Weaver RJ, Marris GC, Bell J, Edwards JP. Identity and mode of action of host endocrine disrupters from the venom of parasitoid wasps. Comp Biochem Physiol -- Part B Biochem Mol Biol. 2000;(Supplement 1):S101.
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(3):207–15.
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.
Rivers DB, Rocco MM, Frayha AR. Venom from the ectoparasitic wasp Nasonia vitripennis increases Na+ influx and activates phospholipase C and phospholipase A2 dependent signal transduction pathways in cultured insect cells. Toxicon. 2002;40:9–21.
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.
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.
Qian C, Fang Q, Wang L, Ye G-Y. Molecular cloning and functional studies of two Kazal-type serine protease inhibitors specifically expressed by Nasonia vitripennis venom apparatus. Toxins. 2015;7:2888–905.
Qian C, Liang D, Liu Y, Wang P, Kausar S, Wei G, et al. Identification of a small pacifastin protease inhibitor from Nasonia vitripennis venom that inhibits humoral immunity of host (Musca domestica). Toxicon Off J Int Soc Toxinology. 2017;131:54–62.
Martinson EO, Martinson VG, Edwards R, Null M, Werren JH. Laterally transferred gene recruited as a venom in parasitoid wasps. Mol Biol Evol. 2016;33:1042–52.
Martinson EO, Siebert AL, He M, Kelkar YD, Doucette LA, Werren JH. Evaluating the evolution and function of the dynamic Venom Y protein in ectoparasitoid wasps. Insect Mol Biol. 2019.
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.
Perkin LC, Friesen KS, Flinn PW, Oppert B. Venom gland components of the ectoparasitoid wasp, Anisopteromalus calandrae. J Venom Res. 2015;6:19–37.
Zhu J-Y. Deciphering the main venom components of the ectoparasitic ant-like bethylid wasp, Scleroderma guani. Toxicon. 2016;113:32–40.
Manzoor A, UlAbdin Z, Webb BA, Arif MJ, Jamil A. De novo sequencing and transcriptome analysis of female venom glands of ectoparasitoid Bracon hebetor (say.) (Hymenoptera: Braconidae). Comp Biochem Physiol Part D Genomics Proteomics. 2016;20:101–10.
Beard RL. Toxicology of Habrobracon venom - a study of a natural insecticide. New Haven: The Connecticut Agricultural Experiment Station; 1952.
Drenth D. Susceptibility of different species of insects to an extract of the venom gland of the wasp Microbracon hebetor (say). Toxicon Off J Int Soc Toxinology. 1974;12:189–92.
Piek T, Veenendaal RL, Mantel P. The pharmacology of Microbracon venom. Comp Biochem Physiol Part C Comp Pharmacol. 1982;72:303–9.
Walther C, Rathmayer W. The effect of Habrobracon venom on excitatory neuromuscular transmission in insects. J Comp Physiol. 1974;89:23–38.
Biondi A, Desneux N, Amiens-Desneux E, Siscaro G, Zappalà L. Biology and developmental strategies of the Palaearctic parasitoid Bracon nigricans (Hymenoptera: Braconidae) on the Neotropical moth Tuta absoluta (Lepidoptera: Gelechiidae). J Econ Entomol. 2013;106:1638–47.
Zappalà L, Bernardo U, Biondi A, Cocco A, Deliperi S, Delrio G, et al. Recruitment of native parasitoids by the exotic pest Tuta absoluta in southern Italy. Bull Insectol. 2012;65:51–61.
Becchimanzi A, Avolio M, Di Lelio I, Marinelli A, Varricchio P, Grimaldi A, et al. Host regulation by the ectophagous parasitoid wasp Bracon nigricans. J Insect Physiol. 2017;101:73–81.
Asgari S, Rivers DB. Venom proteins from endoparasitoid wasps and their role in host-parasite interactions. Annu Rev Entomol. 2011;56:313–35.
Kaas Q, Craik DJ. Bioinformatics-aided venomics. Toxins. 2015;7:2159–87.
Jonckheere W, Dermauw W, Zhurov V, Wybouw N, Bulcke JVd, Villarroel CA, et al. The salivary protein repertoire of the polyphagous spider mite Tetranychus urticae: a quest for effectors. Mol Cell Proteomics. 2016; mcp.M116.058081.
Villarroel CA, Jonckheere W, Alba JM, Glas JJ, Dermauw W, Haring MA, et al. Salivary proteins of spider mites suppress defenses in Nicotiana benthamiana and promote mite reproduction. Plant J. 2016;86:119–31.
Zhao C, Escalante LN, Chen H, Benatti TR, Qu J, Chellapilla S, et al. A massive expansion of effector genes underlies gall-formation in the wheat pest Mayetiola destructor. Curr Biol. 2015;25:613–20.
Aird SD, Aggarwal S, Villar-Briones A, Tin MM-Y, Terada K, Mikheyev AS. Snake venoms are integrated systems, but abundant venom proteins evolve more rapidly. BMC Genomics. 2015;16:647.
Martinson EO, Mrinalini, Kelkar YD, Chang C-H, Werren JH. The evolution of venom by co-option of single copy genes. Curr Biol. 2017;27:2007–2013.e8.
Cusumano A, Duvic B, Jouan V, Ravallec M, Legeai F, Peri E, et al. First extensive characterization of the venom gland from an egg parasitoid: structure, transcriptome and functional role. J Insect Physiol. 2018;107:68–80.
Moreau SJM, Guillot S. Advances and prospects on biosynthesis, structures and functions of venom proteins from parasitic wasps. Insect Biochem Mol Biol. 2005;35:1209–23.
Owen MD, Pfaff LA, Reisman RE, Wypych J. Phospholipase A2 in venom extracts from honey bees (Apis mellifera L.) of different ages. Toxicon. 1990;28:813–20.
Windass JD, Duncan RE, Baule VJ, Christian PD. Toxins from the wasp Bracon hebetor. 1996. https://patents.google.com/patent/WO1996016171A1/en. Accessed 28 Mar 2019.
Mingarro I, Pérez-Payá E, Pinilla C, Appel JR, Houghten RA, Blondelle SE. Activation of bee venom phospholipase A2 through a peptide-enzyme complex. FEBS Lett. 1995;372:131–4.
Guillaume C, Calzada C, Lagarde M, Schrével J, Deregnaucourt C. Interplay between lipoproteins and bee venom phospholipase A2 in relation to their anti-plasmodium toxicity. J Lipid Res. 2006;47:1493–506.
Cavalcante WLG, Noronha-Matos JB, Timóteo MA, Fontes MRM, Gallacci M, Correia-de-Sá P. Neuromuscular paralysis by the basic phospholipase A2 subunit of crotoxin from Crotalus durissus terrificus snake venom needs its acid chaperone to concurrently inhibit acetylcholine release and produce muscle blockage. Toxicol Appl Pharmacol. 2017;334:8–17.
Démares F, Coquerel Q, Richoux G, Linthicum K, Bloomquist J. Fatty acid and related potassium Kv2 channel blockers: toxicity and physiological actions on mosquitoes. Insects. 2018;9:155.
Posada Arias S, Rey-Suárez P, Pereáñez JA, Acosta C, Rojas M, Delazari Dos Santos L, et al. Isolation and functional characterization of an acidic myotoxic phospholipase A2 from colombian Bothrops asper venom. Toxins. 2017;9.
Nielsen VG. Carbon monoxide inhibits the anticoagulant activity of phospholipase A2 purified from Crotalus adamanteus venom. J Thromb Thrombolysis. 2018.
Perumal Samy R, Gopalakrishnakone P, Thwin MM, Chow TKV, Bow H, Yap EH, et al. Antibacterial activity of snake, scorpion and bee venoms: a comparison with purified venom phospholipase A2 enzymes. J Appl Microbiol. 2007;102:650–9.
Annand RR, Kontoyianni M, Penzotti JE, Dudler T, Lybrand TP, Gelb MH. Active site of bee venom phospholipase A2: the role of histidine-34, aspartate-64 and tyrosine-87. Biochemistry. 1996;35:4591–601.
Kini RM. Excitement ahead: structure, function and mechanism of snake venom phospholipase A2 enzymes. Toxicon Off J Int Soc Toxinology. 2003;42:827–40.
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.
Dani MP, Edwards JP, Richards EH. Hydrolase activity in the venom of the pupal endoparasitic wasp, Pimpla hypochondriaca. Comp Biochem Physiol B Biochem Mol Biol. 2005;141:373–81.
Vincent B, Kaeslin M, Roth T, Heller M, Poulain J, Cousserans F, et al. The venom composition of the parasitic wasp Chelonus inanitus resolved by combined expressed sequence tags analysis and proteomic approach. BMC Genomics. 2010;11:693.
Crawford AM, Brauning R, Smolenski G, Ferguson C, Barton D, Wheeler TT, et al. The constituents of Microctonus sp. parasitoid venoms. Insect Mol Biol. 2008;17:313–24.
Gopalapillai R, Kadono-Okuda K, Okuda T. Molecular cloning and analysis of a novel teratocyte-specific carboxylesterase from the parasitic wasp, Dinocampus coccinellae. Insect Biochem Mol Biol. 2005;35:1171–80.
Nakamatsu Y, Tanaka T. Venom of Euplectrus separatae causes hyperlipidemia by lysis of host fat body cells. J Insect Physiol. 2004;50:267–75.
Dorémus T, Urbach S, Jouan V, Cousserans F, Ravallec M, Demettre E, et al. Venom gland extract is not required for successful parasitism in the polydnavirus-associated endoparasitoid Hyposoter didymator (Hym. Ichneumonidae) despite the presence of numerous novel and conserved venom proteins. Insect Biochem Mol Biol. 2013;43:292–307.
Visser B, Lann CL, den Blanken FJ, Harvey JA, van Alphen JJM, Ellers J. Loss of lipid synthesis as an evolutionary consequence of a parasitic lifestyle. Proc Natl Acad Sci. 2010;107:8677–82.
Sim AD, Wheeler D. The venom gland transcriptome of the parasitoid wasp Nasonia vitripennis highlights the importance of novel genes in venom function. BMC Genomics. 2016;17:571.
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.
Asgari S, Zhang G, Zareie R, Schmidt O. A serine proteinase homolog venom protein from an endoparasitoid wasp inhibits melanization of the host hemolymph. Insect Biochem Mol Biol. 2003;33:1017–24.
Cerenius L, Kawabata S, Lee BL, Nonaka M, Söderhäll K. Proteolytic cascades and their involvement in invertebrate immunity. Trends Biochem Sci. 2010;35:575–83.
Colinet D, Anselme C, Deleury E, Mancini D, Poulain J, Azéma-Dossat C, et al. Identification of the main venom protein components of Aphidius ervi, a parasitoid wasp of the aphid model Acyrthosiphon pisum. BMC Genomics. 2014;15:342.
Parkinson N, Richards EH, Conyers C, Smith I, Edwards JP. Analysis of venom constituents from the parasitoid wasp Pimpla hypochondriaca and cloning of a cDNA encoding a venom protein. Insect Biochem Mol Biol. 2002;32:729–35.
Taylor A. Aminopeptidases: structure and function. FASEB J Off Publ Fed Am Soc Exp Biol. 1993;7:290–8.
Dani MP, Richards EH, Isaac RE, Edwards JP. Antibacterial and proteolytic activity in venom from the endoparasitic wasp Pimpla hypochondriaca (Hymenoptera: Ichneumonidae). J Insect Physiol. 2003;49:945–54.
Ragheb D, Dalal S, Bompiani KM, Ray WK, Klemba M. Distribution and biochemical properties of an M1-family aminopeptidase in Plasmodium falciparum indicate a role in vacuolar hemoglobin catabolism. J Biol Chem. 2011;286:27255–65.
Santos PP, Games PD, Azevedo DO, Barros E, de Oliveira LL, Ramos HJ d O, et al. Proteomic analysis of the venom of the predatory ant Pachycondyla striata (Hymenoptera: Formicidae). Arch Insect Biochem Physiol. 2017;96:e21424.
Vaiyapuri S, Wagstaff SC, Watson KA, Harrison RA, Gibbins JM, Hutchinson EG. Purification and functional characterisation of Rhiminopeptidase a, a novel aminopeptidase from the venom of Bitis gabonica rhinoceros. PLoS Negl Trop Dis. 2010;4:e796.
Paciotti S, Codini M, Tasegian A, Ceccarini MR, Cataldi S, Arcuri C, et al. Lysosomal alpha-mannosidase and alpha-mannosidosis. Front Biosci Landmark Ed. 2017;22:157–67.
Fuzita FJ, Pinkse MWH, Patane JSL, Juliano MA, Verhaert PDEM, Lopes AR. Biochemical, transcriptomic and proteomic analyses of digestion in the scorpion Tityus serrulatus: insights into function and evolution of digestion in an ancient arthropod. PLoS One. 2015;10. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0123841.
Fuzita FJ, Pinkse MWH, Patane JSL, Verhaert PDEM, Lopes AR. High throughput techniques to reveal the molecular physiology and evolution of digestion in spiders. BMC Genomics. 2016;17:716.
Walter A, Bechsgaard J, Scavenius C, Dyrlund TS, Sanggaard KW, Enghild JJ, et al. Characterisation of protein families in spider digestive fluids and their role in extra-oral digestion. BMC Genomics. 2017;18:600.
Cohen AC. Extra-oral digestion in predaceous terrestrial arthropoda. Annu Rev Entomol. 1995;40:85–103.
Cohen AC. Solid-to-liquid feeding: the inside(s) story of extra-oral digestion in predaceous arthropoda. Am Entomol. 1998;44:103–17.
Jarjees E, Merritt DJ, Gordh G. Anatomy of the mouthparts and digestive tract during feeding in larvae of the parasitoid wasp Trichogramma australicum Girault (Hymenoptera : Trichogrammatidae). Int J Insect Morphol Embryol. 1998;27:103–10.
Falabella P, Tremblay E, Pennacchio F. Host regulation by the aphid parasitoid Aphidius ervi: the role of teratocytes. Entomol Exp Appl. 2000;97:1–9.
Lavine MD, Strand MR. Insect hemocytes and their role in immunity. Insect Biochem Mol Biol. 2002;32:1295–309.
Pelosi P, Maida R. Odorant-binding proteins in insects. Comp Biochem Physiol B Biochem Mol Biol. 1995;111:503–14.
Ahmed T, Zhang T, Wang Z, He K, Bai S. Molecular cloning, expression profile, odorant affinity, and stability of two odorant-binding proteins in Macrocentrus cingulum Brischke (Hymenoptera: Braconidae). Arch Insect Biochem Physiol. 2017;94:e21374.
Leal WS, Nikonova L, Peng G. Disulfide structure of the pheromone binding protein from the silkworm moth, Bombyx mori. FEBS Lett. 1999;464:85–90.
Yang K, Liu Y, Niu D-J, Wei D, Li F, Wang G-R, et al. Identification of novel odorant binding protein genes and functional characterization of OBP8 in Chilo suppressalis (Walker). Gene. 2016;591:425–32.
Kucharski R, Maleszka J, Maleszka R. A possible role of DNA methylation in functional divergence of a fast evolving duplicate gene encoding odorant binding protein 11 in the honeybee. Proc R Soc B Biol Sci. 2016;283:20160558.
Heavner ME, Gueguen G, Rajwani R, Pagan PE, Small C, Govind S. Partial venom gland transcriptome of a Drosophila parasitoid wasp, Leptopilina heterotoma, reveals novel and shared bioactive profiles with stinging Hymenoptera. Gene. 2013;526:195–204.
Wang L, Zhu J-Y, Qian C, Fang Q, Ye G-Y. Venom of the parasitoid wasp Pteromalus puparum contains an odorant binding protein. Arch Insect Biochem Physiol. 2015;88:101–10.
Wang T, Zhao M, Rotgans BA, Ni G, Dean JFD, Nahrung HF, et al. Proteomic analysis of the venom and venom sac of the woodwasp, Sirex noctilio - towards understanding its biological impact. J Proteome. 2016;146:195–206.
Ishida Y, Ishibashi J, Leal WS. Fatty acid solubilizer from the oral disk of the blowfly. PLoS One. 2013;8:e51779.
Pelosi P, Iovinella I, Zhu J, Wang G, Dani FR. Beyond chemoreception: diverse tasks of soluble olfactory proteins in insects. Biol Rev. 2018;93:184–200.
Noiva R. Protein disulfide isomerase: the multifunctional redox chaperone of the endoplasmic reticulum. Semin Cell Dev Biol. 1999;10:481–93.
Wilkinson B, Gilbert HF. Protein disulfide isomerase. Biochim Biophys Acta BBA - Proteins Proteomics. 1699;2004:35–44.
Safavi-Hemami H, Li Q, Jackson RL, Song AS, Boomsma W, Bandyopadhyay PK, et al. Rapid expansion of the protein disulfide isomerase gene family facilitates the folding of venom peptides. Proc Natl Acad Sci U S A. 2016;113:3227–32.
Violette A, Leonardi A, Piquemal D, Terrat Y, Biass D, Dutertre S, et al. Recruitment of glycosyl hydrolase proteins in a cone snail venomous arsenal: further insights into biomolecular features of Conus venoms. Mar Drugs. 2012;10:258.
Turano C, Coppari S, Altieri F, Ferraro A. Proteins of the PDI family: unpredicted non-ER locations and functions. J Cell Physiol. 2002;193:154–63.
Munro S, Pelham HR. A C-terminal signal prevents secretion of luminal ER proteins. Cell. 1987;48:899–907.
Pennacchio F, Giordana B, Rao R. Chapter 22 - applications of parasitoid virus and venom research in agriculture. In: Beckage NE, Drezen J-M, editors. Parasitoid Viruses. San Diego: Academic Press; 2012. p. 269–83.
Di Lelio I, Illiano A, Astarita F, Gianfranceschi L, Horner D, Varricchio P, et al. Evolution of an insect immune barrier through horizontal gene transfer mediated by a parasitic wasp. Plos Genet Press. 2019.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal. 2011;17:10–2.
Andrews S. FastQC: a quality control tool for high throughput sequence data; 2010.
Krueger F Trim Galore!: A wrapper tool around Cutadapt and FastQC to consistently apply quality and adapter trimming to FastQ files. 0.4; 2015.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat Biotechnol. 2011;29:644–52.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.
Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinforma Oxf Engl. 2012;28:3150–2.
UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2017;45(Database issue):D158–69.
Chen T-W, Gan R-CR WTH, Huang P-J, Lee C-Y, Chen Y-YM, et al. FastAnnotator- an efficient transcript annotation web tool. BMC Genomics. 2012;13(Suppl 7):S9.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785.
Jungo F, Bougueleret L, Xenarios I, Poux S. The UniProtKB/Swiss-Prot Tox-Prot program: a central hub of integrated venom protein data. Toxicon. 2012;60:551–7.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Pringle JWS. Proprioception in insects: II. The action of the campaniform sensilla on the legs. J Exp Biol. 1938;15:114–31.
Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72:248–54.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods San Diego Calif. 2001;25:402–8.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ. Jalview version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25:1189–91.
de Castro E, Sigrist CJA, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, et al. ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006;34(suppl_2):W362–5.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.
Thanks are due to Ilaria Di Lelio, University of Napoli Federico II (Italy), for her valuable help in the rearing and maintenance of B. nigricans and S. littoralis, and to Lucia Zappalà, University of Catania (Italy), for providing the insect material used to start the colony of B. nigricans reared in our laboratory.
This work was supported by the European Union’s Horizon 2020 Research and Innovation Programme, under grant agreement no. 773554 (EcoStack). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
Not applicable, no human participants, human data or human tissue. Research was conducted in compliance with institutional, national and international ethical guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Distribution of transcripts abundance expressed in RPKM
Table S1. List of the most abundant (RPKM>100) annotated transcripts encoding putatively secreted proteins (i.e., positive to SignalP analysis)
Table S2. Venom proteome annotated by BLAST search
Table S3. Homologs of B. nigricans venom proteins occurring in other parasitoid species
Table S4. One-way ANOVA of ΔCt values recorded in venom glands, females devoid of venom glands and males
Figures S2-S9. Multiple sequence alignments of most representative B. nigricans venom proteins
Figures S10-S17. Phylogenetic trees of most representative B. nigricans venom proteins
Table S5. Primers used for qRT-PCR analysis of selected venom components
About this article
Cite this article
Becchimanzi, A., Avolio, M., Bostan, H. et al. Venomics of the ectoparasitoid wasp Bracon nigricans. BMC Genomics 21, 34 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-019-6396-4
- Idiobiont parasitoids
- Host-parasitoid interactions
- Venom proteome
- Venom transcriptome
- Host regulation