The honey bee epigenomes: differential methylation of brain DNA in queens and workers

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.

Reminiscences on the honeybee genome project and the rise of epigenetic concepts in insect science

The sequencing of the honeybee genome in 2006 was an important technological and logistic achievement experience. But what benefits have flown from the honeybee genome project? What does the annotated genomic assembly mean for the study of behavioural complexity and organismal function in honeybees? Here, I discuss several lines of research that have arisen from this project and highlight the rapidly expanding studies on insect epigenomics, emergent properties of royal jelly, the mechanism of nutritional control of development and the contribution of epigenomic regulation to the evolution of sociality. I also argue that the term 'insect epigenetics' needs to be carefully redefined to reflect the diversity of epigenomic toolkits in insects and the impact of lineage-specific innovations on organismal outcomes. The honeybee genome project helped pioneer advances in social insect molecular biology, and fuelled breakthrough research into the role of flexible epigenomic control systems in linking genotype to phenotype.

A review of environmental epigenetics in aquatic invertebrates

Aquatic ecosystems face significant challenges due to increasing human-induced environmental stressors. Recent studies emphasize the role of epigenetic mechanisms in the stress responses and adaptations of organisms to those stressors. Epigenetics influences gene expression, enabling phenotypic plasticity and transgenerational effects. Therefore, understanding the epigenetic responses of aquatic invertebrates to environmental stressors is imperative for aquatic ecosystem research. In this study, we organize the mechanisms of epigenetics in aquatic invertebrates and explore their roles in the responses of aquatic invertebrates to environmental stressors. Furthermore, we discuss the inheritance of epigenetic changes and their influence across generations in aquatic invertebrates. A comprehensive understanding of epigenetic responses is crucial for long-term ecosystem management and conservation strategies in the face of irreversible climate change in aquatic environments. In this review, we synthesize existing knowledge about environmental epigenetics in aquatic invertebrates to provide insights and suggest directions for future research.

Epigenetics and seasonal timing in animals: a concise review

Seasonal adaptation in animals is a complex process that involves genetic, epigenetic, and environmental factors. The present review explores recent studies on epigenetic mechanisms implicated in seasonal adaptation in animals. The review is divided into three main sections, each focusing on a different epigenetic mechanism: DNA methylation, histone modifications, and non-coding RNA. Additionally, the review delves into the current understanding of how these epigenetic factors contribute to the regulation of circadian and seasonal cycles. Understanding these molecular mechanisms provides the first step in deciphering the complex interplay between genetics, epigenetics, and the environment in driving seasonal adaptation in animals. By exploring these mechanisms, a better understanding of how animals adapt to changing environmental conditions can be achieved.

Genome-wide profiling of DNA methylome and transcriptome reveals epigenetic regulation of Urechis unicinctus response to sulfide stress

Sulfide is a well-known environmental pollutant that can have detrimental effects on most organisms. However, few metazoans living in sulfide-rich environments have developed mechanisms to tolerate and adapt to sulfide stress. Epigenetic mechanisms, including DNA methylation, have been shown to play a vital role in environmental stress adaptation. Nevertheless, the precise function of DNA methylation in biological sulfide adaptation remains unclear. Urechis unicinctus, a benthic organism inhabiting sulfide-rich intertidal environments, is an ideal model organism for studying adaptation to sulfide environments. In this study, we conducted a comprehensive analysis of the DNA methylome and transcriptome of U. unicinctus after exposure to 50 μM sulfide. The results revealed dynamic changes in the DNA methylation (5-methylcytosine) landscape in response to sulfide stress, with U. unicinctus exhibiting elevated DNA methylation levels following stress exposure. Integrating differentially expressed genes (DEGs) and differentially methylated regions (DMRs), we identified a crucial role of gene body methylation in predicting gene expression. Furthermore, using a DNA methyltransferase inhibitor, we validated the involvement of DNA methylation in the sulfide stress response and the gene regulatory network influenced by DNA methylation. The results indicated that by modulating DNA methylation levels during sulfide stress, the expression of glutathione S-transferase, glutamyl aminopeptidase, and cytochrome c oxidase could be up-regulated, thereby facilitating the metabolism and detoxification of exogenous sulfides. Moreover, DNA methylation was found to regulate and enhance the oxidative phosphorylation pathway, including NADH dehydrogenase, isocitrate dehydrogenase, and ATP synthase. Additionally, DNA methylation influenced the regulation of Cytochrome P450 and macrophage migration inhibitory factor, both of which are closely associated with oxidative stress and stress resistance. Our findings not only emphasize the role of DNA methylation in sulfide adaptation but also provide novel insights into the potential mechanisms through which marine organisms adapt to environmental changes.

Pervasive relaxed selection in termite genomes

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

Transcriptome analysis of Apis mellifera antennae reveals molecular divergence underlying the division of labour in worker bees

The olfactory system plays a fundamental role in mediating insect behaviour. Worker bees exhibit an age-dependent division of labour, performing discrete sets of behaviours throughout their lifespan. The behavioural states of bees rely on their sense of the environment and chemical communication via their olfactory system, the antennae. However, the olfactory adaptation mechanism of worker bees during their behavioural development remains unclear. In this study, we conducted a comprehensive and quantitative analysis of antennal gene expression in the Apis mellifera of newly emerged workers, nurses, foragers and defenders using RNA-seq. We found that the antenna tissues of honey bees continued developing after transformation from newly emerged workers to adults. Additionally, we identified differentially expressed genes associated with bee development and division of labour. We validated that major royal jelly protein genes are highly and specifically expressed in nurse honey bee workers. Furthermore, we identified and validated significant alternative splicing events correlated with the development and division of labour. These findings provide a comprehensive transcriptome profile and a new perspective on the molecular mechanisms that may underlie the worker honey bee division of labour.

Mixed Patterns of Intergenerational DNA Methylation Inheritance in Acropora

For sessile organisms at high risk from climate change, phenotypic plasticity can be critical to rapid acclimation. Epigenetic markers like DNA methylation are hypothesized as mediators of plasticity; methylation is associated with the regulation of gene expression, can change in response to ecological cues, and is a proposed basis for the inheritance of acquired traits. Within reef-building corals, gene-body methylation (gbM) can change in response to ecological stressors. If coral DNA methylation is transmissible across generations, this could potentially facilitate rapid acclimation to environmental change. We investigated methylation heritability in Acropora, a stony reef-building coral. Two Acropora millepora and two Acropora selago adults were crossed, producing eight offspring crosses (four hybrid, two of each species). We used whole-genome bisulfite sequencing to identify methylated loci and allele-specific alignments to quantify per-locus inheritance. If methylation is heritable, differential methylation (DM) between the parents should equal DM between paired offspring alleles at a given locus. We found a mixture of heritable and nonheritable loci, with heritable portions ranging from 44% to 90% among crosses. gBM was more heritable than intergenic methylation, and most loci had a consistent degree of heritability between crosses (i.e. the deviation between parental and offspring DM were of similar magnitude and direction). Our results provide evidence that coral methylation can be inherited but that heritability is heterogenous throughout the genome. Future investigations into this heterogeneity and its phenotypic implications will be important to understanding the potential capability of intergenerational environmental acclimation in reef building corals.

Metabolic and immune costs balance during natural acclimation of corals in fluctuating environments

Epigenetic modifications based on DNA methylation can rapidly improve the potential of corals to adapt to environmental pressures by increasing their phenotypic plasticity, a factor important for scleractinian corals to adapt to future global warming. However, the extent to which corals develop similar adaptive mechanisms and their specific adaptation processes remain unclear. Here, to reveal the regulatory mechanism by which DNA methylation improves thermal tolerance in Pocillopora damicornis under fluctuating environments, we analyzed genome-wide DNA methylation signatures in P. damicornis and compared the differences in the methylation and transcriptional responses of P. damicornis from fluctuating and stable environments using whole-genome bisulfite sequencing and nanopore-based RNA sequencingtranscriptome sequencing. We discovered low methylation levels in P. damicornis (average methylation 4.14%), with CpG accounting for 74.88%, CHH for 13.27%, and CHG for 11.85% of this methylation. However, methylation levels did not change between coral samples from the fluctuating and stable environments. The varied methylation levels in different regions of the gene revealed that the overall methylation level of the gene body was relatively high and showed a bimodal methylation pattern. Methylation occurs primarily in exons rather than introns within the gene body In P. damicornis, there was only a weak correlation between methylation and transcriptional changes at the individual gene level, and the methylation and gene expression levels generally exhibited a bell-shaped relationship, which we speculate may be due to the specificity of cnidarian species. Correlation analysis between methylation levels and the transcriptome revealed that the highest proportion of the top 20 enriched KEGG pathways was related to immunity. Additionally, P. damicornis collected from a high-temperature pool had a lower metabolic rate than those collected from a low-temperature pool. We hypothesize that the dynamic balance of energy-expenditure costs between immunity and metabolism is an important strategy for increasing P. damicornis tolerance. The fluctuating environment of high-temperature pools may increase the heat tolerance in corals by increasing their immunity and thus lowering their metabolism.

(Epi)genomic adaptation driven by fine geographical scale environmental heterogeneity after recent biological invasions

Elucidating processes and mechanisms involved in rapid local adaptation to varied environments is a poorly understood but crucial component in management of invasive species. Recent studies have proposed that genetic and epigenetic variation could both contribute to ecological adaptation, yet it remains unclear on the interplay between these two components underpinning rapid adaptation in wild animal populations. To assess their respective contributions to local adaptation, we explored epigenomic and genomic responses to environmental heterogeneity in eight recently colonized ascidian (Ciona intestinalis) populations at a relatively fine geographical scale. Based on MethylRADseq data, we detected strong patterns of local environment-driven DNA methylation divergence among populations, significant epigenetic isolation by environment (IBE), and a large number of local environment-associated epigenetic loci. Meanwhile, multiple genetic analyses based on single nucleotide polymorphisms (SNPs) showed genomic footprints of divergent selection. In addition, for five genetically similar populations, we detected significant methylation divergence and local environment-driven methylation patterns, indicating the strong effects of local environments on epigenetic variation. From a functional perspective, a majority of functional genes, Gene Ontology (GO) terms, and biological pathways were largely specific to one of these two types of variation, suggesting partial independence between epigenetic and genetic adaptation. The methylation quantitative trait loci (mQTL) analysis showed that the genetic variation explained only 18.67% of methylation variation, further confirming the autonomous relationship between these two types of variation. Altogether, we highlight the complementary interplay of genetic and epigenetic variation involved in local adaptation, which may jointly promote populations' rapid adaptive capacity and successful invasions in different environments. The findings here provide valuable insights into interactions between invaders and local environments to allow invasive species to rapidly spread, thus contributing to better prediction of invasion success and development of management strategies.

Whole-Genome Bisulfite Sequencing Protocol for the Analysis of Genome-Wide DNA Methylation and Hydroxymethylation Patterns at Single-Nucleotide Resolution

The analysis of genome-wide epigenomic alterations including DNA methylation and hydroxymethylation has become a subject of intensive research for many biological and clinical questions. DNA methylation analysis bears the particular promise to supplement or replace biochemical and imaging-based tests for the next generation of personalized medicine. Whole-genome bisulfite sequencing (WGBS) using next-generation sequencing technologies is currently considered the gold standard for a comprehensive and quantitative analysis of DNA methylation throughout the genome. However, bisulfite conversion does not allow distinguishing between cytosine methylation and hydroxymethylation requiring an additional chemical or enzymatic step to identify hydroxymethylated cytosines. Here, we provide a detailed protocol based on a commercial kit for the preparation of sequencing libraries for the comprehensive whole-genome analysis of DNA methylation and/or hydroxymethylation. The protocol is based on the construction of sequencing libraries from limited amounts of input DNA by ligation of methylated adaptors to the fragmented DNA prior to bisulfite conversion. For analyses requiring a quantitative distinction between 5-methylcytosine and 5-hydroxymethylcytosines levels, an oxidation step is included in the same workflow to perform oxidative bisulfite sequencing (OxBs-Seq). In this case, two sequencing libraries will be generated and sequenced: a classic methylome following bisulfite conversion and analyzing modified cytosines (not distinguishing between methylated and hydroxymethylated cytosines) and a methylome analyzing only methylated cytosines, respectively. Hydroxymethylation levels are deduced from the differences between the two reactions. We also provide a step-by-step description of the data analysis using publicly available bioinformatic tools. The described protocol has been successfully applied to different human and plant samples and yields robust and reproducible results.

DNA Methylation in the Hypothalamic Feeding Center and Obesity

Obesity rates have been increasing worldwide for decades, mainly due to environmental factors, such as diet, nutrition, and exercise. However, the molecular mechanisms through which environmental factors induce obesity remain unclear. Several mechanisms underlie the body's response to environmental factors, and one of the main mechanisms involves epigenetic modifications, such as DNA methylation. The pattern of DNA methylation is influenced by environmental factors, and altered DNA methylation patterns can affect gene expression profiles and phenotypes. DNA methylation may mediate the development of obesity caused by environmental factors. Similar to the factors governing obesity, DNA methylation is influenced by nutrients and metabolites. Notably, DNA methylation is associated with body size and weight programming. The DNA methylation levels of proopiomelanocortin (Pomc) and neuropeptide Y (Npy) in the hypothalamic feeding center, a key region controlling systemic energy balance, are affected by diet. Conditional knockout mouse studies of epigenetic enzymes have shown that DNA methylation in the hypothalamic feeding center plays an indispensable role in energy homeostasis. In this review, we discuss the role of DNA methylation in the hypothalamic feeding center as a potential mechanism underlying the development of obesity induced by environmental factors.

DNA methylation extends lifespan in the bumblebee Bombus terrestris

Epigenetic alterations are a primary hallmark of ageing. In mammals, age-related epigenetic changes alter gene expression profiles, disrupt cellular homeostasis and physiological functions and, therefore, promote ageing. It remains unclear whether ageing is also driven by epigenetic mechanisms in invertebrates. Here, we used a pharmacological hypomethylating agent (RG108) to evaluate the effects of DNA methylation (DNAme) on lifespan in an insect-the bumblebee Bombus terrestris. RG108 extended mean lifespan by 43% and induced the differential methylation of genes involved in hallmarks of ageing, including DNA damage repair and chromatin organization. Furthermore, the longevity gene sirt1 was overexpressed following the treatment. Functional experiments demonstrated that SIRT1 protein activity was positively associated with lifespan. Overall, our study indicates that epigenetic mechanisms are conserved regulators of lifespan in both vertebrates and invertebrates and provides new insights into how DNAme is involved in the ageing process in insects.

Methyl-CpG binding domain (MBD)2/3 specifically recognizes and binds to the genomic mCpG site with a β-sheet in the MBD to affect embryonic development in Bombyx mori

Methyl-CpG (mCpG) binding domain (MBD) proteins especially bind with methylated DNA, and are involved in many important biological processes; however, the binding mechanism between insect MBD2/3 and mCpG remains unclear. In this study, we identified 2 isoforms of the MBD2/3 gene in Bombyx mori, MBD2/3-S and MBD2/3-L. Binding analysis of MBD2/3-L, MBD2/3-S, and 7 mutant MBD2/3-L proteins deficient in β1-β6 or α1 in the MBD showed that β2-β3-turns in the β-sheet of the MBD are necessary for the formation of the MBD2/3-mCpG complex; furthermore, other secondary structures, namely, β4-β6 and an α-helix, play a role in stabilizing the β-sheet structure to ensure that the MBD is able to bind mCpG. In addition, sequence alignment and binding analyses of different insect MBD2/3s indicated that insect MBD2/3s have an intact and conserved MBD that binds to the mCpG of target genes. Furthermore, MBD2/3 RNA interference results showed that MBD2/3-L plays a role in regulating B. mori embryonic development, similar to that of DNA methylation; however, MBD2/3-S without β4-β6 and α-helix does not alter embryonic development. These results suggest that MBD2/3-L recognizes and binds to mCpG through the intact β-sheet structure in its MBD, thus ensuring silkworm embryonic development.

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.

Mechanistic Insights into Alternative Gene Splicing in Oxidative Stress and Tissue Injury

Significance: Oxidative stress (OS) and inflammation are inducers of tissue injury. Alternative splicing (AS) is an essential regulatory step for diversifying the eukaryotic proteome. Human diseases link AS to OS; however, the underlying mechanisms must be better understood. Recent Advances: Genome‑wide profiling studies identify new differentially expressed genes induced by OS-dependent ischemia/reperfusion injury. Overexpression of RNA-binding protein RBFOX1 protects against inflammation. Hypoxia-inducible factor-1α directs polypyrimidine tract binding protein 1 to regulate mouse carcinoembryonic antigen-related cell adhesion molecule 1 (Ceacam1) AS under OS conditions. Heterogeneous nuclear ribonucleoprotein L variant 1 contains an RGG/RG motif that coordinates with transcription factors to influence human CEACAM1 AS. Hypoxia intervention involving short interfering RNAs directed to long-noncoding RNA 260 polarizes M2 macrophages toward an anti-inflammatory phenotype and alleviates OS by inhibiting IL-28RA gene AS. Critical Issues: Protective mechanisms that eliminate reactive oxygen species (ROS) are important for resolving imbalances that lead to chronic inflammation. Defects in AS can cause ROS generation, cell death regulation, and the activation of innate and adaptive immune factors. We propose that AS pathways link redox regulation to the activation or suppression of the inflammatory response during cellular stress. Future Directions: Emergent studies using molecule-mediated RNA splicing are being conducted to exploit the immunogenicity of AS protein products. Deciphering the mechanisms that connect misspliced OS and pathologies should remain a priority. Controlled release of RNA directly into cells with clinical applications is needed as the demand for innovative nucleic acid delivery systems continues to be demonstrated.

Transgenerational Epigenetic DNA Methylation Editing and Human Disease

During gestation, maternal (F0), embryonic (F1), and migrating primordial germ cell (F2) genomes can be simultaneously exposed to environmental influences. Accumulating evidence suggests that operating epi- or above the genetic DNA sequence, covalent DNA methylation (DNAme) can be recorded onto DNA in response to environmental insults, some sites which escape normal germline erasure. These appear to intrinsically regulate future disease propensity, even transgenerationally. Thus, an organism's genome can undergo epigenetic adjustment based on environmental influences experienced by prior generations. During the earliest stages of mammalian development, the three-dimensional presentation of the genome is dramatically changed, and DNAme is removed genome wide. Why, then, do some pathological DNAme patterns appear to be heritable? Are these correctable? In the following sections, I review concepts of transgenerational epigenetics and recent work towards programming transgenerational DNAme. A framework for editing heritable DNAme and challenges are discussed, and ethics in human research is introduced.

The PWWP domain and the evolution of unique DNA methylation toolkits in Hymenoptera

DNMT3 in Hymenoptera has a unique duplication of the essential PWWP domain. Using GST-tagged PWWP fusion proteins and histone arrays we show that these domains have gained new properties and represent the first case of PWWP domains binding to H3K27 chromatin modifications, including H3K27me3, a key modification that is important during development. Phylogenetic analyses of 107 genomes indicate that the duplicated PWWP domains separated into two sister clades, and their distinct binding capacities are supported by 3D modeling. Other features of this unique DNA methylation system include variable copies, losses, and duplications of DNMT1 and DNMT3, and combinatorial generations of DNMT3 isoforms including variants missing the catalytic domain. Some of these losses and duplications of are found only in parasitic wasps. We discuss our findings in the context of the crosstalk between DNA methylation and histone methylation, and the expanded potential of epigenomic modifications in Hymenoptera to drive evolutionary novelties.

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

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

Third-Generation Sequencing Reveals the Adaptive Role of the Epigenome in Three Deep-Sea Polychaetes

The roles of DNA methylation in invertebrates are poorly characterized, and critical data are missing for the phylum Annelida. We fill this knowledge gap by conducting the first genome-wide survey of DNA methylation in the deep-sea polychaetes dominant in deep-sea vents and seeps: Paraescarpia echinospica, Ridgeia piscesae, and Paralvinella palmiformis. DNA methylation calls were inferred from Oxford Nanopore sequencing after assembling high-quality genomes of these animals. The genomes of these worms encode all the key enzymes of the DNA methylation metabolism and possess a mosaic methylome similar to that of other invertebrates. Transcriptomic data of these polychaetes support the hypotheses that gene body methylation strengthens the expression of housekeeping genes and that promoter methylation acts as a silencing mechanism but not the hypothesis that DNA methylation suppresses the activity of transposable elements. The conserved epigenetic profiles of genes responsible for maintaining homeostasis under extreme hydrostatic pressure suggest DNA methylation plays an important adaptive role in these worms.

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.

Integrated Analysis of DNA Methylome and Transcriptome Reveals Epigenetic Regulation of Cold Tolerance in Litopenaeus vannamei

DNA methylation is an important epigenetic modification that has been shown to be associated with responses to non-biological stressors. However, there is currently no research on DNA methylation in response to environmental signals in shrimp. In this study, we conducted a comprehensive comparative analysis of DNA methylation profiles and differentially expressed genes between two strains of Litopenaeus vannamei with significantly different cold tolerance through whole genome bisulfite sequencing (WGBS) and transcriptome sequencing. Between Lv-C and Lv-T (constant temperature of 28 °C and low temperatures of 18 °C and 10 °C) under cytosine-guanine (CG) environments, 39,100 differentially methylated regions (DMRs) were identified, corresponding to 9302 DMR-related genes (DMRGs). The DMRs were mainly located in the gene body (exons and introns). Gene Ontology (GO) analysis showed that these DMRGs were significantly enriched in cell parts, catalytic activity, and metabolic processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed significant enrichment of these DMRGs in pathways such as proteasome (ko03050), oxidative phosphorylation (ko00190), mTOR signaling pathway (ko04150), fatty acid metabolism (ko01212), and fatty acid degradation (ko00071). The comprehensive results suggested that L. vannamei mainly regulates gene expression in response to low temperatures through hypermethylation or demethylation of some genes involved in thermogenesis, glycolysis, the autophagy pathway, the peroxisome, and drug metabolism pathways. These results provide important clues for studying DNA methylation patterns and identifying cold tolerance genes in shrimp.

Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector

BACKGROUND

The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes.

RESULTS

A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups.

CONCLUSIONS

The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.

Potential role of N6-adenine DNA methylation in alternative splicing and endosymbiosis in Paramecium bursaria

N⁶-adenine DNA methylation (6mA), a rediscovered epigenetic mark in eukaryotic organisms, diversifies in abundance, distribution, and function across species, necessitating its study in more taxa. Paramecium bursaria is a typical model organism with endosymbiotic algae of the species Chlorella variabilis. This consortium therefore serves as a valuable system to investigate the functional role of 6mA in endosymbiosis, as well as the evolutionary importance of 6mA among eukaryotes. In this study, we report the first genome-wide, base pair-resolution map of 6mA in P. bursaria and identify its methyltransferase PbAMT1. Functionally, 6mA exhibits a bimodal distribution at the 5' end of RNA polymerase II-transcribed genes and possibly participates in transcription by facilitating alternative splicing. Evolutionarily, 6mA co-evolves with gene age and likely serves as a reverse mark of endosymbiosis-related genes. Our results offer new insights for the functional diversification of 6mA in eukaryotes as an important epigenetic mark.

The Epigenetics of Animal Personality

Animal personality, consistent individual differences in behaviour, is an important concept for understanding how individuals vary in how they cope with environmental challenges. In order to understand the evolutionary significance of animal personality, it is crucial to understand the underlying regulatory mechanisms. Epigenetic marks such as DNA methylation are hypothesised to play a major role in explaining variation in phenotypic changes in response to environmental alterations. Several characteristics of DNA methylation also align well with the concept of animal personality. In this review paper, we summarise the current literature on the role that molecular epigenetic mechanisms may have in explaining personality variation. We elaborate on the potential for epigenetic mechanisms to explain behavioural variation, behavioural development and temporal consistency in behaviour. We then suggest future routes for this emerging field and point to potential pitfalls that may be encountered. We conclude that a more inclusive approach is needed for studying the epigenetics of animal personality and that epigenetic mechanisms cannot be studied without considering the genetic background.

Phenotypic dimorphism between honeybee queen and worker is regulated by complicated epigenetic modifications

Phenotypic dimorphism between queens and workers is an important biological characteristic of honeybees that has been the subject of intensive research. The enormous differences in morphology, lifespan, physiology, and behavior between queens and workers are caused by a complicated set of factors. Epigenetic modifications are considered to play an important role in this process. In this study, we analyzed the differences in chromosome interactions and H3K27ac and H3K4me1 modifications between the queens and workers using high-throughput chromosome conformation capture (Hi-C) and Chromatin immunoprecipitation followed by sequencing (ChIP-Seq) technologies. We found that the queens contain more chromosome interactions and more unique H3K27ac modifications than workers; in contrast, workers have more H3K4me1 modifications than queens. Moreover, we identified Map3k15 as a potential caste gene in queen-worker differentiation. Our results suggest that chromosomal conformation and H3K27ac and H3K4me1 modifications are involved in regulating queen-worker differentiation, which reveals that the queen-worker phenotypic dimorphism is regulated by multiple epigenetic modifications.

DNMT1 mutant ants develop normally but have disrupted oogenesis

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

DNA methylation is associated with codon degeneracy in a species of bumblebee

Social insects display extreme phenotypic differences between sexes and castes even though the underlying genome can be almost identical. Epigenetic processes have been proposed as a possible mechanism for mediating these phenotypic differences. Using whole genome bisulfite sequencing of queens, males, and reproductive female workers we have characterised the sex- and caste-specific methylome of the bumblebee Bombus terrestris. We have identified a potential role for DNA methylation in histone modification processes which may influence sex and caste phenotypic differences. We also find differentially methylated genes generally show low levels of DNA methylation which may suggest a separate function for lowly methylated genes in mediating transcriptional plasticity, unlike highly methylated genes which are usually involved in housekeeping functions. We also examined the relationship between the underlying genome and the methylome using whole genome re-sequencing of the same queens and males. We find DNA methylation is enriched at zero-fold degenerate sites. We suggest DNA methylation may be acting as a targeted mutagen at these sites, providing substrate for selection via non-synonymous changes in the underlying genome. However, we did not see any relationship between DNA methylation and rates of positive selection in our samples. In order to fully assess a possible role for DNA methylation in adaptive processes a specifically designed study using natural population data is needed.

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.

Single nucleotide polymorphism (SNP) in the doublesex ( dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex ( dsx) gene, involved in mosquitos' somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene ( Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

Epigenetic and Genetic Population Structure is Coupled in a Marine Invertebrate

Delineating the relative influence of genotype and the environment on DNA methylation is critical for characterizing the spectrum of organism fitness as driven by adaptation and phenotypic plasticity. In this study, we integrated genomic and DNA methylation data for two distinct Olympia oyster (Ostrea lurida) populations while controlling for within-generation environmental influences. In addition to providing the first characterization of genome-wide DNA methylation patterns in the oyster genus Ostrea, we identified 3,963 differentially methylated loci between populations. Our results show a clear coupling between genetic and epigenetic patterns of variation, with 27% of variation in interindividual methylation differences explained by genotype. Underlying this association are both direct genetic changes in CpGs (CpG-SNPs) and genetic variation with indirect influence on methylation (mQTLs). When comparing measures of genetic and epigenetic population divergence at specific genomic regions this relationship surprisingly breaks down, which has implications for the methods commonly used to study epigenetic and genetic coupling in marine invertebrates.

Autophagy Is Required to Sustain Increased Intestinal Cell Proliferation during Phenotypic Plasticity Changes in Honey Bee (Apis mellifera)

Tissue phenotypic plasticity facilitates rapid adaptation of organisms to biotic and/or abiotic pressure. The reproductive capacity of honey bee workers (Apis mellifera) is plastic and responsive to pheromones produced by broods and the queen. Egg laying workers (ELWs), which could reactivate their ovaries and lay haploid eggs upon queen lost, have been commonly discussed from many aspects. However, it remains unclear whether midgut homeostasis in ELWs is affected during plastic changes. Here, we found that the expression of nutrition- and autophagy-related genes was up-regulated in the midguts of ELWs, compared with that in nurse workers (NWs) by RNA-sequencing. Furthermore, the area and number of autophagosomes were increased, along with significantly increased cell death in the midguts of ELWs. Moreover, cell cycle progression in the midguts of ELWs was increased compared with that in NWs. Consistent with the up-regulation of nutrition-related genes, the body and midgut sizes, and the number of intestinal proliferation cells of larvae reared with royal jelly (RJ) obviously increased more than those reared without RJ in vitro. Finally, cell proliferation was dramatically suppressed in the midguts of ELWs when autophagy was inhibited. Altogether, our data suggested that autophagy was induced and required to sustain cell proliferation in ELWs' midguts, thereby revealing the critical role of autophagy played in the intestines during phenotypic plasticity changes.

Epigenetic biomarkers for animal welfare monitoring

Biomarkers for holistic animal welfare monitoring represent a considerable unmet need in veterinary medicine. Epigenetic modifications, like DNA methylation, provide important information about cellular states and environments, which makes them highly attractive for biomarker development. Up until now, much of the corresponding research has been focused on human cancers. However, the increasing availability of animal genomes and epigenomes has greatly improved our capacity for epigenetic biomarker development. In this review, we provide an overview about animal DNA methylation patterns and the technologies that enable the analysis of these patterns. We also describe the key frameworks for compound DNA methylation biomarkers, DNA methylation clocks and environment-specific DNA methylation signatures, that allow complex, context-dependent readouts about animal health and disease. Finally, we provide practical examples for how these biomarkers could be applied for health and environmental exposure monitoring, two key aspects of animal welfare assessments. Taken together, our article provides an overview about the molecular and biological foundations for the development of epigenetic biomarkers in veterinary science and their application potential in animal welfare monitoring.

Databases and prospects of dynamic gene regulation in eukaryotes: A mini review

In eukaryotes, dynamic regulation enables DNA polymerases to catalyze a variety of RNA products in spatial and temporal patterns. Dynamic gene expression is regulated by transcription factors (TFs) and epigenetics (DNA methylation and histone modification). The applications of biochemical technology and high-throughput sequencing enhance the understanding of mechanisms of these regulations and affected genomic regions. To provide a searchable platform for retrieving such metadata, numerous databases have been developed based on the integration of genome-wide maps (e.g., ChIP-seq, whole-genome bisulfite sequencing, RNA-seq, ATAC-seq, DNase-seq, and MNase-seq data) and functionally genomic annotation. In this mini review, we summarize the main functions of TF-related databases and outline the prevalent approaches used in inferring epigenetic regulations, their associated genes, and functions. We review the literature on crosstalk between TF and epigenetic regulation and the properties of non-coding RNA regulation, which are challenging topics that promise to pave the way for advances in database development.

ChIP-seq profiling of H3K4me3 and H3K27me3 in an invasive insect, Bactrocera dorsalis

Introduction: While it has been suggested that histone modifications can facilitate animal responses to rapidly changing environments, few studies have profiled whole-genome histone modification patterns in invasive species, leaving the regulatory landscape of histone modifications in invasive species unclear. Methods: Here, we screen genome-wide patterns of two important histone modifications, trimethylated Histone H3 Lysine 4 (H3K4me3) and trimethylated Histone H3 Lysine 27 (H3K27me3), in adult thorax muscles of a notorious invasive pest, the Oriental fruit fly Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), using Chromatin Immunoprecipitation with high-throughput sequencing (ChIP-seq). Results: We identified promoters featured by the occupancy of H3K4me3, H3K27me3 or bivalent histone modifications that were respectively annotated with unique genes key to muscle development and structure maintenance. In addition, we found H3K27me3 occupied the entire body of genes, where the average enrichment was almost constant. Transcriptomic analysis indicated that H3K4me3 is associated with active gene transcription, and H3K27me3 is mostly associated with transcriptional repression. Importantly, we identified genes and putative motifs modified by distinct histone modification patterns that may possibly regulate flight activity. Discussion: These findings provide the first evidence of histone modification signature in B. dorsalis, and will be useful for future studies of epigenetic signature in other invasive insect species.

Single nucleotide polymorphism (SNP) in the doublesex (dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex (dsx) gene, involved in mosquitos’ somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene (Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

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.

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

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

Comparative Methylome Analysis Reveals Epigenetic Signatures Associated with Growth and Shell Color in the Pacific Oyster, Crassostrea gigas

Fast growth is one of the most important breeding goals for all economic species such as the Pacific oyster (Crassostrea gigas), an aquaculture mollusk with top global production. Although the genetic basis and molecular mechanisms of growth-related traits have been widely investigated in the oyster, the role of DNA methylation involved in growth regulation remains largely unexplored. In this study, we performed a comparative DNA methylome analysis of two selectively bred C. gigas strains with contrasted phenotypes in growth and shell color based on whole-genome bisulfite sequencing (WGBS). Genome-wide profiling of DNA methylation at the single-base resolution revealed that DNA methylations were widely spread across the genome with obvious hotspots, coinciding with the distribution of genes and repetitive elements. Higher methylation levels were observed within genic regions compared with intergenic and promoter regions. Comparative analysis of DNA methylation allowed the identification of 339,604 differentially methylated CpG sites (DMCs) clustering in 27,600 differentially methylated regions (DMRs). Functional annotation analysis identified 11,033 genes from DMRs which were enriched in biological processes including cytoskeleton system, cell cycle, signal transduction, and protein biosynthesis. Integrative analysis of methylome and transcriptome profiles revealed a positive correlation between gene expression and DNA methylation within gene-body regions. Protein-protein interaction (PPI) analysis of differentially expressed and methylated genes allowed for the identification of integrin beta-6 (homolog of human ITGB3) as a hub modulator of the PI3K/Akt signaling pathway that was involved in various growth-related processes. This work provided insights into epigenetic regulation of growth in oysters and will be valuable resources for studying DNA methylation in invertebrates.

Environmental Epigenetics in Soil Ecosystems: Earthworms as Model Organisms

One of the major emerging concerns within ecotoxicology is the effect of environmental pollutants on epigenetic changes, including DNA methylation, histone modifications, and non-coding RNAs. Epigenetic mechanisms regulate gene expression, meaning that the alterations of epigenetic marks can induce long-term physiological effects that can even be inherited across generations. Many invertebrate species have been used as models in environmental epigenetics, with a special focus on DNA methylation changes caused by environmental perturbations (e.g., pollution). Among soil organisms, earthworms are considered the most relevant sentinel organisms for anthropogenic stress assessment and are widely used as standard models in ecotoxicological testing of soil toxicity. In the last decade, several research groups have focused on assessing the impact of environmental stress on earthworm epigenetic mechanisms and tried to link these mechanisms to the physiological effects. The aim of this review is to give an overview and to critically examine the available literature covering this topic. The high level of earthworm genome methylation for an invertebrate species, responsiveness of epigenome to environmental stimuli, availability of molecular resources, and the possibility to study epigenetic inheritance make earthworms adequate models in environmental epigenomics. However, there are still many knowledge gaps that need to be filled in, before we can fully explore earthworms as models in this field. These include detailed characterization of the methylome using next-generation sequencing tools, exploration of multigenerational and transgenerational effects of pollutants, and information about other epigenetic mechanisms apart from DNA methylation. Moreover, the connection between epigenetic effects and phenotype has to be further explored.

Single-cell transcriptomic analysis of honeybee brains identifies vitellogenin as caste differentiation-related factor

The honeybee (Apis mellifera) is a well-known eusocial insect. In honeybee colonies, thousands of sterile workers, including nurse and forager bees, perform various tasks within or outside the hive, respectively. The queen is the only fertile female and is responsible for reproduction. The queen and workers share similar genomes but occupy different caste statuses. We established single-cell transcriptomic atlases of brains from queens and worker subcastes and identified five major cell groups: Kenyon, optic lobe, olfactory projection, glial, and hemocyte cells. By dividing Kenyon and glial cells into multiple subtypes based on credible markers, we observed that vitellogenin (vg) was highly expressed in specific glial-cell subtypes in brains of queens. Knockdown of vg at the early larval stage significantly suppressed the development into adult queens. We demonstrate vg expression as a "molecular signature" for the queen caste and suggest involvement of vg in regulating caste differentiation.

•

scRNA-seq revealed distinct gene expression in the brains of queens and workers

•

Vitellogenin (vg) may represent a "molecular signature" of the queen caste

•

Knockdown of vg at early larval stage suppressed development into adult queens

•

Vg may be involved in regulating caste differentiation

Complex regulatory role of DNA methylation in caste- and age-specific expression of a termite

The reproductive castes of eusocial insects are often characterized by extreme lifespans and reproductive output, indicating an absence of the fecundity/longevity trade-off. The role of DNA methylation in the regulation of caste- and age-specific gene expression in eusocial insects is controversial. While some studies find a clear link to caste formation in honeybees and ants, others find no correlation when replication is increased across independent colonies. Although recent studies have identified transcription patterns involved in the maintenance of high reproduction throughout the long lives of queens, the role of DNA methylation in the regulation of these genes is unknown. We carried out a comparative analysis of DNA methylation in the regulation of caste-specific transcription and its importance for the regulation of fertility and longevity in queens of the higher termite Macrotermes natalensis. We found evidence for significant, well-regulated changes in DNA methylation in mature compared to young queens, especially in several genes related to ageing and fecundity in mature queens. We also found a strong link between methylation and caste-specific alternative splicing. This study reveals a complex regulatory role of fat body DNA methylation both in the division of labour in termites, and during the reproductive maturation of queens.

Encyclopaedia of eukaryotic DNA methylation: from patterns to mechanisms and functions

DNA methylation is an epigenetic modification with a very long evolutionary history. However, DNA methylation evolves surprisingly rapidly across eukaryotes. The genome-wide distribution of methylation diversifies rapidly in different lineages, and DNA methylation is lost altogether surprisingly frequently. The growing availability of genomic and epigenomic sequencing across organisms highlights this diversity but also illuminates potential factors that could explain why both the DNA methylation machinery and its genome-wide distribution evolve so rapidly. Key to this are new discoveries about the fitness costs associated with DNA methylation, and new theories about how the fundamental biochemical mechanisms of DNA methylation introduction and maintenance could explain how new genome-wide patterns of methylation evolve.

The Tsetse Metabolic Gambit: Living on Blood by Relying on Symbionts Demands Synchronization

Tsetse flies have socioeconomic significance as the obligate vector of multiple Trypanosoma parasites, the causative agents of Human and Animal African Trypanosomiases. Like many animals subsisting on a limited diet, microbial symbiosis is key to supplementing nutrient deficiencies necessary for metabolic, reproductive, and immune functions. Extensive studies on the microbiota in parallel to tsetse biology have unraveled the many dependencies partners have for one another. But far less is known mechanistically on how products are swapped between partners and how these metabolic exchanges are regulated, especially to address changing physiological needs. More specifically, how do metabolites contributed by one partner get to the right place at the right time and in the right amounts to the other partner? Epigenetics is the study of molecules and mechanisms that regulate the inheritance, gene activity and expression of traits that are not due to DNA sequence alone. The roles that epigenetics provide as a mechanistic link between host phenotype, metabolism and microbiota (both in composition and activity) is relatively unknown and represents a frontier of exploration. Here, we take a closer look at blood feeding insects with emphasis on the tsetse fly, to specifically propose roles for microRNAs (miRNA) and DNA methylation, in maintaining insect-microbiota functional homeostasis. We provide empirical details to addressing these hypotheses and advancing these studies. Deciphering how microbiota and host activity are harmonized may foster multiple applications toward manipulating host health, including identifying novel targets for innovative vector control strategies to counter insidious pests such as tsetse.

DNA Methylation Difference between Female and Male Ussuri Catfish (Pseudobagrus ussuriensis) in Brain and Gonad Tissues

DNA methylation has been found to be involved in sex determination and differentiation in many aquaculture species. The Ussuri catfish (Pseudobagrus ussuriensis) is a popular aquaculture fish in China with high economic value in which male-biased sex dimorphism was observed in terms of body size and body weight. In this study, DNA methylation-sensitive RAD sequencing (Methyl-RAD) was used to explore the epigenetic difference between adult male and female samples in brain and gonad tissues. In brain tissues, 5,442,496 methylated cytosine sites were found and 9.94% of these sites were from symmetric CCGG or CCWGG sites. Among these sites, 321 differential DNA methylation sites (DMSs) in 171 genes were identified, while in gonad tissues, 4,043,053 methylated cytosines sites were found in total and 11.70% of them were from CCGG or CCWGG. Among these sites, 78 differential DNA methylation sites were found which were located in 64 genes. We also found several sex-determination genes among these differential methylated genes, such as amh, gsdf and hsd11b2 in brain tissues and slco3a1, socs2 and trim47 in gonad tissues. These results provided evidence for understanding the function of DNA methylation in the sex differentiation in Pseudobagrus ussuriensis, which further deepens the relationship between gene regulation and epigenetics.

Paradigm shifts in animal epigenetics: Research on non-model species leads to new insights into dependencies, functions and inheritance of DNA methylation

Recent investigations with non-model species and whole-genome approaches have challenged several paradigms in animal epigenetics. They revealed that epigenetic variation in populations is not the mere consequence of genetic variation, but is a semi-independent or independent source of phenotypic variation, depending on mode of reproduction. DNA methylation is not positively correlated with genome size and phylogenetic position as earlier believed, but has evolved differently between and within higher taxa. Epigenetic marks are usually not completely erased in the zygote and germ cells as generalized from mouse, but often persist and can be transgenerationally inherited, making them evolutionarily relevant. Gene body methylation and promoter methylation are similar in vertebrates and invertebrates with well methylated genomes but transposon silencing through methylation is variable. The new data also suggest that animals use epigenetic mechanisms to cope with rapid environmental changes and to adapt to new environments. The main benefiters are asexual populations, invaders, sessile taxa and long-lived species.

Dnmt1a is essential for gene body methylation and the regulation of the zygotic genome in a wasp

Gene body methylation (GBM) is an ancestral mode of DNA methylation whose role in development has been obscured by the more prominent roles of promoter and CpG island methylation. The wasp Nasonia vitripennis has little promoter and CpG island methylation, yet retains strong GBM, making it an excellent model for elucidating the roles of GBM. Here we show that N. vitripennis DNA methyltransferase 1a (Nv-Dnmt1a) knockdown leads to failures in cellularization and gastrulation of the embryo. Both of these disrupted events are hallmarks of the maternal-zygotic transition (MZT) in insects. Analysis of the embryonic transcriptome and methylome revealed strong reduction of GBM and widespread disruption of gene expression during embryogenesis after Nv-Dnmt1a knockdown. Strikingly, there was a strong correlation between loss of GBM and reduced gene expression in thousands of methylated loci, consistent with the hypothesis that GBM directly facilitates high levels of transcription. We propose that lower expression levels of methylated genes due to reduced GBM is the crucial direct effect of Nv-Dnmt1 knockdown. Subsequently, the disruption of methylated genes leads to downstream dysregulation of the MZT, culminating in developmental failure at gastrulation.

DNA methylation differences between male and female gonads of the oyster reveal the role of epigenetics in sex determination

DNA methylation involved in sex determination mechanism by regulating gene expression related to sex determination networks are common in vertebrates. However, the mechanism linking epigenetics in invertebrates and sex determination has remained elusive. Here, methylome of the male and female gonads in the oyster Crassostrea gigas were conducted to explore the role of epigenetics in invertebrate sex determination. Comparative analysis of gonadal DNA methylation of females and males revealed that male gonads displayed a higher level of DNA methylation and a greater number of hypermethylated genes. Luxury genes presented hypomethylation, while housekeeping genes got hypermethylation. Genes in the conserved signaling pathways, rather than the key master genes in the sex determination pathway, were the major targets of substantial DNA methylation modification. The negative correlation of expression and promoter methylation in the diacylglycerol kinase delta gene (Dgkd) - a ubiquitously expressed gene - indicated DNA methylation may fine turn the expression of Dgkd and be involved in the process of sex determination. Dgkd can be used as an epigenetic marker to distinguish male C. gigas based on the different methylation regions in the promoter region. The results suggest that DNA methylation mechanisms played potential functional impacts in the sex determination in oysters, which is helpful to deepen the understanding of sex determination in invertebrate.

Common mechanisms cannot explain time- and dose-dependent DNA methylation changes in earthworms exposed to cadmium

DNA hypermethylation caused by environmental pollutants like cadmium (Cd) has already been demonstrated in many invertebrates, including earthworms. However, the exact epigenetic mechanisms that drive this hypermethylation are largely unknown and even basic DNA methylation and demethylation processes are hardly characterized. Therefore, we used an important bioindicator, the earthworm Lumbricus terrestris, as a model organism to determine time- and dose-dependent effects of Cd on global and gene-specific DNA methylation and its underlying mechanisms. We revealed Cd-induced adenine and cytosine hypermethylation using specific antibodies in dot blots and found that the methylation level of adenine compared to cytosine changed even to a bigger extent. However, the levels of hydroxymethylated cytosine did not differ between treatment groups. General methylation and demethylation components like methyltransferases (DNMT1 and 3), and ten-eleven translocation (TET) genes were confirmed in L. terrestris by quantitative RealTime PCR. However, neither gene expression, nor DNMT and TET enzyme activity showed significant differences in the Cd exposure groups. Using bisulfite conversion and sequencing, gene body methylation (gbm) of metallothionein 2 (MT2), one of the most important detoxification proteins, was characterized. Cd-dependent changes in MT2 gbm could, however, not be correlated to MT2 gene activity evaluated by quantitative RealTime PCR. Future directions as well as missing links are discussed in the present study hinting towards the importance of studying epigenetic marks and mechanistic insights in a broad variety of species to deepen our knowledge on the effects of changing environmental conditions.