Social evolution. Genomic signatures of evolutionary transitions from solitary to group living

Eusocial Transition in Blattodea: Transposable Elements and Shifts of Gene Expression

(1) Unravelling the molecular basis underlying major evolutionary transitions can shed light on how complex phenotypes arise. The evolution of eusociality, a major evolutionary transition, has been demonstrated to be accompanied by enhanced gene regulation. Numerous pieces of evidence suggest the major impact of transposon insertion on gene regulation and its role in adaptive evolution. Transposons have been shown to be play a role in gene duplication involved in the eusocial transition in termites. However, evidence of the molecular basis underlying the eusocial transition in Blattodea remains scarce. Could transposons have facilitated the eusocial transition in termites through shifts of gene expression? (2) Using available cockroach and termite genomes and transcriptomes, we investigated if transposons insert more frequently in genes with differential expression in queens and workers and if those genes could be linked to specific functions essential for eusocial transition. (3) The insertion rate of transposons differs among differentially expressed genes and displays opposite trends between termites and cockroaches. The functions of termite transposon-rich queen- and worker-biased genes are related to reproduction and ageing and behaviour and gene expression, respectively. (4) Our study provides further evidence on the role of transposons in the evolution of eusociality, potentially through shifts in gene expression.

Social complexity and brain evolution: insights from ant neuroarchitecture and genomics

Brain evolution is hypothesized to be driven by requirements to adaptively respond to environmental cues and social signals. Diverse models describe how sociality may have influenced eusocial insect-brain evolution, but specific impacts of social organization and other selective forces on brain architecture have been difficult to distinguish. Here, we evaluate predictions derived from and/or inferences made by models of social organization concerning the effects of individual and collective behavior on brain size, structure, and function using results of neuroanatomical and genomic studies. In contrast to the predictions of some models, we find that worker brains in socially complex species have great behavioral and cognitive capacity. We also find that colony size, the evolution of worker physical castes, and task specialization affect brain size and mosaicism, supporting the idea that sensory, processing and motor requirements for behavioral performance select for adaptive allometries of functionally specialized brain centers. We review available transcriptomic and comparative genomic studies seeking to elucidate the molecular pathways functionally associated with social life and the genetic changes that occurred during the evolution of social complexity. We discuss ways forward, using comparative neuroanatomy, transcriptomics, and comparative genomics, to distinguish among multiple alternative explanations for the relationship between the evolution of neural systems and social complexity.

Differential Gene Expression Correlates with Behavioural Polymorphism during Collective Behaviour in Cockroaches

Simple Summary

It is currently well accepted that animals differ from one another in their behaviour and tendency to perform actions, a property we refer to as animal personality. In group-living animals, variation in animal personality can be important to determine group survival, as it determines how individuals interact with each other and with their environment. However, we have little knowledge of the proximal mechanisms underlying personality, particularly in group-living organisms. Here, we investigate the relationship between gene expression and two behavioural types (bold and shy) in a gregarious species: the American cockroach. Our results show that bold individuals have upregulated genes with functions associated with sensory activity (phototaxis and odour detection) and aggressive/dominant behaviour, and suggest that social context can modulate gene expression related to bold/shy characteristics. This work could help identify genes important in the earliest stages of group living and social life, and provides a first step toward establishing cockroaches as a focal group for the study of the evolution of sociality.

Abstract

Consistent inter-individual variation in the propensity to perform different tasks (animal personality) can contribute significantly to the success of group-living organisms. The distribution of different personalities in a group influences collective actions and therefore how these organisms interact with their environment. However, we have little understanding of the proximate mechanisms underlying animal personality in animal groups, and research on this theme has often been biased towards organisms with advanced social systems. The goal of this study is to investigate the mechanistic basis for personality variation during collective behaviour in a species with rudimentary societies: the American cockroach. We thus use an approach which combines experimental classification of individuals into behavioural phenotypes (‘bold’ and ‘shy’ individuals) with comparative gene expression. Our analyses reveal differences in gene expression between behavioural phenotypes and suggest that social context may modulate gene expression related to bold/shy characteristics. We also discuss how cockroaches could be a valuable model for the study of genetic mechanisms underlying the early steps in the evolution of social behaviour and social complexity. This study provides a first step towards a better understanding of the molecular mechanisms associated with differences in boldness and behavioural plasticity in these organisms.

Genome of the bee Holcopasites calliopsidis—a species showing the common apid trait of brood parasitism

Brood parasites represent a substantial but often poorly studied fraction of the wider diversity of bees. Brood parasitic bees complete their life cycles by infiltrating the nests of solitary host bees thereby enabling their offspring to exploit the food provisions intended for the host’s offspring. Here, we present the draft assembly of the bee Holcopasites calliopsidis, the first brood parasitic species to be the subject of detailed genomic analysis. Consistent with previous findings on the genomic signatures of parasitism more broadly, we find that H. calliopsidis has the smallest genome currently known among bees (179 Mb). This small genome does not appear to be the result of purging of repetitive DNA, with some indications of novel repetitive elements which may show signs of recent expansion. Nor does H. calliopsidis demonstrate any apparent net loss of genic content in comparison with nonparasitic species, though many individual gene families do show significant contractions. Although the basis of the small genome size of this species remains unclear, the identification of over 12,000 putative genes—with functional annotation for nearly 10,000 of these—is an important step in investigating the genomic basis of brood parasitism and provides a valuable dataset to be compared against new genomes that remain to be sequenced.

Behavioral genetics and genomics: Mendel's peas, mice, and bees

The question of the heritability of behavior has been of long fascination to scientists and the broader public. It is now widely accepted that most behavioral variation has a genetic component, although the degree of genetic influence differs widely across behaviors. Starting with Mendel's remarkable discovery of "inheritance factors," it has become increasingly clear that specific genetic variants that influence behavior can be identified. This goal is not without its challenges: Unlike pea morphology, most natural behavioral variation has a complex genetic architecture. However, we can now apply powerful genome-wide approaches to connect variation in DNA to variation in behavior as well as analyses of behaviorally related variation in brain gene expression, which together have provided insights into both the genetic mechanisms underlying behavior and the dynamic relationship between genes and behavior, respectively, in a wide range of species and for a diversity of behaviors. Here, we focus on two systems to illustrate both of these approaches: the genetic basis of burrowing in deer mice and transcriptomic analyses of division of labor in honey bees. Finally, we discuss the troubled relationship between the field of behavioral genetics and eugenics, which reminds us that we must be cautious about how we discuss and contextualize the connections between genes and behavior, especially in humans.

Extreme reproductive skew at the dawn of sociality is consistent with inclusive fitness theory but problematic for routes to eusociality

To understand the earliest stages of social evolution, we need to identify species that are undergoing the initial steps into sociality. Amphylaeus morosus is the only unambiguously known social species in the bee family Colletidae and represents an independent origin of sociality within the Apoidea. This allows us to investigate the selective factors promoting the transition from solitary to social nesting. Using genome-wide SNP genotyping, we infer robust pedigree relationships to identify maternity of brood and intracolony relatedness for colonies at the end of the reproductive season. We show that A. morosus forms both matrifilial and full-sibling colonies, both involving complete or almost complete monopolization over reproduction. In social colonies, the reproductive primary was also the primary forager with the secondary female remaining in the nest, presumably as a guard. Social nesting provided significant protection against parasitism and increased brood survivorship in general. We show that secondary females gain large indirect fitness benefits from defensive outcomes, enough to satisfy the conditions of inclusive fitness theory, despite an over-production of males in social colonies. These results suggest an avenue to sociality that involves high relatedness and, very surprisingly, extreme reproductive skew in its earliest stages and raises important questions about the evolutionary steps in pathways to eusociality.

Annotation and Analysis of 3902 Odorant Receptor Protein Sequences from 21 Insect Species Provide Insights into the Evolution of Odorant Receptor Gene Families in Solitary and Social Insects

The gene family of insect olfactory receptors (ORs) has expanded greatly over the course of evolution. ORs enable insects to detect volatile chemicals and therefore play an important role in social interactions, enemy and prey recognition, and foraging. The sequences of several thousand ORs are known, but their specific function or their ligands have only been identified for very few of them. To advance the functional characterization of ORs, we have assembled, curated, and aligned the sequences of 3902 ORs from 21 insect species, which we provide as an annotated online resource. Using functionally characterized proteins from the fly Drosophila melanogaster, the mosquito Anopheles gambiae and the ant Harpegnathos saltator, we identified amino acid positions that best predict response to ligands. We examined the conservation of these predicted relevant residues in all OR subfamilies; the results showed that the subfamilies that expanded strongly in social insects had a high degree of conservation in their binding sites. This suggests that the ORs of social insect families are typically finely tuned and exhibit sensitivity to very similar odorants. Our novel approach provides a powerful tool to exploit functional information from a limited number of genes to study the functional evolution of large gene families.

Gene expression and chromatin conformation differs between worker bees performing different tasks

BACKGROUND

Revealing the effect of transcriptomic regulation on behavioral differences is a fundamental goal in biology, but the relationship between gene regulatory networks and individual behavior differences remains largely unknown. Honey bees are considered as good models for studying the mechanisms underlying gene expression changes and behavioral differences since they exhibit strong and obvious differences in tasks between individuals. The cis-regulatory regions usually contain the binding sites of diverse transcription factor (TFs) influencing bee behavior. Thus, the identification of cis-regulatory elements in the brains across different behavioral states is important for understanding how genomic and transcriptomic variations affect different tasks in honeybees.

METHODS

In this study, we employed transcriptome and genome-wide chromatin accessibility assays to analyze brain tissues of honey bees in different behavioral states for examining the relationship between individual behavior differences and brain gene expression changes. We also used the obtained open chromatin regions to identify cis-motifs associating differentially expressed TFs and genes in order to reveal the transcriptional regulatory mechanism related to the different tasks.

RESULTS

We identified genetic regulatory modules regulating different tasks that contained key TFs (CTCF, Trl and schlank) associated with open chromatin regions enriched for DNA sequence motifs belonging to the family of the corresponding TFs. The most prominent transcriptomic changes, which correlated with chromatin accessibility modifications within their proximal promoter regions, occurred in nervous system development, and were associated with behavior switch.

CONCLUSIONS

Our results revealed the regulatory landscape among three behavioral stages in honeybees and identified interactive molecular networks regulating different tasks. These results provide a comprehensive insight into behavioral differences of honeybees, which offers reference for future study.

Termite sociogenomics: evolution and regulation of caste-specific expressed genes

Termite genomes have been sequenced in at least five species from four different families. Genome-based transcriptome analyses have identified large numbers of protein-coding genes with caste-specific expression patterns. These genes include those involved in caste-specific morphologies and roles, for example high fecundity and longevity in reproductives. Some caste-specific expressed genes belong to multi-gene families, and their genetic architecture and expression profiles indicate they have evolved via tandem gene duplication. Candidate regulatory mechanisms of caste-specific expression include epigenetic regulation (e.g. histone modification and non-coding RNA) and diversification of transcription factors and cis-regulatory elements. We review current knowledge in the area of termite sociogenomics, focussing on the evolution and regulation of caste-specific expressed genes, and discuss future research directions.

Convergent Loss of Prothoracicotropic Hormone, A Canonical Regulator of Development, in Social Bee Evolution

The evolution of insect sociality has repeatedly involved changes in developmental events and their timing. Here, we propose the hypothesis that loss of a canonical regulator of moulting and metamorphosis, prothoracicotropic hormone (PTTH), and its receptor, Torso, is associated with the evolution of sociality in bees. Specifically, we posit that the increasing importance of social influences on early developmental timing in social bees has led to their decreased reliance on PTTH, which connects developmental timing with abiotic cues in solitary insects. At present, the evidence to support this hypothesis includes the absence of genes encoding PTTH and Torso from all fully-sequenced social bee genomes and its presence in all available genomes of solitary bees. Based on the bee phylogeny, the most parsimonious reconstruction of evolutionary events is that this hormone and its receptor have been lost multiple times, across independently social bee lineages. These gene losses shed light on possible molecular and cellular mechanisms that are associated with the evolution of social behavior in bees. We outline the available evidence for our hypothesis, and then contextualize it in light of what is known about developmental cues in social and solitary bees, and the multiple precedences of major developmental changes in social insects.

Burrowing deeper to uncover the genetic architecture of behavioral evolution

Allele-specific expression analysis of hybrid mice provides new insights into the genetic substrates of behavioral evolution. As a complement to QTL mapping, this approach, described by Hu et al. in this issue of Cell Reports, holds promise for identifying causative regulatory loci that influence species-specific behavior.

cis-Regulatory changes in locomotor genes are associated with the evolution of burrowing behavior

How evolution modifies complex, innate behaviors is largely unknown. Divergence in many morphological traits, and some behaviors, is linked to cis-regulatory changes in gene expression. Given this, we compare brain gene expression of two interfertile sister species of Peromyscus mice that show large and heritable differences in burrowing behavior. Species-level differential expression and allele-specific expression in F₁ hybrids indicate a preponderance of cis-regulatory divergence, including many genes whose cis-regulation is affected by burrowing behavior. Genes related to locomotor coordination show the strongest signals of lineage-specific selection on burrowing-induced cis-regulatory changes. Furthermore, genetic markers closest to these candidate genes associate with variation in burrow shape in a genetic cross, suggesting an enrichment for loci affecting burrowing behavior near these candidate locomotor genes. Our results provide insight into how cis-regulated gene expression can depend on behavioral context and how this dynamic regulatory divergence between species may contribute to behavioral evolution.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Convergent Loss of Chemoreceptors across Independent Origins of Slave-Making in Ants

The evolution of an obligate parasitic lifestyle often leads to the reduction of morphological and physiological traits, which may be accompanied by loss of genes and functions. Slave-making ants are social parasites that exploit the work force of closely related ant species for social behaviors such as brood care and foraging. Recent divergence between these social parasites and their hosts enables comparative studies of gene family evolution. We sequenced the genomes of eight ant species, representing three independent origins of ant slavery. During the evolution of eusociality, chemoreceptor genes multiplied due to the importance of chemical communication in insect societies. We investigated the evolutionary fate of these chemoreceptors and found that slave-making ant genomes harbored only half as many gustatory receptors as their hosts', potentially mirroring the outsourcing of foraging tasks to host workers. In addition, parasites had fewer odorant receptors and their loss shows striking patterns of convergence across independent origins of parasitism, in particular in orthologs often implicated in sociality like the 9-exon odorant receptors. These convergent losses represent a rare case of convergent molecular evolution at the level of individual genes. Thus, evolution can operate in a way that is both repeatable and reversible when independent ant lineages lose important social traits during the transition to a parasitic lifestyle.

OUP accepted manuscript

Many species have separate haploid and diploid phases. Theory predicts that each phase should experience the effects of evolutionary forces (like selection) differently. In the haploid phase, all fitness-affecting alleles are exposed to selection, whereas in the diploid phase, those same alleles can be masked by homologous alleles. This predicts that selection acting on genes expressed in haploids should be more effective than diploid-biased genes. Unfortunately, in arrhenotokous species, this prediction can be confounded with the effects of sex-specific expression, as haploids are usually reproductive males. Theory posits that, when accounting for ploidal- and sex-specific expression, selection should be equally efficient on haploid- and diploid-biased genes relative to constitutive genes. Here, we used a multiomic approach in honey bees to quantify the evolutionary rates of haploid-biased genes and test the relative effects of sexual- and haploid-expression on molecular evolution. We found that 16% of the honey bee’s protein-coding genome is highly expressed in haploid tissue. When accounting for ploidy and sex, haploid- and diploid-biased genes evolve at a lower rate than expected, indicating that they experience strong negative selection. However, the rate of molecular evolution of haploid-biased genes was higher than diploid-based genes. Genes associated with sperm storage are a clear exception to this trend with evidence of strong positive selection. Our results provide an important empirical test of theory outlining how selection acts on genes expressed in arrhenotokous species. We propose the haploid life history stage affects genome-wide patterns of diversity and divergence because of both sexual and haploid selection.

Complementary evolution of coding and noncoding sequence underlies mammalian hairlessness

Body hair is a defining mammalian characteristic, but several mammals, such as whales, naked mole-rats, and humans, have notably less hair. To find the genetic basis of reduced hair quantity, we used our evolutionary-rates-based method, RERconverge, to identify coding and noncoding sequences that evolve at significantly different rates in so-called hairless mammals compared to hairy mammals. Using RERconverge, we performed a genome-wide scan over 62 mammal species using 19,149 genes and 343,598 conserved noncoding regions. In addition to detecting known and potential novel hair-related genes, we also discovered hundreds of putative hair-related regulatory elements. Computational investigation revealed that genes and their associated noncoding regions show different evolutionary patterns and influence different aspects of hair growth and development. Many genes under accelerated evolution are associated with the structure of the hair shaft itself, while evolutionary rate shifts in noncoding regions also included the dermal papilla and matrix regions of the hair follicle that contribute to hair growth and cycling. Genes that were top ranked for coding sequence acceleration included known hair and skin genes KRT2, KRT35, PKP1, and PTPRM that surprisingly showed no signals of evolutionary rate shifts in nearby noncoding regions. Conversely, accelerated noncoding regions are most strongly enriched near regulatory hair-related genes and microRNAs, such as mir205, ELF3, and FOXC1, that themselves do not show rate shifts in their protein-coding sequences. Such dichotomy highlights the interplay between the evolution of protein sequence and regulatory sequence to contribute to the emergence of a convergent phenotype.

DNA methylation across the tree of life, from micro to macro-organism

DNA methylation is a process in which methyl (CH3) groups are added to the DNA molecule. The DNA segment does not change in the sequence, but DNA methylation could alter the action of DNA. Different enzymes like DNA methyltransferases (DNMTs) take part in methylation of cytosine/adenine nucleosides in DNA. In prokaryotes, DNA methylation is performed to prevent the attack of phage and also plays a role in the chromosome replication and repair. In fungi, DNA methylation is studied to see the transcriptional changes, as in insects, the DNA methylation is not that well-known, it plays a different role like other organisms. In mammals, the DNA methylation is related to different types of cancers and plays the most important role in the placental development and abnormal DNA methylation connected with diseases like cancer, autoimmune diseases, and rheumatoid arthritis.

Active genic machinery for epigenetic RNA modifications in bees

Epitranscriptomics is an emerging field of investigation dedicated to the study of post-transcriptional RNA modifications. RNA methylations regulate RNA metabolism and processing, including changes in response to environmental cues. Although RNA modifications are conserved from bacteria to eukaryotes, there is little evidence of an epitranscriptomic pathway in insects. Here we identified genes related to RNA m⁶ A (N6-methyladenine) and m⁵ C (5-methylcytosine) methylation machinery in seven bee genomes (Apis mellifera, Melipona quadrifasciata, Frieseomelitta varia, Eufriesea mexicana, Bombus terrestris, Megachile rotundata and Dufourea novaeangliae). In A. mellifera, we validated the expression of methyltransferase genes and found that the global levels of m⁶ A and m⁵ C measured in the fat body and brain of adult workers differ significantly. Also, m⁶ A levels were differed significantly mainly between the fourth larval instar of queens and workers. Moreover, we found a conserved m⁵ C site in the honeybee 28S rRNA. Taken together, we confirm the existence of epitranscriptomic machinery acting in bees and open avenues for future investigations on RNA epigenetics in a wide spectrum of hymenopteran species.

The Birth-and-Death Evolution of Cytochrome P450 Genes in Bees

The birth-and-death model of multigene family evolution describes how gene families evolve and diversify through duplication and deletion. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Extensive studies of insect P450 genes have demonstrated their role in insecticide resistance. Bees are thought to experience toxin exposure through their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some species. Here, we describe the repertoire of P450 genes in the orchid bee Euglossa dilemma. Male orchid bees form perfume bouquets used in courtship displays by collecting volatile compounds, resulting in exposure to compounds known to be toxic. In addition, we conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We find that social behavior and resin collection are not correlated with the repertoire of P450 present in a bee species. However, our analyses revealed that P450 clades can be classified as stable and unstable, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures and exhibit signals of adaptive evolution. This work highlights the complexity of multigene family evolution, revealing that multiple factors contribute to the diversification, stability, and dynamics of this gene family. Furthermore, we provide a resource for future detailed studies investigating the function of different P450s in economically important bee species.

Hymenoptera Genome Database: new genomes and annotation datasets for improved go enrichment and orthologue analyses

We report an update of the Hymenoptera Genome Database (HGD; http://HymenopteraGenome.org), a genomic database of hymenopteran insect species. The number of species represented in HGD has nearly tripled, with fifty-eight hymenopteran species, including twenty bees, twenty-three ants, eleven wasps and four sawflies. With a reorganized website, HGD continues to provide the HymenopteraMine genomic data mining warehouse and JBrowse/Apollo genome browsers integrated with BLAST. We have computed Gene Ontology (GO) annotations for all species, greatly enhancing the GO annotation data gathered from UniProt with more than a ten-fold increase in the number of GO-annotated genes. We have also generated orthology datasets that encompass all HGD species and provide orthologue clusters for fourteen taxonomic groups. The new GO annotation and orthology data are available for searching in HymenopteraMine, and as bulk file downloads.

Sociality and parasite transmission

Parasites and their social hosts form many different relationships. But what kind of selection regimes are important? A look at the parameters that determine fitness of the two parties suggests that social hosts differ from solitary ones primarily in the structure of transmission pathways. Because transmission is, both, the physical encounter of a new host and infecting it, several different elements determine parasite transmission success. These include spatial distance, genetic distance, or the temporal and ecological niche overlaps. Combing these elements into a ‘generalized transmission distance’ that determines parasite fitness aids in the identification of the critical steps. For example, short-distance transmission to genetically similar hosts within the social group is the most frequent process under sociality. Therefore, spatio-genetical distances are the main driver of parasite fitness. Vice versa, the generalized distance identifies the critical host defences. In this case, host defences should be primarily selected to defend against the within-group spread of an infection, especially among closely related group members.

Evolution of sociability by artificial selection

There has been extensive research on the ecology and evolution of social life in animals that live in groups. Less attention, however, has been devoted to apparently solitary species, even though recent research indicates that they also possess complex social behaviors. To address this knowledge gap, we artificially selected on sociability, defined as the tendency to engage in nonaggressive activities with others, in fruit flies. Our goal was to quantify the factors that determine the level of sociability and the traits correlated with this feature. After 25 generations of selection, the high-sociability lineages showed sociability scores about 50% higher than did the low-sociability lineages. Experiments using the evolved lineages indicated that there were no differences in mating success between flies from the low and high lineages. Both males and females from the low lineages, however, were more aggressive than males and females from the high lineages. Finally, the evolved lineages maintained their sociability scores after 10 generations of relaxed selection, suggesting no costs to maintaining low and high sociability, at least under our settings. Sociability is a complex trait, which we currently assess through genomic work on the evolved lineages.

Comparative Genomics Reveals Recent Adaptive Evolution in Himalayan Giant Honeybee Apis laboriosa

The Himalayan giant honeybee, Apis laboriosa, is the largest individual honeybee with major ecological and economic importance in high-latitude environments. However, our understanding of its environmental adaptations is circumscribed by the paucity of genomic data for this species. Here, we provide a draft genome of wild A. laboriosa, along with a comparison to its closely related species, Apis dorsata. The draft genome of A. laboriosa based on the de novo assembly is 226.1 Mbp in length with a scaffold N50 size of 3.34 Mbp, a GC content of 32.2%, a repeat content of 6.86%, and a gene family number of 8,404. Comparative genomics analysis revealed that the genes in A. laboriosa genome have undergone stronger positive selection (2.5 times more genes) and more recent duplication/loss events (6.1 times more events) than those in the A. dorsata genome. Our study implies the potential molecular mechanisms underlying the high-altitude adaptation of A. laboriosa and will catalyze future comparative studies to understand the environmental adaptation of modern honeybees.

Expansion and Accelerated Evolution of 9-Exon Odorant Receptors in Polistes Paper Wasps

Independent origins of sociality in bees and ants are associated with independent expansions of particular odorant receptor (OR) gene subfamilies. In ants, one clade within the OR gene family, the 9-exon subfamily, has dramatically expanded. These receptors detect cuticular hydrocarbons (CHCs), key social signaling molecules in insects. It is unclear to what extent 9-exon OR subfamily expansion is associated with the independent evolution of sociality across Hymenoptera, warranting studies of taxa with independently derived social behavior. Here, we describe OR gene family evolution in the northern paper wasp, Polistes fuscatus, and compare it to four additional paper wasp species spanning ∼40 million years of evolutionary divergence. We find 200 putatively functional OR genes in P. fuscatus, matching predictions from neuroanatomy, and more than half of these are in the 9-exon subfamily. Most OR gene expansions are tandemly arrayed at orthologous loci in Polistes genomes, and microsynteny analysis shows species-specific gain and loss of 9-exon ORs within tandem arrays. There is evidence of episodic positive diversifying selection shaping ORs in expanded subfamilies. Values of omega (dN/dS) are higher among 9-exon ORs compared to other OR subfamilies. Within the Polistes OR gene tree, branches in the 9-exon OR clade experience relaxed negative (relaxed purifying) selection relative to other branches in the tree. Patterns of OR evolution within Polistes are consistent with 9-exon OR function in CHC perception by combinatorial coding, with both natural selection and neutral drift contributing to interspecies differences in gene copy number and sequence.

Expression of methyl farnesoate epoxidase (mfe) and juvenile hormone esterase (jhe) genes and their relation to social organization in the stingless bee Melipona interrupta (Hymenoptera: Apidae)

Social organization in highly eusocial bees relies upon two important processes: caste differentiation in female larvae, and age polyethism in adult workers. Juvenile Hormone (JH) is a key regulator of both processes. Here we investigated the expression of two genes involved in JH metabolism - mfe (biosynthesis) and jhe (degradation) - in the context of social organization in the stingless bee Melipona interrupta. We found evidence that the expression of mfe and jhe genes is related to changes in JH levels during late larval development, where caste determination occurs. Also, both mfe and jhe were upregulated when workers engage in intranidal tasks, but only jhe expression was downregulated at the transition from nursing to foraging activities. This relation is different than expected, considering recent reports of lower JH levels in foragers than nurses in the closely related species Melipona scutellaris. Our findings suggest that highly eusocial bees have different mechanisms to regulate JH and, thus, to maintain their level of social organization.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO₂ detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

The foraging gene as a modulator of division of labour in social insects

The social ants, bees, wasps, and termites include some of the most ecologically-successful groups of animal species. Their dominance in most terrestrial environments is attributed to their social lifestyle, which enable their colonies to exploit environmental resources with remarkable efficiency. One key attribute of social insect colonies is the division of labour that emerges among the sterile workers, which represent the majority of colony members. Studies of the mechanisms that drive division of labour systems across diverse social species have provided fundamental insights into the developmental, physiological, molecular, and genomic processes that regulate sociality, and the possible genetic routes that may have led to its evolution from a solitary ancestor. Here we specifically discuss the conserved role of the foraging gene, which encodes a cGMP-dependent protein kinase (PKG). Originally identified as a behaviourally polymorphic gene that drives alternative foraging strategies in the fruit fly Drosophila melanogaster, changes in foraging expression and kinase activity were later shown to play a key role in the division of labour in diverse social insect species as well. In particular, foraging appears to regulate worker transitions between behavioural tasks and specific behavioural traits associated with morphological castes. Although the specific neuroethological role of foraging in the insect brain remains mostly unknown, studies in genetically tractable insect species indicate that PKG signalling plays a conserved role in the neuronal plasticity of sensory, cognitive and motor functions, which underlie behaviours relevant to division of labour, including appetitive learning, aggression, stress response, phototaxis, and the response to pheromones.

Eusociality in snapping shrimps is associated with larger genomes and an accumulation of transposable elements

Despite progress uncovering the genomic underpinnings of sociality, much less is known about how social living affects the genome. In different insect lineages, for example, eusocial species show both positive and negative associations between genome size and structure, highlighting the dynamic nature of the genome. Here, we explore the relationship between sociality and genome architecture in Synalpheus snapping shrimps that exhibit multiple origins of eusociality and extreme interspecific variation in genome size. Our goal is to determine whether eusociality leads to an accumulation of repetitive elements and an increase in genome size, presumably due to reduced effective population sizes resulting from a reproductive division of labor, or whether an initial accumulation of repetitive elements leads to larger genomes and independently promotes the evolution of eusociality through adaptive evolution. Using phylogenetically informed analyses, we find that eusocial species have larger genomes with more transposable elements (TEs) and microsatellite repeats than noneusocial species. Interestingly, different TE subclasses contribute to the accumulation in different species. Phylogenetic path analysis testing alternative causal relationships between sociality and genome architecture is most consistent with the hypothesis that TEs modulate the relationship between sociality and genome architecture. Although eusociality appears to influence TE accumulation, ancestral state reconstruction suggests moderate TE abundances in ancestral species could have fueled the initial transitions to eusociality. Ultimately, we highlight a complex and dynamic relationship between genome and social evolution, demonstrating that sociality can influence the evolution of the genome, likely through changes in demography related to patterns of reproductive skew.

The role of epigenetics, particularly DNA methylation, in the evolution of caste in insect societies

Eusocial insects can be defined as those that live in colonies and have distinct queens and workers. For most species, queens and workers arise from a common genome, and so caste-specific developmental trajectories must arise from epigenetic processes. In this review, we examine the epigenetic mechanisms that may be involved in the regulation of caste dimorphism. Early work on honeybees suggested that DNA methylation plays a causal role in the divergent development of queen and worker castes. This view has now been challenged by studies that did not find consistent associations between methylation and caste in honeybees and other species. Evidence for the involvement of methylation in modulating behaviour of adult workers is also inconsistent. Thus, the functional significance of DNA methylation in social insects remains equivocal. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Conservation insights from wild bee genetic studies: Geographic differences, susceptibility to inbreeding, and signs of local adaptation

Conserving bees are critical both ecologically and economically. Genetic tools are valuable for monitoring these vital pollinators since tracking these small, fast-flying insects by traditional means is difficult. By surveying the current state of the literature, this review discusses how recent advances in landscape genetic and genomic research are elucidating how wild bees respond to anthropogenic threats. Current literature suggests that there may be geographic differences in the vulnerability of bee species to landscape changes. Populations of temperate bee species are becoming more isolated and more genetically depauperate as their landscape becomes more fragmented, but tropical bee species appear unaffected. These differences may be an artifact of historical differences in land-use, or it suggests that different management plans are needed for temperate and tropical bee species. Encouragingly, genetic studies on invasive bee species indicate that low levels of genetic diversity may not lead to rapid extinction in bees as once predicted. Additionally, next-generation sequencing has given researchers the power to identify potential genes under selection, which are likely critical to species' survival in their rapidly changing environment. While genetic studies provide insights into wild bee biology, more studies focusing on a greater phylogenetic and life-history breadth of species are needed. Therefore, caution should be taken when making broad conservation decisions based on the currently few species examined.

Reciprocity and alignment: quantifying coupling in dynamic interactions

Recent accounts of social cognition focus on how we do things together, suggesting that becoming aligned relies on a reciprocal exchange of information. The next step is to develop richer computational methods that quantify the degree of coupling and describe the nature of the information exchange. We put forward a definition of coupling, comparing it to related terminology and detail, available computational methods and the level of organization to which they pertain, presenting them as a hierarchy from weakest to richest forms of coupling. The rationale is that a temporally coherent link between two dynamical systems at the lowest level of organization sustains mutual adaptation and alignment at the highest level. Postulating that when we do things together, we do so dynamically over time and we argue that to determine and measure instances of true reciprocity in social exchanges is key. Along with this computationally rich definition of coupling, we present challenges for the field to be tackled by a diverse community working towards a dynamic account of social cognition.

Transcriptomic Signatures of Ageing Vary in Solitary and Social Forms of an Orchid Bee

Eusocial insect queens are remarkable in their ability to maximize both fecundity and longevity, thus escaping the typical trade-off between these two traits. Several mechanisms have been proposed to underlie the remolding of the trade-off, such as reshaping of the juvenile hormone (JH) pathway, or caste-specific susceptibility to oxidative stress. However, it remains a challenge to disentangle the molecular mechanisms underlying the remolding of the trade-off in eusocial insects from caste-specific physiological attributes that have subsequently arisen. The socially polymorphic orchid bee Euglossa viridissima represents an excellent model to address the role of sociality per se in longevity as it allows direct comparisons of solitary and social individuals within a common genetic background. We investigated gene expression and JH levels in young and old bees from both solitary and social nests. We found 902 genes to be differentially expressed with age in solitary females, including genes involved in oxidative stress, versus only 100 genes in social dominant females, and 13 genes in subordinate females. A weighted gene coexpression network analysis further highlights pathways related to ageing in this species, including the target of rapamycin pathway. Eleven genes involved in translation, apoptosis, and DNA repair show concurrent age-related expression changes in solitary but not in social females, representing potential differences based on social status. JH titers did not vary with age or social status. Our results represent an important step in understanding the proximate mechanisms underlying the remodeling of the fecundity/longevity trade-off that accompanies the evolutionary transition from solitary life to eusociality.

Insects Provide Unique Systems to Investigate How Early-Life Experience Alters the Brain and Behavior

Early-life experiences have strong and long-lasting consequences for behavior in a surprising diversity of animals. Determining which environmental inputs cause behavioral change, how this information becomes neurobiologically encoded, and the functional consequences of these changes remain fundamental puzzles relevant to diverse fields from evolutionary biology to the health sciences. Here we explore how insects provide unique opportunities for comparative study of developmental behavioral plasticity. Insects have sophisticated behavior and cognitive abilities, and they are frequently studied in their natural environments, which provides an ecological and adaptive perspective that is often more limited in lab-based vertebrate models. A range of cues, from relatively simple cues like temperature to complex social information, influence insect behavior. This variety provides experimentally tractable opportunities to study diverse neural plasticity mechanisms. Insects also have a wide range of neurodevelopmental trajectories while sharing many developmental plasticity mechanisms with vertebrates. In addition, some insects retain only subsets of their juvenile neuronal population in adulthood, narrowing the targets for detailed study of cellular plasticity mechanisms. Insects and vertebrates share many of the same knowledge gaps pertaining to developmental behavioral plasticity. Combined with the extensive study of insect behavior under natural conditions and their experimental tractability, insect systems may be uniquely qualified to address some of the biggest unanswered questions in this field.

Eusociality Shapes Convergent Patterns of Molecular Evolution across Mitochondrial Genomes of Snapping Shrimps

Eusociality is a highly conspicuous and ecologically impactful behavioral syndrome that has evolved independently across multiple animal lineages. So far, comparative genomic analyses of advanced sociality have been mostly limited to insects. Here, we study the only clade of animals known to exhibit eusociality in the marine realm-lineages of socially diverse snapping shrimps in the genus Synalpheus. To investigate the molecular impact of sociality, we assembled the mitochondrial genomes of eight Synalpheus species that represent three independent origins of eusociality and analyzed patterns of molecular evolution in protein-coding genes. Synonymous substitution rates are lower and potential signals of relaxed purifying selection are higher in eusocial relative to noneusocial taxa. Our results suggest that mitochondrial genome evolution was shaped by eusociality-linked traits-extended generation times and reduced effective population sizes that are hallmarks of advanced animal societies. This is the first direct evidence of eusociality impacting genome evolution in marine taxa. Our results also strongly support the idea that eusociality can shape genome evolution through profound changes in life history and demography.

Comparative transcriptomics highlights convergent evolution of energy metabolic pathways in group-living spiders

Although widely thought to be aggressive, solitary, and potentially cannibalistic, some spider species have evolved group-living behaviors. The distinct transition provides the framework to uncover group-living evolution. Here, we conducted a comparative transcriptomic study and examined patterns of molecular evolution in two independently evolved group-living spiders and twelve solitary species. We report that positively selected genes among group-living spider lineages are significantly enriched in nutrient metabolism and autophagy pathways. We also show that nutrient-related genes of group-living spiders convergently experience amino acid substitutions and accelerated relative evolutionary rates. These results indicate adaptive convergence of nutrient metabolism that may ensure energy supply in group-living spiders. The decelerated evolutionary rate of autophagy-related genes in group-living lineages is consistent with an increased constraint on energy homeostasis as would be required in a group-living environment. Together, the results show that energy metabolic pathways play an important role in the transition to group-living in spiders.

Genome Architecture Facilitates Phenotypic Plasticity in the Honeybee (Apis mellifera)

Phenotypic plasticity, the ability of an organism to alter its phenotype in response to an environmental cue, facilitates rapid adaptation to changing environments. Plastic changes in morphology and behavior are underpinned by widespread gene expression changes. However, it is unknown if, or how, genomes are structured to ensure these robust responses. Here, we use repression of honeybee worker ovaries as a model of plasticity. We show that the honeybee genome is structured with respect to plasticity; genes that respond to an environmental trigger are colocated in the honeybee genome in a series of gene clusters, many of which have been assembled in the last 80 My during the evolution of the Apidae. These clusters are marked by histone modifications that prefigure the gene expression changes that occur as the ovary activates, suggesting that these genomic regions are poised to respond plastically. That the linear sequence of the honeybee genome is organized to coordinate widespread gene expression changes in response to environmental influences and that the chromatin organization in these regions is prefigured to respond to these influences is perhaps unexpected and has implications for other examples of plasticity in physiology, evolution, and human disease.

Neurodevelopmental and transcriptomic effects of CRISPR/Cas9-induced somatic orco mutation in honey bees

In insects, odorant receptors facilitate olfactory communication and require the functionality of the highly conserved co-receptor gene orco. Genome editing studies in a few species of ants and moths have revealed that orco can also have a neurodevelopmental function, in addition to its canonical role in adult olfaction, discovered first in Drosophila melanogaster. To extend this analysis, we determined whether orco mutations also affect the development of the adult brain of the honey bee Apis mellifera, an important model system for social behavior and chemical communication. We used CRISPR/Cas9 to knock out orco and examined anatomical and molecular consequences. To increase efficiency, we coupled embryo microinjection with a laboratory egg collection and in vitro rearing system. This new workflow advances genomic engineering technologies in honey bees by overcoming restrictions associated with field studies. We used Sanger sequencing to quickly select individuals with complete orco knockout for neuroanatomical analyses and later validated and described the mutations with amplicon sequencing. Mutant bees had significantly fewer glomeruli, smaller total volume of all the glomeruli, and higher mean individual glomerulus volume in the antennal lobe compared to wild-type controls. RNA-Sequencing revealed that orco knockout also caused differential expression of hundreds of genes in the antenna, including genes related to neural development and genes encoding odorant receptors. The expression of other types of chemoreceptor genes was generally unaffected, reflecting specificity of CRISPR activity in this study. These results suggest that neurodevelopmental effects of orco are related to specific insect life histories.

Systematic detection of brain protein-coding genes under positive selection during primate evolution and their roles in cognition

The human brain differs from that of other primates, but the genetic basis of these differences remains unclear. We investigated the evolutionary pressures acting on almost all human protein-coding genes (N = 11,667; 1:1 orthologs in primates) based on their divergence from those of early hominins, such as Neanderthals, and non-human primates. We confirm that genes encoding brain-related proteins are among the most strongly conserved protein-coding genes in the human genome. Combining our evolutionary pressure metrics for the protein-coding genome with recent data sets, we found that this conservation applied to genes functionally associated with the synapse and expressed in brain structures such as the prefrontal cortex and the cerebellum. Conversely, several genes presenting signatures commonly associated with positive selection appear as causing brain diseases or conditions, such as micro/macrocephaly, Joubert syndrome, dyslexia, and autism. Among those, a number of DNA damage response genes associated with microcephaly in humans such as BRCA1, NHEJ1, TOP3A, and RNF168 show strong signs of positive selection and might have played a role in human brain size expansion during primate evolution. We also showed that cerebellum granule neurons express a set of genes also presenting signatures of positive selection and that may have contributed to the emergence of fine motor skills and social cognition in humans. This resource is available online and can be used to estimate evolutionary constraints acting on a set of genes and to explore their relative contributions to human traits.

Sociality sculpts similar patterns of molecular evolution in two independently evolved lineages of eusocial bees

While it is well known that the genome can affect social behavior, recent models posit that social lifestyles can, in turn, influence genome evolution. Here, we perform the most phylogenetically comprehensive comparative analysis of 16 bee genomes to date: incorporating two published and four new carpenter bee genomes (Apidae: Xylocopinae) for a first-ever genomic comparison with a monophyletic clade containing solitary through advanced eusocial taxa. We find that eusocial lineages have undergone more gene family expansions, feature more signatures of positive selection, and have higher counts of taxonomically restricted genes than solitary and weakly social lineages. Transcriptomic data reveal that caste-affiliated genes are deeply-conserved; gene regulatory and functional elements are more closely tied to social phenotype than phylogenetic lineage; and regulatory complexity increases steadily with social complexity. Overall, our study provides robust empirical evidence that social evolution can act as a major and surprisingly consistent driver of macroevolutionary genomic change.

Gene expression and epigenetics reveal species-specific mechanisms acting upon common molecular pathways in the evolution of task division in bees

A striking feature of advanced insect societies is the existence of workers that forgo reproduction. Two broad types of workers exist in eusocial bees: nurses who care for their young siblings and the queen, and foragers who guard the nest and forage for food. Comparisons between these two worker subcastes have been performed in honeybees, but data from other bees are scarce. To understand whether similar molecular mechanisms are involved in nurse-forager differences across distinct species, we compared gene expression and DNA methylation profiles between nurses and foragers of the buff-tailed bumblebee Bombus terrestris and the stingless bee Tetragonisca angustula. These datasets were then compared to previous findings from honeybees. Our analyses revealed that although the expression pattern of genes is often species-specific, many of the biological processes and molecular pathways involved are common. Moreover, the correlation between gene expression and DNA methylation was dependent on the nucleotide context, and non-CG methylation appeared to be a relevant factor in the behavioral changes of the workers. In summary, task specialization in worker bees is characterized by a plastic and mosaic molecular pattern, with species-specific mechanisms acting upon broad common pathways across species.

Honey as a Functional Food for Apis mellifera

Although nectar is consumed, primarily as a supplemental food, by a broad range of insects spanning at least five orders, it is processed and stored by only a small number of species, most of which are bees and wasps in the superfamily Apoidea. Within this group, Apis mellifera has evolved remarkable adaptations facilitating nectar processing and storage; in doing so, this species utilizes the end product, honey, for diverse functions with few if any equivalents in other phytophagous insects. Honey and its phytochemical constituents, some of which likely derive from propolis, have functional significance in protecting honey bees against microbial pathogens, toxins, and cold stress, as well as in regulating development and adult longevity. The distinctive properties of A. mellifera honey appear to have arisen in multiple ways, including genome modification; partnerships with microbial symbionts; and evolution of specialized behaviors, including foraging for substances other than nectar. That honey making by A. mellifera involves incorporation of exogenous material other than nectar, as well as endogenous products such as antimicrobial peptides and royal jelly, suggests that regarding honey as little more than a source of carbohydrates for bees is a concept in need of revision.

Transposable Elements and the Evolution of Insects

Insects are major contributors to our understanding of the interaction between transposable elements (TEs) and their hosts, owing to seminal discoveries, as well as to the growing number of sequenced insect genomes and population genomics and functional studies. Insect TE landscapes are highly variable both within and across insect orders, although phylogenetic relatedness appears to correlate with similarity in insect TE content. This correlation is unlikely to be solely due to inheritance of TEs from shared ancestors and may partly reflect preferential horizontal transfer of TEs between closely related species. The influence of insect traits on TE landscapes, however, remains unclear. Recent findings indicate that, in addition to being involved in insect adaptations and aging, TEs are seemingly at the cornerstone of insect antiviral immunity. Thus, TEs are emerging as essential insect symbionts that may have deleterious or beneficial consequences on their hosts, depending on context.

Gene Tree Estimation Error with Ultraconserved Elements: An Empirical Study on Pseudapis Bees

Summarizing individual gene trees to species phylogenies using two-step coalescent methods is now a standard strategy in the field of phylogenomics. However, practical implementations of summary methods suffer from gene tree estimation error, which is caused by various biological and analytical factors. Greatly understudied is the choice of gene tree inference method and downstream effects on species tree estimation for empirical data sets. To better understand the impact of this method choice on gene and species tree accuracy, we compare gene trees estimated through four widely used programs under different model-selection criteria: PhyloBayes, MrBayes, IQ-Tree, and RAxML. We study their performance in the phylogenomic framework of $>$800 ultraconserved elements from the bee subfamily Nomiinae (Halictidae). Our taxon sampling focuses on the genus Pseudapis, a distinct lineage with diverse morphological features, but contentious morphology-based taxonomic classifications and no molecular phylogenetic guidance. We approximate topological accuracy of gene trees by assessing their ability to recover two uncontroversial, monophyletic groups, and compare branch lengths of individual trees using the stemminess metric (the relative length of internal branches). We further examine different strategies of removing uninformative loci and the collapsing of weakly supported nodes into polytomies. We then summarize gene trees with ASTRAL and compare resulting species phylogenies, including comparisons to concatenation-based estimates. Gene trees obtained with the reversible jump model search in MrBayes were most concordant on average and all Bayesian methods yielded gene trees with better stemminess values. The only gene tree estimation approach whose ASTRAL summary trees consistently produced the most likely correct topology, however, was IQ-Tree with automated model designation (ModelFinder program). We discuss these findings and provide practical advice on gene tree estimation for summary methods. Lastly, we establish the first phylogeny-informed classification for Pseudapis s. l. and map the distribution of distinct morphological features of the group. [ASTRAL; Bees; concordance; gene tree estimation error; IQ-Tree; MrBayes, Nomiinae; PhyloBayes; RAxML; phylogenomics; stemminess].

Individual differences in honey bee behavior enabled by plasticity in brain gene regulatory networks

Understanding the regulatory architecture of phenotypic variation is a fundamental goal in biology, but connections between gene regulatory network (GRN) activity and individual differences in behavior are poorly understood. We characterized the molecular basis of behavioral plasticity in queenless honey bee (Apis mellifera) colonies, where individuals engage in both reproductive and non-reproductive behaviors. Using high-throughput behavioral tracking, we discovered these colonies contain a continuum of phenotypes, with some individuals specialized for either egg-laying or foraging and 'generalists' that perform both. Brain gene expression and chromatin accessibility profiles were correlated with behavioral variation, with generalists intermediate in behavior and molecular profiles. Models of brain GRNs constructed for individuals revealed that transcription factor (TF) activity was highly predictive of behavior, and behavior-associated regulatory regions had more TF motifs. These results provide new insights into the important role played by brain GRN plasticity in the regulation of behavior, with implications for social evolution.

No obvious transcriptome-wide signature of indirect selection in termites

The evolution of sterile helper castes in social insects implies selection on genes that underlie variation in this nonreproductive phenotype. These focal genes confer no direct fitness and are presumed to evolve through indirect fitness effects on the helper's reproducing relatives. This separation of a gene's phenotypic effect on one caste and its fitness effect on another suggests that genes for this and other forms of reproductive altruism are buffered from selection and will thus evolve closer to the neutral rate than genes directly selected for selfish reproduction. We test this hypothesis by comparing the strength of selection at loci associated in their expression with reproductive versus sterile castes in termites. Specifically, we gather caste-biased gene expression data from four termite transcriptomes and measure the global dN/dS ratio across gene sets and phylogenetic lineages. We find that the majority of examined orthologous gene groups show patterns of nucleotide substitution that are consistent with strong purifying selection and display little evidence for distinct signatures of direct versus indirect selection in reproductive and sterile castes. For one particular species (Reticulitermes flavipes), the strength of purifying selection is relaxed in a reproductive nymph-biased gene set, which opposes the nearly neutral idea. In other species, the synonymous rate (dS) alone was often found to be the highest in the sterile worker caste, suggesting a more subtle signature of indirect selection or an altogether different relationship between caste-biased expression and rates of molecular evolution.