Epigenetics in social insects: a new direction for understanding the evolution of castes

Genome-wide DNA methylation patterns in bumble bee (Bombus vosnesenskii) populations from spatial-environmental range extremes

Unraveling molecular mechanisms of adaptation to complex environments is crucial to understanding tolerance of abiotic pressures and responses to climatic change. Epigenetic variation is increasingly recognized as a mechanism that can facilitate rapid responses to changing environmental cues. To investigate variation in genetic and epigenetic diversity at spatial and thermal extremes, we use whole genome and methylome sequencing to generate a high-resolution map of DNA methylation in the bumble bee Bombus vosnesenskii. We sample two populations representing spatial and environmental range extremes (a warm southern low-elevation site and a cold northern high-elevation site) previously shown to exhibit differences in thermal tolerance and determine positions in the genome that are consistently and variably methylated across samples. Bisulfite sequencing reveals methylation characteristics similar to other arthropods, with low global CpG methylation but high methylation concentrated in gene bodies and in genome regions with low nucleotide diversity. Differentially methylated sites (n = 2066) were largely hypomethylated in the northern high-elevation population but not related to local sequence differentiation. The concentration of methylated and differentially methylated sites in exons and putative promoter regions suggests a possible role in gene regulation, and this high-resolution analysis of intraspecific epigenetic variation in wild Bombus suggests that the function of methylation in niche adaptation would be worth further investigation.

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.

Genome-Wide Identification and Expression Pattern of Cytochrome P450 Genes in the Social Aphid Pseudoregma bambucicola

Cytochrome P450 monooxygenases (P450s) have a variety of functions, including involvement in the metabolism of exogenous substances and the synthesis and degradation of endogenous substances, which are important for the growth and development of insects. Pseudoregma bambucicola is a social aphid that produces genetically identical but morphologically and behaviorally distinct first-instar soldiers and normal nymphs within colonies. In this study, we identified 43 P450 genes based on P. bambucicola genome data. Phylogenetic analysis showed that these genes were classified into 4 clans, 13 families, and 23 subfamilies. The CYP3 and CYP4 clans had a somewhat decreased number of genes. In addition, differential gene expression analysis based on transcriptome data showed that several P450 genes, including CYP18A1, CYP4G332, and CYP4G333, showed higher expression levels in soldiers compared to normal nymphs and adult aphids. These genes may be candidates for causing epidermal hardening and developmental arrest in soldiers. This study provides valuable data and lays the foundation for the study of functions of P450 genes in the social aphid P. bambucicola.

Gut microbiota is a potential factor in shaping phenotypic variation in larvae and adults of female bumble bees

Host symbionts are often considered an essential part of the host phenotype, influencing host growth and development. Bumble bee is an ideal model for investigating the relationship between microbiota and phenotypes. Variations in life history across bumble bees may influence the community composition of gut microbiota, which in turn influences phenotypes. In this study, we explored gut microbiota from four development stages (early-instar larvae, 1st instar; mid-instar larvae, 5th instar; late-instar larvae, 9th instar; and adults) of workers and queens in the bumble bee Bombus terrestris using the full-length 16S rRNA sequencing technology. The results showed that morphological indices (weight and head capsule) were significantly different between workers and queens from 5th instar larvae (p < 0.01). The alpha and beta diversities of gut microbiota were similar between workers and queens in two groups: early instar and mid instar larvae. However, the alpha diversity was significantly different in late instar larvae or adults. The relative abundance of three main phyla of bacteria (Cyanobacteria, Proteobacteria, and Firmicutes) and two genera (Snodgrassella and Lactobacillus) were significantly different (p < 0.01) between workers and queens in late instar larvae or adults. Also, we found that age significantly affected the microbial alpha diversity as the Shannon and ASVs indices differed significantly among the four development stages. Our study suggests that the 5th instar larval stage can be used to judge the morphology of workers or queens in bumble bees. The key microbes differing in phenotypes may be involved in regulating phenotypic variations.

Computational modelling of plasticity-led evolution

Plasticity-led evolution is a form of evolution where a change in the environment induces novel traits via phenotypic plasticity, after which the novel traits are genetically accommodated over generations under the novel environment. This mode of evolution is expected to resolve the problem of gradualism (i.e., evolution by the slow accumulation of mutations that induce phenotypic variation) implied by the Modern Evolutionary Synthesis, in the face of a large environmental change. While experimental works are essential for validating that plasticity-led evolution indeed happened, we need computational models to gain insight into its underlying mechanisms and make qualitative predictions. Such computational models should include the developmental process and gene-environment interactions in addition to genetics and natural selection. We point out that gene regulatory network models can incorporate all the above notions. In this review, we highlight results from computational modelling of gene regulatory networks that consolidate the criteria of plasticity-led evolution. Since gene regulatory networks are mathematically equivalent to artificial recurrent neural networks, we also discuss their analogies and discrepancies, which may help further understand the mechanisms underlying plasticity-led evolution.

Association between parasite microbiomes and caste development and colony structure in a social trematode

Division of labour through the formation of morphologically and functionally distinct castes is a recurring theme in the evolution of animal sociality. The mechanisms driving the differentiation of individuals into distinct castes remain poorly understood, especially for animals forming clonal colonies. We test the association between microbiomes and caste formation within the social trematode Philophthalmus attenuatus, using a metabarcoding approach targeting the bacterial 16S SSU rRNA gene. Clonal colonies of this trematode within snail hosts comprise large reproductive individuals which produce dispersal stages, and small, non-reproducing soldiers which defend the colony against invaders. In colonies extracted directly from field-collected snails, reproductives harboured more diverse bacterial communities than soldiers, and reproductives and soldiers harboured distinct bacterial communities, at all taxonomic levels considered. No single bacterial taxon showed high enough prevalence in either soldiers or reproductives to be singled out as a key driver, indicating that the whole microbial community contributes to these differences. Other colonies were experimentally exposed to antibiotics to alter their bacterial communities, and sampled shortly after treatment and weeks later after allowing for turnover of colony members. At those time points, bacterial communities of the two castes still differed across all antibiotic treatments; however, the caste ratio within colonies changed: after antibiotic disruption and turnover of individuals, new individuals were more likely to become reproductives than in undisturbed control colonies. Our results reveal that each caste has a distinct microbiome; whether the social context affects the microbiota, or whether microbes contribute to modulating the phenotype of individuals, remains to be determined.

Dynamic changes in genomic 5-hydroxymethyluracil and N6-methyladenine levels in the Drosophila melanogaster life cycle and in response to different temperature conditions

In this study, the level of DNA modifications was investigated in three developmental stages of Drosophila melanogaster (larvae, pupae, imago) and in an in vitro model (Schneider 2 cells). Analysis was carried out using two-dimensional ultra-performance liquid chromatography with tandem mass spectrometry. Our method made it possible, for the first time, to analyze a broad spectrum of DNA modifications in the three stages of Drosophila. Each stage was characterized by a specific modification pattern, and the levels of these compounds fluctuated throughout the D. melanogaster life cycle. The level of DNA modification was also compared between insects bred at 25 °C (optimal temperature) and at 18 °C, and the groups differed significantly. The profound changes in N6-methyladenine and 5-hydroxymethyluracil levels during the Drosophila life cycle and as a result of breeding temperature changes indicate that these DNA modifications can play important regulatory roles in response to environmental changes and/or biological conditions. Moreover, the supplementation of Schneider 2 cells with 1 mM L-ascorbic acid caused a time-dependent increase in the level of 5-(hydroxymethyl)-2'-deoxyuridine. These data suggest that a certain pool of this compound may arise from the enzymatic activity of the dTET protein.

Conservation Biology and Reproduction in a Time of Developmental Plasticity

The objective of this review is to ask whether, and how, principles in conservation biology may need to be revisited in light of new knowledge about the power of epigenetics to alter developmental pathways. Importantly, conservation breeding programmes, used widely by zoological parks and aquariums, may appear in some cases to reduce fitness by decreasing animals' abilities to cope when confronted with the 'wild side' of their natural habitats. Would less comfortable captive conditions lead to the selection of individuals that, despite being adapted to life in a captive environment, be better able to thrive if relocated to a more natural environment? While threatened populations may benefit from advanced reproductive technologies, these may actually induce undesirable epigenetic changes. Thus, there may be inherent risks to the health and welfare of offspring (as is suspected in humans). Advanced breeding technologies, especially those that aim to regenerate the rarest species using stem cell reprogramming and artificial gametes, may also lead to unwanted epigenetic modifications. Current knowledge is still incomplete, and therefore ethical decisions about novel breeding methods remain controversial and difficult to resolve.

DNA methylation across the tree of life, from micro to macro-organism

DNA methylation is a process in which methyl (CH3) groups are added to the DNA molecule. The DNA segment does not change in the sequence, but DNA methylation could alter the action of DNA. Different enzymes like DNA methyltransferases (DNMTs) take part in methylation of cytosine/adenine nucleosides in DNA. In prokaryotes, DNA methylation is performed to prevent the attack of phage and also plays a role in the chromosome replication and repair. In fungi, DNA methylation is studied to see the transcriptional changes, as in insects, the DNA methylation is not that well-known, it plays a different role like other organisms. In mammals, the DNA methylation is related to different types of cancers and plays the most important role in the placental development and abnormal DNA methylation connected with diseases like cancer, autoimmune diseases, and rheumatoid arthritis.

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.

Parent-of-origin effects, allele-specific expression, genomic imprinting and paternal manipulation in social insects

Haplo-diploidy and the relatedness asymmetries it generates mean that social insects are prime candidates for the evolution of genomic imprinting. In single-mating social insect species, some genes may be selected to evolve genomic mechanisms that enhance reproduction by workers when they are inherited from a female. This situation reverses in multiple mating species, where genes inherited from fathers can be under selection to enhance the reproductive success of daughters. Reciprocal crosses between subspecies of honeybees have shown strong parent-of-origin effects on worker reproductive phenotypes, and this could be evidence of such genomic imprinting affecting genes related to worker reproduction. It is also possible that social insect fathers directly affect gene expression in their daughters, for example, by placing small interfering RNA molecules in semen. Gene expression studies have repeatedly found evidence of parent-specific gene expression in social insects, but it is unclear at this time whether this arises from genomic imprinting, paternal manipulation, an artefact of cyto-nuclear interactions, or all of these. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

High-quality chromosome-level genome assembly and full-length transcriptome analysis of the pharaoh ant Monomorium pharaonis

Background

Ants with complex societies have fascinated scientists for centuries. Comparative genomic and transcriptomic analyses across ant species and castes have revealed important insights into the molecular mechanisms underlying ant caste differentiation. However, most current ant genomes and transcriptomes are highly fragmented and incomplete, which hinders our understanding of the molecular basis for complex ant societies.

Findings

By hybridizing Illumina, Pacific Biosciences, and Hi-C sequencing technologies, we de novo assembled a chromosome-level genome for Monomorium pharaonis, with a scaffold N50 of 27.2 Mb. Our new assembly provides better resolution for the discovery of genome rearrangement events at the chromosome level. Analysis of full-length isoform sequencing (ISO-seq) suggested that ∼15 Gb of ISO-seq data were sufficient to cover most expressed genes, but the number of transcript isoforms steadily increased with sequencing data coverage. Our high-depth ISO-seq data greatly improved the quality of gene annotation and enabled the accurate detection of alternative splicing isoforms in different castes of M. pharaonis. Comparative transcriptome analysis across castes based on the ISO-seq data revealed an unprecedented number of transcript isoforms, including many caste-specific isoforms. We also identified a number of conserved long non-coding RNAs that evolved specifically in ant lineages and several that were conserved across insect lineages.

Conclusions

We produced a high-quality chromosome-level genome for M. pharaonis, which significantly improved previous short-read assemblies. Together with full-length transcriptomes for all castes, we generated a highly accurate annotation for this ant species. These long-read sequencing results provide a useful resource for future functional studies on the genetic mechanisms underlying the evolution of social behaviors and organization in ants.

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.

Evolutionary relationship of fat body endoreduplication and queen fecundity in termites

Endoreduplication or nuclear genome replication without cell division is widely observed in the metabolically active tissues of plants and animals. The fat body cells of adult female insects produce abundant yolk proteins and become polyploid, which is assumed to accelerate egg production. Recently, it was reported that in termites, endopolyploidy in the fat body occurs only in queens but not in the other females; however, the relationship between the fecundity and ploidy level in the fat body remains unclear. Termite queens exhibit a huge variation in their egg-producing capacity among different species; queens in the species with a foraging lifestyle, in which workers leave the nest to forage outside, are much more fecund than those in the species living in a single piece of wood. In this study, we conducted ploidy analyses on three foragings and three wood-dwelling termites via flow cytometry. In all the species, the fat body of queens contained significantly more polyploid cells than that of other nonreproductive females, considering their body size effect. However, the male fat body, which is not involved in yolk production, did not show consistency in polyploid cell numbers among the species studied. Moreover, highly fecund queens in foraging termites exhibit higher levels of endopolyploidy in their fat body than those with less fecundity in wood-dwelling termites. These results suggest that endopolyploidy in the fat body of termite queens can boost their egg production, and the level of endopolyploidy in their fat body is linked to their fecundity. Our study provides a novel insight into the evolutionary relationship between endoreduplication and caste specialization in social insects.

Unique DNA Methylation Profiles Are Associated with cis-Variation in Honey Bees

DNA methylation is an important epigenetic modification that mediates diverse processes such as cellular differentiation, phenotypic plasticity, and genomic imprinting. Mounting evidence suggests that local DNA sequence variation can be associated with particular DNA methylation states, indicating that the interplay between genetic and epigenetic factors may contribute synergistically to the phenotypic complexity of organisms. Social insects such as ants, bees, and wasps have extensive phenotypic plasticity manifested in their different castes, and this plasticity has been associated with variation in DNA methylation. Yet, the influence of genetic variation on DNA methylation state remains mostly unknown. Here we examine the importance of sequence-specific methylation at the genome-wide level, using whole-genome bisulfite sequencing of the semen of individual honey bee males. We find that individual males harbor unique DNA methylation patterns in their semen, and that genes that are more variable at the epigenetic level are also more likely to be variable at the genetic level. DNA sequence variation can affect DNA methylation by modifying CG sites directly, but can also be associated with local variation in cis that is not CG-site specific. We show that covariation in sequence polymorphism and DNA methylation state contributes to the individual-specificity of epigenetic marks in social insects, which likely promotes their retention across generations, and their capacity to influence evolutionary adaptation.

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.

Evolution, phylogenetic distribution and functional ecology of division of labour in trematodes

Division of labour has evolved in many social animals where colonies consist of clones or close kin. It involves the performance of different tasks by morphologically distinct castes, leading to increased colony fitness. Recently, a form of division of labour has been discovered in trematodes: clonal rediae inside the snail intermediate host belong either to a large-bodied reproductive caste, or to a much smaller and morphologically distinct ‘soldier’ caste which defends the colony against co-infecting trematodes. We review recent research on this phenomenon, focusing on its phylogenetic distribution, its possible evolutionary origins, and how division of labour functions to allow trematode colonies within their snail host to adjust to threats and changing conditions. To date, division of labour has been documented in 15 species from three families: Himasthlidae, Philophthalmidae and Heterophyidae. Although this list of species is certainly incomplete, the evidence suggests that division of labour has arisen independently more than once in the evolutionary history of trematodes. We propose a simple scenario for the gradual evolution of division of labour in trematodes facing a high risk of competition in a long-lived snail host. Starting with initial conditions prior to the origin of castes (size variation among rediae within a colony, size-dependent production of cercariae by rediae, and a trade-off between cercarial production and other functions, such as defence), maximising colony fitness (R₀) can lead to caste formation or the age-structured division of labour observed in some trematodes. Finally, we summarise recent research showing that caste ratios, i.e. relative numbers of reproductive and soldier rediae per colony, become more soldier-biased in colonies exposed to competition from another trematode species sharing the same snail, and also respond to other stressors threatening the host’s survival or the colony itself. In addition, there is evidence of asymmetrical phenotypic plasticity among individual caste members: reproductives can assume defensive functions against competitors in the absence of soldiers, whereas soldiers are incapable of growing into reproductives if the latter’s numbers are reduced. We conclude by highlighting future research directions, and the advantages of trematodes as model systems to study social evolution.

Caste-Specific Expression of Na+/K+-ATPase in the Asian Weaver Ant, Oecophylla smaragdina (Fabricius, 1775)

Social insect colonies adopt different levels of survival strategies and exhibit well-defined reproductive division of labour. Oecophylla smaragdina (Fabricius, 1775) has physiological and behavioral adaptations that enable them to forage at extreme environmental conditions and are lethal to most other insects. Ion homeostasis is the key process in an organism's survival mechanism. Among ion pumps, the ATP-dependent sodium-potassium ion pump is essential for maintaining the Na⁺ and K⁺ ionic balance and is well known as the primary consumer of energy. Oecophylla smaragdina plays pivotal role as a model among social insects for understanding ion homeostasis at the organization level of the castes. We have evaluated the expression and activity of Na⁺/K⁺-ATPase among various castes of O. smaragdina (worker subcastes, queen and male). Real-time PCR and immunoblotting analyses revealed the differential expression of Na⁺/K⁺-ATPase in the castes. Significantly higher expression of Na⁺/K⁺-ATPase mRNA and protein were observed in the minor workers, queen, major workers and males respectively. These results suggest that in the weaver ant colony, the castes might have variously adapted and evolved with a well-developed ion transport mechanism which allows them to perform allocated tasks within the nest and could be a key to their adaptive benefits towards division of labour.

Divergent methylation pattern in adult stage between two forms of Tetranychus urticae (Acari: Tetranychidae)

The two-spotted spider mite, Tetranychus urticae Koch has two forms: green form and red form. Understanding the molecular basis of how these two forms established without divergent genetic background is an intriguing area. As a well-known epigenetic process, DNA methylation has particularly important roles in gene regulation and developmental variation across diverse organisms that do not alter genetic background. Here, to investigate whether DNA methylation could be associated with different phenotypic consequences in the two forms of T. urticae, we surveyed the genome-wide cytosine methylation status and expression level of DNA methyltransferase 3 (Tudnmt3) throughout their entire life cycle. Methylation-sensitive amplification polymorphism (MSAP) analyses of 585 loci revealed variable methylation patterns in the different developmental stages. In particular, principal coordinates analysis (PCoA) indicates a significant epigenetic differentiation between female adults of the two forms. The gene expression of Tudnmt3 was detected in all examined developmental stages, which was significantly different in the adult stage of the two forms. Together, our results reveal the epigenetic distance between the two forms of T. urticae, suggesting that DNA methylation might be implicated in different developmental demands, and contribute to different phenotypes in the adult stage of these two forms.

Honey bees consider larval nutritional status rather than genetic relatedness when selecting larvae for emergency queen rearing

In honey bees and many other social insects, production of queens is a vital task, as colony fitness is dependent on queens. The factors considered by honey bee workers in selecting larvae to rear new queens during emergency queen rearing are poorly understood. Identifying these parameters is critical, both in an evolutionary and apicultural context. As female caste development in honey bees is dependent on larval diet (i.e. nutrition), we hypothesized that larval nutritional state is meticulously assessed and used by workers in selection of larvae for queen rearing. To test this hypothesis, we conducted a series of experiments manipulating the nutritional status of one day old larvae by depriving them of brood food for a four-hour period, and then allowing workers to choose larvae for rearing queens from nutritionally deprived and non-deprived larvae. We simultaneously investigated the role of genetic relatedness in selection of larvae for queen rearing. In all the experiments, significantly greater numbers of non-deprived larvae than deprived larvae were selected for queen rearing irrespective of genetic relatedness. Our results demonstrate that honey bees perceive the nutritional state of larvae and use that information when selecting larvae for rearing queens in the natural emergency queen replacement process.

Optimists or realists? How ants allocate resources in making reproductive investments

Parents often face an investment trade-off between either producing many small or fewer large offspring. When environments vary predictably, the fittest parental solution matches available resources by varying only number of offspring and never optimal individual size. However when mismatches occur often between parental expectations and true resource levels, dynamic models like multifaceted parental investment (MFPI) and parental optimism (PO) both predict offspring size can vary significantly. MFPI is a "realist" strategy: parents assume future environments of average richness. When resources exceed expectations and it is too late to add more offspring, the best-case solution increases investment per individual. Brood size distributions therefore track the degree of mismatch from right-skewed around an optimal size (slight underestimation of resources) to left-skewed around a maximal size (gross underestimation). Conversely, PO is an "optimist" strategy: parents assume maximally good resource futures and match numbers to that situation. Normal or lean years do not affect "core" brood as costs primarily fall on excess "marginal" siblings who die or experience stunted growth (producing left-skewed distributions). Investment patterns supportive of both MFPI and PO models have been observed in nature, but studies that directly manipulate food resources to test predictions are lacking. Ant colonies produce many offspring per reproductive cycle and are amenable to experimental manipulation in ways that can differentiate between MFPI and PO investment strategies. Colonies in a natural population of a harvester ant (Pogonomyrmex salinus) were protein-supplemented over 2 years, and mature sexual offspring were collected annually prior to their nuptial flight. Several results support either MFPI or PO in terms of patterns in offspring size distributions and how protein differentially affected male and female production. Unpredicted by either model, however, is that supplementation affected distributions more strongly across years than within (e.g., small females are significantly rarer in the year after colonies receive protein). Parental investment strategies in P. salinus vary dynamically across years and conditions. Finding that past conditions can more strongly affect reproductive decisions than current ones, however, is not addressed by models of parental investment.

The transcriptomic changes associated with the development of social parasitism in the honeybee Apis mellifera capensis

Social insects are characterized by the division of labor. Queens usually dominate reproduction, whereas workers fulfill non-reproductive age-dependent tasks to maintain the colony. Although workers are typically sterile, they can activate their ovaries to produce their own offspring. In the extreme, worker reproduction can turn into social parasitism as in Apis mellifera capensis. These intraspecific parasites occupy a host colony, kill the resident queen, and take over the reproductive monopoly. Because they exhibit a queenlike behavior and are also treated like queens by the fellow workers, they are so-called pseudoqueens. Here, we compare the development of parasitic pseudoqueens and social workers at different time points using fat body transcriptome data. Two complementary analysis methods-a principal component analysis and a time course analysis-led to the identification of a core set of genes involved in the transition from a social worker into a highly fecund parasitic pseudoqueen. Comparing our results on pseudoqueens with gene expression data of honeybee queens revealed many similarities. In addition, there was a set of specific transcriptomic changes in the parasitic pseudoqueens that differed from both, queens and social workers, which may be typical for the development of the social parasitism in A. m. capensis.

Mechanism of evolution by genetic assimilation : Equivalence and independence of genetic mutation and epigenetic modulation in phenotypic expression

Conrad H. Waddington discovered the phenomenon of genetic assimilation through a series of experiments on fruit flies. In those experiments, artificially exerted environmental stress induced plastic phenotypic changes in the fruit flies, but after some generations, the same phenotypic variant started to appear without the environmental stress. Both the initial state (where the phenotypic changes were environmentally induced and plastic) and the final state (where the phenotypic changes were genetically fixed and constitutive) are experimental facts. However, it remains unclear how the environmentally induced phenotypic change in the first generation becomes genetically fixed in the central process of genetic assimilation itself. We have argued that the key to understanding the mechanism of genetic assimilation lies in epigenetics, and proposed the "cooperative model" in which the evolutionary process depends on both genetic and epigenetic factors. Evolutionary simulations based on the cooperative model reproduced the process of genetic assimilation. Detailed analysis of the trajectories has revealed genetic assimilation is a process in which epigenetically induced phenotypic changes are incrementally and statistically replaced with multiple minor genetic mutations through natural selection. In this scenario, epigenetic and genetic changes may be considered as mutually independent but equivalent in terms of their effects on phenotypic changes. This finding rejects the common (and confused) hypothesis that epigenetically induced phenotypic changes depend on genetic mutations. Furthermore, we argue that transgenerational epigenetic inheritance is not required for evolution by genetic assimilation.

Identification and caste-dependent expression patterns of DNA methylation associated genes in Bombus terrestris

DNA methylation has been proposed to play critical roles in caste fate and behavioral plasticity in bumblebees, however, there is little information on its regulatory mechanisms. Here, we identified six important genes mediating the modification of DNA methylation and determined their expression patterns in the bumblebee Bombus terrestris. There is a complete functional DNA methylation system, including four DNA methyltransferases (DNMT1a, DNMT1b, DNMT2, and DNMT3), a DNA demethylase (Ten-eleven translocation), and a methyl-CpG-binding domain protein in B. terrestris. Most of these genes were highly expressed in fat bodies and gonads but lowly expressed in antennae and brains of bumblebee adults. Besides, these genes exhibited caste-specific expression patterns in bumblebees, with higher transcription levels in queens than workers and drones. Whereas their expression levels showed no remarkable difference in queenright and queenless workers. These results suggested potential roles of DNA methylation-related genes in caste differentiation in bumblebees.

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.

Allele specific expression and methylation in the bumblebee, Bombus terrestris

The social hymenoptera are emerging as models for epigenetics. DNA methylation, the addition of a methyl group, is a common epigenetic marker. In mammals and flowering plants methylation affects allele specific expression. There is contradictory evidence for the role of methylation on allele specific expression in social insects. The aim of this paper is to investigate allele specific expression and monoallelic methylation in the bumblebee, Bombus terrestris. We found nineteen genes that were both monoallelically methylated and monoallelically expressed in a single bee. Fourteen of these genes express the hypermethylated allele, while the other five express the hypomethylated allele. We also searched for allele specific expression in twenty-nine published RNA-seq libraries. We found 555 loci with allele-specific expression. We discuss our results with reference to the functional role of methylation in gene expression in insects and in the as yet unquantified role of genetic cis effects in insect allele specific methylation and expression.

Insights into the biochemical defence and methylation of the solitary bee Osmia rufa L: A foundation for examining eusociality development

We examined age-related biochemical and histological changes in the fat bodies and hemolymph of Osmia rufa males and females. We analysed solitary bees during diapause, in October and in April; as well as the flying insects following diapause, in May and June. The trophocyte sizes, as well as the numbers of lipid droplets were the greatest at the beginning of diapause. Subsequently, they decreased along with age. Triglyceride and glucose concentrations systematically decreased in fat body cells but increased in the hemolymph from October to June. Concentrations/activities of (enzymatic and non-enzymatic) antioxidant and proteolytic systems, as well as phenoloxidase, aspartate aminotransferase, alanine aminotransferase and alkaline phosphatase levels were constant during the diapause, usually lower in the males than the females. Prior to the diapause/overwintering, the concentrations/activities of all the compounds were higher in the fat bodies than in the hemolymph. Later in the spring and in the summer, they increased in the hemolymph and on the body surfaces, while decreasing in the fat bodies. The global DNA methylation levels increased with age. Higher levels were always observed in the males than in the females. The study will promote better understanding of bee evolution and will be useful for the protection and management of solitary bees, with benefits to the environment and agriculture.

Genome, transcriptome and methylome sequencing of a primitively eusocial wasp reveal a greatly reduced DNA methylation system in a social insect

Comparative genomics of social insects has been intensely pursued in recent years with the goal of providing insights into the evolution of social behaviour and its underlying genomic and epigenomic basis. However, the comparative approach has been hampered by a paucity of data on some of the most informative social forms (e.g. incipiently and primitively social) and taxa (especially members of the wasp family Vespidae) for studying social evolution. Here, we provide a draft genome of the primitively eusocial model insect Polistes dominula, accompanied by analysis of caste-related transcriptome and methylome sequence data for adult queens and workers. Polistes dominula possesses a fairly typical hymenopteran genome, but shows very low genomewide GC content and some evidence of reduced genome size. We found numerous caste-related differences in gene expression, with evidence that both conserved and novel genes are related to caste differences. Most strikingly, these -omics data reveal a major reduction in one of the major epigenetic mechanisms that has been previously suggested to be important for caste differences in social insects: DNA methylation. Along with a conspicuous loss of a key gene associated with environmentally responsive DNA methylation (the de novo DNA methyltransferase Dnmt3), these wasps have greatly reduced genomewide methylation to almost zero. In addition to providing a valuable resource for comparative analysis of social insect evolution, our integrative -omics data for this important behavioural and evolutionary model system call into question the general importance of DNA methylation in caste differences and evolution in social insects.

Robust DNA Methylation in the Clonal Raider Ant Brain

Social insects are promising model systems for epigenetics due to their immense morphological and behavioral plasticity. Reports that DNA methylation differs between the queen and worker castes in social insects [1-4] have implied a role for DNA methylation in regulating division of labor. To better understand the function of DNA methylation in social insects, we performed whole-genome bisulfite sequencing on brains of the clonal raider ant Cerapachys biroi, whose colonies alternate between reproductive (queen-like) and brood care (worker-like) phases [5]. Many cytosines were methylated in all replicates (on average 29.5% of the methylated cytosines in a given replicate), indicating that a large proportion of the C. biroi brain methylome is robust. Robust DNA methylation occurred preferentially in exonic CpGs of highly and stably expressed genes involved in core functions. Our analyses did not detect any differences in DNA methylation between the queen-like and worker-like phases, suggesting that DNA methylation is not associated with changes in reproduction and behavior in C. biroi. Finally, many cytosines were methylated in one sample only, due to either biological or experimental variation. By applying the statistical methods used in previous studies [1-4, 6] to our data, we show that such sample-specific DNA methylation may underlie the previous findings of queen- and worker-specific methylation. We argue that there is currently no evidence that genome-wide variation in DNA methylation is associated with the queen and worker castes in social insects, and we call for a more careful interpretation of the available data.

Allele specific expression in worker reproduction genes in the bumblebee Bombus terrestris

Methylation has previously been associated with allele specific expression in ants. Recently, we found methylation is important in worker reproduction in the bumblebee Bombus terrestris. Here we searched for allele specific expression in twelve genes associated with worker reproduction in bees. We found allele specific expression in Ecdysone 20 monooxygenase and IMP-L2-like. Although we were unable to confirm a genetic or epigenetic cause for this allele specific expression, the expression patterns of the two genes match those predicted for imprinted genes.

Nourishment level affects caste-related gene expression in Polistes wasps

Background

Social insects exhibit striking phenotypic plasticity in the form of distinct reproductive (queen) and non-reproductive (worker) castes, which are typically driven by differences in the environment during early development. Nutritional environment and nourishment during development has been shown to be broadly associated with caste determination across social insect taxa such as bees, wasps, and termites. In primitively social insects such as Polistes paper wasps, caste remains flexible throughout adulthood, but there is evidence that nourishment inequalities can bias caste development with workers receiving limited nourishment compared to queens. Dominance and vibrational signaling are behaviors that have also been linked to caste differences in paper wasps, suggesting that a combination of nourishment and social factors may drive caste determination. To better understand the molecular basis of nutritional effects on caste determination, we used RNA-sequencing to investigate the gene expression changes in response to proteinaceous nourishment deprivation in Polistes metricus larvae.

Results

We identified 285 nourishment-responsive transcripts, many of which are related to lipid metabolism and oxidation-reduction activity. Via comparisons to previously identified caste-related genes, we found that nourishment restriction only partially biased wasp gene expression patterns toward worker caste-like traits, which supports the notion that nourishment, in conjunction with social environment, is a determinant of developmental caste bias. In addition, we conducted cross-species comparisons of nourishment-responsive genes, and uncovered largely lineage-specific gene expression changes, suggesting few shared nourishment-responsive genes across taxa.

Conclusion

Overall, the results from this study highlight the complex and multifactorial nature of environmental effects on the gene expression patterns underlying plastic phenotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1410-y) contains supplementary material, which is available to authorized users.

Evolution of an epigenetic gene ensemble within the genus Anopheles

Epigenetic control of gene expression has important implications for the regulation of developmental processes, for mediating homeostasis and responses to the external environment, and for transgenerational inheritance of gene expression patterns. Genes that mediate epigenetic control have been well-characterized in Drosophila melanogaster, and we have identified and analyzed an orthologous gene ensemble in Anopheles gambiae that comprises 169 orthologs related to a 215-member epigenetic gene ensemble in D. melanogaster. We find that this ensemble is highly conserved among anopheline mosquitoes, as we identify only seven gene family expansion/contraction events within the ensemble among 12 mosquito species we have studied within the genus Anopheles. Comparative analyses of the epigenetic gene expression across the genera Drosophila and Anopheles reveal distinct tissue-associated expression patterns in the two genera, but similar temporal expression patterns. The A. gambiae complex and D. melanogaster subgroup epigenetic gene ensembles exhibit similar evolutionary rates, as assessed by their respective dN/dS values. These differences in tissue-associated expression patterns, in contrast to similarities in evolutionary rates and temporal expression patterns, may imply that some members of the epigenetic gene ensemble have been redeployed within one or both genera, in comparison to the most recent common ancestor of these two clades. Members of this epigenetic gene ensemble may constitute another set of potential targets for vector control and enable further reductions in the burden of human malaria, by analogy to recent success in development of small molecule antagonists for mammalian epigenetic machinery.

Eusocial insects as emerging models for behavioural epigenetics

Understanding the molecular basis of how behavioural states are established, maintained and altered by environmental cues is an area of considerable and growing interest. Epigenetic processes, including methylation of DNA and post-translational modification of histones, dynamically modulate activity-dependent gene expression in neurons and can therefore have important regulatory roles in shaping behavioural responses to environmental cues. Several eusocial insect species - with their unique displays of behavioural plasticity due to age, morphology and social context - have emerged as models to investigate the genetic and epigenetic underpinnings of animal social behaviour. This Review summarizes recent studies in the epigenetics of social behaviour and offers perspectives on emerging trends and prospects for establishing genetic tools in eusocial insects.

A survey of DNA methylation across social insect species, life stages, and castes reveals abundant and caste-associated methylation in a primitively social wasp

DNA methylation plays an important role in the epigenetic control of developmental and behavioral plasticity, with connections to the generation of striking phenotypic differences between castes (larger, reproductive queens and smaller, non-reproductive workers) in honeybees and ants. Here, we provide the first comparative investigation of caste- and life stage-associated DNA methylation in several species of bees and vespid wasps displaying different levels of social organization. Our results reveal moderate levels of DNA methylation in most bees and wasps, with no clear relationship to the level of sociality. Strikingly, primitively social Polistes dominula paper wasps show unusually high overall DNA methylation and caste-related differences in site-specific methylation. These results suggest DNA methylation may play a role in the regulation of behavioral and physiological differences in primitively social species with more flexible caste differences.

Behavioural syndromes and social insects: personality at multiple levels

Animal personalities or behavioural syndromes are consistent and/or correlated behaviours across two or more situations within a population. Social insect biologists have measured consistent individual variation in behaviour within and across colonies for decades. The goal of this review is to illustrate the ways in which both the study of social insects and of behavioural syndromes has overlapped, and to highlight ways in which both fields can move forward through the synergy of knowledge from each. Here we, (i) review work to date on behavioural syndromes (though not always referred to as such) in social insects, and discuss mechanisms and fitness effects of maintaining individual behavioural variation within and between colonies; (ii) summarise approaches and principles from studies of behavioural syndromes, such as trade-offs, feedback, and statistical methods developed specifically to study behavioural consistencies and correlations, and discuss how they might be applied specifically to the study of social insects; (iii) discuss how the study of social insects can enhance our understanding of behavioural syndromes-research in behavioural syndromes is beginning to explore the role of sociality in maintaining or developing behavioural types, and work on social insects can provide new insights in this area; and (iv) suggest future directions for study, with an emphasis on examining behavioural types at multiple levels of organisation (genes, individuals, colonies, or groups of individuals).

Genomic and epigenomic insights into nutrition and brain disorders

Considerable evidence links many neuropsychiatric, neurodevelopmental and neurodegenerative disorders with multiple complex interactions between genetics and environmental factors such as nutrition. Mental health problems, autism, eating disorders, Alzheimer's disease, schizophrenia, Parkinson's disease and brain tumours are related to individual variability in numerous protein-coding and non-coding regions of the genome. However, genotype does not necessarily determine neurological phenotype because the epigenome modulates gene expression in response to endogenous and exogenous regulators, throughout the life-cycle. Studies using both genome-wide analysis of multiple genes and comprehensive analysis of specific genes are providing new insights into genetic and epigenetic mechanisms underlying nutrition and neuroscience. This review provides a critical evaluation of the following related areas: (1) recent advances in genomic and epigenomic technologies, and their relevance to brain disorders; (2) the emerging role of non-coding RNAs as key regulators of transcription, epigenetic processes and gene silencing; (3) novel approaches to nutrition, epigenetics and neuroscience; (4) gene-environment interactions, especially in the serotonergic system, as a paradigm of the multiple signalling pathways affected in neuropsychiatric and neurological disorders. Current and future advances in these four areas should contribute significantly to the prevention, amelioration and treatment of multiple devastating brain disorders.

Genomewide analysis indicates that queen larvae have lower methylation levels in the honey bee (Apis mellifera)

The honey bee is a social insect characterized by caste differentiation, by which a young larva can develop into either a queen or a worker. Despite possessing the same genome, queen and workers display marked differences in reproductive capacity, physiology, and behavior. Recent studies have shown that DNA methylation plays important roles in caste differentiation. To further explore the roles of DNA methylation in this process, we analyzed DNA methylome profiles of both queen larvae (QL) and worker larvae (WL) of different ages (2, 4, and 6 day old), by using methylated DNA immunoprecipitation-sequencing (meDIP-seq) technique. The global DNA methylation levels varied between the larvae of two castes. DNA methylation increased from 2-day- to 4-day-old QL and then decreased in 6-day-old larvae. In WL, methylation levels increased with age. The methylcytosines in both larvae were enriched in introns, followed by coding sequence (CDS) regions, CpG islands, 2 kbp downstream and upstream of genes, and 5' and 3' untranslated regions (UTRs). The number of differentially methylated genes (DMGs) in 2-, 4-, and 6-day-old QL and WL was 725, 3,013, and 5,049, respectively. Compared to 4- and 6-day-old WL, a large number of genes in QL were downmethylated, which were involved in many processes including development, reproduction, and metabolic regulation. In addition, some DMGs were concerned with caste differentiation.