Differential enrichment analysis (fold change > 2 and adjusted P value < 0.05) was performed using DESeq2 to gain a quantitative picture of H3K27Ac dynamics upon butyrate treatment. Distribution of hyper-acetylated (in green), hypo-acetylated (in red) and unchanged (in purple) baseline (0 h) H3K27Ac+ peaks in response to 3, 12 and 24 hours of butyrate treatment.