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.

Plant microRNAs (miRNAs) have been shown to play critical roles in plant development. In this study, we employed small RNA combined with degradome sequencing to survey development-related miRNAs and their validated targets during wheat grain development. A total of 186 known miRNAs and 37 novel miRNAs were identified in four small RNA libraries. Moreover, a miRNA-like long hairpin locus was first identified to produce 21~22-nt phased siRNAs that act in trans to cleave target mRNAs. A comparison of the miRNAomes revealed that 55 miRNA families were differentially expressed during the grain development. Predicted and validated targets of these development-related miRNAs are involved in different cellular responses and metabolic processes including cell proliferation, auxin signaling, nutrient metabolism and gene expression. This study provides insight into the complex roles of miRNAs and their targets in regulating wheat grain development.