DNA Methylation and Histone H1 Jointly Repress Transposable Elements and Aberrant Intragenic Transcripts.

DNA methylation and histone H1 mediate transcriptional silencing of genes and transposable elements, but how they interact is unclear. In plants and animals with mosaic genomic methylation, functionally mysterious methylation is also common within constitutively active housekeeping genes. Here, we show that H1 is enriched in methylated sequences, including genes, of Arabidopsis thaliana, yet this enrichment is independent of DNA methylation. Loss of H1 disperses heterochromatin, globally alters nucleosome organization, and activates H1-bound genes, but only weakly de-represses transposable elements. However, H1 loss strongly activates transposable elements hypomethylated through mutation of DNA methyltransferase MET1. Hypomethylation of genes also activates antisense transcription, which is modestly enhanced by H1 loss. Our results demonstrate that H1 and DNA methylation jointly maintain transcriptional homeostasis by silencing transposable elements and aberrant intragenic transcripts. Such functionality plausibly explains why DNA methylation, a well-known mutagen, has been maintained within coding sequences of crucial plant and animal genes.

Genome wide inherited modifications of the tomato epigenome by trans-activated bacterial CG methyltransferase.

BACKGROUND: Epigenetic variation is mediated by epigenetic marks such as DNA methylation occurring in all cytosine contexts in plants. CG methylation plays a critical role in silencing transposable elements and regulating gene expression. The establishment of CG methylation occurs via the RNA-directed DNA methylation pathway and CG methylation maintenance relies on METHYLTRANSFERASE1, the homologue of the mammalian DNMT1. PURPOSE: Here, we examined the capacity to stably alter the tomato genome methylome by a bacterial CG-specific M.SssI methyltransferase expressed through the LhG4/pOP transactivation system. RESULTS: Methylome analysis of M.SssI expressing plants revealed that their euchromatic genome regions are specifically hypermethylated in the CG context, and so are most of their genes. However, changes in gene expression were observed only with a set of genes exhibiting a greater susceptibility to CG hypermethylation near their transcription start site. Unlike gene rich genomic regions, our analysis revealed that heterochromatic regions are slightly hypomethylated at CGs only. Notably, some M.SssI-induced hypermethylation persisted even without the methylase or transgenes, indicating inheritable epigenetic modification. CONCLUSION: Collectively our findings suggest that heterologous expression of M.SssI can create new inherited epigenetic variations and changes in the methylation profiles on a genome wide scale. This open avenues for the conception of epigenetic recombinant inbred line populations with the potential to unveil agriculturally valuable tomato epialleles.

BSXplorer: analytical framework for exploratory analysis of BS-seq data.

BACKGROUND: Bisulfite sequencing detects and quantifies DNA methylation patterns, contributing to our understanding of gene expression regulation, genome stability maintenance, conservation of epigenetic mechanisms across divergent taxa, epigenetic inheritance and, eventually, phenotypic variation. Graphical representation of methylation data is crucial in exploring epigenetic regulation on a genome-wide scale in both plants and animals. This is especially relevant for non-model organisms with poorly annotated genomes and/or organisms where genome sequences are not yet assembled on chromosome level. Despite being a technology of choice to profile DNA methylation for many years now there are surprisingly few lightweight and robust standalone tools available for efficient graphical analysis of data in non-model systems. This significantly limits evolutionary studies and agrigenomics research. BSXplorer is a tool specifically developed to fill this gap and assist researchers in explorative data analysis and in visualising and interpreting bisulfite sequencing data more easily. RESULTS: BSXplorer provides in-depth graphical analysis of sequencing data encompassing (a) profiling of methylation levels in metagenes or in user-defined regions using line plots and heatmaps, generation of summary statistics charts, (b) enabling comparative analyses of methylation patterns across experimental samples, methylation contexts and species, and (c) identification of modules sharing similar methylation signatures at functional genomic elements. The tool processes methylation data quickly and offers API and CLI capabilities, along with the ability to create high-quality figures suitable for publication. CONCLUSIONS: BSXplorer facilitates efficient methylation data mining, contrasting and visualization, making it an easy-to-use package that is highly useful for epigenetic research.

H3K4/H3K9me3 Bivalent Chromatin Domains Targeted by Lineage-Specific DNA Methylation Pauses Adipocyte Differentiation.

Bivalent H3K4me3 and H3K27me3 chromatin domains in embryonic stem cells keep active developmental regulatory genes expressed at very low levels and poised for activation. Here, we show an alternative and previously unknown bivalent modified histone signature in lineage-committed mesenchymal stem cells and preadipocytes that pairs H3K4me3 with H3K9me3 to maintain adipogenic master regulatory genes (Cebpa and Pparg) expressed at low levels yet poised for activation when differentiation is required. We show lineage-specific gene-body DNA methylation recruits H3K9 methyltransferase SETDB1, which methylates H3K9 immediately downstream of transcription start sites marked with H3K4me3 to establish the bivalent domain. At the Cebpa locus, this prevents transcription factor C/EBPß binding, histone acetylation, and further H3K4me3 deposition and is associated with pausing of RNA polymerase II, which limits Cebpa gene expression and adipogenesis.

Gene body methylation safeguards ribosomal DNA transcription by preventing PHF6-mediated enrichment of repressive histone mark H4K20me3.

DNA methylation shows complex correlations with gene expression, and the role of promoter hypermethylation in repressing gene transcription has been well addressed. Emerging evidence indicates that gene body methylation promotes transcription; however, the underlying mechanisms remain to be further investigated. Here, using methylated DNA immunoprecipitation sequencing (MeDIP-seq), bisulfite genomic sequencing, and immunofluorescent labeling, we show that gene body methylation is indeed positively correlated with rRNA gene (rDNA) transcription. Mechanistically, gene body methylation is largely maintained by DNA methyltransferase 1 (DNMT1), deficiency or downregulation of which during myoblast differentiation or nutrient deprivation results in decreased gene body methylation levels, leading to increased gene body occupancy of plant homeodomain (PHD) finger protein 6 (PHF6). PHF6 binds to hypomethylated rDNA gene bodies where it recruits histone methyltransferase SUV4-20H2 to establish the repressive histone modification, H4K20me3, ultimately inhibiting rDNA transcription. These findings demonstrate that DNMT1-mediated gene body methylation safeguards rDNA transcription by preventing enrichment of repressive histone modifications, suggesting that gene body methylation serves to maintain gene expression in response to developmental and/or environmental stresses.

Changes in gene body methylation do not correlate with changes in gene expression in Anthozoa or Hexapoda.

BACKGROUND: As human activity alters the planet, there is a pressing need to understand how organisms adapt to environmental change. Of growing interest in this area is the role of epigenetic modifications, such as DNA methylation, in tailoring gene expression to fit novel conditions. Here, we reanalyzed nine invertebrate (Anthozoa and Hexapoda) datasets to validate a key prediction of this hypothesis: changes in DNA methylation in response to some condition correlate with changes in gene expression. RESULTS: In accord with previous observations, baseline levels of gene body methylation (GBM) positively correlated with transcription, and negatively correlated with transcriptional variation between conditions. Correlations between changes in GBM and transcription, however, were negligible. There was also no consistent negative correlation between methylation and transcription at the level of gene body methylation class (either highly- or lowly-methylated), anticipated under the previously described "seesaw hypothesis". CONCLUSION: Our results do not support the direct involvement of GBM in regulating dynamic transcriptional responses in invertebrates. If changes in DNA methylation regulate invertebrate transcription, the mechanism must involve additional factors or regulatory influences.

Phylogenetic Shifts in Gene Body Methylation Correlate with Gene Expression and Reflect Trait Conservation.

A subset of genes in plant genomes are labeled with DNA methylation specifically at CG residues. These genes, known as gene-body methylated (gbM), have a number of associated characteristics. They tend to have longer sequences, to be enriched for intermediate expression levels, and to be associated with slower rates of molecular evolution. Most importantly, gbM genes tend to maintain their level of DNA methylation between species, suggesting that this trait is under evolutionary constraint. Given the degree of conservation in gbM, we still know surprisingly little about its function in plant genomes or whether gbM is itself a target of selection. To address these questions, we surveyed DNA methylation across eight grass (Poaceae) species that span a gradient of genome sizes. We first established that genome size correlates with genome-wide DNA methylation levels, but less so for genic levels. We then leveraged genomic data to identify a set of 2,982 putative orthologs among the eight species and examined shifts of methylation status for each ortholog in a phylogenetic context. A total of 55% of orthologs exhibited a shift in gbM, but these shifts occurred predominantly on terminal branches, indicating that shifts in gbM are rarely conveyed over time. Finally, we found that the degree of conservation of gbM across species is associated with increased gene length, reduced rates of molecular evolution, and increased gene expression level, but reduced gene expression variation across species. Overall, these observations suggest a basis for evolutionary pressure to maintain gbM status over evolutionary time.

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

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

Genetic, epigenetic and genomic effects on variation of gene expression among grape varieties.

The transcriptional regulatory structure of plant genomes is still relatively unexplored, and little is known about factors that influence expression variation in plants. We used a genetic system consisting of 10 heterozygous grape varieties with high consanguinity and high haplotypic diversity to: (i) identify regions of haplotype sharing through whole-genome resequencing and single-nucleotide polymorphism (SNP) genotyping; (ii) analyse gene expression through RNA-seq in four stages of berry development; and (iii) associate gene expression variation with genetic and epigenetic properties. We found that haplotype sharing in and around genes was positively correlated with similarity in expression and was negatively correlated with the fraction of differentially expressed genes. Genetic and epigenetic properties of the gene and the surrounding region showed significant effects on the extent of expression variation, with negative associations for the level of gene body methylation and mean expression level, and with positive associations for nucleotide diversity, structural diversity and ratio of non-synonymous to synonymous nucleotide diversity. We also observed a spatial dependency of covariation of gene expression among varieties. These results highlight relevant roles for cis-acting factors, selective constraints and epigenetic features of the gene, and the regional context in which the gene is located, in the determination of expression variation. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA385116; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA392287; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA373967 (released upon publication); https://www.ncbi.nlm.nih.gov/bioproject/PRJNA490160 (released upon publication); https://www.ncbi.nlm.nih.gov/bioproject/PRJNA265039; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA265040.

Genomic analysis reveals major determinants of cis-regulatory variation in Capsella grandiflora.

Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants, and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation (i) are under weaker purifying selection and (ii) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby transposable elements (TEs) and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage sensitivity of body-methylated genes. Given the extensive conservation of gbM in flowering plants, this suggests that gbM could be an important predictor of cis-regulatory variation in a wide range of plant species.

Roles of ASYMMETRIC LEAVES2 (AS2) and Nucleolar Proteins in the Adaxial-Abaxial Polarity Specification at the Perinucleolar Region in Arabidopsis.

Leaves of Arabidopsis develop from a shoot apical meristem grow along three (proximal-distal, adaxial-abaxial, and medial-lateral) axes and form a flat symmetric architecture. ASYMMETRIC LEAVES2 (AS2), a key regulator for leaf adaxial-abaxial partitioning, encodes a plant-specific nuclear protein and directly represses the abaxial-determining gene ETTIN/AUXIN RESPONSE FACTOR3 (ETT/ARF3). How AS2 could act as a critical regulator, however, has yet to be demonstrated, although it might play an epigenetic role. Here, we summarize the current understandings of the genetic, molecular, and cellular functions of AS2. A characteristic genetic feature of AS2 is the presence of a number of (about 60) modifier genes, mutations of which enhance the leaf abnormalities of as2. Although genes for proteins that are involved in diverse cellular processes are known as modifiers, it has recently become clear that many modifier proteins, such as NUCLEOLIN1 (NUC1) and RNA HELICASE10 (RH10), are localized in the nucleolus. Some modifiers including ribosomal proteins are also members of the small subunit processome (SSUP). In addition, AS2 forms perinucleolar bodies partially colocalizing with chromocenters that include the condensed inactive 45S ribosomal RNA genes. AS2 participates in maintaining CpG methylation in specific exons of ETT/ARF3. NUC1 and RH10 genes are also involved in maintaining the CpG methylation levels and repressing ETT/ARF3 transcript levels. AS2 and nucleolus-localizing modifiers might cooperatively repress ETT/ARF3 to develop symmetric flat leaves. These results raise the possibility of a nucleolus-related epigenetic repression system operating for developmental genes unique to plants and predict that AS2 could be a molecule with novel functions that cannot be explained by the conventional concept of transcription factors.

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.

H3K23me1 is an evolutionarily conserved histone modification associated with CG DNA methylation in Arabidopsis.

Amino-terminal tails of histones are targets for diverse post-translational modifications whose combinatorial action may constitute a code that will be read and interpreted by cellular proteins to define particular transcriptional states. Here, we describe monomethylation of histone H3 lysine 23 (H3K23me1) as a histone modification not previously described in plants. H3K23me1 is an evolutionarily conserved mark in diverse species of flowering plants. Chromatin immunoprecipitation followed by high-throughput sequencing in Arabidopsis thaliana showed that H3K23me1 was highly enriched in pericentromeric regions and depleted from chromosome arms. In transposable elements it co-localized with CG, CHG and CHH DNA methylation as well as with the heterochromatic histone mark H3K9me2. Transposable elements are often rich in H3K23me1 but different families vary in their enrichment: LTR-Gypsy elements are most enriched and RC/Helitron elements are least enriched. The histone methyltransferase KRYPTONITE and normal DNA methylation were required for normal levels of H3K23me1 on transposable elements. Immunostaining experiments confirmed the pericentromeric localization and also showed mild enrichment in less condensed regions. Accordingly, gene bodies of protein-coding genes had intermediate H3K23me1 levels, which coexisted with CG DNA methylation. Enrichment of H3K23me1 along gene bodies did not correlate with transcription levels. Together, this work establishes H3K23me1 as a so far undescribed component of the plant histone code.

The Evolutionary Dynamics of Orthologs That Shift in Gene Body Methylation between Arabidopsis Species.

DNA methylation labels a specific subset of genes in plant genomes. Recent work has shown that this gene-body methylation (gbM) is a conserved feature of orthologs, because highly methylated genes in one species tend to be highly methylated in another. In this study, we examined the exceptions to that rule by identifying genes that differ in gbM status between two plant species-Arabidopsis thaliana and Arabidopsis lyrata. Using Capsella grandiflora as an outgroup, we polarized the loss and gain of gbM for orthologs in the Arabidopsis lineage. Our survey identified a few hundred genes that differed between ingroup species, out of ∼9,000 orthologs. The estimated rate of gbM gain was ∼2 × 10-9 per gene per year for both ingroup taxa and was similar to the loss rate in A. lyrata. In contrast, A. thaliana had a ∼3-fold higher estimated rate of gbM loss per gene, consistent with a recent diminishment of genome size. As in previous studies, we found that body-methylated genes were expressed broadly across tissues and were also longer than other genic sets. Genes that differed in gbM status exhibited higher variance in expression between species than genes that were body-methylated in both species. Moreover, the gain of gbM in one lineage tended to be associated with an increase of expression in that lineage. The genes that varied in gbM status between species varied more significantly in length between species than other sets of genes; we hypothesize that length is a key feature in the transition between body-methylated and nonmethylated genes.

Evolutionary Consequences of DNA Methylation in a Basal Metazoan.

Gene body methylation (gbM) is an ancestral and widespread feature in Eukarya, yet its adaptive value and evolutionary implications remain unresolved. The occurrence of gbM within protein-coding sequences is particularly puzzling, because methylation causes cytosine hypermutability and hence is likely to produce deleterious amino acid substitutions. We investigate this enigma using an evolutionarily basal group of Metazoa, the stony corals (order Scleractinia, class Anthozoa, phylum Cnidaria). We show that patterns of coral gbM are similar to other invertebrate species, predicting wide and active transcription and slower sequence evolution. We also find a strong correlation between gbM and codon bias, resulting from systematic replacement of CpG bearing codons. We conclude that gbM has strong effects on codon evolution and speculate that this may influence establishment of optimal codons.

Arabidopsis EDM2 promotes IBM1 distal polyadenylation and regulates genome DNA methylation patterns.

DNA methylation is important for the silencing of transposons and other repetitive elements in many higher eukaryotes. However, plant and mammalian genomes have evolved to contain repetitive elements near or inside their genes. How these genes are kept from being silenced by DNA methylation is not well understood. A forward genetics screen led to the identification of the putative chromatin regulator Enhanced Downy Mildew 2 (EDM2) as a cellular antisilencing factor and regulator of genome DNA methylation patterns. EDM2 contains a composite Plant Homeo Domain that recognizes both active and repressive histone methylation marks at the intronic repeat elements in genes such as the Histone 3 lysine 9 demethylase gene Increase in BONSAI Methylation 1 (IBM1) and is necessary for maintaining the expression of these genes by promoting mRNA distal polyadenylation. Because of its role in maintaining IBM1 expression, EDM2 is required for preventing CHG methylation in the bodies of thousands of genes. Our results thus increase the understanding of antisilencing, genome methylation patterns, and regulation of alternative RNA processing by intronic heterochromatin.

DNA methylation in social insects: how epigenetics can control behavior and longevity.

In eusocial insects, genetically identical individuals can exhibit striking differences in behavior and longevity. The molecular basis of such phenotypic plasticity is of great interest to the scientific community. DNA methylation, as well as other epigenetic signals, plays an important role in modulating gene expression and can therefore establish, sustain, and alter organism-level phenotypes, including behavior and life span. Unlike DNA methylation in mammals, DNA methylation in insects, including eusocial insects, is enriched in gene bodies of actively expressed genes. Recent investigations have revealed the important role of gene body methylation in regulating gene expression in response to intrinsic and environmental factors. In this review, we summarize recent advances in DNA methylation research and discuss its significance in our understanding of the epigenetic underpinnings of behavior and longevity.

Virus-induced gene silencing in transgenic plants: transgene silencing and reactivation associate with two patterns of transgene body methylation.

We used bisulfite sequencing to study the methylation of a viral transgene whose expression was silenced upon plum pox virus infection of the transgenic plant and its subsequent recovery as a consequence of so-called virus-induced gene silencing (VIGS). VIGS was associated with a general increase in the accumulation of small RNAs corresponding to the coding region of the viral transgene. After VIGS, the transgene promoter was not methylated and the coding region showed uneven methylation, with the 5' end being mostly unmethylated in the recovered tissue or mainly methylated at CG sites in regenerated silenced plants. The methylation increased towards the 3' end, which showed dense methylation in all three contexts (CG, CHG and CHH). This methylation pattern and the corresponding silenced status were maintained after plant regeneration from recovered silenced tissue and did not spread into the promoter region, but were not inherited in the sexual offspring. Instead, a new pattern of methylation was observed in the progeny plants consisting of disappearance of the CHH methylation, similar CHG methylation at the 3' end, and an overall increase in CG methylation in the 5' end. The latter epigenetic state was inherited over several generations and did not correlate with transgene silencing and hence virus resistance. These results suggest that the widespread CG methylation pattern found in body gene bodies located in euchromatic regions of plant genomes may reflect an older silencing event, and most likely these genes are no longer silenced.

DDM1-mediated gene body DNA methylation is associated with inducible activation of defense-related genes in Arabidopsis.

BACKGROUND: Plants memorize previous pathogen attacks and are "primed" to produce a faster and stronger defense response, which is critical for defense against pathogens. In plants, cytosines in transposons and gene bodies are reported to be frequently methylated. Demethylation of transposons can affect disease resistance by regulating the transcription of nearby genes during defense response, but the role of gene body methylation (GBM) in defense responses remains unclear. RESULTS: Here, we find that loss of the chromatin remodeler decrease in DNA methylation 1 (ddm1) synergistically enhances resistance to a biotrophic pathogen under mild chemical priming. DDM1 mediates gene body methylation at a subset of stress-responsive genes with distinct chromatin properties from conventional gene body methylated genes. Decreased gene body methylation in loss of ddm1 mutant is associated with hyperactivation of these gene body methylated genes. Knockout of glyoxysomal protein kinase 1 (gpk1), a hypomethylated gene in ddm1 loss-of-function mutant, impairs priming of defense response to pathogen infection in Arabidopsis. We also find that DDM1-mediated gene body methylation is prone to epigenetic variation among natural Arabidopsis populations, and GPK1 expression is hyperactivated in natural variants with demethylated GPK1. CONCLUSIONS: Based on our collective results, we propose that DDM1-mediated GBM provides a possible regulatory axis for plants to modulate the inducibility of the immune response.

Dynamic DNA methylation turnover in gene bodies is associated with enhanced gene expression plasticity in plants.

BACKGROUND: In several eukaryotes, DNA methylation occurs within the coding regions of many genes, termed gene body methylation (GbM). Whereas the role of DNA methylation on the silencing of transposons and repetitive DNA is well understood, gene body methylation is not associated with transcriptional repression, and its biological importance remains unclear. RESULTS: We report a newly discovered type of GbM in plants, which is under constitutive addition and removal by dynamic methylation modifiers in all cells, including the germline. Methylation at Dynamic GbM genes is removed by the DRDD demethylation pathway and added by an unknown source of de novo methylation, most likely the maintenance methyltransferase MET1. We show that the Dynamic GbM state is present at homologous genes across divergent lineages spanning over 100 million years, indicating evolutionary conservation. We demonstrate that Dynamic GbM is tightly associated with the presence of a promoter or regulatory chromatin state within the gene body, in contrast to other gene body methylated genes. We find Dynamic GbM is associated with enhanced gene expression plasticity across development and diverse physiological conditions, whereas stably methylated GbM genes exhibit reduced plasticity. Dynamic GbM genes exhibit reduced dynamic range in drdd mutants, indicating a causal link between DNA demethylation and enhanced gene expression plasticity. CONCLUSIONS: We propose a new model for GbM in regulating gene expression plasticity, including a novel type of GbM in which increased gene expression plasticity is associated with the activity of DNA methylation writers and erasers and the enrichment of a regulatory chromatin state.