Next Article in Journal
Reconstituting the Mammalian Apoptotic Switch in Yeast
Next Article in Special Issue
Modulation of Sex Pheromone Discrimination by a UDP-Glycosyltransferase in Drosophila melanogaster
Previous Article in Journal
An Experimental Hut Evaluation of PBO-Based and Pyrethroid-Only Nets against the Malaria Vector Anopheles funestus Reveals a Loss of Bed Nets Efficacy Associated with GSTe2 Metabolic Resistance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

De Novo Transcriptome Identifies Olfactory Genes in Diachasmimorpha longicaudata (Ashmead)

1
Key Laboratory of Integrated Pest Management on Tropical Crops, Ministry of Agriculture and Rural Affairs, Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571101, China
2
Guangxi Key Laboratory for Biology of Crop Diseases and Insect Pests, Institute of Plant Protection, Guangxi Academy of Agricultural Sciences, Nanning 530007, China
3
Department of Entomology, College of Tobacco Science, Guizhou University, Guiyang 550025, China
4
Guangxi Academy of Agricultural Sciences, Nanning 530007, China
5
College of Agriculture, Guangxi University, Nanning 530007, China
*
Authors to whom correspondence should be addressed.
Authors contribute equally.
Genes 2020, 11(2), 144; https://doi.org/10.3390/genes11020144
Submission received: 3 January 2020 / Revised: 17 January 2020 / Accepted: 22 January 2020 / Published: 29 January 2020
(This article belongs to the Special Issue Olfaction: From Genes to Behavior)

Abstract

:
Diachasmimoorpha longicaudata (Ashmead, D. longicaudata) (Hymenoptera: Braconidae) is a solitary species of parasitoid wasp and widely used in integrated pest management (IPM) programs as a biological control agent in order to suppress tephritid fruit flies of economic importance. Although many studies have investigated the behaviors in the detection of their hosts, little is known of the molecular information of their chemosensory system. We assembled the first transcriptome of D. longgicaudata using transcriptome sequencing and identified 162,621 unigenes for the Ashmead insects in response to fruit flies fed with different fruits (guava, mango, and carambola). We annotated these transcripts on both the gene and protein levels by aligning them to databases (e.g., NR, NT, KEGG, GO, PFAM, UniProt/SwissProt) and prediction software (e.g., SignalP, RNAMMER, TMHMM Sever). CPC2 and MIREAP were used to predict the potential noncoding RNAs and microRNAs, respectively. Based on these annotations, we found 43, 69, 60, 689, 26 and 14 transcripts encoding odorant-binding protein (OBP), chemosensory proteins (CSPs), gustatory receptor (GR), odorant receptor (OR), odorant ionotropic receptor (IR), and sensory neuron membrane protein (SNMP), respectively. Sequence analysis identified the conserved six Cys in OBP sequences and phylogenetic analysis further supported the identification of OBPs and CSPs. Furthermore, 9 OBPs, 13 CSPs, 3 GRs, 4IRs, 25 ORs, and 4 SNMPs were differentially expressed in the insects in response to fruit flies with different scents. These results support that the olfactory genes of the parasitoid wasps were specifically expressed in response to their hosts with different scents. Our findings improve our understanding of the behaviors of insects in the detection of their hosts on the molecular level. More importantly, it provides a valuable resource for D. longicaudata research and will benefit the IPM programs and other researchers in this filed.

1. Introduction

Diachasmimoorpha longicaudata (Ashmead, D. longicaudata) is a solitary species of parasitoid wasp of several fruit fly species and has been introduced to many countries as a biological control agent. Its host, Bactrocera dorsalis Hendel, can attack many fruit species and some other plants, such as Caricaceae, Moraceae, Myrtaceae, Rosaceae, and Solaneaceae [1]. It is said that the female Ashmead insects can detect the fly larvae by sound in rotting fruit and that the attractant could be the fungal fermentation products rather than the chemical substances produced by the fly larvae [2,3]. Carrasco and colleagues reported that the presence of fly larvae was essential for the orientation of wasps [4]. Further, it was proposed that cues from the fruit can be used by the wasps directly and that the presence of the host can enhance the attraction towards a patch [5]. Interestingly, chemical compounds produced by the larvae can be detected by wasps to locate the host [6,7]. Once the female parasitoid is on the fruit, a specific chemical compound released by some Tephritidae species can be used to enhance the host search [7]. However, much is unknown about the chemosensory system of parasitoid wasps in response to their hosts.
Animals can use the chemosensory system to detect and discriminate chemical cues in the environment [8]. The chemical sensors of insects are mainly from the antennae system, which is a highly specific and extremely sensitive chemical detector, and olfactory proteins in the antennae can be used by the insects to detect very low-abundance odorants from thousands of odors and further to guide their behaviors, such as forage, mate hunting, host plant location, shelter, and selection of spawning sites [9]. Since the first odorant-binding protein (OBP) was identified in Antheraea polyphemus [10], research on insect OBP has become a hot spot in the field of entomology. However, OBP-related studies have mainly been demonstrated in some important pests [11], and very few have been reported in parasitoid wasps, which are very important natural enemies of the pests.
The chemoreception of insects involves three important events: (i) The uptake of signal molecules from the external environment; (ii) transport (diffusion) through the sensory hair; and (iii) interaction with the chemoreceptor, which in turn activates the cascade of events leading to spike activity in sensory neurons [12]. Some important protein families have been reported to participate in these events, such as OBP, sensory neuron membrane protein (SNMP), chemosensory protein (CSP), odorant receptor (OR), gustatory receptor (GR), and odorant ionotropic receptor (IR) [12]. Insect OBPs are small globular proteins (~135 to 220 amino acids) and are characterized by a specific domain that comprises six α-helices joined by three disulphide bonds [12]. They can be categorized into two subgroups: Pheromone-binding proteins, which are mainly distributed in the male antenna, and general OBP (GOBP), which can be found in multiple tissues of male and female insects and function in the recognition of odorants from plants and other animals [13]. There are about 300 OBPs in the NCBI database and not many studies have been demonstrated to identify the OBPs in parasitoid wasps. Xu and colleagues used transcriptome sequencing and identified 1 CSP, 21 OBPs, 53 ORs, 29 IRs, and 4 SNMPs in Bactrocers minax [14], an oligophagous tephritid insect whose host selection, and oviposition behavior largely depend on the perception of chemical cues. Zhu et al. reported Sgua-OBP1 and Sgua-OBP2 in Scleroderma guani (Hymenoptera: Bethylidae) and NvitOBP in Nasonia vitripennis [15]. Zhang identified 10 OBPs in Microplitis mediator (Halidag) [16]. Zhao et al. identified 25 OBPs, 80 ORs, 10 IRs, 11 CSP, 1 SNMPs, and 17 GRs in adult male and female Chouioia cunea antennae [17]. However, little is known about the olfactory proteins in D. longicaudata.
Harbi and colleagues demonstrated a multistep assay (e.g., olfactory, laboratory, and semi-field trials) and reported the preference of medfly-infected fruits, including apple, orange, peach, and clementine mandarins [18]. This experiment supports that different olfactory genes are expressed in response to different fruit scents. In this study, we used transcriptome sequencing to study the olfactory genes in D. longicaudata. By similarity, we identified a number of OBPs, CSPs, ORs, IRs, SNMPs, and GRs expressed in the Ashmead insects. Our results also showed that different olfactory genes were expressed in the search of their host with different fruit scents. This is the first time to identify the gene and protein sequences for the olfactory products in this species. Our findings will provide a basis for future molecular studies and improve our understanding of the chemosensory system of parasitoid wasps.

2. Materials and Methods

2.1. Insect Rearing

The Ashmead and fruit fly larvae were obtained from the Institute of Plant Protection (IPP), Hainan Academy of Agricultural Sciences and maintained in the experimental fields. The fruit flies were fed with guava (G), mango (M), and carambola (C) separately. A mixture of yeast and sucrose (1:1) was used as supplementary nutrition for the fruit flies in the adult stage. Then, late-second and early-third instar fruit flies were used as hosts for the Ashmead insects. The Ashmead insects were fed with 15% honey water and clean water; the fifth, sixth and seventh generations of the Ashmead adults, which were parasitic to the fruit flies, maintained with G (G1~G3), M (M1~M3), and C (C1~C3), were used as biological replicates. The antenna, head, breast, abdomen, and feet tissues of the one male and two female wasps were mixed together for RNA extraction.

2.2. RNA Isolation, Library Construction, and Deep Sequencing

Total RNA was extracted from the insect tissues using the TRIzol reagent, as previously described [19,20]. The quality and quantity of total RNA were determined by multiple instruments, including a Nanodrop 2000 (Thermo Scientific, MA, USA), Qubit 4 Fluorometer (Invitrogen, CA, USA), and Agilent 2100 Bioanalyzer. Then, the total RNA (1 μg) of each sample was used to build the cDNA library using the TruSeq RNA Library Preparation Kit v2 protocol (Illumina, CA, USA), as described [20]. After the cDNA libraries were quality controlled by the Agilent 2100 Bioanalyzer and qRT-PCR, they were sequenced on the Illumina HiSeqXTEN platform with the paired-end 150 strategy. Raw sequencing reads of these samples can be accessed from the NCBI SRA platform under the accession numbers SRR10766480~SRR10766488.

2.3. Transcriptome Assembly

Raw sequencing reads were cleaned using the trim_galore v0.5.0 and the clean data was quality controlled by FASTQC v0.11.7 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Then, clean reads were used to assemble the transcriptome for each sample using Trinity software with default parameters, as described [21]. In detail, high quality RNA-Seq reads were used to generate overlapping k-mers (25) and Inchworm was used to assemble sorted k-mers into transcript contigs based on the (k-1)-mer overlaps. Next, Chrysalis was used to cluster related Inchworm contigs into components by using grouped raw reads and paired read links. Then, a de Bruijn graph for each cluster was built by Chrysalis and reads were partitioned among the clusters. Finally, Butterfly was used to process the individual graphs and ultimately report the full-length transcripts. To remove redundant sequences, CD-HIT was used to cluster the assembled highly similar transcripts into Unigenes [22], which can be accessed in the NCBI TAS platform under the accession number GIF00000000. BUSCO v4 was used to evaluate the completeness of the assembled Unigenes [23] using the eukaryota_odb10 dataset.

2.4. Annotation for the Transcriptome

We annotated the assembled transcriptome by aligning them to different databases. Initially, the transcriptome was searched against the NR (Non-Redundant Protein Sequence Database), NT (Nucleotide Sequence Database), KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway, gene ontology (GO), Pfam, and UniProt/SwissProt databases using BLAST software, and hits with an e-value > 1 × 10−5 were filtered. Then, BLAST2GO was used to retrieve the GO annotation in terms of the biological process, cellular component, and molecular function [24]. Using the enzyme commission numbers produced by BLAST2GO, we mapped the assembled transcriptome to the KEGG pathway database and obtained the pathway annotation. rRNA transcripts were predicted using RNAMMER [25].

2.5. Likely Protein Identification and Annotation

We next extracted the likely proteins from the assembled transcriptome using TransDecoder. Then, the likely proteins were searched against the UniProtKB/Swiss-Prot database to identify known proteins, functional PFAM domains were identified using HMMER [26], signal peptides were predicted using SignalP [27], and transmembrane domains were predicted using TMHMM Sever v2.0 [28]. The EggNOG database v4.1 [29] was searched against to identify proteins in EuKaryotic Orthologous Groups (KOG), Clusters of Orthologous Groups (COGs), and non-supervised orthologous groups (NOGs). All the annotation for the assembled genes and likely proteins were subjected to the Trinotate v3.1.1 (http://trinotate.github.io) for combination.
Based on the gene annotation and likely protein annotation, we obtained the Unigenes, which were annotated into olfactory gene families, such as OBP, OR, IR, GR, and SNMP, using their names as key words. For the CSP transcripts, we aligned the Unigenes to all the CSP transcripts from NCBI GenBank. Hits with an e-value > 1 × 10−5 were filtered.

2.6. Noncoding Genes and microRNA Genes

Unannotated Ashmead genes were processed by the Coding Potential Calculator (CPC v2) with default parameters to identify potential long noncoding genes [30]. To identify potential microRNA (miRNA) genes, we first mapped all the animal mature miRNAs to the noncoding genes with a maximal of two mismatches [31]. Then, MIREAP was used to predict the miRNA precursors and MIRANDA was used to predict the target genes of these miRNAs [32].

2.7. Differential Expression Analysis

We aligned the clean reads of each sample to the Unigenes using Bowtie2 and profiled the gene expression using RSEM [33]. The trimmed mean of the M-values (TMM) method was used for normalization and edgeR was used to identify differentially expressed genes with the following cut-offs [34]: Count > 5, log2 fold change (log2 fc) >1 or log2 fc < −1, p-value < 0.05, and false discovery rate (FDR) < 0.05.

2.8. Function Enrichment Analysis

We calculated the p-value (calculated by Fisher’s exact test) and q-value (calculated by the R package ‘qvalue’) for each GO term and KEGG pathway involved in the differentially expressed genes. Enriched terms should satisfy the following criteria: p-value < 0.05 and q-value < 0.05.

2.9. Phylogenetic Analysis

Phylogenetic trees were reconstructed for OBPs and CSPs using MEGA7 software [35]. We obtained the likely protein sequences for the top 5 highly expressed OBP and CSP transcripts. These sequences together with the homology protein sequences, obtained from NCBI, were subjected to MEGA7 to create phylogenetic trees using the neighbor-joining method. The bootstrap procedure based on 1000 replicates was used to assess node support, and the node support values < 50% were not shown. Figtree v1.4.3 (https://github.com/rambaut/figtree/) was used to visualize the results.

2.10. qRT-PCR Verification

We used real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) to validate the expression levels of three randomly selected transcripts (TRINITY_DN1020_c0_g1_i3, TRINITY_DN1284_c0_g1_i11, and TRINITY_DN500_c0_g1_i2). Forward and reverse primers of the three transcripts and the internal control (β-actin) were predicted using Prime3 and synthesized at BGI-Shenzhen. The procedure of the qRT-PCR experiment was the same as the previous study [21]. Each transcript was measured three times in every sample and three independent repeats were performed (n = 9). The Delta cycle threshold (ΔCt) was used to present the expression of a transcript in the sample and ΔΔCt was used to show the expression difference between two samples. We used the relative normalized expression (RNE) to show the expression changes: R N E = 2 Ct .

3. Results

3.1. Animal and Transcriptome Sequencing

After the Ashmead animals (three females and one male) were maintained with oriental fruit flies, which were fed guava (G), mango (M), and carambola ©, the antenna, head, breast, abdomen, and feet tissues were mixed together for RNA extraction and transcriptome sequencing. After data cleaning, we obtained a total of ~622.32 million reads (~69.15 million reads on average) and assembled 24,201 to 34,302 genes using Trinity for all the samples (Table 1). After similar genes/transcripts were clustered and merged, we finally obtained 162,621 Unigenes for the Ashmead transcriptome, with an average length of 1425.14 bp. The N50, GC content, and size of the transcriptome were calculated as 3572, 41.88%, and ~231 M, respectively (Table 1). Length distribution analysis showed that 53.05% of the total transcripts were longer than 500 bp and 11,101 transcripts (6.83% of the total transcripts) were longer than 5000 bp (Figure 1A). Last, we used BUSCO to evaluate the completeness of the assembled Unigenes and the results showed 99.6% of the assembled Unigenes were complete. In detail, out of the 255 evaluated BUSCOs in the dataset, 254 were complete, including 228 duplicated and 26 single-copy BUSCOs.

3.2. Transcriptome Annotation

We first annotated the assembled Ashmead transcriptome on the transcript level. All the transcripts were aligned to public databases for full annotation and Figure 1B showed 74,264, 68,185, 54,536, 58,109, 61,794, 58,113, and 22,546 transcripts were aligned to the NR, NT, UniProt/SwissProt, KEGG pathway, KOG, Pfam, and GO databases, respectively. Tthee NR mapping results (Figure 1C) showed the top 10 species aligned by the assembled Ashmead transcripts and the majority of the transcripts were aligned to Diachasma alloeum (45,569 transcripts), Fopius arisanus (5217 transcripts), and Rhinolophus sinicus (4920 transcripts). Unsurprisingly, the top two species together with Ashmead were all from the Braconidae family. GO annotation revealed that 11,131, 7722, and 9078 transcripts were involved in “binding”, “membrane”, and “cellular process”, respectively (Figure 1C). Then, we categorized the KEGG pathway annotation (Figure 1E) into six groups: Cellular progresses, environmental information processing, genetic information processing, human diseases, metabolism, and organism systems. Among them, “signal transduction” is the most significant pathway, which involved 9231 transcripts. KOG annotation also revealed that 14,221 transcripts were involved in the signal transduction mechanisms (Figure 1F).
Next, TranDecoder predicted 88,215 likely proteins encoded by the assembled Ashmead transcriptome and we annotated them using the Trinotate pipeline. It was shown that 59,887, 26,402, 53,725, and 52,572 were aligned to the UniProt/SwissProt, KOG, Pfam, and GO databases, respectively (Figure 1B). In addition, 6642 and 14,672 likely proteins were predicted to contain signal peptides and transmembrane helices (Figure 1B).
In addition, we predicted 74,963 of the unannotated transcripts using CPC2 had low coding probability. Interestingly, 2128 of the unannotated transcripts were predicted to be possible coding transcripts, which might be specific to the Ashmead and require further experiments to be verified. Next, we predicted 87 transcripts have the potential of producing miRNAs (Table S1). Notably, 84 of them were probably derived from the intron regions of coding genes while 1 and 3 derived from the Ashmead specific coding gene and long noncoding genes, respectively.

3.3. Olfactory Genes

We next identified genes encoding the olfactory gene from five families, including OBP, CSP, OR, GR, IR, and SNMP. In the Ashmead transcriptome, we found 43 transcripts encoding OBPs (Table 2, Table S2) by similarity and 35 of them were predicted to have complete ORFs by TransDecoder. The Ashmead OBPs were categorized into four sub-families: OBP-56, -69, -72, and -83 (Table S2). Multiple sequence alignment (Figure S1) showed these OBPs have 6 conserved cysteine residues (Cys) and SignalP predicted that 26 of these OBPs have the signal peptides located in the first 23 amino acids (aa) and 2 located in the first 37 aa. We identified that 69 transcripts had the potential of encoding CSPs in the Ashmead transcriptome by aligning the likely proteins to the known CSPs (Table 2, Table S2) and 63 of them had intact ORFs. SignalP identified 58 CSPs with the signal peptides in the first 28 aa while transmembrane domains were found in 13 CSPs. Surprisingly, 689 Ashmead transcripts had the capacity of encoding ORs and 115 transcripts were found to encode OR-13 (Table 2, Table S2). We also identified nine transcripts encoding OR coreceptors (OR-co). Out of the 347 OR transcripts that had intact ORFs, 276 were predicted to have transmembrane domains (Table 2). In the Ashmead transcriptome, there were 26 transcripts encoding IRs, including 3 IR-21, 18 IR-25, 4 IR-68, and 2 IR-93 (Table 2). Two thirds of the IR transcripts that had intact ORFs were predicted to have transmembrane domains. In addition, we identified 60 GR and 14 SNMP transcripts (Table 2). GR-28 and SNMP-1 were the largest group, which corresponded 19 and 11 transcripts, respectively. Further, we found 55 transcripts that had either the 7tm chemosensory receptor (PF02949) or GOBP (PF01395) PFAM domain in their protein sequences (Table 2, Table S2).

3.4. Gene Expression Profile

We next profiled the gene expression in the Ashmead insects maintained with the fruit flies fed with the three kinds of fruits. After lowly expressed genes (count < 5) were filtered, RSEM identified a total of 63,627 transcripts in the Ashmead animals, of which 46,607, 49,253, and 53,558 transcripts were distributed in G, M, and C, respectively (Table S3). The Venn diagram (Figure 2A) revealed 38,009 transcripts commonly expressed in all samples while 3115, 4571, and 8159 were specifically detected in G, M, and C, respectively. Interestingly, not all the olfactory genes were expressed in the insects and we found 39 OBPs, 65 CSPs, 29 GRs, 382 ORs, 15 IRs, and 13 SNMPs. Figure 2B,C showed the expression levels of these olfactory transcripts in the insects, and revealed that different olfactory genes of Ashmead insects are responsible for the fruit flies with different fruits. According to the average expression levels, we showed the top five highly expressed olfactory transcripts identified in this study (Table 3) and it was revealed that the identities of highly expressed olfactory transcripts were shared by the parasitoid wasps of the fruit flies fed with different fruits. Notably, OBP56 and OBP69 were highly expressed in the insects; IR25a and SNMP1 were the only highly expressed transcript for the IR and SNMP groups, respectively.

3.5. Phylogenetic Analysis

We next compared the sequences of olfactory proteins identified in Ashmead and some other homology organisms and analyzed their phylogenetic relationship. First, we constructed the phylogenetic tree for the top 5 OBPs (Table 3) and 28 other OBP sequences obtained from NCBI. Detailed information, including accession numbers and species, can be accessed in Table S3. The phylogenetic tree of OBPs (Figure 3A) showed high similarity (68.9% to 96.5%) between Ashmead OBPs and other species, such as Aethina tumida (Atumi), Aphidius gifuensis (Agifu), Cephus cinctus (Ccinc), Cotesia chilonis (Cchil). Diachasma alloeum (Dallo), Fopius arisanus (Faris), Megachile rotundata (Mrotu), Meteorus pulchricornis (Mpulc), Microplitis demolitor (Mdemo), and Microplitis mediator (Mmedi). We next performed the phylogenetic analysis for CSPs. The top 5 highly expressed CSPs (Table 3) were compared with 22 CSPs from other species, like Vespa velutina (Vvelu), Ccinc, Sclerodermus, Mpulc, and Yemma signatus (Ysign). Detailed accession numbers of these CSPs can be accessed in Table S3. It is clear that the Ashmead CSPs can be clustered with other known CSPs. The phylogenetic analysis supported the identification and characterization of OBP and CSP transcripts in this study.

3.6. Differential Expression Analysis

Another important goal of this study was to identify genes in the insects in response to their parasitic hosts, which had difference fruit scents. Using edgeR, we identified a total of 2650 transcripts differentially expressed in the Ashmead insects in response to the fruit flies with different scents (Table S4), and the number of differentially expressed transcripts can be seen in Figure 4A. Compared to G, there were 1466 upregulated and 53 downregulated transcripts identified in both C and M (Figure 4B). Some transcripts were specifically expressed in the insects when they were parasitic to the fruit flies with one fruit scent (Figure 4B). We next analyzed the differentially expressed transcripts encoding olfactory proteins in the Ashmead insects in response to the three fruit flies. A total of 58 transcripts encoding olfactory proteins were found, including 9 OBPs, 13 CSPs, 3 GRs, 4IRs, 25 ORs, and 4 SNMPs (Figure 4C). This evidence further supports the existence of multiple pathways of Ashmead insects in response to different fruits scents. In addition to olfactory proteins, some other protein families were differentially expressed in the wasps maintained with fruit flies supplied with different fruits, such as 2 LOC107047718 (Putative 7 transmembrane sweet-taste receptor of 3 gcpr), 88 ribosomal proteins, 109 transcription factors, 99 histones, 11 heat shock proteins, and 43 G-protein coupled receptors/regulators (Table S4). The differential expression of transcripts from different families indicated the complicated regulation mechanisms of parasitoid wasps in response to their hosts with different fruit scents. More experiments are required to explore their functions in this process.

3.7. qRT-PCR

We further used qRT-PCR to validate the expression levels of three randomly selected transcripts in the parasitoid wasps. The primers for the three transcripts and internal control (actin) were predicted using Prime3 and can be accessed in Table 4. We used Log2FC and RNE to show the expression changes of the transcripts in C, M, and G identified by RNA-Seq and qRT-PCR, respectively. Overall, the expression patterns of these transcripts in most comparisons were consistent by both RNA-Seq and qRT-PCR except three (TRINITY_DN1020_c0_g1_i3, TRINITY_DN500_c0_g1_i2 in C_vs_M, and TRINITY_DN1020_c0_g1_i3 in M_vs_G). The high agreement of the gene expression patterns in RNA-Seq and qRT-PCR indicates that the transcripts identified in this study might be functionally expressed in the parasitoid wasps maintained with fruit flies with different scents, which requires future functional experiments.

4. Discussion

This project was initiated in an effort to identify the olfactory proteins of parasitoid wasps. In recent years, many aspects about the perception of pheromones and other odorants have been elucidated [36]. Olfaction is used by insects to recognize volatile cues that allow the detection of food, predators, and mates [37]. We identified 43 OBPs, 69 CSPs, 60 GRs, 689 ORs, 26 IRs, and 14 SNMPs in D. longicaudata (Table 2, Table S2) and some of them were differentially expressed when they were maintained with fruit flies fed with different fruits (Figure 4C, Table S4). This is the first time the chemosensory genes in D. longicaudata have been investigated and our findings provide a basis for elucidating the molecular mechanisms of the olfactory-related behaviors of parasitoid wasps.
OBPs are a group of proteins that specialize in the transport of lipids. In this study, we identified 43 transcripts encoding OBPs and this number is similar to the number of genes encoding OBPs in the Drosophila genome [38]. Addittionally, classic OBPs have been reported to contain six Cys in their sequences [38,39]. Because all the identified OBPs in this study were general OBPs (GOBPs) (Table S2), we found conserved Cys in their sequences (Figure S1), which increased the confidence of using the transcriptome to identify the OBPs in D. longicaudata. In detail, the D. longicaudata OBPs were grouped into four subfamilies: OBP-56, -69, -72, and -83 (Table S2). However, no studies have been demonstrated to distinguish their functions. In general, according to their names, they are called GOBPs because they bind general odorants that are likely to be represented by the same volatiles for most of the species [40]. We also showed the differential expression of nine transcripts encoding OBPs in parasitoid wasps in response to the fruit flies with different scents (Figure 4, Table S4). This might indicate that some specific OBPs are expressed to discriminate different scents.
Similar to OBPs, CSPs are another group of proteins that mediate the olfactory recognition in insects [41]. They are thought to be expressed during nearly the whole life circle of insects [42,43]. The number of CSP genes varies in species. For example, only 4~8 CSP genes in ants, flies, bees, wasps, and anopheline mosquitoes [44]; 19~20 in butterfly, moth, and beetle [45,46,47]; and 27 to 83 in Culex mosquito species [48]. We identified 69 CSP transcripts produced by 58 genes (Table S2), the number of which is similar to Culex mosquito species. CSPs function not only in the chemical communication between insects and the environment but also in some other cellular processes, such as lipid transport, general immunity, insecticide resistance, and xenobiotic degradation [45,49]. In ants, CSPs have been proposed to mediate the recognition of chemical signatures composed of cuticular lipids [50]. The differential expression of CSP transcripts identified in this study (Table S4) may support the ability to recognize different scents from their hosts.
In the present study, we also identified some other olfactory gene families, such as OR, IR, GR, and SNMP, which were differentially expressed in response to the fruit fly with different scents (Table S2, Table S4). OR is the name for all molecules that are expressed in the cell membranes of olfactory receptor neurons and are responsible for the detection of odorants. The ORs form a multigene family consisting of around 800 genes in humans and 1400 genes in mice [51]. We identified 689 OR transcripts derived from 637 genes (Table 2, Table S2). The diversity of ORs might help insects to discriminate as many different odors as possible. GRs are found be expressed exclusively in gustatory receptor neurons [52]. However, many GRs are not related to taste receptors but function in the detection of sugars, bitter compounds, and non-volatile pheromones [53]. Interestingly, GR28B represents a new class of thermosensor and is required for thermotaxis [54]. We found GR28B differentially expressed in C and M (Figure 4C, Table S2). This might be evidence of its new role in the detection of different scents. IR genes are expressed in coeloconic sensilla of the antenna and respond, among others, to water and amines [55]. IRs are not related to insect ORs but rather have evolved from ionotropic glutamate receptors (iGluRs), a conserved family of synaptic ligand-gated ion channels [56]. In this study, we identified 26 IR transcripts, and 18 of them encode IR25A (Table 2, Table S2). It is conceivable that the IR25A ancestor initially evolved as a sensory detector for external glutamate, analogous to the synaptic function of iGluRs, and that it only later acquired a co-receptor function after duplication and diversification of the IR repertoire [56]. We also identified 14 SNMP transcripts in the parasitoid wasps, including 11 encoding SNMP1 and 3 encoding SNMP2 (Table 2, Table S2). The SNMP1 has been shown to be antenna specific and play an important role in pheromone detection [57]. While SNMP2, which acts as a second lepidoperan and also associates with pheromone-sensitive sensilla, has been shown to be expressed in sensilla support cells rather than neurons [58,59]. The identification and differential expression of olfactory-related transcripts revealed the complex chemosensory system of D. longicaudata and supported a diverse function of olfactory genes in discriminating different chemical cues.
The limitations of this project may include the use of the tissue mixture of insects. It is said that some OBPs are expressed in insects with a sex preference. For example, MsepOBP5 exhibited female-biased expression in 0- and 5-day-old adults; MsepOBP22 displayed female-biased expression in 0- and 5-day-old adults but was male-biased in 3-day-old adults [60]. Due to the difficulty of the sample preparation, it is hard to get enough material for sequencing with the same sex. Additionally, it is difficult to determine the tissue-specific olfactory genes. However, our findings provide a basis of future studies about the olfactory system in D. longicaudata.

5. Conclusions

In conclusion, we assembled the first transcriptome for D. longicaudata using transcriptome sequencing and identified 43 OBPs, 69 CSPs, 60 GRs, 689 ORs, 26 IRs, and 14 SNMPs. Further, 9 OBPs, 13 CSPs, 3 GRs, 4IRs, 25 ORs, and 4 SNMPs were differentially expressed in the insects in response to fruit flies with different scents. Our findings provide a basis towards understanding the molecular mechanisms of D. longicaudata in the detection of chemosensory cues. Additionally, the sequences including olfactory genes, noncoding genes, and miRNAs identified in this study can be used in the future and benefit other researchers in this field.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/2/144/s1, Figure S1: Multiple sequence alignment of OBP sequences identified in this study. It showed the conserved six Cys (C1~C6) in these sequences (highlighted in yellow and marked with *). Table S1: miRNA precursors identified in the parasitoid wasp transcriptome. Table S2: Expression profile of transcripts encoding olfactory proteins. Table S3: OBP and CSP sequences for the phylogenetic analysis. Table S4: Differentially expressed transcripts in the parasitoid wasps of fruit flies fed with different fruits.

Author Contributions

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

Funding

This research was funded by the National Natural Science Foundation of China (31560531) and Guangxi Innovation-driven Projects (AA17202017-2).

Acknowledgments

We acknowledge all the members from the Key Laboratory of Intergrated Pest Management on Tropical Crops, Ministry of Agriculture, Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences for their helpful discussions and kindness opinions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. White, I.M.; Elson-Harris, M.M. Fruit Flies of Economic Significance: Their Identification and Bionomics; CAB International: Wallingford, UK, 1992. [Google Scholar]
  2. Lawrence, P.O. Host vibration—A cue to host location by the parasite, Biosteres longicaudatus. Oecologia 1981, 48, 249–251. [Google Scholar] [CrossRef]
  3. Greany, P.D.; Tumlinson, J.H.; Chambers, D.L.; Boush, G.M. Chemically mediated host finding byBiosteres (Opius) longicaudatus, a parasitoid of tephritid fruit fly larvae. J. Chem. Ecol. 1977, 3, 189–195. [Google Scholar] [CrossRef]
  4. Carrasco, M.; Montoya, P.; Cruz-lopez, L.; Rojas, J.C. Response of the Fruit Fly Parasitoid Diachasmimorpha longicaudata (Hymenoptera: Braconidae) to Mango Fruit Volatiles. Environ. Entomol. 2005, 34, 576–583. [Google Scholar] [CrossRef]
  5. Silva, J.W.P.; Bento, J.M.S.; Zucchi, R.A. Olfactory response of three parasitoid species (Hymenoptera: Braconidae) to volatiles of guavas infested or not with fruit fly larvae (Diptera: Tephritidae). Biol. Control. 2007, 41, 304–311. [Google Scholar] [CrossRef]
  6. Duan, J.J.; Messing, R.H. Effects of host substrate and vibration cues on ovipositor-probing behavior in two larval parasitoids of tephritid fruit flies. J. Insect Behav. 2000, 13, 175–186. [Google Scholar] [CrossRef]
  7. Stuhl, C.; Sivinski, J.; Teal, P.; Paranhos, B.; Aluja, M. A compound produced by fruigivorous Tephritidae (Diptera) larvae promotes oviposition behavior by the biological control agent Diachasmimorpha longicaudata (Hymenoptera: Braconidae). Environ. Entomol. 2011, 40, 727–736. [Google Scholar] [CrossRef] [Green Version]
  8. De Bruyne, M.; Warr, C.G. Molecular and cellular organization of insect chemosensory neurons. Bioessays 2006, 28, 23–34. [Google Scholar] [CrossRef]
  9. Zubkov, S.; Gronenborn, A.M.; Byeon, I.J.; Mohanty, S. Structural consequences of the pH-induced conformational switch in A. polyphemus pheromone-binding protein: Mechanisms of ligand release. J. Mol. Biol. 2005, 354, 1081–1090. [Google Scholar] [CrossRef]
  10. Vogt, R.G.; Riddiford, L.M.; Prestwich, G.D. Kinetic properties of a sex pheromone-degrading enzyme: the sensillar esterase of Antheraea polyphemus. Proc. Natl. Acad. Sci. USA 1985, 82, 8827–8831. [Google Scholar] [CrossRef] [Green Version]
  11. Hull, J.J.; Perera, O.P.; Snodgrass, G.L. Cloning and expression profiling of odorant-binding proteins in the tarnished plant bug, Lygus lineolaris. Insect Mol. Biol. 2014, 23, 78–97. [Google Scholar] [CrossRef]
  12. Sanchez-Gracia, A.; Vieira, F.G.; Rozas, J. Molecular evolution of the major chemosensory gene families in insects. Heredity (Edinb.) 2009, 103, 208–216. [Google Scholar] [CrossRef] [PubMed]
  13. Jin, F.; Dong, X.; Xu, X.; Ren, S. cDNA cloning and recombinant expression of the general odorant binding protein II from Spodoptera litura. Sci. China C Life Sci. 2009, 52, 80–87. [Google Scholar] [CrossRef] [PubMed]
  14. Xu, P.; Wang, Y.; Akami, M.; Niu, C.Y. Identification of olfactory genes and functional analysis of BminCSP and BminOBP21 in Bactrocera minax. PLoS ONE 2019, 14, e0222193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Zhu, J.; Yang, P.; Wu, G.; He, Q.; Yin, L.; Xiong, Z. Clone and sequence analysis of the odorant-binding protein genes in Nasonia vitripennis. J. Environ. Entomol. 2010, 32, 476–482. [Google Scholar]
  16. Zhang, K. Characterization of MmedOBP8 and Its Binding Capacity in Microplitis mediator (Halidag); Huazhong Agricultural University: Wuhan, China, 2011. [Google Scholar]
  17. Zhao, Y.; Wang, F.; Zhang, X.; Zhang, S.; Guo, S.; Zhu, G.; Liu, Q.; Li, M. Transcriptome and Expression Patterns of Chemosensory Genes in Antennae of the Parasitoid Wasp Chouioia cunea. PLoS ONE 2016, 11, e0148159. [Google Scholar] [CrossRef] [Green Version]
  18. Harbi, A.; de Pedro, L.; Ferrara, F.A.A.; Tormos, J.; Chermiti, B.; Beitia, F.; Sabater-Munoz, B. Diachasmimorpha longicaudata Parasitism Response to Medfly Host Fruit and Fruit Infestation Age. Insects 2019, 10, 211. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, M.; Xu, R.; Rai, A.; Suwakulsiri, W.; Izumikawa, K.; Ishikawa, H.; Greening, D.W.; Takahashi, N.; Simpson, R.J. Distinct shed microvesicle and exosome microRNA signatures reveal diagnostic markers for colorectal cancer. PLoS ONE 2019, 14, e0210003. [Google Scholar] [CrossRef] [Green Version]
  20. Chen, M.; Xu, R.; Ji, H.; Greening, D.W.; Rai, A.; Izumikawa, K.; Ishikawa, H.; Takahashi, N.; Simpson, R.J. Transcriptome and long noncoding RNA sequencing of three extracellular vesicle subtypes released from the human colon cancer LIM1863 cell line. Sci. Rep. 2016, 6, 38397. [Google Scholar] [CrossRef] [Green Version]
  21. Wei, S.; Ma, X.; Pan, L.; Miao, J.; Fu, J.; Bai, L.; Zhang, Z.; Guan, Y.; Mo, C.; Huang, H.; et al. Transcriptome Analysis of Taxillusi chinensis (DC.) Danser Seeds in Response to Water Loss. PLoS ONE 2017, 12, e0169177. [Google Scholar] [CrossRef] [Green Version]
  22. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 2012, 28, 3150–3152. [Google Scholar] [CrossRef]
  23. Seppey, M.; Manni, M.; Zdobnov, E.M. BUSCO: Assessing Genome Assembly and Annotation Completeness. Methods Mol. Biol. 2019, 1962, 227–245. [Google Scholar] [CrossRef] [PubMed]
  24. Conesa, A.; Gotz, S.; Garcia-Gomez, J.M.; Terol, J.; Talon, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [Green Version]
  25. Lagesen, K.; Hallin, P.; Rodland, E.A.; Staerfeldt, H.H.; Rognes, T.; Ussery, D.W. RNAmmer: Consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35, 3100–3108. [Google Scholar] [CrossRef]
  26. Finn, R.D.; Clements, J.; Eddy, S.R. HMMER web server: Interactive sequence similarity searching. Nucleic Acids Res. 2011, 39, W29–W37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Petersen, T.N.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. Nat. Methods 2011, 8, 785–786. [Google Scholar] [CrossRef] [PubMed]
  28. Krogh, A.; Larsson, B.; von Heijne, G.; Sonnhammer, E.L. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J. Mol. Biol. 2001, 305, 567–580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Powell, S.; Forslund, K.; Szklarczyk, D.; Trachana, K.; Roth, A.; Huerta-Cepas, J.; Gabaldon, T.; Rattei, T.; Creevey, C.; Kuhn, M.; et al. eggNOG v4.0: Nested orthology inference across 3686 organisms. Nucleic Acids Res. 2014, 42, D231–D239. [Google Scholar] [CrossRef]
  30. Kong, L.; Zhang, Y.; Ye, Z.Q.; Liu, X.Q.; Zhao, S.Q.; Wei, L.; Gao, G. CPC: Assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007, 35, W345–349. [Google Scholar] [CrossRef]
  31. Li, R.; Yu, C.; Li, Y.; Lam, T.W.; Yiu, S.M.; Kristiansen, K.; Wang, J. SOAP2: An improved ultrafast tool for short read alignment. Bioinformatics 2009, 25, 1966–1967. [Google Scholar] [CrossRef] [Green Version]
  32. John, B.; Enright, A.J.; Aravin, A.; Tuschl, T.; Sander, C.; Marks, D.S. Human MicroRNA targets. PLoS Biol 2004, 2, e363. [Google Scholar] [CrossRef] [Green Version]
  33. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  36. Zhou, J.J. Odorant-binding proteins in insects. Vitam. Horm. 2010, 83, 241–272. [Google Scholar] [CrossRef] [PubMed]
  37. Zhou, S.S.; Sun, Z.; Ma, W.; Chen, W.; Wang, M.Q. De novo analysis of the Nilaparvata lugens (Stal) antenna transcriptome and expression patterns of olfactory genes. Comp. Biochem. Physiol. Part D Genomics Proteomics 2014, 9, 31–39. [Google Scholar] [CrossRef] [PubMed]
  38. Hekmat-Scafe, D.S.; Scafe, C.R.; McKinney, A.J.; Tanouye, M.A. Genome-wide analysis of the odorant-binding protein gene family in Drosophila melanogaster. Genome Res. 2002, 12, 1357–1369. [Google Scholar] [CrossRef] [Green Version]
  39. Gong, M.; Jiao, L.; Ma, W. Orthogonal immune algorithm with diversity-based selection for numerical optimization. In Proceedings of the First ACM/SIGEVO Summit on Genetic and Evolutionary Computation, Shanghai, China, 12–14 June 2009; pp. 141–148. [Google Scholar]
  40. Vogt, R.G.; Prestwich, G.D.; Lerner, M.R. Odorant-binding-protein subfamilies associate with distinct classes of olfactory receptor neurons in insects. J. Neurobiol. 1991, 22, 74–84. [Google Scholar] [CrossRef]
  41. Vogt, R.G.; Riddiford, L.M. Pheromone binding and inactivation by moth antennae. Nature 1981, 293, 161–163. [Google Scholar] [CrossRef]
  42. Picimbon, J.F.; Dietrich, K.; Angeli, S.; Scaloni, A.; Krieger, J.; Breer, H.; Pelosi, P. Purification and molecular cloning of chemosensory proteins from Bombyx mori. Arch. Insect. Biochem. Physiol. 2000, 44, 120–129. [Google Scholar] [CrossRef]
  43. Ma, C.; Cui, S.; Tian, Z.; Zhang, Y.; Chen, G.; Gao, X.; Tian, Z.; Chen, H.; Guo, J.; Zhou, Z. OcomCSP12, a Chemosensory Protein Expressed Specifically by Ovary, Mediates Reproduction in Ophraella communa (Coleoptera: Chrysomelidae). Front. Physiol. 2019, 10, 1290. [Google Scholar] [CrossRef] [Green Version]
  44. Liu, G.; Yue, S.; Rajashekar, B.; Picimbon, J.-F. Expression of chemosensory protein (CSP) structures in Pediculus humanis corporis and Acinetobacter A. baumannii. SOJ Microbiol. Infect. Dis. 2019, 7, 1–17. [Google Scholar]
  45. Xuan, N.; Guo, X.; Xie, H.Y.; Lou, Q.N.; Lu, X.B.; Liu, G.X.; Picimbon, J.F. Increased expression of CSP and CYP genes in adult silkworm females exposed to avermectins. Insect. Sci. 2015, 22, 203–219. [Google Scholar] [CrossRef] [PubMed]
  46. Ozaki, K.; Utoguchi, A.; Yamada, A.; Yoshikawa, H. Identification and genomic structure of chemosensory proteins (CSP) and odorant binding proteins (OBP) genes expressed in foreleg tarsi of the swallowtail butterfly Papilio xuthus. Insect. Biochem. Mol. Biol. 2008, 38, 969–976. [Google Scholar] [CrossRef] [PubMed]
  47. Liu, G.; Arnaud, P.; Offmann, B.; Picimbon, J.F. Genotyping and Bio-Sensing Chemosensory Proteins in Insects. Sensors (Basel) 2017, 17, 1801. [Google Scholar] [CrossRef] [Green Version]
  48. Mei, T.; Fu, W.B.; Li, B.; He, Z.B.; Chen, B. Comparative genomics of chemosensory protein genes (CSPs) in twenty-two mosquito species (Diptera: Culicidae): Identification, characterization, and evolution. PLoS ONE 2018, 13, e0190412. [Google Scholar] [CrossRef] [Green Version]
  49. Liu, G.; Ma, H.; Xie, H.; Xuan, N.; Guo, X.; Fan, Z.; Rajashekar, B.; Arnaud, P.; Offmann, B.; Picimbon, J.F. Biotype Characterization, Developmental Profiling, Insecticide Response and Binding Property of Bemisia tabaci Chemosensory Proteins: Role of CSP in Insect Defense. PLoS ONE 2016, 11, e0154706. [Google Scholar] [CrossRef]
  50. Ozaki, M.; Wada-Katsumata, A.; Fujikawa, K.; Iwasaki, M.; Yokohari, F.; Satoji, Y.; Nisimura, T.; Yamaoka, R. Ant nestmate and non-nestmate discrimination by a chemosensory sensillum. Science 2005, 309, 311–314. [Google Scholar] [CrossRef] [Green Version]
  51. Niimura, Y. Evolutionary dynamics of olfactory receptor genes in chordates: Interaction between environments and genomic contents. Hum. Genomics 2009, 4, 107–118. [Google Scholar] [CrossRef] [Green Version]
  52. Montell, C. Gustatory receptors: Not just for good taste. Curr. Biol. 2013, 23, R929–R932. [Google Scholar] [CrossRef] [Green Version]
  53. Montell, C. A taste of the Drosophila gustatory receptors. Curr. Opin. Neurobiol. 2009, 19, 345–353. [Google Scholar] [CrossRef] [Green Version]
  54. Ni, L.; Bronk, P.; Chang, E.C.; Lowell, A.M.; Flam, J.O.; Panzano, V.C.; Theobald, D.L.; Griffith, L.C.; Garrity, P.A. A gustatory receptor paralogue controls rapid warmth avoidance in Drosophila. Nature 2013, 500, 580–584. [Google Scholar] [CrossRef] [Green Version]
  55. Benton, R.; Vannice, K.S.; Gomez-Diaz, C.; Vosshall, L.B. Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell 2009, 136, 149–162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Rytz, R.; Croset, V.; Benton, R. Ionotropic receptors (IRs): Chemosensory ionotropic glutamate receptors in Drosophila and beyond. Insect Biochem. Mol. Biol. 2013, 43, 888–897. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Vogt, R.G.; Miller, N.E.; Litvack, R.; Fandino, R.A.; Sparks, J.; Staples, J.; Friedman, R.; Dickens, J.C. The insect SNMP gene family. Insect. Biochem. Mol. Biol. 2009, 39, 448–456. [Google Scholar] [CrossRef] [PubMed]
  58. Rogers, M.E.; Krieger, J.; Vogt, R.G. Antennal SNMPs (sensory neuron membrane proteins) of Lepidoptera define a unique family of invertebrate CD36-like proteins. J. Neurobiol. 2001, 49, 47–61. [Google Scholar] [CrossRef] [PubMed]
  59. Forstner, M.; Gohl, T.; Gondesen, I.; Raming, K.; Breer, H.; Krieger, J. Differential expression of SNMP-1 and SNMP-2 proteins in pheromone-sensitive hairs of moths. Chem. Senses 2008, 33, 291–299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Chang, X.Q.; Nie, X.P.; Zhang, Z.; Zeng, F.F.; Lv, L.; Zhang, S.; Wang, M.Q. De novo analysis of the oriental armyworm Mythimna separata antennal transcriptome and expression patterns of odorant-binding proteins. Comp. Biochem. Physiol Part. D Genom. Proteom. 2017, 22, 120–130. [Google Scholar] [CrossRef]
Figure 1. Overview of the assembled Ashmead transcriptome and annotation. (A) Length distribution of the assembled Ashmead transcriptome. (B) Number of transcripts and likely proteins aligned to different databases. (C) Number of the transcripts aligned to other species in the NR mapping results. (D) GO annotation for the assembled Ashmead transcriptome. (E) KEGG pathway annotation for the assembled Ashmead transcriptome. (F) KOG annotation for the assembled transcriptome.
Figure 1. Overview of the assembled Ashmead transcriptome and annotation. (A) Length distribution of the assembled Ashmead transcriptome. (B) Number of transcripts and likely proteins aligned to different databases. (C) Number of the transcripts aligned to other species in the NR mapping results. (D) GO annotation for the assembled Ashmead transcriptome. (E) KEGG pathway annotation for the assembled Ashmead transcriptome. (F) KOG annotation for the assembled transcriptome.
Genes 11 00144 g001
Figure 2. Expression levels of olfactory transcripts in the Ashmead insects. (A) Number of transcripts identified in the insects stimulated by three fruits. (B) Heat maps of the expression levels of OBP, CSP, and OR transcripts. (C) Heat maps showing the expression levels of transcripts encoding GRs, IRs, and SNMPs in the Ashmead insects.
Figure 2. Expression levels of olfactory transcripts in the Ashmead insects. (A) Number of transcripts identified in the insects stimulated by three fruits. (B) Heat maps of the expression levels of OBP, CSP, and OR transcripts. (C) Heat maps showing the expression levels of transcripts encoding GRs, IRs, and SNMPs in the Ashmead insects.
Genes 11 00144 g002
Figure 3. Phylogenetic trees for OBPs (A) and CSPs (B). Ashmead proteins are highlighted in red, the percentage represents the bootstrap value, and the scale bar represents the evolutionary distance.
Figure 3. Phylogenetic trees for OBPs (A) and CSPs (B). Ashmead proteins are highlighted in red, the percentage represents the bootstrap value, and the scale bar represents the evolutionary distance.
Genes 11 00144 g003
Figure 4. Differentially expressed transcripts in the parasitoid wasps of fruit flies fed with different fruits. (A) Number of differentially expressed transcripts in the insects in response to two fruit scents. (B) Venn diagram of up- (upper panel) and downregulated (lower panel) transcripts identified in C and M, compared to G. (C) A heat map showed the differential expression of 9 OBPs, 13 CSPs, 3 GRs, 4IRs, 25 ORs, and 4 SNMPs in the parasitoid wasps of fruit flies fed with different fruits.
Figure 4. Differentially expressed transcripts in the parasitoid wasps of fruit flies fed with different fruits. (A) Number of differentially expressed transcripts in the insects in response to two fruit scents. (B) Venn diagram of up- (upper panel) and downregulated (lower panel) transcripts identified in C and M, compared to G. (C) A heat map showed the differential expression of 9 OBPs, 13 CSPs, 3 GRs, 4IRs, 25 ORs, and 4 SNMPs in the parasitoid wasps of fruit flies fed with different fruits.
Genes 11 00144 g004
Table 1. Overview of the de novo transcriptome of Ashmead.
Table 1. Overview of the de novo transcriptome of Ashmead.
G1G2G3M1M2M3C1C2C3
Clean reads67,492,57674,983,47670,212,51878,796,85861,848,87672,961,79671,005,59254,703,28670,316,132
Assembled genes34,30226,03225,39824,20133,37533,08133,30824,77224,707
Assembled transcripts49,91440,91239,52838,08950,51751,70450,53038,85639,637
Unigenes162,621
Mean Length (bp)1425.14
N503572
GC (%)41.88
Total bases231,757,796
Expressed transcripts10,73610,32010,645986010,86311,69516,82711,76514,728
Table 2. Olfactory genes identified in the Ashmead transcriptome.
Table 2. Olfactory genes identified in the Ashmead transcriptome.
TypeTranscriptsSub-FamilyIntact ORFSignalPTMHMM
OBP4321 OBP-56, 7 OBP-69, 4 OBP-72, 11 OBP-83352815
CSP695 CSP-1, 3 CSP-3, 2 CSP-4, 5 CSP-5, 10 CSP-6, 3 CSP-7, 4 CSP-8, 37 CSP *635013
OR68945 OR-1, 6 OR-1F12, 21 OR-10, 115 OR-13, 42 OR-2, 46 OR-22, 2 OR-23, 9 OR-24, 2 OR-245, 1 OR-260, 2 OR-266, 1 OR-277, 15 OR-30, 6 OR-33, 31 OR-4, 5 OR-42, 14 OR-43, 2 OR-45, 24 OR-46, 47 OR-47, 26 OR-49, 2 OR-59, 1 OR-5, 1 OR-63, 48 OR-67, 1 OR-69, 6 OR-71, 13 OR-7, 1 OR-81, 62 OR-82, 37 OR-85, 11 OR-92, 11 OR-94, 5 OR-98, 4 OR-9, 9 OR-co, 1 OR-142, 14 OR *34725276
IR263 IR-21, 18 IR-25, 4 IR-68, 2 IR-9318412
GR606 GR-2, 5 GR-107, 3 GR-15, 1 GR-23, 19 GR-28, 4 GR-2, 8 GR-43, 2 GR-64, 3 GR-66, 8 GR *26125
SNMP1411 SNMP-1, 3 SNMP-31019
Others a5416 CR, 38 GOBP101222
a Transcripts containing the chemosensory receptor or the GOBP domain, by PFAM annotation. * Products are not specified.
Table 3. Top five transcripts encoding olfactory proteins in Ashmead insects.
Table 3. Top five transcripts encoding olfactory proteins in Ashmead insects.
TranscriptIDG aM aC aProtein/GeneDescription
OBP
TRINITY_DN22_c0_g3_i141,518.8955,018.9939,228.41OB56DGeneral odorant-binding protein 56d
cluster_contig467020,459.1127,957.5983,228.40OB56HGeneral odorant-binding protein 56h
cluster_contig373815,552.0413,752.5820,997.13OB69AGeneral odorant-binding protein 69a
TRINITY_DN3641_c0_g2_i15,010.5013,074.194,838.31OB69AGeneral odorant-binding protein 69a
TRINITY_DN23528_c0_g1_i13853.935367.956228.73OB56DGeneral odorant-binding protein 56d
CSP
TRINITY_DN1018_c0_g1_i511,065.7415,620.3110,030.62CSPchemosensory protein
TRINITY_DN661_c1_g1_i19447.6710,194.337724.67THK33221.1chemosensory protein 4
TRINITY_DN4258_c0_g1_i14538.674738.335976.67AZQ24964.1chemosensory protein, partial
TRINITY_DN1848_c0_g1_i65121.013247.274919.50CSPchemosensory protein
cluster_contig147483679.273944.222395.00THK33222.1chemosensory protein 5
OR
cluster_contig633041,364.8568,243.0524,104.78LOC107043577Odorant receptor
TRINITY_DN1357_c0_g1_i518,042.9723,565.9311,831.46OR43AOdorant receptor 43a
TRINITY_DN1357_c0_g1_i67765.7927,596.1913,421.25OR43AOdorant receptor 43a
cluster_contig1246816,698.8023,528.668172.54LOC107043576Odorant receptor
cluster_contig97592835.242248.793479.03OR43AOdorant receptor 43a
GR
TRINITY_DN9119_c0_g1_i2455.68489.30545.92GR107gustatory receptor Gr107
TRINITY_DN3525_c0_g1_i10288.94333.98230.17GR43Agustatory receptor for sugar taste 43a-like
cluster_contig36290.1043.8230.95GR43agustatory receptor for sugar taste 43a-like
TRINITY_DN14061_c0_g1_i391.79116.0578.74GR43agustatory receptor for sugar taste 43a-like
TRINITY_DN4338_c0_g1_i490.2277.2264.01GR107gustatory receptor Gr107
IR
TRINITY_DN11328_c0_g1_i141.0939.2745.12IR25AIonotropic receptor 25a
cluster_contig850539.7734.6539.58IR25AIonotropic receptor 25a
TRINITY_DN4424_c0_g1_i216.2612.9247.84IR25AIonotropic receptor 25a
cluster_contig1175219.2125.5718.80IR25AIonotropic receptor 25a
TRINITY_DN13912_c1_g2_i120.8810.7120.37IR25AIonotropic receptor 25a
SNMP
TRINITY_DN1934_c0_g1_i13567.613632.315117.35SNMP1Sensory neuron membrane protein 1
TRINITY_DN4669_c1_g1_i1176.96204.92206.30SNMP1Sensory neuron membrane protein 1
TRINITY_DN1867_c0_g2_i186.7283.9850.41XP_015114627.1Sensory neuron membrane protein 1-like
cluster_contig426470.2073.2346.97SNMP1Sensory neuron membrane protein 1
cluster_contig23744.5490.1034.07SNMP1Sensory neuron membrane protein 1
a Normalized expression, TMM.
Table 4. qRT-PCR validation. Log2 FC represents the log2 values of the fold change of a transcript identified by the RNA-Seq while RNE represents the relative normalized expression of a transcript identified by qRT-PCR.
Table 4. qRT-PCR validation. Log2 FC represents the log2 values of the fold change of a transcript identified by the RNA-Seq while RNE represents the relative normalized expression of a transcript identified by qRT-PCR.
TRINITY_DN1020_c0_g1_i3TRINITY_DN1284_c0_g1_i11TRINITY_DN500_c0_g1_i2
PrimersForwardCAACTTCAAGAACAATCCGACAACACTTATAGACGCATGCCAAGACCGACGTCGCTATGAACGCTTG
ReverseCCACAGCCAGAGACACAGCGGGCTGGAGAACGGGGATGGATTCTGATTTCCAGTACGAATACG
C_vs_GLog2 FC10.009.840
p-value4.12 × 10−286.02 × 10−271
RNE1.132.320.42
C_vs_MLog2 FC−9.471.436.43
p-value2.28 × 10−230.00131.98 × 10−33
RNE1.680.660.00
M_vs_GLog2 FC−0.538.326.66
p-value0.38561.13 × 10−157.18 × 10−34
RNE0.673.5279.36

Share and Cite

MDPI and ACS Style

Tang, L.; Liu, J.; Liu, L.; Yu, Y.; Zhao, H.; Lu, W. De Novo Transcriptome Identifies Olfactory Genes in Diachasmimorpha longicaudata (Ashmead). Genes 2020, 11, 144. https://doi.org/10.3390/genes11020144

AMA Style

Tang L, Liu J, Liu L, Yu Y, Zhao H, Lu W. De Novo Transcriptome Identifies Olfactory Genes in Diachasmimorpha longicaudata (Ashmead). Genes. 2020; 11(2):144. https://doi.org/10.3390/genes11020144

Chicago/Turabian Style

Tang, Liangde, Jimin Liu, Lihui Liu, Yonghao Yu, Haiyan Zhao, and Wen Lu. 2020. "De Novo Transcriptome Identifies Olfactory Genes in Diachasmimorpha longicaudata (Ashmead)" Genes 11, no. 2: 144. https://doi.org/10.3390/genes11020144

APA Style

Tang, L., Liu, J., Liu, L., Yu, Y., Zhao, H., & Lu, W. (2020). De Novo Transcriptome Identifies Olfactory Genes in Diachasmimorpha longicaudata (Ashmead). Genes, 11(2), 144. https://doi.org/10.3390/genes11020144

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