Abstract
Drought stress remains one of the most detrimental environmental cues affecting plant growth and survival. In this work, the DNA methylome changes in mulberry leaves under drought stress (EG) and control (CK) and their impact on gene regulation were investigated by MethylRAD sequencing. The results show 138,464 (37.37%) and 56,241 (28.81%) methylation at the CG and CWG sites (W = A or T), respectively, in the mulberry genome between drought stress and control. The distribution of the methylome was prevalent in the intergenic, exonic, intronic and downstream regions of the mulberry plant genome. In addition, we discovered 170 DMGs (129 in CG sites and 41 in CWG sites) and 581 DMS (413 in CG sites and 168 in CWG sites). Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicates that phenylpropanoid biosynthesis, spliceosome, amino acid biosynthesis, carbon metabolism, RNA transport, plant hormone, signal transduction pathways, and quorum sensing play a crucial role in mulberry response to drought stress. Furthermore, the qRT-PCR analysis indicates that the selected 23 genes enriched in the KEGG pathways are differentially expressed, and 86.96% of the genes share downregulated methylation and 13.04% share upregulation methylation status, indicating the complex link between DNA methylation and gene regulation. This study serves as fundamentals in discovering the epigenomic status and the pathways that will significantly enhance mulberry breeding for adaptation to a wide range of environments.