Genomic dissection of enhancers uncovers principles of combinatorial regulation and cell type-specific wiring of enhancer-promoter contacts

Probabilistic association of differentially expressed genes with cis-regulatory elements

Differential gene expression in response to perturbations is mediated at least in part by changes in binding of transcription factors (TFs) and other proteins at specific genomic regions. Association of these cis-regulatory elements (CREs) with their target genes is a challenging task that is essential to address many biological and mechanistic questions. Many current approaches rely on chromatin conformation capture techniques or single-cell correlational methods to establish CRE-to-gene associations. These methods can be effective but have limitations, including resolution, gaps in detectable association distances, and cost. As an alternative, we have developed DegCre, a nonparametric method that evaluates correlations between measurements of perturbation-induced differential gene expression and differential regulatory signal at CREs to score possible CRE-to-gene associations. It has several unique features, including the ability to use any type of CRE activity measurement, yield probabilistic scores for CRE-to-gene pairs, and assess CRE-to-gene pairings across a wide range of sequence distances. We apply DegCre to six data sets, each using different perturbations and containing a variety of regulatory signal measurements, including chromatin openness, histone modifications, and TF occupancy. To test their efficacy, we compare DegCre associations to Hi-C loop calls and CRISPR-validated CRE-to-gene associations, establishing good performance by DegCre that is comparable or superior to competing methods. DegCre is a novel approach to the association of CREs to genes from a perturbation-differential perspective, with strengths that are complementary to existing approaches and allow for new insights into gene regulation.

Glucocorticoid stimulation induces regionalized gene responses within topologically associating domains

Transcription-factor binding to cis-regulatory regions regulates the gene expression program of a cell, but occupancy is often a poor predictor of the gene response. Here, we show that glucocorticoid stimulation led to the reorganization of transcriptional coregulators MED1 and BRD4 within topologically associating domains (TADs), resulting in active or repressive gene environments. Indeed, we observed a bias toward the activation or repression of a TAD when their activities were defined by the number of regions gaining and losing MED1 and BRD4 following dexamethasone (Dex) stimulation. Variations in Dex-responsive genes at the RNA levels were consistent with the redistribution of MED1 and BRD4 at the associated cis-regulatory regions. Interestingly, Dex-responsive genes without the differential recruitment of MED1 and BRD4 or binding by the glucocorticoid receptor were found within TADs, which gained or lost MED1 and BRD4, suggesting a role of the surrounding environment in gene regulation. However, the amplitude of the response of Dex-regulated genes was higher when the differential recruitment of the glucocorticoid receptor and transcriptional coregulators was observed, reaffirming the role of transcription factor-driven gene regulation and attributing a lesser role to the TAD environment. These results support a model where a signal-induced transcription factor induces a regionalized effect throughout the TAD, redefining the notion of direct and indirect effects of transcription factors on target genes.

A non-coding GWAS variant impacts anthracycline-induced cardiotoxic phenotypes in human iPSC-derived cardiomyocytes

Anthracyclines, widely used to treat breast cancer, have the potential for cardiotoxicity. We have previously identified and validated a germline single nucleotide polymorphism, rs28714259, associated with an increased risk of anthracycline-induced heart failure. We now provide insights into the mechanism by which rs28714259 might confer increased risk of cardiac damage. Using hiPSC-derived cardiomyocyte cell lines with either intrinsic polymorphism or CRISPR-Cas9-mediated deletion of rs28714259 locus, we demonstrate that glucocorticoid receptor signaling activated by dexamethasone pretreatment prior to doxorubicin exposure preserves cardiomyocyte viability and contractility in cardiomyocytes containing the major allele. Homozygous loss of the rs28714259 major allele diminishes dexamethasone's protective effect. We further demonstrate that the risk allele of rs28714259 disrupts glucocorticoid receptor and rs28714259 binding affinity. Finally, we highlight the activation of genes and pathways involved in cardiac hypertrophy signaling that are blocked by the risk allele, suggesting a decreased adaptive survival response to doxorubicin-related stress.

Crosstalk between p38 MAPK and GR Signaling

The p38 MAPK is a signaling pathway important for cells to respond to environmental and intracellular stress. Upon activation, the p38 kinase phosphorylates downstream effectors, which control the inflammatory response and coordinate fundamental cellular processes such as proliferation, apoptosis, and differentiation. Dysregulation of this signaling pathway has been linked to inflammatory diseases and cancer. Secretion of glucocorticoids (GCs) is a classical endocrine response to stress. The glucocorticoid receptor (GR) is the primary effector of GCs and plays an important role in the regulation of cell metabolism and immune response by influencing gene expression in response to hormone-dependent activation. Its ligands, the GCs or steroids, in natural or synthetic variation, are used as standard therapy for anti-inflammatory treatment, severe asthma, autoimmune diseases, and several types of cancer. Several years ago, the GR was identified as one of the downstream targets of p38, and, at the same time, it was shown that glucocorticoids could influence p38 signaling. In this review, we discuss the role of the crosstalk between the p38 and GR in the regulation of gene expression in response to steroids and comprehend the importance and potential of this interplay in future clinical applications.

Multimodal regulatory elements within a hormone-specific super enhancer control a heterogeneous transcriptional response

The hormone-stimulated glucocorticoid receptor (GR) modulates transcription by interacting with thousands of enhancers and GR binding sites (GBSs) throughout the genome. Here, we examined the effects of GR binding on enhancer dynamics and investigated the contributions of individual GBSs to the hormone response. Hormone treatment resulted in genome-wide reorganization of the enhancer landscape in breast cancer cells. Upstream of the DDIT4 oncogene, GR bound to four sites constituting a hormone-dependent super enhancer. Three GBSs were required as hormone-dependent enhancers that differentially promoted histone acetylation, transcription frequency, and burst size. Conversely, the fourth site suppressed transcription and hormone treatment alleviated this suppression. GR binding within the super enhancer promoted a loop-switching mechanism that allowed interaction of the DDIT4 TSS with the active GBSs. The unique functions of each GR binding site contribute to hormone-induced transcriptional heterogeneity and demonstrate the potential for targeted modulation of oncogene expression.

Assessing genome-wide dynamic changes in enhancer activity during early mESC differentiation by FAIRE-STARR-seq

Embryonic stem cells (ESCs) can differentiate into any given cell type and therefore represent a versatile model to study the link between gene regulation and differentiation. To quantitatively assess the dynamics of enhancer activity during the early stages of murine ESC differentiation, we analyzed accessible genomic regions using STARR-seq, a massively parallel reporter assay. This resulted in a genome-wide quantitative map of active mESC enhancers, in pluripotency and during the early stages of differentiation. We find that only a minority of accessible regions is active and that such regions are enriched near promoters, characterized by specific chromatin marks, enriched for distinct sequence motifs, and modeling shows that active regions can be predicted from sequence alone. Regions that change their activity upon retinoic acid-induced differentiation are more prevalent at distal intergenic regions when compared to constitutively active enhancers. Further, analysis of differentially active enhancers verified the contribution of individual TF motifs toward activity and inducibility as well as their role in regulating endogenous genes. Notably, the activity of retinoic acid receptor alpha (RARα) occupied regions can either increase or decrease upon the addition of its ligand, retinoic acid, with the direction of the change correlating with spacing and orientation of the RARα consensus motif and the co-occurrence of additional sequence motifs. Together, our genome-wide enhancer activity map elucidates features associated with enhancer activity levels, identifies regulatory regions disregarded by computational prediction tools, and provides a resource for future studies into regulatory elements in mESCs.

Glucocorticoid receptor wields chromatin interactions to tune transcription for cytoskeleton stabilization in podocytes

Elucidating transcription mediated by the glucocorticoid receptor (GR) is crucial for understanding the role of glucocorticoids (GCs) in the treatment of diseases. Podocyte is a useful model for studying GR regulation because GCs are the primary medication for podocytopathy. In this study, we integrated data from transcriptome, transcription factor binding, histone modification, and genome topology. Our data reveals that the GR binds and activates selective regulatory elements in podocyte. The 3D interactome captured by HiChIP facilitates the identification of remote targets of GR. We found that GR in podocyte is enriched at transcriptional interaction hubs and super-enhancers. We further demonstrate that the target gene of the top GR-associated super-enhancer is indispensable to the effective functioning of GC in podocyte. Our findings provided insights into the mechanisms underlying the protective effect of GCs on podocyte, and demonstrate the importance of considering transcriptional interactions in order to fine-map regulatory networks of GR.

Genome-wide binding potential and regulatory activity of the glucocorticoid receptor's monomeric and dimeric forms

A widely regarded model for glucocorticoid receptor (GR) action postulates that dimeric binding to DNA regulates unfavorable metabolic pathways while monomeric receptor binding promotes repressive gene responses related to its anti-inflammatory effects. This model has been built upon the characterization of the GRdim mutant, reported to be incapable of DNA binding and dimerization. Although quantitative live-cell imaging data shows GRdim as mostly dimeric, genomic studies based on recovery of enriched half-site response elements suggest monomeric engagement on DNA. Here, we perform genome-wide studies on GRdim and a constitutively monomeric mutant. Our results show that impairing dimerization affects binding even to open chromatin. We also find that GRdim does not exclusively bind half-response elements. Our results do not support a physiological role for monomeric GR and are consistent with a common mode of receptor binding via higher order structures that drives both the activating and repressive actions of glucocorticoids.

Characterization of the regulatory 5'-flanking region of bovine mucin 2 (MUC2) gene

Throughout the intestinal epithelium surface there is an intricate polymer network composed by gel-forming mucins, which plays a protective role due to the formation of a physical, chemical and immunological barrier between the organism and the environment. Mucin 2 (MUC2) is the main mucin in the small and large intestine, and it is expressed specifically in the gastrointestinal tract (GIT), which makes its promoter region an important candidate for expression of heterologous genes of biotechnological interest in the GIT of bovine and other ruminants. In order to characterize the bovine MUC2 promoter we designed primers to amplify and isolate a candidate region for this promoter. The amplified sequence was confirmed by sequencing and cloned into a plasmid vector containing the luciferase (LUC) reporter gene. The regulatory sites of the MUC2 promoter already described in the literature were used to find the putative regulatory sites in the bovine MUC2 promoter region. With these data, some deletions were performed in order to find the promoter sequence with greatest expression capacity and specificity. The constructions were tested by transient transfection assays in LoVo cells (human colorectal adenocarcinoma) and bovine fibroblasts. The quantification of the relative expression of the promoter was measured using dual-luciferase assays. Real-time PCR was performed to analyze the expression of endogenous MUC2. The results presented herein prove that the isolated sequence corresponds to the promoter of bovine MUC2 gene, since it was able to induce expression of a reporter gene in an in vitro cell culture experimental platform.

Promoter-interacting expression quantitative trait loci are enriched for functional genetic variants

Expression quantitative trait loci (eQTLs) studies provide associations of genetic variants with gene expression but fall short of pinpointing functionally important eQTLs. Here, using H3K27ac HiChIP assays, we mapped eQTLs overlapping active cis-regulatory elements that interact with their target gene promoters (promoter-interacting eQTLs, pieQTLs) in five common immune cell types (Database of Immune Cell Expression, Expression quantitative trait loci and Epigenomics (DICE) cis-interactome project). This approach allowed us to identify functionally important eQTLs and show mechanisms that explain their cell-type restriction. We also devised an approach to eQTL discovery that relies on HiChIP-based promoter interaction maps as a structural framework for deciding which SNPs to test for association with gene expression, and observe ultra-long-distance pieQTLs (>1 megabase away), including several disease-risk variants. We validated the functional role of pieQTLs using reporter assays, CRISPRi, dCas9-tiling guides and Cas9-mediated base-pair editing. In this article we present a method for functional eQTL discovery and provide insights into relevance of noncoding variants for cell-specific gene regulation and for disease association beyond conventional eQTL mapping.

Mechanisms and Clinical Applications of Glucocorticoid Steroids in Muscular Dystrophy

Glucocorticoid steroids are widely used as immunomodulatory agents in acute and chronic conditions. Glucocorticoid steroids such as prednisone and deflazacort are recommended for treating Duchenne Muscular Dystrophy where their use prolongs ambulation and life expectancy. Despite this benefit, glucocorticoid use in Duchenne Muscular Dystrophy is also associated with significant adverse consequences including adrenal suppression, growth impairment, poor bone health and metabolic syndrome. For other forms of muscular dystrophy like the limb girdle dystrophies, glucocorticoids are not typically used. Here we review the experimental evidence supporting multiple mechanisms of glucocorticoid action in dystrophic muscle including their role in dampening inflammation and myofiber injury. We also discuss alternative dosing strategies as well as novel steroid agents that are in development and testing, with the goal to reduce adverse consequences of prolonged glucocorticoid exposure while maximizing beneficial outcomes.

A preliminary computational outputs versus experimental results: Application of sTRAP, a biophysical tool for the analysis of SNPs of transcription factor-binding sites

BACKGROUND

In the human genome, the transcription factors (TFs) and transcription factor-binding sites (TFBSs) network has a great regulatory function in the biological pathways. Such crosstalk might be affected by the single-nucleotide polymorphisms (SNPs), which could create or disrupt a TFBS, leading to either a disease or a phenotypic defect. Many computational resources have been introduced to predict the TFs binding variations due to SNPs inside TFBSs, sTRAP being one of them.

METHODS

A literature review was performed and the experimental data for 18 TFBSs located in 12 genes was provided. The sequences of TFBS motifs were extracted using two different strategies; in the size similar with synthetic target sites used in the experimental techniques, and with 60 bp upstream and downstream of the SNPs. The sTRAP (http://trap.molgen.mpg.de/cgi-bin/trap_two_seq_form.cgi) was applied to compute the binding affinity scores of their cognate TFs in the context of reference and mutant sequences of TFBSs. The alternative bioinformatics model used in this study was regulatory analysis of variation in enhancers (RAVEN; http://www.cisreg.ca/cgi-bin/RAVEN/a). The bioinformatics outputs of our study were compared with experimental data, electrophoretic mobility shift assay (EMSA).

RESULTS

In 6 out of 18 TFBSs in the following genes COL1A1, Hb ḉᴪ, TF, FIX, MBL2, NOS2A, the outputs of sTRAP were inconsistent with the results of EMSA. Furthermore, no p value of the difference between the two scores of binding affinity under the wild and mutant conditions of TFBSs was presented. Nor, were any criteria for preference or selection of any of the measurements of different matrices used for the same analysis.

CONCLUSION

Our preliminary study indicated some paradoxical results between sTRAP and experimental data. However, to link the data of sTRAP to the biological functions, its optimization via experimental procedures with the integration of expanded data and applying several other bioinformatics tools might be required.

IL-1 Transcriptional Responses to Lipopolysaccharides Are Regulated by a Complex of RNA Binding Proteins

The IL1A and IL1B genes lie in close proximity on chromosome 2 near the gene for their natural inhibitor, IL1RN Despite diverse functions, they are all three inducible through TLR4 signaling but with distinct kinetics. This study analyzed transcriptional induction kinetics, chromosome looping, and enhancer RNA production to understand the distinct regulation of these three genes in human cells. IL1A, IL1B, and IL1RN were rapidly induced after stimulation with LPS; however, IL1B mRNA production was less inhibitable by iBET151, suggesting it does not use pause-release regulation. Surprisingly, chromatin looping contacts between IL1A and IL1B were highly intermingled, although those of IL1RN were distinct, and we focused on comparing IL1A and IL1B transcriptional pathways. Our studies demonstrated that enhancer RNAs were produced from a subset of the regulatory regions, that they were critical for production of the mRNAs, and that they bound a diverse array of RNA binding proteins, including p300 but not CBP. We, furthermore, demonstrated that recruitment of p300 was dependent on MAPKs. Integrator is another RNA binding protein recruited to the promoters and enhancers, and its recruitment was more dependent on NF-κB than MAPKs. We found that integrator and NELF, an RNA polymerase II pausing protein, were associated with RNA in a manner that facilitated interaction. We conclude that IL1A and IL1B share many regulatory contacts, signaling pathways, and interactions with enhancer RNAs. A complex of protein interactions with enhancer RNAs emphasize the role of enhancer RNAs and the overall structural aspects of transcriptional regulation.

Sufficiency analysis of estrogen responsive enhancers using synthetic activators

Enhancers bound by synthetic activators can recreate a transcriptional response to estrogen, drive different levels of gene expression, and work independently to regulate transcription.

Multiple regulatory regions bound by the same transcription factor have been shown to simultaneously control a single gene’s expression. However, it remains unclear how these regulatory regions combine to regulate transcription. Here, we test the sufficiency of promoter-distal estrogen receptor α-binding sites (ERBSs) for activating gene expression by recruiting synthetic activators in the absence of estrogens. Targeting either dCas9-VP16(10x) or dCas9-p300(core) to ERBS induces H3K27ac and activates nearby expression in a manner similar to an estrogen induction, with dCas9-VP16(10x) acting as a stronger activator. The sufficiency of individual ERBSs is highly correlated with their necessity, indicating an inherent activation potential that is associated with the binding of RNA polymerase II and several transcription factors. By targeting ERBS combinations, we found that ERBSs work independently to control gene expression when bound by synthetic activators. The sufficiency results contrast necessity assays that show synergy between these ERBSs, suggesting that synergy occurs between ERBSs in terms of activator recruitment, whereas directly recruiting activators leads to independent effects on gene expression.

Androgen receptor enhancer usage and the chromatin regulatory landscape in human prostate cancers

The androgen receptor (AR) is commonly known as a key transcription factor in prostate cancer development, progression and therapy resistance. Genome-wide chromatin association studies revealed that transcriptional regulation by AR mainly depends on binding to distal regulatory enhancer elements that control gene expression through chromatin looping to gene promoters. Changes in the chromatin epigenetic landscape and DNA sequence can locally alter AR-DNA-binding capacity and consequently impact transcriptional output and disease outcome. The vast majority of reports describing AR chromatin interactions have been limited to cell lines, identifying numerous other factors and interacting transcription factors that impact AR chromatin interactions. Do these factors also impact AR cistromics - the genome-wide chromatin-binding landscape of AR - in vivo? Recent technological advances now enable researchers to identify AR chromatin-binding sites and their target genes in human specimens. In this review, we provide an overview of the different factors that influence AR chromatin binding in prostate cancer specimens, which is complemented with knowledge from cell line studies. Finally, we discuss novel perspectives on studying AR cistromics in clinical samples.

Applications of ENCODE data to Systematic Analyses via Data Integration

Large-scale genomic data have been utilized to generate unprecedented biological findings and new hypotheses. To delineate functional elements in the human genome, the Encyclopedia of DNA Elements (ENCODE) project has generated an enormous amount of genomic data, yielding around 7,000 data profiles in different cell and tissue types. In this article, we reviewed the systematic analyses that have integrated ENCODE data with other data sources to reveal new biological insights, ranging from human genome annotation to the identification of new candidate drugs. These analyses demonstrate the critical impact of ENCODE data on basic biology and translational research.

Expanding the repertoire of glucocorticoid receptor target genes by engineering genomic response elements

This study shows that addition of a single transcription factor binding site can be sufficient to convert a gene into a glucocorticoid receptor target.

The glucocorticoid receptor (GR), a hormone-activated transcription factor, binds to a myriad of genomic binding sites yet seems to regulate a much smaller number of genes. Genome-wide analysis of GR binding and gene regulation has shown that the likelihood of GR-dependent regulation increases with decreased distance of its binding to the transcriptional start site of a gene. To test if we can adopt this knowledge to expand the repertoire of GR target genes, we used CRISPR/Cas-mediated homology-directed repair to add a single GR-binding site directly upstream of the transcriptional start site of each of four genes. To our surprise, we found that the addition of a single GR-binding site can be enough to convert a gene into a GR target. The gain of GR-dependent regulation was observed for two of four genes analyzed and coincided with acquired GR binding at the introduced binding site. However, the gene-specific gain of GR-dependent regulation could not be explained by obvious differences in chromatin accessibility between converted genes and their non-converted counterparts. Furthermore, by introducing GR-binding sequences with different nucleotide compositions, we show that activation can be facilitated by distinct sequences without obvious differences in activity between the GR-binding sequence variants we tested. The approach to use genome engineering to build genomic response elements facilitates the generation of cell lines with tailored repertoires of GR-responsive genes and a framework to test and refine our understanding of the cis-regulatory logic of gene regulation by testing if engineered response elements behave as predicted.

Extensive epigenomic integration of the glucocorticoid response in primary human monocytes and in vitro derived macrophages

Glucocorticoid receptor is a transcription factor that is ubiquitously expressed. Glucocorticoids are circadian steroids that regulate a wide range of bodily functions, including immunity. Here we report that synthetic glucocorticoids affect 1035 mRNAs in isolated healthy human blood monocytes but only 165 in the respective six day-old monocyte-derived macrophages. The majority of the glucocorticoid response in monocytes concerns genes that are dynamic upon monocyte to macrophage differentiation, whereby macrophage-like mRNA levels are often reached in monocytes within four hours of treatment. Concomitantly, over 5000 chromosomal H3K27ac regions undergo remodelling, of which 60% involve increased H3K27ac signal. We find that chromosomal glucocorticoid receptor binding sites correlate with positive but not with negative local epigenomic effects. To investigate further we assigned our data to topologically associating domains (TADs). This shows that about 10% of macrophage TADs harbour at least one GR binding site and that half of all the glucocorticoid-induced H3K27ac regions are confined to these TADs. Our analyses are therefore consistent with the notion that TADs naturally accommodate information from sets of distal glucocorticoid response elements.