High-throughput screening of glucocorticoid-induced enhancer activity reveals mechanisms of stress-related psychiatric disorders.

Exposure to stressful life events increases the risk for psychiatric disorders. Mechanistic insight into the genetic factors moderating the impact of stress can increase our understanding of disease processes. Here, we test 3,662 single nucleotide polymorphisms (SNPs) from preselected expression quantitative trait loci in massively parallel reporter assays to identify genetic variants that modulate the activity of regulatory elements sensitive to glucocorticoids, important mediators of the stress response. Of the tested SNP sequences, 547 were located in glucocorticoid-responsive regulatory elements of which 233 showed allele-dependent activity. Transcripts regulated by these functional variants were enriched for those differentially expressed in psychiatric disorders in the postmortem brain. Phenome-wide Mendelian randomization analysis in 4,439 phenotypes revealed potentially causal associations specifically in neurobehavioral traits, including major depression and other psychiatric disorders. Finally, a functional gene score derived from these variants was significantly associated with differences in the physiological stress response, suggesting that these variants may alter disease risk by moderating the individual set point of the stress response.

Expression of Lineage Transcription Factors Identifies Differences in Transition States of Induced Human Oligodendrocyte Differentiation.

Oligodendrocytes (OLs) are critical for myelination and are implicated in several brain disorders. Directed differentiation of human-induced OLs (iOLs) from pluripotent stem cells can be achieved by forced expression of different combinations of the transcription factors SOX10 (S), OLIG2 (O), and NKX6.2 (N). Here, we applied quantitative image analysis and single-cell transcriptomics to compare different transcription factor (TF) combinations for their efficacy towards robust OL lineage conversion. Compared with S alone, the combination of SON increases the number of iOLs and generates iOLs with a more complex morphology and higher expression levels of myelin-marker genes. RNA velocity analysis of individual cells reveals that S generates a population of oligodendrocyte-precursor cells (OPCs) that appear to be more immature than those generated by SON and to display distinct molecular properties. Our work highlights that TFs for generating iOPCs or iOLs should be chosen depending on the intended application or research question, and that SON might be beneficial to study more mature iOLs while S might be better suited to investigate iOPC biology.

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.

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.

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.

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.

EBF1 binds to EBNA2 and promotes the assembly of EBNA2 chromatin complexes in B cells.

Epstein-Barr virus (EBV) infection converts resting human B cells into permanently proliferating lymphoblastoid cell lines (LCLs). The Epstein-Barr virus nuclear antigen 2 (EBNA2) plays a key role in this process. It preferentially binds to B cell enhancers and establishes a specific viral and cellular gene expression program in LCLs. The cellular DNA binding factor CBF1/CSL serves as a sequence specific chromatin anchor for EBNA2. The ubiquitous expression of this highly conserved protein raises the question whether additional cellular factors might determine EBNA2 chromatin binding selectively in B cells. Here we used CBF1 deficient B cells to identify cellular genes up or downregulated by EBNA2 as well as CBF1 independent EBNA2 chromatin binding sites. Apparently, CBF1 independent EBNA2 target genes and chromatin binding sites can be identified but are less frequent than CBF1 dependent EBNA2 functions. CBF1 independent EBNA2 binding sites are highly enriched for EBF1 binding motifs. We show that EBNA2 binds to EBF1 via its N-terminal domain. CBF1 proficient and deficient B cells require EBF1 to bind to CBF1 independent binding sites. Our results identify EBF1 as a co-factor of EBNA2 which conveys B cell specificity to EBNA2.

Inactivation of intergenic enhancers by EBNA3A initiates and maintains polycomb signatures across a chromatin domain encoding CXCL10 and CXCL9.

Epstein-Barr virus (EBV) causes a persistent infection in human B cells by establishing specific transcription programs to control B cell activation and differentiation. Transcriptional reprogramming of EBV infected B cells is predominantly driven by the action of EBV nuclear antigens, among them the transcriptional repressor EBNA3A. By comparing gene expression profiles of wt and EBNA3A negative EBV infected B cells, we have previously identified a broad array of cellular genes controlled by EBNA3A. We now find that genes repressed by EBNA3A in these cells are significantly enriched for the repressive histone mark H3K27me3, which is installed by Polycomb group (PcG) proteins. This PcG-controlled subset of genes also carries H3K27me3 marks in a variety of other tissues, suggesting that the commitment to PcG silencing is an intrinsic feature of these gene loci that can be used by EBNA3A. In addition, EBNA3A targets frequently reside in co-regulated gene clusters. To study the mechanism of gene repression by EBNA3A and to evaluate the relative contribution of PcG proteins during this process, we have selected the genomic neighbors CXCL10 and CXCL9 as a model for co-repressed and PcG-controlled genes. We show that EBNA3A binds to CBF1 occupied intergenic enhancers located between CXCL10 and CXCL9 and displaces the transactivator EBNA2. This impairs enhancer activity, resulting in a rapid transcriptional shut-down of both genes in a CBF1-dependent manner and initiation of a delayed gain of H3K27me3 marks covering an extended chromatin domain. H3K27me3 marks increase gradually and are maintained by EBNA3A. Our study provides direct evidence that repression by EBNA3A requires CBF1 and that EBNA3A and EBNA2 compete for access to CBF1 at identical genomic sites. Most importantly, our results demonstrate that transcriptional silencing by EBNA3A precedes the appearance of repressive PcG marks and indicate that both events are triggered by loss of enhancer activity.