Mosharrof Mondal
Megan Carver
Judith K. Brown- 1School of Plant Sciences, The University of Arizona, Tucson, AZ, United States
- 2RNAissance Ag LLC, St. Louis, MO, United States
RNA interference (RNAi) has potential to become a major tool for integrated management of insect pests of agricultural crops based on sequence-specificity and low doses of rapidly biodegradable dsRNA. Deploying ‘environmental RNAi’ for control of insect vectors of plant pathogens is of increasing interest for combatting emerging plant diseases. Hemipteran insect vectors, including psyllids, are vascular feeders, making their development difficult to control specifically by targeting with pesticidal chemistries. Psyllids transmit “Candidatus Liberibacter solanacearum” the causal organism of potato zebra chip and tomato vein greening diseases, transmitted, respectively, by the potato or tomato psyllid (PoP). Until now, the optimal effective concentration(s) of double-stranded RNA (dsRNA) required for significant gene knockdown and RNAi persistence in PoP have not been determined. The objective of this study was to optimize RNAi in young PoP adults and 3rd instars for screening by oral delivery of dsRNAs. The minimal effective dsRNA concentrations required for robust knockdown and persistence were evaluated by delivering seven concentrations spanning 0.1 ng/μL to 500 ng/μL over post ingestion-access periods (IAP) ranging from 48 h to 12 days. The PoP gene candidates evaluated as targets were vacuolar ATPase subunit A, clathrin heavy chain, and non-fermenting protein 7, which were evaluated for knockdown by qPCR amplification. The minimum and/or the second most effective dsRNA concentration resulting in effective levels of gene knockdown was 100 ng/μL for all three targets. Higher concentrations did not yield further knockdown, indicating potential RISC saturation at the higher doses. Gene silencing post-IAP of 100 ng/μL dsRNA persisted for 3–5 days in adults and nymphs, with the PoP 3rd instar, followed by teneral and mature adults, respectively, exhibiting the most robust RNAi-response.
Introduction
RNA interference (RNAi) is a biological process that is conserved among eukaryotic organisms. The exploitation of RNAi as a biopesticide technology has potential to expand and possibly reform the repertoire of pesticide regimes available for management of insect pests that damage plants by feeding or serve as vectors of plant pathogens. Though inherently an anti-viral defense system, RNAi also contributes regulatory functions that modulate endogenous gene expression, transcriptionally or post-transcriptionally (Shabalina and Koonin, 2008; Dowling et al., 2017). Small RNAs are the primary effector molecules of the RNAi pathways. These 20–30 nucleotide (nt) long RNAs complementarily bind to the target transcript and trigger degradation of the transcript or blocking of protein synthesis using the transcript as template. Consequently, gene silencing, or RNAi, can be exploited to develop the next generation biopesticides, referred to as ‘environmental RNAi’ (Ivashuta et al., 2015; Mondal et al., 2021). In RNAi technology, long double-stranded RNAs (dsRNA) or single-stranded RNAs (ssRNA) introduced to the target organism are scavenged to the RNAi machinery, usually immediately after ingestion, which is responsible for producing small RNAs that trigger destruction of the target transcript or the inhibition of translation (Agrawal et al., 2003; Baum et al., 2007; Mondal et al., 2020).
The efficacy of RNAi has been shown to result in mortality or other phenotypes similar to that of some traditional chemical pesticides. However, the superior species-specificity and ease of biodegradation of RNAi renders it safer for use in a variety of environments (Christiaens et al., 2022). Also, with respect to human safety, RNAi is considered a superior alternative to chemical pesticides (Baum and Roberts, 2014). Given the desirable features of RNAi, commercial applications that exploit this technology are going to be available for pest management within the next couple of years (Head et al., 2017; Rodrigues et al., 2021). Nonetheless, the success of environmental RNAi biopesticides largely depends on the particular target organism because of variability in responsiveness, or “penetrance”, of the RNAi triggering molecule(s) (Baum and Roberts, 2014). Insects classified in the orders: Blattodea, Coleoptera, and Orthoptera are known to exhibit a robust response, or ‘effective RNAi penetrance’, while dipterans and lepidopterans have been shown to exhibit some extent of recalcitrance to RNAi (Cooper et al., 2019; Zhu and Palli, 2020). Among hemipteran insects, RNAi efficacy has been demonstrated for the genera: Bemisia, Cimex, Halyomorpha, and Oncopeltus, whereas many Sternorrhyncha, including aphids and psyllids, exhibit relatively low RNAi penetrance (Moriyama et al., 2016; Ghosh et al., 2017; Mondal et al., 2020; Zhu and Palli, 2020). Several possible biological barriers are recognized that can influence RNAi efficacy among different organisms. Among these are gut pH, RNase activity in the gut, salivary glands, and other organs and tissues, endosomal entrapment of dsRNA, and variable extents of systemic spread of RNAi, once triggered (Parsons et al., 2018; Cooper et al., 2019; Zhu and Palli, 2020).
Biopesticides that exploit RNAi technology are of interest for management of hemipteran pests and vectors of plant viruses, particularly aphids, mealybugs, psyllids, scales, thrips, and whiteflies. The potato or tomato psyllid (PoP), Bactericera cockerelli (Sulc, 1909) (Hemiptera: Triozidae) is the insect vector of the causal agent of the zebra chip of potato and tomato vein-greening diseases, which can cause yield loss and reduced tuber and fruit quality, respectively (Insect Pests of Potato, 2013). Both diseases are caused by the psyllid-transmitted fastidious bacterium “Candidatus” (“Ca”) Liberibacter solanacearum” (CLso) (Brown et al., 2010; Munyaneza, 2012). A distant relative of PoP, the Asian citrus psyllid (ACP) Diaphorina citri Kuwayama (Hemiptera: Liviidae), is the insect vector of “Ca. Liberibacter asiaticus” (CLas), the causal agent of citrus greening disease, which is recognized as a major threat to the citrus industry worldwide (Halbert and Manjunath, 2004). The internal anatomies of PoP and ACP have been shown to be quite similar (Cicero et al., 2009; Brown, 2016) and render analogous pathophysiological responses to infection by CLso and CLas, respectively (Fisher et al., 2014). Close inspection of psyllid genes involved in the RNAi machinery has revealed similar characteristics of the gene silencing pathways of PoP and ACP, respectively (Fisher et al., 2014; Vyas et al., 2015). Thus, in addition to the development of RNAi biopesticides for management of PoP-transmitted bacterial plant pathogens, PoP has been instrumental in advancing ACP RNAi research. It has been deployed as a surrogate study system for pre-screening dsRNAs in PoP prior to evaluation in the more recalcitrant ACP-CLas pathosystem (Wuriyanghan and Falk, 2013; Fisher et al., 2014). Even so, definitive studies of optimal parameters required for RNAi bioassays in PoP have been quite limited, rendering the relevant characteristics fundamental to effective environmental RNAi in psyllids unknown. Environmental RNAi parameters and paradigms have been established for model organisms, including fruit fly (Diptera), red flour beetle (Lepidoptera), and several well-studied non-model insects. However, such parameters may or may not be applicable to evolutionarily distant relatives because of the propensity of different insects to respond variably to RNAi, and because of the fluid nature of RNAi biology itself. Thus, fundamental studies are required to establish optimal RNAi parameters for each insect species of interest. To now, RNAi in psyllids has not been systematically characterized by the evaluation of psyllid-specific dose effects, persistence of gene silencing, or with respect to the developmental stage(s) that are most RNAi-susceptible to the RNAi pathways that can be best exploited for field-ready biopesticide technologies.
Previous RNAi studies involving PoP have used different delivery methods and range of dsRNA concentrations over variable ingestion-access periods (IAP) (Wuriyanghan et al., 2011; Wuriyanghan and Falk, 2013; Galdeano et al., 2017), making comparisons between the different parameters difficult or impossible. To further RNA biopesticide development it has become essential to establish the environmental RNAi characteristics of psyllids, develop standardized approaches for dsRNA target selection, design, and screening, and determine the optimal dsRNA concentration(s) required for optimal knockdown among the different instars. A standardized delivery approach, together with dsRNA reference doses for knocking down candidate gene targets, will enable comparisons between independent laboratories and expedite translation of RNAi to field application for effective psyllid management. Among the most relevant of the unstudied or poorly-studied features of environmental RNAi in PoP and ACP are the optimal knockdown parameters for the different instars, the minimum and maximum effective dsRNA concentration(s) for optimal gene knockdown, and defined points at which maximum silencing can be achieved, or at which silencing reaches a plateau, and at which gene expression returns to baseline after the initial dsRNA-ingestion.
The objective of this study was to determine the optimal parameters for environmental RNAi in psyllids using the PoP-CLso study system and translate the results downstream to the difficult-to-manipulate ACP-CLas study system. In preliminary experiments, ten PoP genes were selected as potential candidates for gene knockdown, identified in a literature survey of model and non-model arthropod targets previously validated by reverse genetics or RNAi knockdown (Royle, 2006; Forgac, 2007; Weiss et al., 2009; Pinheiro et al., 2018), and therefore were considered likely to result in knockdown in psyllids, including PoP. The laboratory PoP transcriptome databases were mined to confirm the presence/expression of orthologous sequences in libraries constructed from adult and immature PoP instars (Fisher et al., 2014). In preliminary experiments, the potential for gene knockdown in psyllids was demonstrated for approximately ten candidate genes, post-dsRNA ingestion by PoP adults and third instar nymphs (0%–40%), adopting the approach that resulted in significant knockdown in adult whiteflies (Vyas et al., 2017).
Based on preliminary results, vacuolar ATPase A subunit (vATPase-A), clathrin heavy chain (CHC), and non-fermenting protein 7 (Snf7) were selected for the study. Gene knockdown was evaluated by quantitative RT-PCR amplification of the respective transcript in PoP adults, post-oral delivery of seven different concentrations of dsRNA, ranging from 0.1 ng/μL to 500 ng/μL. Second, the extent of RNAi persistence in adult PoP was determined for the minimal effective concentration resulting in optimal gene silencing, determined as 100 ng/μL. Finally, to evaluate the possible effects of psyllid developmental stage on PoP RNAi responsiveness (Araujo et al., 2006; Terenius et al., 2011; Ulrich et al., 2015; Vogel et al., 2019; McCusker et al., 2020), RNAi penetrance was evaluated for PoP third instar nymphs, immature or “teneral” adults, and mature or “reproductive stage” adults.
Materials and methods
Psyllid colony
The “Ca’ Liberibacter solanacearum” (CLso)-infected potato psyllid colony was established in 2004 from adult psyllids collected from greenhouse tomato plants in southern Arizona (Brown et al., 2010). The colony has been maintained by serial transfer to 10–12 leaf stage “Roma” tomato plants at approximately every 4–5 wk. Potato psyllids were identified as the “central haplotype” based on genotype-specific SNPs in the cytochrome oxidase I gene (Fisher et al., 2014). Positive infection of the PoP colony by CLso was demonstrated by PCR amplification of the 16S rRNA gene with CLso-specific primers (Liefting et al., 2009), and confirmatory sequencing of amplicons. Ten individual adult psyllids were collected from each colony (cage), year-round for 12 months and assayed for CLso presence to determine percent infectivity. Routinely, the infection rate for the PoP sub-colonies was 90%–100%. The three developmental stages collected from CLso-positive PoP were identified as the 3rd instar (immature), teneral adult, and mature adult by viewing under a dissecting microscope.
Gene selection for evaluation in the knockdown bioassay
Based on preliminary results (data not shown), the PoP study system was used to evaluate dsRNA-mediated knockdown of PoP vacuolar ATPase A subunit (vATPase-A), clathrin heavy chain (CHC), and non-fermenting protein 7 (Snf7). These targets were selected in part because they were among those exhibiting the highest levels of gene silencing compared to others evaluated in our laboratory. For example, the wings apart-like protein homolog (WAPL), adaptor protein (AP), and Wnt-1 were among the candidates screened that exhibited extremely low to no knockdown (data not shown). Second, the aim was to avoid gene targets for which relative knockdown was less than 40%, to maximize the likelihood that results were reflective of psyllid responsiveness to each increase in dsRNA concentration in the series, for example, 1 ng/μl to 10 ng/μl, and yield a measurable knockdown phenotype for each. Further, the three genes selected are well-known RNAi targets in insects (Baum et al., 2007; Baum and Roberts, 2014; Ulrich et al., 2015; Pereira et al., 2019; Angelotti-Mendonça et al., 2020), and are expressed constitutively in most organisms including insects (Royle, 2006; Forgac, 2007; Weiss et al., 2009; Pinheiro et al., 2018). Finally, expression of each gene has been confirmed by scanning our laboratory transcriptome libraries constructed from PoP adult and nymphal stages, and in the PoP gut and salivary glands of both CLso-infected and -uninfected PoP (data not shown).
Double-stranded RNA region selection, and in-vitro transcription
From the coding region of each gene target, a 200–250 base pair (bp) region was selected as the template for dsRNA synthesis using an in vitro transcription reaction. The dsRNA regions were selected bioinformatically using a custom pipeline to ensure that siRNAs produced from the dsRNA target only the gene of interest, e.g., vATPase-A. Briefly, from an initially selected 200–250 bp coding sequence, 21-mers were obtained from the sense and anti-sense strand. Using Bowtie (Langmead et al., 2009), the 21-mers were mapped to the laboratory databases of PoP transcriptome libraries constructed from whole body adults and nymphs and the alimentary canals (guts) and salivary glands (Fisher et al., 2014), allowing one mismatch per read. The number of mapping events against each transcript was counted using BEDtools (Quinlan and Hall, 2010). When a non-target gene was identified by the mapping, a different region of 200–250 bp was selected and the pipeline was re-run. Finally, a region having the minimum predictable off-target matches was selected for each target gene and cloned into the pGEM®-T Easy vector (Promega, Cat. No. A3600) (Data Sheet S2; Supplementary Tables S2–S5). From the psiCHECK-2 vector (Promega, Cat. No. C8021), a 231-nucleotide (nt) region was selected from the firefly luciferase gene sequence available for the firefly Photinus pyralis (L.) (Genbank Accession no. AY535007) following pipeline screening as described above, and the PCR amplicons were cloned into the pGEM®-T Easy plasmid vector. The sequence for each amplicon was verified by bi-directional DNA (Sanger) sequencing. The dsRNA was synthesized from the cloned fragment of each target through an in vitro transcription reaction using the MEGAscript T7 Transcription Kit (Thermo Fisher Scientific, Cat. No. AM1334), according to the manufacturer’s protocol.
Potato psyllid ingestion-access to double-stranded RNA
The dsRNAs were delivered to psyllids in a sterile 20% sucrose (w/v) solution in an artificial feeding chamber by pipetting the dsRNA-sucrose solution between two layers of thinly stretched Parafilm M® (Denville Scientific Inc. USA) as previously described in studies of the whitefly (Vyas et al., 2017) and potato psyllid (Paredes-Montero et al., 2022) carried out in this laboratory. The Parafilm M® layers were then attached to a 50 ml polypropylene conical centrifuge tube (Fisher Scientific, catalog # 14-432-22). Ten mL of sterile sucrose solution was prepared for each ingestion-access experiment. The sucrose solution contained one drop (∼50 μL) of green food-coloring (Wilton® Icing Color) used to monitor psyllid ingestion, which was indicated by the green-colored honeydew. Experiments were carried to completion only for those in which psyllids were observed to excrete green honeydew. Three to five independent biological replicates were carried out for each treatment/experiment.
For adult experiments, each replicate consisted of a cohort of ten PoP teneral adults, given a 48-h IAP on 200 μL of sucrose-dsRNA solution. Adult psyllids were collected from an infested, CLso-infected tomato plant using a hand-held aspirator, and transferred to a petri dish containing filter paper for a 4-h starvation period. Psyllids were transferred to the artificial feeding chamber and allowed a 2-day IAP on the dsRNA-containing sucrose solution, after which the dsRNA-sucrose solution was replaced with 20% sucrose. Psyllids were held in the chamber for the duration of the experiment in which dose and RNAi persistence were assessed. Integrity of the dsRNAs were evaluated by agarose gel electrophoresis, post-synthesis, pre-IAP to confirm they were of the expected size, and similarly, post-IAP by electrophoresis of 10 μL volume of the sucrose solution. Only experiments for which dsRNAs were found to be intact on the gel (undegraded), were considered in the analysis reported here. The dsRNA was not replenished or replaced with freshly prepared dsRNA in the “long feeding period” evaluated for RNAi persistence assays.
For immature psyllid experiments, a cohort of fifty 3rd instar nymphs were used in each biological replicate. Third instar nymphs were collected from CLso-infected tomato plants and placed into a (covered) petri dish using a fine paintbrush. Nymphs were starved for 6 h to reduce their vigor and encourage subsequent immediate feeding to discourage their roaming away from the sucrose solution and dying. All live 3rd instar psyllids (at least 30 per replicate) were transferred to 3–4 leaf stage “Roma” tomato seedlings, post-IAP.
Total RNA isolation, cDNA synthesis, and quantitative RT-PCR
Tubes containing psyllids were submerged in liquid nitrogen for 5 min and transferred to and stored at −80°C. Psyllid cohorts were homogenized in 1 ml Tri Reagent (Zymo Research Cat. No. R2050-1-200) containing 0.5 mm Zirconium beads in a bead-beater (BioSpec Products) (RPI Research Products). The adult and immature psyllids were homogenized in the bead-beater for 8 and 6 min, respectively. Total RNA was isolated using a widely-used Trizol RNA isolation protocol (Rio et al., 2010). Residual DNA was removed from total RNA using the DNA-free™ DNA Removal Kit (Invitrogen, Lot No. 00522653). Using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Lot No. 00692533), cDNA was synthesized from 2 μg total RNA per sample, according to the manufacturer’s instructions. The qRT-PCR amplification was carried out using a CFX96TM Real-Time PCR System (Bio-Rad Laboratories) and TaqMan RT-PCR master mix (Applied Biosystems; Universal PCR Master Mix, Lot No. 1908161). The psyllid ribosomal protein L5 (RPL5) gene was used as the internal reference gene to normalize candidate gene expression with the ΔΔCt method (Lü et al., 2018). The efficiency of each of the primers/probe combinations was determined based on a 10-fold dilution standard curve (0–10–6). Only reaction efficiencies higher than 90% were considered acceptable for qRT-PCR amplification carried out to quantify differential gene expression. The CFX Maestro software v1.1 (Bio-Rad Laboratories) was used to carry out statistical analysis based on the student’s t-test to determine the differences between the means, and Tukey’s HSD was used to identify means that were significantly different. The primers and probe sequences and gene target sequences are provided in the Data Sheet S1.
Results
Minimal effective concentration of double-stranded RNA for optimal gene knockdown
The aim of the study was to evaluate the effectiveness of gene knockdown in the potato psyllid using a range of concentrations of orally-administered synthetic dsRNA using a 48-h post ingestion-access period (IAP). Seven dsRNA concentrations of 0.1, 1, 10, 100, 150, 200, and 500 ng/μL were evaluated for vATPase-A and CHC, and six concentrations were evaluated for Snf7 (all concentrations mentioned above except the 150 ng/μL concentration). The concentrations of dsRNA were selected based on the range of previously reported doses of dsRNA capable of eliciting an RNAi response in PoP or ACP of 0.3 ng/μL to 1,000 ng/μL (Wuriyanghan et al., 2011; El-Shesheny et al., 2013; Wuriyanghan and Falk, 2013; Andrade and Hunter, 2017; Galdeano et al., 2017; Ghosh et al., 2017; Kishk et al., 2017; Yu et al., 2022). Teneral adults were allowed a 2-day post-IAP on 200 μL sterile 20% sucrose solution containing synthetic dsRNA, sandwiched between two layers of parafilm stretched across an artificial feeding chamber, as described previously (Vyas et al., 2017; Paredes-Montero et al., 2022). Psyllids have been reported to show lower sensitivity to RNAi compared to several insects in certain other insect species or orders (Wuriyanghan and Falk, 2013). However, information regarding potential saturation of the RNA-induced silencing complex (RISC) is lacking. Based on the results of a preliminary study conducted in our laboratory (data not shown), a 48-h IAP was selected for the experiments reported here. Gene knockdown was quantified for surviving PoP adults 4 days post-IAP by RT-qPCR amplification. Gene expression was normalized to gene knockdown levels determined by qPCR amplification for the non-target luciferase dsRNA control at 200 ng/μL (Figure 1).
FIGURE 1. Relative normalized expression of potato psyllid (PoP) genes following ingestion of dsRNAs of the different concentrations (A–C). Relative normalized expression of vATPase-A subunit (A), Snf7 (B), and CHC (C) genes in “Candidatus Liberibacter solanacearum”-infected PoP teneral adults. A 2-day ingestion-access period (IAP) to dsRNAs (100 ng/μL, 200 μL) was provided to psyllids and real time quantitative reverse transcriptase-polymerase chain reaction amplification was carried out 6 days post-IAP. Expression was normalized to the psyllid cohorts receiving 200 ng/μL of the dsRNA luciferase dsRNA non-target control. The Student’s t-test was used to carry out the statistical analysis. Significant groups (p-value <0.0.5) are indicated by an asterisk (*) when the results for the experimental and non-target control groups differed. Lower case letters a-c indicate significantly different groups, determined by Tukey’s HSD test (p-value 0.05). Each experiment consisted of three to five independent biological replicates, with 10 PoP teneral adults per experiment.
Significant knockdown of vATPase-A gene expression was not observed for the dsRNA concentrations of 0.1 and 1 ng/μL. However, increasing the concentration to 10 ng/μL and greater yielded significant gene knockdown. Even so, at a dose of 10 ng/μL only 20% knockdown was observed. Concentrations ranging from 100 ng/μL to 500 ng/μL yielded 42%–50% knockdown, with insignificant variation among the four concentrations (Tukey’s HSD test). Results indicated that the minimal concentration of dsRNA required for optimal knockdown of vATPase-A gene was 100 ng/μL. Increasing the concentration to 500 ng/μL dsRNA did not appreciably increase knockdown, which is potentially suggestive of RISC saturation. Although there is no direct evidence for RISC saturation resulting from treatment with high concentrations of siRNA/dsRNA from any insect studies, the phenomenon is well established in mammalian systems (Grimm et al., 2006; Tan et al., 2012; Varley and Desaulniers, 2021). The results of the knockdown experiments for Snf7 and CHC indicated a similar trend with the 100 ng/μL to 500 ng/μL dsRNA treatments resulting in respective gene silencing levels like those observed for the vATPase-A gene. These results indicated that RNAi penetrance is lower for psyllids than the whitefly B. tabaci and for most coleopteran insects (Mondal et al., 2020; Mehlhorn et al., 2021). Further, a five-fold increase in dsRNA concentration from 100 ng/μL and greater did not result in increased gene knockdown.
Persistence of RNA interference-induced gene silencing in potato psyllids
The persistence of RNAi-induced gene silencing in psyllids has not been well studied. The experiments were carried out to evaluate the persistence of gene silencing in PoP at the pre-determined minimal effective dsRNA concentration of 100 ng/μL (see Results, above). To this end, the endogenous expression of the three selected target genes, vATPase-A (VAT-A), Snf7, and CHC, was quantified in time-course experiments. On day 0, newly emerged teneral adults were collected from PoP colonies and gene expression was quantified on days 0, 4, 7 and 10 by qRT-PCR amplification (Figure 2A). For all experiments, psyllids were held in an artificial feeding chamber for the duration of the study, as indicated. Methods were as described above for dsRNA ingestion-access bioassays. Relatively consistent gene expression levels, based on knockdown, were observed for the three genes evaluated for both the PoP teneral and mature adult stage psyllids (Figure 2A).
FIGURE 2. Persistence of RNAi mediated gene silencing in potato psyllids (A–D). (A) Endogenous expression of “Candidatus Liberibacter solanacearum”-infected PoP vATPase-A, Snf, and CHC genes at day 0, 4, 7 and 10; (B–D) Gene expression for candidate genes was normalized to the respective expression level corresponding to that on day zero (0). Time-course study of relative normalized expression of vATPase subunit-A (B), Snf7 (C), and CHC (D) genes for CLso-infected PoP teneral adults 10-day post-IAP. On day zero (0), psyllids were given an ingestion-access period (IAP) on dsRNA for each respective gene target or the luciferase non-target control (Luc) at 100 ng/μL, 200 μL. On day two (48-h), the dsRNA-sucrose solution was replaced with 20% (w/v) sucrose solution (dsRNA-free). Gene expression was determined by real-time, quantitative reverse transcriptase-polymerase chain reaction amplification, and the results were compared to gene expression levels corresponding to day 0. The asterisk (*) indicates differences in gene expression between day 0 and each respective day for which expression was analyzed, based on the Student’s t-test (p-value <0.0.5). By convention, Tukey’s HSD test was applied to groups exhibiting significant differences relative to day 0 expression. The lower case letters a-d indicate different significant groups determined by Tukey’s HSD test (p-value 0.05). Each experiment consisted of three to five independent biological replicates for each “dsRNA ingestion”- group consisting of 10 PoP teneral adults per experiment.
Next, psyllids were allowed a 48-h IAP to the dsRNA in a 20% sucrose (w/v) solution, designating day 0 as the initial or ‘first day’ of dsRNA ingestion-access. Forty-eight hr post-IAP (day 2), the dsRNA-sucrose solution was replaced with 20% (w/v) sucrose-only solution. Gene knockdown was determined by qRT-PCR amplification over the 12-day time course for PoP teneral adults using qRT-PCR. Gene expression was monitored on days 1, 2, 3, 4, 5, 6, 7, 8, 10, and 12 and compared to that for day 0 for vATPase-A and Chc, while knockdown of the Snf7 gene expression was measured up to day 10. Knockdown among the three genes was variable, as might be expected with respect to initiating significant level knockdown compared to the day 0 starting point, as was the persistence of silencing (Figure 2). Knockdown of the vATPase-A gene expression was significant after 48-h (day 2), and persisted until day10 (8 days post-IAP), with maximum knockdown documented during days 4–7 (Figure 2B). Interestingly, silencing began to decline by day 8 and by day 12 (end of experiment), knockdown of vATPase-A was not significant. In contrast, significant knockdown of Snf7 gene expression was observed on days 4–7, with no significant knockdown detectable by 6–8 days (Figure 2C). The patterns observed for CHC and Snf7 gene knockdown were similar, with significant gene knockdown occurring on days 4–6 and exhibiting a gradual decline to the day 0 levels (Figure 2D). Thus, dsRNA-induced gene silencing was diminished entirely by 4–5 days, post-IAP.
Comparison of RNAi efficiency for different psyllid developmental stages
In insects, the developmental stage has been found to influence RNAi penetrance, a phenomenon resulting from differences in the response to exogenous RNAi triggers (Schmitt-Engel et al., 2015). In general, younger insect developmental stages have been found to be more receptive, or sensitive, to knockdown, compared to more advanced life stages (instars) (Guo et al., 2021; Salvador et al., 2021; Willow and Veromann, 2021). The psyllid life cycle consists of five nymphal instar stages, the newly emerged adult, also referred to as a teneral adult, and the mature adult (Butler and Trumble, 2012). In this study, the responses to environmental RNAi were compared for PoP 3rd instar nymphs, teneral adults (light brown or yellow body color), and mature adults (dark brown body color). Psyllids were allowed an IAP on dsRNA solution, v-ATPase-A, Snf7, or CHC at 100 ng/μL, in 20% (w/v) sucrose solution for 48 h. After the 48-h IAP, psyllids were transferred to a 3–4 leaf stage “Roma” tomato seedling. A similar pattern of gene knockdown was observed for all three psyllid life stages, with the nymphal stage exhibiting the greatest sensitivity to environmental RNAi compared to the teneral and mature adult stages, indicating at least some difference in responsiveness by nymphs and adults (Figure 3). A similar trend was observed for the three candidate gene targets, except for Snf7, for which knockdown in adult psyllids was not significantly different from the non-target control, even though significant knockdown was observed for the 3rd instar nymphs and teneral adults (Figure 3B). Knockdown of vATPase-A and CHC targets for teneral adults was not significantly different from that observed for mature adults. These observations indicate that RNAi penetrance in the 3rd instar psyllid nymphs was significantly higher than for the two older (adult) life stages, indicating that among the developmental stages (instars) examined in this study, potato psyllid nymphs were the most sensitive to environmental RNAi-mediated gene knockdown, compared to the teneral and mature adults.
FIGURE 3. Comparison of gene knockdown levels across the developmental stages in potato psyllid (A–C). Relative normalized expression of vATPase-A subunit (A), Snf7 (B), and CHC (C) genes in “Candidatus Liberibacter solanacearum”-infected mature adults, teneral adults, and 3rd instar nymphs. Psyllids were provided a 2-day inoculation-access period (IAP) on dsRNA for each respective gene target or the luciferase non-target control (Luc) at 100 ng/μL, 200 μL. Gene expression was determined by real time, quantitative reverse-transcriptase polymerase chain reaction amplification, 6 days post-IAP. Relative gene expression was calculated by comparing expression for dsRNA-treated psyllids of each analogous instar cohort, post-IAP on dsRNA for the gene target and non-target Luc control. An asterisk (*) denotes the results of the Student’s t-test (p-value <0.0.5) for experimental and non-target control group comparisons. Each experiment consisted of three biological replicates per instar cohort. Eight to ten mature and teneral adults and fifteen to 20 3rd instar nymphs, respectively, were collected 4-days post-IAP for gene expression analysis.
Discussion
Three key features of environmental RNAi were characterized for the potato psyllid based on ingestion-access bioassay method in this study. Several other delivery approaches have been evaluated for dsRNA in RNAi bioassays, including injection into the abdomen (Dzitoyeva et al., 2001; Ulrich et al., 2015), topical application to external body parts or imbibition of topically-applied dsRNA (Killiny and Kishk, 2017; Vaschetto, 2022, pp. 253–277), ingestion-access to dsRNA applied to leaf disks (Ghosh et al., 2017), transgenic plant-mediated expression to facilitate phloem delivery with “continuous” exposure by over-expression of dsRNA hairpins (Eakteiman et al., 2018), stem cutting- or root- uptake of dsRNA (Ghosh et al., 2017; Hunter and Wintermantel, 2021), and/or dsRNA ingestion using an artificial feeding system (Yu et al., 2013; Ghosh et al., 2017; Zhu and Palli, 2020). Among these approaches, ingestion of dsRNA in sucrose feeding systems is convenient and affords precise control over ingestion of dsRNA that arrives directly in the gut (Vyas et al., 2017). For this study, the latter approach was implemented to define the following three key features of environmental RNAi for PoP, the minimal effective concentration of dsRNA for effective gene knockdown, persistence of the knockdown, and differential penetrance of RNAi in three psyllid developmental stages, 3rd instar, and teneral and mature adults. Gene knockdown was quantified for each concentration by real time qRT-PCR amplification. Based on the minimal effective concentration established, the persistence of gene silencing was evaluated over a 12-day period.
The PoP targets, vATPase-A, Snf7, and CHC, exhibited relatively high knockdown at 40%–60%. Because these genes are involved in insect cellular processes their functions are expected to be critical in most tissues and organs during all or most insect developmental stages. The vATPases are involved in metabolic signaling, oxidative phosphorylation, phagosome, endosome, and lysosome formation, the synaptic vesicle cycle and in other pathways that require pH regulation. The Snf7 and CHC genes are involved in endocytosis, and endocytosis and lysosome pathways, respectively. All three genes are conserved across insect orders and have been evaluated in RNAi research for multiple insect species (Baum et al., 2007; Angelotti-Mendonça et al., 2020), making them amenable for comparative analysis of RNAi characteristics with other better-studied insects. Further, optimal RNAi efficacy for biopesticide development rests on the consistent expression of the genes across all insect developmental stages.
For the three candidate genes studied here, the minimal effective dose of the dsRNA was 100 ng/μL (100 parts per million, ppm), which is at least three orders of magnitude higher than the optimal concentration reported for coleopteran environmental RNAi (Baum and Roberts, 2014). These observations underscore the challenges encountered in developing RNAi as a biopesticide for controlling psyllids, and likely for other hemipteran pests. Psyllids have piercing-sucking mouthparts, and they specialize in phloem-feeding. Delivering sufficient, accessible dsRNA into the plant phloem has been challenging for methods such as foliar spray, trunk injection, and laser-etching. Recently, foliar application of dsRNA was shown to be successful when delivered in a matrix of layered double hydroxide (LDH), referred to as BioClay, for whitefly control (Jain et al., 2022). Although this material has not been evaluated for dsRNA delivery in psyllids, the RNAi knockdown in whitefly has been shown to be superior to RNAi effectiveness in psyllids (Mondal et al., 2020; Yyas et al., 2017). Thus, alternative approaches will be required to deliver concentrations of dsRNA sufficient for achieving gene silencing and persistence in the long term, especially for woody plants, such as Citrus spp.
The persistence of RNAi-induced gene silencing, a crucial characteristic of RNAi, was investigated, with a mindfulness to applications in the field. This key RNAi feature is not well-studied for psyllids or most other insect or arthropod pests of interest. In general, persistence data are scarce in the literature, the community has not established or implemented consistent methodologies and/or concentration regimes, and at the same time diverse parameters and delivery systems have been adopted across other insect study systems, making comparisons between experimental results difficult. Ingestion-access durations ranging from 2 to 8 days have been used in studies of other hemipterans without knowledge of whether exposure/ingestion times could be too short or too lengthy. Consequently, the interpretation of the results depend on many factors, including accurate functional annotations of gene targets, relevance of the gene to the insect life stage, tissue-specific expression of the gene, and optimal concentration of the dsRNA (Andrade and Hunter, 2017; Galdeano et al., 2017; Ghosh et al., 2018; Kishk et al., 2017; Tanning et al., 2016a; Wuriyanghan and Falk 2013; Wuriyanghan, Rosa, and Falk 2011; Yu and Killiny 2020; Yu et al., 2017). To address at least several of the latter caveats, the persistence of gene silencing at the pre-determined, minimal effective dsRNA concentration of 100 ng/μL (see Results, above) was determined. The RNAi persistence results suggest that in PoP, and potentially in other psyllids, the duration of the dsRNA-induced gene silencing is short-term and persists 4–5 days following removal of the dsRNA source. This indicates that the RNAi triggering molecule, i.e., dsRNA, must be available continuously to perpetuate the silencing effect, which is expected to be essential for achieving a discernible phenotype, such as psyllid mortality or reduced fecundity, among others. Therefore, delivery of dsRNA by viral-induced gene silencing vectors (VIGS vectors) expressing dsRNA hairpins, or transgenic plant-mediated dsRNA expression either from the chloroplast or nuclear genome approaches may represent the most viable options available for dsRNA delivery for psyllid pest and vector control (Hajeri et al., 2014; Dong et al., 2022; He, 2022; Wu et al., 2022).
The low RNAi persistence documented in this study was not unexpected for psyllids. Among invertebrates, RNAi persistent has been shown to be robust in C. elegans, and is due primarily to the initially strong systemic RNAi followed by the amplifying effect of the RNAi trigger by the RNA dependent RNA polymerase (RdRP), which further has been shown to propagate gene silencing through several generations (Spracklin et al., 2017). Insects do not encode RdRP and systemic RNAi may not operate in psyllids and/or in other hemipterans. In the pea aphid, the delivery of dsRNA by expression in transgenic plants resulted in knockdown that reached a maximum between day 4 and day 8 (Coleman et al., 2015). Aphids provided longer ingestion exposures to plants did not exhibit increased knockdown, indicating that the RNAi effect reached a plateau after day 8, potentially due to saturation of capacity of RISC. Further, 6 days after the dsRNA source was removed, gene expression levels fell to baseline levels. In another study, expression of pea aphid vATPase subunit E (VAT-E) was silenced by injection-delivered VAT-E dsRNA, and the duration of silencing was diminished after 72 h. This transient RNAi effect was attributed to a rapid degradation of the dsRNA by exo- and endo- nucleases in the pea aphid (Cao et al., 2018). This phenomenon was also observed for psyllids, which also express transcripts of several endonucleases in the gut and salivary glands (data not shown). Also, because the pH of the psyllid gut is basic, dsRNA may not withstand the environment for long periods of time, sufficient to initiate and sustain silencing. The latter barriers potentially limit delivery of any significant amount of dsRNA beyond the gut, resulting in a transient RNAi phenotype. Interestingly, silencing of the vATPase-A gene was initiated sooner and persisted longer compared to CHC and Snf7 (Figure 2). Further, a lower dsRNA concentration (10 ng/μL) resulted in significant knockdown of vATPase-A expression even though only 20% silencing was achieved at this dose (Figure 1). Together, the results indicate that the dsRNA concentration required to achieve an analogous level of gene silencing is lower for vATPase-A than for CHC and Snf7, or that other biological factors such as systemic properties of RNAi could be responsible for the differences in persistence of their silencing. The higher endogenous expression of the vATPase-A gene in psyllid gut compared to the other two genes evaluated (Supplementary Figure S1) might influence knockdown persistence. One possible explanation could be that RNAi more effectively silences genes in the gut epithelium than those expressed in other tissues. Generally, targeting gut-expressed genes results in higher RNAi efficiency in organisms when systemic RNAi is not robust (Vyas et al., 2017). Previously Wuriyanghan et al. (2011), Wuriyanghan and Falk (2013) reported significant knockdown of an ATPase gene in the gut, compared to that for whole body or abdomen samples (Wuriyanghan et al., 2011; Wuriyanghan and Falk, 2013). Although the robustness of systemic RNAi in psyllids is unknown, the greater persistence of knockdown of the vATPase-A gene, compared to results for Snf7 and CHC, has demonstrated that gut-localized genes are effectively knocked down by RNAi in potato psyllid.
In C. elegans, systemic RNA interference defective 2 (Sid2) protein is responsible for passive transport of the dsRNAs from the intestinal lumen to the epithelial cells, whereas the Sid1 protein carries out passive bidirectional transport of the dsRNA between cells (Tomoyasu et al., 2008; Zhu and Palli, 2020). Though Sid2 orthologs are absent in all insects, there are multiple Sid1-like proteins present in insects including psyllids (Taning et al., 2016b). The red flour beetle Tribolium castaneum (Herbst, 1797) has systemic RNAi properties, but Sid1 genes are not required for cell-to-cell transport of the dsRNAs in T. castaneum; thus, among different insects the function of these genes are not completely understood (Tomoyasu et al., 2008). The cellular uptake and systemic properties of RNAi appear to be complex and versatile in insects, and there are potentially multiple pathways involved in this process. Thus, the biological mechanisms underlying differences observed in silencing characteristics of PoP vATPase-A gene and those of Snf7 and CHC cannot yet be explained.
Finally, the 3rd instar (immature) psyllids were more responsive to environmental RNAi than their counterpart teneral or mature adults. Although RNAi penetrance among the different instars of the same psyllid species, or between species, has not been systematically evaluated, Wuriyanghan and Falk (2013) reported that immature psyllids ingested and accumulated higher levels of viral-derived RNA than PoP teneral adults do (Wuriyanghan and Falk, 2013). Other studies have reported similar trends for the different insect developmental stages. For example, in the nymphal stage of the Sunn pest Eurygaster integriceps Puton (Hemiptera: Scutelleridae), gene silencing levels were greater in the nymph than for adults (Amiri and Bandani, 2020). The spotted wing fruit fly Drosophila suzukii (Matsumuru, 1937) larval stage also exhibited higher levels of knockdown across multiple genes, where gene silencing was documented at 32%–42% and 19%–24%, respectively, in larval and adult stages, and mortality was 10%–23% in adults, compared to 22%–42% in larvae (Taning et al., 2016b). By comparison, in two coleopteran pests, the southern and western corn rootworms, v-ATPase-A and Snf7 genes showed higher levels of gene silencing in larvae compared to the adults, at 55%–89% and 48%–83% knockdown, respectively. Similarly, gene silencing experiments with the 2nd and 4th early larval stages of the hemipteran R. prolixus resulted in a frequency of 42% for the 2nd instar and little or no detection in the 4th instar (Araujo et al., 2006).
The distinctive RNAi responses observed across the potato psyllid life stages are likely attributable to several factors. The overall expression of the three gene targets evaluated here did not differ dramatically between the adult and nymphal stages. However, the feeding habits of the psyllid stages are different, with immatures being more voracious feeders than the adults, which is potentially attributed to physiology of younger instars that experience multiple cycles of rapid growing and molting (George et al., 2017; George et al., 2018). Consequently, the immature stages must be more metabolically active, a scenario that might lead to more rapid processing of ingested dsRNAs, thereby avoiding inactivation or degradation in the basic (high pH) environment of the gut lumen. Whether the pH of the adult gut is more basic than the immatures’ is not known. Also, at least two dsRNA-specific endonucleases are encoded by the psyllids, and transcriptomic evidence indicate these are expressed in the gut and salivary glands. These enzymes are expressed to a greater extent in adults than nymphs (Supplementary Figure S2), which could result in greater or more rapid degradation of dsRNAs ingested by adults than by immature stages. Finally, the RNAi core machinery proteins could possibly be involved in the observed differential expression in psyllid adults and nymphs thereby influencing the differential processing of the exogenous dsRNAs in different or all psyllid life stages. To investigate this possibility, relative expression of psyllid RNAi machinery proteins was calculated from data reported in a previously published study, with manual curation of the datasets (Fisher et al., 2014) (Data Sheet S2, Supplementary Table S1). Interestingly, Dicer2 expression was higher in the nymphal compared to adult instars, while Ago2 expression was similar for the two life stages. These observations may possibility suggest that the nymphs harbored more efficient RNAi machinery; however, the two Dicer2 associated proteins, Loqs and R2D2, could not be identified from the PoP transcriptome, precluding the possibility of a satisfactory explanation for the differential processivity of dsRNAs observed for adults and nymphs. The different levels of RNAi efficiency among the different psyllid life stages may have important ramifications to disease management. If nymphs can be controlled effectively before reaching adulthood, the spread of “Ca. Liberibacter” spp. to the plant host could be abated to some extent by reducing the population size prior to reproduction and/or before adult dispersal and colonization of the host plants.
The candidate PoP genes evaluated in this study are conserved among members of the Class: Insecta and their primary biological roles are to carry out fundamental physiological functions. The characteristics of RNAi may vary for genes that encode proteins involved in adaptive or other non-constitutive functions such as those involved in immunity and defenses, or highly specific functions that are essential for one or more life stage. In summary, the RNAi characteristics and parameters established in this study are expected to inform experimental design for RNAi screening and selection of specific phenotypes, e.g., mating interference, mortality, reduced fecundity, and others, by providing baseline results over a range of optimal dsRNA doses to standardize laboratory studies, potentially applicable to PoP and to other psyllid species. The results will also be valuable for the development and optimization of biopesticide studies at the field level, toward effective management of psyllids as pests and vectors of psyllid-transmitted plant pathogens.
Finally, certain environmental and regulatory aspects have been considered in the design of the dsRNA molecules selected for silencing three genes evaluated in this study. The approach involves a priori comparisons of regions of the dsRNA with analogous regions of gene orthologs in genome and/or transcriptome databases for beneficial and organisms of interest in different agroecosystems, in the case of PoP, lacewing, lady beetle, springtail, honeybee, and also humans, to avoid predicted gene sequences (threshold of 0–1 identical nucleotides) with potential off-target effects as siRNAs (data not shown) (Alyokhin et al., 2022). Sequence specificity is an advantage of RNAi-biopesticides over most traditional chemical pesticides, which together with rapid degradation of the dsRNAs in the environment, are promising incentives for advancing this new type of technology and its acceptance in integrated pest and vector management programs (Christiaens et al., 2022; Vaschetto, 2022).
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.
Author contributions
MM performed the experiments, was involved in conceptualization, experimental design, data curation, and formal analysis. MC maintained the insect colony, collected the insects for all bioassays, and edited the manuscript. JB was involved in conceptualization, funding acquisition, experimental design, methodology, project administration, experimental supervision, and editing of the manuscript. MM wrote the original draft and editing of the manuscript.
Funding
The research was supported by USDA-NIFA CDRE grant number 2018-70016-27411.
Conflict of interest
Author MM was employed by RNAissance Ag LLC.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.931951/full#supplementary-material
References
Agrawal N., Dasaradhi P. V. N., Mohmmed A., Malhotra P., Bhatnagar R. K., Mukherjee S. K. (2003). RNA interference: Biology, mechanism, and applications. Microbiol. Mol. Biol. Rev. 67, 657–685. doi:10.1128/MMBR.67.4.657-685.2003
Alyokhin A., Rondon S. I., Gao Y. (2022). Insect pests of potato: Global perspectives on biology and management.
Amiri A., Bandani A. R. (2020). Gluten hydrolase gene silencing using RNAi and its effect on the Sunn pest growth and development. Phytoparasitica 48, 575–587. doi:10.1007/s12600-020-00821-8
Andrade E. C., Hunter W. B. (2017). RNAi feeding bioassay: Development of a non-transgenic approach to control asian citrus psyllid and other hemipterans. Entomol. Exp. Appl. 162, 389–396. doi:10.1111/eea.12544
Angelotti-Mendonça J., Bassan M. M., Marques J. P. R., Yamamoto P. T., Figueira A., Piedade S. M. D. S., et al. (2020). Knockdown of calreticulin, laccase , and Snf7 genes through RNAi is not effective to control the asian citrus psyllid (Hemiptera: Livideae). J. Econ. Entomol. 113, 2931–2940. doi:10.1093/jee/toaa240
Araujo R. N., Santos A., Pinto F. S., Gontijo N. F., Lehane M. J., Pereira M. H. (2006). RNA interference of the salivary gland nitrophorin 2 in the triatomine bug Rhodnius prolixus (Hemiptera: Reduviidae) by dsRNA ingestion or injection. Insect biochem. Mol. Biol. 36, 683–693. doi:10.1016/j.ibmb.2006.05.012
Baum J. A., Bogaert T., Clinton W., Heck G. R., Feldmann P., Ilagan O., et al. (2007). Control of coleopteran insect pests through RNA interference. Nat. Biotechnol. 25, 1322–1326. doi:10.1038/nbt1359
Baum J. A., Roberts J. K. (2014). “Progress towards RNAi-mediated insect pest management,” in Advances in insect physiology (Elsevier), 249–295. doi:10.1016/B978-0-12-800197-4.00005-1
J. K. Brown (Editor) (2016). Vector-mediated transmission of plant pathogens (The American Phytopathological Society). doi:10.1094/9780890545355
Brown J. K., Rehman M., Rogan D., Martin R. R., Idris A. M. (2010). First report of “ candidatus liberibacter psyllaurous” (synonym “ Ca . L. Solanacearum”) associated with ‘tomato vein-greening’ and ‘tomato psyllid yellows’ diseases in commercial greenhouses in Arizona. Plant Dis. 94, 376. doi:10.1094/PDIS-94-3-0376B
Butler C. D., Trumble J. T. (2012). The potato psyllid, Bactericera cockerelli (sulc) (Hemiptera: Triozidae): Life history, relationship to plant diseases, and management strategies. Terr. Arthropod Rev. 5, 87–111. doi:10.1163/187498312X634266
Cao M., Gatehouse J., Fitches E. (2018). A systematic study of RNAi effects and dsRNA stability in Tribolium castaneum and Acyrthosiphon pisum, following injection and ingestion of analogous dsRNAs. Int. J. Mol. Sci. 19, 1079. doi:10.3390/ijms19041079
Christiaens O., Sweet J., Dzhambazova T., Urru I., Smagghe G., Kostov K., et al. (2022). Implementation of RNAi-based arthropod pest control: Environmental risks, potential for resistance and regulatory considerations. J. Pest Sci. 95, 1–15. doi:10.1007/s10340-021-01439-3
Cicero J. M., Brown J. K., Roberts P. D., Stansly P. A. (2009). The digestive system of Diaphorina citri and Bactericera cockerelli (Hemiptera: Psyllidae). Ann. Entomol. Soc. Am. 102, 650–665. doi:10.1603/008.102.0410
Coleman A. D., Wouters R. H. M., Mugford S. T., Hogenhout S. A. (2015). Persistence and transgenerational effect of plant-mediated RNAi in aphids. J. Exp. Bot. 66, 541–548. doi:10.1093/jxb/eru450
Cooper A. M., Silver K., Zhang J., Park Y., Zhu K. Y. (2019). Molecular mechanisms influencing efficiency of RNA interference in insects: Molecular mechanisms influencing the efficiency of RNAi in insects. Pest Manag. Sci. 75, 18–28. doi:10.1002/ps.5126
Dong Y., Wu M., Zhang Q., Fu J., Loiacono F. V., Yang Y., et al. (2022). Control of a sap-sucking insect pest by plastid-mediated RNA interference. Mol. Plant 15, 1176–1191. doi:10.1016/j.molp.2022.05.008
Dowling D., Pauli T., Donath A., Meusemann K., Podsiadlowski L., Petersen M., et al. (2017). Phylogenetic origin and diversification of RNAi pathway genes in insects. Genome Biol. Evol. 8, 3784–3793. doi:10.1093/gbe/evw281
Dzitoyeva S., Dimitrijevic N., Manev H. (2001). Intra-abdominal injection of double-stranded RNA into anesthetized adult Drosophila triggers RNA interference in the central nervous system. Mol. Psychiatry 6, 665–670. doi:10.1038/sj.mp.4000955
Eakteiman G., Moses-Koch R., Moshitzky P., Mestre-Rincon N., Vassão D. G., Luck K., et al. (2018). Targeting detoxification genes by phloem-mediated RNAi: A new approach for controlling phloem-feeding insect pests. Insect biochem. Mol. Biol. 100, 10–21. doi:10.1016/j.ibmb.2018.05.008
El-Shesheny I., Hajeri S., El-Hawary I., Gowda S., Killiny N. (2013). Silencing abnormal wing disc gene of the Asian citrus psyllid, Diaphorina citri disrupts adult wing development and increases nymph mortality. PloS One 8, e65392. doi:10.1371/journal.pone.0065392
Fisher T., Vyas M., He R., Nelson W., Cicero J., Willer M., et al. (2014). Comparison of potato and asian citrus psyllid adult and nymph transcriptomes identified vector transcripts with potential involvement in circulative, propagative liberibacter transmission. Pathogens 3, 875–907. doi:10.3390/pathogens3040875
Forgac M. (2007). Vacuolar ATPases: Rotary proton pumps in physiology and pathophysiology. Nat. Rev. Mol. Cell Biol. 8, 917–929. doi:10.1038/nrm2272
Galdeano D. M., Breton M. C., Lopes J. R. S., Falk B. W., Machado M. A. (2017). Oral delivery of double-stranded RNAs induces mortality in nymphs and adults of the Asian citrus psyllid, Diaphorina citri. PloS One 12, e0171847. doi:10.1371/journal.pone.0171847
George J., Ammar E.-D., Hall D. G., Lapointe S. L. (2017). Sclerenchymatous ring as a barrier to phloem feeding by Asian citrus psyllid: Evidence from electrical penetration graph and visualization of stylet pathways. PLOS ONE 12, e0173520. doi:10.1371/journal.pone.0173520
George J., Ammar E.-D., Hall D. G., Shatters R. G., Lapointe S. L. (2018). Prolonged phloem ingestion by Diaphorina citri nymphs compared to adults is correlated with increased acquisition of citrus greening pathogen. Sci. Rep. 8, 10352. doi:10.1038/s41598-018-28442-6
Ghosh S. K. B., Hunter W. B., Park A. L., Gundersen-Rindal D. E. (2017). Double strand RNA delivery system for plant-sap-feeding insects. PLOS ONE 12, e0171861. doi:10.1371/journal.pone.0171861
Ghosh S. K. B., Hunter W. B., Park A. L., GGundersen-Rindal D. E. (2018). Double-stranded RNA oral delivery methods to induce RNA interference in phloem and plant-sap-feeding Hemipteran insects. J. Vis. Exp. JoVE.. doi:10.3791/57390
Grimm D., Streetz K. L., Jopling C. L., Storm T. A., Pandey K., Davis C. R., et al. (2006). Fatality in mice due to oversaturation of cellular microRNA/short hairpin RNA pathways. Nature 441, 537–541. doi:10.1038/nature04791
Guo W., Guo M., Yang C., Liu Z., Chen S., Lü J., et al. (2021). RNA interference‐mediated silencing of vATPase subunits A and E affect survival and development of the 28‐spotted ladybeetle, Henosepilachna vigintioctopunctata. Insect Sci. 28, 1664–1676. doi:10.1111/1744-7917.12899
Hajeri S., Killiny N., El-Mohtar C., Dawson W. O., Gowda S. (2014). Citrus tristeza virus-based RNAi in citrus plants induces gene silencing in Diaphorina citri, a phloem-sap sucking insect vector of citrus greening disease (Huanglongbing). J. Biotechnol. 176, 42–49. doi:10.1016/j.jbiotec.2014.02.010
Halbert S. E., Manjunath K. L. (2004). Asian citrus psyllids (sternorrhyncha: PSYLLIDAE) and greening disease of citrus: A literature review and assessment of risk in Florida. Fla. Entomol. 87, 330–353. doi:10.1653/0015-4040(2004)087[0330:ACPSPA]2.0.CO;2
He G. (2022). Engineering chloroplasts for insect pest control. Proc. Natl. Acad. Sci. U. S. A. 119, e2205125119. doi:10.1073/pnas.2205125119
Head G. P., Carroll M. W., Evans S. P., Rule D. M., Willse A. R., Clark T. L., et al. (2017). Evaluation of SmartStax and SmartStax PRO maize against Western corn rootworm and northern corn rootworm: Efficacy and resistance management. Pest Manag. Sci. 73, 1883–1899. doi:10.1002/ps.4554
Hunter W. B., Wintermantel W. M. (2021). Optimizing efficient RNAi-mediated control of Hemipteran pests (psyllids, leafhoppers, whitefly): Modified pyrimidines in dsRNA triggers. Plants 10, 1782. doi:10.3390/plants10091782
Ivashuta S., Zhang Y., Wiggins B. E., Ramaseshadri P., Segers G. C., Johnson S., et al. (2015). Environmental RNAi in herbivorous insects. RNA 21, 840–850. doi:10.1261/rna.048116.114
Jain R. G., Fletcher S. J., Manzie N., Robinson K. E., Li P., Lu E., et al. (2022). Foliar application of clay-delivered RNA interference for whitefly control. Nat. Plants 8, 535–548. doi:10.1038/s41477-022-01152-8
Killiny N., Kishk A. (2017). Delivery of dsRNA through topical feeding for RNA interference in the citrus sap piercing-sucking hemipteran, Diaphorina citri. Arch. Insect Biochem. Physiol. 95, 1394. doi:10.1002/arch.21394
Kishk A., Anber H. A. I., AbdEl‐Raof T. K., El‐Sherbeni A.-H. D., Hamed S., Gowda S., et al. (2017). RNA interference of carboxyesterases causes nymph mortality in the Asian citrus psyllid, Diaphorina citri. Arch. Insect Biochem. Physiol. 94, e21377. doi:10.1002/arch.21377
Langmead B., Trapnell C., Pop M., Salzberg S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25. doi:10.1186/gb-2009-10-3-r25
Liefting L. W., Sutherland P. W., Ward L. I., Paice K. L., Weir B. S., Clover G. R. G. (2009). A new ‘ candidatus liberibacter’ species associated with diseases of solanaceous crops. Plant Dis. 93, 208–214. doi:10.1094/PDIS-93-3-0208
Lü J., Yang C., Zhang Y., Pan H. (2018). Selection of reference genes for the normalization of rt-qPCR data in gene expression studies in insects: A systematic review. Front. Physiol. 9, 1560. doi:10.3389/fphys.2018.01560
McCusker P., Hussain W., McVeigh P., McCammick E., Clarke N. G., Robb E., et al. (2020). RNA interference dynamics in juvenile Fasciola hepatica are altered during in vitro growth and development. Int. J. Parasitol. Drugs Drug Resist. 14, 46–55. doi:10.1016/j.ijpddr.2020.08.004
Mehlhorn S., Hunnekuhl V. S., Geibel S., Nauen R., Bucher G. (2021). Establishing RNAi for basic research and pest control and identification of the most efficient target genes for pest control: A brief guide. Front. Zool. 18, 60. doi:10.1186/s12983-021-00444-7
Mondal M., Brown J. K., Flynt A. (2020). Exploiting somatic piRNAs in Bemisia tabac i enables novel gene silencing through RNA feeding. Life Sci. Alliance 3, e202000731. doi:10.26508/lsa.202000731
Mondal M., Peter J., Scarbrough O., Flynt A. (2021). Environmental RNAi pathways in the two-spotted spider mite. BMC Genomics 22, 42. doi:10.1186/s12864-020-07322-2
Moriyama M., Hosokawa T., Tanahashi M., Nikoh N., Fukatsu T. (2016). Suppression of bedbug’s reproduction by RNA interference of vitellogenin. PLOS ONE 11, e0153984. doi:10.1371/journal.pone.0153984
Munyaneza J. E. (2012). Zebra chip disease of potato: Biology, epidemiology, and management. Am. J. Pot. Res. 89, 329–350. doi:10.1007/s12230-012-9262-3
Paredes‐Montero J. R., Arif U., Brown J. K. (2022). Knockdown of ecdysteroid synthesis genes results in impaired molting and high mortality in Bactericera cockerelli (Hemiptera: Triozidae). Pest Manag. Sci. 6848, 2204–2214. doi:10.1002/ps.6848
Parsons K. H., Mondal M. H., McCormick C. L., Flynt A. S. (2018). Guanidinium-functionalized interpolyelectrolyte complexes enabling RNAi in resistant insect pests. Biomacromolecules 19, 1111–1117. doi:10.1021/acs.biomac.7b01717
Pereira A. E., Tenhumberg B., Meinke L. J., Siegfried B. D. (2019). Southern corn rootworm (Coleoptera: Chrysomelidae) adult emergence and population growth assessment after selection with vacuolar ATPase-A double-stranded RNA over multiple generations. J. Econ. Entomol. 112, 1354–1364. doi:10.1093/jee/toz008
Pinheiro D. H., Vélez A. M., Fishilevich E., Wang H., Carneiro N. P., Valencia-Jiménez A., et al. (2018). Clathrin-dependent endocytosis is associated with RNAi response in the Western corn rootworm, Diabrotica virgifera virgifera LeConte. PLOS ONE 13, e0201849. doi:10.1371/journal.pone.0201849
Quinlan A. R., Hall I. M. (2010). BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. doi:10.1093/bioinformatics/btq033
Rio D. C., Ares M., Hannon G. J., Nilsen T. W. (2010). Purification of RNA using TRIzol (TRI reagent). Cold Spring Harb. Protoc. 2010, pdb.prot5439. doi:10.1101/pdb.prot5439
Rodrigues T. B., Mishra S. K., Sridharan K., Barnes E. R., Alyokhin A., Tuttle R., et al. (2021). First sprayable double-stranded RNA-based biopesticide product targets proteasome subunit beta type-5 in Colorado potato beetle (Leptinotarsa decemlineata). Front. Plant Sci. 12, 728652. doi:10.3389/fpls.2021.728652
Royle S. J. (2006). The cellular functions of clathrin. Cell. Mol. Life Sci. 63, 1823–1832. doi:10.1007/s00018-005-5587-0
Salvador R., Niz J. M., Nakaya P. A., Pedarros A., Hopp H. E. (2021). Midgut genes knockdown by oral dsRNA administration produces a lethal effect on cotton boll weevil. Neotrop. Entomol. 50, 121–128. doi:10.1007/s13744-020-00819-1
Schmitt-Engel C., Schultheis D., Schwirz J., Ströhlein N., Troelenberg N., Majumdar U., et al. (2015). The iBeetle large-scale RNAi screen reveals gene functions for insect development and physiology. Nat. Commun. 6, 7822. doi:10.1038/ncomms8822
Shabalina S., Koonin E. (2008). Origins and evolution of eukaryotic RNA interference. Trends Ecol. Evol. 23, 578–587. doi:10.1016/j.tree.2008.06.005
Spracklin G., Fields B., Wan G., Becker D., Wallig A., Shukla A., et al. (2017). The RNAi inheritance machinery of Caenorhabditis elegans. Genetics 206, 1403–1416. doi:10.1534/genetics.116.198812
Tan G. S., Chiu C.-H., Garchow B. G., Metzler D., Diamond S. L., Kiriakidou M. (2012). Small molecule inhibition of RISC loading. ACS Chem. Biol. 7, 403–410. doi:10.1021/cb200253h
Taning C. N. T., Andrade E. C., Hunter W. B., Christiaens O., Smagghe G. (2016a). Asian citrus psyllid RNAi pathway – RNAi evidence. Sci. Rep. 6, 38082. doi:10.1038/srep38082
Taning C. N. T., Christiaens O., Berkvens N., Casteels H., Maes M., Smagghe G. (2016b). Oral RNAi to control Drosophila suzukii: Laboratory testing against larval and adult stages. J. Pest Sci. 89, 803–814. doi:10.1007/s10340-016-0736-9
Terenius O., Papanicolaou A., Garbutt J. S., Eleftherianos I., Huvenne H., Kanginakudru S., et al. (2011). RNA interference in Lepidoptera: An overview of successful and unsuccessful studies and implications for experimental design. J. Insect Physiol. 57, 231–245. doi:10.1016/j.jinsphys.2010.11.006
Tomoyasu Y., Miller S. C., Tomita S., Schoppmeier M., Grossmann D., Bucher G. (2008). Exploring systemic RNA interference in insects: A genome-wide survey for RNAi genes in Tribolium. Genome Biol. 9, R10. doi:10.1186/gb-2008-9-1-r10
Ulrich J., Dao V. A., Majumdar U., Schmitt-Engel C., Schwirz J., Schultheis D., et al. (2015). Large scale RNAi screen in Tribolium reveals novel target genes for pest control and the proteasome as prime target. BMC Genomics 16, 674. doi:10.1186/s12864-015-1880-y
Varley A. J., Desaulniers J.-P. (2021). Chemical strategies for strand selection in short-interfering RNAs. RSC Adv. 11, 2415–2426. doi:10.1039/D0RA07747J
L. M. Vaschetto (Editor) (2022). “RNAi strategies for pest management: Methods and protocols,” Methods in molecular biology (New York, NY: Springer US). doi:10.1007/978-1-0716-1633-8
Vogel E., Santos D., Mingels L., Verdonckt T.-W., Broeck J. V. (2019). RNA interference in insects: Protecting beneficials and controlling pests. Front. Physiol. 9, 1912. doi:10.3389/fphys.2018.01912
Vyas M., Fisher T. W., He R., Nelson W., Yin G., Cicero J. M., et al. (2015). Asian citrus psyllid expression profiles suggest candidatus liberibacter asiaticus-mediated alteration of adult nutrition and metabolism, and of nymphal development and immunity. PLOS ONE 10, e0130328. doi:10.1371/journal.pone.0130328
Vyas M., Raza A., Ali M. Y., Ashraf M. A., Mansoor S., Shahid A. A., et al. (2017). Knock down of whitefly gut gene expression and mortality by orally delivered gut gene-specific dsRNAs. PLOS ONE 12, e0168921. doi:10.1371/journal.pone.0168921
Weiss P., Huppert S., Kölling R. (2009). Analysis of the dual function of the ESCRT-III protein Snf7 in endocytic trafficking and in gene expression. Biochem. J. 424, 89–97. doi:10.1042/BJ20090957
Willow J., Veromann E. (2021). Highly variable dietary RNAi sensitivity among Coleoptera. Front. Plant Sci. 12, 790816. doi:10.3389/fpls.2021.790816
Wu M., Dong Y., Zhang Q., Li S., Chang L., Loiacono F. V., et al. (2022). Efficient control of Western flower thrips by plastid-mediated RNA interference. Proc. Natl. Acad. Sci. U. S. A. 119, e2120081119. doi:10.1073/pnas.2120081119
Wuriyanghan H., Falk B. W. (2013). RNA interference towards the potato psyllid, Bactericera cockerelli, is induced in plants infected with recombinant tobacco mosaic virus (TMV). PLOS ONE 8, e66050. doi:10.1371/journal.pone.0066050
Wuriyanghan H., Rosa C., Falk B. W. (2011). Oral delivery of double-stranded RNAs and siRNAs induces RNAi effects in the potato/tomato psyllid, Bactericerca cockerelli. PloS One 6, e27736. doi:10.1371/journal.pone.0027736
Yu N., Christiaens O., Liu J., Niu J., Cappelle K., Caccia S., et al. (2013). Delivery of dsRNA for RNAi in insects: An overview and future directions: Delivery of dsRNA for RNAi in insects. Insect Sci. 20, 4–14. doi:10.1111/j.1744-7917.2012.01534.x
Yu H., Yi L., Lu Z. (2022). Silencing of chitin-binding protein with PYPV-rich domain impairs cuticle and wing development in the asian citrus psyllid, Diaphorina citri. Insects 13, 353. doi:10.3390/insects13040353
Yu X., Gowda S., Killiny N. (2017). Double-stranded RNA delivery through soaking mediates silencing of the muscle protein 20 and increases mortality to the Asian citrus psyllid, Diaphorina citri. Pest Manag. Sci. 73, 1846–1853. doi:10.1002/ps.4549
Yu X., Killiny N. (2020). RNA interference-mediated control of Asian citrus psyllid, the vector of the huanglongbing bacterial pathogen. Trop. Plant Pathol. 45, 298–305. doi:10.1007/s40858-020-00356-7
Keywords: biopesticide, dsRNA, huanglongbing, pest control, RNA-interference
Citation: Mondal M, Carver M and Brown JK (2022) Characteristics of environmental RNAi in potato psyllid, Bactericera cockerelli (Sulc) (Hemiptera: Psylloidea: Triozidae). Front. Physiol. 13:931951. doi: 10.3389/fphys.2022.931951
Received: 29 April 2022; Accepted: 20 September 2022;
Published: 18 October 2022.
Edited by:
Jin-Jun Wang, Southwest University, ChinaReviewed by:
Xiudao Yu, Gannan Normal University, ChinaRuobing Guan, Henan Agricultural University, China
Amit Roy, Czech University of Life Sciences Prague, Czechia
Mohammad Mehrabadi, Tarbiat Modares University, Iran
Jaba Jagdish, International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), India
Copyright © 2022 Mondal, Carver and Brown. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Judith K. Brown, jbrown@ag.arizona.edu