Differential gene expression underlying ovarian phenotype determination in honey bee, Apis mellifera L., caste development

Abstract

Adult honey bee queens and workers drastically differ in ovary size. This adult ovary phenotype difference becomes established during the final larval instar, when massive programmed cell death leads to the degeneration of 95–99% of the ovariole anlagen in workers. The higher juvenile hormone (JH) levels in queen larvae protect the ovaries against such degeneration. To gain insights into the molecular architecture underlying this divergence critical for adult caste fate and worker sterility, we performed a microarray analysis on fourth and early fifth instar queen and worker ovaries. For the fourth instar we found nine differentially expressed genes (DEGs) with log2FC > 1.0, but this number increased to 56 in early fifth-instar ovaries. We selected 15 DEGs for quantitative PCR (RT-qPCR) analysis. Nine differed significantly by the variables caste and/or development. Interestingly, genes with enzyme functions were higher expressed in workers, while those related to transcription and signaling had higher transcript levels in queens. For the RT-qPCR confirmed genes we analyzed their response to JH. This revealed a significant up-regulation for two genes, a short chain dehydrogenase reductase (sdr) and a heat shock protein 90 (hsp90). Five other genes, including hsp60 and hexamerin 70b (hex70b), were significantly down-regulated by JH. The sdr gene had previously come up as differentially expressed in other transcriptome analyses on honey bee larvae and heat shock proteins are frequently involved in insect hormone responses, this making them interesting candidates for further functional assays.