Transcriptome analysis of an endoparasitoid wasp Cotesia chilonis (Hymenoptera: Braconidae) reveals genes involved in successful parasitism

Decoding bee cleptoparasitism through comparative transcriptomics of Coelioxoides waltheriae and its host Tetrapedia diversipes

Cleptoparasitism, also known as brood parasitism, is a widespread strategy among bee species in which the parasite lays eggs into the nests of the host species. Even though this behavior has significant ecological implications for the dynamics of several species, little is known about the molecular pathways associated with cleptoparasitism. To shed some light on this issue, we used gene expression data to perform a comparative analysis between two solitary neotropical bees: Coelioxoides waltheriae, an obligate parasite, and their specific host Tetrapedia diversipes. We found that ortholog genes involved in signal transduction, sensory perception, learning, and memory formation were differentially expressed between the cleptoparasite and the host. We hypothesize that these genes and their associated molecular pathways are engaged in cleptoparasitism-related processes and, hence, are appealing subjects for further investigation into functional and evolutionary aspects of cleptoparasitism in bees.

Antennal transcriptome analysis of Psyttalia incisi (silvestri) (Hymenoptera: Braconidae): identification and tissue expression profiling of candidate odorant-binding protein genes

BACKGROUND

Olfaction plays an important role in host-seeking by parasitoids, as they can sense chemical signals using sensitive chemosensory systems. Psyttalia incisi (Silvestri) (Hymenoptera: Braconidae) is the dominant parasitoid of Bactrocera dorsalis (Hendel) in fruit-producing regions of southern China. The olfactory behavior of P. incisi has been extensively studied; however, the chemosensory mechanisms of this species are not fully understood.

RESULTS

Bioinformatics analysis of 64,515 unigenes from the antennal transcriptome of both male and female adults P. incisi identified 87 candidate chemosensory genes. These included 13 odorant-binding proteins (OBPs), seven gustatory receptors (GRs), 55 odorant receptors (ORs), 10 ionotropic receptors (IRs), and two sensory neuron membrane proteins (SNMPs). Phylogenetic trees were constructed to predict evolutionary relationships between these chemosensory genes in hymenopterans. Moreover, the tissue expression profiles of 13 OBPs were analyzed by quantitative real-time PCR, revealing high expression of seven OBPs (1, 3, 6, 7, 8, 12, and 13) in the antennae.

CONCLUSION

This study represents the first identification of chemosensory genes and the determination of their expression patterns in different tissues of P. incisi. These results contribute to a better understanding of the function of the chemosensory system of this parasitoid species.

Identification of chemosensory genes by antennal transcriptome analysis and expression profiles of odorant-binding proteins in parasitoid wasp Aulacocentrum confusum

The endoparasitoid wasp, Aulacocentrum confusum (Hymenoptera: Braconidae), is a preponderant natural enemy of the larvae of Glyphodes pyloalis Walker (Lepidoptera: Pyralidae), which is a destructive pest of mulberry trees. We first constructed the antennal transcriptome database of A. confusum. In total, we obtained 48,262,304 clean reads from the dataset and assembled 24,324 unigenes. A total of 12,690 (52.17%) unigenes indicated significant similarity (E-value < 10^-5) compared to known protein sequences of other species from the NCBI non-redundant protein database. Gene ontology (GO) and cluster of orthologous groups (COG) analyses were used to determine the functional categories of these genes. A total of 84 putative chemosensory genes were identified from the antennal transcriptome of A. confusum, including 11 putative odorant-binding protein (OBP) genes, six chemosensory protein (CSP) genes, 44 olfactory receptor (OR) genes (including one olfactory co-receptor, Orco), 19 ionotropic receptor (IR) genes, and four sensory neuron membrane protein (SNMP) genes. Results of qPCR assays indicated that among of 11 AconOBPs, nine AconOBP genes were significantly expressed in the antennae of A. confusum adults. AconOBP8 was significantly expressed in the abdomen and AconOBP10 was highly expressed in the thorax. These findings can build a basis for further study on the processes of chemosensory perception in A. confusum at the molecular level.

Sensory gene identification in the transcriptome of the ectoparasitoid Quadrastichus mendeli

Sensory genes play a key role in the host location of parasitoids. To date, the sensory genes that regulate parasitoids to locate gall-inducing insects have not been uncovered. An obligate ectoparasitoid, Quadrastichus mendeli Kim & La Salle (Hymenoptera: Eulophidae: Tetrastichinae), is one of the most important parasitoids of Leptocybe invasa, which is a global gall-making pest in eucalyptus plantations. Interestingly, Q. mendeli can precisely locate the larva of L. invasa, which induces tumor-like growth on the eucalyptus leaves and stems. Therefore, Q. mendeli-L. invasa provides an ideal system to study the way that parasitoids use sensory genes in gall-making pests. In this study, we present the transcriptome of Q. mendeli using high-throughput sequencing. In total, 31,820 transcripts were obtained and assembled into 26,925 unigenes in Q. mendeli. Then, the major sensory genes were identified, and phylogenetic analyses were performed with these genes from Q. mendeli and other model insect species. Three chemosensory proteins (CSPs), 10 gustatory receptors (GRs), 21 ionotropic receptors (IRs), 58 odorant binding proteins (OBPs), 30 odorant receptors (ORs) and 2 sensory neuron membrane proteins (SNMPs) were identified in Q. mendeli by bioinformatics analysis. Our report is the first to obtain abundant biological information on the transcriptome of Q. mendeli that provided valuable information regarding the molecular basis of Q. mendeli perception, and it may help to understand the host location of parasitoids of gall-making pests.

An Ovarian Protein Involved in Passive Avoidance of an Endoparasitoid To Evade Its Host Immune Response

Through a combination of transcriptomic and proteomic analyses, we identified 817 secreted ovarian proteins from an endoparasitoid wasp, Cotesia chilonis, of which five proteins are probably involved in passive evasion. The results of an encapsulation assay revealed that one of these passive evasion-associated proteins (Crp32B), a homologue of a 32-kDa protein (Crp32) from C. rubecula, could protect resin beads from being encapsulated by host hemocytes in a dose-dependent manner. Crp32B is transcribed in ovarian cells, nurse cells, follicular cells, and oocytes, and the protein is located throughout the ovary and on the egg surface. Moreover, Crp32B has antigenic similarity to several host components. These results indicate that C. chilonis may use molecular mimicry as a mechanism to avoid host cellular immune response.

Analysis of the Gynaephora qinghaiensis pupae immune transcriptome in response to parasitization by Thektogaster sp

As a pest on the Qinghai-Tibet Plateau, Gynaephora qinghaiensis causes severe damage to grassland vegetation and its pupae are also natural hosts of Thektogaster sp. To successfully parasitize, endoparasitoids generally introduce or secrete multiple parasitic factors into the host body during the spawning stage to suppress the host immune response. To study the parasitic effects of Thektogaster sp. on G. qinghaiensis, a transcriptome analysis of immune-related genes in parasitized and nonparasitized G. qinghaiensis pupae was performed. A total of 371,260,704 clean reads were assembled into 118,144 unigenes with an average length of 884.33 base pairs. Of these, 23,660 unigenes were annotated in at least one database and 94,484 unigenes were not annotated in any databases. These findings indicated that the majority of the genetic resources (79.97% of all unigenes) in Gynaephora should be further explored. Parasitization significantly affected the transcriptional profile of G. qinghaiensis pupae. The present study identified 12,322 differentially expressed genes and 57 immune-related genes were identified in parasitized G. qinghaiensis pupae. Most immune-related genes were downregulated, potentially resulting from the inhibitory effect of Thektogaster sp. on G. qinghaiensis pupae after parasitization. Overall, the transcriptome analysis sheds valuable light on the molecular mechanisms of G. qinghaiensis parasitization by Thektogaster sp. and promotes the development of novel biocontrol strategies for Gynaephora based on immune defense.

Expression differences in Aphidius ervi (Hymenoptera: Braconidae) females reared on different aphid host species

The molecular mechanisms that allow generalist parasitoids to exploit many, often very distinct hosts are practically unknown. The wasp Aphidius ervi, a generalist koinobiont parasitoid of aphids, was introduced from Europe into Chile in the late 1970s to control agriculturally important aphid species. A recent study showed significant differences in host preference and host acceptance (infectivity) depending on the host A. ervi were reared on. In contrast, no genetic differentiation between A. ervi populations parasitizing different aphid species and aphids of the same species reared on different host plants was found in Chile. Additionally, the same study did not find any fitness effects in A. ervi if offspring were reared on a different host as their mothers. Here, we determined the effect of aphid host species (Sitobion avenae versus Acyrthosiphon pisum reared on two different host plants alfalfa and pea) on the transcriptome of adult A. ervi females. We found a large number of differentially expressed genes (between host species: head: 2,765; body: 1,216; within the same aphid host species reared on different host plants: alfalfa versus pea: head 593; body 222). As expected, the transcriptomes from parasitoids reared on the same host species (pea aphid) but originating from different host plants (pea versus alfalfa) were more similar to each other than the transcriptomes of parasitoids reared on a different aphid host and host plant (head: 648 and 1,524 transcripts; body: 566 and 428 transcripts). We found several differentially expressed odorant binding proteins and olfactory receptor proteins in particular, when we compared parasitoids from different host species. Additionally, we found differentially expressed genes involved in neuronal growth and development as well as signaling pathways. These results point towards a significant rewiring of the transcriptome of A. ervi depending on aphid-plant complex where parasitoids develop, even if different biotypes of a certain aphid host species (A. pisum) are reared on the same host plant. This difference seems to persist even after the different wasp populations were reared on the same aphid host in the laboratory for more than 50 generations. This indicates that either the imprinting process is very persistent or there is enough genetic/allelic variation between A. ervi populations. The role of distinct molecular mechanisms is discussed in terms of the formation of host fidelity.

Characterization of joining sites of a viral histone H4 on host insect chromosomes

A viral histone H4 (CpBV-H4) is encoded in a polydnavirus, Cotesia plutellae bracovirus (CpBV). It plays a crucial role in parasitism of an endoparasitoid wasp, C. plutellae, against diamondback moth, Plutella xylostella, by altering host gene expression in an epigenetic mode by its N-terminal tail after joining host nucleosomes. Comparative transcriptomic analysis between parasitized and nonparasitized P. xylostella by RNA-Seq indicated that 1,858 genes were altered at more than two folds in expression levels at late parasitic stage, including 877 up-regulated genes and 981 down-regulated genes. Among parasitic factors altering host gene expression, CpBV-H4 alone explained 16.3% of these expressional changes. To characterize the joining sites of CpBV-H4 on host chromosomes, ChIP-Seq (chromatin immunoprecipitation followed by deep sequencing) was applied to chromatins extracted from parasitized larvae. It identified specific 538 ChIP targets. Joining sites were rich (60.2%) in AT sequence. Almost 40% of ChIP targets included short nucleotide repeat sequences presumably recognizable by transcriptional factors and chromatin remodeling factors. To further validate these CpBV-H4 targets, CpBV-H4 was transiently expressed in nonparasitized host at late larval stage and subjected to ChIP-Seq. Two kinds of ChIP-Seqs shared 51 core joining sites. Common targets were close (within 1 kb) to genes regulated at expression levels by CpBV-H4. However, other host genes not close to CpBV-H4 joining sites were also regulated by CpBV-H4. These results indicate that CpBV-H4 joins specific chromatin regions of P. xylostella and controls about one sixth of the total host genes that were regulated by C. plutellae parasitism in an epigenetic mode.

Protein Discovery: Combined Transcriptomic and Proteomic Analyses of Venom from the Endoparasitoid Cotesia chilonis (Hymenoptera: Braconidae)

Many species of endoparasitoid wasps provide biological control services in agroecosystems. Although there is a great deal of information on the ecology and physiology of host/parasitoid interactions, relatively little is known about the protein composition of venom and how specific venom proteins influence physiological systems within host insects. This is a crucial gap in our knowledge because venom proteins act in modulating host physiology in ways that favor parasitoid development. Here, we identified 37 possible venom proteins from the polydnavirus-carrying endoparasitoid Cotesia chilonis by combining transcriptomic and proteomic analyses. The most abundant proteins were hydrolases, such as proteases, peptidases, esterases, glycosyl hydrolase, and endonucleases. Some components are classical parasitoid venom proteins with known functions, including extracellular superoxide dismutase 3, serine protease inhibitor and calreticulin. The venom contains novel proteins, not recorded from any other parasitoid species, including tolloid-like proteins, chitooligosaccharidolytic β-N-acetylglucosaminidase, FK506-binding protein 14, corticotropin-releasing factor-binding protein and vascular endothelial growth factor receptor 2. These new data generate hypotheses and provide a platform for functional analysis of venom components.

Suppressive activity of a viral histone H4 against two host chromatin remodelling factors: lysine demethylase and SWI/SNF

Histone H4, a nucleosome subunit in eukaryotes, plays crucial roles in DNA package and regulation of gene expression through covalent modification. A viral histone H4 encoded in Cotesia plutellae bracovirus (CpBV), a polydnavirus, is called CpBV-H4. It is highly homologous to other histone H4 proteins excepting 38 extra amino acid residues in the N terminus. CpBV-H4 can form octamer with other histone subunits and alter host gene expression. In this study, CpBV-H4 was transiently expressed in a natural host (Plutella xylostella) and its suppressive activity on host gene expression was evaluated by the suppressive subtractive hybridization (SSH) technique. The SSH targets down-regulated by CpBV-H4 were read with the 454 pyrosequencing platform and annotated using the genome of P. xylostella. The down-regulated genes (610 contigs) were annotated in most functional categories based on gene ontology. Among these SSH targets, 115 genes were functionally distinct, including two chromatin remodelling factors: a lysine-specific demethylase (Px-KDM) and a chromatin remodelling complex [Px-SWI/SNF (SWItch/Sucrose Non-Fermentable)]. Px-KDM was highly expressed in all tested tissues during the entire larval period. Suppression of Px-KDM expression by specific RNA interference (RNAi) significantly (P<0.05) reduced haemocyte nodule formation in response to immune challenge and impaired both larval and pupal development. Px-SWI/SNF was expressed in all developmental stages. Suppression of Px-SWI/SNF expression by RNAi reduced cellular immune response and interfered with adult metamorphosis. These results suggest that CpBV-H4 can alter host gene expression by interfering with chromatin modification and remodelling factors in addition to its direct epigenetic control activity.

Sensory and immune genes identification and analysis in a widely used parasitoid wasp Trichogramma dendrolimi (Hymenoptera: Trichogrammatidae)

Trichogramma dendrolimi Matsumura (Hymenoptera: Trichogrammatidae) is one of the preponderant egg parasitoids of Dendrolimus spp., which are important defoliators of coniferous forests. This parasitoid wasp has been widely released to control pine caterpillar and other lepidopteran pests, but its control efficiency needs to be improved. Sensory systems are crucial for T. dendrolimi to locate hosts, and immunity is probably involved after egg deposition in the host cavity; however, few reports have focused on the molecular mechanism of olfactory detection and survival of T. dendrolimi. It is necessary to identify these genes before further functional research is conducted. In this study, we assembled and analyzed the transcriptome of T. dendrolimi using next-generation sequencing technology. The sequencing and assembly resulted in 38 565 contigs with N50 of 3422 bp. Sequence comparison indicate that T. dendrolimi sequences are very similar to those of another parasitoid Nasonia vitripennis. Then the olfactory, vision, and immune-related gene families were identified, and phylogenetic analyses were performed with these genes from T. dendrolimi and other model insect species. Furthermore, phylogenetic tree with odorant binding proteins of T. dendrolimi and their host Dendrolimus was constructed to determine whether convergent evolution exists. These genes can be valid targets for further gene function research. The present study may help us to understand host location and survival mechanisms of T. dendrolimi and to use them more efficiently for pest control in the future.

Effects of the endoparasitoid Cotesia chilonis (Hymenoptera: Braconidae) parasitism, venom, and calyx fluid on cellular and humoral immunity of its host Chilo suppressalis (Lepidoptera: Crambidae) larvae

The larval endoparasitoid Cotesia chilonis injects venom and bracoviruses into its host Chilo suppressalis during oviposition. Here we study the effects of the polydnavirus (PDV)-carrying endoparasitoid C. chilonis (Hymenoptera: Braconidae) parasitism, venom and calyx fluid on host cellular and humoral immunity, specifically hemocyte composition, cellular spreading, encapsulation and melanization. Total hemocyte counts (THCs) were higher in parasitized larvae than in unparasitized larvae in the late stages following parasitization. While both plasmatocyte and granulocyte fractions and hemocyte mortality did not differ between parasitized and unparasitized hosts, in vitro spreading behavior of hemocytes was inhibited significantly by parasitism throughout the course of parasitoid development. C. chilonis parasitism suppressed the encapsulation response and melanization in the early stages. Venom alone did not alter cellular immune responses, including effects on THCs, mortality, hemocyte composition, cell spreading and encapsulation, but venom did inhibit humoral immunity by reducing melanization within 6h after injection. In contrast to venom, calyx fluid had a significant effect on cell spreading, encapsulation and melanization from 6h after injection. Dose-response injection studies indicated the effects of venom and calyx fluid synergized, showing a stronger and more persistent reduction in immune system responses than the effect of either injected alone.

Transcriptome Profile of the Asian Giant Hornet (Vespa mandarinia) Using Illumina HiSeq 4000 Sequencing: De Novo Assembly, Functional Annotation, and Discovery of SSR Markers

Vespa mandarinia found in the forests of East Asia, including Korea, occupies the highest rank in the arthropod food web within its geographical range. It serves as a source of nutrition in the form of Vespa amino acid mixture and is listed as a threatened species, although no conservation measures have been implemented. Here, we performed de novo assembly of the V. mandarinia transcriptome by Illumina HiSeq 4000 sequencing. Over 60 million raw reads and 59,184,811 clean reads were obtained. After assembly, a total of 66,837 unigenes were clustered, 40,887, 44,455, and 22,390 of which showed homologous matches against the PANM, Unigene, and KOG databases, respectively. A total of 15,675 unigenes were assigned to Gene Ontology terms, and 5,132 unigenes were mapped to 115 KEGG pathways. The zinc finger domain (C2H2-like), serine/threonine/dual specificity protein kinase domain, and RNA recognition motif domain were among the top InterProScan domains predicted for V. mandarinia sequences. Among the unigenes, we identified 534,922 cDNA simple sequence repeats as potential markers. This is the first transcriptomic analysis of the wasp V. mandarinia using Illumina HiSeq 4000. The obtained datasets should promote the search for new genes to understand the physiological attributes of this wasp.