Selection and Evaluation of Reference Genes for miRNA Expression Analysis in Bemisia tabaci Under Insecticide Tolerance

A growing number of studies have focused on the microRNA (miRNA) expression in Bemisia tabaci, one devastating agricultural insect pest of the tropical and subtropical areas for which the primary means of control are insecticides. In studying the genetic underpinnings of insecticide resistance, the choice of stable reference genes for normalizing data plays a key role to acquire unbiased expression profile results from quantitative real-time PCR (qPCR) analysis. Expression profiles of 11 selected reference genes were determined systematically in B. tabaci exposure to 11 insecticides. Furthermore, we assessed the stability of all the selected candidates in relation to other variables including sex, tissue type, and developmental stage. Candidate reference gene validation was conducted by analyzing the let-7-5p expression under various experimental treatments. Five programs BestKeeper, NormFinder, geNorm, △Ct, and RefFinder were applied to verify the accuracy of the selected candidates. Our results displayed that the best choices of the selected candidates for pymetrozine, sulfoxaflor, flonicamid, cyantraniliprole, afidopyropen, and deltamethrin treatment were miR-1-3p and miR-100-5p, U6 and miR-100-5p were best for chlorpyrifos and imidacloprid treatments, and U6 and miR-1-3p were best for flupyradifurone and β-cypermethrin treatments. The reference genes miR-624, miR-252, and miR-275 worked best in adult tissues, miR-100-5p and miR-1-3p worked best in either sex, and miR-624 and miR-11 were best to use across developmental stages. Not even one reference gene was found to be suitable for all experimental conditions. Our results contributed to the growing body of the literature on qPCR reference gene selection under various experimental conditions and facilitate further investigation on gene expression changes in B. tabaci, resulting from pesticide exposure.

Introduction

MicroRNAs (miRNAs) are a type of non-coding small RNA, about 18–25 nucleotides (nt) in length, which originate from stem-loop RNA precursors and have been reported to be produced by most living organisms including plants, animals, and bacteria (Bartel, 2004; Lee et al., 2004). miRNAs play key roles in post-transcriptional gene regulation and affect various biological processes in arthropods such as development, host–pathogen interactions, immunity, and metamorphosis (Asgari, 2013; Li et al., 2019). Additionally, miRNAs participate in insecticide resistance by regulating the target gene expression, and recently, a growing number of regulatory functions have involved in insecticide resistance demonstrated in various insect pests (Zhang et al., 2021).

Expression profiling of miRNAs is essential for studying their functions in insect pests. Currently, quantitative RT-PCR (qPCR) has been proven to be one popular method for analyzing the miRNA expression due to the accuracy, sensitivity, specificity, and low costs (Shi and Chiang, 2005). However, because a variety of factors such as template RNA integrity, contamination of genomic DNA, complementary DNA (cDNA) template quality, and reverse transcription reaction efficiency impact the accuracy of qPCR, suitable reference genes must be selected and validated in order to standardize and correct the data to reduce variability introduced by these factors (Huggett et al., 2005). 5SrRNA and U6 snRNA are commonly utilized as reference genes for analyzing the expression of miRNA in insects; however, these genes cannot be consistently appropriate with all treatments (Lim et al., 2011; Xiang et al., 2014). In recent years, screening of reference genes for studies of miRNA in arthropods including Grapholita molesta, Helicoverpa armigera, Bombus lantschouensis, Plutella xylostella, and Galeruca daurica has been performed and found that the reference genes appropriate for assessing the expression of miRNA rely on the species of insect, developmental stage, and environmental conditions, and sometimes, more than one reference gene is needed to normalize the data on gene expression under various abiotic and biotic conditions (Wang X.-W. et al., 2017; Yang et al., 2017; Dong et al., 2019; Zhang et al., 2019; Wang et al., 2021).

Bemisia tabaci, the tobacco whitefly, is known as one devastating piercing-sucking pest around the world, having recently been reported to more than 700 species of plants, including important economical and horticultural crops (Wang X. et al., 2017; Horowitz et al., 2020). Apart from inflicting mechanical damage on plants, B. tabaci can transmit over 200 plant viruses during feeding (Wei et al., 2017). Since control of whiteflies was largely dependent on long-term application of chemical agents, B. tabaci has developed high or extremely high levels of resistance to most commercialized insecticides (Horowitz et al., 2020). Previously, the screening of suitable reference genes for mRNA in B. tabaci had been conducted (Li et al., 2013), and the results were contributed to the studies on mechanisms of insecticide resistance, plant virus–host interactions, chemical sensing, and host preference in whitefly (Wei et al., 2017; Xu et al., 2019; Zeng et al., 2019; He et al., 2020; Wang et al., 2020a,b; Yang et al., 2020). In recent years, several microRNA profiles of B. tabaci under different treatments have been published (Wang et al., 2016; Li et al., 2018; Hasegawa et al., 2020), but as an essential precondition, suitable reference genes of B. tabaci for analyzing the expression of miRNA have not been selected and evaluated.

In our current work, 11 selected reference genes, including 10 miRNAs (miR-9a, miR-1-3p, miR-11, miR-184, miR-275, miR-100-5p, miR-252a, miR-277, miR-279d, and miR-624) from the B. tabaci databases of small RNA sequencing, have been used for the study of selection and evaluation in previous reports (Wang X.-W. et al., 2017; Yang et al., 2017; Dong et al., 2019; Zhang et al., 2019; Wang et al., 2021). With one commonly utilized reference gene, U6 snRNA, the aforementioned 11 candidates were selected to normalize gene expression data for B. tabaci with insecticide tolerance to 11 popular insecticides used against whiteflies in the field, including pymetrozine, chlorpyrifos, imidacloprid, flupyradifurone, sulfoxaflor, flonicamid, cyantraniliprole, afidopyropen, avermectin, deltamethrin, and β-cypermethrin, and most of them have been recorded as cases of resistance in B. tabaci (Wang et al., 2018a; Wang et al., 2020b,c; Zheng et al., 2021; Wang et al., 2022). The influence of tissue type, sex, and developmental stage was also assessed. The results provided several units of appropriate reference genes for more research on the insecticide-associated gene expression of the whitefly.

Materials and Methods

Insects and Chemicals

The susceptible strain of B. tabaci (MED-S) was sampled in 2009 from damaged poinsettias in the city of Beijing, China (Pan et al., 2011). All populations were reared on the Gossypium hirsutum cotton plants, maintained at temperature 27 ± 1°C, relative humidity 60 ± 5%, and under a photoperiod with 14:10 h light and dark cycle. Insecticides utilized were analytically standardized. Afidopyropen and sulfoxaflor were purchased from Dr. Ehrenstorfer, Germany. Flonicamid, pymetrozine, imidacloprid, β-cypermethrin, cyantraniliprole, deltamethrin, avermectin, flupyradifurone, chlorpyrifos, dimethyl sulfoxide, and Triton X-100 were purchased from Sigma-Aldrich, Shanghai, China.

Determining LC₅₀ Values of Insecticides

All whitefly bioassays were conducted using published methods (Wang et al., 2018b), and statistical analysis of insecticide bioassays to determine the LC₅₀ values of the tested insecticides was performed using Polo Plus. (2002). Adult whiteflies were exposed to the LC₅₀ level concentrations of pymetrozine, chlorpyrifos, imidacloprid, flupyradifurone, sulfoxaflor, flonicamid, cyantraniliprole, afidopyropen, avermectin, deltamethrin, or β-cypermethrin. The tested whitefly adults were exposed to distilled water as the control. After 48 h, all survivors of the tested whiteflies gathered for the next step of work about insecticide tolerance. All adults used in bioassays were less than 7 days old. Males and females were used in a 1:1 ratio.

Sample Collection for Reference Gene Validation in Relation to the Developmental Stage, Tissue, and Sex

Five developmental stages of the whitefly (eggs, first and second, third, fourth stage nymphs, and adults) were sampled and put into a tube, respectively. Three replicates with a total of 2500 eggs, 2000 first and second instar nymphs, 500 third instar nymphs, 400 fourth instar nymphs, and 200 adults were sampled. Three kinds of adult tissues including the abdomen, thorax, and head were collected by the dissection of whitefly adults, and three replicates with a total of 1500 heads, 900 thoraces, and 600 abdomens were sampled. For samples of females and males, the sex of whitefly adults was identified under one stereomicroscope, and gender-identical ones were gathered together, and 200 males or females were collected as one sample. In total, 200 survival adults of B. tabaci after the LC₅₀ treatment of each insecticide were randomly collected, and all the aforementioned samples, after collection, were treated with liquid nitrogen for the extraction of RNA.

Extraction of RNA, Synthesis of cDNA, and Design of Primers

Total RNA samples were extracted from treatments in the use of TRIzol reagent kit (Invitrogen, United States), according to the instruction of the manufacturer, and the quality of all the RNA samples was assessed via agarose gel electrophoresis, and also, an ultraviolet spectrophotometer (Nano-Drop-2000, Thermo Scientific) was utilized to check the concentration and purity of the RNA samples, according to the published method (Yang et al., 2017). Then, with the use of RNA as a template, the steps of reverse transcription were carried out following the specification of the Mir-X miRNA First-Strand Synthesis Kit (TaKaRa, China). A total of 10 miRNAs (miR-1-3p, miR-9a, miR-11, miR-184, miR-100-5p, miR-275, miR-252a, miR-279d, miR-277, and miR-624) and U6 snRNA, one commonly used reference gene, were set as the candidates for evaluation, according to the published database of B. tabaci miRNA-seq and our database (Wang et al., 2016; Li et al., 2018; Hasegawa et al., 2020). Primers of the candidate miRNA were designed in utilizing the software of Primer 5.0 based on the specification of the Mir-X miRNA First-Strand Synthesis Kit (TaKaRa, China), and the reverse primer of all miRNAs was a universal reverse primer (mRQ 3′Primer) from the kit.

Quantitative PCR and Validation of Reference Gene Selection

According to the published method, for measuring the expression of miRNAs (Li et al., 2018), qPCR was conducted using the TB Green^® Premix Ex Taq™ II (Tli RNase H Plus) kit (TaKaRa, China). The total mixture of reaction was 20 µL containing 2 µL of cDNA, 10 µL of TB Green Premix Ex Taq II (Tli RNase H Plus), 6.4 µL of distilled water, 0.8 µL reverse primer, and 0.8 µL forward primer. The conditions of qPCR were 95°C for 30 s, followed by 40 cycles of 95°C for 5 s, and 60°C for 34 s, and the method of 2^-△△Ct was applied to assess the expression of miRNA. A series of 10ⁿ-fold diluted cDNA were utilized to establish the standard curves of five points to assess the amplification efficiency of qPCR (E) and the determination coefficient (R²) of each primer pair. E values were measured on the basis of the equation: E = (10^[−1/slope]-1) × 100 (Radonić et al., 2004). As one target gene, let-7-5p was selected to conduct validation of the candidate reference genes with various treatments. In comparison with the expression of let-7-5p, the most stable gene (NF1), the best combination of genes (NF1–2)/(NF1–3), and the most unstable gene (NF11) were set as reference genes, and then, the expression levels of let-7-5p were measured by the published approach (Livak and Schmittgen, 2001).

Data Analysis

The stability of expression in 11 selected reference genes was assessed using Internet-based calculating software RefFinder, which includes four programs geNorm, △Ct, NormFinder, and BestKeeper. (Vandesompele et al., 2002; Andersen et al., 2004; Pfaffl et al., 2004; Silver et al., 2006; Xie et al., 2012). In accordance with ranking each program, the software of RefFinder was used to conduct calculation of the qPCR results for the comprehensive order of rank (Xie et al., 2012). Considering that the order of expression stability could not indicate the suitable number to choose for normalizing data, the program of geNorm was utilized to confirm the suitable number that need to be utilized within specific treatments, depending on the values of pairwise variation (V). When V_{n/n + 1} < 0.15, the suitable number is N. Analysis of statistical differences was conducted by the use of student’s t test between the most unstable gene and the best combination (SPSS, 2011).

Results

Expression Profiles of Candidate Reference Genes

In order to obtain standard curves of selected candidates, qPCR tests were carried out. As shown in Table 1, the determination coefficient (R²) of the 11 selected reference genes was measured (from 0.9564 to 0.9988), demonstrating that the linear relationship between specific working concentrations of the template and their values of Ct was solid. Values of amplification efficiency (E) were also calculated (from 94.51% to 113.59%), demonstrating the primers’ effectiveness. As shown in Figure 1 and Table 2, the mean Ct values of the candidates ranged from 15.64 to 29.95, and the Ct SD values ranged from 1.04 to 4.35, indicating that miR-277 and miR-11 displayed the smallest and largest variation, respectively, among all the selected candidates.

TABLE 1

TABLE 1. Sequences and amplification characteristics of candidate reference genes.

FIGURE 1

FIGURE 1. Cycle threshold (Ct) values of candidate reference genes in Bemisia tabaci.

TABLE 2

TABLE 2. Cycle threshold (Ct) values of candidate reference genes in B. tabaci under different treatments.

Expression Stability of the Selected Reference Genes Under Insecticide Tolerance and Other Experimental Conditions

As shown in Figure 2, in relation to effects of insecticide tolerance, △Ct, NormFinder, and geNorm analyses indicated that miR-1-3p, miR-100-5p, and miR-277 were the most stable reference genes, yet BestKeeper evaluated miR-100-5p, miR-184, and miR-1-3p rather than miR-100-5p, miR-1-3p, and miR-277. The results of RefFinder showed the order of rank for stability of genes was from the least to most stable, miR-11, miR-624, miR-252a, miR-275, miR-9a, miR-279d, U6, miR-184, miR-277, miR-1-3p, and miR-100-5p. As shown in Figure 3, for other treatment settings, the most stable reference genes were miR-624 for developmental stages and adult tissues and miR-100-5p for gender, respectively. As shown in Supplementary Table S1, all data on the order of rank for stability of genes are presented in the part of Supplementary Material.

FIGURE 2

FIGURE 2. Stability of candidate reference genes in B. tabaci under various tolerance for insecticides. In a RefFinder analysis, increasing Geomean values correspond to decreasing gene expression stability. The Geomean values for the following B. tabaci samples are presented. (A) Pymetrozine: samples treated with pymetrozine; (B) chlorpyrifos: samples treated with chlorpyrifos; (C) imidacloprid: samples treated with imidacloprid; (D) flupyradifurone: samples treated with flupyradifurone; (E) sulfoxaflor: samples treated with sulfoxaflor; (F) flonicamid: samples treated with flonicamid; (G) cyantraniliprole: samples treated with cyantraniliprole; (H) afidopyropen: samples treated with afidopyropen; (I) avermectin: samples treated with avermectin; (J) deltamethrin: samples treated with deltamethrin; (K) β-cypermethrin: samples treated with β-cypermethrin; and (L) insecticides: samples treated with all tested insecticides. The candidate reference genes are as follows: miR-1-3p, miR-9a, miR-11, miR-184, miR-100-5p, miR-275, miR-252a, miR-279d, miR-277, miR-624, and U6 snRNA.

FIGURE 3

FIGURE 3. Stability of candidate reference genes in B. tabaci under various conditions. In a RefFinder analysis, increasing Geomean values correspond to decreasing gene expression stability. The Geomean values for the following B. tabaci samples are presented. (A) Developmental stage: samples for all developmental stages; (B) tissue: samples for different tissues of adults; and (C) adult samples for different genders.

Optimal Number of Reference Genes

Based on the values of the paired variation (Figure 4), the values of V2/3 for all insecticide tolerance, developmental stage, and sex were <0.15 except for in adult tissues, in which the V3/4 value was <0.15. According to the results of overall ranking, the most suitable combinations of reference genes were miR-1-3p and miR-100-5p for insecticide tolerance, miR-100-5p and miR-1-5p for sex, miR-9a and miR-624 for the developmental stage, and miR-624, miR-252, and miR-275 for the tissue type.

FIGURE 4

FIGURE 4. Pairwise variation analysis for an accurate normalization in B. tabaci. (A) Insecticide tolerance of 11 chemical agents and (B) other experimental conditions such as developmental stages, adult tissues, and gender.

Validation of the Selected Reference Genes

In order to better verify the reliability of the 11 candidates, let-7-5p was set as one target gene for validating them (Figure 5). Expression profiles of let-7-5p under tolerance for various insecticides were stable, while normalized using NF1 (miR-100-5p) and NF (1–2) (miR-100-5p and miR-1-3p). Moreover, across developmental stages, while by the use of NF1 (miR-624) and NF(1–2) (miR-624 and miR-9a), expression profiles of let-7-5p increased from third instar nymphs to adults, relative to the first two stages, of which a significant difference with the use of NF1 and NF (1–2) was not observed. Similarly, the expression level of let-7-5p in females was higher than that in males while normalized by the use of NF1 (miR-100-5p) and NF (1–2) (miR-100-5p and miR-1-5p), and no significant difference was observed when using NF1 and NF (1–2). Let-7-5p was most highly expressed in the thorax tissue while normalized with the use of NF1 (miR-624) and NF (1-3) (miR-624, miR-252, and miR-275), and the highest expression occurred in the abdomen using NF11 (miR-279d) in comparison with the normalization using NF1 (miR-624) and NF (1-3) (miR-624, miR-252, and miR-275).

FIGURE 5

FIGURE 5. Validation of reference gene selection. Relative expression levels of the target gene let-7-5p under different treatments. (A) Insecticide tolerance, (B) developmental stages, (C) adult tissues, and (D) genders. The results are displayed as the mean ± based on three independent biological replicates. The t-test was used for statistical difference analysis between the best combination of reference genes and the worst, and statistical differences were denoted by * (p < 0.05).

Discussion

A series of published reports have indicated that miRNAs exert regulatory functions in a variety of insects by regulating the expression of target genes (Li et al., 2019; Zhang et al., 2021). Considering that miRNAs account for a very low percentage of the total RNA profile, it is necessary to use accurate and sophisticated methods for analyzing the expression of less abundant miRNAs and exploring their functions (Sokol, 2008). Thus, qPCR has been universally applied to calculate the miRNA expression, and assessment of appropriate reference genes for normalizing results of qPCR is vital to acquire accurate and consistent results. In the current work, 10 miRNAs and one commonly utilized reference gene were selected as candidates and assessed for their ability to reliably normalize the miRNA expression. We utilized the online program RefFinder which includes the four programs Normfinder, BestKeeper, △Ct, and geNorm to identify and rank the best B. tabaci reference genes to use with various treatments. Based on the data, we found that different programs gave disparate rankings of candidate reference gene stability under identical treatments and distinct rankings of gene stability under different treatments in an identical program. In Galeruca daurica, the best combinations of reference genes were miR-9a-5p and U6 snRNA for different stages of development, miR-305-5p and miR-100-5p for different tissue types, and miR-100-5p and miR-275-3p for various temperature settings (Wang et al., 2021). The best reference gene combinations in Helicoverpa armigera were miR-9 and U6 snRNA for stages of development, miR-305 and miR-100 for tissue types, and miR-100 and miR-92a for temperature treatments, and in Grapholita molesta, the most suitable combinations were miR-16 and U6 snRNA for the stages of development, and miR-281 and U6 snRNA for different settings of temperature (Wang X. et al., 2017).

Our current work assessed the expression stability of 11 candidates for evaluating reference genes in B. tabaci with exposure to 11 pesticides, which are widely used for the management of this destructive insect pest. Furthermore, the stability of the 11 candidates was assessed in various types of adult tissue, stages of development, and sex of B. tabaci. Our results indicated that the miR-1-3p expression was highly stable across the insecticides pymetrozine, sulfoxaflor, flonicamid, afidopyropen, deltamethrin, and β-cypermethrin; U6 snRNA was highly stable across the insecticides imidacloprid, flupyradifurone, and chlorpyrifos, and miR-100-5p was highly stable across the insecticides cyantraniliprole and avermectin. Similarly in H. armigera, appropriate reference genes were miR-9 and U6 snRNA for insecticide treatment (Yang et al., 2017). For the other tested variables, miR-624 was the most suitable for the stages of development and tissue types and miR-100-5p for sex. In G. daurica, miR-100-5p was the most suitable for the tissue type (Wang et al., 2021). Also, U6 snRNA can be utilized as one suitable reference gene for the quantification of miRNA for different stages of development (Feng et al., 2014). Thus, our data further demonstrated that not even one reference gene appropriate for whole treatment of experiments can be found out.

MicroRNA let-7 was initially identified in the nematode Caenorhabditis elegans, and it has been found to function in a variety of biological processes in nematodes, flies, and mammals (Lee et al., 2016). To validate the 11 reference gene candidates, the expression levels of let-7-5p in B. tabaci were measured and compared with, the most suitable combination and the least suitable reference gene, respectively, with various treatments. Our data indicated that let-7-5p shows higher levels of expression in adults than any other developmental stage and lower levels in eggs. Similarly, in G. molesta and G. daurica, a higher expression in adults and lower expression in eggs were observed (Wang X.-W. et al., 2017; Wang et al., 2021). Furthermore, it has been demonstrated that inappropriate reference genes can greatly affect miRNA quantification and give rise to inaccurate results, and several publications report similar phenomena in various insect pests, emphasizing the need to assess stable reference genes for reliable determination of the miRNA expression (Wang X. et al., 2017; Yang et al., 2017; Zhang et al., 2019; Wang et al., 2021). In the current work, the expression levels of let-7-5p normalized with the best reference gene, and the best combination of reference genes was similar in each of the tested treatment. Hence, it can be concluded that although the best reference gene is commonly used for normalization in qPCR under most circumstances, it should still be evaluated under new experimental conditions.

Conclusion

In summary, we assessed the choice of stable reference genes for normalizing data, which plays a key role to acquire unbiased expression profile results from qPCR analysis. Expression profiles of 11 selected reference genes were determined systematically in B. tabaci tolerance to 11 insecticides. Furthermore, we assessed the stability of all the selected candidates in relation to other variables including sex, tissue type, and developmental stage. The results contribute to the growing body of the literature on qPCR reference gene selection of miRNA under various experimental conditions and facilitate further investigation on miRNA regulation of gene expression changes in B. tabaci associated with insecticide tolerance.

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 authors.

Author Contributions

QZ, BG, and RW conceived and designed the study. QZ and BG performed the experiments and analyzed the data with the help of CQ, CL, JW, and RW. RW wrote the first draft of the manuscript. QZ, BG, CQ, JW, and RW participated in manuscript drafting and modification.

Funding

This work was supported by the National Natural Science Foundation of China (32172438 and 31972266).

Conflict of Interest

The 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/fgene.2022.899756/full#supplementary-material

References

Andersen, C. L., Jensen, J. L., and Ørntoft, T. F. (2004). Normalization of Real-Time Quantitative Reverse Transcription-PCR Data: a Model-Based Variance Estimation Approach to Identify Genes Suited for Normalization, Applied to Bladder and colon Cancer Data Sets. Cancer Res. 64 (15), 5245–5250. doi:10.1158/0008-5472.CAN-04-0496

PubMed Abstract | CrossRef Full Text | Google Scholar

Asgari, S. (2013). MicroRNA Functions in Insects. Insect Biochem. Mol. Biol. 43 (4), 388–397. doi:10.1016/j.ibmb.2012.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartel, D. P. (2004). MicroRNAs. Cell 116 (2), 281–297. doi:10.1016/s0092-8674(04)00045-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, J., Li, J., Huang, J., and Wu, J. (2019). Identification of Suitable Reference Genes for miRNA Quantitation in Bumblebee (Hymenoptera: Apidae) Response to Reproduction. Apidologie 50, 40–50. doi:10.1007/s13592-018-0616-9

CrossRef Full Text | Google Scholar

Feng, B., Liang, P., and Gao, X. W. (2014). Cloning and Evaluation of U6 snRNA Gene as a Reference Gene for miRNA Quantification in Plutella Xylostella (Lepidoptera: Plutellidae). Acta Ent. Sin. 57 (3), 286–292. doi:10.1093/jisesa/iex007

CrossRef Full Text | Google Scholar

Hasegawa, D. K., Shamimuzzaman, M., Chen, W., Simmons, A. M., Fei, Z., and Ling, K.-S. (2020). Deep Sequencing of Small RNAs in the Whitefly Bemisia Tabaci Reveals Novel MicroRNAs Potentially Associated with Begomovirus Acquisition and Transmission. Insects 11 (9), 562. doi:10.3390/insects11090562

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Y.-Z., Wang, Y.-M., Yin, T.-Y., Fiallo-Olivé, E., Liu, Y.-Q., Hanley-Bowdoin, L., et al. (2020). A Plant DNA Virus Replicates in the Salivary Glands of its Insect Vector via Recruitment of Host DNA Synthesis Machinery. Proc. Natl. Acad. Sci. U.S.A. 117 (29), 16928–16937. doi:10.1073/pnas.1820132117

PubMed Abstract | CrossRef Full Text | Google Scholar

Horowitz, A. R., Ghanim, M., Roditakis, E., Nauen, R., and Ishaaya, I. (2020). Insecticide Resistance and its Management in Bemisia Tabaci Species. J. Pest Sci. 93, 893–910. doi:10.1007/s10340-020-01210-0

CrossRef Full Text | Google Scholar

Huggett, J., Dheda, K., Bustin, S., and Zumla, A. (2005). Real-time RT-PCR Normalisation; Strategies and Considerations. Genes Immun. 6 (4), 279–284. doi:10.1038/sj.gene.6364190

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H., Han, S., Kwon, C. S., and Lee, D. (2016). Biogenesis and Regulation of the Let-7 miRNAs and Their Functional Implications. Protein Cell 7 (2), 100–113. doi:10.1007/s13238-015-0212-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, Y. S., Nakahara, K., Pham, J. W., Kim, K., He, Z., Sontheimer, E. J., et al. (2004). Distinct Roles for Drosophila Dicer-1 and Dicer-2 in the siRNA/miRNA Silencing Pathways. Cell 117 (1), 69–81. doi:10.1016/s0092-8674(04)00261-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Wei, X., Ding, T., and Chu, D. (2018). Genome-wide Profiling of Cardinium-Responsive microRNAs in the Exotic Whitefly, Bemisia Tabaci (Gennadius) Biotype Q. Front. Physiol. 9, 1580. doi:10.3389/fphys.2018.01580

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M.-z., Xiao, H.-m., He, K., and Li, F. (2019). Progress and Prospects of Noncoding RNAs in Insects. J. Integr. Agric. 18 (4), 729–747. doi:10.1016/S2095-3119(18)61976-8

CrossRef Full Text | Google Scholar

Li, R., Xie, W., Wang, S., Wu, Q., Yang, N., Yang, X., et al. (2013). Reference Gene Selection for qRT-PCR Analysis in the Sweetpotato Whitefly, Bemisia Tabaci (Hemiptera: Aleyrodidae). Plos One 8 (1), e53006. doi:10.1371/journal.pone.0053006

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, Q. E., Zhou, L., Ho, Y. K., Wan, G., and Too, H. P. (2011). snoU6 and 5S RNAs Are Not Reliable miRNA Reference Genes in Neuronal Differentiation. Neuroscience 199, 32–43. doi:10.1016/j.neuroscience.2011.10.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 25 (4), 402–408. doi:10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, H., Chu, D., Ge, D., Wang, S., Wu, Q., Xie, W., et al. (2011). Further Spread of and Domination by Bemisia Tabaci (Hemiptera: Aleyrodidae) Biotype Q on Field Crops in China. Jnl. Econ. Entom. 104 (3), 978–985. doi:10.1603/ec11009

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaffl, M. W., Tichopad, A., Prgomet, C., and Neuvians, T. P. (2004). Determination of Stable Housekeeping Genes, Differentially Regulated Target Genes and Sample Integrity: BestKeeper - Excel-based Tool Using Pair-wise Correlations. Biotechnol. Lett. 26 (6), 509–515. doi:10.1023/b:bile.0000019559.84305.47

PubMed Abstract | CrossRef Full Text | Google Scholar

Polo Plus (2002). A User’s Guide to Probit or Logit Analysis. Berkeley, CA, USA: LeOra Software.

Google Scholar

Radonić, A., Thulke, S., Mackay, I. M., Landt, O., Siegert, W., and Nitsche, A. (2004). Guideline to Reference Gene Selection for Quantitative Real-Time PCR. Biochem. Biophysical Res. Commun. 313 (4), 856–862. doi:10.1016/j.bbrc.2003.11.177

CrossRef Full Text | Google Scholar

Shi, R., and Chiang, V. L. (2005). Facile Means for Quantifying microRNA Expression by Real-Time PCR. Biotechniques 39 (4), 519–525. doi:10.2144/000112010

PubMed Abstract | CrossRef Full Text | Google Scholar

Silver, N., Best, S., Jiang, J., and Thein, S. L. (2006). Selection of Housekeeping Genes for Gene Expression Studies in Human Reticulocytes Using Real-Time PCR. BMC Mol. Biol 7, 33. doi:10.1186/1471-2199-7-33

PubMed Abstract | CrossRef Full Text | Google Scholar

Sokol, N. S. (2008). An Overview of the Identification, Detection, and Functional Analysis of Drosophila MicroRNAs. Methods Mol. Biol. 420, 319–334. doi:10.1007/978-1-59745-583-1_20

PubMed Abstract | CrossRef Full Text | Google Scholar

SPSS (2011). Software package version 13.0 for windows. Chicago, IL.

Google Scholar

Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., et al. (2002). Accurate Normalization of Real-Time Quantitative RT-PCR Data by Geometric Averaging of Multiple Internal Control Genes. Genome Biol. 3 (7), research00341. doi:10.1186/gb-2002-3-7-research0034

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, B., Wang, L., Chen, F., Yang, X., Ding, M., Zhang, Z., et al. (2016). MicroRNA Profiling of the Whitefly Bemisia Tabaci Middle East-Aisa Minor I Following the Acquisition of Tomato Yellow Leaf Curl China Virus. Virol. J. 13, 20. doi:10.1186/s12985-016-0469-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H. C., Li, L., Li, Y. Y., Tan, Y., and Pang, B. P. (2021). Evaluation of Reference Genes for miRNA Expression Analysis in Galeruca Daurica (Coleoptera: Chrysomelidae) Using qRT-PCR. Entomol. Res. 51 (8), 393–402. doi:10.1111/1748-5967.12537

CrossRef Full Text | Google Scholar

Wang, R., Che, W., Wang, J., and Luo, C. (2020c). Monitoring Insecticide Resistance and Diagnostics of Resistance Mechanisms in Bemisia Tabaci Mediterranean (Q Biotype) in China. Pestic. Biochem. Physiol. 163, 117–122. doi:10.1016/j.pestbp.2019.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., Fang, Y., Mu, C., Qu, C., Li, F., Wang, Z., et al. (2018b). Baseline Susceptibility and Cross-Resistance of Cycloxaprid, a Novel Cis-Nitromethylene Neonicotinoid Insecticide, in Bemisia Tabaci MED from China. Crop Prot. 110, 283–287. doi:10.1016/j.cropro.2017.02.012

CrossRef Full Text | Google Scholar

Wang, R., Gao, B., Che, W., Qu, C., Zhou, X., and Luo, C. (2022). First Report of Field Resistance to Afidopyropen, the Novel Pyropene Insecticide, on Bemisia Tabaci Mediterranean (Q Biotype) from China. Agronomy 12, 724. doi:10.3390/agronomy12030724

CrossRef Full Text | Google Scholar

Wang, R., Wang, J.-D., Che, W.-N., and Luo, C. (2018a). First Report of Field Resistance to Cyantraniliprole, a New Anthranilic Diamide Insecticide, on Bemisia Tabaci MED in China. J. Integr. Agric. 17 (1), 158–163. doi:10.1016/S2095-3119(16)61613-1

CrossRef Full Text | Google Scholar

Wang, R., Wang, J., Che, W., Fang, Y., and Luo, C. (2020b). Baseline Susceptibility and Biochemical Mechanism of Resistance to Flupyradifurone in Bemisia Tabaci. Crop Prot. 132, 105132. doi:10.1016/j.cropro.2020.105132

CrossRef Full Text | Google Scholar

Wang, R., Wang, J., Zhang, J., Che, W., Feng, H., and Luo, C. (2020a). Characterization of Flupyradifurone Resistance in the Whitefly Bemisia Tabaci Mediterranean (Q Biotype). Pest Manag. Sci. 76 (12), 4286–4292. doi:10.1002/ps.5995

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X.-W., Li, P., and Liu, S.-S. (2017a). Whitefly Interactions with Plants. Curr. Opin. Insect Sci. 19, 70–75. doi:10.1016/j.cois.2017.02.001

CrossRef Full Text | Google Scholar

Wang, X., Li, Y., Zhang, J., Zhang, Q., Liu, X., and Li, Z. (2017b). De Novo characterization of microRNAs in oriental Fruit Moth Grapholita Molesta and Selection of Reference Genes for Normalization of microRNA Expression. Plos One 12 (2), e0171120. doi:10.1371/journal.pone.0171120

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, J., He, Y.-Z., Guo, Q., Guo, T., Liu, Y.-Q., Zhou, X.-P., et al. (2017). Vector Development and Vitellogenin Determine the Transovarial Transmission of Begomoviruses. Proc. Natl. Acad. Sci. U.S.A. 114 (26), 6746–6751. doi:10.1073/pnas.1701720114.Xiang

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, M., Zeng, Y., Yang, R. R., Xu, H. F., Chen, Z., Zhong, J., et al. (2014). U6 Is Not a Suitable Endogenous Control for the Quantification of Circulating microRNAs. Biochem. Biophys. Res. Commun. 454 (1), 210–214. doi:10.1016/j.bbrc.2014.10.064)

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, F., Xiao, P., Chen, D., Xu, L., and Zhang, B. (2012). miRDeepFinder: a miRNA Analysis Tool for Deep Sequencing of Plant Small RNAs. Plant Mol. Biol. 80, 75–84. doi:10.1007/s11103-012-9885-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, H.-X., Qian, L.-X., Wang, X.-W., Shao, R.-X., Hong, Y., Liu, S.-S., et al. (2019). A Salivary Effector Enables Whitefly to Feed on Host Plants by Eliciting Salicylic Acid-Signaling Pathway. Proc. Natl. Acad. Sci. U.S.A. 116 (2), 490–495. doi:10.1073/pnas.1714990116

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Deng, S., Wei, X., Yang, J., Zhao, Q., Yin, C., et al. (2020). MAPK-directed Activation of the Whitefly Transcription Factor CREB Leads to P450-Mediated Imidacloprid Resistance. Proc. Natl. Acad. Sci. U.S.A. 117 (19), 10246–10253. doi:10.1073/pnas.1913603117

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Li, Z., Cao, J., Li, Y., Li, H., Yang, Q., et al. (2017). Identification and Evaluation of Suitable Reference Genes for Normalization of microRNA Expression in Helicoverpa Armigera (Lepidoptera: Noctuidae) Using Quantitative Real-Time PCR. J. Insect Sci. 17 (2), 33. doi:10.1093/jisesa/iex007

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, Y., Yang, Y. T., Wu, Q. J., Wang, S. L., Xie, W., and Zhang, Y. J. (2019). Genome‐wide Analysis of Odorant‐binding Proteins and Chemosensory Proteins in the Sweet Potato Whitefly, Bemisia Tabaci. Insect Sci. 26 (4), 620–634. doi:10.1111/1744-7917.12576

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Jing, X., Chen, W., Bai, J., Vasseur, L., He, W., et al. (2019). Selection of Reference Genes for Expression Analysis of Plant-Derived microRNAs in Plutella Xylostella Using qRT-PCR and ddPCR. Plos One 14 (8), e0220475. doi:10.1371/journal.pone.0220475

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Dou, W., Taning, C. N. T., Smagghe, G., and Wang, J.-J. (2021). Regulatory Roles of microRNAs in Insect Pests: Prospective Targets for Insect Pest Control. Curr. Opin. Biotechnol. 70, 158–166. doi:10.1016/j.copbio.2021.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, H., Xie, W., Fu, B., Xiao, S., Tan, X., Ji, Y., et al. (2021). Annual Analysis of Field‐evolved Insecticide Resistance in Bemisia Tabaci across China. Pest Manag. Sci. 77 (6), 2990–3001. doi:10.1002/ps.6338

PubMed Abstract | CrossRef Full Text | Google Scholar