DNA methylation is maintained with high fidelity in the honey bee germline and exhibits global non-functional fluctuations during somatic development

Multiomics Reveal Associations Between CpG Methylation, Histone Modifications and Transcription in a Species That has Lost DNMT3, the Colorado Potato Beetle

Insects display exceptional phenotypic plasticity, which can be mediated by epigenetic modifications, including CpG methylation and histone modifications. In vertebrates, both are interlinked and CpG methylation is associated with gene repression. However, little is known about these regulatory systems in invertebrates, where CpG methylation is mainly restricted to gene bodies of transcriptionally active genes. A widely conserved mechanism involves the co-transcriptional deposition of H3K36 trimethylation and the targeted methylation of unmethylated CpGs by the de novo DNA methyltransferase DNMT3. However, DNMT3 has been lost multiple times in invertebrate lineages raising the question of how the links between CpG methylation, histone modifications and gene expression are affected by its loss. Here, we report the epigenetic landscape of Leptinotarsa decemlineata, a beetle species that has lost DNMT3 but retained CpG methylation. We combine RNA-seq, enzymatic methyl-seq and CUT&Tag to study gene expression, CpG methylation and patterns of H3K36me3 and H3K27ac histone modifications on a genome-wide scale. Despite the loss of DNMT3, H3K36me3 mirrors CpG methylation patterns. Together, they give rise to signature profiles for expressed and not expressed genes. H3K27ac patterns show a prominent peak at the transcription start site that is predictive of expressed genes irrespective of their methylation status. Our study provides new insights into the evolutionary flexibility of epigenetic modification systems that urge caution when generalizing across species.

Genome-wide DNA methylation patterns in Daphnia magna are not significantly associated with age

BACKGROUND

DNA methylation plays a crucial role in gene regulation and epigenetic inheritance across diverse organisms. Daphnia magna, a model organism in ecological and evolutionary research, has been widely used to study environmental responses, pharmaceutical toxicity, and developmental plasticity. However, its DNA methylation landscape and age-related epigenetic changes remain incompletely understood.

RESULTS

In this study, we characterized DNA methyltransferases (DNMTs) and mapped DNA methylation across the D. magna genome using whole-genome bisulfite sequencing. Our analysis identified three DNMTs: a highly expressed but nonfunctional de novo methyltransferase (DNMT3.1), alongside lowly expressed yet functional de novo methyltransferase (DNMT3.2) and maintenance methyltransferase (DNMT1). D. magna exhibits overall low DNA methylation, targeting primarily CpG dinucleotides. Methylation is sparse at promoters but elevated in the first exons downstream of transcription start sites, with these exons showing hypermethylation relative to adjacent introns. To examine age-associated DNA methylation changes, we analyzed D. magna individuals across multiple life stages. Our results showed no significant global differences in DNA methylation levels between young, mature, and old individuals, nor any age-related clustering in dimensionality reduction analyses. Attempts to construct an epigenetic clock using machine learning models did not yield accurate age predictions, likely due to the overall low DNA methylation levels and lack of robust age-associated methylation changes.

CONCLUSIONS

This study provides a comprehensive characterization of D. magna's DNA methylation landscape and DNMT enzymes, highlighting a distinct pattern of exon-biased CpG methylation. Contrary to prior studies, we found no strong evidence supporting age-associated epigenetic changes, suggesting that DNA methylation may have a limited role in aging in D. magna. These findings enhance our understanding of invertebrate epigenetics and emphasize the need for further research into the interplay between DNA methylation, environmental factors, and gene regulation in D. magna.

Characterization of DNA methylation reader proteins in Arabidopsis thaliana

In plants, cytosine DNA methylation (mC) is largely associated with transcriptional repression of transposable elements, but it can also be found in the body of expressed genes, referred to as gene body methylation (gbM). gbM is correlated with ubiquitously expressed genes; however, its function, or absence thereof, is highly debated. The different outputs that mC can have raise questions as to how it is interpreted-or read-differently in these sequence and genomic contexts. To screen for potential mC-binding proteins, we performed an unbiased DNA affinity pull-down assay combined with quantitative mass spectrometry using methylated DNA probes for each DNA sequence context. All mC readers known to date preferentially bind to the methylated probes, along with a range of new mC-binding protein candidates. Functional characterization of these mC readers, focused on the MBD and SUVH families, was undertaken by ChIP-seq mapping of genome-wide binding sites, their protein interactors, and the impact of high-order mutations on transcriptomic and epigenomic profiles. Together, these results highlight specific context preferences for these proteins, and in particular the ability of MBD2 to bind predominantly to gbM. This comprehensive analysis of Arabidopsis mC readers emphasizes the complexity and interconnectivity between DNA methylation and chromatin remodeling processes in plants.

The functions of DNA methyltransferases during the feeding and development of Haemaphysalis longicornis are potentially associated with lysosome pathways

BACKGROUND

DNA methylation is an epigenetic modification that plays an important role in animal and plant development. Among the diverse types of DNA methylation modifications, methylation of cytosines catalyzed by DNA cytosine methyltransferases (DNMTs) is the most common. Recently, we characterized DNA methyltransferase genes including HlDnmt1 and HlDnmt from the Asian longhorned tick, Haemaphysalis longicornis. However, the dynamic expression and functions of these DNMTs at different developmental stages and feeding statuses of the important vector tick H. longicornis remain unknown.

RESULTS

The expression levels of HlDnmt1 and HlDnmt were significantly different at the four developmental stages: eggs, larvae, nymphs, and adults, with the highest expression levels observed in the larval stage. HlDnmt1 and HlDnmt showed different expression trends in the midguts, ovary, Malpighian tubules, and salivary glands of engorged adults, with the highest expression of HlDnmt1 observed in the ovary and the lowest in the midguts; HlDnmt expression was the highest in the midguts and the lowest in the Malpighian tubules. After RNA interference, the relative expression of HlDnmt1 and HlDnmt in H. longicornis decreased significantly, resulting in a significant decrease in the biting rate of H. longicornis. RNA-seq revealed that the differentially expressed genes were mainly enriched in the biological processes of peptide biosynthesis and the cell components of ribosomes. Molecular functions were mainly concentrated on oxidoreductase activity, ribosome structure composition, serine-type endopeptidase activity, molecular function regulators, and endopeptidase inhibitor activity. KEGG enrichment analysis showed that the differentially expressed genes were mainly enriched in autophagy and lysosome pathways, amino sugar and nucleotide sugar metabolism, glyceride metabolism, ribosomes, and other pathways.

CONCLUSIONS

HlDnmt1 and HlDnmt played an important role during development and feeding of H. longicornis, and their functions were potentially associated with lysosome pathways. These results provide basic knowledge for understanding the epigenetic regulation of the development of the tick H. longicornis, which sheds light on control strategies for ticks and tick-borne diseases.

Non-CG DNA methylation marks the transition from pupa to adult in Helicoverpa armigera

DNA methylation in insects is generally low in abundance, and its role is not well understood. It is often localised in protein coding regions and associated with the expression of 'housekeeping' genes. Few studies have explored DNA methylation dynamics during lifecycle stage transitions in holometabolous (metamorphosing) insects. Using targeted mass spectrometry, we have found a significant difference in global DNA methylation levels between larvae, pupae and adults of Helicoverpa armigera (Lepidoptera: Noctuidae) Hübner, a polyphagous pest of agricultural importance. Whole-genome bisulfite sequencing confirmed these observations and pointed to non-CG context being the primary explanation for the difference observed between pupa and adult. Non-CG methylation was enriched in genes specific to various signalling pathways (Hippo signalling, Hedgehog signalling and mitogen-activated protein kinase (MAPK) signalling) and ATP-dependent chromatin remodelling. Understanding the function of this epigenetic mark could be a target in future studies focusing on integrated pest management.

A role for DNA methylation in bumblebee morphogenesis hints at female-specific developmental erasure

Epigenetic mechanisms, such as DNA methylation, are crucial factors in animal development. In some mammals, almost all DNA methylation is erased during embryo development and re-established in a sex- and cell-specific manner. This erasure and re-establishment is thought to primarily be a vertebrate-specific trait. Insects are particularly interesting in terms of development as many species often undergo remarkable morphological changes en route to maturity, that is, morphogenesis. However, little is known about the role of epigenetic mechanisms in this process across species. We have used whole-genome bisulfite sequencing to track genome-wide DNA methylation changes through the development of an economically and environmentally important pollinator species, the bumblebee Bombus terrestris (Hymenoptera:Apidae Linnaeus). We find overall levels of DNA methylation vary throughout development, and we find developmentally relevant differentially methylated genes throughout. Intriguingly, we have identified a depletion of DNA methylation in ovaries/eggs and an enrichment of highly methylated genes in sperm. We suggest this could represent a sex-specific DNA methylation erasure event. To our knowledge, this is the first suggestion of possible developmental DNA methylation erasure in an insect species. This study lays the required groundwork for functional experimental work to determine if there is a causal nature to the DNA methylation differences identified. Additionally, the application of single-cell methylation sequencing to this system will enable more accurate identification of if or when DNA methylation is erased during development.

Dissecting Mechanisms of Epigenetic Memory Through Computational Modeling

Understanding the mechanistic basis of epigenetic memory has proven to be a difficult task due to the underlying complexity of the systems involved in its establishment and maintenance. Here, we review the role of computational modeling in helping to unlock this complexity, allowing the dissection of intricate feedback dynamics. We focus on three forms of epigenetic memory encoded in gene regulatory networks, DNA methylation, and histone modifications and discuss the important advantages offered by plant systems in their dissection. We summarize the main modeling approaches involved and highlight the principal conceptual advances that the modeling has enabled through iterative cycles of predictive modeling and experiments. Lastly, we discuss remaining gaps in our understanding and how intertwined theory and experimental approaches might help in their resolution.

Unravelling epigenetic mechanisms in Cerastoderma edule genome: a comparison of healthy and neoplastic cockles

Cancer is a multifaceted genetic disease characterized by the acquisition of several essential hallmarks. Notably, certain cancers exhibit horizontal transmissibility, observed across mammalian species and diverse bivalves, the latter referred to as hemic neoplasia. Within this complex landscape, epigenetic mechanisms such as histone modifications and cytosine methylation emerge as fundamental contributors to the pathogenesis of these transmissible cancers. Our study delves into the epigenetic landscape of Cerastoderma edule, focusing on whole-genome methylation and hydroxymethylation profiles in heathy specimens and transmissible neoplasias by means of Nanopore long-read sequencing. Our results unveiled a global hypomethylation in the neoplastic specimens compared to their healthy counterparts, emphasizing the role of DNA methylation in these tumorigenic processes. Furthermore, we verified that intragenic CpG methylation positively correlated with gene expression, emphasizing its role in modulating transcription in healthy and neoplastic cockles, as also highlighted by some up-methylated oncogenic genes. Hydroxymethylation levels were significantly more elevated in the neoplastic samples, particularly within satellites and complex repeats, likely related to structural functions. Additionally, our analysis also revealed distinct methylation and activity patterns in retrotransposons, providing additional insights into bivalve neoplastic processes. Altogether, these findings contribute to understanding the epigenetic dynamics of bivalve neoplasias and shed light on the roles of DNA methylation and hydroxymethylation in tumorigenesis. Understanding these epigenetic alterations holds promise for advancing our broader understanding of cancer epigenetics.

Life-stage specificity and temporal variations in transcriptomes and DNA methylomes of the reef coral Pocillopora damicornis in response to thermal acclimation

Understanding the acclimation capacity of reef corals across generations to thermal stress and its underlying molecular underpinnings could provide insights into their resilience and adaptive responses to future climate change. Here, we acclimated adult brooding coral Pocillopora damicornis to high temperature (32 °C vs. 29 °C) for three weeks and analyzed the changes in phenotypes, transcriptomes and DNA methylomes of adult corals and their brooded larvae. Results showed that although adult corals did not show noticeable bleaching after thermal exposure, they released fewer but larger larvae. Interestingly, larval cohorts from two consecutive lunar days exhibited contrasting physiological resistance to thermal stress, as evidenced by the divergent responses of area-normalized symbiont densities and photochemical efficiency to thermal stress. RNA-seq and whole-genome bisulfite sequencing revealed that adult and larval corals mounted distinct transcriptional and DNA methylation changes in response to thermal stress. Remarkably, larval transcriptomes and DNA methylomes also varied greatly among lunar days and thermal treatments, aligning well with their physiological metrics. Overall, our study shows that changes in transcriptomes and DNA methylomes in response to thermal acclimation can be highly life stage-specific. More importantly, thermally-acclimated adult corals could produce larval offspring with temporally contrasting photochemical performance and thermal resilience, and such variations in larval phenotypes are associated with differential transcriptomes and DNA methylomes, and are likely to increase the likelihood of reproductive success and plasticity of larval propagules under thermal stress.

Modifications to gene body methylation do not alter gene expression plasticity in a reef-building coral

As coral reefs continue to decline due to climate change, the role of coral epigenetics (specifically, gene body methylation, GBM) in coral acclimatization warrants investigation. The evidence is currently conflicting. In diverse animal phyla, the baseline GBM level is associated with gene function: continuously expressed "housekeeping" genes are typically highly methylated, while inducible context-dependent genes have low or no methylation at all. Some authors report an association between GBM and the environment and interpret this observation as evidence of the GBM's role in acclimatization. Yet, others argue that the correlation between GBM change and gene expression change is typically absent or negligible. Here, we used the reef-building coral, Acropora millepora, to test whether environmentally driven changes in GBM are associated with a gene's ability to respond to environmental changes (plasticity) rather than expression level. We analyzed two cases of modified gene expression plasticity observed in a 3-week-long heat acclimatization experiment. The first one was a group of heat-induced genes that failed to revert their expression after the coral was translocated back to the control tank. The second case involved genes that changed the magnitude of their response to the daily temperature fluctuations over the course of the experiment. In both cases, we found negligible or no association with GBM change. We conclude that although both gene expression plasticity and GBM can change during acclimatization, there is no direct association between the two. This adds to the increasing volume of evidence that the function of GBM in invertebrates is unrelated to acclimatization on physiological timescales.

Epigenetics Mechanisms of Honeybees: Secrets of Royal Jelly

Early diets in honeybees have effects on epigenome with consequences on their phenotype. Depending on the early larval diet, either royal jelly (RJ) or royal worker, 2 different female castes are generated from identical genomes, a long-lived queen with fully developed ovaries and a short-lived functionally sterile worker. To generate these prominent physiological and morphological differences between queen and worker, honeybees utilize epigenetic mechanisms which are controlled by nutritional input. These mechanisms include DNA methylation and histone post-translational modifications, mainly histone acetylation. In honeybee larvae, DNA methylation and histone acetylation may be differentially altered by RJ. This diet has biologically active ingredients with inhibitory effects on the de novo methyltransferase DNMT3A or the histone deacetylase 3 HDAC3 to create and maintain the epigenetic state necessary for developing larvae to generate a queen. DNMT and HDAC enzymes work together to induce the formation of a compacted chromatin structure, repressing transcription. Such dialog could be coordinated by their association with other epigenetic factors including the ubiquitin-like containing plant homeodomain (PHD) and really interesting new gene (RING) finger domains 1 (UHRF1). Through its multiple functional domains, UHRF1 acts as an epigenetic reader of both DNA methylation patterns and histone marks. The present review discusses the epigenetic regulation of honeybee's chromatin and how the early diets in honeybees can affect the DNA/histone modifying types of machinery that are necessary to stimulate the larvae to turn into either queen or worker. The review also looks at future directions in epigenetics mechanisms of honeybees, mainly the potential role of UHRF1 in these mechanisms.

Epigenetic regulation and division of labor in social insects

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

Long-term sperm storage in eusocial Hymenoptera

In internally fertilizing species, sperm transfer is not always immediately followed by egg fertilization, and female sperm storage (FSS) may occur. FSS is a phenomenon in which females store sperm in a specialized organ for periods lasting from a few hours to several years, depending on the species. Eusocial hymenopterans (ants, social bees, and social wasps) hold the record for FSS duration. In these species, mating takes place during a single nuptial flight that occurs early in adult life for both sexes; they never mate again. Males die quickly after copulation but survive posthumously as sperm stored in their mates' spermathecae. Reproductive females, also known as queens, have a much longer life expectancy, up to 20 years in some species. Here, we review what is currently known about the molecular adaptations underlying the remarkable FSS capacities in eusocial hymenopterans. Because sperm quality is crucial to the reproductive success of both sexes, we also discuss the mechanisms involved in sperm storage and preservation in the male seminal vesicles prior to ejaculation. Finally, we propose future research directions that should broaden our understanding of this unique biological phenomenon.

Honeybee queen exposure to a widely used fungicide disrupts reproduction and colony dynamic

Pollinators have to cope with a wide range of stressful, not necessarily lethal factors limiting their performance and the ecological services they provide. Among these stressors are pesticides, chemicals that are originally designed to target crop-harming organisms but that also disrupt various functions in pollinators, including flight, communication, orientation and memory. Although all these functions are crucial for reproductive individuals when searching for mates or nesting places, it remains poorly understood how pesticides affect reproduction in pollinators. In this study, we investigated how a widely used fungicide, boscalid, affects reproduction in honey bees (Apis mellifera), an eusocial insect in which a single individual, the queen, fulfills the reproductive functions of the whole colony. Boscalid is a succinate dehydrogenase inhibitor (SDHI) fungicide mainly used on rapeseed flowers to target mitochondrial respiration in fungi but it is also suspected to disrupt foraging-linked functions in bees. We found that immature queen exposure to sublethal, field relevant doses of boscalid disrupted reproduction, as indicated by a dramatic increase in queen mortality during and shortly after the nuptial flights period and a decreased number of spermatozoa stored in the spermatheca of surviving queens. However, we did not observe a decreased paternity frequency in exposed queens that successfully established a colony. Queen exposure to boscalid had detrimental consequences on the colonies they later established regarding brood production, Varroa destructor infection and pollen storage but not nectar storage and population size. These perturbations at the colony-level correspond to nutritional stress conditions, and may have resulted from queen reduced energy provisioning to the eggs. Accordingly, we found that exposed queens had decreased gene expression levels of vitellogenin, a protein involved in egg-yolk formation. Overall, our results indicate that boscalid decreases honey bee queen reproductive quality, thus supporting the need to include reproduction in the traits measured during pesticide risk assessment procedures.

The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L.) response to repeated drought stresses

BACKGROUND

Drought stress has significantly hampered agricultural productivity worldwide and can also result in modifications to DNA methylation levels. However, the dynamics of DNA methylation and its association with the changes in gene transcription and alternative splicing (AS) under drought stress are unknown in linseed, which is frequently cultivated in arid and semiarid regions.

RESULTS

We analysed AS events and DNA methylation patterns in drought-tolerant (Z141) and drought-sensitive (NY-17) linseed under drought stress (DS) and repeated drought stress (RD) treatments. We found that the number of intron-retention (IR) and alternative 3' splice site (Alt3'SS) events were significantly higher in Z141 and NY-17 under drought stress. We found that the linseed response to the DS treatment was mainly regulated by transcription, while the response to the RD treatment was coregulated by transcription and AS. Whole genome-wide DNA methylation analysis revealed that drought stress caused an increase in the overall methylation level of linseed. Although we did not observe any correlation between differentially methylated genes (DMGs) and differentially spliced genes (DSGs) in this study, we found that the DSGs whose gene body region was hypermethylated in Z141 and hypomethylated in NY-17 were enriched in abiotic stress response Gene Ontology (GO) terms. This finding implies that gene body methylation plays an important role in AS regulation in some specific genes.

CONCLUSION

Our study is the first comprehensive genome-wide analysis of the relationship between linseed methylation changes and AS under drought and repeated drought stress. Our study revealed different interaction patterns between differentially expressed genes (DEGs) and DSGs under DS and RD treatments and differences between methylation and AS regulation in drought-tolerant and drought-sensitive linseed varieties. The findings will probably be of interest in the future. Our results provide interesting insights into the association between gene expression, AS, and DNA methylation in linseed under drought stress. Differences in these associations may account for the differences in linseed drought tolerance.

DNA methylation and expression of the egfr gene are associated with worker size in monomorphic ants

The reproductive division of labour is a hallmark of eusocial Hymenoptera. Females are either reproductive queens or non-reproductive workers. In ants, workers often display further task specialisation that is associated with variation in size and/or morphology. Because female polyphenism is typically under environmental control, it is thought epigenetic mechanisms (such as DNA methylation) play a central role since they mediate gene-by-environment interactions. Methylation of the growth-promoting gene epidermal growth factor receptor (egfr) was indeed shown to control worker size in a highly polymorphic ant. However, it remains unknown if egfr methylation could also regulate worker size in monomorphic species. By combining experimental pharmacology and molecular biology, we show that worker size is associated with egfr methylation in two monomorphic ants. Furthermore, we functionally demonstrate that EGFR signalling affects worker size. These results indicate that worker size regulation by egfr methylation has been mechanistically conserved in ants but remains unexploited in monomorphic species.

Comparative studies of genomic and epigenetic factors influencing transcriptional variation in two insect species

Different genes show different levels of expression variability. For example, highly expressed genes tend to exhibit less expression variability. Genes whose promoters have TATA box and initiator motifs tend to have increased expression variability. On the other hand, DNA methylation of transcriptional units, or gene body DNA methylation, is associated with reduced gene expression variability in many species. Interestingly, some insect lineages, most notably Diptera including the canonical model insect Drosophila melanogaster, have lost DNA methylation. Therefore, it is of interest to determine whether genomic features similarly influence gene expression variability in lineages with and without DNA methylation. We analyzed recently generated large-scale data sets in D. melanogaster and honey bee (Apis mellifera) to investigate these questions. Our analysis shows that increased gene expression levels are consistently associated with reduced expression variability in both species, while the presence of TATA box is consistently associated with increased gene expression variability. In contrast, initiator motifs and gene lengths have weak effects limited to some data sets. Importantly, we show that a sequence characteristics indicative of gene body DNA methylation is strongly and negatively associate with gene expression variability in honey bees, while it shows no such association in D. melanogaster. These results suggest the evolutionary loss of DNA methylation in some insect lineages has reshaped the molecular mechanisms concerning the regulation of gene expression variability.

Changes in gene body methylation do not correlate with changes in gene expression in Anthozoa or Hexapoda

Background

As human activity alters the planet, there is a pressing need to understand how organisms adapt to environmental change. Of growing interest in this area is the role of epigenetic modifications, such as DNA methylation, in tailoring gene expression to fit novel conditions. Here, we reanalyzed nine invertebrate (Anthozoa and Hexapoda) datasets to validate a key prediction of this hypothesis: changes in DNA methylation in response to some condition correlate with changes in gene expression.

Results

In accord with previous observations, baseline levels of gene body methylation (GBM) positively correlated with transcription, and negatively correlated with transcriptional variation between conditions. Correlations between changes in GBM and transcription, however, were negligible. There was also no consistent negative correlation between methylation and transcription at the level of gene body methylation class (either highly- or lowly-methylated), anticipated under the previously described “seesaw hypothesis”.

Conclusion

Our results do not support the direct involvement of GBM in regulating dynamic transcriptional responses in invertebrates. If changes in DNA methylation regulate invertebrate transcription, the mechanism must involve additional factors or regulatory influences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08474-z.

Phenotypic Plasticity: What Has DNA Methylation Got to Do with It?

How does one genome give rise to multiple, often markedly different, phenotypes in response to an environmental cue? This phenomenon, known as phenotypic plasticity, is common amongst plants and animals, but arguably the most striking examples are seen in insects. Well-known insect examples include seasonal morphs of butterfly wing patterns, sexual and asexual reproduction in aphids, and queen and worker castes of eusocial insects. Ultimately, we need to understand how phenotypic plasticity works at a mechanistic level; how do environmental signals alter gene expression, and how are changes in gene expression translated into novel morphology, physiology and behaviour? Understanding how plasticity works is of major interest in evolutionary-developmental biology and may have implications for understanding how insects respond to global change. It has been proposed that epigenetic mechanisms, specifically DNA methylation, are the key link between environmental cues and changes in gene expression. Here, we review the available evidence on the function of DNA methylation of insects, the possible role(s) for DNA methylation in phenotypic plasticity and also highlight key outstanding questions in this field as well as new experimental approaches to address these questions.

Differential DNA methylation across environments has no effect on gene expression in the eastern oyster

It has been hypothesized that environmentally induced changes to gene body methylation could facilitate adaptive transgenerational responses to changing environments. We compared patterns of global gene expression (Tag-seq) and gene body methylation (reduced representation bisulfite sequencing) in 80 eastern oysters Crassostrea virginica from six full-sib families, common gardened for 14 months at two sites in the northern Gulf of Mexico that differed in mean salinity. At the time of sampling, oysters from the two sites differed in mass by 60% and in parasite loads by nearly two orders of magnitude. They also differentially expressed 35% of measured transcripts. However, we observed differential methylation at only 1.4% of potentially methylated loci in comparisons between individuals from these different environments, and little correspondence between differential methylation and differential gene expression. Instead, methylation patterns were largely driven by genetic differences among families, with a PERMANOVA analysis indicating nearly a two orders of magnitude greater number of genes differentially methylated between families than between environments. An analysis of CpG observed/expected values (CpG O/E) across the C. virginica genome showed a distinct bimodal distribution, with genes from the first cluster showing the lower CpG O/E values, greater methylation and higher and more stable gene expression, while genes from the second cluster showed lower methylation, and lower and more variable gene expression. Taken together, the differential methylation results suggest that only a small portion of the C. virginica genome is affected by environmentally induced changes in methylation. At this point, there is little evidence to suggest that environmentally induced methylation states would play a leading role in regulating gene expression responses to new environments.

Tools and applications for integrative analysis of DNA methylation in social insects

DNA methylation is a common epigenetic signalling tool and an important biological process which is widely studied in a large array of species. The presence, level and function of DNA methylation vary greatly across species. In some insects, DNA methylation systems are minimal, and overall methylation rates tend to be low in all studied insect species. Low methylation levels probed by whole-genome bisulphite sequencing require great care with respect to data quality control and interpretation. Here, we introduce BWASP/R, a complete workflow that allows efficient, scalable and entirely reproducible analyses of raw DNA methylation sequencing data. Consistent application of quality control filters and analysis parameters provides fair comparisons among different studies and an integrated view of all experiments on one species. We describe the capabilities of the BWASP/R workflow by re-analysing several publicly available social insect WGBS data sets, comprising 70 samples and cumulatively 147 replicates from four different species. We show that the CpG methylome comprises only about 1.5% of CpG sites in the honeybee genome and that the cumulative data are consistent with genetic signatures of site accessibility and physiological control of methylation levels.

Developmental series of gene expression clarifies maternal mRNA provisioning and maternal-to-zygotic transition in a reef-building coral

BACKGROUND

Maternal mRNA provisioning of oocytes regulates early embryogenesis. Maternal transcripts are degraded as zygotic genome activation (ZGA) intensifies, a phenomenon known as the maternal-to-zygotic transition (MZT). Here, we examine gene expression over nine developmental stages in the Pacific rice coral, Montipora capitata, from eggs and embryos at 1, 4, 9, 14, 22, and 36 h-post-fertilization (hpf), as well as swimming larvae (9d), and adult colonies.

RESULTS

Weighted Gene Coexpression Network Analysis revealed four expression peaks, identifying the maternal complement, two waves of the MZT, and adult expression. Gene ontology enrichment revealed maternal mRNAs are dominated by cell division, methylation, biosynthesis, metabolism, and protein/RNA processing and transport functions. The first MZT wave occurs from ~4-14 hpf and is enriched in terms related to biosynthesis, methylation, cell division, and transcription. In contrast, functional enrichment in the second MZT wave, or ZGA, from 22 hpf-9dpf, includes ion/peptide transport and cell signaling. Finally, adult expression is enriched for functions related to signaling, metabolism, and ion/peptide transport. Our proposed MZT timing is further supported by expression of enzymes involved in zygotic transcriptional repression (Kaiso) and activation (Sox2), which peak at 14 hpf and 22 hpf, respectively. Further, DNA methylation writing (DNMT3a) and removing (TET1) enzymes peak and remain stable past ~4 hpf, suggesting that methylome programming occurs before 4 hpf.

CONCLUSIONS

Our high-resolution insight into the coral maternal mRNA and MZT provides essential baseline information to understand parental carryover effects and the sensitivity of developmental success under increasing environmental stress.

Ecotoxicological epigenetics in invertebrates: Emerging tool for the evaluation of present and past pollution burden

The effect of environmental pollution on epigenetic changes and their heredity in affected organisms is of major concern as such changes can play a significant role in adaptation to changing environmental conditions. Changes of epigenetic marks including DNA methylation, histone modifications, and non-coding RNA's can induce changes in gene transcription leading to physiological long-term changes or even transgenerational inheritance. Such mechanisms have until recently been scarcely studied in invertebrate organisms, mainly focusing on model species including Caenorhabditis elegans and Daphnia magna. However, more data are becoming available, particularly focused on DNA methylation changes caused by anthropogenic pollutants in a wide range of invertebrates. This review examines the literature from field and laboratory studies utilising invertebrate species exposed to environmental pollutants and their effect on DNA methylation. Possible mechanisms of epigenetic modifications and their role on physiology and adaptation as well as the incidence of intergenerational and transgenerational inheritance are discussed. Furthermore, critical research challenges are defined and the way forward is proposed. Future studies should focus on the use of next generation sequencing tools to define invertebrate methylomes under environmental stress in higher resolution, those data should further be linked to gene expression patterns and phenotypes and detailed studies focusing on transgenerational effects are encouraged. Moreover, studies of other epigenetic mechanisms in various invertebrate species, apart from DNA methylation would provide better understanding of interconnected cross-talk between epigenetic marks. Taken together incorporating epigenetic studies in ecotoxicology context presents a promising tool for development of sensitive biomarkers for environmental stress assessment.

The pattern of DNA methylation alteration, and its association with the changes of gene expression and alternative splicing during phosphate starvation in tomato

DNA methylation is changed and associates with gene expression alterations in plant response to phosphate starvation (Pi-), a common stress that impacts plant growth and productivity. However, in the horticultural model species Solanum lycopersicum (tomato), the dynamics of DNA methylation and its relationship with changes in gene transcription and alternative splicing (AS) under Pi- are unknown. Here, we performed integrative methylome and transcriptome analyses of tomato seedlings under Pi-deficient and -sufficient conditions. We found Pi- caused a slight increase in the overall methylation level, with millions of differentially methylated cytosines (DmCs) and a few hundred differentially methylated regions (DMRs). We also identified thousands of differentially expressed (DE) and differential AS (DAS) genes induced by Pi-, and found that DmCs were more abundant in non-expressed genes than in DE or DAS genes. Moreover, DNA methylation alterations weakly correlated with transcription changes but not with DAS events, and hyper-CHH-DMRs overlapping with transposable elements (TEs) were enriched in a subset of Pi starvation response (PSR) genes. We propose that changes in DNA methylation may be associated with the differential expression of some PSR genes, but that most of these changes probably control the expression of nearby TEs, rather than directly affecting the transcription or AS of PSR genes. Besides, the pattern of methylation changes upon Pi- may largely be shaped by TE distributions. Together, our study provides comprehensive insights into the association of DNA methylation with gene transcription and AS under Pi- in tomato and may contribute to unveiling novel roles of epigenetic mechanisms in plant stress response.

Origins of human disease: the chrono-epigenetic perspective

Epigenetics has enriched human disease studies by adding new interpretations to disease features that cannot be explained by genetic and environmental factors. However, identifying causal mechanisms of epigenetic origin has been challenging. New opportunities have risen from recent findings in intra-individual and cyclical epigenetic variation, which includes circadian epigenetic oscillations. Cytosine modifications display deterministic temporal rhythms, which may drive ageing and complex disease. Temporality in the epigenome, or the 'chrono' dimension, may help the integration of epigenetic, environmental and genetic disease studies, and reconcile several disparities stemming from the arbitrarily delimited research fields. The ultimate goal of chrono-epigenetics is to predict disease risk, age of onset and disease dynamics from within individual-specific temporal dynamics of epigenomes.

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

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

The role of epigenetics, particularly DNA methylation, in the evolution of caste in insect societies

Eusocial insects can be defined as those that live in colonies and have distinct queens and workers. For most species, queens and workers arise from a common genome, and so caste-specific developmental trajectories must arise from epigenetic processes. In this review, we examine the epigenetic mechanisms that may be involved in the regulation of caste dimorphism. Early work on honeybees suggested that DNA methylation plays a causal role in the divergent development of queen and worker castes. This view has now been challenged by studies that did not find consistent associations between methylation and caste in honeybees and other species. Evidence for the involvement of methylation in modulating behaviour of adult workers is also inconsistent. Thus, the functional significance of DNA methylation in social insects remains equivocal. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

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

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

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

The emergence of the brain non-CpG methylation system in vertebrates

Mammalian brains feature exceptionally high levels of non-CpG DNA methylation alongside the canonical form of CpG methylation. Non-CpG methylation plays a critical regulatory role in cognitive function, which is mediated by the binding of MeCP2, the transcriptional regulator that when mutated causes Rett syndrome. However, it is unclear whether the non-CpG neural methylation system is restricted to mammalian species with complex cognitive abilities or has deeper evolutionary origins. To test this, we investigated brain DNA methylation across 12 distantly related animal lineages, revealing that non-CpG methylation is restricted to vertebrates. We discovered that in vertebrates, non-CpG methylation is enriched within a highly conserved set of developmental genes transcriptionally repressed in adult brains, indicating that it demarcates a deeply conserved regulatory program. We also found that the writer of non-CpG methylation, DNMT3A, and the reader, MeCP2, originated at the onset of vertebrates as a result of the ancestral vertebrate whole-genome duplication. Together, we demonstrate how this novel layer of epigenetic information assembled at the root of vertebrates and gained new regulatory roles independent of the ancestral form of the canonical CpG methylation. This suggests that the emergence of non-CpG methylation may have fostered the evolution of sophisticated cognitive abilities found in the vertebrate lineage.

Impacts of Epigenetic Processes on the Health and Productivity of Livestock

The dynamic changes in the epigenome resulting from the intricate interactions of genetic and environmental factors play crucial roles in individual growth and development. Numerous studies in plants, rodents, and humans have provided evidence of the regulatory roles of epigenetic processes in health and disease. There is increasing pressure to increase livestock production in light of increasing food needs of an expanding human population and environment challenges, but there is limited related epigenetic data on livestock to complement genomic information and support advances in improvement breeding and health management. This review examines the recent discoveries on epigenetic processes due to DNA methylation, histone modification, and chromatin remodeling and their impacts on health and production traits in farm animals, including bovine, swine, sheep, goat, and poultry species. Most of the reports focused on epigenome profiling at the genome-wide or specific genic regions in response to developmental processes, environmental stressors, nutrition, and disease pathogens. The bulk of available data mainly characterized the epigenetic markers in tissues/organs or in relation to traits and detection of epigenetic regulatory mechanisms underlying livestock phenotype diversity. However, available data is inadequate to support gainful exploitation of epigenetic processes for improved animal health and productivity management. Increased research effort, which is vital to elucidate how epigenetic mechanisms affect the health and productivity of livestock, is currently limited due to several factors including lack of adequate analytical tools. In this review, we (1) summarize available evidence of the impacts of epigenetic processes on livestock production and health traits, (2) discuss the application of epigenetics data in livestock production, and (3) present gaps in livestock epigenetics research. Knowledge of the epigenetic factors influencing livestock health and productivity is vital for the management and improvement of livestock productivity.

The effect of DNA methylation on bumblebee colony development

BACKGROUND

Although around 1% of cytosines in bees' genomes are known to be methylated, less is known about methylation's effect on bee behavior and fitness. Chemically altered DNA methylation levels have shown clear changes in the dominance and reproductive behavior of workers in queen-less colonies, but the global effect of DNA methylation on caste determination and colony development remains unclear, mainly because of difficulties in controlling for genetic differences among experimental subjects in the parental line. Here, we investigated the effect of the methylation altering agent decitabine on the developmental rate of full bumblebee colonies. Whole genome bisulfite sequencing was used to assess differences in methylation status.

RESULTS

Our results showed fewer methylated loci in the control group. A total of 22 CpG loci were identified as significantly differentially methylated between treated and control workers with a change in methylation levels of 10% or more. Loci that were methylated differentially between groups participated in pathways including neuron function, oocyte regulation and metabolic processes. Treated colonies tended to develop faster, and therefore more workers were found at a given developmental stage. However, male production followed the opposite trend and it tended to be higher in control colonies.

CONCLUSION

Overall, our results indicate that altered methylation patterns resulted in an improved cooperation between workers, while there were no signs of abnormal worker dominance or caste determination.

Intergenerational transfer of DNA methylation marks in the honey bee

The evolutionary significance of epigenetic inheritance is controversial. While epigenetic marks such as DNA methylation can affect gene function and change in response to environmental conditions, their role as carriers of heritable information is often considered anecdotal. Indeed, near-complete DNA methylation reprogramming, as occurs during mammalian embryogenesis, is a major hindrance for the transmission of nongenetic information between generations. Yet it remains unclear how general DNA methylation reprogramming is across the tree of life. Here we investigate the existence of epigenetic inheritance in the honey bee. We studied whether fathers can transfer epigenetic information to their daughters through DNA methylation. We performed instrumental inseminations of queens, each with four different males, retaining half of each male's semen for whole genome bisulfite sequencing. We then compared the methylation profile of each father's somatic tissue and semen with the methylation profile of his daughters. We found that DNA methylation patterns were highly conserved between tissues and generations. There was a much greater similarity of methylomes within patrilines (i.e., father-daughter subfamilies) than between patrilines in each colony. Indeed, the samples' methylomes consistently clustered by patriline within colony. Samples from the same patriline had twice as many shared methylated sites and four times fewer differentially methylated regions compared to samples from different patrilines. Our findings indicate that there is no DNA methylation reprogramming in bees and, consequently, that DNA methylation marks are stably transferred between generations. This points to a greater evolutionary potential of the epigenome in invertebrates than there is in mammals.

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

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

Benchmarking DNA methylation assays in a reef-building coral

Interrogation of chromatin modifications, such as DNA methylation, has the potential to improve forecasting and conservation of marine ecosystems. The standard method for assaying DNA methylation (whole genome bisulphite sequencing), however, is currently too costly to apply at the scales required for ecological research. Here, we evaluate different methods for measuring DNA methylation for ecological epigenetics. We compare whole genome bisulphite sequencing (WGBS) with methylated CpG binding domain sequencing (MBD-seq), and a modified version of MethylRAD we term methylation-dependent restriction site-associated DNA sequencing (mdRAD). We evaluate these three assays in measuring variation in methylation across the genome, between genotypes, and between polyp types in the reef-building coral Acropora millepora. We find that all three assays measure absolute methylation levels similarly for gene bodies (gbM), as well as exons and 1 Kb windows with a minimum Pearson correlation 0.66. Differential gbM estimates were less correlated, but still concurrent across assays. We conclude that MBD-seq and mdRAD are reliable and cost-effective alternatives to WGBS. The considerably lower sequencing effort required for mdRAD to produce comparable methylation estimates makes it particularly useful for ecological epigenetics.

Contiguous and stochastic CHH methylation patterns of plant DRM2 and CMT2 revealed by single-read methylome analysis

Cytosine methylome data is commonly generated through next-generation sequencing, with analyses averaging methylation states of individual reads. We propose an alternative method of analysing single-read methylome data. Using this method, we identify patterns relating to the mechanism of two plant non-CG-methylating enzymes, CMT2 and DRM2. CMT2-methylated regions show higher stochasticity, while DRM2-methylated regions have higher variation among cells. Based on these patterns, we develop a classifier that predicts enzyme activity in different species and tissues. To facilitate further single-read analyses, we develop a genome browser, SRBrowse, optimised for visualising and analysing sequencing data at single-read resolution.

DNA methylation patterns in the tobacco budworm, Chloridea virescens

DNA methylation is an important epigenetic modification that is prone to stochastic variation and is responsive to environmental factors. Yet changes in DNA methylation could persist across generations and thus play an important role in evolution. In this study, we used methylation-sensitive amplified fragment length polymorphisms (MS-AFLP) to evaluate whether DNA methylation could contribute to the evolution of the sexual communication signal in the noctuid moth Chloridea virescens. We found that most DNA methylation was consistent across tissues, although some methylation sites were specifically found in pheromone glands. We also found significant DNA methylation differences among families and two pheromone phenotype selection lines, and these differences correlated with genetic variation. Most DNA methylation patterns were inherited, although some sites were subject to spontaneous de novo DNA methylation across generations. Thus, DNA methylation likely plays a role in a wide range of processes in moths. Together, our results present an important initial step towards understanding the potential role of DNA methylation in the evolution of sexual communication signals in moths.

Functional genomics of parental care of insects

Parental care was likely the first step most lineages made towards sociality. However, the molecular mechanisms that generate parental care are not broadly characterized. Insects are important as an evolutionary independent group from classic models of parental care, such as, house mice. They provide an opportunity to test the generality of our understanding. With this review, I survey the functional genomics of parental care of insects, summarize several recent advances in the broader framework for studying and understanding parental care, and finish with suggested priorities for further research. Although there are too few studies to draw definitive conclusions, I argue that natural selection appears to be rewiring existing gene networks to produce parental care, that the epigenetic mechanisms influencing parental care are not well understood, and, as an interesting early consensus, that genes strongly associated with carer/offspring interactions appear biased towards proteins that are secreted. I summarize the studies that have functionally validate candidate genes and highlight the increasing need to perform this work. I finish with arguments for both conceptual and practical changes moving forward. I argue that future work can increase the use of predictive frameworks, broaden its definition of conservation of mechanism to gene networks rather than single genes, and increase the use of more established comparative methods. I further highlight the practical considerations of standardizing analyses and reporting, increasing the sampling of both carers and offspring, better characterizing gene regulatory networks, better characterizing taxonomically restricted genes and any consistent role they have underpinning parental care, and using factorial designs to disentangle the influence of multiple variables on the expression of parental care.

Evolution of DNA Methylome Diversity in Eukaryotes

Cytosine DNA methylation (5mC) is a widespread base modification in eukaryotic genomes with critical roles in transcriptional regulation. In recent years, our understanding of 5mC has changed because of advances in 5mC detection techniques that allow mapping of this mark on the whole genome scale. Profiling DNA methylomes from organisms across the eukaryotic tree of life has reshaped our views on the evolution of 5mC. In this review, we explore the macroevolution of 5mC in major eukaryotic groups, and then focus on recent advances made in animals. Genomic 5mC patterns as well as the mechanisms of 5mC deposition tend to be evolutionary labile across large phylogenetic distances; however, some common patterns are starting to emerge. Within the animal kingdom, 5mC diversity has proven to be much greater than anticipated. For example, a previously held common view that genome hypermethylation is a trait exclusive to vertebrates has recently been challenged. Also, data from genome-wide studies are starting to yield insights into the potential roles of 5mC in invertebrate cis regulation. Here we provide an evolutionary perspective of both the well-known and enigmatic roles of 5mC across the eukaryotic tree of life.