MicroRNAs (miRNAs) play a vital role in the nerve regulation of honey bees (Apis mellifera). This study aims to investigate the differences in expression of miRNAs in a honey bee’s brain for olfactory learning tasks and to explore their potential role in a honey bee’s olfactory learning and memory. In this study, 12 day old honey bees with strong and weak olfactory performances were utilized to investigate the influence of miRNAs on olfactory learning behavior. The honey bee brains were dissected, and a small RNA-seq technique was used for high-throughput sequencing. The data analysis of the miRNA sequences revealed that 14 differentially expressed miRNAs (DEmiRNAs) between the two groups, strong (S) and weak (W), for olfactory performance in honey bees were identified, which included seven up-regulated and seven down-regulated. The qPCR verification results of the 14 miRNAs showed that four miRNAs (miR-184-3p, miR-276-3p, miR-87-3p, and miR-124-3p) were significantly associated with olfactory learning and memory. The target genes of these DEmiRNAs were subjected to the GO database annotation and KEGG pathway enrichment analyses. The functional annotation and pathway analysis showed that the neuroactive ligand-receptor interaction pathway, oxidative phosphorylation, biosynthesis of amino acids, pentose phosphate pathway, carbon metabolism, and terpenoid backbone biosynthesis may be a great important pathway related to olfactory learning and memory in honey bees. Our findings together further explained the relationship between olfactory performance and the brain function of honey bees at the molecular level and provides a basis for further study on miRNAs related to olfactory learning and memory in honey bees.

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.

Stingers, evolved from ovipositors, are an important defense organ for the Apidae, Vespidae, and Formicidae species. However, the molecular mechanism of stinger development remains unclear. Here, we show that the earliest time point for the appearance of stingers in Apis mellifera is at the 1-day-old worker pupal stage based on morphological observations and anatomy from the pre-pupal to adult stages. To discover the genes related to stinger development, we first comprehensively compared the stinger transcriptome at different stages and screened 1282, 186, and 166 highly expressed genes in the stingers of 1- and 5-day-old worker pupae and newly emerged worker bees (NEBs), respectively, then identified 25 DEGs involved in the early stage of stinger development. We found that Dll was a key candidate gene in the early development of A. mellifera stingers by combining analyses of the protein–protein interaction network and spatiotemporal expression patterns. An RNAi experiment showed that about 20% of individuals exhibited tip bending in the piercing parts of their stingers in the Dll-dsRNA-treated group, with the morphology presenting as side–side or front–back tip bending. This indicates that Dll plays a vital role in the early development of A. mellifera stingers. Together, our study provides insight into the molecular mechanism of Hymenoptera stinger development and an inspiration for the molecular breeding of gentle honeybee species with stinger abnormalities.

Honeybees are an indispensable pollinator in nature with pivotal ecological, economic, and scientific value. However, a full-length transcriptome for Apis mellifera, assembled with the advanced third-generation nanopore sequencing technology, has yet to be reported. Here, nanopore sequencing of the midgut tissues of uninoculated and Nosema ceranae-inoculated A. mellifera workers was conducted, and the full-length transcriptome was then constructed and annotated based on high-quality long reads. Next followed improvement of sequences and annotations of the current reference genome of A. mellifera. A total of 5,942,745 and 6,664,923 raw reads were produced from midguts of workers at 7 days post-inoculation (dpi) with N. ceranae and 10 dpi, while 7,100,161 and 6,506,665 raw reads were generated from the midguts of corresponding uninoculated workers. After strict quality control, 6,928,170, 6,353,066, 5,745,048, and 6,416,987 clean reads were obtained, with a length distribution ranging from 1 kb to 10 kb. Additionally, 16,824, 17,708, 15,744, and 18,246 full-length transcripts were respectively detected, including 28,019 nonredundant ones. Among these, 43,666, 30,945, 41,771, 26,442, and 24,532 full-length transcripts could be annotated to the Nr, KOG, eggNOG, GO, and KEGG databases, respectively. Additionally, 501 novel genes (20,326 novel transcripts) were identified for the first time, among which 401 (20,255), 193 (13,365), 414 (19,186), 228 (12,093), and 202 (11,703) were respectively annotated to each of the aforementioned five databases. The expression and sequences of three randomly selected novel transcripts were confirmed by RT-PCR and Sanger sequencing. The 5′ UTR of 2082 genes, the 3′ UTR of 2029 genes, and both the 5′ and 3′ UTRs of 730 genes were extended. Moreover, 17,345 SSRs, 14,789 complete ORFs, 1224 long non-coding RNAs (lncRNAs), and 650 transcription factors (TFs) from 37 families were detected. Findings from this work not only refine the annotation of the A. mellifera reference genome, but also provide a valuable resource and basis for relevant molecular and -omics studies.