Abstract
DNA methylation is one of the epigenetic mechanisms that govern gene regulation in response to abiotic stress in plants. Here, we analyzed the role of epigenetic variations by exploring global DNA methylation and integrating it with differential gene expression in response to salinity stress in tolerant and sensitive chickpea genotypes. Genome-wide DNA methylation profiles showed higher CG methylation in the gene body regions and higher CHH methylation in the TE body regions. The analysis of differentially methylated regions (DMRs) suggested more hyper-methylation in response to stress in the tolerant genotype compared to the sensitive genotype. We observed higher enrichment of CG DMRs in genes and CHH DMRs in transposable elements (TEs). A positive correlation of gene expression with CG gene body methylation was observed. The enrichment analysis of DMR-associated differentially expressed genes revealed they are involved in biological processes, such as lateral root development, transmembrane transporter activity, GTPase activity, and regulation of gene expression. Further, a high correlation of CG methylation with CHG and CHH methylation under salinity stress was revealed, suggesting crosstalk among the methylation contexts. Further, we observed small RNA-mediated CHH hypermethylation in TEs. Overall, the interplay between DNA methylation, small RNAs, and gene expression provides new insights into the regulatory mechanism underlying salinity stress response in chickpeas.