Enhancer selectivity in space and time: from enhancer-promoter interactions to promoter activation.

The primary regulators of metazoan gene expression are enhancers, originally functionally defined as DNA sequences that can activate transcription at promoters in an orientation-independent and distance-independent manner. Despite being crucial for gene regulation in animals, what mechanisms underlie enhancer selectivity for promoters, and more fundamentally, how enhancers interact with promoters and activate transcription, remain poorly understood. In this Review, we first discuss current models of enhancer-promoter interactions in space and time and how enhancers affect transcription activation. Next, we discuss different mechanisms that mediate enhancer selectivity, including repression, biochemical compatibility and regulation of 3D genome structure. Through 3D polymer simulations, we illustrate how the ability of 3D genome folding mechanisms to mediate enhancer selectivity strongly varies for different enhancer-promoter interaction mechanisms. Finally, we discuss how recent technical advances may provide new insights into mechanisms of enhancer-promoter interactions and how technical biases in methods such as Hi-C and Micro-C and imaging techniques may affect their interpretation.

3D Genome of macaque fetal brain reveals evolutionary innovations during primate corticogenesis.

Elucidating the regulatory mechanisms of human brain evolution is essential to understanding human cognition and mental disorders. We generated multi-omics profiles and constructed a high-resolution map of 3D genome architecture of rhesus macaque during corticogenesis. By comparing the 3D genomes of human, macaque, and mouse brains, we identified many human-specific chromatin structure changes, including 499 topologically associating domains (TADs) and 1,266 chromatin loops. The human-specific loops are significantly enriched in enhancer-enhancer interactions, and the regulated genes show human-specific expression changes in the subplate, a transient zone of the developing brain critical for neural circuit formation and plasticity. Notably, many human-specific sequence changes are located in the human-specific TAD boundaries and loop anchors, which may generate new transcription factor binding sites and chromatin structures in human. Collectively, the presented data highlight the value of comparative 3D genome analyses in dissecting the regulatory mechanisms of brain development and evolution.

SnapShot: The Regulatory Genome.

Enhancers switch genes on and off in response to a variety of intrinsic and external cellular signals. They are the cornerstone of gene regulation and the most pervasive constituents of the regulatory genome. Sequence polymorphisms in enhancer DNAs are a major source of population diversity and predilection to disease. To view this SnapShot, open or download the PDF.

Chromatin Potential Identified by Shared Single-Cell Profiling of RNA and Chromatin.

Cell differentiation and function are regulated across multiple layers of gene regulation, including modulation of gene expression by changes in chromatin accessibility. However, differentiation is an asynchronous process precluding a temporal understanding of regulatory events leading to cell fate commitment. Here we developed simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq), a highly scalable approach for measurement of chromatin accessibility and gene expression in the same single cell, applicable to different tissues. Using 34,774 joint profiles from mouse skin, we develop a computational strategy to identify cis-regulatory interactions and define domains of regulatory chromatin (DORCs) that significantly overlap with super-enhancers. During lineage commitment, chromatin accessibility at DORCs precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for lineage commitment. We computationally infer chromatin potential as a quantitative measure of chromatin lineage-priming and use it to predict cell fate outcomes. SHARE-seq is an extensible platform to study regulatory circuitry across diverse cells in tissues.

Functional Long Non-coding RNAs Evolve from Junk Transcripts.

Transcriptome studies reveal pervasive transcription of complex genomes, such as those of mammals. Despite popular arguments for functionality of most, if not all, of these transcripts, genome-wide analysis of selective constraints indicates that most of the produced RNA are junk. However, junk is not garbage. On the contrary, junk transcripts provide the raw material for the evolution of diverse long non-coding (lnc) RNAs by non-adaptive mechanisms, such as constructive neutral evolution. The generation of many novel functional entities, such as lncRNAs, that fuels organismal complexity does not seem to be driven by strong positive selection. Rather, the weak selection regime that dominates the evolution of most multicellular eukaryotes provides ample material for functional innovation with relatively little adaptation involved.

Comprehensive In Vivo Interrogation Reveals Phenotypic Impact of Human Enhancer Variants.

Establishing causal links between non-coding variants and human phenotypes is an increasing challenge. Here, we introduce a high-throughput mouse reporter assay for assessing the pathogenic potential of human enhancer variants in vivo and examine nearly a thousand variants in an enhancer repeatedly linked to polydactyly. We show that 71% of all rare non-coding variants previously proposed as causal lead to reporter gene expression in a pattern consistent with their pathogenic role. Variants observed to alter enhancer activity were further confirmed to cause polydactyly in knockin mice. We also used combinatorial and single-nucleotide mutagenesis to evaluate the in vivo impact of mutations affecting all positions of the enhancer and identified additional functional substitutions, including potentially pathogenic variants hitherto not observed in humans. Our results uncover the functional consequences of hundreds of mutations in a phenotype-associated enhancer and establish a widely applicable strategy for systematic in vivo evaluation of human enhancer variants.

TCF-1 promotes chromatin interactions across topologically associating domains in T cell progenitors.

The high mobility group (HMG) transcription factor TCF-1 is essential for early T cell development. Although in vitro biochemical assays suggest that HMG proteins can serve as architectural elements in the assembly of higher-order nuclear organization, the contribution of TCF-1 on the control of three-dimensional (3D) genome structures during T cell development remains unknown. Here, we investigated the role of TCF-1 in 3D genome reconfiguration. Using gain- and loss-of-function experiments, we discovered that the co-occupancy of TCF-1 and the architectural protein CTCF altered the structure of topologically associating domains in T cell progenitors, leading to interactions between previously insulated regulatory elements and target genes at late stages of T cell development. The TCF-1-dependent gain in long-range interactions was linked to deposition of active enhancer mark H3K27ac and recruitment of the cohesin-loading factor NIPBL at active enhancers. These data indicate that TCF-1 has a role in controlling global genome organization during T cell development.

A Genome-wide Framework for Mapping Gene Regulation via Cellular Genetic Screens.

Over one million candidate regulatory elements have been identified across the human genome, but nearly all are unvalidated and their target genes uncertain. Approaches based on human genetics are limited in scope to common variants and in resolution by linkage disequilibrium. We present a multiplex, expression quantitative trait locus (eQTL)-inspired framework for mapping enhancer-gene pairs by introducing random combinations of CRISPR/Cas9-mediated perturbations to each of many cells, followed by single-cell RNA sequencing (RNA-seq). Across two experiments, we used dCas9-KRAB to perturb 5,920 candidate enhancers with no strong a priori hypothesis as to their target gene(s), measuring effects by profiling 254,974 single-cell transcriptomes. We identified 664 (470 high-confidence) cis enhancer-gene pairs, which were enriched for specific transcription factors, non-housekeeping status, and genomic and 3D conformational proximity to their target genes. This framework will facilitate the large-scale mapping of enhancer-gene regulatory interactions, a critical yet largely uncharted component of the cis-regulatory landscape of the human genome.

The cis-Regulatory Atlas of the Mouse Immune System.

A complete chart of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of differentiation and function of an organ system. We generated matched epigenome and transcriptome measurements in 86 primary cell types that span the mouse immune system and its differentiation cascades. This breadth of data enable variance components analysis that suggests that genes fall into two distinct classes, controlled by either enhancer- or promoter-driven logic, and multiple regression that connects genes to the enhancers that regulate them. Relating transcription factor (TF) expression to the genome-wide accessibility of their binding motifs classifies them as predominantly openers or closers of local chromatin accessibility, pinpointing specific cis-regulatory elements where binding of given TFs is likely functionally relevant, validated by chromatin immunoprecipitation sequencing (ChIP-seq). Overall, this cis-regulatory atlas provides a trove of information on transcriptional regulation through immune differentiation and a foundational scaffold to define key regulatory events throughout the immunological genome.

Identifying cis Elements for Spatiotemporal Control of Mammalian DNA Replication.

The temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide early-to-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and the necessity of these "early replication control elements" (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

Recurrent SMARCB1 Mutations Reveal a Nucleosome Acidic Patch Interaction Site That Potentiates mSWI/SNF Complex Chromatin Remodeling.

Mammalian switch/sucrose non-fermentable (mSWI/SNF) complexes are multi-component machines that remodel chromatin architecture. Dissection of the subunit- and domain-specific contributions to complex activities is needed to advance mechanistic understanding. Here, we examine the molecular, structural, and genome-wide regulatory consequences of recurrent, single-residue mutations in the putative coiled-coil C-terminal domain (CTD) of the SMARCB1 (BAF47) subunit, which cause the intellectual disability disorder Coffin-Siris syndrome (CSS), and are recurrently found in cancers. We find that the SMARCB1 CTD contains a basic α helix that binds directly to the nucleosome acidic patch and that all CSS-associated mutations disrupt this binding. Furthermore, these mutations abrogate mSWI/SNF-mediated nucleosome remodeling activity and enhancer DNA accessibility without changes in genome-wide complex localization. Finally, heterozygous CSS-associated SMARCB1 mutations result in dominant gene regulatory and morphologic changes during iPSC-neuronal differentiation. These studies unmask an evolutionarily conserved structural role for the SMARCB1 CTD that is perturbed in human disease.

Collaboration between distinct SWI/SNF chromatin remodeling complexes directs enhancer selection and activation of macrophage inflammatory genes.

Macrophages elicit immune responses to pathogens through induction of inflammatory genes. Here, we examined the role of three variants of the SWI/SNF nucleosome remodeling complex-cBAF, ncBAF, and PBAF-in the macrophage response to bacterial endotoxin (lipid A). All three SWI/SNF variants were prebound in macrophages and retargeted to genomic sites undergoing changes in chromatin accessibility following stimulation. Cooperative binding of all three variants associated with de novo chromatin opening and latent enhancer activation. Isolated binding of ncBAF and PBAF, in contrast, associated with activation and repression of active enhancers, respectively. Chemical and genetic perturbations of variant-specific subunits revealed pathway-specific regulation in the activation of lipid A response genes, corresponding to requirement for cBAF and ncBAF in inflammatory and interferon-stimulated gene (ISG) activation, respectively, consistent with differential engagement of SWI/SNF variants by signal-responsive transcription factors. Thus, functional diversity among SWI/SNF variants enables increased regulatory control of innate immune transcriptional programs, with potential for specific therapeutic targeting.

A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification.

Cytokine expression during T cell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17 bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both in vitro and in vivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

A Pan-Cancer Analysis of Enhancer Expression in Nearly 9000 Patient Samples.

The role of enhancers, a key class of non-coding regulatory DNA elements, in cancer development has increasingly been appreciated. Here, we present the detection and characterization of a large number of expressed enhancers in a genome-wide analysis of 8928 tumor samples across 33 cancer types using TCGA RNA-seq data. Compared with matched normal tissues, global enhancer activation was observed in most cancers. Across cancer types, global enhancer activity was positively associated with aneuploidy, but not mutation load, suggesting a hypothesis centered on "chromatin-state" to explain their interplay. Integrating eQTL, mRNA co-expression, and Hi-C data analysis, we developed a computational method to infer causal enhancer-gene interactions, revealing enhancers of clinically actionable genes. Having identified an enhancer ∼140 kb downstream of PD-L1, a major immunotherapy target, we validated it experimentally. This study provides a systematic view of enhancer activity in diverse tumor contexts and suggests the clinical implications of enhancers.

Risk SNP-Mediated Promoter-Enhancer Switching Drives Prostate Cancer through lncRNA PCAT19.

The prostate cancer (PCa) risk-associated SNP rs11672691 is positively associated with aggressive disease at diagnosis. We showed that rs11672691 maps to the promoter of a short isoform of long noncoding RNA PCAT19 (PCAT19-short), which is in the third intron of the long isoform (PCAT19-long). The risk variant is associated with decreased and increased levels of PCAT19-short and PCAT19-long, respectively. Mechanistically, the risk SNP region is bifunctional with both promoter and enhancer activity. The risk variants of rs11672691 and its LD SNP rs887391 decrease binding of transcription factors NKX3.1 and YY1 to the promoter of PCAT19-short, resulting in weaker promoter but stronger enhancer activity that subsequently activates PCAT19-long. PCAT19-long interacts with HNRNPAB to activate a subset of cell-cycle genes associated with PCa progression, thereby promoting PCa tumor growth and metastasis. Taken together, these findings reveal a risk SNP-mediated promoter-enhancer switching mechanism underlying both initiation and progression of aggressive PCa.

Structural Alterations Driving Castration-Resistant Prostate Cancer Revealed by Linked-Read Genome Sequencing.

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

A Somatically Acquired Enhancer of the Androgen Receptor Is a Noncoding Driver in Advanced Prostate Cancer.

Increased androgen receptor (AR) activity drives therapeutic resistance in advanced prostate cancer. The most common resistance mechanism is amplification of this locus presumably targeting the AR gene. Here, we identify and characterize a somatically acquired AR enhancer located 650 kb centromeric to the AR. Systematic perturbation of this enhancer using genome editing decreased proliferation by suppressing AR levels. Insertion of an additional copy of this region sufficed to increase proliferation under low androgen conditions and to decrease sensitivity to enzalutamide. Epigenetic data generated in localized prostate tumors and benign specimens support the notion that this region is a developmental enhancer. Collectively, these observations underscore the importance of epigenomic profiling in primary specimens and the value of deploying genome editing to functionally characterize noncoding elements. More broadly, this work identifies a therapeutic vulnerability for targeting the AR and emphasizes the importance of regulatory elements as highly recurrent oncogenic drivers.

Analysis of Genetically Diverse Macrophages Reveals Local and Domain-wide Mechanisms that Control Transcription Factor Binding and Function.

Non-coding genetic variation is a major driver of phenotypic diversity and allows the investigation of mechanisms that control gene expression. Here, we systematically investigated the effects of >50 million variations from five strains of mice on mRNA, nascent transcription, transcription start sites, and transcription factor binding in resting and activated macrophages. We observed substantial differences associated with distinct molecular pathways. Evaluating genetic variation provided evidence for roles of ∼100 TFs in shaping lineage-determining factor binding. Unexpectedly, a substantial fraction of strain-specific factor binding could not be explained by local mutations. Integration of genomic features with chromatin interaction data provided evidence for hundreds of connected cis-regulatory domains associated with differences in transcription factor binding and gene expression. This system and the >250 datasets establish a substantial new resource for investigation of how genetic variation affects cellular phenotypes.

Small Molecules Co-targeting CKIα and the Transcriptional Kinases CDK7/9 Control AML in Preclinical Models.

CKIα ablation induces p53 activation, and CKIα degradation underlies the therapeutic effect of lenalidomide in a pre-leukemia syndrome. Here we describe the development of CKIα inhibitors, which co-target the transcriptional kinases CDK7 and CDK9, thereby augmenting CKIα-induced p53 activation and its anti-leukemic activity. Oncogene-driving super-enhancers (SEs) are highly sensitive to CDK7/9 inhibition. We identified multiple newly gained SEs in primary mouse acute myeloid leukemia (AML) cells and demonstrate that the inhibitors abolish many SEs and preferentially suppress the transcription elongation of SE-driven oncogenes. We show that blocking CKIα together with CDK7 and/or CDK9 synergistically stabilize p53, deprive leukemia cells of survival and proliferation-maintaining SE-driven oncogenes, and induce apoptosis. Leukemia progenitors are selectively eliminated by the inhibitors, explaining their therapeutic efficacy with preserved hematopoiesis and leukemia cure potential; they eradicate leukemia in MLL-AF9 and Tet2-/-;Flt3ITD AML mouse models and in several patient-derived AML xenograft models, supporting their potential efficacy in curing human leukemia.

Quantitative control of Ets1 dosage by a multi-enhancer hub promotes Th1 cell differentiation and protects from allergic inflammation.

Multi-enhancer hubs are spatial clusters of enhancers present across numerous developmental programs. Here, we studied the functional relevance of these three-dimensional structures in T cell biology. Mathematical modeling identified a highly connected multi-enhancer hub at the Ets1 locus, comprising a noncoding regulatory element that was a hotspot for sequence variation associated with allergic disease in humans. Deletion of this regulatory element in mice revealed that the multi-enhancer connectivity was dispensable for T cell development but required for CD4+ T helper 1 (Th1) differentiation. These mice were protected from Th1-mediated colitis but exhibited overt allergic responses. Mechanistically, the multi-enhancer hub controlled the dosage of Ets1 that was required for CTCF recruitment and assembly of Th1-specific genome topology. Our findings establish a paradigm wherein multi-enhancer hubs control cellular competence to respond to an inductive cue through quantitative control of gene dosage and provide insight into how sequence variation within noncoding elements at the Ets1 locus predisposes individuals to allergic responses.