Expression of subunits of an insecticide target receptor varies across tissues, life stages, castes, and species of social bees.

Global losses of insects jeopardize ecosystem stability and crop pollination. Robust evidence indicates that insecticides have contributed to these losses. Notably, insecticides targeting nicotinic acetylcholine receptors (nAChRs) have neurotoxic effects on beneficial insects. Because each nAChR consists of five subunits, the alternative arrangements of subunits could create a multitude of receptors differing in structure and function. Therefore, understanding whether the use of subunits varies is essential for evaluating and predicting the effects of insecticides targeting such receptors. To better understand how the use and composition of nAChRs differ within and between insect pollinators, we analysed RNA-seq gene expression data from tissues and castes of Apis mellifera honey bees and life stages and castes of the Bombus terrestris bumble bees. We reveal that all analysed tissues express nAChRs and that relative expression levels of nAChR subunits vary widely across almost all comparisons. Our work thus shows fine-tuned spatial and temporal expression of nAChRs. Given that coexpression of subunits underpins the compositional diversity of functional receptors and that the affinities of insecticides depend on nAChR composition, our findings provide a likely mechanism for the various damaging effects of nAChR-targeting insecticides on insects. Furthermore, our results indicate that the appraisal of insecticide risks should carefully consider variation in molecular targets.