Whole genome bisulfite sequencing methylome analysis of mulberry (Morus alba) reveals epigenome modifications in response to drought stress

DNA methylation plays a significant role in many biological processes. Although some studies of DNA methylation have been performed in woody plant, none is known about the methylation patterns of mulberry (Morus alba). In this study, we performed whole genome bisulfite sequencing under drought stress to generate a methylated cytosines map and assessed the effects of the changes on gene expression combined with transcriptomics. We found that the percentage of methylated cytosines varied depending on the local sequence context (CG, CHG and CHH) and external treatment (control, CK; drought stress, DS). The methylation levels under DS were 8.64% higher than that of CK, and differences that were mainly due to the contribution of mCG (6.24%). Additionally, there were 3,243 different methylation and expression associated genes. In addition, methylated genes were enriched within GO subcategories including catalytic activity, cellular process, metabolic process, response to stimulus and regulation of biological process. This is the first study to comprehensively present methylation patterns in mulberry and reveal widespread DNA methylation changes in response to drought stress, which has the potential to enhance our understanding of links between DNA methylation and the modulation of gene expression in plants subjected to abiotic stresses.