Gene duplication and the evolution of phenotypic diversity in insect societies

Anciently duplicated genes continuously recruited to heart expression in vertebrate evolution are associated with heart chamber increase

Although gene/genome duplications in the early stage of vertebrates have been thought to provide major resources of raw genetic materials for evolutionary innovations, it is unclear whether they continuously contribute to the evolution of morphological complexity during the course of vertebrate evolution, such as the evolution from two heart chambers (fishes) to four heart chambers (mammals and birds). We addressed this issue by our heart RNA-Seq experiments combined with published data, using 13 vertebrates and one invertebrate (sea squirt, as an outgroup). Our evolutionary transcriptome analysis showed that number of ancient paralogous genes expressed in heart tends to increase with the increase of heart chamber number along the vertebrate phylogeny, in spite that most of them were duplicated at the time near to the origin of vertebrates or even more ancient. Moreover, those paralogs expressed in heart exert considerably different functions from heart-expressed singletons: the former are functionally enriched in cardiac muscle and muscle contraction-related categories, whereas the latter play more basic functions of energy generation like aerobic respiration. These findings together support the notion that recruiting anciently paralogous genes that are expressed in heart is associated with the increase of chamber number in vertebrate evolution.

Models for the retention of duplicate genes and their biological underpinnings

Gene content in genomes changes through several different processes, with gene duplication being an important contributor to such changes. Gene duplication occurs over a range of scales from individual genes to whole genomes, and the dynamics of this process can be context dependent. Still, there are rules by which genes are retained or lost from genomes after duplication, and probabilistic modeling has enabled characterization of these rules, including their context-dependence. Here, we describe the biology and corresponding mathematical models that are used to understand duplicate gene retention and its contribution to the set of biochemical functions encoded in a genome.

Lectin diversity and their positive roles in WSSV replication through regulation of calreticulin expression and inhibiting ALFs expression

In biological evolution, gene duplication (GD) generates new genes to facilitate new functions. C-type lectins (CTLs) in crayfish have been extended by GD to expand their family members. In this study, four CTL genes generated by GD were identified from Procambarus clarkii (PcLec1-4). Among these four genes, PcLec1 can also generate new isoforms with different numbers of tandem repeats through DNA slip mispairing. PcLec1-4 was widely expressed in multiple tissues. The expression levels of PcLec1-4 were upregulated in the intestine of P. clarkii upon white spot syndrome virus (WSSV) challenge at multiple time points. Further analysis indicated that GATA transcription factor regulated PcLec1-4 expression. RNA interference and recombinant PcLec1-4 protein injection experiments suggested that PcLec1-4 promoted the expression of calreticulin (PcCRT) and negatively regulated the expression of antimicrobial peptides, thereby promoting WSSV replication. This study contributes to the understanding of the function of CTLs produced by GD during WSSV invasion in crustaceans.

Predicting evolutionary targets and parameters of gene deletion from expression data

Motivation

Gene deletion is traditionally thought of as a nonadaptive process that removes functional redundancy from genomes, such that it generally receives less attention than duplication in evolutionary turnover studies. Yet, mounting evidence suggests that deletion may promote adaptation via the "less-is-more" evolutionary hypothesis, as it often targets genes harboring unique sequences, expression profiles, and molecular functions. Hence, predicting the relative prevalence of redundant and unique functions among genes targeted by deletion, as well as the parameters underlying their evolution, can shed light on the role of gene deletion in adaptation.

Results

Here, we present CLOUDe, a suite of machine learning methods for predicting evolutionary targets of gene deletion events from expression data. Specifically, CLOUDe models expression evolution as an Ornstein-Uhlenbeck process, and uses multi-layer neural network, extreme gradient boosting, random forest, and support vector machine architectures to predict whether deleted genes are "redundant" or "unique", as well as several parameters underlying their evolution. We show that CLOUDe boasts high power and accuracy in differentiating between classes, and high accuracy and precision in estimating evolutionary parameters, with optimal performance achieved by its neural network architecture. Application of CLOUDe to empirical data from Drosophila suggests that deletion primarily targets genes with unique functions, with further analysis showing these functions to be enriched for protein deubiquitination. Thus, CLOUDe represents a key advance in learning about the role of gene deletion in functional evolution and adaptation.

Availability and implementation

CLOUDe is freely available on GitHub (https://github.com/anddssan/CLOUDe).

Four ciliate-specific expansion events occurred during actin gene family evolution of eukaryotes

Actin gene family is a divergent and ancient eukaryotic cellular cytoskeletal gene family, and participates in many essential cellular processes. Ciliated protists offer us an excellent opportunity to investigate gene family evolution, since their gene families evolved faster in ciliates than in other eukaryotes. Nonetheless, actin gene family is well studied in few model ciliate species but little is known about its evolutionary patterns in ciliates. Here, we analyzed the evolutionary pattern of eukaryotic actin gene family based on genomes/transcriptomes of 36 species covering ten ciliate classes, as well as those of nine non-ciliate eukaryotic species. Results showed: (1) Except for conventional actins and actin-related proteins (Arps) shared by various eukaryotes, at least four ciliate-specific subfamilies occurred during evolution of actin gene family. Expansions of Act2 and ArpC were supposed to have happen in the ciliate common ancestor, while expansions of ActI and ActII may have occurred in the ancestor of Armophorea, Muranotrichea, and Spirotrichea. (2) The number of actin isoforms varied greatly among ciliate species. Environmental adaptability, whole genome duplication (WGD) or segmental duplication events, distinct spatial and temporal patterns of expression might play driving forces for the increasement of isoform numbers. (3) The 'birth and death' model of evolution could explain the evolution of actin gene family in ciliates. And actin genes have been generally under strong negative selection to maintain protein structures and physiological functions. Collectively, we provided meaningful information for understanding the evolution of eukaryotic actin gene family.

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.

Expression of fatty acid synthase genes and their role in development and arboviral infection of Aedes aegypti

BACKGROUND

Fatty acids are the building blocks of complex lipids essential for living organisms. In mosquitoes, fatty acids are involved in cell membrane production, energy conservation and expenditure, innate immunity, development and reproduction. Fatty acids are synthesized by a multifunctional enzyme complex called fatty acid synthase (FAS). Several paralogues of FAS were found in the Aedes aegypti mosquito. However, the molecular characteristics and expression of some of these paralogues have not been investigated.

METHODS

Genome assemblies of Ae. aegypti were analyzed, and orthologues of human FAS was identified. Phylogenetic analysis and in silico molecular characterization were performed to identify the functional domains of the Ae. aegypti FAS (AaFAS). Quantitative analysis and loss-of-function experiments were performed to determine the significance of different AaFAS transcripts in various stages of development, expression following different diets and the impact of AaFAS on dengue virus, serotype 2 (DENV2) infection and transmission.

RESULTS

We identified seven putative FAS genes in the Ae. aegypti genome assembly, based on nucleotide similarity to the FAS proteins (tBLASTn) of humans, other mosquitoes and invertebrates. Bioinformatics and molecular analyses suggested that only five of the AaFAS genes produce mRNA and therefore represent complete gene models. Expression levels of AaFAS varied among developmental stages and between male and female Ae. aegypti. Quantitative analyses revealed that expression of AaFAS1, the putative orthologue of the human FAS, was highest in adult females. Transient knockdown (KD) of AaFAS1 did not induce a complete compensation by other AaFAS genes but limited DENV2 infection of Aag2 cells in culture and the midgut of the mosquito.

CONCLUSION

AaFAS1 is the predominant AaFAS in adult mosquitoes. It has the highest amino acid similarity to human FAS and contains all enzymatic domains typical of human FAS. AaFAS1 also facilitated DENV2 replication in both cell culture and in mosquito midguts. Our data suggest that AaFAS1 may play a role in transmission of dengue viruses and could represent a target for intervention strategies.

The Molecular Properties and Roles of Pannier in Harmonia axyridis's Metamorphosis and Melanin Synthesis

The GATA transcription factor Pannier is identified as the major regulatory gene in color pattern formation in the Asian multi-colored ladybird beetle (Harmonia axyridis). however, the mechanisms of Pannier in regulating melanin synthesis and development in H. axyridis remain elusive. In this study, we identified and characterized Pannier in H. axyridis (HaPnr) and showed it to have two alternative spliced variants named HaPnr-α and HaPnr-β. Analyses of developmental stage expression patterns revealed that HaPnr, HaPnr-α and HaPnr-β were constitutively expressed throughout all developmental stages. To examine the role of HaPnr in H. axyridis development, RNA interference was performed in late larvae (the fourth instar) and early pupae (the first day of pupa stage). The transcript levels of HaPnr were effectively suppressed after the injection of double-stranded RNA of HaPnr (dsHaPnr). The fourth instar larvae injected with dsHaPnr reduced the pupation rates to only 61.50%, compared with 88.5% in the dsGFP-injected group. The un-pupated larvae gradually died after 1 week, and visually unaffected pupae emerged into abnormal adults with malformed hind wings and melanin absent from the cuticle. These abnormal adults gradually died 10 days after eclosion. However, when early pupae were injected with dsHaPnr, the normal eclosion rate was achieved at 88.41% on day 6 after the injection. In addition, these successful eclosion adults also showed an absence of melanin in the cuticle, but they could mate normally and have normal fecundity as compared with the control. We further demonstrated that the suppression of HaPnr-α or HaPnr-β individually did not affect the pupation and eclosion process. The suppression of HaPnr-α expression resulted in elytra melanin decreasing in both the conspicua and the succinea subgroup in H. axyridis. Even though the suppression of HaPnr-β expression only affected the melanin synthesis in the succinea subgroup, it significantly prolonged the time taken for melanin synthesis to occur in the conspicua subgroup in H. axyridis. These results indicate that HaPnr plays an essential role in insect development, especially during their metamorphosis, and also support our hypothesis that HaPnr could regulate melanin synthesis in H. axyridis under the combined action with its two splicing variants, HaPnr-α and HaPnr-β.

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.

Genomic and transcriptomic analyses of the subterranean termite Reticulitermes speratus: Gene duplication facilitates social evolution

Gene duplication is a major source of evolutionary innovation and is associated with the increases in biological complexity and adaptive radiation. Termites are model social organisms characterized by a sophisticated caste system. We analyzed the genome of the Japanese subterranean termite, an ecologically and economically important insect acting as a destructive pest. The analyses revealed the significance of gene duplication in social evolution. Gene duplication associated with caste-biased gene expression was prevalent in the termite genome. Many of the duplicated genes were related to social functions, such as chemical communication, social immunity, and defense, and they were often expressed in caste-specific organs. We propose that gene duplication facilitates social evolution through regulatory diversification leading to caste-biased expression and functional specialization.

Termites are model social organisms characterized by a polyphenic caste system. Subterranean termites (Rhinotermitidae) are ecologically and economically important species, including acting as destructive pests. Rhinotermitidae occupies an important evolutionary position within the clade representing a transitional taxon between the higher (Termitidae) and lower (other families) termites. Here, we report the genome, transcriptome, and methylome of the Japanese subterranean termite Reticulitermes speratus. Our analyses highlight the significance of gene duplication in social evolution in this termite. Gene duplication associated with caste-biased gene expression was prevalent in the R. speratus genome. The duplicated genes comprised diverse categories related to social functions, including lipocalins (chemical communication), cellulases (wood digestion and social interaction), lysozymes (social immunity), geranylgeranyl diphosphate synthase (social defense), and a novel class of termite lineage–specific genes with unknown functions. Paralogous genes were often observed in tandem in the genome, but their expression patterns were highly variable, exhibiting caste biases. Some of the assayed duplicated genes were expressed in caste-specific organs, such as the accessory glands of the queen ovary and the frontal glands of soldier heads. We propose that gene duplication facilitates social evolution through regulatory diversification, leading to caste-biased expression and subfunctionalization and/or neofunctionalization conferring caste-specialized functions.

BLAST from the Past: Impacts of Evolving Approaches on Studies of Evolution by Gene Duplication

In 1970, Susumu Ohno hypothesized that gene duplication was a major reservoir of adaptive innovation. However, it was not until over two decades later that DNA sequencing studies uncovered the ubiquity of gene duplication across all domains of life, highlighting its global importance in the evolution of phenotypic complexity and species diversification. Today, it seems that there are no limits to the study of evolution by gene duplication, as it has rapidly coevolved with numerous experimental and computational advances in genomics. In this perspective, we examine word stem usage in PubMed abstracts to infer how evolving discoveries and technologies have shaped the landscape of studying evolution by gene duplication, leading to a more refined understanding of its role in the emergence of novel phenotypes.

Learning Retention Mechanisms and Evolutionary Parameters of Duplicate Genes from Their Expression Data

Learning about the roles that duplicate genes play in the origins of novel phenotypes requires an understanding of how their functions evolve. A previous method for achieving this goal, CDROM, employs gene expression distances as proxies for functional divergence and then classifies the evolutionary mechanisms retaining duplicate genes from comparisons of these distances in a decision tree framework. However, CDROM does not account for stochastic shifts in gene expression or leverage advances in contemporary statistical learning for performing classification, nor is it capable of predicting the parameters driving duplicate gene evolution. Thus, here we develop CLOUD, a multi-layer neural network built on a model of gene expression evolution that can both classify duplicate gene retention mechanisms and predict their underlying evolutionary parameters. We show that not only is the CLOUD classifier substantially more powerful and accurate than CDROM, but that it also yields accurate parameter predictions, enabling a better understanding of the specific forces driving the evolution and long-term retention of duplicate genes. Further, application of the CLOUD classifier and predictor to empirical data from Drosophila recapitulates many previous findings about gene duplication in this lineage, showing that new functions often emerge rapidly and asymmetrically in younger duplicate gene copies, and that functional divergence is driven by strong natural selection. Hence, CLOUD represents a major advancement in classifying retention mechanisms and predicting evolutionary parameters of duplicate genes, thereby highlighting the utility of incorporating sophisticated statistical learning techniques to address long-standing questions about evolution after gene duplication.

Gene Duplication in the Honeybee: Patterns of DNA Methylation, Gene Expression, and Genomic Environment

Gene duplication serves a critical role in evolutionary adaptation by providing genetic raw material to the genome. The evolution of duplicated genes may be influenced by epigenetic processes such as DNA methylation, which affects gene function in some taxa. However, the manner in which DNA methylation affects duplicated genes is not well understood. We studied duplicated genes in the honeybee Apis mellifera, an insect with a highly sophisticated social structure, to investigate whether DNA methylation was associated with gene duplication and genic evolution. We found that levels of gene body methylation were significantly lower in duplicate genes than in single-copy genes, implicating a possible role of DNA methylation in postduplication gene maintenance. Additionally, we discovered associations of gene body methylation with the location, length, and time since divergence of paralogous genes. We also found that divergence in DNA methylation was associated with divergence in gene expression in paralogs, although the relationship was not completely consistent with a direct link between DNA methylation and gene expression. Overall, our results provide further insight into genic methylation and how its association with duplicate genes might facilitate evolutionary processes and adaptation.

R2R3-MYB genes control petal pigmentation patterning in Clarkia gracilis ssp. sonomensis (Onagraceae)

Petal pigmentation patterning is widespread in flowering plants. The genetics of these pattern elements has been of great interest for understanding the evolution of phenotypic diversification. Here, we investigate the genetic changes responsible for the evolution of an unpigmented petal element on a colored background. We used transcriptome analysis, gene expression assays, cosegregation in F₂ plants and functional tests to identify the gene(s) involved in petal coloration in Clarkia gracilis ssp. sonomensis. We identified an R2R3-MYB transcription factor (CgsMYB12) responsible for anthocyanin pigmentation of the basal region ('cup') in the petal of C. gracilis ssp. sonomensis. A functional mutation in CgsMYB12 creates a white cup on a pink petal background. Additionally, we found that two R2R3-MYB genes (CgsMYB6 and CgsMYB11) are also involved in petal background pigmentation. Each of these three R2R3-MYB genes exhibits a different spatiotemporal expression pattern. The functionality of these R2R3-MYB genes was confirmed through stable transformation of Arabidopsis. Distinct spatial patterns of R2R3-MYB expression have created the possibility that pigmentation in different sections of the petal can evolve independently. This finding suggests that recent gene duplication has been central to the evolution of petal pigmentation patterning in C. gracilis ssp. sonomensis.

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.

From phenotype to genotype: the precursor hypothesis predicts genetic influences that facilitate transitions in social behavior

Parental care is expected to be one of the key evolutionary precursors to advanced social behavior. This suggests that there could be common genetic underpinnings to both parental care and sociality. However, little is known of the genetics underlying care. Here, we suggest that ethological predictions of behavioral precursors to care along with a genetic toolkit for behavior provide testable hypotheses and a defined approach to investigating genetics of sociality. We call this the 'precursor hypothesis'.

Coevolution of Genome Architecture and Social Behavior

Although social behavior can have a strong genetic component, it can also result in selection on genome structure and function, thereby influencing the evolution of the genome itself. Here we explore the bidirectional links between social behavior and genome architecture by considering variation in social and/or mating behavior among populations (social polymorphisms) and across closely related species. We propose that social behavior can influence genome architecture via associated demographic changes due to social living. We establish guidelines to exploit emerging whole-genome sequences using analytical approaches that examine genome structure and function at different levels (regulatory vs structural variation) from the perspective of both molecular biology and population genetics in an ecological context.

Rapid functional divergence after small-scale gene duplication in grasses

BACKGROUND

Gene duplication has played an important role in the evolution and domestication of flowering plants. Yet little is known about how plant duplicate genes evolve and are retained over long timescales, particularly those arising from small-scale duplication (SSD) rather than whole-genome duplication (WGD) events.

RESULTS

We address this question in the Poaceae (grass) family by analyzing gene expression data from nine tissues of Brachypodium distachyon, Oryza sativa japonica (rice), and Sorghum bicolor (sorghum). Consistent with theoretical predictions, expression profiles of most grass genes are conserved after SSD, suggesting that functional conservation is the primary outcome of SSD in grasses. However, we also uncover support for widespread functional divergence, much of which occurs asymmetrically via the process of neofunctionalization. Moreover, neofunctionalization preferentially targets younger (child) duplicate gene copies, is associated with RNA-mediated duplication, and occurs quickly after duplication. Further analysis reveals that functional divergence of SSD-derived genes is positively correlated with both sequence divergence and tissue specificity in all three grass species, and particularly with anther expression in B. distachyon.

CONCLUSIONS

Our results suggest that SSD-derived grass genes often undergo rapid functional divergence that may be driven by natural selection on male-specific phenotypes. These observations are consistent with those in several animal species, suggesting that duplicate genes take similar evolutionary trajectories in plants and animals.

Invasive invertebrates associated with highly duplicated gene content

Invasion of alien species has led to serious problems, including the destruction of native ecosystems. In general, invasive species adapt to new environments rapidly, suggesting that they have high genetic diversity that can directly influence environmental adaptability. However, it is not known how genomic architecture causes genetic diversity that leads to invasiveness. Recent studies have showed that the proportion of duplicated genes (P_D ) in whole animal genomes correlate with environmental variability within a habitat. Here, we show that P_D and propagule size significantly explain the differences in species categories (invasive species, noninvasive species, and parasites). P_D correlated negatively with the propagule size. The residual values of regression of P_D on propagule size revealed that the invasive species had higher P_D values and larger propagule size than those of the noninvasive species, whereas the parasites had lower P_D values and smaller propagule size than those of others. There were no correlations between the invasive species and other genomic factors including the genome size, number of genes, and certain gene families. Our results suggest that the P_D values of a genome might be a potential genomic source causing genetic variations for adaptation to diverse environments. The results also showed that the invasiveness status of a species would be predicted by the residual values of regression of P_D on propagule size. Our innovative approach provides a measure to estimate the environmental adaptability of organisms based on genomic data.

Hormonal control and target genes of ftz-f1 expression in the honeybee Apis mellifera: a positive loop linking juvenile hormone, ftz-f1, and vitellogenin

Ftz-f1 is an orphan member of the nuclear hormone receptor superfamily. A 20-hydroxyecdysone pulse allows ftz-f1 gene expression, which then regulates the activity of downstream genes involved in major developmental progression events. In honeybees, the expression of genes like vitellogenin (vg), prophenoloxidase and juvenile hormone-esterase during late pharate-adult development is known to be hormonally controlled in both queens and workers by increasing juvenile hormone (JH) titres in the presence of declining levels of ecdysteroids. Since Ftz-f1 is known for mediating intracellular JH signalling, we hypothesized that ftz-f1 could mediate JH action during the pharate-adult development of honeybees, thus controlling the expression of these genes. Here, we show that ftz-f1 has caste-specific transcription profiles during this developmental period, with a peak coinciding with the increase in JH titre, and that its expression is upregulated by JH and downregulated by ecdysteroids. RNAi-mediated knock down of ftz-f1 showed that the expression of genes essential for adult development (e.g. vg and cuticular genes) depends on ftz-f1 expression. Finally, a double-repressor hypothesis-inspired vg gene knock-down experiment suggests the existence of a positive molecular loop between JH, ftz-f1 and vg.

Conserved Genes Underlie Phenotypic Plasticity in an Incipiently Social Bee

Despite a strong history of theoretical work on the mechanisms of social evolution, relatively little is known of the molecular genetic changes that accompany transitions from solitary to eusocial forms. Here, we provide the first genome of an incipiently social bee that shows both solitary and social colony organization in sympatry, the Australian carpenter bee Ceratina australensis. Through comparative analysis, we provide support for the role of conserved genes and cis-regulation of gene expression in the phenotypic plasticity observed in nest-sharing, a rudimentary form of sociality. Additionally, we find that these conserved genes are associated with caste differences in advanced eusocial species, suggesting these types of mechanisms could pave the molecular pathway from solitary to eusocial living. Genes associated with social nesting in this species show signatures of being deeply conserved, in contrast to previous studies in other bees showing novel and faster-evolving genes are associated with derived sociality. Our data provide support for the idea that the earliest social transitions are driven by changes in gene regulation of deeply conserved genes.