Next Article in Journal
Role of Shiga Toxins in Cytotoxicity and Immunomodulatory Effects of Escherichia coli O157:H7 during Host-Bacterial Interactions in vitro
Next Article in Special Issue
Diet Breadth Mediates the Prey Specificity of Venom Potency in Snakes
Previous Article in Journal
Cellular Effects of T-2 Toxin on Primary Hepatic Cell Culture Models of Chickens
Previous Article in Special Issue
Coevolution of Snake Venom Toxic Activities and Diet: Evidence that Ecological Generalism Favours Toxicological Diversity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 12
Issue 1
10.3390/toxins12010047
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Venom Components and Their Phylogenetic Properties in Some Aculeate Bumblebees and Wasps

by
Kyungjae Andrew Yoon
1,*,
Kyungmun Kim
2,
Woo-Jin Kim
3,
Woo Young Bang
4,
Neung-Ho Ahn
4,
Chang-Hwan Bae
4,
Joo-Hong Yeo
4 and
Si Hyeock Lee
1,2,*
1
Research Institute of Agriculture and Life Sciences, Seoul National University, Seoul 08826, Korea
2
Department of Agricultural Biology, Seoul National University, Seoul 08826, Korea
3
EntoCode Co., Seoul 06028, Korea
4
National Institute of Biological Resources, Environmental Research Complex, Incheon 22689, Korea
*
Authors to whom correspondence should be addressed.
Toxins 2020, 12(1), 47; https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010047
Submission received: 14 December 2019 / Revised: 9 January 2020 / Accepted: 11 January 2020 / Published: 14 January 2020
(This article belongs to the Special Issue Evolutionary Ecology of Venom)
Browse Figure
Versions Notes

Abstract

:
To identify and compare venom components and expression patterns, venom gland-specific transcriptome analyses were conducted for 14 Aculeate bees and wasps. TPM (transcripts per kilobase million) values were normalized using the average transcription level of a reference housekeeping gene (dimethyladenosine transferase). Orthologous venom component genes across the 14 bee and wasp species were identified, and their relative abundance in each species was determined by comparing normalized TPM values. Based on signal sequences in the transcripts, the genes of novel venom components were identified and characterized to encode potential allergens. Most of the allergens and pain-producing factors (arginine kinase, hyaluronidase, mastoparan, phospholipase A1, phospholipase A2, and venom allergen 5) showed extremely high expression levels in social wasps. Acid phosphatase, neprilysin, and tachykinin, which are known allergens and neurotoxic peptides, were found in the venom glands of solitary wasps more often than in social wasps. In the venom glands of bumblebees, few or no transcripts of major allergens or pain-producing factors were identified. Taken together, these results indicate that differential expression patterns of the venom genes in some Aculeate species imply that some wasps and bumblebee species have unique groups of highly expressed venom components. Some venom components reflected the Aculeate species phylogeny, but others did not. This unique evolution of specific venom components in different groups of some wasps and bumblebee species might have been shaped in response to both ecological and behavioral influences.
Keywords:
venom; social wasp; solitary hunting wasp; bumblebee; venom gland; transcriptome analysis
Key Contribution: Venom components and expression patterns were predicted by parallel analyses of venom gland-specific transcriptomes in 14 Aculeate bumblebee and wasp species. Differential expression patterns of some major venom genes in these species implied that venom components might have been evolved in response to both ecological and behavioral influences.

1. Introduction

Wasps are an extremely diverse group of hymenopteran insects that includes an estimated 70,000 species worldwide. Wasps belong to the suborder Apocrita (Hymenoptera), which is divided into Parasitica and Aculeata [1]. Parasitica comprises most of the parasitoid wasps, whereas Aculeata contains parasitoids and predatory wasps, in which the ovipositor is modified to transmit venom rather than being an egg-positioning device [2]. Aculeate wasps are divided into solitary and social groups [3].
Most of Aculeate wasps, approximately 95%, are solitary wasps, which build their own nests and store their prey [4]. Most solitary wasps paralyze and preserve prey to provision their young using venom composed of bioactive molecules with various functions, including paralysis, antimicrobial activity, and developmental arrest [5].
All social wasps belong to the family Vespidae, which occupies only a small portion of Aculeata. Social wasps form colonies of various sizes, and their sting is generally reserved for defensive use against intruders [6,7,8]. Some social vespid wasps have guards that are more likely than others in the colony to sting intruders, and old individuals of the species Polybia occidentalis are reported to defend their colonies more actively than young individuals [9]. Once attacked or disturbed by intruders, including humans, an attack pheromone can cause fatal mass envenomation by attracting all wasps in the colony [10]. Social wasps capture arthropod prey without stinging and masticate the prey for their brood meal [11]. Because social wasps do not paralyze or preserve prey, their venom appears to have evolved as a defensive tool [5]. Because the venoms of social wasps contain various allergens that can cause an anaphylaxis reaction in humans, their medical and clinical aspects have been studied extensively. Social wasp venoms cause pain, local edema, erythema, cardiovascular disturbances, and anaphylactic shock by increasing the permeability of blood vessels [12]. Many patients exposed to wasp venoms develop systemic reactions and immediate or late reactions such as acute renal failure, hemolysis and coagulopathy, and severe neurological disturbances such as cerebral infarction [13,14,15]. Immunological reactions, such as type I hypersensitivity caused by venom allergens, have also been observed [16]. Those symptoms are caused by biologically active peptides in venom, such as mastoparan, chemotactic peptides and kinins [5,17,18]. Venom allergens include arginine kinase, phospholipase A (A1 and A2), hyaluronidase, cysteine-rich secretory proteins, antigen 5, pathogenesis-related proteins (e.g., CAP), and serine proteases [5].
Bumblebees, which belong to the family Apidae, are increasingly used for pollination of greenhouse plants because of their naive characteristics and lifestyle fit for pollinating. They are not aggressive and do not sting unless disturbed near their nests or intentionally touched [19]. Although the frequency of stinging incidents is expected to increase with increased use of bumblebees for pollination in greenhouses, little information has been available on the comparative toxicity of bumblebee venom.
Venom components and their abundance in bee and wasp venom have mostly been studied in individual species [5,20]. Only a few studies have compared the components in venoms of two closely-related wasp species [21,22]. A comparative analysis of venom gland transcriptomes of two closely-related social hornet species (Vespa crabro and V. analis) revealed similar overall venom components, but most of the major venom component genes were more abundantly transcribed in V. crabro [21]. However, little attempt has been made for the systematic analysis and comparison of the venom components in multiple bee and wasp species. The recent introduction of cost-effective, high-throughput deep-sequencing technologies has enabled rapid identification of genes encoding various venom peptides and proteins from venom gland transcriptomes.
Although the transcriptome profiles of venom gland may not reflect the actual composition profiles of venom proteins and peptides, venom gland transcriptomic approaches have shed some light on the composition and evolution of venom in various venomous animals (reviewed by [23]). The transcriptomic approach is particularly useful when target animals have limited sizes as in case of bees and wasps, and thus venom proteomics are not feasible (reviewed by [23]). In addition, high correlations were also reported between venom gland transcriptomes and venom proteomes in 10 of our 11 comparisons of snake venoms [24].
Therefore, we conducted a comparative transcriptomic characterization of the venom glands of 14 Aculeate bee and wasp species. The TPM (transcripts per kilobase million) values of representative venom components were normalized using the average transcription level of a reference housekeeping gene. The orthologous venom component genes were identified across the 14 bee and wasp species, and their relative abundance was compared across species. Based on signal sequences, genes for putative novel venom components were identified and characterized. Finally, we investigated the evolutionary aspects of the venom components in accordance with ecological and behavioral features. This systematic and comparative study of the venom gland transcriptomes of different bee and wasp species provides new insights into the evolution and phylogeny of venom.

2. Results and Discussion

2.1. Differential Transcriptional Profiles of the Venom Glands of Some Aculeate Species

Total reads of 5.89, 6.39, 7.13, 5.74, 6.22, 6.16, 11.3, 11.4, 8.2, 8.07, 8.1, 7.66, 10.9, and 10.6 Gb were obtained from RNA sequencing of the venom glands of Eumenes decoratus, Sphecidae sp., Anterhynchium flavomarginatum, Sceliphron deforme, V. crabro, V. analis, V. dybowskii, V. simillima, Parapolybia varia, Polistes snelleni, Polistes rothneyi, Bombus ardens, B. consobrinus, and B. ussurensis, respectively. De novo assembly of the trimmed sequence data resulted in 18,062, 31,134, 21,512, 24,483, 11,097, 12,531, 21,227, 29,150, 20,378, 27,742, 25,436, 31,366, 21,625, and 26,425 assembled contigs, respectively. A total of 16,357, 27,423, 19,033, 20,097, 10,321, 11,284, 19,627, 27,065, 18,539, 25,879, 23,475, 26,158, 20,050, and 22,619 genes, respectively, were identified by BLAST search (Table 1).
To verify the transcriptional abundance of major venomic genes in Aculeate bees and wasps, 10 venomic gens of two social wasps (P. rothneyi and P. snelleni), one solitary hunting wasp (E. decoratus) and one bumblebee (B. ardens) were investigated by conducting quantitative real-time PCR (qPCR) and compared with TPM value of their transcriptome data (Figure S1). Relative transcription levels showed significantly similar expression pattern with that of TPM values indicating that relative transcriptional abundances of major venomic genes estimated from TPM values are reliable.

2.1.1. Venom Proteins Highly Expressed in Solitary Hunting Wasps

Some venom component genes, such as acid phosphatase, neprilysin, and tachykinin, were abundantly expressed in solitary hunting wasps and expression levels of allergens and pain-producing factors were low in social wasps and bumblebees. Acid phosphatase was found with a high frequency in the venom glands of solitary wasps (Figure 1). In particular, the solitary hunting wasp S. deforme exhibited an extremely high expression level of acid phosphatase (TPM = 83,708). The overall TPM values for acid phosphatase in other wasps and bumblebee species were less than 100 (except that of V. analis), suggesting that acid phosphatase might have a particularly significant role in the venom of solitary hunting wasps (Figure 1). Acid phosphatase in the venom of the endoparasitoid wasp Pteromalus puparum was suggested to affect the host’s immune response and physiology [25]. However, that of the endoparasitoid wasp Pimpla hypochondriaca plays no known role in hemocyte immunology [26].
Neprilysin showed an extremely high expression level (TPM = 7387) in E. decoratus and low expression levels (TPM < 50) in two solitary hunting wasps, Sphecidae sp. and A. flavomarginatum (Figure 1). The social wasps showed none of the transcript for neprilysin, indicating that neprilysin might be specifically expressed in the venom glands of solitary hunting wasps. Neprilysin is involved in metabolism of regulatory peptides in mammalian nervous systems, and it is known to inactivate peptide transmitters and their modulators by terminating brain neuropeptides at peptidergic synapses [27,28]. This implies that paralysis might be one of the major functions of the venom of solitary hunting wasps.
Tachykinin was highly expressed in the solitary hunting wasp E. decoratus (TPM = 4241), and the other wasps and bumblebee species tested showed TPM values lower than 40 (Figure 1). Tachykinin is a neurotoxic peptide in wasp venom that is known to induce hypokinesia, a “zombie-like” behavioral state caused by direct envenomation into the central nervous system. Synthetic AcVTK mimics the venom tachykinins of the parasitoid jewel wasp Ampulex compressa. Injection of synthetic AcVTK into the subesophageal ganglion of the American cockroach (the jewel wasp host) produced behavioral changes, including suppression of the escape response and reduced spontaneous walking [29]. Thus, solitary wasps might use tachykinin to change a prey’s physiology and behavior, ensuring availability of fresh prey for their larvae.
Defensin 2 was abundantly identified in the venom glands of solitary hunting wasps (E. decoratus, Sphecidae sp., and A. flavomarginatum), with only slight expression in some social wasps (Parapolybia varia and Polistes rothneyi). This finding implies that defensin 2 in solitary hunting wasps might be used to preserve prey.

2.1.2. Venom Proteins Highly Expressed in Social Wasps

According to the hierarchical clustering analysis for the 22 representative venom genes, most of the allergens and pain-producing factors (arginine kinase, dipeptidyl peptidase 4, endocuticle structural glycoprotein, hyaluronidase, mastoparan, phospholipase A1, phospholipase A2, serine protease, serine protease inhibitor, venom allergen 5, and vitellogenin) were clustered in a clade although its bootstrap value (52%) did not strongly support the topology (Figure 1; see red rectangle). These venom genes showed significantly high expression levels in the social wasps examined. This suggests that social wasps might have evolved to use these components as defensive tools against intruding species (Figure 1).
Arginine kinase, which is known to cause an allergic reaction in vertebrates [30], was identified in the venom glands of most Aculeate wasps and bumblebee species, suggesting its ancient origins. The arginine kinase of P. varia exhibited extremely high TPM value among the 14 species tested, and most of the social wasps showed values above 300 (except V. crabro, V. simillima and P. rothneyi). Solitary hunting wasps and bumblebees (except B. ardens and B. ussurensis) also showed considerable levels of arginine kinase expression, and expression levels were low in social wasps (i.e., TPM values < 100). B. consobrinus showed a high expression level of arginine kinase (TPM = 832) (Figure 1). Presence of arginine kinase in both wasps and bumblebees suggests that it has likely been used as an allergen in the lineages of both bees and wasps.
Most solitary and social wasps showed high expression levels of dipeptidyl peptidase 4, with particularly high expression in S. deforme (TPM = 54,119) (Figure 1). Dipeptidyl peptidase 4 is reported to be a major allergen in the social wasp Polistes dominula, with more than 60% of Polistes venom sensitized—patients showing IgE-reactivity [31]. Its relatively low expression levels in bumblebees suggest that dipeptidyl peptidase 4 may function as a wasp-specific allergen.
Endocuticle structural glycoprotein, which is known to cause allergic reactions, was highly expressed in most social wasps except V. crabro (Figure 1). The carbohydrate moiety of insect glycoproteins is reported to play a pivotal role in human IgE binding [32], implying that endocuticle structural glycoprotein may contribute to allergic reactions to social wasp venoms.
Most social wasps, V. crabro and V. dybowskii in particular, showed high expression levels of hyaluronidase (TPM values > 1000), although some of the solitary hunting wasps (S. deforme and Sphecidae sp.) also exhibited a moderate level of expression (TPM = 893 and 556) (Figure 1). This result indicates that hyaluronidase, which is a venom spreading factor, might play a major role in the venom of both social and solitary wasps [5].
Mastoparan was identified in the venom of most social wasps (V. crabro, V. analis, P. snelleni, and P. rothneyi), some solitary hunting wasps (A. flavomarginatum and S. deforme), and a bumblebee species (B. ussurensis). V. crabro mastoparan showed an extremely high TPM value among all venom genes in that species. Mastoparan causes allergic inflammation and mast cell degranulation in vertebrates [21]. Therefore, social wasp venoms containing relatively high levels of mastoparan are likely to be highly allergenic to vertebrates (Figure 1). Identification of mastoparan in both wasp and bumblebee species suggests its ancient origin.
Phospholipase A1 was identified in the venom glands of all social wasps and some of the solitary hunting wasps and bumblebees (Figure 1). Most social wasps exhibited extremely high expression of phospholipase A1 (with the exception of P. snelleni and P. rothneyi). Phospholipase A2 was also highly expressed in most social wasps, exhibiting TPM value > 600, especially in V. crabro, V. dybowskii, and V. simillima (Figure 1). As phospholipase A1 and A2 are major Hymenopteran venom allergens in humans [33,34], their high expression in the venoms social wasps suggests that they are likely to be involved in the defensive function of venom.
Abundant transcripts of serine protease were found in the social wasps, among which V. dybowskii showed the highest expression (TPM = 1916). Few or no transcripts of serine protease were expressed in solitary wasps and bumblebees (Figure 1). Serine protease is an important allergen with significant IgE binding activity [35,36,37], which indicates that it might have defensive role to social wasps.
Most of the Aculeate wasps and bumblebee species tested showed abundant transcripts of serine protease inhibitor (Serpin), except V. crabro, V. simillima, E. decorates, and S. deforme. Social wasps such as V. analis, P. varia, and P. rothneyi showed moderate expression of serpin (TPM = 736, 660, and 604, respectively) (Figure 1). Serpin is known to have IgE-sensitizing potential in honeybee venom-allergic patients, indicating its potency as an allergen [38].
The social wasps V. crabro and V. simillima exhibited extremely high expression of venom allergen 5, which is a major allergen identified in many insect venoms [39]. This finding suggests that the venom of social wasps contains abundant allergens (Figure 1).
Vitellogenin is a newly found honeybee venom allergen that shows IgE-sensitizing reactivity in venom-sensitized patients [40]. V. analis, P. varia, P. rothneyi, and A. flavomarginatum exhibited extremely high expression of vitellogenin, and V. dybowskii, and B. ardens showed considerable expression (Figure 1). Thus, social wasp venoms with high expression of vitellogenin might possess high allergenicity.

2.1.3. Venom Proteins Highly Expressed in Bumblebees

In the venom glands of bumblebees, few transcripts for major allergens or pain-producing factors were found, implying that bumblebee venoms are likely to have low toxicity (Figure 1). Interestingly, defensin 2 was not detected in the venom glands of the three bumblebees tested, but defensin 1 (an antimicrobial peptide) showed extremely high expression in B. ardens and B. ussurensis.
Bumblebees (except B. consobrinus) had high expression level of defensin 1. The social wasps V. analis, V. dybowskii, and P. varia showed TPM values > 200; other wasps (except Sphecidae sp.) exhibited little to no defensin 1 transcript expression (Figure 1). The high bumblebee-specific expression of defensin 1, which is well known as an antimicrobial peptide, suggests that bumblebees may have exploited defensin to sanitize their nests, which are built underground. The small number of transcripts for defensin 1 in B. consobrinus suggest that this species might be more susceptible to infection than the other bumblebees examined in this study.

2.1.4. Expression Patterns of Other Venom Proteins

Carboxylesterase 6 causes allergic reactions in humans and is also known as Api m 8, which is a honeybee venom allergen [41]. The social wasp P. varia showed the highest expression of carboxylesterase 6, and the bumblebees B. consobrinus and B. ussurensis exhibited TPM values higher than 200 (Figure 1). Solitary wasps had few or no carboxylesterase 6 transcripts, suggesting that this protein might have a defensive role in social wasp and bumblebee venoms.
Bumblebees, especially B. consobrinus, showed a high expression level of icarapin (TPM = 8901), but all solitary (except for E. decoratus) and social wasps also exhibited high TPM value (Figure 1). Icarapin is a glycoprotein allergen with low abundance in honeybee venom [42] that nonetheless shows greater than 50% IgE-reactivity in honeybee venom-sensitized patients [43]. The extremely high expression of icarapin in bumblebees and its considerable expression in most social and solitary hunting wasps indicate that some Aculeate wasps and bumblebee species are likely to commonly use icarapin for defense by causing allergic reactions.
Metalloproteinase generally showed low expression in the Aculeate wasps and bumblebee species tested, except the solitary hunting wasp Sphecidae sp., which exhibited a considerable number of transcripts (TPM = 516) (Figure 1). Metalloproteinase is a toxin found abundantly in snake venom and causes hemorrhaging and interferes with the hemostatic system [44]. Metalloproteinase is also known to play a role in tumor invasion and metastasis by degrading extracellular matrix [45]. The generally low abundance of metalloproteinase in the venoms of Aculeate species may indicate its minor role for wasps and bumblebee species.
Phospholipase B showed a low abundance in most Aculeate wasps and bumblebee species examined (Figure 1). Although phospholipase B was reported to have cytotoxicity in snake venom [46], its low expression may imply its minor role as a venom component of Aculeate species, especially compared with phospholipase A1 and A2.
Serine carboxypeptidase is an exopeptidase that is a major allergen in honeybee venom [47]. It is known to participate in protein digestion in the guts of animals [48] and is found among the yolk proteins in mosquito oocytes [49]. Serine carboxypeptidase showed low abundance in solitary wasps and bumblebees (except B. consobrinus, TPM = 161) (Figure 1). The presence of few or no transcripts in social wasps indicates that serine carboxypeptidase might have a minor function in the venom of some Aculeate wasps and bumblebee species.

2.2. Identification of Putatively Novel Venom Components in Some Aculeate Bees and Wasps

To identify the properties of novel venom genes, we conducted signal peptide prediction with the top 100 highly expressed genes (Tables S1–S14) and summarized the resulting secretory venom components in Table 2. A total of 15 genes coding secretory proteins were identified in all 14 Aculeate bee/wasp species. Among them, nine genes were found to encode known venom proteins: defensin 1, endocuticle structural glycoprotein, hyaluronidase, icarapin, neprilysin, phospholipase A1 and A2, serine protease, and venom allergen 5 (Table 2). One of the reasons that not all known venom component genes were identified by signal peptide prediction is that many of the venom component genes lack the 5’ end transcript sequences where the signal peptide sequences are located, probably due to incomplete reverse transcription of the 5’end mRNA by reverse transcriptase [50] or incomplete reconstruction of the full-length transcripts from short cDNA reads [51].
Nevertheless, signal peptide prediction enabled identification of six genes encoding uncharacterized secretory proteins (USPs; USP1–USP6) that were commonly expressed in all examined social wasps, bumblebees, and one solitary hunting wasp, A. flavomarginatum (Table 2). Such common and high expression pattern in most Aculeate wasps and bumblebee species examined suggests that the components are likely putative constituents of venom.

Properties of the Six USPs (USP1–USP6) as Putative Allergenic Proteins

USP1 exhibited high expression in most social wasps, especially in V. simillima (TPM = 29,700) (Table 2). The USP1 gene contained 106 to 112 amino acid sequences and was identified as a putative pilosulin from the Australian ant Myrmecia banksi (56.3%, accession number: BAD36780), implying its allergenic function in social wasp venom (Figure S2).
USP2 exhibited high expression in two social wasps, V. analis and P. varia (TPM > 2600), and a considerable expression level in the solitary hunting wasp A. flavomarginatum (TPM = 865). USP2 P. varia showed 70.3% sequence similarity with putative tropomyosin from the Pacific oyster Crassostrea gigas (BAH10152), suggesting USP2 as a putative allergen (Figure S3).
USP3 was highly expressed in two bumblebees, B. consobrinus and B. ussurensis (TPM = 1805 and 1537, respectively). USP3 encodes 167 to 173 amino acid sequences that were identified as related to the putative IgE-binding protein MnSOD of the rubber tree Hevea brasiliensis (CAC13961), with 59.6% sequence similarity (Figure S4).
USP4 was identified in the venom gland of the bumblebees B. ardens and B. ussurensis and had high expression (TPM = 1800 and 1741, respectively). In these two bumblebees, USP4 encoded the same 140-amino acid sequence and exhibited 58.8% sequence similarity with the allergen alpha amylase inhibitor-like protein from common wheat, Triticum aestivum (CAA35597) (Figure S5). Thus, USP3 and USP4 in bumblebees might have the potential to induce allergic sensitization.
USP5 was identified in the social wasp P. varia and the solitary hunting wasp A. flavomarginatum with high expression (TPM = 2439 and 1322, respectively). Only USP5 in A. flavomarginatum, which encodes a 143-amino acid sequence, was identified as the putative venom allergen Sol g 4.01 precursor from the tropical fire ant Solenopsis geminata (AF230383), with 54.6% sequence similarity (Figure S6).
USP6 showed high expression in the social wasp V. dybowskii and the solitary hunting wasp A. flavomarginatum (TPM = 1584 and 1944, respectively). USP6, composed of 129 to 130-amino acid sequences in two species, could not be identified in a BLAST search or an allergen database search (Figure S7).

2.3. Evolutionary Patterns in the Venom of Some Aculeate Bees and Wasps

Aculeata belongs to Apocrita, Hymenoptera and is divided into the three superfamilies Chrysidoidea, Apoidea, and Vespoidea [5]. Bumblebees belong to the Apidae family, and Sphecidae sp. And S. deforme belong to the Sphecidae family, which is divided from the superfamily Apoidea. The Vespa species (V. analis, V. crabro, V. dybowskii, and V. simillima) belong to the subfamily Vespinae, and the Polistes species (P. rothneyi and P. snelleni) and P. varia belong to the subfamily Polistinae, which is grouped with Vespinae. E. decoratus and A. flavomarginatum belong to the subfamily Eumeninae [5].
In the phylogenetic tree of venom components in some Aculeate bees and wasps, eight venom proteins (carboxylesterase 6, dipeptidyl peptidase 4, endocuticle structural glycoprotein, icarapin, major royal jelly protein, phospholipase B, serine protease inhibitor and vitellogenin) matched the evolutionary pattern of Aculeate species (Figures S8–S15). In phylogenetic analysis, icarapin of social wasps and that of bumblebees were clearly separated into two branches. Sphecidae sp. and S. deforme were grouped and closely related to the bumblebee branch, whereas E. decoratus and A. flavomarginatum were grouped and showed a close relationship with the social wasp branch (Figure S11). Thus, icarapin appears to reflect the phylogenetic relationships among some Aculeate bees and wasps.
However, the number of venom component proteins that do not phylogenetically reflect the Aculeate lineage was larger than the number that do reflect. The results from the phylogenetic analyses of acid phosphatase, defensin 2, hyaluronidase, mastoparan, metalloproteinase, neprilysin, phospholipase A1, serine carboxypeptidase and serine protease revealed that one species was positioned as an out group or grouped with another branch in a way not coincident with the evolutionary pattern of Aculeate species (Figures S16–S24). For example, S. deforme was clustered as an out group in acid phosphatase (Figure S16), and Sphecidae sp. was grouped with Formicidae in the phylogenetic tree of defensin 2 (Figure S17). Thus, these venom proteins only appear to semi-reflect the phylogenetic pattern of Aculeate species, and in the cases when venom proteins exhibiting an isolated phylogeny, they likely have evolved to have unique functions. In the phylogenetic trees of arginine kinase, defensin 1, phospholipase A2, tachykinin, and venom allergen 5, two to five species were grouped with other branches or positioned independently (Figures S25–S29). An amino acid sequence alignment showed a high degree of sequence conservancy, but the phylogeny of venom allergen 5 did not clearly match the species phylogeny (Figure S29). Two solitary hunting wasps that belong to different clades in the Aculeate phylogeny were grouped together. Likewise, social wasps and B. ussurensis were grouped together. Thus, venom allergen 5 is a major allergen in the venom of social wasps, but it can be speculated that its evolutionary history is different across different Hymenopteran species. However, as only a limited number of wasps and bumblebee species were compared in this study, a larger scale observation should be conducted to confirm this notion.
The unique evolutionary patterns of these venom proteins or allergens might be caused by specific ecological behaviors. Adaptive traits of the venom could provide a selective advantage in some species, so they exhibit differential evolutionary patterns through the phylogenetic trees. To confirm this idea, however, a robust phylogenetic comparative analysis with larger numbers of species should be conducted. Further studies to identify unique traits, such as ecological behavior or function of the venom, would provide a clearer understanding of venom evolution.
We conducted venom gland—specific transcriptome analyses in 14 Aculeate bees and wasps and found novel venomic genes that show conserved domains with annotated major allergens. These novel genes are potent allergens to be verified and could provide fundamental data for venom immunotherapy (VIT). VIT is accepted as a safe and effective treatment. However, due to genetic differences in Hymenoptera species, side effects can be frequent, and success rates vary [41,52]. To improve the safety and efficacy of Hymenoptera VIT, use of recombinant venom allergens that have domains highly conserved in bees and wasps could provide a solution for systemic allergic side effects. The venom proteins in some Aculeate bees and wasps possessed average amino acid sequence similarities of 38% to 96%. Arginine kinase showed the highest percentage similarity among all venom proteins, exhibiting 92–100% sequence similarity. Because arginine kinase showed considerable expression in the venom glands of solitary and social wasps and bumblebees (Figure 1), use of an epitope region in conserved domains and a recombinant protein could improve the safety and quality of diagnosis approaches and Hymenoptera VIT.

3. Conclusions

Acid phosphatase, neprilysin, and tachykinin, which possess neurotoxic activity, were most highly expressed in solitary hunting wasps, suggesting that they use their venom to paralyze and preserve prey to ensure fresh provisions for their progeny. Most major venom allergens and pain-producing factors (arginine kinase, dipeptidyl peptidase 4, endocuticle structural glycoprotein, hyaluronidase, mastoparan, phospholipase A1, phospholipase A2, serine protease, serine protease inhibitor, venom allergen 5, and vitellogenin) exhibited high expression in social wasps, suggesting that social wasps may use venom as a defensive tool against intruding species and possesses high toxicity. Bumblebees showed few transcripts for allergens or pain-producing factors but the highest expression of an antimicrobial peptide, defensin 1, suggesting that bumblebees possess venom with low toxicity, and that they use their venom for sanitizing purposes. Taken together, these results indicate that some Aculeate bees and wasps tested exhibit distinct venom properties and differential venomic gene expression patterns depending on species and sociality. Some venom components (carboxylesterase 6, dipeptidyl peptidase 4, hyaluronidase, icarapin, phospholipase B, serine carboxypeptidase, serine protease inhibitor, and vitellogenin) reflected the Aculeate species phylogeny, but others (acid phosphatase, arginine kinase, defensin 1, defensin 2, major royal jelly protein, mastoparan, metalloproteinase, phospholipase A1, phospholipase A2, serine protease, tachykinin, and venom allergen 5) did not. It can be speculated that this unique evolution of some venom components might have been shaped by specific ecological behaviors. To confirm this notion, however, a phylogenetic comparative analysis with larger numbers of species should be conducted. To the best of our knowledge, this is the first systematic study of venom gland transcriptomic data to compare major venom components across various species of some Aculeate bumblebees and wasps. A broader scale of comparison for venom components (including minor components across diverse bee and wasp species) would elucidate how venoms evolved according to ecological and behavioral pressures.

4. Materials and Methods

4.1. Wasp and Bumblebee Collection and Total RNA Purification

Female solitary hunting wasps (E. decoratus, Sphecidae sp., A. flavomarginatum and S. deforme), social wasps (V. analis, V. crabro, V. dybowskii, V. simillima, P. varia, P. snelleni, and P. rothneyi), and bumblebees (B. ardens, B. consobrinus, and B. ussurensis) were collected from several southern regions in Korea. Wasps and bumblebee species were collected and placed in 50 mL tube for 2 h at room temperature to unify the stress conditions of collected specimens. Wasps and three bumblebee species were then anesthetized using low-pressure carbon dioxide. Venom glands were only dissected from the venom apparatus, and total RNA was extracted from 20 dissected venom glands using 200 μL of TRI reagent (Molecular Research Center, Cincinnati, OH, USA) according to the manufacturer’s protocol. As the main goal of this study was to overview the overall expression patterns of putative venom genes across several bumblebees and wasp species, no statistical treatment based on multiple transcriptome analyses were made. Instead, we used the same conditions of venom gland dissection and total RNA extraction for all bumblebee and wasp species to minimize the temporal expression bias between different species.

4.2. Construction of an in Silico cDNA Library and Computational Analysis

RNA purity and total RNA integrity were evaluated using a NanoDrop 8000 spectrophotometer (Thermo Fisher Scientific, Waltham MA, USA) and Agilent Technologies 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). De novo assembly of Illumina reads was performed based on a previously reported method [53]. The raw reads from RNA-Seq were filtered using the FASTP program [54] to remove low-quality reads, which were base called with an error rate higher than 0.1% (Q-score < 30). The quality-filtered clean reads were subjected to the Trinity de novo assembler v2.8.4 program (GitHub, San Francisco, CA, USA) [51] to construct contigs. To streamline the contigs, sequences with high homology were clustered by the CD-HIT-EST program [55]. Subsequently, the clustered sequences with protein coding capability were selected by the standalone version of the TransDecoder v5.5 program, which was introduced with Trinity [51], to obtain the final version of the in silico cDNA library. The cDNA sequences were annotated by the BLASTX program with the GenBank reference peptide nr database. The Illumina short reads were mapped to the cDNA sequences using the Kallisto program [56] to calculate the expression rate in transcripts per kilobase million (TPM). To identify genes encoding secreted venom components, the top 100 highly expressed genes with signal peptide sequences of 14 Aculeate bees and wasps (E. decoratus, Sphecidae sp., A. flavomarginatum, S. deforme, V. dybowskii, V. simillima, P. varia, P. snelleni, P. rothneyi, B. ardens, B. consobrinus, and B. ussurensis) were predicted using the Phobius program [57] and SignalP-5.0 server (www.cbs.dtu.dk/services/SignalP/). For the USPs, an allergen database search (www.allergenonline.org) was conducted to predict potential allergenicity. The raw data of transcriptome sequencing from following species have been uploaded to NCBI’s Short Read Archive; E. decoratus (SRR10286816), Sphecidae sp. (SRR10286815), A. flavomarginatum (SRR10286814), S. deforme (SRR10286825), V. crabro (SRR10286826), V. analis (SRR10286827), V. dybowskii (SRR10286821), V. simillima (SRR10286820), P. varia (SRR10286819), P. snelleni (SRR10286817), P. rothneyi (SRR10286818), B. ardens (SRR10286824), B. consobrinus (SRR10286823), and B. ussurensis (SRR10286822).

4.3. Ortholog Analysis of Venom Component Genes

We used the OrthoMCL program v2.0.9 [58] with the Markov Cluster algorithm to identify the orthologous paralog group of venom component genes across wasp and bumblebee species. To do this, we used the BLAST v2.2.28+ program (NCBI, Bethesda, MD, USA) before self-alignment of all proteins, with an e-value cutoff of 1e−5 and a minimum alignment length of 50 bp. False-positives were removed by filtering at least 50% of the query cover. Domain analysis was performed with the HMMER v3.1b2 [59] program using the Pfam [60] (http://pfam.sanger.ac.uk) database to identify additional functions of each protein.
A reference gene (dimethyladenosine transferase) [61] was selected for verification of transcriptional abundance, and total reads were divided and normalized using the average TPM values for the reference gene in the 14 Aculeate bees and wasps (Table S15). Hierachial clustering analysis for the 22 orthologous venom genes was conducted with UPGMA (unweighted-pair group method with arithmetic mean) clustering method and the genetic distance was calculated by Nei’s DA distance measurement [62] using POPTREE2 software. The percentages of bootstrap values based on 1000 replicates are shown on each node. Heat map was generated using the heat map function in TIBCO Spotfire® software (TIBCO Software Inc., Palo Alto, CA, USA).

4.4. Construction of a Phylogenetic Tree

To investigate structural differences in venom proteins, we aligned the deduced amino acid sequences obtained from the transcriptome data with full-length sequences using CLC Main Workbench 7 (CLC Bio, Waltham, MA, USA). Based on amino acid alignment of the orthologous venom genes, we constructed a phylogenetic tree using maximum likelihood method based on the JTT matrix-based model [63] with 1000 bootstrap replications. The percentages of trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) were shown on each node. The tree was drawn to scale, with branch lengths measured in the number of substitutions per site. Phylogenetic analyses were conducted in MEGA X (Pennsylvania State University, PA, USA). Acid phosphatase (AF321918), defensin 1 (AAA52303), defensin 2 (AAC69554), hyaluronidase (AAC70915), metalloproteinase (AAA58658), phospholipase A2 (AAF09020), and serine protease (CCA61110) of Homo sapiens were used as outgroups. The following were also used as outgroups: arginine kinase (JAV48307), phospholipase A1 (JAV47955), and phospholipase B (JAV47950) of Hadrurus spadix; mastoparan of an uncultured bacterium (AP032449) [64]; venom allergen 5 (XP_022659921) of Varroa destructor; carboxylesterase 6 (KK119399) of Stegodyphus mimosarum; icarapin (FX985504) of Odontomachus monticola; dipeptidyl peptidase 4 (XP_021001451) and major royal jelly protein (XP_021000670) of Parasteatoda tepidariorum; and endocuticle structural glycoprotein (PSN39565), neprilysin (PSN37766), serine carboxypeptidase (PSN48057), serine protease inhibitor (PSN55974), tachykinin (PSN47497), and vitellogenin (CAA06379) of Blattella germanica. Major venomic genes of other Hymenopterans, Dipterans, Hemipterans, Lepidopterans and Dictyopterans were included and used as outgroups. Any robust phylogenetic comparative analysis was not feasible due to the limited number of Aculeate wasps and bumblebee species tested in this study.

4.5. Quantitative Real-Time PCR

DNase I (Takara, Kyoto, Japan) was used to remove DNA contamination of total RNA of two social wasps (P. rothneyi and P. snelleni), one solitary hunting wasp (E. decoratus) and one bumblebee (B. ardens). cDNA was synthesized using SuperScriptIV reverse transcriptase (Invitrogen, Carlsbad, CA, USA). Total 10 venomic genes (arginine kinase, defensin 1, dipeptidyl peptidase 4, hyaluronidase, icarapin, phospholipase A2, serine protease inhibitor, tachykinin, vitellogenin, and neprilysin) and a reference housekeeping gene (dimethyladenosine transferase) were selected for verification of transcriptional abundance. Venomic gene specific primers were designed to have similar lengths and %GC content (Table S16) and qPCR was conducted using a previously described protocol [53].

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/12/1/47/s1, Figure S1: Comparison of the relative transcription levels and TPM values of (A) arginine kinase, (B) defensin 1, (C) dipeptidyl peptidase 4, (D) hyaluronidase, (E) icarapin, (F) phospholipase A2, (G) serine protease inhibitor, (H) tachykinin, (I) vitallogenin and (J) neprilysin from P. rothneyi, P. snelleni, E. decoratus and B. ardens, Figure S2: Amino acid alignments of uncharacterized protein 1 from V. crabro, V. simillima, P. varia, P. snelleni, B. consobrinus and M. banksi, Figure S3: Amino acid alignments of uncharacterized protein 2 from V. analis, V. dybowskii, P. varia, A. flavomarginatum and Crassostrea gigas, Figure S4: Amino acid alignments of uncharacterized protein 3 from V. dybowskii, B. consobrinus, B. ussurensis and Hevea brasiliensis, Figure S5: Amino acid alignments of uncharacterized protein 4 from B. ardens, B. ussurensis and Triticum aestivum, Figure S6: Amino acid alignments of uncharacterized protein 5 from P. varia, A. flavomarginatum and S. geminate, Figure S7: Amino acid alignments of uncharacterized protein 6 from V. dybowskii and A. flavomarginatum, Figure S8: Amino acid alignments of carboxylesterase 6. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. snelleni, P. rothneyi, A. flavomarginatum, Sphecidae sp., B. ussurensis and S. mimosarum. (B) Phylogenetic analysis of carboxylesterase 6, Figure S9: Amino acid alignments of dipeptidyl peptidase 4. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. rothneyi, A. flavomarginatum, S. deforme, Sphecidae sp., B. ardens, B. consobrinus, B. ussurensis and P. tepidariorum. (B) Phylogenetic analysis of dipeptidyl peptidase 4, Figure S10: Amino acid alignments of endocuticle structural glycoprotein. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. snelleni, P. rothneyi, A. flavomarginatum, E. decoratus, S. deforme, B. consobrinus, B. ussurensis and Blattella germanica. (B) Phylogenetic analysis of endocuticle structural glycoprotein, Figure S11: Amino acid alignments of icarapin. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. snelleni, P. rothneyi, A. flavomarginatum, E. decoratus, S. deforme, Sphecidae sp., B. ardens, B. consobrinus, B. ussurensis and O. monticola. (B) Phylogenetic analysis of icarapin, Figure S12: Amino acid alignments of major royal jelly protein. (A) Alignment of amino acid sequences from P. varia, P. rothneyi, E. decoratus, S. deforme, B. ardens, B. consobrinus, B. ussurensis and P. tepidariorum. (B) Phylogenetic analysis of major royal jelly protein, Figure S13: Amino acid alignments of phospholipase B. (A) Alignment of amino acid sequences from V. analis, P. varia, P. snelleni, P. rothneyi, E. decoratus, Sphecidae sp., B. ardens, B. consobrinus, B. ussurensis and H. spadix. B) Phylogenetic analysis of phospholipase B, Figure S14: Amino acid alignments of serine protease inhibitor. A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. snelleni, P. rothneyi, A. flavomarginatum, E. decoratus, S. deforme, Sphecidae sp., B. ardens, B. consobrinus and B. germanica. B() Phylogenetic analysis of serine protease inhibitor, Figure S15: Amino acid alignments of vitellogenin. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. snelleni, P. rothneyi, A. flavomarginatum, Sphecidae sp., B. ardens, B. consobrinus, B. ussurensis and B. germanica. (B) Phylogenetic analysis of vitellogenin, Figure S16: Amino acid alignments of acid phosphatase. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, P. varia, P. snelleni, A. flavomarginatum, S. deforme, Sphecidae sp., B. ardens, B. consobrinus, B. ussurensis and H. sapiens. (B) Phylogenetic analysis of acid phosphatase, Figure S17: Amino acid alignments of defensin 2. (A) Alignment of amino acid sequences from P. varia, P. snelleni, P. rothneyi, A. flavomarginatum, E. decoratus, Sphecidae sp. and H. sapiens. (B) Phylogenetic analysis of defensin 2, Figure S18: Amino acid alignments of hyaluronidase. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. snelleni, P. rothneyi, A. flavomarginatum, S. deforme, Sphecidae sp., B. ardens, B. ussurensis and H. sapiens. (B) Phylogenetic analysis of hyaluronidase, Figure S19: Amino acid alignments of mastoparan. (A) Alignment of amino acid sequences from V. crabro, V. analis, P. snelleni, P. rothneyi, A. flavomarginatum, S. deforme, B. ussurensis and uncultured bacterium. (B) Phylogenetic analysis of mastoparan, Figure S20: Amino acid alignments of metalloproteinase. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. snelleni, P. rothneyi, A. flavomarginatum, E. decoratus, S. deforme, Sphecidae sp., B. ardens, B. ussurensis and H. sapiens. (B) Phylogenetic analysis of metalloproteinase, Figure S21: Amino acid alignments of neprilysin. (A) Alignment of amino acid sequences from A. flavomarginatum, E. decoratus, B. consobrinus and B. germanica. (B) Phylogenetic analysis of neprilysin, Figure S22: Amino acid alignments of phospholipase A1. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. rothneyi, E. decoratus, B. consobrinus, B. ussurensis and H. spadix. (B) Phylogenetic analysis of phospholipase A1, Figure S23: Amino acid alignments of serine carboxypeptidase. (A) Alignment of amino acid sequences from V. crabro, V. analis, P. varia, P. rothneyi, A. flavomarginatum, E. decoratus, S. deforme, Sphecidae sp., B. ardens, B. consobrinus and B. germanica. (B) Phylogenetic analysis of serine carboxypeptidase, Figure S24: Amino acid alignments of serine protease. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, P. varia, P. snelleni, P. rothneyi, A. flavomarginatum, E. decorates, S. deforme, Sphecidae sp., B. ardens, B. ussurensis and H. sapiens. (B) Phylogenetic analysis of serine protease, Figure S25: Amino acid alignments of arginine kinase. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. snelleni, A. flavomarginatum, E. decoratus, S. deforme, Sphecidae sp., B. consobrinus and Hadrurus spadix. (B) Phylogenetic analysis of arginine kinase, Figure S26: Amino acid alignments of defensin 1. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, P. varia, P. rothneyi, A. flavomarginatum, Sphecidae sp., B. ardens, B. consobrinus, B. ussurensis and H. sapiens. (B) Phylogenetic analysis of defensin 1, Figure S27: Amino acid alignments of phospholipase A2. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. dybowskii, V. simillima, P. varia, P. rothneyi, A. flavomarginatum, Sphecidae sp., B. ardens, B. ussurensis and H. sapiens. (B) Phylogenetic analysis of phospholipase A2, Figure S28: Amino acid alignments of tachykinin. (A) Alignment of amino acid sequences from V. analis, V. dybowskii, P. varia, P. snelleni, P. rothneyi, E. decorates, Sphecidae sp., B. ardens, B. consobrinus, B. ussurensis, and B. germanica. (B) Phylogenetic analysis of tachykinin, Figure S29: Amino acid alignments of venom allergen 5. (A) Alignment of amino acid sequences from V. crabro, V. analis, V. simillima, P. snelleni, P. rothneyi, A. flavomarginatum, S. deforme, Sphecidae sp., B. ussurensis and Varroa destructor. (B) Phylogenetic analysis of venom allergen 5, Table S1: Annotation of top 100 highly expressed genes in the venom gland of Vespa analis, Table S2: Annotation of top 100 highly expressed genes in the venom gland of Vespa crabro, Table S3: Annotation of top 100 highly expressed genes in the venom gland of Vespa dybowskii, Table S4: Annotation of top 100 highly expressed genes in the venom gland of Vespa simillima, Table S5: Annotation of top 100 highly expressed genes in the venom gland of Parapolybia varia, Table S6: Annotation of top 100 highly expressed genes in the venom gland of Polistes rothneyi, Table S7: Annotation of top 100 highly expressed genes in the venom gland of Polistes snelleni, Table S8: Annotation of top 100 highly expressed genes in the venom gland of Eumenes decoratus, Table S9: Annotation of top 100 highly expressed genes in the venom gland of Sphecidae sp., Table S10: Annotation of top 100 highly expressed genes in the venom gland of Anterhynchium flavomarginatum, Table S11: Annotation of top 100 highly expressed genes in the venom gland of Sceliphron deforme, Table S12: Annotation of top 100 highly expressed genes in the venom gland of Bombus ardens, Table S13: Annotation of top 100 highly expressed genes in the venom gland of Bombus consobrinus, Table S14: Annotation of top 100 highly expressed genes in the venom gland of Bombus ussurensis, Table S15: TPM values of reference housekeeping gene dimethyladenosine transferase in 14 Aculeate bee and wasp species, Table S16: Primers used in quantitative real-time PCR.

Author Contributions

Conceptualization, K.A.Y. and S.H.L.; methodology, K.A.Y. and S.H.L.; software, W.-J.K.; validation, K.A.Y.; formal analysis, K.A.Y. and W.-J.K.; investigation, K.A.Y.; resources, K.K., W.Y.B., N.-H.A., C.-H.B. and J.-H.Y.; data curation, K.A.Y. and W.-J.K.; writing—original draft preparation, K.A.Y. and S.H.L.; writing—review and editing, K.A.Y. and S.H.L.; visualization, K.A.Y.; supervision, S.H.L.; project administration, W.Y.B., N.-H.A., C.-H.B. and J.-H.Y.; funding acquisition, W.Y.B., N.-H.A., C.-H.B. and J.-H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant from the National Institute of Biological Resources (NIBR), funded by the Ministry of Environment (MOE), Republic of Korea (NIBR201711201, NIBR201909201).

Acknowledgments

Kyungmun Kim was supported in part by the Brain Korea 21 Plus Program.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Piek, T. Venoms of the Hymenoptera: Biochemical, Pharmacological and Behavioural Aspects; Elsevier: Amsterdam, The Netherlands, 2013. [Google Scholar]
  2. Gauld, I.D.; Bolton, B.; Huddleston, T.; Fitton, M.; Shaw, M.R.; Noyes, J.; Day, M.; Else, G.; Ferguson, N.; Ward, S. The hymenoptera. Hymenoptera 1988, 1988, 1–332. [Google Scholar]
  3. Dowton, M. Simultaneous analysis of 16S, 28S, COI and morphology in the Hymenoptera: Apocrita—Evolutionary transitions among parasitic wasps. Biol. J. Linn. Soc. 2001, 74, 87–111. [Google Scholar] [CrossRef]
  4. O’Neill, K.M. Solitary Wasps: Behavior and Natural History; Cornell University Press: Ithaca, NY, USA, 2001. [Google Scholar]
  5. Lee, S.H.; Baek, J.H.; Yoon, K.A. Differential Properties of Venom Peptides and Proteins in Solitary vs. Social Hunting Wasps. Toxins 2016, 8, 32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Evans, H.; Eberhard, M.J. The Wasps; University of Michigan Press: Ann Arbor, MI, USA, 1971. [Google Scholar]
  7. Schmidt, J.O. Insect Defenses: Adaptive Mechanisms and Strategies of Prey and Predators; SUNY Press: Albany, NY, USA, 1990. [Google Scholar]
  8. Akre, R.D. Nesting biology and behavior of the baldfaced hornet, Dolichovespula maculata (L.)(Hymenoptera: Vespidae) in the Pacific Northwest. Melanderia 1992, 48, 1–33. [Google Scholar]
  9. Jeanne, R.L.; Williams, N.M.; Yandell, B.S. Age polyethism and defense in a tropical social wasp (Hymenoptera: Vespidae). J. Insect Behav. 1992, 5, 211–227. [Google Scholar] [CrossRef] [Green Version]
  10. Spradbery, J.P. Wasps. In An Account of the Biology and Natural History of Social and Solitary Wasps, with Particular Reference to Those of the British Isles; Sidgwick & Jackson: London, UK, 1973. [Google Scholar]
  11. Richter, M.R. Social wasp (Hymenoptera: Vespidae) foraging behavior. Annu. Rev. Entomol. 2000, 45, 121–150. [Google Scholar] [CrossRef]
  12. Eno, A.E. Pharmacological investigation of oedema induced by venom from the wasp Polistes fuscatus. Toxicon 1997, 35, 1691–1698. [Google Scholar] [CrossRef]
  13. Fisher, B.A.; Antonios, T.F. Atrial flutter following a wasp sting. J. Postgrad. Med. 2003, 49, 254–255. [Google Scholar]
  14. Ferrari, S.; Pietroiusti, A.; Galanti, A.; Compagnucci, M.; Fontana, L. Paroxysmal atrial fibrillation after insect sting. J. Allergy Clin. Immunol. 1996, 98, 759–761. [Google Scholar] [CrossRef]
  15. Sachdev, A.; Mahapatra, M.; D’Cruz, S.; Kumar, A.; Singh, R.; Lehl, S.S. Wasp sting induced neurological manifestations. Neurol. India 2002, 50, 319–321. [Google Scholar]
  16. Sharmila, R.R.; Chetan, G.; Narayanan, P.; Srinivasan, S. Multiple organ dysfunction syndrome following single wasp sting. Indian J. Pediatr. 2007, 74, 1111–1112. [Google Scholar] [CrossRef]
  17. Argiolas, A.; Pisano, J.J. Facilitation of Phospholipase-A2 Activity by Mastoparans, a New Class of Mast-Cell Degranulating Peptides from Wasp Venom. J. Biol. Chem. 1983, 258, 3697–3702. [Google Scholar]
  18. Krishnakumari, V.; Nagaraj, R. Antimicrobial and hemolytic activities of crabrolin, a 13-residue peptide from the venom of the European hornet, Vespa crabro, and its analogs. J. Pept. Res. 1997, 50, 88–93. [Google Scholar] [CrossRef]
  19. Bucher, C.; Korner, P.; Wüthrich, B. Allergy to bumblebee venom. Curr. Opin. Allergy Clin. Immunol. 2001, 1, 361–365. [Google Scholar] [CrossRef]
  20. Moreau, S.J.M.; Asgari, S. Venom Proteins from Parasitoid Wasps and Their Biological Functions. Toxins 2015, 7, 2385–2412. [Google Scholar] [CrossRef] [Green Version]
  21. Yoon, K.A.; Kim, K.; Nguyen, P.; Seo, J.B.; Park, Y.H.; Kim, K.G.; Seo, H.Y.; Koh, Y.H.; Lee, S.H. Comparative bioactivities of mastoparans from social hornets Vespa crabro and Vespa analis. J. Asia-Pac. Entomol. 2015, 18, 825–829. [Google Scholar] [CrossRef]
  22. Baek, J.H.; Lee, S.H. Identification and characterization of venom proteins of two solitary wasps, Eumenes pomiformis and Orancistrocerus drewseni. Toxicon 2010, 56, 554–562. [Google Scholar] [CrossRef]
  23. Sunagar, K.; Morgenstern, D.; Reitzel, A.M.; Moran, Y. Ecological venomics: How genomics, transcriptomics and proteomics can shed new light on the ecology and evolution of venom. J. Proteom. 2016, 135, 62–72. [Google Scholar] [CrossRef]
  24. Rokyta, D.R.; Margres, M.J.; Calvin, K. Post-transcriptional Mechanisms Contribute Little to Phenotypic Variation in Snake Venoms. G3-Genes Genom. Genet. 2015, 5, 2375–2382. [Google Scholar] [CrossRef] [Green Version]
  25. Zhu, J.Y.; Ye, G.Y.; Hu, C. Molecular cloning and characterization of acid phosphatase in venom of the endoparasitoid wasp Pteromalus puparum (Hymenoptera: Pteromalidae). Toxicon 2008, 51, 1391–1399. [Google Scholar] [CrossRef]
  26. Dani, M.P.; Edwards, J.P.; Richards, E.H. Hydrolase activity in the venom of the pupal endoparasitic wasp, Pimpla hypochondriaca. Comp. Biochem. Phys. B 2005, 141, 373–381. [Google Scholar] [CrossRef]
  27. Turner, A.J.; Isaac, R.E.; Coates, D. The neprilysin (NEP) family of zinc metalloendopeptidases: Genomics and function. Bioessays 2001, 23, 261–269. [Google Scholar] [CrossRef]
  28. Isaac, R.E. Neuropeptide-Degrading Endopeptidase Activity of Locust (Schistocerca-Gregaria) Synaptic-Membranes. Biochem. J. 1988, 255, 843–847. [Google Scholar] [CrossRef] [Green Version]
  29. Arvidson, R.; Kaiser, M.; Pan, S.Q.; Libersat, F.; Adams, M.E. Bioinformatic and Functional Analysis of Venom from the Jewel Wasp Ampulex compressa. Faseb J. 2016, 30, 819–821. [Google Scholar]
  30. Binder, M.; Mahler, V.; Hayek, B.; Sperr, W.R.; Scholler, M.; Prozell, S.; Wiedermann, G.; Valent, P.; Valenta, R.; Duchene, M. Molecular and immunological characterization of arginine kinase from the indianmeal moth, Plodia interpunctella, a novel cross-reactive invertebrate pan-allergen. J. Immunol. 2001, 167, 5470–5477. [Google Scholar] [CrossRef] [Green Version]
  31. Schiener, M.; Hilger, C.; Eberlein, B.; Pascal, M.; Kuehn, A.; Revets, D.; Planchon, S.; Pietsch, G.; Serrano, P.; Moreno-Aguilar, C.; et al. The high molecular weight dipeptidyl peptidase IV Pol d 3 is a major allergen of Polistes dominula venom. Sci. Rep. 2018, 8. [Google Scholar] [CrossRef]
  32. Hino, K.; Yamamoto, S.; Sano, O.; Taniguchi, Y.; Kohno, K.; Usui, M.; Fukuda, S.; Hanzawa, H.; Haruyama, H.; Kurimoto, M. Carbohydrate Structures of the Glycoprotein Allergen Cry-J-I from Japanese Cedar (Cryptomeria-Japonica) Pollen. J. Biochem. 1995, 117, 289–295. [Google Scholar] [CrossRef]
  33. Borodina, I.; Jensen, B.M.; Wagner, T.; Hachem, M.A.; Sondergaard, I.; Poulsen, L.K. Expression of Enzymatically Inactive Wasp Venom Phospholipase A1 in Pichia pastoris. PLoS ONE 2011, 6, e21267. [Google Scholar] [CrossRef] [Green Version]
  34. Corthesy, B.; Jilek, S.; Barbey, C.; Spertini, F. Antigen-independent suppression of the allergic immune response to bee venom phospholipase A2 by DNA vaccination in CBA/J mice. J. Allergy Clin. Immunol. 2001, 107, S325. [Google Scholar]
  35. Hoffman, D.R.; Jacobson, R.S. Allergens in Hymenoptera venom. XXVII: Bumblebee venom allergy and allergens. J. Allergy Clin. Immunol. 1996, 97, 812–821. [Google Scholar] [CrossRef]
  36. Winningham, K.M.; Fitch, C.D.; Schmidt, M.; Hoffman, D.R. Hymenoptera venom protease allergens. J. Allergy Clin. Immunol. 2004, 114, 928–933. [Google Scholar] [CrossRef] [PubMed]
  37. Golden, D.B.; Kelly, D.; Hamilton, R.G.; Craig, T.J. Venom immunotherapy reduces large local reactions to insect stings. J. Allergy Clin. Immunol. 2009, 123, 1371–1375. [Google Scholar] [CrossRef] [PubMed]
  38. Michel, Y.; McIntyre, M.; Ginglinger, H.; Ollert, M.; Cifuentes, L.; Blank, S.; Spillner, E. The Putative Serine Protease Inhibitor Api m 6 From Apis mellifera Venom: Recombinant and Structural Evaluation. J. Investig. Allerg. Clin. 2012, 22, 476–484. [Google Scholar]
  39. King, T.P.; Spangfort, M.D. Structure and biology of stinging insect venom allergens. Int. Arch. Allergy Imm. 2000, 123, 99–106. [Google Scholar] [CrossRef] [PubMed]
  40. Blank, S.; Seismann, H.; McIntyre, M.; Ollert, M.; Wolf, S.; Bantleon, F.I.; Spillner, E. Vitellogenins Are New High Molecular Weight Components and Allergens (Api m 12 and Ves v 6) of Apis mellifera and Vespula vulgaris Venom. PLoS ONE 2013, 8, e62009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Blank, S.; Seismann, H.; Bockisch, B.; Braren, I.; Cifuentes, L.; McIntyre, M.; Ruhl, D.; Ring, J.; Bredehorst, R.; Ollert, M.W.; et al. Identification, Recombinant Expression, and Characterization of the 100 kDa High Molecular Weight Hymenoptera Venom Allergens Api m 5 and Ves v 3. J. Immunol. 2010, 184, 5403–5413. [Google Scholar] [CrossRef] [Green Version]
  42. Peiren, N.; de Graaf, D.C.; Brunain, M.; Bridts, C.H.; Ebo, D.G.; Stevens, W.J.; Jacobs, F.J. Molecular cloning and expression of icarapin, a novel IgE-binding bee venom protein. FEBS Lett. 2006, 580, 4895–4899. [Google Scholar] [CrossRef] [Green Version]
  43. Blank, S.; Seismann, H.; Michel, Y.; McIntyre, M.; Cifuentes, L.; Braren, I.; Grunwald, T.; Darsow, U.; Ring, J.; Bredehorst, R.; et al. Api m 10, a genuine A. mellifera venom allergen, is clinically relevant but underrepresented in therapeutic extracts. Allergy 2011, 66, 1322–1329. [Google Scholar] [CrossRef]
  44. Takeda, S.; Takeya, H.; Iwanaga, S. Snake venom metalloproteinases: Structure, function and relevance to the mammalian ADAM/ADAMTS family proteins. Biochim. Biophys. Acta 2012, 1824, 164–176. [Google Scholar] [CrossRef]
  45. Orsolic, N. Bee venom in cancer therapy. Cancer Metastasis Rev. 2012, 31, 173–194. [Google Scholar] [CrossRef]
  46. Bernheimer, A.W.; Linder, R.; Weinstein, S.A.; Kim, K.S. Isolation and characterization of a phospholipase B from venom of Collett’s snake, Pseudechis colletti. Toxicon 1987, 25, 547–554. [Google Scholar] [CrossRef]
  47. Matysiak, J.; Hajduk, J.; Pietrzak, L.; Schmelzer, C.E.H.; Kokot, Z.J. Shotgun proteome analysis of honeybee venom using targeted enrichment strategies. Toxicon 2014, 90, 255–264. [Google Scholar] [CrossRef] [PubMed]
  48. Bown, D.P.; Wilkinson, H.S.; Gatehouse, J.A. Midgut carboxypeptidase from Helicoverpa armigera (Lepidoptera: Noctuidae) larvae: Enzyme characterisation, cDNA cloning and expression. Insect Biochem. Mol. Biol. 1998, 28, 739–749. [Google Scholar] [CrossRef]
  49. Cho, W.L.; Deitsch, K.W.; Raikhel, A.S. An Extraovarian Protein Accumulated in Mosquito Oocytes Is a Carboxypeptidase Activated in Embryos. Proc. Natl. Acad. Sci. USA 1991, 88, 10821–10824. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Zhu, Y.Y.; Machleder, E.M.; Chenchik, A.; Li, R.; Siebert, P.D. Reverse transcriptase template switching: A SMART (TM) approach for full-length cDNA library construction. Biotechniques 2001, 30, 892–897. [Google Scholar] [CrossRef] [Green Version]
  51. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.D.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [Green Version]
  52. Muller, U.; Helbling, A.; Berchtold, E. Immunotherapy with Honeybee Venom and Yellow Jacket Venom IS Different Regarding Efficacy and Safety. J. Allergy Clin. Immunol. 1992, 89, 529–535. [Google Scholar] [CrossRef]
  53. Yoon, K.A.; Kim, K.; Nguyen, P.; Seo, J.B.; Park, Y.H.; Kim, K.G.; Seo, H.Y.; Koh, Y.H.; Lee, S.H. Comparative functional venomics of social hornets Vespa crabro and Vespa analis. J. Asia-Pac. Entomol. 2015, 18, 815–823. [Google Scholar] [CrossRef]
  54. Chen, S.F.; Zhou, Y.Q.; Chen, Y.R.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, 884–890. [Google Scholar] [CrossRef]
  55. Li, W.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef] [Green Version]
  56. Bray, N.L.; Pimentel, H.; Melsted, P.; Pachter, L. Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 2016, 34, 525–527. [Google Scholar] [CrossRef] [PubMed]
  57. Kall, L.; Krogh, A.; Sonnhammer, E.L.L. A combined transmembrane topology and signal peptide prediction method. J. Mol. Biol. 2004, 338, 1027–1036. [Google Scholar] [CrossRef] [PubMed]
  58. Li, L.; Stoeckert, C.J.; Roos, D.S. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 2003, 13, 2178–2189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Eddy, S.R. Profile hidden Markov models. Bioinformatics 1998, 14, 755–763. [Google Scholar] [CrossRef] [PubMed]
  60. Finn, R.D.; Mistry, J.; Schuster-Bockler, B.; Griffiths-Jones, S.; Hollich, V.; Lassmann, T.; Moxon, S.; Marshall, M.; Khanna, A.; Durbin, R.; et al. Pfam: Clans, web tools and services. Nucleic Acids Res. 2006, 34, D247–D251. [Google Scholar] [CrossRef] [Green Version]
  61. Gao, X.K.; Zhang, S.; Luo, J.Y.; Wang, C.Y.; Lu, L.M.; Zhang, L.J.; Zhu, X.Z.; Wang, L.; Lu, H.; Cui, J.J. Comprehensive evaluation of candidate reference genes for gene expression studies in Lysiphlebia japonica (Hymenoptera: Aphidiidae) using RT-qPCR. Gene 2017, 637, 211–218. [Google Scholar] [CrossRef]
  62. Takezaki, N.; Nei, M.; Tamura, K. POPTREE2: Software for Constructing Population Trees from Allele Frequency Data and Computing Other Population Statistics with Windows Interface. Mol. Biol. Evol. 2010, 27, 747–752. [Google Scholar] [CrossRef] [Green Version]
  63. Jones, D.T.; Taylor, W.R.; Thornton, J.M. The Rapid Generation of Mutation Data Matrices from Protein Sequences. Comput. Appl. Biosci. 1992, 8, 275–282. [Google Scholar] [CrossRef]
  64. Moore, A.M.; Ahmadi, S.; Patel, S.; Gibson, M.K.; Wang, B.; Ndao, I.M.; Deych, E.; Shannon, W.; Tarr, P.I.; Warner, B.B.; et al. Gut resistome development in healthy twin pairs in the first year of life (vol 3, 27, 2015). Microbiome 2015, 3. [Google Scholar] [CrossRef] [Green Version]
Toxins 12 00047 g001 550
Figure 1. Hierarchial clustering of the 22 major orthologous venom genes in the 14 Aculeate bumblebee and wasp species (left panel), and their differential expression profiles represented in a heat map (right panel). In the heat map, each row represents the venom gene analyzed, and each column represents the species tested. The percentages of bootstrap values based on 1000 pseudoreplicates are shown on each node of the dendrogram. Red rectangle indicates the clade clustered with most of the allergens and pain-producing factors, which showed significantly high expression levels in the social wasps.
Figure 1. Hierarchial clustering of the 22 major orthologous venom genes in the 14 Aculeate bumblebee and wasp species (left panel), and their differential expression profiles represented in a heat map (right panel). In the heat map, each row represents the venom gene analyzed, and each column represents the species tested. The percentages of bootstrap values based on 1000 pseudoreplicates are shown on each node of the dendrogram. Red rectangle indicates the clade clustered with most of the allergens and pain-producing factors, which showed significantly high expression levels in the social wasps.
Toxins 12 00047 g001
Table
Table 1. Summary of the venom gland transcriptome cDNA libraries of 14 Aculeate bees/wasps.
Table 1. Summary of the venom gland transcriptome cDNA libraries of 14 Aculeate bees/wasps.
SpeciesTotal Number of ReadsTotal Base PairsTrimmed ReadsTotal Number of
Assembled Transcripts		Total Number of
Annotated Transcripts
E. decorates58,988,4105,957,829,41055,681,36018,06216,357
Sphecidae sp.63,954,1926,459,373,39260,641,23631,13427,423
A. flavomarginatum71,380,8787,209,468,67867,788,99021,51219,033
S. deforme57,404,4685,797,851,26852,911,92024,48320,097
V. crabro62,246,2645,756,774,99756,318,83011,09710,321
V. analis61,648,0785,686,981,51955,419,96812,53111,284
V. dybowskii113,343,13411,447,656,534103,185,33221,22719,627
V. simillima114,592,94611,573,887,54689,279,80229,15027,065
P. varia82,025,3868,284,563,98675,702,93620,37818,539
P. snelleni80,798,2168,160,649,81675,374,54027,74225,879
P. rothneyi81,060,7548,187,136,15476,096,05625,43623,475
B. ardens76,622,4047,738,862,80471,427,16831,36626,158
B. consobrinus109,962,55011,106,217,55072,432,17421,62520,050
B. ussurensis106,973,79410,804,353,19499,138,79426,42522,619
Table
Table 2. Expression levels of major secretory venom components in 14 Aculeate bees/wasps.
Table 2. Expression levels of major secretory venom components in 14 Aculeate bees/wasps.
DescriptionTPM Value of Top 100 Highly Expressed Genes a
V. anV. crV. dyV. siP. vaP. roP. snE. deSphA. flS. deB. arB. coB. us
Defensin 1 --1652------1163--
Endocuticle structural glycoprotein6209 -113244043779--4847-----
Hyaluronidase 10,22945661409-1275-8511------
Icarapin1495 1167725----28572081-203946236943
Neprilysin -----4174--1895---
Phospholipase A1 49,149106,39529,42364401246--------
Phospholipase A2 47221547952---------1157
Serine protease 1301-1737947-22,355------
Venom allergen 5 150,56765,46232,24113,212---------
Uncharacterized protein 1 2677-29,7004756-874-----6955-
Uncharacterized protein 217,769 1199-2611----865----
Uncharacterized protein 3 2125---------18051537
Uncharacterized protein 4 ---------1800-1741
Uncharacterized protein 5 1584------1944----
Uncharacterized protein 6 --2439----1322----
a Wasp species abbreviations: V. an, Vespa analis; V. cr, Vespa crabro; V. dy, Vespa dybowskii; V. si, Vespa simillima; P. va, Parapolybia varia; P. ro, Polistes rothneyi; P. sn, Polistes snelleni; E. de, Eumenes decoratus; Sph, Sphecidae sp.; A. fl, Anterhynchium flavomarginatum; S. de, Sceliphron deforme; B. ar, Bombus ardens; B. co, Bombus consobrinus; B. us, Bombus ussurensis.

Share and Cite

MDPI and ACS Style

Yoon, K.A.; Kim, K.; Kim, W.-J.; Bang, W.Y.; Ahn, N.-H.; Bae, C.-H.; Yeo, J.-H.; Lee, S.H. Characterization of Venom Components and Their Phylogenetic Properties in Some Aculeate Bumblebees and Wasps. Toxins 2020, 12, 47. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010047

AMA Style

Yoon KA, Kim K, Kim W-J, Bang WY, Ahn N-H, Bae C-H, Yeo J-H, Lee SH. Characterization of Venom Components and Their Phylogenetic Properties in Some Aculeate Bumblebees and Wasps. Toxins. 2020; 12(1):47. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010047

Chicago/Turabian Style

Yoon, Kyungjae Andrew, Kyungmun Kim, Woo-Jin Kim, Woo Young Bang, Neung-Ho Ahn, Chang-Hwan Bae, Joo-Hong Yeo, and Si Hyeock Lee. 2020. "Characterization of Venom Components and Their Phylogenetic Properties in Some Aculeate Bumblebees and Wasps" Toxins 12, no. 1: 47. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010047

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop