Identifying potential insecticide resistance markers through genomic‐level comparison of Bemisia tabaci (Gennadius) lines

The invasive whitefly (Bemisia tabaci) MED is one of the most economically damaging plant pests. The extensive use of insecticide over decades has led to that the invasive B. tabaci MED has developed resistance to a wide range of insecticide classes, but little is known about the genetic background associated with resistance. To this end, we conducted a comparative genome-wide analysis of single-base nucleotide polymorphisms between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line collected in 1976. First, low-coverage genome sequencings were conducted on DNA isolated from individual whiteflies. The sequencing results were evaluated using an available B. tabaci MED genome as a reference. Significant genetic differences were discovered between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line based on the principal component analyses. Top GO categories and KEGG pathways that might be involved in insecticide resistance development were identified, and several of them have not been previously associated with resistance. Additionally, we identified several genetic loci with novel variations including Cytochrome P450 monooxygenases (P450s), UDP-glucuronosyltransferases (UGTs), Glutathione S-transferases (GSTs), esterase, carboxyl-esterases (COE), ABC transporters, fatty acyl-CoA reductase, voltage-gated sodium channels, GABA receptor, and cuticle proteins (CPs) that were previously reported to have close associations with pesticide resistance in well-studied insect groups that provide an essential resource for the design of insecticide resistance-linked loci arrays insecticide. Our results was obtained solely on resequencing genome data sets, more pesticide bio-assays combined with omics datasets should be further used to verify the markers identified here.