Second-Generation Sequencing Supply an Effective Way to Screen RNAi Targets in Large Scale for Potential Application in Pest Insect Control

Abstract

The key of RNAi approach success for potential insect pest control is mainly dependent on careful target selection and a convenient delivery system. We adopted second-generation sequencing technology to screen RNAi targets. Illumina's RNA-seq and digital gene expression tag profile (DGE-tag) technologies were used to screen optimal RNAi targets from Ostrinia furnalalis. Total 14690 stage specific genes were obtained which can be considered as potential targets, and 47 were confirmed by qRT-PCR. Ten larval stage specific expression genes were selected for RNAi test. When 50 ng/µl dsRNAs of the genes DS10 and DS28 were directly sprayed on the newly hatched larvae which placed on the filter paper, the larval mortalities were around 40∼50%, while the dsRNAs of ten genes were sprayed on the larvae along with artificial diet, the mortalities reached 73% to 100% at 5 d after treatment. The qRT-PCR analysis verified the correlation between larval mortality and the down-regulation of the target gene expression. Topically applied fluorescent dsRNA confirmed that dsRNA did penetrate the body wall and circulate in the body cavity. It seems likely that the combination of DGE-tag with RNA-seq is a rapid, high-throughput, cost less and an easy way to select the candidate target genes for RNAi. More importantly, it demonstrated that dsRNAs are able to penetrate the integument and cause larval developmental stunt and/or death in a lepidopteron insect. This finding largely broadens the target selection for RNAi from just gut-specific genes to the targets in whole insects and may lead to new strategies for designing RNAi-based technology against insect damage.

About the Authors:

Yubing Wang

Contributed equally to this work with: Yubing Wang, Hao Zhang

Affiliation: Key Laboratory of Insect Developmental and Evolutionary Biology, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, People's Republic of China

Hao Zhang

Contributed equally to this work with: Yubing Wang, Hao Zhang

Affiliation: Key Laboratory of Insect Developmental and Evolutionary Biology, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, People's Republic of China

Haichao Li

Affiliation: Key Laboratory of Insect Developmental and Evolutionary Biology, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, People's Republic of China

Xuexia Miao

* E-mail: [email protected]

Affiliation: Key Laboratory of Insect Developmental and Evolutionary Biology, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, People's Republic of China

Introduction

RNA interference (RNAi) refers to double-stranded RNA (dsRNA) mediated gene silencing [1]. Since its discovery, it has been widely used in insect genetic research [2], [3], [4]. Recently, a new hot point is to find a feasible way to use RNAi as an alternative method for practical application of crop protection [5], [6], [7]. In this aspect, two major technologies are considered as the most important, namely, a simple delivery system and a high throughput screen method of optimal RNAi target genes, especially in those pests lacking any genome or transcriptome sequence information.

The dsRNA is delivered either by injection or by feeding insects. It has been reported that in several insect orders the gene can be knocked down by dsRNA injection [8], [9], [10], [11], but injection is not applicable to control pest insects in the field. Oral delivery of dsRNA was demonstrated in the nematode Caenorhabditis elegans [12], [13]. Since then, it was reported in a few insect species that dsRNA can be integrated in artificial diets and results in target gene knockdown [14], [15]. In addition, two research papers indicate that transgenic crops can produce dsRNAs to protect the plant against insect feeding damage [16], [17]. However, this technique is relatively complex and premature in terms of practical application in pest insect control. Thus, a more convenient way should be explored before dsRNA can be used as insecticides.

Regarding the optimal target gene selection, taking advantage of the known genes and cDNA library screening were the main ways up to now [14], [15], [17], [18]. Herein, the selection of target genes from the considered pest insects is still a major challenge. A high throughput screening system of RNAi target genes will be more appreciated when RNAi-based pest insect control can be taken in consideration. Recent results suggest that the next-generation sequencing (NGS) technologies allow directly sequencing the cDNA generated from messenger RNA (RNA-seq) [19], [20]. This new method for analyzing RNA-seq data enables the de novo reconstruction of the transcriptome for a nonmodel organism [21]. NGS technologies have led to novel opportunities for expression profiling in organisms lacking any genome or transcriptome sequence information [22], [23]. Can these technologies be used in RNAi target genes screening?

Asian corn borer Ostrinia furnalalis (Guenée) (ACB) is a major pest of corn in eastern and south-eastern Asia. At present, its control relies mainly on insecticides. Although the transgenic crops producing Bacillus thuringiensis (Bt) toxins for the pest insect control have been succeeful, their efficacy was gradually reduced when pests evolved resistance [24], [25], [26]. Thus there is an urgent need to find an alternative solution to control pest insects. The RNAi-based technology may be considered as an ideal alternate measure when efficient systems of dsRNA formulation and delivery have been developed [5], [7].

In this study we report that RNAi-induced lethality could be achieved by direct delivery via the exoskeleton; moreover an effective high-throughput NGS method was used to select stage-specific RNAi target genes in the Asian corn borer.

Results

Transcriptome of the ACB by RNA-seq

To obtain a global view of transcriptom in the whole ACB life cycle, high-throughput RNA-seq was performed by Solexa/Illumina sequencing technology in Beijing Genomics Institute at Shenzhen. We obtained more than 33 million paired-end reads of 75 bp in length. The total lengths of the reads were over 2.5 giganbases (Gb). The clean reads were de novo assembled into 124,034 contigs with the assemble software of SOAPdenovo [27] (Table 1). Then all the reads were realigned onto contigs again according to the paired-end reads overlap relationship, and all these contigs were joined into 79,825 scaffold sequences. Finally, the intra-scaffold gaps were filled using the paired-end extracted reads. Those sequences which cannot extend to any direction were called unigene. We obtained 45,750 unigenes for the ACB in this project (Table 2). The unigenes were annotated with the database of Nr, Swiss-Prot, KEGG, COG and Gene Ontology (GO) function analysis. The residual unigenes were predicated using the software of ESTscan [28].

[Figure omitted. See PDF.]

Table 1. The contigs assemble results of O. furnacalis transcriptome.

https://doi.org/10.1371/journal.pone.0018644.t001

[Figure omitted. See PDF.]

Table 2. The scaffold and unigene assemble result of O. furnacalis transcriptome.

https://doi.org/10.1371/journal.pone.0018644.t002

The analysis and statistics of DGE-tag sequences

We sequenced four DGE-tag libraries from four ACB developmental stages of egg, larva, pupa and adult. After eliminating low quality tags (containing Ns), copy numbers less than 2 and adaptor sequences, the numbers of total clean tags were about 3.57, 3.54, 3.55 and 3.45 millions for eggs, larvae, pupae and adults, respectively. The four successive developmental stages possessed 71004, 71228, 89637 and 89116 clean tag sequences, respectively. These tag sequences were composed of the DGE-tag datasets presented in Table 3. When the clean tags were annotated with the ACB RNA-seq dataset, 8415, 7988, 9123 and 10253 annotation unigenes were obtained for eggs, larvae, pupae and adults, respectively. Among them, 14690 stage specific unigene tags are unambiguous and can be considered as potential targets (Table S1). Among all of these unigenes, the differential expression and co-expression of clean tags and unigenes between 2 to 4 developmental stages are shown in Figure 1A and 1B.

[Figure omitted. See PDF.]

Figure 1. Venn diagram showing unique and shared gene expression among the ACB.

(A) All clean tag differential expression and co-expression between 2 to 4 developmental stages in four developmental stages (Egg, Larva, Pupa and Adult); (B) These tag sequences were mapping by RNA-seq dataset. We obtained 8415, 7988, 9123 and 10253 annotation unigenes for egg, larvae, pupal and adult stages, respectively. Among all of these unigenes, the differential expression and co-expression clean tags and unigenes between 2 to 4 developmental stages were shown.

https://doi.org/10.1371/journal.pone.0018644.g001

[Figure omitted. See PDF.]

Table 3. DGE-tag raw data and clean tag amount in the four developmental stages of O. furnacalis.

https://doi.org/10.1371/journal.pone.0018644.t003

To obtain the profile of the differential expression genes, the distribution of tag copy numbers in each of the four stages, is summarized in Table 4.

[Figure omitted. See PDF.]

Table 4. The distribution of clean tag copy number in the four developmental stages of O. furnacalis.

https://doi.org/10.1371/journal.pone.0018644.t004

Analysis of 10 larval stage specific genes and qRT-PCR

At first, we selected 47 larval stage high expression tags to determine their expression levels by qRT-PCR. All of them indeed appear consistent with their DGE-tag copy numbers (Table S2). To screen target genes used for producing RNAi of ACB, 10 larval genes with higher expression levels were selected according to the DGE-tag copy number. The cDNA sequences of these tags and their predicted functions were acquired from the RNA-seq dataset and Genebank (Table 5). The qRT-PCR analysis confirmed that the relative expression levels of all 10 genes in the larval stage were higher than those in other developmental stages except DS2 and DS3 (Figure 2A). Results are consistent with their tag copy numbers as seen in Table 5.

[Figure omitted. See PDF.]

Figure 2. The relative expression level of 10 larval stage specific genes and the larval mortality.

(A) The relative expression levels of 10 genes by qRT-PCR. All 10 genes were showing larval stage specific expression genes, and each of them was selected according to the DGE-tag copy number (Table 5). The expression profile was examined from egg to adult stage in which the larval stage involved 1st to 5th instar larvae and pre-pupa. The results were consistent with their DGE-tag copy numbers as shown in Table 5, indicating all of them were larval stage specific except DS2 and DS3 which had relatively higher expression level at pupal and adult stages. It was also consistent with the tag copy number (Table 5). To normalize threshold cycle (Ct) values obtained for each gene, 18S rRNA expression levels were used. (B) The larval mortality of dsRNA directly sprayed on the larvae along with artificial diet. Ten dsRNAs, 50 ng/µl for each (Table 5) were directly sprayed on the newly hatched larvae on artificial diet respectively. dsEYFP and the blank buffer were used as controls. The larval survival was monitored in a 5 d period. At day 5 the corrected mortality reached to more than 70% except DS5. (The data were also shown in Table 6). (C) The larval mortality of dsRNA directly sprayed on the larvae which were placed on filter paper. 50 ng/µl dsRNA of DS10 and DS28 were directly sprayed on the newly hatched larvae on a filter paper pre-soaked in the buffer, and then the larvae were transferred to artificial diet. The blank buffer and dsEYFP were used as controls. The larval survival was monitored in a 5 d period. (D) The 10 gene expression levels at different time after dsRNA treatment. The expression level of 10 genes increased at 1 d after dsRNA spraying (the red bars). 8 out of the 10 genes began to decrease at 3 d after treatment except DS2 and DS34 (the green bars). However, the expression levels of all 10 genes decreased in different degree 5 d after treatment (the purple bars).

https://doi.org/10.1371/journal.pone.0018644.g002

[Figure omitted. See PDF.]

Table 5. Ten larval stage specific expression gene tag sequences, copy number in each stage and their function prediction.

https://doi.org/10.1371/journal.pone.0018644.t005

[Figure omitted. See PDF.]

Table 6. The eggs hatch ratio and larval mortality (%) after spraying 50 ng/µl dsRNA.

https://doi.org/10.1371/journal.pone.0018644.t006

RNAi in ACB via a direct spray of dsRNA

To explore a convenient way of dsRNA delivery, 10 in vitro synthesized dsRNAs (Table 5) were sprayed uniformly on the newly hatched larvae along with an artificial diet. The blank buffer, dsEYFP-containing buffer and artificial diet were used as controls. When 50 ng/µl was sprayed, 9 out of 10 dsRNAs caused a significant developmental stunting and a high level of mortality. After correction, the 5 day mortalities levels were between 73% and 100%. Actually dsRNAs targeting DS2, DS10, DS28 and DS35 led to more than 90% death (Figure 2B; Table 6). As in this experiment dsRNAs were sprayed on the larvae along with artificial diet, there might be two possible pathways of dsRNA delivery, either by penetration of the body wall or by feeding. To test the first possibility (through the body wall), we sprayed 50 ng/µl dsRNAs targeting DS10 and DS28 to the newly hatched larvae placed on moist filter paper and then the larvae were transferred to artificial diet. As a result, 5 day corrected mortalities levels of 40% and 52% for DS10 and DS28 were obtained respectively (Figure 2C), LC50 were 61.02 (53.67–69.14, P<0.05) ng/µl for DS10 and 60.44 ng/µl (51.84–87.13, P<0.05) for DS28. This indicates that dsRNA can be delivered into haemolymph through the exoskeleton and still causes insect mortality. Another evidence is produced using labeled fluorescent dsRNA (see below).

To investigate the correlation of larval mortality with the target gene expression, the expression levels of all 10 genes were measured by qRT-PCR after spraying 50 ng/µl dsRNA. The expression levels of the 10 genes increased 1 day after dsRNA treatment (Figure 2D, red bar). Three days later, profiles began to decrease except that of DS2 and DS34 which showed a slight increase (Figure 2D, green bar). All gene expression levels were knocked down to a certain degree 5 days after dsRNA treatment (Figure 2D, purple bar).

When 20 ng/µl dsRNA solution was sprayed on the larvae along with the artificial diet, the larval development was remained in the first instar until 7 days after treatment and a certain degree of larval death was seen when the control larvae had reached the third instar stage (Figure 3A). When ACB eggs were soaked in 50 ng/µl dsRNA solution for 2 hours, all 10 genes caused an obvious decrease of hatching ratio compared with controls, although results varied for different genes (Table 6; Figure 3B).

[Figure omitted. See PDF.]

Figure 3. The eggs were soaked in dsRNA solutions or the dsRNA was sprayed on larvae.

(A) When 20 ng/µl dsRNA (DS3) feeding solution was sprayed on the larvae, the development of the treated larvae was stunted remarkably. CK (left side): one of the control larvae reaching third instar 7 d after spray with dsRNA blank buffer, DS3 (right side): the larvae remaining in first instar 7 d after spraying dsRNA of DS3 gene; (B) The egg masses were soaked in 50 ng/µl dsRNA solution for 2h, the embryos could develop, but most of them were not able to hatch.

https://doi.org/10.1371/journal.pone.0018644.g003

Use of fluorescent labeled dsRNA

To further confirm that dsRNA can penetrate the egg shell and the larval body wall of ACB, a fluorescent nucleotide was added during the dsRNA synthesizing of gene DS12. When 0.5 µl of 5 µg/µl fluorescent labeled dsRNA was topically applied on the larval protergum of the matured last instar larvae, fluorescence showed that penetration through the larval cuticle took only a few minutes and in 4 hours the whole body was reached (Figure 4, A–D). After 2 hours soaking the eggs in the fluorescent labeled dsRNA solution in a concentration of 50 ng/µl, the fluorescence was observed in the newly hatched larvae (Figure 4E). The gut (Figure 4F), hemocytes (Figure 4G) and silk fiber (Figure 4H) were also clearly labeled, suggesting that the delivered dsRNA not only penetrated into the body cavity but also circulated in heamolymph. To avoid the possibility of the droplet placed on the prothorax moved to the head and taken up orally by the larvae, we dropped 0.5 µl of 5 µg/µl fluorescent labeled dsRNA on the fifth abdominal tergum of the mature larvae, and a similar pattern was seen (figures not shown).

[Figure omitted. See PDF.]

Figure 4. The fluorescent labeled dsRNA dropping test.

The fluorescent labeled dsRNA (0.5 µl of 5 µg/µl) was topically applied on the protergum of the ACB mature larva. (A) The droplet was seen at the surface of the body wall just after dropping; (B) 30 min later, the droplet had already penetrated the body wall, and spread to the head and other thoracic segments; (C) 2 h later, the fluorescent dsRNA had spread to abdominal segments; (D) 4 h later, the fluorescence was seen in the whole body; (E) When the egg was soaked in the solution with 50 ng/µl fluorescent labeled dsRNA for 2 h, the newly hatched larvae were labeled with fluorescent; The fluorescence was observed in the gut(F), hemocytes (G) and silk fiber (H). From A to H, the left panel of each pair was on bright field (BF) and the right panel of each pair was on 559nm (Cy3). Image magnification: A–D, F, H (10×); E (40×); G (400×).

https://doi.org/10.1371/journal.pone.0018644.g004

Discussion

Taking advantage of known genes is a simple and effective way to verify RNAi targets, but the scope of selection was limited [15], . cDNA library screening approach can identify targets that would not necessarily be predicted from functional consideration but has the drawback of being rather labour-intensive comsuming if a large number of insect bioassays are required [14], [17]. RNA-seq is a recently developed deep-sequencing technology. After comparison with the existing technologies (DNA microarrays, cDNA or EST sequencing), Wang et al. (2009) summarized the advantages of RNA-seq [19]. This sequencing-based method does not include any cloning or amplification step, and requires less RNA sample. Furthermore, RNA-seq is not limited to detecting transcripts corresponding to existing genomic sequences, thus it is especially suitable for non-model organisms [19], [20]. The 3′-tag digital gene expression (DGE-tag) is also a second-generation sequencing-based technology. Asmann et al. (2009) comparing DGE-tag with the microarray technology, emphasized that this method has a high sensitivity and reproducibility for transcriptome profiling [23]. Therefore, we chose this method to screen the target genes for RNAi test in this study. Using this method, among the 14690 stage specific unigene tags, 10 larval stage specific genes were tested, 9 of them were confirmed as effective RNAi targets. It indicated that the combination of DGE-tag with RNA-seq is a rapid, cost less and easy way to select candidate target genes for RNAi in comparison with cDNA library screening.

Successful and unsuccessful RNAi experiments have been reviewed in a number of lepidopteran species by Terenius et al. (2011). As indicated in this review, the majority of RNAi studies in Lpidoptera were done using Bombyx mori, Manduca sexta and some Noctuidae species. Systemic RNAi has been demonstrated in some species, such as Hyalophora cecropia and B. mori, in which injection of dsRNA into the pupa can result in phenotypic effects in developing embryos, indicating dsRNA uptake by the developing oocytes of the pupa. However, a great variation of sensitivity to systemic RNAi has been seen among different lepidopteran species and high or no silencing can occur at very different concentrations of dsRNA [31]. Therefore, more work using different species and alternative delivery protocols seems helpful. It is generally thought that either injection or feeding are feasible ways to deliver dsRNA to produce a RNAi effect and potentially cause the development delay and death in insects [6], [18], [29]. Our results revealed that a direct spray of dsRNA can produce such effect in the ACB larvae, including the down-regulation of dsRNA targeted gene expression, developmental delay and death (Figures 2B and 2C; Table 6). In one experiment, after spraying DS10 and DS28 dsRNA on the newly hatched larvae placed on moist filter paper, their development was stunted and they finally died after 5 days of treatment. Corrected mortalities were 40% for DS10 and 52% for DS28. This indicates that the sprayed dsRNAs did penetrate the exoskeleton and could reach hemolymph and produce RNAi effect. In the previous reports no one had tried this kind of method for crop pest insects as it was believed that dsRNA could only go through insect midgut, the only region not covered by chitin exoskeleton. Huvenne and Smagghe deduced that sufficiently protected dsRNA, e.g. by some sort of coating, should be used to spray [7]. Results of the topically application of the fluorescent labeled DS12 dsRNA (Figure 4) further demonstrated that dsRNA does penetrate the integument, enters into heamolymph and labeled several tissues such as midgut, haemocytes and silk fiber. DS12 encoding protein was identified as a lipid requiring mitochondrial enzyme [32]. This suggests that the introduced dsRNA may circulate in heamolymph of ACB larvae. We still do not know whether the dsRNA is able to directly penetrate the body wall from the integumental membrane or reach the site of action via the integument of the tracheal system. However, this finding provides new insights that RNAi targets can be selected from the whole insect instead of gut-specific targets after feeding with dsRNA. So far, only one report had provided evidence that topically applied AaeIAP1 dsRNA products in acetone dilution was able to kill adult female mosquito Aedes aegypti [33]. In our experiments, two dsRNAs were directly spraying on larvae placed on filter paper. They match the genes encoding chymotrypsin-like serine protease C3 (DS10) and an unknown protein (DS28) (Table 5). Chymotrypsin-like serine protease was reported to be involved in food digestion in some lepidopteron insects [34] and interacted with the chitin synthesis from the midgut of Maduca sexta [35]. Recently a trypsin-like protease gene has been cloned in the integument of Helicopverpa armigera [36]. Thus, it is possible that this family gene encoding proteins may also interact with the chitin synthesis of the cuticle in dsRNA treated ACB larvae. In another experiment, we sprayed 10 larval stage specific genes dsRNA on the newly hatched larvae along with artificial diet. As a result, 9 of them caused a significantly developmental stunting and 5 days corrected mortalities reached high levels (between 73% and 100%) (Figure 2B; Table 5). And even four dsRNAs resulted in more than 90% mortality; matching genes of DS2, DS10, DS28 and DS35, which encode LIM protein, chymotrypsin-like serine protease C3, unknown protein and caboxypeptidase 4. The partial function of DS35 and DS10 is associated with the digestion in the midgut [37], [38], but LIM protein (DS2) is localized in muscle and was suggested to act as adapters capable of mediating different protein-protein interactions [39], [40]. DS35 was detected in the molting fluid possibly involved in the degradation of proteins from the old cuticle in B. mori [41]. It indicates that the observed RNAi effects can be achieved through a variety of targets in insects. The qRT-PCR analysis verified the correlation between larval mortality and the down-regulation of the target gene expression (Figure 2C). It is reasonable to think that dsRNA delivery through both the integument contact and oral ingestion may enhance RNAi effect leading to larval death. Because nine of the ten genes tested had a measureable impact on larval mortality, in order to exclude non-specific effects, we tested two more ACB genes, namely a pupal specific gene encoding the transcription elongation factor B polypeptide 3 binding protein (P18) and an adult gene encoding the hypothetical protein (A4). They show high expression levels only in pupae (P18) and adults (A4), respectively, but no expression in other stages. After dsRNAs of these two genes (50 ng/µl) were sprayed on the newly hatched larvae along with artificial diet, the larval mortality caused was similar to that of DS5, the only less susceptible gene among ten genes tested in terms of the RNAi induced larval death (data not shown). It suggests that dsRNAs of the higher expressed larval genes may achieve more effective gene silencing and higher mortality of ACB larvae.

In conclusion, our results demonstrate that a direct delivery via integument is a feasible way to introduce dsRNA into the haemolymph, which then leads to a down-regulation of gene expressions and finally to the development retarding and/or death. However, some questions still remain to be further investigated, such as in which tissues or cells these genes are silenced and whether the proposed receptor mediated endocytosis or the transmembrane channel-mediated uptake are the mechanisms leading to the production of RNAi effect as well as the pathway of RNAi effect leading to larval developmental stunt and death in ACB larvae. Nevertheless, our findings considerably broaden the selection of targets for RNAi research from gut-specific to whole insects. It may lead to new strategies in designing the RNAi-based technology against insect damage.

Materials and Methods

Insect culture

Asian corn borer Ostrinia furnacalis (Guenée) eggs were obtained from Sun Yat-Sen University and reared in laboratory at 25°C and 75% relative humidity on a 14-h-day/10-h-night cycle. Larvae were fed on a modified artificial diet (maize granule 120.0 g, maize flour 32.0 g, soybean flour 120.0 g, Vc 4.0 g, agar-agar 12.0 g, yeast powder 72.0 g, sorbic acid 4.0 g, glucose 60.0 g, gormaldehyde 1.6 mL, water 1000 ml).

Sample collection and RNA isolation

24 samples were collected from egg to adult stages during the whole life cycle of the ACB. The samples were immediately frozen in liquid nitrogen, and stored at −80°C before RNA extraction. Total RNA was isolated using a Qiagen RNA Extraction kit according to the manufacturer's instructions. It was treated with RNase-free DNase I for 30 min at 37°C (New England BioLabs) to remove residual DNA. Then the equivalent amount of the 24 samples was merged into 4 pools of egg, larva, pupa and adult. The total RNA in the whole life cycle came from the 4 pools in equivalent amounts. mRNA was isolated from DNA-free total RNA using Dynabeads mRNA Purification Kit (Invitrogen).

cDNA synthesis, sequencing and data analysis of RNA-seq

First-strand cDNA was synthesized using Oligo(dT)18 primer and reverse transcriptase (Invitrogen). The second-strand cDNA was synthesized using RNase H (Invitrogen) and DNA polymerase I (New England BioLabs). Then the 100–500 bp paired-end cDNA library was constructed for sequencing using the Illumina/Solexa technology.

Short read sequences were de novo assembled using SOAPdenovo program [27]. The unigene expression level was calculated by RPKM (Reads Per kb per Million reads) method [42]. Short clean reads were loaded into computer and the unambiguous sequence fragments as contigs were output. Then the reads were realigned onto the contigs according to the paired-end information and join the unique contigs into scaffolds. Finally, the intrascaffold gaps were filled using the paired-end extracted reads. Those sequences which can not extend to any direction were unigenes. The unigenes were compared with the database of Nr, Swiss-Prot, KEGG and COG and Gene Ontology (GO) function analysis.

Sequence tag preparation, sequencing and DGE-tag annotation

Sequence tag was prepared with Illumina's Digital Gene Expression Tag Profiling Kit according to the manufacturer's protocol. For the readers who want to have a schematic overview of the procedure, please see the reference [43]. When the DGE-tag library has been constructed, it was sequenced by Illumina Cluster Station and Genome Analyzer.

We filtered out low quality tags (containing Ns), copy number below 2 and adaptor sequences. At last, we obtained 3.4∼3.5 million clean sequence DGE-tags for each developmental stages of egg, larva, pupa and adult. The DGE tags, which consist of the CATG restriction enzyme digested site and an additional 17 bp from each transcript, were aligned onto the transcriptome unigenes of ACB. The expression abundance of transcripts was measured by the number of tags mapped.

dsRNA synthesis

dsRNAs were synthesized using the MEGAscript RNAi kit (Ambion, Huntingdon, UK) according to manufacturer's instructions. T7 promoter sequences were tailed to each 5′ end of DNA template by PCR amplifications. Double stranded EYFP (dsEYFP) was generated using pPigbacA3EYFP as template and was used as a negative control in the experiments. All the primer sequences were listed in Table 7. Template DNA and Single-stranded RNA were removed from the transcription reaction by DNase and RNase treatment. dsRNA was purified using MEGAclearTM columns (Ambion) and eluted in 100 ul nuclease free water. dsRNA concentrations were measured using biophotometer (Eppendorf).

[Figure omitted. See PDF.]

Table 7. The in vitro dsRNA synthesize with T7 tail primers of 10 larval stage specific genes.

https://doi.org/10.1371/journal.pone.0018644.t007

Fluorescent labeled dsRNA and the dropping test

The dsRNA synthesis procedure was the same as described above except that 0.25 µl 10 mM fluorescent dCTP labeled with Cy3 was added. Taking advantage of character, T7 RNA polymerase can take dNTP as substrates for synthesis RNA [44]. Then 0.5 µl of 5 µg/µl dsRNA solution was topically applied on the protergum or the fifth abdominal tergum of the ACB mature larvae that had stopped feeding and did not move. The larvae dripped with Cy3 labeled dsRNA were observed under common SZX16 fluorescent microscopy (Olympus) every 0.5 h or 1 h. The haemocytes were observed under confocal fluorescent microscopy FV1000 (Olympus).

Insect bioassay

Three hundred newly hatched larvae were placed in a culture dish (9 cm in diameter). 300 µl 50 ng/µl dsRNA in dsRNA buffer [14] of different genes was sprayed uniformly on them using the micro-sprayer. The experiments were performed in triplicate. After air dried the treated insects for 15 min, the 300 larvae were equally divided into 3 groups and placed in Petri dishes with artificial diet, then the dish was sealed with Mylar film (Dupont) that was pricked with small holes to provide some aeration without desiccating the food. The survival larvae were counted every day for each treatment. LC50 was calculated by using 5 different concentrations, 20, 40, 60, 80, 100 ng/ul dsRNA. The data were analyzed using one-way ANOVA to look for treatment effects compared to the untreated control (P<0.05). In addition, the egg masses were soaked in 50 ng/µl dsRNA soaking buffer for 2 hours and the effect on embryo development and the hatching rate was examined [45].

Quantitative real-time PCR (qRT-PCR)

Expression levels of 10 larval stage specific genes selected according to the DGE-tag copy number were evaluated by qRT-PCR. Primer sequences were provided in Table 8. A qRT-PCR assay for multiple genes was performed with the SYBR® Premix Ex Taq ™ II (TaKaRa). 18S rRNA expression levels were used to normalize Ct values obtained for each gene. qRT-PCR was carried out using a Mastercycler ep realplex instrument (Eppendorf). Data analysis methods were described as the reference [46].

[Figure omitted. See PDF.]

Table 8. The qRT-PCR primers of 10 larval stage specific genes.

https://doi.org/10.1371/journal.pone.0018644.t008

To assess the extent of RNAi, RNA was extracted from pools of 20∼30 dsRNA-treated and surviving larvae using a Qiagen RNA Extraction kit. Primer sequences were provided in Table 8.The samples were then treated with DNAse I (Invitrogen) to remove any genomic DNA contamination, and were used with Superscript II reverse transcriptase (Invitrogen) to make first strand cDNA using random primers. qRT-PCR reactions and statistic were the same with above.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

DGE-tag mapping by RNA-seq dataset and annotated.

https://doi.org/10.1371/journal.pone.0018644.s001

(XLS)

Table S2.

Larval stage high expression tags copy number confirmation by qRT-PCR.

https://doi.org/10.1371/journal.pone.0018644.s002

(XLS)

Acknowledgments

We sincerely thank Prof. Yongping Huang for his suggestion, Prof. Xiangxiong Zhu and Prof. Jean Marc Jallon for their valuable comment and revision of the manuscript.

Author Contributions

Conceived and designed the experiments: XXM. Performed the experiments: YBW HZ HCL. Analyzed the data: YBW HZ HCL. Contributed reagents/materials/analysis tools: YBW HZ HCL. Wrote the paper: XXM.

Citation: Wang Y, Zhang H, Li H, Miao X (2011) Second-Generation Sequencing Supply an Effective Way to Screen RNAi Targets in Large Scale for Potential Application in Pest Insect Control. PLoS ONE6(4): e18644. https://doi.org/10.1371/journal.pone.0018644

References

1. Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, et al. (1998) Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans. Nature 391: 806–811.A. FireS. XuMK MontgomerySA KostasSE Driver1998Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans.Nature391806811

2. Hannon GJ (2002) RNA interference. Nature 418: 244–251.GJ Hannon2002RNA interference.Nature418244251

3. Kuttenkeuler D, Boutros M (2004) Genome-wide RNAi as a route to gene function in Drosophila. Brief Funct Genomic Proteomic 3: 168–176.D. KuttenkeulerM. Boutros2004Genome-wide RNAi as a route to gene function in Drosophila.Brief Funct Genomic Proteomic3168176

4. Chen M, Du Q, Zhang HY, Wang X, Liang Z (2007) High-throughput screening using siRNA (RNAi) libraries. Expert Rev Mol Diagn 7: 281–291.M. ChenQ. DuHY ZhangX. WangZ. Liang2007High-throughput screening using siRNA (RNAi) libraries.Expert Rev Mol Diagn7281291

5. Belles X (2010) Beyond Drosophila: RNAi in vivo and functional genomics in insects. Annual Review of Entomology 55: 111–128.X. Belles2010Beyond Drosophila: RNAi in vivo and functional genomics in insects.Annual Review of Entomology55111128

6. Price DR, Gatehouse JA (2008) RNAi-mediated crop protection against insects. Trends Biotechnol 26: 393–400.DR PriceJA Gatehouse2008RNAi-mediated crop protection against insects.Trends Biotechnol26393400

7. Huvenne H, Smagghe G (2010) Mechanisms of dsRNA uptake in insects and potential of RNAi for pest control: a review. J Insect Physiol 56: 227–235.H. HuvenneG. Smagghe2010Mechanisms of dsRNA uptake in insects and potential of RNAi for pest control: a review.J Insect Physiol56227235

8. Tomoyasu Y, Denell RE (2004) Larval RNAi in Tribolium (Coleoptera) for analyzing adult development. Dev Genes Evol 214: 575–578.Y. TomoyasuRE Denell2004Larval RNAi in Tribolium (Coleoptera) for analyzing adult development.Dev Genes Evol214575578

9. Dong Y, Friedrich M (2005) Nymphal RNAi: systemic RNAi mediated gene knockdown in juvenile grasshopper. BMC Biotechnol 5: 25.Y. DongM. Friedrich2005Nymphal RNAi: systemic RNAi mediated gene knockdown in juvenile grasshopper.BMC Biotechnol525

10. Martin D, Maestro O, Cruz J, Mane-Padros D, Belles X (2006) RNAi studies reveal a conserved role for RXR in molting in the cockroach Blattella germanica. J Insect Physiol 52: 410–416.D. MartinO. MaestroJ. CruzD. Mane-PadrosX. Belles2006RNAi studies reveal a conserved role for RXR in molting in the cockroach Blattella germanica.J Insect Physiol52410416

11. Mutti NS, Park Y, Reese JC, Reeck GR (2006) RNAi knockdown of a salivary transcript leading to lethality in the pea aphid, Acyrthosiphon pisum. J Insect Sci 6: 1–7.NS MuttiY. ParkJC ReeseGR Reeck2006RNAi knockdown of a salivary transcript leading to lethality in the pea aphid, Acyrthosiphon pisum.J Insect Sci617

12. Timmons L, Fire A (1998) Specific interference by ingested dsRNA. Nature 395: 854.L. TimmonsA. Fire1998Specific interference by ingested dsRNA.Nature395854

13. Timmons L, Court DL, Fire A (2001) Ingestion of bacterially expressed dsRNAs can produce specific and potent genetic interference in Caenorhabditis elegans. Gene 263: 103–112.L. TimmonsDL CourtA. Fire2001Ingestion of bacterially expressed dsRNAs can produce specific and potent genetic interference in Caenorhabditis elegans.Gene263103112

14. Araujo RN, Santos A, Pinto FS, Gontijo NF, Lehane MJ, et al. (2006) RNA interference of the salivary gland nitrophorin 2 in the triatomine bug Rhodnius prolixus (Hemiptera: Reduviidae) by dsRNA ingestion or injection. Insect Biochem Mol Biol 36: 683–693.RN AraujoA. SantosFS PintoNF GontijoMJ Lehane2006RNA interference of the salivary gland nitrophorin 2 in the triatomine bug Rhodnius prolixus (Hemiptera: Reduviidae) by dsRNA ingestion or injection.Insect Biochem Mol Biol36683693

15. Turner CT, Davy MW, MacDiarmid RM, Plummer KM, Birch NP, et al. (2006) RNA interference in the light brown apple moth, Epiphyas postvittana (Walker) induced by double-stranded RNA feeding. Insect Mol Biol 15: 383–391.CT TurnerMW DavyRM MacDiarmidKM PlummerNP Birch2006RNA interference in the light brown apple moth, Epiphyas postvittana (Walker) induced by double-stranded RNA feeding.Insect Mol Biol15383391

16. Mao YB, Cai WJ, Wang JW, Hong GJ, Tao XY, et al. (2007) Silencing a cotton bollworm P450 monooxygenase gene by plant-mediated RNAi impairs larval tolerance of gossypol. Nature Biotechnology 25: 1307–1313.YB MaoWJ CaiJW WangGJ HongXY Tao2007Silencing a cotton bollworm P450 monooxygenase gene by plant-mediated RNAi impairs larval tolerance of gossypol.Nature Biotechnology2513071313

17. Baum JA, Bogaert T, Clinton W, Heck GR, Feldmann P, et al. (2007) Control of coleopteran insect pests through RNA interference. Nature Biotechnology 25: 1322–1326.JA BaumT. BogaertW. ClintonGR HeckP. Feldmann2007Control of coleopteran insect pests through RNA interference.Nature Biotechnology2513221326

18. Whyard S, Singh AD, Wong S (2009) Ingested double-stranded RNAs can act as species-specific insecticides. Insect Biochem Mol Biol 39: 824–832.S. WhyardAD SinghS. Wong2009Ingested double-stranded RNAs can act as species-specific insecticides.Insect Biochem Mol Biol39824832

19. Wang Z, Gerstein M, Snyder M (2009) RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 10: 57–63.Z. WangM. GersteinM. Snyder2009RNA-Seq: a revolutionary tool for transcriptomics.Nat Rev Genet105763

20. Haas BJ, Zody MC (2010) Advancing RNA-Seq analysis. Nature Biotechnology 28: 421–423.BJ HaasMC Zody2010Advancing RNA-Seq analysis.Nature Biotechnology28421423

21. Kahvejian A, Quackenbush J, Thompson JF (2008) What would you do if you could sequence everything? Nature Biotechnology 26: 1125–1133.A. KahvejianJ. QuackenbushJF Thompson2008What would you do if you could sequence everything?Nature Biotechnology2611251133

22. Morrissy AS, Morin RD, Delaney A, Zeng T, McDonald H, et al. (2009) Next-generation tag sequencing for cancer gene expression profiling. Genome Res 19: 1825–1835.AS MorrissyRD MorinA. DelaneyT. ZengH. McDonald2009Next-generation tag sequencing for cancer gene expression profiling.Genome Res1918251835

23. Asmann YW, Klee EW, Thompson EA, Perez EA, Middha S, et al. (2009) 3′ tag digital gene expression profiling of human brain and universal reference RNA using Illumina Genome Analyzer. BMC Genomics 10: 531.YW AsmannEW KleeEA ThompsonEA PerezS. Middha20093′ tag digital gene expression profiling of human brain and universal reference RNA using Illumina Genome Analyzer.BMC Genomics10531

24. Sivakumar S, Rajagopal R, Venkatesh GR, Srivastava A, Bhatnagar RK (2007) Knockdown of aminopeptidase-N from Helicoverpla armigera larvae and in transfected Sf21 cells by RNA interference reveals its functional interaction with Bacillus thuringiensis insecticidal protein Cry1Ac. Journal of Biological Chemistry 282: 7312–7319.S. SivakumarR. RajagopalGR VenkateshA. SrivastavaRK Bhatnagar2007Knockdown of aminopeptidase-N from Helicoverpla armigera larvae and in transfected Sf21 cells by RNA interference reveals its functional interaction with Bacillus thuringiensis insecticidal protein Cry1Ac.Journal of Biological Chemistry28273127319

25. Tabashnik BE, Gassmann AJ, Crowder DW, Carriere Y (2008) Insect resistance to Bt crops: evidence versus theory. Nature Biotechnology 26: 199–202.BE TabashnikAJ GassmannDW CrowderY. Carriere2008Insect resistance to Bt crops: evidence versus theory.Nature Biotechnology26199202

26. Tabashnik BE, Van Rensburg JB, Carriere Y (2009) Field-evolved insect resistance to Bt crops: definition, theory, and data. J Econ Entomol 102: 2011–2025.BE TabashnikJB Van RensburgY. Carriere2009Field-evolved insect resistance to Bt crops: definition, theory, and data.J Econ Entomol10220112025

27. Li R, Li Y, Kristiansen K, Wang J (2008) SOAP: short oligonucleotide alignment program. Bioinformatics 24: 713–714.R. LiY. LiK. KristiansenJ. Wang2008SOAP: short oligonucleotide alignment program.Bioinformatics24713714

28. Iseli C, Jongeneel CV, Bucher P (1999) ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol 138–148.C. IseliCV JongeneelP. Bucher1999ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences.Proc Int Conf Intell Syst Mol Biol138148

29. Walshe DP, Lehane SM, Lehane MJ, Haines LR (2009) Prolonged gene knockdown in the tsetse fly Glossina by feeding double stranded RNA. Insect Mol Biol 18: 11–19.DP WalsheSM LehaneMJ LehaneLR Haines2009Prolonged gene knockdown in the tsetse fly Glossina by feeding double stranded RNA.Insect Mol Biol181119

30. Bautista MA, Miyata T, Miura K, Tanaka T (2009) RNA interference-mediated knockdown of a cytochrome P450, CYP6BG1, from the diamondback moth, Plutella xylostella, reduces larval resistance to permethrin. Insect Biochem Mol Biol 39: 38–46.MA BautistaT. MiyataK. MiuraT. Tanaka2009RNA interference-mediated knockdown of a cytochrome P450, CYP6BG1, from the diamondback moth, Plutella xylostella, reduces larval resistance to permethrin.Insect Biochem Mol Biol393846

31. Terenius O, Papanicolaou A, Garbutt JS, Eleftherianos I, Huvenne H, et al. (2011) RNA interference in Lepidoptera: an overview of successful and unsuccessful studies and implications for experimental design. J Insect Physiol 57: 231–245.O. TereniusA. PapanicolaouJS GarbuttI. EleftherianosH. Huvenne2011RNA interference in Lepidoptera: an overview of successful and unsuccessful studies and implications for experimental design.J Insect Physiol57231245

32. Joanisse DR, Storey KB (1996) Fatty acid content and enzymes of fatty acid metabolism in overwintering cold-hardy gall insects. Physiological Zoology 69: 1079–1095.DR JoanisseKB Storey1996Fatty acid content and enzymes of fatty acid metabolism in overwintering cold-hardy gall insects.Physiological Zoology6910791095

33. Pridgeon JW, Zhao LM, Becnel JJ, Strickman DA, Clark GG, et al. (2008) Topically applied AaeIAP1 double-stranded RNA kills female adults of Aedes aegypti. Journal of Medical Entomology 45: 414–420.JW PridgeonLM ZhaoJJ BecnelDA StrickmanGG Clark2008Topically applied AaeIAP1 double-stranded RNA kills female adults of Aedes aegypti.Journal of Medical Entomology45414420

34. Coates BS, Hellmich RL, Lewis LC (2006) Sequence variation in trypsin- and chymotrypsin-like cDNAs from the midgut of Ostrinia nubilalis: methods for allelic differentiation of candidate Bacillus thuringiensis resistance genes. Insect Mol Biol 15: 13–24.BS CoatesRL HellmichLC Lewis2006Sequence variation in trypsin- and chymotrypsin-like cDNAs from the midgut of Ostrinia nubilalis: methods for allelic differentiation of candidate Bacillus thuringiensis resistance genes.Insect Mol Biol151324

35. Broehan G, Zimoch L, Wessels A, Ertas B, Merzendorfer H (2007) A chymotrypsin-like serine protease interacts with the chitin synthase from the midgut of the tobacco hornworm. J Exp Biol 210: 3636–3643.G. BroehanL. ZimochA. WesselsB. ErtasH. Merzendorfer2007A chymotrypsin-like serine protease interacts with the chitin synthase from the midgut of the tobacco hornworm.J Exp Biol21036363643

36. Liu Y, Sui YP, Wang JX, Zhao XF (2009) Characterization of the Trypsin-Like Protease (Ha-Tlp2) Constitutively Expressed in the Integument of the Cotton Bollworm, Helicoverpa Armigera. Archives of Insect Biochemistry and Physiology 72: 74–87.Y. LiuYP SuiJX WangXF Zhao2009Characterization of the Trypsin-Like Protease (Ha-Tlp2) Constitutively Expressed in the Integument of the Cotton Bollworm, Helicoverpa Armigera.Archives of Insect Biochemistry and Physiology727487

37. Hegedus D, Baldwin D, O'Grady M, Braun L, Gleddie S, et al. (2003) Midgut proteases from Mamestra configurata (Lepidoptera: Noctuidae) larvae: characterization, cDNA cloning, and expressed sequence tag analysis. Arch Insect Biochem Physiol 53: 30–47.D. HegedusD. BaldwinM. O'GradyL. BraunS. Gleddie2003Midgut proteases from Mamestra configurata (Lepidoptera: Noctuidae) larvae: characterization, cDNA cloning, and expressed sequence tag analysis.Arch Insect Biochem Physiol533047

38. Zhang C, Zhou D, Zheng S, Liu L, Tao S, et al. (2010) A chymotrypsin-like serine protease cDNA involved in food protein digestion in the common cutworm, Spodoptera litura: Cloning, characterization, developmental and induced expression patterns, and localization. J Insect Physiol 56: 788–799.C. ZhangD. ZhouS. ZhengL. LiuS. Tao2010A chymotrypsin-like serine protease cDNA involved in food protein digestion in the common cutworm, Spodoptera litura: Cloning, characterization, developmental and induced expression patterns, and localization.J Insect Physiol56788799

39. Pechanova O, Stone WD, Monroe W, Nebeker TE, Klepzig KD, et al. (2008) Global and comparative protein profiles of the pronotum of the southern pine beetle, Dendroctonus frontalis. Insect Mol Biol 17: 261–277.O. PechanovaWD StoneW. MonroeTE NebekerKD Klepzig2008Global and comparative protein profiles of the pronotum of the southern pine beetle, Dendroctonus frontalis.Insect Mol Biol17261277

40. Dawid IB (1998) LIM protein interactions: Drosophila enters the stage. Trends Genet 14: 480–482.IB Dawid1998LIM protein interactions: Drosophila enters the stage.Trends Genet14480482

41. Ote M, Mita K, Kawasaki H, Daimon T, Kobayashi M, et al. (2005) Identification of molting fluid carboxypeptidase A (MF-CPA) in Bombyx mori. Comp Biochem Physiol B Biochem Mol Biol 141: 314–322.M. OteK. MitaH. KawasakiT. DaimonM. Kobayashi2005Identification of molting fluid carboxypeptidase A (MF-CPA) in Bombyx mori.Comp Biochem Physiol B Biochem Mol Biol141314322

42. Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature Methods 5: 621–628.A. MortazaviBA WilliamsK. MccueL. SchaefferB. Wold2008Mapping and quantifying mammalian transcriptomes by RNA-Seq.Nature Methods5621628

43. t Hoen PA, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, et al. (2008) Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res 36: e141.PA t HoenY. AriyurekHH ThygesenE. VreugdenhilRH Vossen2008Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms.Nucleic Acids Res36e141

44. Gaur RK, Krupp G (1993) Enzymatic RNA synthesis with deoxynucleoside 5′-O-(1-thiotriphosphates). FEBS Lett 315: 56–60.RK GaurG. Krupp1993Enzymatic RNA synthesis with deoxynucleoside 5′-O-(1-thiotriphosphates).FEBS Lett3155660

45. Gallagher SR (1992) GUS protocols: using the GUS gene as a reporter of gene expression. Academic Press. SR Gallagher1992GUS protocols: using the GUS gene as a reporter of gene expression.Academic Press

46. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25: 402–408.KJ LivakTD Schmittgen2001Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method.Methods25402408