Transcriptome profiling of venom gland from wasp species: De novo assembly, functional annotation, and discovery of molecular markers

Abstract

Background: Vespa velutina, one of the most aggressive and fearful wasps in China, can cause grievous allergies and toxic reactions, leading to organ failure and even death. However, there is little evidence on molecular data regarding wasps. Therefore, we aimed to provide an insight into the transcripts expressed in the venom gland of wasps. Results: In our study, high-throughput RNA sequencing was performed using the venom glands of four wasp species. First, the mitochondrial cytochrome C oxidase submit I (COI) barcoding and the neighbor joining (NJ) tree were used to validate the unique identity and lineage of each individual species. After sequencing, a total of 127,630 contigs were generated and 98,716 coding domain sequences (CDS) were predicted from the four species. The Gene ontology (GO) enrichment analysis of unigenes revealed their functional role in important biological processes (BP), molecular functions (MF) and cellular components (CC). In addition, c-type, p1 type, p2 type and p3 type were the most commonly found simple sequence repeat (SSR) types in the four species of wasp transcriptome. There were differences in the distribution of SSRs and single nucleotide polymorphisms (SNPs) among the four wasp species. Conclusions: The transcriptome data generated in this study will improve our understanding on bioactive proteins and venom-related genes in wasp venom gland and provide a basis for pests control and other applications. To our knowledge, this is the first study on the identification of large-scale genomic data and the discovery of microsatellite markers from V. tropica ducalis and V. analis fabricius.