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 signaling induces transcriptional memory and universally reversible chromatin changes.

Glucocorticoids are stress hormones that elicit cellular responses by binding to the glucocorticoid receptor, a ligand-activated transcription factor. The exposure of cells to this hormone induces wide-spread changes in the chromatin landscape and gene expression. Previous studies have suggested that some of these changes are reversible whereas others persist even when the hormone is no longer around. However, when we examined chromatin accessibility in human airway epithelial cells after hormone washout, we found that the hormone-induced changes were universally reversed after 1 d. Moreover, priming of cells by a previous exposure to hormone, in general, did not alter the transcriptional response to a subsequent encounter of the same cue except for one gene, ZBTB16, that displays transcriptional memory manifesting itself as a more robust transcriptional response upon repeated hormone stimulation. Single-cell analysis revealed that the more robust response is driven by a higher probability of primed cells to activate ZBTB16 and by a subset of cells that express the gene at levels that are higher than the induction levels observed for naïve cells.

Androgen and glucocorticoid receptor direct distinct transcriptional programs by receptor-specific and shared DNA binding sites.

The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied them in an equivalent cellular context. Analysis of chromatin and sequence suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the result of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.

CRUP: a comprehensive framework to predict condition-specific regulatory units.

We present the software Condition-specific Regulatory Units Prediction (CRUP) to infer from epigenetic marks a list of regulatory units consisting of dynamically changing enhancers with their target genes. The workflow consists of a novel pre-trained enhancer predictor that can be reliably applied across cell types and species, solely based on histone modification ChIP-seq data. Enhancers are subsequently assigned to different conditions and correlated with gene expression to derive regulatory units. We thoroughly test and then apply CRUP to a rheumatoid arthritis model, identifying enhancer-gene pairs comprising known disease genes as well as new candidate genes.

Glucocorticoid resistance of allogeneic T cells alters the gene expression profile in the inflamed small intestine of mice suffering from acute graft-versus-host disease.

Glucocorticoids (GCs) play an important role in controlling acute graft-versus-host disease (aGvHD), a frequent complication of allogeneic hematopoietic stem cell transplantation. The anti-inflammatory activity of GCs is mainly ascribed to the modulation of T cells and macrophages, for which reason a genetically induced GC resistance of either of these cell types causes aggravated aGvHD. Since only a few genes are currently known that are differentially regulated under these conditions, we analyzed the expression of 54 candidate genes in the inflamed small intestine of mice suffering from aGvHD when either allogeneic T cells or host myeloid cells were GC resistant using a microfluidic dynamic array platform for high-throughput quantitative PCR. The majority of genes categorized as cytokines (e.g. Il2, Il6), chemokines (e.g. Ccl2, Cxcl1), cell surface receptors (e.g. Fasl, Ctla4) and intracellular molecules (e.g. Dusp1, Arg1) were upregulated in mice transplanted with GC resistant allogeneic T cells. Moreover, the expression of several genes linked to energy metabolism (e.g. Glut1) was altered. Surprisingly, mice harboring GC resistant myeloid cells showed almost no changes in gene expression despite their fatal disease course after aGvHD induction. To identify additional genes in the inflamed small intestine that were affected by a GC resistance of allogeneic T cells, we performed an RNAseq analysis, which uncovered more than 500 differentially expressed transcripts (e.g. Cxcr6, Glut3, Otc, Aoc1, Il1r1, Sphk1) that were enriched for biological processes associated with inflammation and tissue disassembly. The changes in gene expression could be confirmed during full-blown disease but hardly any of them in the preclinical phase using high-throughput quantitative PCR. Further analysis of some of these genes revealed a highly selective expression pattern in T cells, intestinal epithelial cells and macrophages, which correlated with their regulation during disease progression. Collectively, we identified an altered gene expression profile caused by GC resistance of transplanted allogeneic T cells, which could help to define new targets for aGvHD therapy.

Pluripotency reprogramming by competent and incompetent POU factors uncovers temporal dependency for Oct4 and Sox2.

Oct4, along with Sox2 and Klf4 (SK), can induce pluripotency but structurally similar factors like Oct6 cannot. To decode why Oct4 has this unique ability, we compare Oct4-binding, accessibility patterns and transcriptional waves with Oct6 and an Oct4 mutant defective in the dimerization with Sox2 (Oct4defSox2). We find that initial silencing of the somatic program proceeds indistinguishably with or without Oct4. Oct6 mitigates the mesenchymal-to-epithelial transition and derails reprogramming. These effects are a consequence of differences in genome-wide binding, as the early binding profile of Oct4defSox2 resembles Oct4, whilst Oct6 does not bind pluripotency enhancers. Nevertheless, in the Oct6-SK condition many otherwise Oct4-bound locations become accessible but chromatin opening is compromised when Oct4defSox2 occupies these sites. We find that Sox2 predominantly facilitates chromatin opening, whilst Oct4 serves an accessory role. Formation of Oct4/Sox2 heterodimers is essential for pluripotency establishment; however, reliance on Oct4/Sox2 heterodimers declines during pluripotency maintenance.

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

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.

Synthetic STARR-seq reveals how DNA shape and sequence modulate transcriptional output and noise.

The binding of transcription factors to short recognition sequences plays a pivotal role in controlling the expression of genes. The sequence and shape characteristics of binding sites influence DNA binding specificity and have also been implicated in modulating the activity of transcription factors downstream of binding. To quantitatively assess the transcriptional activity of tens of thousands of designed synthetic sites in parallel, we developed a synthetic version of STARR-seq (synSTARR-seq). We used the approach to systematically analyze how variations in the recognition sequence of the glucocorticoid receptor (GR) affect transcriptional regulation. Our approach resulted in the identification of a novel highly active functional GR binding sequence and revealed that sequence variation both within and flanking GR's core binding site can modulate GR activity without apparent changes in DNA binding affinity. Notably, we found that the sequence composition of variants with similar activity profiles was highly diverse. In contrast, groups of variants with similar activity profiles showed specific DNA shape characteristics indicating that DNA shape may be a better predictor of activity than DNA sequence. Finally, using single cell experiments with individual enhancer variants, we obtained clues indicating that the architecture of the response element can independently tune expression mean and cell-to cell variability in gene expression (noise). Together, our studies establish synSTARR as a powerful method to systematically study how DNA sequence and shape modulate transcriptional output and noise.

Loss of the Hematopoietic Stem Cell Factor GATA2 in the Osteogenic Lineage Impairs Trabecularization and Mechanical Strength of Bone.

The transcription factor GATA2 is required for expansion and differentiation of hematopoietic stem cells (HSCs). In mesenchymal stem cells (MSCs), GATA2 blocks adipogenesis, but its biological relevance and underlying genomic events are unknown. We report a dual function of GATA2 in bone homeostasis. GATA2 in MSCs binds near genes involved in skeletal system development and colocalizes with motifs for FOX and HOX transcription factors, known regulators of skeletal development. Ectopic GATA2 blocks osteoblastogenesis by interfering with SMAD1/5/8 activation. MSC-specific deletion of GATA2 in mice increases the numbers and differentiation capacity of bone-derived precursors, resulting in elevated bone formation. Surprisingly, MSC-specific GATA2 deficiency impairs the trabecularization and mechanical strength of bone, involving reduced MSC expression of the osteoclast inhibitor osteoprotegerin and increased osteoclast numbers. Thus, GATA2 affects bone turnover via MSC-autonomous and indirect effects. By regulating bone trabecularization, GATA2 expression in the osteogenic lineage may contribute to the anatomical and cellular microenvironment of the HSC niche required for hematopoiesis.

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

Genomic binding of transcription factors, like the glucocorticoid receptor (GR), is linked to the regulation of genes. However, as we show here, GR binding is a poor predictor of GR-dependent gene regulation even when taking the 3D organization of the genome into account. To connect GR binding sites to the regulation of genes in the endogenous genomic context, we turned to genome editing. By deleting GR binding sites, individually or in combination, we uncovered how cooperative interactions between binding sites contribute to the regulation of genes. Specifically, for the GR target gene GILZ, we show that the simultaneous presence of a cluster of GR binding sites is required for the activity of an individual enhancer and that the GR-dependent regulation of GILZ depends on multiple GR-bound enhancers. Further, by deleting GR binding sites that are shared between different cell types, we show how cell type-specific genome organization and enhancer-blocking can result in cell type-specific wiring of promoter-enhancer contacts. This rewiring allows an individual GR binding site shared between different cell types to direct the expression of distinct transcripts and thereby contributes to the cell type-specific consequences of glucocorticoid signaling.

Identification and characterization of BATF3 as a context-specific coactivator of the glucocorticoid receptor.

The ability of the glucocorticoid receptor (GR) to regulate the transcriptional output of genes relies on its interactions with transcriptional coregulators. However, which coregulators are required for GR-dependent activation is context-dependent and can be influenced by the sequence of the DNA bound by GR and by the nature of the GR isoform responsible for the regulation of a gene. Here, we screened for GR-interacting proteins for which the interaction signal differed between two GR isoforms GRα and GRγ. These isoforms diverge by a single amino acid insertion in a domain, the lever arm, which adopts DNA sequence-specific conformations. We identify Basic Leucine Zipper ATF-Like Transcription Factor 3 (BATF3), an AP-1 family transcription factor, as a GR coregulator whose interaction with GR is modulated by the lever arm. Further, a combination of experiments uncovered that BATF3 acts as a gene-specific coactivator of GR whose coactivator potency is influenced by the sequence of the GR binding site. Together, our findings suggest that GR isoform and the sequence of GR binding site influence the interaction of GR with BATF3, which might direct the assembly of gene-specific regulatory complexes to fine-tune the expression of individual GR target genes.

Airway Epithelial Cells Are Crucial Targets of Glucocorticoids in a Mouse Model of Allergic Asthma.

Although glucocorticoids (GCs) are a mainstay in the clinical management of asthma, the target cells that mediate their therapeutic effects are unknown. Contrary to our expectation, we found that GC receptor (GR) expression in immune cells was dispensable for successful therapy of allergic airway inflammation (AAI) with dexamethasone. Instead, GC treatment was compromised in mice expressing a defective GR in the nonhematopoietic compartment or selectively lacking the GR in airway epithelial cells. Further, we found that an intact GR dimerization interface was a prerequisite for the suppression of AAI and airway hyperresponsiveness by GCs. Our observation that the ability of dexamethasone to modulate gene expression in airway epithelial cells coincided with its potency to resolve AAI supports a crucial role for transcriptional regulation by the GR in this cell type. Taken together, we identified an unknown mode of GC action in the treatment of allergic asthma that might help to develop more specific therapies in the future.

Role of the chromatin landscape and sequence in determining cell type-specific genomic glucocorticoid receptor binding and gene regulation.

The genomic loci bound by the glucocorticoid receptor (GR), a hormone-activated transcription factor, show little overlap between cell types. To study the role of chromatin and sequence in specifying where GR binds, we used Bayesian modeling within the universe of accessible chromatin. Taken together, our results uncovered that although GR preferentially binds accessible chromatin, its binding is biased against accessible chromatin located at promoter regions. This bias can only be explained partially by the presence of fewer GR recognition sequences, arguing for the existence of additional mechanisms that interfere with GR binding at promoters. Therefore, we tested the role of H3K9ac, the chromatin feature with the strongest negative association with GR binding, but found that this correlation does not reflect a causative link. Finally, we find a higher percentage of promoter-proximal GR binding for genes regulated by GR across cell types than for cell type-specific target genes. Given that GR almost exclusively binds accessible chromatin, we propose that cell type-specific regulation by GR preferentially occurs via distal enhancers, whose chromatin accessibility is typically cell type-specific, whereas ubiquitous target gene regulation is more likely to result from binding to promoter regions, which are often accessible regardless of cell type examined.

Sequences flanking the core-binding site modulate glucocorticoid receptor structure and activity.

The glucocorticoid receptor (GR) binds as a homodimer to genomic response elements, which have particular sequence and shape characteristics. Here we show that the nucleotides directly flanking the core-binding site, differ depending on the strength of GR-dependent activation of nearby genes. Our study indicates that these flanking nucleotides change the three-dimensional structure of the DNA-binding site, the DNA-binding domain of GR and the quaternary structure of the dimeric complex. Functional studies in a defined genomic context show that sequence-induced changes in GR activity cannot be explained by differences in GR occupancy. Rather, mutating the dimerization interface mitigates DNA-induced changes in both activity and structure, arguing for a role of DNA-induced structural changes in modulating GR activity. Together, our study shows that DNA sequence identity of genomic binding sites modulates GR activity downstream of binding, which may play a role in achieving regulatory specificity towards individual target genes.

The caffeine-binding adenosine A2A receptor induces age-like HPA-axis dysfunction by targeting glucocorticoid receptor function.

Caffeine is associated with procognitive effects in humans by counteracting overactivation of the adenosine A2A receptor (A2AR), which is upregulated in the human forebrain of aged and Alzheimer's disease (AD) patients. We have previously shown that an anti-A2AR therapy reverts age-like memory deficits, by reestablishment of the hypothalamic-pituitary-adrenal (HPA) axis feedback and corticosterone circadian levels. These observations suggest that A2AR over-activation and glucocorticoid dysfunction are key events in age-related hippocampal deficits; but their direct connection has never been explored. We now show that inducing A2AR overexpression in an aging-like profile is sufficient to trigger HPA-axis dysfunction, namely loss of plasmatic corticosterone circadian oscillation, and promotes reduction of GR hippocampal levels. The synaptic plasticity and memory deficits triggered by GR in the hippocampus are amplified by A2AR over-activation and were rescued by anti-A2AR therapy; finally, we demonstrate that A2AR act on GR nuclear translocation and GR-dependent transcriptional regulation. We provide the first demonstration that A2AR is a major regulator of GR function and that this functional interconnection may be a trigger to age-related memory deficits. This supports the idea that the procognitive effects of A2AR antagonists, namely caffeine, on Alzheimer's and age-related cognitive impairments may rely on its ability to modulate GR actions.

Identification and characterization of DNA sequences that prevent glucocorticoid receptor binding to nearby response elements.

Out of the myriad of potential DNA binding sites of the glucocorticoid receptor (GR) found in the human genome, only a cell-type specific minority is actually bound, indicating that the presence of a recognition sequence alone is insufficient to specify where GR binds. Cooperative interactions with other transcription factors (TFs) are known to contribute to binding specificity. Here, we reasoned that sequence signals preventing GR recruitment to certain loci provide an alternative means to confer specificity. Motif analyses uncovered candidate Negative Regulatory Sequences (NRSs) that interfere with genomic GR binding. Subsequent functional analyses demonstrated that NRSs indeed prevent GR binding to nearby response elements. We show that NRS activity is conserved across species, found in most tissues and that they also interfere with the genomic binding of other TFs. Interestingly, the effects of NRSs appear not to be a simple consequence of changes in chromatin accessibility. Instead, we find that NRSs interact with proteins found at sub-nuclear structures called paraspeckles and that these proteins might mediate the repressive effects of NRSs. Together, our studies suggest that the joint influence of positive and negative sequence signals partition the genome into regions where GR can bind and those where it cannot.

Mechanisms of Glucocorticoid-Regulated Gene Transcription.

One fascinating aspect of glucocorticoid signaling is their broad range of physiological and pharmacological effects. These effects are at least in part a consequence of transcriptional regulation by the glucocorticoid receptor (GR). Activation of GR by glucocorticoids results in tissue-specific changes in gene expression levels with some genes being activated whereas others are repressed. This raises two questions: First, how does GR regulate different subsets of target genes in different tissues? And second, how can GR both activate and repress the expression of genes?To answer these questions, this chapter will describe the function of the various "components" and how they cooperate to mediate the transcriptional responses to glucocorticoids. The first "component" is GR itself. The second "component" is the chromatin and its role in specifying where in the genome GR binds. Binding to the genome however is just the first step in regulating the expression of genes and transcriptional regulation by GR depends on the recruitment of coregulator proteins that either directly or indirectly influence the recruitment and or activity of RNA polymerase II. Ultimately, the integration of inputs including GR isoform, DNA sequence, chromatin and cooperation with coregulators determines which genes are regulated and the direction of their regulation.

ChIP-exo signal associated with DNA-binding motifs provides insight into the genomic binding of the glucocorticoid receptor and cooperating transcription factors.

The classical DNA recognition sequence of the glucocorticoid receptor (GR) appears to be present at only a fraction of bound genomic regions. To identify sequences responsible for recruitment of this transcription factor (TF) to individual loci, we turned to the high-resolution ChIP-exo approach. We exploited this signal by determining footprint profiles of TF binding at single-base-pair resolution using ExoProfiler, a computational pipeline based on DNA binding motifs. When applied to our GR and the few available public ChIP-exo data sets, we find that ChIP-exo footprints are protein- and recognition sequence-specific signatures of genomic TF association. Furthermore, we show that ChIP-exo captures information about TFs other than the one directly targeted by the antibody in the ChIP procedure. Consequently, the shape of the ChIP-exo footprint can be used to discriminate between direct and indirect (tethering to other DNA-bound proteins) DNA association of GR. Together, our findings indicate that the absence of classical recognition sequences can be explained by direct GR binding to a broader spectrum of sequences than previously known, either as a homodimer or as a heterodimer binding together with a member of the ETS or TEAD families of TFs, or alternatively by indirect recruitment via FOX or STAT proteins. ChIP-exo footprints also bring structural insights and locate DNA:protein cross-link points that are compatible with crystal structures of the studied TFs. Overall, our generically applicable footprint-based approach uncovers new structural and functional insights into the diverse ways of genomic cooperation and association of TFs.