To calculate the association between gene expression and chromatin accessibility differences for each gene, we assigned to each gene with at least 1 DAR (FDR < 0.2) in the promoter region the log2 fold change of the DAR with the highest ATAC signal, as well as the log2 fold change of the gene’s expression.