Fig. 3

Computational pipeline for genome-wide bisulfite sequencing data analysis. Reads from bisulfite sequencing are first aligned to the reference genome. The alignment data may be visualized in different tracks for comparison. After methylation calling, the bulk methylation level and genome-wide methylation level are calculated and plotted, and DMRs are determined. To perform an integrative analysis, DNA methylation data are coupled with gene expression, e.g., differentially expressed genes (DEGs), to delineate the regulatory role of DNA methylation