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