Capturing the Onset of PRC2-Mediated Repressive Domain Formation.

Polycomb repressive complex 2 (PRC2) maintains gene silencing by catalyzing methylation of histone H3 at lysine 27 (H3K27me2/3) within chromatin. By designing a system whereby PRC2-mediated repressive domains were collapsed and then reconstructed in an inducible fashion in vivo, a two-step mechanism of H3K27me2/3 domain formation became evident. First, PRC2 is stably recruited by the actions of JARID2 and MTF2 to a limited number of spatially interacting "nucleation sites," creating H3K27me3-forming Polycomb foci within the nucleus. Second, PRC2 is allosterically activated via its binding to H3K27me3 and rapidly spreads H3K27me2/3 both in cis and in far-cis via long-range contacts. As PRC2 proceeds further from the nucleation sites, its stability on chromatin decreases such that domains of H3K27me3 remain proximal, and those of H3K27me2 distal, to the nucleation sites. This study demonstrates the principles of de novo establishment of PRC2-mediated repressive domains across the genome.

Chromatin structure undergoes global and local reorganization during murine dendritic cell development and activation.

Classical dendritic cells (cDCs) are essential for immune responses and differentiate from hematopoietic stem cells via intermediate progenitors, such as monocyte-DC progenitors (MDPs) and common DC progenitors (CDPs). Upon infection, cDCs are activated and rapidly express host defense-related genes, such as those encoding cytokines and chemokines. Chromatin structures, including nuclear compartments and topologically associating domains (TADs), have been implicated in gene regulation. However, the extent and dynamics of their reorganization during cDC development and activation remain unknown. In this study, we comprehensively determined higher-order chromatin structures by Hi-C in DC progenitors and cDC subpopulations. During cDC differentiation, chromatin activation was initially induced at the MDP stage. Subsequently, a shift from inactive to active nuclear compartments occurred at the cDC gene loci in CDPs, which was followed by increased intra-TAD interactions and loop formation. Mechanistically, the transcription factor IRF8, indispensable for cDC differentiation, mediated chromatin activation and changes into the active compartments in DC progenitors, thereby possibly leading to cDC-specific gene induction. Using an infection model, we found that the chromatin structures of host defense-related gene loci were preestablished in unstimulated cDCs, indicating that the formation of higher-order chromatin structures prior to infection may contribute to the rapid responses to pathogens. Overall, these results suggest that chromatin structure reorganization is closely related to the establishment of cDC-specific gene expression and immune functions. This study advances the fundamental understanding of chromatin reorganization in cDC differentiation and activation.

Alu insertion variants alter gene transcript levels.

Alu are high copy number interspersed repeats that have accumulated near genes during primate and human evolution. They are a pervasive source of structural variation in modern humans. Impacts that Alu insertions may have on gene expression are not well understood, although some have been associated with expression quantitative trait loci (eQTLs). Here, we directly test regulatory effects of polymorphic Alu insertions in isolation of other variants on the same haplotype. To screen insertion variants for those with such effects, we used ectopic luciferase reporter assays and evaluated 110 Alu insertion variants, including more than 40 with a potential role in disease risk. We observed a continuum of effects with significant outliers that up- or down-regulate luciferase activity. Using a series of reporter constructs, which included genomic context surrounding the Alu, we can distinguish between instances in which the Alu disrupts another regulator and those in which the Alu introduces new regulatory sequence. We next focused on three polymorphic Alu loci associated with breast cancer that display significant effects in the reporter assay. We used CRISPR to modify the endogenous sequences, establishing cell lines varying in the Alu genotype. Our findings indicate that Alu genotype can alter expression of genes implicated in cancer risk, including PTHLH, RANBP9, and MYC These data show that commonly occurring polymorphic Alu elements can alter transcript levels and potentially contribute to disease risk.

Repetitive Elements: Different Subtypes Hint at Distinct Functions.

Lu et al. report that the association of different repeat types with distinct gene classes goes far beyond what has previously been shown and suggest that such relationship might be essential for gene function and regulation. As an example, they describe how long interspersed nuclear repeat (LINE1) transcripts are recruited together with associated genes to silent nuclear regions.

Close proximity to Igh is a contributing factor to AID-mediated translocations.

Class switch recombination (CSR) has the potential to generate genomic instability in B cells as activation-induced cytidine deaminase (AID), which mediates this process, is known to target many sites outside Igh. Nonetheless we do not fully understand what factors influence AID targeting genome-wide. Given that errors in CSR can lead to dangerous, oncogenic chromosomal translocations it is important to identify the elements that determine which genes are at risk of being "hit" and could be involved in aberrant rearrangements. Here we have investigated the influence of nuclear organization in determining "off-target" activity and the choice of fusion partners. Our studies indicate that the vast majority of known AID-mediated Igh translocation partners are found in chromosomal domains that contact this locus during class switching. Further, these interaction domains can be used to identify other genes that are hit by AID.

Identification of multi-loci hubs from 4C-seq demonstrates the functional importance of simultaneous interactions.

Use of low resolution single cell DNA FISH and population based high resolution chromosome conformation capture techniques have highlighted the importance of pairwise chromatin interactions in gene regulation. However, it is unlikely that associations involving regulatory elements act in isolation of other interacting partners that also influence their impact. Indeed, the influence of multi-loci interactions remains something of an enigma as beyond low-resolution DNA FISH we do not have the appropriate tools to analyze these. Here we present a method that uses standard 4C-seq data to identify multi-loci interactions from the same cell. We demonstrate the feasibility of our method using 4C-seq data sets that identify known pairwise and novel tri-loci interactions involving the Tcrb and Igk antigen receptor enhancers. We further show that the three Igk enhancers, MiEκ, 3'Eκ and Edκ, interact simultaneously in this super-enhancer cluster, which add to our previous findings showing that loss of one element decreases interactions between all three elements as well as reducing their transcriptional output. These findings underscore the functional importance of simultaneous interactions and provide new insight into the relationship between enhancer elements. Our method opens the door for studying multi-loci interactions and their impact on gene regulation in other biological settings.

4C-ker: A Method to Reproducibly Identify Genome-Wide Interactions Captured by 4C-Seq Experiments.

4C-Seq has proven to be a powerful technique to identify genome-wide interactions with a single locus of interest (or "bait") that can be important for gene regulation. However, analysis of 4C-Seq data is complicated by the many biases inherent to the technique. An important consideration when dealing with 4C-Seq data is the differences in resolution of signal across the genome that result from differences in 3D distance separation from the bait. This leads to the highest signal in the region immediately surrounding the bait and increasingly lower signals in far-cis and trans. Another important aspect of 4C-Seq experiments is the resolution, which is greatly influenced by the choice of restriction enzyme and the frequency at which it can cut the genome. Thus, it is important that a 4C-Seq analysis method is flexible enough to analyze data generated using different enzymes and to identify interactions across the entire genome. Current methods for 4C-Seq analysis only identify interactions in regions near the bait or in regions located in far-cis and trans, but no method comprehensively analyzes 4C signals of different length scales. In addition, some methods also fail in experiments where chromatin fragments are generated using frequent cutter restriction enzymes. Here, we describe 4C-ker, a Hidden-Markov Model based pipeline that identifies regions throughout the genome that interact with the 4C bait locus. In addition, we incorporate methods for the identification of differential interactions in multiple 4C-seq datasets collected from different genotypes or experimental conditions. Adaptive window sizes are used to correct for differences in signal coverage in near-bait regions, far-cis and trans chromosomes. Using several datasets, we demonstrate that 4C-ker outperforms all existing 4C-Seq pipelines in its ability to reproducibly identify interaction domains at all genomic ranges with different resolution enzymes.

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.

Structural perturbation of chromatin domains with multiple developmental regulators can severely impact gene regulation and development.

Chromatin domain boundaries delimited by CTCF motifs can restrict the range of enhancer action. However, disruption of domain structure often results in mild gene dysregulation and thus predicting the impact of boundary rearrangements on animal development remains challenging. Here, we tested whether structural perturbation of a chromatin domain with multiple developmental regulators can result in more acute gene dysregulation and severe developmental phenotypes. We targeted clusters of CTCF motifs in a domain of the mouse genome containing three FGF ligand genes-Fgf3, Fgf4, and Fgf15-that regulate several developmental processes. Deletion of the 23.9kb cluster that defines the centromeric boundary of this domain resulted in ectopic interactions of the FGF genes with enhancers located across the deleted boundary that are active in the developing brain. This caused strong induction of FGF expression and perinatal lethality with encephalocele and orofacial cleft phenotypes. Heterozygous boundary deletion was sufficient to cause these fully penetrant phenotypes, and strikingly, loss of a single CTCF motif within the cluster also recapitulated ectopic FGF expression and caused encephalocele. However, such phenotypic sensitivity to perturbation of domain structure did not extend to all CTCF clusters of this domain, nor to all developmental processes controlled by these three FGF genes-for example, the ability to undergo lineage specification in the blastocyst and pre-implantation development were not affected. By tracing the impact of different chromosomal rearrangements throughout mouse development, we start to uncover the determinants of phenotypic robustness and sensitivity to perturbation of chromatin boundaries. Our data show how small sequence variants at certain domain boundaries can have a surprisingly outsized effect and must be considered as potential sources of gene dysregulation during development and disease.

Loss of Ezh2 in the medial ganglionic eminence alters interneuron fate, cell morphology and gene expression profiles.

Introduction: Enhancer of zeste homolog 2 (Ezh2) is responsible for trimethylation of histone 3 at lysine 27 (H3K27me3), resulting in repression of gene expression. Here, we explore the role of Ezh2 in forebrain GABAergic interneuron development. Methods: We removed Ezh2 in the MGE by generating Nkx2-1Cre;Ezh2 conditional knockout mice. We then characterized changes in MGE-derived interneuron fate and electrophysiological properties in juvenile mice, as well as alterations in gene expression, chromatin accessibility and histone modifications in the MGE. Results: Loss of Ezh2 increases somatostatin-expressing (SST+) and decreases parvalbumin-expressing (PV+) interneurons in the forebrain. We observe fewer MGE-derived interneurons in the first postnatal week, indicating reduced interneuron production. Intrinsic electrophysiological properties in SST+ and PV+ interneurons are normal, but PV+ interneurons display increased axonal complexity in Ezh2 mutant mice. Single nuclei multiome analysis revealed differential gene expression patterns in the embryonic MGE that are predictive of these cell fate changes. Lastly, CUT&Tag analysis revealed that some genomic loci are particularly resistant or susceptible to shifts in H3K27me3 levels in the absence of Ezh2, indicating differential selectivity to epigenetic perturbation. Discussion: Thus, loss of Ezh2 in the MGE alters interneuron fate, morphology, and gene expression and regulation. These findings have important implications for both normal development and potentially in disease etiologies.

Med12 is essential for early mouse development and for canonical Wnt and Wnt/PCP signaling.

The Mediator complex is commonly seen as a molecular bridge that connects DNA-bound transcription factors to the RNA polymerase II (Pol II) machinery. It is a large complex of 30 subunits that is present in all eukaryotes. The Med12 subunit has been implicated not only in the regulation of Pol II activity, but also in the binding of transcription factors to the bulk of the Mediator complex. We targeted Med12 in mouse embryonic stem cells to investigate the in vivo function of this subunit. We report here the developmental defects of Med12 hypomorphic mutants that have a drastic reduction in Med12 protein levels. These mutants fail to develop beyond embryonic day 10 and have severe defects in neural tube closure, axis elongation, somitogenesis and heart formation. We show that in Med12 hypomorphic embryos, the Wnt/planar cell polarity pathway is disrupted and that canonical Wnt/beta-catenin signaling is impaired. In agreement with this, embryos that are incapable of Med12 expression failed to establish the anterior visceral endoderm or activate brachyury expression, and did not complete gastrulation.

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.

Enhancer-promoter interactions can bypass CTCF-mediated boundaries and contribute to phenotypic robustness.

How enhancers activate their distal target promoters remains incompletely understood. Here we dissect how CTCF-mediated loops facilitate and restrict such regulatory interactions. Using an allelic series of mouse mutants, we show that CTCF is neither required for the interaction of the Sox2 gene with distal enhancers, nor for its expression. Insertion of various combinations of CTCF motifs, between Sox2 and its distal enhancers, generated boundaries with varying degrees of insulation that directly correlated with reduced transcriptional output. However, in both epiblast and neural tissues, enhancer contacts and transcriptional induction could not be fully abolished, and insertions failed to disrupt implantation and neurogenesis. In contrast, Sox2 expression was undetectable in the anterior foregut of mutants carrying the strongest boundaries, and these animals fully phenocopied loss of SOX2 in this tissue. We propose that enhancer clusters with a high density of regulatory activity can better overcome physical barriers to maintain faithful gene expression and phenotypic robustness.

MLL3/MLL4 methyltransferase activities control early embryonic development and embryonic stem cell differentiation in a lineage-selective manner.

H3K4me1 methyltransferases MLL3 (KMT2C) and MLL4 (KMT2D) are critical for enhancer activation, cell differentiation and development. However, roles of MLL3/4 enzymatic activities and MLL3/4-mediated enhancer H3K4me1 in these processes remain unclear. Here we report that constitutive elimination of both MLL3 and MLL4 enzymatic activities prevents initiation of gastrulation and leads to early embryonic lethality in mice. However, selective elimination of MLL3/4 enzymatic activities in embryonic, but not extraembryonic, lineages leaves gastrulation largely intact. Consistent with this, embryonic stem cells (ESCs) lacking MLL3/4 enzymatic activities can differentiate toward the three embryonic germ layers but show aberrant differentiation to extraembryonic endoderm (ExEn) and trophectoderm. The failure in ExEn differentiation can be attributed to markedly reduced enhancer-binding of the lineage-determining transcription factor GATA6. Furthermore, we show that MLL3/4-catalyzed H3K4me1 is largely dispensable for enhancer activation during ESC differentiation. Together, our findings suggest a lineage-selective, but enhancer activation-independent, role of MLL3/4 methyltransferase activities in early embryonic development and ESC differentiation.

A Therapeutically Targetable NOTCH1-SIRT1-KAT7 Axis in T-cell Leukemia.

T-cell acute lymphoblastic leukemia (T-ALL) is a NOTCH1-driven disease in need of novel therapies. Here, we identify a NOTCH1-SIRT1-KAT7 link as a therapeutic vulnerability in T-ALL, in which the histone deacetylase SIRT1 is overexpressed downstream of a NOTCH1-bound enhancer. SIRT1 loss impaired leukemia generation, whereas SIRT1 overexpression accelerated leukemia and conferred resistance to NOTCH1 inhibition in a deacetylase-dependent manner. Moreover, pharmacologic or genetic inhibition of SIRT1 resulted in significant antileukemic effects. Global acetyl proteomics upon SIRT1 loss uncovered hyperacetylation of KAT7 and BRD1, subunits of a histone acetyltransferase complex targeting H4K12. Metabolic and gene-expression profiling revealed metabolic changes together with a transcriptional signature resembling KAT7 deletion. Consistently, SIRT1 loss resulted in reduced H4K12ac, and overexpression of a nonacetylatable KAT7-mutant partly rescued SIRT1 loss-induced proliferation defects. Overall, our results uncover therapeutic targets in T-ALL and reveal a circular feedback mechanism balancing deacetylase/acetyltransferase activation with potentially broad relevance in cancer. SIGNIFICANCE: We identify a T-ALL axis whereby NOTCH1 activates SIRT1 through an enhancer region, and SIRT1 deacetylates and activates KAT7. Targeting SIRT1 shows antileukemic effects, partly mediated by KAT7 inactivation. Our results reveal T-ALL therapeutic targets and uncover a rheostat mechanism between deacetylase/acetyltransferase activities with potentially broader cancer relevance. This article is highlighted in the In This Issue feature, p. 1.

Extensive co-binding and rapid redistribution of NANOG and GATA6 during emergence of divergent lineages.

Fate-determining transcription factors (TFs) can promote lineage-restricted transcriptional programs from common progenitor states. The inner cell mass (ICM) of mouse blastocysts co-expresses the TFs NANOG and GATA6, which drive the bifurcation of the ICM into either the epiblast (Epi) or the primitive endoderm (PrE), respectively. Here, we induce GATA6 in embryonic stem cells-that also express NANOG-to characterize how a state of co-expression of opposing TFs resolves into divergent lineages. Surprisingly, we find that GATA6 and NANOG co-bind at the vast majority of Epi and PrE enhancers, a phenomenon we also observe in blastocysts. The co-bound state is followed by eviction and repression of Epi TFs, and quick remodeling of chromatin and enhancer-promoter contacts thus establishing the PrE lineage while repressing the Epi fate. We propose that co-binding of GATA6 and NANOG at shared enhancers maintains ICM plasticity and promotes the rapid establishment of Epi- and PrE-specific transcriptional programs.

An epigenome atlas of neural progenitors within the embryonic mouse forebrain.

A comprehensive characterization of epigenomic organization in the embryonic mouse forebrain will enhance our understanding of neurodevelopment and provide insight into mechanisms of neurological disease. Here we collected single-cell chromatin accessibility profiles from four distinct neurogenic regions of the embryonic mouse forebrain using single nuclei ATAC-Seq (snATAC-Seq). We identified thousands of differentially accessible peaks, many restricted to distinct progenitor cell types or brain regions. We integrated snATAC-Seq and single cell transcriptome data to characterize changes of chromatin accessibility at enhancers and promoters with associated transcript abundance. Multi-modal integration of histone modifications (CUT&Tag and CUT&RUN), promoter-enhancer interactions (Capture-C) and high-order chromatin structure (Hi-C) extended these initial observations. This dataset reveals a diverse chromatin landscape with region-specific regulatory mechanisms and genomic interactions in distinct neurogenic regions of the embryonic mouse brain and represents an extensive public resource of a 'ground truth' epigenomic landscape at this critical stage of neurogenesis.

Cavß1 regulates T cell expansion and apoptosis independently of voltage-gated Ca2+ channel function.

TCR stimulation triggers Ca2+ signals that are critical for T cell function and immunity. Several pore-forming α and auxiliary ß subunits of voltage-gated Ca2+ channels (VGCC) were reported in T cells, but their mechanism of activation remains elusive and their contribution to Ca2+ signaling in T cells is controversial. We here identify CaVß1, encoded by Cacnb1, as a regulator of T cell function. Cacnb1 deletion enhances apoptosis and impairs the clonal expansion of T cells after lymphocytic choriomeningitis virus (LCMV) infection. By contrast, Cacnb1 is dispensable for T cell proliferation, cytokine production and Ca2+ signaling. Using patch clamp electrophysiology and Ca2+ recordings, we are unable to detect voltage-gated Ca2+ currents or Ca2+ influx in human and mouse T cells upon depolarization with or without prior TCR stimulation. mRNAs of several VGCC α1 subunits are detectable in human (CaV3.3, CaV3.2) and mouse (CaV2.1) T cells, but they lack transcription of many 5' exons, likely resulting in N-terminally truncated and non-functional proteins. Our findings demonstrate that although CaVß1 regulates T cell function, these effects are independent of VGCC channel activity.

A Tumor Suppressor Enhancer of PTEN in T-cell development and leukemia.

Long-range oncogenic enhancers play an important role in cancer. Yet, whether similar regulation of tumor suppressor genes is relevant remains unclear. Loss of expression of PTEN is associated with the pathogenesis of various cancers, including T-cell leukemia (T-ALL). Here, we identify a highly conserved distal enhancer (PE) that interacts with the PTEN promoter in multiple hematopoietic populations, including T-cells, and acts as a hub of relevant transcription factors in T-ALL. Consistently, loss of PE leads to reduced PTEN levels in T-ALL cells. Moreover, PE-null mice show reduced Pten levels in thymocytes and accelerated development of NOTCH1-induced T-ALL. Furthermore, secondary loss of PE in established leukemias leads to accelerated progression and a gene expression signature driven by Pten loss. Finally, we uncovered recurrent deletions encompassing PE in T-ALL, which are associated with decreased PTEN levels. Altogether, our results identify PE as the first long-range tumor suppressor enhancer directly implicated in cancer.