# Venom profile of the European carpenter bee Xylocopa violacea: Evolutionary and applied considerations on its toxin components

Björn M von Reumont, Sebastien Dutertre, Ivan Koludarov

► To cite this version:

Björn M von Reumont, Sebastien Dutertre, Ivan Koludarov. Venom profile of the European carpenter bee Xylocopa violacea: Evolutionary and applied considerations on its toxin components. Toxicon: X, 2022, 14, pp.100117. 10.1016/j.toxcx.2022.100117. hal-04185943

HAL Id: hal-04185943

https://hal.science/hal-04185943v1

Submitted on 23 Aug 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

Check for
updates

# Venom profile of the European carpenter bee *Xylocopa violacea*: Evolutionary and applied considerations on its toxin components

Björn M. von Reumont<sup>a, b, c, *</sup>, Sebastien Dutertre<sup>d</sup>, Ivan Koludarov<sup>b, **</sup>

## ARTICLE INFO

Handling Editor: Ray Norton

**Keywords:**

Bee venom

Xylopin

Solitary bees

Melittin

Apamin

Proteo-transcriptomics

## ABSTRACT

Modern venomics is increasing its focus on hymenopterans such as honeybees, bumblebees, parasitoid wasps, ants and true wasps. However solitary bees remain understudied in comparison and the few available venom studies focus on short melittin-like sequences and antimicrobial peptides. Herein we describe the first comprehensive venom profile of a solitary bee, the violet carpenter bee *Xylocopa violacea*, by using proteo-transcriptomics. We reveal a diverse and complex venom profile with 43 different protein families identified from dissected venom gland extracts of which 32 are also detected in the defensively injected venom. Melittin and apamin are the most highly secreted components, followed by Phospholipase A2, Icarapin, Secapin and three novel components. Other components, including eight novel protein families, are rather lowly expressed. We further identify multiple forms of apamin-like peptides. The melittin-like sequences of solitary bees separate into two clades, one comprised most sequences from solitary bees including xylopin (the variant in *Xylocopa*), while sequences from *Lasioglossa* appear closer related to melittin-like peptides from *Bombus* (*Bombolittins*). Our study suggests that more proteo-transcriptomic data from other solitary bees should be complemented with corresponding genome data to fully understand the evolution and complexity of bee venom proteins, and is of a particular need to disentangle the ambiguous phylogenetic relations of short peptides.

## 1. Introduction

Hymenopteran insects are the most species-rich animal group on earth including over one million estimated species (Oeyen et al., 2020; Peters et al., 2017). One of their most iconic evolutionary adaptation is the use of venom, which they employ predominately for defensive and predatory purposes (Kukuk et al., 1989; Piek, 1986; Schmidt, 1982; Schmidt et al., 1986). Scientific reports and observations on eusocial and solitary aculeate hymenopterans (distinct by their characteristic waist and the very maneuverable stinger) have been made since the 17th century and include original studies on their venom system and venom properties (Piek, 1986). These early venom studies focused mostly on species that live in close proximity to humans, in particular honeybees, eusocial wasps as well as parasitoid wasps (Leluk et al., 1989; Piek, 1986; Schmidt et al., 1986). Recently, proteomic, transcriptomic and

proteo-transcriptomic analyses have been published that describe in greater detail venom compositions of several hymenopteran species (Bouzid et al., 2013; Danneels et al., 2015; Liu et al., 2015; Ozbek et al., 2019; Robinson et al., 2018; Tan et al., 2020; Teng et al., 2017; Tou-chard et al., 2015, 2018; Wang et al., 2016; Yoon et al., 2020). However, the majority of these studies are restricted to a few species of parasitoid wasps, true wasps, ants, honeybees and bumblebees, see e.g. (Cerpes et al., 2021; de Graaf et al., 2009; Schmidt et al., 1986; Walker et al., 2018). Except for species of honeybees (Apis spp.) and bumblebees (Bombus spp.), bee venoms remain understudied. The long-standing focus on the eusocial honeybees and their close relatives (genera *Apis* and *Bombus*) overlooks that the solitary, wild bees are more speciose and evolutionary older (Branstetter et al., 2017; Peters et al., 2017; Sann et al., 2018). Despite of their tremendous contribution to plant polli-nation (Klein et al., 2018; Koh et al., 2016), their biology and ecology

<footnote><sup>a</sup> Goethe University Frankfurt, Institute for Cell Biology and Neuroscience, Department for Applied Bioinformatics, 60438, Frankfurt am Main, Germany</footnote>

<footnote><sup>b</sup> Justus Liebig University of Giessen, Institute for Insect Biotechnology, Heinrich-Buff-Ring 58, 35392, Giessen, Germany</footnote>

<footnote><sup>c</sup> LOEWE Centre for Translational Biodiversity Genomics (LOEWE TBG), Senckenberganlage 25, 60325, Frankfurt, Germany</footnote>

<footnote><sup>d</sup> IBMM, Université Montpellier, CNRS, ENSCM, 34095, Montpellier, France</footnote>

<footnote>* Corresponding author. Goethe University Frankfurt, Institute for Cell Biology and Neuroscience, Department for Applied Bioinformatics, 60438, Frankfurt am Main, Germany.</footnote>

<footnote>** Corresponding author.</footnote>

<footnote>E-mail address: bmvr@reumont.net (B.M. von Reumont).</footnote>
are still little studied, including their venoms. This taxonomical imbalance is also reflected in the number of curated venom components in the database UniProt (Jungo et al., 2012). From 355 currently manually reviewed venom peptides and proteins, only 12 come from six solitary bee species (Fig. 1).

The few existing studies on solitary bee venom are mostly application-driven and focus on the identification of single compounds in crude venom and their possible bioactivity. A major attention is paid to new antimicrobial peptides from solitary bees (Cujová et al., 2014; Kazuma et al., 2017; Ko et al., 2020; Monincová et al., 2010, 2012, 2014; Raghuraman and Chattopadhyay, 2007; Stöcklin et al., 2010) for which hymenopterans are particular known for (Moreau, 2013; Slaninova et al., 2012). The only proteomic study using whole solitary bee venom is the one of the Japanese carpenter bee *X. appendiculata circumvolans* (Kazuma et al., 2017), however it also focuses on shorter peptides. It may well be that this underrepresentation is linked to often lower specimen abundances and the small sizes of body and venom systems of most solitary bee species, which make them difficult to study. Nevertheless, to understand the origin and evolution of venom in bees, it is important to analyse venoms from older bee lineages and to include more solitary bee species.

The herein studied carpenter bee *Xylocopa violacea* grows up to 2.5 cm making it one of the largest solitary bees in Europe. With its impressive size and dominant blackish colour this solitary bee is often mistaken for a bumblebee or even a beetle. Characteristic nesting behaviour gives the group its name: carpenter bees carve their tubes and breeding chambers into living or dead wood with their large and powerful mandibles (Fig. 2A–C). Interesting in this context is that *X. violacea* appears to expand its range because of the climate change (Banaszak et al., 2019). Probably due to their size and striking habits, the first report on carpenter bee venom was published in 1865, in which the observation was made that its venom could kill a small bird within a few hours (Bert, 1865). And it has long been known that for humans, stings by *Xylocopa* can be very painful (Hardouin, 1948; Hermann and Mullen, 1974). Despite having this potent weapon, carpenter bees are not very aggressive and females only sting as a last resort, especially when protecting the brood. An ironic sidenote here is that the males behaviourally put on quite a show by darting and buzzing very menacingly towards potential aggressors despite actually being stingless and only mimicking a female’s stinging attack. The only human fatality reported following envenomation by *Xylocopa* was from Sri Lanka, but even that one case was later attributed to an unfortunate anaphylactic shock (Kularatne et al., 2016).

# 2. Material and methods

## 2.1. Collection of specimens and proteo-transcriptomic samples

Specimens of *X. violacea* (15 individuals) were collected in July and August 2019/2020 in the alluvial area of the River Wieseck in Giessen, Germany, with the permission to BMvR from the HNLUG Giessen (IV.2 R28). From each individual, the whole venom system composed of venom reservoir, venom duct and venom glands (Fig. 2 D) was dissected on ice under sterile conditions. Striking for *X. violacea* is that the venom reservoir is rather large and shows a very long and stretched tube-like shape in contrast to rather small and round or sack-like reservoirs seen in bees and wasps (Fig. 2 D).

The large and long Dufour gland was excluded and preserved for separate analyses. Crude venom was extracted from glands and venom reservoirs by squeezing with forceps in sterile ultrapure water (Thermo Fisher Scientific, Waltham, MA, USA) after prewashing twice to minimize haemolymph contamination. The tissue of the venom system was then preserved in RNAlater (Thermo Fisher Scientific) to generate transcriptomic data. RNA extraction, library preparation and short-read genome sequencing were outsourced to Macrogen (Seoul, Korea). RNA was extracted in Trizol and a low input protocol (Illumina Truseq) was

used for library preparation. Libraries were sequenced (150-bp paired-end reads) on an Illumina HiSeq 2500 platform (Macrogen). Complementary data of defensively injected venom was generated by holding individuals ($n = 3$) between large forceps and inducing them to sting into sterile parafilm-sealed centrifuge tubes filled with ultrapure water. Due to difficult seasonal variations in populations during collection period, these individuals were collected in France, near Montpellier.

## 2.2. Transcriptome assembly, protein prediction and annotation

The venom gland transcriptome was assembled with the Oyster River Pipeline v2.2.6 (MacManes, 2018) which operates multiple assemblers (Trinity, RNAspades with 55 and 75 kmer length and Shannon using a kmer-length of 75) applying standard settings in the provided Docker image (MacManes, 2018). For further details see (Koludarov et al., 2022). The assembly was then processed running Transdecoder (Minimum length ≥20 amino acids) to predict open reading frames (ORFs) for peptides (Haas, Release v5.5.0), and Kallisto v0.46 (Bray et al., 2016) to calculate individual transcript abundance. The longest ORFs were used as local BLAST queries against ToxProt and UniProt (the latter limited to apocritans only) with an e-value cutoff of $1 \times 10^{-3}$. For subsequent venom protein identification, we only include transcripts that were identified in our proteome data of dissected and injected venom (Supplementary Table 1). The transcriptome raw data is accessible in NCBI under the BioProject PRJNA733472 with the SRA entry SRR14690757.

## 2.3. Proteomic analysis of crude venom

Prior to shotgun proteomics, the *X. violacea* venom samples (dissected and injected) were denatured, reduced, and alkylated. Each sample (~50 µg) was dissolved in 89 µl 100 mM triethylammonium bicarbonate (TEABC), and cysteine residues were reduced by adding 1 µl 1 M DTT (30 min at 60 °C) and modified by adding 10 µl 0.5 M iodoacetamide (incubation for 30 min in the dark). Then 2 µg trypsin (Promega) were added in 100 mM TEABC and incubated overnight at 30 °C. The resulting peptides were purified and concentrated using OMIX Tips C₁₈ reversed-phase resin (Agilent Technologies, Santa Clara, CA, USA). The peptides were dehydrated in a vacuum centrifuge and analysed by NanoLC-MS/MS. Samples were resuspended in 20 µl buffer A (0.1% formic acid) and 1 µl was loaded onto an analytical 25 cm reversed-phase column (Acclaim Pepmap 100 C₁₈) with a 75 mm inner diameter (Thermo Fisher Scientific) and separated on the Ultimate 3000 RSLC system coupled via a nano-electrospray source to a Q Exactive HF-X mass spectrometer (Thermo Fisher Scientific). Peptides were separated using a 6–40% gradient of buffer B (80% acetonitrile in 0.1% formic acid) over 123 min at a flow rate of 300 nl/min. Using data-dependent acquisition mode, full MS/MS scans (375–1500 m/z) were performed in the Orbitrap mass analyser (Thermo Fisher Scientific) with a 60,000 resolution at 200 m/z. For the full scans, $3 \times 10^6$ ions accumulated within a maximum injection time of 60 ms. The 12 most intense ions with charge states ≥2 were sequentially isolated to a target value of $1 \times 10^5$ with a maximum injection time of 45 ms and were fragmented by higher-energy collisional dissociation in the collision cell (normalized collision energy 28%) and detected in the Orbitrap mass analyser at a resolution of 30,000.

All transcriptome assembly-based predicted ORFs were used as specific databases to identify peptides and proteins detected by mass spectrometry for the two dissected and injected venom samples using PEAKS Studio v8.5 (Bioinformatics Solutions, Waterloo, ON, Canada). Carbamidomethylation was set as a fixed modification, and oxidation of methionine as a variable modification, with a maximum of three missed cleavages for trypsin digestion. Parent and fragment mass error tolerances were set at 5 ppm and 0.015 Da, respectively. A false discovery rate (FDR) of 1% and a unique peptide number ≥2 were used to filter out inaccurate proteins. A -10lgP value > 120 was used to estimate whether detected proteins were identified by a sufficient number of reliable
<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Alluvial chart illustrating the 355 curated hymenopteran venom proteins and peptides. Data was mined in Uniprot using the keywords #Hymenotpera and #venom (December 2021). The number of species per group are given in brackets, the groups are colour-coded (blue = bees and bumblebees, orange = solitary bees, green = ants, red = wasps and hornets, purple = parasitoids and solitary wasps). On the right the protein groups are shown, the number of sequences per group are indicated in brackets, for unique single proteins no number is given. Please note that sequences from melittin-like codesan (from *Colletes daviesanus*) and halictin (from *Halictus sexcinctus*) are not included here. Their sequences are only available as text in the respective manuscripts and are not provided in databases. The species names of solitary bees are underlined and given in brackets.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 2.** The habitus of *X. violacea*. Shown are individuals collecting nectar while pollinating in A.) and B.). The typical habitats are dry and warm biotopes, such as orchard meadows or allotment gardens with sufficient deadwood in range. These pictures were taken in the collection location of allotment gardens in Wieseck, Giessen, Germany. Despite the species being polylactic and pollinating a variety of plants, *X. violacea* seems to favour lip petals (Lamiaceae) in the collection area and was picked from wild pea plants. In the bottom left picture, C.) the strong mandibles are seen which are used to carve the brood tubes into wood. The venom system composed of venom reservoir (also referred to as poison sac), venom duct and the paired venom glands are illustrated in D.).</figcaption></figure>

peptides. In order to identify more relevant sequences, the Spider algorithm (PEAKS Studio) was used to find additional mutations or to correct sequences. This algorithm corrects the sequences stored in transcriptomic databases with *de novo* sequences based on MS/MS spectra, allowing the detection of post-translational modifications (PTMs) and mutations. The minimum ion intensity for PTMs and mutations was set to 5%, and the ALC score was set to ≥90 for *de novo* sequences, leading to low precursor mass errors. Transcripts supported by proteomic data were manually filtered by excluding non-venom-related proteins and peptides, such as house-keeping and structural genes (Supplementary Table 1). The proteomic raw data is submitted to PRIDE with the dataset identifier PXD029823 for dissected venom proteomic data (first published in (Koludarov et al., 2022)) and PXD030997 for injected venom proteome data.

## 2.4. Phylogenetic analyses of major venom components

Sequences of melittin, apamin and MCDP peptides were aligned with the software Mafft (Katoh and Toh, 2008), applying the algorithm mafft-L-INS-I. The rather short and highly similar sequences of the discussed peptides make the use of tree reconstructions difficult, also because substitution models are difficult to fit. We used therefore as alternative neighbour networks that indicate in contrast to strictly bifurcating trees also ambiguous and competing relations of sequences.

The networks were reconstructed using Splitstree 5.3.0 (Huson, 1998; Huson and Bryant, 2006) applying the Protein ML JTT model for underlying distance matrix reconstructions. All other sequences of venom proteins were aligned applying mafft-E-INS-I.

# 3. Results

## 3.1. The general venom profile of *X. violacea*

In total, we matched 199 transcripts with the peptide fragments from the mass spectrometry data of the dissected and injected venom samples. 43 different venom protein groups had TPM values above 1.5 (Fig. 3, sorted from the highest to lowest expression value). Only five of these are composed of short peptides, while other 38 families represent larger peptides or proteins. Interesting, however, is that four of these peptide groups (melittin, apamin, secapin and novel 4) are among the nine most highly expressed components.

By far the most dominant secreted venom components in the dissected and injected venoms are well-known major honeybee toxins (Pucca et al., 2019; Wehbe et al., 2019). They belong to the two peptide families of melittin (TPM 365773.84, ~70%) and apamin (TPM 134017.92, ~25%), and the enzyme family phospholipase A2 (TPM 13212.71, ~2%), see Fig. 1. Other dominant but far less expressed components are three novel protein families 5, 6, and 12 (TPMs of
<figure><img src="image_6.png" /><figcaption>**Fig. 3.** Venom profile of *X. violacea*. For the proteo-transcriptomic data based on dissected venom glands ('gland' proteins) only transcripts were considered that were identified in the proteomic data. Expression values are given in TPM (transcripts per million) and are rounded on the second position after the comma. Transcripts that were also identified in crude venom obtained from agitated stings (injected venom) are marked by a black circle, while proteins that were only identified via dissection of glands are marked in red, see also Supplementary Table 1. Asterisks indicate peptides.</figcaption></figure>

4878.70, 747.75 and 1453.61), icarapin (TPM 2241.09), venom allergen 5 (TPM 1745.63) and secapin (TPM 683.13). All other components of the remaining 34 families are far less abundant with TPM values under 500, except for EF1 alpha (TPM 676.58), which was, however, not identified in the injected venom. Eleven families are represented in the venom from the dissected venom glands only and are not identified in the injected venom (See Fig. 3).

## 3.2. The most dominant venom components are melittin-like peptides

Our data reveals that melittin-like peptides are by far the most expressed venom component. This is not surprising since melittin is known as the most dominant venom component in honeybee venom where it makes 50–60% of dried venom mass (Pucca et al., 2019). The

four melittin-like precursor sequences in our transcriptome feature similar signal and propeptide regions compared to known *Apis* and *Bombus* sequences, see Fig. 4. The mature sequences are very similar to two antimicrobial peptides named xylopin (Xac-1, Xac-2) that previously were identified in crude venom from the Japanese carpenter bee *Xylocopa appendiculata* crimuvolans (Kazuma et al., 2017). Because the short sequences in our alignment make it difficult to fit evolutionary substitution models, we used a distance matrix-based neighbour net reconstruction to visualise phylogenetic relations between known melittin and melittin-like sequences. We generated two separate alignments of available melittin-like sequences, one including the complete precursor sequences available for honeybees, bumblebees and our *X. violacea* transcripts (Fig. 4 B), and one that includes the mature sequences only (Supplementary Fig. 1). Both rooted networks show noise
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 4.** Sequence alignment of available melittin sequences reconstructed with MAFFT-L-INSI (A). Neighbour networks reconstructions based on distance matrices (Protein ML JTT model) are shown for entire sequences; (B). Identifiers in red represent wasp sequences, sequences from *X. violacea* are highlighted in blue. For outgroup rooting we used melittin-like peptides from the European common frog (*Rana temporalis*, P56924) and the moor frog (*Rana arvalis*, P86158).</figcaption></figure>
in the alignments, and alternative possible connections of sequences indicate either ambiguously aligned positions or conflicting signal (for alignments see Supplementary Files 1 and 2). However, both networks reveal distinct groups of which two are the honeybee melittins and the bumblebee melittin-like sequences. Within solitary bee melittin-like sequences we find a grouping of peptides from the solitary bees *Osmia rufa* (Osmin), *Xylocopa violacea/Xylocopa appendiculata circumvolans* (Xylopin), *Macropis fulvipes* (Macropin), *Melecta albifrons* (Melectin), *Colletes daviesanus* (Codesane), *Panurgus calcaratus* (Panurgin1) and *Halictus sexcinctus* (Halictin) for which a relation to bombolittins or melittins remains ambiguous. The melittin-like peptides from *Lasioglossa laticens* (Lasioglossin) appear to be closer related to bombolittin sequences than to other solitary bee peptides. Our analysis also indicates that the melittin-like sequences published for *Polistes*, *Vespa* and *Vespula* wasp species probably are contaminations or mistaken samples being identical to bee melittin sequences, see Fig. 4.

### 3.3. The second dominant venom component are apamin-like peptides

The second highly expressed component in *X. violacea* are peptides that are very similar to apamin, a second-highest expressed toxin in honeybee venom (Pucca et al., 2019), and MCDP, which is a peptide that resembles apamin very closely. In our upcoming study (Koludarov et al., 2022) using genomic data we show that apamin and MCDP are not only closely related to each other but also belong to one large peptide family, anthophilin1, which is restricted to bees. Interestingly, the syntenic information of the genomic regions indicate that MCDP is only occurring in *Apis*. Characteristically, apamin and MCDP feature two disulphide bonds that fold the mature peptide by connecting the 1st with the 3rd and the 2nd with the 4th cysteine residue, see Fig. 5 A. The 13 apamin-like transcripts that we identify for *X. violacea*, appear to be more diverse and split into three separate clades with six transcripts for clade one and clade two, which encode the same peptide sequences. The third clade is represented by a single transcript that lacks the 3rd and 4th cysteine but shows an identical signal and propeptide seen in clade two Fig. 5 A.

To investigate the relationships between the apamin clades in *X. violacea* we again used neighbour network reconstructions of the complete and mature transcript sequences including available sequences of apamin and highly similar mast cell degranulating peptides (see alignments given in Supplementary Files 3 and 4). The network based on the full sequence rather separates the species clades (*Apis*, *Bombus* and *Xylocopa*), see Fig. 5 B. It appears that three groups of apamin-like sequences can be recognized.

The network of the mature sequences (Fig. 5 C) also distinguishes three groups of apamin-like sequences in *X. violacea*. Clade two clusters with MCDP from *Apis*, while peptides from clade one are more similar to the MCDP sequence of *Bombus* (which probably is an Apamin). The single transcript in clade three seems more closely related to apamin from *Apis*. However, it becomes apparent that the resolution of the short sequences is too low and ambiguous relationships are apparent that cannot be resolved based on this dataset.

### 3.4. Dominant larger venom peptides and proteins

Several larger proteins were identified in the venom which are, however, far less abundant compared to melittin and apamin (Fig. 3). Phospholipase A2 is with six transcripts (TPM 13212.71) the third most expressed venom component in *X. violacea* (see Supplementary File 5 for its alignment). In honeybee and bumblebee venoms phospholipase A2 is well known as a highly expressed major toxic component that is enhanced in its effect to attack membrane phospholipids by the pore forming melittin (Pucca et al., 2019).

Icarapin is another major allergen known from honeybees (also referred to as Api m 10). So far it has only been described directly from the venom of *Apis* species mostly in the context of bee venom treatment.

However, matches from genomic sequence data to other solitary bees, wasps and ants suggest that this protein is also present in other bees and hymenopterans, which is confirmed here for the solitary bee *X. violacea* (See Supplementary File 6 for aligned sequences).

Venom allergens are a large protein family of which representatives occur in diverse venoms of all kinds of venomous animals (CAP super-family). Venom allergens 5 sequences that we identified in *X. violacea* are similar to variants from wasps and ants, but also to sequences identified in genomic data from *Apis* and *Bombus* venom (See supplementary File 7 for aligned sequences).

Secapin was described in *Apis* and positively tested for various bio-activities, predominantly antimicrobial or antifibrinolytic (Pucca et al., 2019). Only available non-honeybee crude venom evidence comes from ants and suggests a higher similarity between ant and bee secapins. The two transcripts identified in *X. violacea* are more distinct from *Apis* secapin and also have a different mature peptide, with four instead of two Cysteine residues. We have to note that similar to the picture in melittin, the secapin sequences that derive from wasps are close to be identical to bee sequences (mostly 100%), which again indicates a possible contamination or erroneous database upload (See Supplementary File 8).

### 3.5. Novel venom components

Finally, our analysis reveals eight novel venom components for which only vague annotations can be given. For that reason, we refrain here from a more detailed discussion and refer to Table 1 and Supplementary Table 1. Nevertheless, it is striking that the cleaved products from these novel precursor sequences are predominantly peptides, of which several show Cysteine scaffolds (Table 1). This always hints towards a peptide folding based on disulphide bridges which might indicate interesting bioactivities for applied research. Except for number two and four, all novel proteins are identified in the injected venom as well.

## 4. Discussion

### 4.1. *X. violacea* has a complex venom profile with several unique components

The honeybee (*Apis*) is among the very few hymenopteran taxa for which the venom has been known in a great detail for quite a long time (Piek, 1986; Pucca et al., 2019). This predominant focus on the honeybee facilitates even the misleading use of the term 'bee venom' by referring only to the honeybee venoms instead of the whole clade Anthophila (bees *sensu lato*). It is known by now that honeybee venom composition and production can alter during the year in worker bees and that queens have a different venom composition (Danneels et al., 2015; Scaccabarozzi et al., 2021). A recent study on eusocial bumblebees (*Bombus*), the sister group to *Apis*, showcases further that venom compositions are dynamic and vary according to the habitat altitude (Barkan et al., 2020). To understand these obviously highly dynamic adaptations in bee venoms and their general mechanisms in more detail a broader understanding of bee venoms *sensu lato* is necessary, as well as work comparing honeybee and bumblebee venoms to other evolutionary older solitary (non-eusocial) species of bees. The herein discussed venom composition of *X. violacea* with 32 proteinaceous venom components identified in its injected venom and 43 components in the dissected venom system illustrates that solitary bees might have a very complex venom composition. The components which are absent in the injected venom are likely chaperone proteins and enzymes (such as dismutase, carboxypeptidase) that facilitate the maturation of the venom components. To better understand the complexity and possible adaptations of envenomation processes and immunology in bees (*sensu lato*), more data on injected venom components (which interact synergistically) are important from different bee lineages. Most injected
<figure><img src="image_9.png" /><figcaption>B.)</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>C.)</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>Fig. 5. Sequence alignment of available apamin and MCDP sequences reconstructed with MAFFT-L-INSI (A). Neighbour networks reconstructed based on distance matrices (Protein ML JTT model) are shown for entire sequences (B) and mature sequences (C). Sequences from *X. violacea* are highlighted in blue.</figcaption></figure>
<caption>**Table 1**</caption>

Overview of novel venom components identified in *X. violacea*. See Supplementary Table 1 for further details on annotations. Components that are present in the injected venom are indicated by an asterisk. SP column shows the length of predicted signal peptides in amino acid positions. All sequences are given in Supplementary Files 9 to 16.

<table><thead><tr><th>Name</th><th>Transcript name</th><th>Length</th><th>SP</th><th>Pattern</th></tr></thead><tbody><tr><td>Novel 1*</td><td>R662262.p2</td><td>94 aa</td><td>15</td><td>X (3)-C-X (45)-C-X (5)-C-X (5)-C-X (11)</td></tr><tr><td>Novel 1*</td><td>NODE17211_g13121i0. p1</td><td>83 aa</td><td>15</td><td>X (3)-C-X (45)-C-X (5)-C-X (5)-C- (truncated)</td></tr><tr><td>Novel 2</td><td>DN7443_c0g3i1.p4</td><td>53 aa</td><td>NA</td><td>Leucin, Lysn, Isoleucin rich</td></tr><tr><td>Novel 3*</td><td>NODE12433_g1692_i3. p4</td><td>54 aa</td><td>26</td><td>X (2)-C-X-C-X (2)-PPP-X (2)-P-X (5)-C-X (3)-C-X (4)</td></tr><tr><td>Novel 3*</td><td>NODE4379_g3275i0.p6</td><td>54 aa</td><td>26</td><td>X (2)-C-X-C-X (2)-PPP-X (2)-P-X (5)-C-X (3)-C-X (4)</td></tr><tr><td>Novel 3*</td><td>NODE7033_g159i1.p3</td><td>54 aa</td><td>26</td><td>X (2)-C-X-C-X (2)-PPP-X (2)-P-X (5)-C-X (3)-C-X (4)</td></tr><tr><td>Novel 4</td><td>NODE26737_g22473i0. p1</td><td>118 aa</td><td>NA</td><td>Serine, Isoleucin, Asparagin, Isoleucin rich</td></tr><tr><td>Novel 5*</td><td>DN2824_c0g1i3.p3</td><td>71 aa</td><td>42</td><td>X (4)-C-X (7)-C-X (8)-C-X (3)-C-X (3)</td></tr><tr><td>Novel 6*</td><td>R662790.p8</td><td>42 aa</td><td>NA</td><td>Isoleucin, Lysine, Valine, Arginine rich</td></tr><tr><td>Novel 6*</td><td>DN49_c0g1i1.p7</td><td>42 aa</td><td>NA</td><td>Isoleucin, Lysine, Valine, Arginine rich</td></tr><tr><td>Novel 6*</td><td>DN49_c0g1i9.p7</td><td>42 aa</td><td>NA</td><td>Isoleucin, Lysine, Valine, Arginine rich</td></tr><tr><td>Novel 7*</td><td>J678067.p1</td><td>445 aa</td><td>NA</td><td>NA</td></tr><tr><td>Novel 8*</td><td>DN10104_c0g1i1.p1</td><td>165 aa</td><td>NA</td><td>Contains predicted disorder-regions</td></tr></tbody></table>

venom proteins (Fig. 3) we identified in *X. violacea* belong to groups shared with honeybees and bumblebees (Barkan et al., 2017, 2020; Pucca et al., 2019). However, it is striking that novel peptide classes are found in the 9 most highly expressed components of the carpenter bee venom, for which activity and function remain to be studied. The question if these occur also in venoms of other solitary bees or if they are possibly unique to carpenter bees cannot be addressed without further studies of other carpenter bee species' and other solitary bee species' venoms.

## 4.2. Insights into the evolution of bee venom components

Melittin that was first described in honeybees represents their most dominant venom component and is probably the singular best studied Hymenopteran toxin (Pucca et al., 2019). Melittin-like variants (Bombolitin) are also expressed in venom of Bombus (Argiolas and Pisano, 1985), despite bumblebee venom being dominated by phospholipase A2 in contrast to honeybee venom (Barkan et al., 2017, 2020; Yoon et al., 2020). It was illustrated in recent proteomic studies that most venom components from *Apis* and *Bombus* are generally very similar (Barkan et al., 2017; Van Vaerenbergh et al., 2015), which is probably linked to the close phylogenetic relationship of these two genera (Peters et al., 2017), see also Fig. 6. We show that in *X. violacea* most venom components are also similar to homologs known from *Apis* and that melittin is the most dominant venom component. The network reconstruction (Fig. 4) indicates that most melittin-like peptides from solitary bees form a single clade to the exclusion of their bumblebee and honeybee counterparts. In general, melittins from *Bombus* and *Apis* are more distinct from the ones found in solitary bees and only lasioglossins from *Lasioglossa* bees appear closer related to bombolitins. However, our networks clearly indicate that short but partially more distinct peptides are of limited use for the analysis of phylogenetic relationships due to the higher impact of randomly similar aligned amino acids which subsequently can impede appropriate model fitting. As a consequence, to avoid phylogenetic conclusions based on similarity, alternative methods

become important. In our recent genomic analysis, we revealed that melittin genes are found in several available genomes of solitary bee clades including the carpenter bees (Koludarov et al., 2022), see also Fig. 6. Using micro-syntenic patterns of gene regions that flank the exons of melittin in *Apis* we show that this toxin gene is also present in the carpenter bee *C. calcarata* and in the megachilid genus *Osmia*. Searching available wasp genomes, we did not find any matches (neither primary sequence-based nor microsynteny-based) which highly supports our deduction that the sequences published as wasp-melittins must be erroneously assigned (see Fig. 4). From the current data we can deduce that melittin probably evolved already in the earlier bee lineages. However, most of the peptides that we included in our analyses are highly variable and short, which is reflected in the ambiguous relations of sequences that are revealed in the networks (in contrast to strictly bifurcating trees). To prevent false interpretation by noisy alignments and subsequent misled phylogenetic analyses, corresponding genomic data is important that allows to identify the coding genes for these peptides and to reconstruct phylogenies based on the full coding domain sequences and microsyntenic patterns.

Apamin and MCDP are two highly similar toxic peptides originally identified in *Apis* venom that both belong to the recently described toxin family Anthophilin1 (Koludarov et al., 2022). In our network analysis we identified three different forms of apamin-like sequences in *X. violacea*, for which it is hard to derive definite phylogenetic relationships. Interestingly, we showed previously that the gene for MCDP appears unique to *Apis* and was not identified in any other bee clade based on genomic microsynteny analysis. As a consequence, a closer similarity of the apamin-like peptides from *X. violacea* to the *Bombus* sequence (that is from our perspective falsely annotated as MCDP) would reflect instead another distinct apamin variant of which *X. violacea* might have three. To finally test if MCDP is indeed not present in *X. violacea* and supporting its uniqueness to *Apis*, a high-quality genome of *X. violacea* would be needed. Such a genome would be of course important to draw conclusions about the general bee venom evolution as well.

Most other major components in the venom of *X. violacea*, such as Icarapin, Secapin, Phospholipase A2 and Venom allergens (CAP) are also present in other bee venoms and evolved in the early lineage of hymenopterans (Fig. 6). This is also the case for many of the lowly expressed venom components. From a biological perspective the question remains why a predominantly defensive venom of bees (including the herein described venom of *X. violacea*) shows such a complex composition. To address this further, future studies would require more holistic overview of bee venoms from a variety of anthophilan species combining proteo-transcriptomic with genomic data. Likewise it would be important to analyse in a greater detail seasonal or habitual variations of bee venoms, advancing further the existing data (Barkan et al., 2020; Danneels et al., 2015).

## 4.3. Functional and applied aspects of identified venom components

Surprisingly, not only the reason for the complexity of a defensive bee in bees is still not fully understood, but also the biological activities and possible synergistic effects of most of the venom components in bees remain unknown. In general, most studies focus on possible allergenic activity of components in the context of envenomation and possible mitigation if humans are stung, see e.g. (Pucca et al., 2019). Another strong focus is on the effects of short, melittin-like peptides that exhibit higher potential for a variety of translational research objectives, such as promising antimicrobial, anti-tumor and anti-inflammatory activities (Cujová et al., 2014; Monincová et al., 2010, 2010, 2012, 2014, 2010; Slaninova et al., 2012; Stöcklin et al., 2010; Wehbe et al., 2019). Nevertheless, at least for the two major components in honeybee venom, melittin and phospholipase A2, the general mode of action has been described. The various bioactivities of melittin are predominantly based on its pore forming activity which synergistically enhances the effect of
<figure><img src="image_12.png" /><figcaption>**Fig. 6.** Most expressed venom components of *X. violacea* in the context of bee phylogeny and available proteomic data of bees. Red “X” marks indicate missing or reduced venom system from stingless bees (Meliponini). Phylogeny and divergence times are according to (Peters et al., 2017; Sann et al., 2018).</figcaption></figure>
phospholipase A2 which disrupts cell membranes by attacking phospholipids that are exposed by melittin, see e.g. (Pucca et al., 2019). One difficulty with melittin is that it is so destructive that its application potential has proved challenging to tame for translational research and to modify and to apply it in a way that its effects are more targeted (Gajski et al., 2016). Given that melittin likely evolved earlier within solitary bees, the phylogenetically older melittin variants discovered in this study from *X. violacea* venom are interesting because their bioactivities could differ or be less strong and more applicable compared to melittin tested from *Apis*.

## 5. Conclusions

In the present study we demonstrate that solitary bee venoms are as complex as their more known honeybee and bumblebee counterparts and deserve more attention from the research community. Despite superficial similarities in terms of protein classes, bee venoms are dynamic systems that continuously tinker with their molecular weaponry. The present study highlights an untapped biodiscovery potential of bee venoms by revealing entirely novel venom peptide classes as well as new forms of well-known medicinally important molecules like melittin and apamin as part of the new anthophilin family. Our study also illustrates the importance of more dedicated bee venom research, in particular bioactivity, genomic and proteo-transcriptomic studies of neglected solitary bee species to better understand the evolution and adaptations of their biochemical defence system.

### Ethical statement

The authors adhered to Elsevier's publishing ethics policy and ethical guidelines for journal publication.

### Author contributions

BMvR, SD, IK: conceptualization, BMvR: Proteo-transcriptomic/Phylogenetic analyses, Visualization, Writing -original draft. SD: Proteomic analyses, methodology. IK, SD, BMvR: Writing – review and editing, methodology, BMvR: Project administration.

### Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

### Acknowledgement

BMvR and IK acknowledge funding on aculeate venom research by a grant of the German Science Foundation to BMvR (DFG RE3454/6–1). IK is supported by Alexander von Humboldt Forschungsstipendium. BMvR and IK thank Andreas Vilcinskas for providing work space and logistics via the group Animal Venomics at the Institute for Insectobiotechnology at the Justus Liebig University Giessen. Animal Venomics is funded via the Centre for Translational Biodiversity Genomics (LOEWE-TBG) in the programme "LOEWE – Landes-Offensive zur Entwicklung Wissenschaftlich-ökonomischer Exzellenz" of Hesse’s Ministry of Higher Education, Research, and the Arts. BMvR thanks further Ingo Ebersberger for his support. SD acknowledges the Functional Proteomics Platform of Montpellier, where high resolution mass spectrometry experiments were carried out.

## Appendix A. Supplementary data

Supplementary data to this article can be found online at https://doi.org/10.1016/j.toxcx.2022.100117.

## References

Argiolas, A., Pisano, J.J., 1985. Bombolitins, a new class of mast cell degranulating peptides from the venom of the bumblebee *Megabombus pennsylvanicus*. *J. Biol. Chem.* 260, 1437–1444.

Banaszak, J., Banaszak-Cibicka, W., Twerd, L., 2019. Possible expansion of the range of *Xylocopa violacea* L. (Hymenoptera, apiformes, apidae) in europe. *Turk. J. Zool.* 43, 650–656.

Barkan, N.P., Bayazit, M.B., Ozel Demiralp, D., 2017. Proteomic characterization of the venom of five Bombus (*thoracobombus*) species. *Toxins* 9, 362. https://doi.org/10.3390/toxins9110362.

Barkan, N.P., Chevalier, M., Pradervand, J.-N., Guisan, A., 2020. Alteration of bumblebee venom composition toward higher elevation. *Toxins* 12, 4. https://doi.org/10.3390/toxins12010004.

Bert, P., 1865. Contributions à l'otude des venins. II. Venin d'abeille xylocope. Gaz. Med. 20, 770–771.

Bouzid, W., Klopp, C., Verdenaud, M., Ducancel, F., Vétillard, A., 2013. Profiling the venom gland transcriptome of *Tetramorium bicarinatum* (Hymenoptera: formicidae): the first transcriptome analysis of an ant species. *Toxicon* 70, 70–81. https://doi.org/10.1016/j.toxicon.2013.03.010.

Branstetter, M.G., Danforth, B.N., Pitts, J.P., Faircloth, B.C., Ward, P.S., Buffington, M.L., Gates, M.W., Kula, R.R., Brady, S.G., 2017. Phylogenomic insights into the evolution of stinging wasps and the origins of ants and bees. *Curr. Biol. CB* 27, 1019–1025. https://doi.org/10.1016/j.cub.2017.03.027.

Bray, N.L., Pimentel, H., Melsted, P., Pachter, L., 2016. Near-optimial probabilistic RNA-seq quantification. *Nat. Biotechnol.* 34, 525–527. https://doi.org/10.1038/nbt.3519.

Cerpes, U., Repelnig, M.-L., Legat, F.J., 2021. Itch in Hymenoptera sting reactions. Front. Allergy 2, 46. https://doi.org/10.3389/falgy.2021.727776.

Cujová, S., Bednárová, L., Slaninová, J., Straka, J., CEROVSKÝ, V., 2014. Interaction of a novel antimicrobial peptide isolated from the venom of solitary bee *Colletes daviesanus* with phospholipid vesicles and *Escherichia coli* cells. J. Pept. Sci. 20, 885–895. https://doi.org/10.1002/psc.2681.

Danneels, E., Van Vaerenbergh, M., Debyser, G., Devreese, B., de Graaf, D., 2015. Honeybee venom proteome profile of queens and winter bees as determined by a mass spectrometric approach. *Toxins* 7, 4468–4483. https://doi.org/10.3390/toxins7114468.

de Graaf, D.C., Aerts, M., Danneels, E., Devreese, B., 2009. Bee, wasp and ant venomics pave the way for a component-resolved diagnosis of sting allergy. J. Proteomics 72, 145–154. https://doi.org/10.1016/j.jprot.2009.01.017.

Gajski, G., Domijan, A.-M., Zegura, B., Stern, A., Gerić, M., Novak Jovanović, I., Vrhovac, I., Madunić, J., Breljak, D., Filipić, M., Garaj-Vrhovac, V., 2016. Melittin induced cytogenetic damage, oxidative stress and changes in gene expression in human peripheral blood lymphocytes. *Toxicon* 110, 56–67. https://doi.org/10.1016/j.toxicon.2015.12.005.

Haas, B.Jo, Release v5.5.0. TransDecoder [WWW Document]. TransDecoder. URL https://github.com/TransDecoder/TransDecoder.

Hardouin, R., 1948. La vie des Abeilles Solitaires. Libraire Gallimard, Paris.

Hermann, H., Mullen, M., 1974. The hymenopterous poison apparatus. XI. *Xylocopa virgínica* (Hymenoptera: Xylocopidae. J. Ga. Entomol. Soc. 9, 246–252.

Huson, D.H., 1998. *Splitstree: analyzing and visualizing evolutionary data*. Bioinformatics 14, 68–73.

Huson, D.H., Bryant, D., 2006. Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol. 23, 254–267. https://doi.org/10.1093/molbev/msj030.

Jungo, F., Bougueleret, L., Xenarios, I., Poux, S., 2012. The UniProtKB/Swiss-Prot Tox-Prot program: a central hub of integrated venom protein data. *Toxicon* 60, 551–557. https://doi.org/10.1016/j.toxicon.2012.03.010.

Katoh, K., Toh, H., 2008. Recent developments in the MAFFT multiple sequence alignment program. Briefings Bioinf. 9, 286–298. https://doi.org/10.1093/bib/bbn013.

Kazuma, K., Ando, K., Nihei, K., Wang, X., Rangel, M., Franzolin, M.R., Mori-Yasumoto, K., Sekita, S., Kadowaki, M., Satake, M., Konno, K., 2017. Peptidomic analysis of the venom of the solitary bee *Xylocopa appendiculata circumvolans*. J. Venom. Anim. Toxins Incl. Trop. Dis. 23, 40. https://doi.org/10.1186/s40409-017-0130-y.

Klein, A.-M., Boreux, V., Fornoff, F., Mupepele, A.-C., Pufal, G., 2018. Relevance of wild and managed bees for human well-being. Curr. Opin. Insect Sci. Ecol. Parasite. Parasitoid. Biol. Contr. 26, 82–88. https://doi.org/10.1016/j.cois.2018.02.011.

Ko, S.J., Park, E., Asandei, A., Choi, J.-Y., Lee, S.-C., Seo, C.H., Luchian, T., Park, Y., 2020. Bee venom-derived antimicrobial peptide melectin has broad-spectrum potency, cell selectivity, and salt-resistant properties. Sci. Rep. 10, 10145. https://doi.org/10.1038/s41598-020-06695-7.

Koh, I., Lonsdorf, E.V., Williams, N.M., Brittain, C., Isaacs, R., Gibbs, J., Ricketts, T.H., 2016. Modeling the status, trends, and impacts of wild bee abundance in the United States. Proc. Natl. Acad. Sci. Unit. States Am. 113, 140–145. https://doi.org/10.1073/pnas.1517685113.

Koludarov, I., Velasque, M., Timm, T., Lochnit, G., Heinzinger, M., Vilcinskas, A., Gloag, R., Harpur, B.A., Podsiadlowski, L., Rost, B., Dutertre, S., Stolle, E., Reumont, B.M. von, 2022. Bee Core Venom Genes Predominantly Originated in Ancestors of Sawfly Hymenopterans before Aculeate Stingers Evolved. bioRxiv. https://doi.org/10.1101/2022.01.21.477203.

Kukuk, P.F., Eickwort, G.C., Raveret-Richter, M., Alexander, B., Gibson, R., Morse, R.A., Ratnieks, F., 1989. Importance of the sting in the evolution of sociality in the hymenoptera. Ann. Entomol. Soc. Am. 82, 1–5. https://doi.org/10.1093/aesa/82.1.1.
Kularatne, S.A.M., Raveendran, S., Ediringishe, J., Karunaratne, I., Weerakoon, K., 2016. First reported case of fatal stinging by the large carpenter bee *Xylocopa tranquebarica*. *Wilderness Environ. Med.* 27, 262–265. https://doi.org/10.1016/j.wem.2015.12.018.

Leluk, J., Schmidt, J., Jones, D., 1989. Comparative studies on the protein composition of hymenopteran venom reservoirs. *Toxicon* 27, 105–114.

Liu, Z., Chen, S., Zhou, Y., Xie, C., Zhu, B., Zhu, H., Liu, S., Wang, W., Chen, H., Ji, Y., 2015. Deciphering the venomic transcriptome of killer-wasp *Vespa velutina*. *Sci. Rep.* 5, 9454. https://doi.org/10.1038/srep09454.

MacManes, M.D., 2018. The Oyster River Protocol: a multi-assembler and kmer approach for *de novo* transcriptome assembly. *PeerJ* 6, e5428. https://doi.org/10.7717/peerj.5428.

Monincová, L., Budešínský, M., Slaninová, J., Hovorka, O., Cvačka, J., Voburka, Z., Fučík, V., Borovičková, L., Bednárová, L., Straka, J., Cerovský, V., 2010. Novel antimicrobial peptides from the venom of the eusocial bee *Halictus sexcinctus* (Hymenoptera: halictidae) and their analogs. *Amino Acids* 39, 763–775. https://doi.org/10.1007/s00726-010-0519-1.

Monincová, L., Slaninová, J., Fučík, V., Hovorka, O., Voburka, Z., Bednárová, L., Malon, P., Stokrová, J., Cerovský, V., 2012. Lasiocepsin, a novel cyclic antimicrobial peptide from the venom of eusocial bee *Lasioglossum laticeps* (Hymenoptera: halictidae). *Amino Acids* 43, 751–761. https://doi.org/10.1007/s00726-011-1125-6.

Monincová, L., Veverka, V., Slaninová, J., Budešínský, M., Fučík, V., Bednárová, L., Straka, J., Cerovský, V., 2014. Structure-activity study of macropin, a novel antimicrobial peptide from the venom of solitary bee *Macropis fulvipes* (Hymenoptera: melittidae). *J. Pept. Sci.* 20, 375–384. https://doi.org/10.1002/psc.2625.

Moreau, S.J., 2013. It stings a bit but it cleans well": venoms of Hymenoptera and their antimicrobial potential. *J. Insect Physiol.* 59, 186–204. https://doi.org/10.1016/j.jinsphys.2012.10.005.

Oeyen, J.P., Baa-Puyoulet, P., Benoit, J.B., Beukeboom, L.W., Bornberg-Bauer, E., Buttstedt, A., Calevro, F., Cash, E.I., Chao, H., Charles, H., Chen, M.-J.M., Childers, C., Cridge, A.G., Dearden, P., Dinh, H., Doddapaneni, H.V., Dolan, A., Donath, A., Dowling, D., Dugan, S., Duncan, E., Elpidina, E.N., Friedrich, M., Geuverink, E., Gibson, J.D., Grath, S., Grimmelikhuijzen, C.J.P., Große-Wilde, E., Gudobba, C., Han, Y., Hansson, B.S., Hauser, F., Hughes, D.S.T., Ioannidis, P., Jacquín-Joly, E., Jennings, E.C., Jones, J.W., Klasberg, S., Lee, S.L., Lesný, P., Lovegrove, M., Martin, S., Martynov, A.G., Mayer, C., Montagné, N., Moris, V.C., Munoz-Torres, M., Murali, S.C., Muzny, D.M., Oppert, B., Parisot, N., Pauli, T., Peters, R.S., Petersen, M., Pick, C., Persyn, E., Podsiadlowski, L., Poelchau, M.F., Provataris, P., Qu, J., Reijnders, M.J.M.F., von Reumont, B.M., Rosendale, A.J., Simao, F.A., Skelly, J., Sotiriopoulos, A.G., Stahl, A.L., Sumitani, M., Szuter, E.M., Tidswell, O., Tsitlakidis, E., Vedder, L., Waterhouse, R.M., Werren, J.H., Wilbrandt, J., Worley, K.C., Yamamoto, D.S., van de Zande, L., Zdobnov, E.M., Ziesmann, T., Gibbs, R.A., Richards, S., Hatakeyama, M., Misof, B., Niehuis, O., 2020. Sawfly genomes reveal evolutionary acquisitions that fostered the mega-radiation of parasitoid and eusocial Hymenoptera. *Genome Biol. Evol.* 12, 1099–1188. https://doi.org/10.1093/gbe/evaa106.

Ozbek, R., Wielsch, N., Vogel, H., Lochnit, G., Foerster, F., Vilcinskas, A., von Reumont, B.M., 2019. Proteo-Transcriptomic characterization of the venom from the endoparasitoid wasp *Pimpla turionellae* with aspects on its biology and evolution. *Toxins* 11, 721. https://doi.org/10.1038/nprot.2006.468.

Peters, R.S., Krogmann, L., Mayer, C., Donath, A., Gunkel, S., Meusemann, K., Kozlov, A., Podsiadlowski, L., Petersen, M., Lanfear, R., Diez, P.A., Heraty, J., Kjer, K.M., Klopstein, S., Meier, R., Polidori, C., Schmitt, T., Liu, S., Zhou, X., Wappler, T., Rust, J., Misof, B., Niehuis, O., 2017. Evolutionary history of the hymenoptera. *Curr. Biol.* CB 27, 1013–1018. https://doi.org/10.1016/j.cub.2017.01.027.

Piek, T., 1986. Venoms of the Hymenoptera, Biochemical, Pharmacological and Behavioural Aspects. Elsevier.

Pucca, M.B., Cerni, F.A., Oliveira, I.S., Jenkins, T.P., Argemí, L., Sørensen, C.V., Ahmadi, S., Barbosa, J.E., Laustsen, A.H., 2019. Bee updated: current knowledge on bee venom and bee envenoming therapy. *Front. Immunol.* 10, 1–15. https://doi.org/10.3389/fimmu.2019.02090.

Raghuraman, H., Chattopadhyay, A., 2007. Melittin: a membrane-active peptide with diverse functions. *Biosci. Rep.* 27, 189–223. https://doi.org/10.1007/s10540-006-9030-z.

Robinson, S.D., Mueller, A., Clayton, D., Starobova, H., Hamilton, B.R., Payne, R.J., Vetter, I., King, G.F., Undheim, E.A.B., 2018. A comprehensive portrait of the venom of the giant red bull ant, *Myrmecia gulosa*, reveals a hyperdiverse hymenopteran toxin gene family. *Sci. Adv.* 4, eaau4640. https://doi.org/10.1126/sciadv.aau4640.

Sann, M., Niehuis, O., Peters, R.S., Mayer, C., Kozlov, A., Podsiadlowski, L., Bank, S., Meusemann, K., Misof, B., Bleidorn, C., Ohl, M., 2018. Phylogenomic analysis of Apoidea sheds new light on the sister group of bees. *BMC Evol. Biol.* 18, 71. https://doi.org/10.1186/s12862-018-1155-8.

Scaccabarozzi, D., Dods, K., Le, T.T., Gummer, J.P.A., Lussu, M., Milne, L., Campbell, T., Wafujian, B.P., Priddis, C., 2021. Factors driving the compositional diversity of *Apis mellifera* bee venom from a *Corymbia calophylla* (marri) ecosystem, Southwestern Australia. *PLoS One* 16, e0253838. https://doi.org/10.1371/journal.pone.0253838.

Schmidt, J.O., 1982. Biochemistry of insect venoms. *Annu. Rev. Entomol.* 27, 339–368. https://doi.org/10.1146/annurev.en.27.010182.002011.

Schmidt, J.O., Blum, M.S., Overall, W.L., 1986. Comparative enzymology of venoms from stinging hymenoptera. *Toxicon* 24, 907–921.

Slaninova, J., Mlsova, V., Kroupova, H., Alan, L., Tumova, T., Monincova, L., Borovičková, L., Fucík, V., Cerovský, V., 2012. Toxicity study of antimicrobial peptides from wild bee venom and their analogs toward mammalian normal and cancer cells. *Peptides* 33, 18–26. https://doi.org/10.1016/j.peptides.2011.11.002.

Stöcklin, R., Favreau, P., Thai, R., Pflugfelder, J., Bulet, P., Mebs, D., 2010. Structural identification by mass spectrometry of a novel antimicrobial peptide from the venom of the solitary bee *Osmia rufa* (Hymenoptera: megachilidae). *Toxicon* 55, 20–27. https://doi.org/10.1016/j.toxicon.2008.12.011.

Tan, J., Wang, W., Wu, F., Li, Y., Fan, Q., 2020. Transcriptome profiling of venom gland from wasp species: de novo assembly, functional annotation, and discovery of molecular markers. *BMC Genom.* 21, 427. https://doi.org/10.1186/s12864-020-06851-0.

Teng, Z.-W., Xiong, S.-J., Xu, G., Gan, S.-Y., Chen, X., Stanley, D., Yan, Z.-C., Ye, G.-Y., Fang, Q., 2017. Protein discovery: combined transcriptomic and proteomic analyses of venom from the endoparasitoid *Cotesia chilonis* (Hymenoptera: braconidae). *Toxins* 9, 1–24. https://doi.org/10.3390/toxins9040135.

Touchard, A., Dejean, A., Escoubas, P., Orivel, J., 2015. Intraspecific variations in the venom peptidome of the ant *Odontomachus haematodus* (Formicidae: ponerinae) from French Guiana. *J. Hymenoptera Res.* 47, 87–101. https://doi.org/10.3897/JHR.47.6804.

Touchard, A., Téné, N., Song, P.C.T., Lefranc, B., Leprince, J., Treilhou, M., Bonnafé, E., 2018. Deciphering the molecular diversity of an ant venom peptidome through a venomics approach. *J. Proteome Res.* 17, 3503–3516. https://doi.org/10.1021/acs.jproteome.8b00452.

Van Vaerenbergh, M., Debyser, G., Smaghe, G., Devreese, B., de Graaf, D.C., 2015. Unraveling the venom proteome of the bumblebee (*Bombus terrestris*) by integrating a combinatorial peptide ligand library approach with FT-ICR MS. *Toxicon* 102, 81–88. https://doi.org/10.1016/j.toxicon.2013.10.002.

Walker, A.A., Robinson, S.D., Yeates, D.K., Jin, J., Baumann, K., Dobson, J., Fry, B.G., King, G.F., 2018. Entomo-venomics: the evolution, biology and biochemistry of insect venoms. *Toxicon* 154, 15–27. https://doi.org/10.1016/j.toxicon.2018.09.004.

Wang, T., Zhao, M., Rotgans, B.A., Ni, G., Dean, J.F.D., Nahrung, H.F., Cummins, S.F., 2016. Proteomic analysis of the venom and venom sac of the woodwasp, *Sirex noctilio* - towards understanding its biological impact. *J. Proteomics* 146, 195–206. https://doi.org/10.1016/j.jprot.2016.07.002.

Wehbe, R., Frangieh, J., Rima, M., El Obeid, D., Sabatier, J.-M., Fajloun, Z., 2019. Bee Venom: overview of main compounds and bioactivities for therapeutic interests. *Molecules* 24, 2997. https://doi.org/10.3390/molecules24162997.

Yoon, K.A., Kim, K., Kim, W.-J., Bang, W.Y., Ahn, N.-H., Bae, C.-H., Yeo, J.-H., Lee, S.H., 2020. Characterization of venom components and their phylogenetic properties in some aculeate bumblebees and wasps. *Toxins* 12, 47. https://doi.org/10.1186/s40168-015-0095-4.