Transcriptome-wide m6A methylation profile reveals regulatory networks in roots of barley under cadmium stress

Cadmium (Cd) pollutants restrict crop yield and food security in long-term agricultural activities. Crops have evolved adaptive strategies under Cd condition, however, the transcriptional regulatory mechanism of Cd-tolerant genes remains to be largely illustrated. In this study, barley roots were exposed to 5 µM CdCl2 for physiological response and transcriptome-wide m6A methylation profile. Cd stress inhibited root growth after 7 d Cd treatment, which is mainly associated with inhibited absorption of Mn. After Cd treatment, 8151 significantly modified m6A sites and 3920 differentially expressed genes were identified. Transcriptome-wide m6A hypermethylation was widely induced by Cd stress and enriched near the stop codon and 3′ UTR regions. Among 435 m6A modified DEGs, 319 hypermethylated genes were up-regulated and 84 hypomethylated genes were down-regulated, respectively, indicating a positive correlation of m6A methylation and expression. But well-known Cd transporter genes (HvNramp5, HvIRT1, HvHMA3, etc.) were not modified by m6A methylation, except for ABC transporters. We further found key Cd-responding regulatory genes were positively modulated with m6A methylation, including MAPK, WRKY and MYB members. This study proposed a transcriptional regulatory network of Cd stress response in barley roots, which may provide new insight into gene manipulation of controlling low Cd accumulation for crops.