Conservation and divergence of methylation patterning in plants and animals

Epigenomic profiling indicates a role for DNA methylation in the postnatal liver and pancreas development of giant pandas

Giant pandas are unique within Carnivora with a strict bamboo diet. Here, the epigenomic profiles of giant panda liver and pancreas tissues collected from three important feeding stages were investigated using BS-seq. Few differences in DNA methylation profiles were exhibited between no feeding and suckling groups in both tissues. However, we observed a tendency toward a global loss of DNA methylation in the gene-body and promoter region of metabolism-related genes from newborn to adult. Correlation analysis revealed a significant negative correlation between the changes in methylation levels within gene promoters and gene expression. The majority of genes related to nutrition metabolism had lost DNA methylation with increased mRNA expression in adult giant pandas. The few galactose metabolism and unsaturated fatty acid metabolism related genes that were hypomethylated and highly-expressed at early stages of giant panda development may meet the nutritional requirement of this species' highly altricial neonates.

Deep6mAPred: A CNN and Bi-LSTM-based deep learning method for predicting DNA N6-methyladenosine sites across plant species

DNA N6-methyladenine (6mA) is a key DNA modification, which plays versatile roles in the cellular processes, including regulation of gene expression, DNA repair, and DNA replication. DNA 6mA is closely associated with many diseases in the mammals and with growth as well as development of plants. Precisely detecting DNA 6mA sites is of great importance to exploration of 6mA functions. Although many computational methods have been presented for DNA 6mA prediction, there is still a wide gap in the practical application. We presented a convolution neural network (CNN) and bi-directional long-short term memory (Bi-LSTM)-based deep learning method (Deep6mAPred) for predicting DNA 6mA sites across plant species. The Deep6mAPred stacked the CNNs and the Bi-LSTMs in a paralleling manner instead of a series-connection manner. The Deep6mAPred also employed the attention mechanism for improving the representations of sequences. The Deep6mAPred reached an accuracy of 0.9556 over the independent rice dataset, far outperforming the state-of-the-art methods. The tests across plant species showed that the Deep6mAPred is of a remarkable advantage over the state of the art methods. We developed a user-friendly web application for DNA 6mA prediction, which is freely available at http://106.13.196.152:7001/ for all the scientific researchers. The Deep6mAPred would enrich tools to predict DNA 6mA sites and speed up the exploration of DNA modification.

CYP81A68 confers metabolic resistance to ALS and ACCase-inhibiting herbicides and its epigenetic regulation in Echinochloa crus-galli

Long-term and excessive herbicide use has led to some environmental concerns and especially, herbicide resistance evolution in weeds. Here, we confirmed acetolactate synthase (ALS) inhibiting herbicide penoxsulam resistance and cross resistance to acetyl-coenzyme carboxylase (ACCase) inhibiting herbicides (cyhalofop-butyl and metamifop) in a global weed Echinochloa crus-galli population resistant to these herbicides (R). Penoxsulam metabolism study indicated that degradation rate was significantly higher in R than susceptible E. crus-galli population (S). RNA-sequencing revealed that a cytochrome P450 (P450) gene, CYP81A68, expressed higher in R versus S. Rice seedlings overexpressing this CYP81A68 gene are resistant to penoxsulam, cyhalofop-butyl and metamifop, and penoxsulam resistance is due to enhanced metabolism via O-demethylation. Deletion analysis of the CYP81A68 gene promoter identified an efficient region, in which differential methylation of CpG islands occurred between R and S. Collectively, these results demonstrate that upregulation of E. crus-galli CYP81A68 gene endows generalist metabolic resistance to commonly used ALS- and ACCase-inhibiting herbicides in rice fields and epigenetic regulation may play a role in the resistance evolution. This research could contribute to strategies reducing herbicide environmental impacts by judicious selection of alternative herbicide and non-chemical control tactics.

Emerging Cnidarian Models for the Study of Epithelial Polarity

Epithelial tissues are vital to the function of most organs, providing critical functions such as secretion, protection, and absorption. Cells within an epithelial layer must coordinate to create functionally distinct apical, lateral, and basal surfaces in order to maintain proper organ function and organism viability. This is accomplished through the careful targeting of polarity factors to their respective locations within the cell, as well as the strategic placement of post-mitotic cells within the epithelium during tissue morphogenesis. The process of establishing and maintaining epithelial tissue integrity is conserved across many species, as important polarity factors and spindle orientation mechanisms can be found in many phyla. However, most of the information gathered about these processes and players has been investigated in bilaterian organisms such as C. elegans, Drosophila, and vertebrate species. This review discusses the advances made in the field of epithelial polarity establishment from more basal organisms, and the advantages to utilizing these simpler models. An increasing number of cnidarian model organisms have been sequenced in recent years, such as Hydra vulgaris and Nematostella vectensis. It is now feasible to investigate how polarity is established and maintained in basal organisms to gain an understanding of the most basal requirements for epithelial tissue morphogenesis.

A comparative methylome analysis reveals conservation and divergence of DNA methylation patterns and functions in vertebrates

Background

Cytosine DNA methylation is a heritable epigenetic mark present in most eukaryotic groups. While the patterns and functions of DNA methylation have been extensively studied in mouse and human, their conservation in other vertebrates remains poorly explored. In this study, we interrogated the distribution and function of DNA methylation in primary fibroblasts of seven vertebrate species including bio-medical models and livestock species (human, mouse, rabbit, dog, cow, pig, and chicken).

Results

Our data highlight both divergence and conservation of DNA methylation patterns and functions. We show that the chicken genome is hypomethylated compared to other vertebrates. Furthermore, compared to mouse, other species show a higher frequency of methylation of CpG-rich DNA. We reveal the conservation of large unmethylated valleys and patterns of DNA methylation associated with X-chromosome inactivation through vertebrate evolution and make predictions of conserved sets of imprinted genes across mammals. Finally, using chemical inhibition of DNA methylation, we show that the silencing of germline genes and endogenous retroviruses (ERVs) are conserved functions of DNA methylation in vertebrates.

Conclusions

Our data highlight conserved properties of DNA methylation in vertebrate genomes but at the same time point to differences between mouse and other vertebrate species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01270-x.

Gene-body methylation in plants: mechanisms, functions and important implications for understanding evolutionary processes

Gene body methylation (gbM) is an epigenetic mark where gene exons are methylated in the CG context only, as opposed to CHG and CHH contexts (where H stands for A, C, or T). CG methylation is transmitted transgenerationally in plants, opening the possibility that gbM may be shaped by adaptation. This presupposes, however, that gbM has a function that affects phenotype, which has been a topic of debate in the literature. Here, we review our current knowledge of gbM in plants. We start by presenting the well-elucidated mechanisms of plant gbM establishment and maintenance. We then review more controversial topics: the evolution of gbM and the potential selective pressures that act on it. Finally, we discuss the potential functions of gbM that may affect organismal phenotypes: gene expression stabilization and upregulation, inhibition of aberrant transcription (reverse and internal), prevention of aberrant intron retention, and protection against TE insertions. To bolster the review of these topics, we include novel analyses to assess the effect of gbM on transcripts. Overall, a growing body of literature finds that gbM correlates with levels and patterns of gene expression. It is not clear, however, if this is a causal relationship. Altogether, functional work suggests that the effects of gbM, if any, must be relatively small, but there is nonetheless evidence that it is shaped by natural selection. We conclude by discussing the potential adaptive character of gbM and its implications for an updated view of the mechanisms of adaptation in plants.

Methylation of recombinant mononucleosomes by DNMT3A demonstrates efficient linker DNA methylation and a role of H3K36me3

Recently, the structure of the DNMT3A2/3B3 heterotetramer complex bound to a mononucleosome was reported. Here, we investigate DNA methylation of recombinant unmodified, H3K_C4me3 and H3K_C36me3 containing mononucleosomes by DNMT3A2, DNMT3A catalytic domain (DNMT3AC) and the DNMT3AC/3B3C complex. We show strong protection of the nucleosomal bound DNA against methylation, but efficient linker-DNA methylation next to the nucleosome core. High and low methylation levels of two specific CpG sites next to the nucleosome core agree well with details of the DNMT3A2/3B3-nucleosome structure. Linker DNA methylation next to the nucleosome is increased in the absence of H3K4me3, likely caused by binding of the H3-tail to the ADD domain leading to relief of autoinhibition. Our data demonstrate a strong stimulatory effect of H3K36me3 on linker DNA methylation, which is independent of the DNMT3A-PWWP domain. This observation reveals a direct functional role of H3K36me3 on the stimulation of DNA methylation, which could be explained by hindering the interaction of the H3-tail and the linker DNA. We propose an evolutionary model in which the direct stimulatory effect of H3K36me3 on DNA methylation preceded its signaling function, which could explain the evolutionary origin of the widely distributed "active gene body-H3K36me3-DNA methylation" connection.

Epigenetic regulation of seed-specific gene expression by DNA methylation valleys in castor bean

BACKGROUND

Understanding the processes governing angiosperm seed growth and development is essential both for fundamental plant biology and for agronomic purposes. Master regulators of angiosperm seed development are expressed in a seed-specific manner. However, it is unclear how this seed specificity of transcription is established. In some vertebrates, DNA methylation valleys (DMVs) are highly conserved and strongly associated with key developmental genes, but comparable studies in plants are limited to Arabidopsis and soybean. Castor bean (Ricinus communis) is a valuable model system for the study of seed biology in dicots and source of economically important castor oil. Unlike other dicots such as Arabidopsis and soybean, castor bean seeds have a relatively large and persistent endosperm throughout seed development, representing substantial structural differences in mature seeds. Here, we performed an integrated analysis of RNA-seq, whole-genome bisulfite sequencing, and ChIP-seq for various histone marks in the castor bean.

RESULTS

We present a gene expression atlas covering 16 representative tissues and identified 1162 seed-specific genes in castor bean (Ricinus communis), a valuable model for the study of seed biology in dicots. Upon whole-genome DNA methylation analyses, we detected 32,567 DMVs across five tissues, covering ~33% of the castor bean genome. These DMVs are highly hypomethylated during development and conserved across plant species. We found that DMVs have the potential to activate transcription, especially that of tissue-specific genes. Focusing on seed development, we found that many key developmental regulators of seed/endosperm development, including AGL61, AGL62, LEC1, LEC2, ABI3, and WRI1, were located within DMVs. ChIP-seq for five histone modifications in leaves and seeds clearly showed that the vast majority of histone modification peaks were enriched within DMVs, and their remodeling within DMVs has a critical role in the regulation of seed-specific gene expression. Importantly, further experiment analysis revealed that distal DMVs may act as cis-regulatory elements, like enhancers, to activate downstream gene expression.

CONCLUSIONS

Our results point to the importance of DMVs and special distal DMVs behaving like enhancers, in the regulation of seed-specific genes, via the reprogramming of histone modifications within DMVs. Furthermore, these results provide a comprehensive understanding of the epigenetic regulator roles in seed development in castor bean and other important crops.

Distribution and regulatory roles of oxidized 5-methylcytosines in DNA and RNA of the basidiomycete fungi Laccaria bicolor and Coprinopsis cinerea

The formation of three oxidative DNA 5-methylcytosine (5mC) modifications (oxi-mCs)-5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC)-by the TET/JBP family of dioxygenases prompted intensive studies of their functional roles in mammalian cells. However, the functional interplay of these less abundant modified nucleotides in other eukaryotic lineages remains poorly understood. We carried out a systematic study of the content and distribution of oxi-mCs in the DNA and RNA of the basidiomycetes Laccaria bicolor and Coprinopsis cinerea, which are established models to study DNA methylation and developmental and symbiotic processes. Quantitative liquid chromatography-tandem mass spectrometry revealed persistent but uneven occurrences of 5hmC, 5fC and 5caC in the DNA and RNA of the two organisms, which could be upregulated by vitamin C. 5caC in RNA (5carC) was predominantly found in non-ribosomal RNA, which potentially includes non-coding, messenger and small RNA species. Genome-wide mapping of 5hmC and 5fC using the single CG analysis techniques hmTOP-seq and foTOP-seq pointed at involvement of oxi-mCs in the regulation of gene expression and silencing of transposable elements. The implicated diverse roles of 5mC and oxi-mCs in the two fungi highlight the epigenetic importance of the latter modifications, which are often neglected in standard whole-genome bisulfite analyses.

Stage-differentiated ensemble modeling of DNA methylation landscapes uncovers salient biomarkers and prognostic signatures in colorectal cancer progression

Background

Aberrant DNA methylation acts epigenetically to skew the gene transcription rate up or down, contributing to cancer etiology. A gap in our understanding concerns the epigenomics of stagewise cancer progression. In this study, we have developed a comprehensive computational framework for the stage-differentiated modelling of DNA methylation landscapes in colorectal cancer (CRC).

Methods

The methylation β-matrix was derived from the public-domain TCGA data, converted into M-value matrix, annotated with AJCC stages, and analysed for stage-salient genes using an ensemble of approaches involving stage-differentiated modelling of methylation patterns and/or expression patterns. Differentially methylated genes (DMGs) were identified using a contrast against controls (adjusted p-value <0.001 and |log fold-change of M-value| >2), and then filtered using a series of all possible pairwise stage contrasts (p-value <0.05) to obtain stage-salient DMGs. These were then subjected to a consensus analysis, followed by matching with clinical data and performing Kaplan–Meier survival analysis to evaluate the impact of methylation patterns of consensus stage-salient biomarkers on disease prognosis.

Results

We found significant genome-wide changes in methylation patterns in cancer cases relative to controls agnostic of stage. The stage-differentiated models yielded the following consensus salient genes: one stage-I gene (FBN1), one stage-II gene (FOXG1), one stage-III gene (HCN1) and four stage-IV genes (NELL1, ZNF135, FAM123A, LAMA1). All the biomarkers were significantly hypermethylated in the promoter regions, indicating down-regulation of expression and implying a putative CpG island Methylator Phenotype (CIMP) manifestation. A prognostic signature consisting of FBN1 and FOXG1 survived all the analytical filters, and represents a novel early-stage epigenetic biomarker / target.

Conclusions

We have designed and executed a workflow for stage-differentiated epigenomic analysis of colorectal cancer progression, and identified several stage-salient diagnostic biomarkers, and an early-stage prognostic biomarker panel. The study has led to the discovery of an alternative CIMP-like signature in colorectal cancer, reinforcing the role of CIMP drivers in tumor pathophysiology.

Exploring Breakthroughs in Three Traits Belonging to Seed Life

Based on prior knowledge and with the support of new methodology, solid progress in the understanding of seed life has taken place over the few last years. This update reflects recent advances in three key traits of seed life (i.e., preharvest sprouting, genomic imprinting, and stored-mRNA). The first breakthrough refers to cloning of the mitogen-activated protein kinase-kinase 3 (MKK3) gene in barley and wheat. MKK3, in cooperation with ABA signaling, controls seed dormancy. This advance has been determinant in producing improved varieties that are resistant to preharvest sprouting. The second advance concerns to uniparental gene expression (i.e., imprinting). Genomic imprinting primarily occurs in the endosperm. Although great advances have taken place in the last decade, there is still a long way to go to complete the puzzle regarding the role of genomic imprinting in seed development. This trait is probably one of the most important epigenetic facets of developing endosperm. An example of imprinting regulation is polycomb repressive complex 2 (PRC2). The mechanism of PRC2 recruitment to target endosperm with specific genes is, at present, robustly studied. Further progress in the knowledge of recruitment of PRC2 epigenetic machinery is considered in this review. The third breakthrough referred to in this update involves stored mRNA. The role of the population of this mRNA in germination is far from known. Its relations to seed aging, processing bodies (P bodies), and RNA binding proteins (RBPs), and how the stored mRNA is targeted to monosomes, are aspects considered here. Perhaps this third trait is the one that will require greater experimental dedication in the future. In order to make progress, herein are included some questions that are needed to be answered.

Micronutrient supplementation affects DNA methylation in male gonads with potential intergenerational epigenetic inheritance involving the embryonic development through glutamate receptor-associated genes

BACKGROUND

DNA methylation has an important role in intergenerational inheritance. An increasing number of studies have reported evidence of germline inheritance of DNA methylation induced by nutritional signals in mammals. Vitamins and minerals as micronutrients contribute to growth performance in vertebrates, including Atlantic salmon (Salmo salar), and also have a role in epigenetics as environmental factors that alter DNA methylation status. It is important to understand whether micronutrients in the paternal diet can influence the offspring through alterations of DNA methylation signatures in male germ cells.

RESULTS

Here, we show the effect of micronutrient supplementation on DNA methylation profiles in the male gonad through a whole life cycle feeding trial of Atlantic salmon fed three graded levels of micronutrient components. Our results strongly indicate that micronutrient supplementation affects the DNA methylation status of genes associated with cell signalling, synaptic signalling, and embryonic development. In particular, it substantially affects DNA methylation status in the promoter region of a glutamate receptor gene, glutamate receptor ionotropic, NMDA 3A-like (grin3a-like), when the fish are fed both medium and high doses of micronutrients. Furthermore, two transcription factors, histone deacetylase 2 (hdac2) and a zinc finger protein, bind to the hyper-methylated site in the grin3a-like promoter. An estimated function of hdac2 together with a zinc finger indicates that grin3a-like has a potential role in intergenerational epigenetic inheritance and the regulation of embryonic development affected by paternal diet.

CONCLUSIONS

The present study demonstrates alterations of gene expression patterns and DNA methylation signatures in the male gonad when Atlantic salmon are fed different levels of micronutrients. Alterations of gene expression patterns are of great interest because the gonads are supposed to have limited metabolic activities compared to other organs, whereas alterations of DNA methylation signatures are of great importance in the field of nutritional epigenetics because the signatures affected by nutrition could be transferred to the next generation. We provide extensive data resources for future work in the context of potential intergenerational inheritance through the male germline.

Characterisation of the epigenetic architecture of the major malaria vector Anopheles arabiensis (Diptera: Culicidae) after treatment with epigenetic modulators and heavy metals

Anopheles arabiensis (a member of the An. gambiae species complex) is a major vector of malaria in sub-Saharan Africa. Despite its disease vector status, there is currently a paucity of epigenetic information for this species. The aim this study was therefore to analyse global epigenetic markers and their response to metal exposure in insecticide susceptible and resistant laboratory strains of An. arabiensis. This was done using commercially available epigenetic marker quantification kits. In order to validate the efficacy of the kits, several kits were assessed to determine whether changes induced by known epigenetic modulators were detectable using these platforms. The efficacy of the dosages used were determined by examining the effect of the dosages used on insecticide resistant phenotypes. Upon confirmation that the dosages used were sufficient to induce a phenotypic change, the effect on epigenetic markers was assessed. Commercial kits were used to quantify 5-methylcysteine (5-mC) and 5-hydroxymethylcysteine (5-hmC) methylation in DNA, m6A methylation in mRNA as well as Histone Acetyl Transferase (HAT) activity. There was a marked difference in the phenotypic response in adult mosquitoes of the insecticide susceptible strain compared to that of its' resistant counterpart. For males and females of the resistant strain, exposure to nucleic acid modifying drugs typically increased their tolerance to insecticides. The patterns of changes in 5-mC methylation by epigenetic modulators was congruent with previous studies which quantified by mass spectrometry. The two strains differed in methylation patterns under control conditions and responded differentially to larval metal exposure. In the resistant strain, which previously was demonstrated to show increased detoxification enzyme activity and insecticide tolerance after the same treatment, the potential increase in transcriptional activity appeared to be modulated by reduced methylation and increased HAT activity. This study suggests that the commercial epigenetic quantification kits can be used to characterise phenotypic changes in An. arabiensis, and also shows that epigenetic regulation of the response to metal exposure is regulated at the DNA as opposed to the RNA level.

Studying phenotypic variation and DNA methylation across development, ecology and evolution in the clonal marbled crayfish: a paradigm for investigating epigenotype-phenotype relationships in macro-invertebrates

Animals can produce different phenotypes from the same genome during development, environmental adaptation and evolution, which is mediated by epigenetic mechanisms including DNA methylation. The obligatory parthenogenetic marbled crayfish, Procambarus virginalis, whose genome and methylome are fully established, proved very suitable to study this issue in detail. Comparison between developmental stages and DNA methylation revealed low expression of Dnmt methylation and Tet demethylation enzymes from the spawned oocyte to the 256 cell embryo and considerably increased expression thereafter. The global 5-methylcytosine level was 2.78% at mid-embryonic development and decreased slightly to 2.41% in 2-year-old adults. Genetically identical clutch-mates raised in the same uniform laboratory setting showed broad variation in morphological, behavioural and life history traits and differences in DNA methylation. The invasion of diverse habitats in tropical to cold-temperate biomes in the last 20 years by the marbled crayfish was associated with the expression of significantly different phenotypic traits and DNA methylation patterns, despite extremely low genetic variation on the whole genome scale, suggesting the establishment of epigenetic ecotypes. The evolution of marbled crayfish from its parent species Procambarus fallax by autotriploidy a few decades ago was accompanied by a significant increase in body size, fertility and life span, a 20% reduction of global DNA methylation and alteration of methylation in hundreds of genes, suggesting that epigenetic mechanisms were involved in speciation and fitness enhancement. The combined analysis of phenotypic traits and DNA methylation across multiple biological contexts in the laboratory and field in marbled crayfish may serve as a blueprint for uncovering the role of epigenetic mechanisms in shaping of phenotypes in macro-invertebrates.

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.

Dissecting cricket genomes for the advancement of entomology and entomophagy

Significant advances in biophysical methods such as next-generation sequencing technologies have now opened the way to conduct evolutionary and applied research based on the genomic information of greatly diverse insects. Crickets belonging to Orthoptera (Insecta: Polyneoptera), one of the most flourishing groups of insects, have contributed to the development of multiple scientific fields including developmental biology and neuroscience and have been attractive targets in evolutionary ecology for their diverse ecological niches. In addition, crickets have recently gained recognition as food and feed. However, the genomic information underlying their biological basis and application research toward breeding is currently underrepresented. In this review, we summarize the progress of genomics of crickets. First, we outline the phylogenetic position of crickets in insects and then introduce recent studies on cricket genomics and transcriptomics in a variety of fields. Furthermore, we present findings from our analysis of polyneopteran genomes, with a particular focus on their large genome sizes, chromosome number, and repetitive sequences. Finally, how the cricket genome can be beneficial to the food industry is discussed. This review is expected to enhance greater recognition of how important the cricket genomes are to the multiple biological fields and how basic research based on cricket genome information can contribute to tackling global food security.

Quality control and evaluation of plant epigenomics data

Epigenomics is the study of molecular signatures associated with discrete regions within genomes, many of which are important for a wide range of nuclear processes. The ability to profile the epigenomic landscape associated with genes, repetitive regions, transposons, transcription, differential expression, cis-regulatory elements, and 3D chromatin interactions has vastly improved our understanding of plant genomes. However, many epigenomic and single-cell genomic assays are challenging to perform in plants, leading to a wide range of data quality issues; thus, the data require rigorous evaluation prior to downstream analyses and interpretation. In this commentary, we provide considerations for the evaluation of plant epigenomics and single-cell genomics data quality with the aim of improving the quality and utility of studies using those data across diverse plant species.

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

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

Distinctive aspects of the placental epigenome and theories as to how they arise

The placenta has a methylome dramatically unlike that of any somatic cell type. Among other distinctions, it features low global DNA methylation, extensive “partially methylated domains” packed in dense heterochromatin and methylation of hundreds of CpG islands important in somatic development. These features attract interest in part because a substantial fraction of human cancers feature the exact same phenomena, suggesting parallels between epigenome formation in placentation and cancer. Placenta also features an expanded set of imprinted genes, some of which come about by distinctive developmental pathways. Recent discoveries, some from far outside the placental field, shed new light on how the unusual placental epigenetic state may arise. Nonetheless, key questions remain unresolved.

DNA Methylome and LncRNAome Analysis Provide Insights Into Mechanisms of Genome-Dosage Effects in Autotetraploid Cassava

Whole genome duplication (WGD) increases the dosage of all coding and non-coding genes, yet the molecular implications of genome-dosage effects remain elusive. In this study, we generated integrated maps of the methylomes and lncRNAomes for diploid and artificially generated autotetraploid cassava (Manihot esculenta Crantz). We found that transposable elements (TEs) suppressed adjacent protein coding gene (PCG)-expression levels, while TEs activated the expression of nearby long non-coding RNAs (lncRNAs) in the cassava genome. The hypermethylation of DNA transposons in mCG and mCHH sites may be an effective way to suppress the expression of nearby PCGs in autotetraploid cassava, resulting in similar expression levels for most of PCGs between autotetraploid and diploid cassava. In the autotetraploid, decreased methylation levels of retrotransposons at mCHG and mCHH sites contributed to reduced methylation of Gypsy-neighboring long intergenic non-coding RNAs, potentially preserving diploid-like expression patterns in the major of lncRNAs. Collectively, our study highlighted that WGD-induced DNA methylation variation in DNA transposons and retrotransposons may be as direct adaptive responses to dosage of all coding-genes and lncRNAs, respectively.

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

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

Nanopore Sequencing and Data Analysis for Base-Resolution Genome-Wide 5-Methylcytosine Profiling

Whole-genome bisulfite sequencing (WGBS) is currently the gold standard for DNA methylation (5-methylcytosine, 5mC) profiling; however, the destructive nature of sodium bisulfite results in DNA fragmentation and subsequent biases in sequencing data. Such issues have led to the development of bisulfite-free methods for 5mC detection. Nanopore sequencing is a long read nondestructive approach that directly analyzes DNA and RNA fragments in real time. Recently, computational tools have been developed that enable base-resolution detection of 5mC from Oxford Nanopore sequencing data. In this chapter, we provide a detailed protocol for preparation, sequencing, read assembly, and analysis of genome-wide 5mC using Nanopore sequencing technologies.

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.

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

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

Histone H1 prevents non-CG methylation-mediated small RNA biogenesis in Arabidopsis heterochromatin

Flowering plants utilize small RNA (sRNA) molecules to guide DNA methyltransferases to genomic sequences. This RNA-directed DNA methylation (RdDM) pathway preferentially targets euchromatic transposable elements. However, RdDM is thought to be recruited by methylation of histone H3 at lysine 9 (H3K9me), a hallmark of heterochromatin. How RdDM is targeted to euchromatin despite an affinity for H3K9me is unclear. Here, we show that loss of histone H1 enhances heterochromatic RdDM, preferentially at nucleosome linker DNA. Surprisingly, this does not require SHH1, the RdDM component that binds H3K9me. Furthermore, H3K9me is dispensable for RdDM, as is CG DNA methylation. Instead, we find that non-CG methylation is specifically associated with sRNA biogenesis, and without H1 sRNA production quantitatively expands to non-CG-methylated loci. Our results demonstrate that H1 enforces the separation of euchromatic and heterochromatic DNA methylation pathways by excluding the sRNA-generating branch of RdDM from non-CG-methylated heterochromatin.

Cells adapt to different roles by turning different groups of genes on and off. One way cells control which genes are on or off is by creating regions of active and inactive DNA, which are created and maintained by different groups of proteins. Genes in active DNA regions can be turned on, while genes in inactive regions are switched off or silenced. Silenced DNA regions also turn off ‘transposable elements’: pieces of DNA that can copy themselves and move to other regions of the genome if they become active. Transposons can be dangerous if they are activated, because they can disrupt genes or regulatory sequences when they move.

There are different types of active and inactive DNA, but it is not always clear why these differences exist, or how they are maintained over time. In plants, such as the commonly-studied weed Arabidopsis thaliana, there are two types of inactive DNA, called E and H, that can silence transposons. In both types, DNA has small chemicals called methyl groups attached to it, which help inactivate the DNA. Type E DNA is methylated by a process called RNA-directed DNA methylation (RdDM), but RdDM is rarely seen in type H DNA.

Choi, Lyons and Zilberman showed that RdDM is attracted to E and H regions by previously existing methylated DNA. However, in the H regions, a protein called histone H1 blocks RdDM from attaching methyl groups. This helps focus RdDM onto E regions where it is most needed, because E regions contain the types of transposons RdDM is best suited to silence.

When Choi, Lyons and Zilberman examined genetically modified A. thaliana plants that do not produce histone H1, they found that RdDM happened in both E and H regions. There are many more H regions than E regions, so stretching RdDM across both made it less effective at silencing DNA.

This work shows how different DNA silencing processes are focused onto specific genetic regions, helping explain why there are different types of active and inactive DNA within cells. RdDM has been studied as a way to affect crop growth and yield by altering DNA methylation. These results may help such studies by explaining how RdDM is naturally targeted.

EpiDiverse Toolkit: a pipeline suite for the analysis of bisulfite sequencing data in ecological plant epigenetics

The expanding scope and scale of next generation sequencing experiments in ecological plant epigenetics brings new challenges for computational analysis. Existing tools built for model data may not address the needs of users looking to apply these techniques to non-model species, particularly on a population or community level. Here we present a toolkit suitable for plant ecologists working with whole genome bisulfite sequencing; it includes pipelines for mapping, the calling of methylation values and differential methylation between groups, epigenome-wide association studies, and a novel implementation for both variant calling and discriminating between genetic and epigenetic variation.

Dynamic Changes of DNA Methylation During Wild Strawberry (Fragaria nilgerrensis) Tissue Culture

Tissue culture is an important tool for asexual propagation and genetic transformation of strawberry plants. In plant tissue culture, variation of DNA methylation is a potential source of phenotypic variation in regenerated plants. However, the genome wide dynamic methylation patterns of strawberry tissue culture remain unclear. In this study, we used whole-genome bisulfite sequencing (WGBS) to study genomic DNA methylation changes of a wild strawberry Fragaria nilgerrensis at six stages: from explants of shoot tips to outplanting and acclimation. Global methylation levels showed that CG sites exhibited the highest methylation level in all stages with an average of 49.5%, followed by CHG (33.2%) and CHH (12.4%). Although CHH accounted for the lowest proportion of total cytosine methylation, it showed the most obvious methylation change and the most of these changes occurred in the transposable element regions. The overall methylation levels alternately decreased and increased during the entire tissue culture process and the distribution of DNA methylation was non-uniform among different genetic regions. Furthermore, much more differentially methylated regions (DMRs) were detected in dedifferentiation and redifferentiation stages and most of them were transposable elements, suggesting these processes involved activating or silencing of amounts of transposons. The functional enrichment of the DMR-related genes indicated that genes involved in hormone metabolic processes, plant development and the stress response changed methylation throughout the tissue culture process. Finally, the quantitative real-time PCR (qRT-PCR) was conducted to examine the association of methylation and gene expression of a set of different methylated genes. Our findings give deeper insight into the epigenetic regulation of gene expression during the plant tissue cultures process, which will be useful in the efficient control of somaclonal variations and in crop improvement.

Extraction and high-throughput sequencing of oak heartwood DNA: Assessing the feasibility of genome-wide DNA methylation profiling

Tree ring features are affected by environmental factors and therefore are the basis for dendrochronological studies to reconstruct past environmental conditions. Oak wood often provides the data for these studies because of the durability of oak heartwood and hence the availability of samples spanning long time periods of the distant past. Wood formation is regulated in part by epigenetic mechanisms such as DNA methylation. Studies of the methylation state of DNA preserved in oak heartwood thus could identify epigenetic tree ring features informing on past environmental conditions. In this study, we aimed to establish protocols for the extraction of DNA, the high-throughput sequencing of whole-genome DNA libraries (WGS) and the profiling of DNA methylation by whole-genome bisulfite sequencing (WGBS) for oak (Quercus robur) heartwood drill cores taken from the trunks of living standing trees spanning the AD 1776-2014 time period. Heartwood contains little DNA, and large amounts of phenolic compounds known to hinder the preparation of high-throughput sequencing libraries. Whole-genome and DNA methylome library preparation and sequencing consistently failed for oak heartwood samples more than 100 and 50 years of age, respectively. DNA fragmentation increased with sample age and was exacerbated by the additional bisulfite treatment step during methylome library preparation. Relative coverage of the non-repetitive portion of the oak genome was sparse. These results suggest that quantitative methylome studies of oak hardwood will likely be limited to relatively recent samples and will require a high sequencing depth to achieve sufficient genome coverage.

Exposure to extracellular vesicles from Pseudomonas aeruginosa result in loss of DNA methylation at enhancer and DNase hypersensitive site regions in lung macrophages

Various pathogens use differing strategies to evade host immune response including modulating the host's epigenome. Here, we investigate if EVs secreted from P. aeruginosa alter DNA methylation in human lung macrophages, thereby potentially contributing to a dysfunctional innate immune response. Using a genome-wide DNA methylation approach, we demonstrate that P. aeruginosa EVs alter certain host cell DNA methylation patterns. We identified 1,185 differentially methylated CpGs (FDR < 0.05), which were significantly enriched for distal DNA regulatory elements including enhancer regions and DNase hypersensitive sites. Notably, all but one of the 1,185 differentially methylated CpGs were hypomethylated in association with EV exposure. Significantly hypomethylated CpGs tracked to genes including AXL, CFB and CCL23. Gene expression analysis identified 310 genes exhibiting significantly altered expression 48 hours post P. aeruginosa EV treatment, with 75 different genes upregulated and 235 genes downregulated. Some CpGs associated with cytokines such as CSF3 displayed strong negative correlations between DNA methylation and gene expression. Our infection model illustrates how secreted products (EVs) from bacteria can alter DNA methylation of the host epigenome. Changes in DNA methylation in distal DNA regulatory regions in turn can modulate cellular gene expression and potential downstream cellular processes.

Silencing and anti-silencing mechanisms that shape the epigenome in plants

Epigenome information mediates genome function and maintenance by regulating gene expression and chromatin organization. Because the epigenome pattern can change in response to internal and external environments, it may underlie an adaptive genome response that modulates phenotypes during development and in changing environments. Here I summarize recent progress in our understanding of how epigenome patterns are shaped and modulated by concerted actions of silencing and anti-silencing factors mainly in Arabidopsis thaliana. I discuss the dynamic nature of epigenome regulation, which is realized by cooperation and counteraction among silencing and anti-silencing factors, and how the dynamic epigenome mediates robust and plastic responses of plants to fluctuating environments.

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.

DNA methylation and sexual dimorphism: New insights from mealybugs

DNA methylation is an ancient epigenetic pathway found across eukaryotes. Nevertheless, the targets of DNA methylation within genomes evolve extremely rapidly. Arthropods display many such examples. The mealybug Planococcus citri has evolved methylation at promoter sequences, associated with gene silencing just as in mammals. In this issue of Molecular Ecology, Bain et al. (2021), thoroughly characterise mealybug methylation, exploring its potential functions in gene expression and the spectacular sexual dimorphism that is a characteristic of this species. Their results provide new insights into the complex relationship between DNA methylation and gene expression and highlight how rapidly different methylation systems can evolve.

Symbiont shuffling induces differential DNA methylation responses to thermal stress in the coral Montastraea cavernosa

Algal symbiont shuffling in favour of more thermotolerant species has been shown to enhance coral resistance to heat-stress. Yet, the mechanistic underpinnings and long-term implications of these changes are poorly understood. This work studied the modifications in coral DNA methylation, an epigenetic mechanism involved in coral acclimatization, in response to symbiont manipulation and subsequent heat stress exposure. Symbiont composition was manipulated in the great star coral Montastraea cavernosa through controlled thermal bleaching and recovery, producing paired ramets of three genets dominated by either their native symbionts (genus Cladocopium) or the thermotolerant species (Durusdinium trenchi). Single-base genome-wide analyses showed significant modifications in DNA methylation concentrated in intergenic regions, introns and transposable elements. Remarkably, DNA methylation changes in response to heat stress were dependent on the dominant symbiont, with twice as many differentially methylated regions found in heat-stressed corals hosting different symbionts (Cladocopium vs. D. trenchii) compared to all other comparisons. Interestingly, while differential gene body methylation was not correlated with gene expression, an enrichment in differentially methylated regions was evident in repetitive genome regions. Overall, these results suggest that changes in algal symbionts favouring heat tolerant associations are accompanied by changes in DNA methylation in the coral host. The implications of these results for coral adaptation, along with future avenues of research based on current knowledge gaps, are discussed in the present work.

Genetic Variations and Differential DNA Methylation to Face Contrasted Climates in Small Ruminants: An Analysis on Traditionally-Managed Sheep and Goats

The way in which living organisms mobilize a combination of long-term adaptive mechanisms and short-term phenotypic plasticity to face environmental variations is still largely unknown. In the context of climate change, understanding the genetic and epigenetic bases for adaptation and plasticity is a major stake for preserving genomic resources and the resilience capacity of livestock populations. We characterized both epigenetic and genetic variations by contrasting 22 sheep and 21 goats from both sides of a climate gradient, focusing on free-ranging populations from Morocco. We produced for each individual Whole-Genome Sequence at 12X coverage and MeDIP-Seq data, to identify regions under selection and those differentially methylated. For both species, the analysis of genetic differences (F_ST) along the genome between animals from localities with high vs. low temperature annual variations detected candidate genes under selection in relation to environmental perception (5 genes), immunity (4 genes), reproduction (8 genes) and production (11 genes). Moreover, we found for each species one differentially methylated gene, namely AGPTA4 in goat and SLIT3 in sheep, which were both related, among other functions, to milk production and muscle development. In both sheep and goats, the comparison between genomic regions impacted by genetic and epigenetic variations suggests that climatic variations impacted similar biological pathways but different genes.

Genome-wide DNA methylation profile changes associated with shell colouration in the Yesso scallop (Patinopecten yessoensis) as measured by whole-genome bisulfite sequencing

BACKGROUND

Mollusca, a phylum of highly rich species, possess vivid shell colours, but the underlying molecular mechanism remains to be elucidated. DNA methylation, one of the most common epigenetic modifications in eukaryotes, is believed to play a vital role in various biological processes. However, analysis of the effects of DNA methylation on shell colouration has rarely been performed in molluscs, limiting the current knowledge of the molecular mechanism of shell colour formation.

RESULTS

In the present study, to reveal the role of epigenetic regulation in shell colouration, WGBS, the "gold standard" of DNA methylation analysis, was first performed on the mantle tissues of Yesso scallops (Patinopecten yessoensis) with different shell colours (brown and white), and DNA methylomes at single-base resolution were generated. About 3% of cytosines were methylated in the genome of the Yesso scallop. A slight increase in ^mCG percentage and methylation level was found in brown scallops. Sequence preference of nearby methylated cytosines differed between high and low methylation level sites and between the brown- and white-shelled scallops. DNA methylation levels varied among the different genomic regions; all the detected regions in the brown group exhibited higher methylation levels than the white group. A total of 41,175 DMRs (differentially methylated regions) were detected between brown and white scallops. GO functions and pathways associated with the biosynthesis of melanin and porphyrins were significantly enriched for DMRs, among which several key shell colour-related genes were identified. Further, different correlations between mRNA expression levels and DNA methylation status were found in these genes, suggesting that DNA methylation regulates shell colouration in the Yesso scallop.

CONCLUSIONS

This study provides genome-wide DNA methylation landscapes of Yesso scallops with different shell colours, offering new insights into the epigenetic regulatory mechanism underlying shell colour.

Distribution, expression and methylation analysis of positively selected genes provides insights into the evolution in Brassica rapa

It is believed that positive selection is one of the major evolutionary forces underlying organism phenotypic diversification. Nevertheless, the characteristics of positively selected genes (PSGs), have not been well investigated. In this study, we performed a genome-wide analysis of orthologous genes between Brassica rapa (B. rapa) and Brassica oleracea (B. oleracea), and identified 468 putative PSGs. Our data show that, (1) PSGs are enriched in plant hormone signal transduction pathway and the transcription factor family; (2) PSGs are significantly lower expressed than randomly selected non-PSGs; (3) PSGs with tissue specificity are significantly higher expressed in the callus and reproductive tissues (flower and silique) than in vegetable tissues (root, stem and leaf); (4) the proportion of PSGs is positively correlated with the number of retained triplication gene copies, but the expression level of PSGs decay with the increasing of triplication gene copies; (5) the CG and CHG methylation levels of PSGs are significantly higher in introns and UTRs than in the promoter and exon regions; (6) the percent of transposable element is in proportion to the methylation level, and DNA methylation (especially in the CG content) has the tendency to reduce the expression of PSGs. This study provides insights into the characteristics, evolution, function, expression and methylation of PSGs in B. rapa.

Mapping epigenetic divergence in the massive radiation of Lake Malawi cichlid fishes

Epigenetic variation modulates gene expression and can be heritable. However, knowledge of the contribution of epigenetic divergence to adaptive diversification in nature remains limited. The massive evolutionary radiation of Lake Malawi cichlid fishes displaying extensive phenotypic diversity despite extremely low sequence divergence is an excellent system to study the epigenomic contribution to adaptation. Here, we present a comparative genome-wide methylome and transcriptome study, focussing on liver and muscle tissues in phenotypically divergent cichlid species. In both tissues we find substantial methylome divergence among species. Differentially methylated regions (DMR), enriched in evolutionary young transposons, are associated with transcription changes of ecologically-relevant genes related to energy expenditure and lipid metabolism, pointing to a link between dietary ecology and methylome divergence. Unexpectedly, half of all species-specific DMRs are shared across tissues and are enriched in developmental genes, likely reflecting distinct epigenetic developmental programmes. Our study reveals substantial methylome divergence in closely-related cichlid fishes and represents a resource to study the role of epigenetics in species diversification.

Combined Profiling of Transcriptome and DNA Methylome Reveal Genes Involved in Accumulation of Soluble Sugars and Organic Acid in Apple Fruits

Organic acids and soluble sugars are the major determinants of fruit organoleptic quality. Additionally, DNA methylation has crucial regulatory effects on various processes. However, the epigenetic modifications in the regulation of organic acid and soluble sugar accumulation in apple fruits remain uncharacterized. In this study, DNA methylation and the transcriptome were compared between ‘Honeycrisp’ and ‘Qinguan’ mature fruits, which differ significantly regarding soluble sugar and organic acid contents. In both ‘Honeycrisp’ and ‘Qinguan’ mature fruits, the CG context had the highest level of DNA methylation, and then CHG and CHH contexts. The number and distribution of differentially methylated regions (DMRs) varied among genic regions and transposable elements. The DNA methylation levels in all three contexts in the DMRs were significantly higher in ‘Honeycrisp’ mature fruits than in ‘Qinguan’ mature fruits. A combined methylation and transcriptome analysis revealed a negative correlation between methylation levels and gene expression in DMRs in promoters and gene bodies in the CG and CHG contexts and in gene bodies in the CHH context. Two candidate genes (MdTSTa and MdMa11), which encode tonoplast-localized proteins, potentially associated with fruit soluble sugar contents and acidity were identified based on expression and DNA methylation levels. Overexpression of MdTSTa in tomato increased the fruit soluble sugar content. Moreover, transient expression of MdMa11 in tobacco leaves significantly decreased the pH value. Our results reflect the diversity in epigenetic modifications influencing gene expression and will facilitate further elucidating the complex mechanism underlying fruit soluble sugar and organic acid accumulation.

Evolution of CG Methylation Maintenance Machinery in Plants

Cytosine methylation is an epigenetic mark present in most eukaryotic genomes that contributes to the regulation of gene expression and the maintenance of genome stability. DNA methylation mostly occurs at CG sequences, where it is initially deposited by de novo DNA methyltransferases and propagated by maintenance DNA methyltransferases (DNMT) during DNA replication. In this review, we first summarize the mechanisms maintaining CG methylation in mammals that involve the DNA Methyltransferase 1 (DNMT1) enzyme and its cofactor, UHRF1 (Ubiquitin-like with PHD and RING Finger domain 1). We then discuss the evolutionary conservation and diversification of these two core factors in the plant kingdom and speculate on potential functions of novel homologues typically observed in land plants but not in mammals.

Developmental programming of DNA methylation and gene expression patterns is associated with extreme cardiovascular tolerance to anoxia in the common snapping turtle

Background

Environmental fluctuation during embryonic and fetal development can permanently alter an organism’s morphology, physiology, and behaviour. This phenomenon, known as developmental plasticity, is particularly relevant to reptiles that develop in subterranean nests with variable oxygen tensions. Previous work has shown hypoxia permanently alters the cardiovascular system of snapping turtles and may improve cardiac anoxia tolerance later in life. The mechanisms driving this process are unknown but may involve epigenetic regulation of gene expression via DNA methylation. To test this hypothesis, we assessed in situ cardiac performance during 2 h of acute anoxia in juvenile turtles previously exposed to normoxia (21% oxygen) or hypoxia (10% oxygen) during embryogenesis. Next, we analysed DNA methylation and gene expression patterns in turtles from the same cohorts using whole genome bisulfite sequencing, which represents the first high-resolution investigation of DNA methylation patterns in any reptilian species.

Results

Genome-wide correlations between CpG and CpG island methylation and gene expression patterns in the snapping turtle were consistent with patterns observed in mammals. As hypothesized, developmental hypoxia increased juvenile turtle cardiac anoxia tolerance and programmed DNA methylation and gene expression patterns. Programmed differences in expression of genes such as SCN5A may account for differences in heart rate, while genes such as TNNT2 and TPM3 may underlie differences in calcium sensitivity and contractility of cardiomyocytes and cardiac inotropy. Finally, we identified putative transcription factor-binding sites in promoters and in differentially methylated CpG islands that suggest a model linking programming of DNA methylation during embryogenesis to differential gene expression and cardiovascular physiology later in life. Binding sites for hypoxia inducible factors (HIF1A, ARNT, and EPAS1) and key transcription factors activated by MAPK and BMP signaling (RREB1 and SMAD4) are implicated.

Conclusions

Our data strongly suggests that DNA methylation plays a conserved role in the regulation of gene expression in reptiles. We also show that embryonic hypoxia programs DNA methylation and gene expression patterns and that these changes are associated with enhanced cardiac anoxia tolerance later in life. Programming of cardiac anoxia tolerance has major ecological implications for snapping turtles, because these animals regularly exploit anoxic environments throughout their lifespan.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00414-7.

Molecular Docking and ADMET Prediction of Natural Compounds towards SARS Spike Glycoprotein-Human Angiotensin-Converting Enzyme 2 and SARS-CoV-2 Main Protease

More than a decade ago, a novel coronavirus that infects humans, bats, and certain other mammals termed severe acute respiratory syndrome coronavirus (SARS-CoV) caused an epidemic with ~ 10% case fatality, creating global panic and economic damage. Recently, another strain of the virus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), caused an infectious disease (COVID-19) in humans which was detected for the first time in Wuhan, China. Presently, there is no specific therapy available for the treatment of COVID-19. However, social distancing, patient isolation, and supportive medical care make up the current management for this current infectious disease pandemic. The present in silico study evaluated the binding affinities of some natural products (resveratrol, xylopic acid, ellagic acid, kaempferol, and quercetin) to human angiotensin-converting enzyme 2 and coronavirus (SARS-CoV-2) main protease compared to chloroquine, an inhibitor known to prevent cellular entry and replication of the coronavirus. The respective binding energies of the selected natural compounds and chloroquine towards the proteins were computed using PyRx virtual screening tool. The pharmacodynamic and pharmacokinetic attributes of the selected compounds were predicted using admetSAR. Molecular docking analysis showed that the natural compounds had better scores towards the selected protein compared to chloroquine with polar amino acid residues present at the binding sites. The predicted ADMET properties revealed the lower acute oral toxicity of the natural products compared to chloroquine. The study provides evidence suggesting that the relatively less toxic compounds from the natural sources could be repositioned as anti-viral agents to prevent the entry and replication of SARS-CoV-2.

Regulation of DNA methylation on key parasitism genes of Cysticercus cellulosae revealed by integrative epigenomic-transcriptomic analyses

Background

The life cycle of Taenia solium is characterized by different stages of development, requiring various kinds of hosts that can appropriately harbor the eggs (proglottids), the oncospheres, the larvae and the adults. Similar to other metazoan pathogens, T. solium undergoes transcriptional and developmental regulation via epigenetics during its complex lifecycle and host interactions.

Result

In the present study, we integrated whole-genome bisulfite sequencing and RNA-seq technologies to characterize the genome-wide DNA methylation and its effect on transcription of Cysticercus cellulosae of T. solium. We confirm that the T. solium genome in the cysticercus stage is epigenetically modified by DNA methylation in a pattern similar to that of other invertebrate genomes, i.e., sparsely or moderately methylated. We also observed an enrichment of non-CpG methylation in defined genetic elements of the T. solium genome. Furthermore, an integrative analysis of both the transcriptome and the DNA methylome indicated a strong correlation between these two datasets, suggesting that gene expression might be tightly regulated by DNA methylation. Importantly, our data suggested that DNA methylation might play an important role in repressing key parasitism-related genes, including genes encoding excretion-secretion proteins, thereby raising the possibility of targeting DNA methylation processes as a useful strategy in therapeutics of cysticercosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-021-00195-9.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

CHH Methylation Islands: A Nonconserved Feature of Grass Genomes That Is Positively Associated with Transposable Elements but Negatively Associated with Gene-Body Methylation

Methylated CHH (mCHH) islands are peaks of CHH methylation that occur primarily upstream to genes. These regions are actively targeted by the methylation machinery, occur at boundaries between heterochromatin and euchromatin, and tend to be near highly expressed genes. Here we took an evolutionary perspective by studying upstream mCHH islands across a sample of eight grass species. Using a statistical approach to define mCHH islands as regions that differ from genome-wide background CHH methylation levels, we demonstrated that mCHH islands are common and associate with 39% of genes, on average. We hypothesized that islands should be more frequent in genomes of large size, because they have more heterochromatin and hence more need for defined boundaries. We found, however, that smaller genomes tended to have a higher proportion of genes associated with 5' mCHH islands. Consistent with previous work suggesting that islands reflect the silencing of the edge of transposable elements (TEs), genes with nearby TEs were more likely to have mCHH islands. However, the presence of mCHH islands was not a function solely of TEs, both because the underlying sequences of islands were often not homologous to TEs and because genic properties also predicted the presence of 5' mCHH islands. These genic properties included length and gene-body methylation (gbM); in fact, in three of eight species, the absence of gbM was a stronger predictor of a 5' mCHH island than TE proximity. In contrast, gene expression level was a positive but weak predictor of the presence of an island. Finally, we assessed whether mCHH islands were evolutionarily conserved by focusing on a set of 2,720 orthologs across the eight species. They were generally not conserved across evolutionary time. Overall, our data establish additional genic properties that are associated with mCHH islands and suggest that they are not just a consequence of the TE silencing machinery.

Epigenetic regulation of melanogenesis

Melanogenesis is a complex process in which melanin is synthesized in melanocytes and transported to keratinocytes, which involves multiple genes and signaling pathways. Epigenetics refers to the potential genetic changes that affect gene expression without involving changes in the original sequence of DNA nucleotides. DNA methylation regulates the expression of key genes such as tyrosinase (TYR), tyrosinase-related protein 1 (TYRP1), dopachrome tautomerase (DCT) and microphthalmia-associated transcription factor (MITF), as well as paracrine factors such as stem cell factor (SCF) and endothelin-1 (ET-1) in melanogenesis. Potential DNA methylation sites are present in the genes of melanogenesis-related signaling pathways such as "Wnt", "PI3K/Akt/CREB" and "MAPK". H3K27 acetylation is abundant in melanogenesis-related genes. Both the upstream activation and downstream regulation of MITF depend on histone acetyltransferase CBP/p300, and pH-induced H3K27 acetylation may be the amplifying mechanism of MITF's effect. HDAC1 and HDAC10 catalyze histone deacetylation of melanogenesis-related gene promoters. Chromatin remodelers SWI/SNF complex and ISWI complex use the energy of ATP hydrolysis to rearrange nucleosomes, while their active subunits BRG1, BRM and BPTF, act as activators and cofactors of MITF. MicroRNAs (miRNAs) can directly target a large number of melanogenesis-related genes, while long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) regulate melanogenesis in a variety of ways. Interactions exist among the epigenetic mechanisms of melanogenesis. For example, the methyl CpG binding domain protein 2 (MeCP2) links DNA methylation, histone deacetylation, and histone methylation. Epigenetic-based therapy provides novel opportunities for treating dermatoses that are caused by pigmentation disturbances. This review summarizes the epigenetic regulation mechanisms of melanogenesis, and examines the pathogenesis and treatment of epigenetics in pigmentation disorders.

Repression of CHROMOMETHYLASE 3 prevents epigenetic collateral damage in Arabidopsis

DNA methylation has evolved to silence mutagenic transposable elements (TEs) while typically avoiding the targeting of endogenous genes. Mechanisms that prevent DNA methyltransferases from ectopically methylating genes are expected to be of prime importance during periods of dynamic cell cycle activities including plant embryogenesis. However, virtually nothing is known regarding how DNA methyltransferase activities are precisely regulated during embryogenesis to prevent the induction of potentially deleterious and mitotically stable genic epimutations. Here, we report that microRNA-mediated repression of CHROMOMETHYLASE 3 (CMT3) and the chromatin features that CMT3 prefers help prevent ectopic methylation of thousands of genes during embryogenesis that can persist for weeks afterwards. Our results are also consistent with CMT3-induced ectopic methylation of promoters or bodies of genes undergoing transcriptional activation reducing their expression. Therefore, the repression of CMT3 prevents epigenetic collateral damage on endogenous genes. We also provide a model that may help reconcile conflicting viewpoints regarding the functions of gene-body methylation that occurs in nearly all flowering plants.