Evolutionary Transition of Promoter and Gene Body DNA Methylation across Invertebrate-Vertebrate Boundary

The evolutionary consequences of interactions between the epigenome, the genome and the environment

The epigenome is the suite of interacting chemical marks and molecules that helps to shape patterns of development, phenotypic plasticity and gene regulation, in part due to its responsiveness to environmental stimuli. There is increasing interest in understanding the functional and evolutionary importance of this sensitivity under ecologically realistic conditions. Observations that epigenetic variation abounds in natural populations have prompted speculation that it may facilitate evolutionary responses to rapid environmental perturbations, such as those occurring under climate change. A frequent point of contention is whether epigenetic variants reflect genetic variation or are independent of it. The genome and epigenome often appear tightly linked and interdependent. While many epigenetic changes are genetically determined, the converse is also true, with DNA sequence changes influenced by the presence of epigenetic marks. Understanding how the epigenome, genome and environment interact with one another is therefore an essential step in explaining the broader evolutionary consequences of epigenomic variation. Drawing on results from experimental and comparative studies carried out in diverse plant and animal species, we synthesize our current understanding of how these factors interact to shape phenotypic variation in natural populations, with a focus on identifying similarities and differences between taxonomic groups. We describe the main components of the epigenome and how they vary within and between taxa. We review how variation in the epigenome interacts with genetic features and environmental determinants, with a focus on the role of transposable elements (TEs) in integrating the epigenome, genome and environment. And we look at recent studies investigating the functional and evolutionary consequences of these interactions. Although epigenetic differentiation in nature is likely often a result of drift or selection on stochastic epimutations, there is growing evidence that a significant fraction of it can be stably inherited and could therefore contribute to evolution independently of genetic change.

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.

DNA Methylation Machinery in Gastropod Mollusks

The role of DNA methylation in mollusks is just beginning to be understood. This review synthesizes current knowledge on this potent molecular hallmark of epigenetic control in gastropods-the largest class of mollusks and ubiquitous inhabitants of diverse habitats. Their DNA methylation machinery shows a high degree of conservation in CG maintenance methylation mechanisms, driven mainly by DNMT1 homologues, and the presence of MBD2 and MBD2/3 proteins as DNA methylation readers. The mosaic-like DNA methylation landscape occurs mainly in a CG context and is primarily confined to gene bodies and housekeeping genes. DNA methylation emerges as a critical regulator of reproduction, development, and adaptation, with tissue-specific patterns being observed in gonadal structures. Its dynamics also serve as an important regulatory mechanism underlying learning and memory processes. DNA methylation can be affected by various environmental stimuli, including as pathogens and abiotic stresses, potentially impacting phenotypic variation and population diversity. Overall, the features of DNA methylation in gastropods are complex, being an essential part of their epigenome. However, comprehensive studies integrating developmental stages, tissues, and environmental conditions, functional annotation of methylated regions, and integrated genomic-epigenomic analyses are lacking. Addressing these knowledge gaps will advance our understanding of gastropod biology, ecology, and evolution.

Quantitative proteome dynamics across embryogenesis in a model chordate

The evolution of gene expression programs underlying the development of vertebrates remains poorly characterized. Here, we present a comprehensive proteome atlas of the model chordate Ciona, covering eight developmental stages and ∼7,000 translated genes, accompanied by a multi-omics analysis of co-evolution with the vertebrate Xenopus. Quantitative proteome comparisons argue against the widely held hourglass model, based solely on transcriptomic profiles, whereby peak conservation is observed during mid-developmental stages. Our analysis reveals maximal divergence at these stages, particularly gastrulation and neurulation. Together, our work provides a valuable resource for evaluating conservation and divergence of multi-omics profiles underlying the diversification of vertebrates.

The role of DNA methylation reprogramming during sex determination and sex reversal in the Pacific oyster Crassostrea gigas

DNA methylation is instrumental in vertebrate sex reversal. However, the mechanism of DNA methylation regulation regarding sex reversal in invertebrates is unclear. In this study, we used whole genome bisulfite sequencing (WGBS) to map single-base resolution methylation profiles of the Pacific oyster, including female-to-male (FMa-to-FMb) and male-to-female (MFa-to-MFb) sex reversal, as well as sex non-reversed males and females (MMa-to-MMb and FFa-to-FFb). The results showed that global DNA methylation levels increase during female-to-male sex reversals, with a particular increase in the proportion of high methylation levels (mCGs >0.75) and a decrease in the proportion of intermediate methylation levels (0.25 < mCGs <0.75). This increase in DNA methylation was mainly associated with the elevated expression of DNA methylase genes. Genome-wide methylation patterns of females were accurately remodeled to those of males after sex reversal, while the opposite was true for the male-to-female reversal. Those findings directly indicate that alterations in DNA methylation play a significant role in sex reversal in Pacific oysters. Comparative analysis of the DNA methylomes of pre- and post- sex reversal gonadal tissues (FMb-vs-FMa or MFb-vs-MFa) revealed that differentially methylated genes were mainly involved in the biological processes of sex determination or gonadal development. However critical genes such as Dmrt1, Foxl2 and Sox-like, which are involved in the putative sex determination pathway in Pacific oysters, showed almost an absence of methylation modifications, varying greatly from vertebrates. Additionally, comparative analysis of the DNA methylomes of sexual reversal and sex non-reversal (FMa-vs-FFa or MFa-vs-MMa) revealed that heat shock protein genes, such as Hsp68-like and Hsp70B, were important for the occurrence of sex reversal. These findings shed light on the epigenetic mechanisms underlying the maintenance of gonadal plasticity and the reversal of organ architecture in oysters.

Potential Role of DNA Methylation as a Driver of Plastic Responses to the Environment Across Cells, Organisms, and Populations

There is great interest in exploring epigenetic modifications as drivers of adaptive organismal responses to environmental change. Extending this hypothesis to populations, epigenetically driven plasticity could influence phenotypic changes across environments. The canonical model posits that epigenetic modifications alter gene regulation and subsequently impact phenotypes. We first discuss origins of epigenetic variation in nature, which may arise from genetic variation, spontaneous epimutations, epigenetic drift, or variation in epigenetic capacitors. We then review and synthesize literature addressing three facets of the aforementioned model: (i) causal effects of epigenetic modifications on phenotypic plasticity at the organismal level, (ii) divergence of epigenetic patterns in natural populations distributed across environmental gradients, and (iii) the relationship between environmentally induced epigenetic changes and gene expression at the molecular level. We focus on DNA methylation, the most extensively studied epigenetic modification. We find support for environmentally associated epigenetic structure in populations and selection on stable epigenetic variants, and that inhibition of epigenetic enzymes frequently bears causal effects on plasticity. However, there are pervasive confounding issues in the literature. Effects of chromatin-modifying enzymes on phenotype may be independent of epigenetic marks, alternatively resulting from functions and protein interactions extrinsic of epigenetics. Associations between environmentally induced changes in DNA methylation and expression are strong in plants and mammals but notably absent in invertebrates and nonmammalian vertebrates. Given these challenges, we describe emerging approaches to better investigate how epigenetic modifications affect gene regulation, phenotypic plasticity, and divergence among populations.

(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.

Analysis of DNA Methylation Differences during the JIII Formation of Bursaphelenchus xylophilus

DNA methylation is a pivotal process that regulates gene expression and facilitates rapid adaptation to challenging environments. The pinewood nematode (PWN; Bursaphelenchus xylophilus), the causative agent of pine wilt disease, survives at low temperatures through third-stage dispersal juvenile, making it a major pathogen for pines in Asia. To comprehend the impact of DNA methylation on the formation and environmental adaptation of third-stage dispersal juvenile, we conducted whole-genome bisulfite sequencing and transcriptional sequencing on both the third-stage dispersal juvenile and three other propagative juvenile stages of PWN. Our findings revealed that the average methylation rate of cytosine in the samples ranged from 0.89% to 0.99%. Moreover, we observed significant DNA methylation changes in the third-stage dispersal juvenile and the second-stage propagative juvenile of PWN, including differentially methylated cytosine (DMCs, n = 435) and regions (DMRs, n = 72). In the joint analysis of methylation-associated transcription, we observed that 23 genes exhibited overlap between differentially methylated regions and differential gene expression during the formation of the third-stage dispersal juvenile of PWN. Further functional analysis of these genes revealed enrichment in processes related to lipid metabolism and fatty acid synthesis. These findings emphasize the significance of DNA methylation in the development of third-stage dispersal juvenile of PWN, as it regulates transcription to enhance the probability of rapid expansion in PWN.

Glucocorticoid receptor regulates protein chaperone, circadian clock and affective disorder genes in the zebrafish brain

Glucocorticoid resistance is commonly observed in depression, and has been linked to reduced expression and/or function of the glucocorticoid receptor (NR3C1 in human, hereafter referred to as GR). Previous studies have shown that GR-mutant zebrafish exhibit behavioural abnormalities that are indicative of an affective disorder, suggesting that GR plays a role in brain function. We compared the brain methylomes and brain transcriptomes of adult wild-type and GR-mutant zebrafish, and identified 249 differentially methylated regions (DMRs) that are regulated by GR. These include a cluster of CpG sites within the first intron of fkbp5, the gene encoding the glucocorticoid-inducible heat shock protein co-chaperone Fkbp5. RNA-sequencing analysis revealed that genes associated with chaperone-mediated protein folding, the regulation of circadian rhythm and the regulation of metabolism are particularly sensitive to loss of GR function. In addition, we identified subsets of genes exhibiting GR-regulated transcription that are known to regulate behaviour, and are linked to unipolar depression and anxiety. Taken together, our results identify key biological processes and novel molecular mechanisms through which the GR is likely to mediate responses to stress in the adult zebrafish brain, and they provide further support for the zebrafish GR mutant as a model for the study of affective disorders.

Evolution of termination codons of proteins and the TAG-TGA paradox

In most eukaryotes and prokaryotes TGA is used at a significantly higher frequency than TAG as termination codon of protein-coding genes. Although this phenomenon has been recognized several years ago, there is no generally accepted explanation for the TAG-TGA paradox. Our analyses of human mutation data revealed that out of the eighteen sense codons that can give rise to a nonsense codon by single base substitution, the CGA codon is exceptional: it gives rise to the TGA stop codon at an order of magnitude higher rate than the other codons. Here we propose that the TAG-TGA paradox is due to methylation and hypermutabilty of CpG dinucleotides. In harmony with this explanation, we show that the coding genomes of organisms with strong CpG methylation have a significant bias for TGA whereas those from organisms that lack CpG methylation use TGA and TAG termination codons with similar probability.

High temperature influences DNA methylation and transcriptional profiles in sea urchins (Strongylocentrotus intermedius)

BACKGROUND

DNA methylation plays an important role in life processes by affecting gene expression, but it is still unclear how DNA methylation is controlled and how it regulates gene transcription under high temperature stress conditions in Strongylocentrotus intermedius. The potential link between DNA methylation variation and gene expression changes in response to heat stress in S. intermedius was investigated by MethylRAD-seq and RNA-seq analysis. We screened DNA methylation driver genes in order to comprehensively elucidate the regulatory mechanism of its high temperature adaptation at the DNA/RNA level.

RESULTS

The results revealed that high temperature stress significantly affected not only the DNA methylation and transcriptome levels of S. intermedius (P < 0.05), but also growth. MethylRAD-seq analysis revealed 12,129 CG differential methylation sites and 966 CWG differential methylation sites, and identified a total of 189 differentially CG methylated genes and 148 differentially CWG methylated genes. Based on KEGG enrichment analysis, differentially expressed genes (DEGs) are mostly enriched in energy and cell division, immune, and neurological damage pathways. Further RNA-seq analysis identified a total of 1968 DEGs, of which 813 genes were upregulated and 1155 genes were downregulated. Based on the joint MethylRAD-seq and RNA-seq analysis, metabolic processes such as glycosaminoglycan degradation, oxidative phosphorylation, apoptosis, glutathione metabolism, thermogenesis, and lysosomes are regulated by DNA methylation.

CONCLUSIONS

High temperature affected the DNA methylation and expression levels of genes such as MOAP-1, GGT1 and RDH8, which in turn affects the metabolism of HPSE, Cox, glutathione, and retinol, thereby suppressing the immune, energy metabolism, and antioxidant functions of the organism and finally manifesting as stunted growth. In summary, the observations in the present study improve our understanding of the molecular mechanism of the response to high temperature stress in sea urchin.

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.

Associations between DNA methylation and gene regulation depend on chromatin accessibility during transgenerational plasticity

BACKGROUND

Epigenetic processes are proposed to be a mechanism regulating gene expression during phenotypic plasticity. However, environmentally induced changes in DNA methylation exhibit little-to-no association with differential gene expression in metazoans at a transcriptome-wide level. It remains unexplored whether associations between environmentally induced differential methylation and expression are contingent upon other epigenomic processes such as chromatin accessibility. We quantified methylation and gene expression in larvae of the purple sea urchin Strongylocentrotus purpuratus exposed to different ecologically relevant conditions during gametogenesis (maternal conditioning) and modeled changes in gene expression and splicing resulting from maternal conditioning as functions of differential methylation, incorporating covariates for genomic features and chromatin accessibility. We detected significant interactions between differential methylation, chromatin accessibility, and genic feature type associated with differential expression and splicing.

RESULTS

Differential gene body methylation had significantly stronger effects on expression among genes with poorly accessible transcriptional start sites while baseline transcript abundance influenced the direction of this effect. Transcriptional responses to maternal conditioning were 4-13 × more likely when accounting for interactions between methylation and chromatin accessibility, demonstrating that the relationship between differential methylation and gene regulation is partially explained by chromatin state.

CONCLUSIONS

DNA methylation likely possesses multiple associations with gene regulation during transgenerational plasticity in S. purpuratus and potentially other metazoans, but its effects are dependent on chromatin accessibility and underlying genic features.

Starvation-induced changes in sex ratio involve alterations in sex-related gene expression and methylation in Pacific oyster Crassostrea gigas

Aquatic animals are subject to varying degrees of starvation stress in their natural habitats due to food limitations. Consequently, starvation is a crucial environmental factor for sex determination in many species; however, limited research has been conducted on the effects of starvation on sex determination in shellfish. Here, four full sibling families of Pacific oysters were established and subjected to starvation stress. The results demonstrated that starvation caused the sex ratio (female to male) to change from 1:0.78 to 1:1.44 and resulted in a delay in gonadal development. Further studies revealed that the expression levels of DNA methylation-related genes Dnmt1 (DNA methyltransferase 1), Dnmt3a/b (DNA methyltransferase 3a/b) and Tet3 (tet methylcytosine dioxygenase 3) decreased under starvation stress. Conversely, the upregulation of Dmrt1 (doublesex and mab-3 related transcription factor 1), a gene typically associated with males, in females suggests that the increased proportion of males may be linked to starvation-induced high expression of this particular gene. In addition, the gene Dgkd (diacylglycerol kinase delta), which is involved in the regulation of second messenger protein kinase C, was differentially methylated between males and females, with the methylation level of this gene gradually increasing with male development, while the methylation level of this gene decreased under starvation stress. Correlation analysis of Dgkd methylation levels with expression levels showed a negative correlation between DNA methylation and gene expression. Finally dual fluorescence reporter experiments confirmed that DNA methylation suppressed Dgkd expression in vitro. In summary, the results suggest that starvation may regulate Dgkd gene expression through DNA methylation and thus affect Dmrt1 expression, thereby causing sex reversal in the Pacific oyster. The outcomes resolved how environmental factors are involved in sex reversal from an epigenetic perspective and provided a theoretical basis for sex control breeding in the Pacific oyster.

Evolution of homology: From archetype towards a holistic concept of cell type

The concept of homology lies in the heart of comparative biological science. The distinction between homology as structure and analogy as function has shaped the evolutionary paradigm for a century and formed the axis of comparative anatomy and embryology, which accept the identity of structure as a ground measure of relatedness. The advent of single-cell genomics overturned the classical view of cell homology by establishing a backbone regulatory identity of cell types, the basic biological units bridging the molecular and phenotypic dimensions, to reveal that the cell is the most flexible unit of living matter and that many approaches of classical biology need to be revised to understand evolution and diversity at the cellular level. The emerging theory of cell types explicitly decouples cell identity from phenotype, essentially allowing for the divergence of evolutionarily related morphotypes beyond recognition, as well as it decouples ontogenetic cell lineage from cell-type phylogeny, whereby explicating that cell types can share common descent regardless of their structure, function or developmental origin. The article succinctly summarizes current progress and opinion in this field and formulates a more generalistic view of biological cell types as avatars, transient or terminal cell states deployed in a continuum of states by the developmental programme of one and the same omnipotent cell, capable of changing or combining identities with distinct evolutionary histories or inventing ad hoc identities that never existed in evolution or development. It highlights how the new logic grounded in the regulatory nature of cell identity transforms the concepts of cell homology and phenotypic stability, suggesting that cellular evolution is inherently and massively network-like, with one-to-one homologies being rather uncommon and restricted to shallower levels of the animal tree of life.

Comparison of genome-wide DNA methylation patterns between antler precartilage and cartilage

Deer antlers are the only mammalian organs that can fully regenerate after being lost and provide a valuable model for cartilage development. As one of the best-studied epigenetic mechanisms, DNA methylation is known to engage in organ and tissue development. This study aimed to investigate the role of DNA methylation in antler chondrogenesis by comparing whole-genome DNA methylation between precartilage and cartilage. Quantitative reverse transcription PCR (RT-qPCR) showed significant differences in the expression levels of DNA methyltransferase genes (DNMT1, DNMT3A, and DNMT3B) between precartilage and cartilage. Subsequently, we obtained DNA methylation profiles of antler precartilage and cartilage tissues by whole-genome bisulfite sequencing. Although sequencing data indicated that overall methylation levels at CpG and non-CpG sites were similar between precartilage and cartilage, 140,784 differentially methylated regions (DMRs, P < 0.05) and 3,941 DMR-related genes were identified. Gene ontology (GO) analysis of DMR-related genes demonstrated some significantly enriched GO terms (P < 0.05) related to chondrogenesis, including insulin receptor binding, collage trimer, integrin binding, and extracellular matrix structural constituent. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DMR-related genes uncovered that the PI3K/AKT, cortisol synthesis and secretion, glycosaminoglycan biosynthesis-keratan sulfate, Hippo, and NF-κB signaling pathways might play a pivotal role in the transition of precartilage to cartilage. Moreover, we found that 25 DMR-related genes, including CD44, IGF1, ITGAV, ITGB1, RUNX1, COL2A1, COMP, and TAGLN, were most likely involved in antler chondrogenesis. In conclusion, this study revealed the genome-wide DNA methylation patterns of antler precartilage and cartilage, which may contribute to understanding the epigenetic regulation of antler chondrogenesis.

Synonymous alterations of cancer-associated Trp53 CpG mutational hotspots cause fatal developmental jaw malocclusions but no tumors in knock-in mice

Intragenic CpG dinucleotides are tightly conserved in evolution yet are also vulnerable to methylation-dependent mutation, raising the question as to why these functionally critical sites have not been deselected by more stable coding sequences. We previously showed in cell lines that altered exonic CpG methylation can modify promoter start sites, and hence protein isoform expression, for the human TP53 tumor suppressor gene. Here we extend this work to the in vivo setting by testing whether synonymous germline modifications of exonic CpG sites affect murine development, fertility, longevity, or cancer incidence. We substituted the DNA-binding exons 5-8 of Trp53, the mouse ortholog of human TP53, with variant-CpG (either CpG-depleted or -enriched) sequences predicted to encode the normal p53 amino acid sequence; a control construct was also created in which all non-CpG sites were synonymously substituted. Homozygous Trp53-null mice were the only genotype to develop tumors. Mice with variant-CpG Trp53 sequences remained tumor-free, but were uniquely prone to dental anomalies causing jaw malocclusion (p < .0001). Since the latter phenotype also characterises murine Rett syndrome due to dysfunction of the trans-repressive MeCP2 methyl-CpG-binding protein, we hypothesise that CpG sites may exert non-coding phenotypic effects via pre-translational cis-interactions of 5-methylcytosine with methyl-binding proteins which regulate mRNA transcript initiation, expression or splicing, although direct effects on mRNA structure or translation are also possible.

Epigenetics across the evolutionary tree: New paradigms from non-model animals

All animals have evolved solutions to manage their genomes, enabling the efficient organization of meters of DNA strands in the nucleus and allowing for nuanced regulation of gene expression while keeping transposable elements suppressed. Epigenetic modifications are central to accomplishing all these. Recent advances in sequencing technologies and the development of techniques that profile epigenetic marks and chromatin accessibility using reagents that can be used in any species has catapulted epigenomic studies in diverse animal species, shedding light on the multitude of epigenomic mechanisms utilized across the evolutionary tree. Now, comparative epigenomics is a rapidly growing field that is uncovering mechanistic aspects of epigenetic modifications and chromatin organization in non-model invertebrates, ranging from octopus to sponges. This review puts recent discoveries in the epigenetics of non-model invertebrates in historical context, and describes new insight into the patterning and functions of DNA methylation and other highly conserved epigenetic modifications.

Prognostic value of PNN in prostate cancer and its correlation with therapeutic significance

Prostate cancer (PCa) is the most common malignancy. New biomarkers are in demand to facilitate the management. The role of the pinin protein (encoded by PNN gene) in PCa has not been thoroughly explored yet. Using The Cancer Genome Atlas (TCGA-PCa) dataset validated with Gene Expression Omnibus (GEO) and protein expression data retrieved from the Human Protein Atlas, the prognostic and diagnostic values of PNN were studied. Highly co-expressed genes with PNN (HCEG) were constructed for pathway enrichment analysis and drug prediction. A prognostic signature based on methylation status using HCEG was constructed. Gene set enrichment analysis (GSEA) and the TISIDB database were utilised to analyse the associations between PNN and tumour-infiltrating immune cells. The upregulated PNN expression in PCa at both transcription and protein levels suggests its potential as an independent prognostic factor of PCa. Analyses of the PNN's co-expression network indicated that PNN plays a role in RNA splicing and spliceosomes. The prognostic methylation signature demonstrated good performance for progression-free survival. Finally, our results showed that the PNN gene was involved in splicing-related pathways in PCa and identified as a potential biomarker for PCa.

Adapting to climate with limited genetic diversity: Nucleotide, DNA methylation and microbiome variation among populations of the social spider Stegodyphus dumicola

Understanding the role of genetic and nongenetic variants in modulating phenotypes is central to our knowledge of adaptive responses to local conditions and environmental change, particularly in species with such low population genetic diversity that it is likely to limit their evolutionary potential. A first step towards uncovering the molecular mechanisms underlying population-specific responses to the environment is to carry out environmental association studies. We associated climatic variation with genetic, epigenetic and microbiome variation in populations of a social spider with extremely low standing genetic diversity. We identified genetic variants that are associated strongly with environmental variation, particularly with average temperature, a pattern consistent with local adaptation. Variation in DNA methylation in many genes was strongly correlated with a wide set of climate parameters, thereby revealing a different pattern of associations than that of genetic variants, which show strong correlations to a more restricted range of climate parameters. DNA methylation levels were largely independent of cis-genetic variation and of overall genetic population structure, suggesting that DNA methylation can work as an independent mechanism. Microbiome composition also correlated with environmental variation, but most strong associations were with precipitation-related climatic factors. Our results suggest a role for both genetic and nongenetic mechanisms in shaping phenotypic responses to local environments.

Epigenetic machinery is functionally conserved in cephalopods

Background

Epigenetic regulatory mechanisms are divergent across the animal kingdom, yet these mechanisms are not well studied in non-model organisms. Unique features of cephalopods make them attractive for investigating behavioral, sensory, developmental, and regenerative processes, and recent studies have elucidated novel features of genome organization and gene and transposon regulation in these animals. However, it is not known how epigenetics regulates these interesting cephalopod features. We combined bioinformatic and molecular analysis of Octopus bimaculoides to investigate the presence and pattern of DNA methylation and examined the presence of DNA methylation and 3 histone post-translational modifications across tissues of three cephalopod species.

Results

We report a dynamic expression profile of the genes encoding conserved epigenetic regulators, including DNA methylation maintenance factors in octopus tissues. Levels of 5-methyl-cytosine in multiple tissues of octopus, squid, and bobtail squid were lower compared to vertebrates. Whole genome bisulfite sequencing of two regions of the brain and reduced representation bisulfite sequencing from a hatchling of O. bimaculoides revealed that less than 10% of CpGs are methylated in all samples, with a distinct pattern of 5-methyl-cytosine genome distribution characterized by enrichment in the bodies of a subset of 14,000 genes and absence from transposons. Hypermethylated genes have distinct functions and, strikingly, many showed similar expression levels across tissues while hypomethylated genes were silenced or expressed at low levels. Histone marks H3K27me3, H3K9me3, and H3K4me3 were detected at different levels across tissues of all species.

Conclusions

Our results show that the DNA methylation and histone modification epigenetic machinery is conserved in cephalopods, and that, in octopus, 5-methyl-cytosine does not decorate transposable elements, but is enriched on the gene bodies of highly expressed genes and could cooperate with the histone code to regulate tissue-specific gene expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01404-1.

Sequencing Bait: Nuclear and Mitogenome Assembly of an Abundant Coastal Tropical and Subtropical Fish, Atherinomorus stipes

Genetic data from nonmodel species can inform ecology and physiology, giving insight into a species' distribution and abundance as well as their responses to changing environments, all of which are important for species conservation and management. Moreover, reduced sequencing costs and improved long-read sequencing technology allows researchers to readily generate genomic resources for nonmodel species. Here, we apply Oxford Nanopore long-read sequencing and low-coverage (∼1x) whole genome short-read sequencing technology (Illumina) to assemble a genome and examine population genetics of an abundant tropical and subtropical fish, the hardhead silverside (Atherinomorus stipes). These fish are found in shallow coastal waters and are frequently included in ecological models because they serve as abundant prey for commercially and ecologically important species. Despite their importance in sub-tropical and tropical ecosystems, little is known about their population connectivity and genetic diversity. Our A. stipes genome assembly is about 1.2 Gb with comparable repetitive element content (∼47%), number of protein duplication events, and DNA methylation patterns to other teleost fish species. Among five sampled populations spanning 43 km of South Florida and the Florida Keys, we find little population structure suggesting high population connectivity.

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.

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.

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.

GSDMs are potential therapeutic targets and prognostic biomarkers in clear cell renal cell carcinoma

GSDM family is a group of critical proteins that mediate pyroptosis and plays an important role in cell death and inflammation. However, their specific function in clear cell renal cell carcinoma (ccRCC, KIRC) have not been clarified comprehensively. In this study, we assessed the roles of the GSDM family in expression, prognostic value, functional enrichment analysis, genetic alterations, immune infiltration and DNA methylation in ccRCC patients by using different bioinformatics databases. We found that the expression levels of GADMA-E were significantly higher in ccRCC tissues compared with normal tissues, while the expression level of PJVK was decreased. Moreover, survival analysis indicated that upregulation of GSDME was related to poor overall survival (OS) and recurrence-free survival (RFS) of ccRCC patients. The main function of differentially expressed GSDM homologs was related to ion transport. We also found that the expression profiles of the GSDM family were highly correlated with infiltrating immune cells (i.e., CD8+ T cells, CD4+ T cells, B cells, macrophages, neutrophils and dendritic cells), and there were significant differences in the expression of GSDM family in different ccRCC immune subtypes. Furthermore, DNA methylation analysis indicated that the DNA methylation levels of GSDMA/B/D/E were decreased, while the DNA methylation level of PJVK was increased. In conclusion, this study provides integrated information about abnormal GSDM family members as potential biomarkers for the diagnosis and prognosis of ccRCC. Especially, GSDME was a potential clinical target and prognostic biomarkers for patients with ccRCC.

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.

Hit-and-Run Epigenetic Editing for Vectors of Snail-Borne Parasitic Diseases

Snail-borne parasitic diseases represent an important challenge to human and animal health. Control strategies that target the intermediate snail host has proved very effective. Epigenetic mechanisms are involved in developmental processes and therefore play a fundamental role in developmental variation. DNA methylation is an important epigenetic information carrier in eukaryotes that plays a major role in the control of chromatin structure. Epigenome editing tools have been instrumental to demonstrate functional importance of this mark for gene expression in vertebrates. In invertebrates, such tools are missing, and the role of DNA methylation remains unknown. Here we demonstrate that methylome engineering can be used to modify in vivo the CpG methylation level of a target gene in the freshwater snail Biomphalaria glabrata, intermediate host of the human parasite Schistosoma mansoni. We used a dCas9-SunTag-DNMT3A complex and synthetic sgRNA to transfect B. glabrata embryos and observed an increase of CpG methylation at the target site in 50% of the hatching snails.

Pollution induces epigenetic effects that are stably transmitted across multiple generations

It has been hypothesized that the effects of pollutants on phenotypes can be passed to subsequent generations through epigenetic inheritance, affecting populations long after the removal of a pollutant. But there is still little evidence that pollutants can induce persistent epigenetic effects in animals. Here, we show that low doses of commonly used pollutants induce genome‐wide differences in cytosine methylation in the freshwater crustacean Daphnia pulex. Uniclonal populations were either continually exposed to pollutants or switched to clean water, and methylation was compared to control populations that did not experience pollutant exposure. Although some direct changes to methylation were only present in the continually exposed populations, others were present in both the continually exposed and switched to clean water treatments, suggesting that these modifications had persisted for 7 months (>15 generations). We also identified modifications that were only present in the populations that had switched to clean water, indicating a long‐term legacy of pollutant exposure distinct from the persistent effects. Pollutant‐induced differential methylation tended to occur at sites that were highly methylated in controls. Modifications that were observed in both continually and switched treatments were highly methylated in controls and showed reduced methylation in the treatments. On the other hand, modifications found just in the switched treatment tended to have lower levels of methylation in the controls and showed increase methylation in the switched treatment. In a second experiment, we confirmed that sublethal doses of the same pollutants generate effects on life histories for at least three generations following the removal of the pollutant. Our results demonstrate that even low doses of pollutants can induce transgenerational epigenetic effects that are stably transmitted over many generations. Persistent effects are likely to influence phenotypic development, which could contribute to the rapid adaptation, or extinction, of populations confronted by anthropogenic stressors.

The role of DNA methylation in genome defense in Cnidaria and other invertebrates

Considerable attention has recently been focused on the potential involvement of DNA methylation in regulating gene expression in cnidarians. Much of this work has been centered on corals, in the context of changes in methylation perhaps facilitating adaptation to higher seawater temperatures and other stressful conditions. Although first proposed more than 30 years ago, the possibility that DNA methylation systems function in protecting animal genomes against the harmful effects of transposon activity has largely been ignored since that time. Here, we show that transposons are specifically targeted by the DNA methylation system in cnidarians, and that the youngest transposons (i.e., those most likely to be active) are most highly methylated. Transposons in longer and highly active genes were preferentially methylated and, as transposons aged, methylation levels declined, reducing the potentially harmful side effects of CpG methylation. In Cnidaria and a range of other invertebrates, correlation between the overall extent of methylation and transposon content was strongly supported. Present transposon burden is the dominant factor in determining overall level of genomic methylation in a range of animals that diverged in or before the early Cambrian, suggesting that genome defense represents the ancestral role of CpG methylation.

BTBD10 is a Prognostic Biomarker Correlated With Immune Infiltration in Hepatocellular Carcinoma

Background: BTBD10 serves as an activator of Akt family members through decreasing the protein phosphatase 2A-mediated dephosphorylation. The present study attempted to investigate the prognostic value of BTBD10 in hepatocellular carcinoma (HCC), specially, its relationship with tumor-infiltrating lymphocytes (TILs).

Methods: BTBD10 expression was evaluated in HCC using The Cancer Genome Atlas (TCGA) and Xijing Hospital database, and verified in HCC cell lines. Cox analyses were performed to analyze independent prognostic risk factors for HCC. The optimal cut-off value of BTBD10 was calculated, by which all patients were divided into two groups to compare the overall survival (OS). The signaling pathways were predicted, by which BTBD10 may affect the progression of HCC. To investigate the impact of BTBD10 on HCC immunotherapy, correlations between BTBD10 and TILs, immune checkpoints, m6A methylation-related genes and ferroptosis-related genes were assessed. The distribution of half-maximal inhibitory concentration (IC50) of diverse targeted drugs was observed based on the differential expression of BTBD10.

Results: BTBD10 expression was higher in HCC tissues and cell lines than that of normal liver tissues and cells. The patients with high expression of BTBD10 showed a worse OS, as compared to that of BTBD10 low-expressing group. Cox analyses indicated that BTBD10 was an independent prognostic risk factor for HCC. Several molecular pathways of immune responses were activated in HCC patients with high-expressing of BTBD10. Furthermore, BTBD10 expression was demonstrated to be positively correlated with tumor-infiltrating B cells, T cells, macrophages, neutrophils and dendritic cells. Meanwhile, the expression of BTBD10 was synchronized with that of several m6A methylation-related genes, ferroptosis-related genes and immune checkpoints. The IC50 scores of Sorafenib, Navitoclax, Veliparib, Luminespib, and Imatinib were found to be lower in BTBD10 high-expressing HCC group.

Conclusion: BTBD10 negatively regulates tumor immunity in HCC and exhibits adverse effect on the prognosis of HCC, which could be a potential target for immunotherapy.

OUP accepted manuscript

DNA methylation, an important component of eukaryotic epigenetics, varies in pattern and function across Metazoa. Notably, bilaterian vertebrates and invertebrates differ dramatically in gene body methylation (GbM). Using the frequency of cytosine-phospho-guanines (CpGs), which are lost through mutation when methylated, we report the first broad survey of DNA methylation in Cnidaria, the ancient sister group to Bilateria. We find that: 1) GbM differentially relates to expression categories as it does in most bilaterian invertebrates, but distributions of GbM are less discretely bimodal. 2) Cnidarians generally have lower CpG frequencies on gene bodies than bilaterian invertebrates potentially suggesting a compensatory mechanism to replace CpG lost to mutation in Bilateria that is lacking in Cnidaria. 3) GbM patterns show some consistency within taxonomic groups such as the Scleractinian corals; however, GbM patterns variation across a range of taxonomic ranks in Cnidaria suggests active evolutionary change in GbM within Cnidaria. 4) Some but not all GbM variation is associated with life history change and genome expansion, whereas GbM loss is evident in endoparasitic cnidarians. 5) Cnidarian repetitive elements are less methylated than gene bodies, and methylation of both correlate with genome repeat content. 6) These observations reinforce claims that GbM evolved in stem Metazoa. Thus, this work supports overlap between DNA methylation processes in Cnidaria and Bilateria, provides a framework to compare methylation within and between Cnidaria and Bilateria, and demonstrates the previously unknown rapid evolution of cnidarian methylation.

Effects of Cadmium Sulfate on the Brown Garden Snail Cornu aspersum: Implications for DNA Methylation

An extensive literature exists regarding the cellular, physiological, and genetic effects of cadmium (Cd)—A highly toxic, but commonly used trace metal in modern industry. However, limited data are available on its epigenetic effects, especially for terrestrial sentinel invertebrates. We determined Cd retention, total DNA methylation, and the methylation status of 5′ end of the Cd-MT gene in the hepatopancreas of the brown garden snail, Cornu aspersum, fed Cd sulfate for four weeks. Bodyweight changes and survival were also measured. Hepatopancreas cadmium increased in a dose-dependent manner from the third-lowest dose onward, with very large amounts being found for the highest treatment group. However, no mortalities occurred, irrespective of dietary Cd dose. We identified significant genome-wide hypermethylation in specimens given the highest dose, which overlapped with a significant bodyweight decrease. The Cd-MT gene showed an unmethylated 5′ end of the Cd-MT gene and this status was not affected by cadmium exposure. Hepatopancreas DNA methylation is as sensitive as bodyweight to non-lethal concentrations of dietary Cd given as cadmium sulfate but less responsive than tissue accumulation. Such an exposure event, by contrast, does not affect the methylation status of the Cd-MT gene 5′ end.

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 methylome of Biomphalaria glabrata and other mollusks: enduring modification of epigenetic landscape and phenotypic traits by a new DNA methylation inhibitor

BACKGROUND

5-Methylcytosine (5mC) is an important epigenetic mark in eukaryotes. Little information about its role exists for invertebrates. To investigate the contribution of 5mC to phenotypic variation in invertebrates, alteration of methylation patterns needs to be produced. Here, we apply new non-nucleoside DNA methyltransferase inhibitors (DNMTi) to introduce aleatory changes into the methylome of mollusk species.

RESULTS

Flavanone inhibitor Flv1 was efficient in reducing 5mC in the freshwater snails Biomphalaria glabrata and Physa acuta, and to a lesser degree, probably due to lower stability in sea water, in the oyster Crassostrea gigas. Flv1 has no toxic effects and significantly decreased the 5mC level in the treated B. glabrata and in its offspring. Drug treatment triggers significant variation in the shell height in both generations. A reduced representation bisulfite-sequencing method called epiGBS corroborates hypomethylation effect of Flv1 in both B. glabrata generations and identifies seven Differential Methylated Regions (DMR) out of 32 found both in Flv1-exposed snails and its progeny, from which 5 were hypomethylated, demonstrating a multigenerational effect. By targeted bisulfite sequencing, we confirmed hypomethylation in a locus and show that it is associated with reduced gene expression.

CONCLUSIONS

Flv1 is a new and efficient DNMTi that can be used to induce transient and heritable modifications of the epigenetic landscape and phenotypic traits in mollusks, a phylum of the invertebrates in which epigenetics is understudied.

Evolutionary History of DNA Methylation Related Genes in Bivalvia: New Insights From Mytilus galloprovincialis

DNA methylation is an essential epigenetic mechanism influencing gene expression in all organisms. In metazoans, the pattern of DNA methylation changes during embryogenesis and adult life. Consequently, differentiated cells develop a stable and unique DNA methylation pattern that finely regulates mRNA transcription during development and determines tissue-specific gene expression. Currently, DNA methylation remains poorly investigated in mollusks and completely unexplored in Mytilus galloprovincialis. To shed light on this process in this ecologically and economically important bivalve, we screened its genome, detecting sequences homologous to DNA methyltransferases (DNMTs), methyl-CpG-binding domain (MBD) proteins and Ten-eleven translocation methylcytosine dioxygenase (TET) previously described in other organisms. We characterized the gene architecture and protein domains of the mussel sequences and studied their phylogenetic relationships with the ortholog sequences from other bivalve species. We then comparatively investigated their expression levels across different adult tissues in mussel and other bivalves, using previously published transcriptome datasets. This study provides the first insights on DNA methylation regulators in M. galloprovincialis, which may provide fundamental information to better understand the complex role played by this mechanism in regulating genome activity in bivalves.

Effect of DNA methylation on gene transcription is associated with the distribution of methylation sites across the genome in osteoarthritis

Genetics and epigenetics are important subjects in the field of osteoarthritis (OA) research. DNA methylation may affect gene transcription, but the specific mechanisms have remained to be fully elucidated. In the present study, the ChAMP methylation analysis package was used to identify differentially methylated genes (DMGs) from the dataset GSE63695 from the Gene Expression Omnibus (GEO) database. The distribution of differentially methylated sites (DMS) and the total array sites across the genome were analyzed by enrichment analysis. Subsequently, two mRNA expression profiling datasets, GSE114007 and GSE113825, were obtained from the GEO database and common differentially expressed genes (DEGs) were identified using the Limma package. Key genes were screened by analyzing the distribution of DMS across the genome consisting of DEGs and DMGs. A total of 1,662 and 1,986 DEGs were identified between OA and normal human cartilage from the GSE113825 and GSE114007 dataset, respectively. A further screening revealed 292 genes with common differences between the two datasets. A total of 574 DMS containing 394 DMGs were observed between OA and normal cartilage. Integrative analysis revealed a corresponding subset of 15 genes. Of these, 6 genes were verified by reverse transcription-quantitative PCR, confirming that the mRNA expression of 5 genes (MAP1B, FNDC1, ANLN, SCNN1A and STC2) in OA cartilage was consistent with the mRNA expression from the analysis of the datasets. Upon treatment with the DNA methylation inhibitor 5-aza-2'-deoxycytidine, the mRNA levels of FNDC1 and SCNN1A were decreased, and no significant alteration in the mRNA levels of MAP1B, ANLN, KCNN4 and STC2 was observed. The incidence of differential methylation varied in subregions of the genome and the effects on transcription were associated with the distribution of DEGs across the genome. The regulation of this appears more complex than initially postulated. Combining the data on epigenetic differences of OA with the genome or transcriptome data for analysis may improve the understanding of the pathophysiological processes of OA. FNDC1 and SCNN1A may potentially be valuable biomarkers for OA.

The impact of epigenetic information on genome evolution

Epigenetic information affects gene function by interacting with chromatin, while not changing the DNA sequence itself. However, it has become apparent that the interactions between epigenetic information and chromatin can, in fact, indirectly lead to DNA mutations and ultimately influence genome evolution. This review evaluates the ways in which epigenetic information affects genome sequence and evolution. We discuss how DNA methylation has strong and pervasive effects on DNA sequence evolution in eukaryotic organisms. We also review how the physical interactions arising from the connections between histone proteins and DNA affect DNA mutation and repair. We then discuss how a variety of epigenetic mechanisms exert substantial effects on genome evolution by suppressing the movement of transposable elements. Finally, we examine how genome expansion through gene duplication is also partially controlled by epigenetic information. Overall, we conclude that epigenetic information has widespread indirect effects on DNA sequences in eukaryotes and represents a potent cause and constraint of genome evolution. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Alcohol consumption, depression, overweight and cortisol levels as determining factors for NR3C1 gene methylation

The NR3C1 glucocorticoid receptor (GR) gene is a component of the stress response system, which can be regulated by epigenetic mechanisms. NR3C1 methylation has been associated with trauma and mental issues, including depression, post-traumatic stress, anxiety, and personality disorders. Previous studies have reported that stressful events are involved in NR3C1 gene methylation, suggesting that its regulation under environmental effects is complex. The present study aimed to analyze associations involving stressors such as socioeconomic status, health conditions, and lifestyle in relation to NR3C1 methylation in adults. This study included 386 individual users of the Brazilian Public Unified Health System (SUS), and evaluated socioeconomic and health conditions, body mass index, cortisol levels, and lifestyle. Data were correlated with NR3C1 methylation, determined using DNA pyrosequencing. The results showed that alcohol consumption, overweight, and high cortisol levels were related to NR3C1 demethylation, while depression was related to its methylation. Habits, lifestyle, and health status may influence NR3C1 gene regulation via methylation, revealing the complexity of environmental impacts on NR3C1 methylation.

EXT1 methylation promotes proliferation and migration and predicts the clinical outcome of non-small cell lung carcinoma via WNT signalling pathway

DNA methylation is important for lung cancer prognosis. In this work, it is aimed to seek novel biomarkers with DNA methylation‐expression‐pathway pattern and explore its underlying mechanism. Prognostic DNA methylation sites and mRNAs were screened in NSCLC data set from TCGA, and further validated using the samples retrospectively collected, and EXT1 was identified as a potential target. Gene body methylation of three CpG sites (cg03276982, cg11592677, cg16286281) on EXT1 was significantly associated with clinical outcome, and the EXT1 gene expression also predicted prognosis. The expression level of EXT1 was also correlated with its DNA methylation level. This observation was further validated in a new data set consist of 170 samples. Knocking down of EXT1 resulted in decreased proliferation and migration. EXT1 targets were analysed using GSEA. It is found that the WNT signalling is the potential downstream target of EXT1. Further analyses revealed that the EXT1 targets the beta‐catenin and effect migration rate of NSCLC cell lines. The WNT signalling inhibitor, XAV‐939, effectively disrupted the migration promotion effect induced by EXT1. In summary, EXT1 methylation regulates the gene expression, effects the proliferation and migration via WNT pathway and predicted a poor prognosis for NSCLC.

Gene body methylation facilitates the transcription of CTSG via antisense lncRNA AL136018.1 in dermatomyositic myoideum

Dermatomyositis (DM) is characterized as a chronic autoimmune disorder with multiple organ involvement. Our previous study has revealed that Cathepsin G (CTSG) highly expressed in dermatomyositic in vivo is regulated by DNMT3a through DNA methylation of 5'-C-phosphate-G-3' loci at exons and introns. However, the mechanism of gene body methylation on regulating CTSG transcription remains unknown. In this study, we studied quadriceps femoris tissues of six DM patients, and observed that antisense long noncoding RNA AL136018.1 contiguous to CTSG was highly expressed in skeletal muscle tissues of DM and positively correlated with the transcription level and DNA methylation level in gene body of CTSG in vivo. Moreover, we observed that the longer transcript of AL136018.1 (AL136018.1-201) could bind to third and fourth exons and third intron of CTSG via the 3'-end. Finally, AL136018.1-201 could recruit DNMT3a towards gene body via 5'-terminal for adding DNA methylation and facilitating transcription of CTSG. Taken together, our data uncovered a novel epigenetic mechanism behind the gene body methylation for transcriptional regulation of CTSG in DM.

Identification of DNA (de)methylation-related genes and their transcriptional response to environmental challenges in an invasive model ascidian

Marine invasive species are constantly challenged by acute or recurring environmental stresses during their range expansions. DNA methylation-mediated stress memory has been proposed to effectively affect species' response and enhance their overall performance in recurring environmental challenges. In order to further test this proposal in marine invasive species, we identified genes in the DNA methylation and demethylation processes in the highly invasive model species, Ciona robusta, and subsequently investigated the expression patterns of these genes under recurring salinity stresses. After a genome-wide comprehensive survey, we found a total of six genes, including two genes of DNA methyltransferase 3a (DNMT3a1 and DNMT3a2), and one gene of DNA methyltransferase 1 (DNMT1), methyl-CpG-binding domain protein 2 (MBD2), methyl-CpG-binding domain protein 4 (MBD4) and ten-eleven-translocation protein 1 (TET1). Phylogenetic reconstruction and domain arrangement analyses showed that the deduced proteins of the identified genes were evolutionarily conserved and functionally similar with their orthologs. All genes were constitutively expressed in all four tested tissues. Interestingly, we found time-dependent and stress-specific gene expression patterns under high and low salinity stresses. Under the recurring high salinity stresses, DNMT3a1 and TET1 conformed to the definition of memory genes, while under the recurring low salinity stresses, two DNMT3a paralogues were identified as the memory genes. Altogether, our results clearly showed that the transcriptional patterns of (de)methylation-related genes were significantly influenced by environmental stresses, and the transcriptional memory of some (de)methylation-related genes should play crucial roles in DNA methylation-mediated stress memory during the process of biological invasions.

------Widespread conservation and lineage-specific diversification of genome-wide DNA methylation patterns across arthropods

Cytosine methylation is an ancient epigenetic modification yet its function and extent within genomes is highly variable across eukaryotes. In mammals, methylation controls transposable elements and regulates the promoters of genes. In insects, DNA methylation is generally restricted to a small subset of transcribed genes, with both intergenic regions and transposable elements (TEs) depleted of methylation. The evolutionary origin and the function of these methylation patterns are poorly understood. Here we characterise the evolution of DNA methylation across the arthropod phylum. While the common ancestor of the arthropods had low levels of TE methylation and did not methylate promoters, both of these functions have evolved independently in centipedes and mealybugs. In contrast, methylation of the exons of a subset of transcribed genes is ancestral and widely conserved across the phylum, but has been lost in specific lineages. A similar set of genes is methylated in all species that retained exon-enriched methylation. We show that these genes have characteristic patterns of expression correlating to broad transcription initiation sites and well-positioned nucleosomes, providing new insights into potential mechanisms driving methylation patterns over hundreds of millions of years.

MethGET: web-based bioinformatics software for correlating genome-wide DNA methylation and gene expression

Background

DNA methylation is a major epigenetic modification involved in regulating gene expression. The effects of DNA methylation on gene expression differ by genomic location and vary across kingdoms, species and environmental conditions. To identify the functional regulatory roles of DNA methylation, the correlation between DNA methylation changes and alterations in gene expression is crucial. With the advance of next-generation sequencing, genome-wide methylation and gene expression profiling have become feasible. Current bioinformatics tools for investigating such correlation are designed to the assessment of DNA methylation at CG sites. The correlation of non-CG methylation and gene expression is very limited. Some bioinformatics databases allow correlation analysis, but they are limited to specific genomes such as that of humans and do not allow user-provided data.

Results

Here, we developed a bioinformatics web tool, MethGET (Methylation and Gene Expression Teller), that is specialized to analyse the association between genome-wide DNA methylation and gene expression. MethGET is the first web tool to which users can supply their own data from any genome. It is also the tool that correlates gene expression with CG, CHG, and CHH methylation based on whole-genome bisulfite sequencing data. MethGET not only reveals the correlation within an individual sample (single-methylome) but also performs comparisons between two groups of samples (multiple-methylomes). For single-methylome analyses, MethGET provides Pearson correlations and ordinal associations to investigate the relationship between DNA methylation and gene expression. It also groups genes with different gene expression levels to view the methylation distribution at specific genomic regions. Multiple-methylome analyses include comparative analyses and heatmap representations between two groups. These functions enable the detailed investigation of the role of DNA methylation in gene regulation. Additionally, we applied MethGET to rice regeneration data and discovered that CHH methylation in the gene body region may play a role in the tissue culture process, which demonstrates the capability of MethGET for use in epigenomic research.

Conclusions

MethGET is a Python software that correlates DNA methylation and gene expression. Its web interface is publicly available at https://paoyang.ipmb.sinica.edu.tw/Software.html. The stand-alone version and source codes are available on GitHub at https://github.com/Jason-Teng/MethGET.

Population epigenetic divergence exceeds genetic divergence in the Eastern oyster Crassostrea virginica in the Northern Gulf of Mexico

Populations may respond to environmental heterogeneity via evolutionary divergence or phenotypic plasticity. While evolutionary divergence occurs through DNA sequence differences among populations, plastic divergence among populations may be generated by changes in the epigenome. Here, we present the results of a genome-wide comparison of DNA methylation patterns and genetic structure among four populations of Eastern oyster (Crassostrea virginica) in the northern Gulf of Mexico. We used a combination of restriction site-associated DNA sequencing (RADseq) and reduced representation bisulfite sequencing (RRBS) to explore population structure, gene-wide averages of F ST, and DNA methylation differences between oysters inhabiting four estuaries with unique salinity profiles. This approach identified significant population structure despite a moderately low F ST (0.02) across the freshwater boundary of the Mississippi river, a finding that may reflect recent efforts to restore oyster stock populations. Divergence between populations in CpG methylation was greater than for divergence in F ST, likely reflecting environmental effects on DNA methylation patterns. Assessment of CpG methylation patterns across all populations identified that only 26% of methylated DNA was intergenic; and, only 17% of all differentially methylated regions (DMRs) were within these same regions. DMRs within gene bodies between sites were associated with genes known to be involved in DNA damage repair, ion transport, and reproductive timing. Finally, when assessing the correlation between genomic variation and DNA methylation between these populations, we observed population-specific DNA methylation profiles that were not directly associated with single nucleotide polymorphisms or broader gene-body mean F ST trends. Our results suggest that C. virginica may use DNA methylation to generate environmentally responsive plastic phenotypes and that there is more divergence in methylation than divergence in allele frequencies.

Transgenerational acclimation to changes in ocean acidification in marine invertebrates

The rapid pace of increasing oceanic acidity poses a major threat to the fitness of the marine ecosystem, as well as the buffering capacity of the oceans. Disruption in chemical equilibrium in the ocean leads to decreased carbonate ion precipitation, resulting in calcium carbonate saturation. If these trends continue, calcifying invertebrates will experience difficultly maintaining their calcium carbonate exoskeleton and shells. Because malfunction of exoskeleton formation by calcifiers in response to ocean acidification (OA) will have non-canonical biological cascading results in the marine ecosystem, many studies have investigated the direct and indirect consequences of OA on ecosystem- and physiology-related traits of marine invertebrates. Considering that evolutionary adaptation to OA depends on the duration of OA effects, long-term exposure to OA stress over multi-generations may result in adaptive mechanisms that increase the potential fitness of marine invertebrates in response to OA. Transgenerational studies have the potential to elucidate the roles of acclimation, carryover effects, and evolutionary adaptation within and over generations in response to OA. In particular, understanding mechanisms of transgenerational responses (e.g., antioxidant responses, metabolic changes, epigenetic reprogramming) to changes in OA will enhance our understanding of marine invertebrate in response to rapid climate change.

Distinct epigenomic and transcriptomic modifications associated with Wolbachia-mediated asexuality

Wolbachia are maternally transmitted intracellular bacteria that induce a range of pathogenic and fitness-altering effects on insect and nematode hosts. In parasitoid wasps of the genus Trichogramma, Wolbachia infection induces asexual production of females, thus increasing transmission of Wolbachia. It has been hypothesized that Wolbachia infection accompanies a modification of the host epigenome. However, to date, data on genome-wide epigenomic changes associated with Wolbachia are limited, and are often confounded by background genetic differences. Here, we took sexually reproducing Trichogramma free of Wolbachia and introgressed their genome into a Wolbachia-infected cytoplasm, converting them to Wolbachia-mediated asexuality. Wolbachia was then cured from replicates of these introgressed lines, allowing us to examine the genome-wide effects of wasps newly converted to asexual reproduction while controlling for genetic background. We thus identified gene expression and DNA methylation changes associated with Wolbachia-infection. We found no overlaps between differentially expressed genes and differentially methylated genes, indicating that Wolbachia-infection associated DNA methylation change does not directly modulate levels of gene expression. Furthermore, genes affected by these mechanisms exhibit distinct evolutionary histories. Genes differentially methylated due to the infection tended to be evolutionarily conserved. In contrast, differentially expressed genes were significantly more likely to be unique to the Trichogramma lineage, suggesting host-specific transcriptomic responses to infection. Nevertheless, we identified several novel aspects of Wolbachia-associated DNA methylation changes. Differentially methylated genes included those involved in oocyte development and chromosome segregation. Interestingly, Wolbachia-infection was associated with higher levels of DNA methylation. Additionally, Wolbachia infection reduced overall variability in gene expression, even after accounting for the effect of DNA methylation. We also identified specific cases where alternative exon usage was associated with DNA methylation changes due to Wolbachia infection. These results begin to reveal distinct genes and molecular pathways subject to Wolbachia induced epigenetic modification and/or host responses to Wolbachia-infection.

Wolbachia is an extremely common endosymbiotic infection of arthropods and nematodes. One of the reasons why Wolbachia can so successfully infect diverse species is the bacterium’s ability to profoundly alter the reproductive behavior of its host. It has been proposed that Wolbachia may modify host’s epigenetic programs to alter its reproductive behavior. However, it has been difficult to study how epigenetic programs change with Wolbachia infection, due to the confounding effects of genetic backgrounds. Here, we studied host transcriptome and epigenome changes associated with Wolbachia infection in a homogenous genetic background, by carrying out an innovative introgression scheme. By doing so, we show, for the first time, high-resolution molecular consequences of intracellular infection and offer insights into epigenetic and transcriptomic regulation of invertebrates.

Genome-wide analysis of the role of DNA methylation in inbreeding depression of reproduction in Langshan chicken

Inbreeding depression of chicken reproduction is a major concern in the conservation of chicken genetic resources. To investigate the potential DNA methylation sites involved in the inbreeding depression of chicken reproduction, we carried out whole-genome bisulfite sequencing (WGBS) of hypothalamus and ovary tissues from the strongly and weakly inbred Langshan chickens, respectively. 5948 and 4593 differentially methylated regions (DMRs) were identified in the hypothalamus and ovary between the strongly and weakly inbred Langshan chickens, respectively. Large numbers of DMR-related genes (DMGs) were enriched in reproduction-related pathways. By combining the WGBS and transcriptome data, two DMRs in SRD5A1 and CDC27 genes were inferred as the most likely biomarkers of inbreeding depression of reproduction in Langshan chicken. Our study provides the first systematic investigation of the DNA methylation changes in strongly inbred chickens, and extends our understanding of the regulatory mechanisms underlying inbreeding depression in chicken reproduction.