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

Pervasive relaxed selection in termite genomes

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

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.

A postmeiotically bifurcated roadmap of honeybee spermatogenesis marked by phylogenetically restricted genes

Haploid males of hymenopteran species produce gametes through an abortive meiosis I followed by meiosis II that can either be symmetric or asymmetric in different species. Thus, one spermatocyte could give rise to two spermatids with either equal or unequal amounts of cytoplasm. It is currently unknown what molecular features accompany these postmeiotic sperm cells especially in species with asymmetric meiosis II such as bees. Here we present testis single-cell RNA sequencing datasets from the honeybee (Apis mellifera) drones of 3 and 14 days after emergence (3d and 14d). We show that, while 3d testes exhibit active, ongoing spermatogenesis, 14d testes only have late-stage spermatids. We identify a postmeiotic bifurcation in the transcriptional roadmap during spermatogenesis, with cells progressing toward the annotated spermatids (SPT) and small spermatids (sSPT), respectively. Despite an overall similarity in their transcriptomic profiles, sSPTs express the fewest genes and the least RNA content among all the sperm cell types. Intriguingly, sSPTs exhibit a relatively high expression level for Hymenoptera-restricted genes and a high mutation load, suggesting that the special meiosis II during spermatogenesis in the honeybee is accompanied by phylogenetically young gene activities.

Hierarchical network structure as the source of hierarchical dynamics (power-law frequency spectra) in living and non-living systems: How state-trait continua (body plans, personalities) emerge from first principles in biophysics

Living systems are hierarchical control systems that display a small world network structure. In such structures, many smaller clusters are nested within fewer larger ones, producing a fractal-like structure with a 'power-law' cluster size distribution (a mereology). Just like their structure, the dynamics of living systems shows fractal-like qualities: the timeseries of inner message passing and overt behavior contain high frequencies or 'states' (treble) that are nested within lower frequencies or 'traits' (bass), producing a power-law frequency spectrum that is known as a 'state-trait continuum' in the behavioral sciences. Here, we argue that the power-law dynamics of living systems results from their power-law network structure: organisms 'vertically encode' the deep spatiotemporal structure of their (anticipated) environments, to the effect that many small clusters near the base of the hierarchy produce high frequency signal changes and fewer larger clusters at its top produce ultra-low frequencies. Such ultra-low frequencies exert a tonic regulatory pressure that produces morphological as well as behavioral traits (i.e., body plans and personalities). Nested-modular structure causes higher frequencies to be embedded within lower frequencies, producing a power-law state-trait continuum. At the heart of such dynamics lies the need for efficient energy dissipation through networks of coupled oscillators, which also governs the dynamics of non-living systems (e.q., earthquakes, stock market fluctuations). Since hierarchical structure produces hierarchical dynamics, the development and collapse of hierarchical structure (e.g., during maturation and disease) should leave specific traces in system dynamics (shifts in lower frequencies, i.e. morphological and behavioral traits) that may serve as early warning signs to system failure. The applications of this idea range from (bio)physics and phylogenesis to ontogenesis and clinical medicine.

Evolution: How sweat bees gained and lost eusociality

Eusocial insects divide labor between reproductive and non-reproductive individuals. The molecular mechanisms underlying the evolution of these castes have remained mysterious. A comparative genomic study of sweat bees points to a familiar factor as a regulator of behavioral specialization: juvenile hormone.

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.

Social insect transcriptomics and the molecular basis of caste diversity

Studies of gene expression provide fundamentally important information on the molecular mechanisms underlying variation in phenotype. Recent technological advances have allowed for the robust study of gene expression through analysis of whole transcriptomes. Here, we review current advances in social insect transcriptomics and discuss their implications in understanding phenotypic diversity. Recent transcriptomic studies provide detailed inventories of the genes involved in producing distinct phenotypes in social species. These investigations have identified key genes and networks involved in producing distinct social insect castes. Nevertheless, questions concerning the evolution of gene expression patterns remain. We suggest a path forward for studying gene expression in future studies of biological systems.

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.

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.

Sexual dimorphism as a facilitator of worker caste evolution in ants

Ant societies are primarily composed of females, whereby labor is divided into reproductive, or queen, and non-reproductive, or worker, castes. Workers and reproductive queens can differ greatly in behavior, longevity, physiology, and morphology, but queen-worker differences are usually modest relative to the differences in males. Males are short-lived, typically do not provide the colony with labor, often look like a different species, and only occur seasonally. It is these differences that have historically led to their neglect in social insect research, but also why they may facilitate novel phenotypic variation - by increasing the phenotypic variability that is available for selection. In this study, worker variation in multivariate size-shape space paralleled male-queen variation. As worker variation increased within species, so did sexual variation. Across species in two independent genera, using head width as a proxy for body size, sexual size dimorphism correlated with worker polymorphism regardless of whether the ancestral condition was large or small worker/sexual dimorphism. Mounting molecular data support the hypothesis that queen-worker caste determination has co-opted many genes/pathways from sex determination. The molecular evidence, coupled with the observations from this study, leads to the hypothesis that sexual selection and selection on colony-level traits are non-independent, and that sexual dimorphism may even have facilitated the evolution of the distinct worker caste.

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.

Comparative transcriptomics reveals commonalities and differences in the genetic underpinnings of a floral dimorphism

Distyly, a floral dimorphism associated with heteromorphic self-incompatibility and controlled by the S-locus supergene, evolved independently multiple times. Comparative analyses of the first transcriptome atlas for the main distyly model, Primula veris, with other distylous species produced the following findings. A set of 53 constitutively expressed genes in P. veris did not include any of the housekeeping genes commonly used to normalize gene expression in qPCR experiments. The S-locus gene CYP^T acquired its role in controlling style elongation via a change in expression profile. Comparison of genes differentially expressed between floral morphs revealed that brassinosteroids and auxin are the main hormones controlling style elongation in P. veris and Fagopyrum esculentum, respectively. Furthermore, shared biochemical pathways might underlie the expression of distyly in the distantly related P. veris, F. esculentum and Turnera subulata, suggesting a degree of correspondence between evolutionary convergence at phenotypic and molecular levels. Finally, we provide the first evidence supporting the previously proposed hypothesis that distyly supergenes of distantly related species evolved via the recruitment of genes related to the phytochrome-interacting factor (PIF) signaling network. To conclude, this is the first study that discovered homologous genes involved in the control of distyly in distantly related taxa.

Gene transcriptional profiles in gonads of Bacillus taxa (Phasmida) with different cytological mechanisms of automictic parthenogenesis

The evolution of automixis - i.e., meiotic parthenogenesis - requires several features, including ploidy restoration after meiosis and maintenance of fertility. Characterizing the relative contribution of novel versus pre-existing genes and the similarities in their expression and sequence evolution is fundamental to understand the evolution of reproductive novelties. Here we identify gonads-biased genes in two Bacillus automictic stick-insects and compare their expression profile and sequence evolution with a bisexual congeneric species. The two parthenogens restore ploidy through different cytological mechanisms: in Bacillus atticus, nuclei derived from the first meiotic division fuse to restore a diploid egg nucleus, while in Bacillus rossius, diploidization occurs in some cells of the haploid blastula through anaphase restitution. Parthenogens' gonads transcriptional program is found to be largely assembled from genes that were already present before the establishment of automixis. The three species transcriptional profiles largely reflect their phyletic relationships, yet we identify a shared core of genes with gonad-biased patterns of expression in parthenogens which are either male gonads-biased in the sexual species or are not differentially expressed there. At the sequence level, just a handful of gonads-biased genes were inferred to have undergone instances of positive selection exclusively in the parthenogen species. This work is the first to explore the molecular underpinnings of automixis in a comparative framework: it delineates how reproductive novelties can be sustained by genes whose origin precedes the establishment of the novelty itself and shows that different meiotic mechanisms of reproduction can be associated with a shared molecular ground plan.

Genomic signatures of recent convergent transitions to social life in spiders

The transition from solitary to social life is a major phenotypic innovation, but its genetic underpinnings are largely unknown. To identify genomic changes associated with this transition, we compare the genomes of 22 spider species representing eight recent and independent origins of sociality. Hundreds of genes tend to experience shifts in selection during the repeated transition to social life. These genes are associated with several key functions, such as neurogenesis, behavior, and metabolism, and include genes that previously have been implicated in animal social behavior and human behavioral disorders. In addition, social species have elevated genome-wide rates of molecular evolution associated with relaxed selection caused by reduced effective population size. Altogether, our study provides unprecedented insights into the genomic signatures of social evolution and the specific genetic changes that repeatedly underpin the evolution of sociality. Our study also highlights the heretofore unappreciated potential of transcriptomics using ethanol-preserved specimens for comparative genomics and phylotranscriptomics.

Evolutionary endocrinology and the problem of Darwin's tangled bank

Like Darwin's tangled bank of biodiversity, the endocrine mechanisms that give rise to phenotypic diversity also exhibit nearly endless forms. This tangled bank of mechanistic diversity can prove problematic as we seek general principles on the role of endocrine mechanisms in phenotypic evolution. A key unresolved question is therefore: to what degree are specific endocrine mechanisms re-used to bring about replicated phenotypic evolution? Related areas of inquiry are booming in molecular ecology, but behavioral traits are underrepresented in this literature. Here, I leverage the rich comparative tradition in evolutionary endocrinology to evaluate whether and how certain mechanisms may be repeated hotspots of behavioral evolutionary change. At one extreme, mechanisms may be parallel, such that evolution repeatedly uses the same gene or pathway to arrive at multiple independent (or, convergent) origins of a particular behavioral trait. At the other extreme, the building blocks of behavior may be unique, such that outwardly similar phenotypes are generated via lineage-specific mechanisms. This review synthesizes existing case studies, phylogenetic analyses, and experimental evolutionary research on mechanistic parallelism in animal behavior. These examples show that the endocrine building blocks of behavior have some elements of parallelism across replicated evolutionary events. However, support for parallelism is variable among studies, at least some of which relates to the level of complexity at which we consider sameness (i.e. pathway vs. gene level). Moving forward, we need continued experimentation and better testing of neutral models to understand whether, how - and critically, why - mechanism A is used in one lineage and mechanism B is used in another. We also need continued growth of large-scale comparative analyses, especially those that can evaluate which endocrine parameters are more or less likely to undergo parallel evolution alongside specific behavioral traits. These efforts will ultimately deepen understanding of how and why hormone-mediated behaviors are constructed the way that they are.

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.

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.

Do the different life history strategies of ants and honeybees determine fat body morphology?

The separation of two sister groups such as ants and bees in the Cretaceous involved the development of distinctive characteristics to occupy separate ecological niches. From the point of view of biology and ecology, it is important to see how different life history strategies affect the physiology of these insects. The fat body is the most metabolically important tissue in the organism of each insect. Therefore, we conducted a comparative analysis of the morphological image of the subcuticular fat body in different localisation/segments in Formica (Seviformica) cinerea and Apis mellifera mellifera foragers, because of the similarity of their functions in colonies. We observed that the fat bodies of ants and bees were composed of the same cell types: trophocytes and oenocytes. However, in each of the segments, the fat body cells in ants were bigger and there were fewer of them in comparison with bees. The dorsal part of the fat body of ants had a bilayer structure, where the outer layer was formed by binucleated oenocytes. Binucleated oenocytes were also found in the inner layer near the heart and tracheole. In bees, the fat body was unilayered and the trophocytes and oenocytes were present side by side. The similarities and, in particular, the differences in the structure of the fat body are the adaptation of these sister groups to life in a diverse environment.

Foraging Payoffs Change With Group Size in Kin and Non-kin Groups of an Argyrodinae Kleptoparasitic Spider, Argyrodes miniaceus

Evolutionary transitions from solitary to group-living are ubiquitous in animal systems. While the fitness consequences of group size changes are often investigated, the long-standing debate on whether kinship is a prerequisite of sociality is still ongoing. In the current study, we used kleptoparasitic spiders Argyrodes miniaceus (subfamily Argyrodinae, Theridiidae) as a model system to assess the role of group size on the foraging payoffs of kin and non-kin groups. We set up laboratory-manipulated kin and non-kin foraging groups and used feeding occurrence and duration as proxies for foraging benefits and feeding latency and the number of host attacks as estimates of foraging costs. Compared to solitary individuals, feeding durations of successfully fed individuals in groups was not significantly different from that of solitary foragers in both kin and non-kin groups. The occurrences of feeding decreased significantly in group sizes two and above, in non-kin groups, and in group sizes three and above, in kin groups. In kin groups, groups size two had significantly shorter feeding latencies compared to other group sizes, even though feeding duration did not change systematically with group size. Similarly, the number of attacks from the hosts were highest in non-kin groups with more than two individuals and in kin groups with more than three individuals. The juxtaposition of kin and non-kin group showed that A. miniaceus enjoyed the highest foraging payoffs when being solitary or in small groups (group size two). However, host attacks appeared to hamper feeding occurrences in kin groups, which was not observed in non-kin groups. Our results contrast sharply with the feeding benefits of kinship recorded in kin-based groups of sub-social species present in related subfamilies in the Theridiidae.

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.

Time-course RNASeq of Camponotus floridanus forager and nurse ant brains indicate links between plasticity in the biological clock and behavioral division of labor

Background

Circadian clocks allow organisms to anticipate daily fluctuations in their environment by driving rhythms in physiology and behavior. Inter-organismal differences in daily rhythms, called chronotypes, exist and can shift with age. In ants, age, caste-related behavior and chronotype appear to be linked. Brood-tending nurse ants are usually younger individuals and show “around-the-clock” activity. With age or in the absence of brood, nurses transition into foraging ants that show daily rhythms in activity. Ants can adaptively shift between these behavioral castes and caste-associated chronotypes depending on social context. We investigated how changes in daily gene expression could be contributing to such behavioral plasticity in Camponotus floridanus carpenter ants by combining time-course behavioral assays and RNA-Sequencing of forager and nurse brains.

Results

We found that nurse brains have three times fewer 24 h oscillating genes than foragers. However, several hundred genes that oscillated every 24 h in forager brains showed robust 8 h oscillations in nurses, including the core clock genes Period and Shaggy. These differentially rhythmic genes consisted of several components of the circadian entrainment and output pathway, including genes said to be involved in regulating insect locomotory behavior. We also found that Vitellogenin, known to regulate division of labor in social insects, showed robust 24 h oscillations in nurse brains but not in foragers. Finally, we found significant overlap between genes differentially expressed between the two ant castes and genes that show ultradian rhythms in daily expression.

Conclusion

This study provides a first look at the chronobiological differences in gene expression between forager and nurse ant brains. This endeavor allowed us to identify a putative molecular mechanism underlying plastic timekeeping: several components of the ant circadian clock and its output can seemingly oscillate at different harmonics of the circadian rhythm. We propose that such chronobiological plasticity has evolved to allow for distinct regulatory networks that underlie behavioral castes, while supporting swift caste transitions in response to colony demands. Behavioral division of labor is common among social insects. The links between chronobiological and behavioral plasticity that we found in C. floridanus, thus, likely represent a more general phenomenon that warrants further investigation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08282-x.

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.

Segmentation of the subcuticular fat body in Apis mellifera females with different reproductive potentials

Evolution has created different castes of females in eusocial haplodiploids. The difference between them lies in their functions and vulnerability but above all in their reproductive potentials. Honeybee queens are highly fertile. On the other hand, the workers are facultatively sterile. However, rebel workers, i.e. workers that develop in a queenless colony, reproduce more often than normal workers. As a result, the fat body of these bees, which apart from acting as the energy reserve, is also the site of numerous metabolic processes, had to specialize in different functions perfected over millions of years of eusocial evolution. Assuming that the variety of functions manifests itself in the pleomorphic structure of the fat body cells, we predicted that also different parts of the fat body, e.g. from different segments of the abdomen, contain different sets of cells. Such differences could be expected between queens, rebels and normal workers, i.e. females with dramatically different reproductive potentials. We confirmed all these expectations. Although all bees had the same types of cells, their proportion and segmental character corresponded with the caste reproductive potential and physiological characteristics shaped in the evolutionary process. The females with an increased reproductive potential were characterized by the presence of oenocytes in the third tergite and high concentrations of compounds responsible for energy reserves, like glucose, glycogen and triglycerides. Queens had very large trophocytes, especially in the third tergite. Only in workers did we observe intercellular spaces in all the segments of the fat body, as well as high protein concentrations-especially in the sternite. As expected, the rebels combined many features of the queens and normal workers, what with other findings can help understand the ways that led to the origin of different castes in females of eusocial Hymenoptera.

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.

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.

A chromosome-level genome assembly of the parasitoid wasp, Cotesia glomerata (Hymenoptera: Braconidae)

Hymenopterans make up about 20% of all animal species, but most are poorly known and lack high-quality genomic resources. One group of important, yet under-studied hymenopterans, are parasitoid wasps in the family Braconidae. Among this under-studied group are braconid wasps in the genus Cotesia; a clade of ~1,000 species routinely used in studies of physiology, ecology, biological control, and genetics. However, our ability to understand these organisms has been hindered by a lack of genomic resources. We helped bridge this gap by generating a high-quality genome assembly for the parasitoid wasp, Cotesia glomerata (Braconidae; Microgastrinae). We generated this assembly using multiple sequencing technologies, including Oxford Nanopore, whole-genome shotgun sequencing, and 3-D chromatin contact information (Hi-C). Our assembly is one of the most contiguous, complete, and publicly available hymenopteran genomes, represented by 3,355 scaffolds with a scaffold N50 of ~28Mb and a BUSCO score of ~99%. Given the genome sizes found in closely related species, our genome assembly was ~50% larger than expected, which was apparently induced by runaway amplification of three types of repetitive elements: simple repeats, Long Terminal Repeats (LTRs), and Long Interspersed Nuclear Elements (LINEs). This assembly is another step forward for genomics across this hyper-diverse, yet understudied, order of insects. The assembled genomic data and metadata files are publicly available via Figshare (https://doi.org/10.6084/m9.figshare.13010549).

Disentangling the aging gene expression network of termite queens

Background

Most insects are relatively short-lived, with a maximum lifespan of a few weeks, like the aging model organism, the fruit-fly Drosophila melanogaster. By contrast, the queens of many social insects (termites, ants and some bees) can live from a few years to decades. This makes social insects promising models in aging research providing insights into how a long reproductive life can be achieved. Yet, aging studies on social insect reproductives are hampered by a lack of quantitative data on age-dependent survival and time series analyses that cover the whole lifespan of such long-lived individuals. We studied aging in queens of the drywood termite Cryptotermes secundus by determining survival probabilities over a period of 15 years and performed transcriptome analyses for queens of known age that covered their whole lifespan.

Results

The maximum lifespan of C. secundus queens was 13 years, with a median maximum longevity of 11.0 years. Time course and co-expression network analyses of gene expression patterns over time indicated a non-gradual aging pattern. It was characterized by networks of genes that became differentially expressed only late in life, namely after ten years, which associates well with the median maximum lifespan for queens. These old-age gene networks reflect processes of physiological upheaval. We detected strong signs of stress, decline, defense and repair at the transcriptional level of epigenetic control as well as at the post-transcriptional level with changes in transposable element activity and the proteostasis network. The latter depicts an upregulation of protein degradation, together with protein synthesis and protein folding, processes which are often down-regulated in old animals. The simultaneous upregulation of protein synthesis and autophagy is indicative of a stress-response mediated by the transcription factor cnc, a homolog of human nrf genes.

Conclusions

Our results show non-linear senescence with a rather sudden physiological upheaval at old-age. Most importantly, they point to a re-wiring in the proteostasis network and stress as part of the aging process of social insect queens, shortly before queens die.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07649-4.

Ageing and sociality: why, when and how does sociality change ageing patterns?

Individual lifespans vary tremendously between and also within species, but the proximate and ultimate causes of different ageing speeds are still not well understood. Sociality appears to be associated with the evolution of greater longevity and probably also with a larger plasticity of the shape and pace of ageing. For example, reproductives of several termites and ants reach lifespans that surpass those of their non-reproductive nestmates by one or two decades. In this issue, 15 papers explore the interrelations between sociality and individual longevity in both, group-living vertebrates and social insects. Here, we briefly give an overview of the contents of the various contributions, including theoretical and comparative studies, and we explore the similarities and dissimilarities in proximate mechanisms underlying ageing among taxa, with particular emphasis on nutrient-sensing pathways and, in insects, juvenile hormone. These studies point to an underestimated role of more downstream processes. We highlight the need for reliable transcriptomic markers of ageing and a comprehensive ageing theory of social animals, which includes the reproductive potential of workers, and considers the fact that social insect queens reach maturity only after a prolonged period of producing non-reproductive workers. This article is part of the theme issue 'Ageing and sociality: why, when and how does sociality change ageing patterns?'

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?'

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 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.

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.

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.

Brain microRNAs among social and solitary bees

Evolutionary transitions to a social lifestyle in insects are associated with lineage-specific changes in gene expression, but the key nodes that drive these regulatory changes are unknown. We examined the relationship between social organization and lineage-specific microRNAs (miRNAs). Genome scans across 12 bee species showed that miRNA copy-number is mostly conserved and not associated with sociality. However, deep sequencing of small RNAs in six bee species revealed a substantial proportion (20-35%) of detected miRNAs had lineage-specific expression in the brain, 24-72% of which did not have homologues in other species. Lineage-specific miRNAs disproportionately target lineage-specific genes, and have lower expression levels than shared miRNAs. The predicted targets of lineage-specific miRNAs are not enriched for genes with caste-biased expression or genes under positive selection in social species. Together, these results suggest that novel miRNAs may coevolve with novel genes, and thus contribute to lineage-specific patterns of evolution in bees, but do not appear to have significant influence on social evolution. Our analyses also support the hypothesis that many new miRNAs are purged by selection due to deleterious effects on mRNA targets, and suggest genome structure is not as influential in regulating bee miRNA evolution as has been shown for mammalian miRNAs.

Functional genomics of parental care of insects

Parental care was likely the first step most lineages made towards sociality. However, the molecular mechanisms that generate parental care are not broadly characterized. Insects are important as an evolutionary independent group from classic models of parental care, such as, house mice. They provide an opportunity to test the generality of our understanding. With this review, I survey the functional genomics of parental care of insects, summarize several recent advances in the broader framework for studying and understanding parental care, and finish with suggested priorities for further research. Although there are too few studies to draw definitive conclusions, I argue that natural selection appears to be rewiring existing gene networks to produce parental care, that the epigenetic mechanisms influencing parental care are not well understood, and, as an interesting early consensus, that genes strongly associated with carer/offspring interactions appear biased towards proteins that are secreted. I summarize the studies that have functionally validate candidate genes and highlight the increasing need to perform this work. I finish with arguments for both conceptual and practical changes moving forward. I argue that future work can increase the use of predictive frameworks, broaden its definition of conservation of mechanism to gene networks rather than single genes, and increase the use of more established comparative methods. I further highlight the practical considerations of standardizing analyses and reporting, increasing the sampling of both carers and offspring, better characterizing gene regulatory networks, better characterizing taxonomically restricted genes and any consistent role they have underpinning parental care, and using factorial designs to disentangle the influence of multiple variables on the expression of parental care.

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.

Developmental plasticity shapes social traits and selection in a facultatively eusocial bee

Developmental processes are an important source of phenotypic variation, but the extent to which this variation contributes to evolutionary change is unknown. We used integrative genomic analyses to explore the relationship between developmental and social plasticity in a bee species that can adopt either a social or solitary lifestyle. We find genes regulating this social flexibility also regulate development, and positive selection on these genes is influenced by their function during development. This suggests that developmental plasticity may influence the evolution of sociality. Our additional finding of genetic variants linked to differences in social behavior sheds light on how phenotypic variation derived from development may become encoded into the genome, and thus contribute to evolutionary change.

Developmental plasticity generates phenotypic variation, but how it contributes to evolutionary change is unclear. Phenotypes of individuals in caste-based (eusocial) societies are particularly sensitive to developmental processes, and the evolutionary origins of eusociality may be rooted in developmental plasticity of ancestral forms. We used an integrative genomics approach to evaluate the relationships among developmental plasticity, molecular evolution, and social behavior in a bee species (Megalopta genalis) that expresses flexible sociality, and thus provides a window into the factors that may have been important at the evolutionary origins of eusociality. We find that differences in social behavior are derived from genes that also regulate sex differentiation and metamorphosis. Positive selection on social traits is influenced by the function of these genes in development. We further identify evidence that social polyphenisms may become encoded in the genome via genetic changes in regulatory regions, specifically in transcription factor binding sites. Taken together, our results provide evidence that developmental plasticity provides the substrate for evolutionary novelty and shapes the selective landscape for molecular evolution in a major evolutionary innovation: Eusociality.

Comparative Genomics Identifies Putative Signatures of Sociality in Spiders

Comparative genomics has begun to elucidate the genomic basis of social life in insects, but insight into the genomic basis of spider sociality has lagged behind. To begin, to characterize genomic signatures associated with the evolution of social life in spiders, we performed one of the first spider comparative genomics studies including five solitary species and two social species, representing two independent origins of sociality in the genus Stegodyphus. We found that the two social spider species had a large expansion of gene families associated with transport and metabolic processes and an elevated genome-wide rate of molecular evolution compared with the five solitary spider species. Genes that were rapidly evolving in the two social species relative to the five solitary species were enriched for transport, behavior, and immune functions, whereas genes that were rapidly evolving in the solitary species were enriched for energy metabolism processes. Most rapidly evolving genes in the social species Stegodyphus dumicola were broadly expressed across four tissues and enriched for transport functions, but 12 rapidly evolving genes showed brain-specific expression and were enriched for social behavioral processes. Altogether, our study identifies putative genomic signatures and potential candidate genes associated with spider sociality. These results indicate that future spider comparative genomic studies, including broader sampling and additional independent origins of sociality, can further clarify the genomic causes and consequences of social life.

Studying convergent evolution to relate genotype to behavioral phenotype

Neuroscience has a long, rich history in embracing unusual animals for research. Over the past several decades, there has been a technology-driven bottleneck in the species used for neuroscience research. However, an oncoming wave of technologies applicable to many animals hold promise for enabling researchers to address challenging scientific questions that cannot be solved using traditional laboratory animals. Here, we discuss how leveraging the convergent evolution of physiological or behavioral phenotypes can empower research mapping genotype to phenotype interactions. We present two case studies using electric fish and poison frogs and discuss how comparative work can teach us about evolutionary constraint and flexibility at various levels of biological organization. We also offer advice on the potential and pitfalls of establishing novel model systems in neuroscience research. Finally, we end with a discussion on the use of charismatic animals in neuroscience research and their utility in public outreach. Overall, we argue that convergent evolution frameworks can help identify generalizable principles of neuroscience.

Re-thinking the social ladder approach for elucidating the evolution and molecular basis of insect societies

The evolution of large insect societies is a major evolutionary transition that occurred in the long-extinct ancestors of termites, ants, corbiculate bees, and vespid wasps. Researchers have long used 'social ladder thinking': assuming progressive stepwise phenotypic evolution and asserting that extant species with simple societies (e.g. some halictid bees) represent the ancestors of species with complex societies, and thus provide insight into general early steps of eusocial evolution. We discuss how this is inconsistent with data and modern evolutionary 'tree thinking'. Phylogenetic comparative methods with broad sampling provide the best means to make rigorous inferences about ancestral traits and evolutionary transitions that occurred within each lineage, and to determine whether consistent phenotypic and genomic changes occurred across independent lineages.