Regulation of gene expression by DNA methylation is crucial for defining cellular identities and coordinating organism-wide developmental programs in many organisms. In plants, modulation of DNA methylation in response to environmental conditions represents a potentially robust mechanism to regulate gene expression networks; however, examples of dynamic DNA methylation are largely limited to gene imprinting. Here we report an unexpected role for DNA methylation in regulation of the Arabidopsis thaliana immune system. Profiling the DNA methylomes of plants exposed to bacterial pathogen, avirulent bacteria, or salicylic acid (SA) hormone revealed numerous stress-induced differentially methylated regions, many of which were intimately associated with differentially expressed genes. In response to SA, transposon-associated differentially methylated regions, which were accompanied by up-regulation of 21-nt siRNAs, were often coupled to transcriptional changes of the transposon and/or the proximal gene. Thus, dynamic DNA methylation changes within repetitive sequences or transposons can regulate neighboring genes in response to SA stress.

Natural epigenetic variation provides a source for the generation of phenotypic diversity, but to understand its contribution to such diversity, its interaction with genetic variation requires further investigation. Here we report population-wide DNA sequencing of genomes, transcriptomes and methylomes of wild Arabidopsis thaliana accessions. Single cytosine methylation polymorphisms are not linked to genotype. However, the rate of linkage disequilibrium decay amongst differentially methylated regions targeted by RNA-directed DNA methylation is similar to the rate for single nucleotide polymorphisms. Association analyses of these RNA-directed DNA methylation regions with genetic variants identified thousands of methylation quantitative trait loci, which revealed the population estimate of genetically dependent methylation variation. Analysis of invariably methylated transposons and genes across this population indicates that loci targeted by RNA-directed DNA methylation are epigenetically activated in pollen and seeds, which facilitates proper development of these structures.

Methods for profiling DNA methylation differ in the physical principles used to detect modified cytosines. Harris et al. compare the performances of four sequencing-based technologies for genome-wide analysis of DNA methylation and combine two methods to enable detection of allelic differences in epigenetic marks. Analysis of DNA methylation patterns relies increasingly on sequencing-based profiling methods. The four most frequently used sequencing-based technologies are the bisulfite-based methods MethylC-seq and reduced representation bisulfite sequencing (RRBS), and the enrichment-based techniques methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA binding domain sequencing (MBD-seq). We applied all four methods to biological replicates of human embryonic stem cells to assess their genome-wide CpG coverage, resolution, cost, concordance and the influence of CpG density and genomic context. The methylation levels assessed by the two bisulfite methods were concordant (their difference did not exceed a given threshold) for 82% for CpGs and 99% of the non-CpG cytosines. Using binary methylation calls, the two enrichment methods were 99% concordant and regions assessed by all four methods were 97% concordant. We combined MeDIP-seq with methylation-sensitive restriction enzyme (MRE-seq) sequencing for comprehensive methylome coverage at lower cost. This, along with RNA-seq and ChIP-seq of the ES cells enabled us to detect regions with allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression.