Background: Reproductive division of labor in eusocial insects is a striking example of a shared genetic background giving rise to alternative phenotypes, namely queen and worker castes. Queen and worker phenotypes play major roles in the evolution of eusocial insects. Their behavior, morphology and physiology underpin many ecologically relevant colony-level traits, which evolved in parallel in multiple species. Results: Using queen and worker transcriptomic data from 16 ant species we tested the hypothesis that conserved sets of genes are involved in ant reproductive division of labor. We further hypothesized that such sets of genes should also be involved in the parallel evolution of other key traits. We applied weighted gene co-expression network analysis, which clusters co-expressed genes into modules, whose expression levels can be summarized by their ‘eigengenes’. Eigengenes of most modules were correlated with phenotypic differentiation between queens and workers. Furthermore, eigengenes of some modules were correlated with repeated evolution of key phenotypes such as complete worker sterility, the number of queens per colony, and even invasiveness. Finally, connectivity and expression levels of genes within the co-expressed network were strongly associated with the strength of selection. Although caste-associated sets of genes evolve faster than non-caste-associated, we found no evidence for queen- or worker-associated co-expressed genes evolving faster than one another. Conclusions: These results identify conserved functionally important genomic units that likely serve as building blocks of phenotypic innovation, and allow the remarkable breadth of parallel evolution seen in ants, and possibly other eusocial insects as well.

Vitellogenin (Vg), a storage protein, has been extensively studied for its egg-yolk precursor role, and it has been suggested to be fundamentally involved in caste differences in social insects. More than one Vg copy has been reported in several oviparous species, including ants. However, the number and function of different Vgs, their phylogenetic relatedness, and their role in reproductive queens and nonreproductive workers have been studied in few species only. We studied caste-biased expression of Vgs in seven Formica ant species. Only one copy of conventional Vg was identified in Formica species, and three Vg homologs, derived from ancient duplications, which represent yet undiscovered Vg-like genes. We show that each of these Vg-like genes is present in all studied Hymenoptera and some of them in other insects as well. We show that after each major duplication event, at least one of the Vg-like genes has experienced a period of positive selection. This, combined with the observation that the Vg-like genes have acquired or lost specific protein domains suggests sub- or neofunctionalization between Vg and the duplicated genes. In contrast to earlier studies, Vg was not consistently queen biased in its expression, and the caste bias of the three Vg-like genes was highly variable among species. Furthermore, a truncated and Hymenoptera-specific Vg-like gene, Vg-like-C, was consistently worker biased. Multispecies comparisons are essential for Vg expression studies, and for gene expression studies in general, as we show that expression and also, putative functions cannot be generalized even among closely related species.

Social insect colonies are characterized by division of labour, and extensive morphological, physiological and behavioural differences between queens and workers. The storage protein vitellogenin (Vg) affects multiple aspects of social insect life histories, and has been suggested as a key player for caste differentiation and maintenance. Recently, three genes homologous to Vg have been described in the antFormica exsecta. Their role is currently unclear but their structural variation suggests variable functions. We examined the expression patterns of the conventional Vg and the three Vg-like genes using qRT-PCR in the common black antFormica fuscabetween queens and workers, between nurse and foragers workers, and across social contexts (queenless vs. queenright nests), and sampling time. As expected, we found a significant queen caste and nurse task-related increase for the conventional Vg, while Vg-like-C displayed a consistent forager-biased expression pattern. Task (forager vs. nurse) was the only factor that explained expression variation among workers in any of the studied genes. The removal of the queen did not affect expression, although the proportion of fertile nurses increased in queenless nests. The observed expression biases suggest that inFormica fusca,the ancestral duplication has led to alternative social functions for Vg-like genes across castes and tasks. To get a broader picture of the role of gene duplications in social evolution and the roles of Vg-like genes in caste differentiation and maintenance, how these genes achieve these roles at a molecular level need to be investigated further.