MicroRNAs in Honey Bee Caste Determination

Uncovering the Changing Gene Expression Profile of Honeybee (Apis mellifera) Worker Larvae Transplanted to Queen Cells

The reproductive division of labor, based on caste differentiation in social insects, is of great significance in evolution. Generally, a healthy bee colony consists of a queen and numerous workers and drones. Despite being genetically identical, the queen and workers exhibit striking differences in morphology, behavior, and lifespan. The fertilized eggs and larvae selectively develop into queen and worker bees depending on the local nutrition and environment. Bee worker larvae that are transplanted within 3 days of age to queen cells of a bee colony can develop into queens with mature ovaries. This phenomenon is important to understand the regulatory mechanisms of caste differentiation. In this study, we transplanted worker larvae (Apis mellifera) at the age of 1 (L1), 2 (L2), and 3 days (L3) into queen cells until the age of 4 days. Subsequently, genetic changes in these larvae were evaluated. The results revealed that the number of differentially expressed genes (DEGs) in L1 vs. L3 was more than that in L1 vs. L2. Furthermore, many of the genes that were downregulated are mostly involved in metabolism, body development, reproductive ability, and longevity, indicating that these functions decreased with the age of transplantation of the larvae. Moreover, these functions may be critical for worker larvae to undergo the developmental path to become queens. We also found that the DEGs of L1 vs. L2 and L1 vs. L3 were enriched in the MAPK, FoxO, mTOR, Wnt, TGF-beta Hedgehog Toll and Imd, and Hippo signaling pathways. Gene ontology analysis indicated that some genes are simultaneously involved in different biological pathways; through these genes, the pathways formed a mutual regulatory network. Casein kinase 1 (CK 1) was predicted to participate in the FoxO, Wnt, Hedgehog, and Hippo signaling pathways. The results suggest that these pathways cross talked through the network to modify the development of larvae and that CK 1 is an important liaison. The results provide valuable information regarding the regulatory mechanism of environmental factors affecting queen development, thus, amplifying the understanding of caste differentiation in bees.

Phenotypically distinct female castes in honey bees are defined by alternative chromatin states during larval development

The capacity of the honey bee to produce three phenotypically distinct organisms (two female castes; queens and sterile workers, and haploid male drones) from one genotype represents one of the most remarkable examples of developmental plasticity in any phylum. The queen-worker morphological and reproductive divide is environmentally controlled during post-embryonic development by differential feeding. Previous studies implicated metabolic flux acting via epigenetic regulation, in particular DNA methylation and microRNAs, in establishing distinct patterns of gene expression underlying caste-specific developmental trajectories. We produce the first genome-wide maps of chromatin structure in the honey bee at a key larval stage in which developmental canalization into queen or worker is virtually irreversible. We find extensive genome-wide differences in H3K4me3, H3K27ac, and H3K36me3, many of which correlate with caste-specific transcription. Furthermore, we identify H3K27ac as a key chromatin modification, with caste-specific regions of intronic H3K27ac directing the worker caste. These regions may harbor the first examples of caste-specific enhancer elements in the honey bee. Our results demonstrate a key role for chromatin modifications in the establishment and maintenance of caste-specific transcriptional programs in the honey bee. We show that at 96 h of larval growth, the queen-specific chromatin pattern is already established, whereas the worker determination is not, thus providing experimental support for the perceived timing of this critical point in developmental heterochrony in two types of honey bee females. In a broader context, our study provides novel data on environmentally regulated organismal plasticity and the molecular foundation of the evolutionary origins of eusociality.

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.

Plant MicroRNAs in Cross-Kingdom Regulation of Gene Expression

MicroRNAs (miRNAs) are a class of noncoding small RNAs, which play a crucial role in post-transcriptional gene regulation. Recently, various reports revealed that miRNAs could be transmitted between species to mediate cross-kingdom regulation by integrating into a specific target gene-mediated regulatory pathway to exert relevant biological functions. Some scholars and researchers have observed this as an attractive hypothesis that may provide a foundation for novel approaches in the diagnosis, prognosis, and treatment of disease. Meanwhile, others deem the mentioned results were obtained from a “false positive effect” of performed experiments. Here, we focus on several current studies concerning plant miRNA-mediated cross-kingdom regulation (from both fronts) and discuss the existing issues that need further consideration. We also discuss possible miRNA horizontal transfer mechanisms from one species to another and analyze the relationship between miRNA-mediated cross-kingdom regulation and coevolution during a long-term specific host–pathogen interaction.

The evolution of posttranscriptional regulation

"DNA makes RNA makes protein." After transcription, mRNAs undergo a series of intertwining processes to be finally translated into functional proteins. The "posttranscriptional" regulation (PTR) provides cells an extended option to fine-tune their proteomes. To meet the demands of complex organism development and the appropriate response to environmental stimuli, every step in these processes needs to be finely regulated. Moreover, changes in these regulatory processes are important driving forces underlying the evolution of phenotypic differences across different species. The major PTR mechanisms discussed in this review include the regulation of splicing, polyadenylation, decay, and translation. For alternative splicing and polyadenylation, we mainly discuss their evolutionary dynamics and the genetic changes underlying the regulatory differences in cis-elements versus trans-factors. For mRNA decay and translation, which, together with transcription, determine the cellular RNA or protein abundance, we focus our discussion on how their divergence coordinates with transcriptional changes to shape the evolution of gene expression. Then to highlight the importance of PTR in the evolution of higher complexity, we focus on their roles in two major phenomena during eukaryotic evolution: the evolution of multicellularity and the division of labor between different cell types and tissues; and the emergence of diverse, often highly specialized individual phenotypes, especially those concerning behavior in eusocial insects. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution Translation > Translation Regulation RNA Processing > Splicing Regulation/Alternative Splicing.

A comparative analysis of corpora allata-corpora cardiaca microRNA repertoires revealed significant changes during mosquito metamorphosis

The corpora allata (CA) are a pair of endocrine glands with neural connections to the brain and close association with another neuroendocrine organ, the corpora cardiaca (CC). The CA from adult female Aedes aegypti mosquitoes synthesize fluctuating levels of juvenile hormone (JH), which have been linked to the ovarian development and are influenced by nutritional signals. In this study, we investigated the potential involvement of microRNAs (miRNAs), a type of small non-coding RNAs, in the regulation of gene expression in CA-CC complexes during mosquito reproductive development, at stages with distinct JH biosynthesis patterns. We analyzed the miRNA repertoires expressed in the CA-CC of pupae, sugar-fed and blood-fed female Ae. aegypti. In total, 156 mature miRNAs were detected in the CA-CC, with 84 displaying significant differences in expression among the three CA-CC developmental stages. There were more miRNAs that were expressed in pupae, and decreased or were absent after adult emergence, when compared with changes between CA-CC of sugar and blood-fed females. Analysis of the genes identified as potential targets for the CA-CC miRNA repertoires classified them into the broad categories of metabolism, information storage and processing, and cellular processes and signaling; with genes involved in cellular processes and signaling representing the largest portion. Among them, the signal-transduction mechanisms and intracellular trafficking, secretion and vesicular transport contained almost 55% of the genes' targets. A substantial number of miRNAs were differentially abundant in the libraries of the three developmental stages, and those changes were much more notable when pupae and adult stages were compared. We detected putative binding sites for some of the most abundant miRNAs on genes encoding JH biosynthetic enzymes and CC neuropeptides. These studies should help us to gain a better understanding of the regulation of CA-CC activity mediated by miRNAs during major developmental stages in mosquitoes.

A molecular concept of caste in insect societies

The term 'caste' is used to describe the division of reproductive labour that defines eusocial insect societies. The definition of 'caste' has been debated over the last 50 years, specifically with respect to the simplest insect societies; this raises the question of whether a simple categorisation of social behaviour by reproductive state alone is helpful. Gene-level analyses of behaviours of individuals in hymenopteran social insect societies now provide a new empirical base-line for defining caste and understanding the evolution and maintenance of a reproductive division of labour. We review this literature to identify a set of potential molecular signatures that, combined with behavioural, morphological and physiological data, help define caste more precisely; these signatures vary with the type of society, and are likely to be influenced by ecology, life-history, and stage in the colony cycle. We conclude that genomic approaches provide us with additional ways to help quantify and categorise caste, and behaviour in general.

Beyond Royalactin and a master inducer explanation of phenotypic plasticity in honey bees

Distinct female castes produced from one genotype are the trademark of a successful evolutionary invention in eusocial insects known as reproductive division of labour. In honey bees, fertile queens develop from larvae fed a complex diet called royal jelly. Recently, one protein in royal jelly, dubbed Royalactin, was deemed to be the exclusive driver of queen bee determination. However, this notion has not been universally accepted. Here I critically evaluate this line of research and argue that the sheer complexity of creating alternate phenotypes from one genotype cannot be reduced to a single dietary component. An acceptable model of environmentally driven caste differentiation should include the facets of dynamic thinking, such as the concepts of attractor states and genetic hierarchical networks.

In honeybees, genotypically identical females develop into queens or sterile workers, depending on their diets. In this review, Ryszard Maleszka discusses the controversial role of the royal jelly protein Royalactin in caste determination and provides a framework for moving beyond the master inducer concept.

Contrasting Sex-and Caste-Dependent piRNA Profiles in the Transposon Depleted Haplodiploid Honeybee Apis mellifera

Protecting genome integrity against transposable elements is achieved by intricate molecular mechanisms involving PIWI proteins, their associated small RNAs (piRNAs), and epigenetic modifiers such as DNA methylation. Eusocial bees, in particular the Western honeybee, Apis mellifera, have one of the lowest contents of transposable elements in the animal kingdom, and, unlike other animals with a functional DNA methylation system, appear not to methylate their transposons. This raises the question of whether the PIWI machinery has been retained in this species. Using comparative genomics, mass spectrometry, and expressional profiling, we present seminal evidence that the piRNA system is conserved in honeybees. We show that honey bee piRNAs contain a 2'-O-methyl modification at the 3' end, and have a bias towards a 5' terminal U, which are signature features of their biogenesis. Both piRNA repertoire and expression levels are greater in reproductive individuals than in sterile workers. Haploid males, where the detrimental effects of transposons are dominant, have the greatest piRNA levels, but surprisingly, the highest expression of transposons. These results show that even in a transposon-depleted species, the piRNA system is required to guard the vulnerable haploid genome and reproductive castes against transposon-associated genomic instability. This also suggests that dosage plays an important role in the regulation of transposons and piRNAs expression in haplo-diploid systems.

A possible role of DNA methylation in functional divergence of a fast evolving duplicate gene encoding odorant binding protein 11 in the honeybee

Although gene duplication is seen as the main path to evolution of new functions, molecular mechanisms by which selection favours the gain versus loss of newly duplicated genes and minimizes the fixation of pseudo-genes are not well understood. Here, we investigate in detail a duplicate honeybee gene obp11 belonging to a fast evolving insect gene family encoding odorant binding proteins (OBPs). We report that obp11 is expressed only in female bees in rare antennal sensilla basiconica in contrast to its tandem partner obp10 that is expressed in the brain in both females and males (drones). Unlike all other obp genes in the honeybee, obp11 is methylated suggesting that functional diversification of obp11 and obp10 may have been driven by an epigenetic mechanism. We also show that increased methylation in drones near one donor splice site that correlates with higher abundance of a transcript variant encoding a truncated OBP11 protein is one way of controlling its contrasting expression. Our data suggest that like in mammals and plants, DNA methylation in insects may contribute to functional diversification of proteins produced from duplicated genes, in particular to their subfunctionalization by generating complementary patterns of expression.

Transcriptome sequencing reveals high isoform diversity in the ant Formica exsecta

Transcriptome resources for social insects have the potential to provide new insight into polyphenism, i.e., how divergent phenotypes arise from the same genome. Here we present a transcriptome based on paired-end RNA sequencing data for the ant Formica exsecta (Formicidae, Hymenoptera). The RNA sequencing libraries were constructed from samples of several life stages of both sexes and female castes of queens and workers, in order to maximize representation of expressed genes. We first compare the performance of common assembly and scaffolding software (Trinity, Velvet-Oases, and SOAPdenovo-trans), in producing de novo assemblies. Second, we annotate the resulting expressed contigs to the currently published genomes of ants, and other insects, including the honeybee, to filter genes that have annotation evidence of being true genes. Our pipeline resulted in a final assembly of altogether 39,262 mRNA transcripts, with an average coverage of >300X, belonging to 17,496 unique genes with annotation in the related ant species. From these genes, 536 genes were unique to one caste or sex only, highlighting the importance of comprehensive sampling. Our final assembly also showed expression of several splice variants in 6,975 genes, and we show that accounting for splice variants affects the outcome of downstream analyses such as gene ontologies. Our transcriptome provides an outstanding resource for future genetic studies on F. exsecta and other ant species, and the presented transcriptome assembly can be adapted to any non-model species that has genomic resources available from a related taxon.

Integration of lncRNA-miRNA-mRNA reveals novel insights into oviposition regulation in honey bees

BACKGROUND

The honey bee (Apis mellifera) is a highly diverse species commonly used for honey production and pollination services. The oviposition of the honey bee queen affects the development and overall performance of the colony. To investigate the ovary activation and oviposition processes on a molecular level, a genome-wide analysis of lncRNAs, miRNAs and mRNA expression in the ovaries of the queens was performed to screen for differentially expressed coding and noncoding RNAs. Further analysis identified relevant candidate genes or RNAs.

RESULTS

The analysis of the RNA profiles in different oviposition phase of the queens revealed that 740 lncRNAs, 81 miRNAs and 5,481 mRNAs were differently expressed during the ovary activation; 88 lncRNAs, 13 miRNAs and 338 mRNAs were differently expressed during the oviposition inhibition process; and finally, 100 lncRNAs, four miRNAs and 497 mRNAs were differently expressed during the oviposition recovery process. In addition, functional annotation of differentially expressed RNAs revealed several pathways that are closely related to oviposition, including hippo, MAPK, notch, Wnt, mTOR, TGF-beta and FoxO signaling pathways. Furthermore, in the QTL region for ovary size, 73 differentially expressed genes and 14 differentially expressed lncRNAs were located, which are considered as candidate genes affecting ovary size and oviposition. Moreover, a core set of genes served as bridges among different miRNAs were identified through the integrated analysis of lncRNA-miRNA-mRNA network.

CONCLUSION

The observed dramatic expression changes of coding and noncoding RNAs suggest that they may play a critical role in honey bee queens' oviposition. The identified candidate genes for oviposition activation and regulation could serve as a resource for further studies of genetic markers of oviposition in honey bees.

Gene duplication and the evolution of phenotypic diversity in insect societies

Gene duplication is an important evolutionary process thought to facilitate the evolution of phenotypic diversity. We investigated if gene duplication was associated with the evolution of phenotypic differences in a highly social insect, the honeybee Apis mellifera. We hypothesized that the genetic redundancy provided by gene duplication could promote the evolution of social and sexual phenotypes associated with advanced societies. We found a positive correlation between sociality and rate of gene duplications across the Apoidea, indicating that gene duplication may be associated with sociality. We also discovered that genes showing biased expression between A. mellifera alternative phenotypes tended to be found more frequently than expected among duplicated genes than singletons. Moreover, duplicated genes had higher levels of caste-, sex-, behavior-, and tissue-biased expression compared to singletons, as expected if gene duplication facilitated phenotypic differentiation. We also found that duplicated genes were maintained in the A. mellifera genome through the processes of conservation, neofunctionalization, and specialization, but not subfunctionalization. Overall, we conclude that gene duplication may have facilitated the evolution of social and sexual phenotypes, as well as tissue differentiation. Thus this study further supports the idea that gene duplication allows species to evolve an increased range of phenotypic diversity.

Detecting rare asymmetrically methylated cytosines and decoding methylation patterns in the honeybee genome

Context-dependent gene expression in eukaryotes is controlled by several mechanisms including cytosine methylation that primarily occurs in the CG dinucleotides (CpGs). However, less frequent non-CpG asymmetric methylation has been found in various cell types, such as mammalian neurons, and recent results suggest that these sites can repress transcription independently of CpG contexts. In addition, an emerging view is that CpG hemimethylation may arise not only from deregulation of cellular processes but also be a standard feature of the methylome. Here, we have applied a novel approach to examine whether asymmetric CpG methylation is present in a sparsely methylated genome of the honeybee, a social insect with a high level of epigenetically driven phenotypic plasticity. By combining strand-specific ultra-deep amplicon sequencing of illustrator genes with whole-genome methylomics and bioinformatics, we show that rare asymmetrically methylated CpGs can be unambiguously detected in the honeybee genome. Additionally, we confirm differential methylation between two phenotypically and reproductively distinct castes, queens and workers, and offer new insight into the heterogeneity of brain methylation patterns. In particular, we challenge the assumption that symmetrical methylation levels reflect symmetry in the underlying methylation patterns and conclude that hemimethylation may occur more frequently than indicated by methylation levels. Finally, we question the validity of a prior study in which most of cytosine methylation in this species was reported to be asymmetric.

Detecting rare asymmetrically methylated cytosines and decoding methylation patterns in the honeybee genome

Context-dependent gene expression in eukaryotes is controlled by several mechanisms including cytosine methylation that primarily occurs in the CG dinucleotides (CpGs). However, less frequent non-CpG asymmetric methylation has been found in various cell types, such as mammalian neurons, and recent results suggest that these sites can repress transcription independently of CpG contexts. In addition, an emerging view is that CpG hemimethylation may arise not only from deregulation of cellular processes but also be a standard feature of the methylome. Here, we have applied a novel approach to examine whether asymmetric CpG methylation is present in a sparsely methylated genome of the honeybee, a social insect with a high level of epigenetically driven phenotypic plasticity. By combining strand-specific ultra-deep amplicon sequencing of illustrator genes with whole-genome methylomics and bioinformatics, we show that rare asymmetrically methylated CpGs can be unambiguously detected in the honeybee genome. Additionally, we confirm differential methylation between two phenotypically and reproductively distinct castes, queens and workers, and offer new insight into the heterogeneity of brain methylation patterns. In particular, we challenge the assumption that symmetrical methylation levels reflect symmetry in the underlying methylation patterns and conclude that hemimethylation may occur more frequently than indicated by methylation levels. Finally, we question the validity of a prior study in which most of cytosine methylation in this species was reported to be asymmetric.

Plant microRNAs in larval food regulate honeybee caste development

The major environmental determinants of honeybee caste development come from larval nutrients: royal jelly stimulates the differentiation of larvae into queens, whereas beebread leads to worker bee fate. However, these determinants are not fully characterized. Here we report that plant RNAs, particularly miRNAs, which are more enriched in beebread than in royal jelly, delay development and decrease body and ovary size in honeybees, thereby preventing larval differentiation into queens and inducing development into worker bees. Mechanistic studies reveal that amTOR, a stimulatory gene in caste differentiation, is the direct target of miR162a. Interestingly, the same effect also exists in non-social Drosophila. When such plant RNAs and miRNAs are fed to Drosophila larvae, they cause extended developmental times and reductions in body weight and length, ovary size and fecundity. This study identifies an uncharacterized function of plant miRNAs that fine-tunes honeybee caste development, offering hints for understanding cross-kingdom interaction and co-evolution.

How caste has formed in honeybees is an enduring puzzle. The prevailing view is that royal jelly stimulates the differentiation of larvae into queen. Here, we uncover a new mechanism that plant miRNAs in worker bee’s food postpone larval development, thereby inducing sterile worker bees. Thus, the theories about honeybee caste formation need to be re-examined from a new angle besides the traditional focus on royal jelly and its components. Furthermore, since miRNAs are transmitted between species of different kingdoms and can contribute to the phenotype regulation, this new model of horizontal miRNA transfer may open up a new avenue to further study the molecular mechanisms underlying cross-kingdom interaction and co-evolution.

MicroRNAs Associated with Caste Determination and Differentiation in a Primitively Eusocial Insect

In eusocial Hymenoptera (ants, bees and wasps), queen and worker adult castes typically arise via environmental influences. A fundamental challenge is to understand how a single genome can thereby produce alternative phenotypes. A powerful approach is to compare the molecular basis of caste determination and differentiation along the evolutionary trajectory between primitively and advanced eusocial species, which have, respectively, relatively undifferentiated and strongly differentiated adult castes. In the advanced eusocial honeybee, Apis mellifera, studies suggest that microRNAs (miRNAs) play an important role in the molecular basis of caste determination and differentiation. To investigate how miRNAs affect caste in eusocial evolution, we used deep sequencing and Northern blots to isolate caste-associated miRNAs in the primitively eusocial bumblebee Bombus terrestris. We found that the miRNAs Bte-miR-6001-5p and -3p are more highly expressed in queen- than in worker-destined late-instar larvae. These are the first caste-associated miRNAs from outside advanced eusocial Hymenoptera, so providing evidence for caste-associated miRNAs occurring relatively early in eusocial evolution. Moreover, we found little evidence that miRNAs previously shown to be associated with caste in A. mellifera were differentially expressed across caste pathways in B. terrestris, suggesting that, in eusocial evolution, the caste-associated role of individual miRNAs is not conserved.

Age-dependent transcriptional and epigenomic responses to light exposure in the honey bee brain

Light is a powerful environmental stimulus of special importance in social honey bees that undergo a behavioral transition from in-hive to outdoor foraging duties. Our previous work has shown that light exposure induces structural neuronal plasticity in the mushroom bodies (MBs), a brain center implicated in processing inputs from sensory modalities. Here, we extended these analyses to the molecular level to unravel light-induced transcriptomic and epigenomic changes in the honey bee brain. We have compared gene expression in brain compartments of 1- and 7-day-old light-exposed honey bees with age-matched dark-kept individuals. We have found a number of differentially expressed genes (DEGs), both novel and conserved, including several genes with reported roles in neuronal plasticity. Most of the DEGs show age-related changes in the amplitude of light-induced expression and are likely to be both developmentally and environmentally regulated. Some of the DEGs are either known to be methylated or are implicated in epigenetic processes suggesting that responses to light exposure are at least partly regulated at the epigenome level. Consistent with this idea light alters the DNA methylation pattern of bgm, one of the DEGs affected by light exposure, and the expression of microRNA miR-932. This confirms the usefulness of our approach to identify candidate genes for neuronal plasticity and provides evidence for the role of epigenetic processes in driving the molecular responses to visual stimulation.

Mitochondrial structure and dynamics as critical factors in honey bee (Apis mellifera L.) caste development

The relationship between nutrition and phenotype is an especially challenging question in cases of facultative polyphenism, like the castes of social insects. In the honey bee, Apis mellifera, unexpected modifications in conserved signaling pathways revealed the hypoxia response as a possible mechanism underlying the regulation of body size and organ growth. Hence, the current study was designed to investigate possible causes of why the three hypoxia core genes are overexpressed in worker larvae. Parting from the hypothesis that this has an endogenous cause and is not due to differences in external oxygen levels we investigated mitochondrial numbers and distribution, as well as mitochondrial oxygen consumption rates in fat body cells of queen and worker larvae during the caste fate-critical larval stages. By immunofluorescence and electron microscopy we found higher densities of mitochondria in queen larval fat body, a finding further confirmed by a citrate synthase assay quantifying mitochondrial functional units. Oxygen consumption measurements by high-resolution respirometry revealed that queen larvae have higher maximum capacities of ATP production at lower physiological demand. Finally, the expression analysis of mitogenesis-related factors showed that the honey bee TFB1 and TFB2 homologs, and a nutritional regulator, ERR, are overexpressed in queen larvae. These results are strong evidence that the differential nutrition of queen and worker larvae by nurse bees affects mitochondrial dynamics and functionality in the fat body of these larvae, hence explaining their differential hypoxia response.

Epigenetic code and insect behavioural plasticity

Although the nature of the genetic control of adaptive behaviours in insects is a major unresolved problem it is now understood that epigenetic mechanisms, bound by genetic constraints, are prime drivers of brain plasticity arising from both developmental and experience-dependent events. With the recent advancements in methylomics and emerging analyses of histones and non-protein-coding RNAs, insect epigenetics is well positioned to ask more direct questions and importantly, address them experimentally. To achieve rapid progress, insect epigenetics needs to focus on mechanistic explanations of epigenomic dynamics and move beyond low-depth genome-wide analyses to cell-type specific epigenomics. One topic of a high priority is the impact of sequence variants on generating differential methylation patterns and their contribution to behavioural plasticity.