5-Methylcytosine (mC) is a well-characterized DNA modification, and is also predominantly reported in abundant non-coding RNAs in both prokaryotes and eukaryotes. However, the distribution and biological functions of mC in plant mRNAs remain largely unknown. Here, we report transcriptome-wide profiling of RNA mC in Arabidopsis thaliana by applying mC RNA immunoprecipitation followed by a deep-sequencing approach (mC-RIP-seq). LC-MS/MS and dot blot analyses reveal a dynamic pattern of mC mRNA modification in various tissues and at different developmental stages. mC-RIP-seq analysis identified 6045 mC peaks in 4465 expressed genes in young seedlings. We found that mC is enriched in coding sequences with two peaks located immediately after start codons and before stop codons, and is associated with mRNAs with low translation activity. We further demonstrated that an RNA (cytosine-5)-methyltransferase, tRNA-specific methyltransferase 4B (TRM4B), exhibits mC RNA methyltransferase activity. Mutations in TRM4B display defects in root development and decreased mC peaks. TRM4B affects the transcript levels of the genes involved in root development, which is positively correlated with their mRNA stability and mC levels. Our results suggest that mC in mRNA is a new epitranscriptome marker inArabidopsis, and that regulation of this modification is an integral part of gene regulatory networks underlying plant development.

DNA methylation on N-adenine (6mA) has recently been found to be a potentially epigenetic mark in several unicellular and multicellular eukaryotes. However, its distribution patterns and potential functions in land plants, which are primary producers for most ecosystems, remain largely unknown. Here we report global profiling of 6mA sites at single-nucleotide resolution in the genome of Arabidopsis thaliana at different developmental stages using single-molecule real-time sequencing. 6mA sites are widely distributed across the Arabidopsis genome and enriched over the pericentromeric heterochromatin regions. 6mA occurs more frequently in gene bodies than intergenic regions. Analysis of 6mA methylomes and RNA sequencing data demonstrates that 6mA frequency positively correlates with the gene expression level and the transition from vegetative to reproductive growth in Arabidopsis. Our results uncover 6mA as a DNA mark associated with actively expressed genes in Arabidopsis, suggesting that 6mA serves as a hitherto unknown epigenetic mark in land plants.

N 6 -Methyladenine (6mA) DNA methylation has recently been implicated as a potential new epigenetic marker in eukaryotes, including the dicot model Arabidopsis thaliana. However, the conservation and divergence of 6mA distribution patterns and functions in plants remain elusive. Here we report high-quality 6mA methylomes at single-nucleotide resolution in rice based on substantially improved genome sequences of two rice cultivars, Nipponbare (Nip; Japonica) and 93-11 (Indica). Analysis of 6mA genomic distribution and its association with transcription suggest that 6mA distribution and function is rather conserved between rice and Arabidopsis. We found that 6mA levels are positively correlated with the expression of key stressrelated genes, which may be responsible for the difference in stress tolerance between Nip and 93-11. Moreover, we showed that mutations in DDM1 cause defects in plant growth and decreased 6mA level. Our results reveal that 6mA is a conserved DNA modification that is positively associated with gene expression and contributes to key agronomic traits in plants.