Coordination of m6A mRNA methylation and gene transcriptome in rice response to cadmium stress
Abstract N6-methyladenosine (m6A) is the most prevalent internal modification present in mRNAs of all higher eukaryotes. However, the role of the m6A methylomes in rice is still poorly understood. With the development of MeRIP-seq technique, in-depth identification of mRNAs with m6A modification becomes feasible. We investigated the m6A methylomes in roots of cadmium (Cd) group and compared that with the roots in the control (CK) group by m6A sequencing, in 9311 and Nipponbare (NIP), respectively. The results indicated that Cd leads to altered modification profile in 3,406 differential m6A peaks in 9311, and 2,065 differential m6A peaks in NIP. KEGG pathway analysis of genes with differentially modified m6A peaks indicates that the “phenylalanine”, “tyrosine and tryptophan biosynthesis”, “glycine”, “adherens junctions”, “glycerophospholipid metabolism” and “threonine metabolism” signaling pathways may be associated with abnormal roots development of rice due to exposure to cadmium in 9311. “Arginine”, “proline metabolism”, “glycerolipid”, “protein processing in endoplasmic reticulum”, metabolism pathways were significantly enriched in genes with differentially modified m6A peaks in NIP. Different from that in Arabidopsis, the m6A peak (m6A-modified nucleotide position on mRNAs) distribution exhibits preference toward both the stop codon and 3′UTRs region. These findings provide a resource for plant RNA epi-transcriptomics studies and further enlarge our knowledge on the function of RNA m6A modification in plants.