Skip to main content
  • Research article
  • Open access
  • Published:

Genome wide screening of RNAi factors of Sf21 cells reveal several novel pathway associated proteins

Abstract

Background

RNA interference (RNAi) leads to sequence specific knock-down of gene expression and has emerged as an important tool to analyse gene functions, pathway analysis and gene therapy. Although RNAi is a conserved cellular process involving common elements and factors, species-specific differences have been observed among different eukaryotes. Identification of components for RNAi pathway is pursued intensively and successful genome-wide screens have been performed for components of RNAi pathways in various organisms. Functional comparative genomics analysis offers evolutionary insight that forms basis of discoveries of novel RNAi-factors within related organisms. Keeping in view the academic and commercial utility of insect derived cell-line from Spodoptera frugiperda, we pursued the identification and functional analysis of components of RNAi-machinery of Sf21 cell-line using genome-wide application.

Results

The genome and transcriptome of Sf21 was assembled and annotated. In silico application of comparative genome analysis among insects allowed us to identify several RNAi factors in Sf21 line. The candidate RNAi factors from assembled genome were validated by knockdown analysis of candidate factors using the siRNA screens on the Sf21-gfp reporter cell-line. Forty two (42) potential factors were identified using the cell based assay. These include core RNAi elements including Dicer-2, Argonaute-1, Drosha, Aubergine and auxiliary modules like chromatin factors, RNA helicases, RNA processing module, signalling allied proteins and others. Phylogenetic analyses and domain architecture revealed that Spodoptera frugiperda homologs retained identity with Lepidoptera (Bombyx mori) or Coleoptera (Tribolium castaneum) sustaining an evolutionary conserved scaffold in post-transcriptional gene silencing paradigm within insects.

Conclusion

The database of RNAi-factors generated by whole genome association survey offers comprehensive outlook about conservation as well as specific differences of the proteins of RNAi machinery. Understanding the interior involved in different phases of gene silencing also offers impending tool for RNAi-based applications.

Background

RNA interference is an evolutionarily conserved gene-regulatory mechanism for host defence against invading genetic elements. RNAi activity is present in a wide variety of eukaryotic organisms [1–4] and is triggered by the presence of double-stranded RNA (dsRNA) or small interfering RNAs (siRNAs), resulting in the degradation or inactivation of cognate messenger RNA [5, 6]. RNAi thus plays important roles in a broad variety of biological pathways such as antiviral defence, epigenetic regulation, DNA elimination and heterochromatin formation etc [7–9]. Biochemical and structural biology studies have provided insights into the molecular details of RNAi and have shown that RNAi is a two-step process in simplistic terms. The first step, referred to as initiation, involves cleavage of dsRNA into small interfering RNAs of ~21-25nt size with two-nucleotide, 3′ overhanging ends [10, 11]. This step is governed in ATP dependent manner by RNase III class of protein(s) called Dicer(s) and few cognate RNA binding protein(s) such as TRBP or R2D2 [12, 13]. The second step, referred to as ‘Slicing’, is the main pivotal degradation step where siRNAs join a multi-nuclease complex called RISC complex, where the endonuclease degrades the single stranded messenger RNA in a site-specific manner. Slicer activity is guided by Argonaute-2, an RNAase H class of enzyme that has been shown to be associated with FMR, Vig and Tudor-SN proteins [14–16]. Besides these basic factors, other accessory factors also have been shown to participate in the RNAi process [17–19]. However, little is known about RNAi intermediates, RNA-protein complexes and mechanisms of formation of different complexes.

Upon completion of whole genome sequences for a number of organisms, a limited number of reports are available for genome-wide screening of components of RNAi pathways. The first genome-wide screen for the identification of genetic components of RNAi was performed using an engineered RNAi sensor strain of Caenorhabditis elegans[20]. The RNAi screen identified about 90 genes which included 11 well known components of RNAi machinery, RNA binding/processing factors, chromatin associated factors, DNA recombination proteins and nuclear import/export proteins etc. The protein-protein interaction maps of these proteins suggested links between germ line and somatic gene silencing as well as RNA-dependent gene regulatory pathways. Later, a couple of studies carried-out genome-wide RNAi screen in Drosophila cell lines and identified seven known RNAi genes (Ago-2, Tis11, Hsc70-3, Hsc70-4 and hdc) and two annotated genes (CG17265 and CG10883) [21]. RNAi has been successfully used in functional genomic studies in many living insect groups such as Coleoptera, Diptera, Dictyoptera, Hemiptera, Hymenoptera, Isoptera, Lepidoptera, Neuroptera and Orthoptera [22–28]. Silencing efficiency has been shown to vary between various insect species and few of these insects are refractory to robust systemic RNAi effects. To understand the differences in the processes of RNAi among various insect species, it is important to carry-out genome-wide RNAi analysis among different insect species. Up till now such a genome-wide RNAi analysis has been performed in vivo in Tribolium castaneum, Nilaparvata lugens and in silico in Bombyx mori. Many of these RNAi factors have been revealed functionally by administration of respective dsRNAs [27, 29]. In the present study; we carried-out genome-wide analysis of the RNAi components in insect Spodoptera frugiperda (Sf21) cell line using both the bioinformatic and the experimental validation approaches.

Sf21 cell line is originally derived from one of the most agronomically important polyphagous pest Spodoptera and is permissive to multiple virus infection. We have successfully used Sf21cells for the functional validation of number of RNAi suppressor proteins [30–32]. Recently, we have assembled the complete genome and transcriptome sequence of Sf21 cell line. Using information from both the genome and transcriptome data, we comprehensively investigated the repertoire of genes involved in RNAi in Sf21 cell line by comparing the sequence data with orthologues from Bombyx mori, Drosophila melanogaster, Tribolium castaneum and C. elegans. Role of these putative RNAi genes was further confirmed using an insect cell-line expressing gfp reporter, developed in our laboratory [30]. An evolutionary conservation of core RNAi gene set was observed. However, a few new RNAi effector components specific to Spodoptera sp. were also identified. We applied available database of the RNAi-factors from different insects to generate comparative profile of Sf21 based RNAi components with other phylogenetically distinct insects, thus providing insight into diversity of RNAi factors.

Results

In silicoidentification of RNAi factors

To perform a genome-wide analysis of RNAi components in Sf21, a genomic library was generated for Sf21 cells and sequenced on Illumina platform. Additionally, genome-wide transcriptome analysis for Sf21 cell line was also accomplished. The Sf21 genome was assembled using Velvet and gene prediction analysis was performed using Augustus. The transcriptome assembly was done using both Velvet and Oases and corresponding ORF sequences were identified using EMBOSS with default parameters [33] (SUB620801). The ORF sequences corresponding to RNAi factors were predicted with the help of UniProt data set for invertebrates and also using sequences of RNAi factors identified in Caenorhabditis elegans, Tribolium castaneum and Drosophila melanogaster genomes as the guide orthologs. By Blast search, we identified a total of 80 potential RNAi factors from the mined genome and transcriptome dataset of Sf21 cells and these factors were evaluated for validation using three factor-specific siRNAs that were transfected in the Sf21-gfp reporter cell line for gfp reversion [30].

In vivo validation of putative Sf21RNAi candidates by reporter based siRNA screen

We have been using gfp expressing Sf21 cell line for the functional genomic studies as well as to understand host-parasite interactions [30, 34]. The RNAi screen for the putative eighty Sf21 RNAi factors was carried out using gfp fluorescent Sf21 cell line. At least three siRNAs were designed and tested for each of the eighty Sf21 RNAi factors (Additional file 1). Each of these siRNAs was co-transfected with gfp siRNA in the stably gfp expressing Sf21 cell line. Gfp fluorescence was monitored by FACS analysis as well as by microscopic examination. The putative siRNAs that were able to restore the gfp fluorescence of the silenced line were analysed and their corresponding genes/proteins were considered as the true RNAi factors (Table 1). The knock down efficiency of each siRNAs specific to putative candidates has been determined a-priori by performing quantitative Real-Time PCR experiment before using these for gfp-reversion experiment. We show the efficacies of a few representative siRNAs in Additional file 2. These siRNAs targeted three genes, namely, Dcr-1, Ago-1 and Drosha for which gfp reversion was scored well and also another three genes, namely, Loquacious, Tudor and Sil-2, which failed to show low or no reversion.

Table 1 Putative RNAi factors of Sf21 quantified by gfp reversion assay

The schematic of gfp reversion assay for identification of putative RNAi candidates in Sf21 has been shown in Figure 1. Figure 2A shows the reversion in gfp expression with siRNA corresponding to putative Dcr-2 gene by microscopic examination. Figure 2B shows quantitative measurement of gfp fluorescence by FACS analysis in lines transfected with siRNAs corresponding to putative Dcr-2 as well as Ago-1 genes. Each of the siRNA transfection experiments were carried out in triplicate and the number of fluorescent cells was recorded from the FACS data. The average number of gfp expressing cells measured in this way has been displayed in Figure 2C. Figure 2C shows the bar graph with ± SD values showing the reversion in gfp expression for few core and accessory RNAi factors. Following identical regimen and protocol, in total forty two candidate RNAi factors were validated from a pool of 80 potential candidates.

Figure 1
figure 1

A genome-wide screen for identification of Sf21 RNAi factors. Schematic representation of experimental design of the siRNA mediated reversion assay in Sf21cell line for RNAi screening.

Figure 2
figure 2

Reversion Assay for putative RNAi factor. (A) (i) Fluorescence imaging of Sf21-gfp reporter cell line (ii) Sf21-gfp reporter cell line transfected with gfp siRNA (iii) Sf21-gfp reporter cell line co-transfected with putative RNAi factor Sf-Dcr-2 specific siRNA with gfp siRNA. (B) Histogram overlay plot shows different gfp expression in transfected Sf21 cell line. Plots depict number of cells (counts) on y axis vs. expression of the gfp reporter (FLH-1) on x axis. Black trace, normal Sf21 cells; green trace, Sf21-gfp reporter cell line; blue trace, siRNA transfected gfp silenced Sf21-gfp cells; red trace, gfp reverted Sf21-gfp cells with putative RNAi factor transfection (a) Sf-Ago-1 (b) Sf-Dcr-2. (C) The bar graph representation of the FACS results for (i) putative core RNAi factors Dicer-2, Argonaute-1, Aubergine, Drosha and (ii) auxiliary RNAi candidates RNA helicase DDX18/HAS1 subfamily, Integrator complex subunit (Int11), Multi Drug Resistance MDR1A transporter and Histone-3 Lysine-4 N-Methyltransferase in Sf21-gfp reporter cell line represented on the x axis, with the percentage of cells expressing gfp on the y axis. Data shown are mean ± SD of three independent experiments.

The experiments were carried out in several replicates so that the data could be statistically valid. However, the variations amongst the replicates were statistically insignificant. For calculating the gfp-reversion values, we have used the value for the particular siRNA that showed maximum reversion within the set of three siRNAs. The particular siRNA was then transfected three times independently for the reversion experiments and the average value of these replicates was reported accordingly. Additional file 3 shows % of gfp quantification from post transfection FACS result of the functional assay for all three sets of siRNAs from each of a few selected representative candidate genes. These genes include core RNAi factors like Dcr-2, Ago-1, Drosha, Pasha, Aubergine, Loquacious which have shown reversion of gfp and others including Auxilliary RNAi factors, like DDX18/HAS1 subfamily RNA helicase, Multi Drug Resistance (MDR1A), Isocitrate Dehydrogenase, Tudor, Sil-2. The table also contains some genes from 80 putative candidates, namely, Serine/threonine protein phosphatase 2A, Myosin-XV -like and Splicing factor 3 subunit 1. Negligible or mild range of gfp reversion was scored with the latter genes.

These genes were further classified based on their perspective role as Core and Auxiliary RNAi components. Dicer, Argonaute, microprocessor complex (Drosha/Pasha), Aubergine, R2D2 and Loquacious were regarded as integral units of RNAi machinery. The auxiliary components comprised of various subcellular regulatory moieties such as DEAD RNA helicase, RNA processing module, chromatin factors, translational machinery, cell division/cell cycle associated elements, protein kinases, signalling allied components, Sid-like proteins and others as part of RNAi-effector complex [20, 29]. These sequences were further analysed to study the domain organization and ancestral relation with their counterparts from other insects, like B. mori, T. castenum and D. melanogaster.

Spodoptera frugiperdaCore RNAi factors

Ribonuclease III family and dsRBM-containing proteins

Dicer proteins belong to RNase III family of bidentate ribonuclease that selectively cleaves double-stranded RNA precursors with 5′-terminal phosphate and a 2-nucleotide 3′ overhang to produce ~21 nucleotides small interfering RNA [10, 35]. These proteins have several conserved domains: two amino terminal DexH-Box helicase domains, a PAZ (Piwi/Zwille) domain, two tandem RNase III domains and carboxy-terminal dsRNA binding domain(s) [36, 37]. Number of Dicer family of proteins varies in different organisms. A single Dicer protein is involved in both the siRNA and miRNA pathways in C. elegans and mammals [20, 38], while Drosophila, the prime insect model organism and Coleoptera Tribolium; have two Dicer enzymes each; Dcr-1 and Dcr-2. Dm-Dcr-1, which lacks an amino-terminal helicase domain but has a PAZ domain, is involved in the miRNA pathway, while Dm-Dcr-2 lacks PAZ domain but has amino-terminal helicase domain and is involved in siRNA pathway [39]. However, only Bm-Dcr-2 has been annotated yet from Silk worm Bombyx mori genome [40].

We identified two Dicer genes Sf-Dcr-1 and Sf-Dcr-2 from Sf21 genome data. Domain architecture of Sf-Dcr-2 exhibits the unique pattern with ATP-dependent DEAD RNA Helicase/ PAZ/ Ribonuclease III/ dsRBM, resembling other insect models. The search for conserved domains revealed that, in contrast to Dm-Dcr-2, Sf-Dcr-2 retains the PAZ domain and unlike Tc-Dcr-2, it holds a carboxy terminal dsRNA binding motif (Figure 3(i)). The clustering array in phylogenetic study reveals that it shares more similarities with Bm-Dcr-2 in terms of domain conservation (Figure 3B).

Figure 3
figure 3

Domain analysis and phylogenetic tree of Dicer and Ribonuclease III protein family. (A) (i) Domain organisation clearly indicates the presence of an amino terminal DEAD domain with significant similarity with other Dicer homologs of Bombyx, Tribolium and Drosophila that helps to classify conservation of Sf21 RNAi factor Dicer-2 gene. (ii) The other Ribonuclease III component Drosha of Tribolium is devoid of N-terminal dsRBM. While Sf21-Drosha shows similarity with Drosophila and Bombyx. (B) Neighbor-joining tree is based on CLUSTALW alignment of full length proteins Dicer-1, Dicer-2 and Drosha as an out group. The phylogenetic tree distributes three different clusters where Sf-Dicer-2 is in the same clade with Bm-Dicer-2 and Sf-Drosha with a Bm-Ribonuclease III. (*) indicates a Bm-Ribonuclease III as a putative Drosha.

Another variant of RNaseIII endonuclease known as Drosha is essential for miRNA biogenesis. Drosha along with its cofactor crops out ~70 nucleotide pri-miRNA product from the precursor pri-miRNA [41–43]. Two tandem RNaseIII domains of Drosha consecutively cleaves the 3′ strand and 5′ strand of double stranded RNA segment of pri-mRNA, generating a 2- nucleotide (nt) 3′ overhang. Sf-Drosha was annotated from S. frugiperda assembled genome based on the presence of a pair of RNA-specific endonuclease III domain and a dsRBM domain similar to Bombyx mori like protein (Figure 3A(ii)). To assess the involvement of Drosha in Sf21 RNAi pathway, we performed the siRNA screening against Drosha in gfp based reversion assay.

Every Dicer or Dicer-like protein is assisted in its function by its cognate double-stranded RNA-binding domain proteins [12, 44, 45]. For example, Drosha functions in association with Pasha, human Drosha works along with DGCR8; Arabidopsis DCL1 is assisted by HYL1 protein etc [46]. Accordingly, we expected that, along with Sf21-Dcr-2, we should be able to locate Dm-R2D2 type of proteins. We found orthologs of Bombyx R2D2 in the Sf21 genome (Figure 4B). Interestingly, a Loquacious homolog in S. frugiperda, which contains three domains of double-stranded RNA binding motif superfamily, was identified too (Figure 4A). The Sf-Loq dsRBM is highly analogous with Bm-loq in comparison to Tribolium or Drosophila homologs (Figure 4C). Dcr-1 mediates miRNA biogenesis involving functional pre-miRNA processing activity in the presence of its cognate partner Loquacious. Three RNA binding motifs of Sf-Loquacious may provide a strong interaction between Dcr-1 and dsRNA. Taken together, these findings suggested that Sf21, which have the same number of Dicer as well as Dicer associated proteins with Drosophila and Tribolium, with similar domain conformation and phylogenetic proximity with Bombyx mori homologs, might be involved in both RNAi and miRNA pathways (Figures 3 and 4), persisting with their conservation of dsRNA processing components in the passage of evolution.

Figure 4
figure 4

Domain analysis and phylogenetic tree of dsRBM effector proteins. (A) Tc, Bm and Sf21 has three tandem dsRBM in Loquacious. (B) The domain architecture of R2D2 is much conserved while all referring insects have two RNA binding domains. (C) Phylogenecic analysis revealed the dsRBM effector components of Sf21 Pasha, R2D2 and Loquacious are closest with respective Bombyx homologs.

Analysis of available transcriptome data set revealed an additional dsRBM which aligns with Pasha/DGCR8 homolog. Correlative phylogenetic study with domain frame characterizes high resemblance of Bm-DGCR8-like protein with Sf-dsRBM, implying it might be Pasha-like protein for Spodoptera. Hence it is tempting to speculate that the pair of Sf-Drosha and Sf-Pasha is present to process the miRNA biogenesis in the standard manner.

Argonaute proteins

The Argonaute proteins are RNaseH family of proteins and are key components of RISC or miRNP complexes that are involved in post-transcriptional silencing, causing either mRNA-cleavage or transcriptional silencing [47]. Argonaute proteins are characterized by PAZ and Piwi domains; PAZ domain is involved in siRNA binding, while the Piwi domain possesses the RNase activity [15, 48] (Figure 5B). Number of Argonaute proteins varies in different organisms: Eight Argonaute genes exist in humans (4 Piwi subfamily and 4 eIF2C/AGO subfamily) [49, 50]. Tribolium and Drosophila have five members, while brown plant hopper Nilaparvata lugens has two orthologs of Drosophila Argonaute genes [27, 29, 51–53].

Figure 5
figure 5

Alignment Domain analysis and phylogenetic tree of Argonaute family of proteins. (A) Amino acid alignment of Argonaute-1 protein of Bombyx mori, Tribolium castaneum, Drosophila melanogaster and Sf21. Indicated sequences were aligned using CLUSTALX. The conserved amino acid residues Asp, Asp and His (DDH) triad residues are marked with red. (B) Domain architecture of Argonaute proteins. The similarity in the domain architectures of Bm-Ago-1 and Sf21-Ago-1 suggests that they are orthologous and the organisation of PAZ/Piwi domain is highly conserved within the insects. (C) The phylogenetic tree was based on consideration of two different class of Argonaute proteins, miRNA class (Argonaute-1) and siRNA class (Argonaute-2). Sf21 genome contains Ago-1 while Tribolium, Bombyx and Drosophila have both Ago-1 and Ago-2.

We identified Ago-1 and Ago-3 from the annotated genome and assembled transcriptome of S. frugiperda (Figures 5B and 6A). The amino acid sequence alignment of Piwi domain in Sf21 revealed that the catalytic motif (aspartate, aspartate and histidine; DDH) is well conserved with other insect species like silkworm, fly and beetle Argonaute members (Figure 5A). However, transient siRNA knockdown assay monitoring gfp expression validated the role of Sf-Ago-1 as potent miRNA class Argonaute. In view of evolutionary lineage and considering different orthologs of Argonaute proteins, S. frugiperda Argonaute was classified as Ago-1 subclass highly analogous to Bm-Ago-1 (Figure 5C). Subsequent comparison of consensus sequences with homologous proteins shows that Sf-Ago-1 appears in the Ago subfamily, whereas Sf-Ago-3 appears to be clustered in the Piwi subfamily (Figure 6C). In view of structural similarity with Bombyx mori Ago-3 component, we assume functionally, Sf-Ago-3 with its sequence of signature PAZ and Piwi domains might be a possible partner for its piRNA-like small RNAs [54].

Figure 6
figure 6

Domain structure and phylogeny of piRNA class components. (A) and (B) domain search finds both Sf21 Argonaute-3 and Aubergine like protein have conserved lineage of PAZ and Piwi domains. (C) A phylogenetic tree was created using the full-length sequences of piRNA class proteins. Siwi represents silkmoth Piwi.

Aubergine

Piwi (P element-induced wimpy testes) are germ line specific Argonaute family proteins that associate with piRNAs and promote cleavage of expressed transposon RNA targets leading to silencing of transposition [55–57]. Two Piwi subclass of RNA interference (RNAi) accessory proteins; Aubergine and Piwi were identified from Drosophila ovaries [52, 58] and genome wide survey of Tribolium also identified a Piwi/Aub class Argonaute involved in transcriptional silencing along with Tc-Ago-3 [29]. We could identify a single Aubergine protein homolog related to the Piwi/Aubergine subfamily proteins in Sf21 cells (Table 1) (Figure 6B, 6C). As Sf21 cells are developed from ovaries of pupal stage Spodoptera frugiperda, it might possess an active piRNA biogenesis pathway and Sf-Aubergine could induce Piwi-associated RNA silencing that is reminiscent of Drosophila and Bombyx mori.

Spodoptera frugiperdaauxilliary RNAi factors

‘DEAD’ box containing RNA helicase domain proteins

DExD/H-box RNA helicase family of enzymes functions as RNA chaperons because of the presence of a unique RNA recognition motifs over an ATP dependent component [59, 60]. Conserved DExH-box helicase, predominantly Dicer and Ago-related helicase genes have been shown to be required for RNAi in C. elegans and human [61, 62]. ATP-dependent subunits p68 and p72 of DEAD box helicase distinguishes pri-miRNAs and induces their processing through Drosha complex in microRNA biogenesis pathway [63].

We identified three distinct classes of DEAD box RNA helicases related to siRNA pathway in Sf21 genome. One of the probable ATP-dependent RNA helicase resembles Dpb45A (DDX49/DBP8 subfamily) of Drosophila melanogaster and DDX49 of Homo sapiens that also has a close homology with another DEAD box helicase family protein of Bombyx mori in turn (Figure 7A, 7B). Another member from the DEAD helicase that belongs to Vasa/DDX4 RNA binding subfamily was identified. Vasa RNA helicases are known to recruit Aubergine which implicates transposable element silencing in Drosophila[64]. A third Helicase from the set of putative RNAi factors was found to be analogous with RNA helicase Pitchoune of Drosophila and Tribolium and DDX18 of Bombyx mori and Homo sapiens, which belongs to DDX18/HAS1 DEAD box helicase subfamily. Presence of a C-terminal DUF4217 (domain of unknown function) besides the ATP-binding and Helicase domains confers a typical feature, which characterizes Sf-DDX18-like protein under DDX18/HAS1 subfamily (Figure 7). Of the three RNAi related helicases identified in Sf21 genome, knocking down the DDX18/HAS1 mRNA provided the highest reversion level of gfp, thereby.

Figure 7
figure 7

Domain analysis and phylogenetic tree of DEAD RNA helicase domain proteins. (A) Conserved Domains Database search shows both Dbp45A and Vasa RNA helicases has DEAD domain and a C-terminal Helicase domain while DDX18/Has1 subfamily of helicase has an unknown domain which is present in all concerned helicase proteins of this family indicating clear orthology. (B) The neighbour-joining tree is based on the alignment of the conserved DEAD domain containing RNA helicases. The DEAD helicases forms three distinct clusters categorized into three subfamilies Dbp45A, Vasa and DDX18/Has1. Sf-DDX18 like proteins is closest with Bombyx mori but similar with Pitchoune of Drosophila and Tribolium too. Vasa subgroup family of proteins being important for a piRNA pathway suggests Sf-Vasa might participate in Piwi associated RNA silencing.

RNA processing module

A number of candidate RNAi genes encode proteins with known or predicted features in RNA binding or processing. Using the reversion of RNAi sensor assay, we identified five factors that participate in RNAi process and possess RNA binding as well as processing modules (Table 1); U1 small nuclear ribonucleoprotein (U1snRNP), SmG, smg-2 like, Integrator complex subunit (Int11) and Zn finger protein. U1 small nuclear ribonucleoprotein (U1snRNP) is known to inhibit the gene expression by preventing poly(A) polymerase (PAP) mediated 3′ polyadenylation. This factor thus couples the splicing and silencing machineries [65, 66]. In the reversion assay, Sf21 U1snRNP showed a moderate level of reversion of RNAi, thus confirming its role in RNAi in Sf21 cell line.

Small nuclear ribonucleoprotein G (SmG) of both Drosophila and Bombyx have been implicated to participate in RNA binding and spliceosomal machinery [67] but not yet analysed as RNAi factor. However, C. elegans SmG protein actively coordinates nematode RNAi [20]. Gene predictions from Spodoptera frugiperda genome assemblies found SmG with a distinct LSM domain. SiRNA guided knockdown assay of Sf-SmG suggested a conserved function of this protein in Sf21 cells gene regulation.

Genome wide RNAi screen in C. elegans classified precursor proteins associated to post-transcriptional pre-mRNA processing like cleavage and polyadenylation associated factors that play roles in RNAi pathway [20]. The integrator complex pairs with carboxyl-terminal domain (CTD) of the largest subunit of RNA polymerase II for co-transcriptional 3′-end processing of U1 and U2 small nuclear RNAs (snRNAs) [68]. We identified a putative factor namely Integrator complex subunit with one domain of Lactamase B superfamily and a C terminal beta-CASP domain in Sf21 genome. Comparative genomic analysis showed that Sf-Int11 is orthologous to the integrator complex subunit 11 of Bombyx mori. Surprisingly, the functional assay for Sf-Int11 showed a very wide range of gfp reversion in Sf21 cell line, therefore requires further investigation for Int11 involvement in siRNA pathway.

An array of Zn finger proteins has been an integral effector of both C. elegans RNA-silencing pathway and miRNA pathway [20, 69]. Lately a chromatin associated RNA processor zinc finger protein has also been identified as a potent component in Piwi-interacting RNA driven transcriptional silencing paradigm in Drosophila[70, 71]. We identified PHD superfamily Sf-Zn finger protein homologous to Dm-INTS12 with defined PHD-finger [72].

Nonsense-mediated decay (NMD) associated proteins like Smg are essential for RNAi mediated mRNA degradation [73]. Studies suggest that series of phosphorylation and de-phosphorylation of NMD factors of mRNA surveillance complex initiates sequestration of nonsense messenger RNA to P bodies [74]. We found one of the core-proteins of the NMD complex Regulator of nonsense transcripts-1 homolog i.e. UPF1 in yeast or smg-2 like protein of C. elegans in Sf21 genome. Sf-Smg-2 has a conserved N-terminal RNA helicase UPF2 interacting domain. However, we observed low levels of gfp reversion with the three siRNAs used in the present study (Table 1).

Protein Kinase family of signalling and allied components

Many components of RNAi machinery have been shown to undergo phosphorylation and are regulated by EGFR/MAPK signalling pathway [19, 75–77]. The genome wide search for the RNAi-genes predicted a number of kinase family of enzymes in Sf21 cell line that have distinct roles in signal transduction pathways. The comparative genome analysis for S. frugiperda RNAi components followed by functional RNAi reversion assay identified eight signalling related components; calcium/calmodulin-dependent protein kinase I (CaM Kinase), two protein kinases (PKC), Serine/threonine p21-activated kinase (PAK) mbt like protein, an inhibitor of nuclear factor kappa-B kinase subunit beta (IKK-beta), MAP Kinase Kinase 4 and STE20/Fray protein kinase (Table 1).

Multi drug resistance cassette transporter

Mechanistic details and functional analysis identified about ten ATP-binding Cassette (ABC) Transporter genes participating in active RNAi function in C. elegans[17]. A half-transporter member HAF-6 of ABCRNAi subclass and its interacting partners MUT-7/RDE-2 have been shown to induce RNAi [78, 79]. Although ATP-binding cassette (ABC) transporter superfamily genes have been widely annotated in Tribolium and Bombyx[80, 81], their potential roles in mounting RNAi response in insects have not been reported. We identified a multi-drug resistance protein 1A-like (Sf-MDR1A) of S. frugiperda with two conserved ABC transporter transmembrane region that produces a robust response in the cell based RNAi-reversion assay (Table 1).

Tudor

Tudor-SN, a staphylococcal/micrococcal nuclease, which harbours a Tudor domain, is a canonical component of active RISC complex that potentiates dsRNA-mediated transcriptional silencing in Caenorhabditis elegans, Drosophila and mammals [16]. A number of proteins with Tudor domain interact with components of gene silencing machinery like Piwi, Armitage, Yb and Zucchini in Drosophila somatic and germ line piRNA pathway [82]. Recent studies in C. elegans ERI endo-RNAi shows that tandem Tudor domain protein ERI-5 stably loads RNA-dependent RNA polymerase (RdRP) complexes to DCR-1-dependent and DCR-1-independent pathways [83]. Surprisingly, RNAi efficiency in Bombyx mori BmN4 cells remains unchanged despite the knockdown of Bm-Tudor-SN [84]. We identified a Sf21-Tudor domain containing protein with two tandem staphylococcal/micrococcal nuclease and one Tudor domain. Importantly, knockdown of Spodoptera Tudor did not affect RNAi-reversion efficiency in our cell based assay (Table 1).

Sid-1 like protein

In C. elegans, systemic RNAi involves Sid-1 mediated intracellular transport of dsRNA, while in human cells, siRNA internalization requires Sid-1 gated transportation [18, 85]. Among the insect species, Drosophila, the premier model system of RNAi research lacks systemic RNAi and Sid-1 as well [86], while Tribolium possess three sets of Sid-1-related genes (Tc-sirA-C) [29] and also, Sid-1 promotes RNAi effect in BmN4 cell line of silkmoth with regards to dsRNA uptake [87]. We identified Sid-1-like channel protein gene homolog Sf-Sil in Sf21 genome and gfp reversion analysis revealed that Sf-Sil might be a putative RNAi effector component. Phylogenetic profile of vertebrate Sid gene family showed separate clades for Ce-Sid, Bm-Sil, Hs-SidT or Tc-Sir. It appears that S. frugiperda Sil belong to insect clades forming a distinct sub-cluster within Bm-Sil branch depicting a significant convergent sequence similarity with the silk-moth as Sf-Sil-2 protein (Figure 8).

Figure 8
figure 8

Sil protein phylogenetic analysis. The neighbour-joining tree is based on the alignment of sil gene family of proteins in different organisms. Sf-Sil-2 sub-clusters with the silkmoth Bm-Sil-2, while Tc-SirA and Tc-SirB compose a distinct cluster. Orthology of these insect is clear from this phylogenetic analysis that Lepidoptera Sid like proteins are evolutionary conserved. Ce-Sid1-3 and Hs-SidT1 and SidT2 diverge out separately and Tc-SirC joins Hs cluster.

Chromatin factors

Recent studies indicate that heterochromatin is not completely inert; it is also transcribed, and rapidly silenced, by a sequence of events directed by RNA interference components [88]. In addition, heterochromatic modifications are also crucial for epigenetic regulation: gene and transposon silencing, X chromosome inactivation (XCI) and imprinting [89–91]. RNAi regulates the assembly and spreading of heterochromatin at the peri-centromeric dg and dh repeats in Schizosaccharomyces pombe[92]. Genome-wide RNAi screens have shown that chromatin associated proteins, histone deacetylase complex (HDAC), and methyltransferases are essential for C. elegans RNAi [20]. We identified Sf-Histone-3 Lysine-4 N-methyltransferase with conserved PHD-like zinc-binding domain and Sf-Histone deacetylase-3-like protein in genome-wide in-silico RNAi screen of Spodoptera genome and these proteins showed moderate levels of reversion in the RNAi reversion assay, suggesting their participatory role in RNAi (Table 1).

Besides methyltransferase, cognate histone acetyltransferase (HAT) component Gas41, which ushers in transcriptional activation [93], have been found to play a secondary role in canonical RNAi pathway [94]. Gas41 thus provides a link between nuclear structure and gene silencing in C. elegans RNAi. Though RNAi sensitivity of transcription factors has not been reported in insect system, we found Sf-Gas41 in our RNAi screens. Sf-Gas41 has YEATS family domain resembling Drosophila counterpart, thus reflecting conservation within insects for the RNAi machinery.

Translational unit

Homology based functional analysis and evolutionary conservation showed that Argonaute family of proteins are probably the primitive extension of components of translation initiation apparatus [95]. Genome scale RNAi screening recognized eIF2B and 60S ribosomal subunit L23a as the RNAi pathway candidates, along with eIF3D as miRNA pathway candidate gene in case of C. elegans[20, 69]. A recent study exhibited that eIF4AII activity is critical for miRNA mediated gene silencing by restraining translation initiation via recruitment into miRISC [96, 97].We identified Sf-eIF2B-gamma, DEAD-box ATPase and ATP-dependent RNA helicase eIF4A paralog proteins Sf-eIF4AII and Sf-eIF4AIII in our screens. In addition, Sf-RPL23P a Ribosomal L23 super family of protein was recognized as siRNA-guided silencing factor. Identification of these proteins raises the intriguing prospect of a link between RNAi pathway and translational components.

Cell division associated elements

In different model organisms, it has been demonstrated that RNA mediated gene regulatory pathway is a modular scaffold, which can coordinate with fundamental units of cell cycle to attune cellular physiology [88]. In Drosophila, functional activity of Ago-1 requires cytokinesis regulator sticky/citron kinase protein for silencing at transcriptional level [98] whereas in S. pombe, Ago-1 interacts with 14-3-3 proteins and this interaction is required for functioning of Argonaute protein in cell cycle and/or gene-silencing pathways [99, 100]. Further, Stoica et al. [100] demonstrated the microtubule walker Kinesin Motor Protein Cut7 as an interacting partner of Ago-1-RISC component [101]. In the present study, we identified two Kinesin like proteins; KIF18A-like and KIF3A-like along with a Cyclin-dependent kinase 5 homolog in Sf RNAi screen. Among these proteins, Kinesin like motor proteins has been first identified in insect to participate in RNAi machinery. However, further experiments will be required to identify their functional role.

Metabolic factor

Little is known about the metabolic processes that drive the RNA silencing process. In C. elegans, investigators identified mitochondrial isocitrate dehydrogenase as a metabolic factor that participates as a nematode RNAi candidate [20]. We also identified Sf-isocitrate dehydrogenase [NAD] subunit alpha as an important factor in RNAi screen using gfp expressing Sf21 cell line. As shown in Table 1, significant reversion was observed when siRNA corresponding to Sf21 isocitrate dehydrogenase gene were used.

Miscellaneous RNAi factors

Based on robust expression in the gfp reporter system, we were able to detect few other RNAi related proteins that did not cluster with previously known functional RNAi factors. These include Myosin VIIa-like protein, Nucleolar complex protein 2 homolog, WD 40-like repeat domain and S-phase kinase-associated protein (SkpA). Rik1 the crucial element in heterochromatic RNA silencing in fission yeast has C-terminal WD-repeat domain [102]. Another screen for C. elegans microRNA candidate identified ubiquitin-dependent catabolic protein Skp-1 homolog responsible for protein turnover during RNA silencing [69]. It will be interesting to functionally characterize the role(s) of these proteins in RNAi pathway in the Sf21 line.

Discussion

With of the advent of NGS technologies, the presence of a host of small RNAs of diverse types along with their relative abundance has been highlighted across various eukaryotic organisms. The biogenesis and functions of these small RNAs are important research aspects as these RNAs control a gamut of biological functions. A vast majority of these pathways are conserved but also diverse not only across various eukaryotic species but also within tissues of organism. In recent years, such diversities have been reported from a variety of organisms including ciliates, fungi, nematodes, plants and animals. These pathways have also highlighted the diversities in the factors that are responsible for RNAi competence and their interacting partners [103–107]. Thus it is not surprising to investigate the varied nature and abundance of the RNAi factors in different organisms. For example, genome-wide approach has identified ~90 proteins involved in RNAi in C. elegans, while in Drosophila melanogaster only ~20 components have been shown to be associated with RNAi process.

There has been a recent splurge in the research pertaining to insect genomics. This is significant not only from crop protection and food security points of view but also from the angle of application of comparative genomics for discovery of novel insecticidal molecule [108]. Among various insect species, genome-wide approaches have been applied on Tribolium castaneum and Nilaparvata lugens and seventeen to fourteen RNAi associated genes have been described respectively. RNAi is being used in our laboratory to carry-out functional genomic studies in Sf21 cell line as well as in intact insects; Spodoptera litura and Helicoverpa armigera[109, 110]. We have recently assembled Sf21 cell transcriptome and genome data. Using these data set and by setting up a functional screen for the putative RNAi factors, in the present study we explored the components that make the Sf21 insect cells competent for the RNAi. To the best of our knowledge, this is the first genome wide approach for a Lepidopteran RNAi factors with in sillico identification and their validation.

To identify components of RNAi machinery in Sf21 cell a homology search was carried out in Sf21 assembled genome based on previously known RNAi factors from C. elegans, B. mori, T. castaneum and D. melanogaster. These searches lead to the identification of about eighty genes and many of these RNAi candidate homologs were close to B. mori gene inventories. The minimal branch divergence of Sf proteins from the cluster of Bombyx mori homologs in the appropriate phylogenetic tree certainly justifies the notion that the Lepidopteron insects evolved the well conserved functions in the evolution. The molecular phylogeny of the insects also highlights the tremendous degree of sequence similarities of the overall RNAi proteins across the insect kingdom.

To confirm the participation of in silico identified eighty S. frugiperda genes in RNAi pathways, a systematic knock down of each of these 80 genes was performed in reversion of silencing assay using Sf21-gfp reporter cell line earlier reported by Singh et al.[30]. Forty one (Table 1) of the eighty putative RNAi factors (Additional file 1) considerably showed block in RNAi in these cell lines, thereby indicating a role(s) of these factors in RNAi pathway. This number of RNAi factors is higher than any other insect species examined so far and it is close to the number (~90) of RNAi factors described for C. elegans[20]. Thus, our results together with available data of C. elegans provide a platform for appreciating diversity in RNAi composition across eukaryotes.

Among the major RNAi factors identified in Sf21 are Dcr-1, Dcr-2, Drosha, Aubergine, Pasha, Ago-1, Ago-3, Loquacious and R2D2. Although a close similarity was observed between these main RNAi factors and the corresponding Bombyx mori factors, however subtle differences were also observed for the factors and such differences have been existent among related insect species. For example, only a single Dicer homolog, Bm-Dicer-2 has been identified in Bombyx mori genome while we identified two Dicers; Dcr-1 and-2 in Sf21 genome.

It is interesting to note that our biological screen also revealed the Dicer gene, albeit in an unexpected manner. Here, we can only speculate about this surprising finding. As reported earlier, Dicer acts as an intergral siRNA-duplex-binding complex by recognising 2-nt 3′ overhangs. Dicer itself binds siRNAs that mimic the yields of Dicer cleavage (19 bp + 2-nt 3′ overhangs) and serve as a prime senor to harbour effective siRNA duplex which further leads to Argonaute loading [111]. Thereby, in our study, we propose the down regulation of dicer by incorporation of Dicer specific siRNA might interfere in the steps of RISC loading. Moreover, down regulation of Dicer might have some indirect effect on the abundance of other RNAi factors for maintenance of homeostatic equilibrium.

Sf21 Drosha protein show similar domain organization as found in B. mori, T. castaneum and D. melanogaster homologs, however Tc-Drosha lacks a c-terminus dsRBM domain. We also identified other pre-Dicer component; Pasha, thereby suggesting a cross-talk between the miRNA and siRNA pathways in these Sf21 cells. It has been proposed that siRNA and miRNA pathways are partially overlapping and few factors participate in both the pathways. In Drosophila, Dcr-2/R2D2 complex is involved in the processing of siRNA pathway from exogenous origin, while endogenous siRNA processing requires Dcr-2/loquacious complex [112]. Since we identified two Dicer genes; Dcr-1 and Dcr-2 in Sf21 genome along with R2D2 and Loquacious, it is possible that Sf-Dcr-2/R2D2 complex is participating in siRNA dependent silencing mechanism and the complex of Sf-Dcr-1/loquacious could be involved in the miRNA pathway of S. frugiperda. It is also possible that there is synergy between these two pathways. Emerging evidence indicates that both Drosophila and Tribolium requires R2D2 for RNAi initiation [13, 29]. Surprisingly, the Silk moth Bm5 cell line is devoid of R2D2 expression [26]. In Sf21, we identified the dsRNA-binding domains of R2D2, which in association with Dcr-2 might be critical for siRNA binding and triggers the assembly of si-RISC complexes.

Among the prime catalytic RNA-endonucleases, Argonaute-1 and Argonaute-3 were identified in Sf21 genome, whereas Drosophila, Tribolium, Bombyx and Nilaparvata have all three different classes of Argonauts; Ago-1, Ago-2 & Ago-3 proteins. It is possible that Sf-Ago-1 is required for siRNA routed RNAi activity in Spodoptera, which is coherent with the report that Dm-Ago-1 functions downstream of siRNA generation [113]. Taking in view the finding of a recent study, we performed amino acid alignment of Sf-Ago-1 that identifies MC motif in the MID domain with two conserved phenylalanine residues at F522 and F557 like Bm-Ago-1. Such correlation establishes the fact that besides translation repression, Sf-Ago-1 being homologous to Bm-Ago-1 might influence P-body localization [114].

Along with the classical RNAi pathway components, we also identified few bona fide piRNA pathway components such as Ago-3, Aubergine, Vasa RNA helicase and Tudor domain protein. Silkworm Piwi (SIWI) and Bm-Ago-3 endogenous expression in ovary-derived BmN4 cell line has been exclusively found during pupation to execute both primary and secondary steps of piRNA biogenesis [54, 115]. In addition to Piwi/Aubergine, Vasa helicase and Tudor proteins have been reported for canonical ping-pong cycle in Drosophila ovary. Since Sf21 cell line have been developed from pupal ovarian origin and from the outcomes of well-known piRNA pathway tissue-specific factors in high throughput assay, we propose an interface between the classical piRNA and siRNA pathways in Spodoptera. Phylogenetic and protein domain analyses of these factors also indicate a convergent similarity with silkmoth piRNA unit. Recently, transgenic RNAi screen in ovarian germ cells led to the discovery of 74 piRNA biogenesis factors in Drosophila[116]. Therefore, a further functional analysis might unravel the involvement of genes like Ago-3, Aubergine, Armitage or Spindle-E for Sf piRNA gene family.

Based on the results of functional RNAi assay, a number of accessory RNAi factors were discovered in the Sf21 genome. These were classified into eleven groups based on their functional domains and reported functions. Similar classification has also been done for C. elegans RNAi factors [20]. Prominent among the factors, which showed maximum RNAi inhibition, are the members of DEAD-box family of RNA helicases, ABC transporter family, factor(s) possessing RNA binding module such as Integrator complex subunit (Int11) and a metabolic factor, isocitrate dehydrogenase. Few of these factors have been reported to interact with the main components of RNAi machinery. For example, conserved dicer-related DExH-box helicase genes drh-1 and drh-2 interact with RDE-4 protein which is required for direct RNAi in C. elegans[62]. In human, DEAH family of RNA helicase named RNA helicase A (RHA) interacts with Ago-2, Dicer and TRBP through its dsRBD. RHA acts as an activation component of siRISC with a functional role in RNAi-mediated silencing [61]. Recent study in HeLa S3 cells indicates DEAD-box helicase 3 (DDX3) as potent factor with sub-cellular localization coupled to Ago-2 in RNAi pathway [117]. Surprisingly, the Sf-ABC transporter family, present in C. elegans, is a lonesome candidate as it is also found only in Sf21 amongst the insects so far [17, 78, 79]. Hence, our data indicate that RNA interference in Sf cells might require ABC transporter dependent trafficking mechanism. It is quite significant that, in spite of multiple duplication in the lineage from nematodes to arthropods, ABC membrane transporters are non-redundantly required for RNAi function.

The ability of dsRNA import and induction of systemic RNAi is been highly conserved in C. elegans and has been correlated with the expression of sid-1 gene. Tribolium and Nilaparvata lugens, both of the insects in which systemic RNA interference is prominent, have Sid-1 like proteins which have higher sequence similarity with insect Sil than with that of C. elegans Tag-130[27, 29]. However, absence of Sid-1 homolog in Drosophila correlates with lack of systemic spreading of RNAi response in fly. The aphids are RNAi-competent but there are controversies about the systemic response in aphids despite the presence of Sid-1 homologs [118]. The Sid-1 gene has been found in highly diverse insect lineage such as Orthoptera, Phthiraptera, Hemiptera, Coleoptera and Hymenoptera but is absent in Diptera [119]. As far as Spodoptera frugiperda is concerned, there are both structural and functional evidence in favour of systemic RNAi-response. We have previously shown the spreading and retention of silencing at the midgut of the closest insect Spodoptera litura causing systemic RNAi [109]. In the present study, we confirm the existence of Sf-sil-2 gene. This coupled with the observation of RNAi transmission in related organism S. litura, suggests the presence of systemic RNAi in Spodoptera frugiperda. As our functional assay was cell-line based, the presence of Sf-Sil-2 gene also leads us to speculate the intra-cellular amplification of the siRNAs.

A distinct number of protein kinase family members are found to be associated with Sf21 RNAi process, raising an interest for a provisional link between RNAi pathway and signalling mechanism. Both of these two processes might have dynamic & overlapping functions within a cell. A number of reports have shown the role of phosphorylation in Ago protein turn over [19] and Ago-2/RISC localization into distinct centres called cytoplasmic bodies or mammalian P-bodies, guiding efficient cleavage and degradation for siRNA-mediated mRNA decay or translational repression [120, 121]. Although more functional studies are required in this regards, we predict similar role(s) for these kinases in Sf21 RNAi processes. Besides these, phosphorylation of small RNAs might also be required for efficient RNAi activities and a subset of these kinases might be involved in such phosphorylation.

The availability of known and predicted protein interactions of closely related insect species like Tribolium helped us to construct a protein-protein interaction map to facilitate a molecular network consisting both core & auxiliary Sf21 RNAi candidates (Figure 9). A search of interolog prediction of Sf21 consensus sequences present in Tribolium included the above mentioned forty two Sf-RNAi factors. The map provides three distinct clusters with ten small networks of RNAi interactors. The cluster comprising core RNAi candidates postulates connections among new auxiliary factors which includes DEAD helicase DDX18/Has1, three of the RNA processing factors SmG, U1AsnRNP, Regulators of nonsense transcript 1 homolog (smg-2 like protein), Translation component eIF4AII, eIF4AIII, cell division elements KIF18A and KIF3A like proteins and Nucleolar complex protein 2 homolog. Another cluster links the other translational unit eIF2B-gamma, Signalling component PAK mbt like protein, PKC, CaM kinase, Chromatin factor Histone deacetylase 3, Gas41 and CDK5 as well. Integrator complex subunit (Int 11), Histone3 Lysin4 N-methyltransferase, WD 40 repeat domain and Zn finger proteins club to a separate group too. However, the interrelatedness between these two separate clusters and the same with core RNAi group are matters of further investigation. The clustering of other small individual networks like Sil-2, isocytrate dehydrogenase, Multidrug resistance MRD1A proteins also challenges the identification of provisional link of the RNAi mediated processes of Sf21 cells.

Figure 9
figure 9

Interactome Map of Sf21 RNAi factors. Protein-protein interaction data for 42 functionally validated RNAi factors identified in the siRNA mediated high through-put screening were constructed by homology based search in Tribolium. The protein interaction data shows 345 nodes with 1257 edges. The nodes in blue highlighted cluster identifies core RNAi components. The red nodes indicate Sf21 RNAi factor homologous to corresponding Tribolium protein. The pink node identifies the interactors of putative candidates postulated from Tribolium network.

Spodoptera frugiperda derived cell line Sf21 is one of the most common insect cell lines used for expression of proteins via baculovirus expression system. In Sf21 cells RNAi technology successfully allows loss of function phenotypes with stable expression. Understanding the participation of RNAi components might be useful information in S. frugiperda cells for more effective and persistent means of transgenic studies. Involvement of such a huge number of accessory elements in small interfering RNA networking gives a global view about organism specific repertoire of RNAi machinery (Table 2). This justifies the rationale where auxiliary RNAi factors are also integral part of gene silencing apparatus. It is possible to target these secondary elements to somehow regulate the strength of the silencing process. Hence this report will help expand the domain RNAi technology as a larvicidal tool and advance the frontiers of RNA interference biology of the major quarantine pests.

Table 2 Comparative analysis of Sf21 RNAi genes

Conclusions

To understand the dogma as well as the thematic variations of gene silencing, functional and comparative genomics serves as a great tool that attempts to describe gene functions and interactions of various components involved in RNA interference mechanism. The study illuminates the integration of small RNA related pathways with other basic cellular activities that govern gene expression for different cellular actions. The conservation of putative core RNAi factors found in insects highlights the biological characteristics of the insect phyla. Despite conservation the molecular evolution of the prime and accessory RNAi candidates of Sf21 cells is quite appreciable. The significance of the unique RNAi factors of Sf cells needs to be assimilated and thus the new insights in the biology of this specific insect will be explored.

Methods

Annotation of RNAi factors from assembled genome and transcriptome data

To identify and characterize putative RNAi candidates in S. frugiperda, we considered predicted ORF’s from Sf21 whole genome sequencing, mRNAs from whole transcriptome sequences and ESTs extracted from SPODOBASE [122] as well. A search for homologs of RNAi factors present in Spodoptera frugiperda genome was carried out using protein sequences of the known RNAi factors of Bombyx mori, Tribolium castaneum, Drosophila melanogaster and C. elegans using Blastx program. Query sequences over 70% coverage were carefully chosen as putative homologs of the corresponding RNAi genes in S. frugiperda genome followed by phylogenetic analysis, domain architecture and CDS generation of consensus of translated proteins [33] (SUB620801).

siRNA designing

A total of eighty putative factors for RNAi were identified in S. frugiperda and for knockout analysis siRNAs for each of the gene were designed by Dharmacon Inc, siRNA designing tool. Multiple siRNAs were predicted for each candidate gene across the gene length with a varying score. In view of the varying in vivo efficacies of different siRNAs of the same gene, 3-5 siRNAs from different regions of the gene with high scores were selected. siRNA list of RNAi factors validated by functional assay with gfp reversion are summarized in Additional file 4.

Cell culture and transfection

To perform siRNA knockdown assay in the Sf21 cell line, a previously developed Sf21 control line (constitutive gfp expressing cells in which gfp was integrated inside the Sf21 genome) was employed [30]. Sf21 cells (Invitrogen) and Sf21-gfp reporter cell line were cultured in TNM-FH Insect Medium (BD BaculoGold) supplemented with Grace’s medium including trace metals, lactalbumin hydrolysate, yeastolate, and 10% heat inactivated fetal bovine serum along with 300 μg/ml Zeocin (100 mg/mL) (Invitrogen). One hour before transfection, cells (80–90% confluent) were sloughed and cell viability was determined by treating the cells with Trypan Blue. Cells with > 95% viability were uniformly transferred to 24-well plate containing BD Baculogold TNM-FH insect medium with cell density 2.5 × 105/well and allow cells to attach. Prior to the addition of transfection mix, cells settled in the plate were washed three times with max-XP serum-free Insect cell medium (BD Baculogold). For reversion assay Sf21-gfp reporter cells were co-transfected with gfp siRNA (GGU UAU GUA CAG GAA CGC AUU) and test siRNA in the presence of Cellfectin II reagent (Invitrogen) incubated in 200 μl BD Baculogold max-XP serum-free medium for 20 minutes. The transfection mixture for Sf21-gfp reporter cells was prepared with gfp siRNA and Cellfectin II reagent which was treated as control gfp silenced cells. Sf21cells and Sf21-gfp reporter cells were also separately incubated only with serum-free medium. Besides, scrambled siRNA (UUG UCU UGC AUU CGA CUA AUdT) with gfp siRNA was used to transfect reporter cell line as negative control. Four hours after transfection, serum medium was added to the culture plate. 48 hours after transfection cells were processed for FACS analysis. All candidate siRNAs were retested in triplicate.

Fluorescence microscopy and flow cytometric analysis

The level of gfp expression in the Sf21-gfp reporter cell lines and siRNA transfected cells were monitored using a Nikon Eclipse TE2000-U Fluorescence microscope (Nikon, Tokyo, Japan) followed by quantification with fluorescence activated cell sorting (FACS). For the FACS analysis, Sf21 cells were washed with FACS-grade phosphate buffered saline (PBS; BD Biosciences) and re-suspended in 400 μl of FACS-grade PBS. To determine the gfp fluorescence of the cells FACSCalibur flow cytometry (Becton-Dickinson) was used. Fluorescence analysis was performed using BD CellQuest Pro Software (Becton-Dickinson). The % of the parent gate detecting gfp fluorescence of cell population for each set of transfection for the selected candidate genes has been listed in Additional file 3. From the FACS data, the following quantitative parameters were determined. These were (A) the number of Sf21-gfp expressing cells co-transfected with test siRNA and gfp siRNA, (B) the same for Sf21-gfp expressing cells transfected only with gfp siRNA, and (C) Sf21-gfp expressing cells. The % gfp reversion was calculated as {(A − B)/(C − B)} × 100.

Validation and efficiency of siRNA transfection by real time PCR

Transfection efficiency of siRNAs was validated by performing a quantitative Real-Time PCR. Two sets of siRNAs of each of the putative candidates were transfected individually to Sf21-gfp cell line using standard protocol. 48 hours after transfection cells were proceed for RNA extraction using Trizol (Invitrogen). Quality and quantity of RNA was checked by using Nanodrop 2000 Spectrophotometer (Thermo Scientific). 100 ng of total RNA was quantified in a single step RT-qPCR assay with two biological replicates each in triplicates using Verso SYBR Green ROX Mix (Thermo Scientific) following manufacturer’s instructions in PIKOREAl 96 detection system (Thermo Scientific). Beta Actin was used as an endogenous control for RNA expression profiling. siRNA induced down regulation of selected genes were then calculated against Sf21-gfp cell line as a calibrator using 2−ΔΔCT method. Relative expression and efficacy of siRNA induced knock down of six selected candidates (Dcr-1, Ago-1, Drosha, Loquacious, Tudor and Sil-2) were represented in Additional file 2. siRNAs used for knock down of selected genes were listed in Additional file 5 and Primers used for their amplification were listed in Additional file 6.

Domain analysis

Domain architecture of proteins was analysed using NCBI conserved domain [123]. The query sequence was searched using NCBI Conserved Domains Database (CDD) v3.11-45746 PSSMs database. Only significant matches of each of the protein were taken where the bits score was greater than or equal to the gathering threshold for the Pfam domain.

Phylogenetic analysis

Sequence manipulation and evolutionary analyses were conducted in MEGA5. Phylogenetic analyses were followed by CLUSTALW and using Neighbour-joining method with bootstrap test (1500 replicates) along with tree visualization was accomplished using MEGA 5 [124]. P-distance method was used to compute the evolutionary distance and was in the units of the number of amino acid differences per site. All positions containing gaps and missing data in amino acid sequences were eliminated.

Search for putative RNAi factors orthologs

ORF of putative RNAi factors were used as the query in a Uniprot Blastx against arthropods, human and nematodes. Homology search was also crosschecked with Blastx against Tribolium castaneum, Drosophila melanogaster, Bombyx mori, Caenorhabditis elegans, Homo sapiens in the NCBI non-redundant protein sequence database to identify their respective orthologs. Most significant alignment with lowest E-value from each of the organisms were considered with the threshold of >50% in terms of either maximum identity or query coverage.

Interactome analysis of Sf21RNAi factors

Known and Predicted Protein Interactions of Tribolium homologs of Sf21 RNAi factors have been accumulated from STRING 9.1 database. Visualization of molecular interaction networks have been constructed using Cytoscape_v2.8.3 [125].

Accession numbers

The Accession numbers used for identity search are summarized in Additional file 7. The nucleotide sequences that have been used to design siRNA reported in this study have been described in the manuscript [33] (SUB620801).

References

  1. Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC: Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans. Nature. 1998, 391: 806-811.

    CAS  PubMed  Google Scholar 

  2. Hammond SM, Bernstein E, Beach D, Hannon GJ: An RNA-directed nuclease mediates post-transcriptional gene silencing in Drosophila cells. Nature. 2000, 404: 293-296.

    CAS  PubMed  Google Scholar 

  3. Chuang CF, Meyerowitz EM: Specific and heritable genetic interference by double-stranded RNA in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2000, 97: 4985-4990.

    CAS  PubMed Central  PubMed  Google Scholar 

  4. Li YX, Farrell MJ, Liu R, Mohanty N, Kirby ML: Double-stranded RNA injection produces null phenotypes in zebrafish. Dev Biol. 2000, 217: 394-405.

    CAS  PubMed  Google Scholar 

  5. Meister G, Tuschl T: Mechanisms of gene silencing by double-stranded RNA. Nature. 2004, 431: 343-349.

    CAS  PubMed  Google Scholar 

  6. Zamore PD, Tuschl T, Sharp PA, Bartel DP: RNAi: double-stranded RNA directs the ATP-dependent cleavage of mRNA at 21 to 23 nucleotide intervals. Cell. 2000, 101: 25-33.

    CAS  PubMed  Google Scholar 

  7. Matzke MA, Birchler JA: RNAi-mediated pathways in the nucleus. Nat Rev Genet. 2005, 6: 24-35.

    CAS  PubMed  Google Scholar 

  8. Verdel A, Jia S, Gerber S, Sugiyama T, Gygi S, Grewal SI, Moazed D: RNAi-mediated targeting of heterochromatin by the RITS complex. Science. 2004, 303: 672-676.

    CAS  PubMed Central  PubMed  Google Scholar 

  9. Wilkins C, Dishongh R, Moore SC, Whitt MA, Chow M, Machaca K: RNA interference is an antiviral defence mechanism in Caenorhabditis elegans. Nature. 2005, 436: 1044-1047.

    CAS  PubMed  Google Scholar 

  10. Bernstein E, Caudy AA, Hammond SM, Hannon GJ: Role for a bidentate ribonuclease in the initiation step of RNA interference. Nature. 2001, 409: 363-366.

    CAS  PubMed  Google Scholar 

  11. Elbashir SM, Lendeckel W, Tuschl T: RNA interference is mediated by 21- and 22-nucleotide RNAs. Genes Dev. 2001, 15: 188-200.

    CAS  PubMed Central  PubMed  Google Scholar 

  12. Chendrimada TP, Gregory RI, Kumaraswamy E, Norman J, Cooch N, Nishikura K, Shiekhattar R: TRBP recruits the Dicer complex to Ago2 for microRNA processing and gene silencing. Nature. 2005, 436: 740-744.

    CAS  PubMed Central  PubMed  Google Scholar 

  13. Liu Q, Rand TA, Kalidas S, Du F, Kim HE, Smith DP, Wang X: R2D2, a bridge between the initiation and effector steps of the Drosophila RNAi pathway. Science. 2003, 301: 1921-1925.

    CAS  PubMed  Google Scholar 

  14. Filipowicz W: RNAi: the nuts and bolts of the RISC machine. Cell. 2005, 122: 17-20.

    CAS  PubMed  Google Scholar 

  15. Song JJ, Smith SK, Hannon GJ, Joshua-Tor L: Crystal structure of Argonaute and its implications for RISC slicer activity. Science. 2004, 305: 1434-1437.

    CAS  PubMed  Google Scholar 

  16. Caudy AA, Ketting RF, Hammond SM, Denli AM, Bathoorn AM, Tops BB, Silva JM, Myers MM, Hannon GJ, Plasterk RH: A micrococcal nuclease homologue in RNAi effector complexes. Nature. 2003, 425: 411-414.

    CAS  PubMed  Google Scholar 

  17. Sundaram P, Echalier B, Han W, Hull D, Timmons L: ATP-binding cassette transporters are required for efficient RNA interference in Caenorhabditis elegans. Mol Biol Cell. 2006, 17: 3678-3688.

    CAS  PubMed Central  PubMed  Google Scholar 

  18. Winston WM, Molodowitch C, Hunter CP: Systemic RNAi in C. elegans requires the putative transmembrane protein SID-1. Science. 2002, 295: 2456-2459.

    CAS  PubMed  Google Scholar 

  19. Zeng Y, Sankala H, Zhang X, Graves PR: Phosphorylation of Argonaute 2 at serine-387 facilitates its localization to processing bodies. Biochem J. 2008, 413: 429-436.

    CAS  PubMed  Google Scholar 

  20. Kim JK, Gabel HW, Kamath RS, Tewari M, Pasquinelli A, Rual JF, Kennedy S, Dybbs M, Bertin N, Kaplan JM, Vidal M, Ruvkun G: Functional genomic analysis of RNA interference in C. elegans. Science. 2005, 308: 1164-1167.

    CAS  PubMed  Google Scholar 

  21. Dorner S, Lum L, Kim M, Paro R, Beachy PA, Green R: A genomewide screen for components of the RNAi pathway in Drosophila cultured cells. Proc Natl Acad Sci U S A. 2006, 103: 11880-11885.

    CAS  PubMed Central  PubMed  Google Scholar 

  22. Bucher G, Scholten J, Klingler M: Parental RNAi in Tribolium (Coleoptera). Curr Biol. 2002, 12: R85-86.

    CAS  PubMed  Google Scholar 

  23. Dong Y, Friedrich M: Nymphal RNAi: systemic RNAi mediated gene knockdown in juvenile grasshopper. BMC Biotechnol. 2005, 5: 25-

    PubMed Central  PubMed  Google Scholar 

  24. Jarosch A, Moritz RF: Systemic RNA-interference in the honeybee Apis mellifera: tissue dependent uptake of fluorescent siRNA after intra-abdominal application observed by laser-scanning microscopy. J Insect Physiol. 2011, 57: 851-857.

    CAS  PubMed  Google Scholar 

  25. Li J, Wang X, Wang M, Ma W, Hua H: Advance of RNA interference technique in Hemipteran insects. Insect Sci. 2012, 20: 31-39.

    PubMed  Google Scholar 

  26. Swevers L, Liu J, Huvenne H, Smagghe G: Search for limiting factors in the RNAi pathway in silkmoth tissues and the Bm5 cell line: the RNA-binding proteins R2D2 and Translin. PLoS ONE. 2011, 6: e20250-

    CAS  PubMed Central  PubMed  Google Scholar 

  27. Xu HJ, Chen T, Ma XF, Xue J, Pan PL, Zhang XC, Cheng JA, Zhang CX: Genome-wide screening for components of small interfering RNA (siRNA) and micro-RNA (miRNA) pathways in the brown planthopper, Nilaparvata lugens (Hemiptera: Delphacidae). Insect Mol Biol. 2013, 22: 635-647.

    CAS  PubMed  Google Scholar 

  28. Zhou X, Wheeler MM, Oi FM, Scharf ME: RNA interference in the termite Reticulitermes flavipes through ingestion of double-stranded RNA. Insect Biochem Mol Biol. 2008, 38: 805-815.

    CAS  PubMed  Google Scholar 

  29. Tomoyasu Y, Miller SC, Tomita S, Schoppmeier M, Grossmann D, Bucher G: Exploring systemic RNA interference in insects: a genome-wide survey for RNAi genes in Tribolium. Genome Biol. 2008, 9: R10-

    PubMed Central  PubMed  Google Scholar 

  30. Singh G, Popli S, Hari Y, Malhotra P, Mukherjee S, Bhatnagar RK: Suppression of RNA silencing by Flock house virus B2 protein is mediated through its interaction with the PAZ domain of Dicer. FASEB J. 2009, 23: 1845-1857.

    CAS  PubMed  Google Scholar 

  31. Chinnappan M, Singh AK, Kakumani PK, Kumar G, Rooge SB, Kumari A, Varshney A, Rastogi A, Sarin SK, Malhotra P, Mukherjee SK, Bhatnagar RK: Key elements of the RNAi pathway are regulated by hepatitis B virus replication and HBx acts as a viral suppressor of RNA silencing. Biochem J. 2014, 462: 347-358.

    CAS  PubMed  Google Scholar 

  32. Kakumani PK, Ponia SS, S RK, Sood V, Chinnappan M, Banerjea AC, Medigeshi GR, Malhotra P, Mukherjee SK, Bhatnagar RK: Role of RNA interference (RNAi) in dengue virus replication and identification of NS4B as an RNAi suppressor. J Virol. 2013, 87: 8870-8883.

    CAS  PubMed Central  PubMed  Google Scholar 

  33. Kakumani PK, Malhotra P, Mukherjee SK, Bhatnagar RK: A draft genome assembly of the army worm, Spodoptera frugiperda. Genomics. 2014, 104: 134-143.

    CAS  PubMed  Google Scholar 

  34. Singh G, Sachdev B, Sharma N, Seth R, Bhatnagar RK: Interaction of Bacillus thuringiensis vegetative insecticidal protein with ribosomal S2 protein triggers larvicidal activity in Spodoptera frugiperda. Appl Environ Microbiol. 2010, 76: 7202-7209.

    CAS  PubMed Central  PubMed  Google Scholar 

  35. Park JE, Heo I, Tian Y, Simanshu DK, Chang H, Jee D, Patel DJ, Kim VN: Dicer recognizes the 5′ end of RNA for efficient and accurate processing. Nature. 2011, 475: 201-205.

    CAS  PubMed  Google Scholar 

  36. Knight SW, Bass BL: A role for the RNase III enzyme DCR-1 in RNA interference and germ line development in Caenorhabditis elegans. Science. 2001, 293: 2269-2271.

    CAS  PubMed Central  PubMed  Google Scholar 

  37. MacRae IJ, Doudna JA: Ribonuclease revisited: structural insights into ribonuclease III family enzymes. Curr Opin Struct Biol. 2007, 17: 138-145.

    CAS  PubMed  Google Scholar 

  38. Zhang H, Kolb FA, Jaskiewicz L, Westhof E, Filipowicz W: Single processing center models for human Dicer and bacterial RNase III. Cell. 2004, 118: 57-68.

    CAS  PubMed  Google Scholar 

  39. Lee YS, Nakahara K, Pham JW, Kim K, He Z, Sontheimer EJ, Carthew RW: Distinct roles for Drosophila Dicer-1 and Dicer-2 in the siRNA/miRNA silencing pathways. Cell. 2004, 117: 69-81.

    CAS  PubMed  Google Scholar 

  40. Consortium ISG: The genome of a lepidopteran model insect, the silkworm Bombyx mori. Insect Biochem Mol Biol. 2008, 38: 1036-1045.

    Google Scholar 

  41. Landthaler M, Yalcin A, Tuschl T: The human DiGeorge syndrome critical region gene 8 and Its D. melanogaster homolog are required for miRNA biogenesis. Curr Biol. 2004, 14: 2162-2167.

    CAS  PubMed  Google Scholar 

  42. Lee Y, Ahn C, Han J, Choi H, Kim J, Yim J, Lee J, Provost P, Radmark O, Kim S, Kim VN: The nuclear RNase III Drosha initiates microRNA processing. Nature. 2003, 425: 415-419.

    CAS  PubMed  Google Scholar 

  43. Yeom KH, Lee Y, Han J, Suh MR, Kim VN: Characterization of DGCR8/Pasha, the essential cofactor for Drosha in primary miRNA processing. Nucleic Acids Res. 2006, 34: 4622-4629.

    CAS  PubMed Central  PubMed  Google Scholar 

  44. Saito K, Ishizuka A, Siomi H, Siomi MC: Processing of pre-microRNAs by the Dicer-1-Loquacious complex in Drosophila cells. PLoS Biol. 2005, 3: e235-

    PubMed Central  PubMed  Google Scholar 

  45. Liu X, Jiang F, Kalidas S, Smith D, Liu Q: Dicer-2 and R2D2 coordinately bind siRNA to promote assembly of the siRISC complexes. RNA. 2006, 12: 1514-1520.

    CAS  PubMed Central  PubMed  Google Scholar 

  46. Kurihara Y, Takashi Y, Watanabe Y: The interaction between DCL1 and HYL1 is important for efficient and precise processing of pri-miRNA in plant microRNA biogenesis. RNA. 2006, 12: 206-212.

    CAS  PubMed Central  PubMed  Google Scholar 

  47. Hammond SM, Boettcher S, Caudy AA, Kobayashi R, Hannon GJ: Argonaute2, a link between genetic and biochemical analyses of RNAi. Science. 2001, 293: 1146-1150.

    CAS  PubMed  Google Scholar 

  48. Song JJ, Liu J, Tolia NH, Schneiderman J, Smith SK, Martienssen RA, Hannon GJ, Joshua-Tor L: The crystal structure of the Argonaute2 PAZ domain reveals an RNA binding motif in RNAi effector complexes. Nat Struct Biol. 2003, 10: 1026-1032.

    CAS  PubMed  Google Scholar 

  49. Meister G, Landthaler M, Patkaniowska A, Dorsett Y, Teng G, Tuschl T: Human Argonaute2 mediates RNA cleavage targeted by miRNAs and siRNAs. Mol Cell. 2004, 15: 185-197.

    CAS  PubMed  Google Scholar 

  50. Sasaki T, Shiohama A, Minoshima S, Shimizu N: Identification of eight members of the Argonaute family in the human genome small star, filled. Genomics. 2003, 82: 323-330.

    CAS  PubMed  Google Scholar 

  51. Okamura K, Ishizuka A, Siomi H, Siomi MC: Distinct roles for Argonaute proteins in small RNA-directed RNA cleavage pathways. Genes Dev. 2004, 18: 1655-1666.

    CAS  PubMed Central  PubMed  Google Scholar 

  52. Vagin VV, Sigova A, Li C, Seitz H, Gvozdev V, Zamore PD: A distinct small RNA pathway silences selfish genetic elements in the germline. Science. 2006, 313: 320-324.

    CAS  PubMed  Google Scholar 

  53. Brennecke J, Aravin AA, Stark A, Dus M, Kellis M, Sachidanandam R, Hannon GJ: Discrete small RNA-generating loci as master regulators of transposon activity in Drosophila. Cell. 2007, 128: 1089-1103.

    CAS  PubMed  Google Scholar 

  54. Kawaoka S, Minami K, Katsuma S, Mita K, Shimada T: Developmentally synchronized expression of two Bombyx mori Piwi subfamily genes, SIWI and BmAGO3 in germ-line cells. Biochem Biophys Res Commun. 2008, 367: 755-760.

    CAS  PubMed  Google Scholar 

  55. Gunawardane LS, Saito K, Nishida KM, Miyoshi K, Kawamura Y, Nagami T, Siomi H, Siomi MC: A slicer-mediated mechanism for repeat-associated siRNA 5′ end formation in Drosophila. Science. 2007, 315: 1587-1590.

    CAS  PubMed  Google Scholar 

  56. O’Donnell KA, Boeke JD: Mighty Piwis defend the germline against genome intruders. Cell. 2007, 129: 37-44.

    PubMed Central  PubMed  Google Scholar 

  57. Zamore PD: RNA silencing: genomic defence with a slice of pi. Nature. 2007, 446: 864-865.

    CAS  PubMed  Google Scholar 

  58. Saito K, Nishida KM, Mori T, Kawamura Y, Miyoshi K, Nagami T, Siomi H, Siomi MC: Specific association of Piwi with rasiRNAs derived from retrotransposon and heterochromatic regions in the Drosophila genome. Genes Dev. 2006, 20: 2214-2222.

    CAS  PubMed Central  PubMed  Google Scholar 

  59. Cordin O, Banroques J, Tanner NK, Linder P: The DEAD-box protein family of RNA helicases. Gene. 2006, 367: 17-37.

    CAS  PubMed  Google Scholar 

  60. Cordin O, Tanner NK, Doere M, Linder P, Banroques J: The newly discovered Q motif of DEAD-box RNA helicases regulates RNA-binding and helicase activity. EMBO J. 2004, 23: 2478-2487.

    CAS  PubMed Central  PubMed  Google Scholar 

  61. Robb GB, Rana TM: RNA helicase A interacts with RISC in human cells and functions in RISC loading. Mol Cell. 2007, 26: 523-537.

    CAS  PubMed  Google Scholar 

  62. Tabara H, Yigit E, Siomi H, Mello CC: The dsRNA binding protein RDE-4 interacts with RDE-1, DCR-1, and a DExH-box helicase to direct RNAi in C. elegans. Cell. 2002, 109: 861-871.

    CAS  PubMed  Google Scholar 

  63. Fukuda T, Yamagata K, Fujiyama S, Matsumoto T, Koshida I, Yoshimura K, Mihara M, Naitou M, Endoh H, Nakamura T, Akimoto C, Yamamoto Y, Katagiri T, Foulds C, Takezawa S, Kitagawa H, Takeyama K, O'Malley BW, Kato S: DEAD-box RNA helicase subunits of the Drosha complex are required for processing of rRNA and a subset of microRNAs. Nat Cell Biol. 2007, 9: 604-611.

    CAS  PubMed  Google Scholar 

  64. Vagin VV, Klenov MS, Kalmykova AI, Stolyarenko AD, Kotelnikov RN, Gvozdev VA: The RNA interference proteins and vasa locus are involved in the silencing of retrotransposons in the female germline of Drosophila melanogaster. RNA Biol. 2004, 1: 54-58.

    CAS  PubMed  Google Scholar 

  65. Furth PA, Choe WT, Rex JH, Byrne JC, Baker CC: Sequences homologous to 5′ splice sites are required for the inhibitory activity of papillomavirus late 3′ untranslated regions. Mol Cell Biol. 1994, 14: 5278-5289.

    CAS  PubMed Central  PubMed  Google Scholar 

  66. Gunderson SI, Polycarpou-Schwarz M, Mattaj IW: U1 snRNP inhibits pre-mRNA polyadenylation through a direct interaction between U1 70K and poly(A) polymerase. Mol Cell. 1998, 1: 255-264.

    CAS  PubMed  Google Scholar 

  67. Mount SM, Salz HK: Pre-messenger RNA processing factors in the Drosophila genome. J Cell Biol. 2000, 150: F37-44.

    CAS  PubMed  Google Scholar 

  68. Baillat D, Hakimi MA, Naar AM, Shilatifard A, Cooch N, Shiekhattar R: Integrator, a multiprotein mediator of small nuclear RNA processing, associates with the C-terminal repeat of RNA polymerase II. Cell. 2005, 123: 265-276.

    CAS  PubMed  Google Scholar 

  69. Parry DH, Xu J, Ruvkun G: A whole-genome RNAi Screen for C. elegans miRNA pathway genes. Curr Biol. 2007, 17: 2013-2022.

    CAS  PubMed Central  PubMed  Google Scholar 

  70. Donertas D, Sienski G, Brennecke J: Drosophila Gtsf1 is an essential component of the Piwi-mediated transcriptional silencing complex. Genes Dev. 2013, 27: 1693-1705.

    PubMed Central  PubMed  Google Scholar 

  71. Ohtani H, Iwasaki YW, Shibuya A, Siomi H, Siomi MC, Saito K: DmGTSF1 is necessary for Piwi-piRISC-mediated transcriptional transposon silencing in the Drosophila ovary. Genes Dev. 2013, 27: 1656-1661.

    CAS  PubMed Central  PubMed  Google Scholar 

  72. Chen J, Waltenspiel B, Warren WD, Wagner EJ: Functional analysis of the integrator subunit 12 identifies a microdomain that mediates activation of the Drosophila integrator complex. J Biol Chem. 2013, 288: 4867-4877.

    CAS  PubMed Central  PubMed  Google Scholar 

  73. Domeier ME, Morse DP, Knight SW, Portereiko M, Bass BL, Mango SE: A link between RNA interference and nonsense-mediated decay in Caenorhabditis elegans. Science. 2000, 289: 1928-1931.

    CAS  PubMed  Google Scholar 

  74. Eulalio A, Behm-Ansmant I, Izaurralde E: P bodies: at the crossroads of post-transcriptional pathways. Nat Rev Mol Cell Biol. 2007, 8: 9-22.

    CAS  PubMed  Google Scholar 

  75. Adams BD, Claffey KP, White BA: Argonaute-2 expression is regulated by epidermal growth factor receptor and mitogen-activated protein kinase signaling and correlates with a transformed phenotype in breast cancer cells. Endocrinology. 2009, 150: 14-23.

    CAS  PubMed Central  PubMed  Google Scholar 

  76. Paroo Z, Ye X, Chen S, Liu Q: Phosphorylation of the human microRNA-generating complex mediates MAPK/Erk signaling. Cell. 2009, 139: 112-122.

    CAS  PubMed Central  PubMed  Google Scholar 

  77. Rudel S, Wang Y, Lenobel R, Korner R, Hsiao HH, Urlaub H, Patel D, Meister G: Phosphorylation of human Argonaute proteins affects small RNA binding. Nucleic Acids Res. 2011, 39: 2330-2343.

    PubMed Central  PubMed  Google Scholar 

  78. Sundaram P, Han W, Cohen N, Echalier B, Albin J, Timmons L: Caenorhabditis elegans ABCRNAi transporters interact genetically with rde-2 and mut-7. Genetics. 2008, 178: 801-814.

    CAS  PubMed Central  PubMed  Google Scholar 

  79. Timmons LD: ABC transporters and RNAi in Caenorhabditis elegans. J Bioenerg Biomembr. 2007, 39: 459-463.

    CAS  PubMed  Google Scholar 

  80. Broehan G, Kroeger T, Lorenzen M, Merzendorfer H: Functional analysis of the ATP-binding cassette (ABC) transporter gene family of Tribolium castaneum. BMC Genomics. 2013, 14: 6-

    CAS  PubMed Central  PubMed  Google Scholar 

  81. Liu S, Zhou S, Tian L, Guo E, Luan Y, Zhang J, Li S: Genome-wide identification and characterization of ATP-binding cassette transporters in the silkworm, Bombyx mori. BMC Genomics. 2011, 12: 491-

    CAS  PubMed Central  PubMed  Google Scholar 

  82. Handler D, Olivieri D, Novatchkova M, Gruber FS, Meixner K, Mechtler K, Stark A, Sachidanandam R, Brennecke J: A systematic analysis of Drosophila TUDOR domain-containing proteins identifies Vreteno and the Tdrd12 family as essential primary piRNA pathway factors. EMBO J. 2011, 30: 3977-3993.

    CAS  PubMed Central  PubMed  Google Scholar 

  83. Thivierge C, Makil N, Flamand M, Vasale JJ, Mello CC, Wohlschlegel J, Conte D, Duchaine TF: Tudor domain ERI-5 tethers an RNA-dependent RNA polymerase to DCR-1 to potentiate endo-RNAi. Nat Struct Mol Biol. 2011, 19: 90-97.

    PubMed Central  PubMed  Google Scholar 

  84. Zhu L, Tatsuke T, Li Z, Mon H, Xu J, Lee JM, Kusakabe T: Molecular cloning of BmTUDOR-SN and analysis of its role in the RNAi pathway in the silkworm, Bombyx mori (Lepidoptera: Bombycidae). Appl Entomol Zool. 2012, 47: 207-215.

    CAS  Google Scholar 

  85. Duxbury MS, Ashley SW, Whang EE: RNA interference: a mammalian SID-1 homologue enhances siRNA uptake and gene silencing efficacy in human cells. Biochem Biophys Res Commun. 2005, 331: 459-463.

    CAS  PubMed  Google Scholar 

  86. Roignant JY, Carre C, Mugat B, Szymczak D, Lepesant JA, Antoniewski C: Absence of transitive and systemic pathways allows cell-specific and isoform-specific RNAi in Drosophila. RNA. 2003, 9: 299-308.

    CAS  PubMed Central  PubMed  Google Scholar 

  87. Kobayashi I, Tsukioka H, Komoto N, Uchino K, Sezutsu H, Tamura T, Kusakabe T, Tomita S: SID-1 protein of Caenorhabditis elegans mediates uptake of dsRNA into Bombyx cells. Insect Biochem Mol Biol. 2012, 42: 148-154.

    CAS  PubMed  Google Scholar 

  88. Kloc A, Martienssen R: RNAi, heterochromatin and the cell cycle. Trends Genet. 2008, 24: 511-517.

    CAS  PubMed  Google Scholar 

  89. Grewal SI, Moazed D: Heterochromatin and epigenetic control of gene expression. Science. 2003, 301: 798-802.

    CAS  PubMed  Google Scholar 

  90. Lippman Z, Martienssen R: The role of RNA interference in heterochromatic silencing. Nature. 2004, 431: 364-370.

    CAS  PubMed  Google Scholar 

  91. Ogawa Y, Sun BK, Lee JT: Intersection of the RNA interference and X-inactivation pathways. Science. 2008, 320: 1336-1341.

    CAS  PubMed Central  PubMed  Google Scholar 

  92. Martienssen RA, Zaratiegui M, Goto DB: RNA interference and heterochromatin in the fission yeast Schizosaccharomyces pombe. Trends Genet. 2005, 21: 450-456.

    CAS  PubMed  Google Scholar 

  93. Doyon Y, Selleck W, Lane WS, Tan S, Cote J: Structural and functional conservation of the NuA4 histone acetyltransferase complex from yeast to humans. Mol Cell Biol. 2004, 24: 1884-1896.

    CAS  PubMed Central  PubMed  Google Scholar 

  94. Dudley NR, Labbe JC, Goldstein B: Using RNA interference to identify genes required for RNA interference. Proc Natl Acad Sci U S A. 2002, 99: 4191-4196.

    CAS  PubMed Central  PubMed  Google Scholar 

  95. Anantharaman V, Koonin EV, Aravind L: Comparative genomics and evolution of proteins involved in RNA metabolism. Nucleic Acids Res. 2002, 30: 1427-1464.

    CAS  PubMed Central  PubMed  Google Scholar 

  96. Izaurralde E: A role for eIF4AII in microRNA-mediated mRNA silencing. Nat Struct Mol Biol. 2013, 20: 543-545.

    CAS  PubMed  Google Scholar 

  97. Meijer HA, Kong YW, Lu WT, Wilczynska A, Spriggs RV, Robinson SW, Godfrey JD, Willis AE, Bushell M: Translational repression and eIF4A2 activity are critical for microRNA-mediated gene regulation. Science. 2013, 340: 82-85.

    CAS  PubMed  Google Scholar 

  98. Sweeney SJ, Campbell P, Bosco G: Drosophila sticky/citron kinase is a regulator of cell-cycle progression, genetically interacts with Argonaute 1 and modulates epigenetic gene silencing. Genetics. 2008, 178: 1311-1325.

    CAS  PubMed Central  PubMed  Google Scholar 

  99. Kloc A, Zaratiegui M, Nora E, Martienssen R: RNA interference guides histone modification during the S phase of chromosomal replication. Curr Biol. 2008, 18: 490-495.

    CAS  PubMed Central  PubMed  Google Scholar 

  100. Stoica C, Carmichael JB, Parker H, Pare J, Hobman TC: Interactions between the RNA interference effector protein Ago1 and 14-3-3 proteins: consequences for cell cycle progression. J Biol Chem. 2006, 281: 37646-37651.

    CAS  PubMed  Google Scholar 

  101. Stoica C, Park J, Pare JM, Willows S, Hobman TC: The Kinesin motor protein Cut7 regulates biogenesis and function of Ago1-complexes. Traffic. 2010, 11: 25-36.

    CAS  PubMed  Google Scholar 

  102. Li F, Goto DB, Zaratiegui M, Tang X, Martienssen R, Cande WZ: Two novel proteins, dos1 and dos2, interact with rik1 to regulate heterochromatic RNA interference and histone modification. Curr Biol. 2005, 15: 1448-1457.

    CAS  PubMed  Google Scholar 

  103. Bologna NG, Voinnet O: The diversity, biogenesis, and activities of endogenous silencing small RNAs in arabidopsis. Annu Rev Plant Biol. 2014, 65: 473-503.

    CAS  PubMed  Google Scholar 

  104. Dalzell JJ, McVeigh P, Warnock ND, Mitreva M, Bird DM, Abad P, Fleming CC, Day TA, Mousley A, Marks NJ, Maule AG: RNAi effector diversity in nematodes. PLoS Negl Trop Dis. 2011, 5: e1176-

    CAS  PubMed Central  PubMed  Google Scholar 

  105. Nicolas FE, Torres-Martinez S, Ruiz-Vazquez RM: Loss and retention of RNA interference in fungi and parasites. PLoS Pathog. 2013, 9: e1003089-

    CAS  PubMed Central  PubMed  Google Scholar 

  106. Sandoval PY, Swart EC, Arambasic M, Nowacki M: Functional diversification of Dicer-like proteins and small RNAs required for genome sculpting. Dev Cell. 2014, 28: 174-188.

    CAS  PubMed  Google Scholar 

  107. Cecere G, Hoersch S, O’Keeffe S, Sachidanandam R, Grishok A: Global effects of the CSR-1 RNA interference pathway on the transcriptional landscape. Nat Struct Mol Biol. 2014, 21: 358-365.

    CAS  PubMed Central  PubMed  Google Scholar 

  108. Keeling CI, Chiu CC, Aw T, Li M, Henderson H, Tittiger C, Weng HB, Blomquist GJ, Bohlmann J: Frontalin pheromone biosynthesis in the mountain pine beetle, Dendroctonus ponderosae, and the role of isoprenyl diphosphate synthases. Proc Natl Acad Sci U S A. 2013, 110: 18838-18843.

    CAS  PubMed Central  PubMed  Google Scholar 

  109. Rajagopal R, Sivakumar S, Agrawal N, Malhotra P, Bhatnagar RK: Silencing of midgut aminopeptidase N of Spodoptera litura by double-stranded RNA establishes its role as Bacillus thuringiensis toxin receptor. J Biol Chem. 2002, 277: 46849-46851.

    CAS  PubMed  Google Scholar 

  110. Sivakumar S, Rajagopal R, Venkatesh GR, Srivastava A, Bhatnagar RK: Knockdown of aminopeptidase-N from Helicoverpa armigera larvae and in transfected Sf21 cells by RNA interference reveals its functional interaction with Bacillus thuringiensis insecticidal protein Cry1Ac. J Biol Chem. 2007, 282: 7312-7319.

    CAS  PubMed  Google Scholar 

  111. Sakurai K, Amarzguioui M, Kim DH, Alluin J, Heale B, Song MS, Gatignol A, Behlke MA, Rossi JJ: A role for human Dicer in pre-RISC loading of siRNAs. Nucleic Acids Res. 2011, 39: 1510-1525.

    CAS  PubMed Central  PubMed  Google Scholar 

  112. Marques JT, Kim K, Wu PH, Alleyne TM, Jafari N, Carthew RW: Loqs and R2D2 act sequentially in the siRNA pathway in Drosophila. Nat Struct Mol Biol. 2010, 17: 24-30.

    CAS  PubMed Central  PubMed  Google Scholar 

  113. Williams RW, Rubin GM: ARGONAUTE1 is required for efficient RNA interference in Drosophila embryos. Proc Natl Acad Sci U S A. 2002, 99: 6889-6894.

    CAS  PubMed Central  PubMed  Google Scholar 

  114. Zhu L, Masaki Y, Tatsuke T, Li Z, Mon H, Xu J, Lee JM, Kusakabe T: A MC motif in silkworm Argonaute 1 is indispensible for translation repression. Insect Mol Biol. 2013, 22: 320-330.

    CAS  PubMed  Google Scholar 

  115. Kawaoka S, Hayashi N, Suzuki Y, Abe H, Sugano S, Tomari Y, Shimada T, Katsuma S: The Bombyx ovary-derived cell line endogenously expresses PIWI/PIWI-interacting RNA complexes. RNA. 2009, 15: 1258-1264.

    CAS  PubMed Central  PubMed  Google Scholar 

  116. Czech B, Preall JB, McGinn J, Hannon GJ: A transcriptome-wide RNAi screen in the Drosophila ovary reveals factors of the germline piRNA pathway. Mol Cell. 2013, 50: 749-761.

    CAS  PubMed Central  PubMed  Google Scholar 

  117. Kasim V, Wu S, Taira K, Miyagishi M: Determination of the role of DDX3 a factor involved in mammalian RNAi pathway using an shRNA-expression library. PLoS ONE. 2013, 8: e59445-

    CAS  PubMed Central  PubMed  Google Scholar 

  118. Pitino M, Coleman AD, Maffei ME, Ridout CJ, Hogenhout SA: Silencing of aphid genes by dsRNA feeding from plants. PLoS ONE. 2011, 6: e25709-

    CAS  PubMed Central  PubMed  Google Scholar 

  119. Bansal R, Michel AP: Core RNAi machinery and Sid1, a component for systemic RNAi, in the Hemipteran insect, Aphis glycines. Int J Mol Sci. 2013, 14: 3786-3801.

    CAS  PubMed Central  PubMed  Google Scholar 

  120. Pillai RS, Bhattacharyya SN, Artus CG, Zoller T, Cougot N, Basyuk E, Bertrand E, Filipowicz W: Inhibition of translational initiation by Let-7 MicroRNA in human cells. Science. 2005, 309: 1573-1576.

    CAS  PubMed  Google Scholar 

  121. Sen GL, Blau HM: Argonaute 2/RISC resides in sites of mammalian mRNA decay known as cytoplasmic bodies. Nat Cell Biol. 2005, 7: 633-636.

    CAS  PubMed  Google Scholar 

  122. Negre V, Hotelier T, Volkoff AN, Gimenez S, Cousserans F, Mita K, Sabau X, Rocher J, Lopez-Ferber M, d’Alencon E, Audant P, Sabourault C, Bidegainberry V, Hilliou F, Fournier P: SPODOBASE: an EST database for the lepidopteran crop pest Spodoptera. BMC Bioinformatics. 2006, 7: 322-

    PubMed Central  PubMed  Google Scholar 

  123. Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH: CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39: D225-229.

    CAS  PubMed Central  PubMed  Google Scholar 

  124. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739.

    CAS  PubMed Central  PubMed  Google Scholar 

  125. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504.

    CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

We are thankful to Bionivid Technology Pvt. Ltd., Bangalore, India; especially Mr. Rohit Shukla, for his support during in silico analysis. This work was supported by a financial grant (BT/PR10673/AGR/36/579/2008) from the Department of Biotechnology, Govt. of India, India.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Pawan Malhotra, Sunil K Mukherjee or Raj K Bhatnagar.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RKB, SKM and PM proposed and designed the experiments. SG and PKK performed the experiments and AK analysed the bioinformatics data and phylogenetic tree construction. SG, PM, SKM and RKB wrote the paper. All authors read and approved the final manuscript.

Pawan Malhotra, Sunil K Mukherjee and Raj K Bhatnagar contributed equally to this work.

Electronic supplementary material

Additional file 1: In silico predicted gene list of eighty putative candidates for Sf21 RNAi pathway.(DOCX 28 KB)

12864_2014_6685_MOESM2_ESM.docx

Additional file 2: Real-Time analysis to validate gene knockdown by measuring efficiency of siRNA of putative candidates.(DOCX 18 KB)

12864_2014_6685_MOESM3_ESM.docx

Additional file 3: Table containing % gfp expression obtained from post transfected FACS analysis in the functional assay for siRNA treated both gfp reverted and non-reverted genes.(DOCX 21 KB)

Additional file 4: List of siRNAs used for gfp reversion assay.(DOCX 19 KB)

Additional file 5: List of siRNAs used to validate transfection efficiency of selected genes.(DOCX 15 KB)

Additional file 6: List of primers used for qRT-PCR.(DOCX 15 KB)

Additional file 7: Accession numbers used for identity search of RNAi gene homologs with other insects.(DOCX 20 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ghosh, S., Kakumani, P.K., Kumar, A. et al. Genome wide screening of RNAi factors of Sf21 cells reveal several novel pathway associated proteins. BMC Genomics 15, 775 (2014). https://doi.org/10.1186/1471-2164-15-775

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-775

Keywords