Direct observation of a condensate effect on super-enhancer controlled gene bursting.

Enhancers are distal DNA elements believed to loop and contact promoters to control gene expression. Recently, we found diffraction-sized transcriptional condensates at genes controlled by clusters of enhancers (super-enhancers). However, a direct function of endogenous condensates in controlling gene expression remains elusive. Here, we develop live-cell super-resolution and multi-color 3D-imaging approaches to investigate putative roles of endogenous condensates in the regulation of super-enhancer controlled gene Sox2. In contrast to enhancer distance, we find instead that the condensate's positional dynamics are a better predictor of gene expression. A basal gene bursting occurs when the condensate is far (>1 µm), but burst size and frequency are enhanced when the condensate moves in proximity (<1 µm). Perturbations of cohesin and local DNA elements do not prevent basal bursting but affect the condensate and its burst enhancement. We propose a three-way kissing model whereby the condensate interacts transiently with gene locus and regulatory DNA elements to control gene bursting.

Linking chromatin acylation mark-defined proteome and genome in living cells.

A generalizable strategy with programmable site specificity for in situ profiling of histone modifications on unperturbed chromatin remains highly desirable but challenging. We herein developed a single-site-resolved multi-omics (SiTomics) strategy for systematic mapping of dynamic modifications and subsequent profiling of chromatinized proteome and genome defined by specific chromatin acylations in living cells. By leveraging the genetic code expansion strategy, our SiTomics toolkit revealed distinct crotonylation (e.g., H3K56cr) and ß-hydroxybutyrylation (e.g., H3K56bhb) upon short chain fatty acids stimulation and established linkages for chromatin acylation mark-defined proteome, genome, and functions. This led to the identification of GLYR1 as a distinct interacting protein in modulating H3K56cr's gene body localization as well as the discovery of an elevated super-enhancer repertoire underlying bhb-mediated chromatin modulations. SiTomics offers a platform technology for elucidating the "metabolites-modification-regulation" axis, which is widely applicable for multi-omics profiling and functional dissection of modifications beyond acylations and proteins beyond histones.

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.

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.

A Phase Separation Model for Transcriptional Control.

Phase-separated multi-molecular assemblies provide a general regulatory mechanism to compartmentalize biochemical reactions within cells. We propose that a phase separation model explains established and recently described features of transcriptional control. These features include the formation of super-enhancers, the sensitivity of super-enhancers to perturbation, the transcriptional bursting patterns of enhancers, and the ability of an enhancer to produce simultaneous activation at multiple genes. This model provides a conceptual framework to further explore principles of gene control in mammals.

Super-Enhancer-Mediated RNA Processing Revealed by Integrative MicroRNA Network Analysis.

Super-enhancers are an emerging subclass of regulatory regions controlling cell identity and disease genes. However, their biological function and impact on miRNA networks are unclear. Here, we report that super-enhancers drive the biogenesis of master miRNAs crucial for cell identity by enhancing both transcription and Drosha/DGCR8-mediated primary miRNA (pri-miRNA) processing. Super-enhancers, together with broad H3K4me3 domains, shape a tissue-specific and evolutionarily conserved atlas of miRNA expression and function. CRISPR/Cas9 genomics revealed that super-enhancer constituents act cooperatively and facilitate Drosha/DGCR8 recruitment and pri-miRNA processing to boost cell-specific miRNA production. The BET-bromodomain inhibitor JQ1 preferentially inhibits super-enhancer-directed cotranscriptional pri-miRNA processing. Furthermore, super-enhancers are characterized by pervasive interaction with DGCR8/Drosha and DGCR8/Drosha-regulated mRNA stability control, suggesting unique RNA regulation at super-enhancers. Finally, super-enhancers mark multiple miRNAs associated with cancer hallmarks. This study presents principles underlying miRNA biology in health and disease and an unrecognized higher-order property of super-enhancers in RNA processing beyond transcription.

Cis interactions in the Irf8 locus regulate stage-dependent enhancer activation.

Individual elements within a superenhancer can act in a cooperative or temporal manner, but the underlying mechanisms remain obscure. We recently identified an Irf8 superenhancer, within which different elements act at distinct stages of type 1 classical dendritic cell (cDC1) development. The +41-kb Irf8 enhancer is required for pre-cDC1 specification, while the +32-kb Irf8 enhancer acts to support subsequent cDC1 maturation. Here, we found that compound heterozygous Δ32/Δ41 mice, lacking the +32- and +41-kb enhancers on different chromosomes, show normal pre-cDC1 specification but, surprisingly, completely lack mature cDC1 development, suggesting cis dependence of the +32-kb enhancer on the +41-kb enhancer. Transcription of the +32-kb Irf8 enhancer-associated long noncoding RNA (lncRNA) Gm39266 is also dependent on the +41-kb enhancer. However, cDC1 development in mice remained intact when Gm39266 transcripts were eliminated by CRISPR/Cas9-mediated deletion of lncRNA promoters and when transcription across the +32-kb enhancer was blocked by premature polyadenylation. We showed that chromatin accessibility and BATF3 binding at the +32-kb enhancer were dependent on a functional +41-kb enhancer located in cis Thus, the +41-kb Irf8 enhancer controls the subsequent activation of the +32-kb Irf8 enhancer in a manner that is independent of associated lncRNA transcription.

Regulatory T Cell-Specific Epigenomic Region Variants Are a Key Determinant of Susceptibility to Common Autoimmune Diseases.

The contribution of FOXP3-expressing naturally occurring regulatory T (Treg) cells to common polygenic autoimmune diseases remains ambiguous. Here, we characterized genome-wide epigenetic profiles (CpG methylation and histone modifications) of human Treg and conventional T (Tconv) cells in naive and activated states. We found that single-nucleotide polymorphisms (SNPs) associated with common autoimmune diseases were predominantly enriched in CpG demethylated regions (DRs) specifically present in naive Treg cells but much less enriched in activation-induced DRs common in Tconv and Treg cells. Naive Treg cell-specific DRs were largely included in Treg cell-specific super-enhancers and closely associated with transcription and other epigenetic changes in naive and effector Treg cells. Thus, naive Treg cell-specific CpG hypomethylation had a key role in controlling Treg cell-specific gene transcription and epigenetic modification. The results suggest possible contribution of altered function or development of natural Treg cells to the susceptibility to common autoimmune diseases.

Temporal dissection of an enhancer cluster reveals distinct temporal and functional contributions of individual elements.

Many genes are regulated by multiple enhancers that often simultaneously activate their target gene. However, how individual enhancers collaborate to activate transcription is not well understood. Here, we dissect the functions and interdependencies of five enhancer elements that together activate Fgf5 expression during exit from naive murine pluripotency. Four intergenic elements form a super-enhancer, and most of the elements contribute to Fgf5 induction at distinct time points. A fifth, poised enhancer located in the first intron contributes to Fgf5 expression at every time point by amplifying overall Fgf5 expression levels. Despite low individual enhancer activity, together these elements strongly induce Fgf5 expression in a super-additive fashion that involves strong accumulation of RNA polymerase II at the intronic enhancer. Finally, we observe a strong anti-correlation between RNA polymerase II levels at enhancers and their distance to the closest promoter, and we identify candidate elements with properties similar to the intronic enhancer.

Enhancers are activated by p300/CBP activity-dependent PIC assembly, RNAPII recruitment, and pause release.

The metazoan-specific acetyltransferase p300/CBP is involved in activating signal-induced, enhancer-mediated transcription of cell-type-specific genes. However, the global kinetics and mechanisms of p300/CBP activity-dependent transcription activation remain poorly understood. We performed genome-wide, time-resolved analyses to show that enhancers and super-enhancers are dynamically activated through p300/CBP-catalyzed acetylation, deactivated by the opposing deacetylase activity, and kinetic acetylation directly contributes to maintaining cell identity at very rapid (minutes) timescales. The acetyltransferase activity is dispensable for the recruitment of p300/CBP and transcription factors but essential for promoting the recruitment of TFIID and RNAPII at virtually all enhancers and enhancer-regulated genes. This identifies pre-initiation complex assembly as a dynamically controlled step in the transcription cycle and reveals p300/CBP-catalyzed acetylation as the signal that specifically promotes transcription initiation at enhancer-regulated genes. We propose that p300/CBP activity uses a "recruit-and-release" mechanism to simultaneously promote RNAPII recruitment and pause release and thereby enables kinetic activation of enhancer-mediated transcription.

Systematic functional characterization of antisense eRNA of protocadherin α composite enhancer.

Enhancers generate bidirectional noncoding enhancer RNAs (eRNAs) that may regulate gene expression. At present, the eRNA function remains enigmatic. Here, we report a 5' capped antisense eRNA PEARL (Pcdh eRNA associated with R-loop formation) that is transcribed from the protocadherin (Pcdh) α HS5-1 enhancer region. Through loss- and gain-of-function experiments with CRISPR/Cas9 DNA fragment editing, CRISPRi, and CRISPRa, as well as locked nucleic acid strategies, in conjunction with ChIRP, MeDIP, DRIP, QHR-4C, and HiChIP experiments, we found that PEARL regulates Pcdhα gene expression by forming local RNA-DNA duplexes (R-loops) in situ within the HS5-1 enhancer region to promote long-distance chromatin interactions between distal enhancers and target promoters. In particular, increased levels of eRNA PEARL via perturbing transcription elongation factor SPT6 lead to strengthened local three-dimensional chromatin organization within the Pcdh superTAD. These findings have important implications regarding molecular mechanisms by which the HS5-1 enhancer regulates stochastic Pcdhα promoter choice in single cells in the brain.

Evolutionary conservation of VSX2 super-enhancer modules in retinal development.

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. We recently identified developmental stage- and cell type-specific modules within the murine Vsx2 SE. Here, we show that the human VSX2 SE modules have similar developmental stage- and cell type-specific activity in reporter gene assays. By inserting the human sequence of one VSX2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the function of these SE modules during human retinal development, we deleted individual modules in human embryonic stem cells and generated retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module led to disruptions in bipolar neuron development. This prototypical SE serves as a model for understanding developmental stage- and cell type-specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

SP1 undergoes phase separation and activates RGS20 expression through super-enhancers to promote lung adenocarcinoma progression.

Lung adenocarcinoma (LUAD) is the leading cause of cancer-related death worldwide, but the underlying molecular mechanisms remain largely unclear. The transcription factor (TF) specificity protein 1 (SP1) plays a crucial role in the development of various cancers, including LUAD. Recent studies have indicated that master TFs may form phase-separated macromolecular condensates to promote super-enhancer (SE) assembly and oncogene expression. In this study, we demonstrated that SP1 undergoes phase separation and that its zinc finger 3 in the DNA-binding domain is essential for this process. Through Cleavage Under Targets & Release Using Nuclease (CUT&RUN) using antibodies against SP1 and H3K27ac, we found a significant correlation between SP1 enrichment and SE elements, identified the regulator of the G protein signaling 20 (RGS20) gene as the most likely target regulated by SP1 through SE mechanisms, and verified this finding using different approaches. The oncogenic activity of SP1 relies on its phase separation ability and RGS20 gene activation, which can be abolished by glycogen synthase kinase J4 (GSK-J4), a demethylase inhibitor. Together, our findings provide evidence that SP1 regulates its target oncogene expression through phase separation and SE mechanisms, thereby promoting LUAD cell progression. This study also revealed an innovative target for LUAD therapies through intervening in SP1-mediated SE formation.

Dedifferentiation by adenovirus E1A due to inactivation of Hippo pathway effectors YAP and TAZ.

Adenovirus transformed cells have a dedifferentiated phenotype. Eliminating E1A in transformed human embryonic kidney cells derepressed ∼2600 genes, generating a gene expression profile closely resembling mesenchymal stem cells (MSCs). This was associated with a dramatic change in cell morphology from one with scant cytoplasm and a globular nucleus to one with increased cytoplasm, extensive actin stress fibers, and actomyosin-dependent flattening against the substratum. E1A-induced hypoacetylation at histone H3 Lys27 and Lys18 (H3K27/18) was reversed. Most of the increase in H3K27/18ac was in enhancers near TEAD transcription factors bound by Hippo signaling-regulated coactivators YAP and TAZ. E1A causes YAP/TAZ cytoplasmic sequestration. After eliminating E1A, YAP/TAZ were transported into nuclei, where they associated with poised enhancers with DNA-bound TEAD4 and H3K4me1. This activation of YAP/TAZ required RHO family GTPase signaling and caused histone acetylation by p300/CBP, chromatin remodeling, and cohesin loading to establish MSC-associated enhancers and then superenhancers. Consistent results were also observed in primary rat embryo kidney cells, human fibroblasts, and human respiratory tract epithelial cells. These results together with earlier studies suggest that YAP/TAZ function in a developmental checkpoint controlled by signaling from the actin cytoskeleton that prevents differentiation of a progenitor cell until it is in the correct cellular and tissue environment.

The Chromatin Reader ZMYND8 Regulates Igh Enhancers to Promote Immunoglobulin Class Switch Recombination.

Class switch recombination (CSR) is a DNA recombination reaction that diversifies the effector component of antibody responses. CSR is initiated by activation-induced cytidine deaminase (AID), which targets transcriptionally active immunoglobulin heavy chain (Igh) switch donor and acceptor DNA. The 3' Igh super-enhancer, 3' regulatory region (3'RR), is essential for acceptor region transcription, but how this function is regulated is unknown. Here, we identify the chromatin reader ZMYND8 as an essential regulator of the 3'RR. In B cells, ZMYND8 binds promoters and super-enhancers, including the Igh enhancers. ZMYND8 controls the 3'RR activity by modulating the enhancer transcriptional status. In its absence, there is increased 3'RR polymerase loading and decreased acceptor region transcription and CSR. In addition to CSR, ZMYND8 deficiency impairs somatic hypermutation (SHM) of Igh, which is also dependent on the 3'RR. Thus, ZMYND8 controls Igh diversification in mature B lymphocytes by regulating the activity of the 3' Igh super-enhancer.

Dissecting super-enhancer heterogeneity: time to re-examine cancer subtypes?

The heterogeneity of transcriptional regulations by super-enhancers (SEs) is poorly understood in human cancers. Herein, we summarize a bioinformatics workflow for genome-wide SE profiling and identification of subtype-specific SEs and regulatory networks. Dissecting SE heterogeneity provides new insights into cancer biology and alternative therapeutic strategies for cancer precision medicine.

Analysis of super-enhancer using machine learning and its application to medical biology.

The analysis of super-enhancers (SEs) has recently attracted attention in elucidating the molecular mechanisms of cancer and other diseases. SEs are genomic structures that strongly induce gene expression and have been reported to contribute to the overexpression of oncogenes. Because the analysis of SEs and integrated analysis with other data are performed using large amounts of genome-wide data, artificial intelligence technology, with machine learning at its core, has recently begun to be utilized. In promoting precision medicine, it is important to consider information from SEs in addition to genomic data; therefore, machine learning technology is expected to be introduced appropriately in terms of building a robust analysis platform with a high generalization performance. In this review, we explain the history and principles of SE, and the results of SE analysis using state-of-the-art machine learning and integrated analysis with other data are presented to provide a comprehensive understanding of the current status of SE analysis in the field of medical biology. Additionally, we compared the accuracy between existing machine learning methods on the benchmark dataset and attempted to explore the kind of data preprocessing and integration work needed to make the existing algorithms work on the benchmark dataset. Furthermore, we discuss the issues and future directions of current SE analysis.

Multiple genetic variants at the SLC30A8 locus affect local super-enhancer activity and influence pancreatic ß-cell survival and function.

Variants at the SLC30A8 locus are associated with type 2 diabetes (T2D) risk. The lead variant, rs13266634, encodes an amino acid change, Arg325Trp (R325W), at the C-terminus of the secretory granule-enriched zinc transporter, ZnT8. Although this protein-coding variant was previously thought to be the sole driver of T2D risk at this locus, recent studies have provided evidence for lowered expression of SLC30A8 mRNA in protective allele carriers. In the present study, we examined multiple variants that influence SLC30A8 allele-specific expression. Epigenomic mapping has previously identified an islet-selective enhancer cluster at the SLC30A8 locus, hosting multiple T2D risk and cASE associations, which is spatially associated with the SLC30A8 promoter and additional neighboring genes. Here, we show that deletion of variant-bearing enhancer regions using CRISPR-Cas9 in human-derived EndoC-ßH3 cells lowers the expression of SLC30A8 and several neighboring genes and improves glucose-stimulated insulin secretion. While downregulation of SLC30A8 had no effect on beta cell survival, loss of UTP23, RAD21, or MED30 markedly reduced cell viability. Although eQTL or cASE analyses in human islets did not support the association between these additional genes and diabetes risk, the transcriptional regulator JQ1 lowered the expression of multiple genes at the SLC30A8 locus and enhanced stimulated insulin secretion.

NONO enhances mRNA processing of super-enhancer-associated GATA2 and HAND2 genes in neuroblastoma.

High-risk neuroblastoma patients have poor survival rates and require better therapeutic options. High expression of a multifunctional DNA and RNA-binding protein, NONO, in neuroblastoma is associated with poor patient outcome; however, there is little understanding of the mechanism of NONO-dependent oncogenic gene regulatory activity in neuroblastoma. Here, we used cell imaging, biochemical and genome-wide molecular analysis to reveal complex NONO-dependent regulation of gene expression. NONO forms RNA- and DNA-tethered condensates throughout the nucleus and undergoes phase separation in vitro, modulated by nucleic acid binding. CLIP analyses show that NONO mainly binds to the 5' end of pre-mRNAs and modulates pre-mRNA processing, dependent on its RNA-binding activity. NONO regulates super-enhancer-associated genes, including HAND2 and GATA2. Abrogating NONO RNA binding, or phase separation activity, results in decreased expression of HAND2 and GATA2. Thus, future development of agents that target RNA-binding activity of NONO may have therapeutic potential in this cancer context.

Metabolic reprogramming directed by super-enhancers in tumors: An emerging landscape.

Metabolic reprogramming is an essential hallmark of tumors, and metabolic abnormalities are strongly associated with the malignant phenotype of tumor cells. This is closely related to transcriptional dysregulation. Super-enhancers are extremely active cis-regulatory regions in the genome, and can amalgamate a complex set of transcriptional regulatory components that are crucial for establishing tumor cell identity, promoting tumorigenesis, and enhancing aggressiveness. In addition, alterations in metabolic signaling pathways are often accompanied by changes in super-enhancers. Presently, there is a surge in interest in the potential pathogenesis of various tumors through the transcriptional regulation of super-enhancers and oncogenic mutations in super-enhancers. In this review, we summarize the functions of super-enhancers, oncogenic signaling pathways, and tumor metabolic reprogramming. In particular, we focus on the role of the super-enhancer in tumor metabolism and its impact on metabolic reprogramming. This review also discusses the prospects and directions in the field of super-enhancer and metabolic reprogramming.