Role of Non-Coding Regulatory Elements in the Control of GR-Dependent Gene Expression

The glucocorticoid receptor (GR, also known as NR3C1) coordinates molecular responses to stress. It is a potent transcription activator and repressor that influences hundreds of genes. Enhancers are non-coding DNA regions outside of the core promoters that increase transcriptional activity via long-distance interactions. Active GR binds to pre-existing enhancer sites and recruits further factors, including EP300, a known transcriptional coactivator. However, it is not known how the timing of GR-binding-induced enhancer remodeling relates to transcriptional changes. Here we analyze data from the ENCODE project that provides ChIP-Seq and RNA-Seq data at distinct time points after dexamethasone exposure of human A549 epithelial-like cell line. This study aimed to investigate the temporal interplay between GR binding, enhancer remodeling, and gene expression. By investigating a single distal GR-binding site for each differentially upregulated gene, we show that transcriptional changes follow GR binding, and that the largest enhancer remodeling coincides in time with the highest gene expression changes. A detailed analysis of the time course showed that for upregulated genes, enhancer activation persists after gene expression changes settle. Moreover, genes with the largest change in EP300 binding showed the highest expression dynamics before the peak of EP300 recruitment. Overall, our results show that enhancer remodeling may not directly be driving gene expression dynamics but rather be a consequence of expression activation.