Robust DNA Methylation in the Clonal Raider Ant Brain

Detecting differential DNA methylation from sequencing of bisulfite converted DNA of diverse species

DNA methylation is one of the most extensively studied epigenetic modifications of genomic DNA. In recent years, sequencing of bisulfite-converted DNA, particularly via next-generation sequencing technologies, has become a widely popular method to study DNA methylation. This method can be readily applied to a variety of species, dramatically expanding the scope of DNA methylation studies beyond the traditionally studied human and mouse systems. In parallel to the increasing wealth of genomic methylation profiles, many statistical tools have been developed to detect differentially methylated loci (DMLs) or differentially methylated regions (DMRs) between biological conditions. We discuss and summarize several key properties of currently available tools to detect DMLs and DMRs from sequencing of bisulfite-converted DNA. However, the majority of the statistical tools developed for DML/DMR analyses have been validated using only mammalian data sets, and less priority has been placed on the analyses of invertebrate or plant DNA methylation data. We demonstrate that genomic methylation profiles of non-mammalian species are often highly distinct from those of mammalian species using examples of honey bees and humans. We then discuss how such differences in data properties may affect statistical analyses. Based on these differences, we provide three specific recommendations to improve the power and accuracy of DML and DMR analyses of invertebrate data when using currently available statistical tools. These considerations should facilitate systematic and robust analyses of DNA methylation from diverse species, thus advancing our understanding of DNA methylation.

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.

Capture of a functionally active methyl-CpG binding domain by an arthropod retrotransposon family

The repressive capacity of cytosine DNA methylation is mediated by recruitment of silencing complexes by methyl-CpG binding domain (MBD) proteins. Despite MBD proteins being associated with silencing, we discovered that a family of arthropod Copia retrotransposons have incorporated a host-derived MBD. We functionally show how retrotransposon-encoded MBDs preferentially bind to CpG-dense methylated regions, which correspond to transposable element regions of the host genome, in the myriapod Strigamia maritima Consistently, young MBD-encoding Copia retrotransposons (CopiaMBD) accumulate in regions with higher CpG densities than other LTR-retrotransposons also present in the genome. This would suggest that retrotransposons use MBDs to integrate into heterochromatic regions in Strigamia, avoiding potentially harmful insertions into host genes. In contrast, CopiaMBD insertions in the spider Stegodyphus dumicola genome disproportionately accumulate in methylated gene bodies compared with other spider LTR-retrotransposons. Given that transposons are not actively targeted by DNA methylation in the spider genome, this distribution bias would also support a role for MBDs in the integration process. Together, these data show that retrotransposons can co-opt host-derived epigenome readers, potentially harnessing the host epigenome landscape to advantageously tune the retrotransposition process.

Changes in gene DNA methylation and expression networks accompany caste specialization and age-related physiological changes in a social insect

Social insects provide systems for studying epigenetic regulation of phenotypes, particularly with respect to differentiation of reproductive and worker castes, which typically arise from a common genetic background. The role of gene expression in caste specialization has been extensively studied, but the role of DNA methylation remains controversial. Here, we perform well replicated, integrated analyses of DNA methylation and gene expression in brains of an ant (Formica exsecta) with distinct female castes using traditional approaches (tests of differential methylation) combined with a novel approach (analysis of co-expression and co-methylation networks). We found differences in expression and methylation profiles between workers and queens at different life stages, as well as some overlap between DNA methylation and expression at the functional level. Large portions of the transcriptome and methylome are organized into "modules" of genes, some significantly associated with phenotypic traits of castes and developmental stages. Several gene co-expression modules are preserved in co-methylation networks, consistent with possible regulation of caste-specific gene expression by DNA methylation. Surprisingly, brain co-expression modules were highly preserved when compared with a previous study that examined whole-body co-expression patterns in 16 ant species, suggesting that these modules are evolutionarily conserved and for specific functions in various tissues. Altogether, these results suggest that DNA methylation participates in regulation of caste specialization and age-related physiological changes in social insects.

Ants as Object of Gerontological Research

Social insects with identical genotype that form castes with radically different lifespans are a promising model system for studying the mechanisms underlying longevity. The main direction of progressive evolution of social insects, in particular, ants, is the development of the social way of life inextricably linked with the increase in the colony size. Only in a large colony, it is possible to have a developed polyethism, create large food reserves, and actively regulate the nest microclimate. The lifespan of ants hugely varies among genetically similar queens, workers (unproductive females), and males. The main advantage of studies on insects is the determinism of ontogenetic processes, with a single genome leading to completely different lifespans in different castes. This high degree of determinacy is precisely the reason why some researchers (incorrectly) call a colony of ants the "superorganism", emphasizing the fact that during the development, depending on the community needs, ants can switch their ontogenetic programs, which influences their social roles, ability to learn (i.e., the brain [mushroom-like body] plasticity), and, respectively, the spectrum of tasks performed by a given individual. It has been shown that in many types of food behavior, older ants surpass young ones in both performing the tasks and transferring the experience. The balance between the need to reduce the "cost" of non-breeding individuals (short lifespan and small size of workers) and the benefit from experienced long-lived workers possessing useful skills (large size and "non-aging") apparently determines the differences in the lifespan and aging rate of workers in different species of ants. A large spectrum of rigidly determined ontogenetic trajectories in different castes with identical genomes and the possibility of comparison between "evolutionarily advanced" and "primitive" subfamilies (e.g., Formicinae and Ponerinae) make ants an attractive object in the studies of both normal aging and effects of anti-aging drugs.

Changes of gene expression but not cytosine methylation are associated with male parental care reflecting behavioural state, social context and individual flexibility

Behaviour is often a front line response to changing environments. Recent studies show behavioural changes are associated with changes of gene expression; however, these studies have primarily focused on discrete behavioural states. We build on these studies by addressing additional contexts that produce qualitatively similar behavioural changes. We measured levels of gene expression and cytosine methylation, which is hypothesized to regulate the transcriptional architecture of behavioural transitions, within the brain during male parental care of the burying beetle Nicrophorus vespilloides in a factorial design. Male parenting is a suitably plastic behaviour because although male N. vespilloides typically do not provide direct care (i.e. feed offspring) when females are present, levels of feeding by a male equivalent to the female can be induced by removing the female. We examined three different factors: behavioural state (caring versus non-caring), social context (with or without a female mate) and individual flexibility (if a male switched to direct care after his mate was removed). The greatest number of differentially expressed genes were associated with behavioural state, followed by social context and individual flexibility. Cytosine methylation was not associated with changes of gene expression in any of the factors. Our results suggest a hierarchical association between gene expression and the different factors, but that this process is not controlled by cytosine methylation. Our results further suggest that the extent a behaviour is transient plays an underappreciated role in determining its underpinning molecular mechanisms.

Dnmt1 is essential for egg production and embryo viability in the large milkweed bug, Oncopeltus fasciatus

BACKGROUND

The function of cytosine (DNA) methylation in insects remains inconclusive due to a lack of mutant and/or genetic studies.

RESULTS

Here, we provide evidence for the functional role of the maintenance DNA methyltransferase 1 (Dnmt1) in an insect using experimental manipulation. Through RNA interference (RNAi), we successfully posttranscriptionally knocked down Dnmt1 in ovarian tissue of the hemipteran Oncopeltus fasciatus (the large milkweed bug). Individuals depleted for Dnmt1, and subsequently DNA methylation, failed to reproduce. Eggs were inviable and declined in number, and nuclei structure of follicular epithelium was aberrant. Erasure of DNA methylation from gene or transposon element bodies did not reveal a direct causal link to steady-state mRNA levels in somatic cells. These results reveal an important function of Dnmt1 seemingly not contingent on directly controlling gene expression.

CONCLUSIONS

This study provides direct experimental evidence for a functional role of Dnmt1 in egg production and embryo viability and uncovers a trivial role, if any, for DNA methylation in control of gene expression in O. fasciatus.

Epigenetics in Insects: Genome Regulation and the Generation of Phenotypic Diversity

Epigenetic inheritance is fundamentally important to cellular differentiation and developmental plasticity. In this review, we provide an introduction to the field of molecular epigenetics in insects. Epigenetic information is passed across cell divisions through the methylation of DNA, the modification of histone proteins, and the activity of noncoding RNAs. Much of our knowledge of insect epigenetics has been gleaned from a few model species. However, more studies of epigenetic information in traditionally nonmodel taxa will help advance our understanding of the developmental and evolutionary significance of epigenetic inheritance in insects. To this end, we also provide a brief overview of techniques for profiling and perturbing individual facets of the epigenome. Doing so in diverse cellular, developmental, and taxonomic contexts will collectively help shed new light on how genome regulation results in the generation of diversity in insect form and function.

The genomic architecture and molecular evolution of ant odorant receptors

The massive expansions of odorant receptor (OR) genes in ant genomes are notable examples of rapid genome evolution and adaptive gene duplication. However, the molecular mechanisms leading to gene family expansion remain poorly understood, partly because available ant genomes are fragmentary. Here, we present a highly contiguous, chromosome-level assembly of the clonal raider ant genome, revealing the largest known OR repertoire in an insect. While most ant ORs originate via local tandem duplication, we also observe several cases of dispersed duplication followed by tandem duplication in the most rapidly evolving OR clades. We found that areas of unusually high transposable element density (TE islands) were depauperate in ORs in the clonal raider ant, and found no evidence for retrotransposition of ORs. However, OR loci were enriched for transposons relative to the genome as a whole, potentially facilitating tandem duplication by unequal crossing over. We also found that ant OR genes are highly AT-rich compared to other genes. In contrast, in flies, OR genes are dispersed and largely isolated within the genome, and we find that fly ORs are not AT-rich. The genomic architecture and composition of ant ORs thus show convergence with the unrelated vertebrate ORs rather than the related fly ORs. This might be related to the greater gene numbers and/or potential similarities in gene regulation between ants and vertebrates as compared to flies.

Towards reconstructing the ancestral brain gene-network regulating caste differentiation in ants

Specialized queens and life-time unmated workers evolved once in the common ancestor of all ants, but whether caste development across ants continues to be at least partly regulated by a single core set of genes remains obscure. We analysed brain transcriptomes from five ant species (three subfamilies) and reconstructed the origins of genes with caste-biased expression. Ancient genes predating the Neoptera were more likely to regulate gyne (virgin queen) phenotypes, while caste differentiation roles of younger, ant-lineage-specific genes varied. Transcriptome profiling showed that the ancestral network for caste-specific gene-regulation has been maintained, but that signatures of common ancestry are obscured by later modifications. Adjusting for such differences, we identified a core gene-set that: 1. consistently displayed similar directions and degrees of caste-differentiated expression, and 2. have mostly not been reported as being involved in caste differentiation. These core regulatory genes exist in the genomes of ant species that secondarily lost the queen caste, but expression differences for reproductive and sterile workers are minor and similar to social paper wasps that lack differentiated castes. Many caste-biased ant genes have caste-differentiated expression in honeybees, but directions of caste bias were uncorrelated, as expected when permanent castes evolved independently in both lineages.

Signatures of DNA Methylation across Insects Suggest Reduced DNA Methylation Levels in Holometabola

It has been experimentally shown that DNA methylation is involved in the regulation of gene expression and the silencing of transposable element activity in eukaryotes. The variable levels of DNA methylation among different insect species indicate an evolutionarily flexible role of DNA methylation in insects, which due to a lack of comparative data is not yet well-substantiated. Here, we use computational methods to trace signatures of DNA methylation across insects by analyzing transcriptomic and genomic sequence data from all currently recognized insect orders. We conclude that: 1) a functional methylation system relying exclusively on DNA methyltransferase 1 is widespread across insects. 2) DNA methylation has potentially been lost or extremely reduced in species belonging to springtails (Collembola), flies and relatives (Diptera), and twisted-winged parasites (Strepsiptera). 3) Holometabolous insects display signs of reduced DNA methylation levels in protein-coding sequences compared with hemimetabolous insects. 4) Evolutionarily conserved insect genes associated with housekeeping functions tend to display signs of heavier DNA methylation in comparison to the genomic/transcriptomic background. With this comparative study, we provide the much needed basis for experimental and detailed comparative analyses required to gain a deeper understanding on the evolution and function of DNA methylation in insects.

Recent Advances in Behavioral (Epi)Genetics in Eusocial Insects

Eusocial insects live in societies in which distinct family members serve specific roles in maintaining the colony and advancing the reproductive ability of a few select individuals. Given the genetic similarity of all colony members, the diversity of morphologies and behaviors is surprising. Social communication relies on pheromones and olfaction, as shown by mutants of orco, the universal odorant receptor coreceptor, and through electrophysiological analysis of neuronal responses to pheromones. Additionally, neurohormonal factors and epigenetic regulators play a key role in caste-specific behavior, such as foraging and caste switching. These studies start to allow an understanding of the molecular mechanisms underlying social behavior and provide a technological foundation for future studies of eusocial insects. In this review, we highlight recent findings in eusocial insects that advance our understanding of genetic and epigenetic regulations of social behavior and provide perspectives on future studies using cutting-edge technologies.

The effect of maternal care on gene expression and DNA methylation in a subsocial bee

Developmental plasticity describes the influence of environmental factors on phenotypic variation. An important mediator of developmental plasticity in many animals is parental care. Here, we examine the consequences of maternal care on offspring after the initial mass provisioning of brood in the small carpenter bee, Ceratina calcarata. Removal of the mother during larval development leads to increased aggression and avoidance in adulthood. This corresponds with changes in expression of over one thousand genes, alternative splicing of hundreds of genes, and significant changes to DNA methylation. We identify genes related to metabolic and neuronal functions that may influence developmental plasticity and aggression. We observe no genome-wide association between differential DNA methylation and differential gene expression or splicing, though indirect relationships may exist between these factors. Our results provide insight into the gene regulatory context of DNA methylation in insects and the molecular avenues through which variation in maternal care influences developmental plasticity.

Clonal raider ant brain transcriptomics identifies candidate molecular mechanisms for reproductive division of labor

BACKGROUND

Division of labor between reproductive queens and workers that perform brood care is a hallmark of insect societies. However, studies of the molecular basis of this fundamental dichotomy are limited by the fact that the caste of an individual cannot typically be experimentally manipulated at the adult stage. Here we take advantage of the unique biology of the clonal raider ant, Ooceraea biroi, to study brain gene expression dynamics during experimentally induced transitions between reproductive and brood care behavior.

RESULTS

Introducing larvae that inhibit reproduction and induce brood care behavior causes much faster changes in adult gene expression than removing larvae. In addition, the general patterns of gene expression differ depending on whether ants transition from reproduction to brood care or vice versa, indicating that gene expression changes between phases are cyclic rather than pendular. Finally, we identify genes that could play upstream roles in regulating reproduction and behavior because they show large and early expression changes in one or both transitions.

CONCLUSIONS

Our analyses reveal that the nature and timing of gene expression changes differ substantially depending on the direction of the transition, and identify a suite of promising candidate molecular regulators of reproductive division of labor that can now be characterized further in both social and solitary animal models. This study contributes to understanding the molecular regulation of reproduction and behavior, as well as the organization and evolution of insect societies.

Back to the roots: the importance of using simple insect societies to understand the molecular basis of complex social life

The evolutionary trajectories toward insect eusociality come in two broad forms. In species like wasps, bees, and ants, the first helpers remained at the nest primarily to help with brood care. In species like aphids and termites, on the other hand, nest defense was initially the primary ecological driving force. To better understand the molecular basis of these two alternative evolutionary trajectories, it is therefore important to study the mechanistic basis of brood care and nest defense behavior. So far, most studies have compared morphologically distinct castes in advanced eusocial species of ants, bees, wasps, and termites. However, the interpretation of such comparisons is limited by multiple confounding factors and the fact that castes are typically fixed and cannot be manipulated at the adult stage. In this review, we argue that conducting molecular studies of brood care and nest defense in simpler, more flexible insect societies may complement studies of advanced eusocial insects and provide avenues toward more functional analyses. We review the available literature and propose candidate study systems for future molecular investigations of brood care and nest defense in social insects.

Do social insects support Haig's kin theory for the evolution of genomic imprinting?

Although numerous imprinted genes have been described in several lineages, the phenomenon of genomic imprinting presents a peculiar evolutionary problem. Several hypotheses have been proposed to explain gene imprinting, the most supported being Haig's kinship theory. This theory explains the observed pattern of imprinting and the resulting phenotypes as a competition for resources between related individuals, but despite its relevance it has not been independently tested. Haig's theory predicts that gene imprinting should be present in eusocial insects in many social scenarios. These lineages are therefore ideal for testing both the theory's predictions and the mechanism of gene imprinting. Here we review the behavioral evidence of genomic imprinting in eusocial insects, the evidence of a mechanism for genomic imprinting and finally we evaluate recent results showing parent of origin allele specific expression in honeybees in the light of Haig's theory.

Social regulation of insulin signaling and the evolution of eusociality in ants

Queens and workers of eusocial Hymenoptera are considered homologous to the reproductive and brood care phases of an ancestral subsocial life cycle. However, the molecular mechanisms underlying the evolution of reproductive division of labor remain obscure. Using a brain transcriptomics screen, we identified a single gene, insulin-like peptide 2 (ilp2), which is always up-regulated in ant reproductives, likely because they are better nourished than their nonreproductive nestmates. In clonal raider ants (Ooceraea biroi), larval signals inhibit adult reproduction by suppressing ilp2, thus producing a colony reproductive cycle reminiscent of ancestral subsociality. However, increasing ILP2 peptide levels overrides larval suppression, thereby breaking the colony cycle and inducing a stable division of labor. These findings suggest a simple model for the origin of ant eusociality via nutritionally determined reproductive asymmetries potentially amplified by larval signals.

Variation in DNA Methylation Is Not Consistently Reflected by Sociality in Hymenoptera

Changes in gene regulation that underlie phenotypic evolution can be encoded directly in the DNA sequence or mediated by chromatin modifications such as DNA methylation. It has been hypothesized that the evolution of eusocial division of labor is associated with enhanced gene regulatory potential, which may include expansions in DNA methylation in the genomes of Hymenoptera (bees, ants, wasps, and sawflies). Recently, this hypothesis garnered support from analyses of a commonly used metric to estimate DNA methylation in silico, CpG content. Here, we test this hypothesis using direct, nucleotide-level measures of DNA methylation across nine species of Hymenoptera. In doing so, we generated new DNA methylomes for three species of interest, including one solitary and one facultatively eusocial halictid bee and a sawfly. We demonstrate that the strength of correlation between CpG content and DNA methylation varies widely among hymenopteran taxa, highlighting shortcomings in the utility of CpG content as a proxy for DNA methylation in comparative studies of taxa with sparse DNA methylomes. We observed strikingly high levels of DNA methylation in the sawfly relative to other investigated hymenopterans. Analyses of molecular evolution suggest the relatively distinct sawfly DNA methylome may be associated with positive selection on functional DNMT3 domains. Sawflies are an outgroup to all ants, bees, and wasps, and no sawfly species are eusocial. We find no evidence that either global expansions or variation within individual ortholog groups in DNA methylation are consistently associated with the evolution of social behavior.

Genome-Wide Characterization of DNA Methylation in an Invasive Lepidopteran Pest, the Cotton Bollworm Helicoverpa armigera

The genes and genomes of insect pests are shaped by the wide array of selective forces encountered in their environments. While the molecular adaptations that evolve are beginning to be understood at the genomic and transcriptomic level, they have been less well characterized at an epigenetic level. Here, we present a genome-wide map of DNA methylation at single-nucleotide resolution for the cotton bollworm moth, Helicoverpa armigera, a globally invasive pest of agriculture. We show that methylation is almost identical in the larvae and adults of H. armigera and that, through whole-genome bisulfite sequencing (WGBS), at the most ∼0.9% of CpG sites in this species are methylated. We find that DNA methylation occurs primarily in exons, is positively correlated with gene expression, and that methylated genes are enriched for cellular "housekeeping" roles. H. armigera has an exceptional capacity for long-range migration. To explore the role of methylation in influencing the migratory phenotype of H. armigera, we performed targeted bisulfite sequencing on selected loci from 16 genes that were differentially expressed between adult moths exhibiting distinct flight performance in behavioral assays. While most CpG sites in these genes were not methylated between flight phenotypes, we identified hypermethylation in a demethylase (KDM4) that targets lysine-specific histone modifications, which are strongly associated with transcription and methylation. The H. armigera methylome provides new insights into the role of DNA methylation in a noctuid moth and is a valuable resource for further research into the epigenetic control of adaptive traits in this important pest.

Hemimetabolous genomes reveal molecular basis of termite eusociality

Around 150 million years ago, eusocial termites evolved from within the cockroaches, 50 million years before eusocial Hymenoptera, such as bees and ants, appeared. Here, we report the 2-Gb genome of the German cockroach, Blattella germanica, and the 1.3-Gb genome of the drywood termite Cryptotermes secundus. We show evolutionary signatures of termite eusociality by comparing the genomes and transcriptomes of three termites and the cockroach against the background of 16 other eusocial and non-eusocial insects. Dramatic adaptive changes in genes underlying the production and perception of pheromones confirm the importance of chemical communication in the termites. These are accompanied by major changes in gene regulation and the molecular evolution of caste determination. Many of these results parallel molecular mechanisms of eusocial evolution in Hymenoptera. However, the specific solutions are remarkably different, thus revealing a striking case of convergence in one of the major evolutionary transitions in biological complexity.

Developmental DNA methyltransferase expression in the fire ant Solenopsis invicta

DNA methylation is accomplished in animals by 2 classes of enzymes known as DNA methyltransferases, DNMT3 and DNMT1, which perform de novo methylation and maintenance methylation, respectively. Several studies of hymenopteran eusocial insects suggest that DNA methylation is capable of influencing developmental plasticity. However, fundamental questions remain about the patterning of DNA methylation during the course of insect development. In this study, we performed quantitative real-time PCR (qPCR) on transcripts from the single-copy orthologs of DNMT1 and DNMT3 in the red imported fire ant, Solenopsis invicta. In particular, we assessed the expression of S. invicta Dnmt1 and Dnmt3 mRNA during 7 stages of worker development, among behaviorally distinct adults, and among male and female gonads. Dnmt3 was most highly expressed during embryonic development, whereas Dnmt1 was similarly expressed throughout the course of development. Moreover, Dnmt1 and Dnmt3 were highly expressed in testes and ovaries. Neither Dnmt was significantly differentially expressed among heads of behaviorally distinct adult castes. Our results support the hypothesis that extensive patterning of DNA methylation occurs during gametogenesis and embryogenesis in the insect order Hymenoptera.

Differential gene expression is not required for facultative sex allocation: a transcriptome analysis of brain tissue in the parasitoid wasp Nasonia vitripennis

Whole-transcriptome technologies have been widely used in behavioural genetics to identify genes associated with the performance of a behaviour and provide clues to its mechanistic basis. Here, we consider the genetic basis of sex allocation behaviour in the parasitoid wasp Nasonia vitripennis. Female Nasonia facultatively vary their offspring sex ratio in line with Hamilton's theory of local mate competition (LMC). A single female or 'foundress' laying eggs on a patch will lay just enough sons to fertilize her daughters. As the number of 'foundresses' laying eggs on a patch increases (and LMC declines), females produce increasingly male-biased sex ratios. Phenotypic studies have revealed the cues females use to estimate the level of LMC their sons will experience, but our understanding of the genetics underlying sex allocation is limited. Here, we exposed females to three foundress number conditions, i.e. three LMC conditions, and allowed them to oviposit. mRNA was extracted from only the heads of these females to target the brain tissue. The subsequent RNA-seq experiment confirmed that differential gene expression is not associated with the response to sex allocation cues and that we must instead turn to the underlying neuroscience to reveal the underpinnings of this impressive behavioural plasticity.

Epigenetics and developmental plasticity in orthopteroid insects

Developmental plasticity is a key driver of the extraordinary ecological success of insects. Epigenetic mechanisms provide an important link between the external stimuli that initiate polyphenisms, and the stable changes in gene expression that govern alternative insect morphs. We review the epigenetics of orthopteroid insects, focussing on recent research on locusts and termites, two groups which display high levels of phenotypic plasticity, and for which genome sequences have become available in recent years. We examine research on the potential role of DNA methylation, histone modifications, and non-coding RNAs in the regulation of gene expression in these insects. DNA methylation patterns in orthopteroids share a number of characteristics with those of hymenopteran insects, although methylation levels are much higher, and extend to introns and repeat elements. Future examinations of epigenetic mechanisms in these insects will benefit from comparison of tissues from aged-matched individuals from alternative morphs, and adequate biological replication.

Genes and genomic processes underpinning the social lives of ants

The >15000 ant species are all highly social and show great variation in colony organization, complexity and behavior. The mechanisms by which such sociality evolved, as well as those underpinning the elaboration of ant societies since their ∼140 million year old common ancestor, have long been pondered. Here, we review recent insights generated using various genomic approaches. This includes understanding the molecular mechanisms underlying caste differentiation and the diversity of social structures, studying the impact of eusociality on genomic evolutionary rates, and investigating gene expression changes associated with differences in lifespan between castes. Furthermore, functional studies involving RNAi and CRISPR have recently been successfully applied to ants, opening the door to exciting research that promises to revolutionize the understanding of the evolution and diversification of social living.

Building a new research framework for social evolution: intralocus caste antagonism

The breeding and non‐breeding ‘castes’ of eusocial insects provide a striking example of role‐specific selection, where each caste maximises fitness through different morphological, behavioural and physiological trait values. Typically, queens are long‐lived egg‐layers, while workers are short‐lived, largely sterile foragers. Remarkably, the two castes are nevertheless produced by the same genome. The existence of inter‐caste genetic correlations is a neglected consequence of this shared genome, potentially hindering the evolution of caste dimorphism: alleles that increase the productivity of queens may decrease the productivity of workers and vice versa, such that each caste is prevented from reaching optimal trait values. A likely consequence of this ‘intralocus caste antagonism’ should be the maintenance of genetic variation for fitness and maladaptation within castes (termed ‘caste load’), analogous to the result of intralocus sexual antagonism. The aim of this review is to create a research framework for understanding caste antagonism, drawing in part upon conceptual similarities with sexual antagonism. By reviewing both the social insect and sexual antagonism literature, we highlight the current empirical evidence for caste antagonism, discuss social systems of interest, how antagonism might be resolved, and challenges for future research. We also introduce the idea that sexual and caste antagonism could interact, creating a three‐way antagonism over gene expression. This includes unpacking the implications of haplodiploidy for the outcome of this complex interaction.

Origins of Aminergic Regulation of Behavior in Complex Insect Social Systems

Neuromodulators are conserved across insect taxa, but how biogenic amines and their receptors in ancestral solitary forms have been co-opted to control behaviors in derived socially complex species is largely unknown. Here we explore patterns associated with the functions of octopamine (OA), serotonin (5-HT) and dopamine (DA) in solitary ancestral insects and their derived functions in eusocial ants, bees, wasps and termites. Synthesizing current findings that reveal potential ancestral roles of monoamines in insects, we identify physiological processes and conserved behaviors under aminergic control, consider how biogenic amines may have evolved to modulate complex social behavior, and present focal research areas that warrant further study.

Genetics and Evolution of Social Behavior in Insects

The study of insect social behavior has offered tremendous insight into the molecular mechanisms mediating behavioral and phenotypic plasticity. Genomic applications to the study of eusocial insect species, in particular, have led to several hypotheses for the processes underlying the molecular evolution of behavior. Advances in understanding the genetic control of social organization have also been made, suggesting an important role for supergenes in the evolution of divergent behavioral phenotypes. Intensive study of social phenotypes across species has revealed that behavior and caste are controlled by an interaction between genetic and environmentally mediated effects and, further, that gene expression and regulation mediate plastic responses to environmental signals. However, several key methodological flaws that are hindering progress in the study of insect social behavior remain. After reviewing the current state of knowledge, we outline ongoing challenges in experimental design that remain to be overcome in order to advance the field.

Maximizing ecological and evolutionary insight in bisulfite sequencing data sets

Genome-scale bisulfite sequencing approaches have opened the door to ecological and evolutionary studies of DNA methylation in many organisms. These approaches can be powerful. However, they introduce new methodological and statistical considerations, some of which are particularly relevant to non-model systems. Here, we highlight how these considerations influence a study's power to link methylation variation with a predictor variable of interest. Relative to current practice, we argue that sample sizes will need to increase to provide robust insights. We also provide recommendations for overcoming common challenges and an R Shiny app to aid in study design.

Sociomics: Using Omic Approaches to Understand Social Evolution

All of life is social, from genes cooperating to form organisms, to animals cooperating to form societies. Omic approaches offer exceptional opportunities to solve major outstanding problems in the study of how sociality evolves. First, omics can be used to clarify the extent and form of sociality in natural populations. This is especially useful in species where it is difficult to study social traits in natural populations, such as bacteria and other microbes. Second, omics can be used to examine the consequences of sociality for genome evolution and gene expression. This is especially useful in cases where there is clear variation in the level of sociality, such as the social insects. Major tasks for the future are to apply these approaches to a wider range of non-model organisms, and to move from exploratory analyses to the testing of evolutionary theory.

The adaptive significance of phasic colony cycles in army ants

Army ants are top arthropod predators in tropical forests around the world. The colonies of many army ant species undergo stereotypical behavioral and reproductive cycles, alternating between brood care and reproductive phases. In the brood care phase, colonies contain a cohort of larvae that are synchronized in their development and have to be fed. In the reproductive phase larvae are absent and oviposition takes place. Despite these colony cycles being a striking feature of army ant biology, their adaptive significance is unclear. Here we use a modeling approach to show that cyclic reproduction is favored under conditions where per capita foraging costs decrease with the number of larvae in a colony ("High Cost of Entry" scenario), while continuous reproduction is favored under conditions where per capita foraging costs increase with the number of larvae ("Resource Exhaustion" scenario). We argue that the former scenario specifically applies to army ants, because large raiding parties are required to overpower prey colonies. However, once raiding is successful it provides abundant food for a large cohort of larvae. The latter scenario, on the other hand, will apply to non-army ants, because in those species local resource depletion will force workers to forage over larger distances to feed large larval cohorts. Our model provides a quantitative framework for understanding the adaptive value of phasic colony cycles in ants.

Evolution of DNA Methylation across Insects

DNA methylation contributes to gene and transcriptional regulation in eukaryotes, and therefore has been hypothesized to facilitate the evolution of plastic traits such as sociality in insects. However, DNA methylation is sparsely studied in insects. Therefore, we documented patterns of DNA methylation across a wide diversity of insects. We predicted that underlying enzymatic machinery is concordant with patterns of DNA methylation. Finally, given the suggestion that DNA methylation facilitated social evolution in Hymenoptera, we tested the hypothesis that the DNA methylation system will be associated with presence/absence of sociality among other insect orders. We found DNA methylation to be widespread, detected in all orders examined except Diptera (flies). Whole genome bisulfite sequencing showed that orders differed in levels of DNA methylation. Hymenopteran (ants, bees, wasps and sawflies) had some of the lowest levels, including several potential losses. Blattodea (cockroaches and termites) show all possible patterns, including a potential loss of DNA methylation in a eusocial species whereas solitary species had the highest levels. Species with DNA methylation do not always possess the typical enzymatic machinery. We identified a gene duplication event in the maintenance DNA methyltransferase 1 (DNMT1) that is shared by some Hymenoptera, and paralogs have experienced divergent, nonneutral evolution. This diversity and nonneutral evolution of underlying machinery suggests alternative DNA methylation pathways may exist. Phylogenetically corrected comparisons revealed no evidence that supports evolutionary association between sociality and DNA methylation. Future functional studies will be required to advance our understanding of DNA methylation in insects.

Molecular Evolution of Insect Sociality: An Eco-Evo-Devo Perspective

The evolution of eusociality is a perennial issue in evolutionary biology, and genomic advances have fueled steadily growing interest in the genetic changes underlying social evolution. Along with a recent flurry of research on comparative and evolutionary genomics in different eusocial insect groups (bees, ants, wasps, and termites), several mechanistic explanations have emerged to describe the molecular evolution of eusociality from solitary behavior. These include solitary physiological ground plans, genetic toolkits of deeply conserved genes, evolutionary changes in protein-coding genes, cis regulation, and the structure of gene networks, epigenetics, and novel genes. Despite this proliferation of ideas, there has been little synthesis, even though these ideas are not mutually exclusive and may in fact be complementary. We review available data on molecular evolution of insect sociality and highlight key biotic and abiotic factors influencing social insect genomes. We then suggest both phylogenetic and ecological evolutionary developmental biology (eco-evo-devo) perspectives for a more synthetic view of molecular evolution in insect societies.

On the Molecular Basis of Division of Labor in Solenopsis invicta (Hymenoptera: Formicidae) Workers: RNA-seq Analysis

The fire ant Solenopsis invicta Buren is an important invasive pest. Among S. invicta workers behavioral changes depend on age where younger ants are nurses and older ants foragers. To identify potential genes associated with this division of labor, we compared gene expression between foragers and nurses by high-throughput sequencing. In total, we identified 1,618 genes significantly differently expressed between nurses and foragers, of which 542 were upregulated in foragers and 1,076 were upregulated in nurses. Several pathways related to metabolism were significantly enriched, such as lipid storage and fatty acid biosynthesis, which might contribute to the division of labor in S. invicta. Several genes involved in DNA methylation, transcription, and olfactory responses as well as resistance to stress were differentially expressed between nurses and foragers workers. Finally, a comparison between previously published microarray data and our RNA-seq data in S. invicta shows 116 genes overlap, and the GO term myofibril assembly (GO: 0030239) were simultaneously significantly enriched. These results advance knowledge of potentially important genes and molecular pathways associated with worker division of labor in S. invicta. We hope our dataset will provide . candidate target genes to disrupt organization in S. invicta as a control strategy against this invasive pest.

The ecology and evolution of social behavior in microbes

Cooperation has been studied extensively across the tree of life, from eusociality in insects to social behavior in humans, but it is only recently that a social dimension has been recognized and extensively explored for microbes. Research into microbial cooperation has accelerated dramatically and microbes have become a favorite system because of their fast evolution, their convenience as lab study systems and the opportunity for molecular investigations. However, the study of microbes also poses significant challenges, such as a lack of knowledge and an inaccessibility of the ecological context (used here to include both the abiotic and the biotic environment) under which the trait deemed cooperative has evolved and is maintained. I review the experimental and theoretical evidence in support of the limitations of the study of social behavior in microbes in the absence of an ecological context. I discuss both the need and the opportunities for experimental investigations that can inform a theoretical framework able to reframe the general questions of social behavior in a clear ecological context and to account for eco-evolutionary feedback.

The caste- and sex-specific DNA methylome of the termite Zootermopsis nevadensis

Epigenetic inheritance plays an important role in mediating alternative phenotype in highly social species. In order to gain a greater understanding of epigenetic effects in societies, we investigated DNA methylation in the termite Zootermopsis nevadensis. Termites are the most ancient social insects, and developmentally distinct from highly-studied, hymenopteran social insects. We used replicated bisulfite-sequencing to investigate patterns of DNA methylation in both sexes and among castes of Z. nevadensis. We discovered that Z. nevadensis displayed some of the highest levels of DNA methylation found in insects. We also found strong differences in methylation between castes. Methylated genes tended to be uniformly and highly expressed demonstrating the antiquity of associations between intragenic methylation and gene expression. Differentially methylated genes were more likely to be alternatively spliced than not differentially methylated genes, and possessed considerable enrichment for development-associated functions. We further observed strong overrepresentation of multiple transcription factor binding sites and miRNA profiles associated with differential methylation, providing new insights into the possible function of DNA methylation. Overall, our results show that DNA methylation is widespread and associated with caste differences in termites. More generally, this study provides insights into the function of DNA methylation and the success of insect societies.

Genome methylation patterns across castes and generations in a parasitoid wasp

Environmental influences shape phenotypes within and across generations, often through DNA methylations that modify gene expression. Methylations were proposed to mediate caste and task allocation in some eusocial insects, but how an insect's environment affects DNA methylation in its offspring is yet unknown. We characterized parental effects on methylation profiles in the polyembryonic parasitoid wasp Copidosoma koehleri, as well as methylation patterns associated with its simple caste system. We used methylation‐sensitive amplified fragment length polymorphism (MS‐AFLP) to compare methylation patterns, among (1) reproductive and soldier larvae; and (2) offspring (larvae, pupae, and adults) of wasps that were reared at either high or low larval density and mated in the four possible combinations. Methylation frequencies were similar across castes, but the profiles of methylated fragments differed significantly. Parental rearing density did not affect methylation frequencies in the offspring at any developmental stage. Principal coordinate analysis indicated no significant differences in methylation profiles among the four crossbreeding groups and the three developmental stages. Nevertheless, a clustering analysis, performed on a subset of the fragments, revealed similar methylation patterns in larvae, pupae, and adults in two of the four parental crosses. Nine fragments were methylated at two cytosine sites in all larvae, and five others were methylated at two sites in all adults. Thus, DNA methylations correlate with within‐generation phenotypic plasticity due to caste. However, their association with developmental stage and with transgenerational epigenetic effects is not clearly supported.

The genomewide transcriptional response underlying the pea aphid wing polyphenism

Phenotypic plasticity is a key life history strategy used by many plants and animals living in heterogeneous environments. A multitude of studies have investigated the costs and limits of plasticity, as well as the conditions under which it evolves. Much less well understood are the molecular genetic mechanisms that enable an organism to sense its environment and respond in a plastic manner. The pea aphid wing polyphenism is a compelling laboratory model to study these mechanisms. In this polyphenism, environmental stressors like high density cause asexual, viviparous adult female aphids to change the development of their embryos from wingless to winged morphs. The life history trade-offs between the two morphs have been intensively studied, but the molecular mechanisms underlying this process remain largely unknown. We therefore performed a genomewide study of the maternal transcriptome at two time points with and without a crowding stress to discover the maternal molecular changes that lead to the development of winged vs. wingless offspring. We observed significant transcriptional changes in genes associated with odorant binding, neurotransmitter transport, hormonal activity and chromatin remodelling in the maternal transcriptome. We also found that titres of serotonin, dopamine and octopamine were higher in solitary compared to crowded aphids. We use these results to posit a model for how maternal signals inform a developing embryo to be winged or wingless. Our findings add significant insights into the identity of the molecular mechanisms that underlie environmentally induced morph determination and suggest a possible role for biogenic amine regulation in polyphenisms generally.

Chapter 6. Comparative genomic approaches to investigate molecular traits specific to social insects

Ageing is a feature of nearly all known organisms and, by its connection to survival, appears to trade off with fecundity. However, in some organisms such as in queens of social insects, this negative relation appears reversed and individuals live long and reproduce much. Since new experimental techniques, transcriptomes and genomes of many social insects have recently become available, a comparison of these data in a phylogenetic framework becomes feasible. This allows the study of general trends, species specific oddities and evolutionary dynamics of the molecular properties and changes which underlie ageing, fecundity and the reversal of this negative association. In the framework of social insect evolution, we review the most important recent insights, computational methods, their applications and data resources which are available.

Allele-Specific Transcriptome and Methylome Analysis Reveals Stable Inheritance and Cis-Regulation of DNA Methylation in Nasonia

Gene expression divergence between closely related species could be attributed to both cis- and trans- DNA sequence changes during evolution, but it is unclear how the evolutionary dynamics of epigenetic marks are regulated. In eutherian mammals, biparental DNA methylation marks are erased and reset during gametogenesis, resulting in paternal or maternal imprints, which lead to genomic imprinting. Whether DNA methylation reprogramming exists in insects is not known. Wasps of the genus Nasonia are non-social parasitoids that are emerging as a model for studies of epigenetic processes in insects. In this study, we quantified allele-specific expression and methylation genome-wide in Nasonia vitripennis and Nasonia giraulti and their reciprocal F1 hybrids. No parent-of-origin effect in allelic expression was found for >8,000 covered genes, suggesting a lack of genomic imprinting in adult Nasonia. As we expected, both significant cis- and trans- effects are responsible for the expression divergence between N. vitripennis and N. giraulti. Surprisingly, all 178 differentially methylated genes are also differentially methylated between the two alleles in F1 hybrid offspring, recapitulating the parental methylation status with nearly 100% fidelity, indicating the presence of strong cis-elements driving the target of gene body methylation. In addition, we discovered that total and allele-specific expression are positively correlated with allele-specific methylation in a subset of the differentially methylated genes. The 100% cis-regulation in F1 hybrids suggests the methylation machinery is conserved and DNA methylation is targeted by cis features in Nasonia. The lack of genomic imprinting and parent-of-origin differentially methylated regions in Nasonia, together with the stable inheritance of methylation status between generations, suggests either a cis-regulatory motif for methylation at the DNA level or highly stable inheritance of an epigenetic signal in Nasonia.

RNA-sequencing and whole-genome bisulfite sequencing in the hybrid offspring of two Nasonia parasitoid wasp species revealed strong cis-regulation of methylation and allele-specific expression. No gene was found to display genomic imprinting.

The relationship between methylation of genomic DNA and expression of the genes that it encodes—and how this relationship changes during evolution—has been widely studied in mammals, but remains less well understood for insects. Here we analyze the expressed mRNA transcripts and genomic DNA methylation of the hybrid offspring of a pair of Nasonia parasitoid wasp species, producing a wealth of information about the regulation of gene expression. We find that variation in DNA sequence impacts expression on the same strand (called “cis-regulation”), and that cytosine methylation state is also associated in cis with the regulatory consequences of this base alteration. We show that these wasp species lack differential expression dependent on parent-of-origin (called “genomic imprinting”), and that in the hybrids the alleles retain the methylation status of the parental species in a strong cis-regulated fashion. Transcript abundances were also largely driven in a cis-regulated manner, consistent with a correlation between methylation status and expression levels. Despite the many differences between Nasonia and mammals in the impact of genomic DNA methylation, in both groups the use of methylated cytosine has been co-opted in ways that help tune gene expression.

Ant Genetics: Reproductive Physiology, Worker Morphology, and Behavior

Many exciting studies have begun to elucidate the genetics of the morphological and physiological diversity of ants, but as yet few studies have investigated the genetics of ant behavior directly. Ant genomes are marked by extreme rates of gene turnover, especially in gene families related to olfactory communication, such as the synthesis of cuticular hydrocarbons and the perception of environmental semiochemicals. Transcriptomic and epigenetic differences are apparent between reproductive and sterile females, males and females, and workers that differ in body size. Quantitative genetic approaches suggest heritability of task performance, and population genetic studies indicate a genetic association with reproductive status in some species. Gene expression is associated with behavior including foraging, response to queens attempting to join a colony, circadian patterns of task performance, and age-related changes of task. Ant behavioral genetics needs further investigation of the feedback between individual-level physiological changes and socially mediated responses to environmental conditions.

Parent-of-origin effects on genome-wide DNA methylation in the Cape honey bee (Apis mellifera capensis) may be confounded by allele-specific methylation

BACKGROUND

Intersexual genomic conflict sometimes leads to unequal expression of paternal and maternal alleles in offspring, resulting in parent-of-origin effects. In honey bees reciprocal crosses can show strong parent-of-origin effects, supporting theoretical predictions that genomic imprinting occurs in this species. Mechanisms behind imprinting in honey bees are unclear but differential DNA methylation in eggs and sperm suggests that DNA methylation could be involved. Nonetheless, because DNA methylation is multifunctional, it is difficult to separate imprinting from other roles of methylation. Here we use a novel approach to investigate parent-of-origin DNA methylation in honey bees. In the subspecies Apis mellifera capensis, reproduction of females occurs either sexually by fertilization of eggs with sperm, or via thelytokous parthenogenesis, producing female embryos derived from two maternal genomes.

RESULTS

We compared genome-wide methylation patterns of sexually-produced, diploid embryos laid by a queen, with parthenogenetically-produced diploid embryos laid by her daughters. Thelytokous embryos inheriting two maternal genomes had fewer hypermethylated genes compared to fertilized embryos, supporting the prediction that fertilized embryos have increased methylation due to inheritance of a paternal genome. However, bisulfite PCR and sequencing of a differentially methylated gene, Stan (GB18207) showed strong allele-specific methylation that was maintained in both fertilized and thelytokous embryos. For this gene, methylation was associated with haplotype, not parent of origin.

CONCLUSIONS

The results of our study are consistent with predictions from the kin theory of genomic imprinting. However, our demonstration of allele-specific methylation based on sequence shows that genome-wide differential methylation studies can potentially confound imprinting and allele-specific methylation. It further suggests that methylation patterns are heritable or that specific sequence motifs are targets for methylation in some genes.