For the GSEA analyses, the log2 fold changes from the differential expression analysis were used as the ranks, and genes with at least 1 DAR in the promoter region were used as the gene sets, separated into genes with positive log2 fold changes, and those with negative log2 fold changes.