Gene body DNA methylation in plants

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.

Gene duplication and transposition of mobile elements drive evolution of the Rpv3 resistance locus in grapevine

A wild grape haplotype (Rpv3-1) confers resistance to Plasmopara viticola. We mapped the causal factor for resistance to an interval containing a TIR-NB-LRR (TNL) gene pair that originated 1.6-2.6 million years ago by a tandem segmental duplication. Transient coexpression of the TNL pair in Vitis vinifera leaves activated pathogen-induced necrosis and reduced sporulation compared with control leaves. Even though transcripts of the TNL pair from the wild haplotype appear to be partially subject to nonsense-mediated mRNA decay, mature mRNA levels in a homozygous resistant genotype were individually higher than the mRNA trace levels observed for the orthologous single-copy TNL in sensitive genotypes. Allelic expression imbalance in a resistant heterozygote confirmed that cis-acting regulatory variation promotes expression in the wild haplotype. The movement of transposable elements had a major impact on the generation of haplotype diversity, altering the DNA context around similar TNL coding sequences and the GC-content in their proximal 5'-intergenic regions. The wild and domesticated haplotypes also diverged in conserved single-copy intergenic DNA, but the highest divergence was observed in intraspecific and not in interspecific comparisons. In this case, introgression breeding did not transgress the genetic boundaries of the domesticated species, because haplotypes present in modern varieties sometimes predate speciation events between wild and cultivated species.

Epigenetic Footprints of CRISPR/Cas9-Mediated Genome Editing in Plants

CRISPR/Cas9 has been widely applied to various plant species accelerating the pace of plant genome editing and precision breeding in crops. Unintended effects beyond off-target nucleotide mutations are still somewhat unexplored. We investigated the degree and patterns of epigenetic changes after gene editing. We examined changes in DNA methylation in genome-edited promoters of naturally hypermethylated genes (AT1G72350 and AT1G09970) and hypomethylated genes (AT3G17320 and AT5G28770) from Arabidopsis. Transgenic plants were developed via Agrobacterium-mediated floral dip transformation. Homozygous edited lines were selected from segregated T2 plants by an in vitro digestion assay using ribonucleoprotein complex. Bisulfite sequencing comparisons were made between paired groups of edited and non-edited plants to identify changes in DNA methylation of the targeted loci. We found that directed mutagenesis via CRISPR/Cas9 resulted in no unintended morphological or epigenetic alterations. Phenotypes of wild-type, transgenic empty vector, and transgenic edited plants were similar. Epigenetic profiles revealed that methylation patterns of promoter regions flanking target sequences were identical among wild-type, transgenic empty vector, and transgenic edited plants. There was no effect of mutation type on epigenetic status. We also evaluated off-target mutagenesis effects in the edited plants. Potential off-target sites containing up to 4-bp mismatch of each target were sequenced. No off-target mutations were detected in candidate sites. Our results showed that CRISPR/Cas9 did not leave an epigenetic footprint on either the immediate gene-edited DNA and flanking DNA or introduce off-target mutations.

H2A Variants in Arabidopsis: Versatile Regulators of Genome Activity

The eukaryotic nucleosome prevents access to the genome. Convergently evolving histone isoforms, also called histone variants, form diverse families that are enriched over distinct features of plant genomes. Among the diverse families of plant histone variants, H2A.Z exclusively marks genes. Here we review recent research progress on the genome-wide distribution patterns and deposition of H2A.Z in plants as well as its association with histone modifications and roles in plant chromatin regulation. We also discuss some hypotheses that explain the different findings about the roles of H2A.Z in plants.

How Stress Facilitates Phenotypic Innovation Through Epigenetic Diversity

Climate adaptation through phenotypic innovation will become the main challenge for plants during global warming. Plants exhibit a plethora of mechanisms to achieve environmental and developmental plasticity by inducing dynamic alterations of gene regulation and by maximizing natural variation through large population sizes. While successful over long evolutionary time scales, most of these mechanisms lack the short-term adaptive responsiveness that global warming will require. Here, we review our current understanding of the epigenetic regulation of plant genomes, with a focus on stress-response mechanisms and transgenerational inheritance. Field and laboratory-scale experiments on plants exposed to stress have revealed a multitude of temporally controlled, mechanistic strategies integrating both genetic and epigenetic changes on the genome level. We analyze inter- and intra-species population diversity to discuss how methylome differences and transposon activation can be harnessed for short-term adaptive efforts to shape co-evolving traits in response to qualitatively new climate conditions and environmental stress.

Approaches to Whole-Genome Methylome Analysis in Plants

Cytosine methylation as a reversible chromatin mark has been investigated extensively for its influence on gene silencing and the regulation of its dynamic association-disassociation at specific sites within a eukaryotic genome. With the remarkable reductions in cost and time associated with whole-genome DNA sequence analysis, coupled with the high fidelity of bisulfite-treated DNA sequencing, single nucleotide resolution of cytosine methylation repatterning within even very large genomes is increasingly achievable. What remains a challenge is the analysis of genome-wide methylome datasets and, consequently, a clear understanding of the overall influence of methylation repatterning on gene expression or vice versa. Reported data have sometimes been subject to stringent data filtering methods that can serve to skew downstream biological interpretation. These complications derive from methylome analysis procedures that vary widely in method and parameter setting. DNA methylation as a chromatin feature that influences DNA stability can be dynamic and rapidly responsive to environmental change. Consequently, methods to discriminate background "noise" of the system from biological signal in response to specific perturbation is essential in some types of experiments. We describe numerous aspects of whole-genome bisulfite sequence data that must be contemplated as well as the various steps of methylome data analysis which impact the biological interpretation of the final output.

DNA methylation and gene expression profiles characterize epigenetic regulation of lncRNAs in colon adenocarcinoma

The long noncoding RNAs (lncRNAs) are associated with tumorigenesis and progression of cancer. While DNA methylation is a common epigenetic regulator of gene expression, the methylation of lncRNAs was rarely studied. To address this gap, we integrated DNA methylation and RNA-seq data to characterize the landscape of lncRNA methylation in colon adenocarcinoma (COAD). We collected and analyzed the lncRNA expression and methylation data from The Cancer Genome Atlas and Cancer Cell Line Encyclopedia to identify the epigenetically regulated lncRNAs. We further investigated the biological and clinical relevance of the identified lncRNAs via bioinformatics analysis. We identified 20 epigenetically upregulated lncRNAs in COAD, including several well-studied lncRNAs whose methylation regulation were poorly investigated, such as PVT1 and UCA1. We also revealed several novel tumor-associated lncRNAs in COAD, including GATA2-As1 and CYTOR. Next, we explored their biology function using gene set enrichment analysis and competitive endogenous RNA analysis. We characterized the methylation landscape of lncRNA in COAD and identified 20 epigenetically upregulated lncRNAs. Our findings will shed new light on the epigenetic regulation of lncRNA expression by DNA methylation.

Single-base-resolution methylome of giant panda's brain, liver and pancreatic tissue

The giant panda (Ailuropoda melanoleuca) is one of the most endangered mammals, and its conservation has significant ecosystem and cultural service value. Cytosine DNA methylation (5mC) is a stable epigenetic modification to the genome and has multiple functions such as gene regulation. However, DNA methylome of giant panda and its function have not been reported as of yet. Bisulfite sequencing was performed on a 4-day-old male giant panda's brain, liver and pancreatic tissues. We found that the whole genome methylation level was about 0.05% based on reads normalization and mitochondrial DNA was not methylated. Three tissues showed similar methylation tendency in the protein-coding genes of their genomes, but the brain genome had a higher count of methylated genes. We obtained 467 and 1,013 different methylation regions (DMR) genes in brain vs. pancreas and liver, while only 260 DMR genes were obtained in liver vs pancreas. Some lncRNA were also DMR genes, indicating that methylation may affect biological processes by regulating other epigenetic factors. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis indicated that low methylated promoter, high methylated promoter and DMR genes were enriched at some important and tissue-specific items and pathways, like neurogenesis, metabolism and immunity. DNA methylation may drive or maintain tissue specificity and organic functions and it could be a crucial regulating factor for the development of newborn cubs. Our study offers the first insight into giant panda's DNA methylome, laying a foundation for further exploration of the giant panda's epigenetics.

Methylation and gene expression differences between reproductive and sterile bumblebee workers

Phenotypic plasticity is the production of multiple phenotypes from a single genome and is notably observed in social insects. Multiple epigenetic mechanisms have been associated with social insect plasticity, with DNA methylation being explored to the greatest extent. DNA methylation is thought to play a role in caste determination in Apis mellifera, and other social insects, but there is limited knowledge on its role in other bee species. In this study, we analyzed whole genome bisulfite sequencing and RNA-seq data sets from head tissue of reproductive and sterile castes of the eusocial bumblebee Bombus terrestris. We found that genome-wide methylation in B. terrestris is similar to other holometabolous insects and does not differ between reproductive castes. We did, however, find differentially methylated genes between castes, which are enriched for multiple biological processes including reproduction. However, we found no relationship between differential methylation and differential gene expression or differential exon usage between castes. Our results also indicate high intercolony variation in methylation. These findings suggest that methylation is associated with caste differences but may serve an alternate function, other than direct caste determination in this species. This study provides the first insights into the nature of a bumblebee caste-specific methylome as well as its interaction with gene expression and caste-specific alternative splicing, providing greater understanding of the role of methylation in phenotypic plasticity within social bee species. Future experimental work is needed to determine the function of methylation and other epigenetic mechanisms in insects.

Implication for DNA methylation involved in the host transfer of diamondback moth, Plutella xylostella (L.)

DNA methylation exerts extensive impacts on gene expression of various living organisms exposed to environmental variation. However, little is known whether DNA methylation is involved in the host transfer of diamondback moth, Plutella xylostella (L.), a worldwide destructive pest of crucifers. In this study, we found that P. xylostella genome exhibited a relatively low level of DNA methylation on the basis of the CpG O/E prediction and experimental validation. A significant positive linear correlation was observed between the stage-specific expressions of PxDNMT1 and DNA methylation levels (5mC content). Particularly, high levels of DNA methylation and gene expression of PxDNMT1 were observed in eggs and mature females of P. xylostella. After host transfer of P. xylostella from Raphanus sativus to Arabidopsis thaliana, we identified some potential genomic loci that might have changed methylation levels. Using the method of fluorescence-labeled methylation-sensitive amplified polymorphism (F-MSAP), we also found the corresponding genes primarily involved in neural system and signaling. The expressions of six candidate genes were verified by qRT-PCR. One of the genes, Px009600, might be regulated by a DNA methylation-mediated mechanism in response to host transfer. Our study provides evidence for a functional system of DNA methylation in P. xylostella and its possible role in adaptation during host transfer. Further studies should examine methylation as responsive factors to different host plants and environmental cues in insect pests.

EvoChromo: towards a synthesis of chromatin biology and evolution

Over the past few years, interest in chromatin and its evolution has grown. To further advance these interests, we organized a workshop with the support of The Company of Biologists to debate the current state of knowledge regarding the origin and evolution of chromatin. This workshop led to prospective views on the development of a new field of research that we term 'EvoChromo'. In this short Spotlight article, we define the breadth and expected impact of this new area of scientific inquiry on our understanding of both chromatin and evolution.

Extensive profiling in Arabidopsis reveals abundant polysome-associated 24-nt small RNAs including AGO5-dependent pseudogene-derived siRNAs

In a reductionist perspective, plant silencing small (s)RNAs are often classified as mediating nuclear transcriptional gene silencing (TGS) or cytosolic posttranscriptional gene silencing (PTGS). Among the PTGS diagnostics is the association of AGOs and their sRNA cargos with the translation apparatus. In Arabidopsis, this is observed for AGO1 loaded with micro(mi)RNAs and, accordingly, translational-repression (TR) is one layer of plant miRNA action. Using AGO1:miRNA-mediated TR as a paradigm, we explored, with two unrelated polysome-isolation methods, which, among the ten Arabidopsis AGOs and numerous sRNA classes, interact with translation. We found that representatives of all three AGO-clades associate with polysomes, including the TGS-effector AGO4 and stereotypical 24-nt sRNAs that normally mediate TGS of transposons/repeats. Strikingly, approximately half of these annotated 24-nt siRNAs displayed unique matches in coding regions/introns of genes, and in pseudogenes, but not in transposons/repeats commonly found in their vicinity. Protein-coding gene-derived 24-nt sRNAs correlate with gene-body methylation. Those derived from pseudogenes belong to two main clusters defined by their parental-gene expression patterns, and are vastly enriched in AGO5, itself found on polysomes. Based on their tight expression pattern in developing and mature siliques, their biogenesis, and genomic/epigenomic features of their loci-of-origin, we discuss potential roles for these hitherto unknown polysome-enriched, pseudogene-derived siRNAs.

The catalytic core of DEMETER guides active DNA demethylation in Arabidopsis

Flowering plants reproduce via a unique double-fertilization event, producing the zygote and the nutritive endosperm. The genome of the central cell, the precursor of the endosperm, undergoes extensive demethylation prior to fertilization. This epigenetic reconfiguration, directed by the DEMETER (DME) glycosylase at thousands of loci in Arabidopsis, differentiates the epigenetic landscapes of parental genomes and establishes parent of origin-specific expression of many imprinted genes in endosperm essential for seed development. However, how DME is targeted to various locations remains unknown. Here we show that the multidomain DME is organized into 2 functional regions: the C-terminal region, which guides localization and catalysis, and the N-terminal region, which likely recruits chromatin remodelers to facilitate demethylation within heterochromatin.

The Arabidopsis DEMETER (DME) DNA glycosylase demethylates the maternal genome in the central cell prior to fertilization and is essential for seed viability. DME preferentially targets small transposons that flank coding genes, influencing their expression and initiating plant gene imprinting. DME also targets intergenic and heterochromatic regions, but how it is recruited to these differing chromatin landscapes is unknown. The C-terminal half of DME consists of 3 conserved regions required for catalysis in vitro. We show that this catalytic core guides active demethylation at endogenous targets, rescuing dme developmental and genomic hypermethylation phenotypes. However, without the N terminus, heterochromatin demethylation is significantly impeded, and abundant CG-methylated genic sequences are ectopically demethylated. Comparative analysis revealed that the conserved DME N-terminal domains are present only in flowering plants, whereas the domain architecture of DME-like proteins in nonvascular plants mainly resembles the catalytic core, suggesting that it might represent the ancestral form of the 5mC DNA glycosylase found in plant lineages. We propose a bipartite model for DME protein action and suggest that the DME N terminus was acquired late during land plant evolution to improve specificity and facilitate demethylation at heterochromatin targets.

CGT-seq: epigenome-guided de novo assembly of the core genome for divergent populations with large genome

Genetic diversity in plants is remarkably high. Recent whole genome sequencing (WGS) of 67 rice accessions recovered 10,872 novel genes. Comparison of the genetic architecture among divergent populations or between crops and wild relatives is essential for obtaining functional components determining crucial traits. However, many major crops have gigabase-scale genomes, which are not well-suited to WGS. Existing cost-effective sequencing approaches including re-sequencing, exome-sequencing and restriction enzyme-based methods all have difficulty in obtaining long novel genomic sequences from highly divergent population with large genome size. The present study presented a reference-independent core genome targeted sequencing approach, CGT-seq, which employed epigenomic information from both active and repressive epigenetic marks to guide the assembly of the core genome mainly composed of promoter and intragenic regions. This method was relatively easily implemented, and displayed high sensitivity and specificity for capturing the core genome of bread wheat. 95% intragenic and 89% promoter region from wheat were covered by CGT-seq read. We further demonstrated in rice that CGT-seq captured hundreds of novel genes and regulatory sequences from a previously unsequenced ecotype. Together, with specific enrichment and sequencing of regions within and nearby genes, CGT-seq is a time- and resource-effective approach to profiling functionally relevant regions in sequenced and non-sequenced populations with large genomes.

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.

Changes in DNA methylation pattern of apple long-term in vitro shoot culture and acclimatized plants

DNA methylation is a process of epigenetic modification that can alter the functionality of a genome. Using whole-genome bisulfite sequencing, this study quantify the level of DNA methylation in the epigenomes of two diploid apple (Malus x domestica) scion cultivars ('McIntosh' and 'Húsvéti rozmaring') derived from three environmental conditions: in vivo mother plants in an orchard, in vitro culture, and acclimatized in vitro plants. The global DNA methylation levels were not dependent on the source of plant material, and the average level of DNA methylation was 49.77%, 34.65% and 8.77% in CpG, CHG and CHH contexts, respectively. Significant differences in DNA methylation were identified in 586 (specifically 334, 201 and 131 in CpG, CHG and CHH contexts, respectively) out of 45,116 genes, including promoter and coding sequences. These were classified as differentially methylated genes (DMGs). This is a 1.3% difference in the level of DNA methylation of genes in response to a change in the environment. Differential methylation was visualised by MA plots and functional genomic maps were established for biological processes, molecular functions and cellular components. When the DMGs were considered, in vitro tissue culture resulted in the highest level of methylation, but it was lower in acclimatized in vitro plants which was similar to that in the mother tree. Methylation patterns of the two scions differed, indicating cultivar-specific epigenetic regulation of gene expression during adaptation to various environments. After selecting genes that displayed differences larger than ±10% in CpG and CHG contexts, or larger than ±1.35% in the CHH context from among the DMGs, they were annotated in Blast2 GO v5.1.12 for Gene Ontology. DMGs identified as MD07G1113000 (protein transport), MD08G1041600 (extracellular space), MD09G1054800 (phosphatidic acid binding), and MD10G1265800 (not annotated) were methylated in all three contexts in in vitro shoots. These DNA methylation results suggest that epigenetic changes may contribute to the adaptation of apple to environmental changes by modifying the epigenome and thereby gene expression.

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.

Cadmium stress alters cytosine methylation status and expression of a select set of genes in Nicotiana benthamiana

In this paper, we evaluated the genotoxicity of cadmium (Cd) in plants by performing a methylation-sensitive amplification polymorphism (MSAP) on the model plant Nicotiana benthamiana. Among 255 loci examined, 14 genes were found to show altered cytosine methylation patterns in response to Cd stress. Four of those genes (NbMORC3, NbHGSNAT, NbMUT, and NbBG) were selected for further analysis due to their predicted roles in plant development. Cd-induced changes of cytosine methylation status in MSAP fragments of selected genes were confirmed using bisulfite sequencing polymerase chain reaction (BSP). In addition, the expression levels of these genes were found to correlate with cadmium dosage, and a knock-down of these four genes via virus-induced genes silencing (VIGS) led to abnormal development and elevated sensitivity to cadmium stress. Silencing of these four genes resulted in altered cadmium accumulation in different parts of the experimental plants. Our data indicate that cadmium exposure causes dramatic changes in the cytosine methylation status of the plant genome, thus affecting the expression of many genes that are vital for plant growth and are involved in cadmium stress response.

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.

Cytosine Methylation of Isoflavone Synthase Gene in the Genic Region Positively Regulates Its Expression and Isoflavone Biosynthesis in Soybean Seeds

Plants, being sessile organisms, have evolved several dynamic mechanisms of gene regulation. Epigenetic modification especially cytosine methylation and demethylation actively regulates the expression of genes. To understand the role of cytosine methylation during isoflavonoid biosynthesis and accumulation, we performed cytosine methylation analysis in the coding region of two isoforms IFS1 and IFS2 gene, in two contrasting soybean genotypes differing in total isoflavone content (NRC37: high isoflavone; and NRC7: low isoflavone). The results indicated increased 5-mC in both the isoforms in NRC37 (∼20.51% in IFS2 and ∼85% in IFS1) compared with NRC7 (∼7.8% in IFS2 and ∼2.5% in IFS1) genotype, which signifies the positive role of 5-mC in the coding region of the gene leading to enhanced expression. In addition, temporal expression profiling [35 days after flowering (DAF), 45, 55, and 65 DAF] of both the isoforms showed increasing trend of accumulation in both the genotypes with maximum in NRC37 at 65 DAF. To further establish a correlation between methylation and expression of transcripts, we quantified the different isoforms of isoflavone in both the genotypes across all the stages. Therefore, the finding of this study would certainly increase our understanding of epigenetic regulation of isoflavone biosynthetic pathway mediated by the cytosine methylation that would assist molecular breeders to get high-performing soybean genotypes with better isoflavone yield.

Methylome Dynamics of Bovine Gametes and in vivo Early Embryos

DNA methylation undergoes drastic fluctuation during early mammalian embryogenesis. The dynamics of global DNA methylation in bovine embryos, however, have mostly been studied by immunostaining. We adopted the whole genome bisulfite sequencing (WGBS) method to characterize stage-specific genome-wide DNA methylation in bovine sperm, immature oocytes, oocytes matured in vivo and in vitro, as well as in vivo developed single embryos at the 2-, 4-, 8-, and 16-cell stages. We found that the major wave of genome-wide DNA demethylation was complete by the 8-cell stage when de novo methylation became prominent. Sperm and oocytes were differentially methylated in numerous regions (DMRs), which were primarily intergenic, suggesting that these non-coding regions may play important roles in gamete specification. DMRs were also identified between in vivo and in vitro matured oocytes, suggesting environmental effects on epigenetic modifications. In addition, virtually no (less than 1.5%) DNA methylation was found in mitochondrial DNA. Finally, by using RNA-seq data generated from embryos at the same developmental stages, we revealed a weak inverse correlation between gene expression and promoter methylation. This comprehensive analysis provides insight into the critical features of the bovine embryo methylome, and serves as an important reference for embryos produced in vitro, such as by in vitro fertilization and cloning. Lastly, these data can also provide a model for the epigenetic dynamics in human early embryos.

DNA Methylation Analysis of the Citrullus lanatus Response to Cucumber Green Mottle Mosaic Virus Infection by Whole-Genome Bisulfite Sequencing

DNA methylation is an important epigenetic mark associated with plant immunity, butlittle is known about its roles in viral infection of watermelon. We carried out whole-genomebisulfite sequencing of watermelon leaves at 0 h (ck), 48 h, and 25 days post-inoculation withCucumber green mottle mosaic virus (CGMMV). The number of differentially methylated regions(DMRs) increased during CGMMV infection and 2788 DMR-associated genes (DMGs) werescreened out among three libraries. Most DMRs and DMGs were obtained under the CHH context.These DMGs were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG)pathways of secondary biosynthesis and metabolism, plant-pathogen interactions, Toll-likereceptor signaling, and ABC transporters. Additionally, DMGs encoding PR1a, CaMs, calciumbindingprotein, RIN4, BAK1, WRKYs, RBOHs, STKs, and RLPs/RLKs were involved in thewatermelon-CGMMV interaction and signaling. The association between DNA methylation andgene expression was analyzed by RNA-seq and no clear relationship was detected. Moreover,downregulation of genes in the RdDM pathway suggested the reduced RdDM-directed CHHmethylation plays an important role in antiviral defense in watermelon. Our findings providegenome-wide DNA methylation profiles of watermelon and will aid in revealing the molecularmechanism in response to CGMMV infection at the methylation level.

Combining a Simple Method for DNA/RNA/Protein Co-Purification and Arabidopsis Protoplast Assay to Facilitate Viroid Research

Plant–viroid interactions represent a valuable model for delineating structure–function relationships of noncoding RNAs. For various functional studies, it is desirable to minimize sample variations by using DNA, RNA, and proteins co-purified from the same samples. Currently, most of the co-purification protocols rely on TRI Reagent (Trizol as a common representative) and require protein precipitation and dissolving steps, which render difficulties in experimental handling and high-throughput analyses. Here, we established a simple and robust method to minimize the precipitation steps and yield ready-to-use RNA and protein in solutions. This method can be applied to samples in small quantities, such as protoplasts. Given the ease and the robustness of this new method, it will have broad applications in virology and other disciplines in molecular biology.

Analysis of DNA Methylation Patterns Associated with In Vitro Propagated Globe Artichoke Plants Using an EpiRADseq-Based Approach

Globe artichoke represents one of the main horticultural species of the Mediterranean basin, and 'Spinoso sardo' is the most widespread and economically relevant varietal type in Sardinia, Italy. In the last decades, in vitro culture of meristematic apices has increased the frequency of aberrant plants in open-field production. These off-type phenotypes showed highly pinnate-parted leaves and late inflorescence budding, and emerged from some branches of the true-to-type 'Spinoso sardo' plants. This phenomenon cannot be foreseen and is reversible through generations, suggesting the occurrence of epigenetic alterations. Here, we report an exploratory study on DNA methylation patterns in off-type/true-to-type globe artichoke plants, using a modified EpiRADseq technology, which allowed the identification of 2,897 differentially methylated loci (DML): 1,998 in CG, 458 in CHH, and 441 in CHG methylation contexts of which 720, 88, and 152, respectively, were in coding regions. Most of them appeared involved in primary metabolic processes, mostly linked to photosynthesis, regulation of flower development, and regulation of reproductive processes, coherently with the observed phenotype. Differences in the methylation status of some candidate genes were integrated with transcriptional analysis to test whether these two regulation levels might interplay in the emergence and spread of the 'Spinoso sardo' non-conventional phenotype.

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.

Universality of the DNA methylation codes in Eucaryotes

Genetics and epigenetics are tightly linked heritable information classes. Question arises if epigenetics provides just a set of environment dependent instructions, or whether it is integral part of an inheritance system. We argued that in the latter case the epigenetic code should share the universality quality of the genetic code. We focused on DNA methylation. Since availability of DNA methylation data is biased towards model organisms we developed a method that uses kernel density estimations of CpG observed/expected ratios to infer DNA methylation types in any genome. We show here that our method allows for robust prediction of mosaic and full gene body methylation with a PPV of 1 and 0.87, respectively. We used this prediction to complement experimental data, and applied hierarchical clustering to identify methylation types in ~150 eucaryotic species covering different body plans, reproduction types and living conditions. Our analysis indicates that there are only four gene body methylation types. These types do not follow phylogeny (i.e. phylogenetically distant clades can have identical methylation types) but they are consistent within clades. We conclude that the gene body DNA methylation codes have universality similar to the universality of the genetic code and should consequently be considered as part of the inheritance system.

Single-cell expression noise and gene-body methylation in Arabidopsis thaliana

Gene-body methylation (gbM) refers to an increased level of methylated cytosines specifically in a CG sequence context within genes. gbM is found in plant genes with intermediate expression level, which evolve slowly, and is often broadly conserved across millions of years of evolution. Intriguingly however, some plants lack gbM, and thus it remains unclear whether gbM has a function. In animals, there is support for a role of gbM in reducing erroneous transcription and transcription noise, but so far most studies in plants have tested for an effect of gbM on expression level, not noise. Here, we therefore tested whether gbM was associated with reduced expression noise in Arabidopsis thaliana, using single-cell transcriptome sequencing data from root quiescent centre cells. We find that gbM genes have lower expression noise levels than unmethylated genes. However, an analysis of covariance revealed that, if other genomic features are taken into account, this association disappears. Nonetheless, gbM genes were more consistently expressed across single-cell samples, supporting previous inference that gbM genes are constitutively expressed. Finally, we observed that fewer RNAseq reads map to introns of gbM genes than to introns of unmethylated genes, which indicates that gbM might be involved in reducing erroneous transcription by reducing intron retention.

Opportunities and limitations of reduced representation bisulfite sequencing in plant ecological epigenomics

Contents Summary 738 I. Introduction 738 II. RRBS loci as genome-wide epigenetic markers 739 III. Exploiting functional annotation of RRBS loci 739 IV. Limitations of RRBS methods for nonmodel species 740 V. Maximising the impact of RRBS in plants 741 VI. Conclusions 741 Acknowledgements 741 SUMMARY: Investigating the features and implications of epigenetic mechanisms across the breadth of organisms and ecosystems is important for understanding the ecological relevance of epigenetics. Several cost-effective reduced representation bisulfite sequencing approaches (RRBS) have been recently developed and applied to different organisms that lack a well annotated reference genome. These new approaches improve the assessment of epigenetic diversity in ecological settings and may provide functional insights. We assess here the opportunities and limitations of RRBS in nonmodel plant species. Well thought out experimental designs that include complementary gene expression studies, and the improvement of genomics resources for the target group, promise to maximize the effect of future RRBS studies.

Variation in ICE1 Methylation Primarily Determines Phenotypic Variation in Freezing Tolerance in Arabidopsis thaliana

Cold stress is a major abiotic factor plants face during their life cycle. Although plants often exhibit phenotypic variation in cold tolerance, the underlying mechanism remains poorly understood. In the present study, the 50% lethal temperature (LT50) values of 37 Arabidopsis thaliana accessions at latitudes from 15° to 58° ranged from -13.2°C to -4.9°C and were closely correlated with the cold climates of the collection sites. According to a methylation analysis of all C-repeat (CRT)-binding factor (CBF) pathway genes, the coding and promoter regions of AtICE1, a regulator of CBF genes, exhibited the greatest variability in methylation levels among the accessions and included 5-122 methylated cytosine residues. In contrast, unmethylated or only slightly methylated genes in the CBF pathway showed little variation among the accessions. According to a gene expression analysis of four selected A. thaliana populations with distinct methylation patterns, except for the down-regulated gene AtCBF2, the expression levels of all members of the CBF pathway were negatively correlated with AtICE1 gene methylation levels. Treatment of the four A. thaliana populations with the DNA methylation inhibitory reagent 5-azacytidine resulted in a 30.0-78.3% enhancement of freezing tolerance and decreases in LT50 values of approximately 1.9-3.6°C. Similar effects were observed in drm2 mutants, including 30.0-48.3% increases in freezing tolerance and decreases in LT50 values of approximately 0.7-3.4°C. Thus, the AtICE1 methylation-regulated transcription of CBF pathway genes is responsible for the phenotypic variation in the freezing tolerance observed in A. thaliana.

Role of gene body methylation in acclimatization and adaptation in a basal metazoan

Gene body methylation (GBM) has been hypothesized to modulate responses to environmental change, including transgenerational plasticity, but the evidence thus far has been lacking. Here we show that coral fragments reciprocally transplanted between two distant reefs respond predominantly by increase or decrease in genome-wide GBM disparity: The range of methylation levels between lowly and highly methylated genes becomes either wider or narrower. Remarkably, at a broad functional level this simple adjustment correlated very well with gene expression change, reflecting a shifting balance between expressions of environmentally responsive and housekeeping genes. In our experiment, corals in a lower-quality habitat up-regulated genes involved in environmental responses, while corals in a higher-quality habitat invested more in housekeeping genes. Transplanted fragments showing closer GBM match to local corals attained higher fitness characteristics, which supports GBM's role in acclimatization. Fixed differences in GBM between populations did not align with plastic GBM changes and were mostly observed in genes with elevated F _ST, which suggests that they arose predominantly through genetic divergence. However, we cannot completely rule out transgenerational inheritance of acquired GBM states.

Parental experience modifies the Mimulus methylome

BACKGROUND

Transgenerational plasticity occurs when the environmental experience of an organism modifies the growth and development of its progeny. Leaf damage in Mimulus guttatus exhibits transgenerational plasticity mediated through differential expression of hundreds of genes. The epigenetic mechanisms that facilitate this response have yet to be described.

RESULTS

We performed whole genome bisulfite sequencing in the progeny of genetically identical damaged and control plants and developed a pipeline to compare differences in the mean and variance of methylation between treatment groups. We find that parental damage increases the variability of CG and CHG methylation among progeny, but does not alter the overall mean methylation. Instead it has positive effects in some regions and negative in others. We find 3,396 CHH, 203 CG, and 54 CHG Differentially Methylated Regions (DMRs) ranging from tens to thousands of base pairs scattered across the genome. CHG and CHH DMRs tended to overlap with transposable elements. CG DMRs tended to overlap with gene coding regions, many of which were previously found to be differentially expressed.

CONCLUSIONS

Genome-wide increases in methylome variation suggest that parental conditions can increase epigenetic diversity in response to stress. Additionally, the potential association between CG DMRs and differentially expressed genes supports the hypothesis that differential methylation is a mechanistic component of transgenerational plasticity in M. guttatus.

The methylome of the marbled crayfish links gene body methylation to stable expression of poorly accessible genes

Background

The parthenogenetic marbled crayfish (Procambarus virginalis) is a novel species that has rapidly invaded and colonized various different habitats. Adaptation to different environments appears to be independent of the selection of genetic variants, but epigenetic programming of the marbled crayfish genome remains to be understood.

Results

Here, we provide a comprehensive analysis of DNA methylation in marbled crayfish. Whole-genome bisulfite sequencing of multiple replicates and different tissues revealed a methylation pattern that is characterized by gene body methylation of housekeeping genes. Interestingly, this pattern was largely tissue invariant, suggesting a function that is unrelated to cell fate specification. Indeed, integrative analysis of DNA methylation, chromatin accessibility and mRNA expression patterns revealed that gene body methylation correlated with limited chromatin accessibility and stable gene expression, while low-methylated genes often resided in chromatin with higher accessibility and showed increased expression variation. Interestingly, marbled crayfish also showed reduced gene body methylation and higher gene expression variability when compared with their noninvasive mother species, Procambarus fallax.

Conclusions

Our results provide novel insights into invertebrate gene body methylation and its potential role in adaptive gene regulation.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0229-6) contains supplementary material, which is available to authorized users.

The relationship between transgenerational acquired resistance and global DNA methylation in Arabidopsis

Progeny of heavily diseased plants develop transgenerational acquired resistance (TAR). In Arabidopsis, TAR can be transmitted over one stress-free generation. Although DNA methylation has been implicated in the regulation of TAR, the relationship between TAR and global DNA methylation remains unknown. Here, we characterised the methylome of TAR-expressing Arabidopsis at different generations after disease exposure. Global clustering of cytosine methylation revealed TAR-related patterns in the F3 generation, but not in the F1 generation. The majority of differentially methylated positions (DMPs) occurred at CG context in gene bodies. TAR in F3 progeny after one initial generation of disease, followed by two stress-free generations, was lower than TAR in F3 progeny after three successive generations of disease. This difference in TAR effectiveness was proportional to the intensity of differential methylation at a sub-set of cytosine positions. Comparison of TAR-related DMPs with previously characterised cytosine methylation in mutation accumulation lines revealed that ancestral disease stress preferentially acts on methylation-labile cytosine positions, but also extends to methylation-stable positions. Thus, the TAR-related impact of ancestral disease extends beyond stochastic variation in DNA methylation. Our study has shown that the Arabidopsis epigenome responds globally to disease in previous generations and we discuss its contribution to TAR.

Genome-Wide DNA Methylation Comparison between Brassica napus Genic Male Sterile Line and Restorer Line

DNA methylation is an essential epigenetic modification that dynamically regulates gene expression during plant development. However, few studies have determined the DNA methylation profiles of male-sterile rapeseed. Here, we conducted a global comparison of DNA methylation patterns between the rapeseed genic male sterile line 7365A and its near-isogenic fertile line 7365B by whole-genome bisulfite sequencing (WGBS). Profiling of the genome-wide DNA methylation showed that the methylation level in floral buds was lower than that in leaves and roots. Besides, a total of 410 differentially methylated region-associated genes (DMGs) were identified in 7365A relative to 7365B. Traditional bisulfite sequencing polymerase chain reaction (PCR) was performed to validate the WGBS data. Eleven DMGs were found to be involved in anther and pollen development, which were analyzed by quantitative PCR. In particular, Bnams4 was hypo-methylated in 7365A, and its expression was up-regulated, which might affect other DMGs and thus control the male sterility. This study provided genome-wide DNA methylation profiles of floral buds and important clues for revealing the molecular mechanism of genic male sterility in rapeseed.

DNA methylation footprints during soybean domestication and improvement

BACKGROUND

In addition to genetic variation, epigenetic variation plays an important role in determining various biological processes. The importance of natural genetic variation to crop domestication and improvement has been widely investigated. However, the contribution of epigenetic variation in crop domestication at population level has rarely been explored.

RESULTS

To understand the impact of epigenetics on crop domestication, we investigate the variation of DNA methylation during soybean domestication and improvement by whole-genome bisulfite sequencing of 45 soybean accessions, including wild soybeans, landraces, and cultivars. Through methylomic analysis, we identify 5412 differentially methylated regions (DMRs). These DMRs exhibit characters distinct from those of genetically selected regions. In particular, they have significantly higher genetic diversity. Association analyses suggest only 22.54% of DMRs can be explained by local genetic variations. Intriguingly, genes in the DMRs that are not associated with any genetic variation are enriched in carbohydrate metabolism pathways.

CONCLUSIONS

This study provides a valuable map of DNA methylation across diverse accessions and dissects the relationship between DNA methylation variation and genetic variation during soybean domestication, thus expanding our understanding of soybean domestication and improvement.

Expression and DNA methylation of SERK, BBM, LEC2 and WUS genes in in vitro cultures of Boesenbergia rotunda (L.) Mansf

The process of somatic embryogenesis and plant regeneration involve changes in gene expression and have been associated with changes in DNA methylation. Here, we report the expression and DNA methylation patterns of SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK), BABY BOOM (BBM), LEAFY COTYLEDON 2 (LEC2) and WUSCHEL (WUS) in meristematic block of newly emerged shoots from rhizome, embryogenic and non-embryogenic calli, prolonged cell suspension culture, ex vitro leaf, and in vitro leaf of regenerated plants of Boesenbergia rotunda. Among all seven samples, based on qRT-PCR, the highest level of expression of SERK, BBM and LEC2 was in embryogenic callus, while WUS was most highly expressed in meristematic block tissue followed by embryogenic callus. Relatively lower expression was observed in cell suspension culture and watery callus for SERK, LEC2 and WUS and in in vitro leaf for BBM. For gene specific methylation determined by bisulfite sequencing data, embryogenic callus samples had the lowest levels of DNA methylation at CG, CHG and CHH contexts of SERK, LEC2 and WUS. We observed negative correlation between DNA methylation at the CG and CHG contexts and the expression levels of SERK, BBM, LEC2 and WUS. Based on our results, we suggest that relatively higher expression and lower level of DNA methylation of SERK, BBM, LEC2 and WUS are associated with somatic embryogenesis and plant regeneration in B. rotunda.

Genetic and epigenetic divergence of duplicate genes in two legume species

Soybean (Glycine max) and common bean (Phaseolus vulgaris) share a polyploidy event ~59 MYA, followed by a Glycine-specific whole genome duplication (WGD) ~8-13 MYA. Duplicated genes were classified into five categories: singletons, dispersed, proximal, tandem, or WGD/segmental and found strong correlations between gene category and functional annotation. Photosynthesis and transcriptional regulation-related Gene Ontology terms were significantly over-represented in singletons and WGD genes, respectively, aligning with the gene balance hypothesis. We found that the divergence of gene expression and DNA methylation between WGD-derived paralogs increased with age and that WGD genes, initially retained via dosage constraints, subsequently underwent expression divergence, associated with other factors such as DNA methylation. Genes derived from different modes of duplication differed in breadth, level, and specificity of expression in both species. Orthologous genes and ungrouped genes (genes not in an ortholog group) differed in expression patterns. The protein divergence rates of WGD paralog pairs containing an ungrouped gene were higher than those for which both copies had orthologs. We propose that many ungrouped genes are derived from divergent and redundant gene copies, concordant with the neofunctionalization hypothesis. Tandemly duplicated genes were distinct from WGD-derived genes, indicating that mode of duplication contributes to the evolutionary fate of duplicated genes.

The transcriptional landscape of polyploid wheat

The coordinated expression of highly related homoeologous genes in polyploid species underlies the phenotypes of many of the world's major crops. Here we combine extensive gene expression datasets to produce a comprehensive, genome-wide analysis of homoeolog expression patterns in hexaploid bread wheat. Bias in homoeolog expression varies between tissues, with ~30% of wheat homoeologs showing nonbalanced expression. We found expression asymmetries along wheat chromosomes, with homoeologs showing the largest inter-tissue, inter-cultivar, and coding sequence variation, most often located in high-recombination distal ends of chromosomes. These transcriptionally dynamic genes potentially represent the first steps toward neo- or subfunctionalization of wheat homoeologs. Coexpression networks reveal extensive coordination of homoeologs throughout development and, alongside a detailed expression atlas, provide a framework to target candidate genes underpinning agronomic traits in wheat.

Heritable Epigenomic Changes to the Maize Methylome Resulting from Tissue Culture

DNA methylation can contribute to the maintenance of genome integrity and regulation of gene expression. In most situations, DNA methylation patterns are inherited quite stably. However, changes in DNA methylation can occur at some loci as a result of tissue culture resulting in somaclonal variation. To investigate heritable epigenetic changes as a consequence of tissue culture, a sequence-capture bisulfite sequencing approach was implemented to monitor context-specific DNA methylation patterns in ∼15 Mb of the maize genome for a population of plants that had been regenerated from tissue culture. Plants that have been regenerated from tissue culture exhibit gains and losses of DNA methylation at a subset of genomic regions. There was evidence for a high rate of homozygous changes to DNA methylation levels that occur consistently in multiple independent tissue culture lines, suggesting that some loci are either targeted or hotspots for epigenetic variation. The consistent changes inherited following tissue culture include both gains and losses of DNA methylation and can affect CG, CHG, or both contexts within a region. Only a subset of the tissue culture changes observed in callus plants are observed in the primary regenerants, but the majority of DNA methylation changes present in primary regenerants are passed onto offspring. This study provides insights into the susceptibility of some loci and potential mechanisms that could contribute to altered DNA methylation and epigenetic state that occur during tissue culture in plant species.

Locus-specific control of the de novo DNA methylation pathway in Arabidopsis by the CLASSY family

DNA methylation is essential for gene regulation, transposon silencing, and imprinting. Although the generation of specific DNA methylation patterns is critical for these processes, how methylation is regulated at individual loci remains unclear. Here we show that a family of four putative chromatin remodeling factors, CLASSY (CLSY) 1–4, are required for both locus-specific and global regulation of DNA methylation in Arabidopsis. Mechanistically, these factors act in connection with RNA polymerase-IV (Pol-IV) to control the production of 24-nucleotide small interfering RNAs (24nt-siRNAs), which guide DNA methylation. Individually, the CLSYs regulate Pol-IV-chromatin association and 24nt-siRNA production at thousands of distinct loci, and together, they regulate essentially all 24nt-siRNAs. Depending on the CLSYs involved, this regulation relies on different repressive chromatin modifications to facilitate locus-specific control of DNA methylation. Given the conservation between methylation systems in plants and mammals, analogous pathways likely operate in a broad range of organisms.

FACT complex is required for DNA demethylation at heterochromatin during reproduction in Arabidopsis

The DEMETER (DME) DNA glycosylase catalyzes genome-wide DNA demethylation and is required for endosperm genomic imprinting and embryo viability. Targets of DME-mediated DNA demethylation reside in small, euchromatic, AT-rich transposons and at the boundaries of large transposons, but how DME interacts with these diverse chromatin states is unknown. The STRUCTURE SPECIFIC RECOGNITION PROTEIN 1 (SSRP1) subunit of the chromatin remodeler FACT (facilitates chromatin transactions), was previously shown to be involved in the DME-dependent regulation of genomic imprinting in Arabidopsis endosperm. Therefore, to investigate the interaction between DME and chromatin, we focused on the activity of the two FACT subunits, SSRP1 and SUPPRESSOR of TY16 (SPT16), during reproduction in Arabidopsis We found that FACT colocalizes with nuclear DME in vivo, and that DME has two classes of target sites, the first being euchromatic and accessible to DME, but the second, representing over half of DME targets, requiring the action of FACT for DME-mediated DNA demethylation genome-wide. Our results show that the FACT-dependent DME targets are GC-rich heterochromatin domains with high nucleosome occupancy enriched with H3K9me2 and H3K27me1. Further, we demonstrate that heterochromatin-associated linker histone H1 specifically mediates the requirement for FACT at a subset of DME-target loci. Overall, our results demonstrate that FACT is required for DME targeting by facilitating its access to heterochromatin.

Chilling-induced DNA Demethylation is associated with the cold tolerance of Hevea brasiliensis

BACKGROUND

Low temperature influences the development and latex production of rubber trees (Hevea brasiliensis) when extension to suboptimal high-latitude areas. The successful extension of Hevea brasiliensis cultivation to high-latitude areas has long believed to benefit from the breeding of cold-tolerant cultivars. A puzzling incongruity is the variation in cold tolerance among the cultivated clones despite their similar genetic make-up.

RESULTS

To investigate this, we first transferred cultivar Reyan 7-33-97 to short-term cold treatment, and showed that cold-related genes (such as HbICE1 and HbCBF2), cold-responsive (COR) genes, and DNA-methylation related genes (such as HbMET1) were induced by cold treatment. Furthermore, long-term cold treatment not only elevated the transcriptional activities of the HbICE1, HbCBF2, and HbMET, but also induced DNA demethylation of their promoters. Cold treatment increased the transcriptional activities of demethylation-related genes such as the HbDME, HbROS, and HbDML genes, but did not alter the promoter methylation status. Furthermore, the HbICE1 and HbMET promoters showed hypomethylation status in samples collected at the end of winter from 12 different cultivars grown in four geographical locations, but switched to hypermethylation status at the end of summer. Expression of COR was correlated with the low temperature. Given that little genetic diversity exists in the HbICE1 and HbMET promoters among different cultivars, the DNA demethylation induced by cold was highly correlated with low temperature, but not with the genetic backgrounds of cultivars.

CONCLUSION

Cold-induced epigenetic modification might play an important role in cold tolerance of H. brasiliensis.

Epigenetic perspectives on the evolution and domestication of polyploid plant and crops

Polyploidy or whole genome duplication (WGD) is a prominent feature for genome evolution of some animals and all flowering plants, including many important crops such as wheat, cotton, and canola. In autopolyploids, genome duplication often perturbs dosage regulation on biological networks. In allopolyploids, interspecific hybridization could induce genetic and epigenetic changes, the effects of which could be amplified by genome doubling (ploidy changes). Albeit the importance of genetic changes, some epigenetic changes can be stabilized and transmitted as epialleles into the progeny, which are subject to natural selection, adaptation, and domestication. Here we review recent advances for general and specific roles of epigenetic changes in the evolution of flowering plants and domestication of agricultural crops.

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.