Knockdown of Genes Involved in Transcription and Splicing Reveals Novel RNAi Targets for Pest Control

RNA interference (RNAi) is a promising next generation technology for the development of species-specific pest management. The key to successful RNAi based-plant protection is dependent in part on data-driven target gene selection, a challenging task due to the absence of laboratory strains and the seasonality of most pest species. In this study, we aimed to identify novel target genes by performing a knowledge-based approach in order to expand the spectrum of known potent RNAi targets. Recently, the protein-coding genes ncm, Rop, RPII-140, and dre4 have been identified as sensitive RNAi targets for pest control. Based on these potent RNAi targets, we constructed an interaction network and analyzed a selection of 30 genes in the model beetle Tribolium castaneum via injection of dsRNA synthesized by in vitro transcription. Nineteen of these targets induced significant mortality of over 70%, including six that caused 100% lethality. Orthologs of active T. castaneum RNAi targets were verified in the economically important coleopteran pests Diabrotica virgifera virgifera and Brassicogethes aeneus. Knockdown of D. v. virgifera genes coding for transcription factor Spt5, Spt6, and RNA polymerase II subunit RPII-33 caused over 90% mortality in larval feeding assays. Injection of dsRNA constructs targeting RPII-215 or the pre-mRNA-processing factor Prp19 into adult B. aeneus resulted in high lethality rates of 93 and 87%, respectively. In summary, the demonstrated knowledge-based approaches increased the probability of identifying novel lethal RNAi target genes from 2% (whole genome) to 36% (transcription- and splicing-related genes). In addition, performing RNAi pre-screening in a model insect increased also the probability of the identification essential genes in the difficult-to-work-with pest species D. v. virgifera and B. aeneus.


INTRODUCTION
The world population is expected to reach ∼ 10 billion people in 2050, necessitating a substantial crop yield increase to meet the global food demand (Johnson and Jones, 2017;Rohr et al., 2019). An area where significant productivity gains can be made is the reduction of crop losses associated with insect pests, which is estimated at ca. 15% at a multicrop and worldwide scale (Oerke, 2006;Fletcher et al., 2020). The primary solution for insect control, the use of chemical pesticides, is facing challenges like resistance development and growing concerns of undesirable effects on the environment or non-target organisms. Genetically modified (GM) crops expressing insecticidal Bacillus thuringiensis (Bt) proteins provided a technology improvement on pest management that reduced the dependence on chemical insecticides (Phipps and Park, 2002;James, 2009;Areal and Riesgo, 2015). However, wide adoption of Bt trait technology has resulted in field-evolved resistance (Tabashnik et al., 2009(Tabashnik et al., , 2013. Consequently, new pest control strategies to overcome these obstacles are urgently required. A promising alternative insect pest control strategy is the use of double-stranded RNAs (dsRNAs) to knock-down essential pest genes by triggering the conserved eukaryotic RNA interference (RNAi) pathway (Huvenne and Smagghe, 2010;Xue et al., 2012). Upon entry of gene-specific dsRNA into the cell of the intended pest organism, the expression of the targeted endogenous gene transcript will be knocked down post-transcriptionally and can result in insect mortality. Major advantages of this new approach are the rapid biological degradation of dsRNA in the environment and the sequencespecific nature of RNAi, allowing selective pest control that protects non-target organisms like pollinators or predatory insects (Dubelman et al., 2014;Albright et al., 2017;Parker et al., 2019). Baum et al. (2007) demonstrated that transgenic corn expressing a vacuolar ATPase subunit showed protection against feeding damage by Diabrotica virgifera virgifera LeConte (Coleoptera: Chrysomelidae) infestation. Recent examples of successful host-induced gene silencing (HIGS) applications include potato plants resistant to Colorado potato beetles (Zhang et al., 2015), as well as resistant maize plants silencing additional RNAi targets of D. v. virgifera (Hu et al., 2016) (reviewed in Koch and Kogel, 2014). The first commercialized transgenic event that uses HIGS for pest control, jointly developed by Monsanto and Dow AgroSciences, was approved by the three US regulatory agencies (EPA, FDA, and USDA). The combined-trait product known as SmartStax PRO R expresses both Bt toxins and dsRNA targeting the D. v. virgifera Snf7 gene (Head et al., 2017). As RNAi is a very active research area and an emerging pest management tool with potential for plant protection, it is just a question of time before SmartStax Pro and other applied RNA suppression products will appear on the market.
A key step in RNAi based pest control is the identification of novel bioinsecticide targets, a challenging task, as many pests lack genomic data and are often not easy to rear under laboratory conditions (Vogel et al., 2014). Performing large-scale RNAi prescreens in a model organism was often proposed as a solution for target-site validation in pest insects (Alves et al., 2010;Knorr et al., 2013;Ulrich et al., 2015). Recently, a high throughput RNAi screen functionally analyzed 5000 genes of the model beetle Tribolium castaneum, revealing a panel of 100 lethal target genes that caused ≥ 90% mortality (Dönitz et al., 2014;Ulrich et al., 2015). A selection of these lethal T. castaneum RNAi targets was used to test for efficacy of orthologous genes in D. v. virgifera and the oilseed rape pest Brassicogethes aeneus (previously known as Meligethes aeneus). This approach identified the highly active target genes DNA-directed RNA polymerase II subunit (RPII-140), dre4 (homolog of SPT16), and ncm or nucampholin (homolog of CWC22) (Knorr et al., 2018). These sensitive RNAi targets encode polypeptides with essential functions at different steps of transcription and splicing. More specifically, dre4 is part of the FACT (facilitates chromatin transactions) complex that acts as a general chromatin regulator (Belotserkovskaya and Reinberg, 2004;Yang et al., 2016). The core splicing factor ncm is directing eukaryotic translation initiation factor 4A3 (eIF4AIII) into the exon junction complex (EJC) (Barbosa et al., 2012) and RPII-140 catalyzes the transcription of DNA into RNA (Falkenburg et al., 1987).
The primary goal of the current work was to expand on these potent T. castaneum RNAi targets by silencing functionally connected genes. Target gene selection was based on a knowledge based approach supported by functional Protein-Protein Interaction (PPI) networks generated by the STRING search tool (Mering et al., 2005;Szklarczyk et al., 2015) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway information (Ogata et al., 1999). We screened 30 genes in T. castaneum and identified some potent novel RNAi targets that were active in both larval and adult stage. In addition, we evaluated if the newly discovered RNAi targets could be leveraged to the economically important pest species, the Western Corn Rootworm (WCR) D. v. virgifera, a devastating pest of maize in the US Corn Belt, and the pollen beetle B. aeneus, a key pest of oilseed rape. For both of these difficult-to-work with agricultural pest insects we identified novel RNAi target genes that caused more than 90% mortality, suggesting that these genes may be considered as potential RNAi target for RNA-based management.

RNA Isolation and Complementary DNA Synthesis
Total RNA was extracted from T. castaneum larvae and B. aeneus adults using Direct-zol TM RNA MiniPrep (Zymo Research, Irvine, California, USA) according to the manufacturer's instructions. cDNA was synthesized from 500 ng RNA using the First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA).
Total RNA from D. v. virgifera eggs, larvae, and adults was extracted using TRIzol (Life Technologies, Grand Island, NY) according to the manufacturer's instructions. cDNA was synthesized from 1 µg of total RNA using a SuperScript III reverse transcription kit (Thermo Fisher Scientific, Waltham, MA). Orthologs in D. v. virgifera and B. aeneus The T. castaneum genome sequence ver.Tcas5.2 was used for the identification of ortholog proteins of the transcriptome of B. aeneus generated by Vogel et al. (2014) using NCBI BLASTP with an E-value of 0.01. The output were further filtered by a minimum protein identity of 50% as well as a protein coverage of at least 75% and lower than 125%. Resulting hits were ranked by score; redundant and overlapping sequences were removed. D. v. virgifera homologous were identified by performing TBLASTN searches with candidate protein coding sequences that were run against BLASTable databases comprising the unassembled D. v. virgifera sequence reads or the assembled contigs. Significant hits to a D. v. virgifera sequence (defined as lower than <1 × 10-20 for contig homologies and better than E-value of <1 × 10-10 for unassembled sequence reads homologies) were verified using BLASTX against the NCBI non-redundant database. The results of this BLASTX search confirmed that the D. v. virgifera homolog candidate gene sequences identified in the TBLASTN search indeed comprised D. v. virgifera genes, or were the best hits to the non D. v. virgifera candidate gene sequence present in the D. v. virgifera sequences. In the majority of cases, T. castaneum candidate genes that were annotated as encoding a protein, showed clear sequence homology to a sequence or sequences within the D. v. virgifera transcriptome. Sequencher TM v4.9 (Gene Codes Corporation, Ann Arbor, MI) was used in a few cases to assembled contigs that were partially-overlapping into longer contigs.

dsRNA Synthesis
The selected dsRNA sequences were screened against genomic and transcriptomic data to avoid complete matches of more than 19 nt to a potential off-target genes. Gene-specific primers for PCR (Supplementary Table 1) amplification were designed using VectorNTI (Invitrogen, Carlsbad, CA) or Primer3 software and contained a T7 promoter sequence (5 ′ -TAATACGACTCACTATAGGGAGA-3 ′ ) at their 5 ′ ends. The dsRNAs were synthesized using Ambion MEGAscript T7 (Thermo Fisher Scientific, Waltham, MA) according to the manufacturer's protocol and quantified on a NanoDrop 8000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA).

Injection bioassays
Dorsolateral injection of 150 nl dsRNA [250 ng/µl] into T. castaneum larvae and B. aeneus adults was conducted using a glass capillary and a micromanipulator M3301 (World Precision Instruments, Sarasota, FL) under a dissecting stereomicroscope (n = 12). Negative controls received an equivalent amount of water or dsRNA corresponding to the Galleria mellonella metalloproteinase inhibitor (IMPI) gene (not present in B. aeneus) (Supplementary Table 1). Injected T. castaneum larvae were kept on whole-grain flour, whereas B. aeneus beetles were supplied with dried pollen and wet tissues as water source. Survival rates were monitored and insects counted every 48 h for 10 days.

Feeding bioassays
Brassicogethes aeneus and D. v. virgifera feeding bioassays were performed as described previously (Knorr et al., 2018). In brief, B. aeneus adults ingested dsRNA [500 ng/cm 2 ] using artificial diet to ensure continuous uptake of dsRNA. Negative controls received an equivalent amount of water or dsRNA corresponding to the G. mellonella metalloproteinase inhibitor (IMPI) gene (Supplementary Table 1). The number of dead beetles was counted every 48 h for 10 days.
Neonate D. v. virgifera larvae were fed on artificial diet that was overlaid with dsRNA (500 ng/cm 2 ) (Knorr et al., 2018). Water, 0.1X TE buffer and dsRNA targeting yellow fluorescent protein gene (YFP) were used as negative controls. The number of dead insects, and the weights of surviving insects were recorded after 9 days. Growth Inhibition (GI) was calculated based on the average weights of all controls, as follows: GI = [1 -(TWIT/TNIT)/(TWIBC/TNIBC)], where TWIT is the total weight of insects in the treatment, TNIT is the total number of insects in the treatment, TWIBC is the total weight of insects in the buffer control, and TNIBC is the total number of insects in the buffer control. To estimate the 50% lethal concentration (LC50) and the concentration that causes 50% growth inhibition (GI50) of active RNAi targets, 4-fold serial dilutions of dsRNAs were assayed. Values were calculated using log-logistic regression analysis in JMP R Pro from SAS Institute Inc (Supplementary Table 1).

Statistical Analysis
Bioassays were analyzed using an analyses of variance (ANOVA), followed by a Bonferroni-Holm test with significance threshold of p < 0.05 using Daniel's XL toolbox for Excel, version 7.3.4 68. Each experiment was compared to a control group and all experiments were conducted independently at least three times. LC50 (concentration at which 50% of the insects are dead) and GI50 (concentration that causes 50% growth inhibition or GI) values were calculated as described earlier (Knorr et al., 2018).

Knockdown of Potential Target Genes in T. castaneum via Injection of dsRNA
The three most active RNAi targets RPII-140, Dre4, and ncm described in Knorr et al. (2018) are involved in pre-and posttranscriptional mRNA processing. We mined corresponding String PPI networks and KEGG pathway for potential RNAi targets and selected 30 genes for T. castaneum dsRNA injection assays: 18 genes coding for different RNA polymerase subunits (RP genes) and 12 genes involved in various transcriptional processes (non-RP genes) (Table 1, Figure 1). Homologs of the selected 30 genes were identified via BLAST (tblastn) homology searches using Drosophila melanogaster protein isoforms (see Supplementary Table 1 for GenBank IDs) as the queries. To test whether dsRNA targeting these genes may lead to lethal phenotypes in T. castaneum, we designed dsRNA constructs targeting a segment of the open reading frame and applied the dsRNA via injection into the larval hemocoel. Significant mortality of ≥ 50% was observed onset on the 4th day after injection of ∼150 nl of 250 ng/µl dsRNA targeting the four genes cdk8, Prp8, Spt5, and Spt6 (Figure 2A). The detected mortality rates were significantly higher compared to the non-target IMPI dsRNA. Significant mortality of over 70% was caused by 19 out of the 30 applied dsRNA constructs 10 days post injection (dpi). The highest mortality of 100% was recorded for the six genes Spt5, Spt6, RPII-16, RPIII-18, RPII-33, and RPII-215 (Table 1, Figure 2). In addition, silencing RNA polymerase subunits (RPs) caused over 50% larval mortality in all 18 treatments 10 dpi (Table 1, Figure 2B).
We also analyzed the potential lethal effect of the 18 RP dsRNA constructs in T. castaneum female and male adults. Significant mortality rates of over 50% were observed for 14 female and 16 male dsRNA RP constructs 10 dpi (Table 1,  Figure 3), including six female and eight treatments that caused over 90% mortality, respectively (Table 1, Figure 3). RPII-215 was one of the most efficacious RNA targets, causing over 75% in females and even 97% mortality in males 6 dpi. Only knockdown of RPI-1 and RPIII-32 resulted in considerably lower female mortality rates of ∼7 and ∼18%, respectively. Whereas, in males silencing of RPI-1 led to >50% and RPIII-32 to >90% lethality.
Frontiers in Agronomy | www.frontiersin.org FIGURE 1 | Simplified representation of eukaryotic gene expression and post-transcriptional processing, highlighting the genes that were selected as potential RNAi targets. Gene transcription in eukaryotic cells is carried out by three multisubunit RNA polymerases. RNA polymerase I, II, and III (RNAPI, RNPII, and RNAPIII, respectively). RNAPI synthesizes the larger ribosomal RNA (rRNAs) precursor 28S, 18S, and 5.8S, whereas RNAPIII synthesizes 5S rRNAs, transfer RNAs (tRNAs), and small non-coding RNAs sRNAs (Wild and Cramer, 2012;Yokoyama, 2019). Processed rRNAs form the pre-40S and pre-60S ribosomal subunits, which assemble to mature ribosomes in the cytoplasm. RNAPII transcribes sRNAs and precursor messenger RNA (pre-mRNA). The largest RNAPII subunit, the C-terminal domain (CTD), is flexible in structure and undergoes several conformational changes during gene transcription. These changes are regulated, in part, through alterations of the phosphorylation status that are mediated by phosphatases like SSU72 (Liu et al., 2020). The transcription initiation factor TFIID is a multimeric protein complex and plays a central role in mediating promoter responses to various activators and repressors (Lewis and Reinberg, 2003;Guo et al., 2018). Cyclin-dependent kinase 8 (Cdk8) is part of the mediator complex, a coactivator that is involved in regulating the gene transcription (Leclerc et al., 1996;Yin and Wang, 2014). Transcription elongation factor SPT6 enhances transcription elongation by binding to histone H3 and recruitment of other elongation factors like Spt5 (Ardehali et al., 2009;Duina, 2011). Spt5 is a component of the DRB sensitivity-inducing factor complex (DSIF complex) that both inhibits and stimulates transcription elongation (Jennings et al., 2004). Introns of the pre-mRNA are removed in a two-step transesterification reaction by a large ribonucleoprotein complex called the splicosome. Introns of the pre-mRNA are removed in a two-step transesterification reaction by a large ribonucleoprotein complex called the splicosome. The catalytic core of the spliceosome contains the pre-mRNA processing factor 8 (Prp8), which plays a central role in its molecular interactions during spliceosome formation (Grainger and Beggs, 2005). The pre-mRNA splicing factor Cwc25 is required for the first catalytic reaction within the splicosome (Chiu et al., 2009). The splicing factor 3B subunit 4 (SF3B4) encodes a core subunit of the SF3b complex, a part of the U2-type spliceosome (Xiong and Li, 2020). Structural changes of the spliceosome that allow the second reaction are mediated by the DEAH/H-box ATP-dependent RNA helicase Prp16 (Tseng et al., 2011). The large protein complex associated with the DEAH-box ATPase Prp19, named NTC (NineTeen Complex), stabilizes interactions of the activated splicosome (Tseng et al., 2011). The eukaryotic initiation factor 4A-III (eIF4AIII) is another DEAH-box ATP-dependent RNA helicase and a core component of the splicing-dependent multiprotein exon junction complex (EJC) that forms at the junction of two exons (Barbosa et al., 2012;Le Hir et al., 2015). The mature mRNAs are exported to the cytoplasma, were eukaryotic translation initiation factors (eIFs), like the eIF3 complex, stimulate their recruitment to the 40S ribosomal subunit for translation into protein (Poulin andSonenberg, 2013, Valášek et al., 2017).

Knockdown of Potential Target Genes in D. v. virgifera Larvae Upon Feeding
The four most lethal target genes Spt5, Spt6, RPII-215, and RPII-33, as well as the six target genes Prp8, Prp19, Cwc25, RPI-1, RPII-15, and RPIII-128 with moderate activity, were selected from the T. castaneum screen for D. v. virgifera knockdown experiments. Diet overlay bioassays using 500 ng/cm 2 dsRNA identified seven out of the 10 dsRNA constructs that caused significant mortality of over 50%. Seven of the tested target genes, caused over 70% mortality (Figure 4A), including dsRNAs targeting Spt5, Spt6, and RPII-33 transcripts, which resulted in the highest mortality rates of 99, 90, and 88%, respectively. In addition, RNAi treated groups showed significant growth inhibition in D. v. virgifera ranging from 0.72 to 0.99 in all treatments, except for the Cwc25 treatment ( Figure 4B, Supplementary Table 1).
In addition, the median lethal concentration (LC 50 ) and median growth inhibition (GI 50 ) values for dsRNA treatments that caused > 60% mortality were further investigated. A dilution series of dsRNA tested in diet-overlay bioassays led to LC 50 values that ranged from 4.33 to 54.26 ng/cm 2 ( Table 2). The majority of dsRNAs tested resulted in LC 50

Knockdown of Potential Target Genes in B. aeneus Adults via Injection of dsRNA
The four most lethal target genes Spt5, Spt6, RPII-215, and RPII-33, as well as Prp19 and RPIII-128 as moderately potent target genes were selected from the T. castaneum screen for B. aeneus  Table 1). knockdown experiments. The ortholog sequences were identified from the B. aeneus transcriptome and gene-specific dsRNAs were designed. Injection of ∼150 nl of 250 ng/µl into adult beetles caused significant mortality in five of the six tested dsRNA constructs 10 dpi, including high rates of 87% for the Prp19 and 94% for the RPII-215 dsRNA-treated groups (Figure 5). The most rapid response was observed after silencing RPII-215, which caused significant mortality at four dpi, which increased to over 50% six dpi and reaching 87% mortality on day eight of the bioassay (Figure 5).

DISCUSSION
Insect damage to crops along with climate change and the loss of arable land dedicated to agriculture threaten the ability to meet the food demand of a rapidly growing global population (FAO, 2017). Chemical insecticides are still the major control method for insect pests, but widespread resistance development and environmental concerns like non-target toxicity create the need for innovative and sustainable pest control methods (Tabashnik, 1994;Budzinski and Couderchet, 2018). RNAibased crop protection is a new-generation technology for highly efficient, environmentally soundbiopesticides that target insect pests with great specificity. To date, insecticidal RNAi target genes have been identified in key biological systems of the insect physiology like digestion, defense mechanisms, vesicular trafficking, and detoxification (Mao et al., 2007;Bautista et al., 2009;Kim et al., 2015;Kola et al., 2015;Head et al., 2017;Knorr et al., 2018;Cooper et al., 2019). A major challenge for the development of RNAi applications, however, is the identification of novel sensitive target genes, in particular in non-model organisms that often lack genomic data or laboratory colonies (Vogel et al., 2014). Some researchers take advantage of known RNAi-sensitive genes in model organisms and perform targeted gene mining strategies to verify the orthologous RNAi target in other insects (Pitino et al., 2011;Wang et al., 2011;Zhao et al., 2011;Zhu et al., 2011;Li et al., 2013;Rodrigues et al., 2017;Mogilicherla et al., 2018;Koo et al., 2020). Recently, a knowledgebased approach demonstrated that the spectrum of known lethal RNAi targets could be expanded by silencing genes that are functionally connected (Bingsohn et al., 2017).
In this study, we mined the biological networks of RNA transcription and RNA processing, as silencing of genes with essential functions in these processes have demonstrated high mortality rates in three different coleopteran pests (Knorr et al., 2018). It has been shown that mRNA knockdown efficiency is not necessarily positively related to high mortality rates. Several uncontrollable factors can influence the RNAi response, like protein/mRNA half-life, gene function redundancy or the physiological importance of the targeted gene function (Schwinghammer et al., 2011;Li et al., 2013;Rinkevich and Scott, 2013;Bingsohn et al., 2017;Raje et al., 2018;Howell et al., 2020;Koo et al., 2020). Thus, we focused on organism-level analysis and did not consider gene expression as a key criterion for target selection. Out of 30 genes assayed in T. castaneum, 19 genes caused significant mortality of over 70%, including 11 target genes that lead to >90% mortality at day six post injection. Furthermore, silencing of RNA polymerase subunits in T. castaneum larvae caused significant mortality of at least 50% mortality in all 18 treatments and even 100% when knocking down RPII-215, RPII-33, RPII-18, and RPII-16. RNA polymerases subunits are also sensitive targets in adult beetles, both female and males, as knock down of seven subunits (RPI-1, RPII-10 RPII-11 RPII-15 RPII-16 RPII-215, and RPIII-1) led to 100% mortality 10 days after dsRNA injection. Taken together, 36% of the silenced transcription-and splicing-related genes were identified as highly effective RNAi targets (≥90%. mortality). In comparison, only two percent of 5000 randomly screened Tribolium genes in  Table 1).
iBeetle showed mortality rates of over 90% (Dönitz et al., 2014;Ulrich et al., 2015). These results demonstrate that a knowledgebased approach can results in a higher probability of identifying efficient RNAi targets than high-throughput experiments.
To determine if the T. castaneum RNAi-sensitive targets could have a similar effect in other coleopteran pests, we identified orthologous sequences in D. v. virgifera and B. aeneus using transcriptome data mining. Significant growth inhibition was detected for all ten dsRNA constructs that were tested via artificial diet assays using neonate WCR larvae, out of which seven also caused significant mortality of over 50%. Silencing of the transcription elongation factors Spt5 and Spt6, two of the most active T. castaneum RNAi targets in this study, also caused high mortality rates of 99 and 90% in D. v. virgifera, respectively, and the lowest LC50 values of approximately 4 ng/cm 2 . Due to the lack of B. aeneus lab-colonies, pollen beetle experiments were limited to field collected beetles and focused on injection assays of six active T. castaneum target genes, namely Spt5, Spt6, Prp19, RPII-215, RPII-33, and RIII-128. Knockdown of Prp19 and RPII-215 led to high mortality rates of 87 and 94%, respectively. Silencing of RPII-215, resulted in the most rapid response of 50% mortality 6 days after treatment. Delivering dsRNA via injection is, of course, not a practical technique for RNAi field applications (Huvenne and Smagghe, 2010;Xue et al., 2012). Environmental RNAi response in B. aeneus was previously reported (Knorr et al., 2018;Willow et al., 2020a,b) but might lead to variable knockdown efficiency, which, in this initial study, we wanted to avoid by applying, standardized dsRNA doses to the beetles. Oral uptake of dsRNA might have different dynamics than injection and thus, lethal activity needs to be verified in feeding experiments before implementing results of this study in the greenhouse and the field.
Our study demonstrated a high sensitivity to RNAi of genes involved in RNA transcription and mRNA splicing, providing opportunities to screen for additional targets in these functional classes. In particular, silencing genes that are involved in RNAPII mediated mRNA transcription resulted in severe mortality, revealing a prime source for highly active RNAi targets. The Suppressor of Ty (SPT) genes Spt5 and Spt6, most active targets in T. castaneum and D. v. virgifera, encode polypeptides with essential functions in early transcription elongation. Spt5 forms with Spt4 the conserved heterodimer DSIF (DRB sensitivityinducing factor) that regulates processive RNAPII transcription elongation (Price, 2000;Winston, 2001;Jennings et al., 2004;Decker, 2020). The histone chaperone Spt6 assembles histones into nucleosomes and is able to escalate the RNAPII elongation rate (Ardehali et al., 2009;McDonald et al., 2010).
Synthetic agrochemicals have traditionally been the method of choice for pollen beetle control in oil seed rape (OSR), but their excessive application, particularly pesticides of the pyrethroid class, has led to the development of populations with high level of resistance (Zimmer et al., 2014). The decreased sensitivity led to the introduction of non-cross-resistant alternative insecticides, such as neonicotinoids. However, adverse effects on non-target organisms like bees and other pollinators resulted in the ban of several neonicotinoids (Bass and Field, 2018). RNAi mediated pollen beetle control via spray-induced gene silencing (SIGS), whereby target-specific dsRNA is sprayed onto host plants, represents a promising, ecologically safe alternative. Foliar applications of naked dsRNA have already shown to be active against different Coleopteran species (e.g., de Andrade and Hunter, 2016;San Miguel and Scott, 2016). But, OSR Brassica napus crops have a short vulnerability window of only a few weeks, as agronomic crop loss is only caused by adults feeding from flower buds (Williams et al., 2010;Hervé and Cortesero, 2016). Thus, dose response studies need to identify most sensitive target genes in order to develop fast acting dsRNA applications. Formulation of dsRNA might further improve potential HIGS applications. dsRNA formulated with nanoparticles like BioClay or Guanylate Polymers demonstrated increased RNAi efficiency by overcoming problems like. stability or cell uptake (Mitter et al., 2017;Christiaens et al., 2018;Parsons et al., 2018). Behavioral studies have shown a preference of B. aeneus for turnip rape (Brassica rapa) over oilseed rape, which can be used for a push-pull strategy (Vinatier et al., 2012;Hervé and Cortesero, 2016;Skellern and Cook, 2018). Push-pull strategies involves integration of stimuli that deter insect pests away from their main crop (push), while simultaneously luring them to trap plants (pull) (Pyke et al., 1987;Cook et al., 2007). Combining HIGS with trap cropping would reduce the treated field area, leading to a more effective outcome of the dsRNA product.
RNAi represents a promising alternative for ecologically sustainable agricultural insect pest. Species-specific dsRNA based insecticides might complement other biosafe measures (e.g., conservation biocontrol) for managing agricultural pests and could become a part of integrated pest management (IPM) approaches. RNAi pre-screening in a model insect significantly increased the probability of identifying high levels of efficacious RNAi targets in pest species that lack genomic information. However, our data highlights the need to validate candidate targets in the insect species of concern because there is not a known correlation of gene RNAi sensitivity across insect species. Of the analyzed genes, 30-40% caused over 90% mortality in D. v. virgifera and B. aeneus. In addition, knockdown of RPII-215 resulted in significant mortality in all three coleopteran species: 100% in T. castaneum, 93% in B. aeneus and 79% in D. v. virgifera assays. These findings suggest that the function of these selected genes might be conserved across coleopterans. However, silencing of Prp19 resulted in high mortality rates of 81% in T. castaneum and 87% B. aeneus, but not in D. v. virgifera. Knockdown of Spt5, Spt6, and RPII-33 in T. castaneum was 100% lethal in all three cases and also caused significant mortality rates of 99, 90, and 98%, respectively, in D. v. virgifera. Targeting RPII-33 in B. aeneus did not led to any significant mortality rates, but silencing of Spt5 and Spt6 did, albeit to a lesser extent.
To summarize, we demonstrate that a knowledge-based approach increased the efficiency in discovering potent RNAi targets in Tribolium. Furthermore, identified T. castaneum RNAi targets were also highly efficacious in two additional coleopteran pest species, D. v. virgifera and B. aeneus, representing promising candidates for RNAi-based pest control. Next steps should focus on greenhouse and field demonstration of RNAi candidates that might become a promising component in future IPM programs.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
EK, KN, and AV designed the study with input from other team members. MF, PG, and KA identified the D. v. virgifera target orthologs. AB identified the B. aeneus RNAi target orthologs. EK and LT designed and executed T. castaneum and B. aeneus experiments. MR and WL carried out D. v. virgifera diet bioassay. EK, EF, CG, AB, and KN interpreted the data. EK and KN composed the manuscript. All authors contributed to the article and approved the submitted version.

FUNDING
This work was funded by the Hessen State Ministry of Higher Education, Research and the Arts (HMWK) via the LOEWE Center for Insect Biotechnology and Bioresources. This research was also supported under a Fraunhofer IME-Corteva Agriscience joint research collaboration.