OCT2 pre-positioning facilitates cell fate transition and chromatin architecture changes in humoral immunity.

During the germinal center (GC) reaction, B cells undergo profound transcriptional, epigenetic and genomic architectural changes. How such changes are established remains unknown. Mapping chromatin accessibility during the humoral immune response, we show that OCT2 was the dominant transcription factor linked to differential accessibility of GC regulatory elements. Silent chromatin regions destined to become GC-specific super-enhancers (SEs) contained pre-positioned OCT2-binding sites in naive B cells (NBs). These preloaded SE 'seeds' featured spatial clustering of regulatory elements enriched in OCT2 DNA-binding motifs that became heavily loaded with OCT2 and its GC-specific coactivator OCAB in GC B cells (GCBs). SEs with high abundance of pre-positioned OCT2 binding preferentially formed long-range chromatin contacts in GCs, to support expression of GC-specifying factors. Gain in accessibility and architectural interactivity of these regions were dependent on recruitment of OCAB. Pre-positioning key regulators at SEs may represent a broadly used strategy for facilitating rapid cell fate transitions.

H3K27ac bookmarking promotes rapid post-mitotic activation of the pluripotent stem cell program without impacting 3D chromatin reorganization.

During self-renewal, cell-type-defining features are drastically perturbed in mitosis and must be faithfully reestablished upon G1 entry, a process that remains largely elusive. Here, we characterized at a genome-wide scale the dynamic transcriptional and architectural resetting of mouse pluripotent stem cells (PSCs) upon mitotic exit. We captured distinct waves of transcriptional reactivation with rapid induction of stem cell genes and transient activation of lineage-specific genes. Topological reorganization at different hierarchical levels also occurred in an asynchronous manner and showed partial coordination with transcriptional resetting. Globally, rapid transcriptional and architectural resetting associated with mitotic retention of H3K27 acetylation, supporting a bookmarking function. Indeed, mitotic depletion of H3K27ac impaired the early reactivation of bookmarked, stem-cell-associated genes. However, 3D chromatin reorganization remained largely unaffected, suggesting that these processes are driven by distinct forces upon mitotic exit. This study uncovers principles and mediators of PSC molecular resetting during self-renewal.

A molecular roadmap of reprogramming somatic cells into iPS cells.

Factor-induced reprogramming of somatic cells into induced pluripotent stem cells (iPSCs) is inefficient, complicating mechanistic studies. Here, we examined defined intermediate cell populations poised to becoming iPSCs by genome-wide analyses. We show that induced pluripotency elicits two transcriptional waves, which are driven by c-Myc/Klf4 (first wave) and Oct4/Sox2/Klf4 (second wave). Cells that become refractory to reprogramming activate the first but fail to initiate the second transcriptional wave and can be rescued by elevated expression of all four factors. The establishment of bivalent domains occurs gradually after the first wave, whereas changes in DNA methylation take place after the second wave when cells acquire stable pluripotency. This integrative analysis allowed us to identify genes that act as roadblocks during reprogramming and surface markers that further enrich for cells prone to forming iPSCs. Collectively, our data offer new mechanistic insights into the nature and sequence of molecular events inherent to cellular reprogramming.

Histone H1 loss drives lymphoma by disrupting 3D chromatin architecture.

Linker histone H1 proteins bind to nucleosomes and facilitate chromatin compaction1, although their biological functions are poorly understood. Mutations in the genes that encode H1 isoforms B-E (H1B, H1C, H1D and H1E; also known as H1-5, H1-2, H1-3 and H1-4, respectively) are highly recurrent in B cell lymphomas, but the pathogenic relevance of these mutations to cancer and the mechanisms that are involved are unknown. Here we show that lymphoma-associated H1 alleles are genetic driver mutations in lymphomas. Disruption of H1 function results in a profound architectural remodelling of the genome, which is characterized by large-scale yet focal shifts of chromatin from a compacted to a relaxed state. This decompaction drives distinct changes in epigenetic states, primarily owing to a gain of histone H3 dimethylation at lysine 36 (H3K36me2) and/or loss of repressive H3 trimethylation at lysine 27 (H3K27me3). These changes unlock the expression of stem cell genes that are normally silenced during early development. In mice, loss of H1c and H1e (also known as H1f2 and H1f4, respectively) conferred germinal centre B cells with enhanced fitness and self-renewal properties, ultimately leading to aggressive lymphomas with an increased repopulating potential. Collectively, our data indicate that H1 proteins are normally required to sequester early developmental genes into architecturally inaccessible genomic compartments. We also establish H1 as a bona fide tumour suppressor and show that mutations in H1 drive malignant transformation primarily through three-dimensional genome reorganization, which leads to epigenetic reprogramming and derepression of developmentally silenced genes.

TAF5L and TAF6L Maintain Self-Renewal of Embryonic Stem Cells via the MYC Regulatory Network.

Self-renewal and pluripotency of the embryonic stem cell (ESC) state are established and maintained by multiple regulatory networks that comprise transcription factors and epigenetic regulators. While much has been learned regarding transcription factors, the function of epigenetic regulators in these networks is less well defined. We conducted a CRISPR-Cas9-mediated loss-of-function genetic screen that identified two epigenetic regulators, TAF5L and TAF6L, components or co-activators of the GNAT-HAT complexes for the mouse ESC (mESC) state. Detailed molecular studies demonstrate that TAF5L/TAF6L transcriptionally activate c-Myc and Oct4 and their corresponding MYC and CORE regulatory networks. Besides, TAF5L/TAF6L predominantly regulate their target genes through H3K9ac deposition and c-MYC recruitment that eventually activate the MYC regulatory network for self-renewal of mESCs. Thus, our findings uncover a role of TAF5L/TAF6L in directing the MYC regulatory network that orchestrates gene expression programs to control self-renewal for the maintenance of mESC state.

Virus Infection Induces NF-kappaB-dependent interchromosomal associations mediating monoallelic IFN-beta gene expression.

Transcriptional activation of the IFN-beta gene by virus infection requires the cooperative assembly of an enhanceosome. We report that the stochastic and monoallelic expression of the IFN-beta gene depends on interchromosomal associations with three identified distinct genetic loci that could mediate binding of the limiting transcription factor NF-kappaB to the IFN-beta enhancer, thus triggering enhanceosome assembly and activation of transcription from this allele. The probability of a cell to express IFN-beta is dramatically increased when the cell is transfected with any of these loci. The secreted IFN-beta protein induces high-level expression of the enhanceosome factor IRF-7, which in turn promotes enhanceosome assembly and IFN-beta transcription from the remaining alleles and in other initially nonexpressing cells. Thus, the IFN-beta enhancer functions in a nonlinear fashion by working as a signal amplifier.

Chromatin dynamics during cellular reprogramming.

Induced pluripotency is a powerful tool to derive patient-specific stem cells. In addition, it provides a unique assay to study the interplay between transcription factors and chromatin structure. Here, we review the latest insights into chromatin dynamics that are inherent to induced pluripotency. Moreover, we compare and contrast these events with other physiological and pathological processes that involve changes in chromatin and cell state, including germ cell maturation and tumorigenesis. We propose that an integrated view of these seemingly diverse processes could provide mechanistic insights into cell fate transitions in general and might lead to new approaches in regenerative medicine and cancer treatment.

Small molecules facilitate rapid and synchronous iPSC generation.

The reprogramming of somatic cells into induced pluripotent stem cells (iPSCs) upon overexpression of OCT4, KLF4, SOX2 and c-MYC (OKSM) provides a powerful system to interrogate basic mechanisms of cell fate change. However, iPSC formation with standard methods is typically protracted and inefficient, resulting in heterogeneous cell populations. We show that exposure of OKSM-expressing cells to both ascorbic acid and a GSK3-ß inhibitor (AGi) facilitates more synchronous and rapid iPSC formation from several mouse cell types. AGi treatment restored the ability of refractory cell populations to yield iPSC colonies, and it attenuated the activation of developmental regulators commonly observed during the reprogramming process. Moreover, AGi supplementation gave rise to chimera-competent iPSCs after as little as 48 h of OKSM expression. Our results offer a simple modification to the reprogramming protocol, facilitating iPSC induction at unparalleled efficiencies and enabling dissection of the underlying mechanisms in more homogeneous cell populations.

Aberrant silencing of imprinted genes on chromosome 12qF1 in mouse induced pluripotent stem cells.

Induced pluripotent stem cells (iPSCs) have been generated by enforced expression of defined sets of transcription factors in somatic cells. It remains controversial whether iPSCs are molecularly and functionally equivalent to blastocyst-derived embryonic stem (ES) cells. By comparing genetically identical mouse ES cells and iPSCs, we show here that their overall messenger RNA and microRNA expression patterns are indistinguishable with the exception of a few transcripts encoded within the imprinted Dlk1-Dio3 gene cluster on chromosome 12qF1, which were aberrantly silenced in most of the iPSC clones. Consistent with a developmental role of the Dlk1-Dio3 gene cluster, these iPSC clones contributed poorly to chimaeras and failed to support the development of entirely iPSC-derived animals ('all-iPSC mice'). In contrast, iPSC clones with normal expression of the Dlk1-Dio3 cluster contributed to high-grade chimaeras and generated viable all-iPSC mice. Notably, treatment of an iPSC clone that had silenced Dlk1-Dio3 with a histone deacetylase inhibitor reactivated the locus and rescued its ability to support full-term development of all-iPSC mice. Thus, the expression state of a single imprinted gene cluster seems to distinguish most murine iPSCs from ES cells and allows for the prospective identification of iPSC clones that have the full development potential of ES cells.

3D Enhancer-promoter networks provide predictive features for gene expression and coregulation in early embryonic lineages.

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages: the trophectoderm, the epiblast and the primitive endoderm. Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements through which transcriptional regulators enact these fates remain understudied. Here, we characterize, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observe extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although distinct groups of genes are irresponsive to topological changes. In each lineage, a high degree of connectivity, or 'hubness', positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a predictive model for transcriptional regulation (3D-HiChAT) that outperforms models using only 1D promoter or proximal variables to predict levels and cell-type specificity of gene expression. Using 3D-HiChAT, we identify, in silico, candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments, we validate several enhancers that control gene expression in their respective lineages. Our study identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to comprehensively understand lineage-specific transcriptional behaviors.

CRISPR Screening Uncovers a Long-Range Enhancer for ONECUT1 in Pancreatic Differentiation and Links a Diabetes Risk Variant.

Functional enhancer annotation is a valuable first step for understanding tissue-specific transcriptional regulation and prioritizing disease-associated non-coding variants for investigation. However, unbiased enhancer discovery in physiologically relevant contexts remains a major challenge. To discover regulatory elements pertinent to diabetes, we conducted a CRISPR interference screen in the human pluripotent stem cell (hPSC) pancreatic differentiation system. Among the enhancers uncovered, we focused on a long-range enhancer ∼664 kb from the ONECUT1 promoter, since coding mutations in ONECUT1 cause pancreatic hypoplasia and neonatal diabetes. Homozygous enhancer deletion in hPSCs was associated with a near-complete loss of ONECUT1 gene expression and compromised pancreatic differentiation. This enhancer contains a confidently fine-mapped type 2 diabetes associated variant (rs528350911) which disrupts a GATA motif. Introduction of the risk variant into hPSCs revealed substantially reduced binding of key pancreatic transcription factors (GATA4, GATA6 and FOXA2) on the edited allele, accompanied by a slight reduction of ONECUT1 transcription, supporting a causal role for this risk variant in metabolic disease. This work expands our knowledge about transcriptional regulation in pancreatic development through the characterization of a long-range enhancer and highlights the utility of enhancer discovery in disease-relevant settings for understanding monogenic and complex disease.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer.

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

3D enhancer-promoter interactions and multi-connected hubs: Organizational principles and functional roles.

The spatiotemporal control of gene expression is dependent on the activity of cis-acting regulatory sequences, called enhancers, which regulate target genes over variable genomic distances and, often, by skipping intermediate promoters, suggesting mechanisms that control enhancer-promoter communication. Recent genomics and imaging technologies have revealed highly complex enhancer-promoter interaction networks, whereas advanced functional studies have started interrogating the forces behind the physical and functional communication among multiple enhancers and promoters. In this review, we first summarize our current understanding of the factors involved in enhancer-promoter communication, with a particular focus on recent papers that have revealed new layers of complexities to old questions. In the second part of the review, we focus on a subset of highly connected enhancer-promoter "hubs" and discuss their potential functions in signal integration and gene regulation, as well as the putative factors that might determine their dynamics and assembly.

Single-cell analysis of bidirectional reprogramming between early embryonic states reveals mechanisms of differential lineage plasticities.

Two distinct fates, pluripotent epiblast (EPI) and primitive (extra-embryonic) endoderm (PrE), arise from common progenitor cells, the inner cell mass (ICM), in mammalian embryos. To study how these sister identities are forged, we leveraged embryonic (ES) and eXtraembryonic ENdoderm (XEN) stem cells - in vitro counterparts of the EPI and PrE. Bidirectional reprogramming between ES and XEN coupled with single-cell RNA and ATAC-seq analyses uncovered distinct rates, efficiencies and trajectories of state conversions, identifying drivers and roadblocks of reciprocal conversions. While GATA4-mediated ES-to-iXEN conversion was rapid and nearly deterministic, OCT4, KLF4 and SOX2-induced XEN-to-iPS reprogramming progressed with diminished efficiency and kinetics. The dominant PrE transcriptional program, safeguarded by Gata4, and globally elevated chromatin accessibility of EPI underscored the differential plasticities of the two states. Mapping in vitro trajectories to embryos revealed reprogramming in either direction tracked along, and toggled between, EPI and PrE in vivo states without transitioning through the ICM.

Systematic mapping and modeling of 3D enhancer-promoter interactions in early mouse embryonic lineages reveal regulatory principles that determine the levels and cell-type specificity of gene expression.

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages, the trophectoderm (TE), the epiblast (EPI) and the primitive endoderm (PrE). Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements via which transcriptional regulators enact these fates remain understudied. To address this gap, we have characterized, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observed extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although there are distinct groups of genes that are irresponsive to topological changes. In each lineage, a high degree of connectivity or "hubness" positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages, compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a novel predictive model for transcriptional regulation (3D-HiChAT), which outperformed models that use only 1D promoter or proximal variables in predicting levels and cell-type specificity of gene expression. Using 3D-HiChAT, we performed genome-wide in silico perturbations to nominate candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments we validated several novel enhancers that control expression of one or more genes in their respective lineages. Our study comprehensively identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to understand lineage-specific transcriptional behaviors.

Discovery of Competent Chromatin Regions in Human Embryonic Stem Cells.

The mechanisms underlying the ability of embryonic stem cells (ESCs) to rapidly activate lineage-specific genes during differentiation remain largely unknown. Through multiple CRISPR-activation screens, we discovered human ESCs have pre-established transcriptionally competent chromatin regions (CCRs) that support lineage-specific gene expression at levels comparable to differentiated cells. CCRs reside in the same topological domains as their target genes. They lack typical enhancer-associated histone modifications but show enriched occupancy of pluripotent transcription factors, DNA demethylation factors, and histone deacetylases. TET1 and QSER1 protect CCRs from excessive DNA methylation, while HDAC1 family members prevent premature activation. This "push and pull" feature resembles bivalent domains at developmental gene promoters but involves distinct molecular mechanisms. Our study provides new insights into pluripotency regulation and cellular plasticity in development and disease. One sentence summary: We report a class of distal regulatory regions distinct from enhancers that confer human embryonic stem cells with the competence to rapidly activate the expression of lineage-specific genes.

Dynamic network-guided CRISPRi screen identifies CTCF-loop-constrained nonlinear enhancer gene regulatory activity during cell state transitions.

Comprehensive enhancer discovery is challenging because most enhancers, especially those contributing to complex diseases, have weak effects on gene expression. Our gene regulatory network modeling identified that nonlinear enhancer gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Using human embryonic stem cell definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. We discovered a comprehensive set of enhancers for each of the core endoderm-specifying transcription factors. Many enhancers had strong effects mid-transition but weak effects post-transition, consistent with the nonlinear temporal responses to enhancer perturbation predicted by the modeling. Integrating three-dimensional genomic information, we were able to develop a CTCF-loop-constrained Interaction Activity model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and systematic enhancer discovery in both normal and pathological cell state transitions.

Dynamic network-guided CRISPRi screen reveals CTCF loop-constrained nonlinear enhancer-gene regulatory activity in cell state transitions.

Comprehensive enhancer discovery is challenging because most enhancers, especially those affected in complex diseases, have weak effects on gene expression. Our network modeling revealed that nonlinear enhancer-gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Utilizing hESC definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. The screen discovered a comprehensive set of enhancers (4 to 9 per locus) for each of the core endoderm lineage-specifying transcription factors, and many enhancers had strong effects mid-transition but weak effects post-transition. Through integrating enhancer activity measurements and three-dimensional enhancer-promoter interaction information, we were able to develop a CTCF loop-constrained Interaction Activity (CIA) model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and more comprehensive enhancer discovery in both normal and pathological cell state transitions.

ChromaFold predicts the 3D contact map from single-cell chromatin accessibility.

The identification of cell-type-specific 3D chromatin interactions between regulatory elements can help to decipher gene regulation and to interpret the function of disease-associated non-coding variants. However, current chromosome conformation capture (3C) technologies are unable to resolve interactions at this resolution when only small numbers of cells are available as input. We therefore present ChromaFold, a deep learning model that predicts 3D contact maps and regulatory interactions from single-cell ATAC sequencing (scATAC-seq) data alone. ChromaFold uses pseudobulk chromatin accessibility, co-accessibility profiles across metacells, and predicted CTCF motif tracks as input features and employs a lightweight architecture to enable training on standard GPUs. Once trained on paired scATAC-seq and Hi-C data in human cell lines and tissues, ChromaFold can accurately predict both the 3D contact map and peak-level interactions across diverse human and mouse test cell types. In benchmarking against a recent deep learning method that uses bulk ATAC-seq, DNA sequence, and CTCF ChIP-seq to make cell-type-specific predictions, ChromaFold yields superior prediction performance when including CTCF ChIP-seq data as an input and comparable performance without. Finally, fine-tuning ChromaFold on paired scATAC-seq and Hi-C in a complex tissue enables deconvolution of chromatin interactions across cell subpopulations. ChromaFold thus achieves state-of-the-art prediction of 3D contact maps and regulatory interactions using scATAC-seq alone as input data, enabling accurate inference of cell-type-specific interactions in settings where 3C-based assays are infeasible.

Genomic analysis reveals a novel nuclear factor-κB (NF-κB)-binding site in Alu-repetitive elements.

The transcription factor NF-κB is a critical regulator of immune responses. To determine how NF-κB builds transcriptional control networks, we need to obtain a topographic map of the factor bound to the genome and correlate it with global gene expression. We used a ChIP cloning technique and identified novel NF-κB target genes in response to virus infection. We discovered that most of the NF-κB-bound genomic sites deviate from the consensus and are located away from conventional promoter regions. Remarkably, we identified a novel abundant NF-κB-binding site residing in specialized Alu-repetitive elements having the potential for long range transcription regulation, thus suggesting that in addition to its known role, NF-κB has a primate-specific function and a role in human evolution. By combining these data with global gene expression profiling of virus-infected cells, we found that most of the sites bound by NF-κB in the human genome do not correlate with changes in gene expression of the nearby genes and they do not appear to function in the context of synthetic promoters. These results demonstrate that repetitive elements interspersed in the human genome function as common target sites for transcription factors and may play an important role in expanding the repertoire of binding sites to engage new genes into regulatory networks.