Orphan CpG islands amplify poised enhancer regulatory activity and determine target gene responsiveness

CpG island turnover events predict evolutionary changes in enhancer activity

BACKGROUND

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution.

RESULTS

We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon.

CONCLUSIONS

Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

KSHV genome harbors both constitutive and lytically induced enhancers

Kaposi's sarcoma-associated herpesvirus (KSHV) belongs to the gamma-herpesvirus family and is a well-known human oncogenic virus. In infected cells, the viral genome of 165 kbp is circular DNA wrapped in chromatin. The tight control of gene expression is critical for latency, the transition into the lytic phase, and the development of viral-associated malignancies. Distal cis-regulatory elements, such as enhancers and silencers, can regulate gene expression in a position- and orientation-independent manner. Open chromatin is another characteristic feature of enhancers. To systematically search for enhancers, we cloned all the open chromatin regions in the KSHV genome downstream of the luciferase gene and tested their enhancer activity in infected and uninfected cells. A silencer was detected upstream of the latency-associated nuclear antigen promoter. Two constitutive enhancers were identified in the K12p-OriLyt-R and ORF29 Intron regions, where ORF29 Intron is a tissue-specific enhancer. The following promoters: OriLyt-L, PANp, ALTp, and the terminal repeats (TRs) acted as lytically induced enhancers. The expression of the replication and transcription activator (RTA), the master regulator of the lytic cycle, was sufficient to induce the activity of lytic enhancers in uninfected cells. We propose that the TRs that span about 24 kbp region serve as a "viral super-enhancer" that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The silencer and enhancers described here provide an additional layer to the complex gene regulation of herpesviruses.IMPORTANCEIn this study, we performed a systematic functional assay to identify cis-regulatory elements within the genome of the oncogenic herpesvirus, Kaposi's sarcoma-associated herpesvirus (KSHV). Similar to other herpesviruses, KSHV presents both latent and lytic phases. Therefore, our assays were performed in uninfected cells, during latent infection, and under lytic conditions. We identified two constitutive enhancers, one of which seems to be a tissue-specific enhancer. In addition, four lytically induced enhancers, which are all responsive to the replication and transcription activator (RTA), were identified. Furthermore, a silencer was identified between the major latency promoter and the lytic gene locus. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The terminal repeats, spanning a region of about 24 kbp, seem like a "viral super-enhancer" that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA to regulate latency to lytic transition.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

Crosstalk within and beyond the Polycomb repressive system

The development of multicellular organisms depends on spatiotemporally controlled differentiation of numerous cell types and their maintenance. To generate such diversity based on the invariant genetic information stored in DNA, epigenetic mechanisms, which are heritable changes in gene function that do not involve alterations to the underlying DNA sequence, are required to establish and maintain unique gene expression programs. Polycomb repressive complexes represent a paradigm of epigenetic regulation of developmentally regulated genes, and the roles of these complexes as well as the epigenetic marks they deposit, namely H3K27me3 and H2AK119ub, have been extensively studied. However, an emerging theme from recent studies is that not only the autonomous functions of the Polycomb repressive system, but also crosstalks of Polycomb with other epigenetic modifications, are important for gene regulation. In this review, we summarize how these crosstalk mechanisms have improved our understanding of Polycomb biology and how such knowledge could help with the design of cancer treatments that target the dysregulated epigenome.

Polycomb protein binding and looping in the ON transcriptional state

Polycomb group (PcG) proteins mediate epigenetic silencing of important developmental genes by modifying histones and compacting chromatin through two major protein complexes, PRC1 and PRC2. These complexes are recruited to DNA by CpG islands (CGIs) in mammals and Polycomb response elements (PREs) in Drosophila. When PcG target genes are turned OFF, PcG proteins bind to PREs or CGIs, and PREs serve as anchors that loop together and stabilize gene silencing. Here, we address which PcG proteins bind to PREs and whether PREs mediate looping when their targets are in the ON transcriptional state. While the binding of most PcG proteins decreases at PREs in the ON state, one PRC1 component, Ph, remains bound. Further, PREs can loop to each other and with presumptive enhancers in the ON state and, like CGIs, may act as tethering elements between promoters and enhancers. Overall, our data suggest that PREs are important looping elements for developmental loci in both the ON and OFF states.

Increased enhancer-promoter interactions during developmental enhancer activation in mammals

Remote enhancers are thought to interact with their target promoters via physical proximity, yet the importance of this proximity for enhancer function remains unclear. Here we investigate the three-dimensional (3D) conformation of enhancers during mammalian development by generating high-resolution tissue-resolved contact maps for nearly a thousand enhancers with characterized in vivo activities in ten murine embryonic tissues. Sixty-one percent of developmental enhancers bypass their neighboring genes, which are often marked by promoter CpG methylation. The majority of enhancers display tissue-specific 3D conformations, and both enhancer-promoter and enhancer-enhancer interactions are moderately but consistently increased upon enhancer activation in vivo. Less than 14% of enhancer-promoter interactions form stably across tissues; however, these invariant interactions form in the absence of the enhancer and are likely mediated by adjacent CTCF binding. Our results highlight the general importance of enhancer-promoter physical proximity for developmental gene activation in mammals.

Structural variants in the Epb41l4a locus: TAD disruption and Nrep gene misregulation as hypothetical drivers of neurodevelopmental outcomes

Structural variations are a pervasive feature of human genomes, and there is growing recognition of their role in disease development through their impact on spatial chromatin architecture. This understanding has led us to investigate the clinical significance of CNVs in noncoding regions that influence TAD structures. In this study, we focused on the Epb41l4a locus, which contains a highly conserved TAD boundary present in both human chromosome 5 and mouse chromosome 18, and its association with neurodevelopmental phenotypes. Analysis of human data from the DECIPHER database indicates that CNVs within this locus, including both deletions and duplications, are often observed alongside neurological abnormalities, such as dyslexia and intellectual disability, although there is not enough evidence of a direct correlation or causative relationship. To investigate these possible associations, we generated mouse models with deletion and inversion mutations at this locus and carried out RNA-seq analysis to elucidate gene expression changes. We found that modifications in the Epb41l4a TAD boundary led to dysregulation of the Nrep gene, which plays a crucial role in nervous system development. These findings underscore the potential pathogenicity of these CNVs and highlight the crucial role of spatial genome architecture in gene expression regulation.

Enhancer-promoter interactions can form independently of genomic distance and be functional across TAD boundaries

Topologically Associating Domains (TADs) have been suggested to facilitate and constrain enhancer-promoter interactions. However, the role of TAD boundaries in effectively restricting these interactions remains unclear. Here, we show that a significant proportion of enhancer-promoter interactions are established across TAD boundaries in Drosophila embryos, but that developmental genes are strikingly enriched in intra- but not inter-TAD interactions. We pursued this observation using the twist locus, a master regulator of mesoderm development, and systematically relocated one of its enhancers to various genomic locations. While this developmental gene can establish inter-TAD interactions with its enhancer, the functionality of these interactions remains limited, highlighting the existence of topological constraints. Furthermore, contrary to intra-TAD interactions, the formation of inter-TAD enhancer-promoter interactions is not solely driven by genomic distance, with distal interactions sometimes favored over proximal ones. These observations suggest that other general mechanisms must exist to establish and maintain specific enhancer-promoter interactions across large distances.

What are tethering elements?

High-resolution Micro-C maps identified a specialized class of regulatory DNAs termed 'tethering elements' (TEs) in Drosophila. These 300-500-bp elements facilitate specific long-range genomic associations or loops. The POZ-containing transcription factor GAF (GAGA-associated factor) contributes to loop formation. Tether-tether interactions accelerate Hox gene activation by distal enhancers, and coordinate transcription of duplicated genes (paralogs) through promoter-promoter associations. Some TEs engage in ultra-long-range enhancer-promoter and promoter-promoter interactions (meta-loops) in the Drosophila brain. We discuss the basis for tether-tether specificity and speculate on the occurrence of similar elements in vertebrate genomes.

Spatial chromatin accessibility sequencing resolves high-order spatial interactions of epigenomic markers

As the genome is organized into a three-dimensional structure in intracellular space, epigenomic information also has a complex spatial arrangement. However, most epigenetic studies describe locations of methylation marks, chromatin accessibility regions, and histone modifications in the horizontal dimension. Proper spatial epigenomic information has rarely been obtained. In this study, we designed spatial chromatin accessibility sequencing (SCA-seq) to resolve the genome conformation by capturing the epigenetic information in single-molecular resolution while simultaneously resolving the genome conformation. Using SCA-seq, we are able to examine the spatial interaction of chromatin accessibility (e.g. enhancer-promoter contacts), CpG island methylation, and spatial insulating functions of the CCCTC-binding factor. We demonstrate that SCA-seq paves the way to explore the mechanism of epigenetic interactions and extends our knowledge in 3D packaging of DNA in the nucleus.

Sequence characteristics and an accurate model of abundant hyperactive loci in the human genome

Enhancers and promoters are classically considered to be bound by a small set of TFs in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with seemingly no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1,003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected 5 distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

Topology regulatory elements: From shaping genome architecture to gene regulation

The importance of 3D genome topology in the control of gene expression is becoming increasingly apparent, while regulatory mechanisms remain incompletely understood. Several recent studies have identified architectural elements that influence developmental gene expression by shaping locus topology. We refer to these elements as topological regulatory elements (TopoREs) to reflect their dual roles in genome organisation and gene expression. Importantly, these elements do not harbour autonomous transcriptional activation capacity, and instead appear to facilitate enhancer-promoter interactions, contributing to robust and precise timing of transcription. We discuss examples of TopoREs from two classes that are either dependent or independent of CTCF binding. Importantly, identification and interpretation of TopoRE function may shed light on multiple aspects of gene regulation, including the relationship between enhancer-promoter proximity and transcription, and enhancer-promoter specificity. Ultimately, understanding TopoRE diversity and function will aid in the interpretation of how human sequence variation can impact transcription and contribute to disease phenotypes.

Polycomb protein binding and looping mediated by Polycomb Response Elements in the ON transcriptional state

Polycomb group proteins (PcG) mediate epigenetic silencing of important developmental genes and other targets. In Drosophila, canonical PcG-target genes contain Polycomb Response Elements (PREs) that recruit PcG protein complexes including PRC2 that trimethylates H3K27 forming large H3K27me3 domains. In the OFF transcriptional state, PREs loop with each other and this looping strengthens silencing. Here we address the question of what PcG proteins bind to PREs when canonical PcG target genes are expressed, and whether PREs loop when these genes are ON. Our data show that the answer to this question is PRE-specific but general conclusions can be made. First, within a PcG-target gene, some regulatory DNA can remain covered with H3K27me3 and PcG proteins remain bound to PREs in these regions. Second, when PREs are within H3K27ac domains, PcG-binding decreases, however, this depends on the protein and PRE. The DNA binding protein GAF, and the PcG protein Ph remain at PREs even when other PcG proteins are greatly depleted. In the ON state, PREs can still loop with each other, but also form loops with presumptive enhancers. These data support the model that, in addition to their role in PcG silencing, PREs can act as "promoter-tethering elements" mediating interactions between promoter proximal PREs and distant enhancers.

Epromoters are new players in the regulatory landscape with potential pleiotropic roles

Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.

Chromosome-level organization of the regulatory genome in the Drosophila nervous system

Previous studies have identified topologically associating domains (TADs) as basic units of genome organization. We present evidence of a previously unreported level of genome folding, where distant TAD pairs, megabases apart, interact to form meta-domains. Within meta-domains, gene promoters and structural intergenic elements present in distant TADs are specifically paired. The associated genes encode neuronal determinants, including those engaged in axonal guidance and adhesion. These long-range associations occur in a large fraction of neurons but support transcription in only a subset of neurons. Meta-domains are formed by diverse transcription factors that are able to pair over long and flexible distances. We present evidence that two such factors, GAF and CTCF, play direct roles in this process. The relative simplicity of higher-order meta-domain interactions in Drosophila, compared with those previously described in mammals, allowed the demonstration that genomes can fold into highly specialized cell-type-specific scaffolds that enable megabase-scale regulatory associations.

Ultra-long-range interactions between active regulatory elements

Contacts between enhancers and promoters are thought to relate to their ability to activate transcription. Investigating factors that contribute to such chromatin interactions is therefore important for understanding gene regulation. Here, we have determined contact frequencies between millions of pairs of cis-regulatory elements from chromosome conformation capture data sets and analyzed a collection of hundreds of DNA-binding factors for binding at regions of enriched contacts. This analysis revealed enriched contacts at sites bound by many factors associated with active transcription. We show that active regulatory elements, independent of cohesin and polycomb, interact with each other across distances of tens of megabases in vertebrate and invertebrate genomes and that interactions correlate and change with activity. However, these ultra-long-range interactions are not dependent on RNA polymerase II transcription or individual transcription cofactors. Using simulations, we show that a model of chromatin and multivalent binding factors can give rise to long-range interactions via bridging-induced clustering. We propose that long-range interactions between cis-regulatory elements are driven by at least three distinct processes: cohesin-mediated loop extrusion, polycomb contacts, and clustering of active regions.

Profiling the immune epigenome across global cattle breeds

BACKGROUND

Understanding the variation between well and poorly adapted cattle breeds to local environments and pathogens is essential for breeding cattle with improved climate and disease-resistant phenotypes. Although considerable progress has been made towards identifying genetic differences between breeds, variation at the epigenetic and chromatin levels remains poorly characterized. Here, we generate, sequence and analyse over 150 libraries at base-pair resolution to explore the dynamics of DNA methylation and chromatin accessibility of the bovine immune system across three distinct cattle lineages.

RESULTS

We find extensive epigenetic divergence between the taurine and indicine cattle breeds across immune cell types, which is linked to the levels of local DNA sequence divergence between the two cattle sub-species. The unique cell type profiles enable the deconvolution of complex cellular mixtures using digital cytometry approaches. Finally, we show distinct sub-categories of CpG islands based on their chromatin and methylation profiles that discriminate between classes of distal and gene proximal islands linked to discrete transcriptional states.

CONCLUSIONS

Our study provides a comprehensive resource of DNA methylation, chromatin accessibility and RNA expression profiles of three diverse cattle populations. The findings have important implications, from understanding how genetic editing across breeds, and consequently regulatory backgrounds, may have distinct impacts to designing effective cattle epigenome-wide association studies in non-European breeds.

CpG island turnover events predict evolutionary changes in enhancer activity

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. Here we show that turnover of CpG islands (CGIs), which contribute to enhancer activation, is broadly associated with changes in enhancer activity across mammals, including humans. We integrated maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and found that CGI content in enhancers was strongly associated with increased histone modification levels. CGIs showed widespread turnover across species and species-specific CGIs were strongly enriched for enhancers exhibiting species-specific activity across all tissues and species we examined. Genes associated with enhancers with species-specific CGIs showed concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon. Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

Cis-regulatory atlas of primary human CD4+ T cells

Cis-regulatory elements (CRE) are critical for coordinating gene expression programs that dictate cell-specific differentiation and homeostasis. Recently developed self-transcribing active regulatory region sequencing (STARR-Seq) has allowed for genome-wide annotation of functional CREs. Despite this, STARR-Seq assays are only employed in cell lines, in part, due to difficulties in delivering reporter constructs. Herein, we implemented and validated a STARR-Seq-based screen in human CD4+ T cells using a non-integrating lentiviral transduction system. Lenti-STARR-Seq is the first example of a genome-wide assay of CRE function in human primary cells, identifying thousands of functional enhancers and negative regulatory elements (NREs) in human CD4+ T cells. We find an unexpected difference in nucleosome organization between enhancers and NRE: enhancers are located between nucleosomes, whereas NRE are occupied by nucleosomes in their endogenous locations. We also describe chromatin modification, eRNA production, and transcription factor binding at both enhancers and NREs. Our findings support the idea of silencer repurposing as enhancers in alternate cell types. Collectively, these data suggest that Lenti-STARR-Seq is a successful approach for CRE screening in primary human cell types, and provides an atlas of functional CREs in human CD4+ T 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.

Evolution of 3D chromatin organization at different scales

Most animal genomes fold in 3D chromatin domains called topologically associated domains (TADs) that facilitate interactions between cis-regulatory elements (CREs) and promoters. Owing to their critical role in the control of developmental gene expression, we explore how TADs have shaped animal evolution. In the light of recent studies that profile TADs in disparate animal lineages, we discuss their phylogenetic distribution and the mechanisms that underlie their formation. We present evidence indicating that TADs are plastic entities composed of genomic strata of different ages: ancient cores are combined with newer regions and brought into extant TADs through genomic rearrangements. We highlight that newly incorporated TAD strata enable the establishment of new CRE-promoter interactions and in turn new expression patterns that can drive phenotypical innovation. We further highlight how subtle changes in chromatin folding may fine-tune the expression levels of developmental genes and hold a potential for evolutionary significance.

The histone acetyltransferase KAT6A is recruited to unmethylated CpG islands via a DNA binding winged helix domain

The lysine acetyltransferase KAT6A (MOZ, MYST3) belongs to the MYST family of chromatin regulators, facilitating histone acetylation. Dysregulation of KAT6A has been implicated in developmental syndromes and the onset of acute myeloid leukemia (AML). Previous work suggests that KAT6A is recruited to its genomic targets by a combinatorial function of histone binding PHD fingers, transcription factors and chromatin binding interaction partners. Here, we demonstrate that a winged helix (WH) domain at the very N-terminus of KAT6A specifically interacts with unmethylated CpG motifs. This DNA binding function leads to the association of KAT6A with unmethylated CpG islands (CGIs) genome-wide. Mutation of the essential amino acids for DNA binding completely abrogates the enrichment of KAT6A at CGIs. In contrast, deletion of a second WH domain or the histone tail binding PHD fingers only subtly influences the binding of KAT6A to CGIs. Overexpression of a KAT6A WH1 mutant has a dominant negative effect on H3K9 histone acetylation, which is comparable to the effects upon overexpression of a KAT6A HAT domain mutant. Taken together, our work revealed a previously unrecognized chromatin recruitment mechanism of KAT6A, offering a new perspective on the role of KAT6A in gene regulation and human diseases.

Changes in PRC1 activity during interphase modulate lineage transition in pluripotent cells

The potential of pluripotent cells to respond to developmental cues and trigger cell differentiation is enhanced during the G1 phase of the cell cycle, but the molecular mechanisms involved are poorly understood. Variations in polycomb activity during interphase progression have been hypothesized to regulate the cell-cycle-phase-dependent transcriptional activation of differentiation genes during lineage transition in pluripotent cells. Here, we show that recruitment of Polycomb Repressive Complex 1 (PRC1) and associated molecular functions, ubiquitination of H2AK119 and three-dimensional chromatin interactions, are enhanced during S and G2 phases compared to the G1 phase. In agreement with the accumulation of PRC1 at target promoters upon G1 phase exit, cells in S and G2 phases show firmer transcriptional repression of developmental regulator genes that is drastically perturbed upon genetic ablation of the PRC1 catalytic subunit RING1B. Importantly, depletion of RING1B during retinoic acid stimulation interferes with the preference of mouse embryonic stem cells (mESCs) to induce the transcriptional activation of differentiation genes in G1 phase. We propose that incremental enrolment of polycomb repressive activity during interphase progression reduces the tendency of cells to respond to developmental cues during S and G2 phases, facilitating activation of cell differentiation in the G1 phase of the pluripotent cell cycle.

Sensitive Interrogation of Enhancer Activity in Living Cells on a Nanoelectroporation-Probing Platform

Enhancers involved in the upregulation of multiple oncogenes play a fundamental role in tumorigenesis and immortalization. Exploring the activity of enhancers in living cells has emerged as a critical path to a deep understanding of cancer properties, further providing important clues to targeted therapy. However, identifying enhancer activity in living cells is challenging due to the double biological barriers of a cell cytoplasmic membrane and a nuclear membrane, limiting the sensitivity and responsiveness of conventional probing methods. In this work, we developed a nanoelectroporation-probing (NP) platform, which enables intranuclear probe delivery for sensitive interrogation of enhancer activity in living cells. The nanoelectroporation biochip achieved highly focused perforation of the cell cytoplasmic membrane and brought about additional driving force to expedite the delivery of probes into the nucleus. The probes targeting enhancer activity (named "PH probe") are programmed with a cyclic amplification strategy and enable an increase in the fluorescence signals over 100-fold within 1 h. The platform was leveraged to detect the activity of CCAT1 enhancers (CCAT1, colon cancer-associated transcript-1, a long noncoding RNA that functions in tumor invasion and metastasis) in cell samples from clinical lung cancer patients, as well as reveal the heterogeneity of enhancers among different patients. The observations may extend the linkages between enhancers and cancer cells while validating the robustness and reliability of the platform for the assay of enhancer activity. This platform will be a promising toolbox with wide applicable potential for the intranuclear study of living cells.

RINGs, DUBs and Abnormal Brain Growth-Histone H2A Ubiquitination in Brain Development and Disease

During mammalian neurodevelopment, signaling pathways converge upon transcription factors (TFs) to establish appropriate gene expression programmes leading to the production of distinct neural and glial cell types. This process is partially regulated by the dynamic modulation of chromatin states by epigenetic systems, including the polycomb group (PcG) family of co-repressors. PcG proteins form multi-subunit assemblies that sub-divide into distinct, yet functionally related families. Polycomb repressive complexes 1 and 2 (PRC1 and 2) modify the chemical properties of chromatin by covalently modifying histone tails via H2A ubiquitination (H2AK119ub1) and H3 methylation, respectively. In contrast to the PRCs, the Polycomb repressive deubiquitinase (PR-DUB) complex removes H2AK119ub1 from chromatin through the action of the C-terminal hydrolase BAP1. Genetic screening has identified several PcG mutations that are causally associated with a range of congenital neuropathologies associated with both localised and/or systemic growth abnormalities. As PRC1 and PR-DUB hold opposing functions to control H2AK119ub1 levels across the genome, it is plausible that such neurodevelopmental disorders arise through a common mechanism. In this review, we will focus on advancements regarding the composition and opposing molecular functions of mammalian PRC1 and PR-DUB, and explore how their dysfunction contributes to the emergence of neurodevelopmental disorders.

Cell fate decisions, transcription factors and signaling during early retinal development

The development of the vertebrate eyes is a complex process starting from anterior-posterior and dorso-ventral patterning of the anterior neural tube, resulting in the formation of the eye field. Symmetrical separation of the eye field at the anterior neural plate is followed by two symmetrical evaginations to generate a pair of optic vesicles. Next, reciprocal invagination of the optic vesicles with surface ectoderm-derived lens placodes generates double-layered optic cups. The inner and outer layers of the optic cups develop into the neural retina and retinal pigment epithelium (RPE), respectively. In vitro produced retinal tissues, called retinal organoids, are formed from human pluripotent stem cells, mimicking major steps of retinal differentiation in vivo. This review article summarizes recent progress in our understanding of early eye development, focusing on the formation the eye field, optic vesicles, and early optic cups. Recent single-cell transcriptomic studies are integrated with classical in vivo genetic and functional studies to uncover a range of cellular mechanisms underlying early eye development. The functions of signal transduction pathways and lineage-specific DNA-binding transcription factors are dissected to explain cell-specific regulatory mechanisms underlying cell fate determination during early eye development. The functions of homeodomain (HD) transcription factors Otx2, Pax6, Lhx2, Six3 and Six6, which are required for early eye development, are discussed in detail. Comprehensive understanding of the mechanisms of early eye development provides insight into the molecular and cellular basis of developmental ocular anomalies, such as optic cup coloboma. Lastly, modeling human development and inherited retinal diseases using stem cell-derived retinal organoids generates opportunities to discover novel therapies for retinal diseases.

The Mediator complex as a master regulator of transcription by RNA polymerase II

The Mediator complex, which in humans is 1.4 MDa in size and includes 26 subunits, controls many aspects of RNA polymerase II (Pol II) function. Apart from its size, a defining feature of Mediator is its intrinsic disorder and conformational flexibility, which contributes to its ability to undergo phase separation and to interact with a myriad of regulatory factors. In this Review, we discuss Mediator structure and function, with emphasis on recent cryogenic electron microscopy data of the 4.0-MDa transcription preinitiation complex. We further discuss how Mediator and sequence-specific DNA-binding transcription factors enable enhancer-dependent regulation of Pol II function at distal gene promoters, through the formation of molecular condensates (or transcription hubs) and chromatin loops. Mediator regulation of Pol II reinitiation is also discussed, in the context of transcription bursting. We propose a working model for Mediator function that combines experimental results and theoretical considerations related to enhancer-promoter interactions, which reconciles contradictory data regarding whether enhancer-promoter communication is direct or indirect. We conclude with a discussion of Mediator's potential as a therapeutic target and of future research directions.

Distinct roles for CKM-Mediator in controlling Polycomb-dependent chromosomal interactions and priming genes for induction

Precise control of gene expression underpins normal development. This relies on mechanisms that enable communication between gene promoters and other regulatory elements. In embryonic stem cells (ESCs), the cyclin-dependent kinase module Mediator complex (CKM-Mediator) has been reported to physically link gene regulatory elements to enable gene expression and also prime genes for induction during differentiation. Here, we show that CKM-Mediator contributes little to three-dimensional genome organization in ESCs, but it has a specific and essential role in controlling interactions between inactive gene regulatory elements bound by Polycomb repressive complexes (PRCs). These interactions are established by the canonical PRC1 (cPRC1) complex but rely on CKM-Mediator, which facilitates binding of cPRC1 to its target sites. Importantly, through separation-of-function experiments, we reveal that this collaboration between CKM-Mediator and cPRC1 in creating long-range interactions does not function to prime genes for induction during differentiation. Instead, we discover that priming relies on an interaction-independent mechanism whereby the CKM supports core Mediator engagement with gene promoters during differentiation to enable gene activation.

DNA Methylation: Genomewide Distribution, Regulatory Mechanism and Therapy Target

DNA methylation is the most important epigenetic modification involved in the regulation of transcription, imprinting, establishment of X-inactivation, and the formation of a chromatin structure. DNA methylation in the genome is often associated with transcriptional repression and the formation of closed heterochromatin. However, the results of genome-wide studies of the DNA methylation pattern and transcriptional activity of genes have nudged us toward reconsidering this paradigm, since the promoters of many genes remain active despite their methylation. The differences in the DNA methylation distribution in normal and pathological conditions allow us to consider methylation as a diagnostic marker or a therapy target. In this regard, the need to investigate the factors affecting DNA methylation and those involved in its interpretation becomes pressing. Recently, a large number of protein factors have been uncovered, whose ability to bind to DNA depends on their methylation. Many of these proteins act not only as transcriptional activators or repressors, but also affect the level of DNA methylation. These factors are considered potential therapeutic targets for the treatment of diseases resulting from either a change in DNA methylation or a change in the interpretation of its methylation level. In addition to protein factors, a secondary DNA structure can also affect its methylation and can be considered as a therapy target. In this review, the latest research into the DNA methylation landscape in the genome has been summarized to discuss why some DNA regions avoid methylation and what factors can affect its level or interpretation and, therefore, can be considered a therapy target.

Know when to fold 'em: Polycomb complexes in oncogenic 3D genome regulation

Chromatin is spatially and temporally regulated through a series of orchestrated processes resulting in the formation of 3D chromatin structures such as topologically associating domains (TADs), loops and Polycomb Bodies. These structures are closely linked to transcriptional regulation, with loss of control of these processes a frequent feature of cancer and developmental syndromes. One such oncogenic disruption of the 3D genome is through recurrent dysregulation of Polycomb Group Complex (PcG) functions either through genetic mutations, amplification or deletion of genes that encode for PcG proteins. PcG complexes are evolutionarily conserved epigenetic complexes. They are key for early development and are essential transcriptional repressors. PcG complexes include PRC1, PRC2 and PR-DUB which are responsible for the control of the histone modifications H2AK119ub1 and H3K27me3. The spatial distribution of the complexes within the nuclear environment, and their associated modifications have profound effects on the regulation of gene transcription and the 3D genome. Nevertheless, how PcG complexes regulate 3D chromatin organization is still poorly understood. Here we glean insights into the role of PcG complexes in 3D genome regulation and compaction, how these processes go awry during tumorigenesis and the therapeutic implications that result from our insights into these mechanisms.

Enhancer-gene specificity in development and disease

Enhancers control the establishment of spatiotemporal gene expression patterns throughout development. Over the past decade, the development of new technologies has improved our capacity to link enhancers with their target genes based on their colocalization within the same topological domains. However, the mechanisms that regulate how enhancers specifically activate some genes but not others within a given domain remain unclear. In this Review, we discuss recent insights into the factors controlling enhancer specificity, including the genetic composition of enhancers and promoters, the linear and 3D distance between enhancers and their target genes, and cell-type specific chromatin landscapes. We also discuss how elucidating the molecular principles of enhancer specificity might help us to better understand and predict the pathological consequences of human genetic, epigenetic and structural variants.

Deep-learning microscopy image reconstruction with quality control reveals second-scale rearrangements in RNA polymerase II clusters

Fluorescence microscopy, a central tool of biological research, is subject to inherent trade-offs in experiment design. For instance, image acquisition speed can only be increased in exchange for a lowered signal quality, or for an increased rate of photo-damage to the specimen. Computational denoising can recover some loss of signal, extending the trade-off margin for high-speed imaging. Recently proposed denoising on the basis of neural networks shows exceptional performance but raises concerns of errors typical of neural networks. Here, we present a work-flow that supports an empirically optimized reduction of exposure times, as well as per-image quality control to exclude images with reconstruction errors. We implement this work-flow on the basis of the denoising tool Noise2Void and assess the molecular state and 3D shape of RNA polymerase II (Pol II) clusters in live zebrafish embryos. Image acquisition speed could be tripled, achieving 2-s time resolution and 350-nm lateral image resolution. The obtained data reveal stereotyped events of approximately 10 s duration: initially, the molecular mark for recruited Pol II increases, then the mark for active Pol II increases, and finally Pol II clusters take on a stretched and unfolded shape. An independent analysis based on fixed sample images reproduces this sequence of events, and suggests that they are related to the transient association of genes with Pol II clusters. Our work-flow consists of procedures that can be implemented on commercial fluorescence microscopes without any hardware or software modification, and should, therefore, be transferable to many other applications.

Transcriptional coupling of distant regulatory genes in living embryos

The prevailing view of metazoan gene regulation is that individual genes are independently regulated by their own dedicated sets of transcriptional enhancers. Past studies have reported long-range gene-gene associations^1-3, but their functional importance in regulating transcription remains unclear. Here we used quantitative single-cell live imaging methods to provide a demonstration of co-dependent transcriptional dynamics of genes separated by large genomic distances in living Drosophila embryos. We find extensive physical and functional associations of distant paralogous genes, including co-regulation by shared enhancers and co-transcriptional initiation over distances of nearly 250 kilobases. Regulatory interconnectivity depends on promoter-proximal tethering elements, and perturbations in these elements uncouple transcription and alter the bursting dynamics of distant genes, suggesting a role of genome topology in the formation and stability of co-transcriptional hubs. Transcriptional coupling is detected throughout the fly genome and encompasses a broad spectrum of conserved developmental processes, suggesting a general strategy for long-range integration of gene activity.

DNA binding by polycomb-group proteins: searching for the link to CpG islands

Polycomb group proteins predominantly exist in polycomb repressive complexes (PRCs) that cooperate to maintain the repressed state of thousands of cell-type-specific genes. Targeting PRCs to the correct sites in chromatin is essential for their function. However, the mechanisms by which PRCs are recruited to their target genes in mammals are multifactorial and complex. Here we review DNA binding by polycomb group proteins. There is strong evidence that the DNA-binding subunits of PRCs and their DNA-binding activities are required for chromatin binding and CpG targeting in cells. In vitro, CpG-specific binding was observed for truncated proteins externally to the context of their PRCs. Yet, the mere DNA sequence cannot fully explain the subset of CpG islands that are targeted by PRCs in any given cell type. At this time we find very little structural and biophysical evidence to support a model where sequence-specific DNA-binding activity is required or sufficient for the targeting of CpG-dinucleotide sequences by polycomb group proteins while they are within the context of their respective PRCs, either PRC1 or PRC2. We discuss the current knowledge and open questions on how the DNA-binding activities of polycomb group proteins facilitate the targeting of PRCs to chromatin.

Intragenic CpG Islands and Their Impact on Gene Regulation

The mammalian genome is depleted in CG dinucleotides, except at protected regions where they cluster as CpG islands (CGIs). CGIs are gene regulatory hubs and serve as transcription initiation sites and are as expected, associated with gene promoters. Advances in genomic annotations demonstrate that a quarter of CGIs are found within genes. Such intragenic regions are repressive environments, so it is surprising that CGIs reside here and even more surprising that some resist repression and are transcriptionally active within a gene. Hence, intragenic CGI positioning within genes is not arbitrary and is instead, selected for. As a wealth of recent studies demonstrate, intragenic CGIs are embedded within genes and consequently, influence ‘host’ gene mRNA isoform length and expand transcriptome diversity.

Genome organization controls transcriptional dynamics during development

Past studies offer contradictory claims for the role of genome organization in the regulation of gene activity. Here, we show through high-resolution chromosome conformation analysis that the Drosophila genome is organized by two independent classes of regulatory sequences, tethering elements and insulators. Quantitative live imaging and targeted genome editing demonstrate that this two-tiered organization is critical for the precise temporal dynamics of Hox gene transcription during development. Tethering elements mediate long-range enhancer-promoter interactions and foster fast activation kinetics. Conversely, the boundaries of topologically associating domains (TADs) prevent spurious interactions with enhancers and silencers located in neighboring TADs. These two levels of genome organization operate independently of one another to ensure precision of transcriptional dynamics and the reliability of complex patterning processes.

MYC overexpression leads to increased chromatin interactions at superenhancers and MYC binding sites

The MYC oncogene encodes for the MYC protein and is frequently dysregulated across multiple cancer cell types, making it an attractive target for cancer therapy. MYC overexpression leads to MYC binding at active enhancers, resulting in a global transcriptional amplification of active genes. Because super-enhancers are frequently dysregulated in cancer, we hypothesized that MYC preferentially invades into super-enhancers and alters the cancer genome organization. To that end, we performed ChIP-seq, RNA-seq, circular chromosome conformation capture (4C-seq), and Spike-in Quantitative Hi-C (SIQHiC) on the U2OS osteosarcoma cell line with tetracycline-inducible MYC. MYC overexpression in U2OS cells modulated histone acetylation and increased MYC binding at super-enhancers. SIQHiC analysis revealed increased global chromatin contact frequency, particularly at chromatin interactions connecting MYC binding sites at promoters and enhancers. Immunofluorescence staining showed that MYC molecules formed punctate foci at these transcriptionally active domains after MYC overexpression. These results demonstrate the accumulation of overexpressed MYC at promoter–enhancer hubs and suggest that MYC invades into enhancers through spatial proximity. At the same time, the increased protein–protein interactions may strengthen these chromatin interactions to increase chromatin contact frequency. CTCF siRNA knockdown in MYC-overexpressed U2OS cells demonstrated that removal of architectural proteins can disperse MYC and abrogate the increase in chromatin contacts. By elucidating the chromatin landscape of MYC-driven cancers, we can potentially target MYC-associated chromatin interactions for cancer therapy.

Metabolic and Epigenetic Regulation by Estrogen in Adipocytes

Sex hormones contribute to differences between males and females in body fat distribution and associated disease risk. Higher concentrations of estrogens are associated with a more gynoid body shape and with more fat storage on hips and thighs rather than in visceral depots. Estrogen-mediated protection against visceral adiposity is shown in post-menopausal women with lower levels of estrogens and the reduction in central body fat observed after treatment with hormone-replacement therapy. Estrogen exerts its physiological effects via the estrogen receptors (ERα, ERβ and GPR30) in target cells, including adipocytes. Studies in mice indicate that estrogen protects against adipose inflammation and fibrosis also before the onset of obesity. The mechanisms involved in estrogen-dependent body fat distribution are incompletely understood, but involve, e.g., increased mTOR signaling and suppression of autophagy and adipogenesis/lipid storage. Estrogen plays a key role in epigenetic regulation of adipogenic genes by interacting with enzymes that remodel DNA methylation and histone tail post-translational modifications. However, more studies are needed to map the differential epigenetic effects of ER in different adipocyte subtypes, including those in subcutaneous and visceral adipose tissues. We here review recent discoveries of ER-mediated transcriptional and epigenetic regulation in adipocytes, which may explain sexual dimorphisms in body fat distribution and obesity-related disease risk.

The molecular principles of gene regulation by Polycomb repressive complexes

Precise control of gene expression is fundamental to cell function and development. Although ultimately gene expression relies on DNA-binding transcription factors to guide the activity of the transcription machinery to genes, it has also become clear that chromatin and histone post-translational modification have fundamental roles in gene regulation. Polycomb repressive complexes represent a paradigm of chromatin-based gene regulation in animals. The Polycomb repressive system comprises two central protein complexes, Polycomb repressive complex 1 (PRC1) and PRC2, which are essential for normal gene regulation and development. Our early understanding of Polycomb function relied on studies in simple model organisms, but more recently it has become apparent that this system has expanded and diverged in mammals. Detailed studies are now uncovering the molecular mechanisms that enable mammalian PRC1 and PRC2 to identify their target sites in the genome, communicate through feedback mechanisms to create Polycomb chromatin domains and control transcription to regulate gene expression. In this Review, we discuss and contextualize the emerging principles that define how this fascinating chromatin-based system regulates gene expression in mammals.

A bipartite element with allele-specific functions safeguards DNA methylation imprints at the Dlk1-Dio3 locus

Loss of imprinting (LOI) results in severe developmental defects, but the mechanisms preventing LOI remain incompletely understood. Here, we dissect the functional components of the imprinting control region of the essential Dlk1-Dio3 locus (called IG-DMR) in pluripotent stem cells. We demonstrate that the IG-DMR consists of two antagonistic elements: a paternally methylated CpG island that prevents recruitment of TET dioxygenases and a maternally unmethylated non-canonical enhancer that ensures expression of the Gtl2 lncRNA by counteracting de novo DNA methyltransferases. Genetic or epigenetic editing of these elements leads to distinct LOI phenotypes with characteristic alternations of allele-specific gene expression, DNA methylation, and 3D chromatin topology. Although repression of the Gtl2 promoter results in dysregulated imprinting, the stability of LOI phenotypes depends on the IG-DMR, suggesting a functional hierarchy. These findings establish the IG-DMR as a bipartite control element that maintains imprinting by allele-specific restriction of the DNA (de)methylation machinery.

Transcriptional enhancers and their communication with gene promoters

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

The chromatin, topological and regulatory properties of pluripotency-associated poised enhancers are conserved in vivo

Poised enhancers (PEs) represent a genetically distinct set of distal regulatory elements that control the expression of major developmental genes. Before becoming activated in differentiating cells, PEs are already bookmarked in pluripotent cells with unique chromatin and topological features that could contribute to their privileged regulatory properties. However, since PEs were originally characterized in embryonic stem cells (ESC), it is currently unknown whether PEs are functionally conserved in vivo. Here, we show that the chromatin and 3D structural features of PEs are conserved among mouse pluripotent cells both in vitro and in vivo. We also uncovered that the interactions between PEs and their target genes are globally controlled by the combined action of Polycomb, Trithorax and architectural proteins. Moreover, distal regulatory sequences located close to developmental genes and displaying the typical genetic (i.e. CpG islands) and chromatin (i.e. high accessibility and H3K27me3 levels) features of PEs are commonly found across vertebrates. These putative PEs show high sequence conservation within specific vertebrate clades, with only a few being evolutionary conserved across all vertebrates. Lastly, by genetically disrupting PEs in mouse and chicken embryos, we demonstrate that these regulatory elements play essential roles during the induction of major developmental genes in vivo.