Transcriptome comparison between honey bee queen- and worker-destined larvae
Xuan Chen,Yang Hu,Huoqing Zheng,Lianfei Cao,Defang Niu,Dongliang Yu,Yongqiao Sun,Songnian Hu,Fuliang Hu
DOI: https://doi.org/10.1016/j.ibmb.2012.05.004
Abstract:Caste differentiation in the female honey bee is one of the most intriguing polyphenism phenomena. This developmental switch depends on the differential expression of entire suites of the genes involved in the larval fate between the queens and workers. In this study, we compared the transcriptome differences between full-sister queen- (QL) and worker-destined larvae (WL) using high-throughput RNA-Seq. QL and WL at fourth (L4) and fifth instar (L5) were used to prepare four libraries and to generate 50,191,699 (QL4), 57,628,541 (WL4), 56,613,619 (QL5), and 58,626,829 (WL5) usable reads, which were assembled into groups of 7,952, 7,993, 7,971, and 8,023 genes, respectively. The transcriptome changes were investigated using the DEGs Package (DEGseq), which resulted in more than 4,500 differentially expressed genes (DEGs) between the castes. Eight of the DEGs were verified by quantitative real-time RT-PCR (qRT-PCR), and the results supported our sequencing data. All of the DEGs were analysed using Web Gene Ontology Annotation Plot (WEGO) and then mapped using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. These results suggest that over 70% of the DEGs in each instar were more highly expressed in QL than in WL, possibly suggesting that the QL genes had higher transcriptional activity than the WL genes during differentiation. The same gene set is active (but differentially expressed) in both castes, which in turn result in dimorphic females. The L4 stage is a very active gene expression period for both QL and WL before their pupal stage. The activity of the mTOR (a target of rapamycin) encoding gene in the mTOR signalling pathway is higher in QL4 than in WL4, and this difference was no longer present by the L5 feeding stage. The genes down-stream of mTOR maintained this change at the L5 stage. These results could contribute to an in-depth study of the candidate genes during honey bee caste differentiation and improve our current understanding of the polyphenism phenomenon in insects.