Eusocial insects as emerging models for behavioural epigenetics

Gut microbiota influences onset of foraging-related behavior but not physiological hallmarks of division of labor in honeybees

Gut microbes can impact cognition and behavior, but whether they regulate the division of labor in animal societies is unknown. We addressed this question using honeybees since they exhibit division of labor between nurses and foragers and because their gut microbiota can be manipulated. Using automated behavioral tracking and controlling for co-housing effects, we show that gut microbes influence the age at which bees start expressing foraging-like behaviors in the laboratory but have no effects on the time spent in a foraging arena and number of foraging trips. Moreover, the gut microbiota did not influence hallmarks of behavioral maturation such as body weight, cuticular hydrocarbon profile, hypopharyngeal gland size, gene expression, and the proportion of bees maturing into foragers. Overall, this study shows that the honeybee gut microbiota plays a role in controlling the onset of foraging-related behavior without permanent consequences on colony-level division of labor and several physiological hallmarks of behavioral maturation.

IMPORTANCE

The honeybee is emerging as a model system for studying gut microbiota-host interactions. Previous studies reported gut microbiota effects on multiple worker bee phenotypes, all of which change during behavioral maturation-the transition from nursing to foraging. We tested whether the documented effects may stem from an effect of the microbiota on behavioral maturation. The gut microbiota only subtly affected maturation: it accelerated the onset of foraging without affecting the overall proportion of foragers or their average output. We also found no effect of the microbiota on host weight, cuticular hydrocarbon (CHC) profile, hypopharyngeal gland size, and the expression of behavioral maturation-related genes. These results are inconsistent with previous studies reporting effects of the gut microbiota on bee weight and CHC profile. Our experiments revealed that co-housed bees tend to converge in behavior and physiology, suggesting that spurious associations may emerge when rearing environments are not replicated sufficiently or accounted for analytically.

Caste-biased patterns of brain investment in the subterranean termite Reticulitermes flavipes

Investment into neural tissue is expected to reflect the specific sensory and behavioral capabilities of a particular organism. Termites are eusocial insects that exhibit a caste system in which individuals can develop into one of several morphologically and behaviorally distinct castes. However, it is unclear to what extent these differences between castes are reflected in the anatomy of the brain. To address this question, we used deformation-based morphometry to conduct pairwise comparisons between the brains of different castes in the eastern subterranean termite, Reticulitermes flavipes. Workers exhibited enlargement in the antennal lobes and mushroom bodies, while reproductives showed increased investment into the optic lobes and central body. In addition, caste-specific enlargement was observed in regions that could not be mapped to distinct neuropils, most notably in soldiers. These findings demonstrate a significant influence of caste development on brain anatomy in termites alongside convergence with eusocial hymenopteran systems.

Characterization of the Pristionchus pacificus "epigenetic toolkit" reveals the evolutionary loss of the histone methyltransferase complex PRC2

Comparative approaches have revealed both divergent and convergent paths to achieving shared developmental outcomes. Thus, only through assembling multiple case studies can we understand biological principles. Yet, despite appreciating the conservation-or lack thereof-of developmental networks, the conservation of epigenetic mechanisms regulating these networks is poorly understood. The nematode Pristionchus pacificus has emerged as a model system of plasticity and epigenetic regulation as it exhibits a bacterivorous or omnivorous morph depending on its environment. Here, we determined the "epigenetic toolkit" available to P. pacificus as a resource for future functional work on plasticity, and as a comparison with Caenorhabditis elegans to investigate the conservation of epigenetic mechanisms. Broadly, we observed a similar cast of genes with putative epigenetic function between C. elegans and P. pacificus. However, we also found striking differences. Most notably, the histone methyltransferase complex PRC2 appears to be missing in P. pacificus. We described the deletion/pseudogenization of the PRC2 genes mes-2 and mes-6 and concluded that both were lost in the last common ancestor of P. pacificus and a related species P. arcanus. Interestingly, we observed the enzymatic product of PRC2 (H3K27me3) by mass spectrometry and immunofluorescence, suggesting that a currently unknown methyltransferase has been co-opted for heterochromatin silencing. Altogether, we have provided an inventory of epigenetic genes in P. pacificus to compare with C. elegans. This inventory will enable reverse-genetic experiments related to plasticity and has revealed the first loss of PRC2 in a multicellular organism.

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.

miR-2765 Modulates the Seasonal Polyphenism in Cacopsylla chinensis by Targeting a Novel Cold Rreceptor CcTRPC3

Polyphenism is a beneficial way in organisms to better cope with changing circumstances and is a hot topic in entomology, evolutionary biology, and ecology. Until now, this phenomenon has been proven to be season-, density-, and diet-dependent; however, there are very few reports on temperature regulation. Cacopsylla chinensis showed seasonal polyphenism, namely as summer- and winter-form, with obvious diversity in phenotypic characteristics in response to seasonal variation. Previous studies have found that low temperature in autumn is an extremely important element in inducing summer-form change to winter-form, but the underlying regulatory mechanism is still a mystery. Herein, we provided the initial evidence that the third instar of the summer-form is the critical period for developing to the winter-form, and 10 °C induces this transition by affecting the total pigment, chitin level, and thickness of the cuticle. Second, CcTPRC3 was proven to function as a novel cold receptor to control this seasonal polyphenism. Moreover, miR-2765 was found to mediate seasonal polyphenism by inhibiting CcTRPC3 expression. Last, we found that cuticle binding proteins CcCPR4 and CcCPR9 function as the downstream signals of CcTRPC3 to regulate the seasonal polyphenism in C. chinensis. In conclusion, our results displayed a novel signal pathway of miR-2765 and CcTRPC3 for the regulation of seasonal polyphenism in C. chinensis. These findings provide insights into the comprehensive analysis of insect polyphenism and are useful in developing potential strategies to block the phase transition for the pest control of C. chinensis.

Hormonal gatekeeping via the blood-brain barrier governs caste-specific behavior in ants

Here, we reveal an unanticipated role of the blood-brain barrier (BBB) in regulating complex social behavior in ants. Using scRNA-seq, we find localization in the BBB of a key hormone-degrading enzyme called juvenile hormone esterase (Jhe), and we show that this localization governs the level of juvenile hormone (JH3) entering the brain. Manipulation of the Jhe level reprograms the brain transcriptome between ant castes. Although ant Jhe is retained and functions intracellularly within the BBB, we show that Drosophila Jhe is naturally extracellular. Heterologous expression of ant Jhe into the Drosophila BBB alters behavior in fly to mimic what is seen in ants. Most strikingly, manipulation of Jhe levels in ants reprograms complex behavior between worker castes. Our study thus uncovers a remarkable, potentially conserved role of the BBB serving as a molecular gatekeeper for a neurohormonal pathway that regulates social behavior.

DNMT1 mutant ants develop normally but have disrupted oogenesis

Although DNA methylation is an important gene regulatory mechanism in mammals, its function in arthropods remains poorly understood. Studies in eusocial insects have argued for its role in caste development by regulating gene expression and splicing. However, such findings are not always consistent across studies, and have therefore remained controversial. Here we use CRISPR/Cas9 to mutate the maintenance DNA methyltransferase DNMT1 in the clonal raider ant, Ooceraea biroi. Mutants have greatly reduced DNA methylation, but no obvious developmental phenotypes, demonstrating that, unlike mammals, ants can undergo normal development without DNMT1 or DNA methylation. Additionally, we find no evidence of DNA methylation regulating caste development. However, mutants are sterile, whereas in wild-type ants, DNMT1 is localized to the ovaries and maternally provisioned into nascent oocytes. This supports the idea that DNMT1 plays a crucial but unknown role in the insect germline.

Histone 4 lysine 5/12 acetylation enables developmental plasticity of Pristionchus mouth form

Development can be altered to match phenotypes with the environment, and the genetic mechanisms that direct such alternative phenotypes are beginning to be elucidated. Yet, the rules that govern environmental sensitivity vs. invariant development, and potential epigenetic memory, remain unknown. Here, we show that plasticity of nematode mouth forms is determined by histone 4 lysine 5 and 12 acetylation (H4K5/12ac). Acetylation in early larval stages provides a permissive chromatin state, which is susceptible to induction during the critical window of environmental sensitivity. As development proceeds deacetylation shuts off switch gene expression to end the critical period. Inhibiting deacetylase enzymes leads to fixation of prior developmental trajectories, demonstrating that histone modifications in juveniles can carry environmental information to adults. Finally, we provide evidence that this regulation was derived from an ancient mechanism of licensing developmental speed. Altogether, our results show that H4K5/12ac enables epigenetic regulation of developmental plasticity that can be stored and erased by acetylation and deacetylation, respectively.

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

What is in a model?

After reading contradictory claims of model status for some insect species, we feel a brief discussion of the topic may be useful. Here, we document a few examples where clarity on model status seems to be lacking, briefly review work on widely recognized models, and offer criteria for including any given species as a model organism.

Chromatin accessibility-based characterisation of brain gene regulatory networks in three distinct honey bee polyphenisms

The honey bee genome has the capacity to produce three phenotypically distinct organisms (two diploid female castes: queen and worker, and a haploid male drone). Previous studies have implicated metabolic flux acting via epigenetic regulation in directing nutrition-driven phenotypic plasticity in the honey bee. However, the cis-acting DNA regulatory elements that establish tissue and polyphenism -specific epigenomes and gene expression programmes, remain unclear. Using a high resolution multiomic approach including assay for transposase-accessible chromatin by sequencing (ATAC-seq), RNA-seq and ChIP-seq, we produce the first genome-wide maps of the regulatory landscape across all three adult honey bee phenotypes identifying > 5000 regulatory regions in queen, 7500 in worker and 6500 in drone, with the vast majority of these sites located within intronic regions. These regions are defined by positive enrichment of H3K27ac and depletion of H3K4me3 and show a positive correlation with gene expression. Using ATAC-seq footprinting we determine queen, worker and drone -specific transcription factor occupancy and uncover novel phenotype-specific regulatory networks identifying two key nuclear receptors that have previously been implicated in caste-determination and adult behavioural maturation in honey bees; ecdysone receptor and ultraspiracle. Collectively, this study provides novel insights into key gene regulatory networks that are associated with these distinct polyphenisms in the honey bee.

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

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

Behavior Individuality: A Focus on Drosophila melanogaster

Among individuals, behavioral differences result from the well-known interplay of nature and nurture. Minute differences in the genetic code can lead to differential gene expression and function, dramatically affecting developmental processes and adult behavior. Environmental factors, epigenetic modifications, and gene expression and function are responsible for generating stochastic behaviors. In the last decade, the advent of high-throughput sequencing has facilitated studying the genetic basis of behavior and individuality. We can now study the genomes of multiple individuals and infer which genetic variations might be responsible for the observed behavior. In addition, the development of high-throughput behavioral paradigms, where multiple isogenic animals can be analyzed in various environmental conditions, has again facilitated the study of the influence of genetic and environmental variations in animal personality. Mainly, Drosophila melanogaster has been the focus of a great effort to understand how inter-individual behavioral differences emerge. The possibility of using large numbers of animals, isogenic populations, and the possibility of modifying neuronal function has made it an ideal model to search for the origins of individuality. In the present review, we will focus on the recent findings that try to shed light on the emergence of individuality with a particular interest in D. melanogaster.

Gene expression and epigenetic responses of the marine Cladoceran, Evadne nordmanni, and the copepod, Acartia clausi, to elevated CO2

Characterizing the capacity of marine organisms to adapt to climate change related drivers (e.g., pCO2 and temperature), and the possible rate of this adaptation, is required to assess their resilience (or lack thereof) to these drivers. Several studies have hypothesized that epigenetic markers such as DNA methylation, histone modifications and noncoding RNAs, act as drivers of adaptation in marine organisms, especially corals. However, this hypothesis has not been tested in zooplankton, a keystone organism in marine food webs. The objective of this study is to test the hypothesis that acute ocean acidification (OA) exposure alters DNA methylation in two zooplanktonic species-copepods (Acartia clausii) and cladocerans (Evadne nordmanii). We exposed these two species to near-future OA conditions (400 and 900 ppm pCO2) for 24 h and assessed transcriptional and DNA methylation patterns using RNA sequencing and Reduced Representation Bisulfite Sequencing (RRBS). OA exposure caused differential expression of genes associated with energy metabolism, cytoskeletal and extracellular matrix functions, hypoxia and one-carbon metabolism. Similarly, OA exposure also caused altered DNA methylation patterns in both species but the effect of these changes on gene expression and physiological effects remains to be determined. The results from this study form the basis for studies investigating the potential role of epigenetic mechanisms in OA induced phenotypic plasticity and/or adaptive responses in zooplanktonic organisms.

Tramtrack acts during late pupal development to direct ant caste identity

A key question in the rising field of neuroepigenetics is how behavioral plasticity is established and maintained in the developing CNS of multicellular organisms. Behavior is controlled through systemic changes in hormonal signaling, cell-specific regulation of gene expression, and changes in neuronal connections in the nervous system, however the link between these pathways is unclear. In the ant Camponotus floridanus, the epigenetic corepressor CoREST is a central player in experimentally-induced reprogramming of caste-specific behavior, from soldier (Major worker) to forager (Minor worker). Here, we show this pathway is engaged naturally on a large genomic scale during late pupal development targeting multiple genes differentially expressed between castes, and central to this mechanism is the protein tramtrack (ttk), a DNA binding partner of CoREST. Caste-specific differences in DNA binding of ttk co-binding with CoREST correlate with caste-biased gene expression both in the late pupal stage and immediately after eclosion. However, we find a unique set of exclusive Minor-bound genes that show ttk pre-binding in the late pupal stage preceding CoREST binding, followed by caste-specific gene repression on the first day of eclosion. In addition, we show that ttk binding correlates with neurogenic Notch signaling, and that specific ttk binding between castes is enriched for regulatory sites associated with hormonal function. Overall our findings elucidate a pathway of transcription factor binding leading to a repressive epigenetic axis that lies at the crux of development and hormonal signaling to define worker caste identity in C. floridanus.

DNA methylation is not a driver of gene expression reprogramming in young honey bee workers

The presence of DNA methylation marks within genic intervals, also called gene body methylation, is an evolutionarily-conserved epigenetic hallmark of animal and plant methylomes. In social insects, gene body methylation is thought to contribute to behavioural plasticity, for example between foragers and nurse workers, by modulating gene expression. However, recent studies have suggested that the majority of DNA methylation is sequence-specific, and therefore cannot act as a flexible mediator between environmental cues and gene expression. To address this paradox, we examined whole-genome methylation patterns in the brains and ovaries of young honey bee workers that had been subjected to divergent social contexts: the presence or absence of the queen. Although these social contexts are known to bring about extreme changes in behavioral and reproductive traits through differential gene expression, we found no significant differences between the methylomes of workers from queenright and queenless colonies. In contrast, thousands of regions were differentially methylated between colonies, and these differences were not associated with differential gene expression in the subset of genes examined. Methylation patterns were highly similar between brain and ovary tissues and only differed in nine regions. These results strongly indicate that DNA methylation is not a driver of differential gene expression between tissues or behavioral morphs. Finally, despite the lack of difference in methylation patterns, queen presence affected the expression of all four DNA methyltransferase genes, suggesting that these enzymes have roles beyond DNA methylation. Therefore, the functional role of DNA methylation in social insect genomes remains an open question.

Contribution of Epigenetic Mechanisms in the Regulation of Environmentally-Induced Polyphenism in Insects

Simple Summary

Polyphenism is a widespread phenomenon in insects that allows organisms to produce alternative and discrete phenotypes in response to environmental conditions. Epigenetic mechanisms, including histone post-translational modifications, DNA methylation and non-coding RNAs, are essential mechanisms that can promote rapid and flexible changes in the expression of transcriptional programs associated with the production of alternative phenotypes. This review summarizes knowledge regarding the contribution of those mechanisms in the regulation of the most-studied examples of polyphenism in insects.

Abstract

Many insect species display a remarkable ability to produce discrete phenotypes in response to changes in environmental conditions. Such phenotypic plasticity is referred to as polyphenism. Seasonal, dispersal and caste polyphenisms correspond to the most-studied examples that are environmentally-induced in insects. Cues that induce such dramatic phenotypic changes are very diverse, ranging from seasonal cues, habitat quality changes or differential larval nutrition. Once these signals are perceived, they are transduced by the neuroendocrine system towards their target tissues where gene expression reprogramming underlying phenotypic changes occur. Epigenetic mechanisms are key regulators that allow for genome expression plasticity associated with such developmental switches. These mechanisms include DNA methylation, chromatin remodelling and histone post-transcriptional modifications (PTMs) as well as non-coding RNAs and have been studied to various extents in insect polyphenism. Differential patterns of DNA methylation between phenotypes are usually correlated with changes in gene expression and alternative splicing events, especially in the cases of dispersal and caste polyphenism. Combinatorial patterns of histone PTMs provide phenotype-specific epigenomic landscape associated with the expression of specific transcriptional programs, as revealed during caste determination in honeybees and ants. Alternative phenotypes are also usually associated with specific non-coding RNA profiles. This review will provide a summary of the current knowledge of the epigenetic changes associated with polyphenism in insects and highlights the potential for these mechanisms to be key regulators of developmental transitions triggered by environmental cues.

A review on the DNA methyltransferase family of insects: Aspect and prospects

The DNA methyltransferase family contains a conserved set of DNA-modifying enzymatic proteins. They are responsible for epigenetic gene modulation, such as transcriptional silencing, transcription activation, and post-transcriptional modulation. Recent research has revealed that the canonical DNA methyltransferases (DNMTs) biological roles go beyond their traditional functions of establishing and maintaining DNA methylation patterns. Although a complete DNA methylation toolkit is absent in most insect orders, recent evidence indicates the de novo DNA methylation and maintenance function remain conserved. Studies using various molecular approaches provided evidence that DNMTs are multi-functional proteins. However, still in-depth studies on their biological role lack due to the least studied area in insects. Here, we review the DNA methylation toolkit of insects, focusing on recent research on various insect orders, which exhibit DNA methylation at different levels, and for which DNMTs functional studies have become available in recent years. We survey research on the potential roles of DNMTs in the regulation of gene transcription in insect species. DNMTs participate in different physiological processes by interacting with other epigenetic factors. Future studies on insect's DNMTs will benefit to understand developmental processes, responses to various stimuli, and adaptability of insects to different environmental conditions.

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

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

(Epi)Genetic Mechanisms Underlying the Evolutionary Success of Eusocial Insects

Eusocial insects, such as bees, ants, and wasps of the Hymenoptera and termites of the Blattodea, are able to generate remarkable diversity in morphology and behavior despite being genetically uniform within a colony. Most eusocial insect species display caste structures in which reproductive ability is possessed by a single or a few queens while all other colony members act as workers. However, in some species, caste structure is somewhat plastic, and individuals may switch from one caste or behavioral phenotype to another in response to certain environmental cues. As different castes normally share a common genetic background, it is believed that much of this observed within-colony diversity results from transcriptional differences between individuals. This suggests that epigenetic mechanisms, featured by modified gene expression without changing genes themselves, may play an important role in eusocial insects. Indeed, epigenetic mechanisms such as DNA methylation, histone modifications and non-coding RNAs, have been shown to influence eusocial insects in multiple aspects, along with typical genetic regulation. This review summarizes the most recent findings regarding such mechanisms and their diverse roles in eusocial insects.

10-hydroxy-2E-decenoic acid (10HDA) does not promote caste differentiation in Melipona scutellaris stingless bees

In bees from genus Melipona, differential feeding is not enough to fully explain female polyphenism. In these bees, there is a hypothesis that in addition to the environmental component (food), a genetic component is also involved in caste differentiation. This mechanism has not yet been fully elucidated and may involve epigenetic and metabolic regulation. Here, we verified that the genes encoding histone deacetylases HDAC1 and HDAC4 and histone acetyltransferase KAT2A were expressed at all stages of Melipona scutellaris, with fluctuations between developmental stages and castes. In larvae, the HDAC genes showed the same profile of Juvenile Hormone titers-previous reported-whereas the HAT gene exhibited the opposite profile. We also investigated the larvae and larval food metabolomes, but we did not identify the putative queen-fate inducing compounds, geraniol and 10-hydroxy-2E-decenoic acid (10HDA). Finally, we demonstrated that the histone deacetylase inhibitor 10HDA-the major lipid component of royal jelly and hence a putative regulator of honeybee caste differentiation-was unable to promote differentiation in queens in Melipona scutellaris. Our results suggest that epigenetic and hormonal regulations may act synergistically to drive caste differentiation in Melipona and that 10HDA is not a caste-differentiation factor in Melipona scutellaris.

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

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

Reversible plasticity in brain size, behaviour and physiology characterizes caste transitions in a socially flexible ant (Harpegnathos saltator)

Phenotypic plasticity allows organisms to respond to changing environments throughout their lifetime, but these changes are rarely reversible. Exceptions occur in relatively long-lived vertebrate species that exhibit seasonal plasticity in brain size, although similar changes have not been identified in short-lived species, such as insects. Here, we investigate brain plasticity in reproductive workers of the ant Harpegnathos saltator. Unlike most ant species, workers of H. saltator are capable of sexual reproduction, and they compete in a dominance tournament to establish a group of reproductive workers, termed 'gamergates'. We demonstrated that, compared to foragers, gamergates exhibited a 19% reduction in brain volume in addition to significant differences in behaviour, ovarian status, venom production, cuticular hydrocarbon profile, and expression profiles of related genes. In experimentally manipulated gamergates, 6-8 weeks after being reverted back to non-reproductive status their phenotypes shifted to the forager phenotype across all traits we measured, including brain volume, a trait in which changes were previously shown to be irreversible in honeybees and Drosophila. Brain plasticity in H. saltator is therefore more similar to that found in some long-lived vertebrates that display reversible changes in brain volume throughout their lifetimes.

Genetic basis of chemical communication in eusocial insects

Social behavior is one of the most fascinating and complex behaviors in humans and animals. A fundamental process of social behavior is communication among individuals. It relies on the capability of the nervous system to sense, process, and interpret various signals (e.g., pheromones) and respond with appropriate decisions and actions. Eusocial insects, including ants, some bees, some wasps, and termites, display intriguing cooperative social behavior. Recent advances in genetic and genomic studies have revealed key genes that are involved in pheromone synthesis, chemosensory perception, and physiological and behavioral responses to varied pheromones. In this review, we highlight the genes and pathways that regulate queen pheromone-mediated social communication, discuss the evolutionary changes in genetic systems, and outline prospects of functional studies in sociobiology.

Comparative mitochondrial genomics of sponge-dwelling snapping shrimps in the genus Synalpheus: Exploring differences between eusocial and non-eusocial species and insights into phylogenetic relationships in caridean shrimps

The genus Synalpheus is a cosmopolitan clade of marine shrimps found in most tropical regions. Species in this genus exhibit a range of social organizations, including pair-forming, communal breeding, and eusociality, the latter only known to have evolved within this genus in the marine realm. This study examines the complete mitochondrial genomes of seven species of Synalpheus and explores differences between eusocial and non-eusocial species considering that eusociality has been shown before to affect the strength of purifying selection in mitochondrial protein coding genes. The AT-rich mitochondrial genomes of Synalpheus range from 15,421 bp to 15,782 bp in length and comprise, invariably, 13 protein-coding genes (PCGs), two ribosomal RNA genes, and 22 transfer RNA genes. A 648 bp to 994 bp long intergenic space is assumed to be the D-loop. Mitochondrial gene synteny is identical among the studied shrimps. No major differences occur between eusocial and non-eusocial species in nucleotide composition and codon usage profiles of PCGs and in the secondary structure of tRNA genes. Maximum likelihood phylogenetic analysis of the complete concatenated PCG complement of 90 species supports the monophyly of the genus Synalpheus and its family Alpheidae. Moreover, the monophyletic status of the caridean families Alvinocaridae, Atyidae, Thoridae, Lysmatidae, Palaemonidae, and Pandalidae within caridean shrimps are fully or highly supported by the analysis. We therefore conclude that mitochondrial genomes contain sufficient phylogenetic information to resolve relationships at high taxonomic levels within the Caridea. Our analysis of mitochondrial genomes in the genus Synalpheus contributes to the understanding of the coevolution between genomic architecture and sociality in caridean shrimps and other marine organisms.

Early behavioral and molecular events leading to caste switching in the ant Harpegnathos

Ant societies show a division of labor in which a queen is in charge of reproduction while nonreproductive workers maintain the colony. In Harpegnathos saltator, workers retain reproductive ability, inhibited by the queen pheromones. Following the queen loss, the colony undergoes social unrest with an antennal dueling tournament. Most workers quickly abandon the tournament while a few workers continue the dueling for months and become gamergates (pseudoqueens). However, the temporal dynamics of the social behavior and molecular mechanisms underlining the caste transition and social dominance remain unclear. By tracking behaviors, we show that the gamergate fate is accurately determined 3 d after initiation of the tournament. To identify genetic factors responsible for this commitment, we compared transcriptomes of different tissues between dueling and nondueling workers. We found that juvenile hormone is globally repressed, whereas ecdysone biosynthesis in the ovary is increased in gamergates. We show that molecular changes in the brain serve as earliest caste predictors compared with other tissues. Thus, behavioral and molecular data indicate that despite the prolonged social upheaval, the gamergate fate is rapidly established, suggesting a robust re-establishment of social structure.

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

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

Epigenetic regulation of post-embryonic development

Modifications to DNA and core histones influence chromatin organization and expression of the genome. DNA methylation plays a significant role in the regulation of multiple biological processes that regulate behavior and caste differentiation in social insects. Histone modifications play significant roles in the regulation of development and reproduction in other insects. Genes coding for acetyltransferases, deacetylases, methyltransferases, and demethylases that modify core histones have been identified in genomes of multiple insects. Studies on the function and mechanisms of action of some of these enzymes uncovered their contribution to post-embryonic development. The results from studies on epigenetic modifiers could help in the identification of inhibitors of epigenetic modifiers that could be developed to control pests and disease vectors.

A-to-I RNA editing in honeybees shows signals of adaptation and convergent evolution

Social insects exhibit extensive phenotypic diversities among the genetically similar individuals, suggesting a role for the epigenetic regulations beyond the genome level. The ADAR-mediated adenosine-to-inosine (A-to-I) RNA editing, an evolutionarily conserved mechanism, facilitates adaptive evolution by expanding proteomic diversities. Here, we characterize the A-to-I RNA editome of honeybees (Apis mellifera), identifying 407 high-confidence A-to-I editing sites. Editing is most abundant in the heads and shows signatures for positive selection. Editing behavior differs between foragers and nurses, suggesting a role for editing in caste differentiation. Although only five sites are conserved between bees and flies, an unexpectedly large number of genes exhibit editing in both species, albeit at different locations, including the nonsynonymous auto-editing of Adar. This convergent evolution, where the same target genes independently acquire recoding events in distant diverged clades, together with the signals of adaptation observed in honeybees alone, further supports the notion of recoding being adaptive.

•

Nonsynonymous editing sites in honeybees were under positive selection

•

Differential editing may contribute to the phenotypic diversity between sub-castes

•

Target genes acquire editing in different clades, suggesting convergent evolution

Queen pheromone modulates the expression of epigenetic modifier genes in the brain of honeybee workers

Pheromones are used by many insects to mediate social interactions. In the highly eusocial honeybee (Apis mellifera), queen mandibular pheromone (QMP) is involved in the regulation of the reproductive and other behaviour of workers. The molecular mechanisms by which QMP acts are largely unknown. Here, we investigate how genes responsible for epigenetic modifications to DNA, RNA and histones respond to the presence of QMP in the environment. We show that several of these genes are upregulated in the honeybee brain when workers are exposed to artificial QMP. We propose that pheromonal communication systems, such as those used by social insects, evolved to respond to environmental signals by making use of existing epigenomic machineries.

Comparative genome-wide DNA methylation analysis reveals epigenomic differences in response to heat-humidity stress in Bombyx mori

DNA methylation is an important epigenetic modification and has been shown to be involved in the response to abiotic stress. However, there are few studies on DNA methylation in insect response to environmental signals. In this study, we conducted a comprehensive comparative analysis of DNA methylation profiles between two silkworm strains with significantly different resistance to heat and humidity by whole-genome bisulfite sequencing (WGBS). We identified, in total, 2934 differentially methylated regions (DMRs) between RT_48h (resistant strain with high-temperature/humidity treatment for 48 h) and ST_48h (sensitive strain with high-temperature/humidity treatment for 48 h) under cytosine context (CG), which corresponded to 1230 DMR-related genes (DMGs), and the DMRs were primarily located in the gene body (exon and intron) region. Gene ontology (GO) and KEGG analysis showed that these DMGs were most significantly enriched in binding, cellular metabolic process, and RNA transport pathways. Moreover, 10 DMGs have been revealed to be involved in the heat-humidity stress response in the silkworm. The results of this study indicated that DNA methylation plays crucial roles in silkworm response to environmental stressors and provides important clues to identify key resistance genes in silkworm under high-temperature/humidity stress response.

Identification, expression, and artificial selection of silkworm epigenetic modification enzymes

Background

Understanding the genetic basis of phenotype variations during domestication and breeding is of great interest. Epigenetics and epigenetic modification enzymes (EMEs) may play a role in phenotypic variations; however, no comprehensive study has been performed to date. Domesticated silkworm (Bombyx mori) may be utilized as a model in determining how EMEs influence domestication traits.

Results

We identified 44 EMEs in the genome of silkworm (Bombyx mori) using homology searching. Phylogenetic analysis showed that genes in a subfamily among different animals were well clustered, and the expression pattern of EMEs is constant among Bombyx mori, Drosophila melanogaster, and Mus musculus. These are most highly expressed in brain, early embryo, and internal genitalia. By gene-related selective sweeping, we identified five BmEMEs under artificial selection during the domestication and breeding of silkworm. Among these selected genes, BmSuv4–20 and BmDNMT2 harbor selective mutations in their upstream regions that alter transcription factor-binding sites. Furthermore, these two genes are expressed higher in the testis and ovary of domesticated silkworm compared to wild silkworms, and correlations between their expression pattern and meiosis of the sperm and ova were observed.

Conclusions

The domestication of silkworm has induced artificial selection on epigenetic modification markers that may have led to phenotypic changes during domestication. We present a novel perspective to understand the genetic basis underlying animal domestication and breeding.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07155-z.

Epigenetic Molecular Mechanisms in Insects

Insects are the largest animal group on Earth both in biomass and diversity. Their outstanding success has inspired genetics and developmental research, allowing the discovery of dynamic process explaining extreme phenotypic plasticity and canalization. Epigenetic molecular mechanisms (EMMs) are vital for several housekeeping functions in multicellular organisms, regulating developmental, ontogenetic trajectories and environmental adaptations. In Insecta, EMMs are involved in the development of extreme phenotypic divergences such as polyphenisms and eusocial castes. Here, we review the history of this research field and how the main EMMs found in insects help to understand their biological processes and diversity. EMMs in insects confer them rapid response capacity allowing insect either to change with plastic divergence or to keep constant when facing different stressors or stimuli. EMMs function both at intra as well as transgenerational scales, playing important roles in insect ecology and evolution. We discuss on how EMMs pervasive influences in Insecta require not only the control of gene expression but also the dynamic interplay of EMMs with further regulatory levels, including genetic, physiological, behavioral, and environmental among others, as was earlier proposed by the Probabilistic Epigenesis model and Developmental System Theory.

Transcriptomic and Epigenomic Dynamics of Honey Bees in Response to Lethal Viral Infection

Honey bees (Apis mellifera L.) suffer from many brood pathogens, including viruses. Despite considerable research, the molecular responses and dynamics of honey bee pupae to viral pathogens remain poorly understood. Israeli Acute Paralysis Virus (IAPV) is emerging as a model virus since its association with severe colony losses. Using worker pupae, we studied the transcriptomic and methylomic consequences of IAPV infection over three distinct time points after inoculation. Contrasts of gene expression and 5 mC DNA methylation profiles between IAPV-infected and control individuals at these time points - corresponding to the pre-replicative (5 h), replicative (20 h), and terminal (48 h) phase of infection - indicate that profound immune responses and distinct manipulation of host molecular processes accompany the lethal progression of this virus. We identify the temporal dynamics of the transcriptomic response to with more genes differentially expressed in the replicative and terminal phases than in the pre-replicative phase. However, the number of differentially methylated regions decreased dramatically from the pre-replicative to the replicative and terminal phase. Several cellular pathways experienced hyper- and hypo-methylation in the pre-replicative phase and later dramatically increased in gene expression at the terminal phase, including the MAPK, Jak-STAT, Hippo, mTOR, TGF-beta signaling pathways, ubiquitin mediated proteolysis, and spliceosome. These affected biological functions suggest that adaptive host responses to combat the virus are mixed with viral manipulations of the host to increase its own reproduction, all of which are involved in anti-viral immune response, cell growth, and proliferation. Comparative genomic analyses with other studies of viral infections of honey bees and fruit flies indicated that similar immune pathways are shared. Our results further suggest that dynamic DNA methylation responds to viral infections quickly, regulating subsequent gene activities. Our study provides new insights of molecular mechanisms involved in epigenetic that can serve as foundation for the long-term goal to develop anti-viral strategies for honey bees, the most important commercial pollinator.

The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste

BACKGROUND

Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste.

RESULTS

The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication.

CONCLUSIONS

The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

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

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

Transcriptomic analysis of epigenetic modification genes in the termite Reticulitermes speratus

Eusocial insects display a caste system in which different castes are morphologically and physiologically specialized for different tasks. Recent studies have revealed that epigenetic modifications, including DNA methylation and histone modification, mediate caste determination and differentiation, longevity, and polyethism in eusocial insects. Although there has been a growing interest in the relationship between epigenetic mechanisms and phenotypic plasticity in termites, there is little information about differential expression levels among castes and expression sites for these genes in termites. Here we show royal-tissue-specific expression of epigenetic modification genes in the termite Reticulitermes speratus. Using RNA-seq, we identified 74 genes, including three DNA methyltransferases, seven sirtuins, 48 Trithorax group proteins, and 16 Polycomb group proteins. Among these genes, 15 showed king-specific expression, and 52 showed age-dependent differential expression in kings and queens. Quantitative real-time PCR revealed that DNA methyltransferase 3 is expressed specifically in the king's testis and fat body, whereas some histone modification genes are remarkably expressed in the king's testis and queen's ovary. These findings imply that epigenetic modification plays important roles in the gamete production process in termite kings and queens.

Methoprene-Induced Genes in Workers of Formosan Subterranean Termites (Coptotermes formosanus Shiraki)

Termites have a distinct polyphenism controlled by concise hormonal and molecular mechanisms. Workers undergo double molts to transform into soldiers (worker–presoldier–soldier). Juvenile hormone analogs, such as methoprene, can induce workers to transform into presoldiers. However, the molecular mechanism underlying the worker-to-presoldier transformation in Coptotermes formosanus Shiraki is still not clear. We sequenced the transcriptome of workers four days after they had fed on methoprene-treated filter paper and control group workers, which fed on acetone-treated filter paper. The transcriptome of C. formosanus was assembled using the de novo assembly method. Expression levels of unigenes in the methoprene-treated group and the control group were compared. The differentially expressed genes were further analyzed by Gene Ontology (GO) term enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Tetrapyrrole binding, oxidoreductase activity, and metal ion binding were the only three enriched GO terms. Juvenile hormone synthesis was the first ranked enriched pathway. Carbohydrate, amino acid, and lipid metabolism pathways were also enriched. These three pathways may be related to fat body development, which is critical for presoldier formation. Our results have demonstrated the significance of JH synthesis pathways, and pathways related to fat body development in the artificial induction of presoldiers.

DNA Methylation Inhibits the Expression of CFSH in Mud Crab

Crustacean female sex hormone (CFSH) is a key regulator of crustacean sex differentiation. The expression of Sp-CFSH in the mud crab Scylla paramamosain showed a tissue-specific and gender-variant pattern. To explore the role of DNA methylation in Sp-CFSH expression, the 5'-flanking region of Sp-CFSH was cloned, and one CpG island containing 12 CpG sites was found. Results of sodium bisulfite sequencing and methylated DNA immunoprecipitation showed that CpG island methylation was stable in the eyestalk ganglion during ovarian development of the females, which was significantly lower than that in the muscle of the females and in the eyestalk ganglion of the males. Such results suggested that the involvement of DNA methylation in regulating Sp-CFSH expression followed an eyestalk ganglion-specific and gender-variant pattern. The analysis of CpG dinucleotide site methylation and activity of the site-directed mutation (SDM) reporter vector further demonstrated that methylation inhibited Sp-CFSH expression by blocking the binding of transcription factor Sp1. The finding suggested for the first time the involvement of CpG methylation in the regulation of Sp-CFSH expression.

Social influence and interaction bias can drive emergent behavioural specialization and modular social networks across systems

In social systems ranging from ant colonies to human society, behavioural specialization-consistent individual differences in behaviour-is commonplace: individuals can specialize in the tasks they perform (division of labour (DOL)), the political behaviour they exhibit (political polarization) or the non-task behaviours they exhibit (personalities). Across these contexts, behavioural specialization often co-occurs with modular and assortative social networks, such that individuals tend to associate with others that have the same behavioural specialization. This raises the question of whether a common mechanism could drive co-emergent behavioural specialization and social network structure across contexts. To investigate this question, here we extend a model of self-organized DOL to account for social influence and interaction bias among individuals-social dynamics that have been shown to drive political polarization. We find that these same social dynamics can also drive emergent DOL by forming a feedback loop that reinforces behavioural differences between individuals, a feedback loop that is impacted by group size. Moreover, this feedback loop also results in modular and assortative social network structure, whereby individuals associate strongly with those performing the same task. Our findings suggest that DOL and political polarization-two social phenomena not typically considered together-may actually share a common social mechanism. This mechanism may result in social organization in many contexts beyond task performance and political behaviour.

Epigenetic Regulator CoREST Controls Social Behavior in Ants

Ants acquire distinct morphological and behavioral phenotypes arising from a common genome, underscoring the importance of epigenetic regulation. In Camponotus floridanus, "Major" workers defend the colony, but can be epigenetically reprogrammed to forage for food analogously to "Minor" workers. Here, we utilize reprogramming to investigate natural behavioral specification. Reprogramming of Majors upregulates Minor-biased genes and downregulates Major-biased genes, engaging molecular pathways fundamental to foraging behavior. We discover the neuronal corepressor for element-1-silencing transcription factor (CoREST) is upregulated upon reprogramming and required for the epigenetic switch to foraging. Genome-wide profiling during reprogramming reveals CoREST represses expression of enzymes that degrade juvenile hormone (JH), a hormone elevated upon reprogramming. High CoREST, low JH-degrader expression, and high JH levels are mirrored in natural Minors, revealing parallel mechanisms of natural and reprogrammed foraging. These results unveil chromatin regulation via CoREST as central to programming of ant social behavior, with potential far-reaching implications for behavioral epigenetics.

High-Quality Genome Assemblies Reveal Long Non-coding RNAs Expressed in Ant Brains

Ants are an emerging model system for neuroepigenetics, as embryos with virtually identical genomes develop into different adult castes that display diverse physiology, morphology, and behavior. Although a number of ant genomes have been sequenced to date, their draft quality is an obstacle to sophisticated analyses of epigenetic gene regulation. We reassembled de novo high-quality genomes for two ant species, Camponotus floridanus and Harpegnathos saltator. Using long reads enabled us to span large repetitive regions and improve genome contiguity, leading to comprehensive and accurate protein-coding annotations that facilitated the identification of a Gp-9-like gene as differentially expressed in Harpegnathos castes. The new assemblies also enabled us to annotate long non-coding RNAs in ants, revealing caste-, brain-, and developmental-stage-specific long non-coding RNAs (lncRNAs) in Harpegnathos. These upgraded genomes, along with the new gene annotations, will aid future efforts to identify epigenetic mechanisms of phenotypic and behavioral plasticity in ants.

Sex-specific changes in the aphid DNA methylation landscape

Aphids present an ideal system to study epigenetics as they can produce diverse, but genetically identical, morphs in response to environmental stimuli. Here, using whole genome bisulphite sequencing and transcriptome sequencing of the green peach aphid (Myzus persicae), we present the first detailed analysis of cytosine methylation in an aphid and investigate differences in the methylation and transcriptional landscapes of male and asexual female morphs. We found that methylation primarily occurs in a CG dinucleotide (CpG) context and that exons are highly enriched for methylated CpGs, particularly at the 3' end of genes. Methylation is positively associated with gene expression, and methylated genes are more stably expressed than unmethylated genes. Male and asexual female morphs have distinct methylation profiles. Strikingly, these profiles are divergent between the sex chromosome and the autosomes; autosomal genes are hypomethylated in males compared to asexual females, whereas genes belonging to the sex chromosome, which is haploid in males, are hypermethylated. Overall, we found correlated changes in methylation and gene expression between males and asexual females, and this correlation was particularly strong for genes located on the sex chromosome. Our results suggest that differential methylation of sex-biased genes plays a role in aphid sexual differentiation.

Genome-wide effects of social status on DNA methylation in the brain of a cichlid fish, Astatotilapia burtoni

BACKGROUND

Successful social behavior requires real-time integration of information about the environment, internal physiology, and past experience. The molecular substrates of this integration are poorly understood, but likely modulate neural plasticity and gene regulation. In the cichlid fish species Astatotilapia burtoni, male social status can shift rapidly depending on the environment, causing fast behavioral modifications and a cascade of changes in gene transcription, the brain, and the reproductive system. These changes can be permanent but are also reversible, implying the involvement of a robust but flexible mechanism that regulates plasticity based on internal and external conditions. One candidate mechanism is DNA methylation, which has been linked to social behavior in many species, including A. burtoni. But, the extent of its effects after A. burtoni social change were previously unknown.

RESULTS

We performed the first genome-wide search for DNA methylation patterns associated with social status in the brains of male A. burtoni, identifying hundreds of Differentially Methylated genomic Regions (DMRs) in dominant versus non-dominant fish. Most DMRs were inside genes supporting neural development, synapse function, and other processes relevant to neural plasticity, and DMRs could affect gene expression in multiple ways. DMR genes were more likely to be transcription factors, have a duplicate elsewhere in the genome, have an anti-sense lncRNA, and have more splice variants than other genes. Dozens of genes had multiple DMRs that were often seemingly positioned to regulate specific splice variants.

CONCLUSIONS

Our results revealed genome-wide effects of A. burtoni social status on DNA methylation in the brain and strongly suggest a role for methylation in modulating plasticity across multiple biological levels. They also suggest many novel hypotheses to address in mechanistic follow-up studies, and will be a rich resource for identifying the relationships between behavioral, neural, and transcriptional plasticity in the context of social status.

Patterns of reproductive differentiation and reproductive plasticity in the major evolutionary transition to superorganismality

Major evolutionary transitions in individuality are characterised by the formation of new levels of biological complexity from the cooperation of previously independent lower-level units. The evolution of superorganismality in insects is one such major transition, and is characterised by an extreme division of reproductive labour between ancestrally autonomous units, in the form of queen and worker castes. Here, we discuss the nature of plasticity in the emergence of castes across the major transition to superorganismality in insects. We identify key changes in plasticity which act at different levels of selection: a loss of reproductivity plasticity at the individual level is matched by a gain in plasticity at the colony level. Taking multi-level selection into consideration has important implications for formulating testable hypotheses regarding the nature of plasticity in a major transition from a lower to a higher level of biological complexity.

Investigating the evolution and development of biological complexity under the framework of epigenetics

Biological complexity is a key component of evolvability, yet its study has been hampered by a focus on evolutionary trends of complexification and inconsistent definitions. Here, we demonstrate the utility of bringing complexity into the framework of epigenetics to better investigate its utility as a concept in evolutionary biology. We first analyze the existing metrics of complexity and explore the link between complexity and adaptation. Although recently developed metrics allow for a unified framework, they omit developmental mechanisms. We argue that a better approach to the empirical study of complexity and its evolution includes developmental mechanisms. We then consider epigenetic mechanisms and their role in shaping developmental and evolutionary trajectories, as well as the development and organization of complexity. We argue that epigenetics itself could have emerged from complexity because of a need to self‐regulate. Finally, we explore hybridization complexes and hybrid organisms as potential models for studying the association between epigenetics and complexity. Our goal is not to explain trends in biological complexity but to help develop and elucidate novel questions in the investigation of biological complexity and its evolution.

This manuscript argues that biological complexity is better understood under the framework of epigenetics and that the epigenetic interactions emerge from the self‐regulation of complex systems. Hybrids are offered as models to study these properties.

The effects of the neonicotinoid imidacloprid on gene expression and DNA methylation in the buff-tailed bumblebee Bombus terrestris

Neonicotinoids are effective insecticides used on many important arable and horticultural crops. They are nicotinic acetylcholine receptor agonists which disrupt the function of insect neurons and cause paralysis and death. In addition to direct mortality, there are numerous sublethal effects of low doses of neonicotinoids on bees. We hypothesize that some of these large array of effects could be a consequence of epigenetic changes in bees induced by neonicotinoids. We compared whole methylome (BS-seq) and RNA-seq libraries of the brains of buff-tailed bumblebee Bombus terrestris workers exposed to field-realistic doses of the neonicotinoid imidacloprid to libraries from control workers. We found numerous genes which show differential expression between neonicotinoid-treated bees and control bees, but no differentially methylated cytosines in any context. We found CpG methylation to be focused mainly in exons and associated with highly expressed genes. We discuss the implications of our results for future legislation.

RNA editing is abundant and correlates with task performance in a social bumblebee

Colonies of the bumblebee Bombus terrestris are characterized by wide phenotypic variability among genetically similar full-sister workers, suggesting a major role for epigenetic processes. Here, we report a high level of ADAR-mediated RNA editing in the bumblebee, despite the lack of an ADAR1-homolog. We identify 1.15 million unique genomic sites, and 164 recoding sites residing in 100 protein coding genes, including ion channels, transporters, and receptors predicted to affect brain function and behavior. Some edited sites are similarly edited in other insects, cephalopods and even mammals. The global editing level of protein coding and non-coding transcripts weakly correlates with task performance (brood care vs. foraging), but not affected by dominance rank or juvenile hormone known to influence physiology and behavior. Taken together, our findings show that brain editing levels are high in naturally behaving bees, and may be regulated by relatively short-term effects associated with brood care or foraging activities.

Genomic data reveal high conservation but divergent evolutionary pattern of Polycomb/Trithorax group genes in arthropods

Epigenetic gene control is maintained by chromatin-associated Polycomb group (PcG) and Trithorax group (TrxG) genes, which act antagonistically via the interplay between PcG and TrxG regulation to generate silenced or permissive transcriptional states. In this study, we searched for PcG/TrxG genes in 180 arthropod genomes, covering all the sequenced arthropod genomes at the time of conducting this study, to perform a global investigation of PcG/TrxG genes in a phylogenetic frame. Results of ancestral state reconstruction analysis revealed that the ancestor of arthropod species has an almost complete repertoire of PcG/TrxG genes, and most of these genes were seldom lost above order level. The domain diversity analysis indicated that the PcG/TrxG genes show variable extent of domain structure changes; some of these changes could be associated with lineage-specific events. The likelihood ratio tests for selection pressure detected a number of PcG/TrxG genes which underwent episodic positive selection on the branch leading to the insects with holometabolous development. These results suggest that, despite their high conservation across arthropod species, different members of PcG/TrxG genes showed considerable differences in domain structure and sequence divergence in arthropod evolution. Our cross species comparisons using large-scale genomic data provide insights into divergent evolutionary pattern on highly conserved genes in arthropods.