Abstract
In eusocial insects, the molecular basis of worker reproductivity, including how it changes with eusocial complexity, remains relatively poorly understood. To address this, we used mRNA-seq to isolate genes differentially expressed between ovary-active and ovary-inactive workers in the intermediately eusocial bumblebee Bombus terrestris. By comparisons with data from the advanced eusocial honeybee Apis mellifera, which shows reduced worker reproductivity, we characterized gene expression differences associated with change in worker reproductivity as a function of eusocial complexity. By comparisons with genes associated with queen-worker caste development in B. terrestris larvae, we tested the behavioral-morphological caste homology hypothesis, which proposes co-option of genes influencing reproductive division of labor in adults in morphological caste evolution. We conducted comparisons having isolated genes expressed in B. terrestris worker-laid eggs to remove the potential confound caused by gene expression in eggs. Gene expression differences between the B. terrestris worker phenotypes were mainly in fat body and ovary, not brain. Many genes (86%) more highly expressed in ovary of ovary-active workers were also expressed in worker-laid eggs, confirming egg-expressed genes were potentially confounding. Comparisons across B. terrestris and A. mellifera, and with B. terrestris larvae, returned significant percentage overlaps in differentially expressed genes and/or enriched Gene Ontology terms, suggesting conserved gene functions underpin worker reproductivity as it declines with increasing eusocial complexity and providing support for the behavioral-morphological caste homology hypothesis. Therefore, within bees, both a degree of conserved gene use and gene co-option appear to underlie the molecular basis of worker reproductivity and morphological caste evolution.