Honeybee development consists of four stages: embryo, larva, pupa and adult. Embryogenesis, a key process of cell division and differentiation, takes 3 days in honeybees. However, the embryonic transcriptome and the dynamic regulation of embryonic transcription are still largely uncharacterized in honeybees, especially in the Asian honeybee (Apis cerana). Here, we employed high-quality RNA-seq to explore the transcriptome of Asian honeybee embryos at three ages, approximately 24, 48 and 72 h (referred to as Day1, Day2 and Day3, respectively). Nine embryo samples, three from each age, were collected for RNA-seq. According to the staging scheme of honeybee embryos and the morphological features we observed, our Day1, Day2 and Day3 embryos likely corresponded to the late stage four, stage eight and stage ten development stages, respectively. Hierarchical clustering and principal component analysis showed that same-age samples were grouped together, and the Day2 samples had a closer relationship with the Day3 samples than the Day1 samples. Finally, a total of 18,284 genes harboring 55,646 transcripts were detected in the A. cerana embryos, of which 44.5% consisted of the core transcriptome shared by all three ages of embryos. A total of 4088 upregulated and 3046 downregulated genes were identified among the three embryo ages, of which 2010, 3177 and 1528 genes were upregulated and 2088, 2294 and 303 genes were downregulated from Day1 to Day2, from Day1 to Day3 and from Day2 to Day3, respectively. The downregulated genes were mostly involved in cellular, biosynthetic and metabolic processes, gene expression and protein localization, and macromolecule modification; the upregulated genes mainly participated in cell development and differentiation, tissue, organ and system development, and morphogenesis. Interestingly, several biological processes related to the response to and detection of light stimuli were enriched in the first-day A. cerana embryogenesis but not in the Apis mellifera embryogenesis, which was valuable for further investigations. Our transcriptomic data substantially expand the number of known transcribed elements in the A. cerana genome and provide a high-quality view of the transcriptome dynamics of A. cerana embryonic development.

Apis mellifera and Apis cerana are two sibling species of Apidae. Apis cerana is adept at collecting sporadic nectar in mountain and forest region and exhibits stiffer hardiness and acarid resistance as a result of natural selection, whereas Apis mellifera has the advantage of producing royal jelly. To identify differentially expressed genes (DEGs) that affect the development of hypopharyngeal gland (HG) and/or the secretion of royal jelly between these two honeybee species, we performed a digital gene expression (DGE) analysis of the HGs of these two species at three developmental stages (newly emerged worker, nurse and forager). Twelve DGE-tag libraries were constructed and sequenced using the total RNA extracted from the HGs of newly emerged workers, nurses, and foragers of Apis mellifera and Apis cerana. Finally, a total of 1482 genes in Apis mellifera and 1313 in Apis cerana were found to exhibit an expression difference among the three developmental stages. A total of 1417 DEGs were identified between these two species. Of these, 623, 1072, and 462 genes showed an expression difference at the newly emerged worker, nurse, and forager stages, respectively. The nurse stage exhibited the highest number of DEGs between these two species and most of these were found to be up-regulated in Apis mellifera. These results suggest that the higher yield of royal jelly in Apis mellifera may be due to the higher expression level of these DEGs. In this study, we investigated the DEGs between the HGs of two sibling honeybee species (Apis mellifera and Apis cerana). Our results indicated that the gene expression difference was associated with the difference in the royal jelly yield between these two species. These results provide an important clue for clarifying the mechanisms underlying hypopharyngeal gland development and the production of royal jelly.

Honey bees are important species for the study of epigenetics. Female honey bee larvae with the same genotype can develop into phenotypically distinct organisms (sterile workers and fertile queens) depending on conditions such as diet. Previous studies have shown that DNA methylation and histone modification can establish distinct gene expression patterns, leading to caste differentiation. It is unclear whether the histone methylation modification H3K4me1 can also impact caste differentiation. In this study, we analyzed genome-wide H3K4me1 modifications in both queen and worker larvae and found that H3K4me1 marks are more abundant in worker larvae than in queen larvae at both the second and fourth instars, and many genes associated with caste differentiation are differentially methylated. Notably, caste-specific H3K4me1 in promoter regions can direct worker development. Thus, our results suggest that H3K4me1 modification may act as an important regulatory factor in the establishment and maintenance of caste-specific transcriptional programs in honey bees; however, the potential influence of other epigenetic modifications cannot be excluded.

RNA interference (RNAi) has been used successfully to reduce target gene expression and induce specific phenotypes in several species. It has proved useful as a tool to investigate gene function and has the potential to manage pest populations and reduce disease pathogens. However, it is not known whether different administration methods are equally effective at interfering with genes in bees. Therefore, we compared the effects of feeding and injection of small interfering RNA (siRNA) on the messenger RNA (mRNA) levels of alpha-aminoadipic semialdehyde dehydrogenase (ALDH7A1), 4-coumarate-CoA ligase (4CL), and heat shock protein 70 (HSP70). Both feeding and injection of siRNA successfully knocked down the gene but feeding required more siRNA than the injection. Our results suggest that both feeding and injection of siRNA effectively interfere with brain genes in bees. The appropriateness of each method would depend on the situation.

The honeybee is a model organism for studying learning and memory formation and its underlying molecular mechanisms. While DNA methylation is well studied in caste differentiation, its role in learning and memory is not clear in honeybees. Here, we analyzed genome-wide DNA methylation changes during olfactory learning and memory process in A. mellifera using whole genome bisulfite sequencing (WGBS) method. A total of 853 significantly differentially methylated regions (DMRs) and 963 differentially methylated genes (DMGs) were identified. We discovered that 440 DMRs of 648 genes were hypermethylated and 274 DMRs of 336 genes were hypomethylated in trained group compared to untrained group. Of these DMGs, many are critical genes involved in learning and memory, such as Creb, GABA _B R and Ip3k, indicating extensive involvement of DNA methylation in honeybee olfactory learning and memory process. Furthermore, key enzymes for histone methylation, RNA editing and miRNA processing also showed methylation changes during this process, implying that DNA methylation can affect learning and memory of honeybees by regulating other epigenetic modification processes.