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

Genomic and transcriptomic analyses of a social hemipteran provide new insights into insect sociality

The origin of sociality represents one of the most important evolutionary transitions. Insect sociality evolved in some hemipteran aphids, which can produce soldiers and normal nymphs with distinct morphology and behaviour through parthenogenesis. The lack of genomic data resources has hindered the investigations into molecular mechanisms underlying their social evolution. Herein, we generated the first chromosomal-level genome of a social hemipteran (Pseudoregma bambucicola) with highly specialized soldiers and performed comparative genomic and transcriptomic analyses to elucidate the molecular signatures and regulatory mechanisms of caste differentiation. P. bambucicola has a larger known aphid genome of 582.2 Mb with an N50 length of 11.24 Mb, and about 99.6% of the assembly was anchored to six chromosomes with a scaffold N50 of 98.27 Mb. A total of 14,027 protein-coding genes were predicted and 37.33% of the assembly were identified as repeat sequences. The social evolution is accompanied by a variety of changes in genome organization, including expansion of gene families related to transcription factors, transposable elements, as well as species-specific expansions of certain sugar transporters and UGPases involved in carbohydrate metabolism. We also characterized large candidate gene sets linked to caste differentiation and found evidence of expression regulation and positive selection acting on energy metabolism and muscle structure, explaining the soldier-specific traits including morphological and behavioural specialization, developmental arrest and infertility. Overall, this study offers new insights into the molecular basis of social aphids and the evolution of insect sociality and also provides valuable data resources for further comparative and functional studies.

Repeated Shifts in Sociality Are Associated With Fine-tuning of Highly Conserved and Lineage-Specific Enhancers in a Socially Flexible Bee

Comparative genomic studies of social insects suggest that changes in gene regulation are associated with evolutionary transitions in social behavior, but the activity of predicted regulatory regions has not been tested empirically. We used self-transcribing active regulatory region sequencing, a high-throughput enhancer discovery tool, to identify and measure the activity of enhancers in the socially variable sweat bee, Lasioglossum albipes. We identified over 36,000 enhancers in the L. albipes genome from 3 social and 3 solitary populations. Many enhancers were identified in only a subset of L. albipes populations, revealing rapid divergence in regulatory regions within this species. Population-specific enhancers were often proximal to the same genes across populations, suggesting compensatory gains and losses of regulatory regions may preserve gene activity. We also identified 1,182 enhancers with significant differences in activity between social and solitary populations, some of which are conserved regulatory regions across species of bees. These results indicate that social trait variation in L. albipes is associated with the fine-tuning of ancient enhancers as well as lineage-specific regulatory changes. Combining enhancer activity with population genetic data revealed variants associated with differences in enhancer activity and identified a subset of differential enhancers with signatures of selection associated with social behavior. Together, these results provide the first empirical map of enhancers in a socially flexible bee and highlight links between cis-regulatory variation and the evolution of social behavior.

Dedicated developmental programing for group-supporting behaviors in eusocial honeybees

The evolutionary changes from solitary to eusocial living in vertebrates and invertebrates are associated with the diversification of social interactions and the development of queen and worker castes. Despite strong innate patterns, our understanding of the mechanisms manifesting these sophisticated behaviors is still rudimentary. Here, we show that doublesex (dsx) manifests group-supporting behaviors in the honeybee (Apis mellifera) worker caste. Computer-based individual behavioral tracking of worker bees with biallelic stop mutations in colonies revealed that the dsx gene is required for the rate and duration of group-supporting behavior that scales the relationship between bees and their work. General sensorimotor functions remained unaffected. Unexpectedly, unlike in other insects, the dsx gene is required for the neuronal wiring of the mushroom body in which the gene is spatially restricted expressed. Together, our study establishes dedicated programming for group-supporting behaviors and provides insight into the connection between development in the neuronal circuitry and behaviors regulating the formation of a eusocial society.

Reminiscences on the honeybee genome project and the rise of epigenetic concepts in insect science

The sequencing of the honeybee genome in 2006 was an important technological and logistic achievement experience. But what benefits have flown from the honeybee genome project? What does the annotated genomic assembly mean for the study of behavioural complexity and organismal function in honeybees? Here, I discuss several lines of research that have arisen from this project and highlight the rapidly expanding studies on insect epigenomics, emergent properties of royal jelly, the mechanism of nutritional control of development and the contribution of epigenomic regulation to the evolution of sociality. I also argue that the term 'insect epigenetics' needs to be carefully redefined to reflect the diversity of epigenomic toolkits in insects and the impact of lineage-specific innovations on organismal outcomes. The honeybee genome project helped pioneer advances in social insect molecular biology, and fuelled breakthrough research into the role of flexible epigenomic control systems in linking genotype to phenotype.

Evolution of odorant receptor repertoires across Hymenoptera is not linked to the evolution of eusociality

Communication is essential for social organisms. In eusocial insects, olfaction facilitates communication and recognition between nestmates. The study of certain model organisms has led to the hypothesis that odorant receptors are expanded in eusocial Hymenoptera. This has become a widely mentioned idea in the literature, albeit with conflicting reports, and has not been tested with a broad comparative analysis. Here we combined existing genomic and new neuroanatomical data, including from an approximately 100 Myr old fossil ant, across a phylogenetically broad sample of hymenopteran lineages. We find no evidence that variation in the size and evolutionary tempo of odorant receptor repertoires is related to eusociality. Post hoc exploration of our data hinted at loss of flight as a possible factor shaping some of the variation in OR repertoires in Hymenoptera. Nevertheless, our analyses revealed a complex pattern of evolutionary variation, and raise new questions about the ecological, behavioural and social factors that shape olfactory abilities.

No single rescue recipe: genome complexities modulate insect response to climate change

Declines in insect populations have gained formidable attention. Given their crucial role in the ecosystem, the causes of declining insect populations must be investigated. However, the insect clade has been associated with low extinction and high diversification rates. It is unlikely that insects underwent mass extinctions in the past. However, the pace of current climate change could make insect populations vulnerable to extinction. We propose genome size (GS) and transposable elements (TEs) to be rough estimates to assess extinction risk. Larger GS and/or proliferating TEs have been associated with adaptation in rapid climate change scenarios. We speculate that unstable, stressful environmental conditions are strongly associated with GS and TE expansion, which could be further correlated with adaptations. Alternately, stressful conditions trigger TE bursts that are not purged in smaller populations. GS and TE loads could be indicators of small effective populations in the wild, likely experiencing bottlenecks or drastic climatic perturbations, which calls for an urgent assessment of extinction risk.

Major changes in domain arrangements are associated with the evolution of termites

Domains as functional protein units and their rearrangements along the phylogeny can shed light on the functional changes of proteomes associated with the evolution of complex traits like eusociality. This complex trait is associated with sterile soldiers and workers, and long-lived, highly fecund reproductives. Unlike in Hymenoptera (ants, bees, and wasps), the evolution of eusociality within Blattodea, where termites evolved from within cockroaches, was accompanied by a reduction in proteome size, raising the question of whether functional novelty was achieved with existing rather than novel proteins. To address this, we investigated the role of domain rearrangements during the evolution of termite eusociality. Analysing domain rearrangements in the proteomes of three solitary cockroaches and five eusocial termites, we inferred more than 5,000 rearrangements over the phylogeny of Blattodea. The 90 novel domain arrangements that emerged at the origin of termites were enriched for several functions related to longevity, such as protein homeostasis, DNA repair, mitochondrial activity, and nutrient sensing. Many domain rearrangements were related to changes in developmental pathways, important for the emergence of novel castes. Along with the elaboration of social complexity, including permanently sterile workers and larger, foraging colonies, we found 110 further domain arrangements with functions related to protein glycosylation and ion transport. We found an enrichment of caste-biased expression and splicing within rearranged genes, highlighting their importance for the evolution of castes. Furthermore, we found increased levels of DNA methylation among rearranged compared to non-rearranged genes suggesting fundamental differences in their regulation. Our findings indicate the importance of domain rearrangements in the generation of functional novelty necessary for termite eusociality to evolve.

Identification of chemosensory genes in the stingless bee Tetragonisca fiebrigi

Reception of chemical information from the environment is crucial for insects' survival and reproduction. The chemosensory reception mainly occurs by the antennae and mouth parts of the insect, when the stimulus contacts the chemoreceptors located within the sensilla. Chemosensory receptor genes have been well-studied in some social hymenopterans such as ants, honeybees, and wasps. However, although stingless bees are the most representative group of eusocial bees, little is known about their odorant, gustatory, and ionotropic receptor genes. Here, we analyze the transcriptome of the proboscis and antennae of the stingless bee Tetragonisca fiebrigi. We identified and annotated 9 gustatory and 15 ionotropic receptors. Regarding the odorant receptors, we identified 204, and we were able to annotate 161 of them. In addition, we compared the chemosensory receptor genes of T. fiebrigi with those annotated for other species of Hymenoptera. We found that T. fiebrigi showed the largest number of odorant receptors compared with other bees. Genetic expansions were identified in the subfamilies 9-exon, which was also expanded in ants and paper wasps; in G02A, including receptors potentially mediating social behavior; and in GUnC, which has been related to pollen and nectar scent detection. Our study provides the first report of chemosensory receptor genes in T. fiebrigi and represents a resource for future molecular and physiological research in this and other stingless bee species.

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.

Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae)

Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.

The transcriptomic signature of adaptations associated with perfume collection in orchid bees

Secondary sexual traits can convey information on mate quality with the signal honesty maintained by the costly nature of trait expression. Mating signals are also often underpinned by physiological, morphological, and behavioural adaptations, which may require the evolution of novelty, but the genetic basis in many cases is unknown. In orchid bees, males acquire chemical compounds from the environment that act as pheromone-like bouquets (perfumes) during courtship displays. This process could be costly, potentially due to the cognitive demands of learning and the physiological demands of collecting a mix of extrinsic chemical compounds that may require detoxification. Furthermore, a novel trait, a specialized perfume pouch in the hind leg, is required for compound storage. We studied gene expression in the brain, hind leg, and Malpighian tubules-a tissue involved in detoxification-to investigate changes in gene expression following perfume collection. We detected upregulation of genes enriched in functions related to transcription, odorant binding, and receptor activity in the Malpighian tubules. On the other hand, we did not find any evidence for learning processes following perfume collection, or gene expression changes in the hind leg, perhaps due to constitutive expression, or the age of the sampled bees. We did identify high expression of chemosensory proteins in the hind legs, which we suggest could play a role in perfume collection or storage, with further functional studies necessary to determine their binding properties and potential physiological importance. Los rasgos sexuales secundarios pueden servir como indicadores de calidad de la pareja, y en algunos casos la honestidad de la señal se mantiene por el costo de expresar el rasgo. A menudo las señales sexuales están respaldadas por adaptaciones fisiológicas, morfológicas y de comportamiento por lo tanto pueden requerir la evolución de nuevos rasgos, pero en muchos casos se desconoce la base genética. En las abejas de las orquídeas, los machos recolectan compuestos químicos del medio ambiente, los cuales actúan como feromonas (perfumes) durante el despliegue de cortejo. Este proceso podría ser costoso, posiblemente debido a las demandas cognitivas del aprendizaje y las demandas fisiológicas de recolectar una mezcla de compuestos químicos extrínsecos que pueden requerir desintoxicación. Además, se requiere la evolución de un contenedor para almacenar perfumes en la pata trasera. Para investigar los cambios en la expresión génica después de la recolección de perfume, estudiamos la expresión génica en el cerebro, la pata trasera y los túbulos de Malpighi (tejido involucrado en la desintoxicación). Encontramos varios genes regulados positivamente en los túbulos de Malpighi después de la recolección que están enriquecidos en factores de transcripción, proteínas de fijación de olores, y proteínas con actividad de receptor. Por otro lado, no encontramos ninguna evidencia de procesos de aprendizaje posteriores a la recolección de perfumes, o cambios en la expresión génica en la pata trasera, esto quizás debido a la expresión constitutiva o la edad de las abejas muestreadas. Además, identificamos una alta expresión de proteínas quimio-sensoriales en las patas traseras, que podría desempeñar un papel en la recolección o almacenamiento de perfumes. Más estudios funcionales son necesarios para determinar las propiedades de fijación de las proteínas y su potencial importancia fisiológica.

Genomic signatures of eusocial evolution in insects

The genomes of eusocial insects allow the production and regulation of highly distinct phenotypes, largely independent of genotype. Although rare, eusociality has evolved convergently in at least three insect orders (Hymenoptera, Blattodea and Coleoptera). Despite such disparate origins, eusocial phenotypes show remarkable similarity, exhibiting long-lived reproductives and short-lived sterile workers and soldiers. In this article, we review current knowledge on genomic signatures of eusocial evolution. We confirm that especially an increased regulatory complexity and the adaptive evolution of chemical communication are common to several origins of eusociality. Furthermore, colony life itself can shape genomes of divergent taxa in a similar manner. Future research should be geared towards generating more high-quality genomic resources, especially in hitherto understudied clades, such as ambrosia beetles and termites. The application of more sophisticated tools such as machine learning techniques may allow the detection of more subtle convergent genomic footprints of eusociality.

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.

Comparative analyses of mitogenomes in the social bees with insights into evolution of long inverted repeats in the Meliponini

The insect mitogenome is typically a compact circular molecule with highly conserved gene contents. Nonetheless, mitogenome structural variations have been reported in specific taxa, and gene rearrangements, usually the tRNAs, occur in different lineages. Because synapomorphies of mitogenome organizations can provide information for phylogenetic inferences, comparative analyses of mitogenomes have been given increasing attention. However, most studies use a very few species to represent the whole genus, tribe, family, or even order, overlooking potential variations at lower taxonomic levels, which might lead to some incorrect inferences. To provide new insights into mitogenome organizations and their implications for phylogenetic inference, this study conducted comparative analyses for mitogenomes of three social bee tribes (Meliponini, Bombini, and Apini) based on the phylogenetic framework with denser taxonomic sampling at the species and population levels. Comparative analyses revealed that mitogenomes of Apini and Bombini are the typical type, while those of Meliponini show diverse variations in mitogenome sizes and organizations. Large inverted repeats (IRs) cause significant gene rearrangements of protein coding genes (PCGs) and rRNAs in Indo-Malay/Australian stingless bee species. Molecular evolution analyses showed that the lineage with IRs have lower d N/ d S ratios for PCGs than lineages without IRs, indicating potential effects of IRs on the evolution of mitochondrial genes. The finding of IRs and different patterns of gene rearrangements suggested that Meliponini is a hotspot in mitogenome evolution. Unlike conserved PCGs and rRNAs whose rearrangements were found only in the mentioned lineages within Meliponini, tRNA rearrangements are common across all three tribes of social bees, and are significant even at the species level, indicating that comprehensive sampling is needed to fully understand the patterns of tRNA rearrangements, and their implications for phylogenetic inference.

Dynamic Evolution of Repetitive Elements and Chromatin States in Apis mellifera Subspecies

In this study, we elucidate the contribution of repetitive DNA sequences to the establishment of social structures in honeybees (Apis mellifera). Despite recent advancements in understanding the molecular mechanisms underlying the formation of honeybee castes, primarily associated with Notch signaling, the comprehensive identification of specific genomic cis-regulatory sequences remains elusive. Our objective is to characterize the repetitive landscape within the genomes of two honeybee subspecies, namely A. m. mellifera and A. m. ligustica. An observed recent burst of repeats in A. m. mellifera highlights a notable distinction between the two subspecies. After that, we transitioned to identifying differentially expressed DNA elements that may function as cis-regulatory elements. Nevertheless, the expression of these sequences showed minimal disparity in the transcriptome during caste differentiation, a pivotal process in honeybee eusocial organization. Despite this, chromatin segmentation, facilitated by ATAC-seq, ChIP-seq, and RNA-seq data, revealed a distinct chromatin state associated with repeats. Lastly, an analysis of sequence divergence among elements indicates successive changes in repeat states, correlating with their respective time of origin. Collectively, these findings propose a potential role of repeats in acquiring novel regulatory functions.

Comparative Evolutionary Genomics in Insects

Genome sequencing quality, in terms of both read length and accuracy, is constantly improving. By combining long-read sequencing technologies with various scaffolding techniques, chromosome-level genome assemblies are now achievable at an affordable price for non-model organisms. Insects represent an exciting taxon for studying the genomic underpinnings of evolutionary innovations, due to ancient origins, immense species-richness, and broad phenotypic diversity. Here we summarize some of the most important methods for carrying out a comparative genomics study on insects. We describe available tools and offer concrete tips on all stages of such an endeavor from DNA extraction through genome sequencing, annotation, and several evolutionary analyses. Along the way we describe important insect-specific aspects, such as DNA extraction difficulties or gene families that are particularly difficult to annotate, and offer solutions. We describe results from several examples of comparative genomics analyses on insects to illustrate the fascinating questions that can now be addressed in this new age of genomics research.

Proteinortho6: pseudo-reciprocal best alignment heuristic for graph-based detection of (co-)orthologs

Proteinortho is a widely used tool to predict (co)-orthologous groups of genes for any set of species. It finds application in comparative and functional genomics, phylogenomics, and evolutionary reconstructions. With a rapidly increasing number of available genomes, the demand for large-scale predictions is also growing. In this contribution, we evaluate and implement major algorithmic improvements that significantly enhance the speed of the analysis without reducing precision. Graph-based detection of (co-)orthologs is typically based on a reciprocal best alignment heuristic that requires an all vs. all comparison of proteins from all species under study. The initial identification of similar proteins is accelerated by introducing an alternative search tool along with a revised search strategy-the pseudo-reciprocal best alignment heuristic-that reduces the number of required sequence comparisons by one-half. The clustering algorithm was reworked to efficiently decompose very large clusters and accelerate processing. Proteinortho6 reduces the overall processing time by an order of magnitude compared to its predecessor while maintaining its small memory footprint and good predictive quality.

The PWWP domain and the evolution of unique DNA methylation toolkits in Hymenoptera

DNMT3 in Hymenoptera has a unique duplication of the essential PWWP domain. Using GST-tagged PWWP fusion proteins and histone arrays we show that these domains have gained new properties and represent the first case of PWWP domains binding to H3K27 chromatin modifications, including H3K27me3, a key modification that is important during development. Phylogenetic analyses of 107 genomes indicate that the duplicated PWWP domains separated into two sister clades, and their distinct binding capacities are supported by 3D modeling. Other features of this unique DNA methylation system include variable copies, losses, and duplications of DNMT1 and DNMT3, and combinatorial generations of DNMT3 isoforms including variants missing the catalytic domain. Some of these losses and duplications of are found only in parasitic wasps. We discuss our findings in the context of the crosstalk between DNA methylation and histone methylation, and the expanded potential of epigenomic modifications in Hymenoptera to drive evolutionary novelties.

Influence of social lifestyles on host-microbe symbioses in the bees

Microbiomes are increasingly recognised as critical for the health of an organism. In eusocial insect societies, frequent social interactions allow for high-fidelity transmission of microbes across generations, leading to closer host-microbe coevolution. The microbial communities of bees with other social lifestyles are less studied, and few comparisons have been made between taxa that vary in social structure. To address this gap, we leveraged a cloud-computing resource and publicly available transcriptomic data to conduct a survey of microbial diversity in bee samples from a variety of social lifestyles and taxa. We consistently recover the core microbes of well-studied corbiculate bees, supporting this method's ability to accurately characterise microbial communities. We find that the bacterial communities of bees are influenced by host location, phylogeny and social lifestyle, although no clear effect was found for fungal or viral microbial communities. Bee genera with more complex societies tend to harbour more diverse microbes, with Wolbachia detected more commonly in solitary tribes. We present a description of the microbiota of Euglossine bees and find that they do not share the "corbiculate core" microbiome. Notably, we find that bacteria with known anti-pathogenic properties are present across social bee genera, suggesting that symbioses that enhance host immunity are important with higher sociality. Our approach provides an inexpensive means of exploring microbiomes of a given taxa and identifying avenues for further research. These findings contribute to our understanding of the relationships between bees and their associated microbial communities, highlighting the importance of considering microbiome dynamics in investigations of bee health.

New insights into immune genes and other expanded gene families of the house fly, Musca domestica, from an improved whole genome sequence

The house fly, Musca domestica, is a pest of livestock, transmits pathogens of human diseases, and is a model organism in multiple biological research areas. The first house fly genome assembly was published in 2014 and has been of tremendous use to the community of house fly biologists, but that genome is discontiguous and incomplete by contemporary standards. To improve the house fly reference genome, we sequenced, assembled, and annotated the house fly genome using improved techniques and technologies that were not available at the time of the original genome sequencing project. The new genome assembly is substantially more contiguous and complete than the previous genome. The new genome assembly has a scaffold N50 of 12.46 Mb, which is a 50-fold improvement over the previous assembly. In addition, the new genome assembly is within 1% of the estimated genome size based on flow cytometry, whereas the previous assembly was missing nearly one-third of the predicted genome sequence. The improved genome assembly has much more contiguous scaffolds containing large gene families. To provide an example of the benefit of the new genome, we used it to investigate tandemly arrayed immune gene families. The new contiguous assembly of these loci provides a clearer picture of the regulation of the expression of immune genes, and it leads to new insights into the selection pressures that shape their evolution.

Beyond conflict: Kinship theory of intragenomic conflict predicts individual variation in altruistic behaviour

Behavioural variation is essential for animals to adapt to different social and environmental conditions. The Kinship Theory of Intragenomic Conflict (KTIC) predicts that parent-specific alleles can support different behavioural strategies to maximize allele fitness. Previous studies, including in honey bees (Apis mellifera), supported predictions of the KTIC for parent-specific alleles to promote selfish behaviour. Here, we test the KTIC prediction that for altruism-promoting genes (i.e. those that promote behaviours that support the reproductive fitness of kin), the allele with the higher altruism optimum should be selected to be expressed while the other is silenced. In honey bee colonies, workers act altruistically when tending to the queen by performing a 'retinue' behaviour, distributing the queen's mandibular pheromone (QMP) throughout the hive. Workers exposed to QMP do not activate their ovaries, ensuring they care for the queen's brood instead of competing to lay unfertilized eggs. Due to the haplodiploid genetics of honey bees, the KTIC predicts that response to QMP is favoured by the maternal genome. We report evidence for parent-of-origin effects on the retinue response behaviour, ovarian development and gene expression in brains of worker honey bees exposed to QMP, consistent with the KTIC. Additionally, we show enrichment for genes with parent-of-origin expression bias within gene regulatory networks associated with variation in bees' response to QMP. Our study demonstrates that intragenomic conflict can shape diverse social behaviours and influence expression patterns of single genes as well as gene networks.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

The Impacts of Early-Life Experience on Bee Phenotypes and Fitness

Across diverse animal species, early-life experiences have lifelong impacts on a variety of traits. The scope of these impacts, their implications, and the mechanisms that drive these effects are central research foci for a variety of disciplines in biology, from ecology and evolution to molecular biology and neuroscience. Here, we review the role of early life in shaping adult phenotypes and fitness in bees, emphasizing the possibility that bees are ideal species to investigate variation in early-life experience and its consequences at both individual and population levels. Bee early life includes the larval and pupal stages, critical time periods during which factors like food availability, maternal care, and temperature set the phenotypic trajectory for an individual's lifetime. We discuss how some common traits impacted by these experiences, including development rate and adult body size, influence fitness at the individual level, with possible ramifications at the population level. Finally, we review ways in which human alterations to the landscape may impact bee populations through early-life effects. This review highlights aspects of bees' natural history and behavioral ecology that warrant further investigation with the goal of understanding how environmental disturbances threaten these vulnerable species.

The evolutionary history of bees in time and space

Bees are the most significant pollinators of flowering plants. This partnership began ca. 120 million years ago, but the uncertainty of how and when bees spread across the planet has greatly obscured investigations of this key mutualism. We present a novel analysis of bee biogeography using extensive new genomic and fossil data to demonstrate that bees originated in Western Gondwana (Africa and South America). Bees likely originated in the Early Cretaceous, shortly before the breakup of Western Gondwana, and the early evolution of any major bee lineage is associated with either the South American or African land masses. Subsequently, bees colonized northern continents via a complex history of vicariance and dispersal. The notable early absences from large landmasses, particularly in Australia and India, have important implications for understanding the assembly of local floras and diverse modes of pollination. How bees spread around the world from their hypothesized Southern Hemisphere origin parallels the histories of numerous flowering plant clades, providing an essential step to studying the evolution of angiosperm pollination syndromes in space and time.

Epigenetic regulation and division of labor in social insects

Analogous to multicellular organisms, social insect colonies are characterized by division of labor with queens and workers reflecting germline and soma, respectively. In multicellular organisms, such division is achieved through epigenetic factors regulating cell differentiation during development. Analogously, epigenetic regulation is postulated to regulate caste differences in social insects. We summarize recent findings about the role of epigenetics in social insects, focusing on DNA methylation and histone modifications. We specifically address (i) queen versus worker caste differentiation, (ii) queen versus worker caste differences, and (iii) division of labor among workers. Our review provides an overview of an exciting and controversially discussed field in developmental and molecular biology. It shows that our current understanding about the role of epigenetics in regulating division of labor in social insects is still fragmentary but that refined methods with well-replicated samples and targeted questions offer promising insights into this emerging field of socio-epigenomics.

The Relationship between Cognition and Brain Size or Neuron Number

The comparative approach is a powerful way to explore the relationship between brain structure and cognitive function. Thus far, the field has been dominated by the assumption that a bigger brain somehow means better cognition. Correlations between differences in brain size or neuron number between species and differences in specific cognitive abilities exist, but these correlations are very noisy. Extreme differences exist between clades in the relationship between either brain size or neuron number and specific cognitive abilities. This means that correlations become weaker, not stronger, as the taxonomic diversity of sampled groups increases. Cognition is the outcome of neural networks. Here we propose that considering plausible neural network models will advance our understanding of the complex relationships between neuron number and different aspects of cognition. Computational modelling of networks suggests that adding pathways, or layers, or changing patterns of connectivity in a network can all have different specific consequences for cognition. Consequently, models of computational architecture can help us hypothesise how and why differences in neuron number might be related to differences in cognition. As methods in connectomics continue to improve and more structural information on animal brains becomes available, we are learning more about natural network structures in brains, and we can develop more biologically plausible models of cognitive architecture. Natural animal diversity then becomes a powerful resource to both test the assumptions of these models and explore hypotheses for how neural network structure and network size might delimit cognitive function.

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.

Evolutionarily stable strategy analysis and its links to demography and genetics through invasion fitness

Evolutionarily stable strategy (ESS) analysis pioneered by Maynard Smith and Price took off in part because it often does not require explicit assumptions about the genetics and demography of a population in contrast to population genetic models. Though this simplicity is useful, it obscures the degree to which ESS analysis applies to populations with more realistic genetics and demography: for example, how does ESS analysis handle complexities such as kin selection, group selection and variable environments when phenotypes are affected by multiple genes? In this paper, I review the history of the ESS concept and show how early uncertainty about the method lead to important mathematical theory linking ESS analysis to general population genetic models. I use this theory to emphasize the link between ESS analysis and the concept of invasion fitness. I give examples of how invasion fitness can measure kin selection, group selection and the evolution of linked modifier genes in response to variable environments. The ESSs in these examples depend crucially on demographic and genetic parameters, which highlights how ESS analysis will continue to be an important tool in understanding evolutionary patterns as new models address the increasing abundance of genetic and long-term demographic data in natural populations. This article is part of the theme issue 'Half a century of evolutionary games: a synthesis of theory, application and future directions'.

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.

Positive selection in cytochrome P450 genes is associated with gonad phenotype and mating strategy in social bees

The honey bee, Apis mellifera differs from all other social bees in its gonad phenotype and mating strategy. Honey bee queens and drones have tremendously enlarged gonads, and virgin queens mate with several males. In contrast, in all the other bees, the male and female gonads are small, and the females mate with only one or very few males, thus, suggesting an evolutionary and developmental link between gonad phenotype and mating strategy. RNA-seq comparisons of A. mellifera larval gonads revealed 870 genes as differentially expressed in queens versus workers and drones. Based on Gene Ontology enrichment we selected 45 genes for comparing the expression levels of their orthologs in the larval gonads of the bumble bee Bombus terrestris and the stingless bee, Melipona quadrifasciata, which revealed 24 genes as differentially represented. An evolutionary analysis of their orthologs in 13 solitary and social bee genomes revealed four genes with evidence of positive selection. Two of these encode cytochrome P450 proteins, and their gene trees indicated a lineage-specific evolution in the genus Apis, indicating that cytochrome P450 genes may be involved in the evolutionary association of polyandry and the exaggerated gonad phenotype in social bees.

Molecular patterns and processes in evolving sociality: lessons from insects

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

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

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

Wood feeding and social living: Draft genome of the subterranean termite Reticulitermes lucifugus (Blattodea; Termitoidae)

Termites (Insecta, Blattodea, Termitoidae) are a widespread and diverse group of eusocial insects known for their ability to digest wood matter. Herein, we report the draft genome of the subterranean termite Reticulitermes lucifugus, an economically important species and among the most studied taxa with respect to eusocial organization and mating system. The final assembly (~813 Mb) covered up to 88% of the estimated genome size and, in agreement with the Asexual Queen Succession Mating System, it was found completely homozygous. We predicted 16,349 highly supported gene models and 42% of repetitive DNA content. Transposable elements of R. lucifugus show similar evolutionary dynamics compared to that of other termites, with two main peaks of activity localized at 25% and 8% of Kimura divergence driven by DNA, LINE and SINE elements. Gene family turnover analyses identified multiple instances of gene duplication associated with R. lucifugus diversification, with significant lineage-specific gene family expansions related to development, perception and nutrient metabolism pathways. Finally, we analysed P450 and odourant receptor gene repertoires in detail, highlighting the large diversity and dynamical evolutionary history of these proteins in the R. lucifugus genome. This newly assembled genome will provide a valuable resource for further understanding the molecular basis of termites biology as well as for pest control.

The neuroecology of olfaction in bees

The focus of bee neuroscience has for a long time been on only a handful of social honeybee and bumblebee species, out of thousands of bees species that have been described. On the other hand, information about the chemical ecology of bees is much more abundant. Here we attempted to compile the scarce information about olfactory systems of bees across species. We also review the major categories of intra- and inter-specific olfactory behaviors of bees, with specific focus on recent literature. We finish by discussing the most promising avenues for bee olfactory research in the near future.

Maternal body condition and season influence RNA deposition in the oocytes of alfalfa leafcutting bees (Megachile rotundata)

Maternal effects are an important source of phenotypic variance, whereby females influence offspring developmental trajectory beyond direct genetic contributions, often in response to changing environmental conditions. However, relatively little is known about the mechanisms by which maternal experience is translated into molecular signals that shape offspring development. One such signal may be maternal RNA transcripts (mRNAs and miRNAs) deposited into maturing oocytes. These regulate the earliest stages of development of all animals, but are understudied in most insects. Here we investigated the effects of female internal (body condition) and external (time of season) environmental conditions on maternal RNA in the maturing oocytes and 24-h-old eggs (24-h eggs) of alfalfa leafcutting bees. Using gene expression and WGCNA analysis, we found that females adjust the quantity of mRNAs related to protein phosphorylation, transcriptional regulation, and nuclease activity deposited into maturing oocytes in response to both poor body condition and shorter day lengths that accompany the late season. However, the magnitude of these changes was higher for time of season. Females also adjusted miRNA deposition in response to seasonal changes, but not body condition. We did not observe significant changes in maternal RNAs in response to either body condition or time of season in 24-h eggs, which were past the maternal-to-zygotic transition. Our results suggest that females adjust the RNA transcripts they provide for offspring to regulate development in response to both internal and external environmental cues. Variation in maternal RNAs may, therefore, be important for regulating offspring phenotype in response to environmental change.

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.

Current Knowledge on Bee Innate Immunity Based on Genomics and Transcriptomics

As important pollinators, bees play a critical role in maintaining the balance of the ecosystem and improving the yield and quality of crops. However, in recent years, the bee population has significantly declined due to various pathogens and environmental stressors including viruses, bacteria, parasites, and increased pesticide application. The above threats trigger or suppress the innate immunity of bees, their only immune defense system, which is essential to maintaining individual health and that of the colony. In addition, bees can be divided into solitary and eusocial bees based on their life traits, and eusocial bees possess special social immunities, such as grooming behavior, which cooperate with innate immunity to maintain the health of the colony. The omics approach gives us an opportunity to recognize the distinctive innate immunity of bees. In this regard, we summarize innate bee immunity from a genomic and transcriptomic perspective. The genetic characteristics of innate immunity were revealed by the multiple genomes of bees with different kinds of sociality, including honeybees, bumblebees, wasps, leaf-cutter bees, and so on. Further substantial transcriptomic data of different tissues from diverse bees directly present the activation or suppression of immune genes under the infestation of pathogens or toxicity of pesticides.

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.

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

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

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

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

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

Simple Summary

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

Abstract

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

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

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

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

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

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

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

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

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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