Mating flight causes genome-wide transcriptional changes in sexually mature honeybee queens

Abstract In this study, we analyzed the gene and miRNA expression differences between the courted virgin queen (CVQ) and non-courted virgin queen (NCVQ) of Apis mellifera using a high-throughput sequencing method. Through Digital Gene Expression (DGE) sequencing, 452 genes were differentially expressed, out of which, 90 genes were up-regulated and 362 genes were down-regulated in CVQ compared with NCVQ. Through small RNA sequencing, 27 miRNAs showed significant expression difference between these two samples. Moreover, 9 of the differentially expressed genes are the targets of the 11 differentially expressed miRNAs. Besides, 47 novel miRNA candidates were predicted in these two samples. Our results provided valuable information for understanding the molecular mechanism of the transition to functional queens.