Long-Range Enhancer Interactions Are Prevalent in Mouse Embryonic Stem Cells and Are Reorganized upon Pluripotent State Transition

The SWI/SNF ATP-dependent chromatin remodeling complex in cell lineage priming and early development

ATP dependent chromatin remodelers have pivotal roles in transcription, DNA replication and repair, and maintaining genome integrity. SWI/SNF remodelers were first discovered in yeast genetic screens for factors involved in mating type switching or for using alternative energy sources therefore termed SWI/SNF complex (short for SWItch/Sucrose NonFermentable). The SWI/SNF complexes utilize energy from ATP hydrolysis to disrupt histone-DNA interactions and shift, eject, or reposition nucleosomes making the underlying DNA more accessible to specific transcription factors and other regulatory proteins. In development, SWI/SNF orchestrates the precise activation and repression of genes at different stages, safe guards the formation of specific cell lineages and tissues. Dysregulation of SWI/SNF have been implicated in diseases such as cancer, where they can drive uncontrolled cell proliferation and tumor metastasis. Additionally, SWI/SNF defects are associated with neurodevelopmental disorders, leading to disruption of neural development and function. This review offers insights into recent developments regarding the roles of the SWI/SNF complex in pluripotency and cell lineage primining and the approaches that have helped delineate its importance. Understanding these molecular mechanisms is crucial for unraveling the intricate processes governing embryonic stem cell biology and developmental transitions and may potentially apply to human diseases linked to mutations in the SWI/SNF complex.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BACKGROUND

Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.

RESULTS

Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.

CONCLUSIONS

Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

Genetic pathways regulating the longitudinal acquisition of cocaine self-administration in a panel of inbred and recombinant inbred mice

To identify addiction genes, we evaluate intravenous self-administration of cocaine or saline in 84 inbred and recombinant inbred mouse strains over 10 days. We integrate the behavior data with brain RNA-seq data from 41 strains. The self-administration of cocaine and that of saline are genetically distinct. We maximize power to map loci for cocaine intake by using a linear mixed model to account for this longitudinal phenotype while correcting for population structure. A total of 15 unique significant loci are identified in the genome-wide association study. A transcriptome-wide association study highlights the Trpv2 ion channel as a key locus for cocaine self-administration as well as identifying 17 additional genes, including Arhgef26, Slc18b1, and Slco5a1. We find numerous instances where alternate splice site selection or RNA editing altered transcript abundance. Our work emphasizes the importance of Trpv2, an ionotropic cannabinoid receptor, for the response to cocaine.

Super-enhancers complexes zoom in transcription in cancer

Super-enhancers (SEs) consist of multiple typical enhancers enriched at high density with transcription factors, histone-modifying enzymes and cofactors. Oncogenic SEs promote tumorigenesis and malignancy by altering protein-coding gene expression and noncoding regulatory element function. Therefore, they play central roles in the treatment of cancer. Here, we review the structural characteristics, organization, identification, and functions of SEs and the underlying molecular mechanism by which SEs drive oncogenic transcription in tumor cells. We then summarize abnormal SE complexes, SE-driven coding genes, and noncoding RNAs involved in tumor development. In summary, we believe that SEs show great potential as biomarkers and therapeutic targets.

A genome-wide relay of signalling-responsive enhancers drives hematopoietic specification

Developmental control of gene expression critically depends on distal cis-regulatory elements including enhancers which interact with promoters to activate gene expression. To date no global experiments have been conducted that identify their cell type and cell stage-specific activity within one developmental pathway and in a chromatin context. Here, we describe a high-throughput method that identifies thousands of differentially active cis-elements able to stimulate a minimal promoter at five stages of hematopoietic progenitor development from embryonic stem (ES) cells, which can be adapted to any ES cell derived cell type. We show that blood cell-specific gene expression is controlled by the concerted action of thousands of differentiation stage-specific sets of cis-elements which respond to cytokine signals terminating at signalling responsive transcription factors. Our work provides an important resource for studies of hematopoietic specification and highlights the mechanisms of how and where extrinsic signals program a cell type-specific chromatin landscape driving hematopoietic differentiation.

Application of the 3C Method to Study the Developmental Genes in Drosophila Larvae

A transition from one developmental stage to another is accompanied by activation of developmental programs and corresponding gene ensembles. Changes in the spatial conformation of the corresponding loci are associated with this activation and can be investigated with the help of the Chromosome Conformation Capture (3C) methodology. Application of 3C to specific developmental stages is a sophisticated task. Here, we describe the use of the 3C method to study the spatial organization of developmental loci in Drosophila larvae. We critically analyzed the existing protocols and offered our own solutions and the optimized protocol to overcome limitations. To demonstrate the efficiency of our procedure, we studied the spatial organization of the developmental locus Dad in 3rd instar Drosophila larvae. Differences in locus conformation were found between embryonic cells and living wild-type larvae. We also observed the establishment of novel regulatory interactions in the presence of an adjacent transgene upon activation of its expression in larvae. Our work fills the gap in the application of the 3C method to Drosophila larvae and provides a useful guide for establishing 3C on an animal model.

SOX2 transcription factor binding and function

The transcription factor SOX2 is a vital regulator of stem cell activity in various developing and adult tissues. Mounting evidence has demonstrated the importance of SOX2 in regulating the induction and maintenance of stemness as well as in controlling cell proliferation, lineage decisions and differentiation. Recent studies have revealed that the ability of SOX2 to regulate these stem cell features involves its function as a pioneer factor, with the capacity to target nucleosomal DNA, modulate chromatin accessibility and prepare silent genes for subsequent activation. Moreover, although SOX2 binds to similar DNA motifs in different stem cells, its multifaceted and cell type-specific functions are reliant on context-dependent features. These cell type-specific properties include variations in partner factor availability and SOX2 protein expression levels. In this Primer, we discuss recent findings that have increased our understanding of how SOX2 executes its versatile functions as a master regulator of stem cell activities.

Super enhancers as master gene regulators in the pathogenesis of hematologic malignancies

Transcriptional deregulation of multiple oncogenes, tumor suppressors and survival pathways is a cancer cell hallmark. Super enhancers (SE) are long stretches of active enhancers in close linear proximity that ensure extraordinarily high expression levels of key genes associated with cell lineage, function and survival. SE landscape is intrinsically prone to changes and reorganization during the course of normal cell differentiation. This functional plasticity is typically utilized by cancer cells, which remodel their SE landscapes to ensure oncogenic transcriptional reprogramming. Multiple recent studies highlighted structural genetic mechanisms in non-coding regions that create new SE or hijack already existing ones. In addition, alterations in abundance/activity of certain SE-associated proteins or certain viral infections can elicit new super enhancers and trigger SE-driven transcriptional changes. For these reasons, SE profiling emerged as a powerful tool for discovering the core transcriptional regulatory circuits in tumor cells. This, in turn, provides new insights into cancer cell biology, and identifies main nodes of key cellular pathways to be potentially targeted. Since SEs are susceptible to inhibition, their disruption results in exponentially amassing 'butterfly' effect on gene expression and cell function. Moreover, many of SE elements are druggable, opening new therapeutic opportunities. Indeed, SE targeting drugs have been studied preclinically in various hematologic malignancies with promising effects. Herein, we review the unique features of SEs, present different cis- and trans-acting mechanisms through which hematologic tumor cells acquire SEs, and finally, discuss the potential of SE targeting in the therapy of hematologic malignancies.

Genetic control of the pluripotency epigenome determines differentiation bias in mouse embryonic stem cells

Genetically diverse pluripotent stem cells display varied, heritable responses to differentiation cues. Here, we harnessed these disparities through derivation of mouse embryonic stem cells from the BXD genetic reference panel, along with C57BL/6J (B6) and DBA/2J (D2) parental strains, to identify loci regulating cell state transitions. Upon transition to formative pluripotency, B6 stem cells quickly dissolved naïve networks adopting gene expression modules indicative of neuroectoderm lineages, whereas D2 retained aspects of naïve pluripotency. Spontaneous formation of embryoid bodies identified divergent differentiation where B6 showed a propensity toward neuroectoderm and D2 toward definitive endoderm. Genetic mapping identified major trans-acting loci co-regulating chromatin accessibility and gene expression in both naïve and formative pluripotency. These loci distally modulated occupancy of pluripotency factors at hundreds of regulatory elements. One trans-acting locus on Chr 12 primarily impacted chromatin accessibility in embryonic stem cells, while in epiblast-like cells, the same locus subsequently influenced expression of genes enriched for neurogenesis, suggesting early chromatin priming. These results demonstrate genetically determined biases in lineage commitment and identify major regulators of the pluripotency epigenome.

3D genome organization links non-coding disease-associated variants to genes

Genome sequencing has revealed over 300 million genetic variations in human populations. Over 90% of variants are single nucleotide polymorphisms (SNPs), the remainder include short deletions or insertions, and small numbers of structural variants. Hundreds of thousands of these variants have been associated with specific phenotypic traits and diseases through genome wide association studies which link significant differences in variant frequencies with specific phenotypes among large groups of individuals. Only 5% of disease-associated SNPs are located in gene coding sequences, with the potential to disrupt gene expression or alter of the function of encoded proteins. The remaining 95% of disease-associated SNPs are located in non-coding DNA sequences which make up 98% of the genome. The role of non-coding, disease-associated SNPs, many of which are located at considerable distances from any gene, was at first a mystery until the discovery that gene promoters regularly interact with distal regulatory elements to control gene expression. Disease-associated SNPs are enriched at the millions of gene regulatory elements that are dispersed throughout the non-coding sequences of the genome, suggesting they function as gene regulation variants. Assigning specific regulatory elements to the genes they control is not straightforward since they can be millions of base pairs apart. In this review we describe how understanding 3D genome organization can identify specific interactions between gene promoters and distal regulatory elements and how 3D genomics can link disease-associated SNPs to their target genes. Understanding which gene or genes contribute to a specific disease is the first step in designing rational therapeutic interventions.

OUP accepted manuscript

Super enhancers (SEs) are broad enhancer domains usually containing multiple constituent enhancers that hold elevated activities in gene regulation. Disruption in one or more constituent enhancers causes aberrant SE activities that lead to gene dysregulation in diseases. To quantify SE aberrations, differential analysis is performed to compare SE activities between cell conditions. The state-of-art strategy in estimating differential SEs relies on overall activities and neglect the changes in length and structure of SEs. Here, we propose a novel computational method to identify differential SEs by weighting the combinatorial effects of constituent-enhancer activities and locations (i.e. internal dynamics). In addition to overall activity changes, our method identified four novel classes of differential SEs with distinct enhancer structural alterations. We demonstrate that these structure alterations hold distinct regulatory impact, such as regulating different number of genes and modulating gene expression with different strengths, highlighting the differentiated regulatory roles of these unexplored SE features. When compared to the existing method, our method showed improved identification of differential SEs that were linked to better discernment of cell-type-specific SE activity and functional interpretation.

Distinct transcription kinetics of pluripotent cell states

Mouse embryonic stem cells (mESCs) can adopt naïve, ground, and paused pluripotent states that give rise to unique transcriptomes. Here, we use transient transcriptome sequencing (TT-seq) to define both coding and non-coding transcription units (TUs) in these three pluripotent states and combine TT-seq with RNA polymerase II occupancy profiling to unravel the kinetics of RNA metabolism genome-wide. Compared to the naïve state (serum), RNA synthesis and turnover rates are globally reduced in the ground state (2i) and the paused state (mTORi). The global reduction in RNA synthesis goes along with a genome-wide decrease of polymerase elongation velocity, which is related to epigenomic features and alterations in the Pol II termination window. Our data suggest that transcription activity is the main determinant of steady state mRNA levels in the naïve state and that genome-wide changes in transcription kinetics invoke ground and paused pluripotent states.

Single-cell measurement of higher-order 3D genome organization with scSPRITE

Although three-dimensional (3D) genome organization is central to many aspects of nuclear function, it has been difficult to measure at the single-cell level. To address this, we developed 'single-cell split-pool recognition of interactions by tag extension' (scSPRITE). scSPRITE uses split-and-pool barcoding to tag DNA fragments in the same nucleus and their 3D spatial arrangement. Because scSPRITE measures multiway DNA contacts, it generates higher-resolution maps within an individual cell than can be achieved by proximity ligation. We applied scSPRITE to thousands of mouse embryonic stem cells and detected known genome structures, including chromosome territories, active and inactive compartments, and topologically associating domains (TADs) as well as long-range inter-chromosomal structures organized around various nuclear bodies. We observe that these structures exhibit different levels of heterogeneity across the population, with TADs representing dynamic units of genome organization across cells. We expect that scSPRITE will be a critical tool for studying genome structure within heterogeneous populations.

Enhancer-Promoter Communication: It's Not Just About Contact

Cis-regulatory elements such as enhancers can be located even a million base pairs away from their cognate promoter and yet modulate gene transcription. Indeed, the 3D organisation of chromatin enables the establishment of long-range enhancer-promoter communication. The observation of long-range enhancer-promoter chromatin loops at active genes originally led to a model in which enhancers and promoters form physical contacts between each other to control transcription. Yet, recent microscopy data has challenged this prevailing activity-by-contact model of enhancer-promoter communication in transcriptional activation. Live single-cell imaging approaches do not systematically reveal a correlation between enhancer-proximity and transcriptional activation. We therefore discuss the need to move from a static to a dynamic view of enhancer-promoter relationships. We highlight recent studies that not only reveal considerable chromatin movement in specific cell types, but suggest links between chromatin compaction, chromatin movement and transcription. We describe the interplay between enhancer-promoter proximity within the context of biomolecular condensates and the need to understand how condensate microenvironments influence the chromatin binding kinetics of proteins that bind at cis-regulatory elements to activate transcription. Finally, given the complex multi-scale interplay between regulatory proteins, enhancer-promoter proximity and movement, we propose the need to integrate information from complementary single-cell next-generation sequencing and live-cell imaging approaches to derive unified 3D theoretical models of enhancer-promoter communication that are ultimately predictive of transcriptional output and cell fate. In time, improved models will shed light on how tissues grow and diseases emerge.

Long-range promoter-enhancer contacts are conserved during evolution and contribute to gene expression robustness

Gene expression is regulated through complex molecular interactions, involving cis-acting elements that can be situated far away from their target genes. Data on long-range contacts between promoters and regulatory elements are rapidly accumulating. However, it remains unclear how these regulatory relationships evolve and how they contribute to the establishment of robust gene expression profiles. Here, we address these questions by comparing genome-wide maps of promoter-centered chromatin contacts in mouse and human. We show that there is significant evolutionary conservation of cis-regulatory landscapes, indicating that selective pressures act to preserve not only regulatory element sequences but also their chromatin contacts with target genes. The extent of evolutionary conservation is remarkable for long-range promoter–enhancer contacts, illustrating how the structure of regulatory landscapes constrains large-scale genome evolution. We show that the evolution of cis-regulatory landscapes, measured in terms of distal element sequences, synteny, or contacts with target genes, is significantly associated with gene expression evolution.

Improving Chromatin-Interaction Prediction Using Single-Cell Open-Chromatin Profiles and Making Insight Into the Cis-Regulatory Landscape of the Human Brain

Single-cell open-chromatin profiles have the potential to reveal the pattern of chromatin-interaction in a cell type. However, currently available cis-regulatory network prediction methods using single-cell open-chromatin profiles focus more on local chromatin interactions despite the fact that long-range interactions among genomic sites play a significant role in gene regulation. Here, we propose a method that predicts both short and long-range interactions among genomic sites using single-cell open chromatin profiles. Our method, termed as single-cell epigenome based chromatin-interaction analysis (scEChIA) exploits signal imputation and refined L1 regularization. For a few single-cell open-chromatin profiles, scEChIA outperformed other tools even in terms of accuracy of prediction. Using scEChIA, we predicted almost 0.7 million interactions among genomic sites across seven cell types in the human brain. Further analysis revealed cell type for connection between genes and expression quantitative trait locus (eQTL) in the human brain and making insight about target genes of human-accelerated-elements and disease-associated mutations. Our analysis enabled by scEChIA also hints about the possible action of a few transcription factors (TFs), especially through long-range interaction in brain endothelial cells.

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.

Detecting chromosomal interactions in Capture Hi-C data with CHiCAGO and companion tools

Capture Hi-C is widely used to obtain high-resolution profiles of chromosomal interactions involving, at least on one end, regions of interest such as gene promoters. Signal detection in Capture Hi-C data is challenging and cannot be adequately accomplished with tools developed for other chromosome conformation capture methods, including standard Hi-C. Capture Hi-C Analysis of Genomic Organization (CHiCAGO) is a computational pipeline developed specifically for Capture Hi-C analysis. It implements a statistical model accounting for biological and technical background components, as well as bespoke normalization and multiple testing procedures for this data type. Here we provide a step-by-step guide to the CHiCAGO workflow that is aimed at users with basic experience of the command line and R. We also describe more advanced strategies for tuning the key parameters for custom experiments and provide guidance on data preprocessing and downstream analysis using companion tools. In a typical experiment, CHiCAGO takes ~2-3 h to run, although pre- and postprocessing steps may take much longer.

Chromatin lncRNA Platr10 controls stem cell pluripotency by coordinating an intrachromosomal regulatory network

Background

A specific 3-dimensional intrachromosomal architecture of core stem cell factor genes is required to reprogram a somatic cell into pluripotency. As little is known about the epigenetic readers that orchestrate this architectural remodeling, we used a novel chromatin RNA in situ reverse transcription sequencing (CRIST-seq) approach to profile long noncoding RNAs (lncRNAs) in the Oct4 promoter.

Results

We identify Platr10 as an Oct4 - Sox2 binding lncRNA that is activated in somatic cell reprogramming. Platr10 is essential for the maintenance of pluripotency, and lack of this lncRNA causes stem cells to exit from pluripotency. In fibroblasts, ectopically expressed Platr10 functions in trans to activate core stem cell factor genes and enhance pluripotent reprogramming. Using RNA reverse transcription-associated trap sequencing (RAT-seq), we show that Platr10 interacts with multiple pluripotency-associated genes, including Oct4, Sox2, Klf4, and c-Myc, which have been extensively used to reprogram somatic cells. Mechanistically, we demonstrate that Platr10 helps orchestrate intrachromosomal promoter-enhancer looping and recruits TET1, the enzyme that actively induces DNA demethylation for the initiation of pluripotency. We further show that Platr10 contains an Oct4 binding element that interacts with the Oct4 promoter and a TET1-binding element that recruits TET1. Mutation of either of these two elements abolishes Platr10 activity.

Conclusion

These data suggest that Platr10 functions as a novel chromatin RNA molecule to control pluripotency in trans by modulating chromatin architecture and regulating DNA methylation in the core stem cell factor network.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02444-6.

K-mer Content Changes with Node Degree in Promoter-Enhancer Network of Mouse ES Cells

Maps of Hi-C contacts between promoters and enhancers can be analyzed as networks, with cis-regulatory regions as nodes and their interactions as edges. We checked if in the published promoter-enhancer network of mouse embryonic stem (ES) cells the differences in the node type (promoter or enhancer) and the node degree (number of regions interacting with a given promoter or enhancer) are reflected by sequence composition or sequence similarity of the interacting nodes. We used counts of all k-mers (k = 4) to analyze the sequence composition and the Euclidean distance between the k-mer count vectors (k-mer distance) as the measure of sequence (dis)similarity. The results we obtained with 4-mers are interpretable in terms of dinucleotides. Promoters are GC-rich as compared to enhancers, which is known. Enhancers are enriched in scaffold/matrix attachment regions (S/MARs) patterns and depleted of CpGs. Furthermore, we show that promoters are more similar to their interacting enhancers than vice-versa. Most notably, in both promoters and enhancers, the GC content and the CpG count increase with the node degree. As a consequence, enhancers of higher node degree become more similar to promoters, whereas higher degree promoters become less similar to enhancers. We confirmed the key results also for human keratinocytes.

A Tale of Two States: Pluripotency Regulation of Telomeres

Inside the nucleus, chromatin is functionally organized and maintained as a complex three-dimensional network of structures with different accessibility such as compartments, lamina associated domains, and membraneless bodies. Chromatin is epigenetically and transcriptionally regulated by an intricate and dynamic interplay of molecular processes to ensure genome stability. Phase separation, a process that involves the spontaneous organization of a solution into separate phases, has been proposed as a mechanism for the timely coordination of several cellular processes, including replication, transcription and DNA repair. Telomeres, the repetitive structures at the end of chromosomes, are epigenetically maintained in a repressed heterochromatic state that prevents their recognition as double-strand breaks (DSB), avoiding DNA damage repair and ensuring cell proliferation. In pluripotent embryonic stem cells, telomeres adopt a non-canonical, relaxed epigenetic state, which is characterized by a low density of histone methylation and expression of telomere non-coding transcripts (TERRA). Intriguingly, this telomere non-canonical conformation is usually associated with chromosome instability and aneuploidy in somatic cells, raising the question of how genome stability is maintained in a pluripotent background. In this review, we will explore how emerging technological and conceptual developments in 3D genome architecture can provide novel mechanistic perspectives for the pluripotent epigenetic paradox at telomeres. In particular, as RNA drives the formation of LLPS, we will consider how pluripotency-associated high levels of TERRA could drive and coordinate phase separation of several nuclear processes to ensure genome stability. These conceptual advances will provide a better understanding of telomere regulation and genome stability within the highly dynamic pluripotent background.

Regulation of Chromatin Organization in Cell Stemness: The Emerging Role of Long Non-coding RNAs

Chromatin is organized as chromosome territories in the nucleus of an interphase cell. The cell-type- and cell-state-specific organization of chromatin including the location, volume, compaction level, and spatial arrangement of chromosome territories are the major determinants of genome function. In addition, in response to different signaling stimuli and regulatory cues, it is the dynamic adaptation of chromatin structure that establishes and organizes transcriptional programs. It is known that varying levels of stemness are defined by gene regulatory networks. Accordingly, chromatin is the main milieu to host the transcriptional programs and gene regulatory networks responsible for the stemness status of a cell. In this review, our current understanding of the spatial organization of chromatin and the ways by which it defines stemness are discussed. In particular, the role of lncRNAs that regulate and affect chromatin organization and stemness properties are delineated. These roles can be categorized into the topics of specific binding to and epigenetic regulation of the promoter of pluripotency genes, their interaction with transcription factors, coordinating the intra- and inter-chromosomal looping of pluripotency-related genes, and their RNA-independent functions. This review brings together the results of studies that have begun to clarify the emerging roles of lncRNAs in the regulation of chromatin organization adapted for stemness and cancer plasticity.

Auxin-degron system identifies immediate mechanisms of OCT4

The pluripotency factor OCT4 is essential for the maintenance of naive pluripotent stem cells in vitro and in vivo. However, the specific role of OCT4 in this process remains unknown. Here, we developed a rapid protein-level OCT4 depletion system that demonstrates that the immediate downstream response to loss of OCT4 is reduced expression of key pluripotency factors. Our data show a requirement for OCT4 for the efficient transcription of several key pluripotency factors and suggest that expression of trophectoderm markers is a subsequent event. In addition, we find that NANOG is able to bind to the genome in the absence of OCT4, and this binding is in fact enhanced. Globally, however, the active enhancer-associated histone mark H3K27ac is depleted. Our work establishes that, while OCT4 is required for the maintenance of the naive transcription factor network, at a normal embryonic stem cell levels it antagonizes this network through inhibition of NANOG binding.

3D chromatin architecture and epigenetic regulation in cancer stem cells

Dedifferentiation of cell identity to a progenitor-like or stem cell-like state with increased cellular plasticity is frequently observed in cancer formation. During this process, a subpopulation of cells in tumours acquires a stem cell-like state partially resembling to naturally occurring pluripotent stem cells that are temporarily present during early embryogenesis. Such characteristics allow these cancer stem cells (CSCs) to give rise to the whole tumour with its entire cellular heterogeneity and thereby support metastases formation while being resistant to current cancer therapeutics. Cancer development and progression are demarcated by transcriptional dysregulation. In this article, we explore the epigenetic mechanisms shaping gene expression during tumorigenesis and cancer stem cell formation, with an emphasis on 3D chromatin architecture. Comparing the pluripotent stem cell state and epigenetic reprogramming to dedifferentiation in cellular transformation provides intriguing insight to chromatin dynamics. We suggest that the 3D chromatin architecture could be used as a target for re-sensitizing cancer stem cells to therapeutics.

Transcriptional Regulation by (Super)Enhancers: From Discovery to Mechanisms

Accurate control of gene expression in the right cell at the right moment is of fundamental importance to animal development and homeostasis. At the heart of gene regulation lie the enhancers, a class of gene regulatory elements that ensures precise spatiotemporal activation of gene transcription. Mammalian genomes are littered with enhancers, which are frequently organized in cooperative clusters such as locus control regions and superenhancers. Here, we discuss our current knowledge of enhancer biology, including an overview of the discovery of the various enhancer subsets and the mechanistic models used to explain their gene regulatory function.

Widespread reorganisation of pluripotent factor binding and gene regulatory interactions between human pluripotent states

The transition from naive to primed pluripotency is accompanied by an extensive reorganisation of transcriptional and epigenetic programmes. However, the role of transcriptional enhancers and three-dimensional chromatin organisation in coordinating these developmental programmes remains incompletely understood. Here, we generate a high-resolution atlas of gene regulatory interactions, chromatin profiles and transcription factor occupancy in naive and primed human pluripotent stem cells, and develop a network-graph approach to examine the atlas at multiple spatial scales. We uncover highly connected promoter hubs that change substantially in interaction frequency and in transcriptional co-regulation between pluripotent states. Small hubs frequently merge to form larger networks in primed cells, often linked by newly-formed Polycomb-associated interactions. We identify widespread state-specific differences in enhancer activity and interactivity that correspond with an extensive reconfiguration of OCT4, SOX2 and NANOG binding and target gene expression. These findings provide multilayered insights into the chromatin-based gene regulatory control of human pluripotent states.

Base-resolution models of transcription-factor binding reveal soft motif syntax

The arrangement (syntax) of transcription factor (TF) binding motifs is an important part of the cis-regulatory code, yet remains elusive. We introduce a deep learning model, BPNet, that uses DNA sequence to predict base-resolution chromatin immunoprecipitation (ChIP)-nexus binding profiles of pluripotency TFs. We develop interpretation tools to learn predictive motif representations and identify soft syntax rules for cooperative TF binding interactions. Strikingly, Nanog preferentially binds with helical periodicity, and TFs often cooperate in a directional manner, which we validate using clustered regularly interspaced short palindromic repeat (CRISPR)-induced point mutations. Our model represents a powerful general approach to uncover the motifs and syntax of cis-regulatory sequences in genomics data.

Mechanisms of Enhancer-Promoter Interactions in Higher Eukaryotes

In higher eukaryotes, enhancers determine the activation of developmental gene transcription in specific cell types and stages of embryogenesis. Enhancers transform the signals produced by various transcription factors within a given cell, activating the transcription of the targeted genes. Often, developmental genes can be associated with dozens of enhancers, some of which are located at large distances from the promoters that they regulate. Currently, the mechanisms underlying specific distance interactions between enhancers and promoters remain poorly understood. This review briefly describes the properties of enhancers and discusses the mechanisms of distance interactions and potential proteins involved in this process.

Dynamic 3D Chromatin Reorganization during Establishment and Maintenance of Pluripotency

Higher-order chromatin structure is tightly linked to gene expression and therefore cell identity. In recent years, the chromatin landscape of pluripotent stem cells has become better characterized, and unique features at various architectural levels have been revealed. However, the mechanisms that govern establishment and maintenance of these topological characteristics and the temporal and functional relationships with transcriptional or epigenetic features are still areas of intense study. Here, we will discuss progress and limitations of our current understanding regarding how the 3D chromatin topology of pluripotent stem cells is established during somatic cell reprogramming and maintained during cell division. We will also discuss evidence and theories about the driving forces of topological reorganization and the functional links with key features and properties of pluripotent stem cell identity.

How to Use Online Tools to Generate New Hypotheses for Mammary Gland Biology Research: A Case Study for Wnt7b

An increasing number of '-omics' datasets, generated by labs all across the world, are becoming available. They contain a wealth of data that are largely unexplored. Not every scientist, however, will have access to the required resources and expertise to analyze such data from scratch. Fortunately, a growing number of investigators is dedicating their time and effort to the development of user friendly, online applications that allow researchers to use and investigate these datasets. Here, we will illustrate the usefulness of such an approach. Using regulation of Wnt7b expression as an example, we will highlight a selection of accessible tools and resources that are available to researchers in the area of mammary gland biology. We show how they can be used for in silico analyses of gene regulatory mechanisms, resulting in new hypotheses and providing leads for experimental follow up. We also call out to the mammary gland community to join forces in a coordinated effort to generate and share additional tissue-specific '-omics' datasets and thereby expand the in silico toolbox.

High-resolution three-dimensional chromatin profiling of the Chinese hamster ovary cell genome

Chinese hamster ovary (CHO) cell lines are the pillars of a multibillion‐dollar biopharmaceutical industry producing recombinant therapeutic proteins. The effects of local chromatin organization and epigenetic repression within these cell lines result in unpredictable and unstable transgene expression following random integration. Limited knowledge of the CHO genome and its higher order chromatin organization has thus far impeded functional genomics approaches required to tackle these issues. Here, we present an integrative three‐dimensional (3D) map of genome organization within the CHOK1SV® 10E9 cell line in conjunction with an improved, less fragmented CHOK1SV 10E9 genome assembly. Using our high‐resolution chromatin conformation datasets, we have assigned ≈90% of sequence to a chromosome‐scale genome assembly. Our genome‐wide 3D map identifies higher order chromatin structures such as topologically associated domains, incorporates our chromatin accessibility data to enhance the identification of active cis‐regulatory elements, and importantly links these cis‐regulatory elements to target promoters in a 3D promoter interactome. We demonstrate the power of our improved functional annotation by evaluating the 3D landscape of a transgene integration site and two phenotypically different cell lines. Our work opens up further novel genome engineering targets, has the potential to inform vital improvements for industrial biotherapeutic production, and represents a significant advancement for CHO cell line development.

Navigating the crowd: visualizing coordination between genome dynamics, structure, and transcription

The eukaryotic genome is hierarchically structured yet highly dynamic. Regulating transcription in this environment demands a high level of coordination to permit many proteins to interact with chromatin fiber at appropriate sites in a timely manner. We describe how recent advances in quantitative imaging techniques overcome caveats of sequencing-based methods (Hi-C and related) by enabling direct visualization of transcription factors and chromatin at high resolution, from single genes to the whole nucleus. We discuss the contribution of fluorescence imaging to deciphering the principles underlying this coordination within the crowded nuclear space in living cells and discuss challenges ahead.

LncRNA Platr22 promotes super-enhancer activity and stem cell pluripotency

Super-enhancers (SEs) comprise large clusters of enhancers, which are co-occupied by multiple lineage-specific and master transcription factors, and play pivotal roles in regulating gene expression and cell fate determination. However, it is still largely unknown whether and how SEs are regulated by the noncoding portion of the genome. Here, through genome-wide analysis, we found that long noncoding RNA (lncRNA) genes preferentially lie next to SEs. In mouse embryonic stem cells (mESCs), depletion of SE-associated lncRNA transcripts dysregulated the activity of their nearby SEs. Specifically, we revealed a critical regulatory role of the lncRNA gene Platr22 in modulating the activity of a nearby SE and the expression of the nearby pluripotency regulator ZFP281. Through these regulatory events, Platr22 contributes to pluripotency maintenance and proper differentiation of mESCs. Mechanistically, Platr22 transcripts coat chromatin near the SE region and interact with DDX5 and hnRNP-L. DDX5 further recruits p300 and other factors related to active transcription. We propose that these factors assemble into a transcription hub, thus promoting an open and active epigenetic chromatin state. Our study highlights an unanticipated role for a class of lncRNAs in epigenetically controlling the activity and vulnerability to perturbation of nearby SEs for cell fate determination.

JMJD3 acts in tandem with KLF4 to facilitate reprogramming to pluripotency

The interplay between the Yamanaka factors (OCT4, SOX2, KLF4 and c-MYC) and transcriptional/epigenetic co-regulators in somatic cell reprogramming is incompletely understood. Here, we demonstrate that the histone H3 lysine 27 trimethylation (H3K27me3) demethylase JMJD3 plays conflicting roles in mouse reprogramming. On one side, JMJD3 induces the pro-senescence factor Ink4a and degrades the pluripotency regulator PHF20 in a reprogramming factor-independent manner. On the other side, JMJD3 is specifically recruited by KLF4 to reduce H3K27me3 at both enhancers and promoters of epithelial and pluripotency genes. JMJD3 also promotes enhancer-promoter looping through the cohesin loading factor NIPBL and ultimately transcriptional elongation. This competition of forces can be shifted towards improved reprogramming by using early passage fibroblasts or boosting JMJD3’s catalytic activity with vitamin C. Our work, thus, establishes a multifaceted role for JMJD3, placing it as a key partner of KLF4 and a scaffold that assists chromatin interactions and activates gene transcription.

Previous work suggested that histone demethylase JMJD3 is detrimental to somatic cell reprogramming. Here, the authors show that while JMJD3 has a context-independent detrimental effect on early stages of reprogramming, during late stages it activates epithelial and pluripotency genes together with Klf4.

Global hyperactivation of enhancers stabilizes human and mouse naive pluripotency through inhibition of CDK8/19 Mediator kinases

Pluripotent stem cells (PSCs) transition between cell states in vitro, reflecting developmental changes in the early embryo. PSCs can be stabilized in the naive state by blocking extracellular differentiation stimuli, particularly FGF-MEK signalling. Here, we report that multiple features of the naive state in human and mouse PSCs can be recapitulated without affecting FGF-MEK signalling or global DNA methylation. Mechanistically, chemical inhibition of CDK8 and CDK19 (hereafter CDK8/19) kinases removes their ability to repress the Mediator complex at enhancers. CDK8/19 inhibition therefore increases Mediator-driven recruitment of RNA polymerase II (RNA Pol II) to promoters and enhancers. This efficiently stabilizes the naive transcriptional program and confers resistance to enhancer perturbation by BRD4 inhibition. Moreover, naive pluripotency during embryonic development coincides with a reduction in CDK8/19. We conclude that global hyperactivation of enhancers drives naive pluripotency, and this can be achieved in vitro by inhibiting CDK8/19 kinase activity. These principles may apply to other contexts of cellular plasticity.

Manipulating the Mediator complex to induce naïve pluripotency

Human naïve pluripotent stem cells (PSCs) represent an optimal homogenous starting point for molecular interventions and differentiation strategies. This is in contrast to the standard primed PSCs which fluctuate in identity and are transcriptionally heterogeneous. However, despite many efforts, the maintenance and expansion of human naïve PSCs remains a challenge. Here, we discuss our recent strategy for the stabilization of human PSC in the naïve state based on the use of a single chemical inhibitor of the related kinases CDK8 and CDK19. These kinases phosphorylate and negatively regulate the multiprotein Mediator complex, which is critical for enhancer-driven recruitment of RNA Pol II. The net effect of CDK8/19 inhibition is a global stimulation of enhancers, which in turn reinforces transcriptional programs including those related to cellular identity. In the case of pluripotent cells, the presence of CDK8/19i efficiently stabilizes the naïve state. Importantly, in contrast to previous chemical methods to induced the naïve state based on the inhibition of the FGF-MEK-ERK pathway, CDK8/19i-naïve human PSCs are chromosomally stable and retain developmental potential after long-term expansion. We suggest this could be related to the fact that CDK8/19 inhibition does not induce DNA demethylation. These principles may apply to other fate decisions.

Molecular Basis of the Function of Transcriptional Enhancers

Transcriptional enhancers are major genomic elements that control gene activity in eukaryotes. Recent studies provided deeper insight into the temporal and spatial organization of transcription in the nucleus, the role of non-coding RNAs in the process, and the epigenetic control of gene expression. Thus, multiple molecular details of enhancer functioning were revealed. Here, we describe the recent data and models of molecular organization of enhancer-driven transcription.

Interplay between Metabolites and the Epigenome in Regulating Embryonic and Adult Stem Cell Potency and Maintenance

The environment surrounding stem cells has the ability to elicit profound, heritable epigenetic changes orchestrated by multiple epigenetic mechanisms, which can be modulated by the level of specific metabolites. In this review, we highlight the significance of metabolism in regulating stem cell homeostasis, cell state, and differentiation capacity, using metabolic regulation of embryonic and adult muscle stem cells as examples, and cast light on the interaction between cellular metabolism and epigenetics. These new regulatory networks, based on the dynamic interplay between metabolism and epigenetics in stem cell biology, are important, not only for understanding tissue homeostasis, but to determine in vitro culture conditions which accurately support normal cell physiology.

Skeletal muscle enhancer interactions identify genes controlling whole-body metabolism

Obesity and type 2 diabetes (T2D) are metabolic disorders influenced by lifestyle and genetic factors that are characterized by insulin resistance in skeletal muscle, a prominent site of glucose disposal. Numerous genetic variants have been associated with obesity and T2D, of which the majority are located in non-coding DNA regions. This suggests that most variants mediate their effect by altering the activity of gene-regulatory elements, including enhancers. Here, we map skeletal muscle genomic enhancer elements that are dynamically regulated after exposure to the free fatty acid palmitate or the inflammatory cytokine TNFα. By overlapping enhancer positions with the location of disease-associated genetic variants, and resolving long-range chromatin interactions between enhancers and gene promoters, we identify target genes involved in metabolic dysfunction in skeletal muscle. The majority of these genes also associate with altered whole-body metabolic phenotypes in the murine BXD genetic reference population. Thus, our combined genomic investigations identified genes that are involved in skeletal muscle metabolism.

Obesity and type 2 diabetes (T2D) are metabolic disorders characterized by insulin resistance in skeletal muscle. Here, the authors map skeletal muscle enhancer elements dynamically regulated after exposure to free fatty acid palmitate or inflammatory cytokine TNFα and identify target genes involved in metabolic dysfunction in skeletal muscle.

Perspectives on the use of super-enhancers as a defining feature of cell/tissue-identity genes

Super-enhancers (SE) have become a popular concept and are widely used as a feature defining key identity genes. Here, we provide perspectives on the use of SE to define and identify cell/tissue-identity genes. By mining SE and their associated genes using murine functional genomics data, we highlight and discuss current limitations and open questions regarding both the sensitivity and specificity of identity genes/transcription factors predicted by SE. In this context, we point to cell/tissue-specific promoters as an important additional level of information, which we propose to combine with SE when aiming to define potential identity genes.

Transcriptional network dynamics during the progression of pluripotency revealed by integrative statistical learning

The developmental potential of cells, termed pluripotency, is highly dynamic and progresses through a continuum of naive, formative and primed states. Pluripotency progression of mouse embryonic stem cells (ESCs) from naive to formative and primed state is governed by transcription factors (TFs) and their target genes. Genomic techniques have uncovered a multitude of TF binding sites in ESCs, yet a major challenge lies in identifying target genes from functional binding sites and reconstructing dynamic transcriptional networks underlying pluripotency progression. Here, we integrated time-resolved ‘trans-omic’ datasets together with TF binding profiles and chromatin conformation data to identify target genes of a panel of TFs. Our analyses revealed that naive TF target genes are more likely to be TFs themselves than those of formative TFs, suggesting denser hierarchies among naive TFs. We also discovered that formative TF target genes are marked by permissive epigenomic signatures in the naive state, indicating that they are poised for expression prior to the initiation of pluripotency transition to the formative state. Finally, our reconstructed transcriptional networks pinpointed the precise timing from naive to formative pluripotency progression and enabled the spatiotemporal mapping of differentiating ESCs to their in vivo counterparts in developing embryos.

TADCompare: An R Package for Differential and Temporal Analysis of Topologically Associated Domains

Recent research using chromatin conformation capture technologies, such as Hi-C, has demonstrated the importance of topologically associated domains (TADs) and smaller chromatin loops, collectively referred hereafter as "interacting domains." Many such domains change during development or disease, and exhibit cell- and condition-specific differences. Quantification of the dynamic behavior of interacting domains will help to better understand genome regulation. Methods for comparing interacting domains between cells and conditions are highly limited. We developed TADCompare, a method for differential analysis of boundaries of interacting domains between two or more Hi-C datasets. TADCompare is based on a spectral clustering-derived measure called the eigenvector gap, which enables a loci-by-loci comparison of boundary differences. Using this measure, we introduce methods for identifying differential and consensus boundaries of interacting domains and tracking boundary changes over time. We further propose a novel framework for the systematic classification of boundary changes. Colocalization- and gene enrichment analysis of different types of boundary changes demonstrated distinct biological functionality associated with them. TADCompare is available on https://github.com/dozmorovlab/TADCompare and Bioconductor (submitted).

Dynamic CpG methylation delineates subregions within super-enhancers selectively decommissioned at the exit from naive pluripotency

Clusters of enhancers, referred as to super-enhancers (SEs), control the expression of cell identity genes. The organisation of these clusters, and how they are remodelled upon developmental transitions remain poorly understood. Here, we report the existence of two types of enhancer units within SEs typified by distinctive CpG methylation dynamics in embryonic stem cells (ESCs). We find that these units are either prone for decommissioning or remain constitutively active in epiblast stem cells (EpiSCs), as further established in the peri-implantation epiblast in vivo. Mechanistically, we show a pivotal role for ESRRB in regulating the activity of ESC-specific enhancer units and propose that the developmentally regulated silencing of ESRRB triggers the selective inactivation of these units within SEs. Our study provides insights into the molecular events that follow the loss of ESRRB binding, and offers a mechanism by which the naive pluripotency transcriptional programme can be partially reset upon embryo implantation.

Cohesin controls intestinal stem cell identity by maintaining association of Escargot with target promoters

Intestinal stem cells (ISCs) maintain regenerative capacity of the intestinal epithelium. Their function and activity are regulated by transcriptional changes, yet how such changes are coordinated at the genomic level remains unclear. The Cohesin complex regulates transcription globally by generating topologically-associated DNA domains (TADs) that link promotor regions with distant enhancers. We show here that the Cohesin complex prevents premature differentiation of Drosophila ISCs into enterocytes (ECs). Depletion of the Cohesin subunit Rad21 and the loading factor Nipped-B triggers an ISC to EC differentiation program that is independent of Notch signaling, but can be rescued by over-expression of the ISC-specific escargot (esg) transcription factor. Using damID and transcriptomic analysis, we find that Cohesin regulates Esg binding to promoters of differentiation genes, including a group of Notch target genes involved in ISC differentiation. We propose that Cohesin ensures efficient Esg-dependent gene repression to maintain stemness and intestinal homeostasis.

Chromatin regulation and dynamics in stem cells

In eukaryotes, DNA is highly compacted within the nucleus into a structure known as chromatin. Modulation of chromatin structure allows for precise regulation of gene expression, and thereby controls cell fate decisions. Specific chromatin organization is established and preserved by numerous factors to generate desired cellular outcomes. In embryonic stem (ES) cells, chromatin is precisely regulated to preserve their two defining characteristics: self-renewal and pluripotent state. This action is accomplished by a litany of nucleosome remodelers, histone variants, epigenetic marks, and other chromatin regulatory factors. These highly dynamic regulatory factors come together to precisely define a chromatin state that is conducive to ES cell maintenance and development, where dysregulation threatens the survival and fitness of the developing organism.

CRUP: a comprehensive framework to predict condition-specific regulatory units

We present the software Condition-specific Regulatory Units Prediction (CRUP) to infer from epigenetic marks a list of regulatory units consisting of dynamically changing enhancers with their target genes. The workflow consists of a novel pre-trained enhancer predictor that can be reliably applied across cell types and species, solely based on histone modification ChIP-seq data. Enhancers are subsequently assigned to different conditions and correlated with gene expression to derive regulatory units. We thoroughly test and then apply CRUP to a rheumatoid arthritis model, identifying enhancer-gene pairs comprising known disease genes as well as new candidate genes.

KLF4 is involved in the organization and regulation of pluripotency-associated three-dimensional enhancer networks

Cell fate transitions are accompanied by global transcriptional, epigenetic and topological changes driven by transcription factors, as is exemplified by reprogramming somatic cells to pluripotent stem cells through the expression of OCT4, KLF4, SOX2 and cMYC. How transcription factors orchestrate the complex molecular changes around their target gene loci remains incompletely understood. Here, using KLF4 as a paradigm, we provide a transcription-factor-centric view of chromatin reorganization and its association with three-dimensional enhancer rewiring and transcriptional changes during the reprogramming of mouse embryonic fibroblasts to pluripotent stem cells. Inducible depletion of KLF factors in PSCs caused a genome-wide decrease in enhancer connectivity, whereas disruption of individual KLF4 binding sites within pluripotent-stem-cell-specific enhancers was sufficient to impair enhancer-promoter contacts and reduce the expression of associated genes. Our study provides an integrative view of the complex activities of a lineage-specifying transcription factor and offers novel insights into the nature of the molecular events that follow transcription factor binding.

Visualising three-dimensional genome organisation in two dimensions

The three-dimensional organisation of the genome plays a crucial role in developmental gene regulation. In recent years, techniques to investigate this organisation have become more accessible to labs worldwide due to improvements in protocols and decreases in the cost of high-throughput sequencing. However, the resulting datasets are complex and can be challenging to analyse and interpret. Here, we provide a guide to visualisation approaches that can aid the interpretation of such datasets and the communication of biological results.