The evolution of CHROMOMETHYLASES and gene body DNA methylation in plants

RNA directed DNA methylation and seed plant genome evolution

RNA Directed DNA Methylation (RdDM) is a pathway that mediates de novo DNA methylation, an evolutionary conserved chemical modification of cytosine bases, which exists in living organisms and utilizes small interfering RNA. Plants utilize DNA methylation for transposable element (TE) repression, regulation of gene expression and developmental regulation. TE activity strongly influences genome size and evolution, therefore making DNA methylation a key component in understanding divergence in genome evolution among seed plants. Multiple proteins that have extensively been studied in model plant Arabidopsis thaliana catalyze RNA dependent DNA Methylation pathway along with small interfering RNA. Several developmental functions have also been attributed to DNA methylation. This review will highlight aspects of RdDM pathway dynamics, evolution and functions in seed plants with focus on recent findings on conserved and non-conserved attributes between angiosperms and gymnosperms to potentially explain how methylation has impacted variations in evolutionary and developmental complexity among them and advance current understanding of this crucial epigenetic pathway.

RNA directed DNA methylation and seed plant genome evolution

RNA Directed DNA Methylation (RdDM) is a pathway that mediates de novo DNA methylation, an evolutionary conserved chemical modification of cytosine bases, which exists in living organisms and utilizes small interfering RNA. Plants utilize DNA methylation for transposable element (TE) repression, regulation of gene expression and developmental regulation. TE activity strongly influences genome size and evolution, therefore making DNA methylation a key component in understanding divergence in genome evolution among seed plants. Multiple proteins that have extensively been studied in model plant Arabidopsis thaliana catalyze RNA dependent DNA Methylation pathway along with small interfering RNA. Several developmental functions have also been attributed to DNA methylation. This review will highlight aspects of RdDM pathway dynamics, evolution and functions in seed plants with focus on recent findings on conserved and non-conserved attributes between angiosperms and gymnosperms to potentially explain how methylation has impacted variations in evolutionary and developmental complexity among them and advance current understanding of this crucial epigenetic pathway.

Principles of Epigenetic Homeostasis Shared Between Flowering Plants and Mammals

In diverse eukaryotes, epigenetic information such as DNA methylation is stably propagated over many cell divisions and generations, and can remain the same over thousands or millions of years. However, this stability is the product of dynamic processes that add and remove DNA methylation by specialized enzymatic pathways. The activities of these dynamic pathways must therefore be finely orchestrated in order to ensure that the DNA methylation landscape is maintained with high fidelity - a concept we term epigenetic homeostasis. In this review, we summarize recent insights into epigenetic homeostasis mechanisms in flowering plants and mammals, highlighting analogous mechanisms that have independently evolved to achieve the same goal of stabilizing the epigenetic landscape.

A blind and independent benchmark study for detecting differeally methylated regions in plants

MOTIVATION

Bisulfite sequencing (BS-seq) is a state-of-the-art technique for investigating methylation of the DNA to gain insights into the epigenetic regulation. Several algorithms have been published for identification of differentially methylated regions (DMRs). However, the performances of the individual methods remain unclear and it is difficult to optimally select an algorithm in application settings.

RESULTS

We analyzed BS-seq data from four plants covering three taxonomic groups. We first characterized the data using multiple summary statistics describing methylation levels, coverage, noise, as well as frequencies, magnitudes and lengths of methylated regions. Then, simulated datasets with most similar characteristics to real experimental data were created. Seven different algorithms (metilene, methylKit, MOABS, DMRcate, Defiant, BSmooth, MethylSig) for DMR identification were applied and their performances were assessed. A blind and independent study design was chosen to reduce bias and to derive practical method selection guidelines. Overall, metilene had superior performance in most settings. Data attributes, such as coverage and spread of the DMR lengths, were found to be useful for selecting the best method for DMR detection. A decision tree to select the optimal approach based on these data attributes is provided. The presented procedure might serve as a general strategy for deriving algorithm selection rules tailored to demands in specific application settings.

AVAILABILITY AND IMPLEMENTATION

Scripts that were used for the analyses and that can be used for prediction of the optimal algorithm are provided at https://github.com/kreutz-lab/DMR-DecisionTree. Simulated and experimental data are available at https://doi.org/10.6084/m9.figshare.11619045.

CONTACT

ckreutz@imbi.uni-freiburg.de.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Integrative transcriptomics and proteomics analysis constructs a new molecular model for ovule abortion in the female-sterile line of Pinus tabuliformis Carr

Ovule development is critical to plant reproduction and free nuclear mitosis of megagametophyte (FNMM) is vital for ovule development. However, most results of ovule development were based on the studies in angiosperms, and its molecular regulation remained largely unknown in gymnosperms, particularly, during FNMM. In this context, we studied the genome-wide difference between sterile line (SL) and fertile line (FL) ovules using transcriptomics and proteomics approaches in Pinus tabuliformis Carr. Comparative analyses revealed that genes involved in DNA replication, DNA damage repair, Cell cycle, Apoptosis and Energy metabolism were highlighted. Further results showed the low expressions of MCM 2-7, RRM1, etc. perhaps led to abnormal DNA replication and damage repair, and the significantly different expressions of PARP2, CCs1, CCs3, etc. implied that the accumulated DNA double-stranded breaks were failed to be repaired and the cell cycle was arrested at G2/M in SL ovules, potentially resulting in the occurrence of apoptosis. Moreover, the deficiency of ETF-QO might hinder FNMM. Consequently, FNMM stopped and ovule aborted in SL ovules. Our results suggested a selective regulatory mechanism led to FNMM half-stop and ovule abortion in P. tabuliformis and these insights could be exploited to investigate the molecular regulations of ovule development in woody gymnosperms.

Transposon age and non-CG methylation

Silencing of transposable elements (TEs) is established by small RNA-directed DNA methylation (RdDM). Maintenance of silencing is then based on a combination of RdDM and RNA-independent mechanisms involving DNA methyltransferase MET1 and chromodomain DNA methyltransferases (CMTs). Involvement of RdDM, according to this model should decrease with TE age but here we show a different pattern in tomato and Arabidopsis. In these species the CMTs silence long terminal repeat (LTR) transposons in the distal chromatin that are younger than those affected by RdDM. To account for these findings we propose that, after establishment of primary RdDM as in the original model, there is an RNA-independent maintenance phase involving CMTs followed by secondary RdDM. This progression of epigenetic silencing in the gene-rich distal chromatin is likely to influence the transcriptome either in cis or in trans depending on whether the mechanisms are RNA-dependent or -independent.

RNA-directed DNA methylation (RdDM) is thought to silence newly inserted transposable elements (TEs) with RNA-independent mechanisms becoming more prominent as TEs age. Here, the authors show that RdDM continues to silence the oldest intact distal TEs in tomato and Arabidopsis suggesting a second, later phase of RdDM.

A pseudomolecule-scale genome assembly of the liverwort Marchantia polymorpha

Marchantia polymorpha has recently become a prime model for cellular, evo-devo, synthetic biological, and evolutionary investigations. We present a pseudomolecule-scale assembly of the M. polymorpha genome, making comparative genome structure analysis and classical genetic mapping approaches feasible. We anchored 88% of the M. polymorpha draft genome to a high-density linkage map resulting in eight pseudomolecules. We found that the overall genome structure of M. polymorpha is in some respects different from that of the model moss Physcomitrella patens. Specifically, genome collinearity between the two bryophyte genomes and vascular plants is limited, suggesting extensive rearrangements since divergence. Furthermore, recombination rates are greatest in the middle of the chromosome arms in M. polymorpha like in most vascular plant genomes, which is in contrast with P. patens where recombination rates are evenly distributed along the chromosomes. Nevertheless, some other properties of the genome are shared with P. patens. As in P. patens, DNA methylation in M. polymorpha is spread evenly along the chromosomes, which is in stark contrast with the angiosperm model Arabidopsis thaliana, where DNA methylation is strongly enriched at the centromeres. Nevertheless, DNA methylation and recombination rate are anticorrelated in all three species. Finally, M. polymorpha and P. patens centromeres are of similar structure and marked by high abundance of retroelements unlike in vascular plants. Taken together, the highly contiguous genome assembly we present opens unexplored avenues for M. polymorpha research by linking the physical and genetic maps, making novel genomic and genetic analyses, including map-based cloning, feasible.

Natural variation in DNA methylation homeostasis and the emergence of epialleles

In plants and mammals, DNA methylation plays a critical role in transcriptional silencing by delineating heterochromatin from transcriptionally active euchromatin. A homeostatic balance between heterochromatin and euchromatin is essential to genomic stability. This is evident in many diseases and mutants for heterochromatin maintenance, which are characterized by global losses of DNA methylation coupled with localized ectopic gains of DNA methylation that alter transcription. Furthermore, we have shown that genome-wide methylation patterns in Arabidopsis thaliana are highly stable over generations, with the exception of rare epialleles. However, the extent to which natural variation in the robustness of targeting DNA methylation to heterochromatin exists, and the phenotypic consequences of such variation, remain to be fully explored. Here we describe the finding that heterochromatin and genic DNA methylation are highly variable among 725 A. thaliana accessions. We found that genic DNA methylation is inversely correlated with that in heterochromatin, suggesting that certain methylation pathway(s) may be redirected to genes upon the loss of heterochromatin. This redistribution likely involves a feedback loop involving the DNA methyltransferase, CHROMOMETHYLASE 3 (CMT3), H3K9me2, and histone turnover, as highly expressed, long genes with a high density of CMT3-preferred CWG sites are more likely to be methylated. Importantly, although the presence of CG methylation in genes alone may not affect transcription, genes containing CG methylation are more likely to become methylated at non-CG sites and silenced. These findings are consistent with the hypothesis that natural variation in DNA methylation homeostasis may underlie the evolution of epialleles that alter phenotypes.

Genome-wide investigation and transcriptional analysis of cytosine-5 DNA methyltransferase and DNA demethylase gene families in tea plant (Camellia sinensis) under abiotic stress and withering processing

DNA methylation is a highly conserved epigenetic modification involved in many biological processes, including growth and development, stress response, and secondary metabolism. In the plant kingdom, cytosine-5 DNA methyltransferase (C5-MTase) and DNA demethylase (dMTase) genes have been identified in some plant species. However, to the best of our knowledge, no investigator has focused on the identification and analysis of C5-MTase and dMTase genes in tea plants (Camellia sinensis) based on genome-wide levels. In this study, eight CsC5-MTases and four dMTases were identified in tea plants. These CsC5-MTase genes were divided into four subfamilies, including CsMET, CsCMT, CsDRM and CsDNMT2. The CsdMTase genes can be classified into CsROS, CsDME and CsDML. Based on conserved domain analysis of these genes, the gene loss and duplication events occurred during the evolution of CsC5-MTase and CsdMTase. Furthermore, multiple cis-acting elements were observed in the CsC5-MTase and CsdMTase, including light responsiveness, phytohormone responsiveness, stress responsiveness, and plant growth and development-related elements. Then, we investigated the transcript abundance of CsC5-MTase and CsdMTase under abiotic stress (cold and drought) and withering processing (white tea and oolong tea). Notably, most CsC5-MTases, except for CsCMT1 and CsCMT2, were significantly downregulated under abiotic stress, while the transcript abundance of all four CsdMTase genes was significantly induced. Similarly, the same transcript abundance of CsC5-MTase and CsdMTase was found during withering processing of white tea and oolong tea, respectively. In total, our findings will provide a basis for the roles of CsC5-MTase and CsdMTase in response to abiotic stress and the potential functions of these two gene families in affecting tea flavor during tea withering processing.

DNA Methylation Readers in Plants

In plants, DNA methylation occurs in distinct sequence contexts, including CG, CHG, and CHH. Thus, plants have developed a surprisingly diverse set of DNA methylation readers to cope with an extended repertoire of methylated sites. The Arabidopsis genome contains twelve Methyl-Binding Domain proteins (MBD), and nine SET and RING finger-associated (SRA) domain containing proteins belonging to the SUVH clade, in addition to three homologs of UHRF1, namely VIM1-3, all containing SRA domains. In this review, we will highlight several research questions that remain unresolved with respect to the function of plant DNA methylation readers, which can have both de novo demethylase and maintenance activity. We argue that maintenance of CG methylation in plants likely involved actors not found in their mammalian counterparts, and that new evidence suggests significant reprogramming of DNA methylation during plant reproduction as an important new development in the field.

Mesostigma viride Genome and Transcriptome Provide Insights into the Origin and Evolution of Streptophyta

The Streptophyta include unicellular and multicellular charophyte green algae and land plants. Colonization of the terrestrial habitat by land plants is a major evolutionary event that has transformed the planet. So far, lack of genome information on unicellular charophyte algae hinders the understanding of the origin and the evolution from unicellular to multicellular life in Streptophyta. This work reports the high-quality reference genome and transcriptome of Mesostigma viride, a single-celled charophyte alga with a position at the base of Streptophyta. There are abundant segmental duplications and transposable elements in M. viride, which contribute to a relatively large genome with high gene content compared to other algae and early diverging land plants. This work identifies the origin of genetic tools that multicellular Streptophyta have inherited and key genetic innovations required for the evolution of land plants from unicellular aquatic ancestors. The findings shed light on the age-old questions of the evolution of multicellularity and the origin of land plants.

Establishment, maintenance, and biological roles of non-CG methylation in plants

Cytosine DNA methylation is prevalent throughout eukaryotes and prokaryotes. While most commonly thought of as being localized to dinucleotide CpG sites, non-CG sites can also be modified. Such non-CG methylation is widespread in plants, occurring at trinucleotide CHG and CHH (H = A, T, or C) sequence contexts. The prevalence of non-CG methylation in plants is due to the plant-specific CHROMOMETHYLASE (CMT) and RNA-directed DNA Methylation (RdDM) pathways. These pathways have evolved through multiple rounds of gene duplication and gene loss, generating epigenomic variation both within and between species. They regulate both transposable elements and genes, ensure genome integrity, and ultimately influence development and environmental responses. In these capacities, non-CG methylation influence and shape plant genomes.

Evolution of DNA Methylome Diversity in Eukaryotes

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

Interspecific variation across angiosperms in global DNA methylation: phylogeny, ecology and plant features in tropical and Mediterranean communities

The interspecific range of epigenetic variation and the degree to which differences between angiosperm species are related to geography, evolutionary history, ecological settings or species-specific traits, remain essentially unexplored. Genome-wide global DNA cytosine methylation is a tractable 'epiphenotypic' feature suitable for exploring these relationships. Global cytosine methylation was estimated in 279 species from two distant, ecologically disparate geographical regions: Mediterranean Spain and tropical México. At each region, four distinct plant communities were analyzed. Global methylation spanned a 10-fold range among species (4.8-42.2%). Interspecific differences were related to evolutionary trajectories, as denoted by a strong phylogenetic signal. Genomes of tropical species were on average less methylated than those of Mediterranean ones. Woody plants have genomes with lower methylation than perennial herbs, and genomes of widespread species were less methylated than those of species with restricted geographical distribution. The eight communities studied exhibited broad and overlapping interspecific variances in global cytosine methylation and only two of them differed in average methylation. Altogether, our broad taxonomic survey supported global methylation as a plant 'epiphenotypic' trait largely associated with species evolutionary history, genome size, range size and woodiness. Additional studies are required for better understanding the environmental components underlying local and geographical variation.

Genome-Wide Reinforcement of DNA Methylation Occurs during Somatic Embryogenesis in Soybean

Somatic embryogenesis is an important tissue culture technique that sometimes leads to phenotypic variation via genetic and/or epigenetic changes. To understand the genomic and epigenomic impacts of somatic embryogenesis, we characterized soybean (Glycine max) epigenomes sampled from embryos at 10 different stages ranging from 6 weeks to 13 years of continuous culture. We identified genome-wide increases in DNA methylation from cultured samples, especially at CHH sites. The hypermethylation almost exclusively occurred in regions previously possessing non-CG methylation and was accompanied by increases in the expression of genes encoding the RNA-directed DNA methylation (RdDM) machinery. The epigenomic changes were similar between somatic and zygotic embryogenesis. Following the initial global wave of hypermethylation, rare decay events of maintenance methylation were observed, and the extent of the decay increased with time in culture. These losses in DNA methylation were accompanied by downregulation of genes encoding the RdDM machinery and transcriptome reprogramming reminiscent of transcriptomes during late-stage seed development. These results reveal a process for reinforcing already silenced regions to maintain genome integrity during somatic embryogenesis over the short term, which eventually decays at certain loci over longer time scales.

Sex-specific changes in the aphid DNA methylation landscape

Aphids present an ideal system to study epigenetics as they can produce diverse, but genetically identical, morphs in response to environmental stimuli. Here, using whole genome bisulphite sequencing and transcriptome sequencing of the green peach aphid (Myzus persicae), we present the first detailed analysis of cytosine methylation in an aphid and investigate differences in the methylation and transcriptional landscapes of male and asexual female morphs. We found that methylation primarily occurs in a CG dinucleotide (CpG) context and that exons are highly enriched for methylated CpGs, particularly at the 3' end of genes. Methylation is positively associated with gene expression, and methylated genes are more stably expressed than unmethylated genes. Male and asexual female morphs have distinct methylation profiles. Strikingly, these profiles are divergent between the sex chromosome and the autosomes; autosomal genes are hypomethylated in males compared to asexual females, whereas genes belonging to the sex chromosome, which is haploid in males, are hypermethylated. Overall, we found correlated changes in methylation and gene expression between males and asexual females, and this correlation was particularly strong for genes located on the sex chromosome. Our results suggest that differential methylation of sex-biased genes plays a role in aphid sexual differentiation.

Dynamics and function of DNA methylation in plants

DNA methylation is a conserved epigenetic modification that is important for gene regulation and genome stability. Aberrant patterns of DNA methylation can lead to plant developmental abnormalities. A specific DNA methylation state is an outcome of dynamic regulation by de novo methylation, maintenance of methylation and active demethylation, which are catalysed by various enzymes that are targeted by distinct regulatory pathways. In this Review, we discuss DNA methylation in plants, including methylating and demethylating enzymes and regulatory factors, and the coordination of methylation and demethylation activities by a so-called methylstat mechanism; the functions of DNA methylation in regulating transposon silencing, gene expression and chromosome interactions; the roles of DNA methylation in plant development; and the involvement of DNA methylation in plant responses to biotic and abiotic stress conditions.

Evolutionary and Experimental Loss of Gene Body Methylation and Its Consequence to Gene Expression

In flowering plants, gene body methylation (gbM) is associated with a subset of constitutively expressed genes. It has been proposed that gbM modulates gene expression. Here, we show that there are no consistent and direct differences to expression following the loss of gbM. By comparing expression of gbM genes in Arabidopsis thaliana accessions to orthologous genes in two Eutrema salsugineum genotypes, we identified both positive and negative expression differences associated with gbM loss. However, expression is largely unaffected by gbM loss in E. salsugineum Expression differences between species were within the variation of expression observed within A. thaliana accessions that displayed variation in gbM. Furthermore, experimentally induced loss of gbM did not consistently lead to differences in expression compared to wild type. To date, there is no convincing data to support a direct causal link between the presence/absence of gbM and the modulation of expression in flowering plants.

DNA Methylation: Shared and Divergent Features across Eukaryotes

Chemical modification of nucleotide bases in DNA provides one mechanism for conveying information in addition to the genetic code. 5-methylcytosine (5mC) represents the most common chemically modified base in eukaryotic genomes. Sometimes referred to simply as DNA methylation, in eukaryotes 5mC is most prevalent at CpG dinucleotides and is frequently associated with transcriptional repression of transposable elements. However, 5mC levels and distributions are variable across phylogenies, and emerging evidence suggests that the functions of DNA methylation may be more diverse and complex than was previously appreciated. We summarize the current understanding of DNA methylation profiles and functions in different eukaryotic lineages.

Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation

Comparative genomics can unravel the genetic basis of species differences; however, successful reports on quantitative traits are still scarce. Here we present genome assemblies of 31 so-far unassembled Brassicaceae plant species and combine them with 16 previously published assemblies to establish the Brassicaceae Diversity Panel. Using a new interspecies association strategy for quantitative traits, we found a so-far unknown association between the unexpectedly high variation in CG to TG substitution rates in genes and the absence of CHROMOMETHYLASE3 (CMT3) orthologues. Low substitution rates were associated with the loss of CMT3, while species with conserved CMT3 orthologues showed high substitution rates. Species without CMT3 also lacked gene-body methylation (gbM), suggesting an evolutionary trade-off between the unknown function of gbM and low substitution rates in Brassicaceae, possibly due to low mutability of non-methylated cytosines.

Epimutations are associated with CHROMOMETHYLASE 3-induced de novo DNA methylation

In many plant species, a subset of transcribed genes are characterized by strictly CG-context DNA methylation, referred to as gene body methylation (gbM). The mechanisms that establish gbM are unclear, yet flowering plant species naturally without gbM lack the DNA methyltransferase, CMT3, which maintains CHG (H = A, C, or T) and not CG methylation at constitutive heterochromatin. Here, we identify the mechanistic basis for gbM establishment by expressing CMT3 in a species naturally lacking CMT3. CMT3 expression reconstituted gbM through a progression of de novo CHG methylation on expressed genes, followed by the accumulation of CG methylation that could be inherited even following loss of the CMT3 transgene. Thus, gbM likely originates from the simultaneous targeting of loci by pathways that promote euchromatin and heterochromatin, which primes genes for the formation of stably inherited epimutations in the form of CG DNA methylation.

Epimutations are associated with CHROMOMETHYLASE 3-induced de novo DNA methylation

In many plant species, a subset of transcribed genes are characterized by strictly CG-context DNA methylation, referred to as gene body methylation (gbM). The mechanisms that establish gbM are unclear, yet flowering plant species naturally without gbM lack the DNA methyltransferase, CMT3, which maintains CHG (H = A, C, or T) and not CG methylation at constitutive heterochromatin. Here, we identify the mechanistic basis for gbM establishment by expressing CMT3 in a species naturally lacking CMT3. CMT3 expression reconstituted gbM through a progression of de novo CHG methylation on expressed genes, followed by the accumulation of CG methylation that could be inherited even following loss of the CMT3 transgene. Thus, gbM likely originates from the simultaneous targeting of loci by pathways that promote euchromatin and heterochromatin, which primes genes for the formation of stably inherited epimutations in the form of CG DNA methylation.

Tracing the origin and evolution history of methylation-related genes in plants

BACKGROUND

DNA methylation is a crucial epigenetic modification, which is involved in many biological processes, including gene expression regulation, embryonic development, cell differentiation and genomic imprinting etc. And it also involves many key regulatory genes in eukaryotes. By tracing the evolutionary history of methylation-related genes, we can understand the origin and expansion time of these genes, which helps to understand the evolutionary history of plants, and we can also understand the changes of DNA methylation patterns in different species. However, most studies on the evolution of methylation-related genes failed to be carried out for the whole DNA methylation pathway.

RESULTS

In this study, we conducted a comprehensive identification of 33 methylation-related genes in 77 species, and investigated gene origin and evolution throughout the plant kingdom. We found that the origin of genes responsible for methylation maintenance and demethylation evolved early, while most de novo methylation-related genes appeared late. The methylation-related genes were expanded by whole genome duplication and tandem replication, but were also accompanied by a large number of gene absence events in different species. The gene length and intron length varied a lot in different species, but exon structure and functional domains were relatively conserved. The phylogenetic relationships of methylation-related genes were traced to reveal the evolution history of DNA methylation in different species. The expression patterns of methylation-related genes have changed during the evolution of species, and the expression patterns of these genes in different species can be clustered into four categories.

CONCLUSIONS

The study describes a global characterization of DNA methylation-related genes in the plant kingdom. The similarities and differences in origin time, gene structure and phylogenetic relationship of these genes lead us to understand the evolutionary conservation and dynamics of DNA methylation in plants.

DNA methylation in Marchantia polymorpha

Methylation of DNA is an epigenetic mechanism for the control of gene expression. Alterations in the regulatory pathways involved in the establishment, perpetuation and removal of DNA methylation can lead to severe developmental alterations. Our understanding of the mechanistic aspects and relevance of DNA methylation comes from remarkable studies in well-established angiosperm plant models including maize and Arabidopsis. The study of plant models positioned at basal lineages opens exciting opportunities to expand our knowledge on the function and evolution of the components of DNA methylation. In this Tansley Insight, we summarize current progress in our understanding of the molecular basis and relevance of DNA methylation in the liverwort Marchantia polymorpha.

Loss of Gene Body Methylation in Eutrema salsugineum Is Associated with Reduced Gene Expression

Gene body methylation (gbM) is typically characterized by DNA methylation in the CG context within coding regions and is associated with constitutive genes that have moderate to high expression levels. A recent study discovered the loss of gbM in two plant species (Eutrema salsugineum and Conringia planisiliqua), illustrating that gbM is not necessary for survival and reproduction. The same paper stated there was no detectable effect of gbM loss on gene expression (GE). Here, we reanalyzed the GE data and accounted for experimental variability in expression level estimates. We show that the loss of gbM in E. salsugineum is associated with a small but highly significant decrease in GE relative to the closely related species Arabidospis thaliana. Our results are consistent with various evolutionary analyses that suggest gbM has a function, perhaps as a homeostatic effect on GE.

Decrease in DNA methylation 1 (DDM1) is required for the formation of m CHH islands in maize

DNA methylation plays a crucial role in suppressing mobilization of transposable elements and regulation of gene expression. A number of studies have indicated that DNA methylation pathways and patterns exhibit distinct properties in different species, including Arabidopsis, rice, and maize. Here, we characterized the function of DDM1 in regulating genome-wide DNA methylation in maize. Two homologs of ZmDDM1 are abundantly expressed in the embryo and their simultaneous disruption caused embryo lethality with abnormalities in cell proliferation from the early stage of kernel development. We establish that ZmDDM1 is critical for DNA methylation, at CHG sites, and to a lesser extent at CG sites, in heterochromatic regions, and unexpectedly, it is required for the formation of ^m CHH islands. In addition, ZmDDM1 is indispensable for the presence of 24-nt siRNA, suggesting its involvement in the RdDM pathway. Our results provide novel insight into the role of ZmDDM1 in regulating the formation of ^m CHH islands, via the RdDM pathway maize, suggesting that, in comparison to Arabidopsis, maize may have adopted distinct mechanisms for regulating ^m CHH.

Silencing of NbCMT3s has Pleiotropic Effects on Development by Interfering with Autophagy-Related Genes in Nicotiana benthamiana

DNA methylation is a chromatin mark that has a crucial role in regulating gene expression. The chromomethylase (CMT) protein family is a plant-specific DNA methyltransferase that mediates growth and development. However, the roles of CMT3 in autophagy remain to be elucidated. Here, we identified the potential targets of CMT3 in Nicotiana benthamiana (NbCMT3) during developmental programs. Virus-induced gene silencing of NbCMT3/3-2 in N. benthamiana had pleiotropic effects on plant morphology, which indicates its indispensible role in development. Genome-wide transcriptome analysis of NbCMT3/3-2-silenced plants revealed interference with genes related to autophagy and ubiquitination. The expression of NbBeclin 1 and NbHRD1B was higher in NbCMT3/3-2-silenced than control plants. The formation of autophagosomes and starch degradation was disrupted in NbCMT3/3-2-silenced plants, which implies a perturbed autophagic processes. We further generated transgenic N. benthamiana plants carrying a chimeric promoter-reporter construct linking the NbBeclin 1 promoter region and β-glucuronidase (GUS) reporter (pNbBeclin::GUS). NbBeclin 1 promoter activity was significantly enhanced in NbCMT3/3-2-silenced plants. Thus, NbCMT3/3-2 silencing had pleiotropic effects on development by interfering with NbBeclin 1 expression and autophagy-related processes.

RdDM-independent de novo and heterochromatin DNA methylation by plant CMT and DNMT3 orthologs

To properly regulate the genome, cytosine methylation is established by animal DNA methyltransferase 3 s (DNMT3s). While altered DNMT3 homologs, Domains rearranged methyltransferases (DRMs), have been shown to establish methylation via the RNA directed DNA methylation (RdDM) pathway, the role of true-plant DNMT3 orthologs remains elusive. Here, we profile de novo (RPS transgene) and genomic methylation in the basal plant, Physcomitrella patens, mutated in each of its PpDNMTs. We show that PpDNMT3b mediates CG and CHH de novo methylation, independently of PpDRMs. Complementary de novo CHG methylation is specifically mediated by the CHROMOMETHYLASE, PpCMT. Intragenomically, PpDNMT3b functions preferentially within heterochromatin and is affected by PpCMT. In comparison, PpDRMs target active-euchromatic transposons. Overall, our data resolve how DNA methylation in plants can be established in heterochromatin independently of RdDM; suggest that DRMs have emerged to target euchromatin; and link DNMT3 loss in angiosperms to the initiation of heterochromatic CHH methylation by CMT2.

Whether plants have true DNMT3 orthologs and their role in establishing DNA methylation are still unclear. Here, the authors show that DNMT3s are persistent through plant evolution and mediates both de novo and heterochromatin DNA methylation in the early divergent land plant Physcomitrella patens.

Hypomethylated drm1 drm2 cmt3 mutant phenotype of Arabidopsis thaliana is related to auxin pathway impairment

DNA methylation carried out by different methyltransferase classes is a relevant epigenetic modification of DNA which plays a relevant role in the development of eukaryotic organisms. Accordingly, in Arabidopsis thaliana loss of DNA methylation due to combined mutations in genes encoding for DNA methyltransferases causes several developmental abnormalities. The present study describes novel growth disorders in the drm1 drm2 cmt3 triple mutant of Arabidopsis thaliana, defective both in maintenance and de novo DNA methylation, and highlights the correlation between DNA methylation and the auxin hormone pathway. By using an auxin responsive reporter gene, we discovered that auxin accumulation and distribution were affected in the mutant compared to the wild type, from embryo to adult plant stage. In addition, we demonstrated that the defective methylation status also affected the expression of genes that regulate auxin hormone pathways from synthesis to transport and signalling and a direct relationship between differentially expressed auxin-related genes and altered auxin accumulation and distribution in embryo, leaf and root was observed. Finally, we provided evidence of the direct and organ-specific modulation of auxin-related genes through the DNA methylation process.

Diversity of cytosine methylation across the fungal tree of life

The generation of thousands of fungal genomes is leading to a better understanding of genes and genomic organization within the kingdom. However, the epigenome, which includes DNA and chromatin modifications, remains poorly investigated in fungi. Large comparative studies in animals and plants have deepened our understanding of epigenomic variation, particularly of the modified base 5-methylcytosine (5mC), but taxonomic sampling of disparate groups is needed to develop unifying explanations for 5mC variation. Here we utilize the largest phylogenetic resolution of 5mC methyltransferases (5mC MTases) and genome evolution to better understand levels and patterns of 5mC across fungi. We show that extant 5mC MTase genotypes are descendent from ancestral maintenance and de novo genotypes, whereas the 5mC MTases DIM-2 and RID are more recently derived, and that 5mC levels are correlated with 5mC MTase genotype and transposon content. Our survey also revealed that fungi lack canonical gene body methylation, which distinguishes fungal epigenomes from certain insect and plant species. However, some fungal species possess independently derived clusters of contiguous 5mC encompassing many genes. In some cases, DNA repair pathways and the N6-methyladenine (6mA) DNA modification negatively coevolved with 5mC pathways, which additionally contributed to interspecific epigenomic variation across fungi.

Identification and characterisation of hypomethylated DNA loci controlling quantitative resistance in Arabidopsis

Variation in DNA methylation enables plants to inherit traits independently of changes to DNA sequence. Here, we have screened an Arabidopsis population of epigenetic recombinant inbred lines (epiRILs) for resistance against Hyaloperonospora arabidopsidis (Hpa). These lines share the same genetic background, but show variation in heritable patterns of DNA methylation. We identified four epigenetic quantitative trait loci (epiQTLs) that provide quantitative resistance without reducing plant growth or resistance to other (a)biotic stresses. Phenotypic characterisation and RNA-sequencing analysis revealed that Hpa-resistant epiRILs are primed to activate defence responses at the relatively early stages of infection. Collectively, our results show that hypomethylation at selected pericentromeric regions is sufficient to provide quantitative disease resistance, which is associated with genome-wide priming of defence-related genes. Based on comparisons of global gene expression and DNA methylation between the wild-type and resistant epiRILs, we discuss mechanisms by which the pericentromeric epiQTLs could regulate the defence-related transcriptome.

In plants, animals and microbes genetic information is encoded by DNA, which are made up of sequences of building blocks, called nucleotide bases. These sequences can be separated into sections known as genes that each encode specific traits. It was previously thought that only changes to the sequence of bases in a DNA molecule could alter the traits passed on to future generations. However, it has recently become clear that some traits can also be inherited through modifications to the DNA that do not alter its sequence.

One such modification is to attach a tag, known as a methyl group, to a nucleotide base known as cytosine. These methyl tags can be added to, or removed from, DNA to create different patterns of methylation. Previous studies have shown that plants whose DNA is less methylated than normal (‘hypo-methylated’) are more resistant to plant diseases. However, the location and identity of the hypo-methylated DNA regions controlling this resistance remained unknown. To address this problem, Furci, Jain et al. studied how DNA methylation in a small weed known as Arabidopsis thaliana affects how well the plants can resist a disease known as downy mildew.

Furci, Jain et al. studied a population of over 100 A. thaliana lines that have the same DNA sequences but different patterns of DNA methylation. The experiments identified four DNA locations that were less methylated in lines with enhanced resistance to downy mildew. Importantly, this form of resistance did not appear to reduce how well the plants grew, or make them less able to resist other diseases or environmental stresses.

The results of further experiments suggested that reduced methylation at the four DNA regions prime the plant’s immune system, enabling a faster and stronger activation of a multitude of defence genes across the genome after attack by downy mildew.

The next steps following on from this work are to investigate exactly how the four DNA regions with reduced methylation can prime so many different defence genes in the plant. Further research is also needed to determine whether it is possible to breed crop plants with lower levels of methylation at specific DNA locations to improve disease resistance, but without decreasing the amount and quality of food produced.

Natural variation of DNA methylation and gene expression may determine local adaptations of Scots pine populations

Long-lived conifers are vulnerable to climate change because classical evolutionary processes are slow in developing adaptive responses. Therefore, the capacity of a genotype to adopt different phenotypes is important. Gene expression is the primary mechanism that converts genome-encoded information into phenotypes, and DNA methylation is employed in the epigenetic regulation of gene expression. We investigated variations in global DNA methylation and gene expression between three Scots pine (Pinus sylvestris L.) populations located in northern and southern Finland using mature seeds. Gene expression levels were studied in six DNA methyltransferase (DNMT) genes, which were characterized in this study, and in 19 circadian clock genes regulating adaptive traits. In embryos, expression diversity was found for three DNMT genes, which maintain DNA methylation. The expression of two DNMT genes was strongly correlated with climate variables, which suggests a role for DNA methylation in local adaptation. For adaptation-related genes, expression levels showed between-population variation in 11 genes in megagametophytes and in eight genes in embryos, and many of these genes were linked to climate factors. Altogether, our results suggest that differential DNA methylation and gene expression contribute to local adaptation in Scots pine populations and may enhance the fitness of trees under rapidly changing climatic conditions.

Detection of somatic epigenetic variation in Norway spruce via targeted bisulfite sequencing

Epigenetic mechanisms represent a possible mechanism for achieving a rapid response of long-lived trees to changing environmental conditions. However, our knowledge on plant epigenetics is largely limited to a few model species. With increasing availability of genomic resources for many tree species, it is now possible to adopt approaches from model species that permit to obtain single-base pair resolution data on methylation at a reasonable cost. Here, we used targeted bisulfite sequencing (TBS) to study methylation patterns in the conifer species Norway spruce (Picea abies). To circumvent the challenge of disentangling epigenetic and genetic differences, we focused on four clone pairs, where clone members were growing in different climatic conditions for 24 years. We targeted >26.000 genes using TBS and determined the performance and reproducibility of this approach. We characterized gene body methylation and compared methylation patterns between environments. We found highly comparable capture efficiency and coverage across libraries. Methylation levels were relatively constant across gene bodies, with 21.3 ± 0.3%, 11.0 ± 0.4% and 1.3 ± 0.2% in the CG, CHG, and CHH context, respectively. The variance in methylation profiles did not reveal consistent changes between environments, yet we could identify 334 differentially methylated positions (DMPs) between environments. This supports that changes in methylation patterns are a possible pathway for a plant to respond to environmental change. After this successful application of TBS in Norway spruce, we are confident that this approach can contribute to broaden our knowledge of methylation patterns in natural tree populations.

Long-term evolutionary DNA methylation dynamic of protein-coding genes and its underlying mechanism

DNA methylation is an important type of epigenetic modifications for the maintenance of genome functionality and stability. Although there are many studies on DNA methylation patterns, mechanisms, and functions, no study has focused on the evolutionary dynamic of DNA methylation. Here, we present the first genome-wide pattern of evolutionary DNA methylation dynamic in protein-coding genes, by grouping the Arabidopsis thaliana protein-coding genes into several conservation levels representing different evolutionary ages, and by investigating their DNA methylation features for three methylation contexts in both genic and flanking regions. The main results include: in a long-term evolutionary period, (1) genic CHG and CHH methylation levels tend to be decreased over time, which is mainly due to the reductions in the number of siRNA target sites in genes; (2) genic CG methylation levels are firstly reduced and then increased on average over evolutionary time, which is the interactional result of increased proportion and decreased CG methylation level of CG methylated genes; and (3) increased gene length and the stochastic methylation mechanism in CG context may further account for genic CG methylation trend in evolution. The diverse DNA methylation mechanisms in different contexts, together with altered gene length in evolution, could interpret the methylation dynamic of protein-coding genes over evolutionary time. This evolutionary perspective provides a dynamic understanding of the intrinsic relationship between DNA methylation and its functional and evolutionary effects on the genomes.

DNA methylation and genetic degeneration of the Y chromosome in the dioecious plant Silene latifolia

Background

S. latifolia is a model organism for the study of sex chromosome evolution in plants. Its sex chromosomes include large regions in which recombination became gradually suppressed. The regions tend to expand over time resulting in the formation of evolutionary strata. Non-recombination and later accumulation of repetitive sequences is a putative cause of the size increase in the Y chromosome. Gene decay and accumulation of repetitive DNA are identified as key evolutionary events. Transposons in the X and Y chromosomes are distributed differently and there is a regulation of transposon insertion by DNA methylation of the target sequences, this points to an important role of DNA methylation during sex chromosome evolution in Silene latifolia. The aim of this study was to elucidate whether the reduced expression of the Y allele in S. latifolia is caused by genetic degeneration or if the cause is methylation triggered by transposons and repetitive sequences.

Results

Gene expression analysis in S. latifolia males has shown expression bias in both X and Y alleles. To determine whether these differences are caused by genetic degeneration or methylation spread by transposons and repetitive sequences, we selected several sex-linked genes with varying degrees of degeneration and from different evolutionary strata. Immunoprecipitation of methylated DNA (MeDIP) from promoter, exon and intron regions was used and validated through bisulfite sequencing. We found DNA methylation in males, and only in the promoter of genes of stratum I (older). The Y alleles in genes of stratum I were methylation enriched compared to X alleles. There was also abundant and high percentage methylation in the CHH context in most sequences, indicating de novo methylation through the RdDM pathway.

Conclusions

We speculate that TE accumulation and not gene decay is the cause of DNA methylation in the S. latifolia Y sex chromosome with influence on the process of heterochromatinization.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4936-y) contains supplementary material, which is available to authorized users.

Loss of RNA-Directed DNA Methylation in Maize Chromomethylase and DDM1-Type Nucleosome Remodeler Mutants

Plants make use of distinct types of DNA methylation characterized by their DNA methyltransferases and modes of regulation. One type, RNA-directed DNA methylation (RdDM), is guided by small interfering RNAs (siRNAs) to the edges of transposons that are close to genes, areas called mCHH islands in maize (Zea mays). Another type, chromomethylation, is guided by histone H3 lysine 9 methylation to heterochromatin across the genome. We examined DNA methylation and small RNA expression in plant tissues that were mutant for both copies of the genes encoding chromomethylases as well as mutants for both copies of the genes encoding DECREASED DNA METHYLATION1 (DDM1)-type nucleosome remodelers, which facilitate chromomethylation. Both sets of double mutants were nonviable but produced embryos and endosperm. RdDM was severely compromised in the double mutant embryos, both in terms of DNA methylation and siRNAs. Loss of 24-nucleotide siRNA from mCHH islands was coupled with a gain of 21-, 22-, and 24-nucleotide siRNAs in heterochromatin. These results reveal a requirement for both chromomethylation and DDM1-type nucleosome remodeling for RdDM in mCHH islands, which we hypothesize is due to dilution of RdDM components across the genome when heterochromatin is compromised.

Subtle Perturbations of the Maize Methylome Reveal Genes and Transposons Silenced by Chromomethylase or RNA-Directed DNA Methylation Pathways

DNA methylation is a chromatin modification that can provide epigenetic regulation of gene and transposon expression. Plants utilize several pathways to establish and maintain DNA methylation in specific sequence contexts. The chromomethylase (CMT) genes maintain CHG (where H = A, C or T) methylation. The RNA-directed DNA methylation (RdDM) pathway is important for CHH methylation. Transcriptome analysis was performed in a collection of Zea mays lines carrying mutant alleles for CMT or RdDM-associated genes. While the majority of the transcriptome was not affected, we identified sets of genes and transposon families sensitive to context-specific decreases in DNA methylation in mutant lines. Many of the genes that are up-regulated in CMT mutant lines have high levels of CHG methylation, while genes that are differentially expressed in RdDM mutants are enriched for having nearby mCHH islands, implicating context-specific DNA methylation in the regulation of expression for a small number of genes. Many genes regulated by CMTs exhibit natural variation for DNA methylation and transcript abundance in a panel of diverse inbred lines. Transposon families with differential expression in the mutant genotypes show few defining features, though several families up-regulated in RdDM mutants show enriched expression in endosperm tissue, highlighting the potential importance for this pathway during reproduction. Taken together, our findings suggest that while the number of genes and transposon families whose expression is reproducibly affected by mild perturbations in context-specific methylation is small, there are distinct patterns for loci impacted by RdDM and CMT mutants.

Identification and characterization of the cytosine-5 DNA methyltransferase gene family in Salvia miltiorrhiza

Cytosine DNA methylation is highly conserved epigenetic modification involved in a wide range of biological processes in eukaryotes. It was established and maintained by cytosine-5 DNA methyltransferases (C5-MTases) in plants. Through genome-wide identification, eight putative SmC5-MTase genes were identified from the genome of Salvia miltiorrhiza, a well-known traditional Chinese medicine material and an emerging model medicinal plant. Based on conserved domains and phylogenetic analysis, eight SmC5-MTase genes were divided into four subfamilies, including MET, CMT, DRM and DNMT2. Genome-wide comparative analysis of the C5-MTase gene family in S. miltiorrhiza and Arabidopsis thaliana, including gene structure, sequence features, sequence alignment and conserved motifs, was carried out. The results showed conservation and divergence of the members of each subfamily in plants. The length of SmC5-MTase open reading frames ranges widely from 1,152 (SmDNMT2) to 5,034 bp (SmMET1). The intron number of SmC5-MTases varies between 7 (SmDRM1) and 20 (SmCMT1 and SmCMT2b). These features were similar to their counterparts from Arabidopsis. Sequence alignment and conserved motif analysis showed the existence of highly conserved and subfamily-specific motifs in the C5-MTases analyzed. Differential transcript abundance was detected for SmC5-MTases, implying genome-wide variance of DNA methylation in different organs and tissues. Transcriptome-wide analysis showed that the transcript levels of all SmC5-MTase genes was slightly changed under yeast extract and methyl jasmonate treatments. Six SmC5-MTases, including SmMET1, SmCMT1, SmCMT2a, SmCMT2b, SmCMT3 and SmDRM1, were salicylic acid-responsive, suggesting the involvement of SmC5-MTases in salicylic acid-dependent immunity. These results provide useful information for demonstrating the role of DNA methylation in bioactive compound biosynthesis and Dao-di herb formation in medicinal plants.

Specifications of Targeting Heterochromatin Modifications in Plants

Plants encode a diverse repertoire of DNA methyltransferases that have specialized to target cytosines for methylation in specific sequence contexts. These include the de novo methyltransferase, DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), which methylates cytosines in all sequence contexts through an RNA-guided process, the CHROMOMETHYLASES (CMTs), which methylate CHH and CHG cytosines (where H is A, T, or C), and METHYLTRANSFERASE 1 (MET1), which maintains methylation of symmetrical CG contexts. In this review, we discuss the sequence specificities and targeting of each of these pathways. In particular, we highlight recent studies that indicate CMTs preferentially target CWG or CWA/CAW motifs (where W is A or T), and discuss how self-reinforcing feedback loops between DNA methyltransferases and histone modifications characteristic of heterochromatin specify targeting. Finally, the initiating events that lead to gene body methylation are discussed as a model illustrating how interdependent targeting of different silencing pathways can potentiate the establishment of off-target epialleles.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

DNA Methylation and the Evolution of Developmental Complexity in Plants

All land plants so far examined use DNA methylation to silence transposons (TEs). DNA methylation therefore appears to have been co-opted in evolution from an original function in TE management to a developmental function (gene regulation) in both phenotypic plasticity and in normal development. The significance of DNA methylation to the evolution of developmental complexity in plants lies in its role in the management of developmental pathways. As such it is more important in fine tuning the presence, absence, and placement of organs rather than having a central role in the evolution of new organs. Nevertheless, its importance should not be underestimated as it contributes considerably to the range of phenotypic expression and complexity available to plants: the subject of the emerging field of epi-evodevo. Furthermore, changes in DNA methylation can function as a "soft" mutation that may be important in the early stages of major evolutionary novelty.

DNA Methylation Signatures of the Plant Chromomethyltransferases

DNA methylation in plants is traditionally partitioned into CG, CHG and CHH contexts (with H any nucleotide but G). By investigating DNA methylation patterns in trinucleotide contexts in four angiosperm species, we show that such a representation hides spatial and functional partitioning of different methylation pathways and is incomplete. CG methylation (mCG) is largely context-independent whereas, at CHG motifs, there is under-representation of mCCG in pericentric regions of A. thaliana and tomato and throughout the chromosomes of maize and rice. In A. thaliana the biased representation of mCCG in heterochromatin is related to specificities of H3K9 methyltransferase SUVH family members. At CHH motifs there is an over-representation of different variant forms of mCHH that, similarly to mCCG hypomethylation, is partitioned into the pericentric regions of the two dicots but dispersed in the monocot chromosomes. The over-represented mCHH motifs in A. thaliana associate with specific types of transposon including both class I and II elements. At mCHH the contextual bias is due to the involvement of various chromomethyltransferases whereas the context-independent CHH methylation in A. thaliana and tomato is mediated by the RNA-directed DNA methylation process that is most active in the gene-rich euchromatin. This analysis therefore reveals that the sequence context of the methylome of plant genomes is informative about the mechanisms associated with maintenance of methylation and the overlying chromatin structure.

Dense cytosine DNA methylation (mC) in eukaryotes is associated with closed chromatin and gene silencing. In plants it is well known that the sequence context of the mC (either mCG, mCHG or mCHH) provides a clue as to which of several mechanisms is involved but now, based on detailed analyses of the DNA methylome in wild type and mutants of four plant species, we reveal that there is additional information in the mC sequence context. Low mCCG and over-representation of mCAA and mCTA or mCAT in A. thaliana and tomato differentiates regions of the chromosomes near the centromere where methylation is dominated by chromomethyltransferases from the chromosome arms in which mCHH is context-independent and predominantly RNA-directed. Rice and maize have similar sequence context-dependent DNA methylation but the corresponding chromosome domains are not spatially separate as in the dicots. The discovery of the subcomponents of plant methylomes based on sequence context will allow greater resolution in past and future analyses of plant methylomes.

Methylome evolution in plants

Despite major progress in dissecting the molecular pathways that control DNA methylation patterns in plants, little is known about the mechanisms that shape plant methylomes over evolutionary time. Drawing on recent intra- and interspecific epigenomic studies, we show that methylome evolution over long timescales is largely a byproduct of genomic changes. By contrast, methylome evolution over short timescales appears to be driven mainly by spontaneous epimutational events. We argue that novel methods based on analyses of the methylation site frequency spectrum (mSFS) of natural populations can provide deeper insights into the evolutionary forces that act at each timescale.

On the origin and evolutionary consequences of gene body DNA methylation

In plants, CG DNA methylation is prevalent in the transcribed regions of many constitutively expressed genes (gene body methylation; gbM), but the origin and function of gbM remain unknown. Here we report the discovery that Eutrema salsugineum has lost gbM from its genome, to our knowledge the first instance for an angiosperm. Of all known DNA methyltransferases, only CHROMOMETHYLASE 3 (CMT3) is missing from E. salsugineum Identification of an additional angiosperm, Conringia planisiliqua, which independently lost CMT3 and gbM, supports that CMT3 is required for the establishment of gbM. Detailed analyses of gene expression, the histone variant H2A.Z, and various histone modifications in E. salsugineum and in Arabidopsis thaliana epigenetic recombinant inbred lines found no evidence in support of any role for gbM in regulating transcription or affecting the composition and modification of chromatin over evolutionary timescales.