LAUSR

Abstract Neonicotinoid insecticides are now the most widely used insecticides in the world. Previous studies have indicated that sublethal doses of neonicotinoids impair learning, memory capacity, foraging, and immunocompetence in honey bees (Apis mellifera, Linnaeus) (Hymenoptera: Apidae). Despite these, few studies have been carried out on the molecular effects of neonicotinoids. In this study, we focus on the second-generation neonicotinoid thiamethoxam, which is currently widely used in agriculture to protect crops. Using high-throughput RNA-Seq, we investigated the transcriptome profile of honey bees after subchronic exposure to 10 ppb thiamethoxam over 10 d. In total, 609 differentially expressed genes (DEGs) were identified, of which 225 were upregulated and 384 were downregulated. Several genes, including vitellogenin, CSP3, defensin-1, Mrjp1, and Cyp6as5 were selected and further validated using real-time quantitative polymerase chain reaction assays. The functions of some DEGs were identified, and Gene Ontology–enrichment analysis showed that the enriched DEGs were mainly linked to metabolism, biosynthesis, and translation. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that thiamethoxam affected biological processes including ribosomes, the oxidative phosphorylation pathway, tyrosine metabolism pathway, pentose and glucuronate interconversions, and drug metabolism. Overall, our results provide a basis for understanding the molecular mechanisms of the complex interactions between neonicotinoid insecticides and honey bees.