Differential gene expression responsible for caste determination at both larval and adult stages of Bombus terrestris

The influence of gene expression on female caste determination in social Hymenoptera is one of the best-characterized examples of developmental polyphenism. We used mRNA-seq to profile differentially expressed genes (DEGs) between workers and queens at four development stages (early larvae: 1-day-old larvae; mid larvae: 5-day-old larvae; late larvae: 9-day-old larvae; and 7-day-old adult bee) of Bombus terrestris. When DEGs of different female castes were compared, it was found that there were more DEGs at the mid larval and adult stages than at other stages. Caste differentiation may be linked to candidate genes such as Vitellogenin, Kr-h1, Dnmt3, Hexamerin, Yellow, and Chymotrypsin-2. Additionally, the gene expression profiles were observed differently among the four development stages. Comparing one stage to the next showed that, more number of DEGs were found in queens than workers. In the investigation of how gene expression affects larval caste determination, WGCNA showed that five modules were closely associated with larval caste determination. The functions of genes were mainly enriched in energy metabolism and energy production within mitochondria, suggesting that the energy requirements during the larval development process varied between caste determination. The current study sheds light on the patterns of gene expression in the female caste of Bombus terrestris, which will be useful for future research on polyphenism through differential gene expressions.