Transcriptomic Analysis Reveals Patterns of Expression of Stage-Specific Genes in Early Apis cerana Embryos

Abstract

Background/Objectives: Apis cerana development is described as comprising four stages: embryo, larva, pupa, and adult. There are significant differences between workers and drones in terms of physiological functions and social roles, and the formation of the organ primordia occurs during the embryonic stage. Therefore, the objective of this study is to investigate the differential expression of and alternative splicing of genes in worker and drone embryos and to explain their unique developmental patterns. Methods: Long-read sequencing (PacBio Iso-Seq) and short-read sequencing (Illumina RNA-Seq) were used to investigate worker and drone embryo gene expression differences in A. cerana across five developmental points (12, 24, 36, 48, and 60 h). Results: The study identified 59,254 common isoforms, with 5744 and 5106 isoforms specific to worker and drone embryos, respectively. Additionally, a new transcript of the csd gene was identified. The number of differentially expressed genes (3391) and differential splicing events (470 genes) peaked at the 24-h embryonic stage. Differential splicing events of csd, dsx, and Y-y were observed in the worker and drone embryos. Conclusions: The gene expression results indicated that the 24-h embryonic point is a critical period for the expression of genes related to developmental and behavioral differences between workers and drones. The findings provide a theoretical basis for future research on the developmental differences between workers and drones.