Enhancers as information integration hubs in development: lessons from genomics

miR-137 confers robustness to the territorial restriction of the neural plate border

The neural plate border (NPB) of vertebrate embryos is segregated from the neural plate (NP) and epidermal regions, and comprises an intermingled group of progenitors with multiple fate potential. Recent studies have shown that, during the gastrula stage, TFAP2A acts as a pioneer factor in remodeling the epigenetic landscape required to activate components of the NPB induction program. Here, we show that chick Tfap2a has two highly conserved binding sites for miR-137, and both display a reciprocal expression pattern at the NPB and NP, respectively. In addition, ectopic miR-137 expression reduced TFAP2A, whereas its functional inhibition expanded their territorial distribution overlapping with PAX7. Furthermore, we demonstrate that loss of the de novo DNA methyltransferase DNMT3A expanded miR-137 expression to the NPB. Bisulfite sequencing revealed a markedly elevated presence of non-canonical CpH methylation within the miR-137 promoter region when comparing NPB and NP samples. Our findings show that miR-137 contributes to the robustness of NPB territorial restriction in vertebrate development.

KMT2C and KMT2D aberrations in breast cancer

KMT2C and KMT2D are histone lysine methyltransferases responsible for the monomethylation of histone 3 lysine 4 (H3K4) residues at gene enhancer sites. KMT2C/D are the most frequently mutated histone methyltransferases (HMTs) in breast cancer, occurring at frequencies of 10-20% collectively. Frequent damaging and truncating somatic mutations indicate a tumour-suppressive role of KMT2C/D in breast oncogenesis. Recent studies using cell lines and mouse models to replicate KMT2C/D loss show that these genes contribute to oestrogen receptor (ER)-driven transcription in ER+ breast cancers through the priming of gene enhancer regions. This review provides an overview of the functions of KMT2C/D and outlines the recent clinical and experimental evidence of the roles of KMT2C and KMT2D in breast cancer development.

PU.1 and BCL11B sequentially cooperate with RUNX1 to anchor mSWI/SNF to poise the T cell effector landscape

Adaptive immunity relies on specialized effector functions elicited by lymphocytes, yet how antigen recognition activates appropriate effector responses through nonspecific signaling intermediates is unclear. Here we examined the role of chromatin priming in specifying the functional outputs of effector T cells and found that most of the cis-regulatory landscape active in effector T cells was poised early in development before the expression of the T cell antigen receptor. We identified two principal mechanisms underpinning this poised landscape: the recruitment of the nucleosome remodeler mammalian SWItch/Sucrose Non-Fermentable (mSWI/SNF) by the transcription factors RUNX1 and PU.1 to establish chromatin accessibility at T effector loci; and a 'relay' whereby the transcription factor BCL11B succeeded PU.1 to maintain occupancy of the chromatin remodeling complex mSWI/SNF together with RUNX1, after PU.1 silencing during lineage commitment. These mechanisms define modes by which T cells acquire the potential to elicit specialized effector functions early in their ontogeny and underscore the importance of integrating extrinsic cues to the developmentally specified intrinsic program.

The role of lipids in genome integrity and pluripotency

Pluripotent stem cells (PSCs), comprising embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs), offer immense potential for regenerative medicine due to their ability to differentiate into all cell types of the adult body. A critical aspect of harnessing this potential is understanding their metabolic requirements during derivation, maintenance, and differentiation in vitro. Traditional culture methods using fetal bovine serum often lead to issues such as heterogeneous cell populations and diminished pluripotency. Although the chemically-defined 2i/LIF medium has provided solutions to some of these challenges, prolonged culturing of these cells, especially female ESCs, raises concerns related to genome integrity. This review discusses the pivotal role of lipids in genome stability and pluripotency of stem cells. Notably, the introduction of lipid-rich albumin, AlbuMAX, into the 2i/LIF culture medium offers a promising avenue for enhancing the genomic stability and pluripotency of cultured ESCs. We further explore the unique characteristics of lipid-induced pluripotent stem cells (LIP-ESCs), emphasizing their potential in regenerative medicine and pluripotency research.

MAX controls meiotic entry in sexually undifferentiated germ cells

Meiosis is a specialized type of cell division that occurs physiologically only in germ cells. We previously demonstrated that MYC-associated factor X (MAX) blocks the ectopic onset of meiosis in embryonic and germline stem cells in culture systems. Here, we investigated the Max gene's role in mouse primordial germ cells. Although Max is generally ubiquitously expressed, we revealed that sexually undifferentiated male and female germ cells had abundant MAX protein because of their higher Max gene expression than somatic cells. Moreover, our data revealed that this high MAX protein level in female germ cells declined significantly around physiological meiotic onset. Max disruption in sexually undifferentiated germ cells led to ectopic and precocious expression of meiosis-related genes, including Meiosin, the gatekeeper of meiotic onset, in both male and female germ cells. However, Max-null male and female germ cells did not complete the entire meiotic process, but stalled during its early stages and were eventually eliminated by apoptosis. Additionally, our meta-analyses identified a regulatory region that supports the high Max expression in sexually undifferentiated male and female germ cells. These results indicate the strong connection between the Max gene and physiological onset of meiosis in vivo through dynamic alteration of its expression.

Pioneer and PRDM transcription factors coordinate bivalent epigenetic states to safeguard cell fate

Pioneer transcription factors (TFs) regulate cell fate by establishing transcriptionally primed and active states. However, cell fate control requires the coordination of both lineage-specific gene activation and repression of alternative-lineage programs, a process that is poorly understood. Here, we demonstrate that the pioneer TF FOXA coordinates with PRDM1 TF to recruit nucleosome remodeling and deacetylation (NuRD) complexes and Polycomb repressive complexes (PRCs), which establish highly occupied, accessible nucleosome conformation with bivalent epigenetic states, thereby preventing precocious and alternative-lineage gene expression during human endoderm differentiation. Similarly, the pioneer TF OCT4 coordinates with PRDM14 to form bivalent enhancers and repress cell differentiation programs in human pluripotent stem cells, suggesting that this may be a common and critical function of pioneer TFs. We propose that pioneer and PRDM TFs coordinate to safeguard cell fate through epigenetic repression mechanisms.

MAE-seq refines regulatory elements across the genome

Proper cell fate determination relies on precise spatial and temporal genome-wide cooperation between regulatory elements (REs) and their targeted genes. However, the lengths of REs defined using different methods vary, which indicates that there is sequence redundancy and that the context of the genome may be unintelligible. We developed a method called MAE-seq (Massive Active Enhancers by Sequencing) to experimentally identify functional REs at a 25-bp scale. In this study, MAE-seq was used to identify 626879, 541617 and 554826 25-bp enhancers in mouse embryonic stem cells (mESCs), C2C12 and HEK 293T, respectively. Using ∼1.6 trillion 25 bp DNA fragments and screening 12 billion cells, we identified 626879 as active enhancers in mESCs as an example. Comparative analysis revealed that most of the histone modification datasets were annotated by MAE-Seq loci. Furthermore, 33.85% (212195) of the identified enhancers were identified as de novo ones with no epigenetic modification. Intriguingly, distinct chromatin states dictate the requirement for dissimilar cofactors in governing novel and known enhancers. Validation results show that these 25-bp sequences could act as a functional unit, which shows identical or similar expression patterns as the previously defined larger elements, Enhanced resolution facilitated the identification of numerous cell-specific enhancers and their accurate annotation as super enhancers. Moreover, we characterized novel elements capable of augmenting gene activity. By integrating with high-resolution Hi-C data, over 55.64% of novel elements may have a distal association with different targeted genes. For example, we found that the Cdh1 gene interacts with one novel and two known REs in mESCs. The biological effects of these interactions were investigated using CRISPR-Cas9, revealing their role in coordinating Cdh1 gene expression and mESC proliferation. Our study presents an experimental approach to refine the REs at 25-bp resolution, advancing the precision of genome annotation and unveiling the underlying genome context. This novel approach not only advances our understanding of gene regulation but also opens avenues for comprehensive exploration of the genomic landscape.

A chromatin-regulated biphasic circuit coordinates IL-1β-mediated inflammation

Inflammation is characterized by a biphasic cycle consisting initially of a proinflammatory phase that is subsequently resolved by anti-inflammatory processes. Interleukin-1β (IL-1β) is a master regulator of proinflammation and is encoded within the same topologically associating domain (TAD) as IL-37, which is an anti-inflammatory cytokine that opposes the function of IL-1β. Within this TAD, we identified a long noncoding RNA called AMANZI, which negatively regulates IL-1β expression and trained immunity through the induction of IL37 transcription. We found that the activation of IL37 occurs through the formation of a dynamic long-range chromatin contact that leads to the temporal delay of anti-inflammatory responses. The common variant rs16944 present in AMANZI augments this regulatory circuit, predisposing individuals to enhanced proinflammation or immunosuppression. Our work illuminates a chromatin-mediated biphasic circuit coordinating expression of IL-1β and IL-37, thereby regulating two functionally opposed states of inflammation from within a single TAD.

Sequential deregulation of histone marks, chromatin accessibility and gene expression in response to PROTAC-induced degradation of ASH2L

The trithorax protein ASH2L is essential for organismal and tissue development. As a subunit of COMPASS/KMT2 complexes, ASH2L is necessary for methylation of histone H3 lysine 4 (H3K4). Mono- and tri-methylation at this site mark active enhancers and promoters, respectively, although the functional relevance of H3K4 methylation is only partially understood. ASH2L has a long half-life, which results in a slow decrease upon knockout. This has made it difficult to define direct consequences. To overcome this limitation, we employed a PROTAC system to rapidly degrade ASH2L and address direct effects. ASH2L loss resulted in inhibition of proliferation of mouse embryo fibroblasts. Shortly after ASH2L degradation H3K4me3 decreased with its half-life varying between promoters. Subsequently, H3K4me1 increased at promoters and decreased at some enhancers. H3K27ac and H3K27me3, histone marks closely linked to H3K4 methylation, were affected with considerable delay. In parallel, chromatin compaction increased at promoters. Of note, nascent gene transcription was not affected early but overall RNA expression was deregulated late after ASH2L loss. Together, these findings suggest that downstream effects are ordered but relatively slow, despite the rapid loss of ASH2L and inactivation of KMT2 complexes. It appears that the systems that control gene transcription are well buffered and strong effects are only beginning to unfold after considerable delay.

Uncovering oligodendrocyte enhancers that control Cnp expression

Oligodendrocytes (OLs) produce myelin sheaths around axons in the central nervous system (CNS). Myelin accelerates the propagation of action potentials along axons and supports the integrity of axons. Impaired myelination has been linked to neurological and neuropsychiatric disorders. As a major component of CNS myelin, 2',3'-cyclic nucleotide 3'-phosphodiesterase (CNP) plays an indispensable role in the axon-supportive function of myelin. Notably, this function requires a high-level expression of CNP in OLs, as evidenced by downregulated expression of CNP in mental disorders and animal models. Little is known about how CNP expression is regulated in OLs. Especially, OL enhancers that govern CNP remain elusive. We have recently developed a powerful method that links OL enhancers to target genes in a principled manner. Here, we applied it to Cnp, uncovering two OL enhancers for it (termed Cnp-E1 and Cnp-E2). Epigenome editing analysis revealed that Cnp-E1 and Cnp-E2 are dedicated to Cnp. ATAC-seq and ChIP-seq data show that Cnp-E1 and Cnp-E2 are conserved OL-specific enhancers. Single cell multi-omics data that jointly profile gene expression and chromatin accessibility suggest that Cnp-E2 plays an important role in Cnp expression in the early stage of OL differentiation while Cnp-E1 sustains it in mature OLs.

Defining a cohort of anemia-activated cis elements reveals a mechanism promoting erythroid precursor function

Acute anemia elicits broad transcriptional changes in erythroid progenitors and precursors. We previously discovered a cis-regulatory transcriptional enhancer at the sterile alpha motif domain-14 enhancer locus (S14E), defined by a CANNTG-spacer-AGATAA composite motif and occupied by GATA1 and TAL1 transcription factors, is required for survival in severe anemia. However, S14E is only 1 of dozens of anemia-activated genes containing similar motifs. In a mouse model of acute anemia, we identified populations of expanding erythroid precursors, which increased expression of genes that contain S14E-like cis elements. We reveal that several S14E-like cis elements provide important transcriptional control of newly identified anemia-inducing genes, including the Ssx-2 interacting protein (Ssx2ip). Ssx2ip expression was determined to play an important role in erythroid progenitor/precursor cell activities, cell cycle regulation, and cell proliferation. Over a weeklong course of acute anemia recovery, we observed that erythroid gene activation mediated by S14E-like cis elements occurs during a phase coincident with low hematocrit and high progenitor activities, with distinct transcriptional programs activated at earlier and later time points. Our results define a genome-wide mechanism in which S14E-like enhancers control transcriptional responses during erythroid regeneration. These findings provide a framework to understand anemia-specific transcriptional mechanisms, ineffective erythropoiesis, anemia recovery, and phenotypic variability within human populations.

Role of Epigenetics and Metabolomics in Predicting Endothelial Dysfunction in Type 2 Diabetes

Type 2 diabetes (T2D) is a worldwide health problem and cardiovascular disease (CVD) is a leading cause of morbidity and mortality in T2D patients, making the prevention of CVD onset a major priority. It is therefore crucial to optimize diagnosis and treatment to reduce this burden. Endothelial dysfunction is one of the most important prognostic factors for CVD progression, thus novel approaches to identify the early phase of endothelial dysfunction may lead to specific preventive measures to reduce the occurrence of CVD. Nowadays, multiomics approaches have provided unprecedented opportunities to stratify T2D patients into endotypes, improve therapeutic treatment and outcome and amend the survival prediction. Among omics strategies, epigenetics and metabolomics are gaining increasing interest. Recently, a dynamic correlation between metabolic pathways and gene expression through chromatin remodeling, such as DNA methylation, has emerged, indicating new perspectives on the regulatory networks impacting cellular processes. Thus, a better understanding of epigenetic-metabolite relationships can provide insight into the physiological processes altered early in the endothelium that ultimately head to disease development. Here, recent studies on epigenetics and metabolomics related to CVD prevention potentially useful to identify disease biomarkers, as well as new therapies hopefully targeting the early phase of endothelial dysfunction are highlighted.

Transposons contribute to the functional diversification of the head, gut, and ovary transcriptomes across Drosophila natural strains

Transcriptomes are dynamic, with cells, tissues, and body parts expressing particular sets of transcripts. Transposable elements (TEs) are a known source of transcriptome diversity; however, studies often focus on a particular type of chimeric transcript, analyze single body parts or cell types, or are based on incomplete TE annotations from a single reference genome. In this work, we have implemented a method based on de novo transcriptome assembly that minimizes the potential sources of errors while identifying a comprehensive set of gene-TE chimeras. We applied this method to the head, gut, and ovary dissected from five Drosophila melanogaster natural strains, with individual reference genomes available. We found that ∼19% of body part-specific transcripts are gene-TE chimeras. Overall, chimeric transcripts contribute a mean of 43% to the total gene expression, and they provide protein domains for DNA binding, catalytic activity, and DNA polymerase activity. Our comprehensive data set is a rich resource for follow-up analysis. Moreover, because TEs are present in virtually all species sequenced to date, their role in spatially restricted transcript expression is likely not exclusive to the species analyzed in this work.

Molecular basis of FAAH-OUT-associated human pain insensitivity

Chronic pain affects millions of people worldwide and new treatments are needed urgently. One way to identify novel analgesic strategies is to understand the biological dysfunctions that lead to human inherited pain insensitivity disorders. Here we report how the recently discovered brain and dorsal root ganglia-expressed FAAH-OUT long non-coding RNA (lncRNA) gene, which was found from studying a pain-insensitive patient with reduced anxiety and fast wound healing, regulates the adjacent key endocannabinoid system gene FAAH, which encodes the anandamide-degrading fatty acid amide hydrolase enzyme. We demonstrate that the disruption in FAAH-OUT lncRNA transcription leads to DNMT1-dependent DNA methylation within the FAAH promoter. In addition, FAAH-OUT contains a conserved regulatory element, FAAH-AMP, that acts as an enhancer for FAAH expression. Furthermore, using transcriptomic analyses in patient-derived cells we have uncovered a network of genes that are dysregulated from disruption of the FAAH-FAAH-OUT axis, thus providing a coherent mechanistic basis to understand the human phenotype observed. Given that FAAH is a potential target for the treatment of pain, anxiety, depression and other neurological disorders, this new understanding of the regulatory role of the FAAH-OUT gene provides a platform for the development of future gene and small molecule therapies.

Characterization of an eye field-like state during optic vesicle organoid development

Specification of the eye field (EF) within the neural plate marks the earliest detectable stage of eye development. Experimental evidence, primarily from non-mammalian model systems, indicates that the stable formation of this group of cells requires the activation of a set of key transcription factors. This crucial event is challenging to probe in mammals and, quantitatively, little is known regarding the regulation of the transition of cells to this ocular fate. Using optic vesicle organoids to model the onset of the EF, we generate time-course transcriptomic data allowing us to identify dynamic gene expression programmes that characterize this cellular-state transition. Integrating this with chromatin accessibility data suggests a direct role of canonical EF transcription factors in regulating these gene expression changes, and highlights candidate cis-regulatory elements through which these transcription factors act. Finally, we begin to test a subset of these candidate enhancer elements, within the organoid system, by perturbing the underlying DNA sequence and measuring transcriptomic changes during EF activation.

Histone exchange sensors reveal variant specific dynamics in mouse embryonic stem cells

Eviction of histones from nucleosomes and their exchange with newly synthesized or alternative variants is a central epigenetic determinant. Here, we define the genome-wide occupancy and exchange pattern of canonical and non-canonical histone variants in mouse embryonic stem cells by genetically encoded exchange sensors. While exchange of all measured variants scales with transcription, we describe variant-specific associations with transcription elongation and Polycomb binding. We found considerable exchange of H3.1 and H2B variants in heterochromatin and repeat elements, contrasting the occupancy and little exchange of H3.3 in these regions. This unexpected association between H3.3 occupancy and exchange of canonical variants is also evident in active promoters and enhancers, and further validated by reduced H3.1 dynamics following depletion of H3.3-specific chaperone, HIRA. Finally, analyzing transgenic mice harboring H3.1 or H3.3 sensors demonstrates the vast potential of this system for studying histone exchange and its impact on gene expression regulation in vivo.

Regulatory architecture of housekeeping genes is driven by promoter assemblies

Genes that are key to cell identity are generally regulated by cell-type-specific enhancer elements bound by transcription factors, some of which facilitate looping to distant gene promoters. In contrast, genes that encode housekeeping functions, whose regulation is essential for normal cell metabolism and growth, generally lack interactions with distal enhancers. We find that Ronin (Thap11) assembles multiple promoters of housekeeping and metabolic genes to regulate gene expression. This behavior is analogous to how enhancers are brought together with promoters to regulate cell identity genes. Thus, Ronin-dependent promoter assemblies provide a mechanism to explain why housekeeping genes can forgo distal enhancer elements and why Ronin is important for cellular metabolism and growth control. We propose that clustering of regulatory elements is a mechanism common to cell identity and housekeeping genes but is accomplished by different factors binding distinct control elements to establish enhancer-promoter or promoter-promoter interactions, respectively.

Enhancer activation via TCP and HD-ZIP and repression by Dof transcription factors mediate giant cell-specific expression

Proper cell-type identity relies on highly coordinated regulation of gene expression. Regulatory elements such as enhancers can produce cell type-specific expression patterns, but the mechanisms underlying specificity are not well understood. We previously identified an enhancer region capable of driving specific expression in giant cells, which are large, highly endoreduplicated cells in the Arabidopsis thaliana sepal epidermis. In this study, we use the giant cell enhancer as a model to understand the regulatory logic that promotes cell type-specific expression. Our dissection of the enhancer revealed that giant cell specificity is mediated primarily through the combination of two activators and one repressor. HD-ZIP and TCP transcription factors are involved in the activation of expression throughout the epidermis. High expression of HD-ZIP transcription factor genes in giant cells promoted higher expression driven by the enhancer in giant cells. Dof transcription factors repressed the activity of the enhancer such that only giant cells maintained enhancer activity. Thus, our data are consistent with a conceptual model whereby cell type-specific expression emerges from the combined activities of three transcription factor families activating and repressing expression in epidermal cells.

Key transcription factors influence the epigenetic landscape to regulate retinal cell differentiation

How the diverse neural cell types emerge from multipotent neural progenitor cells during central nervous system development remains poorly understood. Recent scRNA-seq studies have delineated the developmental trajectories of individual neural cell types in many neural systems including the neural retina. Further understanding of the formation of neural cell diversity requires knowledge about how the epigenetic landscape shifts along individual cell lineages and how key transcription factors regulate these changes. In this study, we dissect the changes in the epigenetic landscape during early retinal cell differentiation by scATAC-seq and identify globally the enhancers, enriched motifs, and potential interacting transcription factors underlying the cell state/type specific gene expression in individual lineages. Using CUT&Tag, we further identify the enhancers bound directly by four key transcription factors, Otx2, Atoh7, Pou4f2 and Isl1, including those dependent on Atoh7, and uncover the sequential and combinatorial interactions of these factors with the epigenetic landscape to control gene expression along individual retinal cell lineages such as retinal ganglion cells (RGCs). Our results reveal a general paradigm in which transcription factors collaborate and compete to regulate the emergence of distinct retinal cell types such as RGCs from multipotent retinal progenitor cells (RPCs).

Expression of Human Endogenous Retrovirus Group K (HERV-K) HML-2 Correlates with Immune Activation of Macrophages and Type I Interferon Response

Human endogenous retroviruses (HERVs) comprise about 8.3% of the human genome and are capable of producing RNA molecules that can be sensed by pattern recognition receptors, leading to the activation of innate immune response pathways. The HERV-K (HML-2) subgroup is the youngest HERV clade with the highest degree of coding competence. Its expression is associated with inflammation-related diseases. However, the precise HML-2 loci, stimuli, and signaling pathways involved in these associations are not well understood or defined. To elucidate HML-2 expression on a locus-specific level, we used the retroelement sequencing tools TEcount and Telescope to analyze publicly available transcriptome sequencing (RNA-seq) and chromatin immunoprecipitation (ChIP) sequencing data sets of macrophages treated with a wide range of agonists. We found that macrophage polarization significantly correlates with modulation of the expression of specific HML-2 proviral loci. Further analysis demonstrated that the provirus HERV-K102, located in an intergenic region of locus 1q22, constituted the majority of the HML-2 derived transcripts following pro-inflammatory (M1) polarization and was upregulated explicitly in response to interferon gamma (IFN-γ) signaling. We found that signal transducer and activator of transcription 1 and interferon regulatory factor 1 interact with a solo long terminal repeat (LTR) located upstream of HERV-K102, termed LTR12F, following IFN-γ signaling. Using reporter constructs, we demonstrated that LTR12F is critical for HERV-K102 upregulation by IFN-γ. In THP1-derived macrophages, knockdown of HML-2 or knockout of MAVS, an adaptor of RNA-sensing pathways, significantly downregulated genes containing interferon-stimulated response elements (ISREs) in their promoters, suggesting an intermediate role of HERV-K102 in the switch from IFN-γ signaling to the activation of type I interferon expression and, therefore, in a positive feedback loop to enhance pro-inflammatory signaling. IMPORTANCE The human endogenous retrovirus group K subgroup, HML-2, is known to be elevated in a long list of inflammation-associated diseases. However, a clear mechanism for HML-2 upregulation in response to inflammation has not been defined. In this study, we identify a provirus of the HML-2 subgroup, HERV-K102, which is significantly upregulated and constitutes the majority of the HML-2 derived transcripts in response to pro-inflammatory activation of macrophages. Moreover, we identify the mechanism of HERV-K102 upregulation and demonstrate that HML-2 expression enhances interferon-stimulated response element activation. We also demonstrate that this provirus is elevated in vivo and correlates with interferon gamma signaling activity in cutaneous leishmaniasis patients. This study provides key insights into the HML-2 subgroup and suggests that it may participate in enhancing pro-inflammatory signaling in macrophages and probably other immune cells.

Deciphering the multi-scale, quantitative cis-regulatory code

Uncovering the cis-regulatory code that governs when and how much each gene is transcribed in a given genome and cellular state remains a central goal of biology. Here, we discuss major layers of regulation that influence how transcriptional outputs are encoded by DNA sequence and cellular context. We first discuss how transcription factors bind specific DNA sequences in a dosage-dependent and cooperative manner and then proceed to the cofactors that facilitate transcription factor function and mediate the activity of modular cis-regulatory elements such as enhancers, silencers, and promoters. We then consider the complex and poorly understood interplay of these diverse elements within regulatory landscapes and its relationships with chromatin states and nuclear organization. We propose that a mechanistically informed, quantitative model of transcriptional regulation that integrates these multiple regulatory layers will be the key to ultimately cracking the cis-regulatory code.

Twist1-YY1-p300 complex promotes the malignant progression of HCC through activation of miR-9 by forming phase-separated condensates at super-enhancers and relieved by metformin

Hepatocellular carcinoma (HCC) is one of the leading causes of death, which deserves further study to reveal the underlying molecular mechanisms. Studies have shown that miR-9 in associated with poor prognosis in HCC patients. However, the mechanisms of transcriptional activation regulation of miR-9 and its role in the malignant progression of HCC have been rarely investigated. Some transcriptional coactivators can form phase-separated condensates at super-enhancers that compartmentalize and concentrate the transcription apparatus to drive robust gene expression. Here, we demonstrate that Twist1 and YY1 could form a transcriptional complex with p300, creating local high-concentration phase-separated interaction hubs at the super-enhancers of miR-9 and activate its expression to promote the malignant progression of HCC by stimulating the migration and invasion of hepatocellular carcinoma cells. Twist1-YY1-p300 phase-separated condensates were disrupted by metformin (Met) and thus reduce miR-9 expression, thereby inhibiting the malignant progression of HCC. Our study demonstrates that the Twist1 transcriptional factor complex involved in the malignant progression of HCC can form phase separation condensates at super-enhancers of miR-9 to promote the expression of oncogenes in HCC cells. It provides a potential target for the therapy of HCC and offers insights into the mechanism of Met in HCC inhibition.

E2F1 combined with LINC01004 super-enhancer to promote hepatocellular carcinoma cell proliferation and metastasis

INTRODUCTION

Super-enhancer-associated lncRNAs play important roles in the occurrence and development of malignant tumors, including hepatocellular carcinoma (HCC).

OBJECTIVES

The current work aimed to identify and characterize super-enhancer-associated lncRNAs in the pathogenesis of HCC.

METHODS

H3K27ac ChIP-seq data from HepG2 cell line and two HCC tissues were used to identify super-enhancer-associated lncRNAs in HCC. JQ-1 treatment and CRISPR-dCas9 system were performed to confirm super-enhancer activity. Quantitative real-time PCR (qPCR), ChIP-qPCR, and dual-luciferase reporter system assay demonstrated the regulation of E2F1 on super-enhancer. Functional loss experiment was used to identify the function of LINC01004.

RESULTS

In this study, we identified and characterized LINC01004, a novel super-enhancer-associated lncRNA, as a crucial oncogene in HCC. LINC01004 was upregulated in liver cancer tissues and was associated with poor patient prognosis. Moreover, LINC01004 promoted cell proliferation and metastasis of HCC. The binding of E2F1 to the super-enhancer could promote the transcription of LINC01004, while the inhibition of super-enhancer activity decreased LINC01004 expression.

CONCLUSION

This finding might provide mechanistic insights into the molecular mechanisms underlying hepatocarcinogenesis and the biological function of super-enhancer. LINC01004 can serve as a potential therapeutic target for HCC patient.

Loss of the Ash2l subunit of histone H3K4 methyltransferase complexes reduces chromatin accessibility at promoters

Changes in gene expression programs are intimately linked to cell fate decisions. Post-translational modifications of core histones contribute to control gene expression. Methylation of lysine 4 of histone H3 (H3K4) correlates with active promoters and gene transcription. This modification is catalyzed by KMT2 methyltransferases, which require interaction with 4 core subunits, WDR5, RBBP5, ASH2L and DPY30, for catalytic activity. Ash2l is necessary for organismal development and for tissue homeostasis. In mouse embryo fibroblasts (MEFs), Ash2l loss results in gene repression, provoking a senescence phenotype. We now find that upon knockout of Ash2l both H3K4 mono- and tri-methylation (H3K4me1 and me3, respectively) were deregulated. In particular, loss of H3K4me3 at promoters correlated with gene repression, especially at CpG island promoters. Ash2l loss resulted in increased loading of histone H3 and reduced chromatin accessibility at promoters, accompanied by an increase of repressing and a decrease of activating histone marks. Moreover, we observed altered binding of CTCF upon Ash2l loss. Lost and gained binding was noticed at promoter-associated and intergenic sites, respectively. Thus, Ash2l loss and reduction of H3K4me3 correlate with altered chromatin accessibility and transcription factor binding. These findings contribute to a more detailed understanding of mechanistic consequences of H3K4me3 loss and associated repression of gene transcription and thus of the observed cellular consequences.

Systematic discovery and functional dissection of enhancers needed for cancer cell fitness and proliferation

A scarcity of functionally validated enhancers in the human genome presents a significant hurdle to understanding how these cis-regulatory elements contribute to human diseases. We carry out highly multiplexed CRISPR-based perturbation and sequencing to identify enhancers required for cell proliferation and fitness in 10 human cancer cell lines. Our results suggest that the cell fitness enhancers, unlike their target genes, display high cell-type specificity of chromatin features. They typically adopt a modular structure, comprised of activating elements enriched for motifs of oncogenic transcription factors, surrounded by repressive elements enriched for motifs recognized by transcription factors with tumor suppressor functions. We further identify cell fitness enhancers that are selectively accessible in clinical tumor samples, and the levels of chromatin accessibility are associated with patient survival. These results reveal functional enhancers across multiple cancer cell lines, characterize their context-dependent chromatin organization, and yield insights into altered transcription programs in cancer cells.

Super-enhancers in esophageal carcinoma: Transcriptional addictions and therapeutic strategies

The tumorigenesis of esophageal carcinoma arises from transcriptional dysregulation would become exceptionally dependent on specific regulators of gene expression, which could be preferentially attributed to the larger non-coding cis-regulatory elements, i.e. super-enhancers (SEs). SEs, large genomic regulatory entity in close genomic proximity, are underpinned by control cancer cell identity. As a consequence, the transcriptional addictions driven by SEs could offer an Achilles’ heel for molecular treatments on patients of esophageal carcinoma and other types of cancer as well. In this review, we summarize the recent findings about the oncogenic SEs upon which esophageal cancer cells depend, and discuss why SEs could be seen as the hallmark of cancer, how transcriptional dependencies driven by SEs, and what opportunities could be supplied based on this cancer-specific SEs.

JUN Regulation of Injury-Induced Enhancers in Schwann Cells

Schwann cells play a critical role after peripheral nerve injury by clearing myelin debris, forming axon-guiding bands of Büngner, and remyelinating regenerating axons. Schwann cells undergo epigenomic remodeling to differentiate into a repair state that expresses unique genes, some of which are not expressed at other stages of Schwann cell development. We previously identified a set of enhancers that are activated in Schwann cells after nerve injury, and we determined whether these enhancers are preprogrammed into the Schwann cell epigenome as poised enhancers before injury. Poised enhancers share many attributes of active enhancers, such as open chromatin, but are marked by repressive histone H3 lysine 27 (H3K27) trimethylation rather than H3K27 acetylation. We find that most injury-induced enhancers are not marked as poised enhancers before injury indicating that injury-induced enhancers are not preprogrammed in the Schwann cell epigenome. Injury-induced enhancers are enriched with AP-1 binding motifs, and the c-JUN subunit of AP-1 had been shown to be critical to drive the transcriptional response of Schwann cells after injury. Using in vivo chromatin immunoprecipitation sequencing analysis in rat, we find that c-JUN binds to a subset of injury-induced enhancers. To test the role of specific injury-induced enhancers, we focused on c-JUN-binding enhancers upstream of the Sonic hedgehog (Shh) gene, which is only upregulated in repair Schwann cells compared with other stages of Schwann cell development. We used targeted deletions in male/female mice to show that the enhancers are required for robust induction of the Shh gene after injury.SIGNIFICANCE STATEMENT The proregenerative actions of Schwann cells after nerve injury depends on profound reprogramming of the epigenome. The repair state is directed by injury-induced transcription factors, like JUN, which is uniquely required after nerve injury. In this study, we test whether the injury program is preprogrammed into the epigenome as poised enhancers and define which enhancers bind JUN. Finally, we test the roles of these enhancers by performing clustered regularly interspaced short palindromic repeat (CRISPR)-mediated deletion of JUN-bound injury enhancers in the Sonic hedgehog gene. Although many long-range enhancers drive expression of Sonic hedgehog at different developmental stages of specific tissues, these studies identify an entirely new set of enhancers that are required for Sonic hedgehog induction in Schwann cells after injury.

The tissue-specificity associated region and motif of an emx2 downstream enhancer CNE2.04 in zebrafish

BACKGROUND

Expression level of EMX2 plays an important role in the development of nervous system and cancers. CNE2.04, a conserved enhancer downstream of emx2, drives fluorescent protein expression in the similar pattern of emx2.

METHODS

CNE2.04 truncated or motif-mutated transgenic reporter plasmids were constructed and injected into the zebrafish fertilized egg with Tol2 mRNA at the unicellular stage of zebrafish eggs. The green fluorescence expression patterns were observed at 24, 48, and 72 hpf, and the fluorescence rates of different tissues were counted at 48 hpf.

RESULTS

Compared to CNE2.04, CNE2.04-R400 had comparable enhancer activity, while the tissue specificity of CNE2.04-L400 was obviously changed. Motif CCCCTC mutation obviously changed the enhancer activity, while motif CCGCTC mutations also changed it.

CONCLUSION

Due to their correlation with tissue specificity, CNE2.04-R400 is associated with the tissue-specificity of CNE2.04, and motif CCCCTC plays an important role in the enhancer activity of CNE2.04.

Functional analysis of the Vsx2 super-enhancer uncovers distinct cis-regulatory circuits controlling Vsx2 expression during retinogenesis

Vsx2 is a transcription factor essential for retinal proliferation and bipolar cell differentiation, but the molecular mechanisms underlying its developmental roles are unclear. Here, we have profiled VSX2 genomic occupancy during mouse retinogenesis, revealing extensive retinal genetic programs associated with VSX2 during development. VSX2 binds and transactivates its enhancer in association with the transcription factor PAX6. Mice harboring deletions in the Vsx2 regulatory landscape exhibit specific abnormalities in retinal proliferation and in bipolar cell differentiation. In one of those deletions, a complete loss of bipolar cells is associated with a bias towards photoreceptor production. VSX2 occupies cis-regulatory elements nearby genes associated with photoreceptor differentiation and homeostasis in the adult mouse and human retina, including a conserved region nearby Prdm1, a factor implicated in the specification of rod photoreceptors and suppression of bipolar cell fate. VSX2 interacts with the transcription factor OTX2 and can act to suppress OTX2-dependent enhancer transactivation of the Prdm1 enhancer. Taken together, our analyses indicate that Vsx2 expression can be temporally and spatially uncoupled at the enhancer level, and they illuminate important mechanistic insights into how VSX2 is engaged with gene regulatory networks that are essential for retinal proliferation and cell fate acquisition.

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene's promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle

Background

Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined.

Results

In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers.

Conclusions

This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08737-9.

TET1 dioxygenase is required for FOXA2-associated chromatin remodeling in pancreatic beta-cell differentiation

Existing knowledge of the role of epigenetic modifiers in pancreas development has exponentially increased. However, the function of TET dioxygenases in pancreatic endocrine specification remains obscure. We set out to tackle this issue using a human embryonic stem cell (hESC) differentiation system, in which TET1/TET2/TET3 triple knockout cells display severe defects in pancreatic β-cell specification. The integrative whole-genome analysis identifies unique cell-type-specific hypermethylated regions (hyper-DMRs) displaying reduced chromatin activity and remarkable enrichment of FOXA2, a pioneer transcription factor essential for pancreatic endoderm specification. Intriguingly, TET depletion leads to significant changes in FOXA2 binding at the pancreatic progenitor stage, in which gene loci with decreased FOXA2 binding feature low levels of active chromatin modifications and enriches for bHLH motifs. Transduction of full-length TET1 but not the TET1-catalytic-domain in TET-deficient cells effectively rescues β-cell differentiation accompanied by restoring PAX4 hypomethylation. Taking these findings together with the defective generation of functional β-cells upon TET1-inactivation, our study unveils an essential role of TET1-dependent demethylation in establishing β-cell identity. Moreover, we discover a physical interaction between TET1 and FOXA2 in endodermal lineage intermediates, which provides a mechanistic clue regarding the complex crosstalk between TET dioxygenases and pioneer transcription factors in epigenetic regulation during pancreas specification.

Here the authors show that TET1 is required for the generation of functional insulin-producing cells, FOXA2 physically interacts with TET1 and contributes to specific recruitment of TET1 to mediate chromatin opening at the regulatory elements of pancreatic lineage determinants.

JMJD3 intrinsically disordered region links the 3D-genome structure to TGFβ-dependent transcription activation

Enhancers are key regulatory elements that govern gene expression programs in response to developmental signals. However, how multiple enhancers arrange in the 3D-space to control the activation of a specific promoter remains unclear. To address this question, we exploited our previously characterized TGFβ-response model, the neural stem cells, focusing on a ~374 kb locus where enhancers abound. Our 4C-seq experiments reveal that the TGFβ pathway drives the assembly of an enhancer-cluster and precise gene activation. We discover that the TGFβ pathway coactivator JMJD3 is essential to maintain these structures. Using live-cell imaging techniques, we demonstrate that an intrinsically disordered region contained in JMJD3 is involved in the formation of phase-separated biomolecular condensates, which are found in the enhancer-cluster. Overall, in this work we uncover novel functions for the coactivator JMJD3, and we shed light on the relationships between the 3D-conformation of the chromatin and the TGFβ-driven response during mammalian neurogenesis.

Here the authors demonstrate that TGFβ drives multi-enhancer contacts and ultimately gene activation during neuronal commitment, and that this requires the intrinsically disordered region (IDR) of the histone demethylase JMJD3 likely through its role in promoting phase-separated biomolecular condensates.

Leveraging transcriptome and epigenome landscapes to infer regulatory networks during the onset of sexual maturation

Background

Despite sexual development being ubiquitous to vertebrates, the molecular mechanisms underpinning this fundamental transition remain largely undocumented in many organisms. We designed a time course experiment that successfully sampled the period when Atlantic salmon commence their trajectory towards sexual maturation.

Results

Through deep RNA sequencing, we discovered key genes and pathways associated with maturation in the pituitary-ovarian axis. Analyzing DNA methylomes revealed a bias towards hypermethylation in ovary that implicated maturation-related genes. Co-analysis of DNA methylome and gene expression changes revealed chromatin remodeling genes and key transcription factors were both significantly hypermethylated and upregulated in the ovary during the onset of maturation. We also observed changes in chromatin state landscapes that were strongly correlated with fundamental remodeling of gene expression in liver. Finally, a multiomic integrated analysis revealed regulatory networks and identified hub genes including TRIM25 gene (encoding the estrogen-responsive finger protein) as a putative key regulator in the pituitary that underwent a 60-fold change in connectivity during the transition to maturation.

Conclusion

The study successfully documented transcriptome and epigenome changes that involved key genes and pathways acting in the pituitary – ovarian axis. Using a Systems Biology approach, we identified hub genes and their associated networks deemed crucial for onset of maturation. The results provide a comprehensive view of the spatiotemporal changes involved in a complex trait and opens the door to future efforts aiming to manipulate puberty in an economically important aquaculture species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08514-8.

Genome-wide identification of associations between enhancer and alternative splicing in human and mouse

Background

Alternative splicing (AS) increases the diversity of transcriptome and could fine-tune the function of genes, so that understanding the regulation of AS is vital. AS could be regulated by many different cis-regulatory elements, such as enhancer. Enhancer has been experimentally proved to regulate AS in some genes. However, there is a lack of genome-wide studies on the association between enhancer and AS (enhancer-AS association). To bridge the gap, here we developed an integrative analysis on a genome-wide scale to identify enhancer-AS associations in human and mouse.

Result

We collected enhancer datasets which include 28 human and 24 mouse tissues and cell lines, and RNA-seq datasets which are paired with the selected tissues. Combining with data integration and statistical analysis, we identified 3,242 human and 7,716 mouse genes which have significant enhancer-AS associations in at least one tissue. On average, for each gene, about 6% of enhancers in human (5% in mouse) are associated to AS change and for each enhancer, approximately one gene is identified to have enhancer-AS association in both human and mouse. We found that 52% of the human significant (34% in mouse) enhancer-AS associations are the co-existence of homologous genes and homologous enhancers. We further constructed a user-friendly platform, named Visualization of Enhancer-associated Alternative Splicing (VEnAS, http://venas.iis.sinica.edu.tw/), to provide genomic architecture, intuitive association plot, and contingency table of the significant enhancer-AS associations.

Conclusion

This study provides the first genome-wide identification of enhancer-AS associations in human and mouse. The results suggest that a notable portion of enhancers are playing roles in AS regulations. The analyzed results and the proposed platform VEnAS would provide a further understanding of enhancers on regulating alternative splicing.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08537-1.

Long-Distance Repression by Human Silencers: Chromatin Interactions and Phase Separation in Silencers

Three-dimensional genome organization represents an additional layer in the epigenetic regulation of gene expression. Active transcription controlled by enhancers or super-enhancers has been extensively studied. Enhancers or super-enhancers can recruit activators or co-activators to activate target gene expression through long-range chromatin interactions. Chromatin interactions and phase separation play important roles in terms of enhancer or super-enhancer functioning. Silencers are another major type of cis-regulatory element that can mediate gene regulation by turning off or reducing gene expression. However, compared to active transcription, silencer studies are still in their infancy. This review covers the current knowledge of human silencers, especially the roles of chromatin interactions and phase separation in silencers. This review also proposes future directions for human silencer studies.

MAPK-mediated transcription factor GATAd contributes to Cry1Ac resistance in diamondback moth by reducing PxmALP expression

The benefits of biopesticides and transgenic crops based on the insecticidal Cry-toxins from Bacillus thuringiensis (Bt) are considerably threatened by insect resistance evolution, thus, deciphering the molecular mechanisms underlying insect resistance to Bt products is of great significance to their sustainable utilization. Previously, we have demonstrated that the down-regulation of PxmALP in a strain of Plutella xylostella (L.) highly resistant to the Bt Cry1Ac toxin was due to a hormone-activated MAPK signaling pathway and contributed to the resistance phenotype. However, the underlying transcriptional regulatory mechanism remains enigmatic. Here, we report that the PxGATAd transcription factor (TF) is responsible for the differential expression of PxmALP observed between the Cry1Ac susceptible and resistant strains. We identified that PxGATAd directly activates PxmALP expression via interacting with a non-canonical but specific GATA-like cis-response element (CRE) located in the PxmALP promoter region. A six-nucleotide insertion mutation in this cis-acting element of the PxmALP promoter from the resistant strain resulted in repression of transcriptional activity, affecting the regulatory performance of PxGATAd. Furthermore, silencing of PxGATAd in susceptible larvae reduced the expression of PxmALP and susceptibility to Cry1Ac toxin. Suppressing PxMAP4K4 expression in the resistant larvae transiently recovered both the expression of PxGATAd and PxmALP, indicating that the PxGATAd is a positive responsive factor involved in the activation of PxmALP promoter and negatively regulated by the MAPK signaling pathway. Overall, this study deciphers an intricate regulatory mechanism of PxmALP gene expression and highlights the concurrent involvement of both trans-regulatory factors and cis-acting elements in Cry1Ac resistance development in lepidopteran insects.

Gene expression and regulation are associated with adaptive evolution in living organisms. The rapid evolution of insect resistance to Bt insecticidal Cry toxins is frequently associated with reduced expression of diverse midgut genes that code for Cry-toxin receptors. Nonetheless, our current knowledge about the regulation of gene expression of these pivotal receptor genes in insects is limited. Membrane-bound alkaline phosphatase (mALP) is a known receptor for Cry1Ac toxin in diverse insects and here, we report the transcriptional regulatory mechanism of the PxmALP gene related to Cry1Ac resistance in P. xylostella. We identified a MAPK signaling pathway that negatively regulates the PxGATAd transcriptional factor which is involved in the differential expression of PxmALP via interacting with the PxmALP promoter. Furthermore, a cis-acting element mutation repressing the regulatory activity of PxGATAd for PxmALP expression in the Cry1Ac resistant strain was identified. Our study provides an insight into the precise transcriptional regulatory mechanism that regulates PxmALP expression and is involved in the evolution of Bt Cry1Ac resistance in P. xylostella, which provides a paradigm for decoding the regulation landscape of midgut Cry-toxin receptor genes in insects.

CDK7-dependent transcriptional addiction in bone and soft tissue sarcomas: Present and Future

Cancer arises from genetic alterations that invariably contribute to dysregulated transcriptional programs. These dysregulated programs establish and maintain specific cancer cell states, leading to an intensive dependence on a set of certain regulators of gene expression. The CDK7 functions as the core of transcription, and governs RNA polymerase II and the downstream oncogenes expression in cancers. CDK7 inhibition leads to reduced recruitment of super-enhancers-driven oncogenic transcription factors, and the depression of these associated oncogenes expression, which indicates the dependence of transcriptional addiction of cancers on CDK7. Given that specified oncoproteins of sarcomas commonly function at oncogenic transcription, targeting CDK7-denpendent transcriptional addiction may be of guiding significance for the treatment of sarcomas. In this review, we summarize the advances in mechanism of targeted CDK7-dependent transcriptional addiction and discuss the path ahead to potential application discovery in bone and soft tissue sarcomas, providing theoretical considerations for bio-orthogonal therapeutic strategies.

FoxP3 associates with enhancer-promoter loops to regulate Treg-specific gene expression

Gene expression programs are specified by higher-order chromatin structure and enhancer-promoter loops (EPLs). T regulatory cell (Treg) identity is dominantly specified by the transcription factor (TF) FoxP3, whose mechanism of action is unclear. We applied chromatin conformation capture with immunoprecipitation (HiChIP) in Treg and closely related conventional CD4+ T cells (Tconv). EPLs identified by H3K27Ac HiChIP showed a range of connection intensity, with some superconnected genes. TF-specific HiChIP showed that FoxP3 interacts with EPLs at a large number of genes, including some not differentially expressed in Treg versus Tconv, but enriched at the core Treg signature loci that it up-regulates. FoxP3 association correlated with heightened H3K27Ac looping, as ascertained by analysis of FoxP3-deficient Treg-like cells. There was marked asymmetry in the loci where FoxP3 associated at the enhancer- or the promoter-side of EPLs, with enrichment for different transcriptional cofactors. FoxP3 EPL intensity distinguished gene clusters identified by single-cell ATAC-seq as covarying between individual Tregs, supporting a direct transactivation model for FoxP3 in determining Treg identity.

Exploring the role of eRNA in regulating gene expression

eRNAs as the products of enhancers can regulate gene expression via various possible ways, but which regulation way is more reasonable is debatable in biology, and in particular, how eRNAs impact gene expression remains unclear. Here we introduce a mechanistic model of gene expression to address these issues. This model considers three possible regulation ways of eRNA: Type-I by which eRNA regulates transcriptional activity by facilitating the formation of enhancer-promoter (E-P) loop, Type-II by which eRNA directly promotes the mRNA production rate, and mixed regulation (i.e., the combination of Type-I and Type-II). We show that with the increase of the E-P loop length, mRNA distribution can transition from unimodality to bimodality or vice versa in all the three regulation cases. However, in contrast to the other two regulations, Type-II regulation can lead to the highest mean mRNA level and the lowest mRNA noise, independent of the E-P loop length. These results would not only reveal the essential mechanism of how eRNA regulates gene expression, but also imply a new mechanism for phenotypic switching, namely the E-P loop can induce phenotypic switching.

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.

Olfactory receptor coding sequences cause silencing of episomal constructs in multiple cell lines

The mammalian olfactory system consists of sensory neurons with specialized odorant-binding capability accomplished by mutually exclusive odorant receptor (OR) expression. Mutually exclusive OR expression is a complex multi-step process regulated by a number of cis and trans factors, including pan-silencing of all OR genes preceding the robust and stable expression of the one OR selected in each sensory neuron. We transfected two olfactory-placode-derived cell lines modeling immature odorant sensory neurons, as well as the GD25 fibroblast cell line, with episomes containing CMV-driven GFP and TK-driven hygromycin reporter genes. We inserted various coding sequences, along with an IRES, immediately upstream of the GFP gene to produce bicistronic mRNAs driven from the local CMV promoter. We found that the presence of several OR coding sequences resulted in significantly diminished episomal expression of GFP in all three cell lines. These findings suggest that OR coding sequences have intrinsic self-silencing capability that might facilitate mutually exclusive OR expression in olfactory sensory neurons by making it less likely that multiple ORs acquire an above-threshold level of expression at once.

Identifying oligodendrocyte enhancers governing Plp1 expression

Oligodendrocytes (OLs) produce myelin in the central nervous system (CNS), which accelerates the propagation of action potentials and supports axonal integrity. As a major component of CNS myelin, proteolipid protein 1 (Plp1) is indispensable for the axon-supportive function of myelin. Notably, this function requires the continuous high-level expression of Plp1 in OLs. Equally important is the controlled expression of Plp1, as illustrated by Pelizaeus-Merzbacher disease for which the most common cause is PLP1 overexpression. Despite a decade-long search, promoter-distal OL enhancers that govern Plp1 remain elusive. We have recently developed an innovative method that maps promoter-distal enhancers to genes in a principled manner. Here, we applied it to Plp1, uncovering two OL enhancers for it (termed Plp1-E1 and Plp1-E2). Remarkably, clustered regularly interspaced short palindromic repeats (CRISPR) interference epigenome editing showed that Plp1-E1 and Plp1-E2 do not regulate two genes in their vicinity, highlighting their exquisite specificity to Plp1. Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq) data show that Plp1-E1 and Plp1-E2 are OL-specific enhancers that are conserved among human, mouse and rat. Hi-C data reveal that the physical interactions between Plp1-E1/2 and PLP1 are among the strongest in OLs and specific to OLs. We also show that Myrf, a master regulator of OL development, acts on Plp1-E1 and Plp1-E2 to promote Plp1 expression.

Building a Mammalian Retina: An Eye on Chromatin Structure

Regulation of gene expression by chromatin structure has been under intensive investigation, establishing nuclear organization and genome architecture as a potent and effective means of regulating developmental processes. The substantial growth in our knowledge of the molecular mechanisms underlying retinogenesis has been powered by several genome-wide based tools that mapped chromatin organization at multiple cellular and biochemical levels. Studies profiling the retinal epigenome and transcriptome have allowed the systematic annotation of putative cis-regulatory elements associated with transcriptional programs that drive retinal neural differentiation, laying the groundwork to understand spatiotemporal retinal gene regulation at a mechanistic level. In this review, we outline recent advances in our understanding of the chromatin architecture in the mammalian retina during development and disease. We focus on the emerging roles of non-coding regulatory elements in controlling retinal cell-type specific transcriptional programs, and discuss potential implications in untangling the etiology of eye-related disorders.

Predict long-range enhancer regulation based on protein-protein interactions between transcription factors

Long-range regulation by distal enhancers plays critical roles in cell-type specific transcriptional programs. Computational predictions of genome-wide enhancer-promoter interactions are still challenging due to limited accuracy and the lack of knowledge on the molecular mechanisms. Based on recent biological investigations, the protein-protein interactions (PPIs) between transcription factors (TFs) have been found to participate in the regulation of chromatin loops. Therefore, we developed a novel predictive model for cell-type specific enhancer-promoter interactions by leveraging the information of TF PPI signatures. Evaluated by a series of rigorous performance comparisons, the new model achieves superior performance over other methods. The model also identifies specific TF PPIs that may mediate long-range regulatory interactions, revealing new mechanistic understandings of enhancer regulation. The prioritized TF PPIs are associated with genes in distinct biological pathways, and the predicted enhancer-promoter interactions are strongly enriched with cis-eQTLs. Most interestingly, the model discovers enhancer-mediated trans-regulatory links between TFs and genes, which are significantly enriched with trans-eQTLs. The new predictive model, along with the genome-wide analyses, provides a platform to systematically delineate the complex interplay among TFs, enhancers and genes in long-range regulation. The novel predictions also lead to mechanistic interpretations of eQTLs to decode the genetic associations with gene expression.

Role of Enhancers in Development and Diseases

Enhancers are cis-regulatory elements containing short DNA sequences that serve as binding sites for pioneer/regulatory transcription factors, thus orchestrating the regulation of genes critical for lineage determination. The activity of enhancer elements is believed to be determined by transcription factor binding, thus determining the cell state identity during development. Precise spatio-temporal control of the transcriptome during lineage specification requires the coordinated binding of lineage-specific transcription factors to enhancers. Thus, enhancers are the primary determinants of cell identity. Numerous studies have explored the role and mechanism of enhancers during development and disease, and various basic questions related to the functions and mechanisms of enhancers have not yet been fully answered. In this review, we discuss the recently published literature regarding the roles of enhancers, which are critical for various biological processes governing development. Furthermore, we also highlight that altered enhancer landscapes provide an essential context to understand the etiologies and mechanisms behind numerous complex human diseases, providing new avenues for effective enhancer-based therapeutic interventions.

The role of epigenetic modifications, long-range contacts, enhancers and topologically associating domains in the regulation of glioma grade-specific genes

Genome-wide studies have uncovered specific genetic alterations, transcriptomic patterns and epigenetic profiles associated with different glioma types. We have recently created a unique atlas encompassing genome-wide profiles of open chromatin, histone H3K27ac and H3Kme3 modifications, DNA methylation and transcriptomes of 33 glioma samples of different grades. Here, we intersected genome-wide atlas data with topologically associating domains (TADs) and demonstrated that the chromatin organization and epigenetic landscape of enhancers have a strong impact on genes differentially expressed in WHO low grade versus high grade gliomas. We identified TADs enriched in glioma grade-specific genes and/or epigenetic marks. We found the set of transcription factors, including REST, E2F1 and NFKB1, that are most likely to regulate gene expression in multiple TADs, containing specific glioma-related genes. Moreover, many genes associated with the cell-matrix adhesion Gene Ontology group, in particular 14 PROTOCADHERINs, were found to be regulated by long-range contacts with enhancers. Presented results demonstrate the existence of epigenetic differences associated with chromatin organization driving differential gene expression in gliomas of different malignancy.

Nuclear condensates of p300 formed though the structured catalytic core can act as a storage pool of p300 with reduced HAT activity

The transcriptional co-activator and acetyltransferase p300 is required for fundamental cellular processes, including differentiation and growth. Here, we report that p300 forms phase separated condensates in the cell nucleus. The phase separation ability of p300 is regulated by autoacetylation and relies on its catalytic core components, including the histone acetyltransferase (HAT) domain, the autoinhibition loop, and bromodomain. p300 condensates sequester chromatin components, such as histone H3 tail and DNA, and are amplified through binding of p300 to the nucleosome. The catalytic HAT activity of p300 is decreased due to occlusion of the active site in the phase separated droplets, a large portion of which co-localizes with chromatin regions enriched in H3K27me3. Our findings suggest a model in which p300 condensates can act as a storage pool of the protein with reduced HAT activity, allowing p300 to be compartmentalized and concentrated at poised or repressed chromatin regions.

The histone acetyltransferase p300 mostly localizes to active chromatin; however, some repressed genes marked with H3K27me3 are also bound by p300. Here the authors show p300 is capable of phase separation, which relies on its catalytic core, and that p300 catalytic activity is decreased in phase-separated droplets that co-localize with H3K27me3-marked chromatin.

Interaction of the NRF2 and p63 transcription factors promotes keratinocyte proliferation in the epidermis

Epigenetic regulation of cell and tissue function requires the coordinated action of transcription factors. However, their combinatorial activities during regeneration remain largely unexplored. Here, we discover an unexpected interaction between the cytoprotective transcription factor NRF2 and p63- a key player in epithelial morphogenesis. Chromatin immunoprecipitation combined with sequencing and reporter assays identifies enhancers and promoters that are simultaneously activated by NRF2 and p63 in human keratinocytes. Modeling of p63 and NRF2 binding to nucleosomal DNA suggests their chromatin-assisted interaction. Pharmacological and genetic activation of NRF2 increases NRF2–p63 binding to enhancers and promotes keratinocyte proliferation, which involves the common NRF2–p63 target cyclin-dependent kinase 12. These results unravel a collaborative function of NRF2 and p63 in the control of epidermal renewal and suggest their combined activation as a strategy to promote repair of human skin and other stratified epithelia.