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.

Ant colonies are higher-level organisms consisting of specialized reproductive and non-reproductive individuals that differentiate early in development, similar to germ–soma segregation in bilateral Metazoa. Analogous to diverging cell lines, developmental differentiation of individual ants has often been considered in epigenetic terms but the sets of genes that determine caste phenotypes throughout larval and pupal development remain unknown. Here, we reconstruct the individual developmental trajectories of two ant species, Monomorium pharaonis and Acromyrmex echinatior, after obtaining >1,400 whole-genome transcriptomes. Using a new backward prediction algorithm, we show that caste phenotypes can be accurately predicted by genome-wide transcriptome profiling. We find that caste differentiation is increasingly canalized from early development onwards, particularly in germline individuals (gynes/queens) and that the juvenile hormone signalling pathway plays a key role in this process by regulating body mass divergence between castes. We quantified gene-specific canalization levels and found that canalized genes with gyne/queen-biased expression were enriched for ovary and wing functions while canalized genes with worker-biased expression were enriched in brain and behavioural functions. Suppression in gyne larvae of Freja, a highly canalized gyne-biased ovary gene, disturbed pupal development by inducing non-adaptive intermediate phenotypes between gynes and workers. Our results are consistent with natural selection actively maintaining canalized caste phenotypes while securing robustness in the life cycle ontogeny of ant colonies.

During colony relocation, the selection of a new nest involves exploration and assessment of potential sites followed by colony movement on the basis of a collective decision making process. Hygiene and pathogen load of the potential nest sites are factors worker scouts might evaluate, given the high risk of epidemics in group-living animals. Choosing nest sites free of pathogens is hypothesized to be highly efficient in invasive ants as each of their introduced populations is often an open network of nests exchanging individuals (unicolonial) with frequent relocation into new nest sites and low genetic diversity, likely making these species particularly vulnerable to parasites and diseases. We investigated the nest site preference of the invasive pharaoh ant, Monomorium pharaonis, through binary choice tests between three nest types: nests containing dead nestmates overgrown with sporulating mycelium of the entomopathogenic fungus Metarhizium brunneum (infected nests), nests containing nestmates killed by freezing (uninfected nests), and empty nests. In contrast to the expectation pharaoh ant colonies preferentially (84%) moved into the infected nest when presented with the choice of an infected and an uninfected nest. The ants had an intermediate preference for empty nests. Pharaoh ants display an overall preference for infected nests during colony relocation. While we cannot rule out that the ants are actually manipulated by the pathogen, we propose that this preference might be an adaptive strategy by the host to “immunize” the colony against future exposure to the same pathogenic fungus.