Towards reconstructing the ancestral brain gene-network regulating caste differentiation in ants

Juvenile hormone as a key regulator for asymmetric caste differentiation in ants

Caste differentiation involves many functional traits that diverge during larval growth and metamorphosis to produce adults irreversibly adapted to reproductive division of labor. Investigating developmental differentiation is important for general biological understanding and has increasingly been explored for social phenotypes that diverge in parallel from similar genotypes. Here, we use Monomorium pharaonis ants to investigate the extent to which canalized worker development can be shifted toward gyne (virgin-queen) phenotypes by juvenile hormone (JH) treatment. We show that excess JH can activate gyne-biased development in workers so that wing-buds, ocelli, antennal and genital imaginal discs, flight muscles, and gyne-like fat bodies and brains emerge after pupation. However, ovary development remained unresponsive to JH treatment, indicating that JH-sensitive germline sequestration happens well before somatic differentiation. Our findings reveal important qualitative restrictions in the extent to which JH treatment can redirect larval development and that these constraints are independent of body size. Our findings corroborate that JH is a key hormone for inducing caste differentiation but show that this process can be asymmetric for higher colony-level germline versus somatic caste differentiation in superorganisms as defined a century ago by Wheeler. We quantified gene expression changes in response to JH treatment throughout development and identified a set of JH-sensitive genes responsible for the emergence of gyne-like somatic traits. Our study suggests that the gonadotropic role of JH in ovary maturation has shifted from the individual level in solitary insects to the colony level in an evolutionary-derived and highly polygynous superorganism like the pharaoh ant.

Gene regulation supporting sociality shared across lineages and variation in complexity

Across evolutionary lineages, insects vary in social complexity, from those that exhibit extended parental care to those with elaborate divisions of labor. Here, we synthesize the sociogenomic resources from hundreds of species to describe common gene regulatory mechanisms in insects that regulate social organization across phylogeny and levels of social complexity. Different social phenotypes expressed by insects can be linked to the organization of co-expressing gene networks and features of the epigenetic landscape. Insect sociality also stems from processes like the emergence of parental care and the decoupling of ancestral genetic programs. One underexplored avenue is how variation in a group's social environment affects the gene expression of individuals. Additionally, an experimental reduction of gene expression would demonstrate how the activity of specific genes contributes to insect social phenotypes. While tissue specificity provides greater localization of the gene expression underlying social complexity, emerging transcriptomic analysis of insect brains at the cellular level provides even greater resolution to understand the molecular basis of social insect evolution.

A genetic toolkit underlying the queen phenotype in termites with totipotent workers

Social insect castes (e.g., queens, workers) are prime examples of phenotypic plasticity (i.e., different phenotypes arising from the same genotype). Yet, the mechanisms that give rise to highly fertile, long-lived queens versus non-reproducing, short-lived workers are not well understood. Recently, a module of co-expressed genes has been identified that characterizes queens compared to workers of the termite Cryptotermes secundus (Kalotermitidae): the Queen Central Module (QCM). We tested whether the QCM is shared in termite species, in which queens gradually develop via early larval and late larval instars, the latter functioning as totipotent workers (linear development). Similar as in C. secundus, gene expression profiles revealed an enrichment of QCM genes in Zootermopsis angusticollis queens, a species from another termite family (Archotermopsidae). The expression of these QCM genes became gradually enriched during development from early larval instars via workers to queens. Thus, our results support the hypothesis of a conserved genetic toolkit that characterizes termite queens with gradual linear development. Our data also imply a strong caste-specific tissue specificity with the QCM signal being restricted to head-prothorax tissues in termite queens. This tissue-specific expression of key aging-related genes might have facilitated the evolution of a long lifespan in termite queens.

Molecular signatures of alternative reproductive strategies in a facultatively social hover wasp

Social insect reproductives and non-reproductives represent ideal models with which to understand the expression and regulation of alternative phenotypes. Most research in this area has focused on the developmental regulation of reproductive phenotypes in obligately social taxa such as honey bees, while relatively few studies have addressed the molecular correlates of reproductive differentiation in species in which the division of reproductive labour is established only in plastic dominance hierarchies. To address this knowledge gap, we generate the first genome for any stenogastrine wasp and analyse brain transcriptomic data for non-reproductives and reproductives of the facultatively social species Liostenogaster flavolineata, a representative of one of the simplest forms of social living. By experimentally manipulating the reproductive 'queues' exhibited by social colonies of this species, we show that reproductive division of labour in this species is associated with transcriptomic signatures that are more subtle and variable than those observed in social taxa in which colony living has become obligate; that variation in gene expression among non-reproductives reflects their investment into foraging effort more than their social rank; and that genes associated with reproductive division of labour overlap to some extent with those underlying division of labour in the separate polistine origin of wasp sociality but only explain a small portion of overall variation in this trait. These results indicate that broad patterns of within-colony transcriptomic differentiation in this species are similar to those in Polistinae but offer little support for the existence of a strongly conserved 'toolkit' for sociality.

Neuropeptide regulations on behavioral plasticity in social insects

Social insects demonstrate remarkable behavioral flexibility in response to complex external and social environments. One of the most striking examples of this adaptability is the context-dependent division of labor among workers of bees and ants. Neuropeptides, the brain's most diverse group of messenger molecules, play an essential role in modulating this phenotypic plasticity related to labor division in social insects. Integrated omics research and mass spectrometry imaging technology have greatly accelerated the identification and spatiotemporal analysis of neuropeptides. Moreover, key roles of several neuropeptides in age- and caste-dependent behavioral plasticity have been uncovered. This review summarizes recent advances in the characterization, expression, distribution, and functions of neuropeptides in controlling behavioral plasticity in social insects, particularly bees and ants. The article concludes with a discussion of future directions and challenges in understanding the regulation of social behavior by neuropeptides.

Meta-Analysis of Public RNA Sequencing Data Revealed Potential Key Genes Associated with Reproductive Division of Labor in Social Hymenoptera and Termites

Eusociality in insects has evolved independently many times. One of the most notable characteristics of eusociality is the reproductive division of labor. In social insects, the reproductive division of labor is accomplished by queens and workers. Transcriptome analyses of queens and workers have been conducted for various eusocial species. However, the genes that regulate the reproductive division of labor across all or multiple eusocial species have not yet been fully elucidated. Therefore, we conducted a meta-analysis using publicly available RNA-sequencing data from four major groups of social insects. In this meta-analysis, we collected 258 pairs (queen vs. worker) of RNA-sequencing data from 34 eusocial species. The meta-analysis identified a total of 20 genes that were differentially expressed in queens or workers. Out of these, 12 genes have not previously been reported to be involved in the reproductive division of labor. Functional annotation of these 20 genes in other organisms revealed that they could be regulators of behaviors and physiological states related to the reproductive division of labor. These 20 genes, revealed using massive datasets of numerous eusocial insects, may be key regulators of the reproductive division of labor.

Molecular patterns and processes in evolving sociality: lessons from insects

Social insects have provided some of the clearest insights into the origins and evolution of collective behaviour. Over 20 years ago, Maynard Smith and Szathmáry defined the most complex form of insect social behaviour-superorganismality-among the eight major transitions in evolution that explain the emergence of biological complexity. However, the mechanistic processes underlying the transition from solitary life to superorganismal living in insects remain rather elusive. An overlooked question is whether this major transition arose via incremental or step-wise modes of evolution. We suggest that examination of the molecular processes underpinning different levels of social complexity represented across the major transition from solitary to complex sociality can help address this question. We present a framework for using molecular data to assess to what extent the mechanistic processes that take place in the major transition to complex sociality and superorganismality involve nonlinear (implying step-wise evolution) or linear (implying incremental evolution) changes in the underlying molecular mechanisms. We assess the evidence for these two modes using data from social insects and discuss how this framework can be used to test the generality of molecular patterns and processes across other major transitions. This article is part of a discussion meeting issue 'Collective behaviour through time'.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Social complexity, life-history and lineage influence the molecular basis of castes in vespid wasps

A key mechanistic hypothesis for the evolution of division of labour in social insects is that a shared set of genes co-opted from a common solitary ancestral ground plan (a genetic toolkit for sociality) regulates caste differentiation across levels of social complexity. Using brain transcriptome data from nine species of vespid wasps, we test for overlap in differentially expressed caste genes and use machine learning models to predict castes using different gene sets. We find evidence of a shared genetic toolkit across species representing different levels of social complexity. We also find evidence of additional fine-scale differences in predictive gene sets, functional enrichment and rates of gene evolution that are related to level of social complexity, lineage and of colony founding. These results suggest that the concept of a shared genetic toolkit for sociality may be too simplistic to fully describe the process of the major transition to sociality.

Co-expression Gene Networks and Machine-learning Algorithms Unveil a Core Genetic Toolkit for Reproductive Division of Labour in Rudimentary Insect Societies

The evolution of eusociality requires that individuals forgo some or all their own reproduction to assist the reproduction of others in their group, such as a primary egg-laying queen. A major open question is how genes and genetic pathways sculpt the evolution of eusociality, especially in rudimentary forms of sociality-those with smaller cooperative nests when compared with species such as honeybees that possess large societies. We lack comprehensive comparative studies examining shared patterns and processes across multiple social lineages. Here we examine the mechanisms of molecular convergence across two lineages of bees and wasps exhibiting such rudimentary societies. These societies consist of few individuals and their life histories range from facultative to obligately social. Using six species across four independent origins of sociality, we conduct a comparative meta-analysis of publicly available transcriptomes. Standard methods detected little similarity in patterns of differential gene expression in brain transcriptomes among reproductive and non-reproductive individuals across species. By contrast, both supervised machine learning and consensus co-expression network approaches uncovered sets of genes with conserved expression patterns among reproductive and non-reproductive phenotypes across species. These sets overlap substantially, and may comprise a shared genetic "toolkit" for sociality across the distantly related taxa of bees and wasps and independently evolved lineages of sociality. We also found many lineage-specific genes and co-expression modules associated with social phenotypes and possible signatures of shared life-history traits. These results reveal how taxon-specific molecular mechanisms complement a core toolkit of molecular processes in sculpting traits related to the evolution of eusociality.

Canalized gene expression during development mediates caste differentiation in ants

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.

A single-cell transcriptomic atlas tracking the neural basis of division of labour in an ant superorganism

Ant colonies with permanent division of labour between castes and highly distinct roles of the sexes have been conceptualized to be superorganisms, but the cellular and molecular mechanisms that mediate caste/sex-specific behavioural specialization have remained obscure. Here we characterized the brain cell repertoire of queens, gynes (virgin queens), workers and males of Monomorium pharaonis by obtaining 206,367 single-nucleus transcriptomes. In contrast to Drosophila, the mushroom body Kenyon cells are abundant in ants and display a high diversity with most subtypes being enriched in worker brains, the evolutionarily derived caste. Male brains are as specialized as worker brains but with opposite trends in cell composition with higher abundances of all optic lobe neuronal subtypes, while the composition of gyne and queen brains remained generalized, reminiscent of solitary ancestors. Role differentiation from virgin gynes to inseminated queens induces abundance changes in roughly 35% of cell types, indicating active neurogenesis and/or programmed cell death during this transition. We also identified insemination-induced cell changes probably associated with the longevity and fecundity of the reproductive caste, including increases of ensheathing glia and a population of dopamine-regulated Dh31-expressing neurons. We conclude that permanent caste differentiation and extreme sex-differentiation induced major changes in the neural circuitry of ants.

Using single-cell transcriptomics, the authors generate a brain cell atlas for the pharaoh ant including individuals of different sexes and castes and show changes in cell composition underlying division of labour and reproductive specialization.

Validation of Potential Reference Genes for Real-Time qPCR Analysis in Pharaoh Ant, Monomorium pharaonis (Hymenoptera: Formicidae)

Ants are highly diverse social insects living in colonies consisted of up to millions of individuals with reproductive division of labors. Due to the interests in disclosing the genetic and epigenetic regulation mechanisms underlying the distinct developmental trajectories between castes and division of labor in colonies, many ant species have recently been established as laboratory models for evolutionary development and social behavior studies. These functional studies often request a precise quantification of the relative gene expression level, which relies on a stably expressed reference genes for normalization. A core set of reliable reference genes for this purpose however has not been established yet in ants. In the present study, we tested the expression patterns and amplification efficiencies of 12 abundantly expressed candidate genes in Monomorium pharaonis, one of the few ant species that are suitable for laboratory rearing and experimentation. We quantified the expression levels of these genes by RT-qPCR in seven different conditions: embryo development, sexual development, worker development, adult phenotypes, tissues, and two abiotic manipulative treatments in pharaoh ant. Finally, five genes, elongation factor-1 alpha (EF1A), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), TATA-box-binding protein (TATA), tubulin gamma-2 chain-like (TBLg2), heat shock protein 67B2-like (HSP67) were found to be the most stable reference genes across seven conditions. We also identified the most stable reference genes applicable for each distinct condition and the optimal number of reference genes entailed were evaluated. Our study validates reliable reference genes for RT-qPCR analysis which lays the foundation for future studies in pharaoh ant.

Independent variations in genome-wide expression, alternative splicing, and DNA methylation in brain tissues among castes of the buff-tailed bumblebee, Bombus terrestris

Caste differentiation in social hymenopterans is an intriguing example of phenotypic plasticity. However, the co-ordination among gene regulatory factors to mediate caste differentiation remains inconclusive. In this study, we determined the role of gene regulation and related epigenetic processes in pre-imaginal caste differentiation in the primitively eusocial bumblebee Bombus terrestris. By combining RNA-Seq data from Illumina and PacBio and accurately quantifying methylation at whole-genomic base pair resolution, we found that queens, workers, and drones mainly differentiate in gene expression but not in alternative splicing and DNA methylation. Gynes are the most distinct with the lowest global level of whole-genomic methylation and with the largest number of caste-specific transcripts and alternative splicing events. By contrast, workers exhibit few uniquely expressed or alternatively spliced genes. Moreover, several genes involved in hormone and neurotransmitter metabolism are related to caste differentiation, whereas several neuropeptides are linked with sex differentiation. Despite little genome-wide association among differential gene expression, splicing, and differential DNA methylation, the overlapped gene ontology (GO) terms point to nutrition-related activity. Therefore, variations in gene regulation correlate with the behavioral differences among castes and highlight the specialization of toolkit genes in bumblebee gynes at the beginning of the adult stage.

Ants’ Personality and Its Dependence on Foraging Styles: Research Perspectives

The paper is devoted to analyzing consistent individual differences in behavior, also known as “personalities,” in the context of a vital ant task—the detection and transportation of food. I am trying to elucidate the extent to which collective cognition is individual-based and whether a single individual’s actions can suffice to direct the entire colony or colony units. The review analyzes personalities in various insects with different life cycles and provides new insights into the role of individuals in directing group actions in ants. Although it is widely accepted that, in eusocial insects, colony personality emerges from the workers’ personalities, there are only a few examples of investigations of personality at the individual level. The central question of the review is how the distribution of behavioral types and cognitive responsibilities within ant colonies depends on a species’ foraging style. In the context of how workers’ behavioral traits display during foraging, a crucial question is what makes an ant a scout that discovers a new food source and mobilizes its nestmates. In mass recruiting, tandem-running, and even in group-recruiting species displaying leadership, the division of labor between scouts and recruits appears to be ephemeral. There is only little, if any, evidence of ants’ careers and behavioral consistency as leaders. Personal traits characterize groups of individuals at the colony level but not performers of functional roles during foraging. The leader-scouting seems to be the only known system that is based on a consistent personal difference between scouting and foraging individuals.

Relaxed selection underlies genome erosion in socially parasitic ant species

Inquiline ants are highly specialized and obligate social parasites that infiltrate and exploit colonies of closely related species. They have evolved many times convergently, are often evolutionarily young lineages, and are almost invariably rare. Focusing on the leaf-cutting ant genus Acromyrmex, we compared genomes of three inquiline social parasites with their free-living, closely-related hosts. The social parasite genomes show distinct signatures of erosion compared to the host lineages, as a consequence of relaxed selective constraints on traits associated with cooperative ant colony life and of inquilines having very small effective population sizes. We find parallel gene losses, particularly in olfactory receptors, consistent with inquiline species having highly reduced social behavioral repertoires. Many of the genomic changes that we uncover resemble those observed in the genomes of obligate non-social parasites and intracellular endosymbionts that branched off into highly specialized, host-dependent niches.

Comparative transcriptomic analysis of the mechanisms underpinning ageing and fecundity in social insects

The exceptional longevity of social insect queens despite their lifelong high fecundity remains poorly understood in ageing biology. To gain insights into the mechanisms that might underlie ageing in social insects, we compared gene expression patterns between young and old castes (both queens and workers) across different lineages of social insects (two termite, two bee and two ant species). After global analyses, we paid particular attention to genes of the insulin/insulin-like growth factor 1 signalling (IIS)/target of rapamycin (TOR)/juvenile hormone (JH) network, which is well known to regulate lifespan and the trade-off between reproduction and somatic maintenance in solitary insects. Our results reveal a major role of the downstream components and target genes of this network (e.g. JH signalling, vitellogenins, major royal jelly proteins and immune genes) in affecting ageing and the caste-specific physiology of social insects, but an apparently lesser role of the upstream IIS/TOR signalling components. Together with a growing appreciation of the importance of such downstream targets, this leads us to propose the TI-J-LiFe (TOR/IIS-JH-Lifespan and Fecundity) network as a conceptual framework for understanding the mechanisms of ageing and fecundity in social insects and beyond. This article is part of the theme issue 'Ageing and sociality: why, when and how does sociality change ageing patterns?'

An ancient, conserved gene regulatory network led to the rise of oral venom systems

Although oral venom systems are ecologically important characters, how they originated is still unclear. In this study, we show that oral venom systems likely originated from a gene regulatory network conserved across amniotes. This network, which we term the “metavenom network,” comprises over 3,000 housekeeping genes coexpressed with venom and play a role in protein folding and modification. Comparative transcriptomics revealed that the network is conserved between venom glands of snakes and salivary glands of mammals. This suggests that while these tissues have evolved different functions, they share a common regulatory core, that persisted since their common ancestor. We propose several evolutionary mechanisms that can utilize this common regulatory core to give rise to venomous animals from their nonvenomous ancestors.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Transcriptomic analyses of the termite, Cryptotermes secundus, reveal a gene network underlying a long lifespan and high fecundity

Organisms are typically characterized by a trade-off between fecundity and longevity. Notable exceptions are social insects. In insect colonies, the reproducing caste (queens) outlive their non-reproducing nestmate workers by orders of magnitude and realize fecundities and lifespans unparalleled among insects. How this is achieved is not understood. Here, we identified a single module of co-expressed genes that characterized queens in the termite species Cryptotermes secundus. It encompassed genes from all essential pathways known to be involved in life-history regulation in solitary model organisms. By manipulating its endocrine component, we tested the recent hypothesis that re-wiring along the nutrient-sensing/endocrine/fecundity axis can account for the reversal of the fecundity/longevity trade-off in social insect queens. Our data from termites do not support this hypothesis. However, they revealed striking links to social communication that offer new avenues to understand the re-modelling of the fecundity/longevity trade-off in social insects.

Gene expression during larval caste determination and differentiation in intermediately eusocial bumblebees, and a comparative analysis with advanced eusocial honeybees

The queen‐worker caste system of eusocial insects represents a prime example of developmental polyphenism (environmentally‐induced phenotypic polymorphism) and is intrinsic to the evolution of advanced eusociality. However, the comparative molecular basis of larval caste determination and subsequent differentiation in the eusocial Hymenoptera remains poorly known. To address this issue within bees, we profiled caste‐associated gene expression in female larvae of the intermediately eusocial bumblebee Bombus terrestris. In B. terrestris, female larvae experience a queen‐dependent period during which their caste fate as adults is determined followed by a nutrition‐sensitive period also potentially affecting caste fate but for which the evidence is weaker. We used mRNA‐seq and qRT‐PCR validation to isolate genes differentially expressed between each caste pathway in larvae at developmental stages before and after each of these periods. We show that differences in gene expression between caste pathways are small in totipotent larvae, then peak after the queen‐dependent period. Relatively few novel (i.e., taxonomically‐restricted) genes were differentially expressed between castes, though novel genes were significantly enriched in late‐instar larvae in the worker pathway. We compared sets of caste‐associated genes in B. terrestris with those reported from the advanced eusocial honeybee, Apis mellifera, and found significant but relatively low levels of overlap of gene lists between the two species. These results suggest both the existence of low numbers of shared toolkit genes and substantial divergence in caste‐associated genes between Bombus and the advanced eusocial Apis since their last common eusocial ancestor.

High-quality chromosome-level genome assembly and full-length transcriptome analysis of the pharaoh ant Monomorium pharaonis

Background

Ants with complex societies have fascinated scientists for centuries. Comparative genomic and transcriptomic analyses across ant species and castes have revealed important insights into the molecular mechanisms underlying ant caste differentiation. However, most current ant genomes and transcriptomes are highly fragmented and incomplete, which hinders our understanding of the molecular basis for complex ant societies.

Findings

By hybridizing Illumina, Pacific Biosciences, and Hi-C sequencing technologies, we de novo assembled a chromosome-level genome for Monomorium pharaonis, with a scaffold N50 of 27.2 Mb. Our new assembly provides better resolution for the discovery of genome rearrangement events at the chromosome level. Analysis of full-length isoform sequencing (ISO-seq) suggested that ∼15 Gb of ISO-seq data were sufficient to cover most expressed genes, but the number of transcript isoforms steadily increased with sequencing data coverage. Our high-depth ISO-seq data greatly improved the quality of gene annotation and enabled the accurate detection of alternative splicing isoforms in different castes of M. pharaonis. Comparative transcriptome analysis across castes based on the ISO-seq data revealed an unprecedented number of transcript isoforms, including many caste-specific isoforms. We also identified a number of conserved long non-coding RNAs that evolved specifically in ant lineages and several that were conserved across insect lineages.

Conclusions

We produced a high-quality chromosome-level genome for M. pharaonis, which significantly improved previous short-read assemblies. Together with full-length transcriptomes for all castes, we generated a highly accurate annotation for this ant species. These long-read sequencing results provide a useful resource for future functional studies on the genetic mechanisms underlying the evolution of social behaviors and organization in ants.

The gene expression network regulating queen brain remodeling after insemination and its parallel use in ants with reproductive workers

Caste differentiation happens early in development to produce gynes as future colony germlines and workers as present colony soma. However, gynes need insemination to become functional queens, a transition that initiates reproductive role differentiation relative to unmated gynes. Here, we analyze the anatomy and transcriptomes of brains during this differentiation process within the reproductive caste of Monomorium pharaonis Insemination terminated brain growth, whereas unmated control gynes continued to increase brain volume. Transcriptomes revealed a specific gene regulatory network (GRN) mediating both brain anatomy changes and behavioral modifications. This reproductive role differentiation GRN hardly overlapped with the gyne-worker caste differentiation GRN, but appears to be also used by distantly related ants where workers became germline individuals after the queen caste was entirely or partially lost. The genes corazonin and neuroparsin A in the anterior neurosecretory cells were overexpressed in individuals with reduced or nonreproductive roles across all four ant species investigated.

Chromatin accessibility and transcriptome landscapes of Monomorium pharaonis brain

The emergence of social organization (eusociality) is a major event in insect evolution. Although previous studies have investigated the mechanisms underlying caste differentiation and social behavior of eusocial insects including ants and honeybees, the molecular circuits governing sociality in these insects remain obscure. In this study, we profiled the transcriptome and chromatin accessibility of brain tissues in three Monomorium pharaonis ant castes: queens (including mature and un-mated queens), males and workers. We provide a comprehensive dataset including 16 RNA-sequencing and 16 assay for transposase accessible chromatin (ATAC)-sequencing profiles. We also demonstrate strong reproducibility of the datasets and have identified specific genes and open chromatin regions in the genome that may be associated with the social function of these castes. Our data will be a valuable resource for further studies of insect behaviour, particularly the role of brain in the control of eusociality.

Distributed physiology and the molecular basis of social life in eusocial insects

The traditional focus of physiological and functional genomic research is on molecular processes that play out within a single multicellular organism. In the colonial (eusocial) insects such as ants, bees, and termites, molecular and behavioral responses of interacting nestmates are tightly linked, and key physiological processes are regulated at the scale of the colony. Such colony-level physiological processes regulate nestmate physiology in a distributed fashion, through various social communication mechanisms. As a result of physiological decentralization over evolutionary time, organismal mechanisms, for example related to pheromone detection, hormone signaling, and neural signaling pathways, are deployed in novel contexts to influence nestmate and colony traits. Here we explore how functional genomic, physiological, and behavioral studies can benefit from considering the traits of eusocial insects in this light. We highlight functional genomic work exploring how nestmate-level and colony-level traits arise and are influenced by interactions among physiologically-specialized nestmates of various developmental stages. We also consider similarities and differences between nestmate-level (organismal) and colony-level (superorganismal) physiological processes, and make specific hypotheses regarding the physiology of eusocial taxa. Integrating theoretical models of distributed systems with empirical functional genomics approaches will be useful in addressing fundamental questions related to the evolution of eusociality and collective behavior in natural systems.

Has gene expression neofunctionalization in the fire ant antennae contributed to queen discrimination behavior?

Queen discrimination behavior in the fire ant Solenopsis invicta maintains its two types of societies: colonies with one (monogyne) or many (polygyne) queens, yet the underlying genetic mechanism is poorly understood. This behavior is controlled by two supergene alleles, SB and Sb, with ~600 genes. Polygyne workers, having either the SB/SB or SB/Sb genotype, accept additional SB/Sb queens into their colonies but kill SB/SB queens. In contrast, monogyne workers, all SB/SB, reject all additional queens regardless of genotype. Because the SB and Sb alleles have suppressed recombination, determining which genes within the supergene mediate this differential worker behavior is difficult. We hypothesized that the alternate worker genotypes sense queens differently because of the evolution of differential expression of key genes in their main sensory organ, the antennae. To identify such genes, we sequenced RNA from four replicates of pooled antennae from three classes of workers: monogyne SB/SB, polygyne SB/SB, and polygyne SB/Sb. We identified 81 differentially expressed protein-coding genes with 13 encoding potential chemical metabolism or perception proteins. We focused on the two odorant perception genes: an odorant receptor SiOR463 and an odorant-binding protein SiOBP12. We found that SiOR463 has been lost in the Sb genome. In contrast, SiOBP12 has an Sb-specific duplication, SiOBP12b', which is expressed in the SB/Sb worker antennae, while both paralogs are expressed in the body. Comparisons with another fire ant species revealed that SiOBP12b' antennal expression is specific to S. invicta and suggests that queen discrimination may have evolved, in part, through expression neofunctionalization.

Leveraging technological innovations to investigate evolutionary transitions to eusociality

The study of the major transition to eusociality presents several challenges to researchers, largely resulting from the importance of complex behavioral phenotypes and the shift from individual to group level selection. These challenges are being met with corresponding technological improvements. Advances in resource development for non-model taxa, behavioral tracking, nucleic acid sequencing, and reverse genetics are facilitating studies of hypotheses that were previously intractable. These innovations are resulting in the development of new model systems tailored to the exploration of specific behavioral phenotypes and the querying of underlying molecular mechanisms that drive eusocial behaviors. Here, we present a brief overview of how methodological innovations are advancing our understanding of the evolution of eusociality.

Rate variation in the evolution of non-coding DNA associated with social evolution in bees

The evolutionary origins of eusociality represent increases in complexity from individual to caste-based, group reproduction. These behavioural transitions have been hypothesized to go hand in hand with an increased ability to regulate when and where genes are expressed. Bees have convergently evolved eusociality up to five times, providing a framework to test this hypothesis. To examine potential links between putative gene regulatory elements and social evolution, we compare alignable, non-coding sequences in 11 diverse bee species, encompassing three independent origins of reproductive division of labour and two elaborations of eusocial complexity. We find that rates of evolution in a number of non-coding sequences correlate with key social transitions in bees. Interestingly, while we find little evidence for convergent rate changes associated with independent origins of social behaviour, a number of molecular pathways exhibit convergent rate changes in conjunction with subsequent elaborations of social organization. We also present evidence that many novel non-coding regions may have been recruited alongside the origin of sociality in corbiculate bees; these loci could represent gene regulatory elements associated with division of labour within this group. Thus, our findings are consistent with the hypothesis that gene regulatory innovations are associated with the evolution of eusociality and illustrate how a thorough examination of both coding and non-coding sequence can provide a more complete understanding of the molecular mechanisms underlying behavioural evolution. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Convergent eusocial evolution is based on a shared reproductive groundplan plus lineage-specific plastic genes

Eusociality has convergently evolved multiple times, but the genomic basis of caste-based division of labor and degree to which independent origins of eusociality have utilized common genes remain largely unknown. Here we characterize caste-specific transcriptomic profiles across development and adult body segments from pharaoh ants (Monomorium pharaonis) and honey bees (Apis mellifera), representing two independent origins of eusociality. We identify a substantial shared core of genes upregulated in the abdomens of queen ants and honey bees that also tends to be upregulated in mated female flies, suggesting that these genes are part of a conserved insect reproductive groundplan. Outside of this shared groundplan, few genes are differentially expressed in common. Instead, the majority of the thousands of caste-associated genes are plastically expressed, rapidly evolving, and relatively evolutionarily young. These results emphasize that the recruitment of both highly conserved and lineage-specific genes underlie the convergent evolution of novel traits such as eusociality.

Eusocial caste systems have evolved independently multiple times. Here, Warner et al. investigate the amount of shared vs. lineage-specific genes involved in the evolution of caste in pharaoh ants and honey bees by comparing transcriptomes across tissues, developmental stages, and castes.

Gene expression is more strongly associated with behavioural specialization than with age or fertility in ant workers

The ecological success of social insects is based on division of labour, not only between queens and workers, but also among workers. Whether a worker tends the brood or forages is influenced by age, fertility and nutritional status, with brood carers being younger, more fecund and more corpulent. Here, we experimentally disentangle behavioural specialization from age and fertility in Temnothorax longispinosus ant workers and analyse how these parameters are linked to whole-body gene expression. A total of 3,644 genes were associated with behavioural specialization which is ten times more than associated with age and 50 times more than associated with fertility. Brood carers were characterized by an upregulation of three Vitellogenin (Vg) genes, one of which, Vg-like A, was the most differentially expressed gene that was recently shown experimentally to control the switch from brood to worker care. The expression of Conventional Vg was unlinked to behavioural specialization, age or fertility, which contrasts to studies on bees and some ants. Diversity in Vg/Vg-like copy number and expression bias across ants supports subfunctionalization of Vg genes and indicates that some regulatory mechanisms of division of labour diverged in different ant lineages. Simulations revealed that our experimental dissociation of co-varying factors reduced transcriptomic noise, suggesting that confounding factors could potentially explain inconsistencies across transcriptomic studies of behavioural specialization in ants. Thus, our study reveals that worker gene expression is mainly linked to the worker's function for the colony and provides novel insights into the evolution of sociality in ants.