Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

RNA Interference Depletion of the Halloween Gene Disembodied Implies its Potential Application for Management of Planthopper Sogatella furcifera and Laodelphax striatellus

  • Pin-Jun Wan ,

    Contributed equally to this work with: Pin-Jun Wan, Shuang Jia

    Affiliation Education Ministry Key Laboratory of Integrated Management of Crop Diseases and Pests, College of Plant Protection, Nanjing Agricultural University, Nanjing, China

  • Shuang Jia ,

    Contributed equally to this work with: Pin-Jun Wan, Shuang Jia

    Affiliation Education Ministry Key Laboratory of Integrated Management of Crop Diseases and Pests, College of Plant Protection, Nanjing Agricultural University, Nanjing, China

  • Na Li,

    Affiliation Education Ministry Key Laboratory of Integrated Management of Crop Diseases and Pests, College of Plant Protection, Nanjing Agricultural University, Nanjing, China

  • Jin-Mei Fan,

    Affiliation Education Ministry Key Laboratory of Integrated Management of Crop Diseases and Pests, College of Plant Protection, Nanjing Agricultural University, Nanjing, China

  • Guo-Qing Li

    ligq@njau.edu.cn

    Affiliation Education Ministry Key Laboratory of Integrated Management of Crop Diseases and Pests, College of Plant Protection, Nanjing Agricultural University, Nanjing, China

Abstract

Sogatella furcifera and Laodelphax striatellus are economically important rice pests in China by acting as vectors of several rice viruses, sucking the phloem sap and blocking the phloem vessels. Ecdysteroid hormone 20-hydroxyecdysone regulates insect development and reproduction. A cytochrome P450 monooxygenase CYP302A1 (22-hydroxylase), encoded by the Halloween gene disembodied (dib), plays a critical role in ecdysteroidogenesis. The objective of this study is to test whether dib genes are potential targets for RNA interference-based management of S. furcifera and L. striatellus. We cloned and characterized Sfdib and Lsdib. The open reading frame regions of dib genes were generated and used for designing and constructing dsRNA fragments. Experiments were conducted using oral delivery of dsdib to investigate the effectiveness of RNAi in S. furcifera and L. striatellus nymphs. Real-time quantitative reverse transcriptase-PCR analysis demonstrated that continuous ingestion of dsdib at the concentration of 0.01, 0.05 and 0.50 mg/ml diminished Sfdib expression levels by 35.9%, 45.1% and 66.2%, and ecdysone receptor (SfEcR) gene mRNA levels by 34.0%, 36.2% and 58.5% respectively in S. furcifera, and decreased Lsdib expression level by 18.8%, 35.8% and 56.7%, and LsEcR mRNA levels by 25.2%, 46.8% and 68.8% respectively in L. striatellus. The reduction in dib and EcR transcript abundance resulted in observable phenotypes. The development of nymphs was impaired and the survival was negatively affected. Our data will enable the development of new insect control strategies and functional analysis of vital genes in S. furcifera and L. striatellus nymphs.

Introduction

The white-backed planthopper Sogatella furcifera and the small brown planthopper Laodelphax striatellus are economically important pests of rice in China. They cause serious damage to rice plants by acting as viral vectors, and by sucking the phloem sap and blocking the phloem vessels. In China, S. furcifera and L. striatellus are found in all rice-growing areas, especially in the middle and downstream Yangtze River. The most common management strategy to control these planthoppers is chemical treatments. However, this inevitably leads to the development of insecticide resistance, insect resurgence, and serious environmental pollution [1][5]. Alternative control strategies must be developed.

Among the potential control strategies, RNA interference (RNAi), an RNA-dependent gene silencing process induced by double-stranded RNA (dsRNA), has been proposed as a promising approach to control the damage caused by insect pests [6], [7]. In the brown planthopper Nilaparvata lugens, for example, dietary introduction of dsRNAs can not only effectively silence target genes that are mainly expressed in gut [8][10], or widely expressed in various tissues [11], [12], but also knock down target genes which are expressed in other tissues far from the midgut [13]. Similar results have been documented in S. furcifera and L. striatellus [14][16]. These results demonstrate that RNAi has great potential in planthopper management. Identification of the most attractive candidate genes for RNAi in S. furcifera and L. striatellus is the first step.

Ecdysteroid 20-hydroxyecdysone (20E) regulates insect development and reproduction. A cytochrome P450 monooxygenase CYP302A1 (22-hydroxylase), encoded by the Halloween gene disembodied (dib), has been confirmed to be involved in ecdysteroidogenesis in two dipteran species Drosophila melanogaster and Anopheles gambiae, and two lepidopteran species Bombyx mori and Manduca sexta, by expression of their corresponding genes in S2 cells [17][22]. Moreover, in D. melanogaster, several mutants in dib gene have been found [17]. Among them, dibF8 and dibP3 produce truncated proteins that lack the heme binding site. The dib13A allele is characterized by a premature stop codon that removes the last five amino acids [18]. Analysis by RP-HPLC/RIA demonstrate that the homozygous dib mutant embryos have very low titers of ecdysone (E) and 20E and fail to express IMP-E1 and L1, two 20E-inducible genes, in certain tissues of the embryo. Phenotypic analysis reveals that dib mutant embryos produce little or no cuticle, exhibit severe defects in many late morphogenetic processes such as head involution, dorsal closure and gut development, and finally die at the embryonic stage [18]. Furthermore, mostly based on sequence similarity, dib genes are described or predicted in hemipteran Acyrthosiphon pisum [23] and Bemisia tabaci [24], in dipteran Aedes aegypti [25], in lepidopteran Spodoptera littoralis [26], and in hymenopteran Apis mellifera [27]. In addition to insects, a dib gene has been identified in the water flea Daphnia pulex [28].

In the present paper, we cloned and characterized dib genes in S. furcifera and L. striatellus. We demonstrated the insecticidal action of orally-delivered dsRNA (dsdib) in S. furcifera and L. striatellus nymphs. We found that continuous ingestion of dsdib delayed nymphal growth and caused lethality in both planthopper species. Our data will enable the development of new insect control strategies and functional analysis of vital genes.

Materials and Methods

Insect Culture

S. furcifera and L. striatellus were routinely reared on the rice (Oryza sativa) variety Taichung Native 1, using a protocol described recently [14], [15]. In our laboratory, S. furcifera eggs hatched into nymphs within 7 days. Nymphs went through 5 instars, with the first-, second-, third-, fourth- and fifth-instars lasting 2.5, 2.0, 2.0, 3.0 and 3.0 days, respectively. L. striatellus eggs hatched into nymphs within 13 days. The first-, second-, third-, fourth- and fifth-instar nymphs spent an average of 6.0, 5.0, 5.0, 5.0 and 6.0 days, respectively.

Cloning and Sequencing of Full-length Lsdib cDNA

We obtained Sfdib and Lsdib unigenes from the corresponding transcriptome datasets of S. furcifera and L. striatellus recently constructed in our laboratory [14], [15]. The resulting sequences were authenticated by polymerase chain reaction (PCR) using specific primers (Table 1), the same thermal cycling conditions and PCR reaction system as described recently [14], [15]. The full-length Sfdib and Lsdib cDNAs were obtained by 5′- and 3′-RACE, using antisense and sense gene-specific primers (Table 1) and the universal primers in the SMART™ RACE kit, and the standard manufacturer-recommended components and thermal cycling conditions. The PCR products were then cloned into pGEM-T easy vector (Promega), and sequenced at both strands. After full-length cDNAs were obtained, two primer pairs were designed (Table 1) to verify the sequence. Open reading frames were predicted using the editseq program of DNAStar (http://www.dnastar.com). The nucleotide sequences were finally annotated by Dr. David Nelson in accordance with the CYP nomenclature committee convention (http://drnelson.uthsc.edu/cytochromeP450.html) [29]. The resulting sequences were submitted to GenBank (Sfdib, KC579455; Lsdib, KC579460).

thumbnail
Table 1. Primers used in RT-PCR, 5′ and 3′ RACE, synthesizing dsRNA, and performing qRT-PCR.

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

Multiple Sequence Alignment and Phylogenetic Analysis

The annotated Dib proteins from a spider species Tetranychus urticae and 17 insect species were aligned with the predicted SfDib and LsDib sequences using ClustalW2.1 [30]. The alignment was used to construct the maximum-likelihood (ML) trees using RAxML version 7.26 [31] to select the best-fitting model (LG+γ) after estimation by ProtTest [32]. The reliability of the ML tree topology was evaluated by bootstrapping a sampling of 1000 replicates.

Bioassays Using dsRNA, 20E and their Mixtures

DNA samples for the production of Sfdib, Lsdib and enhanced green fluorescent protein (egfp, control) were synthesized by PCR, using a 458 bp Sfdib cDNA, a 406 bp Lsdib cDNA and a 414 bp egfp fragments, and gene-specific primers (Table 1) incorporating the T7 RNA polymerase promoter sequence (5′-taatacgactcactataggg-3′) [14], [15]. PCR products were purified using the Wizard H SV Gel and PCR Clean-Up System (Promega, Madison, WI, USA) and used for dsRNA synthesis using the T7 Ribomax TM Express RNAi System (Promega, Madison, WI, USA). The synthesized dsRNAs were respectively isopropanol-precipitated, resuspended in nuclease-free water, quantified by a spectrophotometry (NanoDrop TM 1000, Thermo Fisher Scientific, USA) at 260 nm, and kept at −70°C until use.

A dietary dsRNA-introducing procedure previously reported [8], [13] was used to feed the insects, with glass cylinders (12 cm in length and 2.8 cm in internal diameter) as feeding chambers. Twenty S. furcifera first-instar nymphs or L. striatellus third-instar nymphs were carefully transferred into each chamber and pre-reared for one day to the next instar, on liquid artificial diet between two layers of stretched Parafilm M (Pechiney Plastic Packaging Company, Chicago, IL, USA) that were placed at both ends of the chamber. The artificial diet containing one of the dsRNAs was then used to feed the nymphs. The diet was changed and dead nymphs were removed daily.

Two bioassays were carried out for each planthopper species, with its corresponding dsdib (dsSfdib for S. furcifera nymphs, dsLsdib for L. striatellus nymphs). The first experiment had five treatments including non-dsRNA diet (blank control), the diet containing dsegfp at the concentration of 0.50 mg/ml (negative control), and the diets containing dsdib at the concentration of 0.05, 0.01 and 0.50 mg/ml. The second bioassay was a rescue experiment. 20E (Sigma, purified by reverse-phase HPLC before experiments) in ethanol was mixed with the diet to obtain a final concentration of 0.50 mg/ml diet, a dose that has no obvious negative effects on planthopper species [14], [15]. There were four treatments: non-dsRNA+ethanol diet (blank control), 0.50 mg/ml of dsegfp+ethanol diet (negative control), 0.50 mg/ml of dsdib+ethanol diet, and 0.50 mg/ml dsdib+0.5 mg/ml of 20E diet. All treatments in the experiments were replicated 25 times (25 chambers), and a total of 250 nymphs in each treatment were used (100 nymphs for bioassays and 150 nymphs for q-PCR).

The experiments lasted 6 days for S. furcifera, and 5 days for L. striatellus nymphs. Mortality was recorded daily. After each experiment, the surviving nymphs were transferred to rice to measure the development delay. The instars were identified by the same method described recently [14], [15]. Mortalities from the termination of experiment to adults emergence were also calculated.

Real-time Quantitative PCR

For S. furcifera, total RNA samples were prepared from whole bodies of the second-instar (Day 1, 2 in the second instar, I2D1 and I2D2), third-instar (I3D1 and I3D2) and fourth-instar (I4D1, I4D2 and I4D3) nymphs on rice, from the head, thorax and abdomen of I4D3 nymphs, and from nymphs subjected to the 6-day bioassays. For L. striatellus, total RNA samples were prepared from I4D1 to I4D5, and I5D1 nymphs on rice, from the head, thorax and abdomen of I4D5 nymphs, and from whole bodies of nymphs which had been exposed to dsRNA for 5 days. Each sample contained 5–10 nymphs and repeated in biological triplicate. All RNAs were isolated using SV Total RNA Isolation System Kit (Promega). The expression levels of Sfdib and SfEcR in each S. furcifera nymphal sample were measured by qRT-PCR as described proviously [14], [15]. The S. furcifera ADP-ribosylation factor (ARF) and ribosomal protein RPL9 were used as internal controls. The mRNA abundance of Lsdib and LsEcR in each L. striatellus nymphal sample was measured using the internal control genes ADP-ribosylation factor (ARF) and elongation factor-1 (EF-1) (primers for these genes are listed in Table 1). Each sample was repeated in technical triplicate. Data were analyzed by the 2−ΔΔCt method [33], using the geometric mean of the two internal control genes for normalization according to the strategy described previously [33], [34].

Data Analysis

The data were given as means ± SE, and were analyzed by ANOVA or a repeated measures ANOVA followed by the Tukey-Kramer test, using SPSS for Windows (SPSS, Chicago, IL, USA).

Results

Molecular Cloning and Sequence Analysis

Two cDNAs of the putative Halloween gene dib (disembodied, cyp302a1) were obtained. Sfdib from S. furcifera contained a complete coding sequence and encoded a 547-amino-acid protein. Similarly, Lsdib from L. striatellus had a complete open reading frame which encoded for a putative 22-hydroxylase of 550 amino acid residues (Fig. 1).

thumbnail
Figure 1. Amino acid sequence alignment of Dib-like proteins.

ClustalX (2.1) multiple sequence alignment program were used. CYP302A1 (Disembodied, Dib) originates from Spodoptera littoralis (SlDib, ACM46003), Anopheles gambiae (AgDib, NP_524810), Nasonia vitripennis (NvDib, XP_001601675), Sogatella furcifera (SfDib) and Laodelphax striatellus (LsDib) respectively. Amino acids with 100%, >80%, and >60% conservation are shaded in black, dark grey and light grey, respectively. The characteristic P450 structure, a mitochondrial targeting segment, Helix C, Helix I, Helix K, PERF motif and Heme-binding domain are indicated above the alignment. Two positively charged residues, a signature for mitochondrial enzyme, are indicated with triangles.

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

The deduced amino acid sequences of SfDib and LsDib share 95% identity. Moreover, SfDib and LsDib have 30% and 30%, 57% and 55%, 56% and 55%, and 54% and 54% of identities with the homologues from A. pisum, Tribolium castaneum, M. sexta and D. melanogaster. Both SfDib and LsDib have five conserved motifs typical of insect P450s, i.e., WxxxR (Helix-C), GxE/DTT/S (Helix-I), ExxR (Helix-K), PxxFxPE/DRF (PERF motif) and PFxxGxRxCxG/A (heme-binding domain), where ‘x’ means any amino acid [35], [36] (Fig. 1).

Insect Dib belongs to the mitochondrial P450 family. The N-termini of SfDib and LsDib have a typical mitochondrial targeting sequence, consisting of several charged residues. Moreover, both SfDib and LsDib possess two positive charged residues near heme-binding domain, another character of typical mitochondrial P450 enzymes (Fig. 1).

The evolutionary relationship of Dib-like proteins derived from a spider mite species and 19 insect species was evaluated (Fig. 2). Those from 19 insects were distant from that came from the spider mite Tetranychus urticae. In 19 insect species, the Dib-like proteins formed a Hymenoptera clade, a Diptera clade and a Lepidoptera clade. Among the Hymenopteran clade, three sub-clades were formed by proteins derived from Apidae, Formicidae and Pteromalidae, supported by 100% 96% and 98% of bootstrap value respectively. In the Apidae sub-clade, two Dibs from the genus Bombus (B_ter and B_imp) and two from the genus Apis (A_mel and A_flo) clustered together, supported by 100% and 99% of bootstrap values. Among the Dibs of the dipteran clade, two from A. aegypti (A_aeg) and A. gambiae (A_gam) formed a mosquito sub-clade with 100% bootstrap support, and then the mosquito sub-clade and that from D. melanogaster (D_mel) joined together to form the dipteran clade, supported by 100% of bootstrap values. However, four Dib-like proteins from Hemiptera formed two clades. Three putative Dib proteins from S. furcifera (S_fur), L. striatellus (L_str) and A. pisum (A_pis) formed one clade, which was well segregated from a putative Rhodnius prolixus Dib. Moreover, the placements of Dibs from T. castaneum (T_cas) and Pediculus humanus (P_hum) were also not well consistent with the traditional view of insect phylogeny (Fig. 2).

thumbnail
Figure 2. A phylogenetic tree of Dib-like proteins from representative insect species.

An unrooted phylogenetic tree was constructed by the maximum-likelihood method based on the protein sequence alignments. Dib-like proteins from the spider mite Tetranychus urticae (T_urt, JGI_V11_204789), and 19 insect species are shown. Bombus terrestris (B_ter, XP_003394738), B. impatiens (B_imp, XP_003492183), Apis florea (A_flo, XP_003695126), A. mellifera (A_mel, XP_001122832), Camponotus floridanus (C_flo, EFN63712), Acromyrmex echinatior (A_ech, EGI60082), Nasonia vitripennis (N_vit, XP_001601675), Anopheles gambiae (A_gam, ABU42523), Aedes aegypti (A_aeg, AAX85205), Drosophila melanogaster (D_mel, NP_524810), Manduca sexta (M_sex, ABC96069), Bombyx mori (B_mori, NP_001036953), Spodoptera littoralis (S_lit, ACM46003), Tribolium castaneum (T_cas, XP_974252), Pediculus humanus (P_hum, XP_002426284), Acyrthosiphon pisum (A_pis, XP_001948299), Rhodnius prolixus (R_pro, RPRC011595-RA), Sogatella furcifera (S_fur) and Laodelphax striatellus (L_str), respectively. The percentiles of bootstrap values (1,000 replicates) are indicated. The scale bar represents the amino acid divergence.

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

Temporal and Spatial Transcript Profiles

For S. furcifera, second-, third- and fourth-instar nymphs lasted an average of 2.0, 2.0 and 3.0 days. Sfdib showed three expression peaks at day 2 of the second-instar (I2D2), day 2 of the third-instar (I3D2) and day 3 of the fourth-instar nymphs (I4D3). In contrast, the expression levels were lower in the newly-molted second (I2D1)-, third (I3D1)- and fourth (I4D1)-instar nymphs. An ANOVA analysis (P<0.05, df = 6,14) showed that the three peaks were significantly higher than the three troughs (Fig. 3A).

thumbnail
Figure 3. Relative transcript levels of dib genes.

The mRNA levels were measured in S. furcifera second through fourth instars (I2D1, I2D2, I3D1, I3D2, I4D1, I4D2 and I4D3) (A) and in L. striatellus fourth (I4D1 to I4D5) and fifth instars (I5D1) (B) at 24 h intervals, and in the head, thorax and abdomen of S. furcifera I4D3 (C) and L. striatellus I4D5 (D) nymphs. For each sample, 3 independent pools of 5–10 nymphs were measured in technical triplicate using qRT-PCR. The values were calculated using the 2−ΔΔCt method. The relative expression levels were the ratios of relative copy numbers in individuals of specific developmental stage or specific body part to that in S. furcifera I2D1, L. striatellus I4D1 or the head of corresponding species. The bars represent averages with vertical lines indicating SE. Bars with the same lowercase or uppercase letters are not statistically different at P = 0.05 or 0.01.

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

In our rearing protocol, L. striatellus fourth-instar nymphs lasted an average of 5.0 days. Lsdib showed two expression peaks, with the first at day 2 and the second at day 4–5 in fourth-instar nymphs. In contrast, the expression levels were lower and formed two troughs in the newly-molted fourth- and fifth-instar nymphs. An ANOVA analysis (P<0.05, df = 5,12) revealed that the two peaks were significantly higher than the two troughs (Fig. 3B).

Figs. 3C and 3D showed the spatial distribution of Sfdib and Lsdib on day 3 and day 5 of S. furcifera and L. striatellus fourth-instar nymphs. Both Sfdib and Lsdib are expressed at a higher level in the thoraces where prothoracic glands were located. Moreover, trace amounts of transcripts were found in the heads and abdomens. An ANOVA analysis (P<0.01, df = 2,6) demonstrated that the expression levels in thoraces were significantly higher than the mRNA levels in heads and abdomens.

Effect of Dietary Introduction of dsRNA on the Expression of dib and EcR Genes

Dietary introduction of dsdib reduced dib transcript levels. In S. furcifera, a 6 day period of continuous ingestion of the diet containing dsSfdib at the concentration of 0.01, 0.05 and 0.50 mg/ml diminished Sfdib expression level by 35.9%, 45.1% and 66.2% respectively. The three expression levels were significantly lower than those in the blank control (Fig. 4A). Similarly, after 5-day ingestion of dsLsdib at the concentration of 0.01, 0.05 and 0.50 mg/ml, the mRNA levels of Lsdib in treated L. striatellus nymphs were reduced by 18.8%, 35.8% and 56.7% respectively, with the last two significantly lower than that the blank control (Fig. 4B). In contrast, Sfdib and Lsdib expression levels in dsegfp-exposed planthoppers were not changed (Fig. 4A and 4B).

thumbnail
Figure 4. Relative expression levels of dib (A and B) and EcR (C and D) genes.

The mRNA levels were measured in S. furcifera and L. striatellus after continuously exposed to dietary dsRNA at different concentrations. For each sample, 3 independent pools of 5–10 nymphs were measured in technical triplicate using qRT-PCR. The values were calculated using the 2−ΔΔCt method. The relative expression levels were the ratios of relative copy numbers in individuals ingested dsRNA to that in blank control. The bars represent averages with vertical lines indicating SE. Bars with the same lowercase letters are not statistically different at P = 0.05.

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

Since Dib is expected to affect the expression of genes involved in the ecdysteroid-signaling pathway, the effects of Sfdib and Lsdib knockdown were examined on the transcript levels of ecdysone receptor genes SfEcR and LsEcR. EcR constitutes one of the two subunits of the 20E heterodimeric nuclear receptor. EcR is also known to be regulated by 20E through a positive feedback loop directly [37] or indirectly [38] in D. melanogaster. As expected, SfEcR expression levels in S. furcifera nymphs that had ingested dsSfdib at the concentration of 0.01, 0.05 and 0.50 mg/ml significantly decreased by 34.0%, 36.2% and 58.5% respectively, when compared with that in blank control (Fig. 4C). Similarly, LsEcR expression levels in L. striatellus nymphs that had ingested dsLsdib at the concentration of 0.01, 0.05 and 0.50 mg/ml reduced by 25.2%, 46.8% and 68.8% respectively, with the last two concentrations significantly lower than that in blank control (Fig. 4D). In contrast, Sfdib and Lsdib expression level in dsegfp-fed planthoppers was not different than insects in the blank control (Fig. 4C and 4D).

Effects of Ingestion of dsRNA on Nymph Survival

Continuous ingestion of dsegfp-containing diet had no significant impact on nymphal survival relative to blank controls in S. furcifera and L. striatellus. However, when S. furcifera and L. striatellus nymphs had been fed on diets containing dsdib at concentrations of 0.01, 0.05 and 0.50 mg/ml, significant impacts on survival were observed (Fig. 5A and 5B).

thumbnail
Figure 5. Survival (A, B) and duration (C, D) of S. furcifera and L. striatellus.

The nymphs were continuously exposed to dsRNA at different concentrations through the fourth-instar stage. The values represent averages with vertical lines indicating SE. Values with the same lowercase letters are not statistically different at P = 0.05.

https://doi.org/10.1371/journal.pone.0086675.g005

After the experiment, the surviving nymphs were transferred to rice until adult emergence. The mortalities from the termination of experiment to adults emergence were 8.4% (±1.3%, SE), 7.6% (±1.4%), 15.3% (±1.5%), 17.4% (±1.8%), and 21.5% (±2.2%) respectively for S. furcifera on normal, egfp-dsRNA-, and 0.01, 0.05 and 0.50 mg/ml dsdib-containing diets; and were 9.1% (±1.5%), 11.8% (±1.3%), 14.9% (±1.6%), 20.2% (±1.7%), and 28.9% (±2.3%) respectively for L. striatellus nymphs. For both planthopper species, the mortalities at 0.50 mg/ml were significantly higher than those on normal and egfp-dsRNA-containing diets (one-way ANOVA analysis, P<0.05, df = 4,20).

Ingestion of dsRNA on Nymph Development

The development of S. furcifera and L. striatellus nymphs that were fed dsdib-containing diet was significantly delayed. More than 90% surviving planthoppers on normal and egfp-dsRNA-contained diets developed to the fourth (for S. furcifera) or fifth (for L. striatellus) instar. In contrast, only 83.7%, 64.6% and 41.3% of S. furcifera survivors reared on diets containing Sfdib-dsRNA at concentrations of 0.01, 0.05 and 0.50 mg/ml reached the fourth-instar (Fig. 5C). Likewise, 74.3%, 63.2% and 37.8% of L. striatellus individuals reared on diets containing dsLsdib at concentrations of 0.01, 0.05 and 0.50 mg/ml reached the fifth instars (Fig. 5D). In S. furcifera, the time required to develop from the 2nd to the 4th instar on normal diet, dsegfp-diet and diets containing dsSfdib at the concentration of 0.01, 0.05 and 0.50 mg/ml were 7.0 (95% confidence interval 6.2–7.8), 6.7 (5.7–7.7), 8.0 (7.3–8.7), 8.4 (7.5–9.3) and 9.2 (8.3–10.1) days respectively. Nymphs fed the normal or dsegfp-containing diets had significantly shorter developmental times than nymphs fed with 0.50 mg/ml dsDib (P<0.05, df = 4,20) (Fig. 5E). Similarly, the average duration of L. striatellus fourth-instar stages on normal, dsegfp- and the diets containing dsLsdib at the concentration of 0.01, 0.05 and 0.50 mg/ml were 4.9 (4.6–5.2), 5.0 (4.4–5.6), 5.6 (5.0–6.2), 5.8 (4.8–6.8) and 6.5 (5.8–7.2) days respectively. Nymphs fed with diet containing 0.50 mg/ml of dsDib took significantly longer to develop than nymphs fed either normal or dsegfp-diet (P<0.05, df = 4,20) (Fig. 5F).

Rescue Experiment

Ingestion of 20E did not affect the expression level of dib in either S. furcifera or L. striatellus nymphs. In contrast, 20E ingestion restored the mRNA expression levels of SfEcR and LsEcR similar to those of controls (Fig. 6C and 6D). Moreover, 20E ingestion almost completely overcame the negative effects on the nymphal survival (Fig. 6E and 6F) and development in dsdib-treated S. furcifera and L. striatellus nymphs (Fig. 6G to 6J).

thumbnail
Figure 6. Effects of dsdib, 20E, and 20E+dsdib ingestion on S. furcifera (A) and L. striatellus (B).

The nymphs were continuously exposed to dsRNA, 0.5/mL 20E and 0.5 mg/mL 20E+dsRNA. The relative expression levels of dib and EcR genes, mortality and duration were measured. The values represent averages with vertical lines indicating SE. Lines with the same lowercase letters are not statistically different at P = 0.05.

https://doi.org/10.1371/journal.pone.0086675.g006

Discussion

Oral delivery of dsRNA has a potential practical application as an insecticide in the field. Moreover, it is more convenient and easy to manipulate than microinjection [8]. Oral delivery of in vitro-synthesized dsRNA has previously been used to induce gene silencing in a variety of pest insects in Lepidoptera [39][42], Coleoptera [43][47], Hymenoptera [48], Diptera [44], [49][53], Hemiptera [8][13], [16], [44], [54][56], Orthoptera [57] and Isoptera [58]. In the present paper, we cloned and characterized putative disembodied (dib) genes in two planthoppers S. furcifera and L. striatellus, and tested whether oral delivery of dsdib effectively down-regulated the target genes.

Characterization of Sfdib and Lsdib Proteins

We provided here three lines of experimental evidence to support the hypothesis that the putative Sfdib and Lsdib encoded functional Dib proteins. Firstly, predicted SfDib and LsDib showed typical protein domains, and temporal and spatial expression patterns associated with ecdysteroid biosynthesis. SfDib and LsDib have five insect conserved P450 motifs, i.e., Helix-C, Helix-I, Helix-K, PERF and heme-binding motifs. Similar structural characters have been documented in Dib-like proteins from other insect species of diverse orders [17][27]. Moreover, the N-termini of SfDib and LsDib have a mitochondrial targeting sequence and possess two positively charged residues near the heme-binding domain, two important features of mitochondrial P450 enzymes. Consistent with the structural characters, mitochondrial localization of Dib was observed in dib transfected Drosophila S2 cells [59]. In addition, 22-hydroxylase activities are detected in mitochondria prepared from Locusta [60] and Manduca [61].

Our qPCR results revealed that Sfdib and Lsdib showed expression peaks in late instar stages. In contrast, their expression levels were lower in the newly-molted nymphs. In N. lugens [62] and S. littoralis [26], the level of E and 20E showed a peak in the later instar stage. Thus, it can be speculated that the temporal expression pattern of Sfdib and Lsdib coincides with E and 20E titers in the haemolymph. Similarly, the variations of expression levels of Bmdib and Msdib were correlated with the changes in E and 20E titers in the haemolymph of B. mori [20] and M. sexta [21]. Moreover, we found that Sfdib and Lsdib were expressed at higher levels in the thoraces where prothoracic glands were located, compared to the mRNA levels in corresponding heads or abdomens. The spatial expression patterns of Sfdib and Lsdib are consistent with their function in the paired prothoracic glands. Similarly, dib genes were observed to be expressed in the prothoracic glands of various insect species such as D. melanogaster [18], [19], B. mori [20] and M. sexta [21]. In addition, Sfdib and Lsdib were slightly but distinctly expressed in the head and abdomens. Similarly, even though dib had highest expression level in the prothoracic gland of M. sexta fifth-instar larvae, substantial expression was observed in the fat body, midgut, ganglia, Malpighian tubules, and epidermis, in descending order. Moreover, dib was also expressed in several pupal tissues and in the adult ovaries of M. sexta [21]. In S. littoralis, dib was also found in midgut, brain, fat body, epidermis and Malpighian tubules in the sixth instar larvae [26]. These tissues may operate as secondary sources of ecdysteroids by catalyzing terminal hydroxylation from late intermediates, and may provide extra-prothoracic gland ecdysteroids for local effects such as paracrine and autocrine during larval development and metamorphosis.

The second line of experimental evidence was that RNAi-mediated depletion of Sfdib and Lsdib caused phenotypic defects similar to insects whose ecdysteroid synthesis had been disturbed [14], [15], [24] or whose ecdysteroid-mediated signaling had been inhibited [63], [64]. RNAi-mediated knockdown of dib reduced ecdysone receptor gene EcR expression at mRNA levels in both S. furcifera and L. striatellus nymphs. Since EcR is a nuclear receptor of ecdysteroids [65], our results raise the possibility that dietary introduction of Lsdib-RNA inhibit ecdysteroid signaling pathway. Since mutations in and RNAi against EcR cause phenotypic defects and lethality in T. castaneum [63], and in L. striatellus and N. lugens [64], we observed the influence of Lsdib-dsRNA ingestion on nymphal development to test the possibility. As expected, RNAi-mediated depletion of dib genes caused nymphal lethality and developmental delay in both S. furcifera and L. striatellus species.

The third line of evidence was that ingestion of 20E almost completely rescued EcR expression, and overcame the negative effects on the survival and the development in dsdib-treated S. furcifera and L. striatellus nymphs. This suggested that dsdib ingestion by S. furcifera and L. striatellus negatively affects ecdysteroidogenesis, and subsequently down-regulates EcR expression in the two planthoppers. Similarly, knockout and knockdown of Halloween genes caused a decrease in ecdysteroid titers in other insect species [18], [19], [66][70]. Thus, we provided the third line of evidence to support the hypothesis that putative Sfdib and Lsdib genes encode functional Dib proteins that play a critical role in ecdysteroidogenesis in S. furcifera and L. striatellus.

Ingestion of dsRNA Showed Potential Application of Pest Management

In this study, the effects of orally delivered dsRNA on dib and EcR transcript levels in S. furcifera and L. striatellus were dose-dependent, with higher concentrations of dsdib causing higher negative effects. Similar results have also been documented in Aedes aegypti [51] and N. lugens [71]. In previous experiments in D. melanogaster and three other non-dipteran insect species, it was found that oral delivery of dsRNA at doses higher than 0.5 mg/ml was no more effective at inducing RNAi and killing the target insects [44]. Similarly, several studies have also observed that increasing the concentration beyond an optimal dose does not improve the extent of RNAi, although the optimal concentration may vary for the species, mode of delivery, and the life stage [55], [57], [72]. Thus, in the present paper, we tested the negative effects of dietary ingestion of dsdib at the concentration of 0.01, 0.05 and 0.50 mg/ml. We did not find the optimal concentration of dsdib in the present paper, and RNAi-mediated knockdown of dib genes was thus incomplete in S. furcifera and L. striatellus. Even so, our data were comparable to gene silencing levels observed in many other insects that have been fed dsRNA, where the extent of RNAi-induced silencing typically ranged between 40% and 60% [9], [41], [49], [51], [52], [54], [58]. Thus, our results in this study offer two possible applications. Firstly, an oral dsRNA delivery method in pest will facilitate the development of higher throughput RNAi screens in S. furcifera and L. striatellus. Such high throughput screens may identify new targets for control technologies. Secondly, the RNAi-induced effects in planthoppers could serve as a possible method for pest control.

It is widely recognized that two factors are critical for dsRNA-based pesticides, namely effectiveness and dsRNA delivery method [73]. As for effectiveness of dsRNA, our results showed that 5- or 6-day’s ingestion of dsdib by second-instar S. furcifera and fourth-instar L. striatellus nymphs killed approximately 15%–25% planthoppers during the experiments, and caused a further 15%–20% of individual to die after the feeding experiments. RNAi depletion of dib gene should also affect embryonic and nymphal development, and adult reproduction [18], however, here we only tested the negative effect of dsdib ingestion on nymphal development. We need to perform further experiments in future to evaluate the influence of dsdib on adult reproduction and developing embryos.

As for dsRNA delivery, it has been reported that transgenic rice expressing insect-specific dsRNAs mediate RNAi-based knockdown of the target genes in N. lugens insects, although the exact dsRNA concentration in rice was not measured [10]. S. furcifera and L. striatellus deposit their eggs in the leaf sheaths of rice plants, and both nymphs and adults feed on rice phloem sap. Considering that a transgenic rice may provide the insecticidal dsRNA to pests in a stable and potent form [10], [40], [43], [74], it can be expected that the transgenic rice may act as a potential ovicide through planthopper embryonic development stage, as well as a stomach toxin to whole nymphs and adults. Thus, the transgenic rice expressing dsRNAs of Halloween genes such as dib provides a potential application for planthopper control and deserves further investigations.

Author Contributions

Conceived and designed the experiments: GQL PJW. Performed the experiments: PJW SJ NL JMF. Analyzed the data: PJW. Contributed reagents/materials/analysis tools: PJW. Wrote the paper: GQL.

References

  1. 1. Wang L, Zhang Y, Han Z, Liu Y, Fang J (2010) Cross-resistance and possible mechanisms of chlorpyrifos resistance in Laodelphax striatellus (Fallén). Pest Management Science 66: 1096–1100.
  2. 2. Sanada-Morimura S, Sakumoto S, Ohtsu R, Otuka A, Huang S-H, et al. (2011) Current status of insecticide resistance in the small brown planthopper, Laodelphax striatellus, in Japan, Taiwan, and Vietnam. Applied Entomology and Zoology 46: 65–73.
  3. 3. Zhang Y, Li S, Wang L, Zi J, He P, et al. (2013) Overexpression of carboxylesterase-1 and mutation (F439H) of acetylcholinesterase-1 are associated with chlorpyrifos resistance in Laodelphax striatellus. Pesticide Biochemistry and Physiology 106: 8–13.
  4. 4. Zhang Y, Guo H, Yang Q, Li S, Wang L, et al. (2012) Overexpression of a P450 gene (CYP6CW1) in buprofezin-resistant Laodelphax striatellus (Fallén). Pesticide Biochemistry and Physiology 104: 277–282.
  5. 5. Nakao T, Hama M, Kawahara N, Hirase K (2012) Fipronil resistance in Sogatella furcifera: Molecular cloning and functional expression of wild-type and mutant RDL GABA receptor subunits. Journal of Pesticide Science 37: 37–44.
  6. 6. Fire A, Xu SQ, Montgomery MK, Kostas SA, Driver SE, et al. (1998) Potent and specific genetic interference by double-stranded RNA in Caenorhabtidis elegans. Nature 39: 806–811.
  7. 7. Belles X (2010) Beyond Drosophila:RNAi in vitro and functional genomics in insects. Annual Review of Entomology 55: 111–128.
  8. 8. Chen J, Zhang D, Yao Q, Zhang J, Dong X, et al. (2010) Feeding-based RNA interference of a trehalose phosphate synthase gene in the brown planthopper, Nilaparvata lugens. Insect Molecular Biology 19: 777–786.
  9. 9. Li J, Chen QH, Lin YJ, Jiang TR, Wu G, et al. (2011) RNA interference in Nilaparvata lugens (Homoptera: Delphacidae) based on dsRNA ingestion. Pest Management Science 67: 852–859.
  10. 10. Zha W, Peng X, Chen R, Du B, Zhu L, et al. (2011) Knockdown of midgut genes by dsRNA-transgenic plant-mediated RNA interference in the hemipteran insect Nilaparvata lugens. PloS One 6: e20504.
  11. 11. Zhang XY, Lu K, Zhou JL, Zhou Q (2013) Molecular characterization and gene functional analysis of Dicer-2 gene from Nilaparvata lugens (Hemiptera: Geometroidea). Insect Science 20: 61–68.
  12. 12. Sun X-Q, Zhang M-X, Yu J-Y, Jin Y, Ling B, et al. (2013) Glutathione S-transferase of brown planthoppers (Nilaparvata lugens) is essential for their adaptation to gramine-containing host plants. PloS One 8: e64026.
  13. 13. He P, Zhang J, Liu NY, Zhang YN, Yang K, et al. (2011) Distinct expression profiles and different functions of odorant binding proteins in Nilaparvata lugens Stål. PLoS One 6: e28921.
  14. 14. Jia S, Wan P-J, Zhou L-T, Mu L-L, Li G-Q (2013) Knockdown of a putative Halloween gene Shade reveals its role in ecdysteroidogenesis in the small brown planthopper Laodelphax striatellus. Gene 531: 168–174.
  15. 15. Jia S, Wan P-J, Zhou L-T, Mu L-L, Li G-Q (2013) Molecular cloning and RNA interference-mediated functional characterization of a Halloween gene spook in the white-backed planthopper Sogatella furcifera. BMC Molecular Biology 14: 19.
  16. 16. Lu D-H, Wu M, Pu J, Zhang Q, Han Z-J (2013) A functional study of two dsRNA binding protein genes in Laodelphax striatellus. Pest Management Science 69: 1034–1039.
  17. 17. Jürgens G, Wieschaus E, Nüsslein-Volhard C, Kluding H (1984) Mutations affecting the pattern of the larval cuticle in Drosophila melanogaster. II. Zygotic loci on the third chromosome. Roux’s Archives of Developmental Biology 193: 283–295.
  18. 18. Chávez VM, Marqués G, Delbecque JP, Kobayashi K, Hollingsworth M, et al. (2000) The Drosophila disembodied gene controls late embryonic morphogenesis and codes for a cytochrome P450 enzyme that regulates embryonic ecdysone levels. Development 127: 4115–4126.
  19. 19. Warren JT, Petryk A, Marqués G, Jarcho M, Parvy J-P, et al. (2002) Molecular and biochemical characterization of two P450 enzymes in the ecdysteroidogenic pathway of Drosophila melanogaster. Proceedings of the National Academy of Sciences of the United States of America 99: 11043–11048.
  20. 20. Niwa R, Sakudoh T, Namiki T, Saida K, Fujimoto Y, et al. (2005) The ecdysteroidogenic P450 Cyp302a1/disembodied from the silkworm, Bombyx mori, is transcriptionally regulated by prothoracicotropic hormone. Insect Molecular Biology 14: 563–571.
  21. 21. Rewitz KF, Rybczynski R, Warren JT, Gilbert LI (2006) Identification, characterization and developmental expression of Halloween genes encoding P450 enzymes mediating ecdysone biosynthesis in the tobacco hornworm, Manduca sexta. Insect Biochemistry and Molecular Biology 36: 188–199.
  22. 22. Pondeville E, Maria A, Jacques J-C, Bourgouin C, Dauphin-Villemant C (2008) Anopheles gambiae males produce and transfer the vitellogenic steroid hormone 20-hydroxyecdysone to females during mating. Proceedings of the National Academy of Sciences of the United States of America 105: 19631–19636.
  23. 23. Christiaens O, Iga M, Velarde R, Rougé P, Smagghe G (2010) Halloween genes and nuclear receptors in ecdysteroid biosynthesis and signalling in the pea aphid. Insect Molecular Biology 19: 187–200.
  24. 24. Luan J-B, Ghanim M, Liu S-S, Czosnek H (2013) Silencing the ecdysone synthesis and signaling pathway genes disrupts nymphal development in the whitefly. Insect Biochemistry and Molecular Biology 43: 740–746.
  25. 25. Sieglaff DH, Duncan KA, Brown MR (2005) Expression of genes encoding proteins involved in ecdysteroidogenesis in the female mosquito, Aedes aegypti. Insect Biochemistry and Molecular Biology 35: 471–490.
  26. 26. Iga M, Smagghe G (2010) Identification and expression profile of Halloween genes involved in ecdysteroid biosynthesis in Spodoptera littoralis. Peptides 31: 456–467.
  27. 27. Yamazaki Y, Kiuchi M, Takeuchi H, Kubo T (2011) Ecdysteroid biosynthesis in workers of the European honeybee Apis mellifera L. Insect Biochemistry and Molecular Biology. 41: 283–293.
  28. 28. Rewitz K, Gilbert L (2008) Daphnia Halloween genes that encode cytochrome P450s mediating the synthesis of the arthropod molting hormone: Evolutionary implications. BMC Evolutionary Biology 8: 60.
  29. 29. Nelson DR (2009) The cytochrome P450 homepage. Human Genomics 4: 59–65.
  30. 30. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, et al. (2007) Clustal W and Clustal X version 2.0. Bioinformatics 23: 2947–2948.
  31. 31. Stamatakis A (2006) RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690.
  32. 32. Darriba D, Taboada GL, Doallo R, Posada D (2011) ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics 27: 1164–1165.
  33. 33. Pfaffl MW (2001) A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Research 29: e45.
  34. 34. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology 3: research0034.
  35. 35. Feyereisen R (1999) Insect P450 enzymes. Annual Review of Entomology 44: 507–533.
  36. 36. Werck-Reichhart D, Feyereisen R (2000) Cytochromes P450: a success story. Genome Biol 1: researchs 3003.
  37. 37. Karim FD, Thummel C (1992) Temporal coordination of regulatory gene expression by the steroid hormone ecdysone. The EMBO Journal 11: 4083–4093.
  38. 38. Varghese J, Cohen SM (2007) microRNA miR-14 acts to modulate a positive autoregulatory loop controlling steroid hormone signaling in Drosophila. Science 21: 2277–2282.
  39. 39. 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 R NA feeding. Insect Molecular Biology 15: 383–391.
  40. 40. 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.
  41. 41. Bautista MA, Anita M, 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 Biochemistry and Molecular Biology 39: 38–46.
  42. 42. Ren X-L, Chen R-R, Zhang Y, Ma Y, Cui J-J, et al. (2013) A Spodoptera exigua cadherin serves as a putative receptor for Bacillus thuringiensis Cry1Ca toxin and shows differential enhancement to Cry1Ca and Cry1Ac toxicity. Applied and Environmental Microbiology 79: 5576–5583.
  43. 43. 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.
  44. 44. Whyard S, Singh AD, Wong S (2009) Ingested double-stranded RNAs can act as species-specific insecticides. Insect Biochemistry and Molecular Biology 39: 824–832.
  45. 45. Zhu F, Xu J, Palli R, Ferguson J, Palli SR (2011) Ingested RNA interference for managing the populations of the Colorado potato beetle, Leptinotarsa decemlineata. Pest Management Science 67: 175–182.
  46. 46. Wan PJ, Lü D, Guo WC, Ahmat T, Yang L, et al.. (2013) Molecular cloning and characterization of a putative proline dehydrogenase gene in the Colorado potato beetle, Leptinotarsa decemlineata (Say). Insect Science DOI: 10.1111/1744-7917.12034.
  47. 47. Zhou L-T, Jia S, Wan P-J, Kong Y, Guo W-C, et al. (2013) RNA interference of a putative S-adenosyl-L-homocysteine hydrolase gene affects larval performance in Leptinotarsa decemlineata (Say). Journal of Insect Physiology 59: 1049–1056.
  48. 48. Maori E, Paldi N, Shafir S, Kalev H, Tsur E, et al. (2009) IAPV, a bee-affecting virus associated with Colony Collapse Disorder can be silenced by dsRNA ingestion. Insect Molecular Biology 18: 55–60.
  49. 49. Walshe DP, Lehane SM, Lehane MJ, Haines LR (2009) Prolonged gene knockdown in the tsetse fly Glossina by feeding double stranded RNA. Insect Molecular Biology 18: 11–19.
  50. 50. Li X, Zhang M, Zhang H (2011) RNA interference of four genes in adult Bactrocera dorsalis by feeding their dsRNAs. PloS One 6: e17788.
  51. 51. Singh AD, Wong S, Ryan CP, Whyard S (2013) Oral delivery of double-stranded RNA in larvae of the yellow fever mosquito, Aedes aegypti: Implications for pest mosquito control. Journal of Insect Science 13: 69.
  52. 52. Zhang X, Zhang J, Zhu K (2010) Chitosan/double-stranded RNA nanoparticle-mediated RNA interference to silence chitin synthase genes through larval feeding in the African malaria mosquito (Anopheles gambiae). Insect Molecular Biology 19: 683–693.
  53. 53. Cancino-Rodezno A, Alexander C, Villaseñor R, Pacheco S, Porta H, et al. (2010) The mitogen-activated protein kinase p38 is involved in insect defense against Cry toxins from Bacillus thuringiensis. Insect Biochemistry and Molecular Biology 40: 58–63.
  54. 54. Araujo R, Santos A, Pinto F, Gontijo N, Lehane M, 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 Biochemistry and Molecular Biology 36: 683–693.
  55. 55. Shakesby AJ, Wallace IS, Isaacs HV, Pritchard J, Roberts DM, et al. (2009) A water-specific aquaporin involved in aphid osmoregulation. Insect Biochemistry and Molecular Biology 39: 1–10.
  56. 56. Yao J, Rotenberg D, Afsharifar A, Barandoc-Alviar K, Whitfield AE (2013) Development of RNAi methods for Peregrinus maidis, the corn planthopper. PloS One 8: e70243.
  57. 57. Meyering-Vos M, Müller A (2007) RNA interference suggests sulfakinins as satiety effectors in the cricket Gryllus bimaculatus. Journal of Insect Physiology 53: 840–848.
  58. 58. Zhou X, Wheeler MM, Oi FM, Scharf ME (2008) RNA interference in the termite Reticulitermes flavipes through ingestion of double-stranded RNA. Insect Biochemistry and Molecular Biology 38: 805–815.
  59. 59. Petryk A, Warren JT, Marqués G, Jarcho MP, Gilbert LI, et al. (2003) Shade is the Drosophila P450 enzyme that mediates the hydroxylation of ecdysone to the steroid insect molting hormone 20-hydroxyecdysone. Proceedings of the National Academy of Sciences of the United States of America 100: 13773–13778.
  60. 60. Kappler C, Kabbouh M, Hetru C, Durst F, Hoffmann JA (1988) Characterization of three hydroxylases involved in the final steps of biosynthesis of the steroid hormone ecdysone in Locusta migratoria (Insecta Orthoptera). Journal of Steriod Biochemistry 31: 891–898.
  61. 61. Grieneisen ML, Warren JT, Gilbert LI (1993) Early steps in ecdysteroid biosynthesis: evidence for the involvement of cytochrome P-450 enzymes. Insect Biochemistry and Molecular Biology 23: 13–23.
  62. 62. Kobayashi M, Uchida M, Kuriyama K (1989) Elevation of 20-hydroxyecdysone level by buprofezin in Nilaparvata lugens Stål nymphs. Pesticide Biochemistry and Physiology 34: 9–16.
  63. 63. Tan A, Palli SR (2008) Edysone receptor isoforms play distinct roles in controlling molting and metamorphosis in the red flour beetle, Tribolium castaneum. Molecular and Cellular Endocrinology 291: 42–49.
  64. 64. Wu W-J, Wang Y, Huang H-J, Bao Y-Y, Zhang C-X (2012) Ecdysone receptor controls wing morphogenesis and melanization during rice planthopper metamorphosis. Journal of Insect Physiology 58: 420–426.
  65. 65. Yamanaka N, Rewitz KF, O’Connor MB (2013) Ecdysone control of developmental transitions: Lessons from Drosophila research. Annual Review of Entomology 58: 497–516.
  66. 66. Kappler C, Kabbouh M, Hetru C, Durst F, Hoffmann JA (1988) Characterization of three hydroxylases involved in the final steps of biosynthesis of the steroid hormone ecdysone in Locusta migratoria (insecta, orthoptera). Journal of Steroid Biochemistry 31: 891–898.
  67. 67. Marchal E, Badisco L, Verlinden H, Vandersmissen T, Van Soest S, et al. (2011) Role of the Halloween genes, Spook and Phantom in ecdysteroidogenesis in the desert locust, Schistocerca gregaria. Journal of Insect Physiology 57: 1240–1248.
  68. 68. Niwa R, Sakudoh T, Namiki T, Saida K, Fujimoto Y, et al. (2005) The ecdysteroidogenic P450 Cyp302a1/disembodied from the silkworm, Bombyx mori, is transcriptionally regulated by prothoracicotropic hormone. Insect Mol Biol 14: 563–571.
  69. 69. Niwa R, Matsuda T, Yoshiyama T, Namiki T, Mita K, et al. (2004) CYP306A1, a cytochrome P450 enzyme, is essential for ecdysteroid biosynthesis in the prothoracic glands of Bombyx and Drosophila. J Biol Chem 279: 35942–35949.
  70. 70. Warren JT, Petryk A, Marqués G, Parvy J-P, Shinoda T, et al. (2004) Phantom encodes the 25-hydroxylase of Drosophila melanogaster and Bombyx mori: a P450 enzyme critical in ecdysone biosynthesis. Insect Biochemistry and Molecular Biology 34: 991–1010.
  71. 71. Liu S, Ding Z, Zhang C, Yang B, Liu Z (2010) Gene knockdown by intro-thoracic injection of double-stranded RNA in the brown planthopper, Nilaparvata lugens. Insect Biochemistry and Molecular Biology 40: 666–671.
  72. 72. Huvenne H, Smagghe G (2010) Mechanisms of dsRNA uptake in insects and potential of RNAi for pest control: a review. Journal of Insect Physiology 56: 227–235.
  73. 73. Zhang H, Li HC, Miao XX (2013) Feasibility, limitation and possible solutions of RNAi-based technology for insect pest control. Insect Science 20: 15–30.
  74. 74. Pitino M, Coleman AD, Maffei ME, Ridout CJ, Hogenhout SA (2011) Silencing of aphid genes by dsRNA feeding from plants. PLoS One 6: e25709.