Variation of DNA methylation patterns associated with gene expression in rice (Oryza sativa) exposed to cadmium

We report genome-wide single-base resolution maps of methylated cytosines and transcriptome change in Cd-exposed rice. Widespread differences were identified in CG and non-CG methylation marks between Cd-exposed and Cd-free rice genomes. There are 2320 non-redundant differentially methylated regions detected in the genome. RNA sequencing revealed 2092 DNA methylation-modified genes differentially expressed under Cd exposure. More genes were found hypermethylated than those hypomethylated in CG, CHH and CHG (where H is A, C or T) contexts in upstream, gene body and downstream regions. Many of the genes were involved in stress response, metal transport and transcription factors. Most of the DNA methylation-modified genes were transcriptionally altered under Cd stress. A subset of loss of function mutants defective in DNA methylation and histone modification activities was used to identify transcript abundance of selected genes. Compared with wide type, mutation of MET1 and DRM2 resulted in general lower transcript levels of the genes under Cd stress. Transcripts of OsIRO2, OsPR1b and Os09g02214 in drm2 were significantly reduced. A commonly used DNA methylation inhibitor 5-azacytidine was employed to investigate whether DNA demethylation affected physiological consequences. 5-azacytidine provision decreased general DNA methylation levels of selected genes, but promoted growth of rice seedlings and Cd accumulation in rice plant.