The dynamics of DNA methylation in maize roots under Pb stress

Plants adapt to adverse conditions through a series of physiological, cellular, and molecular processes, culminating in stress tolerance. However, little is known about the associated regulatory mechanisms at the epigenetic level in maize under lead (Pb) stress. Therefore, in this study, we aimed to compare DNA methylation profiles during the dynamic development of maize roots following Pb treatment to identify candidate genes involved in the response to Pb stress. Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) was used to investigate the genome-wide DNA methylation patterns in maize roots under normal condition (A1) and 3 mM Pb(NO3)2 stress for 12 h (K2), 24 h (K3) and 48 h (K4). The results showed that the average methylation density was the highest in CpG islands (CGIs), followed by the intergenic regions. Within the gene body, the methylation density of the introns was higher than those of the UTRs and exons. In total, 3857 methylated genes were found in 4 tested samples, including 1805 differentially methylated genes for K2 versus A1, 1508 for K3 versus A1, and 1660 for K4 versus A1. Further analysis showed that 140 genes exhibited altered DNA methylation in all three comparisons, including some well-known stress-responsive transcription factors and proteins, such as MYB, AP2/ERF, bZIP, serine-threonine/tyrosine-proteins, pentatricopeptide repeat proteins, RING zinc finger proteins, F-box proteins, leucine-rich repeat proteins and tetratricopeptide repeat proteins. This study revealed the genome-scale DNA methylation patterns of maize roots in response to Pb exposure and identified candidate genes that potentially regulate root dynamic development under Pb stress at the methylation level.