Abstract
The grooming behavior of honeybees serves as a crucial auto-protective mechanism against Varroa mite infestations. Compared to Apis mellifera, Apis cerana demonstrates more effective grooming behavior in removing Varroa mites from the bodies of infested bees. However, the underlying mechanisms regulating grooming behavior remain elusive. In this study, we evaluated the efficacy of the auto-grooming behavior between A. cerana and A. mellifera and employed RNA-sequencing technology to identify differentially expressed genes (DEGs) in bee brains with varying degrees of grooming behavior intensity. We observed that A. cerana exhibited a higher frequency of mite removal between day 5 and day 15 compared to A. mellifera, with day-9 bees showing the highest frequency of mite removal in A. cerana. RNA-sequencing results revealed the differential expression of the HTR2A and SLC17A8 genes in A. cerana and the CCKAR and TpnC47D genes in A. mellifera. Subsequent homology analysis identified the HTR2A gene and SLC17A8 gene of A. cerana as homologous to the HTR2A gene and SLC17A7 gene of A. mellifera. These DEGs are annotated in the neuroactive ligand-receptor interaction pathway, the glutamatergic synaptic pathway, and the calcium signaling pathway. Moreover, CCKAR, TpnC47D, HTR2A, and SLC17A7 may be closely related to the auto-grooming behavior of A. mellifera, conferring resistance against Varroa infestation. Our results further explain the relationship between honeybee grooming behavior and brain function at the molecular level and provide a reference basis for further studies of the mechanism of honeybee grooming behavior.