Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2016; 11(6): e0157684.
Published online 2016 June 22. doi:  10.1371/journal.pone.0157684
PMCID: PMC4917224

Deep Sequencing-Based Transcriptome Analysis Reveals the Regulatory Mechanism of Bemisia tabaci (Hemiptera: Aleyrodidae) Nymph Parasitized by Encarsia sophia (Hymenoptera: Aphelinidae)

Nicolas Desneux, Editor


The whitefly Bemisia tabaci is a genetically diverse complex with multiple cryptic species, and some are the most destructive invasive pests of many ornamentals and crops worldwide. Encarsia sophia is an autoparasitoid wasp that demonstrated high efficiency as bio-control agent of whiteflies. However, the immune mechanism of B. tabaci parasitization by E. sophia is unknown. In order to investigate immune response of B. tabaci to E. Sophia parasitization, the transcriptome of E. sophia parasitized B. tabaci nymph was sequenced by Illumina sequencing. De novo assembly generated 393,063 unigenes with average length of 616 bp, in which 46,406 unigenes (15.8% of all unigenes) were successfully mapped. Parasitization by E. sophia had significant effects on the transcriptome profile of B. tabaci nymph. A total of 1482 genes were significantly differentially expressed, of which 852 genes were up-regulated and 630 genes were down-regulated. These genes were mainly involved in immune response, development, metabolism and host signaling pathways. At least 52 genes were found to be involved in the host immune response, 33 genes were involved in the development process, and 29 genes were involved in host metabolism. Taken together, the assembled and annotated transcriptome sequences provided a valuable genomic resource for further understanding the molecular mechanism of immune response of B. tabaci parasitization by E. sophia.


The whitefly Bemisia tabaci (Hemiptera: Aleyrodidae), is well known as a worldwide invasive pest and may cause severe damage to various vegetables by feeding on phloem sap and transmitting many viruses [1]. It is a complex species containing at least 30 cryptic species [2]. B and Q-types are two most economically damaging and invasive species [3]. There are many studies focus on biological characterization, resistance, invasive mechanism, and biological control of B. tabaci [412]. Over the past years, B. tabaci has demonstrated a remarkable resistance to many groups of chemical insecticides [1316]. Due to the rapid resistance development, it is necessary to explore an alternative and effective management strategy to control B. tabaci. Parasitoid or parasitoid–produced regulatory molecules can be used to improve conventional pest control strategies.

Endoparasitoids have been identified as very important natural enemies of various arthropods, and could be used as biological control agents[1719]. Hymenopteran endoparasitoids deposit their eggs into the host insect haemocoel, whose larvae feed on the host until its death [2021]. Encarsia sophia is one of the specific parasitoids of Aleyrodidae species and has been used as efficacious classic biological control agents in many regions [22]. It can parasitize all instar nymphs of B. tabaci, especially the third and fourth instar nymphs [23]. The female wasp is generated by a bisexual process, but the male wasp is produced by autoparasitism [24]. Homogeneous E. sophia prefers to lay male eggs in the host parasitized by the heterogeneous wasp. When E. sophia and other kinds of wasps are raised or released together, the antecedent colonizers should inhibit the colonization of followers [25]. Previous studies have shown that E. sophia has strong plasticity adaption abilities[26].

However, the relationships between endoparasitoids and their hosts are complicated and involve long-term co-evolution. Many studies have investigated parasitoid biological characteristics, chemical communication, phylogenetic co-evolution, and physiological responses [27]. An increasing number of researchers have focused on revealing the physiological mechanism underlying the parasite induced immune defensive system and the biological development of hosts in order to estimate the co-evolution process between parasitoids and their hosts [2831]. Although several reports have concentrated on the molecular regulation mechanisms, there have only been a few descriptions of related, functional genes [32,33]. Furthermore, the limitations of previous research methods has led to the development of high-throughput RNA sequencing technology (RNA-Seq)[34].

RNA-Seq is widely used to obtain transcriptomes of the organism, tissue, or organ, to identify genes that were regulated under certain conditions, and to reveal the regulatory mechanisms in different organisms [3539]. In recent years, RNA-Seq has increasingly being applied in the biological agents to reveal the interaction mechanisms in the complex parasitoid-host system. Transcriptome profiling of organism under parasitization helps us to obtain a better understanding of host responses and effect on host’s growth, development. As a model species, Drosophila melanogaster and its parasitoid wasp Asobara tabida (Hymenoptera: Braconidae) is a well-studied system. Most genes associated with insect immunity appeared to be differentially expressed after wasp parasitized [40]. Most transcriptome studies on parasitoid-host systems have focused on Lepidoptera and Coleoptera, such as Plutella xylostella, Chilo suppressalis, Tenebrio molitor and Octodonta nipae [4144]. A previous study showed that another parasitoid, Eretmocerus mundus may parasitize B. tabaci and induce the specific transcription of functional genes related to immune responses in the host [45]. However, the host manipulation by the parasitoid is species-specific, and the molecular mechanism of immune system in B. tabaci parasitization by E. sophia has not yet been explored. In this study, we used deep sequencing to explore B. tabaci response to E. sophia parasitization. Our results demonstrate that immune- and metabolic-related genes that are differentially expressed in parasitized versus non-parasitized B. tabaci nymph.

Materials and Methods

Insects Rearing and Parasitization

The biotype Q of Bemesia tabaci was obtained from the greenhouse at the Beijing Academy of Agriculture and Forestry. All experimental populations were derived from one pairs of newly emerged B. tabaci female and male. In our laboratory, the B. tabaci was reared on cotton plants (Zhong-mian-suo 49) in insect proof cages at 26 ± 1°C, and with a photoperiod of 15L: 9D. The purity of the cultures was monitored every three to five generations using the random amplified polymorphic DNA-polymerase chain reaction technique with COI gene [46]. E. sophia was obtained from the greenhouse at Beijing Academy of Agriculture and Forestry. All whitefly instar nymph stages were provided as hosts to E. sophia. Then approximate fifty E. sophia (female to male ratio of 8:1) individuals were released into cages to breed and the newly emerged female and male as parents for five generations breeding.

Thirty pairs of whiteflies were fed on cotton leaf in a micro insect cage and the fresh cotton leaf were provided every 24 hours. When they had reached later 3rd or early 4th instar, they were transferred in culture dish with a piece of cotton leaf, whose petiol were wrapped into soggy cotton, and then the mated E. sophia was released into B. tabaci rearing cage for parasitization. Sixty paired E. sophia were released into one culture dish. Wasp E. sophia were removed after 2 hours parasitization. The first group of samples was collected at 24-hr after parasitization (24AP). At this time period, the parasitoids were at the egg stage in which the embryo had formed and gradually began to move. The brown substance in the egg began to accumulate and chorion had appeared. In other words, the parasitoid possessed immune regulation ability, but the ability was not strong at the egg stage. Therefore, we could identify the immune defense response of the host against the parasitoid. The second sampling period was 72-hr after parasitization (72AP) when the wasps reach larval stage move around and absorb nutrition from the host. At this time, E. sophia may start to regulate host development and metabolism to finish their own development in whiteflies. Each treatment and control had three replicates.

cDNA Library Construction and Illumina Sequencing

Total RNA was extracted from all nymph samples using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instruction and treated with DNaseI. The concentration and integrity of RNA sample were determined using 2100 Bioanalyzer (Agilent Technologies). The first- and second- strand cDNA synthesis, end reparation, addition of “A” bases to 3' ends, ligation of adapters at the end of DNA fragments, and PCR amplification. The cDNA library was qualified and quantified with an Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-time PCR system, respectively, and then sequenced using the Illumina HiSeq2000 platform at the Beijing Genomics Institute (BGI, Shenzhen, China).

Transcriptome Analysis

In order to obtain clean reads, the low quality and adapter-polluted reads were removed from raw data. The good quality reads were assembled using Trinity[47] and assembled sequences were output as unigenes. All raw sequencing data have been deposited in NCBI Sequence Read Archive (SRA) database ( with the following accession numbers: SRR1909644 (24AP), SRR1909651(72AP), SRR1909652 (CK-24AP), and SRR1909653 (CK-72AP). All the open reading frames (ORF) of unigene in B. tabaci were identified. If a unigene had many ORFs, we selected the longest one.

The unigenes were used for BlastX search and annotation against the NCBI non-redundant (nr) (, Swiss-Prot (, Kyoto Encyclopedia of Genes and Genome (KEGG, databases with an E-value cut-off of 10-5. Gene Ontology (GO) annotation of unigenes was analyzed using the Blast2Go software [48], and GO functional classification for all unigens was performed using the WEGO software [49]. In the absence of B. tabaci and E. sophina genome sequences, we selected eight transcriptome datasets of B. tabaci from the NCBI database, and try to utilize the annotation that were the most closely related to B. tabaci gene in the parasitized library.

Differentially Expressed Gene (DEG) Analysis

In order to find all the differentially expressed genes, the same FPKM (Fragments Per Kilobase per Million fragments) value of unigene was first calculated for the treatment and control groups [50]. The results were displayed as fold changes, p-values and q-values. According to the q-value (p-value’s statistical result after PFR (Positive False Rate) correction), a q-value less than 0.05 or the absolute value of fold change greater than 2 represented a significant difference between the treatment and the control.

Quantitative Real-time PCR (qRT-PCR) Validation

The quantitative real-time PCR technique was used to verify the reliability of the deep sequencing. Nine differentially expressed genes were randomly selected. The β-actin gene was used for normalization. The four RNA samples represented nymphs at 24AP and 72AP, and their respective control (non-parasitized nymphs) at the same developmental stages.

First-strand cDNA was synthesized from the total RNA (1.2 μg) by using PrimeScripTM 1st Strand cDNA Synthesis Kit (TaKaRa) with oligo (dT)18 as primer following the manufacture’s protocols. The reaction system consisted of 10 μl of SYBR Green, 0.4 μl of ROX, 2 μl of diluted cDNA, 0.4 μl of each primer and 6.8 μl of distilled water. The reactions were loaded on the CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA) under the following conditions: 50°C for 2 min; 95°C for 2 min; and 40 cycles of 95°C for 10s, 60°C for 15s, and 72°C for 20s, followed by melting curve generation (68°C to 95°C). Data analysis was performed by one-way ANOVA following by Tukey’s test using SPSS software.

Results and Discussion

Illumina Sequencing and de novo Assembly

In order to know how E. sophia parasitization regulated B. tabaci development, immune-response and the differences in regulatory mechanisms between E. sophia egg- and larvae-stages. Approximately, 35 million and 39 million reads were generated from non-parasitized and parasitized B. tabaci nymphs at 24AP, respectively, and 31million and 37 million reads were from non-parasitized and parasitized B. tabaci nymphs at 72AP, respectively. De novo assembly produced 292,696 B. tabaci unigenes with an average size of 616 bp. Of these unigenes, 35.96% were between 200 and 300bp, 27.43% were between 300 and 500bp, 22.65% were between 500 and 1000bp and 13.96% had nucleotide lengths above 1000bp (Fig 1).

Fig 1
Distribution of unigene lengths in the B. tabaci transcriptome.

Functional Annotation and Classification

For functional annotation, the 292,696 unigenes were aligned to the GenBank protein databases with a cut-off E-value of 10−5 using BLASTx. Using this approach, 46,406 unigenes (15.8% of all unigenes) were successfully mapped. In order to predicate protein function, the unigenes were further given a gene ontology (GO) classification and subjected to KEGG pathway analysis. A total of 35,688 unigenes were annotated and assigned to GO terms, which consisted of three main categories: biological process, cellular component and molecular function. A total of 11,993 unigenes were categorized as cellular components, 12,102 unigenes were grouped under the molecular function, and 11,593 unigenes under biological processes. KEGG pathway analysis indicated that there were 4,721 unigenes assigned to different pathways in which translation, signal transduction, neurodegenerative diseases, infectious diseases, and endocrine system were the main B. tabaci pathways after E. sophia parasitizzation.

Enrichment Analysis of DEGs

A total of 1,482 genes appeared to be significantly differentially expressed in the parasitized and non-parasitized B. tabaci, of which 852 genes were differentially up-regulated and 630 genes were differentially down-regulated (Fig 2A). At 24AP, there were 584 genes differentially expressed, of which 356 genes were up-regulated and 228 genes were down-regulated. At 72AP, there were 1,270 genes differentially expressed, of which 698 genes were up-regulated and 572 genes were down-regulated. Out of all of regulated genes, 202 up- and 170 down-regulated genes were found at both time points (Fig 2B) and more genes were up-regulated than that of the down-regulated genes at both 24AP and 72AP (Fig 2A). Furthermore, there was a significant difference in the numbers of differentially expressed genes at 24 hours than at 72 hours after parasitization. When E. sophia emerge in the larvae stage, more genes seemed to be involved in regulatory responses as compared to the egg stage. During the larvae stage, the parasitoid could move freely and began to feed on the host tissues. The distribution of the regulated genes indicated that their expression levels (>95%) were between two- to six-fold higher than at the egg stage (24 AP). Only a few genes changed more than six-fold (Fig 2).

Fig 2
General information about genes that were differentially expressed in response to parasitization.

GO analysis revealed that the DEGs were mainly categorized in the cellular component cluster, that focus on macromolecular, organelle, and cellular levels. In the molecular function cluster, the DEGs were mainly found in structural molecule, binding, and catalytic activity. In the biological process cluster, the DEGs were mainly categorized in cellular and metabolic processes, and cellular component organization or biogenesis (Fig 3). In addition, more genes were involved in cellular processes, metabolic processes, single-organism processes, response to stimuli, biological regulation, localization, and cellular component organization or biogenesis at 72AP. Translation and signal transduction were the two most important pathways according to the KEGG pathways analysis. For KEGG enrichment analysis, genes involved in the immune system, nervous system, endocrine system, and metabolic activities were differentially expressed. The above results showed that parasitization had a great impact on the normal life activities of the host.

Fig 3
GO annotation of differentially expressed genes at 24APand 72AP (level 2).

Effects of Parasitism on the Transcription of Host Immune-related Genes

Vertebrates have a set of immune defense mechanisms that include innate immunity and adaptive immunity, but invertebrates only have innate immunity protection [51]. Insects will initiate their innate immune response when encounting foreign agents, such as bacteria, fungi, virus, and parasitoid. The immune system of insects can be divided into two categories: 1) humoral defense, including the antimicrobial peptides, reactive intermediates of oxygen, melanin formation and clotting; and 2) cellular defense mainly based on haemocytes, such as phagocytosis, encapsulation, microaggregation and nodulation [5254]. Two defense mechanisms are associated with a wide range of immune-related genes.

Our sequencing results indicated that E. sophia parasitism had a significant impact on the transcription of immune-related genes in B. Tabaci nymph (Table 1). We identified several up-regulated genes with homologs known to be involved in immune responses in insects, such as: defensin, knottin, serpin I2, laminin, spectrin, and apolipophorin. Defensin is an antimicrobial peptide, which acts as an innate immunity effector molecule and provides the first protection from pathogen infection. After parasitization by E. sophia at 24AP, we found that the transcription of defensin was up-regulated in B. tabaci nymph. Our results were consistent with previous studies that the mRNA levels of defensin in D. melanogaster and Phlebotomus duboscqi were significantly increased after parasitization [55, 56]. Although the main action targets of defensin are bacteria and fungi, it also plays a role in the host-parasitoid system. Knottins are mini proteins that are present in many different organisms and have various biological functions [57]. After parasitization by E. sophia at both 24AP and 72AP, four knottins were over-expressed. Like defensin, it is also an important antimicrobial peptide.

Table 1
Immune-related genes differentially expressed in B. tabaci after being parasitized by E. sophia.

Serpin I2 was one of the genes having higher levels of up-regulation (5.33-fold) at 72AP in RNA-seq analysis. Quantitative RT-PCR analysis (Fig 4) also show that it was up-regulated by 6.85 folds. Serine proteases are important immune regulatory proteins which play a significant role in the activation of the prophenoloxidase (PPO) cascade. The cascade activation eventually causes melanization to kill parasitized wasp through choking [58], however, serine protease inhibitor (serpin) can prevent the serine proteases activated melanization and weaken host defense for wasp parasization. Although studies have shown that serpins can be regulated by the parasitoids infestation in many hosts, their transcriptional levels are different in different parasitoid-host systems, and even in the same parasitoid-host system, two opposite situations may occur. Mahadav et al. and Song et al. found that serpins were down-regulated in parasitized B. tabaci nymphs and P. xylostella larvae [26,59], while Etebari et al.[41] discovered that serpins were up-regulated 2- to 7-fold after P. xylostella parasitization by Diadegma semiclausum. In C. chilonis parasitized C. suppressalis, three up-regulated and three down-regulated serpins were identified in the fatbody [42]. Different serpins may play different roles in immune defense.

Fig 4
A qRT-PCR analysis of nine randomly selected genes from B.tabaci that showed relative expressions at 24AP and 72AP.

Cellular immunity is another important component of the insect immune system. Laminin can stimulate cell adhesion and cell movement. Cofilin is an actin-binding protein which promotes cell migration and movement by changing the adhesion between cells and the extracellular matrix. Actin plays a significant role in facilitating cellular activities. The up-regulation of these genes showed that the host enhanced hemocyte encapsulation by reinforcing the extension and adhesion of hemocytes. laminin, cofilin, and actin were identified in our study and they were over-expressed at 24AP and 72AP. Ras3 and Rho1 are related to cellular immunity in D. melanogaster [40,60]. In our study, two genes were also identified significantly differentially expressed after parasization which may also be involved in the immune response of B. tabaci. At 24AP and 72AP, these genes were consistently over-expressed, which indicated that the cellular immunity not only defend parasitoid embryo and larval attacking.

Superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) are three common enzymes in organisms. Organisms produce reactive oxygen species (ROS) under environmental stresses, which is cytotoxic to cells. However, the organism utilizes these protective enzymes to eliminate redundant ROS and protect themselves from damage [61]. When Trichoplusia ni is infected by baculoviruses, the expression of manganese superoxide dismutase (MnSOD) significantly reduces oxidative damage [62]. Zhu et al. also discovered that the transcriptional levels of Tenebrio molitor superoxide dismutases were up-regulated following bacterial infection or parasitization by Scleroderma guani [63]. In our study, the expression of SODC, PERO, and CATA were suppressed two-to four-fold at 72AP, but did not show significant changes at 24AP. At 72AP, the E. sophia reaches the larval stage and the damage to host becomes worse than that of egg stage. The decrease in the transcription of protective enzymes showed the parasitoid immune suppressive strategy.

In our study, some genes involved in insecticide resistance or detoxification were found to be differentially expressed under parasitization. Although these genes have no direct connection with defense against parasitoids attacking, they can be regarded as a stress response caused by parasitoid secretions. Takeda et al.[64] confirmed that the activities of glutathione-S-transferase (GST) and cytochrome P450 (CYP) increased in parasitized P. xylostella larvae. We also found most of the cytochrome P450 genes were highly expressed after being parasitized by E. sophia. Some genes were over-expressed at 24AP and others were over-expressed at 72AP.

Heat shock proteins (HSPs) are recognized as a family of highly conserved chaperones which respond to all kinds of environmental stress factors, such as heat, toxins, UV radiation, and invading pathogens by protecting protein from misfolding and denaturation [65]. We identified three heat shock protein genes, which were homologous with D. melanogaster and Anopheles albimanus, and are involved in B. tabaci development. In addition, heat shock 70 kDa protein cognate 3 was found to participate in the immune response. Therefore, we deduced that heat shock protein families can defend the host from damage by participating in the immune response and B. tabaci development.

Effects of parasitism on the transcription of host development-related genes

The parasitoids complete their development by absorbing the host’s hemolymph and tissues. However, the development of the parasitoid and the host are synchronous. A previous study found that Aphidiu servi parasitized Acyrthosiphon pisum late-stage nymph stopped growth [30]. B. tabaci late nymphs parasitized by Encarsia bimaculata also stop growing [66]. After Encarsia formosa parasitized Trialeurodes vaporariorum Westwood nymph, the wasp didn’t molt to until host nymph reached to last instar [67]. In order to complete development, the parasitoids have to change the host’s development to match their own growth. In some cases, parasitoids suppress host’s development and accelerate the host’s early-maturity [68], while, other parasitoids prolong host’s development to meet their own developmental needs. Previous studies have proposed that the wasp might control host’s development through regulating the juvenile hormone and ecdysone levels [69,70].

Juvenile hormone epoxide hydrolases (JHEHs) have been identified as regulatory proteins in the catabolism of juvenile hormones [71,72]. A previous study showed that JHEH transcript levels were down-regulated more than two-fold in P. xylostella after parasitization by D. semiclausum [41]. However, Wu et al. [42] discovered that JHEH and juvenile hormone esterase (JHE) transcript levels increased in C. suppressalis after C. chiilonis parasitization. Based on our transcriptome data, parasitization by E. sophia led to JHEH1 up-regulation at 72AP (Table 2 & Fig 4). However, up-regulation of larvae cuticle protein and down-regulation of pupal cuticle protein might imply that the parasitoid suppressed the host’s development. Thus, the high concentrations of JH may lead to up-regulation of JHEHs and their activity in order to maintain the balance.

Table 2
Developmental-related genes differentially expressed in B. tabaci after being parasitized by E. sophia.

Effects of parasitism on the transcription of host metabolism-related genes

Stearoyl-CoA desaturase (SCD) is an endoplasmic reticulum enzyme that catalyzes the biosynthesis of monounsaturated FA from saturated FA [73]. SCD inactivation causes obesity and abnormal lipid metabolism and one SCD activity, SCD1, was induced by insulin, but inhibited by leptin [74]. We found that at 24AP and 72AP, the transcript levels of SCD in B. tabaci nymph were up-regulated 6.69 and 4.52 times, respectively (Table 3). Furthermore, genes involved in the insulin signaling pathway were also significantly up-regulated. Our result implied that the wasp regulated the lipid metabolism of the host in order to get more nutrients available in host and meet their own needs. A report showed that the wasp preferred to parasitize late instar larvae because of adequate nutrition [75]. Stearoyl-CoA desaturase is an essential enzyme for the parasitic Trypanosoma brucei, and RNA interference of SCD caused a reduction of the parasitemia and an increase in host survival [76]. Environmental stress can influence the organism’s metabolism, same as parasitoid infestation, which is energetically consumption process [77]. We found a high number of differentially expressed transcripts were related to organism metabolism. Metabolic changes occurred at both time points, but a greater amount and different kinds of genes were affected at 72APthan at24AP.

Table 3
Metabolism-related genes differentially expressed in B. tabaci after being parasitized by E. sophia.

E.sophia infection influenced carbohydrate, lipid, and energy metabolism in the host. Some studies have found that trehalose content changed after parasitization [78,79]. In two treatment groups, maltose was degraded to glucose under the action of maltase. Beside the upregulation of maltase, other genes, in the citrate cycle and glycolysis, were over-expressed at 72AP, such as citrate synthase, aconitate hydratase, and glyceraldehydes 3-phosphate dehydrogenase. Glycolysis and the citrate cycle are carbohydrate metabolisms to produce ATP. Citrate synthase, aconitate hydratase and succinyl-CoA synthetase are three essential enzymes in the TCA cycle. Up-regulation of genes that control the synthesis of these enzymes showed that total ATP decreased in the organism. Therefore, the insect needed to obtain more energy by increasing the reaction rate of the TCA cycle. In addition, we found that the transcriptional level of cytochrome c oxidase and f-type H+-transporting ATPase were significantly enhanced. Cytochrome c oxidase is involved in ATP synthesis as a terminase of the mitochondrial inner membrane respiratory chain [80]. However, whether it can be regarded as evidence of enhanced respiration is not clear. Measurement of respiration rate should be investigated in future studies.

There are three types of ion transporting ATPases: P-type, V-type, and F-type. In organisms, their main function is to synthesize ATP and transport H+ [81]. There were more over-expressed f-type H+-ATPases at 72AP than at 24AP. This suggests that E. sophia parasitization of B. tabaci involved increased energy consumption. The host was regulated to produce more energy to supply to the parasitoid. Visser et al. found most wasps lacked a lipid synthesis mechanism and could not accumulate energy [82]. Therefore, it is reasonable to assume that the parasitoid may continually obtain energy from the host in order to complete its development.


In summary, our study first presented comprehensive transcriptome profiles of B. tabaci in response to E. sophia parasitization using RNAseq. The most of differentially expressed genes of B. tabaci after parasization have potential roles in immunity, development and metabolism to meet parasitoids needs. The transcriptome profiles provided a basis for future research in elucidate the host-parasitoid interaction. In addition, the identified immune-, development and detoxification–related genes may be target for B. tabaci control.


This work was supported by National Basic Research Program of China (Grant Number: 2013CB127600) to ZF. We would like to thank DCS, YCL, JWC, and LY for their bioinformatics analysis of the data.

Funding Statement

This research was supported by the National Key Basic Research and Development Program (973 program 2013CB127600) to ZF.

Data Availability

Data Availability

All raw sequencing data have been deposited in NCBI Sequence Read Archive (SRA) database ( with the following accession numbers: SRR1909644 (24AP), SRR1909651(72AP), SRR1909652 (CK-24AP), and SRR1909653 (CK-72AP).


1. Oliveira M, Henneberry T, Anderson P (2001) History, current status, and collaborative research projects for Bemisia tabaci. Crop Prot 20: 709–723.
2. Li ZX (2006) Molecular phylogenetic analysis reveals at least five genetic races of Bemisia tabaci in China. Phytoparasitica 34:431–440.
3. Costa H, Brown J (1991) Variation in biological characteristics and esterase patterns among populations of Bemisia tabaci, and the association of one population with silverleaf symptom induction. Entomologia experimentalis et applicata 61:211–219.
4. Horowitz AR, Kontsedalov S, Khasdan V, Ishaaya I (2005) Biotypes B and Q of Bemisia tabaci and their relevance to neonicotinoid and pyriproxyfen resistance. Arch Insect Biochem Physiol 58:216–225. [PubMed]
5. Wan F, Zhang G, Liu S, Luo C, Chu D, Zhang Y, et al. (2009) Invasive mechanism and management strategy of Bemisia tabaci (Gennadius) biotype B: Progress report of 973 Program on invasive alien species in China. Science in China Series C: Life Sciences 52:88–95. doi: 10.1007/s11427-008-0135-4 [PubMed]
6. Hoddle M, Van Driesche R, Sanderson J (1998) Biology and use of the whitefly parasitoid Encarsia formosa. Annu Rev Entomol 43:645–669. [PubMed]
7. Xiao Y, Chen J, Cantliffe D, Mckenzie C, Houben K, Osborne LS (2011) Establishment of papaya banker plant system for parasitoid, Encarsia sophia (Hymenoptera: Aphilidae) against Bemisia tabaci (Hemiptera: Aleyrodidae) in greenhouse tomato production. Biol Control 58:239–247.
8. Jaworski CC, Chailleux A, Bearez P, Desneux N (2015) Apparent competition between major pest reduces pest population densities on tomato crop, but not yield loss. J Pest Sci 88: 793–803.
9. Zhang YB, Yang NW, Sun LY, Wan FH (2015) Host instar suitability in two invasive whiteflies for the naturally occurring parasitoid Eretmocerus hayati in China. J Pest Sci 88: 225–234.
10. Lin KJ, Lu YH, Wan P, Yang YZ, Wyckhuys KAG, Wu KM (2015) Simultaneous reduction in incidence of Bemisia tabaci (Hemiptera: Aleyrodidae) and Sylepta derogate (Lepidoptera:Pyralidae) using velvetleaf, Abutilon theophrasti as a trap crop. J Pest Sci. 88:49–56.
11. Cruz PL, Baldin ELL, de Castro MJP (2014) Characterization of antibiosis to the silverleaf whitefly Bemisia tabaci biotype B (Hemiptera: Aleyrodidae) in cowpea entries. J Pest Sci. 87:639–645
12. Moreno-Ripoll R, Gabarra R, Symondson WOC, King RA, Agusti N (2014) Do the interations among natural enemies compromise the biological control of the whitefly Bemisia tabaci? J Pest Sci 87:133–141.
13. Cahill M, Gorman K, Day S, Denholm I (1996) Baseline determination and detection of resistance to imidacloprid in Bemisia tabaci (Homoptera: Aleyrodidae). Bull Entomol Res 86: 343–349.
14. Horowitz AR, Kontsedalov S, Ishaaya I (2004) Dynamics of resistance to the neonicotinoids acetamiprid and thiamethoxamthiamethoxam in Bemisia tabaci (Homoptera: Aleyrodidae). J Econ Entomol 97: 2051–2056 [PubMed]
15. Nauen R, Denholm I (2005) Resistance of insect pests to neonicotinoid insecticides: current status and future prospects. Arch Insect Biochem Physiol 58: 200–215. [PubMed]
16. Liang P, Tian YA, Biondi A, Desneux N, Gao XW (2012) Short-term and transgenerational effects of the neonicotinoid nitenpyram on susceptibility to insecticides in two whitefly species. Ecotoxicology 21:1889–1898. doi: 10.1007/s10646-012-0922-3 [PubMed]
17. Turlings TC, Tumlinson JH, Lewis WJ (1990) Exploitation of herbivore-induced plant odors by host-seeking parasitic wasps. Science 250: 1251–1253. [PubMed]
18. Han P, Niu CY, Desneux N (2014) Identification of top-down forces regulating cotton aphid population growth in transgenic Bt cotton in centra China. PLoS ONE. 9:e102980 doi: 10.1371/journal.pone.0102980 [PMC free article] [PubMed]
19. Ali A, Desneux N, Lu YH, Liu B, Wu KM (2016) Characterization of the natural enemy community attacking cotton aphid in the Bt cotton ecosystem in northern China. Sci Rep 6: 24273 doi: 10.1038/srep24273 [PMC free article] [PubMed]
20. Askew RR, Shaw MR (1986) Parasitoid communities: their size, structure and development In: Waage J, Greathead D, editors. Insect Parasitoids. London: Academic; pp. 225–263.
21. Desneux N, Barta RJ, Hoelmer KA, Hopper KR, Heimpel GE (2009) Multifaceted determination of host specificity in an aphid parasitoid. Oecologia 160:387–398. doi: 10.1007/s00442-009-1289-x [PubMed]
22. Simmons AM, Abd-Rabou S (2005) Parasitism of Bemisia tabaci (Homoptera: Aleyrodidae) after multiple releases of Encarsia sophia (Hymenoptera: Aphelinidae) in three vegetable crops. J Agr Urban Entomol 22:73–77.
23. Luo C, Liu TX (2011) Fitness of Encarsia sophia (Hymenoptera: Aphelinidae) parasitizing Trialeurodes vaporariorum and Bemisia tabaci (Hemiptera: Aleyrodidae). Insect Sci 18:84–91.
24. Zang LS, Liu TX, Wan FH (2011) Reevaluation of the value of autoparasitoids in biological control. PLoS ONE 6: e20324 doi: 10.1371/journal.pone.0020324 [PMC free article] [PubMed]
25. Collier TR, Hunter MS (2001) Lethal interference competition in the whitefly parasitoids Eretmocerus eremicus and Encarsia sophia. Oecologia 129:147–154.
26. Mahadav A, Gerling D, Gottlieb Y, Czosnek H, Ghanim M (2008) Parasitization by the wasp Eretmocerus mundus induces transcription of genes related to immune response and symbiotic bacteria proliferation in the whitefly Bemisia tabaci. BMC Genomics 9. [PMC free article] [PubMed]
27. Chen XX, He JH, Shi ZH, Ma Y, Lou YG, Zhu ZR (2000) Research review of interaction between parasitoids and their host. China Entomology of 21th Century——The Entomological Society of China academic annual conference paper highlights, 242–246.
28. Zhang QQ, Wang F, Fang Q, Ye GY (2011) In vitro cellular response of Pieris rapae(Lepidoptera: Pieridae) hemocytes and the effects of Pteromalus puparum venom. Acta Entomologica Sinica 54: 1264–1273.
29. Stettler P, Trenczek T, Wyler T, Pfister-Wilhelm R, Lanzrein B (1998) Overview of parasitism associated effects on host haemocytes in larval parasitoids and comparison with effects of the egg-larval parasitoid Chelonus inanitus on its host Spodoptera littoralis. J Insect Physiol 44:817–831. [PubMed]
30. Pennacchio F, Digilio M, Tremblay E (1995) Biochemical and metabolic alterations in Acyrthosiphon pisum parasitized by Aphidius ervi. Arch Insect Biochem Physiol 30: 351–367.
31. Gan M, Miao XX, Ding DC (2003) Biochemical and metabolic alterations in Aphis craccivora Koch parasitized by Lysiphlebus japonicas Ashmead. Acta Phytophylacica Sinica 30: 255–260.
32. Crozatier M, Ubeda J-M, Vincent A, Meister M (2004) Cellular immune response to parasitization in Drosophila requires the EBF orthologue collier. PLoS Biol 2: e196 [PMC free article] [PubMed]
33. Williams MJ, Ando I, Hultmark D (2005) Drosophila melanogaster Rac2 is necessary for a proper cellular immune response. Genes to Cells 10: 813–823. [PubMed]
34. Ozsolak F, Milos PM (2011) RNA sequencing: advances, challenges and opportunities. Nat Rev Genet 12: 87–98. doi: 10.1038/nrg2934 [PMC free article] [PubMed]
35. Qiu P, Benbow L, Liu SX, Greene JR, Wang LQ (2002) Analysis of a human brain transcriptome map. BMC Genomics 3:1–6. [PubMed]
36. Zhang F, Guo H, Zheng H, Zhou T, Zhou Y, Wang S, et al. (2010) Massively parallel pyrosequencing-based transcriptome analyses of small brown planthopper (Laodelphax striatellus), a vector insect transmitting rice stripe virus (RSV). BMC Genomics 11:303 doi: 10.1186/1471-2164-11-303 [PMC free article] [PubMed]
37. Zeng V, Ewen-Campen B, Horch HW, Roth S, Mito T, Extavour CG (2013) Developmental gene discovery in a hemimetabolous insect: De novo assembly and annotation of a transcriptome for the Cricket Gryllus bimaculatus. PLoS ONE 8(5). [PMC free article] [PubMed]
38. Gandhe AS, Arunkumar KP, John SH, Nagaraju J (2006) Analysis of bacteria-challenged wild silkmoth, Antheraea mylitta (lepidoptera) transcriptome reveals potential immune genes. BMC Genomics 7: 1–10. [PMC free article] [PubMed]
39. Yang N, Xie W, Jones CM, Bass C, Jiao X, Yang X, et al. (2013)Transcriptome profiling of the whitefly Bemisia tabaci reveals stage-specific gene expression signatures for thiamethoxam resistance. Insect Mol Biol 22: 485–496. doi: 10.1111/imb.12038 [PMC free article] [PubMed]
40. Wertheim B, Kraaijeveld AR, Schuster E, Blanc E, Hopkins M, Pletcher SD, et al. (2005) Genome-wide gene expression in response to parasitoid attack in Drosophila. Genome Biol 6: R94 [PMC free article] [PubMed]
41. Etebari K, Palfreyman RW, Schlipalius D, Nielsen LK, Glatz RV, Asgari S (2011) Deep sequencing-based transcriptome analysis of Plutella xylostella larvae parasitized by Diadegma semiclausum. BMC Genomics 12: 266–267. [PMC free article] [PubMed]
42. Wu SF, Sun FD, Qi YX, Yao Y, Fang Q, Huang J, et al. (2013) Parasitization by Cotesia chilonis influences gene expression in fatbody and hemocytes of Chilo suppressalis. PLoS ONE 8: e74309 doi: 10.1371/journal.pone.0074309 [PMC free article] [PubMed]
43. Zhu JY, Yang P, Zhang Z, Wu GX, Yang B (2013) Transcriptomic immune response of Tenebrio molitor pupae to parasitization by Scleroderma guani. PLoS ONE 8: e54411 doi: 10.1371/journal.pone.0054411 [PMC free article] [PubMed]
44. Tang B, Chen J, Hou Y, Meng E (2014) Transcriptome immune analysis of the invasive beetle Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae) parasitized by Tetrastichus brontispae Ferrière (Hymenoptera: Eulophidae). PLoS ONE 9:e91482 doi: 10.1371/journal.pone.0091482 [PMC free article] [PubMed]
45. Hoffmann JA (2003) The immune response of Drosophila. Nature 426:33–38. [PubMed]
46. Chu D, Zhang YJ, Cong B, Xu BY, Wu QJ, Zhu GR (2005) Sequence analysis of mtDNA COI gene and molecular phylogeny of different geographical populations of Bemisia tabaci (Gennadius). Scientia Agricultura Sinica 38: 76–85.
47. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. (2011) Fulllength transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29: 644–652. doi: 10.1038/nbt.1883 [PMC free article] [PubMed]
48. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M (2005) Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21: 3674–3676. [PubMed]
49. Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, et al. (2006) WEGO: a web tool for plotting GO annotations. Nucleic Acids Res 34: W293–W297. [PMC free article] [PubMed]
50. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5: 621–628. doi: 10.1038/nmeth.1226 [PubMed]
51. Hoffmann JA (2003) The immune response of Drosophila. Nature 426:33–38. [PubMed]
52. Hoffmann JA (1995) Innate immunity of insects. Curr Opin Immunol 7: 4–10. [PubMed]
53. Schmid-Hempel P (2005) Evolutionary ecology of insect immune defenses. Annu Rev Entomol 50: 529–551. [PubMed]
54. Feldhaar H, Gross R (2008) Immune reactions of insects on bacterial pathogens and mutualists. Microbes Infect 10: 1082–1088. doi: 10.1016/j.micinf.2008.07.010 [PubMed]
55. Nicolas E, Nappi AJ, Lemaitre B (1996) Expression of antimicrobial peptide genes after infection by parasitoid wasps in Drosophila. Dev Comp Immunol 20:175–181. [PubMed]
56. Boulanger N, Lowenberger C, Volf P, Ursic R, Sigutova L, Sabatiern L, et al. (2004) Characterization of a defensin from the sand fly Phlebotomus duboscqi induced by challenge with bacteria or the protozoan parasite Leishmania major. Infect Immun 72:7140–7146. [PMC free article] [PubMed]
57. Chiche L, Heitz A, Gelly J-C, Gracy J, Chau PT, Ha PT, et al. (2004) Squash inhibitors: from structural motifs to macrocyclic knottins. Curr Protein Pept Sci 5:341–349. [PubMed]
58. Christensen BM, Li J, Chen CC, Nappi AJ (2005): Melanization immune responses in mosquito vectors. Trends Parasitol 21:192–199. [PubMed]
59. Song KH, Jung MK, Eum JH, Hwang IC, Han SS (2008) Proteomic analysis of parasitized Plutella xylostella larvae plasma. J Insect Physiol 54:1271–1280. [PubMed]
60. Bangi E, Pitsouli C, Rahme LG, Cagan R, Apidianakis Y (2012) Immune response to bacteria induces dissemination of Ras-activated Drosophila hindgut cells. Embo Rep 13:569–576. doi: 10.1038/embor.2012.44 [PubMed]
61. Corona M, Robinson G (2006) Genes of the antioxidant system of the honey bee: annotation and phylogeny. Insect mol biol 15:687–701. [PMC free article] [PubMed]
62. Wang P, Oberley LW, Howe D, Jarvis DL, Chauhan G, Murhammer DW (2004) Effect of expression of manganese superoxide dismutase in baculovirus-infected insect cells. Appl Biochem Biotech 119:181–193. [PubMed]
63. Zhu JY, Ze SZ, Stanley DW, Yang B (2014) Parasitization by Scleroderma guani influences expression of superoxide dismutase genes in Tenebrio molitor. Arch Insect Biochem Physiol 87:40–52. doi: 10.1002/arch.21179 [PubMed]
64. Takeda T, Nakamatsu Y, Tanaka T (2006) Parasitization by Cotesia plutellae enhances detoxifying enzyme activity in Plutella xylostella. Pestic Biochem Physiol 86:15–22.
65. Parsell D, Lindquist S (1993) The function of heat-shock proteins in stress tolerance: degradation and reactivation of damaged proteins. Annu Rev Genet 27:437–496. [PubMed]
66. Qian M, Ren S, Hu Q (2005) Effects of endoparasitism by Encarsia bimaculata on the titers of juvenile hormones and 20—hydroxyecdysone in Bemisia tabaci nymphs. Acta entomologica Sinica 49:568–573.
67. Hu JS, Gelman DB, Blackburn MB (2002) Growth and development of Encarsia formosa (Hymenoptera: Aphelinidae) in the greenhouse whitefly, Trialeurodes vaporariorum (Homoptera: Aleyrodidae): effect of host age. Arch Insect Biochem Physiol 49:125–136. [PubMed]
68. Zhu JY, Ye GY, Hu C (2008) Advance in the studies on molecular mechanisms of parasitoid manipulation of host. Acta Phytophylacica Sinica 35:563–568.
69. Rath S, Sinha B (2005) Parasitization of fifth instar tasar silkworm, Antheraea mylitta, by the uzi fly, Blepharipa zebina; a host–parasitoid interaction and its effect on host’s nutritional parameters and parasitoid development. J Invertebr Pathol 88:70–78. [PubMed]
70. Kaeslin M, Pfister-Wilhelm R, Lanzrein B (2005) Influence of the parasitoid Chelonus inanitus and its polydnavirus on host nutritional physiology and implications for parasitoid development. J Insect Physiol 51:1330–1339. [PubMed]
71. Gilbert LI, Granger NA, Roe RM (2000) The juvenile hormones: historical facts and speculations on future research directions. Insect Biochem Mol Biol 30:617–644. [PubMed]
72. Anspaugh DD, Roe RM (2005) Regulation of JH epoxide hydrolase versus JH esterase activity in the cabbage looper, Trichoplusia ni, by juvenile hormone and xenobiotics. J Insect Physiol 51:523–535. [PubMed]
73. Ntambi JM, Miyazaki M, Dobrzyn A (2004) Regulation of stearoyl-CoA desaturase expression. Lipids 39:1061–1065. [PubMed]
74. Zhang R, Zhang YH, Shao D, Wang LD, Gong DQ (2013) The function and regulation of stearoyl-CoA desaturase gene. Chinese Bulletin of Life Sciences 4: 008.
75. Zhang SZ, Wan FH, Zhang F, Hua BZ (2003) Parasitic suitability of two strains of Encarsia formosa on Bemisia tabaci. Chinese J Biol Control 19: 149–153.
76. Alloatti A, Gupta S, Gualdrón-López M, Nguewa PA, Altabe SG, Deumer G, et al. (2011) Stearoyl-CoA desaturase is an essential enzyme for the parasitic protist Trypanosoma brucei. Biochem Bioph Res Co 412:286–290. [PubMed]
77. Holt HL, Aronstein KA, Grozinger CM (2013) Chronic parasitization by Nosema microsporidia causes global expression changes in core nutritional, metabolic and behavioral pathways in honey bee workers (Apis mellifera). BMC Genomics 14:799 doi: 10.1186/1471-2164-14-799 [PMC free article] [PubMed]
78. Dahlman DL, Vinson SB (1977) Effect of calyx fluid from an insect parasitoid on host hemolymph dry weight and trehalose content. J Invertebr Pathol 29:227–229.
79. Nakamatsu Y, Tanaka T (2004) Correlation between concentration of hemolymph nutrients and amount of fat body consumed in lightly and heavily parasitized hosts (Pseudaletia separata). J Insect Physiol 50:135–141. [PubMed]
80. Li LZ, Huang ZX (2001) Progresses in Cytochrome c Oxidase Studies. Chinese J Inorg Chem 17: 761–774.
81. Song LX, Wang RX (1989) Structure and function on ion transporting ATPase. Prog Physiol Sci 20: 334–338.
82. Visser B, Le Lann C, den Blanken FJ, Harvey JA, van Alphen JJ, Ellers J (2010) Loss of lipid synthesis as an evolutionary consequence of a parasitic lifestyle. PNAS 107: 8677–8682. doi: 10.1073/pnas.1001744107 [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science